Genome-wide association to fine-scale ecological heterogeneity within a continuous population of Biscutella laevigata (Brassicaceae)

Population Genetic Structure and Diversity of Cryptic Species of the Plant Genus Macrocarpaea (Gentianaceae) from the Tropical Andes

The Pleistocene climatic oscillations (PCO) that provoked several cycles of glacial-interglacial periods are thought to have profoundly affected species distribution, richness and diversity around the world. While the effect of the PCO on population dynamics at temperate latitudes is well known, considerable questions remain about its impact on the biodiversity of neotropical mountains. Here, we use amplified fragment length polymorphism molecular markers (AFLPs) to investigate the phylogeography and genetic structure of 13 plant species belonging to the gentian genus Macrocarpaea (Gentianaceae) in the tropical Andes. These woody herbs, shrubs or small trees show complex and potentially reticulated relationships, including cryptic species. We show that populations of M. xerantifulva in the dry system of the Rio Marañón in northern Peru have lower levels of genetic diversity compared to other sampled species. We suggest that this is due to a recent demographic bottleneck resulting from the contraction of the montane wet forests into refugia because of the expansion of the dry system into the valley during the glacial cycles of the PCO. This may imply that the ecosystems of different valleys of the Andes might have responded differently to the PCO.

The effects of migration load, selfing, inbreeding depression, and the genetics of adaptation on autotetraploid versus diploid establishment in peripheral habitats

The distribution and abundance of polyploids has intrigued biologists since their discovery in the early 20th century. A pattern in nature that may give insight to processes that shape the distribution and abundance of polyploids is that polyploid populations are sometimes associated with peripheral habitats within the range of a species of mixed ploidy. Here, adaptation and competition of a diploid versus an autotetraploid population in a peripheral habitat are examined theoretically. It is shown that a nascent autotetraploid population adapts to and outcompetes a diploid population in the periphery when the rate of gamete dispersal is high, and when the mode of gene action is recessive for moderate to high rates of selfing. With additive or dominant modes of gene action, the conditions for an autotetraploid to outcompete a diploid in the periphery appear determined more by the rate of selfing and less by gamete dispersal. All of these results are based on empirical work that suggests inbreeding depression is higher in diploids versus autotetraploids. Generally, the results indicate that, although autotetraploids incur minority cytotype exclusion, diploids face burdens themselves. In the case of adaptation to a peripheral habitat, this burden is migration load from gamete and propagule dispersal.

Heritable epigenetic diversity for conservation and utilization of epigenetic germplasm resources of clonal East African Highland banana (EAHB) accessions

KEY MESSAGE

Genetically identical East African Highland banana (EAHB) clones are epigenetically diverse with heritable epialleles that can contribute to morphological diversity. Heritable epigenetic variation can contribute to agronomic traits in crops and should be considered in germplasm conservation. Despite the genetic uniformity arising from a genetic bottleneck of one ancestral clone, followed by subsequent vegetative propagation, East African Highland bananas (EAHBs) display significant phenotypic diversity potentially arising from somatic mutations, heritable epialleles and/or genotype-by-environment interactions. Here, we use DNA methylation profiling across EAHB accessions representing most of the primary EAHB genepool to demonstrate that the genetically uniform EAHB genepool harbours significant epigenetic diversity. By analysing 724 polymorphic DNA methylation sites by methylation-sensitive AFLP across 90 EAHB cultivars, we could differentiate the EAHB varieties according to their regions (Kenya and Uganda). In contrast, there was minimal association of DNA methylation variation with the five morphological groups that are used to classify EAHBs. We further analysed DNA methylation patterns in parent-offspring cohort, which were maintained in offspring generated by sexual (seed) and asexual (vegetative) propagation, with higher levels of altered DNA methylation observed in vegetatively generated offspring. Our results indicate that the phenotypic diversity of near-isogenic EAHBs is mirrored by considerable DNA methylation variation, which is transmitted between generations by both vegetative reproduction and seed reproduction. Genetically uniform vegetatively propagated crops such as EAHBs harbour considerable heritable epigenetic variation, where heritable epialleles could arise in offspring and contribute to functional traits. This study provides a basis for developing strategies for conservation of epigenetic resources and for integration of epimarkers into crop breeding programmes.

Short-distance barriers affect genetic variability of Rhizophora mangle L. in the Yucatan Peninsula

The environmental variability at local scale results in different physiognomic types of mangrove forest. However, this variability has never been considered in studies of mangrove genetic variability. This study analyzed the genetic and morphological variability and structure of Rhizophora mangle at regional and local scales in the Yucatan Peninsula. Thirteen mangrove populations (eight scrub and five tall), located in seven sites, were sampled, and their morphological variability and relationship with the availability of phosphorus and salinity were analyzed. The diversity and genetic structure were estimated at different hierarchical levels with nine microsatellites, also Bayesian inference and Principal Coordinates Analysis were used. We found a great morphological variability of R. mangle that responded to local environmental variability and not to the precipitation gradient of the peninsula. The genetic diversity found in the peninsula was greater than that reported for other populations in Mexico and was grouped into two regions: the Gulf of Mexico and the Caribbean Sea. At a local scale, tall and scrub mangroves had significant genetic differentiation suggesting that ecological barriers promote genetic differentiation within sites. These results need to be considered in future population genetic studies and for mangrove management and conservation.

