基于微卫星分子标记的二斑叶螨和朱砂叶螨种群遗传结构研究
详细信息    本馆镜像全文|  推荐本文 |  |   获取CNKI官网全文
摘要
二斑叶螨(Tetranychus urficae Koch)与朱砂叶螨(T cinnabarinus Boisduval)同属于节肢动物门蛛形纲Arachnida、蜱螨亚纲Acari、真螨目Acariformes、叶螨科Tetranychidae、叶螨属Tetranychus,广泛分布于世界各地,寄主植物广泛,包括果树、花卉、蔬菜、农作物及杂草等多种植物,一般以成螨、若螨群集在叶背面吸食植物汁液,对寄主植物产生机械伤害和毒害作用,引起植物减产和品质下降,给农业生产带来严重的经济损失,这已经成为农业发展中的突出问题。由于两者扩散迅速、世代周期短、发育速度快、繁殖能力强、所以极易在新的环境定殖,形成新的地理种群。本研究利用微卫星标记的方法,揭示中国二斑叶螨和朱砂叶螨的种群遗传结构,为了解其入侵途径和爆发机制,有效地控制其扩散和危害奠定基础。通过比较两者种群遗传结构的差异,为进一步揭示两者的亲缘关系提供可靠的依据。
     一、实验材料的采集和饲养
     在全国范围内选择较有代表性的地区分别对二斑叶螨和朱砂叶螨进行采样,共采集了7个二斑叶螨地理种群和6个朱砂叶螨地理种群,采集的生境主要是果园、花圃和蔬菜地。实验室内采用隔水法在菜豆(Phaseolus vulgaris L.)叶片上进行恒温恒湿隔离饲养,饲养条件为25±1℃,L16:D8,相对湿度60%。
     二、微卫星标记技术在种群遗传结构中的应用
     微卫星标记是种群遗传结构研究中应用最多的一种分子标记方法。因其具有在基因组中分布广泛、多态信息含量高、共显性标记、无损伤取样、一旦开发易于检测等诸多优点,而被越来越多的应用到分子生态、群体遗传学等研究领域。自Navajas等1998年和2002年在二斑叶螨体内成功开发出微卫星位点以来,越来越多的研究者将微卫星标记用于螨的种群遗传结构的研究。
     实验方法:首先需进行单位点扩增,成功扩增物种体内已知的微卫星序列;然后进行各个微卫星位点等位基因和基因型检测,本研究使用的方法是银染—非变性聚丙烯酰胺凝胶电泳法;其次利用相关软件进行数据分析,如种群遗传多样性参数统计,包括等位基因频率、基因型频率、等位基因丰富度、杂合度、基因型多样性等;种群之间分化研究,包括遗传分化指数、基因流、遗传距离等;最后基于种群之间遗传距离的远近,采用距离矩阵法(UPGMA)构建系统发育树,更加直观的反应种群之间的亲缘关系,揭示种群遗传结构,为监控害螨种群涨落趋势,及时有效的对其进行防治具有一定的指导意义。
     三、中国二斑叶螨种群遗传结构的研究
     本研究所选用的3个微卫星位点在二斑叶螨体内成功扩增,扩增结果显示每个位点平均等位基因数3.7个,平均多态信息含量(PIC)0.5247,当PIC>0.5时.该位点为高度多态位点,所以这3个位点作为遗传标记进行二斑叶螨种群遗传结构的研究是切实可行的。表征种群遗传多样性的参数指标(N_A、N_(AE)、H_O和H_E等)在二斑叶螨各个种群中普遍偏低,种群遗传多样性降低。二斑叶螨各个地理种群检测到的私有等位基因占等位基因总数的27.3%,特有基因型占基因型总数的54.5%,单态位点数占检测总数的33%,总样本的基因型多样性远高于各个地理种群的基因型多样性水平,遗传分化指数平均成对F_(ST)值为0.5695,当F_(ST)>0.25时,种群极度分化。以上结果均说明二斑叶螨各个地理种群之间发生了明显的分化。种群之间遗传分化与地理距离呈现一定的相关性(P<0.05,R~2=0.2594)。
     二斑叶螨种群遗传结构表现出两个明显的特点:种群遗传多样性降低和种群之间极度分化。造成以上结果的原因有生物本身的特性(如生殖系统)、分布状态(如栖境破碎、地理隔离等)、自然选择、迁徙(基因交流)、遗传漂变等多种因素综合作用的影响,其中地理隔离引起的小居群范围内的近亲交配和等位基因遗传漂变很可能起到了决定性的作用。
     四、中国朱砂叶螨种群遗传结构的研究
     朱砂叶螨种群遗传结构的研究结果与二斑叶螨有许多类似的地方,如以总体样本为研究对象均能表现出显著的多态性,然而各个地理种群的遗传多样性水平普遍偏低;种群之间极度分化(F_(ST)>>0.25),但是却发现许多种群之间遗传分化与地理距离不成比例的例子;从统计上确定遗传分化程度与地理距离之间的相关关系,结果显示虽然地理距离与种群分化的相关系数很低(二斑叶螨R~2=0.2594;朱砂叶螨R~2=0.277),但是仍然表现出显著的相关性(P<0.05)。所以朱砂叶螨种群遗传结构也表现出种群遗传多样性降低和种群之间极度分化的特点。
     五、我国二斑叶螨和朱砂叶螨种群遗传结构的差异分析
     中国二斑叶螨和朱砂叶螨的种群遗传结构经比较差异显著,主要表现在种群遗传多样性水平不同,朱砂叶螨的总体样本的遗传多样性水平明显高于二斑叶螨,而二斑叶螨各个地理种群的遗传多样性水平却比朱砂叶螨高;种群遗传分化程度不同,朱砂叶螨各个种群之间平均成对F_(ST)值远高于二斑叶螨,说明朱砂叶螨种群遗传分化的程度远高于二斑叶螨;将采自同一个省市的二斑叶螨和朱砂叶螨进行成对比较发现,地理间距缩短并没有降低种群之间的遗传分化水平,说明地理隔离并不是朱砂叶螨和二斑叶螨分化的主要原因,生殖隔离也可能造成同域内种群分化;系统发育分析结果显示7个二斑叶螨地理种群与6个朱砂叶螨地理种群明显聚在不同的分支上。
     综上所述,本研究利用微卫星标记技术,揭示中国二斑叶螨和朱砂叶螨的种群遗传结构及两者的差异,这将有助于我们了解它们的入侵途径和爆发机制,为有效控制其扩散和危害奠定基础,同时,也有助于了解这两种叶螨的分类地位,为现在具有极大争议的同种或异种学说提供科学依据。
