Metabolomic and transcriptomic basis of photoperiodic response regulation in broomcorn millet (Panicum miliaceum L.).

To elucidate the mechanisms underlying photoperiodic responses, we investigated the genomic and metabolomic responses of two broomcorn millet (Panicum miliaceum L.) genotypes. For this purpose, light-insensitive (D32) and light-sensitive (M51) genotypes were exposed to a 16 h photoperiod (long-day (LD) conditions) and an 8 h photoperiod (short-day (SD) conditions), and various transcriptomic and metabolomic changes were investigated. A total of 1664, 2564, 13,017, and 15548 DEGs were identified in the SD-D, LD-D, LD-M, and SD-M groups, respectively. Furthermore, 112 common DEGs were identified as well. Interestingly, most DEGs in the different groups were associated with photosynthesis and phenylpropanoid and carotenoid biosynthesis. In addition, 822 metabolites were identified under different treatments. The main metabolites, including L-malic and fumaric acids, were identified in the negative mode, whereas brucine and loperamide were identified in the positive mode. KEGG analysis revealed that the metabolites in the different groups were enriched in the same metabolic pathway of the TCA cycle. Furthermore, in negative mode, the metabolites of M51 were mainly D-glucose, whereas those of D32 were mainly L-malic and fumaric acids. One photoperiod candidate gene (C2845_PM11G01290), annotated as ATP6B, significantly increased the levels of L-malic and fumaric acids. In conclusion, our study provides a theoretical basis for understanding the molecular mechanisms of photoperiodic response regulation and can be used as a reference for marker development and resource identification in Panicum miliaceum L..

Transcriptome-wide m6A methylation profile reveals tissue specific regulatory networks in switchgrass (Panicum virgatum L.) under cadmium stress.

The heavy metal cadmium (Cd), known for its high toxicity, poses a grave threat to human health through the food chain. N6-methyladenosine (m6A), the most abundant internal modification, regulates plant adaptation to various adversities, yet the panorama of m6A modifications in switchgrass under cadmium stress remains elusive. This study examines the physiological responses of switchgrass roots and shoots exposed to 50 µM CdCl2, alongside an overview of transcriptome-wide m6A methylation patterns. After cadmium treatment, methylation modifications are primarily enriched near stop codons and the 3'UTR region, with a negative correlation between m6A modification and gene expression levels. In shoots, approximately 58 % of DEGs with m6A modifications show upregulation in expression and decrease in m6A peaks, including zinc transporter 4-like (ZIP4). In roots, about 43 % of DEGs with m6A modifications exhibit downregulation in expression and increase in m6A peaks, such as the ABC transporter family member (ABCG25). We further validate the m6A enrichment, gene expression and mRNA stability of ZIP4 in response to Cd treatment. The results suggest that the negative correlation of m6A enrichment and gene expression is due to altered mRNA stability. Our study establishes an m6A regulatory network governing cadmium transport in switchgrass roots and shoots, offering new avenues for candidate gene manipulation in phytoremediation applications of heavy metal pollution.

Disentangling plant- and environment-mediated drivers of active rhizosphere bacterial community dynamics during short-term drought.

Mitigating the effects of climate stress on crops is important for global food security. The microbiome associated with plant roots, the rhizobiome, can harbor beneficial microbes that alleviate stress, but the factors influencing their recruitment are unclear. We conducted a greenhouse experiment using field soil with a legacy of growing switchgrass and common bean to investigate the impact of short-term drought severity on the recruitment of active bacterial rhizobiome members. We applied 16S rRNA and 16S rRNA gene sequencing for both crops and metabolite profiling for switchgrass. We included planted and unplanted conditions to distinguish environment- versus plant-mediated rhizobiome drivers. Differences in community structure were observed between crops and between drought and watered and planted and unplanted treatments within crops. Despite crop-specific communities, drought rhizobiome dynamics were similar across the two crops. The presence of a plant more strongly explained the rhizobiome variation in bean (17%) than in switchgrass (3%), with a small effect of plant mediation during drought observed only for the bean rhizobiome. The switchgrass rhizobiome was stable despite changes in rhizosphere metabolite profiles between planted and unplanted treatments. We conclude that rhizobiome responses to short-term drought are crop-specific, with possible decoupling of plant exudation from rhizobiome responses.

