Transcriptional Regulation of Sorghum Stem Composition: Key Players Identified Through Co-expression Gene Network and Comparative Genomics Analyses

Unveiling the Genomic Symphony: Identification Cultivar-Specific Genes and Enhanced Insights on Sweet Sorghum Genomes Through Comprehensive superTranscriptomic Analysis

Sorghum (Sorghum bicolor (L.) Moench) is a multipurpose crop grown for food, fodder, and bioenergy production. Its cultivated varieties, along with their wild counterparts, contribute to the core genetic pool. Despite the availability of several re-sequenced sorghum genomes, a variable portion of sorghum genomes is not reported during reference genome assembly and annotation. The present analysis used 223 publicly available RNA-seq datasets from seven sweet sorghum cultivars to construct superTranscriptome. This approach yielded 45,864 Representative Transcript Assemblies (RTAs) that showcased intriguing Presence/Absence Variation (PAV) across 15 published sorghum genomes. We found 301 superTranscripts were exclusive to sweet sorghum, including 58 de novo genes encoded core and linker histones, zinc finger domains, glucosyl transferases, cellulose synthase, etc. The superTranscriptome added 2,802 new protein-coding genes to the Sweet Sorghum Reference Genome (SSRG), of which 559 code for different transcription factors (TFs). Our analysis revealed that MULE-like transposases were abundant in the sweet sorghum genome and could play a hidden role in the evolution of sweet sorghum. We observed large deletions in the D locus and terminal deletions in four other NAC encoding loci in the SSRG compared to its wild progenitor (353) suggesting non-functional NAC genes contributed to trait development in sweet sorghum. Moreover, superTranscript-based methods for Differential Exon Usage (DEU) and Differential Gene Expression (DGE) analyses were more accurate than those based on the SSRG. This study demonstrates that the superTranscriptome can enhance our understanding of fundamental sorghum mechanisms, improve genome annotations, and potentially even replace the reference genome.

Bioenergy sorghum nodal root bud development: morphometric, transcriptomic and gene regulatory network analysis

Bioenergy sorghum's large and deep nodal root system and associated microbiome enables uptake of water and nutrients from and deposition of soil organic carbon into soil profiles, key contributors to the crop's resilience and sustainability. The goal of this study was to increase our understanding of bioenergy sorghum nodal root bud development. Sorghum nodal root bud initiation was first observed on the stem node of the 7th phytomer below the shoot apex. Buds were initiated near the upper end of the stem node pulvinus on the side of the stem opposite the tiller bud, then additional buds were added over the next 6-8 days forming a ring of 10-15 nascent nodal root buds around the stem. Later in plant development, a second ring of nodal root buds began forming on the 17th stem node immediately above the first ring of buds. Overall, nodal root bud development can take ~40 days from initiation to onset of nodal root outgrowth. Nodal root buds were initiated in close association with vascular bundles in the rind of the pulvinus. Stem tissue forming nascent nodal root buds expressed sorghum homologs of genes associated with root initiation (WOX4), auxin transport (LAX2, PIN4), meristem activation (NGAL2), and genes involved in cell proliferation. Expression of WOX11 and WOX5, genes involved in root stem niche formation, increased early in nodal root bud development followed by genes encoding PLTs, LBDs (LBD29), LRP1, SMB, RGF1 and root cap LEAs later in development. A nodal root bud gene regulatory network module expressed during nodal root bud initiation predicted connections linking PFA5, SPL9 and WOX4 to genes involved in hormone signaling, meristem activation, and cell proliferation. A network module expressed later in development predicted connections among SOMBRERO, a gene involved in root cap formation, and GATA19, BBM, LBD29 and RITF1/RGF1 signaling. Overall, this study provides a detailed description of bioenergy sorghum nodal root bud development and transcriptome information useful for understanding the regulation of sorghum nodal root bud formation and development.

UV-induced reactive oxygen species and transcriptional control of 3-deoxyanthocyanidin biosynthesis in black sorghum pericarp