Multiscale landscape genomic models to detect signatures of selection in the alpine plant Biscutella laevigata

Plant species are known to adapt locally to their environment, particularly in mountainous areas where conditions can vary drastically over short distances. The climate of such landscapes being largely influenced by topography, using fine-scale models to evaluate environmental heterogeneity may help detecting adaptation to micro-habitats. Here, we applied a multiscale landscape genomic approach to detect evidence of local adaptation in the alpine plant Biscutella laevigata. The two gene pools identified, experiencing limited gene flow along a 1-km ridge, were different in regard to several habitat features derived from a very high resolution (VHR) digital elevation model (DEM). A correlative approach detected signatures of selection along environmental gradients such as altitude, wind exposure, and solar radiation, indicating adaptive pressures likely driven by fine-scale topography. Using a large panel of DEM-derived variables as ecologically relevant proxies, our results highlighted the critical role of spatial resolution. These high-resolution multiscale variables indeed indicate that the robustness of associations between genetic loci and environmental features depends on spatial parameters that are poorly documented. We argue that the scale issue is critical in landscape genomics and that multiscale ecological variables are key to improve our understanding of local adaptation in highly heterogeneous landscapes.

Epigenetic Variance, Performing Cooperative Structure with Genetics, Is Associated with Leaf Shape Traits in Widely Distributed Populations of Ornamental Tree Prunus mume

Increasing evidence shows that epigenetics plays an important role in phenotypic variance. However, little is known about epigenetic variation in the important ornamental tree Prunus mume. We used amplified fragment length polymorphism (AFLP) and methylation-sensitive amplified polymorphism (MSAP) techniques, and association analysis and sequencing to investigate epigenetic variation and its relationships with genetic variance, environment factors, and traits. By performing leaf sampling, the relative total methylation level (29.80%) was detected in 96 accessions of P. mume. And the relative hemi-methylation level (15.77%) was higher than the relative full methylation level (14.03%). The epigenetic diversity (I^∗ = 0.575, h^∗ = 0.393) was higher than the genetic diversity (I = 0.484, h = 0.319). The cultivated population displayed greater epigenetic diversity than the wild populations in both southwest and southeast China. We found that epigenetic variance and genetic variance, and environmental factors performed cooperative structures, respectively. In particular, leaf length, width and area were positively correlated with relative full methylation level and total methylation level, indicating that the DNA methylation level played a role in trait variation. In total, 203 AFLP and 423 MSAP associated markers were detected and 68 of them were sequenced. Homologous analysis and functional prediction suggested that the candidate marker-linked genes were essential for leaf morphology development and metabolism, implying that these markers play critical roles in the establishment of leaf length, width, area, and ratio of length to width.

Ten Years of Landscape Genomics: Challenges and Opportunities

Landscape genomics is a relatively new discipline that aims to reveal the relationship between adaptive genetic imprints in genomes and environmental heterogeneity among natural populations. Although the interest in landscape genomics has increased since this term was coined, studies on this topic remain scarce. Landscape genomics has become a powerful method to scan and determine the genes responsible for the complex adaptive evolution of species at population (mostly) and individual (more rarely) level. This review outlines the sampling strategies, molecular marker types and research categories in 37 articles published during the first 10 years of this field (i.e., 2007-2016). We also address major challenges and future directions for landscape genomics. This review aims to promote interest in conducting additional studies in landscape genomics.

Phylogeography of a tough rock survivor in European dry grasslands

Phylogeographic analyses of plants in Europe have revealed common glacial refugia and migration routes for several trees and herbs with arctic-alpine distributions. The postglacial histories of dry grassland species in central Europe have rarely been analyzed, even though the extremely species-rich habitat is threatened. Sedum album (Crassulaceae) is a common inhabitant of rocky sites in central European dry grasslands. We inferred the phylogeographic history of S. album over its distribution range in Europe. Genetic diversity within and differentiation between 34 S. album populations was examined using AFLP markers. Population isolation was indicated based on the rarity of the fragments and by isolation-by-distance effects. We sequenced the trnL-trnF region in 32 populations and used chloroplast microsatellites to analyze chloroplast haplotype distributions. Two distinct S. album lineages were detected. One lineage was comprised of populations from eastern and central parts of central Europe, and the Apennine Peninsula. A second lineage was comprised of populations from the Iberian Peninsula and western and northern parts of central Europe. Glacial refugia were identified based on the accumulation of ancient chloroplast haplotypes, high diversity of AFLP fragments within populations, and high levels of rare fragments in Liguria, Serbia, the Apennine and Iberian peninsulas. Cryptic refugia were detected in the Czech Republic and Slovakia. Isolation by distance was present all over the distribution range, and it was separately detected in southwestern and central Europe. In western Europe, where a contact zone between the two lineages can be expected, no isolation by distance was detected. Our results suggest migration routes of S. album northeastward from glacial refugia in southern Iberia, northward from the Apennine Peninsula, and northward and westward from the southeastern parts of central Europe. Therefore, central European grasslands were recently colonized by northern cryptic populations and source populations originating in the east and the Apennine Peninsula.

