Molecular population genetics and phenotypic diversification of two populations of the thermophilic cyanobacterium Mastigocladus laminosus

We investigated the distributions of genetic and phenotypic variation for two Yellowstone National Park populations of the heterocyst-forming cyanobacterium Mastigocladus (Fischerella) laminosus that exhibit dramatic phenotypic differences as a result of environmental differences in nitrogen availability. One population develops heterocysts and fixes nitrogen in situ in response to a deficiency of combined nitrogen in its environment, whereas the other population does neither due to the availability of a preferred nitrogen source. Slowly evolving molecular markers, including the 16S rRNA gene and the downstream internal transcribed spacer, are identical among all laboratory isolates from both populations but belie considerable genetic and phenotypic diversity. The total nucleotide diversity at six nitrogen metabolism loci was roughly three times greater than that observed for the human global population. The two populations are genetically differentiated, although variation in performance on different nitrogen sources among genotypes could not be explained by local adaptation to available nitrogen in the respective environments. Population genetic models suggest that local adaptation is mutation limited but also that the populations are expected to continue to diverge due to low migratory gene flow.     

Population differences in Chinook salmon (Oncorhynchus tshawytscha) DNA methylation: Genetic drift and environmental factors

Local adaptation and phenotypic differences among populations have been reported in many species, though most studies focus on either neutral or adaptive genetic differentiation. With the discovery of DNA methylation, questions have arisen about its contribution to individual variation in and among natural populations. Previous studies have identified differences in methylation among populations of organisms, although most to date have been in plants and model animal species. Here we obtained eyed eggs from eight populations of Chinook salmon (Oncorhynchus tshawytscha) and assayed DNA methylation at 23 genes involved in development, immune function, stress response, and metabolism using a gene‐targeted PCR‐based assay for next‐generation sequencing. Evidence for population differences in methylation was found at eight out of 23 gene loci after controlling for developmental timing in each individual. However, we found no correlation between freshwater environmental parameters and methylation variation among populations at those eight genes. A weak correlation was identified between pairwise DNA methylation dissimilarity among populations and pairwise F _ST based on 15 microsatellite loci, indicating weak effects of genetic drift or geographic distance on methylation. The weak correlation was primarily driven by two genes, GTIIBS and Nkef. However, single‐gene Mantel tests comparing methylation and pairwise F _ST were not significant after Bonferroni correction. Thus, population differences in DNA methylation are more likely related to unmeasured oceanic environmental conditions, local adaptation, and/or genetic drift. DNA methylation is an additional mechanism that contributes to among population variation, with potential influences on organism phenotype, adaptive potential, and population resilience.     

Phenotypic diversity,population structure and stress protein-based capacitoring in populations of Xeropicta derbentina,a heat-tolerant land snail species

The shell colour of many pulmonate land snail species is highly diverse. Besides a genetic basis, environmentally triggered epigenetic mechanisms including stress proteins as evolutionary capacitors are thought to influence such phenotypic diversity. In this study, we investigated the relationship of stress protein (Hsp70) levels with temperature stress tolerance, population structure and phenotypic diversity within and among different populations of a xerophilic Mediterranean snail species (Xeropicta derbentina). Hsp70 levels varied considerably among populations, and were significantly associated with shell colour diversity: individuals in populations exhibiting low diversity expressed higher Hsp70 levels both constitutively and under heat stress than those of phenotypically diverse populations. In contrast, population structure (cytochrome c oxidase subunit I gene) did not correlate with phenotypic diversity. However, genetic parameters (both within and among population differences) were able to explain variation in Hsp70 induction at elevated but non-pathologic temperatures. Our observation that (1) population structure had a high explanatory potential for Hsp70 induction and that (2) Hsp70 levels, in turn, correlated with phenotypic diversity while (3) population structure and phenotypic diversity failed to correlate provides empirical evidence for Hsp70 to act as a mediator between genotypic variation and phenotype and thus for chaperone-driven evolutionary capacitance in natural populations.     

