- Research article
- Open Access
Similar patterns of linkage disequilibrium and nucleotide diversity in native and introduced populations of the pea aphid, Acyrthosiphon pisum
© Brisson et al; licensee BioMed Central Ltd. 2009
- Received: 27 June 2008
- Accepted: 26 May 2009
- Published: 26 May 2009
The pea aphid, Acyrthosiphon pisum, is an emerging genomic model system for studies of polyphenisms, bacterial symbioses, host-plant specialization, and the vectoring of plant viruses. Here we provide estimates of nucleotide diversity and linkage disequilibrium (LD) in native (European) and introduced (United States) populations of the pea aphid. Because introductions can cause population bottlenecks, we hypothesized that U.S. populations harbor lower levels of nucleotide diversity and higher levels of LD than native populations.
We sampled four non-coding loci from 24 unique aphid clones from the U. S. (12 from New York and 12 from California) and 24 clones from Europe (12 alfalfa and 12 clover specialists). For each locus, we sequenced approximately 1 kb from two amplicons spaced ~10 kb apart to estimate both short range and longer range LD. We sequenced over 250 kb in total. Nucleotide diversity averaged 0.6% across all loci and all populations. LD decayed slowly within ~1 kb but reached much lower levels over ~10 kb. Contrary to our expectations, neither LD nor nucleotide diversity were significantly different between native and introduced populations.
Both introduced and native populations of pea aphids exhibit low levels of nucleotide diversity and moderate levels of LD. The introduction of pea aphids to North America has not led to a detectable reduction of nucleotide diversity or increase in LD relative to native populations.
- Linkage Disequilibrium
- Nucleotide Diversity
- Effective Population Size
- Linkage Disequilibrium Decay
- Acyrthosiphon Pisum
One important goal of evolutionary biology is to identify the genetic causes of phenotypic variation in wild populations. In theory, in large, randomly mating populations where all loci are in linkage equilibrium, genotype can be connected to phenotype by assaying every DNA polymorphism and associating those polymorphisms with phenotypic variation . In practice, loci are often in linkage disequilibrium (LD) and not all polymorphisms can be surveyed due to cost. LD can be exploited, however, to identify genomic regions, rather than individual nucleotides, that are associated with phenotypic variation [e.g., ]. The scale at which polymorphisms must be sampled depends on the level and decay of LD across the genome. Higher levels of LD allow shallower sampling of polymorphisms, but provide reduced resolution for locating causal variants. Lower levels of LD force denser sampling of polymorphisms, but provide greater precision for linking genotype to phenotype.
Systematic studies of LD have been performed for few species and mainly for plants. Among plants, LD varies widely, with higher levels observed in predominately selfing organisms such as Arabidopsis thaliana [3, 4] and soybeans . Selfing reduces effective heterozygosity and therefore the ability of recombination to reduce LD. Among outcrossers such as sunflowers  and maize [7, 8], LD decays more rapidly, reaching background levels within one kilobase. In humans LD persists from tens to hundreds of kilobases [2, 9, 10]. In a flycatcher (Ficedula albicollis), significant LD extends over 400 to 500 kb . In Caenorhabditis elegans LD generally decays to insignificant levels within a kilobase [12, 13] and in Drosophila melanogaster within even shorter distances .
Here we characterize levels of LD and nucleotide diversity in the pea aphid (Acyrthosiphon pisum). Pea aphids, like most aphids, are cyclical parthenogens . During the spring and summer months, females are parthenogenetic, giving live birth to typically genetically identical daughters via a modified meiosis lacking genetic recombination  (but see  for a recent review of the possibility of genetic variation within an aphid clone). A lineage of asexual females is typically referred to as a clone. Pea aphid males and sexual females are produced asexually in the fall in response to shortened day-length and lower temperatures [18, 19]. Eggs produced by the sexual generation diapause through the cold winter months when adults perish.
The pea aphid genome has been sequenced at 6× coverage by the Human Genome Sequencing Center at Baylor College of Medicine and an assembly has been released http://www.hgsc.bcm.tmc.edu/projects/aphid. Notable aspects of the pea aphid's biology include switching seasonally between asexual and sexual reproduction , harboring obligate and facultative bacterial endosymbionts [20, 21], producing winged and unwinged morphs in response to environmental cues , displaying strong preferences for different host plants , transmitting plant viral diseases , and producing overwintering, diapausing eggs . Many of these traits exhibit variable phenotypes that can be explored using a mapping approach, which would be informed by estimates of levels of nucleotide diversity and LD in this species.
Pea aphids infest field crops of the pea family such as alfalfa (genus Medicago) and clover (genus Trifolium). Although Eurasian in origin, the pea aphid has been introduced to North America, presumably via agriculture, probably many times over the last 200 years . Severe bottleneck events, such as those associated with the introduction of only a few individuals, are expected to dramatically decrease levels of nucleotide diversity and increase levels of LD [27, 28]. Less severe bottlenecks decrease the number of low-frequency single-nucleotide polymorphisms (SNPs) in the introduced population . We therefore hypothesized that introduced populations harbor lower levels of nucleotide diversity and higher levels of LD than native populations. If true, the two populations would have different utility for the coarse-mapping (introduced) and fine-mapping (native) of phenotypic variation.
Aphid clones display increased survivorship and fecundity on their preferred host plants. In New York, pea aphids specialize on alfalfa (Medicago sativa) or red clover (Trifolium pratense) [23, 30] and this specialization has a genetic basis . Specialization on alfalfa, clover and other species of the pea family has also been observed in France , Sweden , and the United Kingdom . In California, pea aphids specialize on white clover (Trifolium repens), but no alfalfa specialists have been reported. Clones found on alfalfa are more generalized in that they are able to feed on either clover or alfalfa . We sequenced four loci from pea aphid clones from both alfalfa and clover in Europe and alfalfa in the U. S. to test whether host plant specialization or population origin was associated with population differentiation at the sampled loci and with different patterns of LD or nucleotide diversity.
Twelve pea aphid clones were collected from alfalfa fields in Yolo and Galt counties of central California by T. Leonardo and S. Nuzhdin . Twelve clones were collected from alfalfa fields in Western and Central France (Domagné and Lusignan) by J. C. Simon. Two clones were collected from clover fields in Germany (coordinates 50°54' W 11°33'N and 50°59' W 11°31'N) and ten from clover fields in France (coordinates 48°06' W 1°47'N, as well as Domagné and Lusignan) by J. C. Simon. Finally, 12 New York clones were collected from alfalfa fields in Tompkins county by M. Caillaud and Otsego county by J. Brisson. From each clone, we extracted DNA from two to six females using the Qiagen (Valencia, CA) DNAeasy kit.
We assembled four contiguous segments of 10 to 20 kb of genomic DNA using sequences from the NCBI trace archive. Two of these loci, CAA and 380, are located on the X chromosome while the other two, 870 and 41, are on autosomes. Pea aphids have three autosomes and one large X chromosome that appears to contain about a third of the nuclear genome content . Small sequence subsets of these 10 to 20 kb segments were previously identified for the CAA, 380 and 870 loci [31, 36], while 41 was developed in this study.
Primers used for PCR and sequencing.
Analyses of LD were performed using the Genetics package  in the R statistical computing environment . Our data are genotypic, not haplotypic data. Thus, for two loci, A and B, AB/ab is indistinguishable from Ab/aB. The frequency of AB, pAB, in this package is estimated via maximum likelihood, and then used to compute the difference in expected versus observed allele pairs, D, where D = pAB - pApB and the correlation coefficient between the pairs, r, where r = -D/(√(pApapBpb)). We thus report LD values between pairs of sites using the squared allele frequency correlation measure, r2. Rare sites (where only one individual differed from the others) were excluded prior to computations. Significance of LD between SNPs was assessed using Fisher's exact test, also implemented using the Genetics package. To estimate the decay of LD over distance, we used nonlinear regression of LD between polymorphic sites versus physical distance . The expected value of r2 using a drift-recombination model is E(r2) = 1/(1 + C), where C equals the population recombination parameter (C = 4Ner). Adding mutation to the model, the expectation becomes:
E(r2) = ((10 + C)/(2+C)(11+C))/(1+(3+C)(12+12C+C2)/(n(2+C)(11+C)),  where n is the number of sequences sampled. We fit this model after replacing C by C*distance in base pairs between sites using PROC NLIN in SAS .
We calculated π, θ, and Tajima's D using R. π measures nucleotide diversity as the average number of nucleotide differences per site between all pairwise comparisons of sequences and θ is the level of heterozygosity determined by the number of variable positions [41, 42]. Tajima's D is a measure of whether the frequency distribution of segregating sites at a locus is consistent with neutral expectations . Levels of population differentiation (Fst) for each of the four loci were estimated using Genepop , which uses a weighted analysis of variance . We also used Genepop to determine whether any of the loci were not in Hardy-Weinberg equilibrium, testing the null hypothesis of random union of gametes using the exact Hardy-Weinberg test .
To measure levels of sequence diversity and LD decay over short distances (less than 1 kb) as well as over longer distances (about 10 kb), we sequenced two amplicons spaced by about 10 kb at four different genomic regions for 12 clones from New York, 12 clones from California, and 24 clones from Europe (22 clones from France and two clones from Germany). The European clones consisted of 12 alfalfa and 12 clover specialists. Two of the genomic regions were from autosomes (41 and 870) and two were from the X chromosome (380 and CAA). A total of approximately 5452 bp was sequenced for each individual, with over 250 kb sequenced in total. None of the regions included coding sequence and all were in intergenic regions as determined by querying the "nr" (non-redundant) database with a blastX (translated versus protein database) search of their associated genomic contigs. We excluded all indel information from our analyses because we used diploid individuals and could not precisely assay the number or size of indels in some of our individuals. Heterozygous indels in some cases led to missing sequence data (see Figure 1).
Measures of nucleotide diversity, Tajima's D, and LD.
Mean r2 short distances
Mean r2 longer distances
All loci possessed indels. We observed two 1 bp indels in locus 41. Locus 870 had a 1 bp, a 2 bp, a 3 bp and a 5 bp indel. Locus CAA showed two 1 bp, a 4 bp and a 6 bp indel. Finally, locus 380 had five 1 bp indels, a 4 bp indel, and a 14 bp indel. There were additional indels that we did not directly observe: multiple individuals had sequencing problems at the same region likely due to a heterozygous indel. These can be inferred from Figure 1.
Tajima's D and Hardy-Weinberg Estimates
Under neutrality, the values of π (based on the number of pairwise differences between individuals) and θ (based on the number of segregating sites) should be equal [42, 47]. Tajima's D test compares these two measures of diversity. D values that are significantly different from zero indicate non-neutral evolution resulting from demographic changes or from selection . We calculated Tajima's D for each locus and each population and observed that across all populations and all loci, with only one exception (locus CAA, California), Tajima's D values were positive (Table 2). Five of 16 values were significantly positive at the 95% level. All loci had at least one significant positive D value. Four of the five significant D values were from European populations.
New York and California populations were in Hardy-Weinberg equilibrium for all loci except the autosomal locus 870 (New York: X2 = 94.2, P = 0.001, California: X2 = 84.2, P = 0.014). Both deviations at locus 870 were towards a deficiency of heterozygotes. All four loci were in Hardy-Weinberg equilibrium within the two European populations.
Fst values for each pairwise population comparison.
We also examined levels of LD across loci. Between the two autosomal loci, 870 and 41 (it is unknown whether they are on the same chromosome), the average r2 was 0.15. Loci CAA and 380 on the X chromosome had an interlocus average r2 of 0.18.
Sequence Diversity and Linkage Disequilibrium
We have investigated the levels of nucleotide diversity and LD across four genomic regions in two introduced pea aphid populations, New York and California, and two host races from native populations in Europe. This is the first estimate of short and longer range LD decay for this species. Previous studies using primarily mtDNA sequence variation, RFLPs, and allozymes have reported low diversity in pea aphid populations [32, 48–50]. Like these studies, we find low diversity, with π averaging 0.0056 across all loci and all four populations, ranging from 0.0008 (locus CAA, California) to 0.0116 (locus 41, Europe, clover).
The two loci on the X chromosome exhibit roughly half the diversity of the two loci on the autosomes. These lower X chromosome values may result, at least in part, from the lower effective population size for the X chromosome. In a population with equal numbers of males and females there are three X chromosomes for every four autosomes because females are XX and males are XO. Under neutrality, nucleotide diversity is directly proportional to effective population size . However, these X chromosome π values remained lower than the autosomal values even after multiplying by 4/3 to compensate for the X to autosomal population size differences. If this pattern remains consistent across autosomal and X chromosome loci, then this implies that other factors, such as selection, may act differentially on the X chromosome and autosomes .
In a previous study , we sequenced portions of the coding regions from 27 different loci from an average of ten individuals from disparate geographical locations. Of these 27 loci, 16 exhibited silent site variation. The average diversity of this silent site polymorphism was 0.0053 (ranging from 0.0006 to 0.0114), very similar to the average of 0.0056 we report here for noncoding variation despite low levels of codon bias in the pea aphid .
The levels of noncoding DNA diversity in the pea aphid are two and a half times lower than the diversity estimates of 0.0153 for intergenic regions of Drosophila simulans  and 0.0142 for noncoding regions of Aedes aegypti . Although it is possible that the pea aphid has a lower mutation rate, the more likely explanation for this comparatively low diversity is small effective population size. The pea aphid has a large absolute population size during its parthenogenetic phase in the spring and summer months, but many individuals are genetic clones and thus effective population size is much lower. Biased sex ratios , high variance in male mating success  and low survival of eggs over the winter might cause additional reductions in effective population size. Finally, as pea aphids are pests on field crops, it is likely that many populations have experienced declines due to pesticides, further reducing the effective population size.
We also examined levels of LD at two scales: within a kilobase, and across approximately 10 kb. Generally, LD is high and decays slowly within a kilobase (Fig. 1). Over ~10 kb, LD has declined to levels only slightly higher than LD averages across loci and therefore what we consider background levels of LD. Thus, this limited survey provides no evidence for long range LD in any of the four populations surveyed. This suggests that LD mapping in the pea aphid will require a relatively high density of molecular markers distributed at the scale of about a kilobase.
Native Versus Introduced Populations
Pea aphids were introduced to North America, probably from Europe . They are mild pests on field crops such as alfalfa and clover, and have likely been introduced to North America several times over the last 100 to 200 years via agricultural routes. If only a few individuals were introduced at a time, the North American populations would likely exhibit decreased levels of variation relative to native populations . We therefore hypothesized that we would observe lower levels of variation in North American populations relative to European populations. This does not appear to be the case, as both π and θ are similar between the two regions. Likewise, LD levels are similar in populations from both regions. We conclude that native and introduced populations will not have differential utility in LD mapping efforts.
It is unknown how many times or where the pea aphid was introduced into North America. It is possible that repeated introductions have occurred, supplementing the molecular diversity provided by earlier introductions. Consistent with this model, the introduced populations are not depauperate for nucleotide diversity. Further, Tajima's D values are largely positive which can be indicative of admixture, a reduction in population size, or balancing selection . It seems unlikely that balancing selection is acting on all four of these non-coding loci. We therefore favor the hypothesis that the positive Tajima D values reflect demographic factors.
Host plant specialization
Pea aphids are specialized to a variety of host plants in Europe, including alfalfa and clover [32, 33]. Similarly, pea aphids in New York demonstrate strong host plant specialization [23, 30], whereas in California both specialists and generalists exist . The specialization in New York either represents introduction of both specialists from the native range, or introduction followed by a secondary event of specialization. All New York individuals were sampled from alfalfa. Fst values between New York and European alfalfa and between New York and European clover populations are similar (average of 0.199 and 0.186, respectively), providing no evidence for the source of the New York populations. Both of these values are higher than for the comparison of European clover and alfalfa specialists (Fst = 0.155). This value is only slightly higher than the level of Fst reported by Simon et al.  for comparisons of clover and alfalfa specialists in France either by allozymes (Fst = 0.098) or microsatellites (Fst = 0.104).
The California pea aphids examined here were sampled from alfalfa, but all are generalists in that they can feed on both hosts (T. Leonardo, pers. comm). Interestingly, population differentiation is high between California and the European clover clones (Fst = 0.228), while lower in comparison with the European alfalfa clones (Fst = 0.168). Although this difference is not statistically significant (one-tailed t-test, P = 0.19) the trend of lower differentiation between California and the alfalfa specialists leads us to hypothesize that the host plant generalist aphids found in California are derived from European alfalfa specialists, but that subsequently there has been interbreeding with other introduced populations. Further tests of this hypothesis will require more extensive population sampling and marker development but will provide insight into the interplay of invasion of a new habitat and the gain or loss of host-plant specialization.
The recent release of the of the pea aphid genome sequence will facilitate mapping efforts aimed at identifying loci underlying variable phenotypes unique to this species, such as polyphenisms and host plant specialization. To inform these future studies, here we provide the first systematic study of nucleotide diversity and LD in both native and introduced populations of the pea aphid. Further, we tested and rejected the hypothesis that bottlenecks associated with the importation of pea aphids into North America have resulted in introduced and native populations being valuable for different mapping needs (coarse versus fine-scale mapping, respectively), finding low levels of nucleotide diversity and moderate levels of LD in all of the populations sampled.
We thank J.-C. Simon, T. Leonardo, and C. Braendle for aphid clones and P. Y. Hui for laboratory assistance. Thanks to R. Bickel, T. Leonardo, A. Tarone, T. Turner and two anonymous reviewers for thoughtful comments on the manuscript. The project was supported by Award Number F32GM076962 from the National Institute of General Medical Sciences to JAB and by grants from National Institutes of Health (GM063622-06A1) and National Science Foundation (IOS-0640339) to DLS and a David & Lucile Packard Foundation Fellowship to DLS.
- Long AD, Langley CH: The power of association studies to detect the contribution of candidate genetic loci to variation in complex traits. Genome Res. 1999, 9: 720-731.PubMed CentralPubMedGoogle Scholar
- Gibbs RA, Belmont JW, Hardenbol P, Willis TD, Yu FL, Yang HM, Ch'ang LY, Huang W, Liu B, Shen Y, et al: The International HapMap Project. Nature. 2003, 426: 789-796. 10.1038/nature02168.View ArticleGoogle Scholar
- Nordborg M, Borevitz JO, Bergelsen J, Berry CC, Chory J, Hagenblad J, Kreitman M, Maloof JN, Noyes T, Oefner PJ, et al: The extent of linkage disequilibirum in Arabidopsis thaliana. Nature Genet. 2002, 30: 190-193. 10.1038/ng813.View ArticlePubMedGoogle Scholar
- Nordborg M, Hu TT, Ishino Y, Jhaveri J, Toomajian C, Zheng H, Bakker E, Calabrese P, Gladstone J, Goyal R, et al: The pattern of polymorphism in Arabidopsis thaliana. PLoS Biol. 2005, 3: e196-10.1371/journal.pbio.0030196.PubMed CentralView ArticlePubMedGoogle Scholar
- Zhu YL, Song QJ, Hyten DL, Van Tassell CP, Matukumalli LK, Grimm DR, Hyatt SM, Fickus EW, Young ND, Cregan PB: Single-nucleotide polymorphisms in soybean. Genetics. 2003, 163 (3): 1123-1134.PubMed CentralPubMedGoogle Scholar
- Liu A, Burke JM: Patterns of nucleotide diversity in wild and cultivated sunflower. Genetics. 2006, 173 (1): 321-330. 10.1534/genetics.105.051110.PubMed CentralView ArticlePubMedGoogle Scholar
- Remington DL, Thornsberry JM, Matsuoka Y, Wilson LM, Whitt SR, Doebley J, Kresovich S, Goodman MM, Buckler ESt: Structure of linkage disequilibrium and phenotypic associations in the maize genome. Proc Natl Acad Sci USA. 2001, 98 (20): 11479-11484. 10.1073/pnas.201394398.PubMed CentralView ArticlePubMedGoogle Scholar
- Tenaillon MI, Sawkins MC, Long AD, Gaut RL, Doebley JF, Gaut BS: Patterns of DNA sequence polymorphism along chromosome 1 of maize (Zea mays ssp. mays L.). Proc Natl Acad Sci USA. 2001, 98 (16): 9161-9166. 10.1073/pnas.151244298.PubMed CentralView ArticlePubMedGoogle Scholar
- Reich DE, Cargill M, Bolk S, Ireland J, Sabeti PC, Richter DJ, Lavery T, Kouyoumjian R, Farhadian SF, Ward R, et al: Linkage disequilibrium in the human genome. Nature. 2001, 411 (6834): 199-204. 10.1038/35075590.View ArticlePubMedGoogle Scholar
- Gabriel SB, Schaffner SF, Nguyen H, Moore JM, Roy J, Blumenstiel B, Higgins J, DeFelice M, Lochner A, Faggart M, et al: The structure of haplotype blocks in the human genome. Science. 2002, 296 (5576): 2225-2229. 10.1126/science.1069424.View ArticlePubMedGoogle Scholar
- Backstrom N, Qvarnstrom A, Gustafsson L, Ellegren H: Levels of linkage disequilibrium in a wild bird population. Biology Letters. 2006, 2: 435-438. 10.1098/rsbl.2006.0507.PubMed CentralView ArticlePubMedGoogle Scholar
- Cutter AD, Baird SE, Charlesworth D: High nucleotide polymorphism and rapid decay of linkage disequilibrium in wild populations of Caenorhabditis remanei. Genetics. 2006, 174: 901-913. 10.1534/genetics.106.061879.PubMed CentralView ArticlePubMedGoogle Scholar
- Cutter AD: Nucleotide plymorphism and linkage disequilibrium in wild populations of the partial selfer Caenorhabditis elegans. Genetics. 2006, 172: 171-184. 10.1534/genetics.105.048207.PubMed CentralView ArticlePubMedGoogle Scholar
- Andolfatto P, Przeworski M: A genome-wide departure from the standard neutral model in natural populations of Drosophila. Genetics. 2000, 156: 257-268.PubMed CentralPubMedGoogle Scholar
- Moran NA: The evolution of aphid life cycles. Ann Rev Ent. 1992, 37: 321-348. 10.1146/annurev.en.37.010192.001541.View ArticleGoogle Scholar
- Blackman RL: Reproduction, cytogenetics and development. Aphids: Their biology, natural enemies & control. Edited by: Minks AK, Harrewijn P. 1987, Amsterdam: Elsevier, 2A: 163-195.Google Scholar
- Loxdale HD: The nature and reality of the aphid clone: genetic variation, adaptation and evolution. Agricult Forest Ent. 2008, 10: 81-90. 10.1111/j.1461-9563.2008.00364.x.View ArticleGoogle Scholar
- Via S: Inducing the sexual forms and hatching the eggs of pea aphids. Ent Exp Appl. 1992, 65: 119-127. 10.1007/BF00221362.View ArticleGoogle Scholar
- MacKay PA, Reeleder DJ, Lamp RJ: Sexual morph production by apterous and alate viviparous Acyrthosiphon pisum (Harris) (Homoptera: Aphididae). Can J Zool. 1983, 61: 952-957.View ArticleGoogle Scholar
- Buchner P: Endosymbiosis of animals with plant microorganisms. 1965, New York: Interscience Publishers, IncGoogle Scholar
- Baumann P, Moran NA, Baumann L: The evolution and genetics of aphid endosymbionts. Bio Science. 1997, 47: 12-20.Google Scholar
- Sutherland ORW: The role of crowding in the production of winged forms by two strains of the pea aphid, Acyrthosiphon pisum. J Insect Phys. 1969, 15: 1385-1410. 10.1016/0022-1910(69)90199-1.View ArticleGoogle Scholar
- Via S: Reproductive isolation between sympatric races of pea aphids. I. Gene flow restriction and habitat choice. Evolution. 1999, 53: 1446-1457. 10.2307/2640891.View ArticleGoogle Scholar
- Nault LR: Arthropod transmission of plant viruses: a new synthesis. Ann Ent Soc Am. 1997, 90: 521-541.View ArticleGoogle Scholar
- Shingleton AW, Sisk GC, Stern DL: Diapause in the pea aphid (Acyrthosiphon pisum) is a slowing but not a cessation of development. BMC Dev Biol. 2003, 3: 7-10.1186/1471-213X-3-7.PubMed CentralView ArticlePubMedGoogle Scholar
- Thomas C: A list of the species of the tribe Aphidini, family Aphidae, found in the United States, which have been heretofore named with descriptions of some new species. State Lab Nat Hist Bull. 1878, 1: 1-16.Google Scholar
- Nei M, Maruyama T, Chakraborty R: The bottleneck effect and genetic variability in populations. Evolution. 1975, 29: 1-10. 10.2307/2407137.View ArticleGoogle Scholar
- Crow JF, Kimura M: An Introduction to Population Genetics Theory. 1970, New York: Harper and RowGoogle Scholar
- Maruyama T, Fuerst PA: Population bottlenecks and nonequilibrium models in population genetics. II. Number of alleles in a small population that was formed by a recent bottleneck. Genetics. 1985, 11: 675-689.Google Scholar
- Via S, Bouck AC, Skillman S: Reproductive isolation between divergent races of pea aphids on two hosts. II. Selection against migrants and hybrids in the parental environments. Evolution. 2000, 54: 1626-1637.View ArticlePubMedGoogle Scholar
- Hawthorne DJ, Via S: Genetic linkage of ecological specialization and reproductive isolation in pea aphids. Nature. 2001, 412: 904-907. 10.1038/35091062.View ArticlePubMedGoogle Scholar
- Simon JC, Carre S, Boutin M, Prunier-Leterme N, Sabater-Munoz B, Latorre A, Bournoville R: Host-based divergence in populations of the pea aphid: insights from nuclear markers and the prevalence of facultative symbionts. Proc R Soc Lond Ser B. 2003, 270: 1703-1712. 10.1098/rspb.2003.2430.View ArticleGoogle Scholar
- Sandstrom J: Temporal changes in host adaptation in the pea aphid, Acyrthosiphon pisum. Ecol Entomol. 1996, 21: 52-56. 10.1111/j.1365-2311.1996.tb00266.x.View ArticleGoogle Scholar
- Ferrari J, Godfray CJ, Faulconbridge AS, Prior K, Via S: Population differentiation and genetic variation in host choice among pea aphids from eight host plant genera. Evolution. 2006, 60: 1574-1584.View ArticlePubMedGoogle Scholar
- Leonardo TE, Muiru GT: Facultative symbionts are associated with host plant specialization in pea aphid populations. Proceedings: Biological Sciences. 2003, 270: S209-S212. 10.1098/rsbl.2003.0064.PubMed CentralGoogle Scholar
- Braendle C, Caillaud MC, Stern DL: Genetic mapping of aphicarus – a sex-linked locus controlling a wing polymorphism in the pea aphid (Acyrthosiphon pisum). Heredity. 2005, 94: 435-442. 10.1038/sj.hdy.6800633.View ArticlePubMedGoogle Scholar
- Warnes GR: The genetics package. R News. 2003, 3: 9-13.Google Scholar
- Team RDC: R: A language and environment for statistical computing. 2006, Vienna, Austria: R Foundation for Statistical Computing, 2.3.1Google Scholar
- Hill WG, Weir BS: Variances and covariances of squared linkage disequilibria in finite populations. Theor Popul Biol. 1988, 33: 54-78. 10.1016/0040-5809(88)90004-4.View ArticlePubMedGoogle Scholar
- Inc. SI: SAS/STAT user's guide. 1989, Cary, NC: SAS Institute, Inc, 2:Google Scholar
- Nei M: Molecular Evolutionary Genetics. 1987, New York: Columbia University PressGoogle Scholar
- Watterson GA: On the number of segregating sites in genetical models without recombination. Theor Popul Biol. 1975, 7: 256-275. 10.1016/0040-5809(75)90020-9.View ArticlePubMedGoogle Scholar
- Tajima F: Statistical method for testing the neutral mutation hypothesis by DNA polymorphism. Genetics. 1989, 123 (3): 585-595.PubMed CentralPubMedGoogle Scholar
- Raymond M, Rousset F: GENEPOP (version 1.2): population genetics software for exact tests and ecumenicism. J Heredity. 1995, 86: 248-249.Google Scholar
- Weir BC, Cockerham CC: Estimating F-statistics for the analysis of population structure. Evolution. 1984, 38: 1358-1370. 10.2307/2408641.View ArticleGoogle Scholar
- Guo SW, Thompson EA: Performing the exact test of Hardy-Weinberg proportion for multiple alleles. Biometrics. 1992, 48: 359-10.2307/2532296.View ArticleGoogle Scholar
- Tajima F: Evolutionary relationship of DNA sequences in finite populations. Genetics. 1983, 105: 437-460.PubMed CentralPubMedGoogle Scholar
- Barrette RJ, Crease TJ, Hebert PDN, Via S: Mitochondrial DNA diversity in the pea aphid, Acyrthosiphon pisum. Genome Biology. 1994, 37: 858-865.View ArticleGoogle Scholar
- Boulding EG: Molecular evidence against phylogenetically distinct host races of the pea aphid (Acyrthosiphon pisum). Genome Biology. 1998, 41: 769-775.View ArticleGoogle Scholar
- Birkle LM, Douglas AE: Low genetic diversity among pea aphid (Acyrthosiphon pisum) biotypes of different plant affiliation. Heredity. 1999, 82: 605-612. 10.1046/j.1365-2540.1999.00509.x.View ArticlePubMedGoogle Scholar
- Hartl DL, Clark AG: Principles of Population Genetics. 1997, Sinauer AssociatesGoogle Scholar
- Begun DJ, Whitley P: Reduced X-linked nucleotide polymorphism in Drosophila simulans. Proc Nat Acad Sci USA. 2000, 97: 5960-5965. 10.1073/pnas.97.11.5960.PubMed CentralView ArticlePubMedGoogle Scholar
- Brisson JA, Nuzhdin SV: Rarity of males in pea aphids results in mutational decay. Science. 2008, 319: 58-10.1126/science.1147919.View ArticlePubMedGoogle Scholar
- Rispe C, Legeai F, Gauthier J-P, Tagu D: Strong heterogeneity in nucleotide composition and codon bias in the pea aphid (Acyrthosiphon pisum) shown by EST-based coding genome reconstruction. J Mol Evol. 2007, 65: 413-424. 10.1007/s00239-007-9023-y.View ArticlePubMedGoogle Scholar
- Begun DJ, Holloway AK, Stevens K, Hilllier LW, Poh Y, Hahn MW, Nista PM, Jones CD, Kern AD, Dewey CN, et al: Population genomics: whole-genome analysis of polymorphism and divergence in Drosophila simulans. PLoS Biol. 2007, 5: e310-10.1371/journal.pbio.0050310.PubMed CentralView ArticlePubMedGoogle Scholar
- Morlais I, Severson DW: Intraspecific DNA variation in nuclear genes of the mosquito Aedes aegypti. Insect Mol Biol. 2003, 12: 631-639. 10.1046/j.1365-2583.2003.00449.x.View ArticlePubMedGoogle Scholar
- Foster WA: Aphid sex ratios. Sex Ratio Handbook: concepts and research methods. Edited by: Hardy I. 2002, Cambridge University Press, 254-265.View ArticleGoogle Scholar
- Sack C, Stern DL: Sex and death in the male pea aphid: Acyrthosiphon pisum: the life-history effects of a wing dimorphism. J Insect Phys. 2007, 7: 45-Google Scholar
- Mayr E: Animal Species and Evolution. 1963, Cambridge: Harvard University PressView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.