Bacterial but Not Fungal Rhizosphere Community Composition Differ among Perennial Grass Ecotypes under Abiotic Environmental Stress.

Environmental change, especially frequent droughts, is predicted to detrimentally impact the North American perennial grasslands. Consistent dry spells will affect plant communities as well as their associated rhizobiomes, possibly altering the plant host performance under environmental stress. Therefore, there is a need to understand the impact of drought on the rhizobiome, and how the rhizobiome may modulate host performance and ameliorate its response to drought stress. In this study, we analyzed bacterial and fungal communities in the rhizospheres of three ecotypes (dry, mesic, and wet) of dominant prairie grass, Andropogon gerardii. The ecotypes were established in 2010 in a common garden design and grown for a decade under persistent dry conditions at the arid margin of the species' range in Colby, Kansas. The experiment aimed to answer whether and to what extent do the different ecotypes maintain or recruit distinct rhizobiomes after 10 years in an arid climate. In order to answer this question, we screened the bacterial and fungal rhizobiome profiles of the ecotypes under the arid conditions of western Kansas as a surrogate for future climate environmental stress using 16S rRNA and ITS2 metabarcoding sequencing. Under these conditions, bacterial communities differed compositionally among the A. gerardii ecotypes, whereas the fungal communities did not. The ecotypes were instrumental in driving the differences among bacterial rhizobiomes, as the ecotypes maintained distinct bacterial rhizobiomes even after 10 years at the edge of the host species range. This study will aid us to optimize plant productivity through the use of different ecotypes under future abiotic environmental stress, especially drought. IMPORTANCE In this study, we used a 10-year long reciprocal garden system, and reports that different ecotypes (dry, mesic, and wet) of dominant prairie grass, Andropogon gerardii can maintain or recruit distinct bacterial but not fungal rhizobiomes after 10 years in an arid environment. We used both 16S rRNA and ITS2 amplicons to analyze the bacterial and fungal communities in the rhizospheres of the respective ecotypes. We showed that A. gerardii might regulate the bacterial community to adapt to the arid environment, in which some ecotypes were not adapted to. Our study also suggested a possible tradeoff between the generalist and the specialist bacterial communities in specific environments, which could benefit the plant host. Our study will provide insights into the plant host regulation of the rhizosphere bacterial and fungal communities, especially during frequent drought conditions anticipated in the future.

Selection of suitable reference genes for quantitive real-time PCR normalization in Miscanthus lutarioriparia.

Miscanthus lutarioriparia, which is found widespread in China, has attracted great attention as a most potential bioenergy plant for years. The quantitative real time PCR (RT-qPCR) has appeared as a sensitive and powerful technique to measure gene expression in living organisms during different development stages. In this study, we evaluated ten candidate genes, including 25S ribosomal RNA gene (25S rRNA), actin1 gene (ACT1), carotenoid-binding protein 20 gene (CBP20), glyceraldehyde-3-phosphate dehydrogenase gene (GAPDH), Ubiquitin gene (UBQ), eukaryotic elongation factor 1-αgene (eEF-1α), α-tubulin gene (α-TUB), ß-tubulin gene (ß-TUB), eukaryotic translation initiation factor 4α-1 gene (eIF-4α) and NAC domain protein gene(NAC) in a series of 30 M. lutarioriparia samples followed by statistical algorithms geNorm and Normfinder to analyze the gene expression stability. The results indicated that eIF-4αand UBQ were the most stable expressed genes while CBP20 showed as the least stable among all the samples. Based on above research, we recommend that at least two top-ranked reference genes should be employed for expression data normalization. The best genes selected in this study will provide a starting point to select reference genes in the future in other tissues and under other experimental conditions in this energy crop candidate.

Interspecific genetic maps in Miscanthus floridulus and M. sacchariflorus accelerate detection of QTLs associated with plant height and inflorescence.

