- Research article
- Open Access
Similar patterns of linkage disequilibrium and nucleotide diversity in native and introduced populations of the pea aphid, Acyrthosiphon pisum
BMC Geneticsvolume 10, Article number: 22 (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.
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.
We wanted to assay LD within short distances (less than 1 kb) and longer distances (over about 10 kb), so we designed two sets of PCR primers for each of the four loci. The resulting amplicons from each set were spaced by about 10 kb. Because we had two amplicons per locus and four loci, in total we utilized eight amplicons (from each of the 48 clones). Specifically, our amplicons (sized between 500 and 1000 bp) were separated by 9249 bp for locus 380, by 11371 bp for locus 870, by 9163 bp for locus 41, and by 10684 bp for locus CAA. After release of the first assembly of the pea aphid genome, Acyr 1.0, we confirmed these distances, with locus 380 on contig 7746 in scaffold 5539, 870 on contig 14350 in scaffold 14350, 41 on contig 6316 in scaffold 12009 and CAA on contig 18522 in scaffold 15751. We also confirmed that all four of the loci are in intergenic regions. Primers are listed in Table 1. 50 μl PCR reactions included 25 μl of 2 × PCR Master Mix (Fermentas, Glen Burnie, MD), 5 μl of each 10 μM primer, 2 μl of genomic DNA, and 13 μl H2O. All PCR reactions were done with negative controls under the following conditions: an initial 95°C anneal, 30 cycles of 95°C for 30 seconds, 55°C for 30 seconds and 72°C for 1 minute, with a final 10 minute extension at 72°C. PCR products were purified using the Invitrogen (Carlsbad, CA) Purelink PCR purification kit and sequenced directly from both ends as well as with the internal primers listed in Table 1.
Sequences were assembled using CodonCode Aligner (CodonCode Corporation, Dedham, MA). Because these were natural pea aphid isolates, all loci had heterozygous positions and some had heterozygous insertions or deletions (indels). The software identified heterozygous positions that were confirmed by visual inspection of the sequence trace file. In some cases we did not sequence across the entire amplicon due to the presence of heterozygous indels. These areas were treated as unknown data. There were also instances where we were unable to amplify a locus in a given genotype. Figure 1 illustrates the gaps in our data for each of the four loci. All indels were excluded from analyses of nucleotide diversity. Sequences have been deposited in GenBank under accession numbers FJ825706–FJ825751 and FJ858381–FJ858708.
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).
Table 2 reports estimates of nucleotide diversity (π) and Watterson's θ. Nucleotide diversity was significantly lower on the X chromosome than the autosomes (1-sided P < 0.001, unpaired t-test comparing the eight autosomal π values to the eight X π values; if we increase the sample size contributing to nucleotide diversity estimates by calculating a π value for the North American samples and a π value for the European samples for each locus and then compare the four autosomal π values to the four X π values, 1-sided P < 0.01). Further, measures of nucleotide diversity were not significantly different between North American and European populations (NY and CA π values averaged versus Europe, alfalfa and clover averaged: 1-sided P = 0.12, paired t-test; again, if we calculate a π value for the North American samples and a π value for the European samples for each locus and compare them, 1-sided P = 0.09, paired t-test).
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 indicated moderate differentiation between most populations, with an average value of 0.188 (Table 3). The lowest value was observed for the comparison of European alfalfa and European clover populations (0.155) and the highest for the comparison of California and clover-feeding populations from Europe (0.228).
We estimated LD levels, as measured by the correlation between variable bases (r2), separately for each population and each locus. The average correlation at both short-range (less than 1 kb) and long-range (about 10 kb) distances is listed in Table 2, and the distribution of r2 within a kilobase for all loci pooled is displayed in Figure 2. Figure 3 illustrates the significance of these pairwise comparisons at both distances. Generally, we observed significant levels of LD within 1 kb distances, with an overall average of 0.38, ranging from 0.19 to 0.72. LD decayed slowly within a kilobase (lines in Figure 2), with the fastest rate of decay observed in the New York population. At the larger scale, average r2 values were low relative to the shorter spatial scale with an average of 0.20, although some of these comparisons were still significant (Figure 3).
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.
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.
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.
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.
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.
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.
Liu A, Burke JM: Patterns of nucleotide diversity in wild and cultivated sunflower. Genetics. 2006, 173 (1): 321-330. 10.1534/genetics.105.051110.
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.
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.
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.
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.
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.
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.
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.
Andolfatto P, Przeworski M: A genome-wide departure from the standard neutral model in natural populations of Drosophila. Genetics. 2000, 156: 257-268.
Moran NA: The evolution of aphid life cycles. Ann Rev Ent. 1992, 37: 321-348. 10.1146/annurev.en.37.010192.001541.
Blackman RL: Reproduction, cytogenetics and development. Aphids: Their biology, natural enemies & control. Edited by: Minks AK, Harrewijn P. 1987, Amsterdam: Elsevier, 2A: 163-195.
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.
Via S: Inducing the sexual forms and hatching the eggs of pea aphids. Ent Exp Appl. 1992, 65: 119-127. 10.1007/BF00221362.
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.
Buchner P: Endosymbiosis of animals with plant microorganisms. 1965, New York: Interscience Publishers, Inc
Baumann P, Moran NA, Baumann L: The evolution and genetics of aphid endosymbionts. Bio Science. 1997, 47: 12-20.
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.
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.
Nault LR: Arthropod transmission of plant viruses: a new synthesis. Ann Ent Soc Am. 1997, 90: 521-541.
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.
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.
Nei M, Maruyama T, Chakraborty R: The bottleneck effect and genetic variability in populations. Evolution. 1975, 29: 1-10. 10.2307/2407137.
Crow JF, Kimura M: An Introduction to Population Genetics Theory. 1970, New York: Harper and Row
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.
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.
Hawthorne DJ, Via S: Genetic linkage of ecological specialization and reproductive isolation in pea aphids. Nature. 2001, 412: 904-907. 10.1038/35091062.
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.
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.
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.
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.
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.
Warnes GR: The genetics package. R News. 2003, 3: 9-13.
Team RDC: R: A language and environment for statistical computing. 2006, Vienna, Austria: R Foundation for Statistical Computing, 2.3.1
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.
Inc. SI: SAS/STAT user's guide. 1989, Cary, NC: SAS Institute, Inc, 2:
Nei M: Molecular Evolutionary Genetics. 1987, New York: Columbia University Press
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.
Tajima F: Statistical method for testing the neutral mutation hypothesis by DNA polymorphism. Genetics. 1989, 123 (3): 585-595.
Raymond M, Rousset F: GENEPOP (version 1.2): population genetics software for exact tests and ecumenicism. J Heredity. 1995, 86: 248-249.
Weir BC, Cockerham CC: Estimating F-statistics for the analysis of population structure. Evolution. 1984, 38: 1358-1370. 10.2307/2408641.
Guo SW, Thompson EA: Performing the exact test of Hardy-Weinberg proportion for multiple alleles. Biometrics. 1992, 48: 359-10.2307/2532296.
Tajima F: Evolutionary relationship of DNA sequences in finite populations. Genetics. 1983, 105: 437-460.
Barrette RJ, Crease TJ, Hebert PDN, Via S: Mitochondrial DNA diversity in the pea aphid, Acyrthosiphon pisum. Genome Biology. 1994, 37: 858-865.
Boulding EG: Molecular evidence against phylogenetically distinct host races of the pea aphid (Acyrthosiphon pisum). Genome Biology. 1998, 41: 769-775.
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.
Hartl DL, Clark AG: Principles of Population Genetics. 1997, Sinauer Associates
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.
Brisson JA, Nuzhdin SV: Rarity of males in pea aphids results in mutational decay. Science. 2008, 319: 58-10.1126/science.1147919.
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.
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.
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.
Foster WA: Aphid sex ratios. Sex Ratio Handbook: concepts and research methods. Edited by: Hardy I. 2002, Cambridge University Press, 254-265.
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-
Mayr E: Animal Species and Evolution. 1963, Cambridge: Harvard University Press
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.
JAB and DLS conceived of the study, JAB carried out the experiments, JAB and SVN performed the analyses, and all three authors participated in the writing of the manuscript. All authors read and approved the final manuscript.