Epigenetic variation associated with responses to different habitats in the context of genetic divergence in Phragmites australis

The mechanisms underlying heritable phenotypic divergence associated with adaptation in response to environmental stresses may involve both genetic and epigenetic variations. Several prior studies have revealed even higher levels of epigenetic variation than genetic variation. However, few population‐level studies have explored the effects of epigenetic variation on species with high levels of genetic diversity distributed across different habitats. Using AFLP and methylation‐sensitive AFLP markers, we tested the hypothesis that epigenetic variation may contribute to differences in plants occupying different habitats when genetic variation alone cannot fully explain adaptation. As a cosmopolitan invasive species, Phragmites australis (common reed) together with high genetic diversity and remarkable adaptability has been suggested as a model for responses to global change and indicators of environmental fluctuations. We found high levels of genetic and epigenetic diversity and significant genetic/epigenetic structure within each of 12 studied populations sampled from four natural habitats of P. australis. Possible adaptive epigenetic variation was suggested by significant correlations between DNA methylation‐based epigenetic differentiation and adaptive genetic divergence in populations across the habitats. Meanwhile, various AMOVAs indicated that some epigenetic differences may respond to various local habitats. A partial Mantel test was used to tease out the correlations between genetic/epigenetic variation and habitat after controlling for the correlation between genetic and epigenetic variations. We found that epigenetic diversity was affected mostly by soil nutrient availability, suggesting that at least some epigenetic differentiation occurred independently of genetic variation. We also found stronger correlations between epigenetic variation and phenotypic traits than between genetic variation and such traits. Overall, our findings indicate that genetically based differentiation correlates with heterogeneous habitats, while epigenetic variation plays an important role in ecological differentiation in natural populations of P. australis. In addition, our results suggest that when assessing global change responses of plant species, intraspecific variation needs to be considered.     

Genetic differentiation in natural populations of a Keystone Bunchgrass (<Emphasis Type="Italic">Aristida stricta</Emphasis>) across its native range

Aristida stricta Michx. (Poaceae) is a perennial bunchgrass native to the Southeastern Coastal Plain of North America where it is a keystone species in the longleaf pine savannas and slash pine flatwoods from southeastern North Carolina to Florida, and westward to the coast of Mississippi. We examined genetic relationships within and among ten populations of A. stricta by using eight inter-simple sequence repeat (ISSR) markers to generate band frequency data for 32 individuals from each sampled population. An analysis of molecular variance showed that 38% of the variation resided among populations while 62% was attributable to variation within populations. Grouping the populations by habitat or by geographic location did not show significant differentiation between the groups. Overall, pair-wise geographic and genetic distances were not correlated. Data indicate that while individuals within each population are genetically diverse, there seemingly are barriers to gene flow across populations leading to their divergence. Each population contains several exclusive loci suggesting that limited gene flow and/or genetic drift are likely leading to this pattern of localization. Our results, coupled with those of the previous studies that presented evidence for local adaptation and phenotypic differences among populations, suggest that there is sufficient differentiation among populations of this species to warrant: (1) maintenance of the existing genetic diversity at individual sites, and (2) use of local seed and plant sources for conservation projects.     

Do different rates of gene flow underlie variation in phenotypic and phenological clines in a montane grasshopper community?

