ABSTRACT Small signaling peptides (SSPs) are important regulators of plant growth, development, and responses to biotic and abiotic stress, yet their role in the C4 grass Sorghum bicolor is largely uncharacterized. To help fill this knowledge gap, 219 sorghum genes that encode SSPs were identified based on SSP sequences previously identified in Arabidopsis thaliana, Zea mays, Oryza sativa, Triticum aestivum, and Brachypodium distachyon. The 219 sorghum SSP‐encoding genes were assigned to 19 gene families, analyzed for the presence of motifs, and aligned with genes that encode SSPs in other plants using phylogenetic analysis. Sorghum genes in 12 of the 19 SSP gene families had not been previously characterized. Expression of the 219 SSP‐encoding genes in sorghum organs, during stem development, and in stem tissues and cell types revealed distinct spatial, temporal, and developmental patterns of expression. Genes associated with the SbCEP and SbRGF families were preferentially expressed in roots, whereas SbEPF genes were expressed in stem epidermal and pith parenchyma cells and panicles. The expression of genes during bioenergy sorghum stem growth and development was investigated because stems account for ~80% of harvested biomass and serve as conduits for water and nutrient transport between leaves and roots. During stem development, 28 SSP genes in several families (CLE, EPF, CEP, GASS, PSY, ES, PSK, CAPE, POE) were expressed at higher levels in zones of cell proliferation. For example, the TDIF homologs SbCLE41 and SbCLE42 were expressed at high levels in nascent stem nodes where they may regulate vascular bundle cambial activity and cell differentiation. A different set of 15 genes in the CIF, POE, CAPE, PSY, CEP, RALF, and CLE families were expressed at higher levels in zones of stem tissue differentiation highlighted by elevated expression of five SbRALFRs in the stem nodal plexus. Cell type–specific expression of many sorghum genes that encode SSPs was observed in fully elongated internodes indicating gene expression is regulated with high spatial resolution. Overall, the results provide a foundation of information for analysis of SSP function in sorghum that can be integrated with knowledge of sorghum gene regulatory networks to modulate traits important for production of sorghum crops.

Although the green revolution adapted a handful of crops to homogeneous and high-input industrialized agriculture, much of the global population still relies on the local production of variable crop cultivars by low-input smallholder farms. This diversity of unhomogenized crops1, like that of the grain and bioenergy crop sorghum2-5, offers raw materials for genetic gain and cultivar improvement. However, breeding efforts can be constrained by highly specialized traits and breeding targets6. Here, to bridge this diversity, we constructed a 33-member pangenome reference and a diversity panel across 1,984 cultivars and landraces. We leveraged these resources to explore the complex interplay among historical contingency, ongoing adaptation and previously uncharacterized structural diversity. Specifically, our analyses conclusively demonstrated multiple nested and deeply diverged structural variants in the domestication gene SHATTERING1, which distinguish the previously established multicentric origin of sorghum. We then applied landscape genomics to reveal how gene flow and secondary contact created the complex genetic mosaic in contemporary breeding networks. As proof of concept for pangenome-accelerated trait discovery, we connected biosynthetic gene cluster structural variation to phenotypic leaf concentration of the cyanogenic glucoside dhurrin. Combined, these approaches will accelerate breeding and trait discovery and provide a framework for similar applications in other crops.

While the green revolution adapted a handful of crops to homogenous and high-input industrialized agriculture, much of the global population still relies on local food production from low-input smallholder farms that grow highly variable crop cultivars. The high diversity of the grain and bioenergy crop sorghum [1][1]–[4][2], and many other crops that were not homogenized during the green revolution [5][3], not only provides the raw materials for breeders to make substantial gains in cultivar improvement, but also constrains breeding efforts due to highly specialized locally adapted plant phenotypes [6][4]. Here, we construct a 33-member pangenome and identify trait-associated variants in 1,988 cultivars and landraces. We then apply these resources to explore the complex interplay between historical contingency, ongoing adaptation, and the potential for future gains through climate-aware genome-enabled breeding. Specifically, our analyses conclusively demonstrate that multiple nested, deeply diverged, and previously uncharacterized structural variants in the domestication gene SHATTERING1 distinguish the previously established multicentric origin of sorghum. We then apply landscape genomics tests to reveal how gene flow, adaptation, and secondary contact created the complex genetic mosaic in current global breeding networks. Further analysis of climate-gene associations highlights candidate loci underlying adaptation, including the biosynthetic gene cluster for the cyanogenic glucoside dhurrin. Combined, the pangenome-informed variants developed here will enable both trait discovery and subsequent marker assays to accelerate breeding and provide a framework for similar applications in other diverse and non-model crops. ### Competing Interest Statement The authors have declared no competing interest. United States Department of Energy, DE-AC02-05CH11231 [1]: #ref-1 [2]: #ref-4 [3]: #ref-5 [4]: #ref-6