Belonging to the Arachinida, Acari, Acariformes, Tetranychidae and Tetranychus, both the two-spotted spider mite Tetranychus urticae Koch and the carmine spider mite T. cinnabarinus (Boisduval) are widely distributed worldwide. The two spider mites have wide range of host plants, including fruit trees, flowers, vegetables, crops and weeds. They live aggregately on the undersurfaces of leaves to suck plant juice, causing damage to host plant and great loss to agriculture. Owing to their rapid development and high reproductive capacities, new geographical populations of these mites are established easily in other areas. In this study, we used the molecular marker of microsatellite to analyze the population genetic structure of T. urticae Koch and T. cinnabarinus (Boisduval) from China. This work will be helpful for understanding their invading approach and outbreak mechanism, laying a solid foundation for controlling their dispersal.
     1. The collection and rearing of experimental material
     Samples of T. urticae Koch and T. cinnabarinus (Boisduval) used in this study were collected from seven and six geographical localities, respectively. Their collections cover an extensive geographical scale representing different regions and environments in China. Mites were then reared separately in the laboratory on leaves of Phaseolus vulgaris L., placed on a water-saturated tissue in a glass dish. All populations were kept in an illumination culture tank at 25±1℃, 60% R.H. and under L16:D8 conditions.
     2. The application of microsatellite markers in the studies of population genetic structure.
     Microsatellite, as a molecular marker, is an important tool in the study of population genetic structure. It has a lot of advantages such as extensive in eukaryotic genome, richness of polymorphism information content, codominant genetic marker, easily to be investigated. Because of many advantages, it has often been used in the fields of molecular ecology, population genetics, and so on. Microsatellite has been used in the population genetic structure of the mite by more and more researchers after successfully exploiting microsatellite locus in Tetranychus urticae Koch by Navajas in 1998 and 2002.
     The method of this experiment was as follows: (1) to amplify known loci successfully; (2) to investigate alleles and allelic genotype at every locus in every sample; (3) PCR products were separated on polyacrylamide gel and silver staining to showing bands; (4) different softwares were used to analyze parameters of population genetic diversity such as allelic frequency, genotypic frequency, allelic richness, heterozygosity, genotypic diversity, and so on. The study of population genetic differentiation includes genetic differentiation index, gene flow, genetic distance, and so on. Finally, the phylogenetic tree based on genetic distances between populations was established. Cluster analysis was performed using the UPGMA algorithm.
     3. The population genetic structure of Tetranychus urticae Koch from China
     The three microsatellite loci selected in this research were amplified successfully. Amplification results showed that a mean number 3.7 of alleles per locus, the average polymorphism information content (PIC) 0.5247. When the value of PIC is bigger than 0.5, the site is highly polymorphic. So these three loci as the genetic marker for the study of T. urticae Koch population genetic structure are practicable. Parameters standing for population genetic diversity is generally low in various geographical populations, these results illustrate the low level genetic diversity of T. urticae Koch populations. The rates of private allele, special genotype and monomorphic locus detected were 27.3%, 54.5% and 33%, respectively. The genotypic diversity of total samples was much higher than respective population. Average pairwise F_(ST) was 0.5695, when F_(ST)>0.25, population was divided extremely. The above results show that extreme division has happened to T. urticae Koch. Genetic differentiation between populations had a direct bearing on the geographical distance.
     The result of population genetic structure of T. urticae Koch showed two distinct characteristics: the lower level of population genetic diversity and extreme differentiation among populations. These results caused by intrinsic biological feature (such as mating system), distribution (such as habitat fragmentation, geographical isolation etc), natural selection, migration (gene flow), genetic drift, etc. Therefore, both inbreeding within populations and genetic drift of alleles have probably had a decisive effect on the population genetic structure of T. urticae Koch.
     4. The population genetic structure of T. cinnabarinus (Boisduval) from China
     The population genetic structure of T. cinnabarinus (Boisduval) displayed many similarities to T. urticae Koch. The overall sample study showed significant polymorphism, but all the geographical populations had generally lower levels of genetic diversity. There existed an extreme differentiation between populations (F_(ST)>> 0.25). However, the genetic differentiations among populations were disproportionate with geographical distance. The correlation analysis between the genetic differentiation and the geographical distance from the point view of statistics showed a low correlation coefficient between the geographic distance and the population division (R~2 = 0.277; R~2 = 0.2594), and a significant correlation (P<0.05). Therefore, the population genetic structure of T. cinnabarinus (Boisduval) also showed a lower level of genetic diversity and the extreme polarization characteristics between populations.
     5. The differences of population genetic structure between Tetranychus urticae Koch and T. cinnabarinus (Boisduval) from China.
     By comparing the population genetic structures of T. urticae Koch and T. cinnabarinus (Boisduval) from China, we found significant differences between them. Firstly, they had different levels of genetic diversity of populations in these two species. The overall level of genetic diversity of T. cinnabarinus (Boisduval) was significantly higher than that of T. urticae Koch, but the various geographical levels of genetic diversity in T. cinnabarinus (Boisduval) were oppositely lower than that of T. urticae Koch. Secondly, the degree of population differentiation between the two spider species was different. Average pairwise F_(ST) between geographical populations of T. cinnabarinus (Boisduval) was far higher than that of T. urticae Koch. This result demonstrated the level of population genetic differentiation of T. cinnabarinus (Boisduval) was much higher than that of T. urticae Koch. Thirdly, by comparing T. cinnabarinus (Boisduval) and T. urticae Koch collected in the same area; we could not find the reduction of the level of genetic differentiation due to the shortening geographical distances between different populations. Geographic isolation between populations wasn't the main reason for genetic differentiation; reproductive isolation could also result in differentiation among populations in the same area. Finally, phylogenetic analysis revealed that seven geographic populations of T. urticae Koch and six geographic populations of T. cinnabarinus (Boisduval) were gathered in different branches.
     In conclusion, by using microsatellite marker, we revealed population genetic stricture of T. urticae Koch and T. cinnabarinus (Boisduval) and differences between them. It will help us understand their ways of invasion and outbreak mechanism, and lay foundations for effectively controlling them. Additionally, these results will also be helpful for understanding their classification status, and offering the scientific evidence to homogeneous or heterogeneous theory.
