AFLP genome scan to detect genetic structure and candidate loci under selection for local adaptation of the invasive weed Mikania micrantha

Epigenetic architecture of Pseudotaxus chienii: Revealing the synergistic effects of climate and soil variables

Whether conifers can withstand environmental changes especially temperature fluctuations has been controversial. Epigenetic analysis may provide new perspectives for solving the issue. Pseudotaxus chienii is an endangered gymnosperm species endemic to China. In this study, we have examined the genetic and epigenetic variations in its natural populations aiming to disentangle the synergistic effects of climate and soil on its population (epi)genetic differentiation by using amplified fragment length polymorphism (AFLP) and methylation-sensitive AFLP (MSAP) techniques. We identified 23 AFLP and 26, 7, and 5 MSAP outliers in P. chienii. Twenty-one of the putative adaptive AFLP loci were found associated with climate and/or soil variables including precipitation, temperature, K, Fe, Zn, and Cu, whereas 21, 7, and 4 MSAP outliers were significantly related to precipitation of wettest month (Bio13), precipitation driest of month (Bio14), percent tree cover (PTC), and soil Fe, Mn, and Cu compositions. Total precipitation and precipitation in the driest seasons were the most influential factors for genetic and epigenetic variation, respectively. In addition, a high full-methylation level and a strong correlation between genetic and epigenetic variation were detected in P. chienii. Climate is found of greater importance than soil in shaping adaptive (epi)genetic differentiation, and the synergistic effects of climate and climate-soil variables were also observed. The identified climate and soil variables should be considered when applying ex situ conservation.

Integrated Analysis of Transcriptome and Small RNAome Reveals the Regulatory Network for Rapid Growth in Mikania micrantha

M. micrantha has caused huge ecological damage and economic losses worldwide due to its rapid growth and serious invasion. However, the underlying molecular mechanisms of its rapid growth and environmental adaption remain unclear. Here, we performed transcriptome and small RNA sequencing with five tissues of M. micrantha to dissect miRNA-mediated regulation in M. micrantha. WGCNA and GO enrichment analysis of transcriptome identified the gene association patterns and potential key regulatory genes for plant growth in each tissue. The genes highly correlated with leaf and stem tissues were mainly involved in the chlorophyll synthesis, response to auxin, the CAM pathway and other photosynthesis-related processes, which promoted the fast growth of M. micrantha. Importantly, we identified 350 conserved and 192 novel miRNAs, many of which displayed differential expression patterns among tissues. PsRNA target prediction analysis uncovered target genes of both conserved and novel miRNAs, including GRFs and TCPs, which were essential for plant growth and development. Further analysis revealed that miRNAs contributed to the regulation of tissue-specific gene expression in M. micrantha, such as mmi-miR396 and mmi-miR319. Taken together, our study uncovered the miRNA-mRNA regulatory networks and the potential vital roles of miRNAs in modulating the rapid growth of M. micrantha.

AFLP-Based Genetic Structure of Lithuanian Populations of Small Balsam (Impatiens parviflora DC.) in Relation to Habitat Characteristics

Currently, there is an increasing focus on understanding the interactions between genetic features of the invader and environmental factors that ensure the success of invasion. The objective of our study was to evaluate the genetic diversity of Lithuanian populations of highly invasive small balsam (Impatiens parviflora) by amplified fragment length polymorphism (AFLP) markers and to relate molecular data to biotope features defined by employing neighboring species of herbaceous plants. Low polymorphism of I. parviflora populations was observed at AFLP loci. Hierarchical analysis of molecular variance did not reveal differentiation of populations depending on biotope, geography, or road types. Bayesian analyses of AFLP data demonstrated many genetic clusters. Our results suggest multiple introductions of I. parviflora into Lithuania. The polymorphism of AFLP loci of populations significantly correlated with the total coverage by herbaceous plants in the sites. Defined by principal component analysis, the variability of study sites was most related to the coverage of herbaceous plants and least related to the molecular features of I. parviflora populations. The sites with I. parviflora were classified into agricultural scrubland, riparian forest, and urban forest biotopes. Of them, urban forest was distinguished by the highest coverage of I. parviflora and the lowest Ellenberg indicatory values for light, soil acidity, and richness in nutrients.

Population structure and adaptive variation of Helichrysum italicum (Roth) G. Don along eastern Adriatic temperature and precipitation gradient

