Conservation genetics of two co-dominant grass species in an endangered grassland ecosystem

Genetic and functional variation across regional and local scales is associated with climate in a foundational prairie grass

Global change forecasts in ecosystems require knowledge of within-species diversity, particularly of dominant species within communities. We assessed site-level diversity and capacity for adaptation in Bouteloua gracilis, the dominant species in the Central US shortgrass steppe biome. We quantified genetic diversity from 17 sites across regional scales, north to south from New Mexico to South Dakota, and local scales in northern Colorado. We also quantified phenotype and plasticity within and among sites and determined the extent to which phenotypic diversity in B. gracilis was correlated with climate. Genome sequencing indicated pronounced population structure at the regional scale, and local differences indicated that gene flow and/or dispersal may also be limited. Within a common environment, we found evidence of genetic divergence in biomass-related phenotypes, plasticity, and phenotypic variance, indicating functional divergence and different adaptive potential. Phenotypes were differentiated according to climate, chiefly median Palmer Hydrological Drought Index and other aridity metrics. Our results indicate conclusive differences in genetic variation, phenotype, and plasticity in this species and suggest a mechanism explaining variation in shortgrass steppe community responses to global change. This analysis of B. gracilis intraspecific diversity across spatial scales will improve conservation and management of the shortgrass steppe ecosystem in the future.

Single nucleotide polymorphism discovery via genotyping by sequencing to assess population genetic structure and recurrent polyploidization in Andropogon gerardii

PREMISE OF THE STUDY

Autopolyploidy, genome duplication within a single lineage, can result in multiple cytotypes within a species. Geographic distributions of cytotypes may reflect the evolutionary history of autopolyploid formation and subsequent population dynamics including stochastic (drift) and deterministic (differential selection among cytotypes) processes. Here, we used a population genomic approach to investigate whether autopolyploidy occurred once or multiple times in Andropogon gerardii, a widespread, North American grass with two predominant cytotypes.

METHODS

Genotyping by sequencing was used to identify single nucleotide polymorphisms (SNPs) in individuals collected from across the geographic range of A. gerardii. Two independent approaches to SNP calling were used: the reference-free UNEAK pipeline and a reference-guided approach based on the sequenced Sorghum bicolor genome. SNPs generated using these pipelines were analyzed independently with genetic distance and clustering.

KEY RESULTS

Analyses of the two SNP data sets showed very similar patterns of population-level clustering of A. gerardii individuals: a cluster of A. gerardii individuals from the southern Plains, a northern Plains cluster, and a western cluster. Groupings of individuals corresponded to geographic localities regardless of cytotype: 6x and 9x individuals from the same geographic area clustered together.

CONCLUSIONS

SNPs generated using reference-guided and reference-free pipelines in A. gerardii yielded unique subsets of genomic data. Both data sets suggest that the 9x cytotype in A. gerardii likely evolved multiple times from 6x progenitors across the range of the species. Genomic approaches like GBS and diverse bioinformatics pipelines used here facilitate evolutionary analyses of complex systems with multiple ploidy levels.

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 sorting of subordinate species in grassland modulated by intraspecific variation in dominant species

Genetic variation in a single species can have predictable and heritable effects on associated communities and ecosystem processes, however little is known about how genetic variation of a dominant species affects plant community assembly. We characterized the genetic structure of a dominant grass (Sorghastrum nutans) and two subordinate species (Chamaecrista fasciculata, Silphium integrifolium), during the third growing season in grassland communities established with genetically distinct (cultivated varieties or local ecotypes) seed sources of the dominant grasses. There were genetic differences between subordinate species growing in the cultivar versus local ecotype communities, indicating that intraspecific genetic variation in the dominant grasses affected the genetic composition of subordinate species during community assembly. A positive association between genetic diversity of S. nutans, C. fasciculata, and S. integrifolium and species diversity established the role of an intraspecific biotic filter during community assembly. Our results show that intraspecific variation in dominant species can significantly modulate the genetic composition of subordinate species.

Convergent and contingent community responses to grass source and dominance during prairie restoration across a longitudinal gradient