Genome-wide association study reveals marker-trait associations for major agronomic traits in proso millet (Panicum miliaceum L.).

MAIN CONCLUSION: The pilot-scale genome-wide association study in the US proso millet identified twenty marker-trait associations for five morpho-agronomic traits identifying genomic regions for future studies (e.g. molecular breeding and map-based cloning). Proso millet (Panicum miliaceum L.) is an ancient grain recognized for its excellent water-use efficiency and short growing season. It is an indispensable part of the winter wheat-based dryland cropping system in the High Plains of the USA. Its grains are endowed with high nutritional and health-promoting properties, making it increasingly popular in the global market for healthy grains. There is a dearth of genomic resources in proso millet for developing molecular tools to complement conventional breeding for developing high-yielding varieties. Genome-wide association study (GWAS) is a widely used method to dissect the genetics of complex traits. In this pilot study of the first-ever GWAS in the US proso millet, 71 globally diverse genotypes of 109 the US proso millet core collection were evaluated for five major morpho-agronomic traits at two locations in western Nebraska, and GWAS was conducted to identify single nucleotide polymorphisms (SNPs) associated with these traits. Analysis of variance showed that there was a significant difference among the genotypes, and all five traits were also found to be highly correlated with each other. Sequence reads from genotyping-by-sequencing (GBS) were used to identify 11,147 high-quality bi-allelic SNPs. Population structure analysis with those SNPs showed stratification within the core collection. The GWAS identified twenty marker-trait associations (MTAs) for the five traits. Twenty-nine putative candidate genes associated with the five traits were also identified. These genomic regions can be used to develop genetic markers for marker-assisted selection in proso millet breeding.

Improvement of little millet (Panicum sumatrense) using novel omics platform and genetic resource integration.

MAIN CONCLUSION: This article explores possible future initiatives, such as the development of targeted breeding and integrated omics approach to boost little millet production, nutritional value, and environmental adaptation. Little millet (P. sumatrense) is a staple grain in many parts of Asia and Africa owing to its abundance in vitamins and minerals and its ability to withstand harsh agro-ecological conditions. Enhancing little millet using natural resources and novel crop improvement strategy is an effective way of boosting nutritional and food security. To understand the genetic makeup of the crop and figure out important characteristics linked to nutritional value, biotic and abiotic resistance, and production, researchers in this field are currently resorting on genomic technology. These realizations have expedited the crop's response to shifting environmental conditions by enabling the production of superior cultivars through targeted breeding. Going forward, further improvements in breeding techniques and genetics may boost the resilience, nutritional content, and production of little millet, which would benefit growers and consumers alike. The research and development on little millet improvement using novel omics platform and the integration of genetic resources are summarized in this review paper. Improved cultivars of little millet that satisfy changing farmer and consumer demands have already been developed through the use of these novel breeding strategies. This article also explores possible future initiatives, such as the development of targeted breeding, genomics, and sustainable agriculture methods. The potential for these measures to boost little millet's overall production, nutritional value, and climate adaptation will be extremely helpful in addressing nutritional security.

A complete reference genome of broomcorn millet.

Broomcorn millet (Panicum miliaceum L.), known for its traits of drought resistance, adaptability to poor soil, short growth period, and high photosynthetic efficiency as a C4 plant, represents one of the earliest domesticated crops globally. This study reports the telomere-to-telomere (T2T) gap-free reference genome for broomcorn millet (AJ8) using PacBio high-fidelity (HiFi) long reads, Oxford Nanopore long-read technologies and high-throughput chromosome conformation capture (Hi-C) sequencing data. The size of AJ8 genome was approximately 834.7 Mb, anchored onto 18 pseudo-chromosomes. Notably, 18 centromeres and 36 telomeres were obtained. The assembled genome showed high quality in terms of completeness (BUSCO score: 99.6%, QV: 61.7, LAI value: 20.4). In addition, 63,678 protein-coding genes and 433.8 Mb (~52.0%) repetitive sequences were identified. The complete reference genome for broomcorn millet provides a valuable resource for genetic studies and breeding of this important cereal crop.

