Are habitat fragmentation, local adaptation and isolation-by-distance driving population divergence in wild riceOryza rufipogon?

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.

Genetic structure of wild rice Zizania latifolia in an expansive heterogeneous landscape along a latitudinal gradient

Global aquatic habitats are undergoing rapid degradation and fragmentation as a result of climate change and changes in land use. Understanding the genetic variability and adaptive potential of aquatic plant species is thus important for conservation purposes. In this study, we investigated the genetic diversity and structure of the extant natural populations of Zizania latifolia from five river basins in China based on 46 microsatellite markers. We tested isolation by environment (IBE), isolation by resistance (IBR), and isolation by distance (IBD) patterns using a reciprocal causal model (RCM). Furthermore, we elucidated the impact of the environment on Z. latifolia genetic diversity using generalized linear models (GLMs) and spatially explicit mixed models. Low genetic diversity (HE = 0.125–0.433) and high genetic differentiation (FST = 0.641, Øpt = 0.654) were found. Higher historical gene flow (MH = 0.212–2.354) than contemporary gene flow (MC = 0.0112–0.0247) and significant bottlenecks in almost all populations were identified, highlighting the negative impact of wetland fragmentation. The IBE model was exclusively supported for all populations and in three river basins. The IBD and IBR models were supported in one river basin each. The maximum temperature of the warmest month and precipitation seasonality were the plausible environmental parameters responsible for the observed pattern of genetic diversity. Local adaptation signatures were found, with nine loci identified as outliers, four of which were gene-linked and associated with environmental variables. Based on these findings, IBE is more important than IBD and IBR in shaping the genetic structure of Z. latifolia.

High-throughput sequencing-based microsatellite genotyping for polyploids to resolve allele dosage uncertainty and improve analyses of genetic diversity, structure and differentiation: A case study of the hexaploid Camellia oleifera

Conventional microsatellite (simple sequence repeat, SSR) genotyping methods cannot accurately identify polyploid genotypes leading to allele dosage uncertainty, introducing biases in population genetic analysis. Here, a new SSR genotyping method was developed to directly infer accurate polyploid genotypes. The frequency distribution of SSR sequences was obtained based on deep-coverage high-throughput sequencing data. Corrections were performed accounting for the "stutter peak" and amplification efficiency of SSR sequences. Perl scripts and an online SSR genotyping tool "SSRSeq" were provided to process the sequencing data and output genotypes with corrected allele dosages. Hexaploid Camellia oleifera is the dominant woody oilseed crop in China. Understanding the geographical pattern of genetic variation in wild C. oleifera is essential for the conservation and utilization of genetic resources. Six wild C. oleifera populations were sampled across geographical ranges in subtropical evergreen broadleaf forests of China. Using 35 SSR markers, the high-throughput sequencing-based SSRSeq method was applied to obtain accurate hexaploid genotypes of wild C. oleifera. The results demonstrated that the new method could resolve allele dosage uncertainty and considerably improve genetic diversity, structure and differentiation analyses for polyploids. The genetic variation patterns of wild C. oleifera across geographical ranges agree with the "central-marginal hypothesis", stating that genetic diversity is high in the central population and declines from the central to the peripheral populations, and genetic differentiation increases from the centre to the periphery. This method and findings can facilitate the utilization of wild C. oleifera genetic resources for the breeding of cultivated C. oleifera.

Genome-Wide Analysis of CCT Transcript Factors to Identify Genes Contributing to Photoperiodic Flowering in Oryza rufipogon

