Distinct Expression and Methylation Patterns for Genes with Different Fates following a Single Whole-Genome Duplication in Flowering Plants

NnSUS1 encodes a sucrose synthase involved in sugar accumulation in lotus seed cotyledons

Fresh lotus seeds are gaining favor with consumers for their crunchy texture and natural sweetness. However, the intricacies of sugar accumulation in lotus seeds remain elusive, which greatly hinders the quality improvement of fresh lotus seeds. This study endeavors to elucidate this mechanism by identifying and characterizing the sucrose synthase (SUS) gene family in lotus. Comprising five distinct members, namely NnSUS1 to NnSUS5, each gene within this family features a C-terminal glycosyl transferase1 (GT1) domain. Among them, NnSUS1 is the predominately expressed gene, showing high transcript abundance in the floral organs and cotyledons. NnSUS1 was continuously up-regulated from 6 to 18 days after pollination (DAP) in lotus cotyledons. Furthermore, NnSUS1 demonstrates co-expression relationships with numerous genes involved in starch and sucrose metabolism. To investigate the function of NnSUS1, a transient overexpression system was established in lotus cotyledons, which confirmed the gene's contribution to sugar accumulation. Specifically, transient overexpression of NnSUS1 in seed cotyledons leads to a significant increase in the levels of total soluble sugar, including sucrose and fructose. These findings provide valuable theoretical insights for improving sugar content in lotus seeds through molecular breeding methods.

Lineage-Specific CYP80 Expansion and Benzylisoquinoline Alkaloid Diversity in Early-Diverging Eudicots

Menispermaceae species, as early-diverging eudicots, can synthesize valuable benzylisoquinoline alkaloids (BIAs) like bisbenzylisoquinoline alkaloids (bisBIAs) and sinomenines with a wide range of structural diversity. However, the evolutionary mechanisms responsible for their chemo-diversity are not well understood. Here, a chromosome-level genome assembly of Menispermum dauricum is presented and demonstrated the occurrence of two whole genome duplication (WGD) events that are shared by Ranunculales and specific to Menispermum, providing a model for understanding chromosomal evolution in early-diverging eudicots. The biosynthetic pathway for diverse BIAs in M. dauricum is reconstructed by analyzing the transcriptome and metabolome. Additionally, five catalytic enzymes - one norcoclaurine synthase (NCS) and four cytochrome P450 monooxygenases (CYP450s) - from M. dauricum are responsible for the formation of the skeleton, hydroxylated modification, and C-O/C-C phenol coupling of BIAs. Notably, a novel leaf-specific MdCYP80G10 enzyme that catalyzes C2'-C4a phenol coupling of (S)-reticuline into sinoacutine, the enantiomer of morphinan compounds, with predictable stereospecificity is discovered. Moreover, it is found that Menispermum-specific CYP80 gene expansion, as well as tissue-specific expression, has driven BIA diversity in Menispermaceae as compared to other Ranunculales species. This study sheds light on WGD occurrences in early-diverging eudicots and the evolution of diverse BIA biosynthesis.

Reshaped DNA methylation cooperating with homoeolog-divergent expression promotes improved root traits in synthesized tetraploid wheat

Polyploidization is a major event driving plant evolution and domestication. However, how reshaped epigenetic modifications coordinate gene transcription to generate phenotypic variations during wheat polyploidization is currently elusive. Here, we profiled transcriptomes and DNA methylomes of two diploid wheat accessions (SlSl and AA) and their synthetic allotetraploid wheat line (SlSlAA), which displayed elongated root hair and improved root capability for nitrate uptake and assimilation after tetraploidization. Globally decreased DNA methylation levels with a reduced difference between subgenomes were observed in the roots of SlSlAA. DNA methylation changes in first exon showed strong connections with altered transcription during tetraploidization. Homoeolog-specific transcription was associated with biased DNA methylation as shaped by homoeologous sequence variation. The hypomethylated promoters showed significantly enriched binding sites for MYB, which may affect gene transcription in response to root hair growth. Two master regulators in root hair elongation pathway, AlCPC and TuRSL4, exhibited upregulated transcription levels accompanied by hypomethylation in promoter, which may contribute to the elongated root hair. The upregulated nitrate transporter genes, including NPFs and NRTs, also are significantly associated with hypomethylation, indicating an epigenetic-incorporated regulation manner in improving nitrogen use efficiency. Collectively, these results provided new insights into epigenetic changes in response to crop polyploidization and underscored the importance of epigenetic regulation in improving crop traits.

Genome-Wide Identification and Co-Expression Networks of WOX Gene Family in Nelumbo nucifera

WUSCHEL-related homeobox (WOX) genes are a class of plant-specific transcription factors, regulating the development of multiple tissues. However, the genomic characterizations and expression patterns of WOX genes have not been analyzed in lotus. In this study, 15 NnWOX genes were identified based on the well-annotated reference genome of lotus. According to the phylogenetic analysis, the NnWOX genes were clustered into three clades, i.e., ancient clade, intermediate clade, and WUS clade. Except for the conserved homeobox motif, we further found specific motifs of NnWOX genes in different clades and divergence gene structures, suggesting their distinct functions. In addition, two NnWOX genes in the ancient clade have conserved expression patterns and other NnWOX genes exhibit different expression patterns in lotus tissues, suggesting a low level of functional redundancy in lotus WOX genes. Furthermore, we constructed the gene co-expression networks for each NnWOX gene. Based on weighted gene co-expression network analysis (WGCNA), ten NnWOX genes and their co-expressed genes were assigned to the modules that were significantly related to the cotyledon and seed coat. We further performed RT-qPCR experiments, validating the expression levels of ten NnWOX genes in the co-expression networks. Our study reveals comprehensive genomic features of NnWOX genes in lotus, providing a solid basis for further function studies.

Genomic evidence for rediploidization and adaptive evolution following the whole-genome triplication

Whole-genome duplication (WGD), or polyploidy, events are widespread and significant in the evolutionary history of angiosperms. However, empirical evidence for rediploidization, the major process where polyploids give rise to diploid descendants, is still lacking at the genomic level. Here we present chromosome-scale genomes of the mangrove tree Sonneratia alba and the related inland plant Lagerstroemia speciosa. Their common ancestor has experienced a whole-genome triplication (WGT) approximately 64 million years ago coinciding with a period of dramatic global climate change. Sonneratia, adapting mangrove habitats, experienced extensive chromosome rearrangements post-WGT. We observe the WGT retentions display sequence and expression divergence, suggesting potential neo- and sub-functionalization. Strong selection acting on three-copy retentions indicates adaptive value in response to new environments. To elucidate the role of ploidy changes in genome evolution, we improve a model of the polyploidization-rediploidization process based on genomic evidence, contributing to the understanding of adaptive evolution during climate change.

FLOWERING LOCUS T genes control floral induction in lotus

The transition to flowering is a vital process in the lotus life cycle that significantly impacts its ornamental value and seed production. However, the molecular basis of floral transition in lotus remains largely unknown. Here, eight homologous FLOWERING LOCUS T (FT) genes were initially characterized in lotus, which were designated as NnFT1-NnFT8. All of these genes were found to possess the conserved PEBP domain and exhibited high transcript levels in both lotus leaves and floral organs. The proNnFT:β-glucuronidase (GUS) assay exhibited GUS staining in the vascular tissues of leaves. Furthermore, subcellular localization revealed that NnFT proteins were present in various cellular organelles, including the nucleus, cytoplasm, and endoplasmic reticulum. Overexpression of two NnFT homologs, NnFT2 and NnFT3, rescued the late flowering phenotype in the Arabidopsis ft-10 mutant, indicating the stimulative roles of NnFTs in floral induction. Moreover, NnFTs demonstrated interactions with a bZIP transcription factor, FLOWERING LOCUS D (NnFD), both in vitro and in vivo. These findings will not only deepen our understanding of the regulatory mechanism underlying lotus floral transition, but also provide valuable genetic resources for creating new lotus varieties with extended blooming periods using molecular strategies in the future.

