Geographic variation in adaptation at the molecular level: a case study of plant immunity genes

Natural selection imposed by interacting species frequently varies among geographic locations and can lead to local adaptation, where alternative phenotypes are found in different populations. Little is known, however, about whether geographically variable selection acting on traits that mediate species interactions is consistent or strong enough to influence patterns of nucleotide variation at individual loci. To investigate this question, we examined patterns of nucleotide diversity and population structure at 16 plant innate immunity genes, with putative functions in defending plants against pathogens or herbivores, from six populations of teosinte (Zea mays ssp. parviglumis). Specifically, we tested whether patterns of population structure and within-population diversity at immunity genes differed from patterns found at nonimmunity (reference) loci and from neutral expectations derived from coalescent simulations of structured populations. For the majority of genes, we detected no strong evidence of geographically variable selection. However, in the wound-induced serine protease inhibitor (wip1), which inhibits the hydrolysis of dietary proteins in insect herbivores, one population showed unusually high levels of genetic differentiation, very low levels of nucleotide polymorphism, and was fixed for a novel replacement substitution in the active site of the protein. Taken together, these data suggest that wip1 experienced a recent selective sweep in one geographic region; this pattern may reflect local adaptation or an ongoing species-wide sweep. Overall, our results indicate that a signature of local adaptation at the molecular level may be uncommon-particularly for traits that are under complex genetic control.     

Using multilocus sequence data to assess population structure, natural selection, and linkage disequilibrium in wild tomatoes

We employed a multilocus approach to examine the effects of population subdivision and natural selection on DNA polymorphism in 2 closely related wild tomato species (Solanum peruvianum and Solanum chilense), using sequence data for 8 nuclear loci from populations across much of the species' range. Both species exhibit substantial levels of nucleotide variation. The species-wide level of silent nucleotide diversity is 18% higher in S. peruvianum (pi(sil) approximately 2.50%) than in S. chilense (pi(sil) approximately 2.12%). One of the loci deviates from neutral expectations, showing a clinal pattern of nucleotide diversity and haplotype structure in S. chilense. This geographic pattern of variation is suggestive of an incomplete (ongoing) selective sweep, but neutral explanations cannot be entirely dismissed. Both wild tomato species exhibit moderate levels of population differentiation (average F(ST) approximately 0.20). Interestingly, the pooled samples (across different demes) exhibit more negative Tajima's D and Fu and Li's D values; this marked excess of low-frequency polymorphism can only be explained by population (or range) expansion and is unlikely to be due to population structure per se. We thus propose that population structure and population/range expansion are among the most important evolutionary forces shaping patterns of nucleotide diversity within and among demes in these wild tomatoes. Patterns of population differentiation may also be impacted by soil seed banks and historical associations mediated by climatic cycles. Intragenic linkage disequilibrium (LD) decays very rapidly with physical distance, suggesting high recombination rates and effective population sizes in both species. The rapid decline of LD seems very promising for future association studies with the purpose of mapping functional variation in wild tomatoes.     

Subdivision and haplotype structure in natural populations of Arabidopsis lyrata

High levels of inbreeding are expected to cause a strong reduction in levels of genetic variability, effective recombination rates and in adaptation compared with related outcrossing populations. We examined patterns of DNA polymorphism at five nuclear loci and one chloroplast locus within and between four populations of the outcrossing plant Arabidopsis lyrata, a close relative of the highly self-fertilizing model species A. thaliana. The observed patterns are compared with species-wide polymorphism at orthologous loci, as well as within- and between-population patterns at other studied loci in A. thaliana. In addition to evidence for much higher average within-population diversity, species-wide levels of silent polymorphism are generally higher in A. lyrata than in A. thaliana, unlike the results from a previous study of the ADH locus. However, polymorphism is also low in the North American A. lyrata subspecies lyrata compared with the European subspecies petraea, suggesting either a population bottleneck in North American populations or recent admixture involving diverged European populations. Differentiation between the two subspecies is strong, although there are few fixed differences, suggesting that their isolation is recent. Estimates of intralocus recombination rates and analysis of haplotype structure in European A. lyrata populations indicate lower recombination than predicted based on the variability together with physical recombination rates estimated from A. thaliana. This may be due to strong population subdivision, or to recent departures from demographic equilibrium such as a bottleneck or population admixture. Alternatively, there may be consistently lower recombination rates in the outcrossing species. In contrast, estimates of recombination rates from species-wide samples of A. thaliana are close to the values expected assuming a high rate of self-fertilization. Complex population histories in both A. thaliana and A. lyrata complicate theoretical predictions and empirical tests of the effects of inbreeding on polymorphism and molecular evolution.     

