Genome Alignment Spanning Major Poaceae Lineages Reveals Heterogeneous Evolutionary Rates and Alters Inferred Dates for Key Evolutionary Events

Deciphering recursive polyploidization in Lamiales and reconstructing their chromosome evolutionary trajectories

Lamiales is an order of core eudicots with abundant diversity, and many Lamiales plants have important medicinal and ornamental values. Here, we comparatively reanalyzed 11 Lamiales species with well-assembled genome sequences and found evidence that Lamiales plants, in addition to a hexaploidization or whole-genome triplication (WGT) shared by core eudicots, experienced further polyploidization events, establishing new groups in the order. Notably, we identified a whole-genome duplication (WGD) occurred just before the split of Scrophulariaceae from the other Lamiales families, such as Acanthaceae, Bignoniaceae, and Lamiaceae, suggesting its likely being the causal reason for the establishment and fast divergence of these families. We also found that a WGT occurred ∼68 to 78 million years ago (Mya), near the split of Oleaceae from the other Lamiales families, implying that it may have caused their fast divergence and the establishment of the Oleaceae family. Then, by exploring and distinguishing intra- and intergenomic chromosomal homology due to recursive polyploidization and speciation, respectively, we inferred that the Lamiales ancestral cell karyotype had 11 proto-chromosomes. We reconstructed the evolutionary trajectories from these proto-chromosomes to form the extant chromosomes in each Lamiales plant under study. We must note that most of the inferred 11 proto-chromosomes, duplicated during a WGD thereafter, have been well preserved in jacaranda (Jacaranda mimosifolia) genome, showing the credibility of the present inference implementing a telomere-centric chromosome repatterning model. These efforts are important to understand genome repatterning after recursive polyploidization, especially shedding light on the origin of new plant groups and angiosperm cell karyotype evolution.

Wild grass-derived alleles represent a genetic architecture for the resilience of modern common wheat to stresses

This review explores the integration of wild grass-derived alleles into modern bread wheat breeding to tackle the challenges of climate change and increasing food demand. With a focus on synthetic hexaploid wheat, this review highlights the potential of genetic variability in wheat wild relatives, particularly Aegilops tauschii, for improving resilience to multifactorial stresses like drought, heat, and salinity. The evolutionary journey of wheat (Triticum spp.) from diploid to hexaploid species is examined, revealing significant genetic contributions from wild grasses. We also emphasize the importance of understanding incomplete lineage sorting in the genomic evolution of wheat. Grasping this information is crucial as it can guide breeders in selecting the appropriate alleles from the gene pool of wild relatives to incorporate into modern wheat varieties. This approach improves the precision of phylogenetic relationships and increases the overall effectiveness of breeding strategies. This review also addresses the challenges in utilizing the wheat wild genetic resources, such as the linkage drag and cross-compatibility issues. Finally, we culminate the review with future perspectives, advocating for a combined approach of high-throughput phenotyping tools and advanced genomic techniques to comprehensively understand the genetic and regulatory architectures of wheat under stress conditions, paving the way for more precise and efficient breeding strategies.

Whole-genome sequencing of allotetraploid bermudagrass reveals the origin of Cynodon and candidate genes for salt tolerance

Bermudagrass (Cynodon dactylon) is a globally distributed, extensively used warm-season turf and forage grass with high tolerance to salinity and drought stress in alkaline environments. However, the origin of the species and genetic mechanisms for salinity tolerance in the species are basically unknown. Accordingly, we set out to study evolution divergence events in the Cynodon genome and to identify genes for salinity tolerance. We developed a 604.0 Mb chromosome-level polyploid genome sequence for bermudagrass 'A12359' (n = 18). The C. dactylon genome comprises 2 complete sets of homoeologous chromosomes, each with approximately 30 000 genes, and most genes are conserved as syntenic pairs. Phylogenetic study showed that the initial Cynodon species diverged from Oropetium thomaeum approximately 19.7-25.4 million years ago (Mya), the A and B subgenomes of C. dactylon diverged approximately 6.3-9.1 Mya, and the bermudagrass polyploidization event occurred 1.5 Mya on the African continent. Moreover, we identified 82 candidate genes associated with seven agronomic traits using a genome-wide association study, and three single-nucleotide polymorphisms were strongly associated with three salt resistance genes: RAP2-2, CNG channels, and F14D7.1. These genes may be associated with enhanced bermudagrass salt tolerance. These bermudagrass genomic resources, when integrated, may provide fundamental insights into evolution of diploid and tetraploid genomes and enhance the efficacy of comparative genomics in studying salt tolerance in Cynodon.

Identification, molecular evolution, codon bias, and expansion analysis of NLP transcription factor family in foxtail millet (Setaria italica L.) and closely related crops

The NODULE-INCEPTION-like protein (NLP) family is a plant-specific transcription factor (TF) family involved in nitrate transport and assimilation in plants, which are essential for improving plant nitrogen use efficiency. Currently, the molecular nature and evolutionary trajectory of NLP genes in the C4 model crop foxtail millet are unknown. Therefore, we performed a comprehensive analysis of NLP and molecular evolution in foxtail millet by scanning the genomes of foxtail millet and representative species of the plant kingdom. We identified seven NLP genes in the foxtail millet genome, all of which are individually and separately distributed on different chromosomes. They were not structurally identical to each other and were mainly expressed on root tissues. We unearthed two key genes (Si5G004100.1 and Si6G248300.1) with a variety of excellent characteristics. Regarding its molecular evolution, we found that NLP genes in Gramineae mainly underwent dispersed duplication, but maize NLP genes were mainly generated via WGD events. Other factors such as base mutations and natural selection have combined to promote the evolution of NLP genes. Intriguingly, the family in plants showed a gradual expansion during evolution with more duplications than losses, contrary to most gene families. In conclusion, this study advances the use of NLP genetic resources and the understanding of molecular evolution in cereals.

The super-pangenome of Populus unveils genomic facets for its adaptation and diversification in widespread forest trees

Understanding the underlying mechanisms and links between genome evolution and adaptive innovations stands as a key goal in evolutionary studies. Poplars, among the world's most widely distributed and cultivated trees, exhibit extensive phenotypic diversity and environmental adaptability. In this study, we present a genus-level super-pangenome comprising 19 Populus genomes, revealing the likely pivotal role of private genes in facilitating local environmental and climate adaptation. Through the integration of pangenomes with transcriptomes, methylomes, and chromatin accessibility mapping, we unveil that the evolutionary trajectories of pangenes and duplicated genes are closely linked to local genomic landscapes of regulatory and epigenetic architectures, notably CG methylation in gene-body regions. Further comparative genomic analyses have enabled the identification of 142 202 structural variants across species that intersect with a significant number of genes and contribute substantially to both phenotypic and adaptive divergence. We have experimentally validated a ∼180-bp presence/absence variant affecting the expression of the CUC2 gene, crucial for leaf serration formation. Finally, we developed a user-friendly web-based tool encompassing the multi-omics resources associated with the Populus super-pangenome (http://www.populus-superpangenome.com). Together, the present pioneering super-pangenome resource in forest trees not only aids in the advancement of breeding efforts of this globally important tree genus but also offers valuable insights into potential avenues for comprehending tree biology.

Genome organization and botanical diversity

The rich diversity of angiosperms, both the planet's dominant flora and the cornerstone of agriculture, is integrally intertwined with a distinctive evolutionary history. Here, we explore the interplay between angiosperm genome organization and botanical diversity, empowered by genomic approaches ranging from genetic linkage mapping to analysis of gene regulation. Commonality in the genetic hardware of plants has enabled robust comparative genomics that has provided a broad picture of angiosperm evolution and implicated both general processes and specific elements in contributing to botanical diversity. We argue that the hardware of plant genomes-both in content and in dynamics-has been shaped by selection for rather substantial differences in gene regulation between plants and animals such as maize and human, organisms of comparable genome size and gene number. Their distinctive genome content and dynamics may reflect in part the indeterminate development of plants that puts strikingly different demands on gene regulation than in animals. Repeated polyploidization of plant genomes and multiplication of individual genes together with extensive rearrangement and differential retention provide rich raw material for selection of morphological and/or physiological variations conferring fitness in specific niches, whether natural or artificial. These findings exemplify the burgeoning information available to employ in increasing knowledge of plant biology and in modifying selected plants to better meet human needs.

