alexa Diversity of Subterranean Termites in South India Based on COI Gene
ISSN: 2376-0214
Journal of Biodiversity, Bioprospecting and Development
Make the best use of Scientific Research and information from our 700+ peer reviewed, Open Access Journals that operates with the help of 50,000+ Editorial Board Members and esteemed reviewers and 1000+ Scientific associations in Medical, Clinical, Pharmaceutical, Engineering, Technology and Management Fields.
Meet Inspiring Speakers and Experts at our 3000+ Global Conferenceseries Events with over 600+ Conferences, 1200+ Symposiums and 1200+ Workshops on Medical, Pharma, Engineering, Science, Technology and Business
All submissions of the EM system will be redirected to Online Manuscript Submission System. Authors are requested to submit articles directly to Online Manuscript Submission System of respective journal.
  • Research Article   
  • J Biodivers Biopros Dev 2017, Vol 4(1): 161
  • DOI: 10.4172/2376-0214.1000161

Diversity of Subterranean Termites in South India Based on COI Gene

Srinivasa Murthy KS*, Yeda Lubna Banu and Ramakrishna P
Division of Molecular Entomology, ICAR-NBAIR, India
*Corresponding Author: Srinivasa Murthy KS, Division of Molecular Entomology, ICAR-NBAIR, Karnataka, 560024, India, Tel: 080 2351 1982, Email: [email protected]

Received Date: Jun 27, 2016 / Accepted Date: Apr 04, 2017 / Published Date: Apr 10, 2017


The diversity of subterranean termites collected from various locations in South India were characterised based on the COI gene using specific primers. Sequence analysis and divergence among the species was assessed. Genbank accession numbers were obtained for the different species. Phylogenetic tree based on neighbour- joining method was drawn on the basis of multiple sequence alignment, which revealed clustering of individuals according to the genera. Among the species, Odontotermes longignathus was more prevalent than others. The utility of COI gene to study the systematics of termites, their evolution and relatedness that would have implication on their management is discussed.

Keywords: Subterranean termites; CO1 gene; Genbank; Phylogenetic tree; Odontotermes longignathus


Termites (Isoptera) represent up to 95% of soil insect biomass [1,2] show an elaborated morphology and complex behaviour. The family Termitidae is considered to be the most evolved group with about 85% of all known genera and nearly 70% known species Ohkuma et al. Termites are considered as the most abundant invertebrates that are classified in about 280 genera, and over 2800 species within 14 subfamilies [1,3-5]. In India about 300 species within seven families have been reported [6] The subterranean termites which are of economic importance Wang et al. [7] to agriculture are mostly found in temperate climates Groves et al. [8] Termites are often referred as “ecosystem engineers” [9,10] as they play a vital role in recycling of plant materials and wood, modifying and improving the soil condition and composition, and providing food for other animals [3,7,11] and are also considered as potent catalyst due to their role involved in converting lignocellulose into biofuels Manjula et al. [12] which is of industrial value.

Termites are the serious pests of agricultural and horticultural crops that mainly destroys the roots and above ground parts and feed on paper, wood and timber Murthy et al. [10]. They are found to considerably damage the artificial structures and commodities due to their large colony size, varying feeding preferences and nesting behaviour Wang et al. [7]. It has been estimated that worldwide the overall cost annually for the control of damage caused by termites is more than $20 billion NY [13].

The understanding of termites in various biological processes is extremely limited regardless of their importance in agriculture and hence it is very important to understand their biology and ecology, which greatly relies on accurate species identification Singla et al. [14]. Termite on systematics was exclusively based on taxonomical identification based morphological characters of individuals belonging to various castes (e.g., soldiers or workers) Kambhampati and Eggleton [3]. Nevertheless, the caste differentiation, eusocial behaviour, varying physiological functions and crypto-biotic structure had contributed to ambiguities in their morphological identification [7,10] and species diagnosis has become a challenging task Kirton et al. [15].

The use of molecular methods which are fast and reliable, complementary to the morphological identification [2,16], are helpful in estimating evolutionary relatedness between the species Singla et al. [14]. Studies on mitochondrial genome sequences such as the ATrich region, 16S rDNA and cytochrome oxidase genes have shown an efficient alternative for species identification and phylogenetic studies [2,17-23]. The mitochondrial DNA is more abundant as the mitochondrial genes evolve more rapidly, than the nuclear genome Wang et al. [7], therefore at species level mitochondrial DNA is more suitable Masters et al. [24] and various other regions also can be sequenced [8,25,26].

