Genomic mechanisms of climate adaptation in polyploid bioenergy switchgrass

Characterization of the m6A gene family in switchgrass and functional analysis of PvALKBH10 during flowering

N6-methyladenosine (m6A), a nucleotide modification that is frequently seen in RNA, plays a crucial role in plant growth, development and stress resistance. However, the m6A regulatory machinery in switchgrass (Panicum virgatum L.), a model plant for cellulose-to-ethanol conversion, remains largely unknown. In this study, we identified 57 candidate genes involved in m6A-regulation in the switchgrass genome, and analyzed their chromosomal distribution, evolutionary relationships, and functions. Notably, we observed distinct gene expression patterns under salt and drought stress, with salt stress inducing writer and eraser genes, alongside drought stress predominantly affecting reader genes. Additionally, we knocked out PvALKBH10, an m6A demethylase gene, via CRISPR/Cas9 and found its potential function in controlling flowering time. This study provides insight into the genomic organization and evolutionary features of m6A-associated putative genes in switchgrass, and therefore serves as the basis for further functional studies.

miR156-PvSPL2 controls culm development by transcriptional repression of switchgrass CYTOKININ OXIDASE/DEHYDROGENASE4

Culm development in grasses can be controlled by both miR156 and cytokinin. However, the crosstalk between the miR156-SPL module and the cytokinin metabolic pathway remains largely unknown. Here, we found CYTOKININ OXIDASE/DEHYDROGENASE4 (PvCKX4) plays a negative regulatory role in culm development of the bioenergy grass Panicum virgatum (switchgrass). Overexpression of PvCKX4 in switchgrass reduced the internode diameter and length without affecting tiller number. Interestingly, we also found that PvCKX4 was always upregulated in miR156 overexpressing (miR156OE) transgenic switchgrass lines. Additionally, upregulation of either miR156 or PvCKX4 in switchgrass reduced the content of isopentenyl adenine (iP) without affecting trans-zeatin (tZ) accumulation. It is consistent with the evidence that the recombinant PvCKX4 protein exhibited much higher catalytic activity against iP than tZ in vitro. Furthermore, our results showed that miR156-targeted SPL2 bound directly to the promoter of PvCKX4 to repress its expression. Thus, alleviating the SPL2-mediated transcriptional repression of PvCKX4 through miR156 overexpression resulted in a significant increase in cytokinin degradation and impaired culm development in switchgrass. On the contrary, suppressing PvCKX4 in miR156OE transgenic plants restored iP content, internode diameter, and length to wild-type levels. Most strikingly, the double transgenic lines retained the same increased tiller numbers as the miR156OE transgenic line, which yielded more biomass than the wild type. These findings indicate that the miR156-SPL module can control culm development through transcriptional repression of PvCKX4 in switchgrass, which provides a promising target for precise design of shoot architecture to yield more biomass from grasses.

Current status of community resources and priorities for weed genomics research

Weeds are attractive models for basic and applied research due to their impacts on agricultural systems and capacity to swiftly adapt in response to anthropogenic selection pressures. Currently, a lack of genomic information precludes research to elucidate the genetic basis of rapid adaptation for important traits like herbicide resistance and stress tolerance and the effect of evolutionary mechanisms on wild populations. The International Weed Genomics Consortium is a collaborative group of scientists focused on developing genomic resources to impact research into sustainable, effective weed control methods and to provide insights about stress tolerance and adaptation to assist crop breeding.

C4 cereal and biofuel crop microbiomes

Microbiomes provide multiple life-support functions for plants, including nutrient acquisition and tolerance to abiotic and biotic stresses. Considering the importance of C4 cereal and biofuel crops for food security under climate change conditions, more attention has been given recently to C4 plant microbiome assembly and functions. Here, we review the current status of C4 cereal and biofuel crop microbiome research with a focus on beneficial microbial traits for crop growth and health. We highlight the importance of environmental factors and plant genetics in C4 crop microbiome assembly and pinpoint current knowledge gaps. Finally, we discuss the potential of foxtail millet as a C4 model species and outline future perspectives of C4 plant microbiome research.

Genetically correlated leaf tensile and morphological traits are driven by growing season length in a widespread perennial grass

PREMISE

Leaf tensile resistance, a leaf's ability to withstand pulling forces, is an important determinant of plant ecological strategies. One potential driver of leaf tensile resistance is growing season length. When growing seasons are long, strong leaves, which often require more time and resources to construct than weak leaves, may be more advantageous than when growing seasons are short. Growing season length and other ecological conditions may also impact the morphological traits that underlie leaf tensile resistance.

METHODS

To understand variation in leaf tensile resistance, we measured size-dependent leaf strength and size-independent leaf toughness in diverse genotypes of the widespread perennial grass Panicum virgatum (switchgrass) in a common garden. We then used quantitative genetic approaches to estimate the heritability of leaf tensile resistance and whether there were genetic correlations between leaf tensile resistance and other morphological traits.

RESULTS

Leaf tensile resistance was positively associated with aboveground biomass (a proxy for fitness). Moreover, both measures of leaf tensile resistance exhibited high heritability and were positively genetically correlated with leaf lamina thickness and leaf mass per area (LMA). Leaf tensile resistance also increased with the growing season length in the habitat of origin, and this effect was mediated by both LMA and leaf thickness.

CONCLUSIONS

Differences in growing season length may promote selection for different leaf lifespans and may explain existing variation in leaf tensile resistance in P. virgatum. In addition, the high heritability of leaf tensile resistance suggests that P. virgatum will be able to respond to climate change as growing seasons lengthen.

N6-methyladenosine analysis unveils key mechanisms underlying long-term salt stress tolerance in switchgrass (Panicum virgatum)

N6-methyladenosine (m6A) RNA modification is critical for plant growth, development, and environmental stress response. While short-term stress impacts on m6A are well-documented, the consequences of prolonged stress remain underexplored. This study conducts a thorough transcriptome-wide analysis of m6A modifications following 28-day exposure to 200 mM NaCl. We detected 11,149 differentially expressed genes (DEGs) and 12,936 differentially methylated m6A peaks, along with a global decrease in m6A levels. Notably, about 62% of m6A-modified DEGs, including demethylase genes like PvALKBH6_N, PvALKBH9_K, and PvALKBH10_N, showed increased expression and reduced m6A peaks, suggesting that decreased m6A methylation may enhance gene expression under salt stress. Consistent expression and methylation patterns were observed in key genes related to ion homeostasis (e.g., H+-ATPase 1, High-affinity K+transporter 5), antioxidant defense (Catalase 1/2, Copper/zinc superoxide dismutase 2, Glutathione synthetase 1), and osmotic regulation (delta 1-pyrroline-5-carboxylate synthase 2, Pyrroline-5-carboxylate reductase). These findings provide insights into the adaptive mechanisms of switchgrass under long-term salt stress and highlight the potential of regulating m6A modifications as a novel approach for crop breeding strategies focused on stress resistance.

Allopolyploid subgenome identification and implications for evolutionary analysis

Whole-genome duplications (WGDs) are widespread genomic events in eukaryotes that are hypothesized to contribute to the evolutionary success of many lineages, including flowering plants, Saccharomyces yeast, and vertebrates. WGDs generally can be classified into autopolyploids (ploidy increase descended from one species) or allopolyploids (ploidy increase descended from multiple species). Assignment of allopolyploid progenitor species (called subgenomes in the polyploid) is important to understanding the biology and evolution of polyploids, including the asymmetric subgenome evolution following hybridization (biased fractionation). Here, I review the different methodologies used to identify the ancestors of allopolyploid subgenomes, discuss the advantages and disadvantages of these methods, and outline the implications of how these methods affect the subsequent evolutionary analysis of these genomes.