The chromosome-level genome of the submerged plant Cryptocoryne crispatula provides insights into the terrestrial-freshwater transition in Araceae

Plant terrestrialization (i.e. the transition to a terrestrial environment) is a significant evolutionary event that has been intensively studied. While certain plant lineages, particularly in angiosperms, have re-adapted to freshwater habitats after colonizing terrene, however, the molecular mechanism of the terrestrial-freshwater (T-F) transition remains limited. Here, the basal monocot Araceae was selected as the study object to explore the T-F transition adaptation mechanism by comparative genomic analysis. Our findings revealed that the substitution rates significantly increased in the lineage of freshwater Araceae, which may promote their adaptation to the freshwater habitat. Additionally, 20 gene sets across all four freshwater species displayed signs of positive selection contributing to tissue development and defense responses in freshwater plants. Comparative synteny analysis showed that genes specific to submerged plants were enriched in cellular respiration and photosynthesis. In contrast, floating plants were involved in regulating gene expression, suggesting that gene and genome duplications may provide the original material for plants to adapt to the freshwater environment. Our study provides valuable insights into the genomic aspects of the transition from terrestrial to aquatic environments in Araceae, laying the groundwork for future research in the angiosperm.

Exploring domestication pattern in lotus: insights from dispensable genome assembly

Lotus (Nelumbo nucifera Gaertn.), an important aquatic plant in horticulture and ecosystems, has been cultivated for more than 7000 years and domesticated into three different subgroups: flower lotus, rhizome lotus, and seed lotus. To explore the domesticated regions of each subgroup, re-sequencing data of 371 lotus accessions collected from the public database were aligned to the genome of 'China-Antique (CA)'. Unmapped reads were used to build the dispensable genome of each subgroup using a metagenome-like assembly strategy. More than 27 Mb of the dispensable genome in these three subgroups and the wild group was assembled, of which 11,761 genes were annotated. Some of the contigs in the dispensable genome were similar to the genomic segments of other lotus accessions other than 'CA'. The annotated genes in each subgroup played essential roles in specific developmental processes. Dissection of selective signals in three cultivated subgroups also demonstrated that subgroup-specific metabolic pathways, such as the brassinosteroids metabolism enrichment in FL, associated with these selected genes in each subgroup and the contigs in dispensable genome nearly located in the domesticated regions of each subgroup, respectively. Our data presented a valuable resource for facilitating lotus genomic studies, complemented the helpful information to the reference genome, and shed light on the selective signals of domesticated subgroups.

Unlocking secrets of nature's chemists: Potential of CRISPR/Cas-based tools in plant metabolic engineering for customized nutraceutical and medicinal profiles

Plant species have evolved diverse metabolic pathways to effectively respond to internal and external signals throughout their life cycle, allowing adaptation to their sessile and phototropic nature. These pathways selectively activate specific metabolic processes, producing plant secondary metabolites (PSMs) governed by genetic and environmental factors. Humans have utilized PSM-enriched plant sources for millennia in medicine and nutraceuticals. Recent technological advances have significantly contributed to discovering metabolic pathways and related genes involved in the biosynthesis of specific PSM in different food crops and medicinal plants. Consequently, there is a growing demand for plant materials rich in nutrients and bioactive compounds, marketed as "superfoods". To meet the industrial demand for superfoods and therapeutic PSMs, modern methods such as system biology, omics, synthetic biology, and genome editing (GE) play a crucial role in identifying the molecular players, limiting steps, and regulatory circuitry involved in PSM production. Among these methods, clustered regularly interspaced short palindromic repeats-CRISPR associated protein (CRISPR/Cas) is the most widely used system for plant GE due to its simple design, flexibility, precision, and multiplexing capabilities. Utilizing the CRISPR-based toolbox for metabolic engineering (ME) offers an ideal solution for developing plants with tailored preventive (nutraceuticals) and curative (therapeutic) metabolic profiles in an ecofriendly way. This review discusses recent advances in understanding the multifactorial regulation of metabolic pathways, the application of CRISPR-based tools for plant ME, and the potential research areas for enhancing plant metabolic profiles.

Transcription factor NnMYB5 controls petal color by regulating GLUTATHIONE S-TRANSFERASE2 in Nelumbo nucifera

Lotus (Nelumbo spp.) is an important aquatic ornamental genus in the family Nelumbonaceae comprising only 2 species: Nelumbo lutea with yellow flowers and Nelumbo nucifera with red or white flowers. The petal color variations between these 2 species have previously been associated with the potential activities of FLAVONOL SYNTHASE (FLS) and MYB5. However, the underlying genetic mechanisms of flower color divergence within the N. nucifera species remain unclear. Here, quantitative trait locus mapping led to the identification of MYB5, a candidate gene controlling petal color in N. nucifera. Genotyping of 213 natural lotus accessions revealed an 80 kb presence/absence variant (PAV) of the NnMYB5 gene that is associated with petal color variation. Transcriptome analysis, dual-luciferase, and yeast 1-hybrid assays showed that NnMYB5 could directly activate the anthocyanin transporter gene GLUTATHIONE S-TRANSFERASE2 (NnGST2). Heterologous expression of NnGST2 in Arabidopsis (Arabidopsis thaliana) and its overexpression in lotus petals induced anthocyanin accumulation. Deletion of the 80 kb PAV within NnMYB5 inactivated NnGST2 expression and blocked anthocyanin accumulation in white N. nucifera petals. In contrast, the anthocyanin deficiency of N. lutea occurred due to pseudogenized NlMYB5 alleles. Our results establish a regulatory link between NnMYB5 and NnGST2 in petal anthocyanin accumulation and demonstrate the independent mechanisms controlling flower coloration in Nelumbo.

Metabolite profiling and screening of callus browning-related genes in lotus (Nelumbo nucifera)

Callus browning is a major drawback to lotus callus proliferation and regeneration. However, the underlying mechanism of its formation remains largely unknown. Herein, we aimed to explore the metabolic and molecular basis of lotus callus browning by combining histological staining, high-throughput metabolomics, and transcriptomic assays for lotus callus at three browning stages. Histological stained brown callus cross sections displayed severe cell death symptoms, accompanied by an obvious accumulation of polyphenols and lignified materials. Widely targeted metabolomics revealed extensively decreased accumulation of most detected flavonoids and benzylisoquinoline alkaloids (BIAs), as well as a few phenolic acids, amino acids and their derivatives in callus with browning symptoms. Conversely, the contents of most detected tannins were significantly increased. Subsequent comparative transcriptomics identified a set of differentially expressed genes (DEGs) associated with the biosynthesis and regulation of flavonoids and BIAs in lotus. Notably, callus browning was coupled with significantly up-regulated expression of two polyphenol oxidase (PPO) and 17 peroxidase (POD) encoding genes, while the expression of ethylene associated genes remained at marginal levels. These results suggest that lotus callus browning is primarily controlled at the level of metabolism, wherein the oxidation of flavonoids and BIAs is crucially decisive.

Chromosome-level genome assembly of watershield (Brasenia schreberi)

Watershield (Brasenia schreberi) is an aquatic plant that belongs to the basal angiosperm family Cabombaceae. This species has been cultivated as an aquatic vegetable for more than 3000 years in East Asia, but the natural populations have greatly declined in recent decades and have become endangered in several countries of East Asia. In this study, by using PacBio long reads, Illumina short reads, and Hi-C sequencing data, we assembled the genome of B. schreberi, which was approximately 1170.4 Mb in size with a contig N50 of 7.1 Mb. Of the total assembled sequences, 93.6% were anchored to 36 pseudochromosomes with a scaffold N50 of 28.9 Mb. A total of 74,699 protein-coding genes were predicted in the B. schreberi genome, and 558 Mb of repetitive elements occupying 47.69% of the genome were identified. BUSCO analysis yielded a completeness score of 95.8%. The assembled high-quality genome of B. schreberi will be a valuable reference for the study of conservation, evolution and molecular breeding in this species.