Miscanthus is recognized as a promising lignocellulosic crop for the production of bioethanol and bioproducts worldwide. To facilitate the identification of agronomical important traits and establish genetics knowledge, two genetic maps were developed from a controlled interspecific cross between M. floridulus and M. sacchariflorus. A total of 650 SSR markers were mapped in M. floridulus, spanning 19 linkage groups and 2053.31 cM with an average interval of 3.25 cM. The map of M. sacchariflorus comprised 495 SSR markers in 19 linkage groups covering 1684.86 cM with an average interval of 3.54 cM. The estimation on genome length indicated that the genome coverage of parental genetic maps were 93.87% and 89.91%, respectively. Eighty-eight bi-parental common markers were allowed to connect the two maps, and six pairs of syntenic linkage groups were recognized. Furthermore, quantitative trait loci (QTL) mapping of three agronomic traits, namely, plant height (PH), heading time (HT), and flowering time (FT), demonstrated that a total of 66 QTLs were identified in four consecutive years using interval mapping and multiple-QTL model. The LOD value of these QTLs ranged from 2.51 to 10.60, and the phenotypic variation explained varied from 9.50 to 37.10%. QTL cluster in syntenic groups MF19/MS7 contained six stable QTLs associated with PH, HT, and FT. In conclusion, we report for the first time the genetic mapping of biomass traits in M. floridulus and M. sacchariflorus. These results will be a valuable genetic resource, facilitating the discovery of essential genes and breeding of Miscanthus.

Coexpression network revealing the plasticity and robustness of population transcriptome during the initial stage of domesticating energy crop Miscanthus lutarioriparius.

KEY MESSAGE: Coexpression network revealing genes with Co-variation Expression pattern (CE) and those with Top rank of Expression fold change (TE) played different roles in responding to new environment of Miscanthus lutarioriparius. Variation in gene expression level, the product of genetic and/or environmental perturbation, determines the robustness-to-plasticity spectrum of a phenotype in plants. Understanding how expression variation of plant population response to a new field is crucial to domesticate energy crops. Weighted Gene Coexpression Network Analysis (WGCNA) was used to explore the patterns of expression variation based on 72 Miscanthus lutarioriparius transcriptomes from two contrasting environments, one near the native habitat and the other in one harsh domesticating region. The 932 genes with Co-variation Expression pattern (CE) and other 932 genes with Top rank of Expression fold change (TE) were identified and the former were strongly associated with the water use efficiency (r ≥ 0.55, P ≤ 10-7). Functional enrichment of CE genes were related to three organelles, which well matched the annotation of twelve motifs identified from their conserved noncoding sequence; while TE genes were mostly related to biotic and/or abiotic stress. The expression robustness of CE genes with high genetic diversity kept relatively stable between environments while the harsh environment reduced the expression robustness of TE genes with low genetic diversity. The expression plasticity of CE genes was increased less than that of TE genes. These results suggested that expression variation of CE genes and TE genes could account for the robustness and plasticity of acclimation ability of Miscanthus, respectively. The patterns of expression variation revealed by transcriptomic network would shed new light on breeding and domestication of energy crops.

NMR-based Metabolomics to Study the Cold-acclimation Strategy of Two Miscanthus Genotypes.

INTRODUCTION: Abiotic stress is a major cause of yield loss in plant culture. Miscanthus, a perennial C4 grass, is now considered a major source of renewable energy, especially for biofuel production. During the first year of planting in Northern Europe, Miscanthus was exposed to frost temperature, which generated high mortality in young plants and large loss of yield. One strategy to avoid such loss is to apply cold-acclimation, which confers on plants a better resistance to low temperature. OBJECTIVES: The aim of this study is to describe the effect of a cold-acclimation period on the metabolome of two Miscanthus genotypes that vary in their frost sensitivity at the juvenile stage. Miscanthus × giganteus (GIG) is particularly sensitive to frost, whereas Miscanthus sinensis August Feder (AUG) is tolerant. MATERIALS AND METHODS: Polar metabolite extraction was performed on Miscanthus, grown in non-acclimated or cold-acclimated conditions. Extracts were analysed by 1 H-NMR followed by multivariate statistical analysis. Discriminant metabolites were identified. RESULTS: More than 40 metabolites were identified in the two Miscanthus genotypes. GIG and AUG showed a different metabolic background before cold treatment, probably related to their genetic background. After cold-acclimation, GIG and AUG metabolomes remained different. The tolerant genotype showed notably higher levels of accumulation in proline, sucrose and maltose when subjected to cold. CONCLUSION: These two genotypes seem to have a different adaptation strategy in cold conditions. The studied change in the metabolome concerns different types of molecules related to the cold-tolerant behaviour of Miscanthus. Copyright © 2016 John Wiley & Sons, Ltd.