Restoring prairie on formerly cultivated land begins by selecting propagule seed sources and the diversity of species to reintroduce. This study examined the effects of dominant grass propagule source (cultivar vs. non-cultivar) and sown propagule diversity (grass:forb sowing ratio) on plant community structure. Two field experiments were established in Kansas and Illinois consisting of identical split plot designs. Dominant grass source was assigned as the whole-plot factor, and sown dominance of grasses (five levels of seeded grass dominance) as the subplot factor. Species density, cover, and diversity were quantified for 5 years. The effect of dominant grass source on the cover of focal grasses, sown species, and volunteer species was contingent upon location, with variation between dominant grass sources observed exclusively in Kansas. Species density and diversity showed regionally convergent patterns in response to dominant grass source. Contrary to our hypotheses, total species density and diversity were not lower in the presence of grass cultivars, the grass source we had predicted would be more competitive. Sown grass dominance effects on the cover of the focal grass species were contingent upon location resulting from establishment corresponding better to the assigned treatments in Illinois. All other cover groups showed regionally convergent patterns, with lower cover of volunteers and higher cover of sown forbs, diversity, and species density in the lowest sown grass dominance treatment in both sites. Thus, decisions regarding the diversity of propagules to reintroduce had more consequence for plant community structure than cultivar or non-cultivar source of dominant grasses.

No effect of seed source on multiple aspects of ecosystem functioning during ecological restoration: cultivars compared to local ecotypes of dominant grasses

Genetic principles underlie recommendations to use local seed, but a paucity of information exists on the genetic distinction and ecological consequences of using different seed sources in restorations. We established a field experiment to test whether cultivars and local ecotypes of dominant prairie grasses were genetically distinct and differentially influenced ecosystem functioning. Whole plots were assigned to cultivar and local ecotype grass sources. Three subplots within each whole plot were seeded to unique pools of subordinate species. The cultivar of the increasingly dominant grass, Sorghastrum nutans, was genetically different than the local ecotype, but genetic diversity was similar between the two sources. There were no differences in aboveground net primary production, soil carbon accrual, and net nitrogen mineralization rate in soil between the grass sources. Comparable productivity of the grass sources among the species pools for four years shows functional equivalence in terms of biomass production. Subordinate species comprised over half the aboveground productivity, which may have diluted the potential for documented trait differences between the grass sources to influence ecosystem processes. Regionally developed cultivars may be a suitable alternative to local ecotypes for restoration in fragmented landscapes with limited gene flow between natural and restored prairie and negligible recruitment by seed.

Environmental metabolomics links genotype to phenotype and predicts genotype abundance in wild plant populations

'The Holy Grail' of plant ecology is to uncover rules that associate species and traits with environmental constraints, community composition and subsequent ecosystem functioning. These aims have been crystallized in recent years within the context of global climate change and environmental pollution, increasing the urgency of the need to predict how vegetation will respond across spatial scales. We investigated whether genetic diversity is associated with the way in which phenotypic plasticity within plant populations is realized and whether this is related to genotype abundance. We used environmental metabolomics to demonstrate biochemical variation between co-occurring genotypes of Carex caryophyllea L. A novel combined metabolomic/functional trait analysis was used to test the functionality of this variation in governing plasticity to variation in edaphic conditions, with particular reference to metabolic pathways that play important roles in growth-related traits. We show that genetic diversity within a wild C. caryophyllea population relates to differences in metabolic composition and functional traits in response to soil nutrient variation, influencing genotype abundance within a community. Our findings highlight the vital role genetic diversity plays within a population in facilitating plant phenotypic plasticity and the potential usefulness of environmental metabolomics to future ecological studies.

Genomic and resistance gene homolog diversity of the dominant tallgrass prairie species across the U.S. Great Plains precipitation gradient

BACKGROUND

Environmental variables such as moisture availability are often important in determining species prevalence and intraspecific diversity. The population genetic structure of dominant plant species in response to a cline of these variables has rarely been addressed. We evaluated the spatial genetic structure and diversity of Andropogon gerardii populations across the U.S. Great Plains precipitation gradient, ranging from approximately 48 cm/year to 105 cm/year.

METHODOLOGY/PRINCIPAL FINDINGS

Genomic diversity was evaluated with AFLP markers and diversity of a disease resistance gene homolog was evaluated by PCR-amplification and digestion with restriction enzymes. We determined the degree of spatial genetic structure using Mantel tests. Genomic and resistance gene homolog diversity were evaluated across prairies using Shannon's index and by averaging haplotype dissimilarity. Trends in diversity across prairies were determined using linear regression of diversity on average precipitation for each prairie. We identified significant spatial genetic structure, with genomic similarity decreasing as a function of distance between samples. However, our data indicated that genome-wide diversity did not vary consistently across the precipitation gradient. In contrast, we found that disease resistance gene homolog diversity was positively correlated with precipitation.