Genome-wide profiling of histone (H3) lysine 4 (K4) tri-methylation (me3) under drought, heat, and combined stresses in switchgrass

BACKGROUND

Switchgrass (Panicum virgatum L.) is a warm-season perennial (C4) grass identified as an important biofuel crop in the United States. It is well adapted to the marginal environment where heat and moisture stresses predominantly affect crop growth. However, the underlying molecular mechanisms associated with heat and drought stress tolerance still need to be fully understood in switchgrass. The methylation of H3K4 is often associated with transcriptional activation of genes, including stress-responsive. Therefore, this study aimed to analyze genome-wide histone H3K4-tri-methylation in switchgrass under heat, drought, and combined stress.

RESULTS

In total, ~ 1.3 million H3K4me3 peaks were identified in this study using SICER. Among them, 7,342; 6,510; and 8,536 peaks responded under drought (DT), drought and heat (DTHT), and heat (HT) stresses, respectively. Most DT and DTHT peaks spanned 0 to + 2000 bases from the transcription start site [TSS]. By comparing differentially marked peaks with RNA-Seq data, we identified peaks associated with genes: 155 DT-responsive peaks with 118 DT-responsive genes, 121 DTHT-responsive peaks with 110 DTHT-responsive genes, and 175 HT-responsive peaks with 136 HT-responsive genes. We have identified various transcription factors involved in DT, DTHT, and HT stresses. Gene Ontology analysis using the AgriGO revealed that most genes belonged to biological processes. Most annotated peaks belonged to metabolite interconversion, RNA metabolism, transporter, protein modifying, defense/immunity, membrane traffic protein, transmembrane signal receptor, and transcriptional regulator protein families. Further, we identified significant peaks associated with TFs, hormones, signaling, fatty acid and carbohydrate metabolism, and secondary metabolites. qRT-PCR analysis revealed the relative expressions of six abiotic stress-responsive genes (transketolase, chromatin remodeling factor-CDH3, fatty-acid desaturase A, transmembrane protein 14C, beta-amylase 1, and integrase-type DNA binding protein genes) that were significantly (P < 0.05) marked during drought, heat, and combined stresses by comparing stress-induced against un-stressed and input controls.

CONCLUSION

Our study provides a comprehensive and reproducible epigenomic analysis of drought, heat, and combined stress responses in switchgrass. Significant enrichment of H3K4me3 peaks downstream of the TSS of protein-coding genes was observed. In addition, the cost-effective experimental design, modified ChIP-Seq approach, and analyses presented here can serve as a prototype for other non-model plant species for conducting stress studies.

Plant pangenomes for crop improvement, biodiversity and evolution

Plant genome sequences catalogue genes and the genetic elements that regulate their expression. Such inventories further research aims as diverse as mapping the molecular basis of trait diversity in domesticated plants or inquiries into the origin of evolutionary innovations in flowering plants millions of years ago. The transformative technological progress of DNA sequencing in the past two decades has enabled researchers to sequence ever more genomes with greater ease. Pangenomes - complete sequences of multiple individuals of a species or higher taxonomic unit - have now entered the geneticists' toolkit. The genomes of crop plants and their wild relatives are being studied with translational applications in breeding in mind. But pangenomes are applicable also in ecological and evolutionary studies, as they help classify and monitor biodiversity across the tree of life, deepen our understanding of how plant species diverged and show how plants adapt to changing environments or new selection pressures exerted by human beings.

Widespread incomplete lineage sorting and introgression shaped adaptive radiation in the Gossypium genus

Cotton (Gossypium) stands as a crucial economic crop, serving as the primary source of natural fiber for the textile sector. However, the evolutionary mechanisms driving speciation within the Gossypium genus remain unresolved. In this investigation, we leveraged 25 Gossypium genomes and introduced four novel assemblies-G. harknessii, G. gossypioides, G. trilobum, and G. klotzschianum (Gklo)-to delve into the speciation history of this genus. Notably, we encountered intricate phylogenies potentially stemming from introgression. These complexities are further compounded by incomplete lineage sorting (ILS), a factor likely to have been instrumental in shaping the swift diversification of cotton. Our focus subsequently shifted to the rapid radiation episode during a concise period in Gossypium evolution. For a recently diverged lineage comprising G. davidsonii, Gklo, and G. raimondii, we constructed a finely detailed ILS map. Intriguingly, this analysis revealed the non-random distribution of ILS regions across the reference Gklo genome. Moreover, we identified signs of robust natural selection influencing specific ILS regions. Noteworthy variations pertaining to speciation emerged between the closely related sister species Gklo and G. davidsonii. Approximately 15.74% of speciation structural variation genes and 12.04% of speciation-associated genes were estimated to intersect with ILS signatures. These findings enrich our understanding of the role of ILS in adaptive radiation, shedding fresh light on the intricate speciation history of the Gossypium genus.

Economic and Sustainability Impacts of Yield and Composition Variation in Bioenergy Crops: Switchgrass (Panicum virgatum L.)

Economically viable production of biobased products and fuels requires high-yielding, high-quality, sustainable process-advantaged crops, developed using bioengineering or advanced breeding approaches. Understanding which crop phenotypic traits have the largest impact on biofuel economics and sustainability outcomes is important for the targeted feedstock crop development. Here, we evaluated biomass yield and cell-wall composition traits across a large natural variant population of switchgrass (Panicum virgatum L.) grown across three common garden sites. Samples from 331 switchgrass genotypes were collected and analyzed for carbohydrate and lignin components. Considering plant survival and biomass after multiple years of growth, we found that 84 of the genotypes analyzed may be suited for commercial production in the southeastern U.S. These genotypes show a range of growth and compositional traits across the population that are apparently independent of each other. We used these data to conduct techno-economic analyses and life cycle assessments evaluating the performance of each switchgrass genotype under a standard cellulosic ethanol process model with pretreatment, added enzymes, and fermentation. We find that switchgrass yield per area is the largest economic driver of the minimum fuel selling price (MSFP), ethanol yield per hectare, global warming potential (GWP), and cumulative energy demand (CED). At any yield, the carbohydrate content is significant but of secondary importance. Water use follows similar trends but has more variability due to an increased dependence on the biorefinery model. Analyses presented here highlight the primary importance of plant yield and the secondary importance of carbohydrate content when selecting a feedstock that is both economical and sustainable.

Plant physical defenses contribute to a latitudinal gradient in resistance to insect herbivory within a widespread perennial grass

PREMISE

Herbivore pressure can vary across the range of a species, resulting in different defensive strategies. If herbivory is greater at lower latitudes, plants may be better defended there, potentially driving a latitudinal gradient in defense. However, relationships that manifest across the entire range of a species may be confounded by differences within genetic subpopulations, which may obscure the drivers of these latitudinal gradients.

METHODS

We grew plants of the widespread perennial grass Panicum virgatum in a common garden that included genotypes from three genetic subpopulations spanning an 18.5° latitudinal gradient. We then assessed defensive strategies of these plants by measuring two physical resistance traits-leaf mass per area (LMA) and leaf ash, a proxy for silica-and multiple measures of herbivory by caterpillars of the generalist herbivore fall armyworm (Spodoptera frugiperda).

RESULTS