Immortelle (Helichrysum italicum (Roth) G. Don; Asteraceae) is a perennial plant species native to the Mediterranean region, known for many properties with wide application mainly in perfume and cosmetic industry. A total of 18 wild H. italicum populations systematically sampled along the eastern Adriatic environmental gradient were studied using AFLP markers to determine genetic diversity and structure and to identify loci potentially responsible for adaptive divergence. Results showed higher levels of intrapopulation diversity than interpopulation diversity. Genetic differentiation among populations was significant but low, indicating extensive gene flow between populations. Bayesian analysis of population structure revealed the existence of two genetic clusters. Combining the results of F_ST - outlier analysis (Mcheza and BayeScan) and genome-environment association analysis (Samβada, LFMM) four AFLP loci strongly associated with the bioclimatic variables Bio03 Isothermality, Bio08 Mean temperature of the wettest quarter, Bio15 Precipitation seasonality, and Bio17 Precipitation of driest quarter were found to be the main variables driving potential adaptive genetic variation in H. italicum along the eastern Adriatic environmental gradient. Redundancy analysis revealed that the partitioning of genetic variation was mainly associated with the adaptation to temperature oscillations. The results of the research may contribute to a clearer understanding of the importance of local adaptations for the genetic differentiation of Mediterranean plants and allow the planning of appropriate conservation strategies. However, considering that the identified outlier loci may be linked to genes under selection rather than being the target of natural selection, future studies must aim at their additional analysis.

High genetic and epigenetic variation of transposable elements: Potential drivers to rapid adaptive evolution for the noxious invasive weed Mikania micrantha

Why invasive species can rapidly adapt to novel environments is a puzzling question known as the genetic paradox of invasive species. This paradox is explainable in terms of transposable elements (TEs) activity, which are theorized to be powerful mutational forces to create genetic variation. Mikania micrantha, a noxious invasive weed, in this sense provides an excellent opportunity to test the explanation. The genetic and epigenetic variation of 21 invasive populations of M. micrantha in southern China have been examined by using transposon display (TD) and transposon methylation display (TMD) techniques to survey 12 TE superfamilies. Our results showed that M. micrantha populations maintained an almost equally high level of TE-based genetic and epigenetic variation and they have been differentiated into subpopulations genetically and epigenetically. A similar positive spatial genetic and epigenetic structure pattern was observed within 300 m. Six and seven TE superfamilies presented significant genetic and epigenetic isolation by distance (IBD) pattern. In total, 59 genetic and 86 epigenetic adaptive TE loci were identified. Of them, 51 genetic and 44 epigenetic loci were found to correlate with 25 environmental variables (including precipitation, temperature, vegetation coverage, and soil metals). Twenty-five transposon-inserted genes were sequenced and homology-based annotated, which are found to be involved in a variety of molecular and cellular functions. Our research consolidates the importance of TE-associated genetic and epigenetic variation in the rapid adaptation and invasion of M. micrantha.

Population Genomics Reveals Gene Flow and Adaptive Signature in Invasive Weed Mikania micrantha

A long-standing and unresolved issue in invasion biology concerns the rapid adaptation of invaders to nonindigenous environments. Mikania micrantha is a notorious invasive weed that causes substantial economic losses and negative ecological consequences in southern China. However, the contributions of gene flow, environmental variables, and functional genes, all generally recognized as important factors driving invasive success, to its successful invasion of southern China are not fully understood. Here, we utilized a genotyping-by-sequencing approach to sequence 306 M. micrantha individuals from 21 invasive populations. Based on the obtained genome-wide single nucleotide polymorphism (SNP) data, we observed that all the populations possessed similar high levels of genetic diversity that were not constrained by longitude and latitude. Mikania micrantha was introduced multiple times and subsequently experienced rapid-range expansion with recurrent high gene flow. Using FST outliers, a latent factor mixed model, and the Bayesian method, we identified 38 outlier SNPs associated with environmental variables. The analysis of these outlier SNPs revealed that soil composition, temperature, precipitation, and ecological variables were important determinants affecting the invasive adaptation of M. micrantha. Candidate genes with outlier signatures were related to abiotic stress response. Gene family clustering analysis revealed 683 gene families unique to M. micrantha which may have significant implications for the growth, metabolism, and defense responses of M. micrantha. Forty-one genes showing significant positive selection signatures were identified. These genes mainly function in binding, DNA replication and repair, signature transduction, transcription, and cellular components. Collectively, these findings highlight the contribution of gene flow to the invasion and spread of M. micrantha and indicate the roles of adaptive loci and functional genes in invasive adaptation.