Photoperiod sensitivity is a dominant determinant for the phase transition in cereal crops. CCT (CONSTANS, CO-like, and TOC1) transcription factors (TFs) are involved in many physiological functions including the regulation of the photoperiodic flowering. However, the functional roles of CCT TFs have not been elucidated in the wild progenitors of crops. In this study, we identified 41 CCT TFs, including 19 CMF, 17 COL, and five PRR TFs in Oryza rufipogon, the presumed wild ancestor of Asian cultivated rice. There are thirty-eight orthologous CCT genes in Oryza sativa, of which ten pairs of duplicated CCT TFs are shared with O. rufipogon. We investigated daily expression patterns, showing that 36 OrCCT genes exhibited circadian rhythmic expression. A total of thirteen OrCCT genes were identified as putative flowering suppressors in O. rufipogon based on rhythmic and developmental expression patterns and transgenic phenotypes. We propose that OrCCT08, OrCCT24, and OrCCT26 are the strong functional alleles of rice DTH2, Ghd7, and OsPRR37, respectively. The SD treatment at 80 DAG stimulated flowering of the LD-grown O. rufipogon plants. Our results further showed that the nine OrCCT genes were significantly downregulated under the treatment. Our findings would provide valuable information for the construction of photoperiodic flowering regulatory network and functional characterization of the CCT TFs in both O. rufipogon and O. sativa.

Conservation recommendations for Oryza rufipogon Griff. in China based on genetic diversity analysis

Over the past 30 years, human disturbance and habitat fragmentation have severely endangered the survival of common wild rice (Oryza rufipogon Griff.) in China. A better understanding of the genetic structure of O. rufipogon populations will therefore be useful for the development of conservation strategies. We examined the diversity and genetic structure of natural O. rufipogon populations at the national, provincial, and local levels using simple sequence repeat (SSR) markers. Twenty representative populations from sites across China showed high levels of genetic variability, and approximately 44% of the total genetic variation was among populations. At the local level, we studied fourteen populations in Guangxi Province and four populations in Jiangxi Province. Populations from similar ecosystems showed less genetic differentiation, and local environmental conditions rather than geographic distance appeared to have influenced gene flow during population genetic evolution. We identified a triangular area, including northern Hainan, southern Guangdong, and southwestern Guangxi, as the genetic diversity center of O. rufipogon in China, and we proposed that this area should be given priority during the development of ex situ and in situ conservation strategies. Populations from less common ecosystem types should also be given priority for in situ conservation.

Population genetic structure and gene flow of rare and endangered Tetraena mongolica Maxim. revealed by reduced representation sequencing

BACKGROUND

Studying population genetic structure and gene flow of plant populations and their influencing factors is of particular significance in the field of conservation biology, especially important for species such as rare and endangered plants. Tetraena mongolica Maxim. (TM), belongs to Zygophyllaceae family, a rare and endangered plant with narrow distribution. However, for the last decade, due to excessive logging, urban expansion, industrial and tourism development, habitat fragmentation and loss of natural habitats have become major threats to the population of endangered plants.

RESULTS

In this study, genetic diversity, population genetic structure and gene flow of TM populations were evaluated by reduced representation sequencing technology, and a total of more than 133.45 GB high-quality clean reads and 38,097 high-quality SNPs were generated. Analysis based on multiple methods, we found that the existing TM populations have moderate levels of genetic diversity, and very low genetic differentiation as well as high levels of gene flow between populations. Population structure and principal coordinates analysis showed that 8 TM populations can be divided into two groups. The Mantel test detected no significant correlation between geographical distances and genetic distance for the whole sampling. Moreover, the migration model indicated that the gene flow is more of a north to south migration pattern in history.

CONCLUSIONS

This study demonstrates that the present genetic structure is mainly due to habitat fragmentation caused by urban sprawl, industrial development and coal mining. Our recommendation with respect to conservation management is that, all 8 populations should be preserved as a whole population, rather than just those in the core area of TM nature reserve. In particular, the populations near the edge of TM distribution in cities and industrial areas deserve our special protection.

Inferring the Origin of Cultivated Zizania latifolia, an Aquatic Vegetable of a Plant-Fungus Complex in the Yangtze River Basin

