Genome assemblies of 11 bamboo species highlight diversification induced by dynamic subgenome dominance.

Polyploidy (genome duplication) is a pivotal force in evolution. However, the interactions between parental genomes in a polyploid nucleus, frequently involving subgenome dominance, are poorly understood. Here we showcase analyses of a bamboo system (Poaceae: Bambusoideae) comprising a series of lineages from diploid (herbaceous) to tetraploid and hexaploid (woody), with 11 chromosome-level de novo genome assemblies and 476 transcriptome samples. We find that woody bamboo subgenomes exhibit stunning karyotype stability, with parallel subgenome dominance in the two tetraploid clades and a gradual shift of dominance in the hexaploid clade. Allopolyploidization and subgenome dominance have shaped the evolution of tree-like lignified culms, rapid growth and synchronous flowering characteristic of woody bamboos as large grasses. Our work provides insights into genome dominance in a remarkable polyploid system, including its dependence on genomic context and its ability to switch which subgenomes are dominant over evolutionary time.

The YABBY gene SHATTERING1 controls activation rather than patterning of the abscission zone in Setaria viridis.

Abscission is predetermined in specialized cell layers called the abscission zone (AZ) and activated by developmental or environmental signals. In the grass family, most identified AZ genes regulate AZ anatomy, which differs among lineages. A YABBY transcription factor, SHATTERING1 (SH1), is a domestication gene regulating abscission in multiple cereals, including rice and Setaria. In rice, SH1 inhibits lignification specifically in the AZ. However, the AZ of Setaria is nonlignified throughout, raising the question of how SH1 functions in species without lignification. Crispr-Cas9 knockout mutants of SH1 were generated in Setaria viridis and characterized with histology, cell wall and auxin immunofluorescence, transmission electron microscopy, hormonal treatment and RNA-Seq analysis. The sh1 mutant lacks shattering, as expected. No differences in cell anatomy or cell wall components including lignin were observed between sh1 and the wild-type (WT) until abscission occurs. Chloroplasts degenerated in the AZ of WT before abscission, but degeneration was suppressed by auxin treatment. Auxin distribution and expression of auxin-related genes differed between WT and sh1, with the signal of an antibody to auxin detected in the sh1 chloroplast. SH1 in Setaria is required for activation of abscission through auxin signaling, which is not reported in other grass species.

A happy accident: a novel turfgrass reference genome.

Poa pratensis, commonly known as Kentucky bluegrass, is a popular cool-season grass species used as turf in lawns and recreation areas globally. Despite its substantial economic value, a reference genome had not previously been assembled due to the genome's relatively large size and biological complexity that includes apomixis, polyploidy, and interspecific hybridization. We report here a fortuitous de novo assembly and annotation of a P. pratensis genome. Instead of sequencing the genome of a C4 grass, we accidentally sampled and sequenced tissue from a weedy P. pratensis whose stolon was intertwined with that of the C4 grass. The draft assembly consists of 6.09 Gbp with an N50 scaffold length of 65.1 Mbp, and a total of 118 scaffolds, generated using PacBio long reads and Bionano optical map technology. We annotated 256K gene models and found 58% of the genome to be composed of transposable elements. To demonstrate the applicability of the reference genome, we evaluated population structure and estimated genetic diversity in P. pratensis collected from three North American prairies, two in Manitoba, Canada and one in Colorado, USA. Our results support previous studies that found high genetic diversity and population structure within the species. The reference genome and annotation will be an important resource for turfgrass breeding and study of bluegrasses.

Grain shattering by cell death and fracture in Eragrostis tef.