Cytochrome c oxidase subunit I gene (CO1) is one of the three mitochondrial DNA (mtDNA) encoded subunits of respiratory complex IV. It is a key enzyme in aerobic metabolism. CO1 gene is the most conservative protein-coding gene in the mitochondrial genome. Since, the mutation rate in CO 1 gene is fast enough it can differentiate precisely the closely related termite species and assess their phylogeny, in understanding the evolutionary relationships, the gene has been extensively used to understand the diversity and genetic relatedness.

In the present study, we characterized twelve species of termites of family termitidae obtained from various locations in South India based on the mitochondrial CO 1 gene and studied their phylogenetic relation.

Materials And Methods

Collection and identification of termite samples

Termite specimens collected from different locations in South India were preserved in absolute alcohol and stored at -80°C at the Division of Molecular Entomology, NBAIR-ICAR Bangalore, India. Taxonomical identification of these specimens was done at the Division of Entomology Indian Agricultural Research Institute, New Delhi, Institute of Wood Science Technology Bangalore, and Centre for Insect taxonomy, University of Agricultural Sciences, Bangalore.

Isolation of genomic DNA

Isolation of genomic DNA was carried out using Qiagen DNeasy blood tissue kit following the established protocols. Termite samples were stored in absolute alcohol, were first washed well with distilled water followed by dissection of the head region. The dissected head was dried and collected in 1.5ml eppendorf tube to which 180 μl of ATL buffer was added and homogenised using micro pestle. 20 μl of proteinase K was added to it and mixed methodically by vortexing and then incubated in water bath at 56°C overnight. 100 μl of AL buffer was added and kept further for incubation at 56°C for 10 min, after incubation 100 μl of 100% ethanol was added to the eppendorf tube and vortexed. The solution was then transferred into the mini spin columns with silica member that binds the genomic DNA. The columns were then centrifuged at 8000 rpm for 5 minutes and the flowthrough collected in a tube was discarded. The column was then further given two washes with two different buffers present in kit i.e., AW1 and AW2 and centrifuged at 8000 rpm for 5 min simultaneously. These columns were then transferred to a new fresh 1.5 μl eppendorf tubes in order to elute out the DNA bound to the Silica member in the column. The elution was carried out by pipetting 100 μl of sterile water to the columns, centrifuged at 8000 rpm for 5 min. The DNA was then checked at 1% agarose gel and stored in 4°C until PCR was done.

Quantification of DNA by nanodrop spectrophotometer

The extracted DNA was quantified using nanodrop spectrophotometer and the amount of DNA present in the samples was identified. Since, all the organic compounds show a characteristic absorption, the nitrogenous bases in DNA show a strong absorption at a wavelength of 260 nm.

PCR amplification of CO1 gene

The DNA obtained was then amplified for a portion of mitochondrial CO1 gene fragment, using the universal primers CO-1 F 5’ GGTCAACAAATCATAAAGATATTGG 3’ and CO-1 R 5’ TAACTTCAGGCTGACCAAAAAATCA 3’, which were obtained from M/S Bioserve biotechnologies (India) Pvt Ltd. Each PCR reaction mixture of 25 μl consisted 2.5 μl of 10x PCR buffer with 15 mM MgCl2, 2.0 μl of dNTP’s mix, 1 μl of each forward and reverse primer, 1 μl of Taq Polymerase (1U/μl), 2.5 μl of template DNA and 15 μl of sterile water. The PCR was carried out in thermal-cycler (BioRad, USA) and the conditions set were initial denaturation at 950 C for 4 minutes, denaturation at 950 C for 30 seconds, annealing and extension at 500ºC and 720ºC for 1 min, respectively and final extension for 7 min at 720ºC. The PCR was carried out for a total of 34 cycles. Amplified DNA was then checked on 1% agarose gel using a DNA ladder of 250 bp and the gel was visualized in gel dock (Figure 1).


Figure 1: PCR amplified products of CO 1 gene of various termite populations, Lane M: 250bp DNA ladder, Lanes L1-L6: amplified products.


Dideoxy method or chain termination method was used for sequencing of DNA, wherein, the modified bases dideoxy bases i.e., ddNTP’s were used. The sequencing of amplified CO1 product was carried out at M/S. Eurofins Pvt Ltd, Bangalore. The sequence data was retrieved in the form of chromatograms which was then submitted to genbank for obtaining the accession numbers.

Sequence analysis and data interpretation