Sorghum bicolor is a deep-rooted, heat- and drought-tolerant crop that thrives on marginal lands and is increasingly valued for its applications in biofuel, bioenergy, and biopolymer production. The sorghum stem, which can reach 4–5 m in length, serves as the primary reservoir of both lignocellulosic biomass and soluble sugars, making it a promising bioenergy feedstock. Although recent advances in genetic, genomic, and transcriptomic resources have improved our understanding of sorghum biology, comprehensive genome-wide analyses of functional dynamics across diverse organ types and developmental stages remain limited. In particular, candidate genes with stem preferred expression pattern or their associated cis-regulatory elements, which may program key stem-related functions and enable organ- or tissue-specific engineering, have not yet been identified. To address these gaps, we reanalyzed a published RNA-seq dataset to identify genes with organ-preferential expression and to infer representative organ functions across major developmental stages. Our analysis revealed that the sorghum stem exhibits distinct temporal functional signatures, which correlate with the developmental dynamics of stem-specific genes and their associated regulatory elements. We further identified a set of genes with ubiquitous stem-specific expression across diverse sorghum genotypes, suggesting their universal importance and broad potential for genetic engineering applications. Among them, SbTALE03 and SbTALE04 emerged as stem hub transcription factors (TF). Both genes were empirically validated for their stem specificity across stages. Gene regulatory network analysis further indicated that these TFs participate in stage-specific transcriptional programs that maintain and regulate stem development. This study presents a genome-wide analysis of organ-specific gene expression, functions, and regulatory networks in sorghum, with a focus on genes preferentially-expressed in stems and their promoter motifs. We identified a set of core stem-specific genes with ubiquitous expression across genotypes and developmental stages, including two experimentally validated transcription factors with potential roles in stem development. These findings offer valuable candidates for further functional characterization and genetic engineering aimed at improving sorghum stem biomass and composition.

Cytochrome P450s of the CYP79 family catalyze two N-hydroxylation reactions, converting a selected number of amino acids into the corresponding oximes. The sorghum genome (Sorghum bicolor) harbours nine CYP79A encoding genes, and here sequence comparisons of the CYP79As along with their substrate recognition sites (SRSs) are provided. The substrate specificity of previously uncharacterized CYP79As was investigated by transient expression in Nicotiana benthamiana and subsequent transformation of the oximes formed into the corresponding stable oxime glucosides catalyzed by endogenous UDPG-glucosyltransferases (UGTs). CYP79A61 uses phenylalanine as a substrate, whereas CYP79A91, CYP79A93, and CYP79A95 use valine and isoleucine as substrates, with CYP79A93 showing the ability also to use phenylalanine. CYP79A94 uses isoleucine as a substrate. Analysis of 249 sorghum transcriptomes from two different sorghum cultivars showed the expression levels and tissue-specific expression of the CYP79As. CYP79A1 is the committed gene in dhurrin formation and was the highest expressed gene in most tissues/organs. CYP79A61 was primarily expressed in fully developed leaf blades and leaf sheaths. CYP79A91 and CYP79A92 were expressed mainly in roots >200 cm below ground, while CYP79A93 and CYP79A94 were most highly expressed in the leaf collar and leaf sheath, respectively. The possible signalling effects of the oximes and their metabolites produced in different sorghum tissues are discussed.

Bioenergy sorghum is a highly productive drought tolerant C4 grass that accumulates 80

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