Geography Plays a More Important Role than Soil Composition on Structuring Genetic Variation of Pseudometallophyte Commelina communis

Pseudometallophytes are excellent models to study microevolution and local adaptation to soil pollution, as they can grow both on metalliferous and contrasting non-metalliferous soils. Although, there has been accumulating evidence for the effects of edaphic conditions and geographical isolation on the genetic structure of pesudometallophytes, it is still a difficult problem in evolutionary biology to assess their relative importance. In this study, we investigated the spatial patterns of genetic variability, population differentiation and genetic groups in pseudometallophyte Commelina communis with 12 microsatellite loci. Eight metallicolous and six non-metallicolous populations of C. communis were sampled from cupriferous sites and surrounding non-contaminated areas in China. Neither significant reduction in genetic diversity nor apparent founder and bottleneck effects were observed in metallicolous populations of C. communis. Based on Bayesian and Neighbor-Joining clustering analyses and a principal coordinates analysis, all sampled populations were found to be mainly separated into three genetic groups, corresponding well to their geographical locations rather than edaphic origins. Moreover, a significant and strong correlation between population genetic divergence and geographical distance were detected by Mantel test (r = 0.33; P < 0.05) and multiple matrix regression with randomization (MMRR; βD = 0.57, P < 0.01). However, the effect of copper concentration on genetic patterns of C. communis was not significant (MMRR; βE = -0.17, P = 0.12). Our study clearly demonstrated that the extreme edaphic conditions in metalliferous areas had limited effects on the genetic variability in C. communis. Geographic distance played a more important role in affecting the genetic structure of C. communis than soil composition did. In C. communis, the geographically disjunctive populations on metalliferous soils had multiple origins and evolved independently from nearby non-metallicolous populations.

AFLP Genome Scanning Reveals Divergent Selection in Natural Populations of Liriodendron chinense (Magnoliaceae) along a Latitudinal Transect

Understanding adaptive genetic variation and its relation to environmental factors are important for understanding how plants adapt to climate change and for managing genetic resources. Genome scans for the loci exhibiting either notably high or low levels of population differentiation (outlier loci) provide one means of identifying genomic regions possibly associated with convergent or divergent selection. In this study, we combined Amplified Fragment Length Polymorphism (AFLP) genome scan and environmental association analysis to test for signals of natural selection in natural populations of Liriodendron chinense (Chinese Tulip Tree; Magnoliaceae) along a latitudinal transect. We genotyped 276 individuals from 11 populations of L. chinense using 987 AFLP markers. Both frequency-based (Dfdist and BayeScan) and correlation-based (MLM) methods were applied to detect outlier loci. Our analyses recovered both neutral and potentially adaptive genetic differentiation among populations of L. chinense. We found moderate genetic diversity within populations and high genetic differentiation among populations with reduced genetic diversity toward the periphery of the species ranges. Nine AFLP marker loci showed evidence of being outliers for population differentiation for both detection methods. Of these, six were strongly associated with at least one climate factor. Temperature, precipitation, and radiation were found to be three important factors influencing local adaptation of L. chinense. The outlier AFLP loci are likely not the target of natural selection, but the neighboring genes of these loci might be involved in local adaptation. Hence, these candidates should be validated by further studies.

Inheritance and Variation of Genomic DNA Methylation in Diploid and Triploid Pacific Oyster (Crassostrea gigas)

DNA methylation is an important epigenetic mechanism that could be responsive to environmental changes indicating a potential role in natural selection and adaption. In order to evaluate an evolutionary role of DNA methylation, it is essential to first gain a better insight into inheritability. To address this question, this study investigated DNA methylation variation from parents to offspring in the Pacific oyster Crassostrea gigas using fluorescent-labeled methylation-sensitive amplified polymorphism (F-MSAP) analysis. Most of parental methylated loci were stably transmitted to offspring segregating following Medelian expectation. However, methylated loci deviated more often than non-methylated loci and offspring showed a few de novo methylated loci indicating DNA methylation changes from parents to offspring. Interestingly, some male-specific methylated loci were found in this study which might help to explore sex determination in oyster. Despite environmental stimuli, genomic stresses such as polyploidization also can induce methylation changes. This study also compared global DNA methylation level and individual methylated loci between diploid and triploid oysters. Results showed no difference in global methylation state but a few ploidy-specific loci were detected. DNA methylation variation during polyploidization was less than autonomous methylation variation from parents to offspring.

Variation in genomic methylation in natural populations of Populus simonii is associated with leaf shape and photosynthetic traits