引文
1.蔡双虎,程立生.二斑叶螨的研究进展.热带农业科学,2003,23(2):68-74.
    2.陈志杰,张淑莲,张美荣.麦棉间套棉叶螨种群动态与生态控制技术研究.陕西师范大学学报(自然科学版),1997,25:137-142.
    3.谌有光.警惕二斑叶螨在我国果产区蔓延为害.植保技术与推广,1997,17(5):24-25.
    4.程立生.二斑叶螨的检疫重要性.植物检疫,1994,8(3):164.
    5.程立生.中国朱砂叶螨各地理种群形态变异研究.热带作物学报,1998,19(1):83-86.
    6.何飞,霍金龙,曾养志.微卫星DNA分子标记技术及其在动物学研究中的应用.动物养殖,2005,6:35-38.
    7.匡海源,程立生.关于区分朱砂叶螨和二斑叶螨两个近似种的研究.昆虫学报,1990,33(1):109-116.
    8.匡海源.农螨学.北京:农业出版社,1986.
    9.蓝明扬,沈伟华,黄超群,等.亚洲璃眼蜱基因组多态性DNA的研究.中国媒介生物学及控制杂志,1996,7(2):101-104.
    10.李明,亓小红.微卫星标记的基因组学研究进展.天津农学院学报,2006,13:47-51.
    11.刘群红,张浩.两种粉螨的酯酶同工酶和相关蛋白质电泳的研究.锦州医学院学报,2001,22(3):8-9.
    12.马淑娥,王家红,徐春荣.二斑叶螨发规律及防治技术.中国果树,1998,2:34-36.
    13.孟和生,王开运,姜兴印,等.二斑叶螨发生危害特点及防治对策.昆虫知识,2001,38(1):52-54.
    14.乔中东,殷国荣,赵瑞君,等.四种硬蜱的随机扩增多态性DNA分析.中国人兽共患病杂志,1997,13(1):54-56.
    15.曲鲁江,李显耀,杜志强等.微卫星PCR产物变性与非变性PAGE-银染检测方法的比较.遗传,2004,26(4):522-524.
    16.王冬生.蔬菜害虫治理技术-第二讲害螨的发生与防治.上海蔬菜,2001(4):41-42.
    17.吴孔明.朱砂叶螨为害对棉花生理活动的影响.植物保护学报,1989,16(2):99-106.
    18.谢霖.中国二斑叶螨和朱砂叶螨种群分子遗传结构的研究.南京农业大学硕士论文,2006.
    19.晏鹏,吴孝兵,史燕,等.微卫星多态性检测技术及其在保护遗传学中的应用.应用生态学报,2003,14(3):461-464.
    20.周玉书,朴春树,刘池林.警惕二斑叶螨在北方果区为害蔓延.植物保护,1996,22(5):S1-52
    21.Ali S,Muller C R,Epplen J H.DNA fmgerprinting by oligonucleotide probe specific for simple repeats. Hum Genet, 1986, 74: 239-246.
    22. Antonovics J. Evolution in closely adjacent plant populations. IV. Manifold effects of gene flow. Heredity, 1968, 23: 507-524.
    23. Avise J C. Phylogeography: The history and formation of species. Harvard University Press, Cambridge, MA, 2000.
    24. Bailly X, Migeon A, Navajas M. Aanalysis of microsatellite variation in the spider mite pest Tetranychus turkestani (Acari: Tetranychidae) reveals population genetic structure and raises questions about related ecological factors. Bio J Linn Soc, 2004 (82): 69 - 78.
    25. Botstein D, White R L, Skolnick M, Davis R W. Construction of a genetic linkage map in man using restriction fragment length polymorphisms. Am J Human Genet, 1980, 32: 314 - 331.
    26. Boudreaux H B. Biological aspects of some phytophagous mites. Ann Rev Enotomol, 1963, 8: 137-154.
    27. Boudreaux H B. Revision of the two-spotted spider mite complex, Tetranychus telarius (Linnaeus). Ann Entomol Soc Am, 1956, 49: 43-48.
    28. Boudreaux H B. The effect of relative humidity on egg-laying, hatching, and survival in various spider mites. J Insect Physiol, 1958, 2: 65-72.
    29. Carbonnelle S, Hance T, Migeon A, Baret P, Cros-Arteil S, Navajas M. MicroSatellite marker reveal spatial genetic structure of Tetranychus urticae (Acari: Tetranychidae) population along a latitudinal gradient in Europe. Exp Appl Acarol, 2007,41: 225 - 241.
    30. Coote T, Bruford M W. Universally applicable microsatellites for analysis of genetic variation in apes and old world monkeys. J Heredity, 1996, 87: 406-410.
    31. Coyne J A, Orr H A. Patterns of speciation in Drosophila. Evolution, 1989, 43: 362-381.
    32. Crow J F, Kimura M. An introduction to population genetics theory. New York: Harper & Row, 1970, 95-145.
    33. de Boer R. Partial hybrid sterility between strains of the arrhenotokous spider mite, Tetranychus urticae complex (Acari: Tetranychidae) . Genetica, 1982, 58: 23-33.
    34. Delaye C, Aeschlimann A, Renaud F, Rosenthal B, De Meeus T. Isolation and characterization of microsatellite markers in the Ixodes ricinus complex (Acari: Ixodidade) .Mol Ecol, 1998, 7: 357-363.
    35. Dupont L M. On gene flow between Tetranychus urticae Koch, 1836 and Tetranychus cinnabarinus (Boisduval) Boudreaux, 1956 (Acari: Tetranychidae): synonymy between the two species. Entomol Exp Appl , 1979, 25: 297-303.
    36. Ellegren H. Microsatellite mutations in the germline: implications for evolutionary inference. Trends Genet, 2000b, 16: 551-558.
    37. Fenton B, Birch A N E, Malloch G, Lanham P G, Brennan R M. Gall mite molecular phylogeny and its relationship to the evolution of plant host specificity. Exp Appl Acarol, 2000, 24: 831 - 861.
    38. Feselstein (2001) PHYLIP: Phylogeny Inference Package, version 3.6, Department of genetics. University of Washington, Seattle.
    39. Goka K, Takafuji A, Toda S, Hamamura T, Osakabe M, Komazaki S. Genetic distinctness between two forms of Tetranychus urticae Koch (Acari: Tetranychidae) detected by electrophoresis. Exp Appl Acarol, 1996, 20: 683-693.
    40. Goka K, Takafuji A. Allozyme variations among populations of the two-spotted spider mite, Tetranychus urticae Koch, in Japan. Appl Entomol Zool, 1995, 30 (4): 567-579.
    41. Gomi K, Gotoh T, Noda H. Wolbachia having no effect on reproductive incompatibility in Tetranychus kanzawai Kishida (Acari: Tetranychidae) .Appl Entomol Zool, 1997, 32: 485-490.
    42. Gotoh T, Tokioka T. Genetic compatibility among diapausing red, non-diapausing red and diapausing green forms of the two-spotted spider mite, Tetranychus urticae Koch (Acari: Tetranychidae) . Jpn J Entomol, 1996, 64: 215-225.
    43. Guo F Y, Zhang Z Q, Zhao Z M. Pesticide resistance of Tetranychus cinnabarinus (Acari: Tetranychidae) in China: a review. Syste Appl Acarol, 1998, 3: 3-7.
    44. Hamada H, Petrino M, Kakunaga T. A novel repeated element with Z-DNA-forming potential is widely found in evolutionarily diverse eukaryotic genomes. Proc Natl Acad Sci, 1982, 79: 646-649.
    45. Hance T, Neuberg P, Noel-Lastelle C. The use of fecundity, lobe biometry and the RAPD-PCR technique in order to compare strains of Tetranychus sp. Exp Appl Acarol, 1998, 22: 649-666.
    46. Harr B, Todorova J, Schlterer C. Mismatch repair-driven mutational bias in D.mela- nogaster. Mol Cell, 2002, 10: 199-205.
    47. Hartl D L, Clark A G. Principles of population genetics 3~(rd) edn. Sinquer Associates, Sunderland, MA, 1997: 283-285.
    48. Hasting J. Thermoregulation in the dog-day cicada Tibicen duryi (Homoptera, Cicadidae). Trans. K. Acad. Sci., 1989, 50: 145-149.
    49. Hedrick P W. Genetics of population, Indiana: Second ed. Jones and Bartlett Publisher, Inc. , 2000, 205 - 207.
    50. Helle W, Van de Bund C F. Crossbreeding experiments with some species of the Tetranychus urticae group. Entomol Exp Appl, 1962, 5: 159-165.
    51. Hillis D M, Moritz C, Mable B K. Molecular systernatics, 2~(nd) edn. Sinauer, Sunderland, MA, 1996.
    52. Hinomoto N, Osakabe M, Gotoh T, Takafuji A. Phylogenetic analysis of green and red forms of the two-spotted spider mite, Tetranychus urticae Koch (Acari: Tetranychidae) , in Japan, based on mitochondrial cytochrome oxidase subunit I sequences. Appl Entomol Zool, 2001a, 36(4) : 459 -464.
    53. Hinomoto N, Takafuji A. Genetic diversity and phylogeny of the Kanzawa spider mite, Tetranychus kanzawai, in Japan. Exp Appl Acarol, 2001b, 25: 355 - 370.
    54. Hinomoto N, Takafuji A. Studies on the population structures of the two-spotted mite Tetranychus urticae Koch, by allozyme variability analysis. Appl Entomol Zool, 1994, 29: 259-266.
    55. Ho C C, Lo C C, Chen W H. Spider mite (Acari: Tetranychidae) on various crops in Taiwan. J Agr Res China, 1997, 46: 333-346.
    56. Jarne P, Lagoda P J L. Microsatellites from molecules to populations and back. Trends Ecol Evol, 1996, 11: 424-429.
    57. KattiMV, Sami-SubbuR, RanjekarPK, et al. Amino acid repeat patterns in protein sequences: their diversity and structural-functional implications. Protein Sci, 2000, 9: 1203 - 1209.
    58. Kephart S R. Starch gel electrophoresis of plant isozymes: A comparative analysis of techniques. Amer J Botany, 1990, 77 (5) : 693-712.
    59. Lan M Y, Sheng W H, et al. RAPD-PCR method for detecting DNA polymorphism in Boophilus microplus (Acari: Ixodidae). System Appl Acarol, 1996, 1: 1-4.
    60. Li S B, Lai J H, Gao S H, et al. STR polymorphism in five Chinese ethnic groups. Acta Genet Sin, 2000, 27 (12): 1035-41.
    61. Love J M, Knight A M, Mcaleer M A, et al. Towards construction of a high-resolution map of mouse genome using PCR-analyzed microsatellite. Nucleic Acids Res. 1990, 18: 4123-4130.
    62. Mixson T R, Fang Q Q, McLain D K, Oliver J H. Population structure of the blacklegged tick Ixodes scapularis revealed by SSCP data using the mitochondrial Cyt b and the nuclear ITS1 markers. Acta Zool Sinica, 2004, 50 (2) : 176-186.
    63. Moxon E R, Wills C. DNA microsatellites : Agents of evolution? Sci Am, 1999, 280 (1) : 94 -100.
    64. Nakamura Y. Variable number of tandem repeats (VNTR) markers of human gene mapping. Science, 1987, 235: 1616-1622.
    65. Navajas M, Gutierrez J, Lagnel J, Boursot P. Mitochondrial cytochrome oxidase I in tetranychid mites: a comparison between molecular phylogeny and changes of morphological and life history traits. Bull Entomol Res, 1996, 86: 407-417.
    66. Navajas M J, Lagnel J, Gutierrez J, Boursot P. Species-wide homogeneity of nuclear ribosomal ITS2 sequences in the spider mite Tetranychus urticae contrasts with extensive mitochondrial COI polymorphism. Heredity, 1998, 80: 742-752.
    67. Navajas M J, Thistlewood H M A, Lagnel J, Hushes C. Microsatellite sequences are under-represented in two mite genomes. Insect Mol Biol, 1998, 7 (3): 249 - 256.
    68. Navajas M J. Host plant associations in the spider mite Tetranychus urticae (Acari: Tetranychidae): insights from molecular phylogeography. Exp App Acarol, 1998 , 22: 201 - 214.
    69. Navajas M, Lagnel J, Gutierrez J, Boursot P. Species-wide homogeneity of nuclear ribosomal ITS2 sequences in the spider mite Tetranychus urticae contrasts with extensive mitochondrial COI polymorphism. Heredity, 1998, 80: 742-752.
    70. Navajas M, Gutierrez J, Lagnel J, Fauvel G, Gotoh T. DNA sequences and cross- breeding experiments in hawthorn spider mite Amphitetranychus viennensis reveal high genetic differentiation between Japanese and French populations. Entomol Exp Appl, 1999, 90: 113-122.
    71. Navajas MJ, Perrot-Minnot MJ. Lagnel J, Migeon A, Bourse T, Cornuet J M. Genetic structure of a greenhouse population of the spider mite Tetranychus urticae: spatio-temporal analysis with microsatellite markers. Insect Mol Biol, 2002, 11 (2): 157-165.
    72. Nei M, kumar S. Molecular evolution and phylogenetics. Oxford University Press, NewYork, 2000.
    73. Nei M. Analysis of gene diversity in subdivided population. PNAS, 1973, 70: 3321 - 3323.
    74. Nei M. Estimation of average heterozygosity and genetic distance from a small number of individuals. Genetics, 1978, 89: 583-590.
    75. Nei M. Genetic distance between populations. American Naturalist. 1972, 106: 283 - 292.
    76. Nei M. Molecular evolutionary genetics. Columbia University Press, New York, 1987.
    77. Osakabe M, Sakagami Y. RFLP analysis of ribosomal DNA in sibling species of spider mite, genus Panonychus (Acari: Tetranychidae) . Insect Mol Biol, 1994, 3: 63-66.
    78. Page R D M. TREEVIEW: An application to display phylogenetic trees on personal computers. Comp Appl Biosci, 1996, 12: 357-358.
    79. Parr W J, Hussey N W. Further studies on the reproductive isolation of geographical strains in the Tetranychus urticae complex. Entomol Exp Appl, 1960, 3: 137-141.
    80. Passos D T, da Silva S S, Richter M F, Ozaki L S. Detection of genomic variability in different populations of the cattle tick Boophilus microplus in southern Brazil. Vet Parasitol, 1999, 87: 83 -92.
    81. Perez T, Albornoz J, Dominguez A. An evaluation of RAPD fragment reproducibility and nature. Mol Ecol, 1998, 7: 1347-1357.
    82. Powell W, Morgente M, Andre C. et al. The comparison of RFLP, RAPD and SSR (microsatellie) markers for germplasm analysis. Mol Breeding, 1996, 2: 225 - 238.
    83. Raymond M, Rousset F. GENEPOP (Version 1.2) : population genetics software for exact tests and ecumenicism. J Heredity, 1995, 86: 248-249.
    84. Reed D H, Frankham R. How closely correlated are molecular and quantitative measures of genetic diversity: a meta-analysis. Evolution, 2001,55:1095-1103.
    85. Rousset F. Genetic differentiation and estimation of gene flow from F-statistics under isolation by distance. Genetics, 1997, 145: 1219-1228.
    86. Saccheri I, Kuussaari M, Kankare M, Vikman P, Fortelius W, Hanski I. Inbreeding and extinction in a butterfly metapopulation. Nature, 1998, 392 : 491-494.
    87. Salomone N, Frati F, Bernini F. Investigation on the taxonomic status of Steganacarus magnus and Steganacarus anomalus (Acari: Oribatida) using mitochondrial DNA sequences. Exp Appl Acarol, 1996, 20: 607-615.
    88. Schaal B A. Measurement of gene flow in Lupinus texensis. Nature, 1980, 284: 450-451.
    89. Schafer R, Zischler H, Birsner U, Becker A, Epplen J T. Optimized oligonucleotide probes for DNA fingerprinting. Electrophoresis, 1988, 9: 369-374.
    90. Schemske D W. Population structure and local selection in Impatiens (Balsaminaceae), a selfing annual. Evolution, 1984, 38: 817-832.
    91. Shannon C E, Weaver W. The mathematical theory of communication. University of Illinois Press, Press, Urbana, IL, 1949, 117-126.
    92. Simon J C, Baumann S, Sunnucks P, et al. Reproductive mode and population genetic structure of the cereal aphid Sitobion avenae studied using phenotypic and microsatellite markers. Mol Ecol, 1999, 8: 531-545.
    93. Skinner D M, Beatie W G, Blattner F R, et al. The Sequence of a Hermit Crab Satellite DNA in (-TAGG-)_n- (-TAGG-)_n. Biochemistry, 1974, 13: 3930-3937.
    
    94. Slakin M. A measure of population subdivision based on microsatellie allele frequencies. Genetics, 1995, 139: 457-462.
    95. Slatkin M. Gene flow and the geographic structure of natural population. Science, 1987, 236: 787 -792.
    96. Sugasawa J, Kitashima Y, Gotoh T. Hybrid affinities between the green and the red forms of the two-spotted spider mite Tetranychus urticae (Acari: Tetranychidae) under laboratory and semi-natural conditions. Appl Entomol Zool, 2002, 37: 127-139.
    97. Sunnucks P. Efficient genetic markers for population biology. Trends Ecol Evol, 2000, 15 (5) : 199-203.
    98.Takafuji A,So P M,Tsuno N.Inter- and intra-population variations in diapause attribute of the two-spotted spider mite,Tetranychus urticae Koch,in Japan.Res Popul Ecol,1991,33:331-344.
    99.Tautz D.,Renz M.Simple sequences are ubiquitous repetitive components of eukaryote genomes.Nucl Acids Res,1984,12:4217-4138.
    100.Travis C G,Ryan S O,Wolfgang S,etal.Microsatellite DNA loci for genetic studies of Cranes.Proc North A m Crane Workshop,1997,7:36-45.
    101.Tsagkarakou A,Guillemaud T,Rousset F,Navajas M.Molecular identification of a Wolbachia endosymbiont in a Tetranychus urticae strain(Acari:Tetranychidae).Insect Mol Biol,1996,5:217-221.
    102.Tsagkarakou A,Navajas M,Lagnel J,Pasteur N.Population structure in the spider mite Tetranychus urticae(Acari:Tetranychidae) from Crete based on multiple allozymes.Heredity,1997,78:84-92.
    103.Tsagkarakou A,Navajas M,Papaioannou-Souliotis P,Pasteur N.Gene flow among Tetranychus urticae(Acari:Tetranychidae) populations in Greece.Mol Ecol,1998,7:71-79.
    104.Tsagkarakou A,Navajas M,Rousset F,Pasteur N.Genetic differentiation in Tetranychus urticae (Acari:Tetranychidae) from greenhouses in France.Exp Appl Acarol,1999,23:365-378.
    105.Van de Bund C F,Helle W.Investigations on the Tetranychus urticae complex in north west Europe (Acari:Tetranychidae).Entomol Exp Appl,1960,3:142-156.
    106.Vanoverbeke J,De Meester L.Among-population genetic differentiation in the cyclical parthenogen Daphnia magna(Crustacea,Anomopoda) and it relation to geographic distance and clonal diversity.Hydrobiologia,1997,360:135-142.
    107.Veerman A.The pigments ofTetranychus cinnabarinus(Acari:Tetranychidae).Comp Biochem Physiol,1970,36:749-763.
    108.Wang Z,Weber J L,Zhong G,et al.Survey of plant short tandem DNA repeats.Nature,1994,359:794-801.
    109.Weber J L,May P E.Abundant class of human DNA polymorphisms which can be typed using the polymerase chain reaction.Am J Hum Genet,1989,44:388-396.
    110.Weber J L,Wong C.Mutation of human short tandem repeats.Hum Mol Genet,1993,(2):1123-1128.
    111.Weeks A R,Van Opijnen T,Breeuwer J A J.AFLP fingerprinting for assessing intraspecific variation and genome mapping in mites.Exp Appl Acarol,2000,24:775-793.
    112.Weir B S,Cokerham C C.Estimating F-statistics for the analysis of population structure. Evolution, 1984, 38: 1358 - 1370.
    113. Wesson D M. McLain D K , Oliver J H, Piesman J, Collins F H. Investigation of the validity of species status of Ixodes dammini (Acari: ixodidae) using rDNA. Proc Natl Acad Sci. USA, 1993, 90: 10221-10225.
    114. Wren J D, Forgacs E, Fondon J W, et al. Garner Repeat polymorphisms within gene regions: phenotypic and evolutionary implications. Am J Hum Genet, 2000, 67: 345-356.
    115. Wright S. Eolution in Mendelian population. Genetics, 1931, 16: 97-159.
    116. Wright S. Evolution and the Genetics of Population, 1969, Vol.2, The Thory of Gene Frequencies. Univ Chicago Press, Chicago, IL.
    117. Wright S. Evolution and the genetics of populations, 1978, vol.4. Variability within and among natural populations. Univ Chicago Press, Chicago, IL.
    118. Wright S. The genetical structure of population. Ann. Eugen., 1951, 15: 323-354.
    119. Xie L, Hong X Y, Xue X F. Population genetic structure of the two-spotted spider mite (Acari: Tetranychidae) from China. Ann Entomol Soc Amer, 2006, 99: 959-965.
    120. Yeh F C, Yang R C, Boyle T. Microsoft window-based freeware for population genetic analysis (POPGENE), version 1.31. Genetics, 1999, 147: 245-251.
    121. Young S S Y, Wrensch D L, Kongchuensin M. Geographic variations and combining abilities in the two-spotted spider mite, Tetranychus urticae. Entomol Exp Appl, 1985, 39: 109-113.
    122. Zhang Z Q, Jacobson R J. Using adult female morphological characters for differentiating Tetranychus urticae complex (Acari: Tetranychidae) from greenhouse tomato crops in UK. Syst Appl Acarol, 2000, 5: 69-76.