Abscission, known as shattering in crop species, is a highly regulated process by which plants shed parts. Although shattering has been studied extensively in cereals and a number of regulatory genes have been identified, much diversity in the process remains to be discovered. Teff (Eragrostis tef) is a crop native to Ethiopia that is potentially highly valuable worldwide for its nutritious grain and drought tolerance. Previous work has suggested that grain shattering in Eragrostis might have little in common with other cereals. In this study, we characterize the anatomy, cellular structure, and gene regulatory control of the abscission zone (AZ) in E. tef. We show that the AZ of E. tef is a narrow stalk below the caryopsis, which is common in Eragrostis species. X-ray microscopy, scanning electron microscopy, transmission electron microscopy, and immunolocalization of cell wall components showed that the AZ cells are thin walled and break open along with programmed cell death (PCD) at seed maturity, rather than separating between cells as in other studied species. Knockout of YABBY2/SHATTERING1, documented to control abscission in several cereals, had no effect on abscission or AZ structure in E. tef. RNA sequencing analysis showed that genes related to PCD and cell wall modification are enriched in the AZ at the early seed maturity stage. These data show that E. tef drops its seeds using a unique mechanism. Our results provide the groundwork for understanding grain shattering in Eragrostis and further improvement of shattering in E. tef.

Diverse ecological functions and the convergent evolution of grass awns.

The awn of grasses is a long, conspicuous outgrowth of the floral bracts in a grass spikelet. It is known to impact agricultural yield, but we know little about its broader ecological function, nor the selective forces that lead to its evolution. Grass awns are phenotypically diverse across the extant ~12,000 species of Poaceae. Awns have been lost and gained repeatedly over evolutionary time, between and within lineages, suggesting that they could be under selection and might provide adaptive benefit in some environments. Despite the phylogenetic context, we know of no studies that have tested whether the origin of awns correlates with putative selective forces on their form and function. Presence or absence of awns is not plastic; rather, heritability is high. The awns of grasses often are suggested as adaptations for dispersal, and most experimental work has been aimed at testing this hypothesis. Proposed dispersal functions include soil burial, epizoochory, and aerial orientation. Awns may also protect the seed from drought, herbivores, or fire by helping it become buried in soil. We do not fully understand the fitness or nutrient costs of awn production, but in some species awns function in photosynthesis, providing carbon to the seed. Here we show that awns likely provide an adaptive advantage, but argue that studies on awn function have lacked critical phylogenetic information to demonstrate adaptive convergent evolution, are taxonomically biased, and often lack clear alternative hypotheses.

Molecular, cellular, and developmental foundations of grass diversity.

Humans have cultivated grasses for food, feed, beverages, and construction materials for millennia. Grasses also dominate the landscape in vast parts of the world, where they have adapted morphologically and physiologically, diversifying to form ~12,000 species. Sequences of hundreds of grass genomes show that they are essentially collinear; nonetheless, not all species have the same complement of genes. Here, we focus on the molecular, cellular, and developmental bases of grain yield and dispersal-traits that are essential for domestication. Distinct genes, networks, and pathways were selected in different crop species, reflecting underlying genomic diversity. With increasing genomic resources becoming available in nondomesticated species, we anticipate advances in coming years that illuminate the ecological and economic success of the grasses.

Pleiotropic and nonredundant effects of an auxin importer in Setaria and maize.

Directional transport of auxin is critical for inflorescence and floral development in flowering plants, but the role of auxin influx carriers (AUX1 proteins) has been largely overlooked. Taking advantage of available AUX1 mutants in green millet (Setaria viridis) and maize (Zea mays), we uncover previously unreported aspects of plant development that are affected by auxin influx, including higher order branches in the inflorescence, stigma branch number, glume (floral bract) development, and plant fertility. However, disruption of auxin flux does not affect all parts of the plant, with little obvious effect on inflorescence meristem size, time to flowering, and anther morphology. In double mutant studies in maize, disruptions of ZmAUX1 also affect vegetative development. A green fluorescent protein (GFP)-tagged construct of the Setaria AUX1 protein Sparse Panicle1 (SPP1) under its native promoter showed that SPP1 localizes to the plasma membrane of outer tissue layers in both roots and inflorescences, and accumulates specifically in inflorescence branch meristems, consistent with the mutant phenotype and expected auxin maxima. RNA-seq analysis indicated that most gene expression modules are conserved between mutant and wild-type plants, with only a few hundred genes differentially expressed in spp1 inflorescences. Using clustered regularly interspaced short palindromic repeats (CRISPR)-Cas9 technology, we disrupted SPP1 and the other four AUX1 homologs in S. viridis. SPP1 has a larger effect on inflorescence development than the others, although all contribute to plant height, tiller formation, and leaf and root development. The AUX1 importers are thus not fully redundant in S. viridis. Our detailed phenotypic characterization plus a stable GFP-tagged line offer tools for future dissection of the function of auxin influx proteins.