DNA methylation, one of the best-studied types of chromatin modification, suppresses the expression of transposable elements, pseudogenes, repetitive sequences, and individual genes. However, the extent and variation of genome-wide DNA methylation in natural populations of plants remain relatively unknown. To investigate variation in DNA methylation and whether this variation associates with important plant traits, including leaf shape and photosynthesis, 20 413 DNA methylation sites were examined in a poplar association population (505 individuals) using methylation-sensitive amplification polymorphism (MSAP) technology. Calculation of epi-population structure and kinships assigned individuals into subsets (K=3), revealing that the natural population of P. simonii consists of three subpopulations. Population epigenetic distance and geographic distance showed a significant correlation (r=0.4688, P<0.001), suggesting that environmental factors may affect epigenetics. Single-marker approaches were also used to identify significant marker-trait associations, which found 1087 high-confidence DNA methylation markers associated with different phenotypic traits explaining ~5-15% of the phenotypic variance. Among these loci, 147 differentially methylated fragments were obtained by sequencing, representing 130 candidate genes. Expression analysis of six candidate genes indicated that these genes might play important roles in leaf development and regulation of photosynthesis. This study provides association analysis to study the effects of DNA methylation on plant development and these data indicate that epigenetics bridges environmental and genetic factors in affecting plant growth and development.

Rapid genetic divergence in response to 15 years of simulated climate change

Genetic diversity may play an important role in allowing individual species to resist climate change, by permitting evolutionary responses. Our understanding of the potential for such responses to climate change remains limited, and very few experimental tests have been carried out within intact ecosystems. Here, we use amplified fragment length polymorphism (AFLP) data to assess genetic divergence and test for signatures of evolutionary change driven by long-term simulated climate change applied to natural grassland at Buxton Climate Change Impacts Laboratory (BCCIL). Experimental climate treatments were applied to grassland plots for 15 years using a replicated and spatially blocked design and included warming, drought and precipitation treatments. We detected significant genetic differentiation between climate change treatments and control plots in two coexisting perennial plant study species (Festuca ovina and Plantago lanceolata). Outlier analyses revealed a consistent signature of selection associated with experimental climate treatments at individual AFLP loci in P. lanceolata, but not in F. ovina. Average background differentiation at putatively neutral AFLP loci was close to zero, and genomewide genetic structure was associated neither with species abundance changes (demography) nor with plant community-level responses to long-term climate treatments. Our results demonstrate genetic divergence in response to a suite of climatic environments in reproductively mature populations of two perennial plant species and are consistent with an evolutionary response to climatic selection in P. lanceolata. These genetic changes have occurred in parallel with impacts on plant community structure and may have contributed to the persistence of individual species through 15 years of simulated climate change at BCCIL.

Variation in Heavy Metal Accumulation and Genetic Diversity at a Regional Scale Among Metallicolous and Non-Metallicolous Populations of the Facultative Metallophyte Biscutella laevigata subsp. laevigata

Biscutella laevigata is a facultative metallophyte, with populations on non-metalliferous and metalliferous soils. Some of its metallicolous populations have been shown to hyperaccumulate thallium or lead in nature. Only Tl hyperaccumulation has been experimentally confirmed. We aimed to compare the patterns of metal (hyper)accumulation and genetic diversity among populations of B. laevigata subsp. laevigata in NE Italy. None of the populations exhibited foliar hyperaccumulation of Cu, Zn, or Pb. The root-to-shoot accumulation rates for these metals were unchanged or decreased rather than enhanced in the metallicolous populations, in comparison with the non-metallicolous ones. Hyperaccumulation of Tl was confined to the population of the Cave del Predil mine. This population was genetically very distinct from the others, as demonstrated by AFLP-based cluster analysis. The two other mine populations did not surpass the threshold for Tl hyperaccumulation, but showed enhanced foliar Tl concentrations and root-to-shoot translocation rates, in comparison with the non-metallicolous populations. Genetic analysis suggested that adaptation to metalliferous soil must have been independently evolved in the metallicolous populations.

Rapid divergence of ecotypes of an invasive plant

Invasive species demonstrate rapid evolution within a very short period of time allowing one to understand the underlying mechanism(s). Lantana camara, a highly invasive plant of the tropics and subtropics, has expanded its range and successfully established itself almost throughout India. In order to uncover the processes governing the invasion dynamics, 218 individuals from various locations across India were characterized with six microsatellites. By integrating genetic data with niche modelling, we examined the effect of drift and environmental selection on genetic divergence. We found multiple genetic clusters that were non-randomly distributed across space. Spatial autocorrelation revealed a strong fine-scale structure, i.e. isolation by distance. In addition, we obtained evidence of inhibitory effects of selection on gene flow, i.e. isolation by environmental distance. Perhaps, local adaptation in response to selection is offsetting gene flow and causing the populations to diverge. Niche models suggested that temperature and precipitation play a major role in the observed spatial distribution of this plant. Based on a non-random distribution of clusters, unequal gene flow among them and different bioclimatic niche requirements, we concluded that the emergence of ecotypes represented by two genetic clusters is underway. They may be locally adapted to specific climatic conditions, and perhaps at the very early stages of ecological divergence.

Genetic diversity and structure of Dalmatian pyrethrum (Tanacetum cinerariifolium Trevir./Sch./Bip., Asteraceae) within the Balkan refugium

