The water lily genome and the early evolution of flowering plants

Genome evolution and diversity of wild and cultivated potatoes

Potato (Solanum tuberosum L.) is the world's most important non-cereal food crop, and the vast majority of commercially grown cultivars are highly heterozygous tetraploids. Advances in diploid hybrid breeding based on true seeds have the potential to revolutionize future potato breeding and production^1-4. So far, relatively few studies have examined the genome evolution and diversity of wild and cultivated landrace potatoes, which limits the application of their diversity in potato breeding. Here we assemble 44 high-quality diploid potato genomes from 24 wild and 20 cultivated accessions that are representative of Solanum section Petota, the tuber-bearing clade, as well as 2 genomes from the neighbouring section, Etuberosum. Extensive discordance of phylogenomic relationships suggests the complexity of potato evolution. We find that the potato genome substantially expanded its repertoire of disease-resistance genes when compared with closely related seed-propagated solanaceous crops, indicative of the effect of tuber-based propagation strategies on the evolution of the potato genome. We discover a transcription factor that determines tuber identity and interacts with the mobile tuberization inductive signal SP6A. We also identify 561,433 high-confidence structural variants and construct a map of large inversions, which provides insights for improving inbred lines and precluding potential linkage drag, as exemplified by a 5.8-Mb inversion that is associated with carotenoid content in tubers. This study will accelerate hybrid potato breeding and enrich our understanding of the evolution and biology of potato as a global staple food crop.

Evolution of coastal forests based on a full set of mangrove genomes

Genomic studies are now poised to explore whole communities of species. The ~70 species of woody plants that anchor the coastal ecosystems of the tropics, collectively referred to as mangroves, are particularly suited to this exploration. In this study, we de novo sequenced the whole genomes of 32 mangroves, which we combined with other sequences of 30 additional species, comprising almost all mangroves globally. These community-wide genomic data will be valuable for ecology, evolution and biodiversity research. While the data revealed 27 independent origins of mangroves, the total phylogeny shows only modest increases in species number, even in coastal areas of active speciation, suggesting that mangrove extinction is common. A possible explanation for common extinction is the frequent sea-level rises and falls (SLRs and SLFs) documented in the geological record. Indeed, near-extinctions of species with extremely small population size (N) often happened during periods of rapid SLR, as revealed by the genome-wide heterozygosity of almost all mangroves. Reduction in N has possibly been further compounded by population fragmentation and the subsequent accumulation of deleterious mutations, thus pushing mangroves even closer to extinction. Crucially, the impact of the next SLR will be exacerbated by human encroachment into these mangrove habitats, potentially altering the ecosystems of tropical coasts irreversibly.

The chromosome-level genome for Toxicodendron vernicifluum provides crucial insights into Anacardiaceae evolution and urushiol biosynthesis

The lacquer tree (Toxicodendron vernicifluum (Stokes) F.A. Barkley) is an important tree with economic, industrial, and medicinal values. Here, we generated the reference genome of T. vernicifluum at the chromosome level with 491.93 Mb in size, in which 98.26% of the assembled contigs were anchored onto 15 pseudochromosomes with the scaffold N50 of 32.97 Mb. Comparative genomic analysis revealed the gene families related to urushiol biosynthesis were expanded, contributing to the ecological fitness and biological adaptability of the lacquer tree. We combined multi-omics data to identify genes that encode key enzymes in the T. vernicifluum urushiol and lignin biosynthetic pathways. Furthermore, the unique active metabolites, such as butin and fisetin, in cultivar lacquers were identified by metabolism profiling. Our work would provide crucial insights into metabolite synthesis such as urushiol and lignin, meanwhile offer a basis for further exploration of the cultivation and breeding of T. vernicifluum and other Anacardiaceae members.

•

We provide the first chromosome-level genome for Toxicodendron vernicifluum

•

We explore the phylogenetic position of lacquer tree

•

We identify genes involved in the urushiol and lignin biosynthetic pathways

•

Analysis of the unique active metabolites in cultivar lacquers by LC-MS/MS is reported

The adaptive evolution of Euryale ferox to the aquatic environment through paleo-hexaploidization

Occupation of living space is one of the main driving forces of adaptive evolution, especially for aquatic plants whose leaves float on the water surface and thus have limited living space. Euryale ferox, from the angiosperm basal family Nymphaeaceae, develops large, rapidly expanding leaves to compete for space on the water surface. Microscopic observation found that the cell proliferation of leaves is almost completed underwater, while the cell expansion occurs rapidly after they grow above water. To explore the mechanism underlying the specific development of leaves, we performed sequences assembly and analyzed the genome and transcriptome dynamics of E. ferox. Through reconstruction of the three sub-genomes generated from the paleo-hexaploidization event in E. ferox, we revealed that one sub-genome was phylogenetically closer to Victoria cruziana, which also exhibits gigantic floating leaves. Further analysis revealed that while all three sub-genomes promoted the evolution of the specific leaf development in E. ferox, the genes from the sub-genome closer to V. cruziana contributed more to this adaptive evolution. Moreover, we found that genes involved in cell proliferation and expansion, photosynthesis, and energy transportation were over-retained and showed strong expression association with the leaf development stages, such as the expression divergence of SWEET orthologs as energy uploaders and unloaders in the sink and source leaf organs of E. ferox. These findings provide novel insights into the genome evolution through polyploidization, as well as the adaptive evolution regarding the leaf development accomplished through biased gene retention and expression sub-functionalization of multi-copy genes in E. ferox.

Genomic Analysis Based on Chromosome-Level Genome Assembly Reveals an Expansion of Terpene Biosynthesis of Azadirachta indica

Azadirachta indica (neem), an evergreen tree of the Meliaceae family, is a source of the potent biopesticide azadirachtin. The lack of a chromosome-level assembly impedes an in-depth understanding of its genome architecture and the comparative genomic analysis of A. indica. Here, a high-quality genome assembly of A. indica was constructed using a combination of data from Illumina, PacBio, and Hi-C technology, which is the first chromosome-scale genome assembly of A. indica. Based on the length of our assembly, the genome size of A. indica is estimated to be 281 Mb anchored to 14 chromosomes (contig N50 = 6 Mb and scaffold N50 = 19 Mb). The genome assembly contained 115 Mb repetitive elements and 25,767 protein-coding genes. Evolutional analysis revealed that A. indica didn't experience any whole-genome duplication (WGD) event after the core eudicot γ event, but some genes and genome segment might likely experienced recent duplications. The secondary metabolite clusters, TPS genes, and CYP genes were also identified. Comparative genomic analysis revealed that most of the A. indica-specific TPS genes and CYP genes were located on the terpene-related clusters on chromosome 13. It is suggested that chromosome 13 may play an important role in the specific terpene biosynthesis of A. indica. The gene duplication events may be responsible for the terpene biosynthesis expansion in A. indica. The genomic dataset and genomic analysis created for A. indica will shed light on terpene biosynthesis in A. indica and facilitate comparative genomic research of the family Meliaceae.

Evolutionary History and Functional Diversification of the JmjC Domain-Containing Histone Demethylase Gene Family in Plants

Histone demethylases containing JumonjiC (JmjC) domains regulate gene transcription and chromatin structure by changing the methylation status of lysine residues and play an important role in plant growth and development. In this study, a total of 332 JmjC family genes were identified from 21 different plant species. The evolutionary analysis results showed that the JmjC gene was detected in each species, that is, the gene has already appeared in algae. The phylogenetic analysis showed that the KDM3/JHDM2 subfamily genes may have appeared when plants transitioned from water to land, but were lost in lycophytes (Selaginella moellendorffii). During the evolutionary process, some subfamily genes may have been lost in individual species. According to the analysis of the conserved domains, all of the plant JmjC genes contained a typical JmjC domain, which was highly conserved during plant evolution. The analysis of cis-acting elements showed that the promoter region of the JmjC gene was rich in phytohormones and biotic and abiotic stress-related elements. The transcriptome data analysis and protein interaction analyses showed that JmjC genes play an important role in plant growth and development. The results clarified the evolutionary history of JmjC family genes in plants and lay the foundation for the analysis of the biological functions of JmjC family genes.