Population Genetics of Odontarrhena (Brassicaceae) from Albania: The Effects of Anthropic Habitat Disturbance, Soil, and Altitude on a Ni-Hyperaccumulator Plant Group from a Major Serpentine Hotspot

Albanian taxa and populations of the genus Odontarrhena are most promising candidates for research on metal tolerance and Ni-agromining, but their genetic structure remains unknown. We investigated phylogenetic relationships and genetic differentiation in relation to distribution and ploidy of the taxa, anthropic site disturbance, elevation, soil type, and trace metals at each population site. After performing DNA sequencing of selected accessions, we applied DNA-fingerprinting to analyze the genetic structure of 32 populations from ultramafic and non-ultramafic outcrops across Albania. Low sequence divergence resulted in poorly resolved phylograms, but supported affinity between the two diploid serpentine endemics O. moravensis and O. rigida. Analysis of molecular variance (AMOVA) revealed significant population differentiation, but no isolation by distance. Among-population variation was higher in polyploids than in diploids, in which genetic distances were lower. Genetic admixing at population and individual level occurred especially in the polyploids O. chalcidica, O. decipiens, and O. smolikana. Admixing increased with site disturbance. Outlier loci were higher in serpentine populations but decreased along altitude with lower drought and heat stress. Genetic variability gained by gene flow and hybridization at contact zones with "resident" species of primary ultramafic habitats promoted expansion of the tetraploid O. chalcidica across anthropogenic sites.

Genetic variability and structure of an important wild steppe grass Psathyrostachys juncea (Triticeae: Poaceae) germplasm collection from north and central Asia

Russian wildrye, Psathyrostachys junceus (Fisch.) Nevski, is an important wild steppe perennial grass, which is characterized by longevity and short robust rhizomes. It also has highly resistance in grazing and abiotic stress. In this study, the genetic diversity of eleven P. juncea wild germplasms from north and central Asia was investigated using AFLP markers. The P. juncea populations were divided into three clades in both UPGMA dendrogram and PCoA clustering corresponding to the three genetic memberships in STRUCTURE analysis. The genetic specificity of Xinjiang (XJ) populations was revealed by the highest N_e (1.5411) and H_j (0.3553) and their dispersion shown in UPGMA. High inbreeding coefficient (Wright’s F statistics, f = 0.496, F_st = 0.128) was observed although a moderate degree of gene flow (N_m = 1.4736) existed, that may ascribe to habitat fragmentation or the low seedling recruitment, which ultimately resulted in decrease of population size and their genetic diversity. The key factors inducing moderate genetic differentiation detected in this study were isolation by distance (IBD), climatic variabilities and geographical barriers. All these results provide insights into the study of genetic status and germplasm collecting of Russian wildrye.

Adaptive fitness of Sapindus emarginatus Vahl populations towards future climatic regimes and the limiting factors of its distribution

Sapindus emarginatus Vahl (Sapindaceae) also known as ‘Indian Soap nut’ is significantly important for saponin content in its fruits. However, its current population in India is heavily fragmented due to a lack of sustainable harvesting practices. Moreover, changing climatic regimes may further limit its distribution and possibly compromise the survival of the species in nature. The aim of the present study was to: predict the future distribution range of S. emarginatus; identify the bioclimatic variables limiting this distribution and to evaluate its adaptive fitness and genomic resilience towards these variables. To determine future species distribution range and identify limiting bioclimatic variables, we applied two different ecological niche models (ENMs; BioClim and MaxEnt) on real occurrence data (n = 88 locations). The adaptive fitness of the species was evaluated by quantifying the genetic variability with AFLP markers and marker-environmental associations, using AFLP-associated Bayesian statistics. We found 77% overlap between the baseline (2030) and predicted (2100) species distribution ranges, which were primarily determined by maximum temperature (T_MAX) and mean annual precipitation (MAP). The T_MAX and MAP contributed 43.1% and 27.1%, respectively to ENM model prediction. Furthermore, AFLP loci significantly associated with bioclimatic variables, and T_MAX and MAP represent the lowest proportion (6.15%), confirming to the severe response of the species genome towards these variables. Nevertheless, the very low Linkage disequilibrium (LD) in these loci (4.54%) suggests that the current sensitivity to T_MAX and MAP is subject to change during recombination. Moreover, a combination of high heterozygosity (0.40–0.43) and high within-population variability (91.63 ± 0.31%) confirmed high adaptive fitness to maintain reproductive success. Therefore, the current populations of S. emarginatus have substantial genomic resilience towards future climate change, albeit significant conservation efforts (including mass multiplication) are warranted to avoid future deleterious impacts of inbreeding depression on the fragmented populations.