Establishment of genome-editing system and assembly of a near-complete genome in broomcorn millet.

The ancient crop broomcorn millet (Panicum miliaceum L.) is an indispensable orphan crop in semi-arid regions due to its short life cycle and excellent abiotic stress tolerance. These advantages make it an important alternative crop to increase food security and achieve the goal of zero hunger, particularly in light of the uncertainty of global climate change. However, functional genomic and biotechnological research in broomcorn millet has been hampered due to a lack of genetic tools such as transformation and genome-editing techniques. Here, we successfully performed genome editing of broomcorn millet. We identified an elite variety, Hongmi, that produces embryogenic callus and has high shoot regeneration ability in in vitro culture. We established an Agrobacterium tumefaciens-mediated genetic transformation protocol and a clustered regularly interspaced short palindromic repeats (CRISPR)/Cas9-mediated genome-editing system for Hongmi. Using these techniques, we produced herbicide-resistant transgenic plants and edited phytoene desaturase (PmPDS), which is involved in chlorophyll biosynthesis. To facilitate the rapid adoption of Hongmi as a model line for broomcorn millet research, we assembled a near-complete genome sequence of Hongmi and comprehensively annotated its genome. Together, our results open the door to improving broomcorn millet using biotechnology.

Application of CRISPR/Cas12i.3 for targeted mutagenesis in broomcorn millet (Panicum miliaceum L.).

A CRISPR/Cas12i.3-based gene editing platform is established in broomcorn millet (Panicum miliaceum) and used to create new elite germplasm for this ancient crop.

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

PREMISE: Leaf tensile resistance, a leaf's ability to withstand pulling forces, is an important determinant of plant ecological strategies. One potential driver of leaf tensile resistance is growing season length. When growing seasons are long, strong leaves, which often require more time and resources to construct than weak leaves, may be more advantageous than when growing seasons are short. Growing season length and other ecological conditions may also impact the morphological traits that underlie leaf tensile resistance. METHODS: To understand variation in leaf tensile resistance, we measured size-dependent leaf strength and size-independent leaf toughness in diverse genotypes of the widespread perennial grass Panicum virgatum (switchgrass) in a common garden. We then used quantitative genetic approaches to estimate the heritability of leaf tensile resistance and whether there were genetic correlations between leaf tensile resistance and other morphological traits. RESULTS: Leaf tensile resistance was positively associated with aboveground biomass (a proxy for fitness). Moreover, both measures of leaf tensile resistance exhibited high heritability and were positively genetically correlated with leaf lamina thickness and leaf mass per area (LMA). Leaf tensile resistance also increased with the growing season length in the habitat of origin, and this effect was mediated by both LMA and leaf thickness. CONCLUSIONS: Differences in growing season length may promote selection for different leaf lifespans and may explain existing variation in leaf tensile resistance in P. virgatum. In addition, the high heritability of leaf tensile resistance suggests that P. virgatum will be able to respond to climate change as growing seasons lengthen.

UDP-glycosyltransferase UGT96C10 functions as a novel detoxification factor for conjugating the activated dinitrotoluene sulfonate in switchgrass.

