Dynamics of gene expression during development and expansion of vegetative stem internodes of bioenergy sorghum

Transcriptome and metabolome analyses reveal regulatory networks associated with nutrition synthesis in sorghum seeds

Cereal seeds are vital for food, feed, and agricultural sustainability because they store and provide essential nutrients to human and animal food and feed systems. Unraveling molecular processes in seed development is crucial for enhancing cereal grain yield and quality. We analyze spatiotemporal transcriptome and metabolome profiles during sorghum seed development in the inbred line 'BTx623'. Morphological and molecular analyses identify the key stages of seed maturation, specifying starch biosynthesis onset at 5 days post-anthesis (dpa) and protein at 10 dpa. Transcriptome profiling from 1 to 25 dpa reveal dynamic gene expression pathways, shifting from cellular growth and embryo development (1-5 dpa) to cell division, fatty acid biosynthesis (5-25 dpa), and seed storage compounds synthesis in the endosperm (5-25 dpa). Network analysis identifies 361 and 207 hub genes linked to starch and protein synthesis in the endosperm, respectively, which will help breeders enhance sorghum grain quality. The availability of this data in the sorghum reference genome line establishes a baseline for future studies as new pangenomes emerge, which will consider copy number and presence-absence variation in functional food traits.

Low-cost and reliable substrate-based phenotyping platform for screening salt tolerance of cutting propagation-dependent grass, paspalum vaginatum

BACKGROUND

Salt tolerance in plants is defined as their ability to grow and complete their life cycle under saline conditions. Staple crops have limited salt tolerance, but forage grass can survive in large unexploited saline areas of costal or desert land. However, due to the restriction of self-incompatible fertilization in many grass species, vegetative propagation via stem cuttings is the dominant practice; this is incompatible with current methodologies of salt-tolerance phenotyping, which have been developed for germination-based seedling growth. Therefore, the performance of seedlings from cuttings under salt stress is still fuzzy. Moreover, the morphological traits involved in salt tolerance are still mostly unknown, especially under experimental conditions with varying levels of stress.

RESULTS

To estimate the salt tolerance of cutting propagation-dependent grasses, a reliable and low-cost workflow was established with multiple saline treatments, using Paspalum vaginatum as the material and substrate as medium, where cold stratification and selection of stem segments were the two variables used to control for experimental errors. Average leaf number (ALN) was designated as the best criterion for evaluating ion-accumulated salt tolerance. The reliability of ALN was revealed by the consistent results among four P. vaginatum genotypes, and three warm-season (pearl millet, sweet sorghum, and wild maize) and four cold-season (barley, oat, rye, and ryegrass) forage cultivars. Dynamic curves simulated by sigmoidal mathematical models were well-depicted for the calculation of the key parameter, Salt50. The reliability of the integrated platform was further validated by screening 48 additional recombinants, which were previously generated from a self-fertile mutant of P. vaginatum. The genotypes displaying extreme ALN-based Salt50 also exhibited variations in biomass and ion content, which not only confirmed the reliability of our phenotyping platform but also the representativeness of the aerial ALN trait for salt tolerance.

CONCLUSIONS

Our phenotyping platform is proved to be compatible with estimations in both germination-based and cutting propagation-dependent seedling tolerance under salt stresses. ALN and its derived parameters are prone to overcome the species barriers when comparing salt tolerance of different species together. The accuracy and reliability of the developed phenotyping platform is expected to benefit breeding programs in saline agriculture.

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/).

Integrated VIS/NIR Spectrum and Genome-Wide Association Study for Genetic Dissection of Cellulose Crystallinity in Wheat Stems

Cellulose crystallinity is a crucial factor influencing stem strength and, consequently, wheat lodging. However, the genetic dissection of cellulose crystallinity is less reported due to the difficulty of its measurement. In this study, VIS/NIR spectra and cellulose crystallinity were measured for a wheat accession panel with diverse genetic backgrounds. We developed a reliable VIS/NIR model for cellulose crystallinity with a high determination coefficient (R2) (0.95) and residual prediction deviation (RPD) (4.04), enabling the rapid screening of wheat samples. A GWAS of the cellulose crystallinity in 326 wheat accessions revealed 14 significant SNPs and 13 QTLs. Two candidate genes, TraesCS4B03G0029800 and TraesCS5B03G1085500, were identified. In summary, this study establishes an efficient method for the measurement of cellulose crystallinity in wheat stems and provides a genetic basis for enhancing lodging resistance in wheat.

The Physiological Basis of Alfalfa Plant Height Establishment