Genetic Diversity and Population Divergence of a Rare, Endemic Grass (Elymus breviaristatus) in the Southeastern Qinghai-Tibetan Plateau

Elymus breviaristatus is a grass species only distributed in the southeast of Qinghai-Tibetan Plateau (QTP), which has suffered from serious habitat fragmentation. Therefore, understanding patterns of genetic diversity within and among natural E. breviaristatus populations could provide insight for future conservation strategies. In this study, sequence-related amplified polymorphism markers were employed to investigate the genetic diversity and hierarchical structure of seven E. breviaristatus populations from QTP, China. Multiple measures of genetic diversity indicated that there is low to moderate genetic variation within E. breviaristatus populations, consistent with its presumed mating system. In spite of its rarity, E. breviaristatus presented high genetic diversity that was equivalent to or even higher than that of widespread species. Bayesian clustering approaches, along with clustering analysis and principal coordinate analysis partitioned the studied populations of E. breviaristatus into five genetic clusters. Differentiation coefficients (Fst, GST, etc.) and AMOVA analysis revealed considerable genetic divergence among different populations. BARRIER analyses indicated that there were two potential barriers to gene flow among the E. breviaristatus populations. Despite these patterns of differentiation, genetic distances between populations were independent of geographic distances (r = 0.2197, p = 0.2534), indicating little isolation by distance. Moreover, despite detecting a common outlier by two methods, bioclimatic factors (altitude, annual mean temperature, and annual mean precipitation) were not related to diversity parameters, indicating little evidence for isolation caused by the environment. These patterns of diversity within and between populations are used to propose a conservation strategy for E. breviaristatus.

Divergent selection and genetic structure of Sideritis scardica populations from southern Balkan Peninsula as revealed by AFLP fingerprinting

Sideritis scardica Giseb. is a subalpine/alpine plant species endemic to the central part of the Balkan Peninsula. In this study, we combined Amplified Fragment Length Polymorphism (AFLP) and environmental data to examine the adaptive genetic variations in S. scardica natural populations sampled in contrasting environments. A total of 226 AFLP loci were genotyped in 166 individuals from nine populations. The results demonstrated low gene diversity, ranging from 0.095 to 0.133 and significant genetic differentiation ranging from 0.115 to 0.408. Seven genetic clusters were revealed by Bayesian clustering methods as well as by Discriminant Analysis of Principal Components and each population formed its respective cluster. The exception were populations P02 Mt. Shara and P07 Mt. Vermio, that were admixed between two clusters. Both landscape genetic methods Mcheza and BayeScan identified a total of seven (3.10%) markers exhibiting higher levels of genetic differentiation among populations. The spatial analysis method Samβada detected 50 individual markers (22.12%) associated with bioclimatic variables, among them seven were identified by both Mcheza and BayeScan as being under directional selection. Four bioclimatic variables associated with five out of seven outliers were related to precipitation, suggesting that this variable is the key factor affecting the adaptive variation of S. scardica.

Salt-tolerant native plants have greater responses to other environments when compared to salt-tolerant invasive plants

The strong expansion potential of invasive plants is often attributed to fast adaptive responses to stress. However, the evolution of tolerance to one stressor may affect the responses to other stressors. Currently, it remains unclear what effect the evolution to one stressor might have on the responses to other single or combined stressors. Moreover, it is unknown how this might differ between invasive and native species.Invasive plants (Mikania micrantha and Bidens pilosa) and native plants (Merremia hederacea and Sida acuta) from low- and high-salinity habitats were grown under control and stressful conditions [salt stress, water stress (drought/waterlogging), and their combinations]. We explored the effects of evolved salt tolerance on the responses to water stress/combined stresses and the underlying trait mechanisms.The high-salinity populations of all species exhibited stronger salt tolerance than the low-salinity populations. As to the tolerance to other stressors, the high-salinity and low-salinity populations of the invasive species were similar, whereas the high-salinity populations of the native species exhibited stronger tolerance than the low-salinity populations under most stress treatments. However, the enhanced salt tolerance in native species was accompanied by reduced total biomass under control condition. The stress tolerance of native species correlated with leaf production rate and allocation to root, while the performance of native species under control condition correlated with leaf morphology and carbon assimilation rate. This suggests a trade-off between salt tolerance and performance in the native but not the invasive species, probably resulting from altered phenotypic/physiological traits.

SYNTHESIS