Species responses to environmental change are likely to depend on existing genetic and phenotypic variation, as well as evolutionary potential. A key challenge is to determine whether gene flow might facilitate or impede genomic divergence among populations responding to environmental change, and if emergent phenotypic variation is dependent on gene flow rates. A general expectation is that patterns of genetic differentiation in a set of codistributed species reflect differences in dispersal ability. In less dispersive species, we predict greater genetic divergence and reduced gene flow. This could lead to covariation in life‐history traits due to local adaptation, although plasticity or drift could mirror these patterns. We compare genome‐wide patterns of genetic structure in four phenotypically variable grasshopper species along a steep elevation gradient near Boulder, Colorado, and test the hypothesis that genomic differentiation is greater in short‐winged grasshopper species, and statistically associated with variation in growth, reproductive, and physiological traits along this gradient. In addition, we estimate rates of gene flow under competing demographic models, as well as potential gene flow through surveys of phenological overlap among populations within a species. All species exhibit genetic structure along the elevation gradient and limited gene flow. The most pronounced genetic divergence appears in short‐winged (less dispersive) species, which also exhibit less phenological overlap among populations. A high‐elevation population of the most widespread species, Melanoplus sanguinipes, appears to be a sink population derived from low elevation populations. While dispersal ability has a clear connection to the genetic structure in different species, genetic distance does not predict growth, reproductive, or physiological trait variation in any species, requiring further investigation to clearly link phenotypic divergence to local adaptation.     

Invasion strategies in clonal aquatic plants: are phenotypic differences caused by phenotypic plasticity or local adaptation?

Background and Aims

The successful spread of invasive plants in new environments is often linked to multiple introductions and a diverse gene pool that facilitates local adaptation to variable environmental conditions. For clonal plants, however, phenotypic plasticity may be equally important. Here the primary adaptive strategy in three non-native, clonally reproducing macrophytes (Egeria densa, Elodea canadensis and Lagarosiphon major) in New Zealand freshwaters were examined and an attempt was made to link observed differences in plant morphology to local variation in habitat conditions.

Methods

Field populations with a large phenotypic variety were sampled in a range of lakes and streams with different chemical and physical properties. The phenotypic plasticity of the species before and after cultivation was studied in a common garden growth experiment, and the genetic diversity of these same populations was also quantified.

Key Results

For all three species, greater variation in plant characteristics was found before they were grown in standardized conditions. Moreover, field populations displayed remarkably little genetic variation and there was little interaction between habitat conditions and plant morphological characteristics.

Conclusions

The results indicate that at the current stage of spread into New Zealand, the primary adaptive strategy of these three invasive macrophytes is phenotypic plasticity. However, while limited, the possibility that genetic diversity between populations may facilitate ecotypic differentiation in the future cannot be excluded. These results thus indicate that invasive clonal aquatic plants adapt to new introduced areas by phenotypic plasticity. Inorganic carbon, nitrogen and phosphorous were important in controlling plant size of E. canadensis and L. major, but no other relationships between plant characteristics and habitat conditions were apparent. This implies that within-species differences in plant size can be explained by local nutrient conditions. All together this strongly suggests that invasive clonal aquatic plants adapt to a wide range of habitats in introduced areas by phenotypic plasticity rather than local adaptation.     

The influence of population structure on gene expression and flowering time variation in the ubiquitous weed Capsella bursa‐pastoris (Brassicaceae)

Population structure is a potential problem when testing for adaptive phenotypic differences among populations. The observed phenotypic differences among populations can simply be due to genetic drift, and if the genetic distance between them is not considered, the differentiation may be falsely interpreted as adaptive. Conversely, adaptive and demographic processes might have been tightly associated and correcting for the population structure may lead to false negatives. Here, we evaluated this problem in the cosmopolitan weed Capsella bursa‐pastoris. We used RNA‐Seq to analyse gene expression differences among 24 accessions, which belonged to a much larger group that had been previously characterized for flowering time and circadian rhythm and were genotyped using genotyping‐by‐sequencing (GBS) technique. We found that clustering of accessions for gene expression retrieved the same three clusters that were obtained with GBS data previously, namely Europe, the Middle East and Asia. Moreover, the three groups were also differentiated for both flowering time and circadian rhythm variation. Correction for population genetic structure when analysing differential gene expression analysis removed all differences among the three groups. This may suggest that most differences are neutral and simply reflect population history. However, geographical variation in flowering time and circadian rhythm indicated that the distribution of adaptive traits might be confounded by population structure. To bypass this confounding effect, we compared gene expression differentiation between flowering ecotypes within the genetic groups. Among the differentially expressed genes, FLOWERING LOCUS C was the strongest candidate for local adaptation in regulation of flowering time.     