Plant height plays an important role in crop yield, product quality, and cultivation management. However, the physiological mechanisms that regulate the establishment of plant height in alfalfa plants remain unclear. Herein, we measured plant height traits, leaf characteristics, photosynthetic physiology, cell wall composition, and endogenous hormone contents of tall- and short-stalked alfalfa materials at different reproductive periods. We analyzed the physiology responsible for differences in plant height. The results demonstrated that the number of internodes in tall- and short-stalked alfalfa materials tended to converge with the advancement of the fertility period. Meanwhile, the average internode length (IL) of tall-stalked materials was significantly higher than that of short-stalked materials at different fertility periods, with internode length identified as the main trait determining the differences in alfalfa plant height. Leaf characteristics, which are closely related to photosynthetic capacity, are crucial energy sources supporting the expression of plant height traits, and we found that an increase in the number of leaves contributed to a proportional increase in plant height. Additionally, a significant positive correlation was observed between plant height and leaf dry weight per plant during the branching and early flowering stages of alfalfa. The leaves of alfalfa affect plant height through photosynthesis, with the budding stage identified as the key period for efficient light energy utilization. Plant height at the budding stage showed a significant positive correlation with soluble sugar (SS) content and a significant negative correlation with intercellular CO2 concentration. Moreover, we found that alfalfa plant height was significantly correlated with the contents of indole-3-acetic acid in stem tips (SIAA), gibberellin A3 in leaves (LGA3), zeatin in stem tips (SZT), and abscisic acid in leaves (LABA). Further investigation revealed that SS, SIAA, and LGA3 contents were important physiological indicators affecting alfalfa plant height. This study provides a theoretical basis for understanding the formation of alfalfa plant height traits and for genetic improvement studies.

Sorghum: A Multipurpose Crop

Sorghum (Sorghum bicolor L.) is one of the top five cereal crops in the world in terms of production and planting area and is widely grown in areas with severe abiotic stresses such as drought and saline-alkali land due to its excellent stress resistance. Moreover, sorghum is a rare multipurpose crop that can be classified into grain sorghum, energy sorghum, and silage sorghum according to its domestication direction and utilization traits, endowing it with broad breeding and economic value. In this review, we mainly discuss the latest research progress and regulatory genes of agronomic traits of sorghum as a grain, energy, and silage crop, as well as the future improvement direction of multipurpose sorghum. We also emphasize the feasibility of cultivating multipurpose sorghum through genetic engineering methods by exploring potential targets using wild sorghum germplasm and genetic resources, as well as genomic resources.

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.

Ten new high-quality genome assemblies for diverse bioenergy sorghum genotypes

INTRODUCTION

Sorghum (Sorghum bicolor (L.) Moench) is an agriculturally and economically important staple crop that has immense potential as a bioenergy feedstock due to its relatively high productivity on marginal lands. To capitalize on and further improve sorghum as a potential source of sustainable biofuel, it is essential to understand the genomic mechanisms underlying complex traits related to yield, composition, and environmental adaptations.

METHODS

Expanding on a recently developed mapping population, we generated de novo genome assemblies for 10 parental genotypes from this population and identified a comprehensive set of over 24 thousand large structural variants (SVs) and over 10.5 million single nucleotide polymorphisms (SNPs).

RESULTS

We show that SVs and nonsynonymous SNPs are enriched in different gene categories, emphasizing the need for long read sequencing in crop species to identify novel variation. Furthermore, we highlight SVs and SNPs occurring in genes and pathways with known associations to critical bioenergy-related phenotypes and characterize the landscape of genetic differences between sweet and cellulosic genotypes.

DISCUSSION

These resources can be integrated into both ongoing and future mapping and trait discovery for sorghum and its myriad uses including food, feed, bioenergy, and increasingly as a carbon dioxide removal mechanism.

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.

Bioenergy sorghum stem growth regulation: intercalary meristem localization, development, and gene regulatory network analysis