Functional differentiation and genetic diversity of rice cation exchanger (CAX) genes and their potential use in rice improvement

Cation exchanger (CAX) genes play an important role in plant growth/development and response to biotic and abiotic stresses. Here, we tried to obtain important information on the functionalities and phenotypic effects of CAX gene family by systematic analyses of their expression patterns, genetic diversity (gene CDS haplotypes, structural variations, gene presence/absence variations) in 3010 rice genomes and nine parents of 496 Huanghuazhan introgression lines, the frequency shifts of the predominant gcHaps at these loci to artificial selection during modern breeding, and their association with tolerances to several abiotic stresses. Significant amounts of variation also exist in the cis-regulatory elements (CREs) of the OsCAX gene promoters in 50 high-quality rice genomes. The functional differentiation of OsCAX gene family were reflected primarily by their tissue and development specific expression patterns and in varied responses to different treatments, by unique sets of CREs in their promoters and their associations with specific agronomic traits/abiotic stress tolerances. Our results indicated that OsCAX1a and OsCAX2 as general signal transporters were in many processes of rice growth/development and responses to diverse environments, but they might be of less value in rice improvement. OsCAX1b, OsCAX1c, OsCAX3 and OsCAX4 was expected to be of potential value in rice improvement because of their associations with specific traits, responsiveness to specific abiotic stresses or phytohormones, and relatively high gcHap and CRE diversity. Our strategy was demonstrated to be highly efficient to obtain important genetic information on genes/alleles of specific gene family and can be used to systematically characterize the other rice gene families.

The global distribution of angiosperm genome size is shaped by climate

Angiosperms, which inhabit diverse environments across all continents, exhibit significant variation in genome sizes, making them an excellent model system for examining hypotheses about the global distribution of genome size. These include the previously proposed large genome constraint, mutational hazard, polyploidy-mediated, and climate-mediated hypotheses. We compiled the largest genome size dataset to date, encompassing 16 017 (> 5% of known) angiosperm species, and analyzed genome size distribution using a comprehensive geographic distribution dataset for all angiosperms. We observed that angiosperms with large range sizes generally had small genomes, supporting the large genome constraint hypothesis. Climate was shown to exert a strong influence on genome size distribution along the global latitudinal gradient, while the frequency of polyploidy and the type of growth form had negligible effects. In contrast to the unimodal patterns along the global latitudinal gradient shown by plant size traits and polyploid proportions, the increase in angiosperm genome size from the equator to 40-50°N/S is probably mediated by different (mostly climatic) mechanisms than the decrease in genome sizes observed from 40 to 50°N northward. Our analysis suggests that the global distribution of genome sizes in angiosperms is mainly shaped by climatically mediated purifying selection, genetic drift, relaxed selection, and environmental filtering.

Sorghum SbGhd7 is a major regulator of floral transition and directly represses genes crucial for flowering activation

Molecular genetic understanding of flowering time regulation is crucial for sorghum development. GRAIN NUMBER, PLANT HEIGHT AND HEADING DATE 7 (SbGhd7) is one of the six classical loci conferring photoperiod sensitivity of sorghum flowering. However, its functions remain poorly studied. The molecular functions of SbGhd7 were characterized. The gene regulatory network controlled by SbGhd7 was constructed and validated. The biological roles of SbGhd7 and its major targets were studied. SbGhd7 overexpression (OE) completely prevented sorghum flowering. Additionally, we show that SbGhd7 is a major negative regulator of flowering, binding to the promoter motif TGAATG(A/T)(A/T/C) and repressing transcription of the major florigen FLOWERING LOCUS T 10 (SbFT10) and floral activators EARLY HEADING DATE (SbEhd1), FLAVIN-BINDING, KELCH REPEAT, F-BOX1 (SbFKF1) and EARLY FLOWERING 3 (SbELF3). Reinforcing the direct effect of SbGhd7, SbEhd1 OE activated the promoters of three functional florigens (SbFT1, SbFT8 and SbFT10), dramatically accelerating flowering. Our studies demonstrate that SbGhd7 is a major repressor of sorghum flowering by directly and indirectly targeting genes for flowering activation. The mechanism appears ancient. Our study extends the current model of floral transition regulation in sorghum and provides a framework for a comprehensive understanding of sorghum photoperiod response.

Assembly and comparative analysis of the initial complete mitochondrial genome of Primulina hunanensis (Gesneriaceae): a cave-dwelling endangered plant

BACKGROUND

Primulina hunanensis, a troglobitic plant within the Primulina genus of Gesneriaceae family, exhibits robust resilience to arid conditions and holds great horticultural potential as an ornamental plant. The work of chloroplast genome (cpDNA) has been recently accomplished, however, the mitochondrial genome (mtDNA) that is crucial for plant evolution has not been reported.

RESULTS

In this study, we sequenced and assembled the P. hunanensis complete mtDNA, and elucidated its evolutionary and phylogenetic relationships. The assembled mtDNA spans 575,242 bp with 43.54% GC content, encompassing 60 genes, including 37 protein-coding genes (PCGs), 20 tRNA genes, and 3 rRNA genes. Notably, high number of repetitive sequences in the mtDNA and substantial sequence translocation from chloroplasts to mitochondria were observed. To determine the evolutionary and taxonomic positioning of P. hunanensis, a phylogenetic tree was constructed using mitochondrial PCGs from P. hunanensis and 32 other taxa. Furthermore, an exploration of PCGs relative synonymous codon usage, identification of RNA editing events, and an investigation of collinearity with closely related species were conducted.

CONCLUSIONS

This study reports the initial assembly and annotation of P. hunanensis mtDNA, contributing to the limited mtDNA repository for Gesneriaceae plants and advancing our understanding of their evolution for improved utilization and conservation.

Striking a growth-defense balance: Stress regulators that function in maize development

Maize (Zea mays) cultivation is strongly affected by both abiotic and biotic stress, leading to reduced growth and productivity. It has recently become clear that regulators of plant stress responses, including the phytohormones abscisic acid (ABA), ethylene (ET), and jasmonic acid (JA), together with reactive oxygen species (ROS), shape plant growth and development. Beyond their well established functions in stress responses, these molecules play crucial roles in balancing growth and defense, which must be finely tuned to achieve high yields in crops while maintaining some level of defense. In this review, we provide an in-depth analysis of recent research on the developmental functions of stress regulators, focusing specifically on maize. By unraveling the contributions of these regulators to maize development, we present new avenues for enhancing maize cultivation and growth while highlighting the potential risks associated with manipulating stress regulators to enhance grain yields in the face of environmental challenges.

Whole Genome Duplication Events Likely Contributed to the Aquatic Adaptive Evolution of Parkerioideae

As the only aquatic lineage of Pteridaceae, Parkerioideae is distinct from many xeric-adapted species of the family and consists of the freshwater Ceratopteris species and the only mangrove ferns from the genus Acrostichum. Previous studies have shown that whole genome duplication (WGD) has occurred in Parkerioideae at least once and may have played a role in their adaptive evolution; however, more in-depth research regarding this is still required. In this study, comparative and evolutionary transcriptomics analyses were carried out to identify WGDs and explore their roles in the environmental adaptation of Parkerioideae. Three putative WGD events were identified within Parkerioideae, two of which were specific to Ceratopteris and Acrostichum, respectively. The functional enrichment analysis indicated that the lineage-specific WGD events have played a role in the adaptation of Parkerioideae to the low oxygen concentrations of aquatic habitats, as well as different aquatic environments of Ceratopteris and Acrostichum, such as the adaptation of Ceratopteris to reduced light levels and the adaptation of Acrostichum to high salinity. Positive selection analysis further provided evidence that the putative WGD events may have facilitated the adaptation of Parkerioideae to changes in habitat. Moreover, the gene family analysis indicated that the plasma membrane H+-ATPase (AHA), vacuolar H+-ATPase (VHA), and suppressor of K+ transport growth defect 1 (SKD1) may have been involved in the high salinity adaptation of Acrostichum. Our study provides new insights into the evolution and adaptations of Parkerioideae in different aquatic environments.

Highly differentiated genomic properties underpin the different cell walls of Poaceae and eudicots