6mA DNA Methylation on Genes in Plants Is Associated with Gene Complexity, Expression and Duplication

N6-methyladenine (6mA) DNA methylation has emerged as an important epigenetic modification in eukaryotes. Nevertheless, the evolution of the 6mA methylation of homologous genes after species and after gene duplications remains unclear in plants. To understand the evolution of 6mA methylation, we detected the genome-wide 6mA methylation patterns of four lotus plants (Nelumbo nucifera) from different geographic origins by nanopore sequencing and compared them to patterns in Arabidopsis and rice. Within lotus, the genomic distributions of 6mA sites are different from the widely studied 5mC methylation sites. Consistently, in lotus, Arabidopsis and rice, 6mA sites are enriched around transcriptional start sites, positively correlated with gene expression levels, and preferentially retained in highly and broadly expressed orthologs with longer gene lengths and more exons. Among different duplicate genes, 6mA methylation is significantly more enriched and conserved in whole-genome duplicates than in local duplicates. Overall, our study reveals the convergent patterns of 6mA methylation evolution based on both lineage and duplicate gene divergence, which underpin their potential role in gene regulatory evolution in plants.

Identification, evolution and expression analyses of the whole genome-wide PEBP gene family in Brassica napus L

BACKGROUND

With the release of genomic data for B.rapa, B.oleracea, and B.napus, research on the genetic and molecular functions of Brassica spp. has entered a new stage. PEBP genes in plants play an important role in the transition to flowering as well as seed development and germination. Molecular evolutionary and functional analyses of the PEBP gene family in B.napus based on molecular biology methods can provide a theoretical basis for subsequent investigations of related regulators.

RESULTS

In this paper, we identified a total of 29 PEBP genes from B.napus that were located on 14 chromosomes and 3 random locations. Most members contained 4 exons and 3 introns; motif 1 and motif 2 were the characteristic motifs of PEBP members. On the basis of intraspecific and interspecific collinearity analyses, it is speculated that fragment replication and genomic replication are the main drivers of for the amplification and evolution of the PEBP gene in the B.napus genome. The results of promoter cis-elements prediction suggest that BnPEBP family genes are inducible promoters, which may directly or indirectly participate in multiple regulatory pathways of plant growth cycle. Furthermore, the tissue-specific expression results show that the expression levels of BnPEBP family genes in different tissues were quite different, but the gene expression organization and patterns of the same subgroup were basically the same. qRT‒PCR revealed certain spatiotemporal patterns in the expression of the PEBP subgroups in roots, stems, leaves, buds, and siliques, was tissue-specific, and related to function.

CONCLUSIONS

A systematic comparative analysis of the B.napus PEBP gene family was carried out at here. The results of gene identification, phylogenetic tree construction, structural analysis, gene duplication analysis, prediction of promoter cis-elements and interacting proteins, and expression analysis provide a reference for exploring the molecular mechanisms of BnPEBP family genes in future research.

Comprehensive genome-wide identification and functional characterization of MAPK cascade gene families in Nelumbo

Mitogen-activated protein kinase (MAPK) cascade signaling pathway plays pivotal roles in various plant biological processes. However, systematic study of MAPK cascade gene families is yet to be conducted in lotus. Herein, 198 putative MAPK genes, including 152 MAP3Ks, 15 MKKs, and 31 MPKs genes were identified in Nelumbo. Segmental duplication was identified as the predominant factor driving MAPK cascade gene family expansion in lotus. MAPK cascade genes in N. nucifera and N. lutea shared high degree of sequence homologies, with 84, 9, and 19 homologous MAP3K, MKK, and MPK gene pairs being detected between the two species, respectively, with most genes predominantly undergoing purifying selection. Gene expression profiling indicated that NnMAPK cascade genes were extensively involved in plant development and submergence stress response. Co-expression analysis revealed potential interaction between transcription factors (TFs) and NnMAPK cascade genes in various biological processes. NnMKK showed predicted interactions with multiple NnMAP3K or NnMPK proteins, which suggested that functional diversity of MAPK cascade genes could be as a result of their complex protein interaction mechanisms. This first systematic analysis of MAPK cascade families in lotus provides deeper insights into their evolutionary dynamics and functional properties, which potentially could be crucial for lotus genetic improvement.

Expression inheritance and constraints on cis- and trans-regulatory mutations underlying lotus color variation

Both cis- and trans-regulatory mutations drive changes in gene expression that underpin plant phenotypic evolution. However, how and why these two major types of regulatory mutations arise in different genes and how gene expression is inherited and associated with these regulatory changes are unclear. Here, by studying allele-specific expression in F1 hybrids of pink-flowered sacred lotus (Nelumbo nucifera) and yellow-flowered American lotus (N. lutea), we reveal the relative contributions of cis- and trans-regulatory changes to interspecific expression rewiring underlying petal color change and how the expression is inherited in hybrids. Although cis-only variants influenced slightly more genes, trans-only variants had a stronger impact on expression differences between species. In F1 hybrids, genes under cis-only and trans-only regulatory effects showed a propensity toward additive and dominant inheritance, respectively, whereas transgressive inheritance was observed in genes carrying both cis- and trans-variants acting in opposite directions. By investigating anthocyanin and carotenoid coexpression networks in petals, we found that the same category of regulatory mutations, particularly trans-variants, tend to rewire hub genes in coexpression modules underpinning flower color differentiation between species; we identified 45 known genes with cis- and trans-regulatory variants significantly correlated with flower coloration, such as ANTHOCYANIN 5-AROMATIC ACYLTRANSFERASE (ACT), GLUTATHIONE S-TRANSFERASE F11 (GSTF11), and LYCOPENE Ε-CYCLASE (LCYE). Notably, the relative abundance of genes in different categories of regulatory divergence was associated with the inferred magnitude of constraints like expression level and breadth. Overall, our study suggests distinct selective constraints and modes of gene expression inheritance among different regulatory mutations underlying lotus petal color divergence.

Tracing 100 million years of grass genome evolutionary plasticity

Grasses derive from a family of monocotyledonous plants that includes crops of major economic importance such as wheat, rice, sorghum and barley, sharing a common ancestor some 100 million years ago. The genomic attributes of plant adaptation remain obscure and the consequences of recurrent whole genome duplications (WGD) or polyploidization events, a major force in plant evolution, remain largely speculative. We conducted a comparative analysis of omics data from ten grass species to unveil structural (inversions, fusions, fissions, duplications, substitutions) and regulatory (expression and methylation) basis of genome plasticity, as possible attributes of plant long lasting evolution and adaptation. The present study demonstrates that diverged polyploid lineages sharing a common WGD event often present the same patterns of structural changes and evolutionary dynamics, but these patterns are difficult to generalize across independent WGD events as a result of non-WGD factors such as selection and domestication of crops. Polyploidy is unequivocally linked to the evolutionary success of grasses during the past 100 million years, although it remains difficult to attribute this success to particular genomic consequences of polyploidization, suggesting that polyploids harness the potential of genome duplication, at least partially, in lineage-specific ways. Overall, the present study clearly demonstrates that post-polyploidization reprogramming is more complex than traditionally reported in investigating single species and calls for a critical and comprehensive comparison across independently polyploidized lineages.

Genome structure-based Juglandaceae phylogenies contradict alignment-based phylogenies and substitution rates vary with DNA repair genes