Genome-Wide Identification of TLP Gene Family and Their Roles in Carya cathayensis Sarg in Response to Botryosphaeria dothidea

Hickory (Carya cathayensis) is a critical tree species of the genus Carya from the Juglandaceae family that contains nutrient-rich nuts. Due to large-scale soil degradation, the pests and diseases of hickory are becoming more and more serious. Thaumatin-like proteins (TLPs) are vital proteins involved in the complex defense process of plant pathogens. In this study, 40 CcTLP genes were identified genome-widely and phylogenetically grouped into three subfamilies. The sequence of CcTLPs had a conservative pattern, such as eight stable disulfide bonds, REDDD, and G-X-[GF]-X-C-X-T-[GA]-D-C-X(1,2)-G-X-(2,3)-C structure. In total, 57 cis-elements related to stress-responsive, light-responsive, phytohormone-responsive, and plant-responsive were discovered. Under salicylate (SA), methyl jasmonate (MeJA), and ethephon (ETH) treatments, the expressions of CcTLP28, CcTLP29, CcTLP30, CcTLP31, CcTLP32, CcTLP33, CcTLP37, CcTLP38, and CcTLP39 had different patterns. This is an indication that most of the TLP genes were upregulated by SA and downregulated by MeJA. Notably, seven TLP genes were significantly upregulated under the Botryosphaeria dothidea inoculation, especially CcTLP31, with an over 20-fold change. Nine genes were shown by subcellular localization analysis to be located at the plasma membrane and cytoplasm. The knowledge of the disease-resistant function of the CcTLP family in hickory is promoted by these results. A foundation reference for the molecular breeding of this plant in the future is provided by our findings.

Diversification of the phytophagous lineages of true bugs (Insecta: Hemiptera: Heteroptera) shortly after that of the flowering plants

More than 95% of phytophagous true bug (Hemiptera: Heteroptera) species belong to four superfamilies: Miroidea (Cimicomorpha), Pentatomoidea, Coreoidea, and Lygaeoidea (all Pentatomomorpha). These iconic groups of highly diverse, overwhelmingly phytophagous insects include several economically prominent agricultural and silvicultural pest species, though their evolutionary history has not yet been well resolved. In particular, superfamily- and family-level phylogenetic relationships of these four lineages have remained controversial, and the divergence times of some crucial nodes for phytophagous true bugs have hitherto been little known, which hampers a better understanding of the evolutionary processes and patterns of phytophagous insects. In the present study, we used 150 species and concatenated nuclear and mitochondrial protein-coding genes and rRNA genes to infer the phylogenetic relationships within the Terheteroptera (Cimicomorpha + Pentatomomorpha) and estimated their divergence times. Our results support the monophyly of Cimicomorpha, Pentatomomorpha, Miroidea, Pentatomoidea, Pyrrhocoroidea, Coreoidea, and Lygaeoidea. The phylogenetic relationships across phytophagous lineages are largely congruent at deep nodes across the analyses based on different datasets and tree-reconstructing methods with just a few exceptions. Estimated divergence times and ancestral state reconstructions for feeding habit indicate that phytophagous true bugs explosively radiated in the Early Cretaceous-shortly after the angiosperm radiation-with the subsequent diversification of the most speciose clades (Mirinae, Pentatomidae, Coreinae, and Rhyparochromidae) in the Late Cretaceous.

Floral Scents and Fruit Aromas: Functions, Compositions, Biosynthesis, and Regulation

Floral scents and fruit aromas are crucial volatile organic compounds (VOCs) in plants. They are used in defense mechanisms, along with mechanisms to attract pollinators and seed dispersers. In addition, they are economically important for the quality of crops, as well as quality in the perfume, cosmetics, food, drink, and pharmaceutical industries. Floral scents and fruit aromas share many volatile organic compounds in flowers and fruits. Volatile compounds are classified as terpenoids, phenylpropanoids/benzenoids, fatty acid derivatives, and amino acid derivatives. Many genes and transcription factors regulating the synthesis of volatiles have been discovered. In this review, we summarize recent progress in volatile function, composition, biosynthetic pathway, and metabolism regulation. We also discuss unresolved issues and research perspectives, providing insight into improvements and applications of plant VOCs.

MADS-box transcription factors determine the duration of temporary winter dormancy in closely related evergreen and deciduous Iris spp

Winter dormancy (WD) is a crucial strategy for plants coping with potentially deadly environments. In recent decades, this process has been extensively studied in economically important perennial eudicots due to changing climate. However, in evergreen monocots with no chilling requirements, dormancy processes are so far a mystery. In this study, we compared the WD process in closely related evergreen (Iris japonica) and deciduous (I. tectorum) iris species across crucial developmental time points. Both iris species exhibit a 'temporary' WD process with distinct durations, and could easily resume growth under warm conditions. To decipher transcriptional changes, full-length sequencing for evergreen iris and short read RNA sequencing for deciduous iris were applied to generate respective reference transcriptomes. Combining results from a multipronged approach, SHORT VEGETATIVE PHASE and FRUITFULL (FUL) from MADS-box was associated with a dormancy- and a growth-related module, respectively. They were co-expressed with genes involved in phytohormone signaling, carbohydrate metabolism, and environmental adaptation. Also, gene expression patterns and physiological changes in the above pathways highlighted potential abscisic acid and jasmonic acid antagonism in coordinating growth and stress responses, whereas differences in carbohydrate metabolism and reactive oxygen species scavenging might lead to species-specific WD durations. Moreover, a detailed analysis of MIKCCMADS-box in irises revealed common features described in eudicots as well as possible new roles for monocots during temporary WD, such as FLOWERING LOCUS C and FUL. In essence, our results not only provide a portrait of temporary WD in perennial monocots but also offer new insights into the regulatory mechanism underlying WD in plants.

Genome and systems biology of Melilotus albus provides insights into coumarins biosynthesis

Melilotus species are used as green manure and rotation crops worldwide and contain abundant pharmacologically active coumarins. However, there is a paucity of information on its genome and coumarin production and function. Here, we reported a chromosome-scale assembly of Melilotus albus genome with 1.04 Gb in eight chromosomes, containing 71.42% repetitive elements. Long terminal repeat retrotransposon bursts coincided with declining of population sizes during the Quaternary glaciation. Resequencing of 94 accessions enabled insights into genetic diversity, population structure, and introgression. Melilotus officinalis had relatively larger genetic diversity than that of M. albus. The introgression existed between M. officinalis group and M. albus group, and gene flows was from M. albus to M. officinalis. Selection sweep analysis identified candidate genes associated with flower colour and coumarin biosynthesis. Combining genomics, BSA, transcriptomics, metabolomics, and biochemistry, we identified a β-glucosidase (BGLU) gene cluster contributing to coumarin biosynthesis. MaBGLU1 function was verified by overexpression in M. albus, heterologous expression in Escherichia coli, and substrate feeding, revealing its role in scopoletin (coumarin derivative) production and showing that nonsynonymous variation drives BGLU enzyme activity divergence in Melilotus. Our work will accelerate the understanding of biologically active coumarins and their biosynthetic pathways, and contribute to genomics-enabled Melilotus breeding.

Comprehensive Analysis of Endogenous Volatile Compounds, Transcriptome, and Enzyme Activity Reveals PmCAD1 Involved in Cinnamyl Alcohol Synthesis in Prunus mume

Floral scent is an important economic and ornamental trait of Prunus mume. The floral volatiles from most cultivars of P. mume in composition exist significant differences. Cinnamyl alcohol was one of the main floral volatile compounds with distinct abundances in different cultivars, namely, 'Zaohua Lve,' 'Zao Yudie,' 'Fenpi Gongfen,' 'Jiangsha Gongfen,' and 'Fenhong Zhusha.' Based on the determination of endogenous volatiles of full-blooming flowers, vital enzyme activity and transcriptomes were comprehensively analyzed to screen the key potential genes involved in cinnamyl alcohol synthesis. Transcriptome combining with enzyme activity level analysis suggested that the expression levels of three PmCADs were highly correlated with the cinnamyl alcohol dehydrogenase (CAD) enzyme activities in six cultivars. Furthermore, phylogenetic tree and transcriptome analysis suggested that PmCAD1 and PmCAD2 might contribute to the cinnamyl alcohol synthesis. Relative expression analyses and enzyme activity assays showed that PmCAD1 played an important role in cinnamyl alcohol biosynthesis in vitro. Overall, this research lays a theoretical foundation for clarifying comprehensively the molecular biosynthesis mechanism of floral volatiles in P. mume.

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.