Plant cell walls of Poaceae and eudicots differ substantially, both in the content and composition of their components. However, the genomic and genetic basis underlying these differences is not fully resolved. In this research, we analyzed multiple genomic properties of 150 cell wall gene families across 169 angiosperm genomes. The properties analyzed include gene presence/absence, copy number, synteny, occurrence of tandem gene clusters, and phylogenetic gene diversity. Results revealed a profound genomic differentiation of cell wall genes between Poaceae and eudicots, often associated with the cell wall diversity between these plant groups. For example, overall patterns of gene copy number variation and synteny were clearly divergent between Poaceae and eudicot species. Moreover, differential Poaceae-eudicot copy number and genomic contexts were observed for all the genes within the BEL1-like HOMEODOMAIN 6 regulatory pathway, which respectively induces and represses secondary cell wall synthesis in Poaceae and eudicots. Similarly, divergent synteny, copy number, and phylogenetic gene diversification were observed for the major biosynthetic genes of xyloglucans, mannans, and xylans, potentially contributing to the differences in content and types of hemicellulosic polysaccharides differences in Poaceae and eudicot cell walls. Additionally, the Poaceae-specific tandem clusters and/or higher copy number of PHENYLALANINE AMMONIA-LYASE, CAFFEIC ACID O-METHYLTRANSFERASE, or PEROXIDASE genes may underly the higher content and larger variety of phenylpropanoid compounds observed in Poaceae cell walls. All these patterns are discussed in detail in this study, along with their evolutionary and biological relevance for cell wall (genomic) diversification between Poaceae and eudicots.

Genome-wide expansion and reorganization during grass evolution: from 30 Mb chromosomes in rice and Brachypodium to 550 Mb in Avena

BACKGROUND

The BOP (Bambusoideae, Oryzoideae, and Pooideae) clade of the Poaceae has a common ancestor, with similarities to the genomes of rice, Oryza sativa (2n = 24; genome size 389 Mb) and Brachypodium, Brachypodium distachyon (2n = 10; 271 Mb). We exploit chromosome-scale genome assemblies to show the nature of genomic expansion, structural variation, and chromosomal rearrangements from rice and Brachypodium, to diploids in the tribe Aveneae (e.g., Avena longiglumis, 2n = 2x = 14; 3,961 Mb assembled to 3,850 Mb in chromosomes).

RESULTS

Most of the Avena chromosome arms show relatively uniform expansion over the 10-fold to 15-fold genome-size increase. Apart from non-coding sequence diversification and accumulation around the centromeres, blocks of genes are not interspersed with blocks of repeats, even in subterminal regions. As in the tribe Triticeae, blocks of conserved synteny are seen between the analyzed species with chromosome fusion, fission, and nesting (insertion) events showing deep evolutionary conservation of chromosome structure during genomic expansion. Unexpectedly, the terminal gene-rich chromosomal segments (representing about 50 Mb) show translocations between chromosomes during speciation, with homogenization of genome-specific repetitive elements within the tribe Aveneae. Newly-formed intergenomic translocations of similar extent are found in the hexaploid A. sativa.

CONCLUSIONS

The study provides insight into evolutionary mechanisms and speciation in the BOP clade, which is valuable for measurement of biodiversity, development of a clade-wide pangenome, and exploitation of genomic diversity through breeding programs in Poaceae.

Chromosome-level genomes of three key Allium crops and their trait evolution

Allium crop breeding remains severely hindered due to the lack of high-quality reference genomes. Here we report high-quality chromosome-level genome assemblies for three key Allium crops (Welsh onion, garlic and onion), which are 11.17 Gb, 15.52 Gb and 15.78 Gb in size with the highest recorded contig N50 of 507.27 Mb, 109.82 Mb and 81.66 Mb, respectively. Beyond revealing the genome evolutionary process of Allium species, our pathogen infection experiments and comparative metabolomic and genomic analyses showed that genes encoding enzymes involved in the metabolic pathway of Allium-specific flavor compounds may have evolved from an ancient uncharacterized plant defense system widely existing in many plant lineages but extensively boosted in alliums. Using in situ hybridization and spatial RNA sequencing, we obtained an overview of cell-type categorization and gene expression changes associated with spongy mesophyll cell expansion during onion bulb formation, thus indicating the functional roles of bulb formation genes.

Diploid and tetraploid genomes of Acorus and the evolution of monocots

Monocots are a major taxon within flowering plants, have unique morphological traits, and show an extraordinary diversity in lifestyle. To improve our understanding of monocot origin and evolution, we generate chromosome-level reference genomes of the diploid Acorus gramineus and the tetraploid Ac. calamus, the only two accepted species from the family Acoraceae, which form a sister lineage to all other monocots. Comparing the genomes of Ac. gramineus and Ac. calamus, we suggest that Ac. gramineus is not a potential diploid progenitor of Ac. calamus, and Ac. calamus is an allotetraploid with two subgenomes A, and B, presenting asymmetric evolution and B subgenome dominance. Both the diploid genome of Ac. gramineus and the subgenomes A and B of Ac. calamus show clear evidence of whole-genome duplication (WGD), but Acoraceae does not seem to share an older WGD that is shared by most other monocots. We reconstruct an ancestral monocot karyotype and gene toolkit, and discuss scenarios that explain the complex history of the Acorus genome. Our analyses show that the ancestors of monocots exhibit mosaic genomic features, likely important for that appeared in early monocot evolution, providing fundamental insights into the origin, evolution, and diversification of monocots.

Whole-genome and dispersed duplication, including transposed duplication, jointly advance the evolution of TLP genes in seven representative Poaceae lineages

BACKGROUND

In the evolutionary study of gene families, exploring the duplication mechanisms of gene families helps researchers understand their evolutionary history. The tubby-like protein (TLP) family is essential for growth and development in plants and animals. Much research has been done on its function; however, limited information is available with regard to the evolution of the TLP gene family. Herein, we systematically investigated the evolution of TLP genes in seven representative Poaceae lineages.

RESULTS

Our research showed that the evolution of TLP genes was influenced not only by whole-genome duplication (WGD) and dispersed duplication (DSD) but also by transposed duplication (TRD), which has been neglected in previous research. For TLP family size, we found an evolutionary pattern of progressive shrinking in the grass family. Furthermore, the evolution of the TLP gene family was at least affected by evolutionary driving forces such as duplication, purifying selection, and base mutations.

CONCLUSIONS

This study presents the first comprehensive evolutionary analysis of the TLP gene family in grasses. We demonstrated that the TLP gene family is also influenced by a transposed duplication mechanism. Several new insights into the evolution of the TLP gene family are presented. This work provides a good reference for studying gene evolution and the origin of duplication.

A pan-grass transcriptome reveals patterns of cellular divergence in crops

Different plant species within the grasses were parallel targets of domestication, giving rise to crops with distinct evolutionary histories and traits1. Key traits that distinguish these species are mediated by specialized cell types2. Here we compare the transcriptomes of root cells in three grass species-Zea mays, Sorghum bicolor and Setaria viridis. We show that single-cell and single-nucleus RNA sequencing provide complementary readouts of cell identity in dicots and monocots, warranting a combined analysis. Cell types were mapped across species to identify robust, orthologous marker genes. The comparative cellular analysis shows that the transcriptomes of some cell types diverged more rapidly than those of others-driven, in part, by recruitment of gene modules from other cell types. The data also show that a recent whole-genome duplication provides a rich source of new, highly localized gene expression domains that favour fast-evolving cell types. Together, the cell-by-cell comparative analysis shows how fine-scale cellular profiling can extract conserved modules from a pan transcriptome and provide insight on the evolution of cells that mediate key functions in crops.

How to survive in the world's third poplar: Insights from the genome of the highest altitude woody plant, Hippophae tibetana (Elaeagnaceae)

Hippophae tibetana (Tibetan sea-buckthorn) is one of the highest distributed woody plants in the world (3,000-5,200 meters a.s.l.). It is characterized by adaptation to extreme environment and important economic values. Here, we combined PacBio Hifi platform and Hi-C technology to assemble a 1,452.75 Mb genome encoding 33,367 genes with a Contig N50 of 74.31 Mb, and inferred its sexual chromosome. Two Hippophae-specific whole-genome duplication events (18.7-21.2 million years ago, Ma; 28.6-32.4 Ma) and long terminal repeats retroelements (LTR-RTs) amplifications were detected. Comparing with related species at lower altitude, Ziziphus jujuba (<1, 700 meters a.s.l.), H. tibetana had some significantly rapid evolving genes involved in adaptation to high altitude habitats. However, comparing with Hippophae rhamnoides (<3, 700 meters a.s.l.), no rapid evolving genes were found except microtubule and microtubule-based process genes, H. tibetana has a larger genome, with extra 2, 503 genes (7.5%) and extra 680.46 Mb transposable elements (TEs) (46.84%). These results suggest that the changes in the copy number and regulatory pattern of genes play a more important role for H. tibetana adapting to more extreme and variable environments at higher altitude by more TEs and more genes increasing genome variability and expression plasticity. This suggestion was supported by two findings: nitrogen-fixing genes of H. tibetana having more copies, and intact TEs being significantly closer genes than fragmentary TEs. This study provided new insights into the evolution of alpine plants.