Chromatograms were edited in order to remove the ambiguous bases. These edited sequences were then aligned using Basic Local Alignment Search Tool (BLAST), with the sequences of same or related genera retrieved from the nucleotide database (PUBMED) of National Centre for Biological Information (NCBI). The CO1 nucleotide sequences of the termite species included in our present study were aligned and compared with the species obtained from PUBMED, using CLUSTAL W alignment Thompson et al. [27].

Phylogenetic analysis

Phylogenetic tree was constructed using neighbour- joining method and the evolutionary distances were computed using pdistance method with a boot strap consensus of 1000 replicates Tamura and Nie [28]. Constructed phylogenetic tree was visualised using tree viewer program.

Results and Discussion

The genomic DNA was collected from 12 populations of termites and the CO1 gene was characterized, using the universal primers CO-1 F 5’ GGTCAACAAATCATAAAGATATTGG 3’ and CO-1 R 5’ TAACTTCAGGCTGACCAAAAAATCA 3’. The amplified CO1 product was sequenced at M/S. Eurofins Pvt Ltd, Bangalore. The sequence data was submitted to NCBI genbank and the accession numbers were obtained for the populations of termites (Table 1). Place of collection Code Latitude Longitude Molecular identity GenBank Accession no.
1 Mysore Mys TE-2 12.296° 76.639° Hypotermesxenotermitis KT224387
2 Attur Attur TE-7 13.099° 77.568° Odontotermeslongignathus KT224388
3 PKP PKP TE-1 12.38° 75.05° Odontotermesescherichi KT224389
4 Sivaganga Sivaganga TE-1 9.843° 78.481° Odontotermeslongignathus  KT224390
5 Marat Marat TE-2 12.95° 77.70° Hypotermes  sp.  KT224391
6 Attur Attur TE-8 13.099° 77.568°  Odontotermeslongignathus KT254244
7 Sivaganga Sivaganga TE-2 9.843° 78.481° Hypotermesmakhamensis  KT274763
8 Sivaganga Sivaganga TE-3 9.843° 78.481° Hypotermesmakhamensis  KT274765 
9 Marat Marat TE-1 12.95° 77.70° Hypotermesxenotermitis  KT274764  
10 KAU,Thrissur BCCP59/2015 10.528° 76.214° Odontotermeslongignathus KT719274
11 Dasarahalli Dast-1 13.096° 77.835° Odontotermeslongignathus  KM015486
12 Ooty Ooty-TE-2 11.00° 78.00° Nasutitermesoctopilis  KM657478

Table 1: Place of collection of Termite species and their GenBank accession numbers.

Nucleotide analysis

The complete gene analysis of nucleotide sequence for each of the collected termite species showed a considerably high percentage of A+T base composition content, with an average composition of A+T=54.88% and G+C=45.11% of various species (Table 2). Identified species Sub-Family A+T (%) G+C (%)
1 Odontotermeslongignathus (DAST-1)  Macrotermitinae  57.28 42.72
2 Hypotermesxenotermitis (Mys TE-2)  Macrotermitinae  55.06 44.94
3 Odontotermeslongignathus (Attur TE-7)  Macrotermitinae  53.68 46.32
4 Odontotermeslongignathus (Sivaganga TE-1)  Macrotermitinae  53.93 46.07
5 Hypotermes  sp. (Marat TE-2)  Macrotermitinae  54.69 45.31
6 Odontotermeslongignathus (Attur TE-8)  Macrotermitinae  54.72 45.28
7 Hypotermesmakhamensis (Sivanganga TE-2)  Macrotermitinae  54.46 45.54
8 Hypotermesmakhamensis (Sivaganga TE-3)  Macrotermitinae  54.17 45.83
9 Odontotermeslongignathus (BCCP 59 Thrissur)  Macrotermitinae  54.52 45.48
10 Odontotermesescherichi  (PKP TE-1)  Macrotermitinae  52.78 47.22
11 Hypotermesxenotermitis (Marat TE-1)  Macrotermitinae  54.89 45.11
12 Nasutitermesoctopilis (Ooty-TE-2) Nasutitermitinae 58.41 41.59

Table 2: A+T and G+C content (%) for partial CO1 of the different termite species.

The variation in A+T (%) among the different populations was 4.63% and 5.68% with respect to G+C (%). The populations from Attur, Sivaganaga, Thrissur and Marthalli had the least variation in the base composition of A+T (%) (54.17 to 54.89) and (45.28 to 45.83), respectively. The estimated Transition/Transversion bias (R) is 1.65. Substitution pattern and rates were estimated under the Tamura-Nei [28] model. The nucleotide frequencies were A=32.20%, T/U=24.96%, C=27.27%, and G=15.57%. A tree topology was automatically computed for estimating the ML values W et al. [29]. The maximum Log likelihood for this computation was -1053.298, which involved analysis of 12 nucleotide sequences.