Genetic control of branching patterns in grass inflorescences.

Inflorescence branching in the grasses controls the number of florets and hence the number of seeds. Recent data on the underlying genetics come primarily from rice and maize, although new data are accumulating in other systems as well. This review focuses on a window in developmental time from the production of primary branches by the inflorescence meristem through to the production of glumes, which indicate the transition to producing a spikelet. Several major developmental regulatory modules appear to be conserved among most or all grasses. Placement and development of primary branches are controlled by conserved auxin regulatory genes. Subtending bracts are repressed by a network including TASSELSHEATH4, and axillary branch meristems are regulated largely by signaling centers that are adjacent to but not within the meristems themselves. Gradients of SQUAMOSA-PROMOTER BINDING-like and APETALA2-like proteins and their microRNA regulators extend along the inflorescence axis and the branches, governing the transition from production of branches to production of spikelets. The relative speed of this transition determines the extent of secondary and higher order branching. This inflorescence regulatory network is modified within individual species, particularly as regards formation of secondary branches. Differences between species are caused both by modifications of gene expression and regulators and by presence or absence of critical genes. The unified networks described here may provide tools for investigating orphan crops and grasses other than the well-studied maize and rice.

The Streptochaeta Genome and the Evolution of the Grasses.

In this work, we sequenced and annotated the genome of Streptochaeta angustifolia, one of two genera in the grass subfamily Anomochlooideae, a lineage sister to all other grasses. The final assembly size is over 99% of the estimated genome size. We find good collinearity with the rice genome and have captured most of the gene space. Streptochaeta is similar to other grasses in the structure of its fruit (a caryopsis or grain) but has peculiar flowers and inflorescences that are distinct from those in the outgroups and in other grasses. To provide tools for investigations of floral structure, we analyzed two large families of transcription factors, AP2-like and R2R3 MYBs, that are known to control floral and spikelet development in rice and maize among other grasses. Many of these are also regulated by small RNAs. Structure of the gene trees showed that the well documented whole genome duplication at the origin of the grasses (ρ) occurred before the divergence of the Anomochlooideae lineage from the lineage leading to the rest of the grasses (the spikelet clade) and thus that the common ancestor of all grasses probably had two copies of the developmental genes. However, Streptochaeta (and by inference other members of Anomochlooideae) has lost one copy of many genes. The peculiar floral morphology of Streptochaeta may thus have derived from an ancestral plant that was morphologically similar to the spikelet-bearing grasses. We further identify 114 loci producing microRNAs and 89 loci generating phased, secondary siRNAs, classes of small RNAs known to be influential in transcriptional and post-transcriptional regulation of several plant functions.

Intraspecific variation in elemental accumulation and its association with salt tolerance in Paspalum vaginatum.