Genome sequencing reveals evidence of adaptive variation in the genus Zea

Maize is a globally valuable commodity and one of the most extensively studied genetic model organisms. However, we know surprisingly little about the extent and potential utility of the genetic variation found in wild relatives of maize. Here, we characterize a high-density genomic variation map from 744 genomes encompassing maize and all wild taxa of the genus Zea, identifying over 70 million single-nucleotide polymorphisms. The variation map reveals evidence of selection within taxa displaying novel adaptations. We focus on adaptive alleles in highland teosinte and temperate maize, highlighting the key role of flowering-time-related pathways in their adaptation. To show the utility of variants in these data, we generate mutant alleles for two flowering-time candidate genes. This work provides an extensive sampling of the genetic diversity of Zea, resolving questions on evolution and identifying adaptive variants for direct use in modern breeding.

Two independent allohexaploidizations and genomic fractionation in Solanales

Solanales, an order of flowering plants, contains the most economically important vegetables among all plant orders. To date, many Solanales genomes have been sequenced. However, the evolutionary processes of polyploidization events in Solanales and the impact of polyploidy on species diversity remain poorly understood. We compared two representative Solanales genomes (Solanum lycopersicum L. and Ipomoea triloba L.) and the Vitis vinifera L. genome and confirmed two independent polyploidization events. Solanaceae common hexaploidization (SCH) and Convolvulaceae common hexaploidization (CCH) occurred ∼43-49 and ∼40-46 million years ago (Mya), respectively. Moreover, we identified homologous genes related to polyploidization and speciation and constructed multiple genomic alignments with V. vinifera genome, providing a genomic homology framework for future Solanales research. Notably, the three polyploidization-produced subgenomes in both S. lycopersicum and I. triloba showed significant genomic fractionation bias, suggesting the allohexaploid nature of the SCH and CCH events. However, we found that the higher genomic fractionation bias of polyploidization-produced subgenomes in Solanaceae was likely responsible for their more abundant species diversity than that in Convolvulaceae. Furthermore, through genomic fractionation and chromosomal structural variation comparisons, we revealed the allohexaploid natures of SCH and CCH, both of which were formed by two-step duplications. In addition, we found that the second step of two paleohexaploidization events promoted the expansion and diversity of β-amylase (BMY) genes in Solanales. These current efforts provide a solid foundation for future genomic and functional exploration of Solanales.

Genome-Wide Investigation and Characterization of SWEET Gene Family with Focus on Their Evolution and Expression during Hormone and Abiotic Stress Response in Maize

The sugar will eventually be exported transporters (SWEET) family is an important group of transport carriers for carbon partitioning in plants and has important functions in growth, development, and abiotic stress tolerance. Although the SWEET family is an important sugar transporter, little is known of the functions of the SWEET family in maize (Zea mays), especially in response to abiotic stresses. To further explore the response pattern of maize SWEET to abiotic stress, a bioinformatics-based approach was used to predict and identify the maize SWEET gene (ZmSWEET) family. Twenty-four ZmSWEET genes were identified using the MaizeGDB database. Phylogenetic analysis resolved these twenty-four genes into four clades. One tandem and five segmental duplication events were identified, which played a major role in ZmSWEET family expansion. Synteny analysis provided insight into the evolutionary characteristics of the ZmSWEET genes with those of three graminaceous crop species. A heatmap showed that most ZmSWEET genes responded to at least one type of abiotic stress. By an abscisic acid signaling pathway, among which five genes were significantly induced under NaCl treatment, eight were obviously up-regulated under PEG treatment and five were up-regulated under Cd stress, revealing their potential functions in response to abiotic stress. These findings will help to explain the evolutionary links of the ZmSWEET family and contribute to future studies on the functional characteristics of ZmSWEET genes, and then improve abiotic stress tolerance in maize through molecular breeding.

Genomic insights into rapid speciation within the world's largest tree genus Syzygium

Species radiations, despite immense phenotypic variation, can be difficult to resolve phylogenetically when genetic change poorly matches the rapidity of diversification. Genomic potential furnished by palaeopolyploidy, and relative roles for adaptation, random drift and hybridisation in the apportionment of genetic variation, remain poorly understood factors. Here, we study these aspects in a model radiation, Syzygium, the most species-rich tree genus worldwide. Genomes of 182 distinct species and 58 unidentified taxa are compared against a chromosome-level reference genome of the sea apple, Syzygium grande. We show that while Syzygium shares an ancient genome doubling event with other Myrtales, little evidence exists for recent polyploidy events. Phylogenomics confirms that Syzygium originated in Australia-New Guinea and diversified in multiple migrations, eastward to the Pacific and westward to India and Africa, in bursts of speciation visible as poorly resolved branches on phylogenies. Furthermore, some sublineages demonstrate genomic clines that recapitulate cladogenetic events, suggesting that stepwise geographic speciation, a neutral process, has been important in Syzygium diversification.

Genome-wide identification and multiple abiotic stress transcript profiling of potassium transport gene homologs in Sorghum bicolor

Potassium (K⁺) is the most abundant cation that plays a crucial role in various cellular processes in plants. Plants have developed an efficient mechanism for the acquisition of K⁺ when grown in K⁺ deficient or saline soils. A total of 47 K⁺ transport gene homologs (27 HAKs, 4 HKTs, 2 KEAs, 9 AKTs, 2 KATs, 2 TPCs, and 1 VDPC) have been identified in Sorghum bicolor. Of 47 homologs, 33 were identified as K⁺ transporters and the remaining 14 as K⁺ channels. Chromosome 2 has been found as the hotspot of K⁺ transporters with 9 genes. Phylogenetic analysis revealed the conservation of sorghum K⁺ transport genes akin to Oryza sativa. Analysis of regulatory elements indicates the key roles that K⁺ transport genes play under different biotic and abiotic stress conditions. Digital expression data of different developmental stages disclosed that expressions were higher in milk, flowering, and tillering stages. Expression levels of the genes SbHAK27 and SbKEA2 were higher during milk, SbHAK17, SbHAK11, SbHAK18, and SbHAK7 during flowering, SbHAK18, SbHAK10, and 23 other gene expressions were elevated during tillering inferring the important role that K⁺ transport genes play during plant growth and development. Differential transcript expression was observed in different tissues like root, stem, and leaf under abiotic stresses such as salt, drought, heat, and cold stresses. Collectively, the in-depth genome-wide analysis and differential transcript profiling of K⁺ transport genes elucidate their role in ion homeostasis and stress tolerance mechanisms.

Wheat genomic study for genetic improvement of traits in China

Bread wheat (Triticum aestivum L.) is a major crop that feeds 40% of the world's population. Over the past several decades, advances in genomics have led to tremendous achievements in understanding the origin and domestication of wheat, and the genetic basis of agronomically important traits, which promote the breeding of elite varieties. In this review, we focus on progress that has been made in genomic research and genetic improvement of traits such as grain yield, end-use traits, flowering regulation, nutrient use efficiency, and biotic and abiotic stress responses, and various breeding strategies that contributed mainly by Chinese scientists. Functional genomic research in wheat is entering a new era with the availability of multiple reference wheat genome assemblies and the development of cutting-edge technologies such as precise genome editing tools, high-throughput phenotyping platforms, sequencing-based cloning strategies, high-efficiency genetic transformation systems, and speed-breeding facilities. These insights will further extend our understanding of the molecular mechanisms and regulatory networks underlying agronomic traits and facilitate the breeding process, ultimately contributing to more sustainable agriculture in China and throughout the world.

Integrative genomics analysis of the ever-shrinking pectin methylesterase (PME) gene family in foxtail millet (Setaria italica)