A Genomic Resource for the Development, Improvement, and Exploitation of Sorghum for Bioenergy.

With high productivity and stress tolerance, numerous grass genera of the Andropogoneae have emerged as candidates for bioenergy production. To optimize these candidates, research examining the genetic architecture of yield, carbon partitioning, and composition is required to advance breeding objectives. Significant progress has been made developing genetic and genomic resources for Andropogoneae, and advances in comparative and computational genomics have enabled research examining the genetic basis of photosynthesis, carbon partitioning, composition, and sink strength. To provide a pivotal resource aimed at developing a comparative understanding of key bioenergy traits in the Andropogoneae, we have established and characterized an association panel of 390 racially, geographically, and phenotypically diverse Sorghum bicolor accessions with 232,303 genetic markers. Sorghum bicolor was selected because of its genomic simplicity, phenotypic diversity, significant genomic tools, and its agricultural productivity and resilience. We have demonstrated the value of sorghum as a functional model for candidate gene discovery for bioenergy Andropogoneae by performing genome-wide association analysis for two contrasting phenotypes representing key components of structural and non-structural carbohydrates. We identified potential genes, including a cellulase enzyme and a vacuolar transporter, associated with increased non-structural carbohydrates that could lead to bioenergy sorghum improvement. Although our analysis identified genes with potentially clear functions, other candidates did not have assigned functions, suggesting novel molecular mechanisms for carbon partitioning traits. These results, combined with our characterization of phenotypic and genetic diversity and the public accessibility of each accession and genomic data, demonstrate the value of this resource and provide a foundation for future improvement of sorghum and related grasses for bioenergy production.

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.

Inferential considerations for low-count RNA-seq transcripts: a case study on the dominant prairie grass Andropogon gerardii.

BACKGROUND: Differential expression (DE) analysis of RNA-seq data still poses inferential challenges, such as handling of transcripts characterized by low expression levels. In this study, we use a plasmode-based approach to assess the relative performance of alternative inferential strategies on RNA-seq transcripts, with special emphasis on transcripts characterized by a small number of read counts, so-called low-count transcripts, as motivated by an ecological application in prairie grasses. Big bluestem (Andropogon gerardii) is a wide-ranging dominant prairie grass of ecological and agricultural importance to the US Midwest while edaphic subspecies sand bluestem (A. gerardii ssp. Hallii) grows exclusively on sand dunes. Relative to big bluestem, sand bluestem exhibits qualitative phenotypic divergence consistent with enhanced drought tolerance, plausibly associated with transcripts of low expression levels. Our dataset consists of RNA-seq read counts for 25,582 transcripts (60% of which are classified as low-count) collected from leaf tissue of individual plants of big bluestem (n = 4) and sand bluestem (n = 4). Focused on low-count transcripts, we compare alternative ad-hoc data filtering techniques commonly used in RNA-seq pipelines and assess the inferential performance of recently developed statistical methods for DE analysis, namely DESeq2 and edgeR robust. These methods attempt to overcome the inherently noisy behavior of low-count transcripts by either shrinkage or differential weighting of observations, respectively. RESULTS: Both DE methods seemed to properly control family-wise type 1 error on low-count transcripts, whereas edgeR robust showed greater power and DESeq2 showed greater precision and accuracy. However, specification of the degree of freedom parameter under edgeR robust had a non-trivial impact on inference and should be handled carefully. When properly specified, both DE methods showed overall promising inferential performance on low-count transcripts, suggesting that ad-hoc data filtering steps at arbitrary expression thresholds may be unnecessary. A note of caution is in order regarding the approximate nature of DE tests under both methods. CONCLUSIONS: Practical recommendations for DE inference are provided when low-count RNA-seq transcripts are of interest, as is the case in the comparison of subspecies of bluestem grasses. Insights from this study may also be relevant to other applications focused on transcripts of low expression levels.