Bioenergy sorghum is a highly productive drought tolerant C₄ grass that accumulates 80% of its harvestable biomass in approximately 4 m length stems. Stem internode growth is regulated by development, shading, and hormones that modulate cell proliferation in intercalary meristems (IMs). In this study, sorghum stem IMs were localized above the pulvinus at the base of elongating internodes using magnetic resonance imaging, microscopy, and transcriptome analysis. A change in cell morphology/organization occurred at the junction between the pulvinus and internode where LATERAL ORGAN BOUNDARIES (SbLOB), a boundary layer gene, was expressed. Inactivation of an AGCVIII kinase in DDYM (dw2) resulted in decreased SbLOB expression, disrupted IM localization, and reduced internode cell proliferation. Transcriptome analysis identified approximately 1000 genes involved in cell proliferation, hormone signaling, and other functions selectively upregulated in the IM compared with a non-meristematic stem tissue. This cohort of genes is expressed in apical dome stem tissues before localization of the IM at the base of elongating internodes. Gene regulatory network analysis identified connections between genes involved in hormone signaling and cell proliferation. The results indicate that gibberellic acid induces accumulation of growth regulatory factors (GRFs) known to interact with ANGUSTIFOLIA (SbAN3), a master regulator of cell proliferation. GRF:AN3 was predicted to induce SbARF3/ETT expression and regulate SbAN3 expression in an auxin-dependent manner. GRFs and ARFs regulate genes involved in cytokinin and brassinosteroid signaling and cell proliferation. The results provide a molecular framework for understanding how hormone signaling regulates the expression of genes involved in cell proliferation in the stem IM.

Identification of Candidate Forage Yield Genes in Sorghum (Sorghum bicolor L.) Using Integrated Genome-Wide Association Studies and RNA-Seq

Genetic dissection of forage yield traits is critical to the development of sorghum as a forage crop. In the present study, association mapping was performed with 85,585 SNP markers on four forage yield traits, namely plant height (PH), tiller number (TN), stem diameter (SD), and fresh weight per plant (FW) among 245 sorghum accessions evaluated in four environments. A total of 338 SNPs or quantitative trait nucleotides (QTNs) were associated with the four traits, and 21 of these QTNs were detected in at least two environments, including four QTNs for PH, ten for TN, six for SD, and one for FW. To identify candidate genes, dynamic transcriptome expression profiling was performed at four stages of sorghum development. One hundred and six differentially expressed genes (DEGs) that were enriched in hormone signal transduction pathways were found in all stages. Weighted gene correlation network analysis for PH and SD indicated that eight modules were significantly correlated with PH and that three modules were significantly correlated with SD. The blue module had the highest positive correlation with PH and SD, and the turquoise module had the highest negative correlation with PH and SD. Eight candidate genes were identified through the integration of genome-wide association studies (GWAS) and RNA sequencing. Sobic.004G143900, an indole-3-glycerol phosphate synthase gene that is involved in indoleacetic acid biosynthesis, was down-regulated as sorghum plants grew in height and was identified in the blue module, and Sobic.003G375100, an SD candidate gene, encoded a DNA repair RAD52-like protein 1 that plays a critical role in DNA repair-linked cell cycle progression. These findings demonstrate that the integrative analysis of omics data is a promising approach to identify candidate genes for complex traits.

Growth and yield performance of sorghum (Sorghum bicolor L.) crop under anthracnose stress in dryland crop-livestock farming system

Dual-purpose sorghum response to anthracnose disease, growth, and yield was undertaken in Derashe and Arba Minch trial sites during March–June 2018 and 2019. Five sorghum varieties and Rara (local check) were arranged in a randomized complete block design with four replications. Variety Chelenko exhibited the tallest main crop plant height (430 cm) while Dishkara was the tallest (196.65 cm) at ratoon crop harvesting. Rara had a higher tiller number (main = 6.73, ratoon = 9.73) among the varieties. Dishkara and Chelenko varieties produced 50 and 10% more dry biomass yield (DBY) than the overall mean DBY, while Konoda produced 40% less. Although the anthracnose infestation was highest on the varieties Konoda (percentage severity index [PSI] = 20.37%) and NTJ_2 (PSI = 32.19%), they produced significantly (p < .001) higher grain yield (3.89 t/ha) than others. Under anthracnose pressure, Chelenko and Dishkara varieties are suggested for dry matter yield while NTJ_2 for grain yield production in the study area and similar agroecology.

Regulation of dhurrin pathway gene expression during Sorghum bicolor development

Developmental and organ-specific expression of genes in dhurrin biosynthesis, bio-activation, and recycling offers dynamic metabolic responses optimizing growth and defence responses in Sorghum. Plant defence models evaluate the costs and benefits of resource investments at different stages in the life cycle. Poor understanding of the molecular regulation of defence deployment and remobilization hampers accuracy of the predictions. Cyanogenic glucosides, such as dhurrin are phytoanticipins that release hydrogen cyanide upon bio-activation. In this study, RNA-seq was used to investigate the expression of genes involved in the biosynthesis, bio-activation and recycling of dhurrin in Sorghum bicolor. Genes involved in dhurrin biosynthesis were highly expressed in all young developing vegetative tissues (leaves, leaf sheath, roots, stems), tiller buds and imbibing seeds and showed gene specific peaks of expression in leaves during diel cycles. Genes involved in dhurrin bio-activation were expressed early in organ development with organ-specific expression patterns. Genes involved in recycling were expressed at similar levels in the different organ during development, although post-floral initiation when nutrients are remobilized for grain filling, expression of GSTL1 decreased > tenfold in leaves and NITB2 increased > tenfold in stems. Results are consistent with the establishment of a pre-emptive defence in young tissues and regulated recycling related to organ senescence and increased demand for nitrogen during grain filling. This detailed characterization of the transcriptional regulation of dhurrin biosynthesis, bioactivation and remobilization genes during organ and plant development will aid elucidation of gene regulatory networks and signalling pathways that modulate gene expression and dhurrin levels. In-depth knowledge of dhurrin metabolism could improve the yield, nitrogen use efficiency and stress resilience of Sorghum.