The colonial green microalga Botryococcus braunii is well known for producing liquid hydrocarbons that can be utilized as biofuel feedstocks. B. braunii is taxonomically classified as a single species made up of three chemical races, A, B, and L, that are mainly distinguished by the hydrocarbons produced. We previously reported a B race draft nuclear genome, and here we report the draft nuclear genomes for the A and L races. A comparative genomic study of the three B. braunii races and 14 other algal species within Chlorophyta revealed significant differences in the genomes of each race of B. braunii. Phylogenomically, there was a clear divergence of the three races with the A race diverging earlier than both the B and L races, and the B and L races diverging from a later common ancestor not shared by the A race. DNA repeat content analysis suggested the B race had more repeat content than the A or L races. Orthogroup analysis revealed the B. braunii races displayed more gene orthogroup diversity than three closely related Chlamydomonas species, with nearly 24-36% of all genes in each B. braunii race being specific to each race. This analysis suggests the three races are distinct species based on sufficient differences in their respective genomes. We propose reclassification of the three chemical races to the following species names: Botryococcus alkenealis (A race), Botryococcus braunii (B race), and Botryococcus lycopadienor (L race).

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.

ABSTRACT A collection of 47 bacteria isolated from the mucilage of aerial roots of energy sorghum is available at the Great Lakes Bioenergy Research Center, Michigan State University, Michigan, USA. We enriched bacteria with putative plant-beneficial phenotypes and included information on phenotypic diversity, taxonomy, and whole genome sequences.