Evidence for a large-scale population structure of Arabidopsis thaliana from genome-wide single nucleotide polymorphism markers

Population-based methods for the genetic mapping of adaptive traits and the analysis of natural selection require that the population structure and demographic history of a species are taken into account. We characterized geographic patterns of genetic variation in the model plant Arabidopsis thaliana by genotyping 115 genome-wide single nucleotide polymorphism (SNP) markers in 351 accessions from the whole species range using a matrix-assisted laser desorption/ionization time-of-flight assay, and by sequencing of nine unlinked short genomic regions in a subset of 64 accessions. The observed frequency distribution of SNPs is not consistent with a constant-size neutral model of sequence polymorphism due to an excess of rare polymorphisms. There is evidence for a significant population structure as indicated by differences in genetic diversity between geographic regions. Accessions from Central Asia have a low level of polymorphism and an increased level of genome-wide linkage disequilibrium (LD) relative to accessions from the Iberian Peninsula and Central Europe. Cluster analysis with the structure program grouped Eurasian accessions into K=6 clusters. Accessions from the Iberian Peninsula and from Central Asia constitute distinct populations, whereas Central and Eastern European accessions represent admixed populations in which genomes were reshuffled by historical recombination events. These patterns likely result from a rapid postglacial recolonization of Eurasia from glacial refugial populations. Our analyses suggest that mapping populations for association or LD mapping should be chosen from regional rather than a species-wide sample or identified genetically as sets of individuals with similar average genetic distances. Electronic Supplementary Material Supplementary material is available for this article at and is accessible for authorized users.     

GENETIC DIVERSITY AND POPULATION STRUCTURE IN CAMELLIA JAPONICA L. (THEACEAE)