Comparative Analysis of the MADS-Box Genes Revealed Their Potential Functions for Flower and Fruit Development in Longan (Dimocarpus longan)

Longan (Dimocarpus longan Lour.) is an important economic crop widely planted in tropical and subtropical regions, and flower and fruit development play decisive effects on the longan yield and fruit quality formation. MCM1, AGAMOUS, DEFICIENS, Serum Response Factor (MADS)-box transcription factor family plays important roles for the flowering time, floral organ identity, and fruit development in plants. However, there is no systematic information of MADS-box family in longan. In this study, 114 MADS-box genes were identified from the longan genome, phylogenetic analysis divided them into type I (Mα, Mβ, Mγ) and type II (MIKC*, MIKC ^C ) groups, and MIKC ^C genes were further clustered into 12 subfamilies. Comparative genomic analysis of 12 representative plant species revealed the conservation of type II in Sapindaceae and analysis of cis-elements revealed that Dof transcription factors might directly regulate the MIKC ^C genes. An ABCDE model was proposed for longan based on the phylogenetic analysis and expression patterns of MADS-box genes. Transcriptome analysis revealed that MIKC ^C genes showed wide expression spectrums, particularly in reproductive organs. From 35 days after KClO₃ treatment, 11 MIKC genes were up-regulated, suggesting a crucial role in off-season flower induction, while DlFLC, DlSOC1, DlSVP, and DlSVP-LIKE may act as the inhibitors. The gene expression patterns of longan fruit development indicated that DlSTK, DlSEP1/2, and DlMADS53 could be involved in fruit growth and ripening. This paper carried out the whole genome identification and analysis of the longan MADS-box family for the first time, which provides new insights for further understanding its function in flowers and fruit.

Characterization of Terpene synthase variation in flowers of wild aquilegia species from Northeastern Asia

There are several causes for the great diversity in floral terpenes. The terpene products are determined by the catalytic fidelity, efficiency and plasticity of the active sites of terpene synthases (TPSs). However, the molecular mechanism of TPS in catalyzing terpene biosynthesis and its evolutionary fate in wild plant species remain largely unknown. In this study, the functionality of terpene synthases and their natural variants were assessed in two Northeastern Asia endemic columbine species and their natural hybrid. Synoptically, TPS7, TPS8, and TPS9 were highly expressed in these Aquilegia species from the Zuojia population. The in vitro and in vivo enzymatic assays revealed that TPS7 and TPS8 mainly produced (+)-limonene and β-sesquiphellandrene, respectively, whereas TPS9 produced pinene, similar to the major components released from Aquilegia flowers. Multiple sequence alignment of Aquilegia TPS7 and TPS8 in the Zuojia population revealed amino acid polymorphisms. Domain swapping and amino acid substitution assays demonstrated that 413A, 503I and 529D had impacts on TPS7 catalytic activity, whereas 420G, 538F and 545 L affected the ratio of β-sesquiphellandrene to β-bisabolene in TPS8. Moreover, these key polymorphic amino acid residues were found in Aquilegia species from the Changbai Mountain population. Interestingly, amino acid polymorphisms in TPSs were present in individuals with low expression levels, and nonsynonymous mutations could impact the catalytic activity or product specificity of these genes. The results of this study will shed new light on the function and evolution of TPS genes in wild plant species and are beneficial to the modification of plant fragrances.

The chrysanthemum lavandulifolium genome and the molecular mechanism underlying diverse capitulum types

Cultivated chrysanthemum (Chrysanthemum × morifolium Ramat.) is a beloved ornamental crop due to the diverse capitula types among varieties, but the molecular mechanism of capitulum development remains unclear. Here, we report a 2.60 Gb chromosome-scale reference genome of C. lavandulifolium, a wild Chrysanthemum species found in China, Korea and Japan. The evolutionary analysis of the genome revealed that only recent tandem duplications occurred in the C. lavandulifolium genome after the shared whole genome triplication (WGT) in Asteraceae. Based on the transcriptomic profiling of six important developmental stages of the radiate capitulum in C. lavandulifolium, we found genes in the MADS-box, TCP, NAC and LOB gene families that were involved in disc and ray floret primordia differentiation. Notably, NAM and LOB30 homologs were specifically expressed in the radiate capitulum, suggesting their pivotal roles in the genetic network of disc and ray floret primordia differentiation in chrysanthemum. The present study not only provides a high-quality reference genome of chrysanthemum but also provides insight into the molecular mechanism underlying the diverse capitulum types in chrysanthemum.

Differential efficacy of water lily cultivars in phytoremediation of eutrophic water contaminated with phosphorus and nitrogen

Excessive inputs of phosphorus (P) and nitrogen (N) trigger eutrophication of the water bodies, which promotes the undesirable growth of algal bloom and deterioration of the water quality, and aquatic biodiversity. Macrophytes provide an environmentally benign and economically viable paradigm for the ecological restoration of eutrophic aquatic ecosystems. Water lily (Nymphaea) is largely used as ornamental plant for landscaping, and it has been documented that water lily possesses the potentiality in mitigating polluted aquatic environments. In the present study, water lily cultivars Nymphaea Texas Dawn (NTD), Nymphaea Colorado (NC), Nymphaea Madame Wilfron Gonnère (NMWG), and Nymphaea Sunshine Princess (NSP) were investigated for their potency in alleviating the eutrophication. The concentrations of total P and total N were significantly higher in the leaves of NC and NSP compared with NTD and NMWG. Therefore, NC and NSP were selected for subsequent studies to decipher their recuperation efficacy on eutrophic waters at different growth stages. NC and NSP significantly reduced the concentrations of eutrophication indicators i.e., total P, NH₄⁺-N, and chemical oxygen demand in different gradients of the simulated eutrophic water in a growth-dependent manner. On the contrary, NC and NSP triggered a significant increase in the concentration of dissolved oxygen particularly at the seedling stage. Notably, the concentrations of total P (shoot and root) and total N (root) were relatively higher in NSP than NC. The study thus revealed a growth-dependent differential efficacy of NSP and NC in mitigating the different eutrophic waters.

RNA-seq and phytohormone analysis reveals the culm color variation of Bambusa oldhamii Munro

BACKGROUND

The clumping bamboo Bambusa oldhamii Munro, known as "green bamboo", is famous for its edible bamboo shoots and fast-growing timber. The green and yellow striped-culm B. oldhamii variety, named B. oldhamii f. revoluta W.T. Lin & J. Y. Lin, is an attractive system for researching the culm color variation of B. oldhamii.

METHODS

Millions of clean reads were generated and assembled into 604,900 transcripts, and 383,278 unigenes were acquired with RNA-seq technology. The quantification of ABA, IAA, JA, GA₁, GA₃, GA₄, and GA₇ was performed using HPLC-MS/MS platforms.

RESULTS

Differential expression analysis showed that 449 unigenes were differentially expressed genes (DEGs), among which 190 DEGs were downregulated and 259 DEGs were upregulated in B. oldhamii f. revoluta. Phytohormone contents, especially GA₁ and GA_7, were higher in B. oldhamii. Approximately 21 transcription factors (TFs) were differentially expressed between the two groups: the bZIP, MYB, and NF-YA transcription factor families had the most DEGs, indicating that those TFs play important roles in B. oldhamii culm color variation. RNA-seq data were confirmed by quantitative RT-PCR analysis of the selected genes; moreover, phytohormone contents, especially those of ABA, GA₁ and GA₇, were differentially accumulated between the groups. Our study provides a basal gene expression and phytohormone analysis of B. oldhamii culm color variation, which could provide a solid fundamental theory for investigating bamboo culm color variation.

