Chromosome-level assembly of the Neolamarckia cadamba genome provides insights into the evolution of cadambine biosynthesis

Identification of the CONSTANS-like family in Cymbidium sinense, and their functional characterization

BACKGROUND

Cymbidium sinense is an orchid that is typically used as a potted plant, given its high-grade ornamental characteristics, and is most frequently distributed in China and SE Asia. The inability to strictly regulate flowering in this economically important potted and cut-flower orchid is a bottleneck that limits its industrial development. Studies on C. sinense flowering time genes would help to elucidate the mechanism regulating flowering. There are very few studies on the genetic regulation of flowering pathways in C. sinense. Photoperiod significantly affects the flowering of C. sinense, but it was unknown how the CONSTANS gene family is involved in regulating flowering.

RESULTS

In this study, eight CONSTANS-like genes were identified and cloned. They were divided into three groups based on a phylogenetic analysis. Five representative CsCOL genes (CsCOL3/4/6/8/9) were selected from the three groups to perform expression characterization and functional study. CsCOL3/4/6/8/9 are nucleus-localized proteins, and all five CsCOL genes were expressed in all organs, mainly in leaves followed by sepals. The expression levels of CsCOL3/4 (group I) were higher in all organs than other CsCOL genes. Developmental stage specific expression revealed that the expression of CsCOL3/4/9 peaked at the initial flowering stage. In contrast, the transcript level of CsCOL6/8 was highest at the pedicel development stage. Photoperiodic experiments demonstrated that the transcripts of the five CsCOL genes exhibited distinct diurnal rhythms. Under LD conditions, the overexpression of CsCOL3/4 promoted early flowering, and CsCOL6 had little effect on flowering time, whereas CsCOL8 delayed flowering of Arabidopsis thaliana. However, under SD conditions, overexpression of CsCOL4/6/8 promoted early flowering and the rosette leaves growth, and CsCOL3 induced flower bud formation in transgenic Arabidopsis.

CONCLUSION

The phylogenetic analysis, temporal and spatial expression patterns, photoperiodic rhythms and functional study indicate that CsCOL family members in C. sinense were involved in growth, development and flowering regulation through different photoperiodic pathway. The results will be useful for future research on mechanisms pertaining to photoperiod-dependent flowering, and will also facilitate genetic engineering-based research that uses Cymbidium flowering time genes.

The Rauvolfia tetraphylla genome suggests multiple distinct biosynthetic routes for yohimbane monoterpene indole alkaloids

Monoterpene indole alkaloids (MIAs) are a structurally diverse family of specialized metabolites mainly produced in Gentianales to cope with environmental challenges. Due to their pharmacological properties, the biosynthetic modalities of several MIA types have been elucidated but not that of the yohimbanes. Here, we combine metabolomics, proteomics, transcriptomics and genome sequencing of Rauvolfia tetraphylla with machine learning to discover the unexpected multiple actors of this natural product synthesis. We identify a medium chain dehydrogenase/reductase (MDR) that produces a mixture of four diastereomers of yohimbanes including the well-known yohimbine and rauwolscine. In addition to this multifunctional yohimbane synthase (YOS), an MDR synthesizing mainly heteroyohimbanes and the short chain dehydrogenase vitrosamine synthase also display a yohimbane synthase side activity. Lastly, we establish that the combination of geissoschizine synthase with at least three other MDRs also produces a yohimbane mixture thus shedding light on the complex mechanisms evolved for the synthesis of these plant bioactives.

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.

Genome-Wide Identification and Expression Analysis of WRKY Gene Family in Neolamarckia cadamba

The WRKY transcription factor family plays important regulatory roles in multiple biological processes in higher plants. They have been identified and functionally characterized in a number of plant species, but very little is known in Neolamarckia cadamba, a 'miracle tree' for its fast growth and potential medicinal resource in Southeast Asia. In this study, a total of 85 WRKY genes were identified in the genome of N. cadamba. They were divided into three groups according to their phylogenetic features, with the support of the characteristics of gene structures and conserved motifs of protein. The NcWRKY genes were unevenly distributed on 22 chromosomes, and there were two pairs of segmentally duplicated events. In addition, a number of putative cis-elements were identified in the promoter regions, of which hormone- and stress-related elements were shared in many NcWRKYs. The transcript levels of NcWRKY were analyzed using the RNA-seq data, revealing distinct expression patterns in various tissues and at different stages of vascular development. Furthermore, 16 and 12 NcWRKY genes were confirmed to respond to various hormone treatments and two different abiotic stress treatments, respectively. Moreover, the content of cadambine, the active metabolite used for the various pharmacological activities found in N. cadamba, significantly increased after Methyl jasmonate treatment. In addition, expression of NcWRKY64/74 was obviously upregulated, suggesting that they may have a potential function of regulating the biosynthesis of cadambine in response to MeJA. Taken together, this study provides clues into the regulatory roles of the WRKY gene family in N. cadamba.