Crop domestication is one of the essential topics in evolutionary biology. Cultivated Zizania latifolia, domesticated as the special form of a plant-fungus (the host Zizania latifolia and the endophyte Ustilago esculenta) complex, is a popular aquatic vegetable endemic in East Asia. The rapid domestication of cultivated Z. latifolia can be traced in the historical literature but still need more evidence. This study focused on deciphering the genetic relationship between wild and cultivated Z. latifolia, as well as the corresponding parasitic U. esculenta. Twelve microsatellites markers were used to study the genetic variations of 32 wild populations and 135 landraces of Z. latifolia. Model simulations based on approximate Bayesian computation (ABC) were then performed to hierarchically infer the population history. We also analyzed the ITS sequences of the smut fungus U. esculenta to reveal its genetic structure. Our results indicated a significant genetic divergence between cultivated Z. latifolia and its wild ancestors. The wild Z. latifolia populations showed significant hierarchical genetic subdivisions, which may be attributed to the joint effect of isolation by distance and hydrological unconnectedness between watersheds. Cultivated Z. latifolia was supposedly domesticated once in the low reaches of the Yangtze River. The genetic structure of U. esculenta also indicated a single domestication event, and the genetic variations in this fungus might be associated with the diversification of cultivars. These findings provided molecular evidence in accordance with the historical literature that addressed the domestication of cultivated Z. latifolia involved adaptive evolution driven by artificial selection in both the plant and fungus.

Seed characteristic variations and genetic structure of wild Zizania latifolia along a latitudinal gradient in China: implications for neo-domestication as a grain crop

Crop wild relatives are not only important genetic resources for crop improvement, but also domestication candidates for selecting new crops. As a close relative of American wild rice Zizania palustris, Z. latifolia is a perennial aquatic grass widely distributed in China. Although Z. latifolia has been domesticated and cultivated as an aquatic vegetable for >1000 years, a neo-domestication for grain production needs to be soundly evaluated. In this study, we investigated the seed characteristic variations and genetic structure of 15 Z. latifolia wild populations along a latitudinal gradient in China. Our results showed that Z. latifolia tended to produce relatively larger seeds with lower moisture content and lower investments in seed pericarp at lower latitudes. The width, size, shape, seed-pericarp ratio and relative water content of seeds were significantly associated with climatic variables. The seeds of Z. latifolia showed a relatively low germination percentage and strong dormancy, which might hinder the neo-domestication. In addition, high genetic differentiation had been found among Z. latifolia populations, which could be attributed to isolation by distance. This study offered preliminary information for the utilization and conservation of wild Z. latifolia. It suggested that the wild populations in the middle and lower reaches of the Yangtze River could be good candidates for grain crop domestication due to appropriate seed traits and high genetic diversity. The neo-domestication of wild Z. latifolia requires further researches on the genetic mechanism of the Domestication Syndrome and more works on artificial breeding.

Genetic differentiation, local adaptation and phenotypic plasticity in fragmented populations of a rare forest herb

Background

Due to habitat loss and fragmentation, numerous forest species are subject to severe population decline. Investigating variation in genetic diversity, phenotypic plasticity and local adaptation should be a prerequisite for implementing conservation actions. This study aimed to explore these aspects in ten fragmented populations of Physospermum cornubiense in view of translocation measures across its Italian range.

Methods

For each population we collected environmental data on landscape (habitat size, quality and fragmentation) and local conditions (slope, presence of alien species, incidence of the herbivorous insect Metcalfa pruinosa and soil parameters). We measured vegetative and reproductive traits in the field and analysed the genetic population structure using ISSR markers (STRUCTURE and AMOVA). We then estimated the neutral (F_ST) and quantitative (P_ST) genetic differentiation of populations.

Results

The populations exhibited moderate phenotypic variation. Population size (range: 16–655 individuals), number of flowering adults (range: 3–420 individuals) and inflorescence size (range: 5.0–8.4 cm) were positively related to Mg soil content. Populations’ gene diversity was moderate (Nei-H = 0.071–0.1316); STRUCTURE analysis identified five different clusters and three main geographic groups: upper, lower, and Apennine/Western Po plain. Fragmentation did not have an influence on the local adaptation of populations, which for all measured traits showed P_ST < F_ST, indicating convergent selection.