Toxic hydrogen sulfide and dark caves: phenotypic and genetic divergence across two abiotic environmental gradients in Poecilia mexicana

Divergent natural selection drives evolutionary diversification. It creates phenotypic diversity by favoring developmental plasticity within populations or genetic differentiation and local adaptation among populations. We investigated phenotypic and genetic divergence in the livebearing fish Poecilia mexicana along two abiotic environmental gradients. These fish typically inhabit nonsulfidic surface rivers, but also colonized sulfidic and cave habitats. We assessed phenotypic variation among a factorial combination of habitat types using geometric and traditional morphometrics, and genetic divergence using quantitative and molecular genetic analyses. Fish in caves (sulfidic or not) exhibited reduced eyes and slender bodies. Fish from sulfidic habitats (surface or cave) exhibited larger heads and longer gill filaments. Common-garden rearing suggested that these morphological differences are partly heritable. Population genetic analyses using microsatellites as well as cytochrome b gene sequences indicate high population differentiation over small spatial scale and very low rates of gene flow, especially among different habitat types. This suggests that divergent environmental conditions constitute barriers to gene flow. Strong molecular divergence over short distances as well as phenotypic and quantitative genetic divergence across habitats in directions classic to fish ecomorphology suggest that divergent selection is structuring phenotypic variation in this system.     

Divergent Selection and Local Adaptation in Disjunct Populations of an Endangered Conifer,Keteleeria davidiana var. formosana (Pinaceae)

The present study investigated the genetic diversity, population structure, F _ST outliers, and extent and pattern of linkage disequilibrium in five populations of Keteleeria davidiana var. formosana, which is listed as a critically endangered species by the Council of Agriculture, Taiwan. Twelve amplified fragment length polymorphism primer pairs generated a total of 465 markers, of which 83.74% on average were polymorphic across populations, with a mean Nei’s genetic diversity of 0.233 and a low level of genetic differentiation (approximately 6%) based on the total dataset. Linkage disequilibrium and HICKORY analyses suggested recent population bottlenecks and inbreeding in K. davidiana var. formosana. Both STRUCTURE and BAPS observed extensive admixture of individual genotypes among populations based on the total dataset in various clustering scenarios, which probably resulted from incomplete lineage sorting of ancestral variation rather than a high rate of recent gene flow. Our results based on outlier analysis revealed generally high levels of genetic differentiation and suggest that divergent selection arising from environmental variation has been driven by differences in temperature, precipitation, and humidity. Identification of ecologically associated outliers among environmentally disparate populations further support divergent selection and potential local adaptation.     

Genetic diversity analysis of Hepatacodium miconioides at different altitude in Tiantai Mountain, Zhejiang Province, China and its relationship with environmental factors

The genetic diversity within and among populations of Hepatacodium miconioides collected at three different altitudes in Tiantai Mountain, Zhejiang Province and its relationships to environmental factors were analyzed by random amplified polymorphic DNA (RAPD) technique. Amplification using 12 random primers of 60 plants and 122 repetitive loci were produced. The percentage of polymorphic loci of three populations ranged from 18.85% to 23.77% with an average of 21.86%, indicating the relatively low genetic diversity of H. miconioides. The average Shannon index of phenotypic diversity (0.1329) and Nei index (0.0925) within populations were relatively low. A distinct genetic differentiation existed among populations of H. miconioides in spite of the relatively small geographical distribution. The average genetic diversity within populations of H. miconioides accounted for 33.58% of the total genetic diversity while the genetic diversity among populations accounted for 66.42% as estimated by the Shannon index of phenotypic diversity, The genetic differentiation among populations of H. miconioides was 0.6546, as estimated by Nei index. The gene flow estimated from G _ST was only 0.2656 and it indicated that gene flow among populations of H. miconioides was relatively low. The mean value of the genetic identity among populations of H. miconioides was 0.7126 and the average of genetic distance of H. miconioides was 0.3412. The genetic identity between populations at the elevation of 990 m and at the elevation of 780 m was the highest. The genetic identity between population at the elevation 500 m and other two populations was relatively low. The correlation analysis showed that the genetic diversity within populations was significantly related with the soil total nitrogen.     