In lineages of allopolyploid origin, sets of homoeologous chromosomes may coexist that differ in gene content and syntenic structure. Presence or absence of genes and microsynteny along chromosomal blocks can serve to differentiate subgenomes and to infer phylogenies. We here apply genome-structural data to infer relationships in an ancient allopolyploid lineage, the walnut family (Juglandaceae), by using seven chromosome-level genomes, two of them newly assembled. Microsynteny and gene-content analyses yield identical topologies that place Platycarya with Engelhardia as did a 1980s morphological-cladistic study. DNA-alignment-based topologies here and in numerous earlier studies instead group Platycarya with Carya and Juglans, perhaps misled by past hybridization. All available data support a hybrid origin of Juglandaceae from extinct or unsampled progenitors nested within, or sister to, Myricaceae. Rhoiptelea chiliantha, sister to all other Juglandaceae, contains proportionally more DNA repair genes and appears to evolve at a rate 2.6- to 3.5-times slower than the remaining species.

MicroRNAs play regulatory roles in genomic balance

Classic genetics studies found that genomic imbalance caused by changing the dosage of part of the genome (aneuploidy) has more detrimental effects than altering the dosage of the whole genome (ploidy). Previous analysis revealed global modulation of gene expression triggered by aneuploidy across various species, including maize (Zea mays), Arabidopsis, yeast, mammals, etc. Plant microRNAs (miRNAs) are a class of 20- to 24-nt endogenous small noncoding RNAs that carry out post-transcriptional gene expression regulation. That miRNAs and their putative targets are preferentially retained as duplicates after whole-genome duplication, as are many transcription factors and signaling components, indicates miRNAs are likely to be dosage-sensitive and potentially involved in genomic balance networks. This review addresses the following questions regarding the role of miRNAs in genomic imbalance. (1) How do aneuploidy and polyploidy impact the expression of miRNAs? (2) Do miRNAs play a regulatory role in modulating the expression of their targets under genomic imbalance?

Evolution of the CBL and CIPK gene families in Medicago: genome-wide characterization, pervasive duplication, and expression pattern under salt and drought stress

BACKGROUND

Calcineurin B-like proteins (CBLs) are ubiquitous Ca²⁺ sensors that mediate plant responses to various stress and developmental processes by interacting with CBL-interacting protein kinases (CIPKs). CBLs and CIPKs play essential roles in acclimatization of crop plants. However, evolution of these two gene families in the genus Medicago is poorly understood.

RESULTS

A total of 68 CBL and 135 CIPK genes have been identified in five genomes from Medicago. Among these genomes, the gene number of CBLs and CIPKs shows no significant difference at the haploid genome level. Phylogenetic and comprehensive characteristic analyses reveal that CBLs and CIPKs are classified into four clades respectively, which is validated by distribution of conserved motifs. The synteny analysis indicates that the whole genome duplication events (WGDs) have contributed to the expansion of both families. Expression analysis demonstrates that two MsCBLs and three MsCIPKs are specifically expressed in roots, mature leaves, developing flowers and nitrogen fixing nodules of Medicago sativa spp. sativa, the widely grown tetraploid species. In particular, the expression of these five genes was highly up-regulated in roots when exposed to salt and drought stress, indicating crucial roles in stress responses.

CONCLUSIONS

Our study leads to a comprehensive understanding of evolution of CBL and CIPK gene families in Medicago, but also provides a rich resource to further address the functions of CBL-CIPK complexes in cultivated species and their closely related wild relatives.

Evolution of isoform-level gene expression patterns across tissues during lotus species divergence

Both gene duplication and alternative splicing (AS) drive the functional diversity of gene products in plants, yet the relative contributions of the two key mechanisms to the evolution of gene function are largely unclear. Here, we studied AS in two closely related lotus plants, Nelumbo lutea and Nelumbo nucifera, and the outgroup Arabidopsis thaliana, for both single-copy and duplicated genes. We show that most splicing events evolved rapidly between orthologs and that the origin of lineage-specific splice variants or isoforms contributed to gene functional changes during species divergence within Nelumbo. Single-copy genes contain more isoforms, have more AS events conserved across species, and show more complex tissue-dependent expression patterns than their duplicated counterparts. This suggests that expression divergence through isoforms is a mechanism to extend the expression breadth of genes with low copy numbers. As compared to isoforms of local, small-scale duplicates, isoforms of whole-genome duplicates are less conserved and display a less conserved tissue bias, pointing towards their contribution to subfunctionalization. Through comparative analysis of isoform expression networks, we identified orthologous genes of which the expression of at least some of their isoforms displays a conserved tissue bias across species, indicating a strong selection pressure for maintaining a stable expression pattern of these isoforms. Overall, our study shows that both AS and gene duplication contributed to the diversity of gene function during the evolution of lotus.

Distinct composition and amplification dynamics of transposable elements in sacred lotus (Nelumbo nucifera Gaertn.)

Sacred lotus (Nelumbo nucifera Gaertn.) is a basal eudicot plant with a unique lifestyle, physiological features, and evolutionary characteristics. Here we report the unique profile of transposable elements (TEs) in the genome, using a manually curated repeat library. TEs account for 59% of the genome, and hAT (Ac/Ds) elements alone represent 8%, more than in any other known plant genome. About 18% of the lotus genome is comprised of Copia LTR retrotransposons, and over 25% of them are associated with non-canonical termini (non-TGCA). Such high abundance of non-canonical LTR retrotransposons has not been reported for any other organism. TEs are very abundant in genic regions, with retrotransposons enriched in introns and DNA transposons primarily in flanking regions of genes. The recent insertion of TEs in introns has led to significant intron size expansion, with a total of 200 Mb in the 28 455 genes. This is accompanied by declining TE activity in intergenic regions, suggesting distinct control efficacy of TE amplification in different genomic compartments. Despite the prevalence of TEs in genic regions, some genes are associated with fewer TEs, such as those involved in fruit ripening and stress responses. Other genes are enriched with TEs, and genes in epigenetic pathways are the most associated with TEs in introns, indicating a dynamic interaction between TEs and the host surveillance machinery. The dramatic differential abundance of TEs with genes involved in different biological processes as well as the variation of target preference of different TEs suggests the composition and activity of TEs influence the path of evolution.

Natural products of medicinal plants: biosynthesis and bioengineering in post-genomic era

Globally, medicinal plant natural products (PNPs) are a major source of substances used in traditional and modern medicine. As we human race face the tremendous public health challenge posed by emerging infectious diseases, antibiotic resistance and surging drug prices etc., harnessing the healing power of medicinal plants gifted from mother nature is more urgent than ever in helping us survive future challenge in a sustainable way. PNP research efforts in the pre-genomic era focus on discovering bioactive molecules with pharmaceutical activities, and identifying individual genes responsible for biosynthesis. Critically, systemic biological, multi- and inter-disciplinary approaches integrating and interrogating all accessible data from genomics, metabolomics, structural biology, and chemical informatics are necessary to accelerate the full characterization of biosynthetic and regulatory circuitry for producing PNPs in medicinal plants. In this review, we attempt to provide a brief update on the current research of PNPs in medicinal plants by focusing on how different state-of-the-art biotechnologies facilitate their discovery, the molecular basis of their biosynthesis, as well as synthetic biology. Finally, we humbly provide a foresight of the research trend for understanding the biology of medicinal plants in the coming decades.

Genome-wide association study of traits in sacred lotus uncovers MITE-associated variants underlying stamen petaloid and petal number variations

