The chromosome-scale genome of Magnolia officinalis provides insight into the evolutionary position of magnoliids

Pan-evolutionary and regulatory genome architecture delineated by an integrated macro- and microsynteny approach

The forthcoming massive genome data generated by the Earth BioGenome Project will open up a new era of comparative genomics, for which genome synteny analysis provides an important framework. Profiling genome synteny represents an essential step in elucidating genome architecture, regulatory blocks/elements and their evolutionary history. Here we describe PanSyn, ( https://github.com/yhw320/PanSyn ), the most comprehensive and up-to-date genome synteny pipeline, providing step-by-step instructions and application examples to demonstrate its usage. PanSyn inherits both basic and advanced functions from existing popular tools, offering a user-friendly, highly customized approach for genome macrosynteny analysis and integrated pan-evolutionary and regulatory analysis of genome architecture, which are not yet available in public synteny software or tools. The advantages of PanSyn include: (i) advanced microsynteny analysis by functional profiling of microsynteny genes and associated regulatory elements; (ii) comprehensive macrosynteny analysis, including the inference of karyotype evolution from ancestors to extant species; and (iii) functional integration of microsynteny and macrosynteny for pan-evolutionary profiling of genome architecture and regulatory blocks, as well as integration with external functional genomics datasets from three- or four-dimensional genome and ENCODE projects. PanSyn requires basic knowledge of the Linux environment and Perl programming language and the ability to access a computer cluster, especially for large-scale genomic comparisons. Our protocol can be easily implemented by a competent graduate student or postdoc and takes several days to weeks to execute for dozens to hundreds of genomes. PanSyn provides yet the most comprehensive and powerful tool for integrated evolutionary and functional genomics.

A high-quality haplotype genome of Michelia alba DC reveals differences in methylation patterns and flower characteristics

Michelia alba DC is a highly valuable ornamental plant of the Magnoliaceae family. This evergreen tropical tree commonly grows in Southeast Asia and is adored for its delightful fragrance. Our study assembled the M. alba haplotype genome MC and MM by utilizing Nanopore ultralong reads, Pacbio Hifi long reads and parental second-generation data. Moreover, the first methylation map of Magnoliaceae was constructed based on the methylation site data obtained using Nanopore data. Metabolomic datasets were generated from the flowers of three different species to assess variations in pigment and volatile compound accumulation. Finally, transcriptome data were generated to link genomic, methylation, and morphological patterns to reveal the reasons underlying the differences between M. alba and its parental lines in petal color, flower shape, and fragrance. We found that the AP1 and AP2 genes are crucial in M. alba petal formation, while the 4CL, PAL, and C4H genes control petal color. The data generated in this study serve as a foundation for future physiological and biochemical research on M. alba, facilitate the targeted improvement of M. alba varieties, and offer a theoretical basis for molecular research on Michelia L.

The American Cherimoya Genome Reveals Insights into the Intra-Specific Divergence, the Evolution of Magnoliales, and a Putative Gene Cluster for Acetogenin Biosynthesis

Annona cherimola (cherimoya) is a species renowned for its delectable fruit and medicinal properties. In this study, we developed a chromosome-level genome assembly for the cherimoya 'Booth' cultivar from the United States. The genome assembly has a size of 794 Mb with a N50 = 97.59 Mb. The seven longest scaffolds account for 87.6% of the total genome length, which corresponds to the seven pseudo-chromosomes. A total of 45,272 protein-coding genes (≥30 aa) were predicted with 92.9% gene content completeness. No recent whole genome duplications were identified by an intra-genome collinearity analysis. Phylogenetic analysis supports that eudicots and magnoliids are more closely related to each other than to monocots. Moreover, the Magnoliales was found to be more closely related to the Laurales than the Piperales. Genome comparison revealed that the 'Booth' cultivar has 200 Mb less repeats than the Spanish cultivar 'Fino de Jete', despite their highly similar (>99%) genome sequence identity and collinearity. These two cultivars were diverged during the early Pleistocene (1.93 Mya), which suggests a different origin and domestication of the cherimoya. Terpene/terpenoid metabolism functions were found to be enriched in Magnoliales, while TNL (Toll/Interleukin-1-NBS-LRR) disease resistance gene has been lost in Magnoliales during evolution. We have also identified a gene cluster that is potentially responsible for the biosynthesis of acetogenins, a class of natural products found exclusively in Annonaceae. The cherimoya genome provides an invaluable resource for supporting characterization, conservation, and utilization of Annona genetic resources.

Identification and Validation of Magnolol Biosynthesis Genes in Magnolia officinalis

Bacterial infections pose a significant risk to human health. Magnolol, derived from Magnolia officinalis, exhibits potent antibacterial properties. Synthetic biology offers a promising approach to manufacture such natural compounds. However, the plant-based biosynthesis of magnolol remains obscure, and the lack of identification of critical genes hampers its synthetic production. In this study, we have proposed a one-step conversion of magnolol from chavicol using laccase. After leveraging 20 transcriptomes from diverse parts of M. officinalis, transcripts were assembled, enriching genome annotation. Upon integrating this dataset with current genomic information, we could identify 30 laccase enzymes. From two potential gene clusters associated with magnolol production, highly expressed genes were subjected to functional analysis. In vitro experiments confirmed MoLAC14 as a pivotal enzyme in magnolol synthesis. Improvements in the thermal stability of MoLAC14 were achieved through selective mutations, where E345P, G377P, H347F, E346C, and E346F notably enhanced stability. By conducting alanine scanning, the essential residues in MoLAC14 were identified, and the L532A mutation further boosted magnolol production to an unprecedented level of 148.83 mg/L. Our findings not only elucidated the key enzymes for chavicol to magnolol conversion, but also laid the groundwork for synthetic biology-driven magnolol production, thereby providing valuable insights into M. officinalis biology and comparative plant science.

The chromosome-scale genome of Magnolia sinica (Magnoliaceae) provides insights into the conservation of plant species with extremely small populations (PSESP)

Magnolia sinica (Magnoliaceae) is a highly threatened tree endemic to southeast Yunnan, China. In this study, we generated for the first time a high-quality chromosome-scale genome sequence from M. sinica, by combining Illumina and ONT data with Hi-C scaffolding methods. The final assembled genome size of M. sinica was 1.84 Gb, with a contig N50 of ca. 45 Mb and scaffold N50 of 92 Mb. Identified repeats constituted approximately 57% of the genome, and 43,473 protein-coding genes were predicted. Phylogenetic analysis shows that the magnolias form a sister clade with the eudicots and the order Ceratophyllales, while the monocots are sister to the other core angiosperms. In our study, a total of 21 individuals from the 5 remnant populations of M. sinica, as well as 22 specimens belonging to 8 related Magnoliaceae species, were resequenced. The results showed that M. sinica had higher genetic diversity (θw = 0.01126 and θπ = 0.01158) than other related species in the Magnoliaceae. However, population structure analysis suggested that the genetic differentiation among the 5 M. sinica populations was very low. Analyses of the demographic history of the species using different models consistently revealed that 2 bottleneck events occurred. The contemporary effective population size of M. sinica was estimated to be 10.9. The different patterns of genetic loads (inbreeding and numbers of deleterious mutations) suggested constructive strategies for the conservation of these 5 different populations of M. sinica. Overall, this high-quality genome will be a valuable genomic resource for conservation of M. sinica.

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.

Potential herb‒drug interactions between anti-COVID-19 drugs and traditional Chinese medicine

Coronavirus disease 2019 (COVID-19), caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), has spread worldwide. Effective treatments against COVID-19 remain urgently in need although vaccination significantly reduces the incidence, hospitalization, and mortality. At present, antiviral drugs including Nirmatrelvir/Ritonavir (Paxlovid^TM), Remdesivir, and Molnupiravir have been authorized to treat COVID-19 and become more globally available. On the other hand, traditional Chinese medicine (TCM) has been used for the treatment of epidemic diseases for a long history. Currently, various TCM formulae against COVID-19 such as Qingfei Paidu decoction, Xuanfei Baidu granule, Huashi Baidu granule, Jinhua Qinggan granule, Lianhua Qingwen capsule, and Xuebijing injection have been widely used in clinical practice in China, which may cause potential herb-drug interactions (HDIs) in patients under treatment with antiviral drugs and affect the efficacy and safety of medicines. However, information on potential HDIs between the above anti-COVID-19 drugs and TCM formulae is lacking, and thus this work seeks to summarize and highlight potential HDIs between antiviral drugs and TCM formulae against COVID-19, and especially pharmacokinetic HDIs mediated by metabolizing enzymes and/or transporters. These well-characterized HDIs could provide useful information on clinical concomitant medicine use to maximize clinical outcomes and minimize adverse and toxic effects.