Our work suggests that the evolution of tolerance to one stressor may have stronger effects on the tolerance to other stressors of the native compared with the invasive species. This may be a new paradigm to explain the greater advantage of invasive vs. native species in highly stressful habitats.

Comparative transcriptome analysis of the invasive weed Mikania micrantha with its native congeners provides insights into genetic basis underlying successful invasion

Background

Mikania micrantha H.B.K. (Asteraceae) is one of the world’s most invasive weeds which has been rapidly expanding in tropical Asia, including China, while its close relative M. cordata, the only Mikania species native to China, shows no harm to the local ecosystems. These two species are very similar in morphology but differ remarkably in several ecological and physiological traits, representing an ideal system for comparative analysis to investigate the genetic basis underlying invasion success. In this study, we performed RNA-sequencing on the invader M. micrantha and its native congener M. cordata in China, to unravel the genetic basis underlying the strong invasiveness of M. micrantha. For a more robust comparison, another non-invasive congener M. cordifolia was also sequenced and compared.

Results

A total of 52,179, 55,835, and 52,983 unigenes were obtained for M. micrantha, M. cordata, and M. cordifolia, respectively. Phylogenetic analyses and divergence time dating revealed a relatively recent split between M. micrantha and M. cordata, i.e., approximately 4.81 million years ago (MYA), after their divergence with M. cordifolia (8.70 MYA). Gene ontology classifications, pathway assignments and differential expression analysis revealed higher representation or significant up-regulation of genes associated with photosynthesis, energy metabolism, protein modification and stress response in M. micrantha than in M. cordata or M. cordifolia. Analysis of accelerated evolution and positive selection also suggested the importance of these related genes and processes to the adaptability and invasiveness of M. micrantha. Particularly, most (77 out of 112, i.e. 68.75%) positively selected genes found in M. micrantha could be classified into four groups, i.e., energy acquisition and utilization (10 genes), growth and reproduction (13 genes), protection and repair (34 genes), and signal transduction and expression regulation (20 genes), which may have contributed to the high adaptability of M. micrantha to various new environments and the capability to occupy a wider niche, reflected in its high invasiveness.

Conclusions

We characterized the transcriptomes of the invasive species M. micrantha and its non-invasive congeners, M. cordata and M. cordifolia. A comparison of their transcriptomes provided insights into the genetic basis of the high invasiveness of M. micrantha.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4784-9) contains supplementary material, which is available to authorized users.

Comparative chloroplast genomics between the invasive weed Mikania micrantha and its indigenous congener Mikania cordata: Structure variation, identification of highly divergent regions, divergence time estimation, and phylogenetic analysis

Mikania micrantha and Mikania cordata are the only two species in genus Mikania (Asteraceae) in China. They share very similar morphological and life-history characteristics but occupy quite different habitats. Most importantly, they generate totally different ecological consequences. While M. micrantha has become an exotic invasive weed, M. cordata exists as an indigenous species with no harmful effects on native plants or habitats. As a continuous study of our previously reported M. micrantha chloroplast (cp) genome, in this study we have further sequenced the M. cordata cp genome to (1) conduct a comparative genome analysis to gain insights into the mechanism of invasiveness; (2) develop cp markers to examine the population genetic adaptation of M. micrantha; and (3) screen variable genome regions of phylogenetic utility. The M. cordata chloroplast genome is 151,984 bp in length and displays a typical quadripartite structure. The number and distribution of protein coding genes, tRNA genes, and rRNA genes of M. cordata are identical to those of M. micrantha. The main difference lays in that the pseudogenization of ndhF and a 118-bp palindromic repeat only arises in M. cordata. Fourteen highly divergent regions, 235 base substitutions, and 58 indels were identified between the two cp genomes. Phylogenetic inferences revealed a sister relationship between M. micrantha and M. cordata whose divergence was estimated to occur around 1.78 million years ago (MYA). Twelve cpSSR loci were detected to be polymorphic and adopted to survey the genetic adaptation of M. micrantha populations. No cpSSR loci were found to undergo selection. Our results build a foundation to examine the invasive mechanism of Mikania weed.

The emergence of the hyperinvasive vine, Mikania micrantha (Asteraceae), via admixture and founder events inferred from population transcriptomics

