Landscape Population Genomics of Forsythia (Forsythia suspensa) Reveal That Ecological Habitats Determine the Adaptive Evolution of Species

Landscape genomics reveals adaptive genetic differentiation driven by multiple environmental variables in naked barley on the Qinghai-Tibetan Plateau

Understanding the local adaptation of crops has long been a concern of evolutionary biologists and molecular ecologists. Identifying the adaptive genetic variability in the genome is crucial not only to provide insights into the genetic mechanism of local adaptation but also to explore the adaptation potential of crops. This study aimed to identify the climatic drivers of naked barley landraces and putative adaptive loci driving local adaptation on the Qinghai-Tibetan Plateau (QTP). To this end, a total of 157 diverse naked barley accessions were genotyped using the genotyping-by-sequencing approach, which yielded 3123 high-quality SNPs for population structure analysis and partial redundancy analysis, and 37,636 SNPs for outlier analysis. The population structure analysis indicated that naked barley landraces could be divided into four groups. We found that the genomic diversity of naked barley landraces could be partly traced back to the geographical and environmental diversity of the landscape. In total, 136 signatures associated with temperature, precipitation, and ultraviolet radiation were identified, of which 13 had pleiotropic effects. We mapped 447 genes, including a known gene HvSs1. Some genes involved in cold stress and regulation of flowering time were detected near eight signatures. Taken together, these results highlight the existence of putative adaptive loci in naked barley on QTP and thus improve our current understanding of the genetic basis of local adaptation.

Genome-wide investigation of the WRKY transcription factor gene family in weeping forsythia: expression profile and cold and drought stress responses

Weeping forsythia is a wide-spread shrub in China with important ornamental, medicinal and ecological values. It is widely distributed in China's warm temperate zone. In plants, WRKY transcription factors play important regulatory roles in seed germination, flower development, fruit ripening and coloring, and biotic and abiotic stress response. To date, WRKY transcription factors have not been systematically studied in weeping forsythia. In this study, we identified 79 WRKY genes in weeping forsythia and classified them according to their naming rules in Arabidopsis thaliana. Phylogenetic tree analysis showed that, except for IIe subfamily, whose clustering was inconsistent with A. thaliana clustering, other subfamily clustering groups were consistent. Cis-element analysis showed that WRKY genes related to pathogen resistance in weeping forsythia might be related to methyl jasmonate and salicylic acid-mediated signaling pathways. Combining cis-element and expression pattern analyses of WRKY genes showed that more than half of WRKY genes were involved in light-dependent development and morphogenesis in different tissues. The gene expression results showed that 13 WRKY genes were involved in drought response, most of which might be related to the abscisic acid signaling pathway, and a few of which might be regulated by MYB transcription factors. The gene expression results under cold stress showed that 17 WRKY genes were involved in low temperature response, and 9 of them had low temperature responsiveness cis-elements. Our study of WRKY family in weeping forsythia provided useful resources for molecular breeding and important clues for their functional verification.

EST-SSR Primer Development and Genetic Structure Analysis of Psathyrostachys juncea Nevski

Psathyrostachys juncea is a perennial forage grass which plays an important role in soil and water conservation and ecological maintenance in cold and dry areas of temperate regions. In P. juncea, a variety of biotic and abiotic stress related genes have been used in crop improvement, indicating its agronomic, economic, forage, and breeding value. To date, there have been few studies on the genetic structure of P. juncea. Here, the genetic diversity and population structure of P. juncea were analyzed by EST-SSR molecular markers to evaluate the genetic differentiation related to tillering traits in P. juncea germplasm resources. The results showed that 400 simple sequence repeat (SSR) loci were detected in 2,020 differentially expressed tillering related genes. A total of 344 scored bands were amplified using 103 primer pairs, out of which 308 (89.53%) were polymorphic. The Nei's gene diversity of 480 individuals was between 0.092 and 0.449, and the genetic similarity coefficient was between 0.5008 and 0.9111, with an average of 0.6618. Analysis of molecular variance analysis showed that 93% of the variance was due to differences within the population, and the remaining 7% was due to differences among populations. Psathyrostachys juncea materials were clustered into five groups based on population genetic structure, principal coordinate analysis and unweighted pair-group method with arithmetic means (UPGMA) analysis. The results were similar between clustering methods, but a few individual plants were distributed differently by the three models. The clustering results, gene diversity and genetic similarity coefficients showed that the overall genetic relationship of P. juncea individuals was relatively close. A Mantel test, UPGMA and structural analysis also showed a significant correlation between genetic relationship and geographical distribution. These results provide references for future breeding programs, genetic improvement and core germplasm collection of P. juncea.

Large- and small-scale geographic structures affecting genetic patterns across populations of an Alpine butterfly