Most plant species, including most crops, perform poorly in salt-affected soils because high sodium levels are cytotoxic and can disrupt the uptake of water and important nutrients. Halophytes are species that have evolved adaptations to overcome these challenges and may be a useful source of knowledge for salt tolerance mechanisms and genes that may be transferable to crop species. The salt content of saline habitats can vary dramatically by location, providing ample opportunity for different populations of halophytic species to adapt to their local salt concentrations; however, the extent of this variation, and the physiology and polymorphisms that drive it, remain poorly understood. Differential accumulation of inorganic elements between genotypes or populations may play an important role in local salinity adaptation. To test this, we investigated the relationships between population structure, tissue ion concentrations, and salt tolerance in 17 "fine-textured" genotypes of the halophytic turfgrass seashore paspalum (Paspalum vaginatum Swartz). A high-throughput ionomics pipeline was used to quantify the shoot concentration of 18 inorganic elements across three salinity treatments. We found a significant relationship between population structure and ion accumulation, with strong correlations between principal components derived from genetic and ionomic data. Additionally, genotypes with higher salt tolerance accumulated more K and Fe and less Ca than less tolerant genotypes. Together these results indicate that differences in ion accumulation between P. vaginatum populations may reflect locally adapted salt stress responses.

Conserved noncoding sequences provide insights into regulatory sequence and loss of gene expression in maize.

Thousands of species will be sequenced in the next few years; however, understanding how their genomes work, without an unlimited budget, requires both molecular and novel evolutionary approaches. We developed a sensitive sequence alignment pipeline to identify conserved noncoding sequences (CNSs) in the Andropogoneae tribe (multiple crop species descended from a common ancestor ∼18 million years ago). The Andropogoneae share similar physiology while being tremendously genomically diverse, harboring a broad range of ploidy levels, structural variation, and transposons. These contribute to the potential of Andropogoneae as a powerful system for studying CNSs and are factors we leverage to understand the function of maize CNSs. We found that 86% of CNSs were comprised of annotated features, including introns, UTRs, putative cis-regulatory elements, chromatin loop anchors, noncoding RNA (ncRNA) genes, and several transposable element superfamilies. CNSs were enriched in active regions of DNA replication in the early S phase of the mitotic cell cycle and showed different DNA methylation ratios compared to the genome-wide background. More than half of putative cis-regulatory sequences (identified via other methods) overlapped with CNSs detected in this study. Variants in CNSs were associated with gene expression levels, and CNS absence contributed to loss of gene expression. Furthermore, the evolution of CNSs was associated with the functional diversification of duplicated genes in the context of maize subgenomes. Our results provide a quantitative understanding of the molecular processes governing the evolution of CNSs in maize.

The CLV3 Homolog in Setaria viridis Selectively Controls Inflorescence Meristem Size.

The CLAVATA pathway controls meristem size during inflorescence development in both eudicots and grasses, and is initiated by peptide ligands encoded by CLV3/ESR-related (CLE) genes. While CLV3 controls all shoot meristems in Arabidopsis, evidence from cereal grasses indicates that different meristem types are regulated by different CLE peptides. The rice peptide FON2 primarily controls the size of the floral meristem, whereas the orthologous peptides CLE7 and CLE14 in maize have their most dramatic effects on inflorescence and branch meristems, hinting at diversification among CLE responses in the grasses. Setaria viridis is more closely related to maize than to rice, so can be used to test whether the maize CLE network can be generalized to all members of subfamily Panicoideae. We used CRISPR-Cas9 in S. viridis to knock out the SvFON2 gene, the closest homolog to CLV3 and FON2. Svfon2 mutants developed larger inflorescence meristems, as in maize, but had normal floral meristems, unlike Osfon2, suggesting a panicoid-specific CLE network. Vegetative traits such as plant height, tiller number and leaf number were not significantly different between mutant and wild type plants, but time to heading was shorter in the mutants. In situ hybridization showed strong expression of Svfon2 in the inflorescence and branch meristems, consistent with the mutant phenotype. Using bioinformatic analysis, we predicted the co-expression network of SvFON2 and its signaling components, which included genes known to control inflorescence architecture in maize as well as genes of unknown function. The similarity between SvFON2 function in Setaria and maize suggests that its developmental specialization in inflorescence meristem control may be shared among panicoid grasses.