Dalmatian pyrethrum (Tanacetum cinerariifolium Trevir. /Sch./ Bip.) is an outcrossing, perennial insecticidal plant, restricted to the eastern Adriatic coast (Mediterranean). Amplified fragment-length polymorphisms (AFLP) were used to investigate the genetic diversity and structure within and among 20 natural plant populations. The highest level of gene diversity, the number of private alleles and the frequency down-weighted marker values (DW) were found in northern Adriatic populations and gradually decreased towards the southern boundary of the species range. Genetic impoverishment of these southern populations is most likely the result of human-related activities. An analysis of molecular variance (AMOVA) indicated that most of the genetic diversity was attributed to differences among individuals within populations (85.78%), which are expected due to the outcrossing nature of the species. A Bayesian analysis of the population structure identified two dominant genetic clusters. A spatial analysis of the genetic diversity indicated that 5.6% of the genetic differentiation resulted from isolation by distance (IBD), while 12.3% of the genetic differentiation among populations followed the pattern of isolation by environmental distance (IBED). Knowledge of the genetic diversity patterns of the natural populations and the mechanism behind these patterns is required for the exploitation and possible conservation management of this endemic and economically important species.

The importance of species traits for species distribution on oceanic islands

Understanding species' ability to colonize new habitats is a key knowledge allowing us to predict species' survival in the changing landscapes. However, most studies exploring this topic observe distribution of species in landscapes which are under strong human influence being fragmented only recently and ignore the fact that the species distribution in these landscapes is far from equilibrium. Oceanic islands seem more appropriate systems for studying the relationship between species traits and its distribution as they are fragmented without human contribution and as they remained unchanged for a long evolutionary time. In our study we compared the values of dispersal as well as persistence traits among 18 species pairs from the Canary Islands differing in their distribution within the archipelago. The data were analyzed both with and without phylogenetic correction. The results demonstrate that no dispersal trait alone can explain the distribution of the species in the system. They, however, also suggest that species with better dispersal compared to their close relatives are better colonizers. Similarly, abundance of species in the archipelago seems to be an important predictor of species colonization ability only when comparing closely related species. This implies that analyses including phylogenetic correction may provide different insights than analyses without such a correction and both types of analyses should be combined to understand the importance of various plant traits for species colonization ability.

Contrasting evolutionary trajectories of multiple retrotransposons following independent allopolyploidy in wild wheats

Transposable elements (TEs) are expectedly central to genome evolution. To assess the impact of TEs in driving genome turnover, we used allopolyploid genomes, showing considerable deviation from the predicted additivity of their diploid progenitors and thus having undergone major restructuring. Genome survey sequencing was used to select 17 putatively active families of long terminal repeat retrotransposons. Genome-wide TE insertions were genotyped with sequence-specific amplified polymorphism (SSAP) in diploid progenitors and their derived polyploids, and compared with changes in random sequences to assess restructuring of four independent Aegilops allotetraploid genomes. Generally, TEs with different evolutionary trajectories from those of random sequences were identified. Thus, TEs presented family-specific and species-specific dynamics following polyploidy, as illustrated by Sabine showing proliferation in particular polyploids, but massive elimination in others. Contrasting with that, only a few families (BARE1 and Romani) showed proliferation in all polyploids. Overall, TE divergence between progenitors was strongly correlated with the degree of restructuring in polyploid TE fractions. TE families present evolutionary trajectories that are decoupled from genome-wide changes after allopolyploidy and have a pervasive impact on their restructuring.

Genetics of cryptic speciation within an Arctic mustard, Draba nivalis

Crossing experiments indicate that hybrid sterility barriers frequently have developed within diploid, circumpolar plant species of the genus Draba. To gain insight into the rapid evolution of postzygotic reproductive isolation in this system, we augmented the linkage map of one of these species, D. nivalis, and searched for quantitative trait loci (QTLs) associated with reproductive isolation. The map adds 63 new dominant markers to a previously published dataset of 31 co-dominant microsatellites. These markers include 52 amplified fragment length polymorphisms (AFLPs) and 11 sequence-specific amplified polymorphisms (SSAPs) based on retrotransposon sequence. 22 markers displaying transmission ratio distortion were further included in the map. We resolved eight linkage groups with a total map length of 894 cM. Significant genotype-trait associations, or quantitative trait loci (QTL), were detected for reproductive phenotypes including pollen fertility (4 QTLs), seed set (3 QTLs), flowering time (3 QTLs) and number of flowers (4 QTLs). Observed patterns of inheritance were consistent with the influence of both nuclear-nuclear interactions and chromosomal changes on these traits. All seed set QTLs and one pollen fertility QTL displayed underdominant effects suggestive of the involvement of chromosomal rearrangements in hybrid sterility. Interestingly, D. nivalis is predominantly self-fertilizing, which may facilitate the establishment of underdominant loci and contribute to reproductive isolation.

Variation in genomic methylation in natural populations of chinese white poplar

Background

It is thought that methylcytosine can be inherited through meiosis and mitosis, and that epigenetic variation may be under genetic control or correlation may be caused by neutral drift. However, DNA methylation also varies with tissue, developmental stage, and environmental factors. Eliminating these factors, we analyzed the levels and patterns, diversity and structure of genomic methylcytosine in the xylem of nine natural populations of Chinese white poplar.

Principal Findings