Black pericarp sorghum has notable value due to the biosynthesis of 3-deoxyanthocyanidins (3-DOAs), a rare class of bioactive polyphenols valued as antioxidant food additives and as bioactive compounds with cytotoxicity to human cancer cells. A metabolic and transcriptomic study was conducted to ascertain the cellular events leading to the activation of 3-DOA biosynthesis in black sorghum pericarp. Prolonged exposure of pericarp during grain maturation to high-fluence ultraviolet (UV) light resulted in elevated levels of reactive oxygen species (ROS) and the activation of 3-DOA biosynthesis in pericarp tissues. In conjunction with 3-DOA biosynthesis was the transcriptional activation of specific family members of early and late flavonoid biosynthesis pathway genes as well as the downstream activation of defense-related pathways. Promoter analysis of genes highly correlated with 3-DOA biosynthesis in black pericarp were enriched in MYB and HHO5/ARR-B motifs. Light microscopy studies of black pericarp tissues suggest that 3-DOAs are predominantly localized in the epicarp and are associated with the cell wall. A working model of UV-induced 3-DOA biosynthesis in black pericarp is proposed that shares features of plant immunity associated with pathogen attack or mechanical wounding. The present model depicts ROS accumulation, the transcriptional activation of receptor kinases and transcription factors (TFs) including NAC, WRKY, bHLH, AP2, and C2H2 Zinc finger domain. This study identified key biosynthetic and regulatory genes of 3-DOA accumulation in black pericarp and provided a deeper understanding of the gene networks and cellular events controlling this tissue-and genotype-specific trait.

Stem lodging Resistance-1 controls stem strength by positively regulating the biosynthesis of cell wall components in Capsicum annuum L

Lodging presents a significant challenge in cultivating high-yield crops with extensive above-ground biomass, yet the molecular mechanisms underlying this phenomenon in the Solanaceae family remain largely unexplored. In this study, we identified a gene, CaSLR1 (Capsicum annuum Stem Lodging Resistance 1), which encodes a MYELOBLASTOSIS (MYB) family transcription factor, from a lodging-affected C. annuum EMS mutant. The suppression of CaSLR1 expression in pepper led to notable stem lodging, reduced thickness of the secondary cell wall, and decreased stem strength. A similar phenotype was observed in tomato with the knockdown of SlMYB61, the orthologous gene to CaSLR1. Further investigations demonstrated that CaNAC6, a gene involved in secondary cell wall (SCW) formation, is co-expressed with CaSLR1 and acts as a positive regulator of its expression, as confirmed through yeast one-hybrid, dual-luciferase reporter assays, and electrophoretic mobility shift assays. These findings elucidate the CaNAC6-CaSLR1 module that contributes to lodging resistance, emphasizing the critical role of CaSLR1 in the lodging resistance regulatory network.

Cell-type-specific transcriptomics uncovers spatial regulatory networks in bioenergy sorghum stems

Bioenergy sorghum is a low-input, drought-resilient, deep-rooting annual crop that has high biomass yield potential enabling the sustainable production of biofuels, biopower, and bioproducts. Bioenergy sorghum's 4-5 m stems account for ~80% of the harvested biomass. Stems accumulate high levels of sucrose that could be used to synthesize bioethanol and useful biopolymers if information about cell-type gene expression and regulation in stems was available to enable engineering. To obtain this information, laser capture microdissection was used to isolate and collect transcriptome profiles from five major cell types that are present in stems of the sweet sorghum Wray. Transcriptome analysis identified genes with cell-type-specific and cell-preferred expression patterns that reflect the distinct metabolic, transport, and regulatory functions of each cell type. Analysis of cell-type-specific gene regulatory networks (GRNs) revealed that unique transcription factor families contribute to distinct regulatory landscapes, where regulation is organized through various modes and identifiable network motifs. Cell-specific transcriptome data was combined with known secondary cell wall (SCW) networks to identify the GRNs that differentially activate SCW formation in vascular sclerenchyma and epidermal cells. The spatial transcriptomic dataset provides a valuable source of information about the function of different sorghum cell types and GRNs that will enable the engineering of bioenergy sorghum stems, and an interactive web application developed during this project will allow easy access and exploration of the data (https://mc-lab.shinyapps.io/lcm-dataset/).

Transcriptome-wide expression landscape and starch synthesis pathway co-expression network in sorghum