Genome-wide identification of the Liriodendron chinense WRKY gene family and its diverse roles in response to multiple abiotic stress

BACKGROUND

Liriodendron chinense (Lchi) is a tree species within the Magnoliaceae family and is considered a basal angiosperm. The too low or high temperature or soil drought will restrict its growth as the adverse environmental conditions, thus improving L. chinense abiotic tolerance was the key issues to study. WRKYs are a major family of plant transcription factors known to often be involved in biotic and abiotic stress responses. So far, it is still largely unknown if and how the LchiWRKY gene family is tied to regulating L. chinense stress responses. Therefore, studying the involvement of the WRKY gene family in abiotic stress regulation in L. chinense could be very informative in showing how this tree deals with such stressful conditions.

RESULTS

In this research, we performed a genome-wide analysis of the Liriodendron chinense (Lchi) WRKY gene family, studying their classification relationships, gene structure, chromosomal locations, gene duplication, cis-element, and response to abiotic stress. The 44 members of the LchiWRKY gene family contain a significant amount of sequence diversity, with their lengths ranging from 525 bp to 40,981 bp. Using classification analysis, we divided the 44 LchiWRKY genes into three phylogenetic groups (I, II, II), with group II then being further divided into five subgroups (IIa, IIb, IIc, IId, IIe). Comparative phylogenetic analysis including the WRKY families from 17 plant species suggested that LchiWRKYs are closely related to the Magnolia Cinnamomum kanehirae WRKY family, and has fewer family members than higher plants. We found the LchiWRKYs to be evenly distributed across 15 chromosomes, with their duplication events suggesting that tandem duplication may have played a major role in LchiWRKY gene expansion model. A Ka/Ks analysis indicated that they mainly underwent purifying selection and distributed in the group IId. Motif analysis showed that LchiWRKYs contained 20 motifs, and different phylogenetic groups contained conserved motif. Gene ontology (GO) analysis showed that LchiWRKYs were mainly enriched in two categories, i.e., biological process and molecular function. Two group IIc members (LchiWRKY10 and LchiWRKY37) contain unique WRKY element sequence variants (WRKYGKK and WRKYGKS). Gene structure analysis showed that most LchiWRKYs possess 3 exons and two different types of introns: the R- and V-type which are both contained within the WRKY domain (WD). Additional promoter cis-element analysis indicated that 12 cis-elements that play different functions in environmental adaptability occur across all LchiWRKY groups. Heat, cold, and drought stress mainly induced the expression of group II and I LchiWRKYs, some of which had undergone gene duplication during evolution, and more than half of which had three exons. LchiWRKY33 mainly responded to cold stress and LchiWRKY25 mainly responded to heat stress, and LchiWRKY18 mainly responded to drought stress, which was almost 4-fold highly expressed, while 5 LchiWRKYs (LchiWRKY5, LchiWRKY23, LchiWRKY14, LchiWRKY27, and LchiWRKY36) responded equally three stresses with more than 6-fold expression. Subcellular localization analysis showed that all LchiWRKYs were localized in the nucleus, and subcellular localization experiments of LchiWRKY18 and 36 also showed that these two transcription factors were expressed in the nucleus.

CONCLUSIONS

This study shows that in Liriodendron chinense, several WRKY genes like LchiWRKY33, LchiWRKY25, and LchiWRKY18, respond to cold or heat or drought stress, suggesting that they may indeed play a role in regulating the tree's response to such conditions. This information will prove a pivotal role in directing further studies on the function of the LchiWRKY gene family in abiotic stress response and provides a theoretical basis for popularizing afforestation in different regions of China.

GreeNC 2.0: a comprehensive database of plant long non-coding RNAs

The Green Non-Coding Database (GreeNC) is one of the reference databases for the study of plant long non-coding RNAs (lncRNAs). Here we present our most recent update where 16 species have been updated, while 78 species have been added, resulting in the annotation of more than 495 000 lncRNAs. Moreover, sequence clustering was applied providing information about sequence conservation and gene families. The current version of the database is available at: http://greenc.sequentiabiotech.com/wiki2/Main_Page.

Two divergent haplotypes from a highly heterozygous lychee genome suggest independent domestication events for early and late-maturing cultivars

Lychee is an exotic tropical fruit with a distinct flavor. The genome of cultivar ‘Feizixiao’ was assembled into 15 pseudochromosomes, totaling ~470 Mb. High heterozygosity (2.27%) resulted in two complete haplotypic assemblies. A total of 13,517 allelic genes (42.4%) were differentially expressed in diverse tissues. Analyses of 72 resequenced lychee accessions revealed two independent domestication events. The extremely early maturing cultivars preferentially aligned to one haplotype were domesticated from a wild population in Yunnan, whereas the late-maturing cultivars that mapped mostly to the second haplotype were domesticated independently from a wild population in Hainan. Early maturing cultivars were probably developed in Guangdong via hybridization between extremely early maturing cultivar and late-maturing cultivar individuals. Variable deletions of a 3.7 kb region encompassed by a pair of CONSTANS-like genes probably regulate fruit maturation differences among lychee cultivars. These genomic resources provide insights into the natural history of lychee domestication and will accelerate the improvement of lychee and related crops.

Two divergent haplotypes from a highly heterozygous lychee genome of the cultivar ‘Feizixiao’ and resequencing of 72 lychee accessions provide insights into the genome evolution and domestication history of lychee.

Genome sequencing and population resequencing provide insights into the genetic basis of domestication and diversity of vegetable soybean

Vegetable soybean is one of the most important vegetables in China, and the demand for this vegetable has markedly increased worldwide over the past two decades. Here, we present a high-quality de novo genome assembly of the vegetable soybean cultivar Zhenong 6 (ZN6), which is one of the most popular cultivars in China. The 20 pseudochromosomes cover 94.57% of the total 1.01 Gb assembly size, with contig N50 of 3.84 Mb and scaffold N50 of 48.41 Mb. A total of 55 517 protein-coding genes were annotated. Approximately 54.85% of the assembled genome was annotated as repetitive sequences, with the most abundant long terminal repeat transposable elements. Comparative genomic and phylogenetic analyses with grain soybean Williams 82, six other Fabaceae species and Arabidopsis thaliana genomes highlight the difference of ZN6 with other species. Furthermore, we resequenced 60 vegetable soybean accessions. Alongside 103 previously resequenced wild soybean and 155 previously resequenced grain soybean accessions, we performed analyses of population structure and selective sweep of vegetable, grain, and wild soybean. They were clearly divided into three clades. We found 1112 and 1047 genes under selection in the vegetable soybean and grain soybean populations compared with the wild soybean population, respectively. Among them, we identified 134 selected genes shared between vegetable soybean and grain soybean populations. Additionally, we report four sucrose synthase genes, one sucrose-phosphate synthase gene, and four sugar transport genes as candidate genes related to important traits such as seed sweetness and seed size in vegetable soybean. This study provides essential genomic resources to promote evolutionary and functional genomics studies and genomically informed breeding for vegetable soybean.

The complete chloroplast genome of Nymphaea thermarum (Nymphaeaceae) from Rwanda, Africa

Nymphaea thermarum is classified in the Nymphaeaceae, and is the smallest water lily in the world. It has been extinct its native environment and needs urgent protection. Here, we report and characterize the complete chloroplast genome of N. thermarum. The total length of the chloroplast genome is 159,849 bp and the GC content is 39.2% (A: 30.1%, C: 20.0%, G: 19.2%, T: 30.8%). The chloroplast genome consists of 8 rRNA, 37 tRNA, and 85 protein-coding genes. Phylogenetic analysis of N. thermarum fully resolved this taxon in a clade with Nymphaea capensis. The chloroplast genome of N. thermarum provides scientific guidance for its conservation genetics and also contributes genome resources for the phylogenetic relationship of Nymphaea.

Gramene: A Resource for Comparative Analysis of Plants Genomes and Pathways