Environmental correlates of cytotype distribution in Andropogon gerardii (Poaceae).

UNLABELLED: • PREMISE OF THE STUDY: Information about geographic distribution of cytotypes can provide insight into the origin and maintenance of autopolyploid complexes and builds a foundation for understanding cytotype differentiation and the dynamics of mixed-ploidy populations. Here, we investigate environmental correlates of the geographic distributions of 6x and 9x individuals in the ecologically dominant grass Andropogon gerardii to examine the role of climate in shaping patterns of cytotype distribution in this species.• METHODS: Flow cytometry was used to estimate ploidy level in 352 individuals from 32 populations across North America. Ecological differentiation of cytotypes was tested by relating BIOCLIM variables to cytotype distribution using principal components analysis and partial linear regression.• KEY RESULTS: Broad geographic sampling confirmed two primary cytotypes-6x (hexaploid) and 9x (enneaploid)-and revealed that 9x plants are more common than previously thought. Enneaploids occur frequently in the southern portions of the range, with hexaploids dominating in northern regions. Mixed-ploidy populations were common (46.9%). Principal components analysis and partial linear regression indicated that reduced summer precipitation and increased variation in diurnal and seasonal temperature range were significant predictors of the frequency of 9x plants in a population.• CONCLUSIONS: Results indicate that (1) geographic distribution of 6x and 9x individuals is nonrandom; (2) environmental variables are associated with cytotype distribution in A. gerardii; and (3) nearly half of populations surveyed include both 6x and 9x individuals. The persistence of mixed-ploidy populations may reflect a combination of recurrent polyploid formation and the prevalence of clonal reproduction.

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.

Environmental and genetic variation in leaf anatomy among populations of Andropogon gerardii (Poaceae) along a precipitation gradient.

PREMISE OF THE STUDY: Phenotypes of two Andropogon gerardii subspecies, big bluestem and sand bluestem, vary throughout the prairie ecosystem of North America. This study sought to determine the role of genetics and environment in driving adaptive variation of leaf structure in big bluestem and sand bluestem. • METHODS: Four populations of big bluestem and one population of sand bluestem were planted in common gardens at four sites across a precipitation gradient from western Kansas to southern Illinois. Internal leaf structure and trichome density of A. gerardii were examined by light microscopy to separate genetic and environmentally controlled traits. Leaf thickness, midrib thickness, bulliform cells, interveinal distance, vein size, and trichome density were quantified. • KEY RESULTS: At all planting sites, sand bluestem and the xeric population of A. gerardii had thicker leaves and fewer bulliform cells compared with mesic populations. Environment and genetic source population were both influential for leaf anatomy. Leaves from plants grown in mesic sites (Carbondale, Illinois and Manhattan, Kansas) had thicker midribs, larger veins, fewer trichomes, and a greater proportion of bulliform cells compared to plants grown in drier sites (Colby and Hays, Kansas). • CONCLUSIONS: Water availability has driven adaptive variation in leaf structure in populations of A. gerardii, particularly between sand bluestem and big bluestem. Genetically based differences in leaves of A. gerardii indicate adaptive variation and evolutionary forces differentiating sand bluestem from big bluestem. Environmental responses of A. gerardii leaves suggest an ability to adjust to drought, even in populations adapted to mesic home environments.

Genetic diversity of a dominant C4 grass is altered with increased precipitation variability.