Biological invasions that involve well-documented rapid adaptations to new environments provide unequalled opportunities for testing evolutionary hypotheses. Mikania micrantha Kunth (Asteraceae), a perennial herbaceous vine native to tropical Central and South America, successfully invaded tropical Asia in the early 20th century. It is regarded as one of the most aggressive weeds in the world. To elucidate the molecular and evolutionary processes underlying this invasion, we extensively sampled this weed throughout its invaded range in South-East and South Asia and surveyed its genetic structure using variants detected from population transcriptomics. Clustering results suggest that more than one source population contributed to this invasion. Computer simulations using genomewide genetic variation support a scenario of admixture and founder events during invasion. The genes differentially expressed between native and invasive populations were found to be involved in oxidative and high light intensity stress responses, pointing to a possible ecological mechanism of adaptation. Our results provide a foundation for further detailed mechanistic and population studies of this ecologically and economically important invasion. This line of research promises to provide new mitigation strategies for invasive species as well as insights into mechanisms of adaptation.

Invasive Chloroplast Population Genetics of Mikania micrantha in China: No Local Adaptation and Negative Correlation between Diversity and Geographic Distance

Two fundamental questions on how invasive species are able to rapidly colonize novel habitat have emerged. One asks whether a negative correlation exists between the genetic diversity of invasive populations and their geographic distance from the origin of introduction. The other is whether selection on the chloroplast genome is important driver of adaptation to novel soil environments. Here, we addressed these questions in a study of the noxious invasive weed, Mikania micrantha, which has rapidly expanded in to southern China after being introduced to Hong Kong in 1884. Seven chloroplast simple sequence repeats (cpSSRs) were used to investigate population genetics in 28 populations of M. micrantha, which produced 39 loci. The soil compositions for these populations, including Mg abundance, were measured. The results showed that M. micrantha possessed relatively high cpSSR variation and differentiation among populations. Multiple diversity indices were quantified, and none was significantly correlated with distance from the origin of introduction. No evidence for "isolation by distance," significant spatial structure, bottlenecks, nor linkage disequilibrium was detected. We also were unable to identify loci on the chloroplast genome that exhibited patterns of differentiation that would suggest adaptive evolution in response to soil attributes. Soil Mg had only a genome-wide effect instead of being a selective factor, which highlighted the association between Mg and the successful invasion. This study characterizes the role of the chloroplast genome of M. micrantha during its recent invasion of southern China.

Understanding invasion history and predicting invasive niches using genetic sequencing technology in Australia: case studies from Cucurbitaceae and Boraginaceae

Part of the challenge in dealing with invasive plant species is that they seldom represent a uniform, static entity. Often, an accurate understanding of the history of plant introduction and knowledge of the real levels of genetic diversity present in species and populations of importance is lacking. Currently, the role of genetic diversity in promoting the successful establishment of invasive plants is not well defined. Genetic profiling of invasive plants should enhance our understanding of the dynamics of colonization in the invaded range. Recent advances in DNA sequencing technology have greatly facilitated the rapid and complete assessment of plant population genetics. Here, we apply our current understanding of the genetics and ecophysiology of plant invasions to recent work on Australian plant invaders from the Cucurbitaceae and Boraginaceae. The Cucurbitaceae study showed that both prickly paddy melon (Cucumis myriocarpus) and camel melon (Citrullus lanatus) were represented by only a single genotype in Australia, implying that each was probably introduced as a single introduction event. In contrast, a third invasive melon, Citrullus colocynthis, possessed a moderate level of genetic diversity in Australia and was potentially introduced to the continent at least twice. The Boraginaceae study demonstrated the value of comparing two similar congeneric species; one, Echium plantagineum, is highly invasive and genetically diverse, whereas the other, Echium vulgare, exhibits less genetic diversity and occupies a more limited ecological niche. Sequence analysis provided precise identification of invasive plant species, as well as information on genetic diversity and phylogeographic history. Improved sequencing technologies will continue to allow greater resolution of genetic relationships among invasive plant populations, thereby potentially improving our ability to predict the impact of these relationships upon future spread and better manage invaders possessing potentially diverse biotypes and exhibiting diverse breeding systems, life histories and invasion histories.

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.

Local adaptation to temperature and precipitation in naturally fragmented populations of Cephalotaxus oliveri, an endangered conifer endemic to China

Cephalotaxus oliveri is an endangered tertiary relict conifer endemic to China. The species survives in a wide range from west to east with heterogeneous climatic conditions. Precipitation and temperature are main restrictive factors for distribution of C. oliveri. In order to comprehend the mechanism of adaptive evolution to climate variation, we employed ISSR markers to detect adaptive evolution loci, to identify the association between variation in temperature and precipitation and adaptive loci, and to investigate the genetic structure for 22 C. oliveri natural populations. In total, 14 outlier loci were identified, of which five were associated with temperature and precipitation. Among outlier loci, linkage disequilibrium (LD) was high (42.86%), which also provided strong evidence for selection. In addition, C. oliveri possessed high genetic variation (93.31%) and population differentiation, which may provide raw material to evolution and accelerate local adaptation, respectively. Ecological niche modeling showed that global warming will cause a shift for populations of C. oliveri from south to north with a shrinkage of southern areas. Our results contribute to understand the potential response of conifers to climatic changes, and provide new insights for conifer resource management and conservation strategies.

