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.

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.

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.

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.

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.

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.

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.

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.

Overexpression of PvBiP2 improved biomass yield and cadmium tolerance in switchgrass (Panicum virgatum L.).

Switchgrass (Panicum virgatum L.), the prime bioenergy feedstock crop, is one ideal candidate for phytoremediation of cadmium (Cd). The absorption of Cd imposes severe endoplasmic reticulum (ER)-stress in plants. ER chaperone binding proteins (BiPs) are important modulators in ER-stress responses. The objective of this study was to characterize one Cd-responsive BiP gene, PvBiP2, in switchgrass for its roles in Cd tolerance and plant growth. PvBiP2 was up-regulated by Cd and the ER-stress inducer, dithiothreitol (DTT) and could be trans-activated by one Cd-responsive heat shock transcription factor PvHsfA4. Overexpression of PvBiP2 in switchgrass significantly increased its plant growth with higher height, stem diameter, leaf width, internode length, and tiller numbers than those of the wildtype (WT) plants under non-stress conditions. After 30 days of Cd treatment, the PvBiP2 over-expression transgenic lines showed 40-45% higher dry biomass accumulation with net photosynthesis rate (Pn), but lower electrolyte leakage (EL), malondialdehyde (MDA), and glutathione (GSH) levels than WT. Moreover, over-expressing PvBiP2 led to ∼90-140% Cd accumulation in plants but 46-57% lower Cd translocation rates to shoots. Together, the PvHSFA4-PvBiP2 module acted as positive regulators in plant Cd tolerance, and over-expressing PvBiP2 promoted plant vegetative growth as well as Cd tolerance making it an ideal molecular target for genetic improvement in switchgrass in the future.

Legacies of precipitation influence primary production in Panicum virgatum.

Precipitation is a key driver of primary production worldwide, but primary production does not always track year-to-year variation in precipitation linearly. Instead, plant responses to changes in precipitation may exhibit time lags, or legacies of past precipitation. Legacies can be driven by multiple mechanisms, including persistent changes in plant physiological and morphological traits and changes to the physical environment, such as plant access to soil water. We used three precipitation manipulation experiments in central Texas, USA to evaluate the magnitude, duration, and potential mechanisms driving precipitation legacies on aboveground primary production of the perennial C4 grass, Panicum virgatum. Specifically, we performed a rainout shelter study, where eight genotypes grew under different precipitation regimes; a transplant study, where plants that had previously grown in a rainout shelter under different precipitation regimes were moved to a common environment; and a mesocosm study, where the effect of swapping precipitation regime was examined with a single genotype. Across these experiments, plants previously grown under wet conditions generally performed better than expected when exposed to drought. Panicum virgatum exhibited stronger productivity legacies of past wet years on current-year responses to drought than of past dry years on current-year responses to wet conditions. Additionally, previous year tiller counts, a proxy for meristem availability, were important in determining legacy effects on aboveground production. As climate changes and precipitation extremes-both dry and wet-become more common, these results suggest that populations of P. virgatum may become less resilient.

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.

Cold acclimation threshold induction temperatures of switchgrass ecotypes grown under a long and short photoperiod.

Plants can cold acclimate to enhance their freezing tolerance by sensing declining temperature and photoperiod cues. However, the factors influencing genotypic variation in the induction of cold acclimation are poorly understood among perennial grasses. We hypothesized that the more northern upland switchgrass (Panicum virgatum L.) ecotype develops a higher degree of freezing tolerance by initiating cold acclimation at higher temperatures as compared with the coastal and southern lowland ecotypes. First, we determined the optimal method for assessing freezing tolerance and the length of exposure to 8/4°C required to induce the maximum level of freezing tolerance in the most northern upland and most southern lowland genotypes. We characterized the maximum freezing tolerance of eight uplands, three coastal and five lowland genotypes grown for 21 days at 8/4°C and a 10 or 16 h photoperiod. Next, we identified the temperature required to induce cold acclimation by exposing the 16 genotypes for 7 days at 20-6°C constant temperatures under a 10 or 16 h photoperiod. Cold acclimation initiated at temperatures 5 and 7°C higher in upland than in coastal and lowland genotypes. Among upland genotypes the shorter photoperiod induced cold acclimation at a 1°C higher temperature. Genotypes originating from a more northern latitude initiate cold acclimation at higher temperatures and develop higher maximum freezing tolerances. An earlier response to declining temperatures may provide the upland ecotype with additional time to prepare for winter and provide an advantage when plants are subjected to the rapid changes in fall temperature associated with injurious frosts.

An Agrobacterium strain auxotrophic for methionine is useful for switchgrass transformation.

Auxotrophic strains of Agrobacterium tumefaciens can contribute to the development of more efficient transformation systems, especially for crops historically considered recalcitrant. Homologous recombination was used to derive methionine auxotrophs of two common A. tumefaciens strains, LBA4404 and EHA105. The EHA105 strains were more efficient for switchgrass transformation, while both the EHA105 and LBA4404 strains worked equally well for the rice control. Event quality, as measured by transgene copy number, was not affected by auxotrophy, but was higher for the LBA4404 strains than the EHA105 strains. Ultimately, the use of auxotrophs reduced bacterial overgrowth during co-cultivation and decreased the need for antibiotics.