On average, the relative total methylation and non-methylation levels were approximately 26.567% and 42.708% (P<0.001), respectively. Also, the relative CNG methylation level was higher than the relative CG methylation level. The relative methylation/non-methylation levels were significantly different among the nine natural populations. Epigenetic diversity ranged from 0.811 (Gansu) to 1.211 (Shaanxi), and the coefficients of epigenetic differentiation (G_ST = 0.159) were assessed by Shannon’s diversity index. Co-inertia analysis indicated that methylation-sensitive polymorphism (MSP) and genomic methylation pattern (CG-CNG) profiles gave similar distributions. Using a between-group eigen analysis, we found that the Hebei and Shanxi populations were independent of each other, but the Henan population intersected with the other populations, to some degree.

Conclusions

Genome methylation in Populus tomentosa presented tissue-specific characteristics and the relative 5′-CCGG methylation level was higher in xylem than in leaves. Meanwhile, the genome methylation in the xylem shows great epigenetic variation and could be fixed and inherited though mitosis. Compared to genetic structure, data suggest that epigenetic and genetic variation do not completely match.

Are adaptive loci transferable across genomes of related species? Outlier and environmental association analyses in Alpine Brassicaceae species

Local adaptation is one possible response of organisms to survive in a changing environment. However, the genetic basis of adaptation is not well understood, especially in nonmodel species. To infer recurrent patterns of local adaptation, we investigated whether the same putative adaptive loci reoccur in related species. We performed genome scans using amplified fragment length polymorphism (AFLP) markers on populations of five Alpine Brassicaceae species sampled across a wide range of environmental conditions. To identify markers potentially under directional selection, we performed outlier and environmental association analyses using a set of topo-climatic variables available as GIS layers. Several AFLP loci showed signatures of adaptation, of which one, found in Cardamine resedifolia (Cre_P1_212.5), was associated with precipitation. We sequence-characterized this candidate locus and genotyped single nucleotide polymorphisms (SNPs) found within this locus for all species. Testing for environmental associations of SNPs revealed the same association of this locus in Arabis alpina but not in other study species. Cumulative statistical evidence indicates that locus Cre_P1_212.5 is environmentally relevant or is linked to a gene under selection in our study range. Furthermore, the locus shows an association to the same potentially selective factor in at least one other related species. These findings help to identify trends in plant adaptation in Alpine ecosystems in response to particular environmental parameters.

Environmental heterogeneity explains the genetic structure of Continental and Mediterranean populations of Fraxinus angustifolia Vahl

Tree species with wide distributions often exhibit different levels of genetic structuring correlated to their environment. However, understanding how environmental heterogeneity influences genetic variation is difficult because the effects of gene flow, drift and selection are confounded. We investigated the genetic variation and its ecological correlates in a wind-pollinated Mediterranean tree species, Fraxinus angustifolia Vahl, within a recognised glacial refugium in Croatia. We sampled 11 populations from environmentally divergent habitats within the Continental and Mediterranean biogeographical regions. We combined genetic data analyses based on nuclear microsatellite loci, multivariate statistics on environmental data and ecological niche modelling (ENM). We identified a geographic structure with a high genetic diversity and low differentiation in the Continental region, which contrasted with the significantly lower genetic diversity and higher population divergence in the Mediterranean region. The positive and significant correlation between environmental and genetic distances after controlling for geographic distance suggests an important influence of ecological divergence of the sites in shaping genetic variation. The ENM provided support for niche differentiation between the populations from the Continental and Mediterranean regions, suggesting that contemporary populations may represent two divergent ecotypes. Ecotype differentiation was also supported by multivariate environmental and genetic distance analyses. Our results suggest that despite extensive gene flow in continental areas, long-term stability of heterogeneous environments have likely promoted genetic divergence of ashes in this region and can explain the present-day genetic variation patterns of these ancient populations.

Landscape genetics of plants

Landscape genetics is the amalgamation of landscape ecology and population genetics to help with understanding microevolutionary processes such as gene flow and adaptation. In this review, we examine why landscape genetics of plants lags behind that of animals, both in number of studies and consideration of landscape elements. The classical landscape distance/resistance approach to study gene flow is challenging in plants, whereas boundary detection and the assessment of contemporary gene flow are more feasible. By contrast, the new field of landscape genetics of adaptive genetic variation, establishing the relationship between adaptive genomic regions and environmental factors in natural populations, is prominent in plant studies. Landscape genetics is ideally suited to study processes such as migration and adaptation under global change.

Tracking genes of ecological relevance using a genome scan in two independent regional population samples of Arabis alpina

Understanding the genetic basis of adaptation in response to environmental variation is fundamental as adaptation plays a key role in the extension of ecological niches to marginal habitats and in ecological speciation. Based on the assumption that some genomic markers are correlated to environmental variables, we aimed to detect loci of ecological relevance in the alpine plant Arabis alpina L. sampled in two regions, the French (99 locations) and the Swiss (109 locations) Alps. We used an unusually large genome scan [825 amplified fragment length polymorphism loci (AFLPs)] and four environmental variables related to temperature, precipitation and topography. We detected linkage disequilibrium among only 3.5% of the considered AFLP loci. A population structure analysis identified no admixture in the study regions, and the French and Swiss Alps were differentiated and therefore could be considered as two independent regions. We applied generalized estimating equations (GEE) to detect ecologically relevant loci separately in the French and Swiss Alps. We identified 78 loci of ecological relevance (9%), which were mainly related to mean annual minimum temperature. Only four of these loci were common across the French and Swiss Alps. Finally, we discuss that the genomic characterization of these ecologically relevant loci, as identified in this study, opens up new perspectives for studying functional ecology in A. alpina, its relatives and other alpine plant species.