Understanding the genetic variants responsible for floral trait diversity is important for the molecular breeding of ornamental flowers. Widely used in water gardening for thousands of years, the sacred lotus exhibits a wide range of diversity in floral organs. Nevertheless, the genetic variations underlying various morphological characteristics in lotus remain largely unclear. Here, we performed a genome-wide association study of sacred lotus for 12 well-recorded ornamental traits. Given a moderate linkage disequilibrium level of 32.9 kb, we successfully identified 149 candidate genes responsible for seven flower traits and plant size variations, including many pleiotropic genes affecting multiple floral-organ-related traits, such as NnKUP2. Notably, we found a 2.75-kb presence-and-absence genomic fragment significantly associated with stamen petaloid and petal number variations, which was further confirmed by re-examining another independent population dataset with petal number records. Intriguingly, this fragment carries MITE transposons bound by siRNAs and is related to the expression differentiation of a nearby candidate gene between few-petalled and double-petalled lotuses. Overall, these genetic variations and candidate genes responsible for diverse lotus traits revealed by our GWAS highlight the role of transposon variations, particularly MITEs, in shaping floral trait diversity.

Identification of QTLs and a putative candidate gene involved in rhizome enlargement of Asian lotus (Nelumbo nucifera)

QTL mapping studies identified three reliable QTLs of rhizome enlargement in lotus. NnBEL6 located within the confidence interval of the major QTL cqREI-LG2 is a key candidate gene enhancing rhizome enlargement. Lotus (Nelumbo) is perennial aquatic plant with nutritional, pharmacological, and ornamental significance. Rhizome is an underground lotus stem that acts as a storage organ and as a reproductive tissue for asexual production. The enlargement of lotus rhizome is an important adaptive strategy for surviving the cold winter. The aims of this study were to identify quantitative trait loci (QTLs) for rhizome enlargement traits including rhizome enlargement index (REI) and number of enlarged rhizome (NER), and to uncover their associated candidate genes. A high-density genetic linkage map was constructed, consisting of 2935 markers binned from 236,840 SNPs. A total of 14 significant QTLs were detected for REI and NER, which explained 6.7-22.3% of trait variance. Three QTL regions were repeatedly identified in at least 2 years, and a major QTL, designated cqREI-LG2, with a rhizome-enlargement effect and about 20% of the phenotypic contribution was identified across the 3 climatic years. A candidate NnBEL6 gene located within the confidence interval of cqREI-LG2 was considered to be putatively involved in lotus rhizome enlargement. The expression of NnBEL6 was exclusively induced by rhizome swelling. Sequence comparison of NnBEL6 among lotus cultivars revealed a functional Indel site in its promoter that likely initiates the rhizome enlargement process. Transgenic potato assay was used to confirm the role of NnBEL6 in inducing tuberization. The successful identification QTLs and functional validation of NnBEL6 gene reported in this study will enrich our knowledge on the genetic basis of rhizome enlargement in lotus.

Diversity and determinants of recombination landscapes in flowering plants

During meiosis, crossover rates are not randomly distributed along the chromosome and their location may have a strong impact on the functioning and evolution of the genome. To date, the broad diversity of recombination landscapes among plants has rarely been investigated and a formal comparative genomic approach is still needed to characterize and assess the determinants of recombination landscapes among species and chromosomes. We gathered genetic maps and genomes for 57 flowering plant species, corresponding to 665 chromosomes, for which we estimated large-scale recombination landscapes. We found that the number of crossover per chromosome spans a limited range (between one to five/six) whatever the genome size, and that there is no single relationship across species between genetic map length and chromosome size. Instead, we found a general relationship between the relative size of chromosomes and recombination rate, while the absolute length constrains the basal recombination rate for each species. At the chromosome level, we identified two main patterns (with a few exceptions) and we proposed a conceptual model explaining the broad-scale distribution of crossovers where both telomeres and centromeres play a role. These patterns correspond globally to the underlying gene distribution, which affects how efficiently genes are shuffled at meiosis. These results raised new questions not only on the evolution of recombination rates but also on their distribution along chromosomes.

Meiotic recombination is a universal feature of sexually reproducing species. During meiosis, crossovers play a fundamental role for the proper segregation of chromosomes during meiosis and reshuffles alleles among chromosomes. How much variation in recombination is expected within a genome and among different species remains a central question for understanding the evolution of recombination. We characterized and compared recombination landscapes in a large set of plant species with a wide range of genome size. We found that the number of crossovers varied little among species, from one mandatory to no more than five or six crossovers per chromosomes, whatever the genome size. However, we identified two main patterns of variation along chromosomes (with a few exceptions) that can be explained by a new conceptual model where chromosome length, chromosome structure and gene density play a role. The strong association between gene density and recombination was already known, but raised new questions not only about the evolution of recombination rates but also on their distribution along chromosomes.

Characterization of Transposon-Derived Accessible Chromatin Regions in Rice (Oryza Sativa)

Growing evidence indicates that transposons or transposable elements (TEs)-derived accessible chromatin regions (ACRs) play essential roles in multiple biological processes by interacting with trans-acting factors. However, the function of TE-derived ACRs in the regulation of gene expression in the rice genome has not been well characterized. In this study, we examined the chromatin dynamics in six types of rice tissues and found that ~8% of ACRs were derived from TEs and exhibited distinct levels of accessibility and conservation as compared to those without TEs. TEs exhibited a TE subtype-dependent impact on ACR formation, which can be mediated by changes in the underlying DNA methylation levels. Moreover, we found that tissue-specific TE-derived ACRs might function in the tissue development through the modulation of nearby gene expression. Interestingly, many genes in domestication sweeps were found to overlap with TE-derived ACRs, suggesting their potential functions in the rice domestication. In addition, we found that the expression divergence of 1070 duplicate gene pairs were associated with TE-derived ACRs and had distinct distributions of TEs and ACRs around the transcription start sites (TSSs), which may experience different selection pressures. Thus, our study provides some insights into the biological implications of TE-derived ACRs in the rice genome. Our results imply that these ACRs are likely involved in the regulation of tissue development, rice domestication and functional divergence of duplicated genes.

The multiple fates of gene duplications: Deletion, hypofunctionalization, subfunctionalization, neofunctionalization, dosage balance constraints, and neutral variation

Gene duplications have long been recognized as a contributor to the evolution of genes with new functions. Multiple copies of genes can result from tandem duplication, from transposition to new chromosomes, or from whole-genome duplication (polyploidy). The most common fate is that one member of the pair is deleted to return the gene to the singleton state. Other paths involve the reduced expression of both copies (hypofunctionalization) that are held in duplicate to maintain sufficient quantity of function. The two copies can split functions (subfunctionalization) or can diverge to generate a new function (neofunctionalization). Retention of duplicates resulting from doubling of the whole genome occurs for genes involved with multicomponent interactions such as transcription factors and signal transduction components. In contrast, these classes of genes are underrepresented in small segmental duplications. This complementary pattern suggests that the balance of interactors affects the fate of the duplicate pair. We discuss the different mechanisms that maintain duplicated genes, which may change over time and intersect.

The slow-evolving Acorus tatarinowii genome sheds light on ancestral monocot evolution

Monocots are one of the most diverse groups of flowering plants, and tracing the evolution of their ancestral genome into modern species is essential for understanding their evolutionary success. Here, we report a high-quality assembly of the Acorus tatarinowii genome, a species that diverged early from all the other monocots. Genome-wide comparisons with a range of representative monocots characterized Acorus as a slowly evolved genome with one whole-genome duplication. Our inference of the ancestral monocot karyotypes provides new insights into the chromosomal evolutionary history assigned to modern species and reveals the probable molecular functions and processes related to the early adaptation of monocots to wetland or aquatic habitats (that is, low levels of inorganic phosphate, parallel leaf venation and ephemeral primary roots). The evolution of ancestral gene order in monocots is constrained by gene structural and functional features. The newly obtained Acorus genome offers crucial evidence for delineating the origin and diversification of monocots, including grasses.