Chia (Salvia hispanica) Gene Expression Atlas Elucidates Dynamic Spatio-Temporal Changes Associated With Plant Growth and Development

Chia (Salvia hispanica L.), now a popular superfood and a pseudocereal, is one of the richest sources of dietary nutrients such as protein, fiber, and polyunsaturated fatty acids (PUFAs). At present, the genomic and genetic information available in the public domain for this crop are scanty, which hinders an understanding of its growth and development and genetic improvement. We report an RNA-sequencing (RNA-Seq)-based comprehensive transcriptome atlas of Chia sampled from 13 tissue types covering vegetative and reproductive growth stages. We used ~355 million high-quality reads of total ~394 million raw reads from transcriptome sequencing to generate de novo reference transcriptome assembly and the tissue-specific transcript assemblies. After the quality assessment of the merged assemblies and implementing redundancy reduction methods, 82,663 reference transcripts were identified. About 65,587 of 82,663 transcripts were translated into 99,307 peptides, and we were successful in assigning InterPro annotations to 45,209 peptides and gene ontology (GO) terms to 32,638 peptides. The assembled transcriptome is estimated to have the complete sequence information for ~86% of the genes found in the Chia genome. Furthermore, the analysis of 53,200 differentially expressed transcripts (DETs) revealed their distinct expression patterns in Chia's vegetative and reproductive tissues; tissue-specific networks and developmental stage-specific networks of transcription factors (TFs); and the regulation of the expression of enzyme-coding genes associated with important metabolic pathways. In addition, we identified 2,411 simple sequence repeats (SSRs) as potential genetic markers from the transcripts. Overall, this study provides a comprehensive transcriptome atlas, and SSRs, contributing to building essential genomic resources to support basic research, genome annotation, functional genomics, and molecular breeding of Chia.

The AGCVIII kinase Dw2 modulates cell proliferation, endomembrane trafficking, and MLG/xylan cell wall localization in elongating stem internodes of Sorghum bicolor

Stems of bioenergy sorghum (Sorghum bicolor L. Moench.), a drought-tolerant C4 grass, contain up to 50 nodes and internodes of varying length that span 4-5 m and account for approximately 84% of harvested biomass. Stem internode growth impacts plant height and biomass accumulation and is regulated by brassinosteroid signaling, auxin transport, and gibberellin biosynthesis. In addition, an AGCVIII kinase (Dw2) regulates sorghum stem internode growth, but the underlying mechanism and signaling network are unknown. Here we provide evidence that mutation of Dw2 reduces cell proliferation in internode intercalary meristems, inhibits endocytosis, and alters the distribution of heteroxylan and mixed linkage glucan in cell walls. Phosphoproteomic analysis showed that Dw2 signaling influences the phosphorylation of proteins involved in lipid signaling (PLDδ), endomembrane trafficking, hormone, light, and receptor signaling, and photosynthesis. Together, our results show that Dw2 modulates endomembrane function and cell division during sorghum internode growth, providing insight into the regulation of monocot stem development.

High planting density induces the expression of GA3-oxidase in leaves and GA mediated stem elongation in bioenergy sorghum

The stems of bioenergy sorghum hybrids at harvest are > 4 m long, contain > 40 internodes and account for ~ 80% of harvested biomass. In this study, bioenergy sorghum hybrids were grown at four planting densities (~ 20,000 to 132,000 plants/ha) under field conditions for 60 days to investigate the impact shading has on stem growth and biomass accumulation. Increased planting density induced a > 2-fold increase in sorghum internode length and a ~ 22% decrease in stem diameter, a typical shade avoidance response. Shade-induced internode elongation was due to an increase in cell length and number of cells spanning the length of internodes. SbGA3ox2 (Sobic.003G045900), a gene encoding the last step in GA biosynthesis, was expressed ~ 20-fold higher in leaf collar tissue of developing phytomers in plants grown at high vs. low density. Application of GA3 to bioenergy sorghum increased plant height, stem internode length, cell length and the number of cells spanning internodes. Prior research showed that sorghum plants lacking phytochrome B, a key photoreceptor involved in shade signaling, accumulated more GA1 and displayed shade avoidance phenotypes. These results are consistent with the hypothesis that increasing planting density induces expression of GA3-oxidase in leaf collar tissue, increasing synthesis of GA that stimulates internode elongation.