Discussion

The variation of phenotypic traits across sites was attributed to plastic response rather than local adaptation. Plant translocation from suitable source populations to endangered ones should particularly take into account provenance according to identified genetic clusters and specific soil factors.

The Tanggula Mountains enhance population divergence in Carex moorcroftii: a dominant sedge on the Qinghai-Tibetan Plateau

High-altitude mountains are often geographic barriers to gene flow and play important roles in shaping population divergence. The central Qinghai-Tibetan Plateau (QTP) stands the location of the Tanggula Mountains (TM). We use the TM as a case, using Carex moorcroftii, a dominant species on the QTP to test the effects of geographic barriers on plant population divergence. We sampled 18 C. moorcroftii populations along a north-south transect crossing the TM to investigate the correlations of genetic variation and morphological traits with climate variables. The results showed this species holds high genetic diversity (H_e = 0.58) and the surveyed populations can be genetically clustered into two groups: populations from the north face of TM, and the other from the south. Gene flow between populations within groups is higher than those between groups. The traits, number and mass of seeds, mass of root and infructescence significantly varied among populations. Mantel-tests detected a weak but significantly positive correlation between genetic and geographic (R² = 0.107, p = 0.032) and climatic distance (R² = 0.162, p = 0.005), indicating both isolation by distance and isolation by environment. These findings together suggest high-altitude mountains of TM interrupt habitat continuity, result in distinct climatic conditions on both sides, increasing population divergence of plant species.

Genetically engineered rice endogenous 5-enolpyruvoylshikimate-3-phosphate synthase (epsps) transgene alters phenology and fitness of crop-wild hybrid offspring

Genetically engineered (GE) rice endogenous epsps (5-enolpyruvoylshikimate-3-phosphate synthase) gene overexpressing EPSPS can increase glyphosate herbicide-resistance of cultivated rice. This type of epsps transgene can enhance the fecundity of rice crop-weed hybrid offspring in the absence of glyphosate, stimulating great concerns over undesired environmental impacts of transgene flow to populations of wild relatives. Here, we report the substantial alteration of phenology and fitness traits in F1-F3 crop-wild hybrid descendants derived from crosses between an epsps GE rice line and two endangered wild rice (Oryza rufipogon) populations, based on the common-garden field experiments. Under the glyphosate-free condition, transgenic hybrid lineages showed significantly earlier tillering and flowering, as well as increased fecundity and overwintering survival/regeneration abilities. In addition, a negative correlation was observed between the contents of endogenous EPSPS of wild, weedy, and cultivated rice parents and fitness differences caused by the incorporation of the epsps transgene. Namely, a lower level of endogenous EPSPS in the transgene-recipient populations displayed a more pronounced enhancement in fitness. The altered phenology and enhanced fitness of crop-wild hybrid offspring by the epsps transgene may cause unwanted environmental consequences when this type of glyphosate-resistance transgene introgressed into wild rice populations through gene flow.

Ecological genomics of local adaptation in Cornus florida L. by genotyping by sequencing