Gramene is an integrated bioinformatics resource for accessing, visualizing, and comparing plant genomes and biological pathways. Originally targeting grasses, Gramene has grown to host annotations for over 90 plant genomes including agronomically important cereals (e.g., maize, sorghum, wheat, teff), fruits and vegetables (e.g., apple, watermelon, clementine, tomato, cassava), specialty crops (e.g., coffee, olive tree, pistachio, almond), and plants of special or emerging interest (e.g., cotton, tobacco, cannabis, or hemp). For some species, the resource includes multiple varieties of the same species, which has paved the road for the creation of species-specific pan-genome browsers. The resource also features plant research models, including Arabidopsis and C4 warm-season grasses and brassicas, as well as other species that fill phylogenetic gaps for plant evolution studies. Its strength derives from the application of a phylogenetic framework for genome comparison and the use of ontologies to integrate structural and functional annotation data. This chapter outlines system requirements for end-users and database hosting, data types and basic navigation within Gramene, and provides examples of how to (1) explore Gramene's search results, (2) explore gene-centric comparative genomics data visualizations in Gramene, and (3) explore genetic variation associated with a gene locus. This is the first publication describing in detail Gramene's integrated search interface-intended to provide a simplified entry portal for the resource's main data categories (genomic location, phylogeny, gene expression, pathways, and external references) to the most complete and up-to-date set of plant genome and pathway annotations.

Plant protein-coding gene families: Their origin and evolution

Steady advances in genome sequencing methods have provided valuable insights into the evolutionary processes of several gene families in plants. At the core of plant biodiversity is an extensive genetic diversity with functional divergence and expansion of genes across gene families, representing unique phenomena. The evolution of gene families underpins the evolutionary history and development of plants and is the subject of this review. We discuss the implications of the molecular evolution of gene families in plants, as well as the potential contributions, challenges, and strategies associated with investigating phenotypic alterations to explain the origin of plants and their tolerance to environmental stresses.

Plant genome sequence assembly in the era of long reads: Progress, challenges and future directions

Third-generation long-read sequencing is transforming plant genomics. Oxford Nanopore Technologies and Pacific Biosciences are offering competing long-read sequencing technologies and enable plant scientists to investigate even large and complex plant genomes. Sequencing projects can be conducted by single research groups and sequences of smaller plant genomes can be completed within days. This also resulted in an increased investigation of genomes from multiple species in large scale to address fundamental questions associated with the origin and evolution of land plants. Increased accessibility of sequencing devices and user-friendly software allows more researchers to get involved in genomics. Current challenges are accurately resolving diploid or polyploid genome sequences and better accounting for the intra-specific diversity by switching from the use of single reference genome sequences to a pangenome graph.

The investigation of intraspecific characteristics and comparative analyses of the complete mitochondrial genome of Stegobium paniceum (Linnaeus, 1758) (Coleoptera: Ptinidae) assembled from public NGS raw reads of the black truffle, Tuber melanosporum

Due to the rapid development of NGS technologies, a huge amount of NGS raw reads have been accumulated in public repositories, such as the Short Read Archive of NCBI. We successfully rescued the complete mitochondrial genome of Stegobium paniceum, a drug store beetle, from public NGS raw reads of truffle generated from the whole genome project. The circular mitogenome of S. paniceum is 15,474 bp long including 13 protein-coding genes, two ribosomal RNA genes, 22 transfer RNAs, and a single large non-coding region of 803 bp. Intraspecific transfer RNAs structure and sequence variations were investigated and simple sequence repeats identified from three S. paniceum mitochondrial genomes were compared showing their diversities as fundamental data to utilize them in various aspects including developing efficient molecular markers in the family, Ptinidae. Phylogenetic analysis of 23 Bostrichoidea mitochondrial genomes presented better species identification based on phylogenetic analyses and the optimal options for constructing phylogenetic trees based on Bostrichoidea mitochondrial genomes. Our results present not only utilization of public NGS raw read sequences but also intraspecific features of S. paniceum mitochondrial genomes and comparative analysis of Bostrichoidea mitochondrial genomes in various aspects.

Combining novel technologies with interdisciplinary basic research to enhance horticultural crops

Horticultural crops mainly include fruits, vegetables, ornamental trees and flowers, and tea trees (Melaleuca alternifolia). They produce a variety of nutrients for the daily human diet in addition to the nutrition provided by staple crops, and some of them additionally possess ornamental and medicinal features. As such, horticultural crops make unique and important contributions to both food security and a colorful lifestyle. Under the current climate change scenario, the growing population and limited arable land means that agriculture, and especially horticulture, has been facing unprecedented challenges to meet the diverse demands of human daily life. Breeding horticultural crops with high quality and adaptability and establishing an effective system that combines cultivation, post-harvest handling, and sales becomes increasingly imperative for horticultural production. This review discusses characteristic and recent research highlights in horticultural crops, focusing on the breeding of quality traits and the mechanisms that underpin them. It additionally addresses challenges and potential solutions in horticultural production and post-harvest practices. Finally, we provide a prospective as to how emerging technologies can be implemented alongside interdisciplinary basic research to enhance our understanding and exploitation of horticultural crops.

Pathogen-driven coevolution across the CBP60 plant immune regulator subfamilies confers resilience on the regulator module

Components of the plant immune signaling network need mechanisms that confer resilience against fast-evolving pathogen effectors that target them. Among eight Arabidopsis CaM-Binding Protein (CBP) 60 family members, AtCBP60g and AtSARD1 are partially functionally redundant, major positive immune regulators, and AtCBP60a is a negative immune regulator. We investigated possible resilience-conferring evolutionary mechanisms among the CBP60a, CBP60g and SARD1 immune regulatory subfamilies. Phylogenetic analysis was used to investigate the times of CBP60 subfamily neofunctionalization. Then, using the pairwise distance rank based on the newly developed analytical platform Protein Evolution Analysis in a Euclidean Space (PEAES), hypotheses of specific coevolutionary mechanisms that could confer resilience on the regulator module were tested. The immune regulator subfamilies diversified around the time of angiosperm divergence and have been evolving very quickly. We detected significant coevolutionary interactions across the immune regulator subfamilies in all of 12 diverse core eudicot species lineages tested. The coevolutionary interactions were consistent with the hypothesized coevolution mechanisms. Despite their unusually fast evolution, members across the CBP60 immune regulator subfamilies have influenced the evolution of each other long after their diversification in a way that could confer resilience on the immune regulator module against fast-evolving pathogen effectors.

The biochemical and molecular investigation of flower color and scent sheds lights on further genetic modification of ornamental traits in Clivia miniata

Clivia miniata is renowned for its evergreen and strap-like leaves, whereas its floral color and scent are lacking diversity. Here, anthocyanin, volatile terpene, and carotenoid metabolisms were integrally investigated in C. miniata flowers. The results showed that pelargonidins and lutein might cooperate to confer orange or yellow color to C. miniata flowers, but only a trace amount of (+)-limonene was detected. The expression levels of CmF3'H and CmDFR appeared to be responsible for the ratio of cyanidin and pelargonidin derivatives in C. miniata, and the low expression of CmF3'H was responsible for the lack of cyanidins in flowers. Moreover, the CmF3'H promoter could not be activated by CmMYBAs, suggesting that it was controlled by novel regulators. Only two CmTPSs were functional, with CmTPS2 responsible for (+)-limonene synthesis, contributing to the monotonous flower volatile terpenes of C. miniata. CmCCD1a and CmCCD1b were able to cleave carotenoids at the 5,6 (5',6'), and 9,10 (9',10') positions to generate volatile apocarotenoids, whereas the substrates found in low-quantities or specific subcellular localizations of CmCCD1s might constrain volatile apocarotenoid release. Consequently, activating F3'H and introducing novel F3'5'H or versatile TPS may be effective ways to modify the floral color and scent, respectively. Alternatively, modifying the carotenoid flux or CCD1 localization might affect floral color and scent simultaneously. Taking these results together, the present study provides a preliminary deciphering of the genetic constraints underlying flower color and scent development, and proposes possible schemes for further genetic modification of ornamental traits in C. miniata and other plants.