Population Genetic Structure and Reproductive Strategy of the Introduced Grass Centotheca lappacea in Tropical Land-Use Systems in Sumatra

Intensive transformation of lowland rainforest into oil palm and rubber monocultures is the most common land-use practice in Sumatra (Indonesia), accompanied by invasion of weeds. In the Jambi province, Centotheca lappacea is one of the most abundant alien grass species in plantations and in jungle rubber (an extensively used agroforest), but largely missing in natural rainforests. Here, we investigated putative genetic differentiation and signatures for adaptation in the introduced area. We studied reproductive mode and ploidy level as putative factors for invasiveness of the species. We sampled 19 populations in oil palm and rubber monocultures and in jungle rubber in two regions (Bukit Duabelas and Harapan). Amplified fragment length polymorphisms (AFLP) revealed a high diversity of individual genotypes and only a weak differentiation among populations (FST = 0.173) and between the two regions (FST = 0.065). There was no significant genetic differentiation between the three land-use systems. The metapopulation of C. lappacea consists of five genetic partitions with high levels of admixture; all partitions appeared in both regions, but with different proportions. Within the Bukit Duabelas region we observed significant isolation-by-distance. Nine AFLP loci (5.3% of all loci) were under natural diversifying selection. All studied populations of C. lappacea were diploid, outcrossing and self-incompatible, without any hints of apomixis. The estimated residence time of c. 100 years coincides with the onset of rubber and oil palm planting in Sumatra. In the colonization process, the species is already in a phase of establishment, which may be enhanced by efficient selection acting on a highly diverse gene pool. In the land-use systems, seed dispersal might be enhanced by adhesive spikelets. At present, the abundance of established populations in intensively managed land-use systems might provide opportunities for rapid dispersal of C. lappacea across rural landscapes in Sumatra, while the invasion potential in rainforest ecosystems appears to be moderate as long as they remain undisturbed.

Small-Scale Habitat-Specific Variation and Adaptive Divergence of Photosynthetic Pigments in Different Alkali Soils in Reed Identified by Common Garden and Genetic Tests

Flexibility of photosynthetic pigment traits is an important adaptive mechanism through which plants can increase mean fitness in a variable environment. Unlike morphological traits in plants, photosythesis has been shown to exhibit phenotypic plasticity, responding rapidly to environmental conditions. Meanwhile, local adaptation at small scales is considered to be rare. Thus, detecting the small-scale adaptive divergence of photosynthetic pigments presents a challenge. Leaf concentrations of photosynthetic pigments under stressful conditions may be reduced or maintained. Concentrations of some pigments and/or ratio of Chlorophyll a (Chla) to Chlorophyll b (Chlb) do not change markedly in some species, such as the common reed, Phragmites australis, a cosmopolitan grass and common invader. Little is known about photosynthetic responses of this plant to varying levels of alkali salt. Few studies have attempted to account for the relationship between pigment accumulation and leaf position in wild plant populations in grasslands. In this study, photosynthetic pigment concentrations and the total Chl(a+b)/Car ratio decreased as the growing season progressed and were shown to be significantly lower in the habitat with a higher soil pH value and less moisture when compared between habitats. The Chla/Chlb ratio did not differ significantly between habitats, although it increased significantly over time. Leaves in the middle position may be functionally important in the response to soil conditions because only pigment concentrations and the Chl(a+b)/Car ratio of those leaves varied between habitats significantly. The outlier loci, used to evaluate molecular signatures of selection, were detected by Arlequin, Bayescan, and Bayenv analyses. In the simulated habitats of common garden, the local genotypes had higher values of Chla, Chlb, Chl(a+b), Car in their home habitat than did genotypes originating from the other habitat. Q_ST-F_ST comparisons provided evidence of divergent selection. It appears likely that soil moisture, pH and electric conductivity drove local adaptation. Combined approaches that utilize information on phenotypes from field and common garden experiments, genome-wide markers, and environmental data will be the most informative for understanding the adaptive nature of the intraspecific divergence.