Climate change has the potential to alter the genetic diversity of plant populations with consequences for community dynamics and ecosystem processes. Recent research focused on changes in climatic means has found evidence of decreased precipitation amounts reducing genetic diversity. However, increased variability in climatic regimes is also predicted with climate change, but the effects of this aspect of climate change on genetic diversity have yet to be investigated. After 10 years of experimentally increased intra-annual variability in growing season precipitation regimes, we report that the number of genotypes of the dominant C(4) grass, Andropogon gerardii Vitman, has been significantly reduced in native tallgrass prairie compared with unmanipulated prairie. However, individuals showed a different pattern of genomic similarity with increased precipitation variability resulting in greater genome dissimilarity among individuals when compared to unmanipulated prairie. Further, we found that genomic dissimilarity was positively correlated with aboveground productivity in this system. The increased genomic dissimilarity among individuals in the altered treatment alongside evidence for a positive correlation of genomic dissimilarity with phenotypic variation suggests ecological sorting of genotypes may be occurring via niche differentiation. Overall, we found effects of more variable precipitation regimes on population-level genetic diversity were complex, emphasizing the need to look beyond genotype numbers for understanding the impacts of climate change on genetic diversity. Recognition that future climate change may alter aspects of genetic diversity in different ways suggests possible mechanisms by which plant populations may be able to retain a diversity of traits in the face of declining biodiversity.

Invasion of an intact plant community: the role of population versus community level diversity.

To improve the understanding of how native plant diversity influences invasion, we examined how population and community diversity may directly and indirectly be related to invasion in a natural field setting. Due to the large impact of the dominant C(4) grass species (Andropogon gerardii) on invasion resistance of tallgrass prairie, we hypothesized that genetic diversity and associated traits within a population of this species would be more strongly related to invasion than diversity or traits of the rest of the community. We added seeds of the exotic invasive C(4) grass, A. bladhii, to 1-m(2) plots in intact tallgrass prairie that varied in genetic diversity of A. gerardii and plant community diversity, but not species richness. We assessed relationships among genetic diversity and traits of A. gerardii, community diversity, community aggregated traits, resource availability, and early season establishment and late-season persistence of the invader using structural equation modeling (SEM). SEM models suggested that community diversity likely enhanced invasion indirectly through increasing community aggregated specific leaf area as a consequence of more favorable microclimatic conditions for seedling establishment. In contrast, neither population nor community diversity was directly or indirectly related to late season survival of invasive seedlings. Our research suggests that while much of diversity-invasion research has separately focused on the direct effects of genetic and species diversity, when taken together, we find that the role of both levels of diversity on invasion resistance may be more complex, whereby effects of diversity may be primarily indirect via traits and vary depending on the stage of invasion.

Genomic and small RNA sequencing of Miscanthus x giganteus shows the utility of sorghum as a reference genome sequence for Andropogoneae grasses.

BACKGROUND: Miscanthus x giganteus (Mxg) is a perennial grass that produces superior biomass yields in temperate environments. The essentially uncharacterized triploid genome (3n = 57, x = 19) of Mxg is likely critical for the rapid growth of this vegetatively propagated interspecific hybrid. RESULTS: A survey of the complex Mxg genome was conducted using 454 pyrosequencing of genomic DNA and Illumina sequencing-by-synthesis of small RNA. We found that the coding fraction of the Mxg genome has a high level of sequence identity to that of other grasses. Highly repetitive sequences representing the great majority of the Mxg genome were predicted using non-cognate assembly for de novo repeat detection. Twelve abundant families of repeat were observed, with those related to either transposons or centromeric repeats likely to comprise over 95% of the genome. Comparisons of abundant repeat sequences to a small RNA survey of three Mxg organs (leaf, rhizome, inflorescence) revealed that the majority of observed 24-nucleotide small RNAs are derived from these repetitive sequences. We show that high-copy-number repeats match more of the small RNA, even when the amount of the repeat sequence in the genome is accounted for. CONCLUSIONS: We show that major repeats are present within the triploid Mxg genome and are actively producing small RNAs. We also confirm the hypothesized origins of Mxg, and suggest that while the repeat content of Mxg differs from sorghum, the sorghum genome is likely to be of utility in the assembly of a gene-space sequence of Mxg.