Studies on Lotus Genomics and the Contribution to Its Breeding

Lotus (Nelumbo nucifera), under the Nelumbonaceae family, is one of the relict plants possessing important scientific research and economic values. Because of this, much attention has been paid to this species on both its biology and breeding among the scientific community. In the last decade, the genome of lotus has been sequenced, and several high-quality genome assemblies are available, which have significantly facilitated functional genomics studies in lotus. Meanwhile, re-sequencing of the natural and genetic populations along with different levels of omics studies have not only helped to classify the germplasm resources but also to identify the domestication of selected regions and genes controlling different horticultural traits. This review summarizes the latest progress of all these studies on lotus and discusses their potential application in lotus breeding.

Genome-Wide Characterization and Comprehensive Analysis of NAC Transcription Factor Family in Nelumbo nucifera

NAC (NAM, ATAF, and CUC) is a ubiquitously expressed plant-specific transcription factor (TF) family which is involved in the regulation of various biological processes. However, a systematic characterization of NAC gene family is yet to be reported in lotus. Here, 82 NnNAC genes which included five predicted membrane-bound NAC proteins were identified in the lotus genome. Phylogenetic analysis revealed seven-subfamily clusters (I–VII) of NnNAC proteins, with homologous gene pairs displaying similar conserved motifs and gene structure characteristics. Transactivation assay of NnNAC proteins revealed an extensive transcriptional activation capacity which is mediated by the highly divergent C-terminal activation domain (AD). Expression analysis of NnNAC genes in lotus tissues showed high transcript levels in root, stamen, petal and seed coat. In addition, 30 and 29 differentially expressed NnNAC candidate genes putatively involved in lotus seed development and response to complete submergence stress, respectively, were identified. Overall, our study provides potentially useful candidate gene resources for future molecular breeding of lotus varieties with novel agronomic traits.

BLSSpeller to discover novel regulatory motifs in maize

With the decreasing cost of sequencing and availability of larger numbers of sequenced genomes, comparative genomics is becoming increasingly attractive to complement experimental techniques for the task of transcription factor (TF) binding site identification. In this study, we redesigned BLSSpeller, a motif discovery algorithm, to cope with larger sequence datasets. BLSSpeller was used to identify novel motifs in Zea mays in a comparative genomics setting with 16 monocot lineages. We discovered 61 motifs of which 20 matched previously described motif models in Arabidopsis. In addition, novel, yet uncharacterized motifs were detected, several of which are supported by available sequence-based and/or functional data. Instances of the predicted motifs were enriched around transcription start sites and contained signatures of selection. Moreover, the enrichment of the predicted motif instances in open chromatin and TF binding sites indicates their functionality, supported by the fact that genes carrying instances of these motifs were often found to be co-expressed and/or enriched in similar GO functions. Overall, our study unveiled several novel candidate motifs that might help our understanding of the genotype to phenotype association in crops.

Comprehensive Analysis for GRF Transcription Factors in Sacred Lotus (Nelumbo nucifera)

Sacred lotus (Nelumbo nucifera) is an aquatic perennial plant with essential food, ornamental, and pharmacological value. Growth-regulating factor (GRF) is a transcription factor (TF) family that plays an important role in regulating the growth and development of plants. In this study, a comprehensive analysis of the GRF family in N. nucifera was performed, and its role in N. nucifera development was studied. A total of eight GRF genes were identified in the N. nucifera genome. Phylogenetic analysis divided the 38 GRF genes into six clades, while the NuGRFs only contained five clades. The analyses of gene structures, motifs, and cis-acting regulatory elements of the GRF gene family were performed. In addition, the chromosome location and collinearity were analyzed. The expression pattern based on transcriptomic data and real-time reverse transcription-quantitative PCR (qRT-PCR) revealed that the GRF genes were expressed in multiple organs and were abundant in actively growing tissues, and the expression levels decreased as the age of N. nucifera increased. Then, 3D structures of the NuGRF proteins were predicted by homology modeling. Finally, the subcellular localization of GRF1 was ascertained in the tobacco leaf through a vector. Therefore, this study provides a comprehensive overview of the GRF TF family in N. nucifera.

Comparative analyses of American and Asian lotus genomes reveal insights into petal color, carpel thermogenesis and domestication

Nelumbo lutea (American lotus), which differs from Nelumbo nucifera (Asian lotus) morphologically, is one of the two remaining species in the basal eudicot family Nelumbonaceae. Here, we assembled the 843-Mb genome of American lotus into eight pseudochromosomes containing 31 382 protein-coding genes. Comparative analyses revealed conserved synteny without large chromosomal rearrangements between the genomes of American and Asian lotus and identified 29 533 structural variants (SVs). Carotenoid and anthocyanin pigments determine the yellow and red petal colors of American and Asian lotus, respectively. The structural genes encoding enzymes of the carotenoid and anthocyanin biosynthesis pathways were conserved between two species but differed in expression. We detected SVs caused by repetitive sequence expansion or contraction among the anthocyanin biosynthesis regulatory MYB genes. Further transient overexpression of candidate NnMYB5 induced anthocyanin accumulation in lotus petals. Alternative oxidase (AOX), uncoupling proteins (UCPs), and sugar metabolism and transportation contributed to carpel thermogenesis. Carpels produce heat with sugars transported from leaves as the main substrates, because there was weak tonoplast sugar transporter (TST) activity, and with SWEETs were highly expressed during thermogenesis. Cell proliferation-related activities were particularly enhanced in the warmer carpels compared with stamens during the cold night before blooming, which suggested that thermogenesis plays an important role in flower protogyny. Population genomic analyses revealed deep divergence between American and Asian lotus, and independent domestication affecting seed, rhizome, and flower traits. Our findings provide a high-quality reference genome of American lotus for exploring the genetic divergence and variation between two species and revealed possible genomic bases for petal color, carpel thermogenesis and domestication in lotus.

Expression rewiring and methylation of non-coding RNAs involved in rhizome phenotypic variations of lotus ecotypes

Non-coding RNAs (ncRNAs), including miRNAs, lncRNAs, and circRNAs, emerge as crucial components for gene regulation. Nelumbo nucifera (lotus), a horticulturally important plant, differentiates into a temperate ecotype of enlarged rhizomes and a tropical ecotype of thin rhizomes. Nevertheless, whether and how ncRNAs can be rewired in expression and differentially methylated contributing to adaptive divergence of this storage organ in lotus ecotypes is unclear. Herein, we study the expression behaviors and DNA methylation patterns of ncRNAs in temperate and tropical lotus rhizomes. By whole transcriptome sequencing, we found both mRNAs and lncRNAs have divergent expression patterns between ecotypes, whereas miRNAs and circRNAs tended to be accession-specific or noisier in expression. The differentially expressed ncRNAs are involved in phenotypic differentiation of lotus rhizome between ecotypes, as the genes that interacted with them in the competing endogenous RNA network are enriched in functions including carbohydrate metabolism and plant hormone signaling, being critical to rhizome enlargement. Intriguingly, ncRNA-targeted genes are less prone to show positive selection or differential expression during ecotypic divergence due to constraints from ncRNA-mRNA interactions. The methylation levels of ncRNAs generally tend to be higher in temperate lotus than in tropical lotus, and differential methylation of lncRNAs also tends to have expression changes. Overall, our study of ncRNAs and their targets highlights the role of ncRNAs in rhizome growth variation between lotus ecotypes through expression rewiring and methylation modification.

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.

The chromosome-level genome of a free-floating aquatic weed Pistia stratiotes provides insights into its rapid invasion