The gene expression landscape across different tissues and developmental stages reflects their biological functions and evolutionary patterns. Integrative and comprehensive analyses of all transcriptomic data in an organism are instrumental to obtaining a comprehensive picture of gene expression landscape. Such studies are still very limited in sorghum, which limits the discovery of the genetic basis underlying complex agricultural traits in sorghum. We characterized the genome-wide expression landscape for sorghum using 873 RNA-sequencing (RNA-seq) datasets representing 19 tissues. Our integrative analysis of these RNA-seq data provides the most comprehensive transcriptomic atlas for sorghum, which will be valuable for the sorghum research community for functional characterizations of sorghum genes. Based on the transcriptome atlas, we identified 595 housekeeping genes (HKGs) and 2080 tissue-specific expression genes (TEGs) for the 19 tissues. We identified different gene features between HKGs and TEGs, and we found that HKGs have experienced stronger selective constraints than TEGs. Furthermore, we built a transcriptome-wide co-expression network (TW-CEN) comprising 35 modules with each module enriched in specific Gene Ontology terms. High-connectivity genes in TW-CEN tend to express at high levels while undergoing intensive selective pressure. We also built global and seed-preferential co-expression networks of starch synthesis pathways, which indicated that photosynthesis and microtubule-based movement play important roles in starch synthesis. The global transcriptome atlas of sorghum generated by this study provides an important functional genomics resource for trait discovery and insight into starch synthesis regulation in sorghum.

Transcriptomic analysis of developing sorghum grains to detect genes related to cell wall biosynthesis and remodelling

OBJECTIVE

Sorghum (Sorghum bicolor (L.) Moench) is the fifth most important grain produced in the world. Interest for cultivating sorghum is increasing all over the world in the context of climate change, due to its low input and water requirements. Like other cultivated cereals, sorghum has significant nutritional value thanks to its protein, carbohydrate and dietary fiber content, these latter mainly consisting of cell wall polysaccharides. This work describes for the first time a transcriptomic analysis dedicated to identify the genes involved in the biosynthesis and remodelling of cell walls both in the endosperm and outer layers of sorghum grain during its development. Further analysis of these transcriptomic data will improve our understanding of cell wall assembly, which is a key component of grain quality.

DATA DESCRIPTION

This research delineates the steps of our analysis, starting with the cultivation conditions and the grain harvest at different stages of development, followed by the laser microdissection applied to separate the endosperm from the outer layers. It also describes the procedures implemented to generate RNA libraries and to obtain a normalized and filtered table of transcript counts, and finally determine the number of putative cell wall-related genes already listed in literature.

Phenylpropanoids Following Wounding and Infection of Sweet Sorghum Lines Differing in Responses to Stalk Pathogens

Sweet sorghum (Sorghum bicolor) lines M81-E and Colman were previously shown to differ in responses to Fusarium thapsinum and Macrophomina phaseolina, stalk rot pathogens that can reduce the yields and quality of biomass and extracted sugars. Inoculated tissues were compared for transcriptomic, phenolic metabolite, and enzymatic activity during disease development 3 and 13 days after inoculation (DAI). At 13 DAI, M81-E had shorter mean lesion lengths than Colman when inoculated with either pathogen. Transcripts encoding monolignol biosynthetic and modification enzymes were associated with transcriptional wound (control) responses of both lines at 3 DAI. Monolignol biosynthetic genes were differentially coexpressed with transcriptional activator SbMyb76 in all Colman inoculations, but only following M. phaseolina inoculation in M81-E, suggesting that SbMyb76 is associated with lignin biosynthesis during pathogen responses. In control inoculations, defense-related genes were expressed at higher levels in M81-E than Colman. Line, treatment, and timepoint differences observed in phenolic metabolite and enzyme activities did not account for observed differences in lesions. However, generalized additive models were able to relate metabolites, but not enzyme activities, to lesion length for quantitatively modeling disease progression: in M81-E, but not Colman, sinapic acid levels positively predicted lesion length at 3 DAI when cell wall-bound syringic acid was low, soluble caffeic acid was high, and lactic acid was high, suggesting that sinapic acid may contribute to responses at 3 DAI. These results provide potential gene targets for development of sweet sorghum varieties with increased stalk rot resistance to ensure biomass and sugar quality.

Epicuticular wax accumulation and regulation of wax pathway gene expression during bioenergy Sorghum stem development