Codon positions included were 1st+2nd+3rd+Noncoding. All positions containing gaps and missing data were eliminated. There were a total of 206 positions in the final dataset. Evolutionary analyses were conducted in MEGA6 Szalanski et al. [30] software.

Divergence and percent identity

The divergence and percent identity was calculated using MEGA 6 software based on the sequence alignment. The overall average was found to be 0.06 for the total 48 nucleotide sequences. The number of base differences per site from between sequences is shown. Standard error estimate(s) are shown above the diagonal Tamura et al. [31]. There were a total of 190 positions in the final dataset.

Odontotermes longignathus population from Attur, Dasarahalli and Marthahalli had negligible divergence with the similarity co-efficient ranging between 0.01-0.03%, indicating the proximity in geographical location and topography. However, the population from Thrissur had low divergence of 0.03% Tamura et al. [31]. Therefore the occurrence of termite populations of species similarity across the geographical barriers did not vary, notwithstanding the habitat. Population from Mysore (Hypotermes xenotermitis) and Marthahalli indicated similar trend (Figure 2). Our observations are in broad conformity with the reports by earlier workers [6,10]. The greater occurrence of Odontotermes longignathus among the collections could be due to widest niche breadth [6,31-33].


Figure 2: Matrix showing the divergence and percent identity among the 12 species under study and 36 species of termites retrieved from the database (out-group taxa).

Phylogenetic analysis

A total of 48 strains were used for the phylogenetic analysis. The phylogenetic tree built using Hypotermes xenotermitis, Hypotermes makhamensis, Odontotermes longignathus, Nasutitermes octopilis, Odontotermes escherichi, Hypoteres sp were divided into two major clusters and three smaller clusters. The phylogenetic tree displayed in (Figure 3), shows that all the species belonging to Hypotermes xenotermitis were clubbed together with bootstrap score of 97, 53, 1 and 68. All the Hypotermes xenotermitis collected from Ratnagiri, India were grouped together with the bootstrap score of 97. The first major cluster was further divided into 6 sub-clusters. In all the 6 sub-clusters only the strains collected from Ratnagiri, India were clubbed together.


Figure 3: Phylogeny of some South Indian termites

The Hypotermes xenotermitis strains collected from Pune, India formed two clusters. In one cluster, most of the Hypotermes xenotermitis species collected from Pune, India were grouped together with bootstrap score of 53. In another cluster, only four sequences with the ID of KT879848, KT879846, KT879847 and KT879845 were grouped together with the bootstrap score of 63. These four gene sequences of Hypotermes xenotermitis species collected from Pune showed divergent from the other gene sequences collected from the same place. Apart from Hypotermes xenotermitis strains, four Hypotermes makhamensis strains were also collected from Pune and Sivaganga. All the four sequences along with the Hypotermes xenotermitis strain collected from Marat were grouped together and formed a cluster. The two strains (KT879848 and KT879846) collected from Pune were closely related and grouped together with bootstrap score of 87, whereas the two strains collected from Sivaganga were also closely related and grouped together to form a sub-cluster with bootstrap score of 98. The Hypotermes xenotermitis strain collected from Marat region was formed a single clade in the cluster. The Odontotermes longignathus strain collected from Thrissur and Attur were grouped together with a good bootstrap score of 70. This cluster was sub-divided into two small sub-clusters, where the Odontotermes longignathus strain collected from Attur were formed a one cluster with a bootstrap score of 71 and Odontotermes longignathus collected from Thrissur formed a single clade with a bootstrap score of 98. Apart from these, the Nasutitermes octopilis, Hypoteres sp, Odontotermes escherichi and only one strain of Odontotermes longignathus collected from DAST were formed a single clade which was appeared to be very divergent from other strains used to construct the phylogenetic tree. Overall, the phylogenetic tree analysis suggested that the Hypotermes xenotermitis strains collected from Ratnagiri were appeared to be closely related with the Hypotermes xenotermitis strains collected from Pune, India.