The Wild Sugarcane and Sorghum Kinomes: Insights Into Expansion, Diversification, and Expression Patterns

The protein kinase (PK) superfamily is one of the largest superfamilies in plants and the core regulator of cellular signaling. Despite this substantial importance, the kinomes of sugarcane and sorghum have not been profiled. Here, we identified and profiled the complete kinomes of the polyploid Saccharum spontaneum (Ssp) and Sorghum bicolor (Sbi), a close diploid relative. The Sbi kinome was composed of 1,210 PKs; for Ssp, we identified 2,919 PKs when disregarding duplications and allelic copies, and these were related to 1,345 representative gene models. The Ssp and Sbi PKs were grouped into 20 groups and 120 subfamilies and exhibited high compositional similarities and evolutionary divergences. By utilizing the collinearity between the species, this study offers insights into Sbi and Ssp speciation, PK differentiation and selection. We assessed the PK subfamily expression profiles via RNA-Seq and identified significant similarities between Sbi and Ssp. Moreover, coexpression networks allowed inference of a core structure of kinase interactions with specific key elements. This study provides the first categorization of the allelic specificity of a kinome and offers a wide reservoir of molecular and genetic information, thereby enhancing the understanding of Sbi and Ssp PK evolutionary history.

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.

Phylogenomic Analysis of R2R3 MYB Transcription Factors in Sorghum and their Role in Conditioning Biofuel Syndrome

Background

Large scale cultivation of sorghum for food, feed, and biofuel requires concerted efforts for engineering multipurpose cultivars with optimised agronomic traits. Due to their vital role in regulating the biosynthesis of phenylpropanoid-derived compounds, biomass composition, biotic, and abiotic stress response, R2R3-MYB family transcription factors are ideal targets for improving environmental resilience and economic value of sorghum.

Methods

We used diverse computational biology tools to survey the sorghum genome to identify R2R3-MYB transcription factors followed by their structural and phylogenomic analysis. We used in-house generated as well as publicly available high throughput expression data to analyse the R2R3 expression patterns in various sorghum tissue types.

Results

We have identified a total of 134 R2R3-MYB genes from sorghum and developed a framework to predict gene functions. Collating information from the physical location, duplication, structural analysis, orthologous sequences, phylogeny, and expression patterns revealed the role of duplications in clade-wise expansion of the R2R3-MYB family as well as intra-clade functional diversification. Using publicly available and in-house generated RNA sequencing data, we provide MYB candidates for conditioning biofuel syndrome by engineering phenylpropanoid biosynthesis and sugar signalling pathways in sorghum.

Conclusion

The results presented here are pivotal to prioritize MYB genes for functional validation and optimize agronomic traits in sorghum.

Shade signals alter the expression of circadian clock genes in newly-formed bioenergy sorghum internodes

Stem internodes of bioenergy sorghum inbred R.07020 are longer at high plant density (shade) than at low plant density (control). Initially, the youngest newly-formed subapical stem internodes of shade-treated and control plants are comparable in length. However, full-length internodes of shade-treated plants are three times longer than the internodes of the control plants. To identify the early molecular events associated with internode elongation in response to shade, we analyzed the transcriptome of the newly-formed internodes of shade-treated and control plants sampled between 4 and 6 hr after the start of the light period (14 hr light/10 hr dark). Sorghum genes homologous to the Arabidopsis shade marker genes ATHB2 and PIL1 were not differentially expressed. The results indicate that shade signals promote internode elongation indirectly because sorghum internodes are not illuminated and grow while enclosed with leaf sheaths. Sorghum genes homologous to the Arabidopsis morning-phased circadian clock genes LHY, RVE, and LNK were downregulated and evening-phased genes such as TOC1, PRR5, and GI were upregulated in young internodes in response to shade. We hypothesize that a change in the function or patterns of expression of the circadian clock genes is the earliest molecular event associated with internode elongation in response to shade in bioenergy sorghum. Increased expression of CycD1, which promotes cell division, and decreased expression of cell wall-loosening and MBF1-like genes, which promote cell expansion, suggest that shade signals promote internode elongation in bioenergy sorghum in part through increasing cell number by delaying transition from cell division to cell expansion.

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