Across all genetic subpopulations, low-latitude plants experienced less herbivory than high-latitude plants. Within genetic subpopulations, however, this relationship was inconsistent-the most widely distributed and phenotypically variable subpopulation (Atlantic) exhibited more consistent latitudinal trends than either of the other two subpopulations. The two physical resistance traits, LMA and leaf ash, were both highly heritable and positively associated with resistance to different measures of herbivory across all subpopulations, indicating their importance in defense against herbivores. Again, however, these relationships were inconsistent within subpopulations.

CONCLUSIONS

Defensive gradients that occur across the entire species range may not arise within localized subpopulations. Thus, identifying the drivers of latitudinal gradients in herbivory defense may depend on adequately sampling the diversity within a species.

Pangenome analysis reveals genomic variations associated with domestication traits in broomcorn millet

Broomcorn millet (Panicum miliaceum L.) is an orphan crop with the potential to improve cereal production and quality, and ensure food security. Here we present the genetic variations, population structure and diversity of a diverse worldwide collection of 516 broomcorn millet genomes. Population analysis indicated that the domesticated broomcorn millet originated from its wild progenitor in China. We then constructed a graph-based pangenome of broomcorn millet based on long-read de novo genome assemblies of 32 representative accessions. Our analysis revealed that the structural variations were highly associated with transposable elements, which influenced gene expression when located in the coding or regulatory regions. We also identified 139 loci associated with 31 key domestication and agronomic traits, including candidate genes and superior haplotypes, such as LG1, for panicle architecture. Thus, the study's findings provide foundational resources for developing genomics-assisted breeding programs in broomcorn millet.

Identification of loci associated with water use efficiency and symbiotic nitrogen fixation in soybean

Soybean (Glycine max) production is greatly affected by persistent and/or intermittent droughts in rainfed soybean-growing regions worldwide. Symbiotic N2 fixation (SNF) in soybean can also be significantly hampered even under moderate drought stress. The objective of this study was to identify genomic regions associated with shoot carbon isotope ratio (δ13C) as a surrogate measure for water use efficiency (WUE), nitrogen isotope ratio (δ15N) to assess relative SNF, N concentration ([N]), and carbon/nitrogen ratio (C/N). Genome-wide association mapping was performed with 105 genotypes and approximately 4 million single-nucleotide polymorphism markers derived from whole-genome resequencing information. A total of 11, 21, 22, and 22 genomic loci associated with δ13C, δ15N, [N], and C/N, respectively, were identified in two environments. Nine of these 76 loci were stable across environments, as they were detected in both environments. In addition to the 62 novel loci identified, 14 loci aligned with previously reported quantitative trait loci for different C and N traits related to drought, WUE, and N2 fixation in soybean. A total of 58 Glyma gene models encoding for different genes related to the four traits were identified in the vicinity of the genomic loci.

Genome evolution and initial breeding of the Triticeae grass Leymus chinensis dominating the Eurasian Steppe

Leymus chinensis, a dominant perennial grass in the Eurasian Steppe, is well known for its remarkable adaptability and forage quality. Hardly any breeding has been done on the grass, limiting its potential in ecological restoration and forage productivity. To enable genetic improvement of the untapped, important species, we obtained a 7.85-Gb high-quality genome of L. chinensis with a particularly long contig N50 (318.49 Mb). Its allotetraploid genome is estimated to originate 5.29 million years ago (MYA) from a cross between the Ns-subgenome relating to Psathyrostachys and the unknown Xm-subgenome. Multiple bursts of transposons during 0.433-1.842 MYA after genome allopolyploidization, which involved predominantly the Tekay and Angela of LTR retrotransposons, contributed to its genome expansion and complexity. With the genome resource available, we successfully developed a genetic transformation system as well as the gene-editing pipeline in L. chinensis. We knocked out the monocot-specific miR528 using CRISPR/Cas9, resulting in the improvement of yield-related traits with increases in the tiller number and growth rate. Our research provides valuable genomic resources for Triticeae evolutionary studies and presents a conceptual framework illustrating the utilization of genomic information and genome editing to accelerate the improvement of wild L. chinensis with features such as polyploidization and self-incompatibility.

Genome-Wide Association Studies of Biluochun Tea Plant Populations in Dongting Mountain and Comprehensive Identification of Candidate Genes Associated with Core Agronomic Traits by Four Analysis Models

The elite germplasm resources are key to the beautiful appearance and pleasant flavor of Biluochun tea. We collected and measured the agronomic traits of 95 tea plants to reveal the trait diversity and breeding value of Biluochun tea plant populations. The results revealed that the agronomic traits of Biluochun tea plant populations were diverse and had high breeding value. Additionally, we resequenced these tea plant populations to reveal genetic diversity, population structure, and selection pressure. The Biluochun tea plant populations contained two groups and were least affected by natural selection based on the results of population structure and selection pressure. More importantly, four non-synonymous single nucleotide polymorphisms (nsSNPs) and candidate genes associated with (-)-gallocatechin gallate (GCG), (-)-gallocatechin (GC), and caffeine (CAF) were detected using at least two GWAS models. The results will promote the development and application of molecular markers and the utilization of elite germplasm from Biluochun populations.

Evolutionary Analyses of Gene Expression Divergence in Panicum hallii: Exploring Constitutive and Plastic Responses Using Reciprocal Transplants

The evolution of gene expression is thought to be an important mechanism of local adaptation and ecological speciation. Gene expression divergence occurs through the evolution of cis- polymorphisms and through more widespread effects driven by trans-regulatory factors. Here, we explore expression and sequence divergence in a large sample of Panicum hallii accessions encompassing the species range using a reciprocal transplantation experiment. We observed widespread genotype and transplant site drivers of expression divergence, with a limited number of genes exhibiting genotype-by-site interactions. We used a modified FST-QST outlier approach (QPC analysis) to detect local adaptation. We identified 514 genes with constitutive expression divergence above and beyond the levels expected under neutral processes. However, no plastic expression responses met our multiple testing correction as QPC outliers. Constitutive QPC outlier genes were involved in a number of developmental processes and responses to abiotic environments. Leveraging earlier expression quantitative trait loci results, we found a strong enrichment of expression divergence, including for QPC outliers, in genes previously identified with cis and cis-environment interactions but found no patterns related to trans-factors. Population genetic analyses detected elevated sequence divergence of promoters and coding sequence of constitutive expression outliers but little evidence for positive selection on these proteins. Our results are consistent with a hypothesis of cis-regulatory divergence as a primary driver of expression divergence in P. hallii.

A single amino acid change led to structural and functional differentiation of PvHd1 to control flowering in switchgrass

Switchgrass, a forage and bioenergy crop, occurs as two main ecotypes with different but overlapping ranges of adaptation. The two ecotypes differ in a range of characteristics, including flowering time. Flowering time determines the duration of vegetative development and therefore biomass accumulation, a key trait in bioenergy crops. No causal variants for flowering time differences between switchgrass ecotypes have, as yet, been identified. In this study, we mapped a robust flowering time quantitative trait locus (QTL) on chromosome 4K in a biparental F2 population and characterized the flowering-associated transcription factor gene PvHd1, an ortholog of CONSTANS in Arabidopsis and Heading date 1 in rice, as the underlying causal gene. Protein modeling predicted that a serine to glycine substitution at position 35 (p.S35G) in B-Box domain 1 greatly altered the global structure of the PvHd1 protein. The predicted variation in protein compactness was supported in vitro by a 4 °C shift in denaturation temperature. Overexpressing the PvHd1-p.35S allele in a late-flowering CONSTANS-null Arabidopsis mutant rescued earlier flowering, whereas PvHd1-p.35G had a reduced ability to promote flowering, demonstrating that the structural variation led to functional divergence. Our findings provide us with a tool to manipulate the timing of floral transition in switchgrass cultivars and, potentially, expand their cultivation range.