Fine‐grained adaptive divergence in an amphibian: genetic basis of phenotypic divergence and the role of nonrandom gene flow in restricting effective migration among wetlands

Adaptive ecological differentiation among sympatric populations is promoted by environmental heterogeneity, strong local selection and restricted gene flow. High gene flow, on the other hand, is expected to homogenize genetic variation among populations and therefore prevent local adaptation. Understanding how local adaptation can persist at the spatial scale at which gene flow occurs has remained an elusive goal, especially for wild vertebrate populations. Here, we explore the roles of natural selection and nonrandom gene flow (isolation by breeding time and habitat choice) in restricting effective migration among local populations and promoting generalized genetic barriers to neutral gene flow. We examined these processes in a network of 17 breeding ponds of the moor frog Rana arvalis, by combining environmental field data, a common garden experiment and data on variation in neutral microsatellite loci and in a thyroid hormone receptor (TRβ) gene putatively under selection. We illustrate the connection between genotype, phenotype and habitat variation and demonstrate that the strong differences in larval life history traits observed in the common garden experiment can result from adaptation to local pond characteristics. Remarkably, we found that haplotype variation in the TRβ gene contributes to variation in larval development time and growth rate, indicating that polymorphism in the TRβ gene is linked with the phenotypic variation among the environments. Genetic distance in neutral markers was correlated with differences in breeding time and environmental differences among the ponds, but not with geographical distance. These results demonstrate that while our study area did not exceed the scale of gene flow, ecological barriers constrained gene flow among contrasting habitats. Our results highlight the roles of strong selection and nonrandom gene flow created by phenological variation and, possibly, habitat preferences, which together maintain genetic and phenotypic divergence at a fine‐grained spatial scale.     

Standing variation boosted by multiple sources of introduction contributes to the success of the introduced species, Lotus corniculatus

Although ecological differences between native and introduced ranges have been considered to drive rapid expansion of invasive species, recent studies suggest that rapid evolutionary responses of invasive species to local environments may also be common. Such expansion across heterogeneous environments by adaptation to local habitats requires genetic variation. In this study, we investigated the source and role of standing variation in successful invasion of heterogeneous abiotic environments in a self-incompatible species, Lotus corniculatus. We compared phenotypic and genetic variation among cultivars, natives, and introduced genotypes, and found substantial genetic variation within both native and introduced populations. Introduced populations possessed genotypes derived from both cultivars and native populations, and had lower population differentiation, indicating multiple sources of introduction and population admixture among the sources in the introduced range. Both cultivars and introduced populations had similarly outperforming phenotypes on average, with increased biomass and earlier flowering compared with native populations, but those phenotypes were within the range of the variation in phenotypes of the native populations. In addition, clinal variation within introduced populations was detected along a climatic gradient. Multiple introductions from different sources, including cultivars, may have contributed to pre-adaptive standing variation in the current introduced populations. We conclude that both introduction of cultivar genotypes and natural selection in local environments contributed to current patterns of genetic and phenotypic variation observed in the introduced populations.     

Strong intraspecific variation in genetic diversity and genetic differentiation in Daphnia magna: the effects of population turnover and population size