The occurrence of different species of termites could be attributed to varied ecology, food availability, geography and the dispersed pattern [20]. Knowledge of their diversity, distribution and understanding evolutionary relationships would provide insights to tailor effective management strategies (Table 3) [10,14,32,34-37]. Organism name
1 KM015486_Odontotermes_longignathus_strain_DAST-1
2 KM657478_Nasutitermes_octopilis_strain_Ooty-TE-2
3 KT224387_Hypotermes_xenotermitis_strain_Mys_TE-2
4 KT724955_Hypotermes_xenotermitis_strain_PPPU-5**
5 KT879849_Hypotermes_xenotermitis_strain_PPPU-20**
6 KT879849_Hypotermes_xenotermitis_strain_PPPU-15**
7 KT879849_Hypotermes_xenotermitis_strain_PPPU-13**
8 KT879849_Hypotermes_xenotermitis_strain_PPPU-12**
9 KT879849_Hypotermes_xenotermitis_strain_PPPU-11**
10 KT879849_Hypotermes_xenotermitis_strain_PPPU-10**
11 KT879849_Hypotermes_xenotermitis_strain_PPPU-9**
12 KT879849_Hypotermes_xenotermitis_strain_PPPU-6**
13 KT879849_Hypotermes_xenotermitis_strain_PPPU-5**
14 KT879849_Hypotermes_xenotermitis_strain_PPPU-19**
15 KT879849_Hypotermes_xenotermitis_strain_PPPU-18**
16 KT879849_Hypotermes_xenotermitis_strain_PPPU-17**
17 KT879849_Hypotermes_xenotermitis_strain_PPPU-16**
18 KT224388_Odontotermes_longignathus_strain_Attur_TE-7
19 KT224389_Odontotermes_escherichi_strain_PKP_TE-1
20 KT224391_Hypotermes_sp._Marat_TE-2
21 KT254244_Odontotermes_longignathus_strain_Attur_TE-8
22 KT274763_Hypotermes_makhamensis_strain_Sivanganga_TE-2
23 KT887717_Hypotermes_xenotermitis_strain_PPRA-20**
24 KT887717_Hypotermes_xenotermitis_strain_PPRA-19**
25 KT887717_Hypotermes_xenotermitis_strain_PPRA-18**
26 KT887717_Hypotermes_xenotermitis_strain_PPRA-17**
27 KT887717_Hypotermes_xenotermitis_strain_PPRA-16**
28 KT887717_Hypotermes_xenotermitis_strain_PPRA-15**
29 KT887717_Hypotermes_xenotermitis_strain_PPRA-14**
30 KT887717_Hypotermes_xenotermitis_strain_PPRA-13**
31 KT887717_Hypotermes_xenotermitis_strain_PPRA-12
32 KT887717_Hypotermes_xenotermitis_strain_PPRA-11**
33 KT274764_Hypotermes_xenotermitis_strain_Marat_TE-1
34 KT879843_Hypotermes_xenotermitis_strain_PPPU-14**
35 KT879843_Hypotermes_xenotermitis_strain_PPPU-3**
36 KT879843_Hypotermes_xenotermitis_strain_PPPU-2**
37 KT887707_Hypotermes_xenotermitis_strain_PPRA-10**
38 KT887707_Hypotermes_xenotermitis_strain_PPRA-9**
39 KT887707_Hypotermes_xenotermitis_strain_PPRA8**
40 KT887707_Hypotermes_xenotermitis_strain_PPRA-7**
41 KT887707_Hypotermes_xenotermitis_strain_PPRA-6**
42 KT887707_Hypotermes_xenotermitis_strain_PPRA-5**
43 KT887707_Hypotermes_xenotermitis_strain_PPRA-4**
44 KT887707_Hypotermes_xenotermitis_strain_PPRA-3**
45 KT887707_Hypotermes_xenotermitis_strain_PPRA-2**
46 KT887707_Hypotermes_xenotermitis_strain_PPRA-1**
47 KT274765_Hypotermes_makhamensis_strain_Sivaganga_TE-3
48 KT719274_Odontotermes_longignathus_strain_BCCP_59_Thrissur
Note:where ** marked indicate sequences from the NCBI

Table 3: Chronological order of termite species as indicated in the similarity matrix.


Citation: Murthy KSS, Banu Yl, Ramakrishna P (2017) Diversity of Subterranean Termites in South India Based on COI Gene. J Biodivers Biopros Dev 4: 161. Doi: 10.4172/2376-0214.1000161

Copyright: ©2015 Murthy KSS, et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Select your language of interest to view the total content in your interested language

Post Your Comment Citation
Share This Article
Relevant Topics
Article Usage
  • Total views: 1973
  • [From(publication date): 0-2017 - Nov 14, 2019]
  • Breakdown by view type
  • HTML page views: 1815
  • PDF downloads: 158
Share This Article