Dinitrotoluene sulfonates (DNTSes) are highly toxic hazards regulated by the Resource Conservation and Recovery Act (RCRA) in the United States. The trinitrotoluene (TNT) red water formed during the TNT purification process consists mainly of DNTSes. Certain plants, including switchgrass, reed and alfalfa, can detoxify low concentrations of DNTS in TNT red water-contaminated soils. However, the precise mechanism by which these plants detoxify DNTS remains unknown. In order to aid in the development of phytoremediation resources with high DNTS removal rates, we identified and characterized 1-hydroxymethyl-2,4-dinitrobenzene sulfonic acid (HMDNBS) and its glycosylated product HMDNBS O-glucoside as the degradation products of 2,4-DNT-3-SO3Na, the major isoform of DNTS in TNT red water-contaminated soils, in switchgrass via LC-MS/MS- and NMR-based metabolite analyses. Transcriptomic analysis revealed that 15 UDP-glycosyltransferase genes were dramatically upregulated in switchgrass plants following 2,4-DNT-3-SO3Na treatment. We expressed, purified and assayed the activity of recombinant UGT proteins in vitro and identified PvUGT96C10 as the enzyme responsible for the glycosylation of HMDNBS in switchgrass. Overexpression of PvUGT96C10 in switchgrass significantly alleviated 2,4-DNT-3-SO3Na-induced plant growth inhibition. Notably, PvUGT96C10-overexpressing transgenic switchgrass plants removed 83.1% of 2,4-DNT-3-SO3Na in liquid medium after 28 days, representing a 3.2-fold higher removal rate than that of control plants. This work clarifies the DNTS detoxification mechanism in plants for the first time, suggesting that PvUGT96C10 is crucial for DNTS degradation. Our results indicate that PvUGT96C10-overexpressing plants may hold great potential for the phytoremediation of TNT red water-contaminated soils.

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

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

Genome-wide identification of GRF transcription factors and their expression profile in stem meristem of broomcorn millet (Panicum miliaceum L.).

To understand the genome-wide information of the GRF family genes in broomcorn millet and their expression profile in the vegetative meristems, bioinformatic methods and transcriptome sequencing were used to analyze the characteristics, physical and chemical properties, phylogenetic relationship, chromosome distribution, gene structure, cis-acting elements and expression profile in stem meristem for the GRF family members. The results showed that the GRF gene family of millet contains 21 members, and the PmGRF gene is unevenly distributed on 12 chromosomes. The lengths of PmGRF proteins vary from 224 to 618 amino acids, and the isoelectric points are between 4.93-9.69. Each member of the family has 1-4 introns and 2-5 exons. The protein PmGRF13 is localized in both the nucleus and chloroplast, and the rest PmGRF proteins are located in the nucleus. Phylogenetic analysis showed that the 21 GRF genes were divided into 4 subfamilies (A,B,C and D) in broomcorn millet. The analysis of cis-acting elements showed that there were many cis-acting elements involved in light response, hormone response, drought induction, low temperature response and other environmental stress responses in the 2000 bp sequence upstream of the GRF genes. Transcriptome sequencing and qRT-PCR analyses showed that the expression levels of PmGRF3 and PmGRF12 in the dwarf variety Zhang778 were significantly higher than those of the tall variety Longmi12 in the internode and node meristems at the jointing stage, while the expression patterns of PmGRF4, PmGRF16 and PmGRF21 were reverse. In addition, the expression levels of PmGRF2 and PmGRF5 in the internode of Zhang778 were significantly higher than Longmi12. The other GRF genes were not or insignificantly expressed. These results indicated that seven genes, PmGRF2, PmGRF3, PmGRF4, PmGRF5, PmGRF12, PmGRF16 and PmGRF21, were related to the formation of plant height in broomcorn millet.

Genomic variation in weedy and cultivated broomcorn millet accessions uncovers the genetic architecture of agronomic traits.

Large-scale genomic variations are fundamental resources for crop genetics and breeding. Here we sequenced 1,904 genomes of broomcorn millet to an average of 40× sequencing depth and constructed a comprehensive variation map of weedy and cultivated accessions. Being one of the oldest cultivated crops, broomcorn millet has extremely low nucleotide diversity and remarkably rapid decay of linkage disequilibrium. Genome-wide association studies identified 186 loci for 12 agronomic traits. Many causative candidate genes, such as PmGW8 for grain size and PmLG1 for panicle shape, showed strong selection signatures during domestication. Weedy accessions contained many beneficial variations for the grain traits that are largely lost in cultivated accessions. Weedy and cultivated broomcorn millet have adopted different loci controlling flowering time for regional adaptation in parallel. Our study uncovers the unique population genomic features of broomcorn millet and provides an agronomically important resource for cereal crops.