Most sorghum biomass accumulates in stem secondary cell walls (SCW). As sorghum stems are used as raw materials for various purposes such as feed, energy and fiber reinforced polymers, identifying the genes responsible for SCW establishment is highly important. Taking advantage of studies performed in model species, most of the structural genes contributing at the molecular level to the SCW biosynthesis in sorghum have been proposed while their regulatory factors have mostly not been determined. Validation of the role of several MYB and NAC transcription factors in SCW regulation in Arabidopsis and a few other species has been provided. In this study, we contributed to the recent efforts made in grasses to uncover the mechanisms underlying SCW establishment. We reported updated phylogenies of NAC and MYB in 9 different species and exploited findings from other species to highlight candidate regulators of SCW in sorghum. We acquired expression data during sorghum internode development and used co-expression analyses to determine groups of co-expressed genes that are likely to be involved in SCW establishment. We were able to identify two groups of co-expressed genes presenting multiple evidences of involvement in SCW building. Gene enrichment analysis of MYB and NAC genes provided evidence that while NAC SECONDARY WALL THICKENING PROMOTING FACTOR NST genes and SECONDARY WALL-ASSOCIATED NAC DOMAIN PROTEIN gene functions appear to be conserved in sorghum, NAC master regulators of SCW in sorghum may not be as tissue compartmentalized as in Arabidopsis. We showed that for every homolog of the key SCW MYB in Arabidopsis, a similar role is expected for sorghum. In addition, we unveiled sorghum MYB and NAC that have not been identified to date as being involved in cell wall regulation. Although specific validation of the MYB and NAC genes uncovered in this study is needed, we provide a network of sorghum genes involved in SCW both at the structural and regulatory levels.

Microscale thermophoresis as a powerful tool for screening glycosyltransferases involved in cell wall biosynthesis

BACKGROUND

Identification and characterization of key enzymes associated with cell wall biosynthesis and modification is fundamental to gain insights into cell wall dynamics. However, it is a challenge that activity assays of glycosyltransferases are very low throughput and acceptor substrates are generally not available.

RESULTS

We optimized and validated microscale thermophoresis (MST) to achieve high throughput screening for glycosyltransferase substrates. MST is a powerful method for the quantitative analysis of protein-ligand interactions with low sample consumption. The technique is based on the motion of molecules along local temperature gradients, measured by fluorescence changes. We expressed glycosyltransferases as YFP-fusion proteins in tobacco and optimized the MST method to allow the determination of substrate binding affinity without purification of the target protein from the cell lysate. The application of this MST method to the β-1,4-galactosyltransferase AtGALS1 validated the capability to screen both nucleotide-sugar donor substrates and acceptor substrates. We also expanded the application to members of glycosyltransferase family GT61 in sorghum for substrate screening and function prediction.

CONCLUSIONS

This method is rapid and sensitive to allow determination of both donor and acceptor substrates of glycosyltransferases. MST enables high throughput screening of glycosyltransferases for likely substrates, which will narrow down their in vivo function and help to select candidates for further studies. Additionally, this method gives insight into biochemical mechanism of glycosyltransferase function.

Common metabolic networks contribute to carbon sink strength of sorghum internodes: implications for bioenergy improvement

BACKGROUND

Sorghum bicolor (L.) is an important bioenergy source. The stems of sweet sorghum function as carbon sinks and accumulate large amounts of sugars and lignocellulosic biomass and considerable amounts of starch, therefore providing a model of carbon allocation and accumulation for other bioenergy crops. While omics data sets for sugar accumulation have been reported in different genotypes, the common features of primary metabolism in sweet genotypes remain unclear. To obtain a cohesive and comparative picture of carbohydrate metabolism between sorghum genotypes, we compared the phenotypes and transcriptome dynamics of sugar-accumulating internodes among three different sweet genotypes (Della, Rio, and SIL-05) and two non-sweet genotypes (BTx406 and R9188).

RESULTS

Field experiments showed that Della and Rio had similar dynamics and internode patterns of sugar concentration, albeit distinct other phenotypes. Interestingly, cellulose synthases for primary cell wall and key genes in starch synthesis and degradation were coordinately upregulated in sweet genotypes. Sweet sorghums maintained active monolignol biosynthesis compared to the non-sweet genotypes. Comparative RNA-seq results support the role of candidate Tonoplast Sugar Transporter gene (TST), but not the Sugars Will Eventually be Exported Transporter genes (SWEETs) in the different sugar accumulations between sweet and non-sweet genotypes.