Plant genome resequencing and population genomics: Current status and future prospects

Advances in DNA sequencing technology have sparked a genomics revolution, driving breakthroughs in plant genetics and crop breeding. Recently, the focus has shifted from cataloging genetic diversity in plants to exploring their functional significance and delivering beneficial alleles for crop improvement. This transformation has been facilitated by the increasing adoption of whole-genome resequencing. In this review, we summarize the current progress of population-based genome resequencing studies and how these studies affect crop breeding. A total of 187 land plants from 163 countries have been resequenced, comprising 54 413 accessions. As part of resequencing efforts 367 traits have been surveyed and 86 genome-wide association studies have been conducted. Economically important crops, particularly cereals, vegetables, and legumes, have dominated the resequencing efforts, leaving a gap in 49 orders, including Lycopodiales, Liliales, Acorales, Austrobaileyales, and Commelinales. The resequenced germplasm is distributed across diverse geographic locations, providing a global perspective on plant genomics. We highlight genes that have been selected during domestication, or associated with agronomic traits, and form a repository of candidate genes for future research and application. Despite the opportunities for cross-species comparative genomics, many population genomic datasets are not accessible, impeding secondary analyses. We call for a more open and collaborative approach to population genomics that promotes data sharing and encourages contribution-based credit policy. The number of plant genome resequencing studies will continue to rise with the decreasing DNA sequencing costs, coupled with advances in analysis and computational technologies. This expansion, in terms of both scale and quality, holds promise for deeper insights into plant trait genetics and breeding design.

Adaptive evolution to the natural and anthropogenic environment in a global invasive crop pest, the cotton bollworm

The cotton bollworm, Helicoverpa armigera, is set to become the most economically devastating crop pest in the world, threatening food security and biosafety as its range expands across the globe. Key to understanding the eco-evolutionary dynamics of H. armigera, and thus its management, is an understanding of population connectivity and the adaptations that allow the pest to establish in unique environments. We assembled a chromosome-scale reference genome and re-sequenced 503 individuals spanning the species range to delineate global patterns of connectivity, uncovering a previously cryptic population structure. Using a genome-wide association study (GWAS) and cell line expression of major effect loci, we show that adaptive changes in a temperature- and light-sensitive developmental pathway enable facultative diapause and that adaptation of trehalose synthesis and transport underlies cold tolerance in extreme environments. Incorporating extensive pesticide resistance monitoring, we also characterize a suite of novel pesticide and Bt resistance alleles under selection in East China. These findings offer avenues for more effective management strategies and provide insight into how insects adapt to variable climatic conditions and newly colonized environments.

A temporal analysis and response to nitrate availability of 3D root system architecture in diverse pennycress (Thlaspi arvense L.) accessions

Introduction

Roots have a central role in plant resource capture and are the interface between the plant and the soil that affect multiple ecosystem processes. Field pennycress (Thlaspi arvense L.) is a diploid annual cover crop species that has potential utility for reducing soil erosion and nutrient losses; and has rich seeds (30-35% oil) amenable to biofuel production and as a protein animal feed. The objective of this research was to (1) precisely characterize root system architecture and development, (2) understand plastic responses of pennycress roots to nitrate nutrition, (3) and determine genotypic variance available in root development and nitrate plasticity.

Methods

Using a root imaging and analysis pipeline, the 4D architecture of the pennycress root system was characterized under four nitrate regimes, ranging from zero to high nitrate concentrations. These measurements were taken at four time points (days 5, 9, 13, and 17 after sowing).

Results

Significant nitrate condition response and genotype interactions were identified for many root traits, with the greatest impact observed on lateral root traits. In trace nitrate conditions, a greater lateral root count, length, density, and a steeper lateral root angle was observed compared to high nitrate conditions. Additionally, genotype-by-nitrate condition interaction was observed for root width, width:depth ratio, mean lateral root length, and lateral root density.

Discussion

These findings illustrate root trait variance among pennycress accessions. These traits could serve as targets for breeding programs aimed at developing improved cover crops that are responsive to nitrate, leading to enhanced productivity, resilience, and ecosystem service.

A graph-based genome and pan-genome variation of the model plant Setaria

Setaria italica (foxtail millet), a founder crop of East Asian agriculture, is a model plant for C4 photosynthesis and developing approaches to adaptive breeding across multiple climates. Here we established the Setaria pan-genome by assembling 110 representative genomes from a worldwide collection. The pan-genome is composed of 73,528 gene families, of which 23.8%, 42.9%, 29.4% and 3.9% are core, soft core, dispensable and private genes, respectively; 202,884 nonredundant structural variants were also detected. The characterization of pan-genomic variants suggests their importance during foxtail millet domestication and improvement, as exemplified by the identification of the yield gene SiGW3, where a 366-bp presence/absence promoter variant accompanies gene expression variation. We developed a graph-based genome and performed large-scale genetic studies for 68 traits across 13 environments, identifying potential genes for millet improvement at different geographic sites. These can be used in marker-assisted breeding, genomic selection and genome editing to accelerate crop improvement under different climatic conditions.

Genome of tetraploid sour cherry (Prunus cerasus L.) 'Montmorency' identifies three distinct ancestral Prunus genomes