Gene functional descriptions offer a crucial line of evidence for candidate genes underlying trait variation. Conversely, plant responses to environmental cues represent important resources to decipher gene function and subsequently provide molecular targets for plant improvement through gene editing. However, biological roles of large proportions of genes across the plant phylogeny are poorly annotated. Here we describe the Joint Genome Institute (JGI) Plant Gene Atlas, an update able data resource consisting of transcript abundance assays spanning 18 diverse species. To integrate across these diverse genotypes, we anal yzed e xpression pr ofiles, b uilt gene c lusters that exhibited tissue / condition specific expression, and tested for transcriptional response to environmental queues. We discovered extensive phylogenetically constrained and condition-specific expression profiles f or genes without an y previously documented functional annotation. Such conserved expression patterns and tightly co-expressed gene clusters let us assign expression derived additional biological information to 64 495 genes with otherwise unknown functions. The ever-expanding Gene Atlas resource is available at JGI Plant Gene Atlas ( https://plantgeneatlas.jgi.doe.gov ) and Phytozome ( https://phytozome .jgi.doe .gov/), providing bulk access to data and user-specified queries of gene sets. Combined, these web interfaces let users access differentiall y e xpressed genes, track orthologs across the Gene Atlas plants, graphically represent co-expressed genes, and visualize gene ontology and pathway enrichments.

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 stem cell-type gene expression and regulation was available to enable engineering. To obtain this information, Laser Capture Microdissection (LCM) 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 TF families contribute to distinct regulatory landscapes, where regulation is organized through various modes and identifiable network motifs. Cell-specific transcriptome data was combined with a stem developmental transcriptome dataset to identify the GRN that differentially activates the secondary cell wall (SCW) formation in stem xylem sclerenchyma and epidermal cells. The cell-type transcriptomic dataset provides a valuable source of information about the function of sorghum stem cell types and GRNs that will enable the engineering of bioenergy sorghum stems.### Competing Interest StatementThe authors have declared no competing interest.

Phyllosphere exudates create specialized microhabitats that shape microbial community diversity. We explored the microbiome associated with two sorghum phyllosphere exudates, the epicuticular wax and aerial root mucilage. We assessed the microbiome associated with the wax from sorghum plants over two growth stages, and the root mucilage additionally from nitrogen-fertilized and nonfertilized plants. In parallel, we isolated and characterized hundreds of bacteria from wax and mucilage, and integrated data from cultivation-independent and cultivation-dependent approaches to gain insights into exudate diversity and bacterial phenotypes. We found that Sphingomonadaceae and Rhizobiaceae families were the major taxa in the wax regardless of water availability and plant developmental stage to plants. The cultivation-independent mucilage-associated bacterial microbiome contained the families Erwiniaceae, Flavobacteriaceae, Rhizobiaceae, Pseudomonadaceae, and Sphingomonadaceae, and its structure was strongly influenced by sorghum development but only modestly influenced by fertilization. In contrast, the fungal community structure of mucilage was strongly affected by the year of sampling but not by fertilization or plant developmental stage, suggesting a decoupling of fungal–bacterial dynamics in the mucilage. Our bacterial isolate collection from wax and mucilage had several isolates that matched 100% to detected amplicon sequence variants, and were enriched on media that selected for phenotypes that included phosphate solubilization, putative diazotrophy, resistance to desiccation, capability to grow on methanol as a carbon source, and ability to grow in the presence of linalool and β-caryophyllene (terpenes in sorghum wax). This work expands our understanding of the microbiome of phyllosphere exudates and supports our long-term goal to translate microbiome research to support sorghum cultivation.

Mixed-linkage glucan (MLG) is a component of the cell wall (CW) of grasses and is composed of glucose monomers linked by β-1,3 and β-1,4 bonds. MLG is believed to have several biological functions, such as the mobilizable storage of carbohydrates and structural support of the CW. The extracellular levels of MLG are largely controlled by rates of synthesis mediated by cellulose synthase-like (CSL) enzymes, and turnover by lichenases. Economically important crops like sorghum accumulate MLG to variable levels during development. While in sorghum, like other grasses, there is one major MLG synthase (CSLF6), the identity of lichenases is yet unknown. To fill this gap, we identified three sorghum lichenases (SbLCH1-3) and characterized them in leaves in relation to the expression of SbCSLF6, and the abundance of MLG and starch. We established that SbLCH1-3 are secreted to the apoplast, consistent with a role of degrading MLG extracellularly. Furthermore, while SbCSLF6 expression was associated with cell development, the SbLCH genes exhibited distinct development-, cell-type-specific and diel-regulated expression. Therefore, our study identifies three functional sorghum MLG lichenases and highlights that MLG accumulation in sorghum leaves is likely controlled by the activity of lichenases that tune MLG levels, possibly to suit distinct cell and developmental needs in planta. These findings have important implications for improving the growth, yield, and composition of sorghum as a feedstock.

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.

Texas A&M University recently completed a set of Automated Precision Phenotyping (APP) Greenhouses that incorporate robotic systems for automated collection of advanced sensor-based plant phenotypes. Transiting the length of a greenhouse is a gantry beam, on which a rolling truck provides a second axis of motion along the gantry. Attached to the truck is a 3.0-m long robotic arm that is controlled to position a sensor head at virtually any position relative to any plant in a greenhouse. The robotic arm can be programmed to operate quickly and safely in complicated scanning patterns to enable data collection on all plants in the greenhouse within a time window of a few hours, ensuring consistent conditions during data collection. The sensor head includes a high-speed multispectral camera and eventually a Raman spectrometer. Relative to phenotyping greenhouses at other institutions, the APP Greenhouses have the advantage of maximum flexibility in configuration of plants in the greenhouses, in positioning of sensors relative to the plants, and in the types of sensors used, making research capabilities in the APP Greenhouses truly unique. Preliminary data have been collected on sorghum and maize plants. Four-band multispectral images have been collected daily, scanning the side of each plant from top to bottom. Preliminary software development is directed at automated image stitching to create a full side-view image of each plant, from which consistent metrics can be automatically calculated, such as plant height, stalk diameter, leaf angle, etc.

Bioenergy sorghum hybrids are being developed with enhanced drought tolerance and high levels of stem sugars. Raffinose family oligosaccharides (RFOs) contribute to plant environmental stress tolerance, sugar storage, transport, and signaling. To better understand the role of RFOs in sorghum, genes involved in myo-inositol and RFO metabolism were identified and relative transcript abundance analyzed during development. Genes involved in RFO biosynthesis (SbMIPS1, SbInsPase, SbGolS1, SbRS) were more highly expressed in leaves compared to stems and roots, with peak expression early in the morning in leaves. SbGolS, SbRS, SbAGA1 and SbAGA2 were also expressed at high levels in the leaf collar and leaf sheath. In leaf blades, genes involved in myo-inositol biosynthesis (SbMIPS1, SbInsPase) were expressed in bundle sheath cells, whereas genes involved in galactinol and raffinose synthesis (SbGolS1, SbRS) were expressed in mesophyll cells. Furthermore, SbAGA1 and SbAGA2, genes that encode neutral-alkaline alpha-galactosidases that hydrolyze raffinose, were differentially expressed in minor vein bundle sheath cells and major vein and mid-rib vascular and xylem parenchyma. This suggests that raffinose synthesized from sucrose and galactinol in mesophyll cells diffuses into vascular bundles where hydrolysis releases sucrose for long distance phloem transport. Increased expression (>20-fold) of SbAGA1 and SbAGA2 in stem storage pith parenchyma of sweet sorghum between floral initiation and grain maturity, and higher expression in sweet sorghum compared to grain sorghum, indicates these genes may play a key role in non-structural carbohydrate accumulation in stems.