SIGNIFICANCE

Prairie remnants differ in the genetic resources they maintain. Selection and evolution in this disease resistance homolog is environmentally dependent. Overall, we found that, though this environmental gradient may not predict genomic diversity, individual traits such as disease resistance genes may vary significantly.

Genetic variation among Stipa grandis P. Smirn populations with different durations of fencing in the Inner Mongolian Steppe

The genetic structure of a population should be carefully considered in ecological restoration because it may play a critical role in maintaining the persistence of a restored ecosystem. In the present study, we examined genetic diversity and genetic structure of Stipa grandis P. Smirn populations from fenced and grazed plots using amplified fragment length polymorphism markers. Molecular genetic variation showed that the genetic diversity of the fenced populations was greater compared with the overgrazed population. There was a significant variation among the populations (Fst = 0.3689, P < 0.001) by AMOVA analysis, and the gene flow was 0.4039 among the populations. The results from a comparison of limited morphological characteristics and from an unweighted pair group method with arithmetic mean cluster analysis and non-metric multi-dimensional scaling analysis suggested that genetic differentiation had occurred between the fenced populations and the grazed populations. The largest genetic diversity was in the moderately grazed population, which might be related to higher population density and greater sexual reproduction due to less disturbances in the plots. The genetic diversity of the long-term (24 years) fenced population was similar to that of a short-term fenced population (fenced for 11 years). These results suggested that the genetic diversity in the overgrazed population might be increased to some extent through fencing, but this effect did not occur beyond 11 years.

Genetic diversity and population differentiation of the dominant species Stipa krylovii in the Inner Mongolia Steppe

Random amplified polymorphic DNA was used to assess the level of genetic diversity and genetic structure of Stipa krylovii (Gramineae), an important dominant species in the northern grasslands of China. Genetic diversity was low within S. krylovii populations, and diversity at the population level was associated with precipitation and cumulative temperature variations. There was much genetic differentiation among populations and among habitats as well. A Mantel test indicated no significant correlation between genetic distance and geographic distance of populations. A nonmetric multidimensional scaling analysis revealed some spatial relationships among the 90 individuals in a two-dimensional plot. Habitat fragmentation and degradation throughout the geographic range of S. krylovii could account for the low genetic diversity and high genetic differentiation of the species. Such information will be useful for conservation managers trying to plan an effective strategy to protect this important species.

Genetic fingerprinting of germplasm accessions as an aid for species conservation: a case study with Borderea chouardii (Dioscoreaceae), one of the most critically endangered Iberian plants

BACKGROUND AND AIMS

Molecular markers have changed previous expectations about germplasm collections of endangered plants, as new perspectives aim at holding a significant representation of all the genetic diversity in the studied species to accomplish further conservation initiatives successfully. Borderea chouardii is a critically endangered allotetraploid dioecious member of Dioscoreaceae, known from a single population in the Iberian pre-Pyrenees. This population was reported to be highly structured into two genetically distinct groups of individuals corresponding to their spatial separation along the vertical cliff where it grows. In 1999, the Spanish Government of Aragón launched the first conservation programme for the ex situ preservation of this species, and since then a seed collection has been conserved at the Germplasm Bank of the Universidad Politécnica de Madrid. However, as some seed samples had not been labelled clearly at the time of collection, their origin was uncertain.

METHODS

Genetic variation in germplasm accessions of B. chouardii was investigated using microsatellite (simple sequence repeat; SSR) markers.

KEY RESULTS

The 17 primer pairs used detected 62 SSR alleles in the 46 samples analysed from five different germplasm stocks. Eight alleles scored from the wild population were not detected in the germplasm samples analysed. The relatedness of the germplasm samples to the wild subpopulations through neighbour-joining clustering, principal coordinates analysis (PCO) and assignment tests revealed a biased higher representation of the genetic diversity of the lower cliff (43 samples) subpopulation than that of the upper cliff (three samples).

CONCLUSIONS

The collection of additional samples from the upper cliff is recommended to achieve a better representation of the genetic diversity of this subpopulation. It is also recommended that these stocks should be managed separately according to their distinct microspatial origin in order to preserve the genetic substructuring of the wild population.