CONCLUSIONS

Comparisons of the expression dynamics of carbon metabolic genes across the RNA-seq data sets identify several candidate genes with contrasting expression patterns between sweet and non-sweet sorghum lines, including genes required for cellulose and monolignol synthesis (CesA, PTAL, and CCR), starch metabolism (AGPase, SS, SBE, and G6P-translocator SbGPT2), and sucrose metabolism and transport (TPP and TST2). The common transcriptome features of primary metabolism identified here suggest the metabolic networks contributing to carbon sink strength in sorghum internodes, prioritize the candidate genes for manipulating carbon allocation with bioenergy purposes, and provide a comparative and cohesive picture of the complexity of carbon sink strength in sorghum stem.

APETALA2 control of barley internode elongation

Many plants dramatically elongate their stems during flowering, yet how this response is coordinated with the reproductive phase is unclear. We demonstrate that microRNA (miRNA) control of APETALA2 (AP2) is required for rapid, complete elongation of stem internodes in barley, especially of the final ‘peduncle’ internode directly underneath the inflorescence. Disrupted miR172 targeting of AP2 in the Zeo1.b barley mutant caused lower mitotic activity, delayed growth dynamics and premature lignification in the peduncle leading to fewer and shorter cells. Stage- and tissue-specific comparative transcriptomics between Zeo1.b and its parent cultivar showed reduced expression of proliferation-associated genes, ectopic expression of maturation-related genes and persistent, elevated expression of genes associated with jasmonate and stress responses. We further show that applying methyl jasmonate (MeJA) phenocopied the stem elongation of Zeo1.b, and that Zeo1.b itself was hypersensitive to inhibition by MeJA but less responsive to promotion by gibberellin. Taken together, we propose that miR172-mediated restriction of AP2 may modulate the jasmonate pathway to facilitate gibberellin-promoted stem growth during flowering.

Summary: Regulation of reproductive stem elongation in barley by APETALA2 suggests a pivotal role for phase change repression of JA-associated responses to promote internode growth.

Multi-omics analysis of the development and fracture resistance for maize internode

The maize stalk is an important mechanical supporting tissue. The stalk fracture resistance is closely related to lodging resistance, and thus the yield. In this study, we showed that the basal zone (BZ) was more fragile than the middle zone (MZ) of the stalk internode before tasseling. In order to clarify the relationship between the different zones and fragile resistance between the internodes, we systematically analyzed the phenotypic, metabolomic and transcriptomic differences. The results indicated that the BZ zone had lower stalk strength, which corresponded to the results of less lignin, cellulose and hemicellulose than that of the MZ. The 27 highly enriched metabolites and 4430 highly expressed genes in the BZ mainly participated in pentose phosphate, and in ribosome and sterol synthesis pathways, respectively. In addition, the BZ had higher vascular bundles density but smaller size compared with the MZ. By contrast, the 28 highly enriched known metabolites and 4438 highly expressed genes in the MZ were mainly involved in lignin synthesis, and secondary metabolites synthesis, respectively, especially the phenylpropanoid synthesis. The results provide a deeper understanding of the relationship between development and fracture differences in stalk, and may facilitate the improvement of field management practice to reduce lodging.

How plants grow up

A plant's lateral structures, such as leaves, branches and flowers, literally hinge on the shoot axis, making its integrity and growth fundamental to plant form. In all plants, subapical proliferation within the shoot tip displaces cells downward to extrude the cylindrical stem. Following the transition to flowering, many plants show extensive axial elongation associated with increased subapical proliferation and expansion. However, the cereal grasses also elongate their stems, called culms, due to activity within detached intercalary meristems which displaces cells upward, elevating the grain-bearing inflorescence. Variation in culm length within species is especially relevant to cereal crops, as demonstrated by the high-yielding semi-dwarfed cereals of the Green Revolution. Although previously understudied, recent renewed interest the regulation of subapical and intercalary growth suggests that control of cell division planes, boundary formation and temporal dynamics of differentiation, are likely critical mechanisms coordinating axial growth and development in plants.

Sorghum stem aerenchyma formation is regulated by SbNAC_D during internode development