Pistia stratiotes (Araceae), commonly referred to as water lettuce, is one of the most notorious weeds that cause severe damage to the economy and natural ecosystems of infested areas. In order to explore the mechanism of its rapid invasion, here, we assembled a high-quality chromosome-level genome for P. stratiotes based on the Illumina sequencing, PacBio sequencing, and Hi-C scaffolding technology. The assembled genome is 311.87 Mb in size with a contig N50 of 1.08 Mb. The contigs were further anchored on 14 pseudochromosomes with a scaffold N50 of 21.21 Mb. A total of 20,356 protein-coding genes were predicted, of which 79.35% were functionally annotated here. Evolutionary analysis showed that P. stratiotes and Colocasia esculenta were clustered together as sister lineages that diverged approximately 61 Mya. The synteny analyses indicated that two whole-genome duplication (WGD) events occurred within a short period in P. stratiotes. Moreover, comparative genome analysis indicated that the expansion of gene families corresponding to disease resistance might contribute to rapid invasion in P. stratiotes. Also, we analyzed the disease-resistance gene family (NBS-LRR) involved in plant defense. A genome-wide search in P. stratiotes genome identified 85 NBS-LRR genes in this study. In conclusion, our present study provides some new insights into the evolution of the invasive aquatic plant P. stratiotes. Our reference genome will also provide valuable resources for future invasion genomics research programs.

Transcriptome-Wide Characterization of Alkaloids and Chlorophyll Biosynthesis in Lotus Plumule

Lotus plumule is a green tissue in the middle of seeds that predominantly accumulates bisbenzylisoquinoline alkaloids (bis-BIAs) and chlorophyll (Chl). However, the biosynthetic mechanisms of these two metabolites remain largely unknown in lotus. This study used physiological and RNA sequencing (RNA-Seq) approaches to characterize the development and molecular mechanisms of bis-BIAs and Chl biosynthesis in lotus plumule. Physiological analysis revealed that exponential plumule growth occurred between 9 and 15 days after pollination (DAP), which coincided with the onset of bis-BIAs biosynthesis and its subsequent rapid accumulation. Transcriptome analysis of lotus plumule identified a total of 8,725 differentially expressed genes (DEGs), representing ~27.7% of all transcripts in the lotus genome. Sixteen structural DEGs, potentially associated with bis-BIAs biosynthesis, were identified. Of these, 12 encoded O-methyltransferases (OMTs) are likely involved in the methylation and bis-BIAs diversity in lotus. In addition, functionally divergent paralogous and redundant homologous gene members of the BIAs biosynthesis pathway, as well as transcription factors co-expressed with bis-BIAs and Chl biosynthesis genes, were identified. Twenty-two genes encoding 16 conserved enzymes of the Chl biosynthesis pathway were identified, with the majority being significantly upregulated by Chl biosynthesis. Photosynthesis and Chl biosynthesis pathways were simultaneously activated during lotus plumule development. Moreover, our results showed that light-driven Pchlide reduction is essential for Chl biosynthesis in the lotus plumule. These results will be useful for enhancing our understanding of alkaloids and Chl biosynthesis in plants.

Chromosome evolution and the genetic basis of agronomically important traits in greater yam

The nutrient-rich tubers of the greater yam, Dioscorea alata L., provide food and income security for millions of people around the world. Despite its global importance, however, greater yam remains an orphan crop. Here, we address this resource gap by presenting a highly contiguous chromosome-scale genome assembly of D. alata combined with a dense genetic map derived from African breeding populations. The genome sequence reveals an ancient allotetraploidization in the Dioscorea lineage, followed by extensive genome-wide reorganization. Using the genomic tools, we find quantitative trait loci for resistance to anthracnose, a damaging fungal pathogen of yam, and several tuber quality traits. Genomic analysis of breeding lines reveals both extensive inbreeding as well as regions of extensive heterozygosity that may represent interspecific introgression during domestication. These tools and insights will enable yam breeders to unlock the potential of this staple crop and take full advantage of its adaptability to varied environments.

Genomic variation reveals demographic history and biological adaptation of the ancient relictual, lotus (Nelumbo Adans)

Lotus (Nelumbo Adans.), a relict plant, is the testimony of long-term sustained ecological success, but the underlying genetic changes related to its survival strategy remains unclear. Here, we assembled the high-quality lotus genome, investigated genome variation of lotus mutation accumulation (MA) lines and reconstructed the demographic history of wild Asian lotus, respectively. We identified and validated 43 base substitutions fixed in MA lines, implying a spontaneous mutation rate of 1.4 × 10-9 base/generation in lotus shoot stem cells. The past history of lotus revealed that the ancestors of lotus in eastern and southern Asia could be traced back ~20 million years ago (Mya) and experienced twice significant bottlenecks and population splits. We further identified the selected genes among three lotus groups in different habitats, suggesting that 453 genes between tropical and temperate group and 410 genes between two subgroups from Northeastern China and the Yangtze River - Yellow River Basin might play important roles in natural selection in lotus's adaptation and resilience. Our findings not only improve an understanding of the lotus evolutionary history and the genetic basis of its survival advantages, but also provide valuable data for addressing various questions in evolution and protection for the relict plants.

Chromosome-level genome assembly of Aristolochia contorta provides insights into the biosynthesis of benzylisoquinoline alkaloids and aristolochic acids

Aristolochic acids (AAs) and their derivatives exist in multiple Aristolochiaceae species which had been or are being used as medicinal materials. During the past decades, AAs have received increasing attention due to their nephrotoxicity and carcinogenecity. Elimination of AAs in medicinal materials using biotechnological approaches is important to improve medication safety. However, it has not been achieved because of the limited information of AA biosynthesis available. Here, we report a high-quality reference-grade genome assembly of the AA-containing vine, Aristolochia contorta. Total size of the assembly is 209.27 Mb, which is assembled into 7 pseudochromosomes. Synteny analysis, Ks distribution and 4DTv suggest absences of whole-genome duplication events in A. contorta after the angiosperm-wide WGD. Based on genomic, transcriptomic and metabolic data, pathways and candidate genes of benzylisoquinoline alkaloid (BIA) and AA biosynthesis in A. contorta were proposed. Five O-methyltransferase genes, including AcOMT1-3, AcOMT5 and AcOMT7, were cloned and functionally characterized. The results provide a high-quality reference genome for AA-containing species of Aristolochiaceae. It lays a solid foundation for further elucidation of AA biosynthesis and regulation and molecular breeding of Aristolochiaceae medicinal materials.

Spatial and Temporal Transcriptomic Heredity and Asymmetry in an Artificially Constructed Allotetraploid Wheat (AADD)

Polyploidy, or whole-genome duplication (WGD), often induces dramatic changes in gene expression due to "transcriptome shock. " However, questions remain about how allopolyploidy (the merging of multiple nuclear genomes in the same nucleus) affects gene expression within and across multiple tissues and developmental stages during the initial foundation of allopolyploid plants. Here, we systematically investigated the immediate effect of allopolyploidy on gene expression variation in an artificial allopolyploidy system consisting of a constructed allotetraploid wheat (AADD genome, accession AT2) and its diploid progenitors Triticum urartu and Aegilops tauschii. We performed comprehensive RNA sequencing of 81 samples from different genotypes, tissues, and developmental stages. First, we found that intrinsic interspecific differences between the diploid parents played a major role in establishing the expression architecture of the allopolyploid. Nonetheless, allopolyploidy per se also induced dramatic and asymmetric patterns of differential gene expression between the subgenomes, and genes from the D subgenome exhibited a more drastic response. Second, analysis of homoeolog expression bias (HEB) revealed that the D subgenome exhibited significant expression bias and that de novo-generated HEB was attributed mainly to asymmetrical differential gene expression. Homoeolog-specific expression (HSE) analyses showed that the cis-only regulatory pattern was predominant in AT2, reflecting significant divergence between the parents. Co-expression network analysis revealed that homoeolog expression connectivity (HEC) was significantly correlated with sequence divergence in cis elements between subgenomes. Interestingly, allopolyploidy-induced reconstruction of network modules was also associated with different HSE patterns. Finally, a transcriptome atlas of spike development demonstrated that the phenotypic similarity of AT2 to T. urartu may be attributed to the combination of relatively stable expression of A-subgenome genes and drastic downregulation of their D-subgenome homoeologs. These findings provide a broad, multidimensional characterization of allopolyploidy-induced transcriptomic responses and suggest that allopolyploidy can have immediate and complex regulatory effects on the expression of nuclear genes.