Bioenergy sorghum is a drought-tolerant high-biomass C4 grass targeted for production on annual cropland marginal for food crops due primarily to abiotic constraints. To better understand the overall contribution of stem wax to bioenergy sorghum's resilience, the current study characterized sorghum stem cuticular wax loads, composition, morphometrics, wax pathway gene expression and regulation using vegetative phase Wray, R07020, and TX08001 genotypes. Wax loads on sorghum stems (~103-215 µg/cm2) were much higher than Arabidopsis stem and leaf wax loads. Wax on developing sorghum stem internodes was enriched in C28/30 primary alcohols (~65%) while stem wax on fully developed stems was enriched in C28/30 aldehydes (~80%). Scanning Electron Microscopy showed minimal wax on internodes prior to the onset of elongation and that wax tubules first appear associated with cork-silica cell complexes when internode cell elongation is complete. Sorghum homologs of genes involved in wax biosynthesis/transport were differentially expressed in the stem epidermis. Expression of many wax pathway genes (i.e., SbKCS6, SbCER3-1, SbWSD1, SbABCG12, SbABCG11) is low in immature apical internodes then increases at the onset of stem wax accumulation. SbCER4 is expressed relatively early in stem development consistent with accumulation of C28/30 primary alcohols on developing apical internodes. High expression of two SbCER3 homologs in fully elongated internodes is consistent with a role in production of C28/30 aldehydes. Gene regulatory network analysis aided the identification of sorghum homologs of transcription factors that regulate wax biosynthesis (i.e., SbSHN1, SbWRI1/3, SbMYB94/96/30/60, MYS1) and other transcription factors that could regulate and specify expression of the wax pathway in epidermal cells during cuticle development.

Harnessing the Genetic Basis of Sorghum Biomass-Related Traits to Facilitate Bioenergy Applications

The extensive use of fossil fuels and global climate change have raised ever-increasing attention to sustainable development, global food security and the replacement of fossil fuels by renewable energy. Several C4 monocot grasses have excellent photosynthetic ability, stress tolerance and may rapidly produce biomass in marginal lands with low agronomic inputs, thus representing an important source of bioenergy. Among these grasses, Sorghum bicolor has been recognized as not only a promising bioenergy crop but also a research model due to its diploidy, simple genome, genetic diversity and clear orthologous relationship with other grass genomes, allowing sorghum research to be easily translated to other grasses. Although sorghum molecular genetic studies have lagged far behind those of major crops (e.g., rice and maize), recent advances have been made in a number of biomass-related traits to dissect the genetic loci and candidate genes, and to discover the functions of key genes. However, molecular and/or targeted breeding toward biomass-related traits in sorghum have not fully benefited from these pieces of genetic knowledge. Thus, to facilitate the breeding and bioenergy applications of sorghum, this perspective summarizes the bioenergy applications of different types of sorghum and outlines the genetic control of the biomass-related traits, ranging from flowering/maturity, plant height, internode morphological traits and metabolic compositions. In particular, we describe the dynamic changes of carbohydrate metabolism in sorghum internodes and highlight the molecular regulators involved in the different stages of internode carbohydrate metabolism, which affects the bioenergy utilization of sorghum biomass. We argue the way forward is to further enhance our understanding of the genetic mechanisms of these biomass-related traits with new technologies, which will lead to future directions toward tailored designing sorghum biomass traits suitable for different bioenergy applications.

Unleashing the power within short-read RNA-seq for plant research: Beyond differential expression analysis and toward regulomics

RNA-seq has become a state-of-the-art technique for transcriptomic studies. Advances in both RNA-seq techniques and the corresponding analysis tools and pipelines have unprecedently shaped our understanding in almost every aspects of plant sciences. Notably, the integration of huge amount of RNA-seq with other omic data sets in the model plants and major crop species have facilitated plant regulomics, while the RNA-seq analysis has still been primarily used for differential expression analysis in many less-studied plant species. To unleash the analytical power of RNA-seq in plant species, especially less-studied species and biomass crops, we summarize recent achievements of RNA-seq analysis in the major plant species and representative tools in the four types of application: (1) transcriptome assembly, (2) construction of expression atlas, (3) network analysis, and (4) structural alteration. We emphasize the importance of expression atlas, coexpression networks and predictions of gene regulatory relationships in moving plant transcriptomes toward regulomics, an omic view of genome-wide transcription regulation. We highlight what can be achieved in plant research with RNA-seq by introducing a list of representative RNA-seq analysis tools and resources that are developed for certain minor species or suitable for the analysis without species limitation. In summary, we provide an updated digest on RNA-seq tools, resources and the diverse applications for plant research, and our perspective on the power and challenges of short-read RNA-seq analysis from a regulomic point view. A full utilization of these fruitful RNA-seq resources will promote plant omic research to a higher level, especially in those less studied species.