The genome of Magnolia hypoleuca provides a new insight into cold tolerance and the evolutionary position of magnoliids

Magnolia hypoleuca Sieb. & Zucc, a member of the Magnoliaceae of magnoliids, is one of the most economically valuable, phylogenetic and ornamental tree species in Eastern China. Here, the 1.64 Gb chromosome-level assembly covers 96.64% of the genome which is anchored to 19 chromosomes, with a contig N50 value of 1.71 Mb and 33,873 protein-coding genes was predicted. Phylogenetic analyses between M. hypoleuca and other 10 representative angiosperms suggested that magnoliids were placed as a sister group to the eudicots, rather than sister to monocots or both monocots and eudicots. In addition, the relative timing of the whole-genome duplication (WGD) events about 115.32 Mya for magnoliid plants. M. hypoleuca was found to have a common ancestor with M. officinalis approximately 23.4 MYA, and the climate change of OMT (Oligocene-Miocene transition) is the main reason for the divergence of M. hypoleuca and M. officinalis, which was along with the division of Japanese islands. Moreover, the TPS gene expansion observed in M. hypoleuca might contribute to the enhancement of flower fragrance. Tandem and proximal duplicates of younger age that have been preserved have experienced more rapid sequence divergence and a more clustered distribution on chromosomes contributing to fragrance accumulation, especially phenylpropanoid, monoterpenes and sesquiterpenes and cold tolerance. The stronger selective pressure drived the evolution of tandem and proximal duplicates toward plant self-defense and adaptation. The reference M. hypoleuca genome will provide insights into the evolutionary process of M. hypoleuca and the relationships between the magnoliids with monocots and eudicots, and enable us to delve into the fragrance and cold tolerance produced by M. hypoleuca and provide more robust and deep insight of how the Magnoliales evolved and diversified.