Evolutionary analysis of angiosperm dehydrin gene family reveals three orthologues groups associated to specific protein domains

Dehydrins (DHNs) are a family of plant proteins that play important roles on abiotic stress tolerance and seed development. They are classified into five structural subgroups: K-, SK-, YK-, YSK-, and KS-DHNs, according to the presence of conserved motifs named K-, Y- and S- segments. We carried out a comparative structural and phylogenetic analysis of these proteins, focusing on the less-studied KS-type DHNs. A search for conserved motifs in DHNs from 56 plant genomes revealed that KS-DHNs possess a unique and highly conserved N-terminal, 15-residue amino acid motif, not previously described. This novel motif, that we named H-segment, is present in DHNs of angiosperms, gymnosperms and lycophytes, suggesting that HKS-DHNs were present in the first vascular plants. Phylogenetic and microsynteny analyses indicate that the five structural subgroups of angiosperm DHNs can be assigned to three groups of orthologue genes, characterized by the presence of the H-, F- or Y- segments. Importantly, the hydrophilin character of DHNs correlate with the phylogenetic origin of the DHNs rather than to the traditional structural subgroups. We propose that angiosperm DHNs can be ultimately subdivided into three orthologous groups, a phylogenetic framework that should help future studies on the evolution and function of this protein family.

The evolution of N6-methyladenosine regulators in plants

As a reversible modification, N⁶-methyladenosine (m⁶A) plays key roles in series of biological processes. Although it has been found that m⁶A modification is regulated by writers, erasers and readers, their evolutionary processes are still not clearly and systematically described. In the present work, we identified 1592 m⁶A modification regulators from 65 representative plant species and performed the phylogenetic relationships, sequence structure, selection pressure, and codon usage analysis across species. The regulators from different species or subfamilies were distinguishable based on the phylogenetic trees. Although the gene structure was structurally and functionally conserved for each kind of regulators, the unique exon/intron structures and motif organizations were observed among different families. The selection pressure analysis demonstrated that the regulators experienced purifying selection. Interestingly, the selection pressure for the regulators in higher plants was more relaxed, indicating that they might have acquired new functions during evolution. In addition, the different codon usage preferences were observed for the different kinds of m⁶A modification regulators. These results will not only facilitate our understanding of the evolution of m⁶A regulators, but also shed light on how the evolutionary differences affect their functional divergence.

The red flower wintersweet genome provides insights into the evolution of magnoliids and the molecular mechanism for tepal color development

Wintersweet (Chimonanthus praecox) is one of the most important ornamental plants. Its color is mainly determined by the middle tepals. However, the molecular mechanisms underlying the intriguing flower color development among different wintersweet groups are still largely unknown. In addition, wintersweet belongs to magnoliids, and the phylogenetic position of magnoliids remains to be determined conclusively. Here, the whole genome of red flower wintersweet, a new wintersweet type, was sequenced and assembled with high quality. The genome comprised 11 super-scaffolds (chromosomes) with a total size of 737.03 Mb. Based on the analyses of the long branch attraction, incomplete lineage sorting, sparse taxon sampling, and other factors, we suggest that a bifurcating tree may not fully represent the complex early diversification of the angiosperms and that magnoliids are most likely sister to the eudicots. The wintersweet genome appears to have undergone two whole-genome duplication (WGD) events: a recent WGD event representing an independent event specific to the Calycanthaceae and an ancient WGD event shared by Laurales. By integrating genomic, transcriptomic, and metabolomic data, CpANS1 and the transcription factor CpMYB1 were found to play key roles in regulating tepal color development, whereas CpMYB1 needs to form a complex with bHLH and WD40 to fully perform its regulatory function. The present study not only provides novel insights into the evolution of magnoliids and the molecular mechanism for flower color development, but also lays the foundation for subsequent functional genomics study and molecular breeding of wintersweet.

The evolution of imprinting in plants: beyond the seed

Genomic imprinting results in the biased expression of alleles depending on if the allele was inherited from the mother or the father. Despite the prevalence of sexual reproduction across eukaryotes, imprinting is only found in placental mammals, flowering plants, and some insects, suggesting independent evolutionary origins. Numerous hypotheses have been proposed to explain the selective pressures that favour the innovation of imprinted gene expression and each differs in their experimental support and predictions. Due to the lack of investigation of imprinting in land plants, other than angiosperms with triploid endosperm, we do not know whether imprinting occurs in species lacking endosperm and with embryos developing on maternal plants. Here, we discuss the potential for uncovering additional examples of imprinting in land plants and how these observations may provide additional support for one or more existing imprinting hypotheses.

High-quality evergreen azalea genome reveals tandem duplication-facilitated low-altitude adaptability and floral scent evolution

Azalea belongs to Rhododendron, which is one of the largest genera of flowering plants and is well known for the diversity and beauty in its more than 1000 woody species. Rhododendron contains two distinct groups: the most high-altitude and a few low-altitude species; however, the former group is difficult to be domesticated for urban landscaping, and their evolution and adaptation are little known. Rhododendron ovatum has broad adaptation in low-altitude regions but possesses evergreen characteristics like high-altitude species, and it has floral fragrance that is deficient in most cultivars. Here we report the chromosome-level genome assembly of R. ovatum, which has a total length of 549 Mb with scaffold N50 of 41 Mb and contains 41 264 predicted genes. Genomic micro-evolutionary analysis of R. ovatum in comparison with two high-altitude Rhododendron species indicated that the expansion genes in R. ovatum were significantly enriched in defence responses, which may account for its adaptability in low altitudes. The R. ovatum genome contains much more terpene synthase genes (TPSs) compared with the species that lost floral fragrance. The subfamily b members of TPS are involved in the synthesis of sesquiterpenes as well as monoterpenes and play a major role in flora scent biosynthesis and defence responses. Tandem duplication is the primary force driving expansion of defence-responsive genes for extensive adaptability to the low-altitude environments. The R. ovatum genome provides insights into low-altitude adaptation and gain or loss of floral fragrance for Rhododendron species, which are valuable for alpine plant domestication and floral scent breeding.

The MADS-box gene PpPI is a key regulator of the double-flower trait in peach

Double flower is an invaluable trait in ornamental peach, but the mechanism underlying its development remains largely unknown. Here, we report the roles of ABCE model genes in double flower development in peach. A total of nine ABCE regulatory genes, including eight MADS-box genes and one AP2/EREBP gene, were identified in the peach genome. Subcellular localization assay showed that all the ABCE proteins were localized in the nucleus. Four genes, PpAP1, PpAP3, PpSEP3, and PpPI, showed a difference in expression levels between single and double flowers. Ectopic overexpression of PpPI increased petal number in Arabidopsis, while transgenic lines overexpressing PpAP3 or PpSEP3 were morphologically similar to wild-type. Ectopic overexpression of PpAP1 resulted in a significant decrease in the number of basal leaves and caused early flowering. These results suggest that PpPI is likely crucial for double flower development in peach. In addition, double flowers have petaloid sepals and stamens, and single flower could occasionally change to be double flower by converting stamens to petals in peach, suggesting that the double-flower trait is likely to have evolved from an ancestral single-flower structure. Our results provide new insights into mechanisms underlying the double-flower trait in peach.

The Cymbidium genome reveals the evolution of unique morphological traits

The marvelously diverse Orchidaceae constitutes the largest family of angiosperms. The genus Cymbidium in Orchidaceae is well known for its unique vegetation, floral morphology, and flower scent traits. Here, a chromosome-scale assembly of the genome of Cymbidium ensifolium (Jianlan) is presented. Comparative genomic analysis showed that C. ensifolium has experienced two whole-genome duplication (WGD) events, the most recent of which was shared by all orchids, while the older event was the τ event shared by most monocots. The results of MADS-box genes analysis provided support for establishing a unique gene model of orchid flower development regulation, and flower shape mutations in C. ensifolium were shown to be associated with the abnormal expression of MADS-box genes. The most abundant floral scent components identified included methyl jasmonate, acacia alcohol and linalool, and the genes involved in the floral scent component network of C. ensifolium were determined. Furthermore, the decreased expression of photosynthesis-antennae and photosynthesis metabolic pathway genes in leaves was shown to result in colorful striped leaves, while the increased expression of MADS-box genes in leaves led to perianth-like leaves. Our results provide fundamental insights into orchid evolution and diversification.