Hybridization, polyploidy and clonality influence geographic patterns of diversity and salt tolerance in the model halophyte seashore paspalum (Paspalum vaginatum).

In plant species, variation in levels of clonality, ploidy and interspecific hybridization can interact to influence geographic patterns of genetic diversity. These factors commonly vary in plants that specialize on saline habitats (halophytes) and may play a role in how they adapt to salinity variation across their range. One such halophyte is the turfgrass and emerging genomic model system seashore paspalum (Paspalum vaginatum Swartz). To investigate how clonal propagation, ploidy variation, and interspecific hybridization vary across ecotypes and local salinity levels in wild P. vaginatum, we employed genotyping-by-sequencing, cpDNA sequencing and flow cytometry in 218 accessions representing > 170 wild collections from throughout the coastal southern United States plus USDA germplasm. We found that the two morphologically distinct ecotypes of P. vaginatum differ in their adaptive strategies. The fine-textured ecotype is diploid and appears to reproduce in the wild both sexually and by clonal propagation; in contrast, the coarse-textured ecotype consists largely of clonally-propagating triploid and diploid genotypes. The coarse-textured ecotype appears to be derived from hybridization between fine-textured P. vaginatum and an unidentified Paspalum species. These clonally propagating hybrid genotypes are more broadly distributed than clonal fine-textured genotypes and may represent a transition to a more generalist adaptive strategy. Additionally, the triploid genotypes vary in whether they carry one or two copies of the P. vaginatum subgenome, indicating multiple evolutionary origins. This variation in subgenome composition shows associations with local ocean salinity levels across the sampled populations and may play a role in local adaptation.

A genome resource for green millet Setaria viridis enables discovery of agronomically valuable loci.

Wild and weedy relatives of domesticated crops harbor genetic variants that can advance agricultural biotechnology. Here we provide a genome resource for the wild plant green millet (Setaria viridis), a model species for studies of C4 grasses, and use the resource to probe domestication genes in the close crop relative foxtail millet (Setaria italica). We produced a platinum-quality genome assembly of S. viridis and de novo assemblies for 598 wild accessions and exploited these assemblies to identify loci underlying three traits: response to climate, a 'loss of shattering' trait that permits mechanical harvest and leaf angle, a predictor of yield in many grass crops. With CRISPR-Cas9 genome editing, we validated Less Shattering1 (SvLes1) as a gene whose product controls seed shattering. In S. italica, this gene was rendered nonfunctional by a retrotransposon insertion in the domesticated loss-of-shattering allele SiLes1-TE (transposable element). This resource will enhance the utility of S. viridis for dissection of complex traits and biotechnological improvement of panicoid crops.

Evolutionary Dynamics of Transposable Elements Following a Shared Polyploidization Event in the Tribe Andropogoneae.

Both polyploidization and transposable element (TE) activity are known to be major drivers of plant genome evolution. Here, we utilize the Zea-Tripsacum clade to investigate TE activity and accumulation after a shared polyploidization event. Comparisons of TE evolutionary dynamics in various Zea and Tripsacum species, in addition to two closely related diploid species, Urelytrum digitatum and Sorghum bicolor, revealed variation in repeat content among all taxa included in the study. The repeat composition of Urelytrum is more similar to that of Zea and Tripsacum compared to Sorghum, despite the similarity in genome size with the latter. Although LTR-retrotransposons were abundant in all species, we observed an expansion of the copia superfamily, specifically in Z. mays and T. dactyloides, species that have adapted to more temperate environments. Additional analyses of the genomic distribution of these retroelements provided evidence of biased insertions near genes involved in various biological processes including plant development, defense, and macromolecule biosynthesis. Specifically, copia insertions in Zea and T. dactyloides were significantly enriched near genes involved in abiotic stress response, suggesting independent evolution post Zea-Tripsacum divergence. The lack of copia insertions near the orthologous genes in S. bicolor suggests that duplicate gene copies generated during polyploidization may offer novel neutral sites for TEs to insert, thereby providing an avenue for subfunctionalization via TE insertional mutagenesis.