Nutri-cereal tissue-specific transcriptome atlas during development: Functional integration of gene expression to identify mineral uptake pathways in little millet (Panicum sumatrense).

Little millet (Panicum sumatrense Roth ex Roem. & Schult.) is an essential minor millet of southeast Asia and Africa's temperate and subtropical regions. The plant is stress-tolerant, has a short life cycle, and has a mineral-rich nutritional profile associated with unique health benefits. We report the developmental gene expression atlas of little millet (genotype JK-8) from ten tissues representing different stages of its life cycle, starting from seed germination and vegetative growth to panicle maturation. The developmental transcriptome atlas led to the identification of 342 827 transcripts. The BUSCO analysis and comparison with the transcriptomes of related species confirm that this study presents high-quality, in-depth coverage of the little millet transcriptome. In addition, the eFP browser generated here has a user-friendly interface, allowing interactive visualizations of tissue-specific gene expression. Using these data, we identified transcripts, the orthologs of which in Arabidopsis and rice are involved in nutrient acquisition, transport, and response pathways. The comparative analysis of the expression levels of these transcripts holds great potential for enhancing the mineral content in crops, particularly zinc and iron, to address the issue of "hidden hunger" and to attain nutritional security, making it a valuable asset for translational research.

Genetic linkage map construction and QTL analysis for plant height in proso millet (Panicum miliaceum L.).

KEY MESSAGE: A genetic linkage map representing proso millet genome was constructed with SSR markers, and a major QTL corresponding to plant height was mapped on chromosome 14 of this map. Proso millet (Panicum miliaceum L.) has the lowest water requirements of all cultivated cereal crops. However, the lack of a genetic map and the paucity of genomic resources for this species have limited the utility of proso millet for detailed genetic studies and hampered genetic improvement programs. In this study, 97,317 simple sequence repeat (SSR) markers were developed based on the genome sequence of the proso millet landrace Longmi 4. Using some of these markers in conjunction with previously identified SSRs, an SSR-based linkage map for proso millet was successfully constructed using a large mapping population (316 F2 offspring). In total, 186 SSR markers were assigned to 18 linkage groups corresponding to the haploid chromosomes. The constructed map had a total length of 3033.42 centimorgan (cM) covering 78.17% of the assembled reference genome. The length of the 18 linkage groups ranged from 88.89 cM (Chr. 15) to 274.82 cM (Chr. 16), with an average size of 168.17 cM. To our knowledge, this is the first genetic linkage map for proso millet based on SSR markers. Plant height is one of the most important traits in crop improvement. A major QTL was repeatedly detected in different environments, explaining 8.70-24.50% of the plant height variations. A candidate gene affecting auxin biosynthesis and transport, and ROS homeostasis regulation was predicted. Thus, the linkage map and QTL analysis provided herein will promote the development of gene mining and molecular breeding in proso millet.

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

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

Blocking miR528 function promotes tillering and regrowth in switchgrass.

MiRNAs have been reported to be the key regulators involving a wide range of biological processes in diverse plant species, but their functions in switchgrass, an important biofuel and forage crop, are largely unknown. Here, we reported the novel function of miR528, which has expanded to four copies in switchgrass, in controlling biomass trait of tillering number and regrowth rate after mowing. Blocking miR528 activity by expressing short tandem target mimic (STTM) increased tiller number and regrowth rate after mowing. The quadruple pvmir528 mutant lines derived from genome editing also showed such improved traits. Degradome and RNA-seq analysis, combined with in situ hybridization assay revealed that up-regulation of two miR528 targets coding for Cu/Zn-SOD enzymes, might be responsible for the improved traits of tillering and regrowth in pvmir528 mutant. Additionally, natural variations in the miR528-SOD interaction exist in C3 and C4 monocot species, implying the distinct regulatory strength of the miR528-SOD module during monocot evolution. Overall, our data illuminated a novel role of miR528 in controlling biomass traits and provided a new target for genetic manipulation-mediated crop improvement.