Patterns of population structure and environmental associations to aridity across the range of loblolly pine (Pinus taeda L., Pinaceae)

Natural populations of forest trees exhibit striking phenotypic adaptations to diverse environmental gradients, thereby making them appealing subjects for the study of genes underlying ecologically relevant phenotypes. Here, we use a genome-wide data set of single nucleotide polymorphisms genotyped across 3059 functional genes to study patterns of population structure and identify loci associated with aridity across the natural range of loblolly pine (Pinus taeda L.). Overall patterns of population structure, as inferred using principal components and Bayesian cluster analyses, were consistent with three genetic clusters likely resulting from expansions out of Pleistocene refugia located in Mexico and Florida. A novel application of association analysis, which removes the confounding effects of shared ancestry on correlations between genetic and environmental variation, identified five loci correlated with aridity. These loci were primarily involved with abiotic stress response to temperature and drought. A unique set of 24 loci was identified as F(ST) outliers on the basis of the genetic clusters identified previously and after accounting for expansions out of Pleistocene refugia. These loci were involved with a diversity of physiological processes. Identification of nonoverlapping sets of loci highlights the fundamental differences implicit in the use of either method and suggests a pluralistic, yet complementary, approach to the identification of genes underlying ecologically relevant phenotypes.

Epigenetic variation in mangrove plants occurring in contrasting natural environment

BACKGROUND

Epigenetic modifications, such as cytosine methylation, are inherited in plant species and may occur in response to biotic or abiotic stress, affecting gene expression without changing genome sequence. Laguncularia racemosa, a mangrove species, occurs in naturally contrasting habitats where it is subjected daily to salinity and nutrient variations leading to morphological differences. This work aims at unraveling how CpG-methylation variation is distributed among individuals from two nearby habitats, at a riverside (RS) or near a salt marsh (SM), with different environmental pressures and how this variation is correlated with the observed morphological variation.

PRINCIPAL FINDINGS

Significant differences were observed in morphological traits such as tree height, tree diameter, leaf width and leaf area between plants from RS and SM locations, resulting in smaller plants and smaller leaf size in SM plants. Methyl-Sensitive Amplified Polymorphism (MSAP) was used to assess genetic and epigenetic (CpG-methylation) variation in L. racemosa genomes from these populations. SM plants were hypomethylated (14.6% of loci had methylated samples) in comparison to RS (32.1% of loci had methylated samples). Within-population diversity was significantly greater for epigenetic than genetic data in both locations, but SM also had less epigenetic diversity than RS. Frequency-based (G(ST)) and multivariate (beta(ST)) methods that estimate population structure showed significantly greater differentiation among locations for epigenetic than genetic data. Co-Inertia analysis, exploring jointly the genetic and epigenetic data, showed that individuals with similar genetic profiles presented divergent epigenetic profiles that were characteristic of the population in a particular environment, suggesting that CpG-methylation changes may be associated with environmental heterogeneity.

CONCLUSIONS

In spite of significant morphological dissimilarities, individuals of L. racemosa from salt marsh and riverside presented little genetic but abundant DNA methylation differentiation, suggesting that epigenetic variation in natural plant populations has an important role in helping individuals to cope with different environments.

Evolutionary consequences of autopolyploidy

Autopolyploidy is more common in plants than traditionally assumed, but has received little attention compared with allopolyploidy. Hence, the advantages and disadvantages of genome doubling per se compared with genome doubling coupled with hybridizations in allopolyploids remain unclear. Autopolyploids are characterized by genomic redundancy and polysomic inheritance, increasing effective population size. To shed light on the evolutionary consequences of autopolyploidy, we review a broad range of studies focusing on both synthetic and natural autopolyploids encompassing levels of biological organization from genes to evolutionary lineages. The limited evidence currently available suggests that autopolyploids neither experience strong genome restructuring nor wide reorganization of gene expression during the first generations following genome doubling, but that these processes may become more important in the longer term. Biogeographic and ecological surveys point to an association between the formation of autopolyploid lineages and environmental change. We thus hypothesize that polysomic inheritance may provide a short-term evolutionary advantage for autopolyploids compared to diploid relatives when environmental change enforces range shifts. In addition, autopolyploids should possess increased genome flexibility, allowing them to adapt and persist across heterogeneous landscapes in the long run.

Divergent selection in trailing- versus leading-edge populations of Biscutella laevigata

BACKGROUND AND AIMS

Knowledge on how climate-induced range shifts might affect natural selection is crucial to understand the evolution of species ranges.

METHODS

