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.

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.

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.

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.

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.

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.

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.

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.

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

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

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

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

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

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

Transcriptome and DNA methylome divergence of inflorescence development between 2 ecotypes in Panicum hallii.

The morphological diversity of the inflorescence determines flower and seed production, which is critical for plant adaptation. Hall's panicgrass (Panicum hallii, P. hallii) is a wild perennial grass that has been developed as a model to study perennial grass biology and adaptive evolution. Highly divergent inflorescences have evolved between the 2 major ecotypes in P. hallii, the upland ecotype (P. hallii var hallii, HAL2 genotype) with compact inflorescence and large seed and the lowland ecotype (P. hallii var filipes, FIL2 genotype) with an open inflorescence and small seed. Here we conducted a comparative analysis of the transcriptome and DNA methylome, an epigenetic mark that influences gene expression regulation, across different stages of inflorescence development using genomic references for each ecotype. Global transcriptome analysis of differentially expressed genes (DEGs) and co-expression modules underlying the inflorescence divergence revealed the potential role of cytokinin signaling in heterochronic changes. Comparing DNA methylome profiles revealed a remarkable level of differential DNA methylation associated with the evolution of P. hallii inflorescence. We found that a large proportion of differentially methylated regions (DMRs) were located in the flanking regulatory regions of genes. Intriguingly, we observed a substantial bias of CHH hypermethylation in the promoters of FIL2 genes. The integration of DEGs, DMRs, and Ka/Ks ratio results characterized the evolutionary features of DMR-associated DEGs that contribute to the divergence of the P. hallii inflorescence. This study provides insights into the transcriptome and epigenetic landscape of inflorescence divergence in P. hallii and a genomic resource for perennial grass biology.

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

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

Genome-wide characterization of LBD transcription factors in switchgrass (Panicum virgatum L.) and the involvement of PvLBD12 in salt tolerance.

KEY MESSAGE: PvLBD12 enhanced the salt tolerance by increasing proline accumulation, improving K+ accumulation, and decreasing reactive oxygen species level in switchgrass. Abiotic stresses are the serious factors which limit plant development and productivity and restrict the agricultural economy. It is important, therefore, to understand the mechanism of abiotic tolerance in plants. Lateral organ boundaries domain (LBD) proteins as plant-specific transcription factors play important function in plant lateral organ development, plant regeneration, and abiotic stress. In our study, we identify 69 LBD members from switchgrass genome-wide sequences and classify them based on their homology with LBD proteins in Arabidopsis. RT-qPCR showed that PvLBD genes had different expression patterns under abiotic stress conditions, indicating that they play important roles in various stress. PvLBD12 was selected as a candidate gene for further functional analysis because it had the highest expression level under salt stress. Overexpression of PvLBD12 enhanced salt tolerance by altering a wide range of physiological responses (like increased proline accumulation, reduced malondialdehyde production, improved K+ accumulation, and reduced Na+ absorption) in switchgrass. Some stress response genes such as proline biosynthesis gene PvP5CS1, vacuolar Na+(K+)/H+ antiporter gene PvNHX1, two key ROS-scavenging enzyme genes PvCAT and PvSOD were all upregulated in PvLBD12 overexpression lines. Taken together, PvLBD12 plays a pivotal role in response to salt stress by increasing proline accumulation, improving K+ accumulation, reducing Na+ absorption, and decreasing reactive oxygen species level. It will be better to understand the potential biological functions of LBD genes in other plants.

Genomic prediction of switchgrass winter survivorship across diverse lowland populations.

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

Transcriptome analysis reveals new insights in the starch biosynthesis of non-waxy and waxy broomcorn millet (Panicum miliaceum L.).

Broomcorn millet is a popular cereal with health benefits, and its grains are rich in starch. However, the differences in the pathway and key genes involved in starch biosynthesis of waxy and non-waxy broomcorn millet grain remain unclear. Therefore, the grain and starch physicochemical index and transcriptomic analyses of two genotypes of broomcorn millet were conducted at 3, 6, 9, 12, 15, 18, and 21 days after pollination. The phenotypic and physiological results indicated that the starch synthetic process of non-waxy and waxy broomcorn millet was significantly different. The amylose, amylopectin, and total starch contents of non-waxy broomcorn millet were 1.99, 4.74, and 6.73 mg/grain, while those of waxy broomcorn millet were 0.34, 5.94, and 6.28 mg/grain, respectively. The transcriptomic analysis revealed that 106 differentially expressed genes were identified, which were mainly enriched in the "amino sugar and nucleotide sugar metabolism", "pyruvate metabolism", "galactose metabolism", and "starch and sucrose metabolism" pathways. The WGCNA suggested that a total of 31 hub genes were correlated with starch biosynthesis. These findings provide a new approach to studying the starch synthesis in broomcorn millet.

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

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