Transcriptional and metabolic changes associated with internode development and reduced cinnamyl alcohol dehydrogenase activity in sorghum

The molecular mechanisms associated with secondary cell wall (SCW) deposition in sorghum remain largely uncharacterized. Here, we employed untargeted metabolomics and large-scale transcriptomics to correlate changes in SCW deposition with variation in global gene expression profiles and metabolite abundance along an elongating internode of sorghum, with a major focus on lignin and phenolic metabolism. To gain deeper insight into the metabolic and transcriptional changes associated with pathway perturbations, a bmr6 mutant [with reduced cinnamyl alcohol dehydrogenase (CAD) activity] was analyzed. In the wild type, internode development was accompanied by an increase in the content of oligolignols, p-hydroxybenzaldehyde, hydroxycinnamate esters, and flavonoid glucosides, including tricin derivatives. We further identified modules of genes whose expression pattern correlated with SCW deposition and the accumulation of these target metabolites. Reduced CAD activity resulted in the accumulation of hexosylated forms of hydroxycinnamates (and their derivatives), hydroxycinnamaldehydes, and benzenoids. The expression of genes belonging to one specific module in our co-expression analysis correlated with the differential accumulation of these compounds and contributed to explaining this metabolic phenotype. Metabolomics and transcriptomics data further suggested that CAD perturbation activates distinct detoxification routes in sorghum internodes. Our systems biology approach provides a landscape of the metabolic and transcriptional changes associated with internode development and with reduced CAD activity in sorghum.

PlantNexus: A Gene Co-expression Network Database and Visualization Tool for Barley and Sorghum

Global gene co-expression networks (GCNs) are powerful tools for functional genomics whereby putative functions and regulatory mechanisms can be inferred by gene co-expression. Cereal crops, such as Hordeum vulgare (barley) and Sorghum bicolor (sorghum), are among the most important plants to civilization. However, co-expression network tools for these plants are lacking. Here, we have constructed global GCNs for barley and sorghum using existing RNA-seq data sets. Meta-information was manually curated and categorized by tissue type to also build tissue-specific GCNs. To enable GCN searching and visualization, we implemented a website and database named PlantNexus. PlantNexus is freely available at https://plantnexus.ohio.edu/.

Profiling Alternative 3' Untranslated Regions in Sorghum using RNA-seq Data

Sorghum is an important crop widely used for food, feed, and fuel. Transcriptome-wide studies of 3′ untranslated regions (3′UTR) using regular RNA-seq remain scarce in sorghum, while transcriptomes have been characterized extensively using Illumina short-read sequencing platforms for many sorghum varieties under various conditions or developmental contexts. 3′UTR is a critical regulatory component of genes, controlling the translation, transport, and stability of messenger RNAs. In the present study, we profiled the alternative 3′UTRs at the transcriptome level in three genetically related but phenotypically contrasting lines of sorghum: Rio, BTx406, and R9188. A total of 1,197 transcripts with alternative 3′UTRs were detected using RNA-seq data. Their categorization identified 612 high-confidence alternative 3′UTRs. Importantly, the high-confidence alternative 3′UTR genes significantly overlapped with the genesets that are associated with RNA N⁶-methyladenosine (m⁶A) modification, suggesting a clear indication between alternative 3′UTR and m⁶A methylation in sorghum. Moreover, taking advantage of sorghum genetics, we provided evidence of genotype specificity of alternative 3′UTR usage. In summary, our work exemplifies a transcriptome-wide profiling of alternative 3′UTRs using regular RNA-seq data in non-model crops and gains insights into alternative 3′UTRs and their genotype specificity.