Phylogenetic studies of magnoliids: Advances and perspectives

Magnoliids are the largest flowering plant clades outside of the eudicots and monocots, which are distributed worldwide and have high economic, ornamental and ecological values. Eudicots, monocots and magnoliids are the three major clades of Mesangiospermae, and their phylogenetic relationship is one of the most interesting issues. In recent years, with the continuous accumulation of genomic information, the evolutionary status of magnoliids has become a hot spot in plant phylogenetic research. Although great efforts have been made to study the evolution of magnoliids using molecular data from several representative species such as nuclear genome, plastid genome, mitochondrial genome, and transcriptome, the results of current studies on the phylogenetic status of magnoliids are inconsistent. Here, we systematically describe the current understanding of the molecular research on magnoliid phylogeny and review the differences in the evolutionary state of magnoliids. Understanding the research approaches and limitations of magnoliid phylogeny can guide research strategies to further improve the study of the phylogenetic evolution of magnoliids.

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

Genomic Data Reveals Population Genetic and Demographic History of Magnolia fistulosa (Magnoliaceae), a Plant Species With Extremely Small Populations in Yunnan Province, China

Elucidating the genetic background of threatened species is fundamental to their management and conservation, and investigating the demographic history of these species is helpful in the determination of the threats facing them. The woody species of the genus Magnolia (Magnoliaceae) have high economic, scientific and ecological values. Although nearly half of all Magnolia species have been evaluated as threatened, to date there has been no population genetic study employing Next Generation Sequencing (NGS) technology in this genus. In the present study, we investigate the conservation genomics of Magnolia fistulosa, a threatened species endemic to the limestone area along the Sino-Vietnamese border, using a double digest restriction-site-associated DNA-sequencing (ddRAD-seq) approach. To increase the reliability of our statistical inferences, we employed two approaches, Stacks and ipyrad, for SNP calling. A total of 15,272 and 18,960, respectively, putatively neutral SNPs were generated by Stacks and ipyrad. Relatively high genetic diversity and large population divergence were detected in M. fistulosa. Although higher absolute values were calculated using the ipyrad data set, the two data sets showed the same trends in genetic diversity (π, H e), population differentiation (F ST) and inbreeding coefficients (F IS). A change in the effective population size of M. fistulosa within the last 1 Ma was detected, including a population decline about 0.5-0.8 Ma ago, a bottleneck event about 0.2-0.3 Ma ago, population fluctuations during the last glacial stage, and the recovery of effective population size after the last glacial maximum. Our findings not only lay the foundation for the future conservation of this species, but also provide new insights into the evolutionary history of the genus Magnolia in southeastern Yunnan, China.