Evolution of RGF/GLV/CLEL Peptide Hormones and Their Roles in Land Plant Growth and Regulation

Rooting is a key innovation during plant terrestrialization. RGFs/GLVs/CLELs are a family of secreted peptides, playing key roles in root stem cell niche maintenance and pattern formation. The origin of this peptide family is not well characterized. RGFs and their receptor genes, RGIs, were investigated comprehensively using phylogenetic and genetic analyses. We identified 203 RGF genes from 24 plant species, representing a variety of land plant lineages. We found that the RGF genes originate from land plants and expand via multiple duplication events. The lineage-specific RGF duplicates are retained due to their regulatory divergence, while a majority of RGFs experienced strong purifying selection in most land plants. Functional analysis indicated that RGFs and their receptor genes, RGIs, isolated from liverwort, tomato, and maize possess similar biological functions with their counterparts from Arabidopsis in root development. RGFs and RGIs are likely coevolved in land plants. Our studies shed light on the origin and functional conservation of this important peptide family in plant root development.

One Hundred Years of Gene Balance: How Stoichiometric Issues Affect Gene Expression, Genome Evolution, and Quantitative Traits

A century ago experiments with the flowering plant Datura stramonium and the fruit fly Drosophila melanogaster revealed that adding an extra chromosome to a karyotype was much more detrimental than adding a whole set of chromosomes. This phenomenon was referred to as gene balance and has been recapitulated across eukaryotic species. Here, we retrace some developments in this field. Molecular studies suggest that the basis of balance involves stoichiometric relationships of multi-component interactions. This concept has implication for the mechanisms controlling gene expression, genome evolution, sex chromosome evolution/dosage compensation, speciation mechanisms, and the underlying genetics of quantitative traits.

Multi-tissue transcriptome analysis of two Begonia species reveals dynamic patterns of evolution in the chalcone synthase gene family

Begonia is an important horticultural plant group, as well as one of the most speciose Angiosperm genera, with over 2000 described species. Genus wide studies of genome size have shown that Begonia has a highly variable genome size, and analysis of paralog pairs has previously suggested that Begonia underwent a whole genome duplication. We address the contribution of gene duplication to the generation of diversity in Begonia using a multi-tissue RNA-seq approach. We chose to focus on chalcone synthase (CHS), a gene family having been shown to be involved in biotic and abiotic stress responses in other plant species, in particular its importance in maximising the use of variable light levels in tropical plants. We used RNA-seq to sample six tissues across two closely related but ecologically and morphologically divergent species, Begonia conchifolia and B. plebeja, yielding 17,012 and 19,969 annotated unigenes respectively. We identified the chalcone synthase gene family members in our Begonia study species, as well as in Hillebrandia sandwicensis, the monotypic sister genus to Begonia, Cucumis sativus, Arabidopsis thaliana, and Zea mays. Phylogenetic analysis suggested the CHS gene family has high duplicate turnover, all members of CHS identified in Begonia arising recently, after the divergence of Begonia and Cucumis. Expression profiles were similar within orthologous pairs, but we saw high inter-ortholog expression variation. Sequence analysis showed relaxed selective constraints on some ortholog pairs, with substitutions at conserved sites. Evidence of pseudogenisation and species specific duplication indicate that lineage specific differences are already beginning to accumulate since the divergence of our study species. We conclude that there is evidence for a role of gene duplication in generating diversity through sequence and expression divergence in Begonia.

Genomic insights into molecular adaptation to intertidal environments in the mangrove Aegiceras corniculatum

Mangroves have colonised extreme intertidal environments characterised by high salinity, hypoxia and other abiotic stresses. Aegiceras corniculatum, a pioneer mangrove species that has evolved two specialised adaptive traits (salt secretion and crypto-vivipary) is an attractive ecological model to investigate molecular mechanisms underlying adaptation to intertidal environments. We assembled de novo a high-quality reference genome of A. corniculatum and performed comparative genomic and transcriptomic analyses to investigate molecular mechanisms underlying adaptation to intertidal environments. We provide evidence that A. corniculatum experienced a whole-genome duplication (WGD) event c. 35 Ma. We infer that maintenance of cellular environmental homeostasis is an important adaptive process in A. corniculatum. The 14-3-3 and H⁺ -ATPase protein-coding genes, essential for the salt homeostasis, were preferentially retained after the recent WGD event. Using comparative transcriptomics, we show that genes upregulated under high-salt conditions are involved in salt transport and ROS scavenging. We also found that all homologues of DELAY OF GERMINATION1 (DOG1) had lost their heme-binding ability in A. corniculatum, and that this may contribute to crypto-vivipary. Our study provides insight into the genomic correlates of phenotypic adaptation to intertidal environments. This could contribute not only within the genomics community, but also to the field of plant evolution.

Epimutations Define a Fast-Ticking Molecular Clock in Plants

Stochastic gains and losses of DNA methylation at CG dinucleotides are a frequent occurrence in plants. These spontaneous 'epimutations' occur at a rate that is 100 000 times higher than the genetic mutation rate, are effectively neutral at the genome-wide scale, and are stably inherited across mitotic and meiotic cell divisions. Mathematical models have been extraordinarily successful at describing how epimutations accumulate in plant genomes over time, making this process one of the most predictable epigenetic phenomena to date. Here, we propose that their high rate and effective neutrality make epimutations a powerful new molecular clock for timing evolutionary events of the recent past and for age dating of long-lived perennials such as trees.

Patterns and Processes of Diploidization in Land Plants

Most land plants are now known to be ancient polyploids that have rediploidized. Diploidization involves many changes in genome organization that ultimately restore bivalent chromosome pairing and disomic inheritance, and resolve dosage and other issues caused by genome duplication. In this review, we discuss the nature of polyploidy and its impact on chromosome pairing behavior. We also provide an overview of two major and largely independent processes of diploidization: cytological diploidization and genic diploidization/fractionation. Finally, we compare variation in gene fractionation across land plants and highlight the differences in diploidization between plants and animals. Altogether, we demonstrate recent advancements in our understanding of variation in the patterns and processes of diploidization in land plants and provide a road map for future research to unlock the mysteries of diploidization and eukaryotic genome evolution.

Analysis of the Coptis chinensis genome reveals the diversification of protoberberine-type alkaloids

Chinese goldthread (Coptis chinensis Franch.), a member of the Ranunculales, represents an important early-diverging eudicot lineage with diverse medicinal applications. Here, we present a high-quality chromosome-scale genome assembly and annotation of C. chinensis. Phylogenetic and comparative genomic analyses reveal the phylogenetic placement of this species and identify a single round of ancient whole-genome duplication (WGD) shared by the Ranunculaceae. We characterize genes involved in the biosynthesis of protoberberine-type alkaloids in C. chinensis. In particular, local genomic tandem duplications contribute to member amplification of a Ranunculales clade-specific gene family of the cytochrome P450 (CYP) 719. The functional versatility of a key CYP719 gene that encodes the (S)-canadine synthase enzyme involved in the berberine biosynthesis pathway may play critical roles in the diversification of other berberine-related alkaloids in C. chinensis. Our study provides insights into the genomic landscape of early-diverging eudicots and provides a valuable model genome for genetic and applied studies of Ranunculales.