Sour cherry (Prunus cerasus L.) is a valuable fruit crop in the Rosaceae family and a hybrid between progenitors closely related to extant Prunus fruticosa (ground cherry) and Prunus avium (sweet cherry). Here we report a chromosome-scale genome assembly for sour cherry cultivar Montmorency, the predominant cultivar grown in the USA. We also generated a draft assembly of P. fruticosa to use alongside a published P. avium sequence for syntelog-based subgenome assignments for 'Montmorency' and provide compelling evidence P. fruticosa is also an allotetraploid. Using hierarchal k-mer clustering and phylogenomics, we show 'Montmorency' is trigenomic, containing two distinct subgenomes inherited from a P. fruticosa-like ancestor (A and A') and two copies of the same subgenome inherited from a P. avium-like ancestor (BB). The genome composition of 'Montmorency' is AA'BB and little-to-no recombination has occurred between progenitor subgenomes (A/A' and B). In Prunus, two known classes of genes are important to breeding strategies: the self-incompatibility loci (S-alleles), which determine compatible crosses, successful fertilization, and fruit set, and the Dormancy Associated MADS-box genes (DAMs), which strongly affect dormancy transitions and flowering time. The S-alleles and DAMs in 'Montmorency' and P. fruticosa were manually annotated and support subgenome assignments. Lastly, the hybridization event 'Montmorency' is descended from was estimated to have occurred less than 1.61 million years ago, making sour cherry a relatively recent allotetraploid. The 'Montmorency' genome highlights the evolutionary complexity of the genus Prunus and will inform future breeding strategies for sour cherry, comparative genomics in the Rosaceae, and questions regarding neopolyploidy.

Transposon signatures of allopolyploid genome evolution

Hybridization brings together chromosome sets from two or more distinct progenitor species. Genome duplication associated with hybridization, or allopolyploidy, allows these chromosome sets to persist as distinct subgenomes during subsequent meioses. Here, we present a general method for identifying the subgenomes of a polyploid based on shared ancestry as revealed by the genomic distribution of repetitive elements that were active in the progenitors. This subgenome-enriched transposable element signal is intrinsic to the polyploid, allowing broader applicability than other approaches that depend on the availability of sequenced diploid relatives. We develop the statistical basis of the method, demonstrate its applicability in the well-studied cases of tobacco, cotton, and Brassica napus, and apply it to several cases: allotetraploid cyprinids, allohexaploid false flax, and allooctoploid strawberry. These analyses provide insight into the origins of these polyploids, revise the subgenome identities of strawberry, and provide perspective on subgenome dominance in higher polyploids.

Molecular evolution, diversification, and expression assessment of MADS gene family in Setaria italica, Setaria viridis, and Panicum virgatum

KEY MESSAGE

This paper sheds light on the evolution and expression patterns of MADS genes in Setaria and Panicum virgatum. SiMADS51 and SiMADS64 maybe involved in the ABA-dependent pathway of drought response. The MADS gene family is a key regulatory factor family that controls growth, reproduction, and response to abiotic stress in plants. However, the molecular evolution of this family is rarely reported. Here, a total of 265 MADS genes were identified in Setaria italica (foxtail millet), Setaria viridis (green millet), and Panicum virgatum (switchgrass) and analyzed by bioinformatics, including physicochemical characteristics, subcellular localization, chromosomal position and duplicate, motif distribution, genetic structure, genetic evolvement, and expression patterns. Phylogenetic analysis was used to categorize these genes into M and MIKC types. The distribution of motifs and gene structure were similar for the corresponding types. According to a collinearity study, the MADS genes have been mostly conserved during evolution. The principal cause of their expansion is segmental duplication. However, the MADS gene family tends to shrink in foxtail millet, green millet, and switchgrass. The MADS genes were subjected to purifying selection, but several positive selection sites were also identified in three species. And most of the promoters of MADS genes contain cis-elements related to stress and hormonal response. RNA-seq and quantitative Real-time PCR (qRT-PCR) analysis also were examined. SiMADS genes expression levels are considerably changed in reaction to various treatments, following qRT-PCR analysis. This sheds fresh light on the evolution and expansion of the MADS family in foxtail millet, green millet, and switchgrass, and lays the foundation for further research on their functions.

Mapping quantitative trait loci for biomass yield and yield-related traits in lowland switchgrass (Panicum virgatum L.) multiple populations

Switchgrass can be used as an alternative source for bioenergy production. Many breeding programs focus on the genetic improvement of switchgrass for increasing biomass yield. Quantitative trait loci (QTL) mapping can help to discover marker-trait associations and accelerate the breeding process through marker-assisted selection. To identify significant QTL, this study mapped 7 hybrid populations and one combined of 2 hybrid populations (30-96 F1s) derived from Alamo and Kanlow genotypes. The populations were evaluated for biomass yield, plant height, and crown size in a simulated-sward plot with 2 replications at 2 locations in Tennessee from 2019 to 2021. The populations showed significant genetic variation for the evaluated traits and exhibited transgressive segregation. The 17,251 single nucleotide polymorphisms (SNPs) generated through genotyping-by-sequencing (GBS) were used to construct a linkage map using a fast algorithm for multiple outbred families. The linkage map spanned 1,941 cM with an average interval of 0.11 cM between SNPs. The QTL analysis was performed on evaluated traits for each and across environments (year and location) that identified 5 QTL for biomass yield (logarithm of the odds, LOD 3.12-4.34), 4 QTL for plant height (LOD 3.01-5.64), and 7 QTL for crown size (LOD 3.0-4.46) (P ≤ 0.05). The major QTL for biomass yield, plant height, and crown size resided on chromosomes 8N, 6N, and 8K explained phenotypic variations of 5.6, 5.1, and 6.6%, respectively. SNPs linked to QTL could be incorporated into marker-assisted breeding to maximize the selection gain in switchgrass breeding.

DNA methylation signatures of duplicate gene evolution in angiosperms

Gene duplication is a source of evolutionary novelty. DNA methylation may play a role in the evolution of duplicate genes (paralogs) through its association with gene expression. While this relationship has been examined to varying extents in a few individual species, the generalizability of these results at either a broad phylogenetic scale with species of differing duplication histories or across a population remains unknown. We applied a comparative epigenomics approach to 43 angiosperm species across the phylogeny and a population of 928 Arabidopsis (Arabidopsis thaliana) accessions, examining the association of DNA methylation with paralog evolution. Genic DNA methylation was differentially associated with duplication type, the age of duplication, sequence evolution, and gene expression. Whole genome duplicates were typically enriched for CG-only gene-body methylated or unmethylated genes, while single-gene duplications were typically enriched for non-CG methylated or unmethylated genes. Non-CG methylation, in particular, was characteristic of more recent single-gene duplicates. Core angiosperm gene families differentiated into those which preferentially retain paralogs and 'duplication-resistant' families, which convergently reverted to singletons following duplication. Duplication-resistant families that still have paralogous copies were, uncharacteristically for core angiosperm genes, enriched for non-CG methylation. Non-CG methylated paralogs had higher rates of sequence evolution, higher frequency of presence-absence variation, and more limited expression. This suggests that silencing by non-CG methylation may be important to maintaining dosage following duplication and be a precursor to fractionation. Our results indicate that genic methylation marks differing evolutionary trajectories and fates between paralogous genes and have a role in maintaining dosage following duplication.

Genetic determinants of switchgrass-root-associated microbiota in field sites spanning its natural range

A fundamental goal in plant microbiome research is to determine the relative impacts of host and environmental effects on root microbiota composition, particularly how host genotype impacts bacterial community composition. Most studies characterizing the effect of plant genotype on root microbiota undersample host genetic diversity and grow plants outside of their native ranges, making the associations between host and microbes difficult to interpret. Here, we characterized the root microbiota of a large diversity panel of switchgrass, a North American native C4 bioenergy crop, in three field locations spanning its native range. Our data, composed of 1,961 samples, suggest that field location is the primary determinant of microbiome composition; however, substantial heritable variation is widespread across bacterial taxa, especially those in the Sphingomonadaceae family. Despite diverse compositions, relatively few highly prevalent taxa make up the majority of the switchgrass root microbiota, a large fraction of which is shared across sites. Local genotypes preferentially recruit/filter for local microbes, supporting the idea of affinity between local plants and their microbiota. Using genome-wide association, we identified loci impacting the abundance of >400 microbial strains and found an enrichment of genes involved in immune responses, signaling pathways, and secondary metabolism. We found loci associated with over half of the core microbiota (i.e., microbes in >80% of samples), regardless of field location. Finally, we show a genetic relationship between a basal plant immunity pathway and relative abundances of root microbiota. This study brings us closer to harnessing and manipulating beneficial microbial associations via host genetics.

Tracing 100 million years of grass genome evolutionary plasticity

Grasses derive from a family of monocotyledonous plants that includes crops of major economic importance such as wheat, rice, sorghum and barley, sharing a common ancestor some 100 million years ago. The genomic attributes of plant adaptation remain obscure and the consequences of recurrent whole genome duplications (WGD) or polyploidization events, a major force in plant evolution, remain largely speculative. We conducted a comparative analysis of omics data from ten grass species to unveil structural (inversions, fusions, fissions, duplications, substitutions) and regulatory (expression and methylation) basis of genome plasticity, as possible attributes of plant long lasting evolution and adaptation. The present study demonstrates that diverged polyploid lineages sharing a common WGD event often present the same patterns of structural changes and evolutionary dynamics, but these patterns are difficult to generalize across independent WGD events as a result of non-WGD factors such as selection and domestication of crops. Polyploidy is unequivocally linked to the evolutionary success of grasses during the past 100 million years, although it remains difficult to attribute this success to particular genomic consequences of polyploidization, suggesting that polyploids harness the potential of genome duplication, at least partially, in lineage-specific ways. Overall, the present study clearly demonstrates that post-polyploidization reprogramming is more complex than traditionally reported in investigating single species and calls for a critical and comprehensive comparison across independently polyploidized lineages.

Genomic prediction of switchgrass winter survivorship across diverse lowland populations

In the North-Central United States, lowland ecotype switchgrass can increase yield by up to 50% compared with locally adapted but early flowering cultivars. However, lowland ecotypes are not winter tolerant. The mechanism for winter damage is unknown but previously has been associated with late flowering time. This study investigated heading date (measured for two years) and winter survivorship (measured for three years) in a multi-generation population generated from two winter-hardy lowland individuals and diverse southern lowland populations. Sequencing data (311,776 markers) from 1,306 individuals were used to evaluate genome-wide trait prediction through cross-validation and progeny prediction (n = 52). Genetic variance for heading date and winter survivorship was additive with high narrow-sense heritability (0.64 and 0.71, respectively) and reliability (0.68 and 0.76, respectively). The initial negative correlation between winter survivorship and heading date degraded across generations (F1r = -0.43, pseudo-F2r = -0.28, pseudo-F2 progeny r = -0.15). Within-family predictive ability was moderately high for heading date and winter survivorship (0.53 and 0.52, respectively). A multi-trait model did not improve predictive ability for either trait. Progeny predictive ability was 0.71 for winter survivorship and 0.53 for heading date. These results suggest that lowland ecotype populations can obtain sufficient survival rates in the northern United States with two or three cycles of effective selection. Despite accurate genomic prediction, naturally occurring winter mortality successfully isolated winter tolerant genotypes and appears to be an efficient method to develop high-yielding, cold-tolerant switchgrass cultivars.

Understanding local plant extinctions before it is too late: bridging evolutionary genomics with global ecology

Understanding evolutionary genomic and population processes within a species range is key to anticipating the extinction of plant species before it is too late. However, most models of biodiversity risk under global change do not account for the genetic variation and local adaptation of different populations. Population diversity is critical to understanding extinction because different populations may be more or less susceptible to global change and, if lost, would reduce the total diversity within a species. Two new modeling frameworks advance our understanding of extinction from a population and evolutionary angle: Rapid climate change-driven disruptions in population adaptation are predicted from associations between genomes and local climates. Furthermore, losses of population diversity from global land-use transformations are estimated by scaling relationships of species' genomic diversity with habitat area. Overall, these global eco-evolutionary methods advance the predictability - and possibly the preventability - of the ongoing extinction of plant species.

Biased mutations and gene losses underlying diploidization of the tetraploid broomcorn millet genome

Broomcorn millet (Panicum miliaceum L.) is one of the earliest domesticated crops, and is a valuable resource to secure food diversity and combat drought stresses under the global warming scenario. However, due to the absence of extant diploid progenitors, the polyploidy genome of broomcorn millet remains poorly understood. Here, we report the chromosome-scale genome assembly of broomcorn millet. We divided the broomcorn millet genome into two subgenomes using the genome sequence of Panicum hallii, a diploid relative of broomcorn millet. Our analyses revealed that the two subgenomes diverged at ~4.8 million years ago (Mya), while the allotetraploidization of broomcorn millet may have occurred about ~0.48 Mya, suggesting that broomcorn millet is a relatively recent allotetraploid. Comparative analyses showed that subgenome B was larger than subgenome A in size, which was caused by the biased accumulation of long terminal repeat retrotransposons in the progenitor of subgenome B before polyploidization. Notably, the accumulation of biased mutations in the transposable element-rich subgenome B led to more gene losses. Although no significant dominance of either subgenome was observed in the expression profiles of broomcorn millet, we found the minimally expressed genes in P. hallii tended to be lost during diploidization of broomcorn millet. These results suggest that broomcorn millet is at the early stage of diploidization and that mutations likely occurred more on genes that were marked with lower expression levels.

Gene-by-environment interactions in plants: Molecular mechanisms, environmental drivers, and adaptive plasticity

Plants demonstrate a broad range of responses to environmental shifts. One of the most remarkable responses is plasticity, which is the ability of a single plant genotype to produce different phenotypes in response to environmental stimuli. As with all traits, the ability of plasticity to evolve depends on the presence of underlying genetic diversity within a population. A common approach for evaluating the role of genetic variation in driving differences in plasticity has been to study genotype-by-environment interactions (G × E). G × E occurs when genotypes produce different phenotypic trait values in response to different environments. In this review, we highlight progress and promising methods for identifying the key environmental and genetic drivers of G × E. Specifically, methodological advances in using algorithmic and multivariate approaches to understand key environmental drivers combined with new genomic innovations can greatly increase our understanding about molecular responses to environmental stimuli. These developing approaches can be applied to proliferating common garden networks that capture broad natural environmental gradients to unravel the underlying mechanisms of G × E. An increased understanding of G × E can be used to enhance the resilience and productivity of agronomic systems.

Genotype-environment associations to reveal the molecular basis of environmental adaptation

A fundamental goal in plant biology is to identify and understand the variation underlying plants' adaptation to their environment. Climate change has given new urgency to this goal, as society aims to accelerate adaptation of ecologically important plant species, endangered plant species, and crops to hotter, less predictable climates. In the pre-genomic era, identifying adaptive alleles was painstaking work, leveraging genetics, molecular biology, physiology, and ecology. Now, the rise of genomics and new computational approaches may facilitate this research. Genotype-environment associations (GEAs) use statistical associations between allele frequency and environment of origin to test the hypothesis that allelic variation at a given gene is adapted to local environments. Researchers may scan the genome for GEAs to generate hypotheses on adaptive genetic variants (environmental genome-wide association studies). Despite the rapid adoption of these methods, many important questions remain about the interpretation of GEA findings, which arise from fundamental unanswered questions on the genetic architecture of adaptation and limitations inherent to association-based analyses. We outline strategies to ground GEAs in the underlying hypotheses of genetic architecture and better test GEA-generated hypotheses using genetics and ecophysiology. We provide recommendations for new users who seek to learn about the molecular basis of adaptation. When combined with a rigorous hypothesis testing framework, GEAs may facilitate our understanding of the molecular basis of climate adaptation for plant improvement.

A Transferable Machine Learning Framework for Predicting Transcriptional Responses of Genes Across Species

Leveraging existing resources in studied species to predict gene functions has the potential to rapidly expand understanding of annotated genes in other, less well-studied, species with assembled genomes. However, orthology is not a reliable predictor for the transcriptional responses of genes to stress. Machine learning methods can quantitatively estimate expression patterns and gene functions using known annotations and collections of features describing each gene. In this chapter, we describe a supervised machine learning framework to predict stress-responsive genes across species using only features derived from nucleotide sequences, using the example of cold stress-responsive genes in different Panicoid grass species.

Legacies of precipitation influence primary production in Panicum virgatum

Precipitation is a key driver of primary production worldwide, but primary production does not always track year-to-year variation in precipitation linearly. Instead, plant responses to changes in precipitation may exhibit time lags, or legacies of past precipitation. Legacies can be driven by multiple mechanisms, including persistent changes in plant physiological and morphological traits and changes to the physical environment, such as plant access to soil water. We used three precipitation manipulation experiments in central Texas, USA to evaluate the magnitude, duration, and potential mechanisms driving precipitation legacies on aboveground primary production of the perennial C₄ grass, Panicum virgatum. Specifically, we performed a rainout shelter study, where eight genotypes grew under different precipitation regimes; a transplant study, where plants that had previously grown in a rainout shelter under different precipitation regimes were moved to a common environment; and a mesocosm study, where the effect of swapping precipitation regime was examined with a single genotype. Across these experiments, plants previously grown under wet conditions generally performed better than expected when exposed to drought. Panicum virgatum exhibited stronger productivity legacies of past wet years on current-year responses to drought than of past dry years on current-year responses to wet conditions. Additionally, previous year tiller counts, a proxy for meristem availability, were important in determining legacy effects on aboveground production. As climate changes and precipitation extremes-both dry and wet-become more common, these results suggest that populations of P. virgatum may become less resilient.

Population genomics and subgenome evolution of the allotetraploid frog Xenopus laevis in southern Africa

Allotetraploid genomes have two distinct genomic components called subgenomes that are derived from separate diploid ancestral species. Many genomic characteristics such as gene function, expression, recombination, and transposable element mobility may differ significantly between subgenomes. To explore the possibility that subgenome population structure and gene flow may differ as well, we examined genetic variation in an allotetraploid frog-the African clawed frog (Xenopus laevis)-over the dynamic and varied habitat of its native range in southern Africa. Using reduced representation genome sequences from 91 samples from 12 localities, we found no strong evidence that population structure and gene flow differed substantially by subgenome. We then compared patterns of population structure in the nuclear genome to the mitochondrial genome using Sanger sequences from 455 samples from 183 localities. Our results provide further resolution to the geographic distribution of mitochondrial and nuclear diversity in this species and illustrate that population structure in both genomes corresponds roughly with variation in seasonal rainfall and with the topography of southern Africa.

On the hidden temporal dynamics of plant adaptation

Adaptation to a wide range of environments is a major driver of plant diversity. It is now possible to catalog millions of potential adaptive genomic differences segregating between environments within a plant species in a single experiment. Understanding which of these changes contributes to adaptive phenotypic divergence between plant populations is a major goal of evolutionary biologists and crop breeders. In this review, we briefly highlight the approaches frequently used to understand the genetic basis of adaptive phenotypes in plants, and we discuss some of the limitations of these methods. We propose that direct observation of the process of adaptation using multigenerational studies and whole genome sequencing is a crucial missing component of recent studies of plant adaptation because it complements several shortcomings of sampling-based techniques.

The leaf economics spectrum of triploid and tetraploid C4 grass Miscanthus x giganteus

The leaf economics spectrum (LES) describes multivariate correlations in leaf structural, physiological and chemical traits, originally based on diverse C₃ species grown under natural ecosystems. However, the specific contribution of C₄ species to the global LES is studied less widely. C₄ species have a CO₂ concentrating mechanism which drives high rates of photosynthesis and improves resource use efficiency, thus potentially pushing them towards the edge of the LES. Here, we measured foliage morphology, structure, photosynthesis, and nutrient content for hundreds of genotypes of the C₄ grass Miscanthus× giganteus grown in two common gardens over two seasons. We show substantial trait variations across M.× giganteus genotypes and robust genotypic trait relationships. Compared to the global LES, M.× giganteus genotypes had higher photosynthetic rates, lower stomatal conductance, and less nitrogen content, indicating greater water and photosynthetic nitrogen use efficiency in the C₄ species. Additionally, tetraploid genotypes produced thicker leaves with greater leaf mass per area and lower leaf density than triploid genotypes. By expanding the LES relationships across C₃ species to include C₄ crops, these findings highlight that M.× giganteus occupies the boundary of the global LES and suggest the potential for ploidy to alter LES traits.

Local adaptation: Causal agents of selection and adaptive trait divergence

Divergent selection across the landscape can favor the evolution of local adaptation in populations experiencing contrasting conditions. Local adaptation is widely observed in a diversity of taxa, yet we have a surprisingly limited understanding of the mechanisms that give rise to it. For instance, few have experimentally confirmed the biotic and abiotic variables that promote local adaptation, and fewer yet have identified the phenotypic targets of selection that mediate local adaptation. Here, we highlight critical gaps in our understanding of the process of local adaptation and discuss insights emerging from in-depth investigations of the agents of selection that drive local adaptation, the phenotypes they target, and the genetic basis of these phenotypes. We review historical and contemporary methods for assessing local adaptation, explore whether local adaptation manifests differently across life history, and evaluate constraints on local adaptation.

The Boechera model system for evolutionary ecology

Model systems in biology expand the research capacity of individuals and the community. Closely related to Arabidopsis, the genus Boechera has emerged as an important ecological model owing to the ability to integrate across molecular, functional, and eco-evolutionary approaches. Boechera species are broadly distributed in relatively undisturbed habitats predominantly in western North America and provide one of the few experimental systems for identification of ecologically important genes through genome-wide association studies and investigations of selection with plants in their native habitats. The ecologically, evolutionarily, and agriculturally important trait of apomixis (asexual reproduction via seeds) is common in the genus, and field experiments suggest that abiotic and biotic environments shape the evolution of sex. To date, population genetic studies have focused on the widespread species B. stricta, detailing population divergence and demographic history. Molecular and ecological studies show that balancing selection maintains genetic variation in ~10% of the genome, and ecological trade-offs contribute to complex trait variation for herbivore resistance, flowering phenology, and drought tolerance. Microbiome analyses have shown that host genotypes influence leaf and root microbiome composition, and the soil microbiome influences flowering phenology and natural selection. Furthermore, Boechera offers numerous opportunities for investigating biological responses to global change. In B. stricta, climate change has induced a shift of >2 weeks in the timing of first flowering since the 1970s, altered patterns of natural selection, generated maladaptation in previously locally-adapted populations, and disrupted life history trade-offs. Here we review resources and results for this eco-evolutionary model system and discuss future research directions.

Cold acclimation threshold induction temperatures of switchgrass ecotypes grown under a long and short photoperiod

Plants can cold acclimate to enhance their freezing tolerance by sensing declining temperature and photoperiod cues. However, the factors influencing genotypic variation in the induction of cold acclimation are poorly understood among perennial grasses. We hypothesized that the more northern upland switchgrass (Panicum virgatum L.) ecotype develops a higher degree of freezing tolerance by initiating cold acclimation at higher temperatures as compared with the coastal and southern lowland ecotypes. First, we determined the optimal method for assessing freezing tolerance and the length of exposure to 8/4°C required to induce the maximum level of freezing tolerance in the most northern upland and most southern lowland genotypes. We characterized the maximum freezing tolerance of eight uplands, three coastal and five lowland genotypes grown for 21 days at 8/4°C and a 10 or 16 h photoperiod. Next, we identified the temperature required to induce cold acclimation by exposing the 16 genotypes for 7 days at 20-6°C constant temperatures under a 10 or 16 h photoperiod. Cold acclimation initiated at temperatures 5 and 7°C higher in upland than in coastal and lowland genotypes. Among upland genotypes the shorter photoperiod induced cold acclimation at a 1°C higher temperature. Genotypes originating from a more northern latitude initiate cold acclimation at higher temperatures and develop higher maximum freezing tolerances. An earlier response to declining temperatures may provide the upland ecotype with additional time to prepare for winter and provide an advantage when plants are subjected to the rapid changes in fall temperature associated with injurious frosts.

Hybridisation and chloroplast capture between distinct Themeda triandra lineages in Australia

Ecotypes are distinct populations within a species that are adapted to specific environmental conditions. Understanding how these ecotypes become established, and how they interact when reunited, is fundamental to elucidating how ecological adaptations are maintained. This study focuses on Themeda triandra, a dominant grassland species across Asia, Africa and Australia. It is the most widespread plant in Australia, where it has distinct ecotypes that are usually restricted to either wetter and cooler coastal regions or the drier and hotter interior. We generate a reference genome for T. triandra and use whole genome sequencing for over 80 Themeda accessions to reconstruct the evolutionary history of T. triandra and related taxa. Organelle phylogenies confirm that Australia was colonized by T. triandra twice, with the division between ecotypes predating their arrival in Australia. The nuclear genome provides evidence of differences in the dominant ploidal level and gene-flow among the ecotypes. In northern Queensland there appears to be a hybrid zone between ecotypes with admixed nuclear genomes and shared chloroplast haplotypes. Conversely, in the cracking claypans of Western Australia, there is cytonuclear discordance with individuals possessing the coastal chloroplast and interior clade nuclear genome. This chloroplast capture is potentially a result of adaptive introgression, with selection detected in the rpoC2 gene which is associated with water use efficiency. The reason that T. triandra is the most widespread plant in Australia appears to be a result of distinct ecotypic genetic variation and genome duplication, with the importance of each depending on the geographic scale considered.

Comparative transcriptomics and metabolomics reveal specialized metabolite drought stress responses in switchgrass (Panicum virgatum)

Switchgrass (Panicum virgatum) is a bioenergy model crop valued for its energy efficiency and drought tolerance. The related monocot species rice (Oryza sativa) and maize (Zea mays) deploy species-specific, specialized metabolites as core stress defenses. By contrast, specialized chemical defenses in switchgrass are largely unknown. To investigate specialized metabolic drought responses in switchgrass, we integrated tissue-specific transcriptome and metabolite analyses of the genotypes Alamo and Cave-in-Rock that feature different drought tolerance. The more drought-susceptible Cave-in-Rock featured an earlier onset of transcriptomic changes and significantly more differentially expressed genes in response to drought compared to Alamo. Specialized pathways showed moderate differential expression compared to pronounced transcriptomic alterations in carbohydrate and amino acid metabolism. However, diterpenoid-biosynthetic genes showed drought-inducible expression in Alamo roots, contrasting largely unaltered triterpenoid and phenylpropanoid pathways. Metabolomic analyses identified common and genotype-specific flavonoids and terpenoids. Consistent with transcriptomic alterations, several root diterpenoids showed significant drought-induced accumulation, whereas triterpenoid abundance remained predominantly unchanged. Structural analysis verified select drought-responsive diterpenoids as oxygenated furanoditerpenoids. Drought-dependent transcriptome and metabolite profiles provide the foundation to understand the molecular mechanisms underlying switchgrass drought responses. Accumulation of specialized root diterpenoids and corresponding pathway transcripts supports a role in drought stress tolerance.

Identification of the genes at S and Z reveals the molecular basis and evolution of grass self-incompatibility

Self-incompatibility (SI) is a feature of many flowering plants, whereby self-pollen is recognized and rejected by the stigma. In grasses (Poaceae), the genes controlling this phenomenon have not been fully elucidated. Grasses have a unique two-locus system, in which two independent genetic loci (S and Z) control self-recognition. S and Z are thought to have arisen from an ancient duplication, common to all grasses. With new chromosome-scale genome data, we examined the genes present at S- and Z-loci, firstly in ryegrass (Lolium perenne), and subsequently in ~20 other grass species. We found that two DUF247 genes and a short unstructured protein (SP/ZP) were present at both S- and Z- in all SI species, while in self-compatible species these genes were often lost or mutated. Expression data suggested that DUF247 genes acted as the male components and SP/ZP were the female components. Consistent with their role in distinguishing self- from non-self, all genes were hypervariable, although key secondary structure features were conserved, including the predicted N-terminal cleavage site of SP/ZP. The evolutionary history of these genes was probed, revealing that specificity groups at the Z-locus arose before the advent of various grass subfamilies/species, while specificity groups at the S-locus arose after the split of Panicoideae, Chloridoideae, Oryzoideae and Pooideae. Finally, we propose a model explaining how the proteins encoded at the S and Z loci might function to specify self-incompatibility.

Temporal regulation of cold transcriptional response in switchgrass

Switchgrass low-land ecotypes have significantly higher biomass but lower cold tolerance compared to up-land ecotypes. Understanding the molecular mechanisms underlying cold response, including the ones at transcriptional level, can contribute to improving tolerance of high-yield switchgrass under chilling and freezing environmental conditions. Here, by analyzing an existing switchgrass transcriptome dataset, the temporal cis-regulatory basis of switchgrass transcriptional response to cold is dissected computationally. We found that the number of cold-responsive genes and enriched Gene Ontology terms increased as duration of cold treatment increased from 30 min to 24 hours, suggesting an amplified response/cascading effect in cold-responsive gene expression. To identify genomic sequences likely important for regulating cold response, machine learning models predictive of cold response were established using k-mer sequences enriched in the genic and flanking regions of cold-responsive genes but not non-responsive genes. These k-mers, referred to as putative cis-regulatory elements (pCREs) are likely regulatory sequences of cold response in switchgrass. There are in total 655 pCREs where 54 are important in all cold treatment time points. Consistent with this, eight of 35 known cold-responsive CREs were similar to top-ranked pCREs in the models and only these eight were important for predicting temporal cold response. More importantly, most of the top-ranked pCREs were novel sequences in cold regulation. Our findings suggest additional sequence elements important for cold-responsive regulation previously not known that warrant further studies.

Identification of microRNAs responsive to arbuscular mycorrhizal fungi in Panicum virgatum (switchgrass)

BACKGROUND

MicroRNAs (miRNAs) are important post-transcriptional regulators involved in the control of a range of processes, including symbiotic interactions in plants. MiRNA involvement in arbuscular mycorrhizae (AM) symbiosis has been mainly studied in model species, and our study is the first to analyze global miRNA expression in the roots of AM colonized switchgrass (Panicum virgatum), an emerging biofuel feedstock. AM symbiosis helps plants gain mineral nutrition from the soil and may enhance switchgrass biomass production on marginal lands. Our goals were to identify miRNAs and their corresponding target genes that are controlling AM symbiosis in switchgrass.

RESULTS

Through genome-wide analysis of next-generation miRNA sequencing reads generated from switchgrass roots, we identified 122 mature miRNAs, including 28 novel miRNAs. By comparing miRNA expression profiles of AM-inoculated and control switchgrass roots, we identified 15 AM-responsive miRNAs across lowland accession "Alamo", upland accession "Dacotah", and two upland/lowland F₁ hybrids. We used degradome sequencing to identify target genes of the AM-responsive miRNAs revealing targets of miRNAs residing on both K and N subgenomes. Notably, genes involved in copper ion binding were targeted by downregulated miRNAs, while upregulated miRNAs mainly targeted GRAS family transcription factors.

CONCLUSION

Through miRNA analysis and degradome sequencing, we revealed that both upland and lowland switchgrass genotypes as well as upland-lowland hybrids respond to AM by altering miRNA expression. We demonstrated complex GRAS transcription factor regulation by the miR171 family, with some miR171 family members being AM responsive while others remained static. Copper miRNA downregulation was common amongst the genotypes tested and we identified superoxide dismutases and laccases as targets, suggesting that these Cu-miRNAs are likely involved in ROS detoxification and lignin deposition, respectively. Other prominent targets of the Cu miRNAs were blue copper proteins. Overall, the potential effect of AM colonization on lignin deposition pathways in this biofuel crop highlights the importance of considering AM and miRNA in future biofuel crop development strategies.