Using historical demographic perspectives gathered from regional-scale phylogeography on the alpine herb Biscutella laevigata, indirect inferences on gene flow and signature of selection based on AFLP genotyping were compared between local populations persisting at the trailing edge and expanding at the leading edge.

KEY RESULTS

Spatial autocorrelation revealed that gene flow was two times more restricted at the trailing edge and genome scans indicated divergent selection in this persisting population. In contrast, no pattern of selection emerged in the expanding population at the leading edge.

CONCLUSIONS

Historical effects may determine different architecture of genetic variation and selective patterns within local populations, what is arguably important to understand evolutionary processes acting across the species ranges.

Pleistocene distribution range shifts were accompanied by breeding system divergence within Hornungia alpina (Brassicaceae) in the Alps

Impact of glacial history on the phylogeography of silicate-dwelling plants of the European Alps has been particularly well studied, whereas virtually no data are available for species growing on different bedrock types, as for Hornungia alpina. Bayesian clustering of AFLP data only partly support the distinction of three subspecies as morphologically defined. Whereas the phylogeographical N-group corresponds to subsp. alpina, the congruence of the SW-group and SE-group with subsp. brevicaulis, and subsp. austroalpina, respectively, is limited. High levels of rarity and genetic diversity in the N-group suggest Pleistocene survival along the outer margin of the Alpine arc. For subsp. brevicaulis we suggest a single origin from a refugium in the Southwestern Alps, whereas subsp. austroalpina might have originated twice in the Southern and Southeastern Alps. Different levels of genetic diversity and partitioning of genetic variation indicate a divergence in breeding system, which is corroborated by pollinator exclusion experiments revealing self-incompatibility in the N-group and autonomous selfing in the SE-group.

Darwin's wind hypothesis: does it work for plant dispersal in fragmented habitats?

Using the wind-dispersed plant Mycelis muralis, we examined how landscape fragmentation affects variation in seed traits contributing to dispersal. Inverse terminal velocity (Vt(-1)) of field-collected achenes was used as a proxy for individual seed dispersal ability. We related this measure to different metrics of landscape connectivity, at two spatial scales: in a detailed analysis of eight landscapes in Spain and along a latitudinal gradient using 29 landscapes across three European regions. In the highly patchy Spanish landscapes, seed Vt(-1)increased significantly with increasing connectivity. A common garden experiment suggested that differences in Vt(-1) may be in part genetically based. The Vt(-1) was also found to increase with landscape occupancy, a coarser measure of connectivity, on a much broader (European) scale. Finally, Vt(-1)was found to increase along a south-north latitudinal gradient. Our results for M. muralis are consistent with 'Darwin's wind dispersal hypothesis' that high cost of dispersal may select for lower dispersal ability in fragmented landscapes, as well as with the 'leading edge hypothesis' that most recently colonized populations harbour more dispersive phenotypes.

Genetic markers in the playground of multivariate analysis

Multivariate analyses such as principal component analysis were among the first statistical methods employed to extract information from genetic markers. From their early applications to current innovations, these approaches have proven to be efficient for the analysis of the genetic variability in various contexts such as human genetics, conservation and adaptation studies. However, because multivariate analysis is a wide and diversified area of statistics, choosing a method appropriate to both the data and to the question being asked can be difficult. Moreover, some particularities of genetic markers need to be taken into account when using multivariate methods. As a consequence, multivariate analyses are often used as black boxes, which results in frequent mistakes in the literature. In this review, we provide a critical analysis of the application of multivariate methods to genetic markers, using a general framework that unifies all these methods for the sake of clarity. First, we focus on some common mistakes in these applications and ways to avoid these pitfalls. We then detail the most critical particularities of allele frequencies that demand adaptations of multivariate methods, and we propose solutions to the subsequent problems. Finally, we tackle several questions of interest in which multivariate analysis has a great role to play, such as the study of the typological coherence of different genetic markers, or the investigation of spatial genetic patterns.

Fine-scale genetic structure and marginal processes in an expanding population of Biscutella laevigata L. (Brassicaceae)

Evolutionary processes acting at the expanding margins of a species' range are still poorly understood. Genetic drift is considered prevalent in marginal populations, and the maintenance of genetic diversity during recolonization might seem puzzling. To investigate such processes, a fine-scale investigation of 219 individuals was performed within a population of Biscutella laevigata (Brassicaceae), located at the leading edge of its range. The survey used amplified fragment length polymorphisms (AFLPs). As commonly reported across the whole species distribution range, individual density and genetic diversity decreased along the local axis of recolonization of this expanding population, highlighting the enduring effect of the historical colonization on present-day diversity. The self-incompatibility system of the plant may have prevented local inbreeding in newly found patches and sustained genetic diversity by ensuring gene flow from established populations. Within the more continuously populated region, spatial analysis of genetic structure revealed restricted gene flow among individuals. The distribution of genotypes formed a mosaic of relatively homogenous patches within the continuous population. This pattern could be explained by a history of expansion by long-distance dispersal followed by fine-scale diffusion (that is, a stratified dispersal combination). The secondary contact among expanding patches apparently led to admixture among differentiated genotypes where they met (that is, a reshuffling effect). This type of dynamics could explain the maintenance of genetic diversity during recolonization.