Understanding factors influencing patterns of genetic diversity and the population genetic structure of species is of particular importance in the current era of global climate change and habitat loss. These factors include the evolutionary history of a species as well as heterogeneity in the environment it occupies, which in turn can change across time. Most studies investigating spatio-temporal genetic patterns have focused on patterns across wide geographic areas rather than local variation, but the latter can nevertheless be important particularly in topographically complex areas. Here, we consider these issues in the Sooty Copper butterfly (Lycaena tityrus) from the European Alps, using genome-wide SNPs identified through RADseq. We found strong genetic differentiation within the Alps with four genetic clusters, indicating western, central, and eastern refuges, and a strong reduction of genetic diversity from west to east. This reduction in diversity may suggest that the southwestern refuge was the largest one in comparison to other refuges. Also, the high genetic diversity in the west may result from (a) admixture of different western refuges, (b) more recent demographic changes, or (c) introgression of lowland L. tityrus populations. At small spatial scales, populations were structured by several landscape features and especially by high mountain ridges and large river valleys. We detected 36 outlier loci likely under altitudinal selection, including several loci related to membranes and cellular processes. We suggest that efforts to preserve alpine L. tityrus should focus on the genetically diverse populations in the western Alps, and that the dolomite populations should be treated as genetically distinct management units, since they appear to be currently more threatened than others. This study demonstrates the usefulness of SNP-based approaches for understanding patterns of genetic diversity, gene flow, and selection in a region that is expected to be particularly vulnerable to climate change.

Non-linear genetic diversity and notable population differentiation caused by low gene flow of bermudagrass [Cynodon dactylon (L.) Pers.] along longitude gradients

Background

Environmental variation related to ecological habitat is the main driver of plant adaptive divergence. Longitude plays an important role in the formation of plant population structure, indicating that environmental differentiation can significantly shape population structure.

Methods

Genetic diversity and population genetic structure were estimated using 105 expressed sequence tag-derived simple sequence repeat (EST-SSR) loci. A total of 249 C. dactylon (L.) Pers. (common bermudagrass) individuals were sampled from 13 geographic sites along the longitude (105°57′34″–119°27′06″E).

Results

There was no obvious linear trend of intra-population genetic diversity along longitude and the intra-population genetic diversity was not related to climate in this study. Low gene flow (Nm = 0.7701) meant a rich genetic differentiation among populations of C. dactylon along longitude gradients. Significantly positive Mantel correlation (r = 0.438, P = 0.001) was found between genetic distance and geographical interval while no significant partial Mantel correlation after controlling the effect of mean annual precipitation, which indicated geographic distance correlated with mean annual precipitation affect genetic distance. The genetic diversity of C. dactylon with higher ploidy level was higher than that with lower ploidy level and groups of individuals with higher ploidy level were separated further away by genetic distance from the lower ploidy levels. Understanding the different genetic bases of local adaptation comparatively between latitude and longitude is one of the core findings in the adaptive evolution of plants.

Transcriptomic responses to drought stress among natural populations provide insights into local adaptation of weeping forsythia

BACKGROUND

Understanding the genetic mechanisms of local adaptation is an important emerging topic in molecular ecology and evolutionary biology.

RESULTS

Here, we identify the physiological changes and differential expression of genes among different weeping forsythia populations under drought stress in common garden experiments. Physiological results showed that HBWZ might have higher drought tolerance among four populations. RNA-seq results showed that significant differential expression in the genes responding to the synthesis of flavonoids, aromatic substances, aromatic amino acids, oxidation-reduction process, and transmembrane transport occured among four populations. By further reanalysis of results of previous studies, sequence differentiation was found in the genes related to the synthesis of aromatic substances among different weeping forsythia populations.

CONCLUSIONS

Overall, our study supports the hypothesis that the dual differentiation in gene efficiency and expression increases among populations in response to heterogeneous environments and is an important evolutionary process of local adaptation. Here, we proposed a new working model of local adaptation of weeping forsythia populations under different intensities of drought stress, which provides new insights for understanding the genetic mechanisms of local adaptation for non-model species.

Genome sequencing and population genomics modeling provide insights into the local adaptation of weeping forsythia

Understanding the genetic basis underlying the local adaptation of nonmodel species is a fundamental goal in evolutionary biology. In this study, we explored the genetic mechanisms of the local adaptation of Forsythia suspensa using genome sequence and population genomics data obtained from specific-locus amplified fragment sequencing. We assembled a high-quality reference genome of weeping forsythia (Scaffold N50 = 7.3 Mb) using ultralong Nanopore reads. Then, genome-wide comparative analysis was performed for 15 natural populations of weeping forsythia across its current distribution range. Our results revealed that candidate genes associated with local adaptation are functionally correlated with solar radiation, temperature and water variables across heterogeneous environmental scenarios. In particular, solar radiation during the period of fruit development and seed drying after ripening, cold, and drought significantly contributed to the adaptive differentiation of F. suspensa. Natural selection exerted by environmental factors contributed substantially to the population genetic structure of F. suspensa. Our results supported the hypothesis that adaptive differentiation should be highly pronounced in the genes involved in signal crosstalk between different environmental variables. Our population genomics study of F. suspensa provides insights into the fundamental genetic mechanisms of the local adaptation of plant species to climatic gradients.

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.

