Genetic diversity and differentiation of populations of Chlorops oryzae (Diptera, Chloropidae)
BMC Ecology volume 20, Article number: 22 (2020)
Chlorops oryzae is an important pest of rice crops. There have been frequent outbreaks of this pest in recent years and it has become the main rice pest in some regions. To elucidate the molecular mechanism of frequent C. oryzae outbreaks, we estimated the genetic diversity and genetic differentiation of 20 geographical populations based on a dataset of ISSR markers and COI sequences.
ISSR data revealed a high level of genetic diversity among the 20 populations as measured by Shannon’s information index (I), Nei’s gene diversity (H), and the percentage of polymorphic bands (PPB). The mean coefficient of gene differentiation (Gst) was 0.0997, which indicates that only 9.97% genetic variation is between populations. The estimated gene flow (Nm) value was 4.5165, indicating a high level of gene flow and low, or medium, genetic differentiation among some populations. The results of a Mantel test revealed no significant correlation between genetic and geographic distance among populations, which means there is no evidence of significant genetic isolation by distance. An UPGMA (unweighted pair-group method with arithmetic averages) dendrogram based on genetic identity, did not indicate any major geographic structure for the 20 populations examined. mtDNA COI data indicates low nucleotide (0.0007) and haplotype diversity (0.36) in all populations. Fst values suggest that the 20 populations have low, or medium, levels of genetic differentiation. And the topology of a Neighbor-Joining tree suggests that there are no independent groups among the populations examined.
Our results suggest that C. oryzae populations have high genetic diversity at the species level. There is evidence of frequent gene flow and low, or medium, levels of genetic differentiation among some populations. There is no significant correlation between genetic and geographic distance among C. oryzae populations, and therefore no significant isolation by distance. All results are consistent with frequent gene exchange between populations, which could increase the genetic diversity, and hence, adaptability of C. oryzae, thereby promoting frequent outbreaks of this pest. Such knowledge may provide a scientific basis for predicting future outbreaks.
Chlorops oryzae (Diptera, Chloropidae) is an important pest of paddy rice plants, which inflicts significant economic damage to rice crops throughout Asia. However, the dispersal ability, potential for long-distance dispersal, pattern of migration, ecological amplitude, and population size of C. oryzae are still uncertain. Most research on this species has focused on the physiology and ecology [1,2,3], and its genetics is relatively unstudied .
Frequent C. oryzae outbreaks in recent years have caused the species to become a major pest in some regions. The propensity for outbreaks may itself play an important role in homogenizing genetic variation and intensifying gene flow between pest populations [5, 6]. We hypothesized that frequent gene flow between populations enhances the species’ overall adaptability, promoting the frequent outbreaks that occur today. In other words, the frequent outbreaks of C. oryzae are associated with the species’ genetic diversity, population demography and high rate of gene flow between populations.
To test this hypothesis, we evaluated the genetic structure of different geographic populations of C. oryzae and the level of gene flow between them. We quantified the genetic diversity and degrees of genetic differentiation of 20 different geographical populations, which may provide a scientific basis for predicting future outbreaks.
We used two effective and promising DNA markers, mitochondrial DNA (mtDNA) and inter-simple sequence repeat (ISSR), to examine between-population differences. Studies of genetic variation between pest populations can not only provide information on their population structure in different geographical regions, but also deduce the demographic history of this species [7,8,9]. Yi et al. used microsatellite and mtDNA loci to investigate the genetic divergence and dispersal ability of Bactrocera dorsalis (Hendel) on six offshore islands in South China, which results indicated that these populations have high genetic diversity, frequent gene flow and low, or medium, levels of genetic differentiation. Thus, the geographic isolation of the six islands is no barrier to the dispersal of B. dorsalis . Research on the genetic diversity and population structure of Leucinodes orbonalis collected from a variety of agro-climatic conditions found almost no genetic diversity and no significant genetic variation among the mitochondrial gene Cytochrome Oxidase I (COI) gene sequences of the populations examined. However, a few genetically distinct populations were associated with some specific habitat requirements . Similarly, genetic differentiation in ISSR markers and the COI gene among Iranian populations of Hishimonus phycitis may have been induced by geographical and ecological isolation and may have an impact on the vectoring capability of this insect .
Our results not only provide information on the genetic structure and phylogeography research of C. oryzae, but also provide a potential scientific basis for monitoring and controlling this pest.
COI gene analysis
Genetic diversity and differentiation
In total, 432 individuals collected from different locations were used to amplify 684 bp of the COI gene sequence, which defined 47 haplotypes. All 47 haplotypes had 43 variable sites, including 26 singleton variable sites and 17 parsimony informative sites. The mean total nucleotide frequencies of A, T, C and G in the nucleotide sequences from the 20 different populations were 29.98%, 36.56%, 16.95% and 16.51%, respectively, which shows an obvious AT bias (66.54%). The transition/transversion rate ratio was observed to be higher with purines (36.158) than pyrimidines (19.381). The overall transition/transversion bias was 12.022.
Genetic diversity parameters of the 20 populations and the results of neutrality tests are shown in Table 1. Haplotype diversity (Hd) for each population ranged from 0 to 0.71739 and the average number of differences (k) ranged from 0 to 1.35145. Nucleotide diversity (Π) for each population ranged from 0 to 0.00198. Tajima’s D and Fu’s Fs test of neutrality of 19 populations showed a negative value. When all samples were calculated as one population, Tajima’s D and Fu’s Fs values were negative and 1‰ significant, which is strong evidence of population expansion.
Based on our sequence database, inter- and intraspecific genetic distances of C. oryzae populations were 0.00012–0.00184 and 0–0.00198 (Table 2) respectively, which indicates no significant genetic differentiation. The Fst values between populations ranged from − 0.0595 to 0.1174 (Table 3), indicating that inter-population differences are relatively low. AMOVA results suggest that 97.28% of all genetic variation is within, and only 2.72% between, populations (Table 4).
Haplotype network and population tree
Evolutionary relationships among the haplotypes were depicted using the median-joining network method (Fig. 1). Among the 47 haplotypes, H1 was the most frequent haplotype, which occupied a central position of the network and was diversified by 46 haplotypes. Especially, H23 and H45 can be derived from two different haplotypes with just one mutational step, respectively. The topology of the C. oryzae population Neighbor-Joining tree suggested that there were no independent groups in all populations (Fig. 2). It is worth mentioning that NX population is most distant on the haplotype network and population tree, however, there is no certain landscape features around the NX location. The reason for preventing emigration of individuals is still unknown.
197 fragments were generated from the 9 primers, ranging in size from 250 to 2000 bp, and with an average of 21.89 fragments per primer. The percentage of polymorphic bands was 100% showing high genetic diversity at the species level. Among the 20 populations, PPB values ranged from 89.85% to 99.49%. The Na and Ne were 2.0000 ± 0.0000 (1.8985 ± 0.3028 to 1.9949 ± 0.0712 among populations) and 1.6460 ± 0.2373 (1.4992 ± 0.3270 to 1.6908 ± 0.2784 among populations), respectively. The H and I were 0.3785 ± 0.0984 (0.2992 ± 0.1588 to 0.3893 ± 0.1196) and 0.5609 ± 0.1151 (0.4547 ± 0.2091 to 0.5692 ± 0.1457), respectively, within all samples (Table 5).
The Gst was 0.0997, which indicates that 90.03% genetic variation is within populations and only 9.97% between populations. The estimated Nm value was 4.5165 (Table 6). These results suggest that genetic differentiation among populations of C. oryzae is impeded by high gene flow. Table 7 lists the genetic identity (above diagonal) and genetic distance (below diagonal) among populations.
The relationship between genetic and geographic distance was shown in Fig. 3. A Mantel test revealed no significant correlation between genetic and geographic distance (r = 0.54675, p = 0.9992) among C. oryzae populations, and there was therefore no evidence of significant isolation by distance.
An UPGMA dendrogram constructed based on genetic identity (Fig. 4) grouped the 20 populations into two major clusters. The dendrogram did not reveal any major geographic structure for these populations.
The rapid development of molecular techniques has made it possible to directly measure genetic differentiation and genetic diversity among populations . The faster mutation rate and relatively conserved sequence of the mtDNA COI gene is ideally suited to species identification via DNA barcoding . Analysis of mtCOI gene variation is regarded as an important and reliable tool for defining cryptic species , evaluating biodiversity , identifying samples , and distinguishing closely related species .
Various surveys have demonstrated the reliability of ISSR markers, which can generate more polymorphisms than either RAPD or RFLP, . For example, Dioscorea hispida were grouped into 10 vital groups based on information that provided by ISSR markers, proving the existence of significant variation among germplasm specimens. D. hispida shows a high level of genetic diversity among accessions, which suggests that ISSR markers have been very effective in detecting polymorphism in this species . ISSR primers have been used to determine the potential for the diversification of cassava crops in Angola, revealing genetic diversity within populations and genetic information sharing among the three main taxa . The ISSR molecular marker technique has also been used to distinguish between citrus rootstock species, and to reveal the broad genetic base and high genetic variability among these .
Genetic diversity, also known as gene diversity, is the foundation of biodiversity that guarantees the evolution of species. The percentage of polymorphism, haplotype diversity and nucleotide diversity are used to evaluate population genetic diversity . High levels of genetic diversity are indicative of the strong viability and adaptability of populations . Results from this study indicate that the genetic diversity of C. oryzae was high between populations that were sampled. This may increase fitness in populations to changing conditions. Crawford and Whitney  also showed that genetic diversity increases the ability of species to colonize on a short-term ecological timescale by increasing the possibility of population survival, growth and reproduction under novel environments. Episyrphus balteatus and Sphaerophoria scripta European populations successfully adapt to changing environmental conditions and have great colonization abilities due to the high genetic diversity . The Tajima’s D and Fu’s Fs test of neutrality values for 19 populations were negative. When all samples were calculated as one population, Tajima’s D and Fu’s Fs values were negative and 1‰ significant, strongly indicating population expansion. This result is consistent with the idea that the negative Tajima’s D and Fu’s Fs test values demonstrated demographic expansion has occurred in these populations.
Fst values between C. oryzae populations ranged from − 0.0595 to 0.1174, and the Gst value was 0.0997, both of which are indicative of low genetic differentiation between populations. This suggests frequent gene flow between populations, which may increase the species’ adaptability to environmental change . Moreover, gene flow can not only demonstrate the probable genetic differentiation and genetic infiltration among populations, but also reduce the genetic differences among populations . In this research, the Nm value of C. oryzae was 4.5165 which indicates high gene flow and low, or medium, genetic differentiation among some populations. High gene flow may impede genetic differentiation in C. oryzae. Gene flow (or lack of gene flow) plays a crucial role in genetic differentiation, affects the overall adaption of entire species and adaptative divergence between populations . It has been traditionally considered as a homogeneous force that limits adaptive differences [29, 30], and recent studies have shown that it can also promote adaptation to local environmental conditions [31, 32]. For example, moderate gene flow increases the adaptation capabilities of Rhagoletis cerasi populations (which occupy different habitates in fragmented landscapes) to local habitates, thus preventing them from becoming extinct due to genetic processes [33, 34].
An AMOVA based on Fst values indicates that most of the genetic variation was resulted from the difference within populations. Furthermore, a Mantel test indicates no significant correlation between genetic and geographic distance. This result is consistent with the findings of Yang et al.  who compared correlation of the symmetric matrix constituted by geographic and genetic distances to analyze the existence of isolation among populations of Odontotermes formosanus in different regions. These authors found no significant correlation between geographic distance and genetic distance and no significant isolation by distance. Overall, we found a high level of genetic diversity and a low degree of population differentiation among populations of C. oryzae, and the gene flow was unaffected by geographic distance. Similarly, geographic distance did not appear to affect gene flow between 10 geographically separated populations Oedaleus infernalis . Fst and Gst values for these populations are low, and the gene flow is high, indicating a low level of genetic differentiation and high gene flow among populations . The correlation between genetic and geographic distance was insignificant .
Our results provide important, new information on the genetic diversity and genetic differentiation of C. oryzae, and suggest that high gene flow between populations contributes to the now frequent outbreaks of this pest. However, further research on both additional geographical populations and different genetic markers are necessary before definitive conclusions can be reached. Furthermore, future work can focus on doing a more comprehensive ecological and behavioural research to understand the natural history of C. oryzae in greater detail.
This study showed that the now frequent outbreaks of C. oryzae may due to high gene flow between populations. We have found that these populations have high genetic diversity at the species level, whereas exhibited low genetic differentiation. High genetic diversity and frequent gene flow between populations may enhance the tolerance of populations to environmental variability and increase the adaptability to novel environmental pressures, leading to frequent outbreaks what had happened and what will happen in a large scale.
Sample collection and DNA extraction
400 specimens of C. oryzae were collected from different parts of Hunan province, China, and an additional 32 specimens from Zhejiang and Guizhou provinces (Fig. 5, Additional file 1: Table S1). Samples were soaked in 100% ethanol and stored at − 20 °C until their genomic DNA (gDNA) was isolated. After removing the residual ethanol, DNA was extracted from each individual using an Ezup Column Animal Genomic DNA Extraction Kit as per the manufacturer’s instructions (Sangon Biotech, Shanghai, China).
COI PCR amplification and sequencing
The COI was amplified with the COIF (5′-CTAGGTGCTCCAGATATAGCATTTC-3′) and COIR (5′-GGCTAAAACAACTCCTGTTAATCC-3′) primers from isolated DNA. PCR was performed in 20 μL volumes comprised of 10 μL PrimeSTAR Max DNA Polymerase (TaKaRa, Tokyo, Japan), 1 μL of each primer (10 mmol/L), 1 μL of template DNA solution (70 ng/μL), and 7 μL double distilled water. Amplifications were conducted as follows: 34 cycles of denaturation at 94 °C for 30 s, annealing at 59 °C for 30 s, and extension at 72 °C for 1 min. All the PCR products were checked by electrophoresis on a 1.2% agarose gel and bidirectional sequencing was completed by TSINGKE (Beijing, China).
ISSR PCR amplification
A total of sixteen primers from the University of British Columbia Biotechnology Laboratory Primer kit No.9 were tested for PCR and nine (Additional file 1: Table S2) that could produce reproducible, clear, polymorphic electrophoretic bands were chosen for further analysis. PCR was performed in 20 μL volumes comprised of 10 μL Premix Taq™ (Ex Taq™ Version 2.0 plus dye) (TaKaRa, Tokyo, Japan), 2 μL of primer (10 mmol/L), 1 μL of template DNA solution (70 ng/μL), and 7 μL double distilled water. Amplifications were carried out as follows: an initial denaturing at 94 °C for 3 min, followed by 34 cycles of denaturing at 94 °C for 30 s, annealing at an optimized temperature for 30 s, and extension at 72 °C for 1 min, with a final extension 7 min at 72 °C. All PCR products were electrophoretically separated on 2% agarose.
COI gene data analysis
COI sequences were edited manually with BioEdit v.7.0.9 to produce consensus sequences of 685 bp for each specimen . All indices for sequence polymorphic sites, DNA polymorphism, genetic differentiation, neutrality tests [Tajima’s D  and Fu’s Fs  ], and haplotype analyses were executed using DnaSp v.5.10 . A haplotype network, which included haplotype frequencies, was calculated using Network v.4.6 . Intra- and inter-specific genetic distances and transition/transversion ratios in each codon were computed based on COI gene sequences using MEGA v.7.0 . A population phylogenetic tree based on genetic distances was constructed using the Neighbor-Joining tree model in MEGA v.7.0. Analysis of molecular variance (AMOVA) was performed with Arlequin software v.3.5.2 .
ISSR data analysis
Amplified ISSR fragments were scored as present (1) or absent (0) according to the molecular weight (bp) and the resulting matrix of binary values was used for further analyses. The observed number of alleles (Na), effective number of alleles (Ne), Nei’s gene diversity (H), Shannon’s information index (I), the percentage of polymorphic bands (PPB), total gene diversity (Ht), genetic diversity within populations (Hs), coefficient of gene differentiation (Gst), and Gene flow (Nm) were calculated using POPOGENE v.1.31 . Cluster analysis was used to construct dendrograms using the UPGMA (unweighted pair-group method with arithmetic averages) in NTSYSpc v.2.1 . To estimate the existence of isolation among 20 C. oryzae populations, the geographic distances between them were calculated using the google-maps-distance-calculator (www. daftlogic.com/projects-google-maps-distance-calculator. htm) (Additional file 1: Table S3). The correlation between geographic and genetic distances was further analyzed using the MXCOMP program and a Mantel test in NTSYSpc v.2.1.
Availability of data and materials
The authors declare that the data supporting the finding of this study are available in the article and its additional files.
Inter-simple sequence repeat
- COI :
Mitochondrial gene Cytochrome Oxidase I
Analysis of molecular variance
- Na :
The observed number of alleles
- Ne :
Effective number of alleles
- H :
Nei’s gene diversity
- I :
Shannon’s information index
The percentage of polymorphic bands
- Ht :
Total gene diversity
- Hs :
Genetic diversity within populations
- Gst :
Coefficient of gene differentiation
- Nm :
Unweighted pair-group method with arithmetic averages
- Hd :
- k :
The average number of differences
- Π :
Takeda M. Genetic basis of photoperiodic control of summer and winter diapause in geographic ecotypes of the rice stem maggot, Chlorops oryzae. Entomol Exp Appl. 1998;86(1):59–70.
Takeda M, Nagata T. Photoperiodic responses during larval development and diapause of two geographic ecotypes of the rice stem maggot, Chlorops oryzae. Entomol Exp Appl. 1992;63(3):273–81.
Takeda M. Effects of photoperiod and temperature on larval development and summer diapause in two geographic ecotypes of the rice stem maggot, Chlorops oryzae Matsumura (Diptera: chloropidae). Appl Entomol Zool. 1997;32(1):63–74.
Qiu L, Tao S, He H, Ding W, Li Y. Transcriptomics reveal the molecular underpinnings of chemosensory proteins in Chlorops oryzae. BMC Genom. 2018;19(1):890.
Chapuis MP, Lecoq M, Michalakis Y, Loiseau A, Sword GA, Piry S, Estoup A. Do outbreaks affect genetic population structure? A worldwide survey in Locusta migratoria, a pest plagued by microsatellite null alleles. Mol Ecol. 2008;17(16):3640–53.
Chapuis MP, Loiseau A, Michalakis Y, Lecoq M, Franc A, Estoup A. Outbreaks, gene flow and effective population size in the migratory locust, Locusta migratoria: a regional-scale comparative survey. Mol Ecol. 2009;18(5):792–800.
Assefa Y, Mitchell A, Conlong DE. Phylogeography of Eldana saccbarina Walker (Lepidoptera: Pyralidae). Ann Soc Entomol Fr. 2006;42(3–4):331–7.
Assefa Y, Conlong DE, Van Den Berg J, Martin LA. Ecological genetics and host range expansion by Busseola fusca (Lepidoptera: Noctuidae). Environ Entomol. 2015;44(4):1265–74.
Assefa Y, Goftishu M, Capdevielle-Dulac C, Ru BL. Clarifying the source of Conicofrontia sesamoides Hampson (Lepidoptera: Noctuidae) population in South African sugarcane using morphological identification and mitochondrial DNA analysis. Phytoparasitica. 2017;45(1):45–55.
Yi C, Zheng C, Zeng L, Xu Y. High genetic diversity in the offshore island populations of the tephritid fruit fly Bactrocera dorsalis. BMC Ecol. 2016;16(1):46.
Palraju M, Paulchamy R, Sundaram J. Population genetic structure and molecular diversity of Leucinodes orbonalis based on mitochondrial COI gene sequences. Mitochondrial DNA A DNA Mapp Seq Anal. 2018;29(8):1231–9.
Hemmati C, Moharramipour S, Seyahooei MA, Bagheri A, Mehrabadi M. Population genetic structure of Hishimonus phycitis (Hem.: Cicadellidae), vector of lime witches’ broom phytoplasma. J Agric Sci Technol. 2018;20(5):999–1012.
Mutun S, Borst DW. Intraspecific mitochondrial DNA variation and historical biogeography of the eastern lubber grasshopper, Romalea microptera. Ann Entomol Soc Am. 2004;97(4):681–96.
Saccone C, De Giorgi C, Gissi C, Pesole G, Reyes A. Evolutionary genomics in Metazoa: the mitochondrial DNA as a model system. Gene. 1999;238(1):195–209.
Bucklin A, Wiebe PH, Smolenack SB, Copley NJ, Beaudet JG, Bonner KG, Färberlorda J, Pierson JJ. DNA barcodes for species identification of Euphausiids (Euphausiacea, Crustacea). J Plankt Res. 2007;29(6):483–93.
Smith MA, Fisher BL, Hebert PD. DNA barcoding for effective biodiversity assessment of a hyperdiverse arthropod group: the ants of Madagascar. Philos Trans R Soc Lond B Biol Sci. 2005;360(1462):1825–34.
Lakra WS, Goswami M, Gopalakrishnan A, Singh DP, Singh A, Nagpure NS. Genetic relatedness among fish species of Genus Channa using mitochondrial DNA genes. Biochem Syst Ecol. 2010;38(6):1212–9.
Hebert PD, Ratnasingham S, deWaard JR. Barcoding animal life: cytochrome c oxidase subunit 1 divergences among closely related species. Proc Biol Sci. 2003;270(Suppl 1):S96–9.
Prevost A, Wilkinson MJ. A new system of comparing PCR primers applied to ISSR fingerprinting of potato cultivars. Theor Appl Genet. 1999;98(1):107–12.
Nudin NFH, Ali AM, Ngah N, Mazlan NZ, Mat N, Ghani MNA, Alias N, Zakaria AJ, Jahan MS. ISSR marker-assisted genetic diversity analysis of Dioscorea hispida and selection of the best variety for sustainable production. C R Biol. 2017;340(8):359–66.
Sandra DJA, Ricardo FCM, Carlos SL, Cláudia FF, Vanderlei SS, Pascoal AM. Genetic structure of cassava populations (Manihot esculenta Crantz) from Angola assessed through (ISSR) markers. Afr J Biotechnol. 2019;18(7):144–54.
Salis C, Papadakis IE, Kintzios S, Hagidimitriou M. In vitro propagation and assessment of genetic relationships of citrus rootstocks using ISSR molecular markers. Not Bot Horti Agrobo. 2017;45(2):383–91.
Barrett RDH, Schluter D. Adaptation from standing genetic variation. Trends Ecol Evol. 2008;23(1):38–44.
Crawford KM, Whitney KD. Population genetic diversity influences colonization success. Mol Ecol. 2010;19:1253–63.
Raymond L, Plantegenes M, Vialatte A. Migration and dispersal may drive to high genetic variation and significant mixing: the case of two agriculturally important, continental hoverflies (Episyrphus balteatus and Sphaerophoria scripta). Mol Ecol. 2013;22:5329–39.
Xu Y, Mai JW, Yu BJ, Hu HX, Yuan L, Jashenko R, Ji R. Study on the genetic differentiation of geographic populations of Calliptamus italicus (Orthoptera: Acrididae) in Sino-Kazakh border areas based on mitochondrial COI and COII genes. J Econ Entomol. 2019;112(4):1912–9.
Millar CI, Libby WJ, Falk DA, Holsinger KE. Strategies for conserving clinal, ecotypic, and disjunct population diversity in widespread species. Genet Conserv Rare Plants. 1991;149:170.
Garant D, Forde SE, Hendry AP. The multifarious effects of dispersal and gene flow on contemporary adaptation. Funct Ecol. 2007;21(3):434–43.
Barton N, Partridge L. Limits to natural selection. BioEssays. 2000;22:1075–84.
Lenormand T. Gene flow and the limits to natural selection. Trends Ecol Evol. 2002;17:183–9.
Gandon S, Nuismer SL. Interactions between genetic drift, gene flow, and selection mosaics drive parasite local adaptation. Am Nat. 2009;173:212–24.
Ribeiro A, Llord P, Bowie RCK. A tight balance between natural selection and gene flow in a southern African arid-zone endemic bird. Evol. 2011;65:3499–514.
Lynch M, Conery J, Burger R. Mutation accumulation and the extinction of small population. Am Nat. 1995;146:485–518.
Crnokrak P, Roff DA. Inbreeding depression in the wild. Heredity. 1999;83:260–70.
Yang YQ, Pu LJ, Wang Q, Wang Z, Pang ZP, Long YH. Population diversity of Odontotermes formosanus (Shiraki) (Termitidae, Macrotermitinae) from different geographic locations in Anhui province. China. Sociobiol. 2018;65(3):497–505.
Sun W, Dong H, Gao YB, Su QF, Qian HT, Bai HY, Zhang ZT, Cong B. Genetic variation and geographic differentiation among populations of the nonmigratory agricultural pest Oedaleus infernalis (Orthoptera: Acridoidea) in China. J Insect Sci. 2015;15:150.
Hall TA. BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucleic Acid Symp Ser. 1999;41:95–8.
Tajima F. Statistical method for testing the neutral mutation hypothesis by DNA polymorphism. Genetics. 1989;123(3):585–95.
Fu YX. Statistical tests of neutrality of mutations against population growth, hitchhiking and background selection. Genetics. 1997;147(2):915–25.
Librado P, Rozas J. DnaSP v5: a software for comprehensive analysis of DNA polymorphism data. Bioinformatics. 2009;25(11):1451–2.
Polzin T, Daneshmand SV. On Steiner trees and minimum spanning trees in hypergraphs. Oper Res Lett. 2003;31:12–20.
Kumar S, Stecher G, Tamura K. MEGA7: molecular evolutionary genetics analysis version 7.0 for bigger datasets. Mol Biol Evol. 2016;33(7):1870–4.
Excoffier L, Lischer HE. Arlequin suite ver 3.5: a new series of programs to perform population genetics analyses under Linux and Windows. Mol Ecol Resour. 2010;10(3):564–7.
Yeh FC, Yang RC, Boyle T, Ye ZH, Mao JX. Popgene, the user friendly shareware for population genetic analysis. Edmonton: Molecular Biology and Biotechnology Center. University of Alberta; 1997.
Rohlf FJ. NTSYS-pc Numerical taxonomy and multivariate analysis system’, Version 2.02e. Setauket: Exeter Software; 1997.
We wish to thank Zhongcai Li, Xinwen Li, Xiaoping Tan and Zhengbing Zhang for their help in collecting the material in this study.
This study was supported by Changsha City Science and Technology Bureau Key Projects, China (kc1701003).
Double first-class construction project of Hunan Agricultural University (SYL2019029)
Ethic approval and consent to participate
No permissions were required for sampling. There is no any regional/national legislature/law require ethical approval for sampling.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Zhou, A., Tian, P., Li, Z. et al. Genetic diversity and differentiation of populations of Chlorops oryzae (Diptera, Chloropidae). BMC Ecol 20, 22 (2020). https://doi.org/10.1186/s12898-020-00293-8
- Chlorops oryzae
- Mitochondrial DNA
- Genetic diversity
- Genetic differentiation