Sorghum bicolor is a drought-resilient C4 grass used for production of grain, forage, sugar, and biomass. Sorghum genotypes capable of accumulating high levels of stem sucrose have solid stems that contain low levels of aerenchyma. The D-locus on SBI06 modulates the extent of aerenchyma formation in sorghum stems and leaf midribs. A QTL aligned with this locus was identified and fine-mapped in populations derived from BTx623*IS320c, BTx623*R07007, and BTx623*Standard broomcorn. Analysis of coding polymorphisms in the fine-mapped D-locus showed that genotypes that accumulate low levels of aerenchyma encode a truncated NAC transcription factor (Sobic.006G147400, SbNAC_d1), whereas parental lines that accumulate higher levels of stem aerenchyma encode full-length NAC TFs (SbNAC-D). During vegetative stem development, aerenchyma levels are low in nonelongated stem internodes, internode growing zones, and nodes. Aerenchyma levels increase in recently elongated internodes starting at the top of the internode near the center of the stem. SbNAC_D was expressed at low levels in nonelongated internodes and internode growing zones and at higher levels in regions of stem internodes that form aerenchyma. SbXCP1, a gene encoding a cysteine protease involved in programmed cell death, was induced in SbNAC_D genotypes in parallel with aerenchyma formation in sorghum stems but not in SbNAC_d1 genotypes. Several sweet sorghum genotypes encode the recessive SbNAC_d1 allele and have low levels of stem aerenchyma. Based on these results, we propose that SbNAC_D is the D-gene identified by Hilton (1916) and that allelic variation in SbNAC_D modulates the extent of aerenchyma formation in sorghum stems.

Current Models for Transcriptional Regulation of Secondary Cell Wall Biosynthesis in Grasses

Secondary cell walls mediate many crucial biological processes in plants including mechanical support, water and nutrient transport and stress management. They also provide an abundant resource of renewable feed, fiber, and fuel. The grass family contains the most important food, forage, and biofuel crops. Understanding the regulatory mechanism of secondary wall formation in grasses is necessary for exploiting these plants for agriculture and industry. Previous research has established a detailed model of the secondary wall regulatory network in the dicot model species Arabidopsis thaliana. Grasses, branching off from the dicot ancestor 140-150 million years ago, display distinct cell wall morphology and composition, suggesting potential for a different secondary wall regulation program from that established for dicots. Recently, combined application of molecular, genetic and bioinformatics approaches have revealed more transcription factors involved in secondary cell wall biosynthesis in grasses. Compared with the dicots, grasses exhibit a relatively conserved but nevertheless divergent transcriptional regulatory program to activate their secondary cell wall development and to coordinate secondary wall biosynthesis with other physiological processes.

Involvement of histone acetylation and deacetylation in regulating auxin responses and associated phenotypic changes in plants

The most recent outcomes about the transcription factors and transcription complexes mediated auxin signaling pathway by the histone acetylation and deacetylation. The phytohormone auxin, is required to regulate its accumulation spatiotemporally and responses to orchestrate various developmental levels in plants. Histone acetylation and deacetylation modulate auxin biosynthesis, its distribution and accumulation. In the absence of auxin, histone deacetylase represses the expression of auxin-responsive genes. Various transcription factors and transcription complexes facilitate the proper regulation of auxin signaling pathway genes. The primary and lateral root development, promotion of flowering and initiation of seed germination are all regulated by auxin-mediated histone acetylation and deacetylation. These findings conclude the auxin mode of action, which is mediated by histone acetylation and deacetylation, and associated phenotypic responses in plants, along with the underlying mechanism of these modifications.

A Growing Stem Inhibits Bud Outgrowth - The Overlooked Theory of Apical Dominance

Three theories of apical dominance, direct, diversion, and indirect, were proposed in the 1930s to explain how auxin synthesized in the shoot apex might inhibit axillary bud outgrowth, and thus shoot branching. The direct and diversion theories of apical dominance have been investigated in detail, and they are replaced with the current auxin transport canalization and second messenger theories, respectively. These two current theories still cannot entirely explain the phenomenon of apical dominance. Although there is ample evidence that the inhibition of bud outgrowth by auxin from the shoot apex is linked to stem elongation and highly branched auxin biosynthesis or signaling mutants are dwarf, the third theory proposed in the 1930s, the indirect theory, that explains apical dominance as auxin-induced stem growth indirectly inhibits bud outgrowth has been overlooked. The indirect theory did not propose how a growing stem might inhibit bud outgrowth. Recent discoveries indicate bud dormancy (syn. quiescence, paradormancy) in response to intrinsic and environmental factors in diverse species is linked to enhanced growth of the main shoot and reduced sugar level in the buds. Since a growing stem is a strong sink for sugars, and sugar is indispensable for shoot branching, the indirect theory of apical dominance might now be explained as auxin-induced stem growth inhibits bud outgrowth by diverting sugars away from buds. Detailed study of the indirect theory and the effect of source-sink status on dormancy and outgrowth of axillary buds will advance our knowledge of apical dominance and shoot branching in plants.