The Chloranthus sessilifolius genome provides insight into early diversification of angiosperms

Most extant angiosperms belong to Mesangiospermae, which comprises eudicots, monocots, magnoliids, Chloranthales and Ceratophyllales. However, phylogenetic relationships between these five lineages remain unclear. Here, we report the high-quality genome of a member of the Chloranthales lineage (Chloranthus sessilifolius). We detect only one whole genome duplication within this species and find that polyploidization events in different Mesangiospermae lineage are mutually independent. We also find that the members of all floral development-related gene lineages are present in C. sessilifolius despite its extremely simplified flower. The AP1 and PI genes, however, show a weak floral tissue-specialized expression. Our phylogenomic analyses suggest that Chloranthales and magnoliids are sister groups, and both are together sister to the clade comprising Ceratophyllales and eudicots, while the monocot lineage is sister to all other Mesangiospermae. Our findings suggest that in addition to hybridization, incomplete lineage sorting may largely account for phylogenetic inconsistencies between the observed gene trees.

Chromosome-Scale Reference Genome of Amphicarpaea edgeworthii: A New Resource for Amphicarpic Plants Research and Complex Flowering Pattern

Amphicarpaea edgeworthii, an annual twining herb, is a widely distributed species and an attractive model for studying complex flowering types and evolutionary mechanisms of species. Herein, we have generated a high-quality assembly of A. edgeworthii by using a combination of PacBio, 10× Genomics libraries, and Hi-C mapping technologies. The final 11 chromosome-level scaffolds covered 90.61% of the estimated genome (343.78Mb), which is a chromosome-scale assembled genome of an amphicarpic plant. Subsequently, we characterized the genetic diversity and population structure of A. edgeworthii species by resequencing individuals collected from their natural area of distribution. Using transcriptome profiling, we observed that specific phenotypes are regulated by a complex network of light, hormones, and MADS-box gene families. These data are beneficial for the discovery of genes that control major agronomic traits and spur genetic improvement of and functional genetic studies in legumes, as well as supply comparative genetic resources for other amphicarpic plants.

High-quality reference genome sequences of two coconut cultivars provide insights into evolution of monocot chromosomes and differentiation of fiber content and plant height

BACKGROUND

Coconut is an important tropical oil and fruit crop whose evolutionary position renders it a fantastic species for the investigation of the evolution of monocot chromosomes and the subsequent differentiation of ancient plants.

RESULTS

Here, we report the assembly and annotation of reference-grade genomes of Cn. tall and Cn. dwarf, whose genome sizes are 2.40 Gb and 2.39 Gb, respectively. The comparative analysis reveals that the two coconut subspecies diverge about 2-8 Mya while the conserved Arecaceae-specific whole-genome duplication (ω WGD) occurs approximately 47-53 Mya. It additionally allows us to reconstruct the ancestral karyotypes of the ten ancient monocot chromosomes and the evolutionary trajectories of the 16 modern coconut chromosomes. Fiber synthesis genes in Cn. tall, related to lignin and cellulose synthesis, are found at a higher copy number and expression level than dwarf coconuts. Integrated multi-omics analysis reveals that the difference in coconut plant height is the result of altered gibberellin metabolism, with both the GA20ox copy number and a single-nucleotide change in the promoter together leading to the difference in plant height between Cn. tall and Cn. dwarf.

CONCLUSION

We provide high-quality coconut genomes and reveal the genetic basis of trait differences between two coconuts through multi-omics analysis. We also reveal that the selection of plant height has been targeted for the same gene for millions of years, not only in natural selection of ancient plant as illustrated in coconut, but also for artificial selection in cultivated crops such as rice and maize.

Underwater CAM photosynthesis elucidated by Isoetes genome

To conserve water in arid environments, numerous plant lineages have independently evolved Crassulacean Acid Metabolism (CAM). Interestingly, Isoetes, an aquatic lycophyte, can also perform CAM as an adaptation to low CO₂ availability underwater. However, little is known about the evolution of CAM in aquatic plants and the lack of genomic data has hindered comparison between aquatic and terrestrial CAM. Here, we investigate underwater CAM in Isoetes taiwanensis by generating a high-quality genome assembly and RNA-seq time course. Despite broad similarities between CAM in Isoetes and terrestrial angiosperms, we identify several key differences. Notably, Isoetes may have recruited the lesser-known 'bacterial-type' PEPC, along with the 'plant-type' exclusively used in other CAM and C4 plants for carboxylation of PEP. Furthermore, we find that circadian control of key CAM pathway genes has diverged considerably in Isoetes relative to flowering plants. This suggests the existence of more evolutionary paths to CAM than previously recognized.

Metabolite analysis in Nymphaea 'Blue Bird' petals reveal the roles of flavonoids in color formation, stress amelioration, and bee orientation

In this study, metabolome of open petals (OP) and closed petals (CP) from Nymphaea 'Blue Bird' was firstly investigated. A total of 455 metabolites was identified in Nymphaea 'Blue Bird' petals, which was mainly composed of 100 flavonoids, 83 phenolic acids, 64 amino acids and derivatives, 60 lipids, 32 alkaloids, 32 organic acids, 24 nucleotides and derivatives, and 12 lignans and coumarins. By differential analysis, 192 metabolites were identified with variable importance in project ≥ 1, among which 83 and 109 metabolites were up- and down-regulated in OP group, respectively. Further analysis (Log₂ fold change ≥ 1) identified 26 and 7 metabolites exhibited significantly lower and higher contents in CP group, relative to OP group. Importantly, KEGG analysis indicated that flavonoid biosynthesis exhibited the most significant enrichment. qRT-PCR analysis indicated that the PAL, CHS, and HCDBR genes showed a significantly higher expression in OP group than in CP group. These data explain the increase of naringenin chalcone and phloretin in OP. However, there was no significant difference of total flavonoids between OP and CP groups. Considering the increase of H₂O₂ content and ultraviolet (UV) absorption peak in OP, our results implied that diurnal stressful conditions induced the degradation of flavonoids, which contributed to environmental stress amelioration. Moreover, a higher absorption peak of 360-380 nm UV was observed in the extract liquor of OP. The sensitivity maximum of the UV-photoreceptor of bees is situated around 340-380 nm UV. This suggested, as noted for the maximum absorption of dihydrokaempferol in 340-370 nm, rhythmic accumulation and loss of these differential flavonoids in Nymphaea 'Blue Bird' petals might enhance UV pattern to some degree, influencing pollinator attraction.

Biosynthesis and emission of methyl hexanoate, the major constituent of floral scent of a night-blooming water lily Victoriacruziana

Among the factors that have made flowering plants the most species-rich lineage of land plants is the interaction between flower and insect pollinators, for which floral scent plays a pivotal role. Water lilies belong to the ANA (Amborellales, Nymphaeales, and Austrobaileyales) grade of basal flowering plants. In this study, Victoria cruziana was investigated as a model night-blooming water lily for floral scent biosynthesis. Four volatile compounds, including three benzenoids and one fatty acid methyl ester methyl hexanoate, were detected from the flowers of V. cruziana during their first bloom, with methyl hexanoate accounting for 45 % of total floral volatile emission. Emission rates were largely constant before significant drop starting at the end of second bloom. To understand the molecular basis of floral scent biosynthesis in V. cruziana, particularly methyl hexanoate, a transcriptome from the whole flowers at the full-bloom stage was created and analyzed. Methyl hexanoate was hypothesized to be biosynthesized by SABATH methyltransferases. From the transcriptome, three full-length SABATH genes designated VcSABATH1-3 were identified. A full-length cDNA for each of the three VcSABATH genes was expressed in Escherichia coli to produce recombinant proteins. When tested in in vitro methyltransferase enzyme assays with different fatty acids, both VcSABATH1 and VcSABATH3 exhibited highest levels of activity with hexanoic acid to produce methyl hexanoate, with the specific activity of VcSABATH1 being about 15 % of that for VcSABATH3. VcSABATH1 and VcSABATH3 showed the highest levels of expression in stamen and pistil, respectively. In phylogenetic analysis, three VcSABATH genes clustered with other water lily SABATH methyltransferase genes including the one known for making other fatty acid methyl esters, implying both a common evolutionary origin and functional divergence. Fatty acid methyl esters are not frequent constituents of floral scents of mesangiosperms, pointing to the importance for the evolution of novel fatty acid methyltransferase for making fatty acid methyl esters in the pollination biology of water lilies.