Molecular Phylogeography Analysis Reveals Population Dynamics and Genetic Divergence of a Widespread Tree Pterocarya stenoptera in China

The geological events, past climatic fluctuations, and river systems played key roles in the spatial distribution, population dynamics, and genetic differentiation of species. In this work, we selected Pterocarya stenoptera, a widespread tree species in China, to test the roles of these factors. Four noncoding spacers, eight microsatellite (simple sequence repeat) markers, and species distribution modeling were used to examine the phylogeographical pattern of P. stenoptera. Based on chloroplast DNA data, populations of P. stenoptera were clearly clustered into three groups. The divergence time of these groups fell into the stage of the Qinghai–Tibet Movement, 1.7–2.6 Ma. For simple sequence repeat data, only one western marginal population YNYB could be separated from other populations, whereas other populations were mixed together. Our results indicated that the environmental heterogeneity resulting from the Qinghai–Tibet movement might be response for this genetic divergence. The climatic fluctuations in the Pleistocene did not cause the substantial range shift of P. stenoptera, while the fluctuations affected its population size. Moreover, we also confirmed the river systems did not act as channels or barrier of dispersal for P. stenoptera.

Landscape genetics reveals that adaptive genetic divergence in Pinus bungeana (Pinaceae) is driven by environmental variables relating to ecological habitats

Background

Understanding the genetic basis of local adaptation has long been the concern of biologists. Identifying these adaptive genetic variabilities is crucial not only to improve our knowledge of the genetic mechanism of local adaptation but also to explore the adaptation potential of species.

Results

Using 10 natural populations and 12 start codon targeted (SCoT) markers, a total of 430 unambiguous loci were yielded. The Bayesian analysis of population structure clearly demonstrated that the 10 populations of P. bungeana could be subdivided into three groups. Redundancy analysis showed that this genetic divergence was caused by divergence selection from environmental variables related to the ecological habitats of “avoidance of flooding” and “avoidance of high temperature and humidity.” LFMM results indicated that Bio1, Bio5, Bio8, Bio12, Bio14, and Bio16, which are related to the ecological habitat of P. bungeana, were correlated with the highest numbers of environment-associated loci (EAL).

Conclusions

The results of EAL characterization in P. bungeana clearly supported the hypothesis that environmental variations related to the ecological habitat of species are the key drivers of species adaptive divergence. Moreover, a method to calculate the species landscape adaptation index and quantify the adaptation potential of species was proposed and verified using ecological niche modeling. This model could estimate climatically suitable areas of species spatial distribution. Taking the results together, this study improves the current understanding on the genetic basis of local adaptation.

Electronic supplementary material

The online version of this article (10.1186/s12862-019-1489-x) contains supplementary material, which is available to authorized users.

Adaptive genetic differentiation in Pterocarya stenoptera (Juglandaceae) driven by multiple environmental variables were revealed by landscape genomics

BACKGROUND

The investigation of the genetic basis of local adaptation in non-model species is an interesting focus of evolutionary biologists and molecular ecologists. Identifying these adaptive genetic variabilities on the genome responsible can provide insight into the genetic mechanism of local adaptation.

RESULTS

We investigated the spatial distribution of genetic variation in 22 natural populations of Pterocarya stenoptera across its distribution area in China to provide insights into the complex interplay between multiple environmental variables and adaptive genetic differentiation. The Bayesian analysis of population structure showed that the 22 populations of P. stenoptera were subdivided into two groups. Redundancy analysis demonstrated that this genetic differentiation was caused by the divergent selection of environmental difference. A total of 44 outlier loci were mutually identified by Arlequin and BayeScan, 43 of which were environment-associated loci (EAL). The results of latent factor mixed model analysis showed that solar radiation in June (Sr6), minimum temperature of the coldest month (Bio6), temperature seasonality (Bio4), and water vapor pressure in January (Wvp1) were associated with the highest numbers of EAL. Sr6 was associated with the ecological habitat of "prefered light", and Bio6 and Wvp1 were associated with the ecological habitat of "warm and humid environment".

CONCLUSIONS

Our results provided empirical evidence that environmental variables related to the ecological habitats of species play key roles in driving adaptive differentiation of species genome.

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.