Discovering local adaptation, its genetic underpinnings, and environmental drivers is important for conserving forest species. Ecological genomic approaches coupled with next-generation sequencing are useful means to detect local adaptation and uncover its underlying genetic basis in nonmodel species. We report results from a study on flowering dogwood trees (Cornus florida L.) using genotyping by sequencing (GBS). This species is ecologically important to eastern US forests but is severely threatened by fungal diseases. We analyzed subpopulations in divergent ecological habitats within North Carolina to uncover loci under local selection and associated with environmental-functional traits or disease infection. At this scale, we tested the effect of incorporating additional sequencing before scaling for a broader examination of the entire range. To test for biases of GBS, we sequenced two similarly sampled libraries independently from six populations of three ecological habitats. We obtained environmental-functional traits for each subpopulation to identify associations with genotypes via latent factor mixed modeling (LFMM) and gradient forests analysis. To test whether heterogeneity of abiotic pressures resulted in genetic differentiation indicative of local adaptation, we evaluated Fst per locus while accounting for genetic differentiation between coastal subpopulations and Piedmont-Mountain subpopulations. Of the 54 candidate loci with sufficient evidence of being under selection among both libraries, 28-39 were Arlequin-BayeScan Fst outliers. For LFMM, 45 candidates were associated with climate (of 54), 30 were associated with soil properties, and four were associated with plant health. Reanalysis of combined libraries showed that 42 candidate loci still showed evidence of being under selection. We conclude environment-driven selection on specific loci has resulted in local adaptation in response to potassium deficiencies, temperature, precipitation, and (to a marginal extent) disease. High allele turnover along ecological gradients further supports the adaptive significance of loci speculated to be under selection.

Population Genetic Structure of Glycyrrhiza inflata B. (Fabaceae) Is Shaped by Habitat Fragmentation, Water Resources and Biological Characteristics

BACKGROUND

Habitat fragmentation, water resources and biological characteristics are important factors that shape the genetic structure and geographical distribution of desert plants. Analysis of the relationships between these factors and population genetic variation should help to determine the evolutionary potential and conservation strategies for genetic resources for desert plant populations. As a traditional Chinese herb, Glycyrrhiza inflata B. (Fabaceae) is restricted to the fragmented desert habitat in China and has undergone a dramatic decline due to long-term over-excavation. Determining the genetic structure of the G. inflata population and identifying a core collection could help with the development of strategies to conserve this species.

RESULTS

We investigated the genetic variation of 25 G. inflata populations based on microsatellite markers. A high level of population genetic divergence (FST = 0.257), population bottlenecks, reduced gene flow and moderate genetic variation (HE = 0.383) were detected. The genetic distances between the populations significantly correlated with the geographical distances, and this suggests that habitat fragmentation has driven a special genetic structure of G. inflata in China through isolation by distance. STRUCTURE analysis showed that G. inflata populations were structured into three clusters and that the populations belonged to multiple water systems, which suggests that water resources were related to the genetic structure of G. inflata. In addition, the biological characteristics of the perennial species G. inflata, such as its long-lived seeds, asexual reproduction, and oasis ecology, may be related to its resistance to habitat fragmentation. A core collection of G. inflata, that included 57 accessions was further identified, which captured the main allelic diversity of G. inflata.

CONCLUSIONS

Recent habitat fragmentation has accelerated genetic divergence. The population genetic structure of G. inflata has been shaped by habitat fragmentation, water resources and biological characteristics. This genetic information and core collection will facilitate the conservation of wild germplasm and breeding of this Chinese medicinal plant.

Great influence of geographic isolation on the genetic differentiation of Myriophyllum spicatum under a steep environmental gradient

Understanding how natural processes affect population genetic structures is an important issue in evolutionary biology. One effective method is to assess the relative importance of environmental and geographical factors in the genetic structure of populations. In this study, we examined the spatial genetic variation of thirteen Myriophyllum spicatum populations from the Qinghai-Tibetan Plateau (QTP) and adjacent highlands (Yunnan-Guizhou Plateau, YGP) by using microsatellite loci and environmental and geographical factors. Bioclim layers, hydrological properties and elevation were considered as environmental variables and reduced by principal component analysis. The genetic isolation by geographic distance (IBD) was tested by Mantel tests and the relative importance of environmental variables on population genetic differentiation was determined by a partial Mantel test and multiple matrix regression with randomization (MMRR). Two genetic clusters corresponding to the QTP and YGP were identified. Both tests and MMRR revealed a significant and strong correlation between genetic divergence and geographic isolation under the influence of environmental heterogeneity at the overall and finer spatial scales. Our findings suggested the dominant role of geography on the evolution of M. spicatum under a steep environmental gradient in the alpine landscape as a result of dispersal limitation and genetic drift.