Sterile Spikelets Contribute to Yield in Sorghum and Related Grasses.

Sorghum (Sorghum bicolor) and its relatives in the grass tribe Andropogoneae bear their flowers in pairs of spikelets in which one spikelet (seed-bearing or sessile spikelet [SS]) of the pair produces a seed and the other is sterile or male (staminate). This division of function does not occur in other major cereals such as wheat (Triticum aestivum) or rice (Oryza sativa). Additionally, one bract of the SS spikelet often produces a long extension, the awn, that is in the same position as, but independently derived from, that of wheat and rice. The function of the sterile spikelet is unknown and that of the awn has not been tested in Andropogoneae. We used radioactive and stable isotopes of carbon, RNA sequencing of metabolically important enzymes, and immunolocalization of ribulose-1,5-bisphosphate carboxylase/oxygenase (Rubisco) to show that the sterile spikelet assimilates carbon, which is translocated to the largely heterotrophic SS. The awn shows no evidence of photosynthesis. These results apply to distantly related species of Andropogoneae. Removal of sterile spikelets in sorghum significantly decreases seed weight (yield) by ∼9%. Thus, the sterile spikelet, but not the awn, affects yield in the cultivated species and fitness in the wild species.

The Genomes of the Allohexaploid Echinochloa crus-galli and Its Progenitors Provide Insights into Polyploidization-Driven Adaptation.

The hexaploid species Echinochloa crus-galli is one of the most detrimental weeds in crop fields, especially in rice paddies. Its evolutionary history is similar to that of bread wheat, arising through polyploidization after hybridization between a tetraploid and a diploid species. In this study, we generated and analyzed high-quality genome sequences of diploid (E. haploclada), tetraploid (E. oryzicola), and hexaploid (E. crus-galli) Echinochloa species. Gene family analysis showed a significant loss of disease-resistance genes such as those encoding NB-ARC domain-containing proteins during Echinochloa polyploidization, contrary to their significant expansionduring wheat polyploidization, suggesting that natural selection might favor reduced investment in resistance in this weed to maximize its growth and reproduction. In contrast to the asymmetric patterns of genome evolution observed in wheat and other crops, no significant differences in selection pressure were detected between the subgenomes in E. oryzicola and E. crus-galli. In addition, distinctive differences in subgenome transcriptome dynamics during hexaploidization were observed between E. crus-galli and bread wheat. Collectively, our study documents genomic mechanisms underlying the adaptation of a major agricultural weed during polyploidization. The genomic and transcriptomic resources of three Echinochloa species and new insights into the polyploidization-driven adaptive evolution would be useful for future breeding cereal crops.

Comprehensive 3D phenotyping reveals continuous morphological variation across genetically diverse sorghum inflorescences.

●Inflorescence architecture in plants is often complex and challenging to quantify, particularly for inflorescences of cereal grasses. Methods for capturing inflorescence architecture and for analyzing the resulting data are limited to a few easily captured parameters that may miss the rich underlying diversity. ●Here, we apply X-ray computed tomography combined with detailed morphometrics, offering new imaging and computational tools to analyze three-dimensional inflorescence architecture. To show the power of this approach, we focus on the panicles of Sorghum bicolor, which vary extensively in numbers, lengths, and angles of primary branches, as well as the three-dimensional shape, size, and distribution of the seed. ●We imaged and comprehensively evaluated the panicle morphology of 55 sorghum accessions that represent the five botanical races in the most common classification system of the species, defined by genetic data. We used our data to determine the reliability of the morphological characters for assigning specimens to race and found that seed features were particularly informative. ●However, the extensive overlap between botanical races in multivariate trait space indicates that the phenotypic range of each group extends well beyond its overall genetic background, indicating unexpectedly weak correlation between morphology, genetic identity, and domestication history.