Pectin methylesterase (PME) plays a vital role in the growth and development of plants. Their genes can be classified into two types, with Type-1 having an extra domain, PMEI. PME genes in foxtail millet (Setaria italica L.) have not been identified, and their sequence features and evolution have not been explored. Here, we identified 41 foxtail millet PME genes. Decoding the pro-region, containing the PMEI domain, revealed its more active nature than the DNA encoding PME domain, easier to be lost to produce Type-2 PME genes. We inferred that the active nature of the pro-region could be related to its harbouring more repetitive DNA sequences. Further, we revealed that though whole-genome duplication and tandem duplication contributed to producing new copies of PME genes, phylogenetic analysis provided clear evidence of ever-shrinking gene family size in foxtail millet and the other grasses in the past 100 million years. Phylogenetic analysis also supports the existence of two gene groups, Group I and Group II, with genes in Group II being more conservative. Our research contributes to understanding how DNA sequence structure affects the functional innovation and evolution of PME genes.

GGDB: A Grameneae genome alignment database of homologous genes hierarchically related to evolutionary events

The genomes of Gramineae plants have been preferentially sequenced owing to their economic value. These genomes are often quite complex, for example harboring many duplicated genes, and are the main source of genetic innovation and often the result of recurrent polyploidization. Deciphering these complex genome structures and linking duplicated genes to specific polyploidization events are important for understanding the biology and evolution of plants. However, efforts have been hampered by the complexity of analyzing these genomes. Here, we analyzed 29 well-assembled and up-to-date Gramineae genome sequences by hierarchically relating duplicated genes in collinear regions to specific polyploidization or speciation events. We separated duplicated genes produced by each event, established lists of paralogous and orthologous genes, and ultimately constructed an online database, GGDB (http://www.grassgenome.com/). Homologous gene lists from each plant and between plants can be displayed, searched, and downloaded from the database. Interactive comparison tools are deployed to demonstrate homology among user-selected plants and to draw genome-scale or local alignment figures and gene-based phylogenetic trees corrected by exploiting gene collinearity. Using these tools and figures, users can easily detect structural changes in genomes and explore the effects of paleo-polyploidy on crop genome structure and function. The GGDB will provide a useful platform for improving our understanding of genome changes and functional innovation in Gramineae plants.

Pleiotropic ZmICE1 Is an Important Transcriptional Regulator of Maize Endosperm Starch Biosynthesis

Starch, the major component of cereal grains, affects crop yield and quality and is widely used in food and industrial applications. The biosynthesis of maize starch is a complex process involving a series of functional enzymes. However, the sophisticated regulatory mechanisms of starch biosynthetic genes have not been fully elaborated. The basic/helix-loop-helix (bHLH) transcription factors are widely distributed in eukaryotes and participate in many physiological processes. In this study, 202 bHLH encoding genes were identified in the maize genome by Blast method. ZmICE1 gene, which belongs to the ICE subfamily of the bHLH family, was obtained and expressed mainly in maize filling endosperm and co-expressed with 14 starch biosynthesis genes. Based on the comparative analyses across different plant species, we revealed that the gene structures and protein domains of the ICE subfamily were conserved between monocots and dicots, suggesting their functional conservation feature. Yeast activation and subcellular localization assays suggested that ZmICE1 had transcriptional activation activity and localized in the nucleus. Yeast one-hybrid assays confirmed that ZmICE1 could directly bind to the promoters of ZmSSIIa and ZmGBSSI. Transient gene expression analysis in maize endosperm revealed that ZmICE1 positively regulated the expression of ZmSSIIa, but inhibited the expression of ZmGBSSI. Our results indicated that ZmICE1 could function as a regulator of maize starch biosynthesis.

Evolution and origin of bread wheat

Bread wheat (Triticum aestivum, genome BBAADD) is a young hexaploid species formed only 8,500-9,000 years ago through hybridization between a domesticated free-threshing tetraploid progenitor, genome BBAA, and Aegilops tauschii, the diploid donor of the D subgenome. Very soon after its formation, it spread globally from its cradle in the fertile crescent into new habitats and climates, to become a staple food of humanity. This extraordinary global expansion was probably enabled by allopolyploidy that accelerated genetic novelty through the acquisition of new traits, new intergenomic interactions, and buffering of mutations, and by the attractiveness of bread wheat's large, tasty, and nutritious grain with high baking quality. New genome sequences suggest that the elusive donor of the B subgenome is a distinct (unknown or extinct) species rather than a mosaic genome. We discuss the origin of the diploid and tetraploid progenitors of bread wheat and the conflicting genetic and archaeological evidence on where it was formed and which species was its free-threshing tetraploid progenitor. Wheat experienced many environmental changes throughout its evolution, therefore, while it might adapt to current climatic changes, efforts are needed to better use and conserve the vast gene pool of wheat biodiversity on which our food security depends.

The long road of functional recruitment-The evolution of a gene duplicate to pyrrolizidine alkaloid biosynthesis in the morning glories (Convolvulaceae)

In plants, homospermidine synthase (HSS) is a pathway-specific enzyme initiating the biosynthesis of pyrrolizidine alkaloids (PAs), which function as a chemical defense against herbivores. In PA-producing Convolvulaceae ("morning glories"), HSS originated from deoxyhypusine synthase at least >50 to 75 million years ago via a gene duplication event and subsequent functional diversification. To study the recruitment of this ancient gene duplicate to PA biosynthesis, the presence of putative hss gene copies in 11 Convolvulaceae species was analyzed. Additionally, various plant parts from seven of these species were screened for the presence of PAs. Although all of these species possess a putative hss copy, PAs could only be detected in roots of Ipomoea neei (Spreng.) O'Donell and Distimake quinquefolius (L.) A.R.Simões & Staples in this study. A precursor of PAs was detected in roots of Ipomoea alba L. Thus, despite sharing high sequence identities, the presence of an hss gene copy does not correlate with PA accumulation in particular species of Convolvulaceae. In vitro activity assays of the encoded enzymes revealed a broad spectrum of enzyme activity, further emphasizing a functional diversity of the hss gene copies. A recently identified HSS specific amino acid motif seems to be important for the loss of the ancestral protein function-the activation of the eukaryotic initiation factor 5A (eIF5A). Thus, the motif might be indicative for a change of function but allows not to predict the new function. This emphasizes the challenges in annotating functions for duplicates, even for duplicates from closely related species.

Conversion between duplicated genes generated by polyploidization contributes to the divergence of poplar and willow

BACKGROUND

Gene conversion has an important effect on duplicate genes produced by polyploidization. Poplar (Populus trichocarpa) and willow (Salix brachista) are leading models and excellent green plants in the Salicaceae. Although much attention has been paid to the evolution of duplicated genes in poplar and willow, the role of conversion between duplicates generated from polyploidization remains poorly understood.

RESULTS

Here, through genomic synteny analyses, we identified duplicate genes generated by the Salicaceae common tetraploidization (SCT) in the poplar and willow genomes. We estimated that at least 0.58% and 0.25% of poplar and willow duplicates were affected by whole-gene conversion after the poplar-willow divergence, with more (5.73% and 2.66%) affected by partial-gene conversion. Moreover, we found that the converted duplicated genes were unevenly distributed on each chromosome in the two genomes, and the well-preserved homoeologous chromosome regions may facilitate the conversion of duplicates. Notably, we found that conversion maintained the similarity of duplicates, likely contributing to the conservation of certain sequences, but is essentially accelerated the rate of evolution and increased species divergence. In addition, we found that converted duplicates tended to have more similar expression patterns than nonconverted duplicates. We found that genes associated with multigene families were preferentially converted. We also found that the genes encoding conserved structural domains associated with specific traits exhibited a high frequency of conversion.

CONCLUSIONS

Extensive conversion between duplicate genes generated from the SCT contributes to the diversification of the family Salicaceae and has had long-lasting effects on those genes with important biological functions.

Exploring the SiCCT Gene Family and Its Role in Heading Date in Foxtail Millet

CCT transcription factors are involved in the regulation of photoperiod and abiotic stress in Arabidopsis and rice. It is not clear that how CCT gene family expand and regulate heading date in foxtail millet. In this study, we conducted a systematic analysis of the CCT gene family in foxtail millet. Thirty-nine CCT genes were identified and divided into four subfamilies based on functional motifs. Analysis showed that dispersed duplication played a predominant role in the expansion of CCT genes during evolution. Nucleotide diversity analysis suggested that genes in CONSTANS (COL)-like, CCT MOTIF FAMILY (CMF)-like, and pseudoresponse response regulator (PRR)-like subfamilies were subjected to selection. Fifteen CCT genes were colocalized with previous heading date quantitative trait loci (QTL) and genome-wide association analysis (GWAS) signals. Transgenic plants were then employed to confirm that overexpression of the CCT gene SiPRR37 delayed the heading date and increased plant height. Our study first investigated the characterization and expansion of the CCT family in foxtail millet and demonstrated the role of SiPRR37. These results lay a significant foundation for further research on the function of CCT genes and provide a cue for the regulation of heading date.

Paleopolyploidies and Genomic Fractionation in Major Eudicot Clades

Eudicots account for ~75% of living angiosperms, containing important food and energy crops. Recently, high-quality genome sequences of several eudicots including Aquilegia coerulea and Nelumbo nucifera have become available, providing an opportunity to investigate the early evolutionary characteristics of eudicots. We performed genomic hierarchical and event-related alignments to infer homology within and between representative species of eudicots. The results provide strong evidence for multiple independent polyploidization events during the early diversification of eudicots, three of which are likely to be allopolyploids: The core eudicot-common hexaploidy (ECH), Nelumbo-specific tetraploidy (NST), and Ranunculales-common tetraploidy (RCT). Using different genomes as references, we constructed genomic alignment to list the orthologous and paralogous genes produced by polyploidization and speciation. This could provide a fundamental framework for studying other eudicot genomes and gene(s) evolution. Further, we revealed significantly divergent evolutionary rates among these species. By performing evolutionary rate correction, we dated RCT to be ~118-134 million years ago (Mya), after Ranunculales diverged with core eudicots at ~123-139 Mya. Moreover, we characterized genomic fractionation resulting from gene loss and retention after polyploidizations. Notably, we revealed a high degree of divergence between subgenomes. In particular, synonymous nucleotide substitutions at synonymous sites (Ks) and phylogenomic analyses implied that A. coerulea might provide the subgenome(s) for the gamma-hexaploid hybridization.

Structural diversity, biosynthesis, and function of plant falcarin-type polyacetylenic lipids

The polyacetylenic lipids falcarinol, falcarindiol, and associated derivatives, termed falcarins, have a widespread taxonomical distribution in the plant kingdom and have received increasing interest for their demonstrated health-promoting properties as anti-cancer and anti-inflammatory agents. These fatty acid-derived compounds are also linked to plant pathogen resistance through their potent antimicrobial properties. Falcarin-type polyacetylenes, which contain two conjugated triple bonds, are derived from structural modifications of the common fatty acid oleic acid. In the past half century, much progress has been made in understanding the structural diversity of falcarins in the plant kingdom, whereas limited progress has been made on elucidating falcarin function in plant-pathogen interactions. More recently, an understanding of the biosynthetic machinery underlying falcarin biosynthesis has emerged. This review provides a concise summary of the current state of knowledge on falcarin structural diversity, biosynthesis, and plant defense properties. We also present major unanswered questions about falcarin biosynthesis and function.

Identification and Characterization of SOG1 (Suppressor of Gamma-Response 1) Homologues in Plants Using Data Mining Resources and Gene Expression Profiling

SOG1 (Suppressor of the Gamma response 1) is the master-regulator of plant DNA damage response (DDR), a highly coordinated network of DNA damage sensors, transducers, mediators, and effectors, with highly coordinated activities. SOG1 transcription factor belongs to the NAC/NAM protein family, containing the well-conserved NAC domain and five serine-glutamine (SQ) motifs, preferential targets for phosphorylation by ATM and ATR. So far, the information gathered for the SOG1 function comes from studies on the model plant Arabidopsis thaliana. To expand the knowledge on plant-specific DDR, it is opportune to gather information on other SOG1 orthologues. The current study identified plants where multiple SOG1 homologues are present and evaluated their functions by leveraging the information contained in publicly available transcriptomics databases. This analysis revealed the presence of multiple SOG1 sequences in thirteen plant species, and four (Medicago truncatula, Glycine max, Kalankoe fedtschenkoi, Populus trichocarpa) were selected for gene expression data mining based on database availability. Additionally, M. truncatula seeds and seedlings exposed to treatments known to activate DDR pathways were used to evaluate the expression profiles of MtSOG1a and MtSOG1b. The experimental workflow confirmed the data retrieved from transcriptomics datasets, suggesting that the SOG1 homologues have redundant functions in different plant species.

Reshuffling of the ancestral core-eudicot genome shaped chromatin topology and epigenetic modification in Panax

All extant core-eudicot plants share a common ancestral genome that has experienced cyclic polyploidizations and (re)diploidizations. Reshuffling of the ancestral core-eudicot genome generates abundant genomic diversity, but the role of this diversity in shaping the hierarchical genome architecture, such as chromatin topology and gene expression, remains poorly understood. Here, we assemble chromosome-level genomes of one diploid and three tetraploid Panax species and conduct in-depth comparative genomic and epigenomic analyses. We show that chromosomal interactions within each duplicated ancestral chromosome largely maintain in extant Panax species, albeit experiencing ca. 100–150 million years of evolution from a shared ancestor. Biased genetic fractionation and epigenetic regulation divergence during polyploidization/(re)diploidization processes generate remarkable biochemical diversity of secondary metabolites in the Panax genus. Our study provides a paleo-polyploidization perspective of how reshuffling of the ancestral core-eudicot genome leads to a highly dynamic genome and to the metabolic diversification of extant eudicot plants.

The role of polyploidization generated genomic diversity in shaping the hierarchical genome architecture remains unclear. Here, the authors show that repatterning of the ancestral eudicot genome has resulted in multi-dimensional genome plasticity and secondary metabolite diversification via comparisons of Panax genomes.

A recent burst of gene duplications in Triticeae

Gene duplication provides raw genetic materials for evolution and potentially novel genes for crop improvement. The two seminal genomic studies of Aegilops tauschii both mentioned the large number of genes independently duplicated in recent years, but the duplication mechanism and the evolutionary significance of these gene duplicates have not yet been investigated. Here, we found that a recent burst of gene duplications (hereafter abbreviated as the RBGD) has probably occurred in all sequenced Triticeae species. Further investigations of the characteristics of the gene duplicates and their flanking sequences suggested that transposable element (TE) activity may have been involved in generating the RBGD. We also characterized the duplication timing, retention pattern, diversification, and expression of the duplicates following the evolution of Triticeae. Multiple subgenome-specific comparisons of the duplicated gene pairs clearly supported extensive differential regulation and related functional diversity among such pairs in the three subgenomes of bread wheat. Moreover, several duplicated genes from the RBGD have evolved into key factors that influence important agronomic traits of wheat. Our results provide insights into a unique source of gene duplicates in Triticeae species, which has increased the gene dosage together with the two polyploidization events in the evolutionary history of wheat.

Identification, Molecular Characteristics, and Evolution of GRF Gene Family in Foxtail Millet (Setaria italica L.)

Growth-regulating factor (GRF) is a multigene family that plays a vital role in the growth and development of plants. In the past, the GRF family of many plants has been studied. However, there is not a report about identification and evolution of GRF in foxtail millet (Setaria italia). Here, we identified 10 GRF genes in foxtail millet. Seven (70.00%) were regulated by Sit-miR396, and there were 19 optimal codons in GRFs of foxtail millet. Additionally, we found that WGD or segmental duplication have affected GRFs in foxtail millet between 15.07 and 45.97 million years ago. Regarding the GRF gene family of land plants, we identified a total of 157 GRF genes in 15 representative land plants. We found that GRF gene family originated from Group E, and the GRF gene family in monocots was gradually shrinking. Also, more loss resulted from the small number of GRF genes in lower plants. Exploring the evolution of GRF and functional analysis in the foxtail millet help us to understand GRF better and make a further study about the mechanism of GRF. These results provide a basis for the genetic improvement of foxtail millet and indicate an improvement of the yield.

The genome sequence provides insights into salt tolerance of Achnatherum splendens (Gramineae), a constructive species of alkaline grassland

Achnatherum splendens Trin. (Gramineae) is a constructive species of the arid grassland ecosystem in Northwest China and is a major forage grass. It has good tolerance of salt and drought stress in alkaline habitats. Here, we report its chromosome-level genome, determined through a combination of Illumina HiSeq sequencing, PacBio sequencing and Hi-C technology. The final assembly of the ~1.17 Gb genome sequence had a super-scaffold N50 of 40.3 Mb. A total of 57 374 protein-coding genes were annotated, of which 54 426 (94.5%) genes have functional protein annotations. Approximately 735 Mb (62.37%) of the assembly were identified as repetitive elements, and among these, LTRs (40.53%) constitute the highest proportion, having made a major contribution to the expansion of genome size in A. splendens. Phylogenetic analysis revealed that A. splendens diverged from the Brachypodium distachyon-Hordeum vulgare-Aegilops tauschii subclade around 37 million years ago (Ma) and that a clade comprising these four species diverged from the Phyllostachys edulis clade ~47 Ma. Genomic synteny indicates that A. splendens underwent an additional species-specific whole-genome duplication (WGD) 18-20 Ma, which further promoted an increase in copies of numerous saline-alkali-related gene families in the A. splendens genome. By transcriptomic analysis, we further found that many of these duplicated genes from this extra WGD exhibited distinct functional divergence in response to salt stress. This WGD, therefore, contributed to the strong resistance to salt stress and widespread arid adaptation of A. splendens.

Chromosome-length genome assemblies of six legume species provide insights into genome organization, evolution, and agronomic traits for crop improvement

INTRODUCTION

Legume crops are an important source of protein and oil for human health and in fixing atmospheric N₂ for soil enrichment. With an objective to accelerate much-needed genetic analyses and breeding applications, draft genome assemblies were generated in several legume crops; many of them are not high quality because they are mainly based on short reads. However, the superior quality of genome assembly is crucial for a detailed understanding of genomic architecture, genome evolution, and crop improvement.

OBJECTIVES

Present study was undertaken with an objective of developing improved chromosome-length genome assemblies in six different legumes followed by their systematic investigation to unravel different aspects of genome organization and legume evolution.

METHODS

We employed in situ Hi-C data to improve the existing draft genomes and performed different evolutionary and comparative analyses using improved genome assemblies.

RESULTS

We have developed chromosome-length genome assemblies in chickpea, pigeonpea, soybean, subterranean clover, and two wild progenitor species of cultivated groundnut (A. duranensis and A. ipaensis). A comprehensive comparative analysis of these genome assemblies offered improved insights into various evolutionary events that shaped the present-day legume species. We highlighted the expansion of gene families contributing to unique traits such as nodulation in legumes, gravitropism in groundnut, and oil biosynthesis in oilseed legume crops such as groundnut and soybean. As examples, we have demonstrated the utility of improved genome assemblies for enhancing the resolution of "QTL-hotspot" identification for drought tolerance in chickpea and marker-trait associations for agronomic traits in pigeonpea through genome-wide association study. Genomic resources developed in this study are publicly available through an online repository, 'Legumepedia'.

CONCLUSION

This study reports chromosome-length genome assemblies of six legume species and demonstrates the utility of these assemblies in crop improvement. The genomic resources developed here will have significant role in accelerating genetic improvement applications of legume crops.

Integrated Genomic Analyses From Low-Depth Sequencing Help Resolve Phylogenetic Incongruence in the Bamboos (Poaceae: Bambusoideae)

The bamboos (Bambusoideae, Poaceae) comprise a major grass lineage with a complex evolutionary history involving ancient hybridization and allopolyploidy. About 1700 described species are classified into three tribes, Olyreae (herbaceous bamboos), Bambuseae (tropical woody bamboos), and Arundinarieae (temperate woody bamboos). Nuclear analyses strongly support monophyly of the woody tribes, whereas plastome analyses strongly support paraphyly, with Bambuseae sister to Olyreae. Our objectives were to clarify the origin(s) of the woody bamboo tribes and resolve the nuclear vs. plastid conflict using genomic tools. For the first time, plastid and nuclear genomic information from the same bamboo species were combined in a single study. We sampled 51 species of bamboos representing the three tribes, estimated their genome sizes and generated low-depth sample sequence data, from which plastomes were assembled and nuclear repeats were analyzed. The distribution of repeat families was found to agree with nuclear gene phylogenies, but also provides novel insights into nuclear evolutionary history. We infer two early, independent hybridization events, one between an Olyreae ancestor and a woody ancestor giving rise to the two Bambuseae lineages, and another between two woody ancestors giving rise to the Arundinarieae. Retention of the Olyreae plastome associated with differential dominance of nuclear genomes and subsequent diploidization in some lineages explains the paraphyly observed in plastome phylogenetic estimations. We confirm ancient hybridization and allopolyploidy in the origins of the extant woody bamboo lineages and propose biased fractionation and diploidization as important factors in their evolution.

Unraveling the Genome of a High Yielding Colombian Sugarcane Hybrid

Recent developments in High Throughput Sequencing (HTS) technologies and bioinformatics, including improved read lengths and genome assemblers allow the reconstruction of complex genomes with unprecedented quality and contiguity. Sugarcane has one of the most complicated genomes among grassess with a haploid length of 1Gbp and a ploidies between 8 and 12. In this work, we present a genome assembly of the Colombian sugarcane hybrid CC 01-1940. Three types of sequencing technologies were combined for this assembly: PacBio long reads, Illumina paired short reads, and Hi-C reads. We achieved a median contig length of 34.94 Mbp and a total genome assembly of 903.2 Mbp. We annotated a total of 63,724 protein coding genes and performed a reconstruction and comparative analysis of the sucrose metabolism pathway. Nucleotide evolution measurements between orthologs with close species suggest that divergence between Saccharum officinarum and Saccharum spontaneum occurred <2 million years ago. Synteny analysis between CC 01-1940 and the S. spontaneum genome confirms the presence of translocation events between the species and a random contribution throughout the entire genome in current sugarcane hybrids. Analysis of RNA-Seq data from leaf and root tissue of contrasting sugarcane genotypes subjected to water stress treatments revealed 17,490 differentially expressed genes, from which 3,633 correspond to genes expressed exclusively in tolerant genotypes. We expect the resources presented here to serve as a source of information to improve the selection processes of new varieties of the breeding programs of sugarcane.

De novo assembly, annotation, and comparative analysis of 26 diverse maize genomes

We report de novo genome assemblies, transcriptomes, annotations, and methylomes for the 26 inbreds that serve as the founders for the maize nested association mapping population. The number of pan-genes in these diverse genomes exceeds 103,000, with approximately a third found across all genotypes. The results demonstrate that the ancient tetraploid character of maize continues to degrade by fractionation to the present day. Excellent contiguity over repeat arrays and complete annotation of centromeres revealed additional variation in major cytological landmarks. We show that combining structural variation with single-nucleotide polymorphisms can improve the power of quantitative mapping studies. We also document variation at the level of DNA methylation and demonstrate that unmethylated regions are enriched for cis-regulatory elements that contribute to phenotypic variation.

Comparisons within the Rice GA 2-Oxidase Gene Family Revealed Three Dominant Paralogs and a Functional Attenuated Gene that Led to the Identification of Four Amino Acid Variants Associated with GA Deactivation Capability

BACKGROUND

GA 2-oxidases (GA2oxs) are involved in regulating GA homeostasis in plants by inactivating bioactive GAs through 2β-hydroxylation. Rice GA2oxs are encoded by a family of 10 genes; some of them have been characterized, but no comprehensive comparisons for all these genes have been conducted.

RESULTS

Rice plants with nine functional GA2oxs were demonstrated in the present study, and these genes not only were differentially expressed but also revealed various capabilities for GA deactivation based on their height-reducing effects in transgenic plants. Compared to that of wild-type plants, the relative plant height (RPH) of transgenic plants was scored to estimate their reducing effects, and 8.3% to 59.5% RPH was observed. Phylogenetic analysis of class I GA2ox genes revealed two functionally distinct clades in the Poaceae. The OsGA2ox3, 4, and 8 genes belonging to clade A showed the most severe effect (8.3% to 8.7% RPH) on plant height reduction, whereas the OsGA2ox7 gene belonging to clade B showed the least severe effect (59.5% RPH). The clade A OsGA2ox3 gene contained two conserved C186/C194 amino acids that were crucial for enzymatic activity. In the present study, these amino acids were replaced with OsGA2ox7-conserved arginine (C186R) and proline (C194P), respectively, or simultaneously (C186R/C194P) to demonstrate their importance in planta. Another two amino acids, Q220 and Y274, conserved in OsGA2ox3 were substituted with glutamic acid (E) and phenylalanine (F), respectively, or simultaneously to show their significance in planta. In addition, through sequence divergence, RNA expression profile and GA deactivation capability analyses, we proposed that OsGA2ox1, OsGA2ox3 and OsGA2ox6 function as the predominant paralogs in each of their respective classes.

CONCLUSIONS

This study demonstrates rice has nine functional GA2oxs and the class I GA2ox genes are divided into two functionally distinct clades. Among them, the OsGA2ox7 of clade B is a functional attenuated gene and the OsGA2ox1, OsGA2ox3 and OsGA2ox6 are the three predominant paralogs in the family.

Chromosome number is key to longevity of polyploid lineages

Polyploidy is ubiquitous and often recursive in plant lineages, most frequently resulting in extinction but occasionally associated with great evolutionary success. However, instead of chromosome numbers exponentially increasing due to recurrent polyploidy, most angiosperm species have fewer than 14 chromosome pairs. Following genome duplication, diploidisation can render one copy of essential genes nonfunctional without fitness cost. In isolated subpopulations, alternate (homoeologous) gene copies can be lost, creating incompatibilities that reduce fitness of hybrids between subpopulations, constraining exchange of favourable genetic changes and reducing species fitness. When multiple sets of incompatible genes are genetically linked, their deleterious effects are not independent. The effective number of independently acting sets of incompatible loci in hybrids is limited by chromosome number and recombination. Therefore, species with many chromosomes are subject to a higher fitness penalty during diploidisation. Karyotypic changes, especially fusions, that reduce gene flow are normally fitness disadvantages, but during the diploidisation process, can increase fitness by reducing mixing of differentially diploidised alleles. Fitness penalties caused by diploidisation favour accelerated karyotypic change, with each change increasing barriers to gene flow, contributing to speciation. Lower chromosome numbers and increased chromosome fusions confer advantages to surviving the diploidisation process following polyploid formation, by independent mechanisms.

Conversion between 100-million-year-old duplicated genes contributes to rice subspecies divergence

BACKGROUND

Duplicated gene pairs produced by ancient polyploidy maintain high sequence similarity over a long period of time and may result from illegitimate recombination between homeologous chromosomes. The genomes of Asian cultivated rice Oryza sativa ssp. indica (XI) and Oryza sativa ssp. japonica (GJ) have recently been updated, providing new opportunities for investigating ongoing gene conversion events and their impact on genome evolution.

RESULTS

Using comparative genomics and phylogenetic analyses, we evaluated gene conversion rates between duplicated genes produced by polyploidization 100 million years ago (mya) in GJ and XI. At least 5.19-5.77% of genes duplicated across the three rice genomes were affected by whole-gene conversion after the divergence of GJ and XI at ~ 0.4 mya, with more (7.77-9.53%) showing conversion of only portions of genes. Independently converted duplicates surviving in the genomes of different subspecies often use the same donor genes. The ongoing gene conversion frequency was higher near chromosome termini, with a single pair of homoeologous chromosomes, 11 and 12, in each rice genome being most affected. Notably, ongoing gene conversion has maintained similarity between very ancient duplicates, provided opportunities for further gene conversion, and accelerated rice divergence. Chromosome rearrangements after polyploidization are associated with ongoing gene conversion events, and they directly restrict recombination and inhibit duplicated gene conversion between homeologous regions. Furthermore, we found that the converted genes tended to have more similar expression patterns than nonconverted duplicates. Gene conversion affects biological functions associated with multiple genes, such as catalytic activity, implying opportunities for interaction among members of large gene families, such as NBS-LRR disease-resistance genes, contributing to the occurrence of the gene conversion.

CONCLUSION

Duplicated genes in rice subspecies generated by grass polyploidization ~ 100 mya remain affected by gene conversion at high frequency, with important implications for the divergence of rice subspecies.

The Pharus latifolius genome bridges the gap of early grass evolution

The grass family (Poaceae) includes all commercial cereal crops and is a major contributor to biomass in various terrestrial ecosystems. The ancestry of all grass genomes includes a shared whole-genome duplication (WGD), named rho (ρ) WGD, but the evolutionary significance of ρ-WGD remains elusive. We sequenced the genome of Pharus latifolius, a grass species (producing a true spikelet) in the subfamily Pharoideae, a sister lineage to the core Poaceae including the (Panicoideae, Arundinoideae, Chloridoideae, Micrairoideae, Aristidoideae, and Danthonioideae (PACMAD) and Bambusoideae, Oryzoideae, and Pooideae (BOP) clades. Our results indicate that the P. latifolius genome has evolved slowly relative to cereal grass genomes, as reflected by moderate rates of molecular evolution, limited chromosome rearrangements and a low rate of gene loss for duplicated genes. We show that the ρ-WGD event occurred approximately 98.2 million years ago (Ma) in a common ancestor of the Pharoideae and the PACMAD and BOP grasses. This was followed by contrasting patterns of diploidization in the Pharus and core Poaceae lineages. The presence of two FRIZZY PANICLE-like genes in P. latifolius, and duplicated MADS-box genes, support the hypothesis that the ρ-WGD may have played a role in the origin and functional diversification of the spikelet, an adaptation in grasses related directly to cereal yields. The P. latifolius genome sheds light on the origin and early evolution of grasses underpinning the biology and breeding of cereals.

The celery genome sequence reveals sequential paleo-polyploidizations, karyotype evolution and resistance gene reduction in apiales

Celery (Apium graveolens L. 2n = 2x = 22), a member of the Apiaceae family, is among the most important and globally grown vegetables. Here, we report a high-quality genome sequence assembly, anchored to 11 chromosomes, with total length of 3.33 Gb and N50 scaffold length of 289.78 Mb. Most (92.91%) of the genome is composed of repetitive sequences, with 62.12% of 31 326 annotated genes confined to the terminal 20% of chromosomes. Simultaneous bursts of shared long-terminal repeats (LTRs) in different Apiaceae plants suggest inter-specific exchanges. Two ancestral polyploidizations were inferred, one shared by Apiales taxa and the other confined to Apiaceae. We reconstructed 8 Apiales proto-chromosomes, inferring their evolutionary trajectories from the eudicot common ancestor to extant plants. Transcriptome sequencing in three tissues (roots, leaves and petioles), and varieties with different-coloured petioles, revealed 4 and 2 key genes in pathways regulating anthocyanin and coumarin biosynthesis, respectively. A remarkable paucity of NBS disease-resistant genes in celery (62) and other Apiales was explained by extensive loss and limited production of these genes during the last ~10 million years, raising questions about their biotic defence mechanisms and motivating research into effects of chemicals, for example coumarins, that give off distinctive odours. Celery genome sequencing and annotation facilitates further research into important gene functions and breeding, and comparative genomic analyses in Apiales.

Quantitative trait mapping of plant architecture in two BC1F2 populations of Sorghum Bicolor × S. halepense and comparisons to two other sorghum populations

Comparing populations derived, respectively, from polyploid Sorghum halepense and its progenitors improved knowledge of plant architecture and showed that S. halepense harbors genetic novelty of potential value for sorghum improvement Vegetative growth and the timing of the vegetative-to-reproductive transition are critical to a plant's fitness, directly and indirectly determining when and how a plant lives, grows and reproduces. We describe quantitative trait analysis of plant height and flowering time in the naturally occurring tetraploid Sorghum halepense, using two novel BC₁F₂ populations totaling 246 genotypes derived from backcrossing two tetraploid Sorghum bicolor x S. halepense F₁ plants to a tetraploidized S. bicolor. Phenotyping for two years each in Bogart, GA and Salina, KS allowed us to dissect variance into narrow-sense genetic (QTLs) and environmental components. In crosses with a common S. bicolor BTx623 parent, comparison of QTLs in S. halepense, its rhizomatous progenitor S. propinquum and S. bicolor race guinea which is highly divergent from BTx623 permit inferences of loci at which new alleles have been associated with improvement of elite sorghums. The relative abundance of QTLs unique to the S. halepense populations may reflect its polyploidy and subsequent 'diploidization' processes often associated with the formation of genetic novelty, a possibility further supported by a high level of QTL polymorphism within sibling lines derived from a common S. halepense parent. An intriguing hypothesis for further investigation is that polyploidy of S. halepense following 96 million years of abstinence, coupled with natural selection during its spread to diverse environments across six continents, may provide a rich collection of novel alleles that offer potential opportunities for sorghum improvement.