Establishing Tetraploid Embryogenic Cell Lines of Magnolia officinalis to Facilitate Tetraploid Plantlet Production and Phenotyping

The production of synthetic polyploids for plant breeding is compromised by high levels of mixoploids and low numbers of solid polyploid regenerants during in vitro induction. Somatic embryogenesis could potentially contribute to the maximization of solid polyploid production due to the single cell origin of regenerants. In the present study, a novel procedure for establishing homogeneous tetraploid embryogenic cell lines in Magnolia officinalis has been established. Embryogenic cell aggregate (ECA) about 100-200 μm across, and consisting of dozens of cells, regenerated into a single colony of new ECAs and somatic embryos following colchicine treatment. Histological analysis indicated that the few cells that survived some colchicine regimes still regenerated to form a colony. In some colonies, 100% tetraploid somatic embryos were obtained without mixoploid formation. New granular ECA from single colonies with 100% tetraploid somatic embryos were isolated and cultured individually to proliferate into cell lines. These cell lines were confirmed to be homogeneous tetraploid by flow cytometry. Many tetraploid somatic embryos and plantlets were differentiated from these cell lines and the stability of ploidy level through the somatic embryogenesis process was confirmed by flow cytometry and chromosome counting. The establishment of homogeneous polyploid cell lines, which were presumed to represent individual polyploidization events, might expand the phenotypic variations of the same duplicated genome and create novel breeding opportunities using newly generated polyploid plantlets.

The bZIP transcription factors in Liriodendron chinense: Genome-wide recognition, characteristics and cold stress response

The basic leucine zipper (bZIP) is a transcription factor family that plays critical roles in abiotic and biotic stress responses as well as plant development and growth. A comprehensive genome-wide study in Liriodendron chinense was conducted to identify 45 bZIP transcription factors (LchibZIPs), which were divided into 13 subgroups according the phylogenetic analysis. Proteins in the same subgroup shared similar gene structures and conserved domains, and a total of 20 conserved motifs were revealed in LchibZIP proteins. Gene localization analysis revealed that LchibZIP genes were unequally distributed across 16 chromosomes, and that 4 pairs of tandem and 9 segmental gene duplications existed. Concluding that segmental duplication events may be strongly associated with the amplification of the L. chinense bZIP gene family. We also assessed the collinearity of LchibZIPs between the Arabidopsis and Oryza and showed that the LchibZIP is evolutionarily closer to O. sativa as compared to the A. thaliana. The cis-regulatory element analysis showed that LchibZIPs clustered in one subfamily are involved in several functions. In addition, we gathered novel research suggestions for further exploration of the new roles of LchibZIPs from protein-protein interactions and gene ontology annotations of the LchibZIP proteins. Using the RNA-seq data and qRT-PCR we analyzed the gene expression patterns of LchibZIP genes, and showed that LchibZIP genes regulate cold stress, especially LchibZIP4 and LchibZIP7; and LchibZIP2 and LchibZIP28 which were up-regulated and down-regulated by cold stress, respectively. Studies of genetic engineering and gene function in L. chinense can benefit greatly from the thorough investigation and characterization of the L. chinense bZIP gene family.