Physiological and transcriptomic responses to magnesium deficiency in Neolamarckia Cadamba

Magnesium (Mg²⁺) is a critical component of chlorophyll and enzymes involved in various physiological and biochemical processes essential for plant growth, biomass accumulation, and photosynthesis. Mg²⁺ deficiency (MgD) is common in hot and rainy subtropical areas due to its easy loss from soil. Neolamarckia cadamba, an important tropical tree in South Asia, faces severe effects of MgD, however, the responses of N. cadamba to MgD stress remain unclear. In here, effects of N. cadamba under MgD stress were investigated. The study revealed that MgD had lower plant biomass, fresh and dry weight, root length, root volume, and surface area compared to CK (normal Mg²⁺). As treatment time increased, the leaves began to yellow, and lesions appeared. Chlorophyll a, chlorophyll b, and total chlorophyll content, along with fluorescence-related parameters and leaf photosynthetic capacity, were significantly reduced in MgD stress compared to CK treatment. Transcriptome analysis showed that transporters as well as transcription factors (TFs) from MYC (v-myc avian myelocytomatosis viral oncogene homolog), MYB (v-myb avian myeloblastosis viral oncogene homolog), bHLH (basic helix-loop-helix) and WRKY families were upregulated in leaves at 10 d of MgD stress, indicating that magnesium signaling transduction might be activated to compensate MgD. In addition, genes including chlorophyll(ide) b reductase (NYC1/NOL) chlorophyll/bacteriochlorophyll synthase (G4) and 7-hydroxymethyl chlorophyll a reductase synthesizing (HCAR) chlorophyll a and chlorophyll b were down-regulated in leaves, while those scavenging reactive oxygen species (ROS) were mainly up-regulated at 10 d of MgD stress. These results shed light on underlying MgD in N. cadamba.

Population Structure and Genetic Diversity in the Natural Distribution of Neolamarckia cadamba in China

Neolamarckia cadamba (Roxb.) Bosser is a fast-growing deciduous tree species and belongs to the Neolamarckia genus of the Rubiaceae family. This species has great economic and medical values in addition to being an important timber species for multiple industrial purposes. However, few studies have examined the genetic diversity and population structure in the natural distribution of this species in China. Here, we applied both the haploid nrDNA ITS (619 bp for aligned sequences) and mtDNA (2 polymorphic loci) markers to investigate 10 natural populations (239 individuals in total) that covered most of the distribution of the species in China. The results showed that the nucleotide diversity was π = 0.1185 ± 0.0242 for the nrDNA ITS markers and π = 0.00038 ± 0.00052 for the mtDNA markers. The haplotype diversity for the mtDNA markers was h = 0.1952 ± 0.2532. The population genetic differentiation was small (Fstn = 0.0294) for the nrDNA ITS markers but large (Fstm = 0.6765) for the mtDNA markers. There were no significant effects of isolation by distance (IBD), by elevation, and by two climatic factors (annual average precipitation and tem perature). A geographic structure among populations (Nst<Gst) was absent. Phylogenetic analysis showed a highly genetic mixture among individuals of the ten populations. Pollen flow was substantially greater than seed flow (mp/ms ≫ 1.0) and played a dominant role in shaping population genetic structure. The nrDNA ITS sequences were neutral and all local populations did not undergo demographic expansion. The overall results provide fundamental information for the genetic conservation and breeding of this miraculous tree.

Gene coexpression networks allow the discovery of two strictosidine synthases underlying monoterpene indole alkaloid biosynthesis in Uncaria rhynchophylla

Plant-derived monoterpene indole alkaloids (MIAs) from Uncaria rhynchophylla (UR) have huge medicinal properties in treating Alzheimer's disease, Parkinson's disease, and depression. Although many bioactive UR-MIA products have been isolated as drugs, their biosynthetic pathway remains largely unexplored. In this study, untargeted metabolome identified 79 MIA features in UR tissues (leaf, branch stem, hook stem, and stem), of which 30 MIAs were differentially accumulated among different tissues. Short time series expression analysis captured 58 pathway genes and 12 hub regulators responsible for UR-MIA biosynthesis and regulation, which were strong links with main UR-MIA features. Coexpression networks further pointed to two strictosidine synthases (UrSTR1/5) that were coregulated with multiple MIA-related genes and highly correlated with UR-MIA features (r > 0.7, P < 0.005). Both UrSTR1/5 catalyzed the formation of strictosidine with tryptamine and secologanin as substrates, highlighting the importance of key residues (UrSTR1: Glu309, Tyr155; UrSTR5: Glu295, Tyr141). Further, overexpression of UrSTR1/5 in UR hairy roots constitutively increased the biosynthesis of bioactive UR-MIAs (rhynchophylline, isorhynchophylline, corynoxeine, etc), whereas RNAi of UrSTR1/5 significantly decreased UR-MIA biosynthesis. Collectively, our work not only provides candidates for reconstituting the biosynthesis of bioactive UR-MIAs in heterologous hosts but also highlights a powerful strategy for mining natural product biosynthesis in medicinal plants.

