Open access journal: http://periodicos.uefs.br/ojs/index.php/sociobiology ISSN: 0361-6525 DOI: 10.13102/sociobiology.v68i1.5911Sociobiology 68(1): e5911 (June, 2021) Introduction The subterranean termite Globitermes sulphureus is commonly found in Malaysia, Singapore, Thailand, and Vietnam (Ahmad, 1965; Bordereau et al., 1997; Kuswanto et al., 2015; Lee et al., 2007; Ngee & Lee, 2002). This termite belongs to higher group termites which possess only bacteria and archaea in their gut (Bujang et al., 2014). As a wood feeder termite, this species has been reported to infest premises’ wood structures (Ab Majid & Ahmad, 2009; Neoh et al., 2011). Moreover, it was also reported as the primary pest in agricultural sectors such as coconut and oil palm plantations (Lee et al., 2003). G. sulphureus is recognized as a pest of significant economic importance in Southeast Asia (Rust & Su, 2012). Abstract The subterranean higher termite Globitermes sulphureus (Blattodea: Termitidae) is a peridomestic forager and regarded as a significant pest in Southeast Asia. In this study, populations of G. sulphureus from the USM main campus area were investigated based on partial sequences of the mitochondrial COII gene. The genetic diversity was determined using DnaSP v5 software, while the phylogenetic relationship was defined using Neighbor-joining (NJ) and maximum likelihood (ML) methods using Molecular Evolutionary Genetics Analysis (MEGA 7) software. A total of 2 haplotypes were detected among 5 sample sequences distinguished through two variable sites. Also, both phylogenetic trees gave similar topology and supporting the results from haplotype diversity. Based on the haplotype diversity and molecular phylogeny, it is proposed that geographic isolation and lack of human activities have contributed to the neutral genetic diversity of G. sulphureus. Sociobiology An international journal on social insects Nurul Akmar Hussin, Abdul Hafiz Ab Majid Article History Edited by Qiuying Huang, HZAU, China Received 09 October 2020 Initial acceptance 13 February 2021 Final acceptance 11 March 2021 Publication date 11 June 2021 Keywords Genetic diversity, phylogenetic, termite, Globitermes sulphureus, COII. Corresponding author Abdul Hafiz Ab Majid Household & Structural Urban Entomology Laboratory Vector Control Research Unit, School of Biological Sciences Universiti Sains Malaysia Penang, 11800 Minden, Malaysia. E-Mail: abdhafiz@usm.my Most of the currently available data for G. sulphureus distribution in Malaysia is based on samples gathered from home infestations, disturbed forest, forest reserve, and rural areas (Ab Majid & Ahmad, 2009; Aiman Hanis et al., 2014; Khizam & Ab Majid, 2019; Lee et al., 2004). However, its invasion biology, geographical distribution, and the pattern of introducing this species are still less understood. The phylogenetic relationships among populations and the impact of introductions to new habitats on genetic diversity and colony structure of G. sulphureus are relatively unknown (Khizam & Ab Majid, 2021). Previous researches focused on the genetic diversity and phylogenetics of the lower termites, especially in the genus Reticulitermes and Coptotermes, and dry wood termites genus Incisitermes, Cryptotermes, Kalotermes, Neotermes, Household & Structural Urban Entomology Laboratory, Vector Control Research Unit, School of Biological Sciences, Universiti Sains Malaysia, Penang, Malaysia RESEARCH ARTICLE - TERMITES Genetic diversity and phylogenetic relationship of higher termite Globitermes sulphureus (Haviland) (Blattodea: Termitidae) Nurul A Hussin, Abdul HA Majid – Genetic diversity and phylogenetic of Globitermes sulphureus2 and Mastotermes isolated from multiple locations (Austin et al., 2012; Husseneder et al., 2012; Leniaud et al., 2010; Pinzon & Houseman, 2009; Szalanski et al., 2008; Thompson et al., 2000; Yeap et al., 2007). The genetic diversity and phylogenetics of higher group termites have also been widely studied, such as genus Odontotermes, Microtermes, Microcerotermes, Macrotermes, Labiotermes, and Nasutitermes (Dupont et al., 2009; Ab Majid et al., 2018; Murthy et al., 2015; Ohkuma et al., 2004; Singla et al., 2016; Singla et al., 2015). Most of the studies mentioned above use COII gene sequences as molecular markers in the genetic diversity and phylogenetic relationship. The COII is a subunit of cytochrome c oxidase in the mitochondria (Frati et al., 1997). COII is considered the fastest evolving gene compared to 12S and 16S genes (Yeap et al., 2007). The rapid evolution of these genes has been proven helpful for deducing phylogenetic relationships between closely related insect species due to the relatively high degree of variation at the 3’ end of this gene (Aly et al., 2012; Singla et al., 2016). G. sulphureus is a peridomestic forager and mound- building termite (Lee et al., 2003). Its mound is easily identified based on a dome shape and a dark brown color (Ahmad, 1965). G. sulphureus is easily identified since the soldiers possess a bright-yellow colored body (Ab Majid & Ahmad, 2011; Hussin & Ab Majid, 2017; Hussin et al., 2018; Khizam & Ab Majid, 2019). The goal of this study is to determine the genetic diversity and phylogenetic relationship among local populations of G. sulphureus in Universiti Sains Malaysia (USM) main campus, Penang by using partial sequences of the mitochondrial COII gene. There may also be differences in the termite genome from several mounds of G. sulphureus in response to local ecological conditions. Materials and Methods Termites Collection Termite specimens were collected from the termite mounds. There are five accessible termite mounds identified in Universiti Sains Malaysia (Fig 1). The five colonies are colony 1: Durian Valley (DV), colony 2: Indah Kembara (IK), colony 3: Padang Minden (PM), colony 4: Tasik Harapan (TH), colony 5: Bakti Permai (BP) (Table 1). The termites were collected in a universal bottle containing 10 mL of 90% ethanol. Fig 1. Map of the USM main campus, Penang showing the sites of G. sulphureus populations being collected. The abbreviation details are listed in Table 1. The map was edited and retrieved from Meng et al., (2002). Table 1. Collection sites, isolated code, source, and GPS coordinate of each G. sulphureus in this study. Isolated Code Collection Site Source Lat_Lon DV Durian Valley Mound 5˚21’34.9” N 100˚18’19.3” E IK Indah Kembara Mound 5˚21’20.0” N 100˚17’35.2” E PM Padang Minden Mound 5˚21’34.6” N 100˚18’29.4” E TH Tasik Harapan Mound 5˚21’14.5” N 100˚18’01.1” E BP Bakti Permai Mound 5˚21’30.7” N 100˚18’04.9” E Sociobiology 68(1): e5911 (June, 2021) 3 DNA Extraction and PCR Amplification Ten worker termites from each collection site were rinsed in sterile distilled water and dried with a paper towel. The workers’ heads were cut off from the body using sterile dissecting scissors for DNA extraction. Genomic DNA extraction was performed using the DNeasy Blood & Tissues Kit (Qiagen, Germany) for COII gene PCR amplification. The COII gene was amplified using a pair of primer, Forward (5’-CATTGCACCCGCAATCATCC-3’) and Reverse (5’-GAATCTGTGGTTTGCTCCTCCGC-3’). The PCR reaction mixture (50 μL) contains 25 μL of 2x TopTaq Master Mix (Qiagen, Germany), producing a final concentration of 1.5 mM MgCl2, 1.25 units TopTaq DNA polymerase, 1× PCR buffer, and 200 µM of each dNTP, together with 1 μL (10 μM) of each primer, 5 μL (10x) CoralLoad Concentrate (as a substitute to loading dye), 10-100 ng of a bulk DNA template, and sterile distilled water. The PCR reaction profile comprises an initial denaturation at 94 °C for 2 min, followed by 40 cycles of denaturation at 94 °C for 45 seconds, primer annealing at 50 °C for 45 seconds, the first extension step at 72 °C for 60 seconds, and the final extension step at 72 °C for 5 min. The PCR product was purified with MEGAquick-spinTM Total Fragment DNA Purification Kit (iNtRON Biotechnology, South Korea). Sequencing and Analysis of COII Gene Sequences The purified PCR products were sequenced using the Sanger sequencing machine at First Base Laboratories Sdn. Bhd., Malaysia. Raw data obtained were extracted using FinchTV 1.4 (www.geospiza.com). Then, the forward and reverse sequences were aligned using T-Coffee (Notredame et al., 2000). The aligned sequences were edited manually to remove the low-quality bases before converting the files to FASTA format. The FASTA files were BLASTN search at the NCBI database (https://blast.ncbi.nlm.nih.gov) to compare with all COII gene sequences available in the GenBank. The species identified are based on ≥ 99% sequences similarity with the NCBI database. The partial COII sequences were submitted to the GenBank. Genetic Diversity and Phylogenetic Analysis All sequences were aligned automatically using the multiple alignment algorithm in ClustalX v2.1 (Larkin et al., 2007; Thompson et al., 1997) with default settings. Genetic characteristics such as haplotype diversity, nucleotide diversity, and the total number of mutations were calculated using DnaSP v5 (Librado & Rozas, 2009). The pairwise genetic distance was calculated using the p-distance method in MEGA 7. The phylogenetic relationship was carried out using the Molecular Evolutionary Genetics Analysis (MEGA) version 7 (Kumar et al., 2016). The phylogenetic trees were constructed using the Neighbor-Joining (NJ) (Saitou & Nei, 1987) and Maximum Likelihood (ML) method based on p-distance (Nei & Kumar, 2000) and Hasegawa- Kishino-Yano model (Hasegawa et al., 1985), respectively. The C. gestroi (lower termite) with an accession number GU931692.1 was used as an outgroup. Bootstrap analysis with 1000 resamplings was used to establish the NJ and ML trees (Fellsenstein, 1985). Results Five COII gene sequences were used with an average amplicon size of 420 base pairs after editing and removing the low-quality bases (bp). BLASTN results confirm the termite species is G. sulphureus, and the accession number for submitted sequences is written in Table 2. The edited sequences are used for genetic diversity and phylogenetic analysis. Isolated code Collection sites Haplotype Hd Pi Accession no. TH Tasik Harapan 1 MF997559.1 IK Indah Kembara 1 MF997563.1 BP Bakti Permai 1 0.400 0.00192 MF997545.1 PM Padang Minden 1 MF997577.1 DV Durian Valley 2 MF997568.1 Nucleotide Analysis The average base frequencies for five COII nucleotide sequences are A = 35.5%, T = 22.9%, G = 16.1%, and C = 25.5%. The total content of A+T is 58.4%, much higher than Table 2. Haplotype, haplotype diversity (Hd), nucleotide diversity (Pi), and accession G. sulphureus from five sample locations. C+G = 41.6%. The nucleotide diversity, Pi, is 0.00192 (Table 2). Two polymorphic or singleton sites are observed in the COII nucleotide sequence of DV samples at positions 252 and 255 (Table 3). The pairwise genetic distance between partial COII gene sequences is 0.000 to 0.005 (Table 4). https://blast.ncbi.nlm.nih.gov Nurul A Hussin, Abdul HA Majid – Genetic diversity and phylogenetic of Globitermes sulphureus4 Haplotype Analysis We observed two haplotypes, with haplotype 1 consisting of four colonies and haplotype 2 consisting of only one colony (Table 2), distinguished by two variable sites (Table 3). The haplotype diversity (Hd) of the five samples is 0.400 (Table 2). The relationship between haplotypes is further confirmed upon employing NJ and ML methods as implemented in MEGA 7 (Fig 2 and Fig 3). Phylogenetic Relationship Inferred From COII Genes The phylogenetic relationship of G. sulphureus from five locations is analyzed using two approaches; distance matrix: NJ (Fig 2) and character-based method: ML (Fig 3). Both trees show a monophyletic group among five samples but with three separate clades corresponding to the outgroup species, C. gestroi CG003TW. The first clade comprises haplotype 1, while the second clade comprises haplotype 2. The third clade shows the separation of Termitidae from the Rhinotermitidae family. The trees are also supported by pairwise genetic distance values (Table 4). Colony DV has a pairwise genetic distance value of 0.005 between other colonies, causing it to be located in the second clade. Meanwhile, all colonies in haplotype 1 have 0.000 value of the pairwise genetic distance between them. Nucleotide position Samples 248 249 250 251 252 253 254 255 256 257 258 TH C A A T C C G G T T A IK . . . . . . . . . . . BP . . . . . . . . . . . PM . . . . . . . . . . . DV . . . . T . . A . . . Table 3. Aligned haplotype sequences showing polymorphic sites (highlighted in green color). 1 2 3 4 5 1. TH 2. IK 0.000 3. BP 0.000 0.000 4. PM 0.000 0.000 0.000 5. DV 0.005 0.005 0.005 0.005 Table 4. Pairwise genetic distance (p-distance) in species under study. Fig 2. Neighbor-Joining tree built on the partial COII gene sequences of G. sulphureus from five different mounds. C. gestroi is used as an outgroup. The number at the branch refers to the bootstrap value. The evolutionary distances are computed using the p-distance method. Bar = 0.02 substitutions per nucleotide position. Sociobiology 68(1): e5911 (June, 2021) 5 Discussion There are five accessible mounds of G. sulphureus identified on the USM main campus. Colony DV is located at Durian Valley, known as natural forest/habitat in the USM main campus area. This location is protected to preserve flora and fauna there (Meng et al., 2002). Colony PM is located at the bamboo trees near the Minden Field, where students actively play football. Colony BP is located closer to the hostel’s area, while colony IK is near the main roadside. Lastly, colony TH is located between hostels and the lake. Globitermus sulphureus is regarded as a peridomestic forager and previously found in a disturbed forest area rather than a natural forest area (Aiman Hanis et al., 2014). A disturbed forest area is a place that had been cleared and developed for eco-tourism activities where various wooden structures and facilities were built. G. sulphureus is also commonly found in rural and urban areas (Lee et al., 2007; Lee et al., 2003; Ab Majid & Ahmad, 2009). Urbanization of plantation areas attracted this termite, causing infestation at the door and window frames of houses (Ngee & Lee, 2002). Therefore, the encounter of this species in the USM area is not surprising. Partial COII gene sequences of G. sulphureus isolated from five colonies are used for genetic diversity and phylogenetic analysis. The composition of different nucleotides in COII gene sequences from the colonies is calculated. From the results, the COII gene sequences show a high percentage of A+T (58.4%) than G+C (41.6%). Adenine and thymine bias in nucleotide sequences is consistent with the data on COII mitochondrial genes of other higher group termites species (genus Microtermes, Microcerotermes, and Odontotermes) and lower termites (genus Coptotermes) (Singla et al., 2016; Yeap et al., 2007). High A+T content is typical in insect mitochondrial DNA (Keller et al., 2007). However, the adenine-thymine percentage recorded in this study is lower than reported by those studies (59.87 - 62.0%). This result might occur because partial sequences of the COII gene are used in this study instead of full sequences (675 - 680 bp). The pairwise genetic distance for this study ranges from 0.000 to 0.005. This value is lower than reported previously regarding the genetic distance comparison between different genus of Odontotermes (ranged from 0.025 to 0.072) and genus Microtermes (ranged from 0.028 to 0.229) (Singla et al., 2016), and higher than the pairwise genetic distance between Macrotermes carbonarius populations (0.003) (Ab Majid et al., 2018). This result suggests that populations within a species have smaller genetic distances than populations within a genus because intraspecific populations are more closely related and have a recent common ancestor. From the haplotype analysis, two haplotypes are formed with DV colony solely in haplotype 2. The DV colony has two variable sites causing the pairwise genetic distance value to be 0.005 between other colonies. This result might occur due to the biased nature of natural forest zones (DV) compared to urban zones (IK, PM, BP, and TH). In a previous study, Szalanski et al. (2008) demonstrated that the genus Reticulitermes isolated from Lake Wedington (National Forest) has a higher frequency of rare haplotypes than Reticulitermes isolated from urban areas. Geographic Fig 3. Maximum likelihood tree inferred from the partial COII gene sequences of G. sulphureus from five different mounds. C. gestroi is used as an outgroup. The number at the branch refers to the bootstrap value. The evolutionary history is computed using the Hasegawa-Kishino-Yano model. Bar = 0.02 substitutions per nucleotide position. Nurul A Hussin, Abdul HA Majid – Genetic diversity and phylogenetic of Globitermes sulphureus6 isolation and the lack of human involvement in forested regions contributed to the rare haplotype of Reticulitermes species. Besides, demographic fluctuations and natural selections can affect any particular species’ neutral genetic diversity (Ellegren & Galtier, 2016). Two variable sites are detected in this study. This is contrary to other higher group termite (M. carbonarius) isolated from USM main campus, which has only one variable site (Ab Majid et al., 2018). Both variable sites are silent mutations since changes in the nucleotide bases do not affect amino acid sequences. Silent mutations in this study can be regarded as synonymous mutations since the nucleotide bases only change at the codon’s third position. Synonymous mutations usually occur by neutral selection and later become fixed (Fouks & Lattorff, 2016). It is long thought to be without phenotypic consequences but is currently recognized as critical in shaping gene expression, protein folding, cellular function, and the organism’s fitness (Plotkin & Kudla, 2011; Zwart et al., 2018). However, the consequences of synonymous mutations in termites fitness remained understood. The phylogenetic trees supported genetic diversity results where two clusters are formed regarding the haplotypes. Both trees (NJ and ML) show monophyletic of all five colonies. Bootstrap values (> 80%) show the significance of both trees. In the phylogenetic tree, closely related organisms are group together according to the order, family, subfamily, genus, and species (Ab Majid et al., 2018; Bourguignon et al., 2014; Singla et al., 2015; Yeap et al., 2007). Since there is no parsimony-informative site detected in this study, parsimony analysis is excluded. More sampling sites and analysis are needed to confirm the geographic distribution and the native and invasive termite G. sulphureus in USM main campus and Malaysia. Genotype analysis may explain the breeding pattern at the microsatellite level and reveal the nature of haplotype variability within G. sulphureus. In conclusion, this study demonstrated the COII gene’s ability to differentiate between G. sulphureus populations from a few different USM main campus locations. The genetic diversity analysis shows nucleotide divergence between isolated populations. Phylogenetic analysis supports the haplotype relationship of G. sulphureus. However, geographical area influences on the species’ genetic diversity require more sampling sites and further analysis such as microsatellite genotyping. Acknowledgement This research is supported by Research University Grant (RUi) (1001 / PBIOLOGI / 8011104) funded by Universiti Sains Malaysia . Author Contributions NAH: conceptualization, methodology, investigation, data curation, visualization, formal analysis, writing-original draft. AHAM: supervision, conceptualization, methodology, visualization, validation, project administration, resources, funding acquisition, writing -review & editing Declaration – Conflict of interest The authors report no conflicts of interest. The authors alone are responsible for the content and writing of the paper. Data availability statement The data that support the findings of this study are openly available in [NCBI] at [https://www.ncbi.nlm. nih.gov/], reference number MF997559.1, MF997563.1, MF997545.1, MF997577.1, MF997568.1]. References Ab Majid, A.H. & Ahmad, A.H. (2009). The status of subterranean termite infestation in Penang, Seberang Perai and Kedah, Malaysia. Malaysia Application Biology, 38: 37-48. Ab Majid, A.H. & Ahmad, A.H. (2011). Foraging population, territory and control of Globitermes sulphureus (Isoptera: Termitidae) with Fipronil in Penang, Malaysia. Malaysian Applied Biology, 40: 61-65. Ab Majid, A.H., Shen, E.Y., Heng, C.Y. & Foong, L.C. (2018). Genetic variation, diversity and molecular phylogenetic of higher group termite Macrotermes carbonarious Hagen (Blattodea: Termitidae). Malaysian Applied Biology, 47: 97-104. Ahmad, M. (1965). Termites (Isoptera) of Thailand. American Museum of Natural History, 131: 1-114. Aiman Hanis, J., Abu Hassan, A., Nurita, A.T. & Che Salmah, M.R. (2014). Community structure of termites in a hill dipterocarp forest of Belum- Temengor Forest Complex, Malaysia: Emergence of pest species. Raffles Bulletin of Zoology, 62: 3-11. Aly, S.M., Wen, J., Wang, X. & Cai, J. (2012). Cytochrome oxidase II gene “short fragment” applicability in identification of forensically important insects. Romanian Journal of Legal Medicine, 20: 231-236. doi: 10.4323/rjlm.2012.231 Austin, J.W., Allen, L.S., Solorzano, C., Magnus, R. & Scheffrahn, R.H. (2012). Mitochondrial DNA genetic diversity of the drywood termites Incisitermes minor and I. snyderi (Isoptera: Kalotermitidae). Florida Entomologist, 95: 75-81. doi: 10.1653/024.095.0112 Bordereau, C., Robert, A., Van Tuyen, V. & Peppuy, A. (1997). Suicidal defensive behaviour by frontal gland dehiscence in Globitermes sulphureus Haviland soldiers (Isoptera). Insectes Sociaux, 44: 289-296. doi: 10.1007/s000400050049 Bourguignon, T., Lo, N., Cameron, S.L., Sobotn, J., Hayashi, Y., Shigenobu, S., Watanabe, D., Roisin, Y., Miura, T. & Evans, T.A. (2014). The evolutionary history of termites as inferred from 66 mitochondrial genomes. Molecular Biology and Evolution, 32: 406-421. doi: 10.1093/molbev/msu308 Sociobiology 68(1): e5911 (June, 2021) 7 Bujang, N., Harrison, N. & Su, N. (2014). A phylogenetic study of endo-beta-1,4-glucanase in higher termites. Insectes Sociaux, 61: 29-40. doi: 10.1007/s00040-013-0321-7 Dupont, L., Roy, V., Bakkali, A. & Harry, M. (2009). Genetic variability of the soil-feeding termite Labiotermes labralis (Termitidae, Nasutitermitinae) in the Amazonian primary forest and remnant patches. Insect Conservation and Diversity, 2: 53-61. doi: 10.1111/j.1752-4598.2008.00040.x Ellegren, H. & Galtier, N. (2016). Determinants of genetic diversity. Nature Publishing Group, 17: 422-433. doi: 10.1038/ nrg.2016.58 Fellsenstein, J. (1985). Confidence limits on phylogenies : An approach using the bootstrap. Evolution, 39: 783-791. doi: 10.3389/fimmu.2015.00048 Fouks, B. & Lattorff, H.M.G. (2016). Contrasting evolutionary rates between social and parasitic bumblebees for three social effect genes. Frontiers in Ecology and Evolution, 4: 1-9. doi: 10.3389/fevo.2016.00064 Frati, F., Simon, C., Sullivan, J. & Swofford, D.L. (1997). Evolution of the mitochondrial cytochrome oxidase II gene in Collembola. Journal of Molecular Evolution, 44: 145-158. doi: 10.1007/PL00006131 Hasegawa, M., Kishino, H. & Yano, T. (1985). Dating of the human-ape splitting by a molecular clock of mitochondrial DNA. Journal of Molecular Evolution, 22: 160-174. doi: 10.1007/BF02101694 Husseneder, C., Simms, D.M., Delatte, J.R., Wang,C., Grace, J.K. & Vargo, E.L. (2012). Genetic diversity and colony breeding structure in native and introduced ranges of the Formosan subterranean termite, Coptotermes formosanus. Biological Invasions, 14: 419-437. doi: 10.1007/s10530-011- 0087-7 Hussin, N.A. & Ab Majid, A.H. (2017). Inter and intra termites colonies comparisons of gut microbial diversity from worker and soldier caste of Globitermes sulphureus (Blattodea: Termitidae) using 16S rRNA gene. Malaysian Journal of Microbiology, 13: 228-234. Hussin, N.A., Zarkasi, K.Z. & Ab Majid, A.H. (2018). Characterization of gut bacterial community associated with worker and soldier castes of Globitermes sulphureus Haviland (Blattodea: Termitidae) using 16S rRNA metagenomic. Journal of Asia-Pacific Entomology, 21: 1268-1274. doi: 10.1016/J.ASPEN.2018.10.002 Keller, I., Bensasson, D. & Nichols, R.A. (2007). Transition- transversion bias is not universal: A counter example from grasshopper pseudogenes. PLoS Genetics, 3: 0185-0191. doi: 10.1371/journal.pgen.0030022 Khizam, N.A.N. & Ab Majid, A.H. (2019). Development and annotation of species-specific microsatellite markers from transcriptome sequencing for a higher group termite, Globitermes sulphureus Haviland (Blattodea: Termitidae). Meta Gene, 20: 1-6. doi: 10.1016/j.mgene.2019.100568 Khizam, N.A.N. & Ab Majid, A.H. (2021). Population genetic structure and breeding pattern of higher group termite Globitermes sulphureus Haviland (Blattodea: Termitidae). Sociobiology, 68: e5772. doi: 10.13102/sociobiology.v68i1.5772 Kumar, S., Stecher, G. & Tamura, K. (2016). MEGA7: Molecular evolutionary genetics analysis version 7.0 for bigger datasets. Molecular Biology and Evolution, 33: 1870- 1874. doi: 10.1093/molbev/msw054 Kuswanto, E., Ahmad, I. & Dungani, R. (2015). Threat of subterranean termites attack in the Asian countries and their control: A review. Asian Journal of Applied Sciences, 8: 227- 239. doi: 10.3923/ajaps.2015.227.239 Larkin, M.A., Blackshields, G., Brown, N.P., Chenna, R., Mcgettigan, P.A., McWilliam, H., Valentin, F., Wallace, I.M., Wilm, A., Lopez, R., Thompson,J.D., Gibson, T.J. & Higgins, D. G. (2007). Clustal W and Clustal X version 2.0. Bioinformatics, 23: 2947-2948. doi: 10.1093/bioinformatics/ btm404 Lee, C.Y., Ngee, P.S., Lee, L.C. & Na, J.P.S. (2004). Survey of termite diversity in Pantai Acheh Forest Reserve, Penang Island, Malaysia. Jurnal Biosains, 15: 91-99. Lee, C.Y., Yap, J., Ngee, P.S. & Jaal, Z. (2003). Foraging colonies of a higher mound-building subterranean termite, Globitermes sulphureus (Haviland) in Malaysia. Japanese Journal of Environmental Entomology and Zoology, 14: 105-112. Lee, C.Y., Vongkaluang, C. & Lenz, M. (2007). Challenges to subterranean termite management of multi-genera faunas in Southeast Asia and Australia. Sociobiology, 50: 213-221. Leniaud, L., Dedeine, F., Pichon, A., Dupont, S. & Bagnères, A.G. (2010). Geographical distribution, genetic diversity and social organization of a new European termite, Reticulitermes urbis (Isoptera: Rhinotermitidae). Biological Invasions, 12: 1389-1402. doi: 10.1007/s10530-009-9555-8 Meng, L.L., Badarulzaman, N., Mui, L.Y., Awang, H. & Ta, T.L. (2002). The university in the garden, policies & guidelines (Vol. 1). Murthy, S., Rajeshwari, K. & Jalali, T. (2015). Genetic diversity among Indian termites based on mitochondrial 12S rRNA gene. European Journal of Zoological Research, 4: 1-6. Nei, M. & Kumar, S. (2000). Molecular evolution and phylogenetics. Oxford University Press, 333 p Neoh, K.B., Jalaludin, N.A. & Lee, C.Y. (2011). Elimination of field colonies of a mound-building termite Globitermes sulphureus (Isoptera: Termitidae) by bistrifluron bait. Journal of Economic Entomology, 104: 607-613. doi: 10.1603/ EC10161 Nurul A Hussin, Abdul HA Majid – Genetic diversity and phylogenetic of Globitermes sulphureus8 Ngee, P.S. & Lee, C.Y. (2002). Colony characterization of a mound-building subterranean termite, Globitermes sulphureus (Isoptera: Termitidae) using modified single-mark recapture technique. Sociobiology, 40: 525-532. Notredame, C., Higgins, D.G. & Heringa, J. (2000). T-Coffee: A novel method for fast and accurate multiple sequence alignment. Journal of Molecular Biology, 302: 205-217. doi: 10.1006/jmbi.2000.4042 Ohkuma, M., Yuzawa, H., Amornsak, W., Sornnuwat, Y., Takematsu, Y., Yamada, A., Vongkaluang, C., Sarnthoy, O., Kirtibutr, N., Noparatnaraporn, N., Kudo, T. & Inoue, T. (2004). Molecular phylogeny of Asian termites (Isoptera) of the families Termitidae and Rhinotermitidae based on mitochondrial COII sequences. Molecular Phylogenetics and Evolution, 31: 701-710. doi: 10.1016/j.ympev.2003.09.009 Pinzon, O.P. & Houseman, R.M. (2009). Species diversity and intraspecific genetic variation of Reticulitermes (Isoptera: Rhinotermitidae) subterranean termites in Woodland and Urban Environments of Missouri. Annals of the Entomological Society of America, 102: 868-880. doi: 10.1603/008.102.0513 Plotkin, J. B. & Kudla, G. (2011). Synonymous but not the same: the causes and consequences of codon bias. National Review of Genetics, 12: 32-42. doi: 10.1038/nrg2899. Synonymous Rust, M.K. & Su, N.Y. (2012). Managing social insects of urban importance. Annual Review of Entomology, 57: 355- 375. doi: 10.1146/annurev-ento-120710-100634 Saitou, N. & Nei, M. (1987). The neighbor-joining method: a new method for reconstructing phylogenetic trees. Molecular Biology and Evolution, 4: 406-425. doi: 10.1093/ oxfordjournals.molbev.a040454 Singla, M., Goyal, N., Sharma, R. & Singla, N. (2016). Reconstructing phylogenetic relationship among Indian termite species inferred from COII gene sequences ( lattodea: Isoptera : Termitidae ). Journal of Entomology and Zoology Studies, 4: 1-7. Singla, M., Goyal, N., Sobti, R.C. & Sharma, V.L. (2015). Estimating molecular phylogeny of some Indian termites combining partial COI sequences, Journal of Entomology and Zoology Studies, 3: 213-218. Szalanski, A.L., Austin, J.W. & McKern, J.A. (2008). Genetic diversity of Reticulitermes termites (Isoptera Rhinotermitidae) from Lake Wedington, Arkansas. Sociobiology, 52: 95-106. Thompson, G.J., Miller, L.R., Lenz, M. & Crozier, R. H. (2000). Phylogenetic analysis and trait evolution in Australian lineages, of drywood termites (Isoptera, Kalotermitidae). Molecular Phylogenetics and Evolution, 17: 419-429. doi: 10.1006/mpev.2000.0852 Thompson, J.D., Gibson, T.J., Plewniak, F., Jeanmougin, F. & Higgins, D.G. (1997). The CLUSTAL X windows interface: Flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucleic Acids Research, 25: 4876- 4882. doi: 10.1093/nar/25.24.4876 Yeap, B., Othman, A.S., Lee, V.S. & Lee, C. (2007). Genetic relationship between Coptotermes gestroi and Coptotermes vastator (Isoptera : Rhinotermitidae ). Journal of Economic Entomology, 100: 467-474. Zwart, M.P., Schenk, M.F., Hwang, S., Koopmanschap, B., de Lange, N., van de Pol, L., Nga, T.T.T., Szendro, I.G., Krug, J. & de Visser, J.A.G.M. (2018). Unraveling the causes of adaptive benefits of synonymous mutations in TEM-1 β-lactamase. Heredity, 121: 406-421. doi: 10.1038/s41437- 018-0104-z _Hlk67661769