Ecotypes of an ecologically dominant prairie grass (Andropogon gerardii) exhibit genetic divergence across the U.S. Midwest grasslands' environmental gradient

Big bluestem (Andropogon gerardii) is an ecologically dominant grass with wide distribution across the environmental gradient of U.S. Midwest grasslands. This system offers an ideal natural laboratory to study population divergence and adaptation in spatially varying climates. Objectives were to: (i) characterize neutral genetic diversity and structure within and among three regional ecotypes derived from 11 prairies across the U.S. Midwest environmental gradient, (ii) distinguish between the relative roles of isolation by distance (IBD) vs. isolation by environment (IBE) on ecotype divergence, (iii) identify outlier loci under selection and (iv) assess the association between outlier loci and climate. Using two primer sets, we genotyped 378 plants at 384 polymorphic AFLP loci across regional ecotypes from central and eastern Kansas and Illinois. Neighbour-joining tree and PCoA revealed strong genetic differentiation between Kansas and Illinois ecotypes, which was better explained by IBE than IBD. We found high genetic variability within prairies (80%) and even fragmented Illinois prairies, surprisingly, contained high within-prairie genetic diversity (92%). Using Bayenv2, 14 top-ranked outlier loci among ecotypes were associated with temperature and precipitation variables. Six of seven BayeScanFST outliers were in common with Bayenv2 outliers. High genetic diversity may enable big bluestem populations to better withstand changing climates; however, population divergence supports the use of local ecotypes in grassland restoration. Knowledge of genetic variation in this ecological dominant and other grassland species will be critical to understanding grassland response and restoration challenges in the face of a changing climate.

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.

Signatures of selection in the Iberian honey bee (Apis mellifera iberiensis) revealed by a genome scan analysis of single nucleotide polymorphisms

Understanding the genetic mechanisms of adaptive population divergence is one of the most fundamental endeavours in evolutionary biology and is becoming increasingly important as it will allow predictions about how organisms will respond to global environmental crisis. This is particularly important for the honey bee, a species of unquestionable ecological and economical importance that has been exposed to increasing human-mediated selection pressures. Here, we conducted a single nucleotide polymorphism (SNP)-based genome scan in honey bees collected across an environmental gradient in Iberia and used four FST -based outlier tests to identify genomic regions exhibiting signatures of selection. Additionally, we analysed associations between genetic and environmental data for the identification of factors that might be correlated or act as selective pressures. With these approaches, 4.4% (17 of 383) of outlier loci were cross-validated by four FST -based methods, and 8.9% (34 of 383) were cross-validated by at least three methods. Of the 34 outliers, 15 were found to be strongly associated with one or more environmental variables. Further support for selection, provided by functional genomic information, was particularly compelling for SNP outliers mapped to different genes putatively involved in the same function such as vision, xenobiotic detoxification and innate immune response. This study enabled a more rigorous consideration of selection as the underlying cause of diversity patterns in Iberian honey bees, representing an important first step towards the identification of polymorphisms implicated in local adaptation and possibly in response to recent human-mediated environmental changes.

Differential tolerance to copper, but no evidence of population-level genetic differences in a widely-dispersing native barnacle

Despite many estuaries having high levels of metal pollution, species are found to persist in these stressful environments. Populations of estuarine invertebrates exposed to toxic concentrations of such metals may be under selection. However, in species with a wide-dispersal potential, any short-term results of localized selection may be counteracted by external recruitment from populations not under selection. The barnacle Amphibalanus variegatus is found in nearshore coastal environments as well as sheltered embayments and estuaries, including metal-impacted estuaries, from New South Wales, Australia to Western Australia. The fertilised eggs of A. variegatus are brooded internally and released as larvae (nauplii), which remain in the water-column for ~14 days before settling. Hence the species has a considerable dispersal capacity. The purpose of this study was to examine whether populations of A. variegatus from metal-impacted sites, displayed a greater tolerance to a toxicant (copper) than reference populations. Adult barnacles where collected from twenty sites within two metal-impacted and fourteen sites within two reference estuaries. Within 24 h, adults were induced to spawn and the offspring were exposed to copper in a laboratory assay. Larvae collected from the metal-impacted estuaries demonstrated a greater tolerance to copper compared to those from reference sites. To determine if selection/localised in the metal impacted sites was occurring, the genetic structure of populations at three sites was examined using an AFLP methodology. No evidence of unique population identity and or selection (outlier loci) was detected suggesting that: (1) the tolerance displayed in the assay was derived from acclimation during development; and/or (2) that the populations are open preventing the fixation of any unique alleles.