Emerging mechanistic insights into the regulation of specialized metabolism in plants

Plants biosynthesize a broad range of natural products through specialized and species-specific metabolic pathways that are fuelled by core metabolism, together forming a metabolic network. Specialized metabolites have important roles in development and adaptation to external cues, and they also have invaluable pharmacological properties. A growing body of evidence has highlighted the impact of translational, transcriptional, epigenetic and chromatin-based regulation and evolution of specialized metabolism genes and metabolic networks. Here we review the forefront of this research field and extrapolate to medicinal plants that synthetize rare molecules. We also discuss how this new knowledge could help in improving strategies to produce useful plant-derived pharmaceuticals.

Evolution and Expression Patterns of the Fructose 1,6-Bisphosptase Gene Family in a Miracle Tree (Neolamarckia cadamba)

Neolamarckia cadamba (N. cadamba) is a fast-growing tree species with tremendous economic and ecological value; the study of the key genes regulating photosynthesis and sugar accumulation is very important for the breeding of N. cadamba. Fructose 1,6-bisphosptase (FBP) gene has been found to play a key role in plant photosynthesis, sugar accumulation and other growth processes. However, no systemic analysis of FBPs has been reported in N. cadamba. A total of six FBP genes were identifed and cloned based on the N. cadamba genome, and these FBP genes were sorted into four groups. The characteristics of the NcFBP gene family were analyzed such as phylogenetic relationships, gene structures, conserved motifs, and expression patterns. A cis-acting element related to circadian control was first found in the promoter region of FBP gene. Phylogenetic and quantitative real-time PCR analyses showed that NcFBP5 and NcFBP6 may be chloroplast type 1 FBP and cytoplasmic FBP, respectively. FBP proteins from N. cadamba and 22 other plant species were used for phylogenetic analyses, indicating that FBP family may have expanded during the evolution of algae to mosses and differentiated cpFBPase1 proteins in mosses. This work analyzes the internal relationship between the evolution and expression of the six NcFBPs, providing a scientific basis for the evolutionary pattern of plant FBPs, and promoting the functional studies of FBP genes.

The Vinca minor genome highlights conserved evolutionary traits in monoterpene indole alkaloid synthesis

Vinca minor, also known as the lesser periwinkle, is a well-known species from the Apocynaceae, native to central and southern Europe. This plant synthesizes monoterpene indole alkaloids, which are a class of specialized metabolites displaying a wide range of bioactive- and pharmacologically important properties. Within the almost 50 monoterpene indole alkaloids it produces, V. minor mainly accumulates vincamine, which is commercially used as a nootropic. Using a combination of Oxford Nanopore Technologies long read- and Illumina short-read sequencing, a 679,098 Mb V. minor genome was assembled into 296 scaffolds with an N50 scaffold length of 6 Mb, and encoding 29,624 genes. These genes were functionally annotated and used in a comparative genomic analysis to establish gene families and to investigate gene family expansion and contraction across the phylogenetic tree. Furthermore, homology-based monoterpene indole alkaloid gene predictions together with a metabolic analysis across 4 different V. minor tissue types guided the identification of candidate monoterpene indole alkaloid genes. These candidates were finally used to identify monoterpene indole alkaloid gene clusters, which combined with synteny analysis allowed for the discovery of a functionally validated vincadifformine-16-hydroxylase, reinforcing the potential of this dataset for monoterpene indole alkaloids gene discovery. It is expected that access to these resources will facilitate the elucidation of unknown monoterpene indole alkaloid biosynthetic routes with the potential of transferring these pathways to heterologous expression systems for large-scale monoterpene indole alkaloid production.

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.

A highly contiguous, scaffold-level nuclear genome assembly for the fever tree (Cinchona pubescens Vahl) as a novel resource for Rubiaceae research

The Andean fever tree (Cinchona L.; Rubiaceae) is a source of bioactive quinine alkaloids used to treat malaria. C. pubescens Vahl is a valuable cash crop within its native range in northwestern South America, however, genomic resources are lacking. Here we provide the first highly contiguous and annotated nuclear and plastid genome assemblies using Oxford Nanopore PromethION-derived long-read and Illumina short-read data. Our nuclear genome assembly comprises 603 scaffolds with a total length of 904 Mbp (∼82% of the full genome based on a genome size of 1.1 Gbp/1C). Using a combination of de novo and reference-based transcriptome assemblies we annotated 72,305 coding sequences comprising 83% of the BUSCO gene set and 4.6% fragmented sequences. Using additional plastid and nuclear datasets we place C. pubescens in the Gentianales order. This first genomic resource for C. pubescens opens new research avenues, including the analysis of alkaloid biosynthesis in the fever tree.