Transcriptomic analysis of ncRNAs and mRNAs interactions during drought stress in switchgrass.

Switchgrass (Panicum virgatum L.) plays a pivotal role as a bioenergy feedstock in the production of cellulosic ethanol and contributes significantly to enhancing ecological grasslands and soil quality. The utilization of non-coding RNAs (ncRNAs) has gained momentum in deciphering the intricate genetic responses to abiotic stress in various plant species. Nevertheless, the current research landscape lacks a comprehensive exploration of the responses of diverse ncRNAs, including long non-coding RNAs (lncRNAs), circular RNAs (circRNAs), and microRNAs (miRNAs), to drought stress in switchgrass. In this study, we employed whole transcriptome sequencing to comprehensively characterize the expression profiles of both mRNA and ncRNAs during episodes of drought stress in switchgrass. Our analysis identified a total of 12,511 mRNAs, 59 miRNAs, 38 circRNAs, and 368 lncRNAs that exhibited significant differential expression between normal and drought-treated switchgrass leaves. Notably, the majority of up-regulated mRNAs displayed pronounced enrichment within the starch and sucrose metabolism pathway, as validated through KEGG analysis. Co-expression analysis illuminated that differentially expressed (DE) lncRNAs conceivably regulated 1308 protein-coding genes in trans and 7110 protein-coding genes in cis. Furthermore, both cis- and trans-target mRNAs of DE lncRNAs exhibited enrichment in four common KEGG pathways. The intricate interplay between lncRNAs and circRNAs with miRNAs via miRNA response elements was explored within the competitive endogenous RNA (ceRNA) network framework. As a result, we constructed elaborate regulatory networks, including lncRNA-novel_miRNA480-mRNA, lncRNA-novel_miRNA304-mRNA, lncRNA/circRNA-novel_miRNA122-PvSS4, and lncRNA/circRNA-novel_miRNA14-PvSS4, and subsequently validated the functionality of the target gene, starch synthase 4 (PvSS4). Furthermore, through the overexpression of PvSS4, we ascertained its capacity to enhance drought tolerance in yeast. However, it is noteworthy that PvSS4 did not exhibit any discernible impact under salt stress conditions. These findings, as presented herein, not only contribute substantively to our understanding of ceRNA networks but also offer a basis for further investigations into their potential functions in response to drought stress in switchgrass.

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

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

Dynamic Reconfiguration of Switchgrass Proteomes in Response to Rust (Puccinia novopanici) Infection.

Switchgrass (Panicum virgatum L.) can be infected by the rust pathogen (Puccinia novopanici) and results in lowering biomass yields and quality. Label-free quantitative proteomics was conducted on leaf extracts harvested from non-infected and infected plants from a susceptible cultivar (Summer) at 7, 11, and 18 days after inoculation (DAI) to follow the progression of disease and evaluate any plant compensatory mechanisms to infection. Some pustules were evident at 7 DAI, and their numbers increased with time. However, fungal DNA loads did not appreciably change over the course of this experiment in the infected plants. In total, 3830 proteins were identified at 1% false discovery rate, with 3632 mapped to the switchgrass proteome and 198 proteins mapped to different Puccinia proteomes. Across all comparisons, 1825 differentially accumulated switchgrass proteins were identified and subjected to a STRING analysis using Arabidopsis (A. thaliana L.) orthologs to deduce switchgrass cellular pathways impacted by rust infection. Proteins associated with plastid functions and primary metabolism were diminished in infected Summer plants at all harvest dates, whereas proteins associated with immunity, chaperone functions, and phenylpropanoid biosynthesis were significantly enriched. At 18 DAI, 1105 and 151 proteins were significantly enriched or diminished, respectively. Many of the enriched proteins were associated with mitigation of cellular stress and defense.