Camellia japonica is a widespread and morphologically diverse tree native to parts of Japan and adjacent islands. Starch gel electrophoresis was used to score allelic variation at 20 loci in seeds collected from 60 populations distributed throughout the species range. In comparison with other plant species, the level of genetic diversity within C. japonica populations is very high: 66.2% of loci were polymorphic on average per population, with a mean number of 2.16 alleles per locus; the mean observed and panmictic heterozygosities were 0.230 and 0.265, respectively. Genotypic proportions at most loci in most populations fit Hardy-Weinberg expectations. However, small heterozygote deficiencies were commonly observed (mean population fixation index = 0.129). It is suggested that the most likely cause of the observed deficiencies is population subdivision into genetically divergent subpopulations. The overall level of population differentiation is greater than is typically observed in out-breeders: The mean genetic distance and identity (Nei's D and I) between pairs of populations were 0.073 and 0.930, respectively, and Wright's Fst was 0.144. Differences among populations appeared to be manifested as variation in gene frequencies at many loci rather than variation in allelic composition per se. However, the patterns of variation were not random. Reciprocal clinal variation of gene frequencies was observed for allele pairs at six loci. In addition, principal components analysis revealed that populations tended to genetically cluster into four regions representing the geographic areas Kyushu, Shikoku, western Honshu, and eastern Honshu. There was a significant relationship between genetic and geographic distance (r = 0.61; P < 0.01). Analysis of variance on allozyme frequencies showed that there was approximately four times as much differentiation among populations within regions, as among regions. It is likely that the observed patterns of population relationships result from the balance between genetic drift in small subpopulations and gene flow between them.     

Nucleotide variation in Quercus crispula Blume

Quercus is attractive for evolutionary studies, primarily for developing the concepts of the species, speciation and adaptation; however, remarkably little is known about levels of nucleotide polymorphism in the nuclear functional genes of this genus. This article provides the first characterization of levels of nucleotide polymorphism in 11 gene fragments in natural populations of a Quercus species, Quercus crispula Blume. Results show that the level of nucleotide variation in this oak is generally higher than that in conifers, as high as that in a European oak, but lower than that in an aspen. The level of population recombination is relatively high. Within-population inbreeding is negligible and between-population differentiation is modest. The decay of linkage disequilibrium is significantly faster in the species-wide samples and the three northernmost populations than in the other populations. Statistical tests support the hypothesis of a recent bottleneck for several populations in the southern part of Japan. The amounts and patterns of nucleotide variation, recombination and linkage disequilibrium, and genetic differentiation observed among populations of this species are contradictory to our expectations, given the recent colonization history of the northern Japan populations.     

Tracing the geographic origins of major avocado cultivars

It has been difficult to infer the genetic history of avocado breeding, owing to the role of hybridization in the origin of contemporary avocado cultivars. To address this difficulty, we used the model-based clustering program, STRUCTURE, and nucleotide polymorphism in 5960 bp of sequence from 4 nuclear loci to examine population structure in 21 wild avocado accessions. The origins of 33 cultivars were inferred relative to the wild sample. Nucleotide sequence diversity in domesticated avocados ranged between 80% and 90% of that observed for the same loci in wild avocado, depending on the diversity statistic used for comparison. Substantial genetic differentiation among 3 geographic groups of wild germplasm corresponded to the classically defined horticultural races of avocado. Previously undetected genetic differentiation was revealed in wild populations from Central Mexico, where 2 subpopulations were distinguished based on elevation and latitude.     

Continent‐wide population genomic structure and phylogeography of North America’s most destructive conifer defoliator,the spruce budworm (Choristoneura fumiferana)

The spruce budworm, Choristoneura fumiferana, is presumed to be panmictic across vast regions of North America. We examined the extent of panmixia by genotyping 3,650 single nucleotide polymorphism (SNP) loci in 1975 individuals from 128 collections across the continent. We found three spatially structured subpopulations: Western (Alaska, Yukon), Central (southeastern Yukon to the Manitoba–Ontario border), and Eastern (Manitoba–Ontario border to the Atlantic). Additionally, the most diagnostic genetic differentiation between the Central and Eastern subpopulations was chromosomally restricted to a single block of SNPs that may constitute an island of differentiation within the species. Geographic differentiation in the spruce budworm parallels that of its principal larval host, white spruce (Picea glauca), providing evidence that spruce budworm and spruce trees survived in the Beringian refugium through the Last Glacial Maximum and that at least two isolated spruce budworm populations diverged with spruce/fir south of the ice sheets. Gene flow in the spruce budworm may also be affected by mountains in western North America, habitat isolation in West Virginia, regional adaptations, factors related to dispersal, and proximity of other species in the spruce budworm species complex. The central and eastern geographic regions contain individuals that assign to Eastern and Central subpopulations, respectively, indicating that these barriers are not complete. Our discovery of previously undetected geographic and genomic structure in the spruce budworm suggests that further population modelling of this ecologically important insect should consider regional differentiation, potentially co‐adapted blocks of genes, and gene flow between subpopulations.     

Sonneratia ovata Backer–A genetically depauperate mangrove species

It is expected that geographically widespread and outcrossing tree species would have high level of genetic variation. Three chloroplast regions and six nuclear genes were sequenced to assess the genetic variation of a species of mangroves, Sonneratia ovata, from China and Thailand. No nucleotide polymorphism was found in the three chloroplast regions and in the six nuclear genes from all of the four populations examined. The depauperation of polymorphism in S. ovata in comparison to moderate polymorphisms of other congeneric species is surprising, particularly considering high level of polymorphism in the past and relatively wide geographic distribution of the species. Since multiple independent loci were surveyed in this study, the most plausible explanation for our observation is that S. ovata has experienced severe demographic bottlenecks in the Pleistocene glaciation, followed by subsequent recolonization by sea currents in China and Southeast Asia regions. The lack of polymorphism may also be attributable to its small population size, despite its wide geographic distribution.     

High nucleotide polymorphism and rapid decay of linkage disequilibrium in wild populations of Caenorhabditis remanei

The common ancestor of the self-fertilizing nematodes Caenorhabditis elegans and C. briggsae must have reproduced by obligate outcrossing, like most species in this genus. However, we have only a limited understanding about how genetic variation is patterned in such male-female (gonochoristic) Caenorhabditis species. Here, we report results from surveying nucleotide variation of six nuclear loci in a broad geographic sample of wild isolates of the gonochoristic C. remanei. We find high levels of diversity in this species, with silent-site diversity averaging 4.7%, implying an effective population size close to 1 million. Additionally, the pattern of polymorphisms reveals little evidence for population structure or deviation from neutral expectations, suggesting that the sampled C. remanei populations approximate panmixis and demographic equilibrium. Combined with the observation that linkage disequilibrium between pairs of polymorphic sites decays rapidly with distance, this suggests that C. remanei will provide an excellent system for identifying the genetic targets of natural selection from deviant patterns of polymorphism and linkage disequilibrium. The patterns revealed in this obligately outcrossing species may provide a useful model of the evolutionary circumstances in C. elegans' gonochoristic progenitor. This will be especially important if self-fertilization evolved recently in C. elegans history, because most of the evolutionary time separating C. elegans from its known relatives would have occurred in a state of obligate outcrossing.     

Spatial heterogeneity of mitochondrial DNA and allozymes among populations of white-tailed deer and mule deer.

A white-tailed deer (Odocoileus virginianus) population in northeastern Minnesota and a mule deer (O. hemionus) population in the Bridger Mountains Montana, have previously been shown to be spatially subdivided into contiguous subpopulations. We assessed the degree of genetic differentiation among subpopulations and tested the hypothesis that differentiation will be greater for mitochondrial DNA (mtDNA) than for nuclear-encoded allozymes. Differentiation of the white-tailed deer subpopulations was significant for two allozyme loci but not for mtDNA, and the overall degree of differentiation was low. Gene flow, recent founding of the subpopulations, and polygamous breeding structure may all have contributed to this pattern. Greater differentiation was evident among disjunct populations than between the contiguous subpopulations of white-tailed deer. The contiguous mule deer subpopulations were significantly differentiated for mtDNA and one allozyme locus. Differentiation was greater for mtDNA than for allozymes. These results are consistent with demographic data that indicate mule deer males disperse more than do females. Disjunct mule deer populations may be similar or dramatically different in mtDNA haplotype frequencies that do not necessarily vary with geographic distance. Current and historical gene flow and breeding structure will influence population genetic patterns.     

Genetic diversity and the evolutionary history of plant immunity genes in two species of Zea

Plant pathogenesis-related genes (PR genes) code for enzymes, enzyme inhibitors, and other peptides that confer resistance to pathogens and herbivores. Although several PR genes have been the subject of molecular population genetic analyses, a general understanding of their long-term evolutionary dynamics remains incomplete. Here we analyze sequence data from 17 PR genes from two closely related teosinte species of central Mexico. In addition to testing whether patterns of diversity at individual loci depart from expectations under a neutral model, we compared patterns of diversity at defense genes, as a class, to nondefense genes. In Zea diploperennis, the majority of defense genes have patterns of diversity consistent with neutral expectations while at least two genes showed evidence of recent positive selection consistent with arms-race models of antagonistic coevolution. In Zea mays ssp. parviglumis, by contrast, analyses of both defense and nondefense genes revealed strong and consistent departures from the neutral model, suggestive of nonequilibrium population dynamics or population structure. Nevertheless, we found a significant excess of replacement polymorphism in defense genes compared to nondefense genes. Although we cannot exclude relaxed selective constraint as an explanation, our results are consistent with temporally variable (transient or episodic) selection or geographically variable selection acting on parviglumis defense genes. The different patterns of diversity found in the two Zea species may be explained by parviglumis' greater distribution and population structure together with geographic variation in selection.     

Genetic variation,population structure and identification of yellow catfish, <Emphasis Type="Italic">Mystus nemurus</Emphasis> (C&;V) in Thailand using RAPD,ISSR and SCAR marker

Random amplified polymorphic DNA (RAPD) and inter-simple sequence repeat (ISSR) markers were used to investigate the genetic structure of four subpopulations of Mystus nemurus in Thailand. The 7 RAPD and 7 ISSR primers were selected. Of 83 total RAPD fragments, 80 (96.39%) were polymorphic loci, and of 81 total ISSR fragments, 75 (92.59%) were polymorphic loci. Genetic variation and genetic differentiation obtained from RAPD fragments or ISSR fragments showed similar results. Percentage of polymorphic loci (%P), observed number of alleles, effective number of alleles, Nei’s gene diversity (H) and Shannon’s information index revealed moderate to high level of genetic variations within each M. nemurus subpopulation and overall population. High levels of genetic differentiations were received from pairwise unbiased genetic distance (D) and coefficient of differentiation. Mantel test between D or gene flow and geographical distance showed a low to moderate correlation. Analysis of molecular variance indicated that variations among subpopulations were higher than those within subpopulations. The UPGMA dendrograms, based on RAPD and ISSR, showing the genetic relationship among subpopulations are grouped into three clusters; Songkhla (SK) subpopulation was separated from the other subpopulations. The candidate species-specific and subpopulation-specific RAPD fragments were sequenced and used to design sequence-characterized amplified region primers which distinguished M. nemurus from other species and divided SK subpopulation from the other subpopulations. The markers used in this study should be useful for breeding programs and future aquacultural development of this species in Thailand.     

THE VARIABLE GENOMIC ARCHITECTURE OF ISOLATION BETWEEN HYBRIDIZING SPECIES OF HOUSE MICE

Studies of the genetics of hybrid zones can provide insight into the genomic architecture of species boundaries. By examining patterns of introgression of multiple loci across a hybrid zone, it may be possible to identify regions of the genome that have experienced selection. Here, we present a comparison of introgression in two replicate transects through the house mouse hybrid zone through central Europe, using data from 41 single nucleotide markers. Using both genomic and geographic clines, we found many differences in patterns of introgression between the two transects, as well as some similarities. We found that many loci may have experienced the effects of selection at linked sites, including selection against hybrid genotypes, as well as positive selection in the form of genotypes introgressed into a foreign genetic background. We also found many positive associations of conspecific alleles among unlinked markers, which could be caused by epistatic interactions. Different patterns of introgression in the two transects highlight the challenge of using hybrid zones to identify genes underlying isolation and raise the possibility that the genetic basis of isolation between these species may be dependent on the local population genetic make-up or the local ecological setting.     

A multilocus sequence survey in Arabidopsis thaliana reveals a genome-wide departure from a neutral model of DNA sequence polymorphism

The simultaneous analysis of multiple genomic loci is a powerful approach to studying the effects of population history and natural selection on patterns of genetic variation of a species. By surveying nucleotide sequence polymorphism at 334 randomly distributed genomic regions in 12 accessions of Arabidopsis thaliana, we examined whether a standard neutral model of nucleotide sequence polymorphism is consistent with observed data. The average nucleotide diversity was 0.0071 for total sites and 0.0083 for silent sites. Although levels of diversity are variable among loci, no correlation with local recombination rate was observed, but polymorphism levels were correlated for physically linked loci (<250 kb). We found that observed distributions of Tajima's D- and D/D(min)- and of Fu and Li's D-, D*- and F-, F*-statistics differed significantly from the expected distributions under a standard neutral model due to an excess of rare polymorphisms and high variances. Observed and expected distributions of Fay and Wu's H were not different, suggesting that demographic processes and not selection at multiple loci are responsible for the deviation from a neutral model. Maximum-likelihood comparisons of alternative demographic models like logistic population growth, glacial refugia, or past bottlenecks did not produce parameter estimates that were more consistent with observed patterns. However, exclusion of highly polymorphic "outlier loci" resulted in a fit to the logistic growth model. Various tests of neutrality revealed a set of candidate loci that may evolve under selection.     

Among population differentiation at nuclear genes in native Scots pine (Pinus sylvestris L.) in Scotland

In the Scottish Highlands, Scots pine is at the north-western extreme of its wide natural distribution. Here, the remaining native populations are patchily distributed in highly variable environments, from the more continental, drier eastern Highlands to the milder, wetter Atlantic Ocean coast. As these pinewoods are the remnants of a naturally established forest, they form a valuable system for analysis of genetic and adaptive variation in heterogeneous environments. Using samples from across the Scottish population, we analysed data from nuclear and mitochondrial genes to assess patterns of within and between population genetic variation. Within population diversity levels were high, and significant genetic differentiation among pairs of Scottish populations at relatively small spatial scales was present at several nuclear loci. At these loci, no differentiation had been found among continental populations, even those separated by large geographic distances. Overall, no clear clustering of Scottish samples was found in population structure analysis suggesting that geographically distant populations with high intra-population nucleotide diversity are not strongly isolated or diverged from each other. Scottish populations lacked a mitotype that is widespread in eastern and north-eastern Europe, indicating that pines from that area may not have participated in the most recent colonisation of the British Isles.     

Range-wide genetic differentiation of Eugenia dysenterica (Myrtaceae) populations in Brazilian Cerrado

The analysis of geographic patterns in genetic variation has been one of the most important current tools to understand ecological and evolutionary processes underlying population structure. However, inferring such processes from population data may be misleading if biased geographic samples are analyzed. Here we expand previous analyses of Eugenia dysenterica population structure in Brazilian Cerrado, analyzing a larger number of populations distributed throughout a broader geographic region covering most of species' range. We provide new estimates of genetic diversity and population structure based on SSR markers from both neutral and genic regions, using several cluster and ordination techniques. These analyzes reveal a continuous northwestern-southeastern gradient in population differentiation, and not two distinct clusters of populations as suggested in some previous studies. This more comprehensive analysis also reinforces that a simple process of stochastic differentiation do not explain the observed patterns. Moreover, we conclude that explanations for population differentiation may focus on why genetic diversity decreases toward southeastern populations and not necessarily on barriers and interruption of gene flow creating regional patterns of population differentiation.     

Population structure and genetic basis of the agronomic traits of upland cotton in China revealed by a genome‐wide association study using high‐density SNPs

Gossypium hirsutum L. represents the largest source of textile fibre, and China is one of the largest cotton‐producing and cotton‐consuming countries in the world. To investigate the genetic architecture of the agronomic traits of upland cotton in China, a diverse and nationwide population containing 503 G. hirsutum accessions was collected for a genome‐wide association study (GWAS) on 16 agronomic traits. The accessions were planted in four places from 2012 to 2013 for phenotyping. The CottonSNP63K array and a published high‐density map based on this array were used for genotyping. The 503 G. hirsutum accessions were divided into three subpopulations based on 11 975 quantified polymorphic single‐nucleotide polymorphisms (SNPs). By comparing the genetic structure and phenotypic variation among three genetic subpopulations, seven geographic distributions and four breeding periods, we found that geographic distribution and breeding period were not the determinants of genetic structure. In addition, no obvious phenotypic differentiations were found among the three subpopulations, even though they had different genetic backgrounds. A total of 324 SNPs and 160 candidate quantitative trait loci (QTL) regions were identified as significantly associated with the 16 agronomic traits. A network was established for multieffects in QTLs and interassociations among traits. Thirty‐eight associated regions had pleiotropic effects controlling more than one trait. One candidate gene, Gh_D08G2376, was speculated to control the lint percentage (LP). This GWAS is the first report using high‐resolution SNPs in upland cotton in China to comprehensively investigate agronomic traits, and it provides a fundamental resource for cotton genetic research and breeding.