Plastid phylogenomic insights into relationships of all flowering plant families

BACKGROUND

Flowering plants (angiosperms) are dominant components of global terrestrial ecosystems, but phylogenetic relationships at the familial level and above remain only partially resolved, greatly impeding our full understanding of their evolution and early diversification. The plastome, typically mapped as a circular genome, has been the most important molecular data source for plant phylogeny reconstruction for decades.

RESULTS

Here, we assembled by far the largest plastid dataset of angiosperms, composed of 80 genes from 4792 plastomes of 4660 species in 2024 genera representing all currently recognized families. Our phylogenetic tree (PPA II) is essentially congruent with those of previous plastid phylogenomic analyses but generally provides greater clade support. In the PPA II tree, 75% of nodes at or above the ordinal level and 78% at or above the familial level were resolved with high bootstrap support (BP ≥ 90). We obtained strong support for many interordinal and interfamilial relationships that were poorly resolved previously within the core eudicots, such as Dilleniales, Saxifragales, and Vitales being resolved as successive sisters to the remaining rosids, and Santalales, Berberidopsidales, and Caryophyllales as successive sisters to the asterids. However, the placement of magnoliids, although resolved as sister to all other Mesangiospermae, is not well supported and disagrees with topologies inferred from nuclear data. Relationships among the five major clades of Mesangiospermae remain intractable despite increased sampling, probably due to an ancient rapid radiation.

CONCLUSIONS

We provide the most comprehensive dataset of plastomes to date and a well-resolved phylogenetic tree, which together provide a strong foundation for future evolutionary studies of flowering plants.

The Melastoma dodecandrum genome and the evolution of Myrtales

Melastomataceae have abundant morphological diversity with high economic and ornamental merit in Myrtales. The phylogenetic position of Myrtales is still contested. Here, we report the first chromosome-level genome assembly of Melastoma dodecandrum in Melastomataceae. The assembled genome size was 299.81 Mb with a contig N50 value of 3.00 Mb. Genome evolution analysis indicated that M. dodecandrum, Eucalyptus grandis and Punica granatum were clustered into a clade of Myrtales and formed a sister group with the ancestor of fabids and malvids. We found that M. dodecandrum experienced four whole-genome polyploidization events: the ancient event was shared with most eudicots, one event was shared with Myrtales, and the other two events were unique to M. dodecandrum. Moreover, we identified MADS-box genes and found that the AP1-like genes expanded, and AP3-like genes might have undergone subfunctionalization. We found that the SUAR63-like genes and AG-like genes showed different expression patterns in stamens, which may be associated with heteranthery. In addition, we found that LAZY1-like genes were involved in the negative regulation of stem branching development, which may be related to its creeping features. Our study sheds new light on the evolution of Melastomataceae and Myrtales, which provides a comprehensive genetic resource for future research.

Phylogeny of gracillariid leaf-mining moths: evolution of larval behaviour inferred from phylogenomic and Sanger data

Gracillariidae is the most taxonomically diverse cosmopolitan leaf-mining moth family, consisting of nearly 2000 named species in 105 described genera, classified into eight extant subfamilies. The majority of gracillariid species are internal plant feeders as larvae, creating mines and galls in plant tissue. Despite their diversity and ecological adaptations, their phylogenetic relationships, especially among subfamilies, remain uncertain. Genomic data (83 taxa, 589 loci) were integrated with Sanger data (130 taxa, 22 loci), to reconstruct a phylogeny of Gracillariidae. Based on analyses of both datasets combined and analyzed separately, monophyly of Gracillariidae and all its subfamilies, monophyly of the clade "LAMPO" (subfamilies: Lithocolletinae, Acrocercopinae, Marmarinae, Phyllocnistinae, and Oecophyllembiinae) and relationships of its subclade "AMO" (subfamilies: Acrocercopinae, Marmarinae, and Oecophyllembiinae) were strongly supported. A sister-group relationship of Ornixolinae to the remainder of the family, and a monophyletic leaf roller lineage (Callicercops Vári + Parornichinae) + Gracillariinae, as sister to the "LAMPO" clade were supported by the most likely tree. Dating analyses indicate a mid-Cretaceous (105.3 Ma) origin of the family, followed by a rapid diversification into the nine subfamilies predating the Cretaceous-Palaeogene extinction. We hypothesize that advanced larval behaviours, such as making keeled or tentiform blotch mines, rolling leaves and galling, allowed gracillariids to better avoid larval parasitoids allowing them to further diversify. Finally, we stabilize the classification by formally re-establishing the subfamily ranks of Marmarinae stat.rev., Oecophyllembiinae stat.rev. and Parornichinae stat.rev., and erect a new subfamily, Callicercopinae Li, Ohshima and Kawahara to accommodate the enigmatic genus Callicercops.

Physalis floridana CRABS CLAW mediates neofunctionalization of GLOBOSA genes in carpel development

Floral B-function MADS-box genes, such as GLOBOSA (GLO), function in corolla and stamen organ identity specification. The functions of these genes outside these floral whorls are rarely reported. DOLL1 is a GLO gene controlling corolla and androecium organ identity. In this study we found that, in Physalis floridana double-layered-lantern 1 (doll1) mutant pollinated with wild-type pollen, fruit set was extremely low, indicating that doll1 females are dysfunctional. Stigma and style structure, stigma receptivity, pollen tube guidance, and embryo sac development were also impaired in doll1. P. floridana CRABS CLAW (PFCRC), predominantly expressed in carpels, was repressed in doll1 native carpels. Loss-of-function of PFCRC altered carpel meristem determinacy, carpel closure, and ovule number, and the resultant 'pistil' consisted of multiple spirally-arranged dorsiventral carpels occasionally with 1-2 naked ovules on the margin and trichomes at each mutated carpel tip, implying an alteration of carpel organ identity. Regulatory and genetic interactions between B-class MADS-box genes and PFCRC were revealed in a context-dependent manner in floral development. Our work reveals a new role for the B-function genes in carpel and ovule development via regulating PFCRC, providing a new understanding of genetic regulatory networks between MADS-domain and CRC transcription factors in mediating carpel organ specification, functionality, and origin.

Taxus yunnanensis genome offers insights into gymnosperm phylogeny and taxol production

Taxol, a natural product derived from Taxus, is one of the most effective natural anticancer drugs and the biosynthetic pathway of Taxol is the basis of heterologous bio-production. Here, we report a high-quality genome assembly and annotation of Taxus yunnanensis based on 10.7 Gb sequences assembled into 12 chromosomes with contig N50 and scaffold N50 of 2.89 Mb and 966.80 Mb, respectively. Phylogenomic analyses show that T. yunnanensis is most closely related to Sequoiadendron giganteum among the sampled taxa, with an estimated divergence time of 133.4-213.0 MYA. As with most gymnosperms, and unlike most angiosperms, there is no evidence of a recent whole-genome duplication in T. yunnanensis. Repetitive sequences, especially long terminal repeat retrotransposons, are prevalent in the T. yunnanensis genome, contributing to its large genome size. We further integrated genomic and transcriptomic data to unveil clusters of genes involved in Taxol synthesis, located on the chromosome 12, while gene families encoding hydroxylase in the Taxol pathway exhibited significant expansion. Our study contributes to the further elucidation of gymnosperm relationships and the Taxol biosynthetic pathway.