Theory predicts that genetic diversity and genetic differentiation may strongly vary among populations of the same species depending on population turnover and local population sizes. Yet, despite the importance of these predictions for evolutionary and conservation issues, empirical studies comparing high‐turnover and low‐turnover populations of the same species are scarce. In this study, we used Daphnia magna, a freshwater crustacean, as a model organism for such a comparison. In the southern/central part of its range, D. magna inhabits medium‐sized, stable ponds, whereas in the north, it occurs in small rock pools with strong population turnover. We found that these northern populations have a significantly lower genetic diversity and higher genetic differentiation compared to the southern/central populations. Total genetic diversity across populations was only about half and average within‐population diversity only about a third of that in southern/central populations. Moreover, an average southern population contains more genetic diversity than the whole metapopulation system in the north. We based our analyses both on silent sites and microsatellites. The similarity of our results despite the contrasting mutation rates of these markers suggests that the differences are caused by contemporary rather than by historical processes. Our findings show that variation in population turnover and population size may have a major impact on the genetic diversity and differentiation of populations, and hence may lead to differences in evolutionary processes like local adaptation, hybrid vigour and breeding system evolution in different parts of a species range.     

Unraveling the Effects of Selection and Demography on Immune Gene Variation in Free-Ranging Plains Zebra (Equus quagga) Populations

Demography, migration and natural selection are predominant processes affecting the distribution of genetic variation among natural populations. Many studies use neutral genetic markers to make inferences about population history. However, the investigation of functional coding loci, which directly reflect fitness, is critical to our understanding of species'' ecology and evolution. Immune genes, such as those of the Major Histocompatibility Complex (MHC), play an important role in pathogen recognition and provide a potent model system for studying selection. We contrasted diversity patterns of neutral data with MHC loci, ELA-DRA and -DQA, in two southern African plains zebra (Equus quagga) populations: Etosha National Park, Namibia, and Kruger National Park, South Africa. Results from neutrality tests, along with observations of elevated diversity and low differentiation across populations, supported previous genus-level evidence for balancing selection at these loci. Despite being low, MHC divergence across populations was significant and may be attributed to drift effects typical of geographically separated populations experiencing little to no gene flow, or alternatively to shifting allele frequency distributions driven by spatially variable and fluctuating pathogen communities. At the DRA, zebra exhibited geographic differentiation concordant with microsatellites and reduced levels of diversity in Etosha due to highly skewed allele frequencies that could not be explained by demography, suggestive of spatially heterogeneous selection and local adaptation. This study highlights the complexity in which selection affects immune gene diversity and warrants the need for further research on the ecological mechanisms shaping patterns of adaptive variation among natural populations.     

Root architecture,plant size and soil nutrient variation in natural populations of <Emphasis Type="Italic">Arabidopsis thaliana</Emphasis>

Phenotypic variation in ecologically important traits may vary at large and small geographic scales, and may be shaped by natural selection. Here our explicit aim is to evaluate phenotypic differentiation among local populations and examine its relationship with ecological edaphic and climatic features that could lead to local adaptation. We characterized six populations of the model plant Arabidopsis thaliana over 3 years in the field in its native range. At each site, we measured edaphic conditions and aboveground and belowground phenotypes. In addition, we grew plants from the six characterized populations in a common greenhouse along with an additional fifteen populations from the Iberian Peninsula to examine evolutionary and ecological differentiation among populations, and relationships between geographic and ecological distance to phenotypic differences among populations. Significant differences in aboveground and belowground traits, population density, and micro- and macro-nutrient soil concentrations were found among the field populations. In particular, root architectural traits differed significantly among field populations. Complex patterns of ecological differences among population and plant phenotypes emerged when examining edaphic conditions in the Extremadura region, and geographic and climate variables at a broader scale of the Iberian Peninsula. We report levels of phenotypic variation at the local scale comparable to those found at broad geographic scales and report that local edaphic conditions contribute to population-level phenotypic variation in root and shoot traits. To our knowledge, these are the first reports of among population root architectural variation from natural field populations for this model organism. We demonstrate how ecological features, such as soil nutrients, can be associated with the phenotypic variation of A. thaliana measured in natural populations and may contribute to adaptive differentiation at a local scale.     

Spatial variation in climate mediates gene flow across an island archipelago

High levels of gene flow among partially isolated populations can overwhelm selection and limit local adaptation. This process, known as “gene swamping,” can homogenize genetic diversity among populations and reduce the capacity of a species to withstand rapid environmental change. We studied brown anole lizards (Anolis sagrei) distributed across seven islands in The Bahamas. We used microsatellite markers to estimate gene flow among islands and then examined the correlation between thermal performance and island temperature. The thermal optimum for sprint performance was correlated with both mean and maximum island temperature, whereas performance breadth was not correlated with any measure of temperature variation. Gene flow between islands decreased as the difference between mean island temperatures increased, even when those islands were adjacent to one another. These data suggest that phenotypic variation is the result of either (1) local genetic adaptation with selection against immigrants maintaining variation in the thermal optimum, (2) irreversible forms of adaptive plasticity such that immigrants have reduced fitness, or (3) an interaction between fixed genetic differences and plasticity. In general, the patterns of gene flow we observed suggest that local thermal environments represent important ecological filters that can mediate gene flow on relatively fine geographic scales.     

Assessment of variability among morphological and molecular characters in wild populations of mint [Mentha longifolia (L.) L.] germplasm

Mentha longifolia is an important medicinal and aromatic perennial herb that exhibits wide distribution range from sub-tropical to temperate regions. In the present study, agro-morphological traits and genetic differences in 19 different populations of M. longifolia were studied to evaluate the level and extent of its diversity. Analysis of variance (ANOVA) showed that the different phenotypic characters show considerable differences among various populations and was significant at p < 0.05. Molecular diversity analysis performed by using arbitrary amplified eleven ISSR primers generated a total of 121 amplicons that range within the size of 200–2500 base pairs (bp). Each primer on average generated 11 amplicons with percentage polymorphism being 100. The analysis of molecular variance (AMOVA) showed more (64%) among population genetic diversity and less (36%) within the populations. Greater genetic differentiation (Gst = 0.6852) among these populations occurs due to low gene flow (Nm = 0.2297) and greater habitat variability. Geographic and genetic distances were positively correlated according to Mantel’s test. In order to remove any kind of biases, we used R software to perform cluster and redundancy analysis to analyse the extent of relatedness among studied populations. In terms of morphological and molecular aspects, the populations were grouped into four and five clusters respectively based on hierarchical clustering method. The results demonstrated that M. longifolia displays a great degree of morphological and genetic variation and can be utilized in breeding, genetic improvement, and gene bank conservation programmes in future.     

A Genome Scan for Genes Underlying Microgeographic-Scale Local Adaptation in a Wild Arabidopsis Species

Adaptive divergence at the microgeographic scale has been generally disregarded because high gene flow is expected to disrupt local adaptation. Yet, growing number of studies reporting adaptive divergence at a small spatial scale highlight the importance of this process in evolutionary biology. To investigate the genetic basis of microgeographic local adaptation, we conducted a genome-wide scan among sets of continuously distributed populations of Arabidopsis halleri subsp. gemmifera that show altitudinal phenotypic divergence despite gene flow. Genomic comparisons were independently conducted in two distinct mountains where similar highland ecotypes are observed, presumably as a result of convergent evolution. Here, we established a de novo reference genome and employed an individual-based resequencing for a total of 56 individuals. Among 527,225 reliable SNP loci, we focused on those showing a unidirectional allele frequency shift across altitudes. Statistical tests on the screened genes showed that our microgeographic population genomic approach successfully retrieve genes with functional annotations that are in line with the known phenotypic and environmental differences between altitudes. Furthermore, comparison between the two distinct mountains enabled us to screen out those genes that are neutral or adaptive only in either mountain, and identify the genes involved in the convergent evolution. Our study demonstrates that the genomic comparison among a set of genetically connected populations, instead of the commonly-performed comparison between two isolated populations, can also offer an effective screening for the genetic basis of local adaptation.