Xue-Chan Tian1†
Jing-Fang Guo1†
Xue-Mei Yan1Tian-Le Shi1
Shuai Nie1Shi-Wei Zhao1Yu-Tao Bao1Zhi-Chao Li1Lei Kong1Guang-Ju Su2
Jian-Feng Mao1,3*Jinxing Lin1*- 1National Engineering Research Center of Tree Breeding and Ecological Restoration, State Key Laboratory of Tree Genetics and Breeding, Key Laboratory of Genetics and Breeding in Forest Trees and Ornamental Plants, Ministry of Education, College of Biological Sciences and Technology, Beijing Forestry University, Beijing, China
- 2National Tree Breeding Station for Nanmu in Zhuxi, Forest Farm of Zhuxi County, Hubei, China
- 3Department of Plant Physiology, Umeå Plant Science Centre, Umeå University, Umeå, Sweden
Wood decay resistance (WDR) is marking the value of wood utilization. Many trees of the Lauraceae have exceptional WDR, as evidenced by their use in ancient royal palace buildings in China. However, the genetics of WDR remain elusive. Here, through comparative genomics, we revealed the unique characteristics related to the high WDR in Lauraceae trees. We present a 1.27-Gb chromosome-level assembly for Lindera megaphylla (Lauraceae). Comparative genomics integrating major groups of angiosperm revealed Lauraceae species have extensively shared gene microsynteny associated with the biosynthesis of specialized metabolites such as isoquinoline alkaloids, flavonoid, lignins and terpenoid, which play significant roles in WDR. In Lauraceae genomes, tandem and proximal duplications (TD/PD) significantly expanded the coding space of key enzymes of biosynthesis pathways related to WDR, which may enhance the decay resistance of wood by increasing the accumulation of these compounds. Among Lauraceae species, genes of WDR-related biosynthesis pathways showed remarkable expansion by TD/PD and conveyed unique and conserved motifs in their promoter and protein sequences, suggesting conserved gene collinearity, gene expansion and gene regulation supporting the high WDR. Our study thus reveals genomic profiles related to biochemical transitions among major plant groups and the genomic basis of WDR in the Lauraceae.
Introduction
Wood is an exceptionally useful biomaterial, with myriad uses in construction, pulp and paper, and as a biofuel. Moreover, wood is a renewable material. One problem with using wood as a renewable biomaterial is that many microbes and insects have evolved to use wood as an energy source, producing enzymes that break down the components of the wood. Some species have evolved mechanisms to resist microbial damage and oxidation; many species with high wood decay resistance (WDR), such as teak (Tectona grandis), redwood (Sequoia sempervirens), and mahogany (Swietenia mahagoni) are rare and extremely valuable. Therefore, understanding the genetic basis and molecular mechanisms of WDR has the potential to provide effective information for improving WDR in commercially grown tree species. Wood is mainly composed of cellulose, hemicellulose, and lignin, which provide structural support for trees and resistance to microbial attack (Nascimento et al., 2013). Generally, lignin, a phenolic compound that is extremely resistant to degradation by certain fungi and plant diseases, acts as the basal component of wood durability by covering and protecting cellulose (Vance et al., 1980; Mounguengui et al., 2016). Further, trees resistant to decay exhibit significant production or accumulation of some bioactive compounds that function as antifungal compounds, antioxidants, or insect antifeedants, and are the main factors contributing to WDR (Nascimento et al., 2013). WDR is influenced by alkaloids such as indols and beta-carboline alkaloids, which have strong antifungal activity, as well as berberine and palmatine, which have shown good antifeedant and antioxidant activities (Kawaguchi et al., 1989; Anouhe et al., 2018; Ekeuku et al., 2020; Imenshahidi and Hosseinzadeh, 2020). Moreover, flavonoids are phenolic compounds with strong fungicidal activity, natural antioxidants and are excellent free radical scavengers, which have a significant effect on improving WDR (Schultz and Nicholas, 2000). In addition, terpenoids, including triterpenoids, diterpenoids, sesquiterpenoids, and monoterpenoids, have important antifungal, antifeeding, and antioxidant abilities, and contribute greatly to WDR (Park et al., 2000; Isman, 2002).
Lauraceae, a family of the order Laurales in the Magnoliids, includes about 67 genera and over 2,500 species (Anouhe et al., 2018). Lauraceae species are economically important, playing important roles in timber production, medicine, spice production, and ecological afforestation (Anouhe et al., 2018). A distinguishing feature of most Lauraceae species is the extremely high decay resistance of wood, including resistance to fungi, insect erosion, and oxidation (Jagels et al., 2005). Nanmu species, a group of tree species belonging to the Lauraceae family, are characterized by their straight trunks, fragrant and dense wood, and most notably by their super WDR (Jiao et al., 2022). Given these valuable traits, Nanmu wood is a precious natural resource that has historically been exploited, for example, for the construction of royal palaces (Xie et al., 2015). Generally, most species of the Phoebe and Machilu genera are recognized as Nanmu (e.g., Phoebe zhennan and Machilu nanmu) (Jiao et al., 2022). Another tree, Lindera megaphylla, has all superior qualities of the generally accepted Nanmu species, and was extensively used for the construction of royal buildings in Beijing in the Qing dynasty (Figure S1). L. megaphylla accumulates a variety of alkaloids (Chou et al., 1994) that promote resistance to microbial infection and herbivore attack, increasing the antifeeding and antioxidant activities of its wood (Kawaguchi et al., 1989; Ekeuku et al., 2020). L. megaphylla also has a wide range of medicinal properties due to alkaloid accumulation (Cao et al., 2016). In addition, the wood of some other Lauraceae species, e.g., Cinnamomum (Zhou et al., 2019) and Litsea species, have good natural durability and are highly valuable in construction, furniture, sculpture, and other building applications. With the development of society, there is increasing demand for naturally durable wood. However, genetic studies on the natural durability of wood, especially of Lauraceae species, are limited. Therefore, it is of great significance to identify the genes of biosynthetic pathways related to WDR, to investigate whether the WDR-related gene families have expanded significantly, and to reveal whether there are unique and conserved characteristics of WDR-related genes in Lauraleae species.
The phylogenetic location of Magnoliids remains to be further clarified. Lindera megaphylla belongs to Lauraceae, which together with Canellales, Piperales, and Magnoliales, constitutes the Magnoliids, including 9,000 species (The Angiosperm Phylogeny Group et al., 2016). Although multiple genomes of Magnoliids have been published, the relationship between magnoliids, eudicots, and monocots remains discordant. For example, the gene sequence-based phylogenomic analyses of Liriodendron chinense (Chen et al., 2019), Piper nigrum (Hu et al., 2019), Persea americana (Rendon-Anaya et al., 2019) and Phoebe bournei (Chen et al., 2020a) supported the Magnoliids as sister to the monocots-eudicots clade, while analyses of Cinnamomum kanehirae (Chaw et al., 2019a), Chimonanthus salicifolius (Lv et al., 2020) and Chimonanthus praecox (Shang et al., 2020a) supported Magnoliids as sister clade of eudicots. In addition, the phylogenomic analyses of Litsea cubeba suggested that the definite evolutionary relationships between Magnoliids, monocots, and eudicots remains to be resolved due to the possibility of incomplete lineage sorting (ILS) (Chen et al., 2020b). Microsynteny, gene colocality or collinearity, is the local conservation of gene order or gene neighborhood. Microsynteny provides valuable information to infer gene and genome evolution (Bowers et al., 2003; Van de Peer, 2004; Dewey, 2011), and is significant in phylogenetic inferences (Zhao and Schranz, 2019; Zhao et al., 2021c).
Here, we generated a chromosome-level genome assembly of L. megaphylla with long-read sequencing and Hi-C scaffolding technologies. The wood of L. megaphylla is dense and durable, making it an ideal material for construction, furniture, and shipbuilding. We conducted phylogenomic reconstruction of main angiosperm groups based on multiple strategies of concatenation, coalescent-based, and network-based microsynteny. Further, through the comparative genomics, especially shared gene microsynteny among major angiosperm lineages, we identified unique gene duplications and conserved microsynteny associated with isoquinoline alkaloids (IA), flavonoids, lignin, and terpenoids biosynthesis in Lauraceae species, which may be associated with outstanding wood durability in Lauraceae trees. The genome resources and findings presented here provide a basis for further evolutionary or functional studies in Lauraceae species, and for additional exploration of Lauraceae wood decay resistance.
Results
L. megaphylla genome sequencing, assembly, and gene annotation
As a first step to understand genomics of WDR in Lauraceae species with significant WDR, we sequenced the genome L. megaphylla. According to k-mer analysis, the genome size of L. megaphylla (Figure S1) was estimated to be ~1.3 Gb, with a 0.5% heterozygosity rate (Figure S2 and Note S1 for details). We generated 178.78 Gb (10.3 million reads, roughly 130× coverage) of Oxford Nanopore Technologies (ONT) long reads (Table S1) for primary assembly, 160.28 Gb (1068 million reads, 120× coverage, PCR-free library) of Illumina paired-end reads for correction and polishing, and 223.23 Gb (1488.194 million reads, 170× coverage) of Hi-C paired-end reads for scaffolding (Figure S3 and Table S1). A final genome assembly of 1.27 Gb was obtained, which consisted of 486 scaffolds, including 12 chromosome-level scaffolds, with a scaffold N50 of 104 Mb (Figure 1; Table 1; Tables S2, 3). The high confidence of the genome assembly was supported by high ten-fold minimum genome coverages of 95.1% (Illumina) and 99.6% (ONT), as well as the high mapping rates of 99.2% (Illumina) and 81.3% (ONT) reads. A 90.7% (1,306 complete genes) Benchmarking Universal Single Copy Orthologs (BUSCO) recovery score (Simão et al., 2015) and a high LTR Assembly Index (LAI) (Ou et al., 2018) score of 12.40 revealed a high completeness in the final assembly (Table S4).
Figure 1 Phylogenomic analysis of three major angiosperm groups. (A) Phylogenetic tree of 18 plant species generated by the concatenation-based method. Pie charts indicate the predicted expansion (red) and contraction (blue) of the gene family. The numbers represent divergence time of each node (Mya, million years ago), and values in brackets are 95% confidence intervals for the time of divergence. The yellow circle shows the WGD events identified in Lauraceae species. (B) Comparison of phylogenetic trees produced by the concatenation- and multi-species coalescent (MSC)-based methods. (C) Comparison of phylogenetic trees generated using concatenation- and microsynteny-based methods. (D) Comparison of phylogenetic trees produced using the microsynteny- and MSC-based methods.
Table 1 Statistics of the Lindera megaphylla genome assembly and annotation.
A total of 32,586 protein-coding genes were predicted from the final assembly (Table S5). The average lengths of total gene regions, transcripts, coding sequences, exons, and introns were 7,693.8, 1,410.1, 1,250.9, 270, and 1,094.7 bp, respectively (Table S5). In addition, we annotated 579 tRNAs, 248 rRNAs (including five 28S, six 18S, and 237 5S rRNAs), and 803 other non-coding RNAs (Table S6). The strongly supported gene annotation was evidenced by a 91.7% complete BUSCO score, as well as by 85.9% of the predicted genes (29,400 genes) with an annotation edit distance (AED) lower than 0.5 (Table S4 and Figure S4). More results of genome annotation are available in Note S4 and Table S7.
We identified 34,888 gene families, of which 6,340 are shared among all 18 species (Table S8) (see “Methods” section for details). And 885 expansion gene families in Lauraceae were enriched in isoquinoline alkaloid biosynthesis, flavonol biosynthesis, phenylpropanoid catabolism, lignin catabolic processes, and sulfur compound transport (Figure S5). All of these processes are tightly associated with resistance to bacteria and fungi, insect attacks, and high wood durability. The expanded gene families in L. megaphylla were also enriched in isoquinoline alkaloid biosynthesis, positive regulation of flavonoid biosynthesis, and isoflavone 7-O-glucosyltransferase activity (Figure S6). Similarly, these processes are all tightly associated with wood decay resistance.
Results of transposable element and other repeat annotation are available in Note S5, Figure S7 and Tables S9, S10.
Phylogenetic placement of Magnoliids
To determine the phylogenetic position of the Magnoliids relative to monocots and eudicots, phylogenetic trees were constructed using three distinct methods (concatenation-, coalescent-, and microsynteny-based approaches). For the concatenation-based approach, we constructed a phylogenetic tree using 885 low-copy orthologs from 18 species, with Amborella trichopoda and Nymphaea colorata as the outgroup (Figure 1A) (see “Methods” section). Results showed that the Maximum likelihood (ML) trees placed the Magnoliids as sister to the eudicots (Figure 1A). Phylogenetic analysis indicated that divergence time between Magnoliids and eudicots was 158.8 million years ago (Mya), with 95% confidence intervals of 143.8-174.5 Mya (Figure 1A), which overlaps with the C. kanehirae genome (136-209 Mya) (Chaw et al., 2019a). Lauraceae divergence was 124.16 Mya (Figure 1A), which was approximately equal to Phoebe bournei (Chen et al., 2020a). In addition, L. megaphylla diverged from C. kanehirae and L. cubeba around 10.52 Mya and 8.45 Mya, respectively (Figure 1A).
To reduce the influence of incomplete lineage sorting (ILS) on the determination of phylogenetic position, we also performed coalescent-based analyses of gene trees from the 855 low-copy gene families with ASTRAL-Pro (version 1.1.2) (Zhang et al., 2020a). The result from the coalescent-based analysis with strongly supported topology was highly consistent with the results of the concatenation-based method, placing Magnoliids as a sister group to eudicots after their divergence from monocots (Figure 1B). In addition, to reduce the interference caused by gene duplication and loss, ancestor hybridization, and lateral gene transfer in the homology assessment of plants, a novel method for phylogenetic tree reconstruction based on genome-wide synteny network data has been proposed (Zhao and Schranz, 2019; Zhao et al., 2021c). This method, microsynteny or gene order conservation, has been considered to be a valuable and alternative phylogenetic character in addition to sequence-based characters (Zhao et al., 2021c). The microsynteny-based analysis results confirmed that Magnoliids and eudicots are sister groups, which was topologically identical to the results of the above two methods (Figures 1C, D). These results strongly support that Magnoliids and eudicots are sister branches of monocots.
Microsynteny sharing and functional implications
To examine the lineage-specific microsynteny profile of major plant groups (Magnoliids, monocots, and eudicots), the genome synteny cluster obtained from microsynteny-based analysis of 16 species excluding N. colorata and A. trichopoda was analyzed. Interestingly, the Lauraceae species L. megaphylla, L. cubeba, and C. kanehirae had the largest number of microsyntenic clusters, with 15,347, 14,879, and 14,830 from each species, respectively (Figure S8A). The number of microsyntenic clusters shared by Magnoliids-eudicots (3,840) was significantly more than that shared by Magnoliids-monocots (871) and eudicots-monocots (491) (Figure S8B). Based on the heatmap of correlation in shared microsynteny, we observed a strong correlation between Magnoliids and eudicots (see “Methods” section) (Figure 2A). In contrast, the monocots showed a weak correlation with the other two clades, especially S. bicolor and O. sativa, which belong to the Poaceae (Figure 2A). These data signified a closer relationship between Magnoliids and eudicots.
Figure 2 Analysis of gene microsyntenic clusters. (A) Heatmap of the number of microsyntenic clusters shared among Magnoliid, eudicot, and monocot species. (B) Venn diagram showing the microsyntenic clusters shared among Magnoliids, eudicots, and monocots. (C) The black solid circle on the left surrounds microsyntenic clusters shared by eudicots and monocots, where blue dots represent eudicots, and red dots represent monocots. The black dotted rectangle on the right highlights an example of a cluster (in a subnetwork) shared between eudicots and monocots. (D) The black solid circle on the right surrounds microsyntenic clusters shared by Magnoliids and monocots, where yellow dots represent Magnoliids and red dots represent monocots. The black dotted rectangle on the left highlights an example of a cluster (in a subnetwork) shared between Magnoliids and monocots. (E) The black solid circle on the left surrounds microsyntenic clusters shared by Magnoliids and eudicots, where yellow dots represent Magnoliids and blue dots represent eudicots. The black dotted rectangle on the right highlights an example of a cluster (in a subnetwork) shared between Magnoliids and eudicots.
Next, we examined the functional implications of the shared or group-specific microsyntenic clusters among the three clades by removing the species-specific cluster (see “Methods” section). We discovered 2,839, 1,758, and 1,208 clusters specific to Magnoliids, eudicots, and monocots, respectively (Figure 2B). The number of synteny clusters common to Magnoliids-eudicots was still the largest (358), followed by Magnoliids-monocots (54), and eudicots-monocots (37) (Figures 2B–E, S9). As revealed in the UpSet plot, the Poaceae species Sorghum bicolor and Oryza sativa shared the largest number of clusters (6,283), followed by Piper nigrum and Musa acuminata with 3,243 and 1,651 species-specific clusters, respectively (Figure S8C). Four Lauraceae species also shared many clusters (1,460) (Figure S8C). Excluding these species-specific and clade-specific clusters, the six Magnoliids and six eudicot species shared the most clusters (39) (Figure S8C). These results further supported that Magnoliids and eudicots may be most closely related.
Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Ontology (GO) enrichment analyses of eudicot-specific microsyntenic clusters showed that they were mainly associated with terms related to a series of signaling pathways (Figure S10). Clusters specific to the Magnoliids were mainly enriched in terms such as isoquinoline alkaloid biosynthesis, ribosome biogenesis, brassinosteroid biosynthetic process, phospholipid biosynthetic process, and secondary metabolite biosynthesis (Figure S10). The microsyntenic clusters in monocots were mainly enriched in terms such as histidine metabolism, chloroalkane limonene and pinene degradation, cell plate assembly, and pyrimidine metabolism (Figure S10). Remarkably, the synteny clusters shared by Magnoliids, eudicots, and monocots were significantly enriched in sesquiterpenoid, diterpenoid, and triterpenoid biosynthesis (Figure S10). This finding indicates that the genes involved in terpenoid biosynthesis are conserved among plant clades, indicating the importance of terpenoids in various plants. In addition, the unique clusters in Lauraceae were mainly enriched in isoquinoline alkaloid biosynthesis, phenylpropanoid metabolic process, secondary metabolic process and lignin metabolic process, revealing potential links to the super WDR of Lauraceae trees (Figure S11).
Inference of whole-genome duplication in Lauraceae species are available in Note S7 and Figures S12-S14.
Tandem duplicate/proximal duplicate gene duplications in Lauraceae
A total of 28,838, 22,618 and 25,951 duplicated genes originating from whole-genome duplicates (WGD), tandem duplicates (TD), proximal duplicates (PD), dispersed duplicates (DSD), and transposed duplicates (TRD) were annotated in L. megaphylla, C. kanehirae, and L. cubeba, respectively (Figure 3A and Table S11). Aside from 18.39% TD/PD genes in Aquilegia coerulea, high TD/PD ratios were found in the Lauraceae species L. megaphylla (19.97%), L. cubeba (18.36%), C. kanehirae (22.45%), and P. bournei (18.44%) (Figure 3A and Table S12). Gene families expanded via TD/PD duplications in Lauraceae were functionally enriched in GO categories significantly associated with wood decay resistance, such as lignin catabolism, isoquinoline alkaloid biosynthesis, flavonol biosynthesis, and phenylpropanoid catabolism (Figure 3B). KEGG enrichment confirmed this pattern, showing that TD/PD duplications were enriched in isoquinoline alkaloid biosynthesis, flavone and flavonol biosynthesis, phenylpropanoid and flavonoid biosynthesis, monoterpenoid biosynthesis, antibiotic biosynthesis, defense response to bacterium, response to oxidative stress, cyanoamino acid metabolism, tropane, piperidine and pyridine alkaloid biosynthesis, and sulfur metabolism (Figure 3B). In addition to these functions, KEGG and GO analyses also revealed significant enrichment of TD/PD duplications in the biosynthesis of various terpenoids in L. megaphylla, including diterpenoid, monoterpenoid, sesquiterpenoid, and triterpenoid biosynthesis (Figure S15) In summary, these results indicate that local gene duplication in Lauraceae contributed to the expansion of secondary metabolite biosynthesis genes related to WDR.
Figure 3 The expansion of duplicated genes. (A) The stacked bar chart shows the proportion of genes derived from five duplication types (WGD whole-genome duplication, TD tandem duplication, PD proximal duplication, TRD transposed duplication and DSD dispersed duplication). (B) GO and KEGG functional enrichment analysis of expanded genes arising from tandem and proximal duplicates (TD/PD) in Lauraceae. The red line represents GO enrichment and the blue line represents KEGG enrichment. Blue letters indicate terms related to wood decay resistance.
Genes involved in benzylisoquinoline alkaloid biosynthesis
Three different benzylisoquinoline alkaloid (BIA) biosynthesis pathways were annotated in four Lauraceae species (L. megaphylla, L. cubeba, C. kanehirae, and P. bournei), including magnoflorine, berberine, and palmatine biosynthesis pathways (Figure 4A and Table S13), all of which were important for improving decay resistance of wood. The termite antifeeding activities of berberine and palmatine have been well demonstrated (Kawaguchi et al., 1989; Park et al., 2000). A total of twelve gene families related to BIA biosynthesis were identified. The enzymes 4OMT, 6OMT, SOMT, and CoOMT belong to the O-methyltransferase (OMT) family, and CYP80G, CYP80B, and CYP719A belong to the cytochrome P450 (CYP) family. These enzymes are mainly found in Magnoliids and A. coerulea, but rarely in monocots and other core eudicots (Figure 4B). In addition, TD/PD duplication contributed the most BIA biosynthesis genes in Lauraceae, especially the 4OMT, 6OMT, CoOMT, CYP80G, CYP80B, CYP719A, and CNMT ((S)-coclaurine-N-methyltransferase) genes (Figures 4C, S16, and Table S14).
Figure 4 Characteristics of benzylisoquinoline alkaloid genes in Lauraceae. (A) The biosynthesis pathway of isoquinoline alkaloids. TyrAT, tyrosine aminotransferase; PDC, 4-hydroxyphenylpyruvate decarboxylase; NCS, (S)-norcoclaurine synthase; 6OMT, (RS)-norcoclaurine 6-O-methyltransferase; CNMT, (S)-coclaurine-N-methyltransferase; CYP80B, N-methylcoclaurine 3′-hydroxylase; 4OMT, 3′-hydroxy-N-methyl-(S)-coclaurine 4′-O-methyltransferase; CYP80G, (S)-corytuberine synthase; RNMT, reticuline N-methyltransferase; BBE, berberine bridge enzyme; SMT, (S)-scoulerine 9-O-methyltransferase; CYP719A, (S)-canadine synthase; THBO, tetrahydroberberine oxidase; CoOMT, columbamine O-methyltransferase. (B) Number of annotated genes in each enzyme gene family (4OMT, 6OMT, SOMT, CoOMT, CYP80G, CYP80B, CYP719A, BBE, NMT, NCS, and TyrAT) for each species. (C) Proportion of tandem (TD) and proximal (PD) duplication genes in each enzyme gene family (4OMT, 6OMT, SOMT, CoOMT, CYP80G, CYP80B, CYP719A, BBE, NMT, NCS, and TyrAT) for each species. (D) Microsyntenic gene clusters associated with subfamilies of the OMT gene family (here, 4OMT, 6OMT, and CoOMT). Circles in dashed red line denote the syntenic clusters (here, C4, C5, C7, C10 and C12) unique to Lauraceae species. (E) Heatmap of 12 microsyntenic clusters in (D), five of which are Lauraceae-specific and highlighted by a red square. The color in the heatmap represents the gene number in each cluster for each species. (F) Phylogenetic analysis of OMT gene families. The red stars represent genes within Lauraceae-specific gene clusters identified in (D).
A total of 12 microsyntenic clusters were identified as related to OMT gene families (here, 4OMT, 6OMT, and CoOMT) (Figure 4D). Five of these 12 microsyntenic clusters were specific to Lauraceae, including C5, C10 and C12 associated with 6OMT, C7 associated with 4OMT, and C4 with CoOMT (Figures 4D–F). These genes on Lauraceae-specific microsyntenic clusters may play an important role in the unique WDR of Lauraceae species.
6OMT is involved in the rate-limiting step of isoquinoline biosynthesis (Robin et al., 2016). Phylogenetic analysis showed that the 6OMT genes in Lauraceae could be divided into five groups (Figure 5A). Genes of the Lauraceae-specific clusters C5, C10 and C12 were located in groups 2,4 and 5 respectively (Figures 5A–C). Protein sequence analysis found a Lauraceae-specific motif (motif 9) among genes in group 2 (C5). Genes in group 4 (C10) and 5 (C12) shared another Lauraceae-specific motif (motif 12) (Figure 5A). In addition, we identified several conserved motifs unique to Lauraceae through sequence analysis of gene promoters. Motif 1 existed in both C5 and C10 genes and overlapped with the predicted binding sites of bHLH transcription factors (TFs) (Figure 5A). Motif 8 was specific to C5 genes and overlapped with the predicted binding sites of ERF TFs (Figure 5A). Interestingly, six motifs (motif 5, 4, 2, 1, 3, and 6) formed a tandem cluster unique to Lauraceae genes in group 4 (C10). These motifs were the predicted binding sites of GATA, B3, ERF, bHLH, Trihelix, and MYB transcription factors (Figure 5A).
Figure 5 Specificity of 6OMT and 4OMT genes in Lauraceae. (A) The different panels illustrate the phylogenetic tree of the 6OMT gene family (left), the distribution of motifs in the promoter sequences and the predicted transcription factor binding sites (TFBS) (middle), and the distribution of motifs in protein sequences (right). Thick squares represent motifs and thin ones represent TFBSs. Dashed boxes highlight genes and promoter motifs unique to Lauraceae species. (B) The syntenic block containing the 6OMT gene family within the Lauraceae-specific microsynteny gene cluster (C5), which was identified in Figures 4D, E. This syntenic block was compared among P. bournei, C. kanehirae, L. cubeba, and L. megaphylla. Chartreuse squares represent the 6OMT genes and aquamarine squares represent other genes on the syntenic block. (C) The syntenic block containing the 6OMT gene family within the Lauraceae-specific microsynteny gene cluster (C10), which was identified in Figures 4D, E. This syntenic block was compared among P. bournei, C. kanehirae, L. cubeba, and L. megaphylla. Blue squares represent 6OMT genes and aquamarine squares represent other genes on the syntenic block. (D) The syntenic block containing the 4OMT gene family within the Lauraceae-specific microsynteny gene cluster (C7), which was identified in Figures 4D, E. This syntenic block was compared among P. bournei, C. kanehirae, L. cubeba, and L. megaphylla. Red squares represent 4OMT genes and aquamarine squares represent other genes on the syntenic block. (E) The different panels show the phylogenetic tree of the 4OMT gene family (left),the distribution of motifs in the promoter sequences and the predicted transcription factor binding sites (TFBS) (middle), and the distribution of motifs in protein sequences (right). Thick squares represent motifs and thin ones represent TFBSs. Dashed boxes highlight genes and promoter motifs unique to Lauraceae species.
In addition to 6OMT, 4OMT is also an important rate-limiting enzyme in BIA biosynthesis (Inui et al., 2012). The 4OMT genes in Lauraceae were divided into two groups, with C7 genes located in group 2 (Figures 5D, E). Four Lauraceae-specific and conserved motifs (motifs 10-14) were identified among the protein sequences of these group 2 (C7) genes (Figure 5E). Although no Lauraceae-specific microsyntenic cluster in group 1, phylogenetic analysis results showed that they were located in Lauraceae-specific clades, and two Lauraceae-specific motifs (motif 9 and motif 10) were identified (Figure 5E). Sequence analysis of gene promoters revealed that both groups (group1 and group 2) shared a common Lauraceae-specific DNA motif (motif 1), but only group 1 contained potential MYB transcription factor binding sites (TFBSs) (Figure 5D). One Lauraceae-specific promoter motif (motif 6) in group 2 (C7) overlapped with predicted WRKY and bHLH TFBSs (Figure 5E). Interestingly, we identified a unique promoter motif (motif 8) among group 1 genes, and found a conserved cluster formed by six motifs (motif 8, 9, 3, 5, 2, and 1) that overlap with WRKY, bHLH, B3, ERF, MYB, and C2H2 TFBSs. These TFBS clusters may play key roles in coordinating specific gene expression as well as efficient activation and regulation of alkaloid biosynthesis.
Columbamine O-methyltransferase (CoOMT) is a vital enzyme that catalyzes the formation of tetrahydropalmatine, an isoquinoline alkaloid. The C4, a Lauraceae-specific microsyntenic cluster contained all CoOMT genes (Figures 4E, S17A). We found that TD/PD duplications occurred before Lauraceae speciation, producing three major CoOMT groups (group 1, 2, and 3) (Figures S17A, 17B). In L. megaphylla, all members of the CoOMT family were found in one TD/PD cluster on chromosome 3 (Figure S17B). Two Lauraceae-specific motifs (motif 9 and motif 10) among CoOMT protein sequences were identified (Figure S17A). In addition, we identified two Lauraceae-specific promoter motifs (motif 3 and motif 4), of which motif 3 is the potential TFBS of WRKY, ERF, and MYB TFs, and motif 4 is the potential TFBS of bHLH TFs (Figures S17A, 17C).
Cytochrome P450 monooxygenases (CYPs) play an important role in the structural and functional diversity of alkaloids. The CYP80B, CYP80G, and CYP719A gene families play key oxidative roles in BIA metabolism (Hagel and Facchini, 2013; Nguyen and Dang, 2021). A total of 20 microsyntenic clusters were identified as related to the CYP gene family, among which three clusters were unique to Lauraceae (C11, C12, and C5) (Figures S18A, 18B, S19). Specifically, microsyntenic cluster C11 is related to the CYP719A family, and C12 and C5 are related to the CYP80G family. TD/PD expansion of genes on C5 cluster occurred in all Lauraceae species, especially in L. megaphylla (Figure S19A).
CYP719A catalyzes the conversion of (S)-tetrahydrocolumbamine to (S)-tetrahydroberberine, and is an essential enzyme in berberine biosynthesis (Ikezawa et al., 2003). According to the phylogenetic tree, CYP719A genes from Lauraceae can be divided into two groups. All members of C11 were classified into group 2, and these genes are located in a species-specific TD/PD cluster found on L. megaphylla chromosome 8 (Figure S20B). A motif unique to Lauraceae (motif 12) was discovered in the protein sequences of these CYP719A genes (Figures S20A, 20C). Further, three Lauraceae-specific motifs (motif 1, motif 4, and motif 8) were found in the promoters (Figures S20A, 20B). Among these motifs, motif 1 contains NAC TFBSs, motif 4 contains bHLH and ERF TFBSs, and motif 8 contains MYB TFBSs (Figure S20C). Notably, TFs such as bHLH, NAC, WRKY, and MYB have been implicated in the regulation of BIA biosynthesis in plants (Yamada et al., 2011; Zhou and Memelink, 2016; Deng et al., 2018). Here, we identified Lauraceae-specific and conserved protein sequences, TFBS motifs, and TFBS clusters among BIA biosynthesis genes. It is found that the genes related to BIA biosynthesis in Lauraceae species are significantly different from those in other species. These findings are valuable in the genetic dissection of BIA biosynthesis in Lauraceae species.
Characterization of genes involved in phenolic compound biosynthesis
We next examined the lignin and flavonoid biosynthesis pathways, which are the downstream branches of phenylpropanoid metabolism related to phenol biosynthesis (Figure 6A and Table S15). Phenolic compounds can protect wood from decaying organisms and improve WDR. The key reactions of general phenylpropanoid biosynthesis involve three enzymes: phenylalanine ammonia-lyase (PAL), cinnamate 4-hydroxylase (C4H), and 4-coumarate coenzyme A ligase (4CL). Among these enzymes, we found that the C4H and 4CL genes underwent remarkable TD/PD duplication events in Lauraceae (Figures 6A, B, S21, and Table S16).
Figure 6 Characteristics of flavonoid and lignin genes in Lauraceae. (A) Biosynthesis pathways of general phenylpropanoids, flavonoids, and lignin. PAL, phenylalanine ammonia-lyase; C4H, cinnamate-4-hydroxylase; 4CL, 4-coumarate CoA ligase 4; CHS, chalcone synthase; CHI, chalcone isomerase; F3H, flavanone 3-hydroxylase; FLS, flavonol synthase; F3′H, flavonoid 3′-hydroxylase; F3′5′H, flavonoid 3′,5′-hydroxylase; DFR, dihydroflavonol 4-reductase; ANS, anthocyanidin synthase; C3′H, p-coumaroyl shikimate 3′-hydroxylase; CCR, cinnamoyl-CoA reductase; CAD, (hydroxy)cinnamyl alcohol dehydrogenase; HCT, hydroxycinnamoyl-CoA:shikimate/quinate hydroxycinnamoyltransferase; CCoAOMT, caffeoyl-CoAO methyltransferase; F5H, coniferaldehyde/ferulate 5-hydroxylase; COMT, caffeicacid/5-hydroxyferulic acid O-methyltransferase. (B) Proportion of tandem and proximal duplication genes in 4CL, C4H, F3H, F3′5′H, C3′H, CAD, CCoAOMT, COMT, F5H and HCT gene families in each species. (C) Heatmap of 13 microsyntenic clusters of C4H gene families in 18 species. The Lauraceae-specific cluster is highlighted by a red square. Colors in the heatmap indicate gene number in each cluster for each species. (D) Phylogenetic analysis of C4H gene families. The gene names of Lauraceae species are shown in red, and red stars represent Lauraceae-specific gene clusters identified in (C). Yellow stars show the tandem and proximal duplication (TD/PD) genes. (E) The phylogenetic tree of the C4H gene family (left), the distribution of motifs in the promoter sequences and the predicted transcription factor binding sites (TFBS) (middle), and the distribution of motifs in protein sequences (right) are shown. Thick squares represent motifs and thin ones represent TFBSs. Dashed boxes highlight genes and promoter motifs unique to Lauraceae species. (F) The syntenic block containing the C4H gene family within the Lauraceae-specific microsynteny gene cluster (C4) identified in (C). Here, this syntenic block was compared among L. chinense, P. bournei, C. kanehirae, L. cubeba, and L. megaphylla. Chartreuse squares represent C4H genes and blue squares represent other genes on the syntenic block.
Our microsynteny analysis of C4H genes revealed a Lauraceae-specific cluster (C4) (Figure 6C). Genes of this C4 cluster were divided into two groups resulted from the Lauraceae-specific TD/PD duplication (Figure 6D). Two Lauraceae-specific protein motifs (motif 12 and motif 13) were identified in these C4 cluster genes (Figures 6E, F, S22). We also found two motifs (motif 2 and motif 4) specific to Lauraceae in the promoter regions of these genes (Figures 6E, S22). Motif 2 overlapped with C2H2 TFBSs and motif 8 with that of MYB and ERF TFs (Figures 6E, S22). Members of all these TF families are involved in the regulation of phenylpropanoid biosynthesis (Ma et al., 2017; Mondal and Roy, 2018; Teng et al., 2018). In addition, among the C4 genes, these two Lauraceae-specific motifs were clustered together with motifs 7, 9, 6, 3, 1, and 8, forming a very distinct cluster of ERF, MYB, bHLH, and ERF TFBSs. This motif cluster was shared among Lauraceae species and L. chinense (Figures 6E, S22).
Sequence analysis of promoter regions revealed two motifs (motif 8 and motif 9) unique to Lauraceae of PAL genes (Figure 7A). Motif 8 overlapped with TCP TFBSs and motif 9 with that of TCP and GATA TFs (Figure 7A). TCP TFs play an important role in plant defense and have been found to enhance flavonoid biosynthesis of Arabidopsis thaliana (Li and Zachgo, 2013; Li, 2014). Moreover, overexpression of a GATA gene can enhance the activity of the phenylpropanoid biosynthesis pathway in Solanum lycopersicum (Zhao et al., 2021b). Similar to the PAL genes, although there was no Lauraceae-specific collinearity cluster found related to 4CL genes, a motif in the promoters unique to Lauraceae (motif 7) was identified and overlapped with C2H2 TFBSs (Figure 7B). Moreover, conserved TFBS clusters were also found among the promoters of PAL and 4CL genes. These TFBSs were of TFs belonging to the TCP, BFR-BPC, C2H2, ERF, MYB, GATA, and GRAS families (Figures 7A, B).
Figure 7 Characterization of PAL and 4CL genes in Lauraceae species. (A) Different panels represent of the phylogenetic tree of the PAL gene family (the left), the distribution of motifs in the promoter sequence and the predicated transcription factor binding sites (TFBS) (the right). Fat squares represent the motifs and thin ones the TFBSs. Red boxes highlight the promoter motifs uniquely found among the Lauraceae species. (B) Different panels represent of the phylogenetic tree of the 4CL gene family (the left), the distribution of motifs in the promoter sequence and the predicated transcription factor binding sites (TFBS) (the right). Fat squares represent the motifs and thin ones the TFBSs. Purple boxes highlight the promoter motifs uniquely found among the Lauraceae species.
The biosynthesis pathways of taxifolin, myricetin, catechin, quercetin, and kaempferol have been annotated in Lauraceae species (Figure 6A). All of these flavonoids have been reported to improve plant WDR (Nascimento et al., 2013). TD/PD duplications accounted for expansions of F3H (flavanone 3-hydroxylase) and F3′5′H (flavonoid 3′,5′-hydroxylase) genes in Lauraceae (Figures 6B, S21, and Table S17). F3H is an important rate-limiting enzyme in flavonoid biosynthesis pathway. Enzymatic gene families of the lignin biosynthesis pathway include C3′H, HCT, CCR, CAD, CCoAOMT, F5H, and COMT, all of which were expanded through TD/PD duplications in Lauraceae (Figures 6B, S23, and Table S18). Microsynteny analysis revealed two Lauraceae-specific conserved gene clusters (C9 and C24) associated with lignin pathway genes (HCT and CCR) (Figures 8A–D). Although no Lauraceae-specific motifs and TFBSs were found among genes of the C9 and C24 clusters, all of these genes showed obvious TD/PD expansion (Figures 8B, D).
Figure 8 Lauraceae-specific CCR and HCT genes in lignin biosynthesis pathway. (A) Heatmap of 26 microsynteny clusters identified to be related with CCR gene family, one of which specific to Lauraceae were highlighted in a red square. Color in the heatmap was determined by the gene number found in each cluster for each species. (B) The syntenic block containing of CCR gene family within the Lauraceae-specific microsynteny gene cluster (C9) identified in (A). Here this syntenic block was compared among P. bournei, C. kanehirae, L. cubeba and L. megaphylla. Yellow squares represent the CCR genes and blue ones represent other genes on the syntenic block. (C) Heatmap of 27 microsynteny clusters identified to be related with HCT gene familiy, one of which specific to Lauraceae were highlighted in a red square. Color in the heatmap was determined by the gene number found in each cluster for each species. (D) The syntenic block containing of HCT gene family inside the Lauraceae-specific microsynteny gene cluster (C24) in (C). Here this syntenic block was compared among P. bournei, C. kanehirae, L. cubeba and L. megaphylla. Red squares represents the HCT genes and blue ones represents other genes on the syntenic block.
Remarkable TD/PD duplications were also found for TPS gene family of Lauraceae species, which may be associated with the super WDR. Details are available in Note S8, Figures S24-26 and Table S19.
Discussion
Our genomic investigation, especially the gene microsynteny profiling, may contribute to resolving the phylogenetic position of Magnoliids relative to eudicots and monocots, the other two major angiosperm groups. Although multiple assemblies of magnoliid genomes have been published, such as C. kanehirae (Chaw et al., 2019a), L. chinense (Chen et al., 2019), P. nigrum (Hu et al., 2019), P. americana (Rendon-Anaya et al., 2019), P. bournei (Chen et al., 2020a), L. cubeba (Chen et al., 2020b), C. salicifolius (Lv et al., 2020), and C. praecox (Shang et al., 2020a), the phylogenetic placement of Magnoliids still remains unclear. Our phylogenetic analyses using three different methods (concatenation-, coalescent-, and microsynteny-based approaches) confirmed that Magnoliids are the sister group of eudicots, which is in line with previous genomic analyses (Chaw et al., 2019a; Lv et al., 2020; Shang et al., 2020a) and phylotranscriptomic analyses of 92 streptophytes (Wickett et al., 2014) and 20 representative angiosperms (Zeng et al., 2014). In addition, the microsyntenic clusters of 16 species in Magnoliids, eudicots, and monocots were further analyzed. There were significantly more shared clusters in Magnoliids-eudicots compared with Magnoliids-monocots and eudicots-monocots, which strongly supports the finding that Magnoliids and eudicots are sister groups. The three clades were enriched in different GO and KEGG terms, indicating their functional divergence. The genes of Lauraceae-specific microsyntenic clusters were significantly enriched in terms including isoquinoline alkaloid biosynthesis, phenylpropanoid metabolic and lignin metabolic processes, suggesting that various Lauraceae-specific biochemical processes may influence its wood decay resistance.
A variety of bioactive compounds, including terpenoids, alkaloids, and phenolic compounds such as flavonoids, have been associated with WDR (Nascimento et al., 2013; Anouhe et al., 2018). In addition to the dual fungicidal and antioxidant effects of bioactive compounds, other factors such as lignin content also impact WDR (Vance et al., 1980; Nascimento et al., 2013; Mounguengui et al., 2016). Apart from annotating enzymes involved in the biosynthesis of isoquinoline alkaloids, flavonoids, terpenoids, and lignins, we characterized genes, gene syntenies, gene expansions, and gene promoter motifs specific to Lauraceae, which help to track genomic characters potentially related with the super wood decay resistance.
The biosynthetic pathways of three benzylisoquinoline alkaloids (BIA), namely magnoflorine, berberine, and palmatine, were annotated in Lauraceae species. Both berberine and palmatine exhibit significant antifeedant activity against termites (Kawaguchi et al., 1989). Magnoflorine is an aporphine-type BIA that has antibacterial and insecticidal effects, and may also play a role in improving WDR (Okon et al., 2020). Our comparative analyses demonstrated that the OMT, CYP, and BBE gene families involved BIA biosynthesis showed specific expansion in Lauraceae. Most members of these gene families originated from TD/PD duplications, which greatly enriched the enzymatic genes of the BIA biosynthesis pathway. These data indicate the significant value of TD/PD duplications in BIA biosynthesis. In the OMT gene family, a total of four Lauraceae-specific microsyntenic clusters were identified, including genes of the 4OMT, 6OMT, and CoOMT subfamilies. Again, TD and PD duplications were associated with significant expansion of the CoOMT gene family in L. megaphylla, which may have contributed to the accumulation of palmatine, thereby further improving WDR.
In addition to the Lauraceae-specific gene microsyntenic clusters uncovered for the biosynthesis of bioactive compounds related to WDR, we also found conserved TFBS clusters in the promoter regions of genes in these conserved clusters. These conserved TFBS clusters suggest conserved transcriptional regulation of secondary metabolite biosynthesis efficiency, which may lead to the high WDR trait shared among many Lauraceae woods. In the OMT gene family, the Lauraceae-specific promoter motifs were mainly TFBSs of bHLH, MYB, ERF and WRKY TFs. In the CYP gene family, the conserved promoter motifs were generally TFBSs of bHLH, MYB, ERF, and NAC TFs, all of which have been reported to be involved in the regulation of BIA biosynthesis (Yamada et al., 2011; Deng et al., 2018). In addition, we found that B3, GATA, Trihelix, and C2H2 TFs may bind these Lauraceae-specific TFBS clusters. However, their involvement in the regulation of alkaloid biosynthesis requires further evaluation. Compared with other species, the unique characteristics of Lauraceae species in BIA biosynthesis suggest that isoquinoline alkaloids may play a large proportion of roles in the decay resistance of Lauraceae.
There are diverse metabolic branches downstream of the general phenylpropanoid biosynthesis. Of these branches, we investigated the lignin and flavonoid pathways in the present study. The C4H and 4CL genes of the general phenylpropanoid pathway, F3H and F3′5′H of the flavonoid pathway, and all gene families of the lignin pathway have undergone significant TD/PD duplication in Lauraceae. C4H is the second key enzyme in the general phenylpropanoid biosynthesis pathway, and belongs to the CYP73A subfamily. C4H directly affects the biosynthesis and yield of flavonoids and lignin in plants (Ryan et al., 2002; Millar et al., 2007). Lauraceae-specific genes were found in the C4H gene family. In addition to carrying motifs in the coding and promoter regions that were different from other species, these C4H genes also had unique TFBS clusters specific to Lauraceae. Such TFBSs in the clusters are adjacent to each other, including binding sites of bHLH, C2H2, ERF, and MYB TFs, which all have important regulatory functions in phenylpropanoid biosynthesis (Ma and Constabel, 2019; Yadav et al., 2020; Meng et al., 2021). Moreover, TD/PD events also occurred in the Lauraceae-specific genes of the C4H gene family, which greatly increased their coding space, and further contributed to the WDR of Lauraceae species. In addition, Lauraceae-specific TFBS clusters were also found in the promoter regions of genes encoding PAL and 4CL. PAL is a rate-limiting enzyme that catalyzes the first step in the phyenylpropanoid biosynthesis pathway. Thus, it plays an important role in phenylpropanoid biosynthesis (Zhao et al., 2021a). 4CL, the third enzyme in the general phenylpropanoid biosynthesis pathway, participates in monolignol biosynthesis through the production of p-coumaroyl-CoA, a precursor for the biosynthesis of lignin, flavonoid compounds, and plant defense compounds (isoflavonoids). Therefore, compared with other plant groups, general phenylpropanoid biosynthesis genes in Lauraceae are highly unique, which affects the biosynthesis of flavonoids and lignin and may improve the natural durability of Lauraceae wood. Studies found that functional disruption of CCR and HCT genes affects lignin content (Thévenin et al., 2011; Wang et al., 2015). Although microsyntenic clusters were notable in the CCR and HCT gene families, no unique motifs were found among the protein sequences and promoter regions of homologous genes. We suspected that these Lauraceae-specific genes may have arisen more recently and have not yet diverged significantly from the original genes, in addition, these genes also showed significant TD/PD expansion.
In summary, we investigated the WDR of Lauraceae species by identifying microsynteny clusters among different angiosperm lineages. The Lauraceae-specific biosynthetic genes related to WDR, the conserved motifs of the encoding proteins, the unique and conserved gene expansion and TFBS clusters may play a vital role in increasing and regulating WDR, which may be the main reason for the super decay resistance of Lauraceae. The present genome resources and investigation lay the foundation for molecular breeding or genetic engineering of Lauraceae, and provide key resources for further exploration of the naturally durable wood of Lauraceae species.
Materials and methods
Plant material
A healthy, fruitful, mature L. megaphylla individual was selected and used for whole genome sequencing. This individual was collected from naturally regenerated forest at the National Tree Breeding Station for Nanmu in Zhuxi, Forest Farm of Zhuxi County, Hubei, China. For RNA sequencing, flower buds, stems, buds, and leaves were sampled from healthy trees in the same location, with three replicates per tissue. Tissues were immediately flash frozen and stored at -80 °C for subsequent nucleic acid extractions.
Genome sequencing
For Nanopore sequencing, PromethION libraries were prepared and sequenced on a Nanopore PromethION platform. For Illumina sequencing, 150-bp paired-end (PE) libraries were prepared for sequencing on an Illumina HiSeq X Ten platform. The Hi-C library prepared with the MboI restriction enzyme was sequenced in an Illumina HiSeq X Ten to generate 1488.194 million reads (~223 Gb, roughly 170x coverage of the assembled genome) from 150-bp PE reads. For RNA sequencing, four tissues (flower buds, stems, buds, and leaves) were used to construct mRNA sequencing libraries, after which 150-bp PE sequencing was performed in an Illumina HiSeq X Ten. RNA sequencing produced 996.020 million raw reads (~145 Gb).
More details regarding genome sequencing are available in Note S2.
De novo genome assembly and quality control
De novo genome assembly involved three steps: primary assembly, Hi-C scaffolding, and polishing. First, we used SMARTdenovo (see “URLs” section), WTDBG (version 2.1) (Ruan and Li, 2020), and Canu (version 1.7) (Koren et al., 2017) to generate four of the primary assemblies from ONT long reads. Then, one primary assembly (v0.3, with reasonably sized assembly, fewest contigs, and highest contig N50) was chosen as the optimal assembly, and further polished with three rounds of pilon (see “URLs” section) with clean Illumina reads to generate assembly v1.0. Based on Hi-C data and assembly v1.0, primary scaffolds were produced with 3D-DNA (version 180922) (see “URLs” section). These scaffolds were inspected and manually corrected using Juicebox (version 1.8) (see “URLs” section) and re-scaffolded by 3D-DNA. Afterwards, we optimized the new scaffolds with gap closing using LR_Gapcloser (version 1.1) (see “URLs” section) followed by four rounds of pilon polishing.
Benchmarking Universal Single Copy Orthologs (BUSCO) and LTR Assembly Index (LAI) were used to assess genome completeness and continuity. To evaluate the completeness of the assembly and uniformity of the sequencing, 178 Gb of ONT reads, 160 Gb of clean Illumina reads, and 90 Gb of RNA sequencing reads were aligned to the assembly genome using BWA-MEM (see “URLs” section), minimap2 (Li, 2018), and HiSat2 (version 2.1.0) (see “URLs” section), respectively.
More details of genome assembly are available in Note S3.
Genome annotation
Protein-coding genes were predicted using the MAKER2 pipeline (Holt and Yandell, 2011) including ab initio, homolog proteins, and EST-based prediction methods. We annotated non-coding RNAs (ncRNAs) with several databases and software including tRNAscan-SE (version 1.3.1) (Lowe and Eddy, 1997), RNAMMER (version 1.2) (Lagesen et al., 2007), Rfam database (version 9.1) (see “URLs” section), and BLASTN (version 2.2.28+).
Functions of predicted genes were annotated using sequence similarity searches by BLAT (version 36) (Kent, 2002) with 30% identity and 1e-05 E-value cutoff, as well as domain similarity annotations using InterProScan (version 5.27-66.0) (see “URLs” section). The completeness of genome annotation was assessed using BUSCO. Centurion (Varoquaux et al., 2015) was used to infer the location of all centromeres in the genome based on corrected Hi-C data.
Repeated elements were annotated using RepeatModeler (version 1.0.10) (see “URLs” section) and RepeatMasker (version 4.0.7, rmblast-2.2.28) (see “URLs” section) with homology-based and de novo approaches. In addition, we examined classification, age distribution, birth, and death of LTR-RTs.
More details of genome annotation are available in Note S4 and S5.
Gene family and phylogenetic inference
To determine the phylogenetic relationships among Magnoliids, we used Orthofinder (version 2.3.1) (Emms and Kelly, 2019) to identify gene families from 6 eudicots including Aquilegia coerulea (Filiault et al., 2018), Populus trichocarpa (Tuskan et al., 2006), Arabidopsis thaliana (Michael et al., 2018), Coffea canephora (Denoeud et al., 2014), Solanum lycopersicum (Sato et al., 2012) and Cajanus cajan (Varshney et al., 2012), 4 monocots including Zostera marina (Olsen et al., 2016), Sorghum bicolor (Deschamps et al., 2018), Musa acuminata (D’Hont et al., 2012) and Oryza sativa (Ouyang et al., 2006), 6 Magnoliids including Piper nigrum (Hu et al., 2019), Liriodendron chinense (Chen et al., 2019), Persea americana (Rendon-Anaya et al., 2019), Cinnamomum kanehirae (Chaw et al., 2019a), Litsea cubeba (Chen et al., 2020b) and Lindera megaphylla and 2 outgroup species including Amborella trichopoda (Albert et al., 2013) and Nymphaea colorata (Zhang et al., 2020b). A total of 34,888 orthogroups, including 112 orthologous single-copy gene families and 885 low-copy orthologs with minimum of 83.3% of species having single-copy genes in any orthogroup. Amino acid sequence alignment was performed on these low-copy genes using MUSCLE (version 3.8.31) (Edgar, 2004).
Phylogenetic trees were constructed using concatenation-, coalescent-, and microsynteny-based approaches (Zhao et al., 2021c). For the concatenation-based approach, the maximum likelihood tree was constructed based on concatenated low-copy amino acid sequences with IQ-TREE (version 1.6.7) (Nguyen et al., 2014), employing the best-fit model (-m JTT+F+R5) with ultrafast bootstrapping (-bb 1000). For the coalescent-based approach, gene trees of 855 low-copy gene families were inferred by IQ-TREE. Next, we removed low bootstrap support branches (less than 50%) using the Newick utilities. Then, gene trees were used to construct species trees with ASTRAL-pro. Quartet support of each node was estimated for this coalescent tree. Finally, the microsynteny-based method included two steps. First, after an all-by-all protein alignment of the whole genome was performed using DIAMOND (Buchfink et al., 2015), pairwise synteny blocks were identified using MCScanX (Wang et al., 2012). Then, microsyntenic clusters were detected using Infomap (see “URLs” section).
The maximum likelihood (ML) phylogenetic tree was generated with IQ-TREE (version 1.6.7), using the Mk+R+FO model and ultrafast bootstrapping (-bb 1000). The ML tree constructed using the coalescent-based approach was used as an input tree to estimate divergence time with the MCMCTree program in the PAML package (version 4.9h) (Yang, 2007). Dating was calibrated according to the TimeTree web service (http://www.timetree.org/) by placing soft bounds at four split nodes as constraints for calibrating tree age: (1) the A. trichopoda node (173-199 Mya), (2) L. chinense (117-130 Mya), (3) O. sativa-S. bicolor (42-52 Mya), and (4) P. trichocarpa-A. thaliana (98-177 Mya). Expansion and contraction of gene families were inferred with CAFÉ (version 4.1) (De Bie et al., 2006).
Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Ontology (GO) enrichment analyses were performed using the R package clusterProfiler (version 3.6.0) (Yu et al., 2012).
Additional details are available in Note S6.
Analysis of microsyntenic clusters
The microsyntenic clusters were identified with a computational pipeline previously setup (Zhao et al., 2021c). Key steps in the process are as follows. After an all-vs-all reciprocal sequence similarity search for all annotated genomes using DIAMOND (Buchfink et al., 2015), pairwise synteny block detection was performed using MCScanX (Wang et al., 2012). Then the synteny network was clustered using the Infomap algorithm (see “URLs” section). After that, a synteny cluster matrix was obtained, and the number of each species in each cluster was noted, in which the rows and columns correspond to the various species and clusters, respectively. The matrix was then converted into a binary matrix for phylogenetic inference, where 1 denoted the presence of a specific cluster for the species and 0 denoted its absence. This matrix was analyzed using the cor function in R tools to obtain the correlation coefficient between species. The correlation matrix was then plotted and visualized using the R package corrplot (Wei and Simko, 2017). The states of synteny clusters of 16 species were visualized using the UpSetR package (Conway et al., 2017). Clusters shared among Magnoliids-eudicots, Magnoliids-monocots, and eudicots-monocots were further visualized in Cytoscape (Shannon et al., 2003). To select representative clusters of Magnoliids, eudicots, and monocots, the following criteria were set: for Magnoliids, microsyntenic clusters present in four or more species were reserved; for eudicots, microsyntenic clusters present in four or more species were reserved; for monocots, microsyntenic clusters present in three or more species were reserved.
Genome duplication
We examined genome-wide gene duplications in L. megaphylla, C. kanehirae, and L. cubeba using DupGen_finder (Qiao et al., 2019) with default parameters. The duplicated genes were annotated into five different gene duplication models, including whole-genome duplication (WGD), tandem duplication (TD), proximal duplication (less than 10 gene distance on the same chromosome: PD), transposed duplications (TRD), or dispersed duplications (DSD).
Secondary metabolite biosynthesis pathways
Protein sequences from sequenced Lindera genomes were processed with the Ensemble Enzyme Prediction Pipeline (E2P2) package (version 3.1) (see “URLs” section) to identify putative enzymes. Based on these enzymatic annotations, we then constructed a metabolic pathway database by querying the Plant Metabolic Network (see “URLs” section). The derived pathway database was then validated using SAVI (version 3.1) (Schlapfer et al., 2017) to remove any false positives and redundant pathways, such as non-plant pathway variants, as well as pathways already included in larger pathways. Gene family trees were constructed using IQ-TREE (version 1.6.7) with 1,000 bootstrap replicates. The sequences spanning 2 kb upstream of genes were used to identify transcription factor binding sites (TFBS) in promoters. Putative TF binding sites for suspected promoter sequences were predicted by PlantRegMap (Tian et al., 2019) with q-value ≤ 0.05.
URLs
SMARTdenovo [https://github.com/ruanjue/smartdenovo];
Pilon [http://github.com/broadinstitute/pilon];
3D-DNA (version 180922) [https://github.com/theaidenlab/3d-dna];
Juicebox (version 1.8) [https://github.com/aidenlab/Juicebox];
LR_Gapcloser (version 1.1) [https://github.com/CAFS-bioinformatics/LR_Gapcloser];
BWA-MEM [https://github.com/lh3/bwa];
HiSat2 (version 2.1.0) [https://github.com/infphilo/hisat2];
RepeatMasker [http://www.repeatmasker.org ];
RepeatModeler [http://www.repeatmasker.org ];
Rfam database (version 9.1) [http://eggnogdb.embl.de ];
InterProScan (version 5.27-66.0) [http://www.ebi.ac.uk/InterProScan];
Infomap algorithm [https://github.com/mapequation/infomap];
Timetree web service (http://www.timetree.org/ );
PMN Ensemble Enzyme Prediction Pipeline (E2P2, version 3.1) (https://gitlab.com/rhee-lab/E2P2);
Plant Metabolic Network (https://www.plantcyc.org).
Data availability statement
The data presented in the study are deposited both in the NCBI repository with the accession number SRP382804 (https://www.ncbi.nlm.nih.gov/), and in the Genome Warehouse in National Genomics Data Center with the accession number GWHBKHA00000000 (https://ngdc.cncb.ac.cn/gwh).
Author contributions
J-XL conceived and designed the study; X-CT, J-FG, X-MY, T-LS, SN, S-WZ, Y-TB, Z-CL, and LK prepared the materials and performed related analysis; G-JS provided the specimens; X-CT, J-FG, and J-FM wrote the manuscript; J-FM involved in structuring and polishing the manuscript. All authors contributed to the article and approved the submitted version.
Acknowledgments
We acknowledge support from the Key Program of the National Natural Science Foundation of China (32030010).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2023.1122549/full#supplementary-material
References
Albert, V. A., Barbazuk, W. B., Depamphilis, C. W., Der, J. P., Leebens-Mack, J., Ma, H., et al. (2013). The Amborella genome and the evolution of flowering plants. Science 342, 1241089. doi: 10.1126/science.1241089
Anouhe, J.-B. S., Niamké, F. B., Faustin, M., Virieux, D., Pirat, J.-L., Adima, A. A., et al. (2018). The role of extractives in the natural durability of the heartwood of Dicorynia guianensis amsh: new insights in antioxydant and antifungal properties. Ann. For. Sci. 75, 1–10. doi: 10.1007/s13595-018-0691-0
Bowers, J. E., Chapman, B. A., Rong, J., Li, C., Chai, X., Tu, P.. (2003). Unravelling angiosperm genome evolution by phylogenetic analysis of chromosomal duplication events. Nature 422, 433–438. doi: 10.1038/nature01521
Buchfink, B., Xie, C., Huson, D. H. (2015). Fast and sensitive protein alignment using DIAMOND. Nat. Methods 12, 59–60. doi: 10.1038/nmeth.3176
Cao, Y., Xuan, B., Peng, B., Wang, H.-Y., Lin, C.-Y. I., Wu, C.-S., et al. (2016). The genus Lindera: a source of structurally diverse molecules having pharmacological significance. Phytochem. Rev. 15, 869–906. doi: 10.1007/s11101-015-9432-2
Chaw, S.-M., Liu, Y.-C., Wu, Y.-W., Zhao, C., Wang, P., Xue, L., et al. (2019a). Stout camphor tree genome fills gaps in understanding of flowering plant genome evolution. Nat. Plants 5, 63–73. doi: 10.1038/s41477-018-0337-0
Chen, J., Hao, Z., Guang, X., Gao, M., Wang, J.-Y., Liu, K.-W., et al. (2019). Liriodendron genome sheds light on angiosperm phylogeny and species–pair differentiation. Nat. Plants 5, 18–25. doi: 10.1038/s41477-018-0323-6
Chen, Y.-C., Li, Z., Zhao, Y.-X., Jiang, Y. T., Liu, X. D., Liao, X. Y., et al. (2020b). The Litsea genome and the evolution of the laurel family. Nat. Commun. 11, 1675. doi: 10.1038/s41467-020-15493-5
Chen, S. P., Sun, W. H., Xiong, Y. F., Chen, C.-F.. (2020a). The Phoebe genome sheds light on the evolution of magnoliids. Horticulture Res. 7, 146. doi: 10.1038/s41438-020-00368-z
Chou, C.-J., Lin, L.-C., Chen, K.-T., Hahn, M. W.. (1994). Northalifoline, a new isoquinolone alkaloid from the pedicels of Lindera megaphylla. J. Natural Products 57, 689–694. doi: 10.1021/np50108a001
Conway, J. R., Lex, A., Gehlenborg, N. (2017). UpSetR: an r package for the visualization of intersecting sets and their properties. Bioinformatics 33, 2938–2940. doi: 10.1093/bioinformatics/btx364
De Bie, T., Cristianini, N., Demuth, J. P., Hahn, M. W.. (2006). CAFE: a computational tool for the study of gene family evolution. Bioinformatics 22, 1269–1271. doi: 10.1093/bioinformatics/btl097
Deng, X., Zhao, L., Fang, T., Xiong, Y., Ogutu, C., Yang, D., et al. (2018). Investigation of benzylisoquinoline alkaloid biosynthetic pathway and its transcriptional regulation in lotus. Horticulture Res. 5, 29. doi: 10.1038/s41438-018-0035-0
Denoeud, F., Carretero-paulet, L., Dereeper, A., Droc, G., Guyot, R., Pietrella, M., et al. (2014). The coffee genome provides insight into the convergent evolution of caffeine biosynthesis. Science 345, 1181–1184. doi: 10.1126/science.1255274
Deschamps, S., Zhang, Y., Llaca, V., Ye, L., Sanyal, A., King, M., et al. (2018). A chromosome-scale assembly of the sorghum genome using nanopore sequencing and optical mapping. Nat. Commun. 9, 4844. doi: 10.1038/s41467-018-07271-1
Dewey, C. N. (2011). Positional orthology: putting genomic evolutionary relationships into context. Briefings Bioinf. 12, 401–412. doi: 10.1093/bib/bbr040
D’Hont, A., Denoeud, F., Aury, J.-M., Baurens, F.-C., Carreel, F., Garsmeur, O., et al. (2012). The banana (Musa acuminata) genome and the evolution of monocotyledonous plants. Nature 488, 213–217. doi: 10.1038/nature11241
Edgar, R. C. (2004). MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 32, 1792–1797. doi: 10.1093/nar/gkh340
Ekeuku, S. O., Pang, K.-L., Chin, K.-Y. (2020). Palmatine as an agent against metabolic syndrome and its related complications: a review. Drug Design Dev. Ther. 14, 4963–4974. doi: 10.2147/DDDT.S280520
Emms, D. M., Kelly, S. (2019). OrthoFinder: phylogenetic orthology inference for comparative genomics. Genome Biol. 20, 1–14. doi: 10.1186/s13059-019-1832-y
Filiault, D. L., Ballerini, E. S., Mandáková, T., Aköz, G., Derieg, N. J., Schmutz, J., et al. (2018). The aquilegia genome provides insight into adaptive radiation and reveals an extraordinarily polymorphic chromosome with a unique history. Elife 7, e36426. doi: 10.7554/eLife.36426.050
Hagel, J. M., Facchini, P. J. (2013). Benzylisoquinoline alkaloid metabolism: a century of discovery and a brave new world. Plant Cell Physiol. 54, 647–672. doi: 10.1093/pcp/pct020
Holt, C., Yandell, M. (2011). MAKER2: an annotation pipeline and genome-database management tool for second-generation genome projects. BMC Bioinf. 12, 491. doi: 10.1186/1471-2105-12-491
Hu, L., Xu, Z., Wang, M., Fan, R., Yuan, D., Wu, B., et al. (2019). The chromosome-scale reference genome of black pepper provides insight into piperine biosynthesis. Nat. Commun. 10, 4702. doi: 10.1038/s41467-019-12607-6
Ikezawa, N., Tanaka, M., Nagayoshi, M., Shinkyo, R., Sakaki, T., Inouye, K., et al. (2003). Molecular cloning and characterization of CYP719, a methylenedioxy bridge-forming enzyme that belongs to a novel P450 family, from cultured Coptis japonica cells*. J. Biol. Chem. 278, 38557–38565. doi: 10.1074/jbc.M302470200
Imenshahidi, M., Hosseinzadeh, H. (2020). Berberine neuroprotection and antioxidant activity. oxidative stress and dietary antioxidants in neurological diseases. Eds. Martin, C. R., Preedy, V. R. (London: Academic Press), 199–216.
Inui, T., Kawano, N., Shitan, N., et al. (2012). Improvement of benzylisoquinoline alkaloid productivity by overexpression of 3′-hydroxy-N-methylcoclaurine 4′-o-methyltransferase in transgenic Coptis japonica plants. Biol. Pharm. Bull. 35, 650–659. doi: 10.1248/bpb.35.650
Jagels, R., Visscher, G., Wheeler, E. (2005). An Eocene high arctic angiosperm wood. IAWA J. 26, 387–392. doi: 10.1163/22941932-02603009
Jiao, L., Lu, Y., Zhang, M., Chen, Y., Wang, Z., Guo, Y., et al. (2022). Ancient plastid genomes solve the tree species mystery of the imperial wood “Nanmu“ in the forbidden city, the largest existing wooden palace complex in the world. PLANTS PEOPLE PLANET 4, 696–709. doi: 10.1002/ppp3.10311
Kawaguchi, H., Kim, M., Ishida, M., Ahn, Y.-J., Yamamoto, T., Yamaoka, R., et al. (1989). Several antifeedants from Phellodendron amurense against Reticulitermes speratus. Agric. Biol. Chem. 53, 2635–2640. doi: 10.1080/00021369.1989.10869702
Kent, W. J. (2002). BLAT-the BLAST-like alignment tool. Genome Res. 12, 656–664. doi: 10.1101/gr.229202
Koren, S., Walenz, B., Berlin, K., Miller, J., Bergman, N., Phillippy, A., et al. (2017). Canu: scalable and accurate long-read assembly via adaptive K-mer weighting and repeat separation. Genome Res. 27, 722–736. doi: 10.1101/gr.215087.116
Lagesen, K., Hallin, P., Rødland, E. A., Stærfeldt, H.-H., Rognes, T., Ussery, D. W., et al. (2007). RNAmmer: consistent and rapid annotation of ribosomal RNA genes. Nucleic Acids Res. 35, 3100–3108. doi: 10.1093/nar/gkm160
Li, S. (2014). Transcriptional control of flavonoid biosynthesis: fine-tuning of the MYB-bHLH-WD40 (MBW) complex. Plant Signaling Behav. 9, e27522. doi: 10.4161/psb.27522
Li, H. (2018). Minimap2: pairwise alignment for nucleotide sequences. Bioinformatics 34, 3094–3100. doi: 10.1093/bioinformatics/bty191
Li, S., Zachgo, S. (2013). TCP 3 interacts with R2R3-MYB proteins, promotes flavonoid biosynthesis and negatively regulates the auxin response in Arabidopsis thaliana. Plant J. 76, 901–913. doi: 10.1111/tpj.12348
Lowe, T. M., Eddy, S. R. (1997). tRNAscan-SE: a program for improved detection of transfer RNA genes in genomic sequence. Nucleic Acids Res. 25, 955–964. doi: 10.1093/nar/25.5.955
Lv, Q., Qiu, J., Liu, J., Li, Z., Zhang, W., Wang, Q., et al. (2020). The Chimonanthus salicifolius genome provides insight into magnoliid evolution and flavonoid biosynthesis. Plant J. 103, 1910–1923. doi: 10.1111/tpj.14874
Ma, D., Constabel, C. P. (2019). MYB repressors as regulators of phenylpropanoid metabolism in plants. Trends Plant Sci. 24, 275–289. doi: 10.1016/j.tplants.2018.12.003
Ma, R., Xiao, Y., Lv, Z., Tan, H., Chen, R., Li, Q., et al. (2017). AP2/ERF transcription factor, Ii049, positively regulates lignan biosynthesis in Isatis indigotica through activating salicylic acid signaling and lignan/lignin pathway genes. Front. Plant Sci. 8, 1361. doi: 10.3389/fpls.2017.01361
Meng, X., Wang, Y., Li, J., Jiao, N., Zhang, X., Zhang, Y., et al. (2021). RNA Sequencing reveals phenylpropanoid biosynthesis genes and transcription factors for Hevea brasiliensis reaction wood formation. Front. Genet. 12, 763841–763841. doi: 10.3389/fgene.2021.763841
Michael, T. P., Jupe, F., Bemm, F., Motley, S. T., Sandoval, J. P., Lanz, C., et al. (2018). High contiguity arabidopsis thaliana genome assembly with a single nanopore flow cell. Nat. Commun. 9, 541. doi: 10.1038/s41467-018-03016-2
Millar, D. J., Long, M., Donovan, G., Fraser, P. D., Boudet, A.-M., Danoun, S., et al. (2007). Introduction of sense constructs of cinnamate 4-hydroxylase (CYP73A24) in transgenic tomato plants shows opposite effects on flux into stem lignin and fruit flavonoids. Phytochemistry 68, 1497–1509. doi: 10.1016/j.phytochem.2007.03.018
Mondal, S. K., Roy, S. (2018). Genome-wide sequential, evolutionary, organizational and expression analyses of phenylpropanoid biosynthesis associated MYB domain transcription factors in Arabidopsis. J. Biomolecular Structure Dynamics 36, 1577–1601. doi: 10.1080/07391102.2017.1329099
Mounguengui, S., Saha Tchinda, J.-B., Ndikontar, M., Dumarçay, S., Attéké, C., Perrin, D., et al. (2016). Total phenolic and lignin contents, phytochemical screening, antioxidant and fungal inhibition properties of the heartwood extractives of ten Congo basin tree species. Ann. For. Sci. 73, 287–296. doi: 10.1007/s13595-015-0514-5
Nascimento, M. D., Santana, A., Maranhão, C., Oliveira, L., Bieber, L.. (2013). Phenolic extractives and natural resistance of wood. Biodegradation-Life Sci. 801, 349–370. doi: 10.5772/56358
Nguyen, T.-D., Dang, T.-T. T. (2021). Cytochrome P450 enzymes as key drivers of alkaloid chemical diversification in plants. Front. Plant Sci. 12, 682181. doi: 10.3389/fpls.2021.682181
Nguyen, L.-T., Schmidt, H. A., von Haeseler, A., Minh, B. Q.. (2014). IQ-TREE: a fast and effective stochastic algorithm for estimating maximum-likelihood phylogenies. Mol. Biol. Evol. 32, 268–274. doi: 10.1093/molbev/msu300
Okon, E., Kukula-Koch, W., Jarzab, A., Halasa, M., Stepulak, A., Wawruszak, A., et al. (2020). Advances in chemistry and bioactivity of magnoflorine and magnoflorine-containing extracts. Int. J. Mol. Sci. 21, 1330. doi: 10.3390/ijms21041330
Olsen, J. L., Rouzé, P., Verhelst, B., Lin, Y.-C., Bayer, T., Collen, J., et al. (2016). The genome of the seagrass zostera marina reveals angiosperm adaptation to the sea. Nature 530, 331–335. doi: 10.1038/nature16548
Ou, S., Chen, J., Jiang, N. (2018). Assessing genome assembly quality using the LTR assembly index (LAI). Nucleic Acids Res. 46, e126. doi: 10.1093/nar/gky730
Ouyang, S., Zhu, W., Hamilton, J., Lin, H., Campbell, M., Childs, K., et al. (2006). The TIGR rice genome annotation resource: improvements and new features. Nucleic Acids Res. 35, D883–D887. doi: 10.1093/nar/gkl976
Park, I.-K., Lee, H.-S., Lee, S.-G., Park, J.-D., Ahn, Y.-J.. (2000). Antifeeding activity of isoquinoline alkaloids identified in Coptis japonica roots against Hyphantria cunea (Lepidoptera: Arctiidae) and Agelastica coerulea (Coleoptera: Galerucinae). J. Economic Entomology 93, 331–335. doi: 10.1603/0022-0493-93.2.331
Qiao, X., Li, Q., Yin, H., Qi, K., Li, L., Wang, R., et al. (2019). Gene duplication and evolution in recurring polyploidization–diploidization cycles in plants. Genome Biol. 20, 1–23. doi: 10.1186/s13059-019-1650-2
Rendon-Anaya, M., Ibarra-Laclette, E., Mendez-Bravo, A. (2019). The avocado genome informs deep angiosperm phylogeny, highlights introgressive hybridization, and reveals pathogen-influenced gene space adaptation. Proc. Natl. Acad. Sci. U.S.A. 116, 17081–17089. doi: 10.1073/pnas.1822129116
Robin, A. Y., Giustini, C., Graindorge, M., Matringe, M., Dumas, R.. (2016). Crystal structure of norcoclaurine-6-O-methyltransferase, a key rate-limiting step in the synthesis of benzylisoquinoline alkaloids. Plant J. 87, 641–653. doi: 10.1111/tpj.13225
Ruan, J., Li, H. (2020). Fast and accurate long-read assembly with wtdbg2. Nat. Methods 17, 155–158. doi: 10.1038/s41592-019-0669-3
Ryan, K. G., Swinny, E. E., Markham, K. R., Winefield, C.. (2002). Flavonoid gene expression and UV photoprotection in transgenic and mutant Petunia leaves. Phytochemistry 59, 23–32. doi: 10.1016/S0031-9422(01)00404-6
Sato, S., Tabata, S., Hirakawa, H., Asamizu, E., Shirasawa, K., Isobe, S., et al. (2012). The tomato genome sequence provides insights into fleshy fruit evolution. Nature 485, 635–641. doi: 10.1038/nature11119
Schlapfer, P., Zhang, P., Wang, C., Kim, T., Banf, M., Chae, L., et al. (2017). Genome-wide prediction of metabolic enzymes, pathways and gene clusters in plants. Plant Physiol. 173, 2041–2059. doi: 10.1104/pp.16.01942
Schultz, T. P., Nicholas, D. D. (2000). Naturally durable heartwood: evidence for a proposed dual defensive function of the extractives. Phytochemistry 54, 47–52. doi: 10.1016/S0031-9422(99)00622-6
Shang, J., Tian, J., Cheng, H., Yan, Q., Li, L., Jamal, A., et al. (2020a). The chromosome-level wintersweet (Chimonanthus praecox) genome provides insights into floral scent biosynthesis and flowering in winter. Genome Biol. 21, 1–28. doi: 10.1186/s13059-020-02088-y
Shannon, P., Markiel, A., Ozier, O., Baliga, N. S., Wang, J. T., Ramage, D., et al. (2003). Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 13, 2498–2504. doi: 10.1101/gr.1239303
Simão, F. A., Waterhouse, R. M., Ioannidis, P., Ioannidis, P., Kriventseva, E. V., Zdobnov, E. M., et al. (2015). BUSCO: assessing genome assembly and annotation completeness with single-copy orthologs. Bioinformatics 31, 3210–3212. doi: 10.1093/bioinformatics/btv351
Teng, K., Tan, P., Guo, W., Yue, Y., Fan, X., Wu, J., et al. (2018). Heterologous expression of a novel Zoysia japonica C2H2 zinc finger gene, ZjZFN1, improved salt tolerance in Arabidopsis. Front. Plant Sci. 9, 1159. doi: 10.3389/fpls.2018.01159
The Angiosperm Phylogeny Group, Chase, M. W., Christenhusz, M. J. M., Fay, M. F., Byng, J. W., Judd, W. S., et al. (2016). An update of the angiosperm phylogeny group classification for the orders and families of flowering plants: APG IV. Botanical J. Linn. Soc. 181, 1–20. doi: 10.1111/boj.12385
Thévenin, J., Pollet, B., Letarnec, B., Saulnier, L., Gissot, L., Maia-Grondard, A., et al. (2011). The simultaneous repression of CCR and CAD, two enzymes of the lignin biosynthetic pathway, results in sterility and dwarfism in Arabidopsis thaliana. Mol. Plant 4, 70–82. doi: 10.1093/mp/ssq045
Tian, F., Yang, D.-C., Meng, Y.-Q., Jin, J., Gao, G.. (2019). PlantRegMap: charting functional regulatory maps in plants. Nucleic Acids Res. 48, D1104–D1113. doi: 10.1093/nar/gkz1020
Tuskan, G. A., Difazio, S., Jansson, S., Bohlmann, J., Grigoriev, I., Hellsten, U., et al. (2006). The genome of black cottonwood, populus trichocarpa (Torr. & Gray). Science 313, 1596–1604. doi: 10.1126/science.1128691
Vance, C., Kirk, T., Sherwood, R. (1980). Lignification as a mechanism of disease resistance. Annu. Rev. Phytopathol. 18, 259–288. doi: 10.1146/annurev.py.18.090180.001355
Van de Peer, Y. (2004). Computational approaches to unveiling ancient genome duplications. Nat. Rev. Genet. 5, 752–763. doi: 10.1038/nrg1449
Varoquaux, N., Liachko, I., Ay, F., Burton, J. N., Shendure, J., Dunham, M. J., et al. (2015). Accurate identification of centromere locations in yeast genomes using Hi-c. Nucleic Acids Res. 43, 5331–5339. doi: 10.1093/nar/gkv424
Varshney, R. K., Chen, W., Li, Y., Bharti, A. K., Saxena, R. K., Schlueter, J. A., et al. (2012). Draft genome sequence of pigeonpea (Cajanus cajan), an orphan legume crop of resource-poor farmers. Nat. Biotechnol. 30, 83. doi: 10.1038/nbt.2022
Wang, G.-F., He, Y., Strauch, R., Olukolu, B. A., Nielsen, D., Li, X., et al. (2015). Maize homologs of hydroxycinnamoyltransferase, a key enzyme in lignin biosynthesis, bind the nucleotide binding leucine-rich repeat Rp1 proteins to modulate the defense response. Plant Physiol. 169, 2230–2243. doi: 10.1104/pp.15.00703
Wang, Y., Tang, H., Debarry, J. D., Tan, X., Li, J., Wang, X., et al. (2012). MCScanX: a toolkit for detection and evolutionary analysis of gene synteny and collinearity. Nucleic Acids Res. 40, e49. doi: 10.1093/nar/gkr1293
Wei, T., Simko, V. (2017). R package “corrplot”: Visualization of a correlation matrix (Version 0.84). Available from: http://CRAN.R-project.org/package=corrplot.
Wickett, N. J., Mirarab, S., Nguyen, N., Warnow, T., Carpenter, E., Matasci, N., et al. (2014). Phylotranscriptomic analysis of the origin and early diversification of land plants. Proc. Natl. Acad. Sci. 111, E4859–E4868. doi: 10.1073/pnas.1323926111
Xie, J., Qi, J., Huang, X., Zhou, N., Hu, Y.. (2015). Comparative analysis of modern and ancient buried Phoebe zhennan wood: surface color, chemical components, infrared spectroscopy, and essential oil composition. J. Forestry Res. 26, 501–507. doi: 10.1007/s11676-015-0034-z
Yadav, V., Wang, Z., Wei, C., Amo, A., Ahmed, B., Yang, X., et al. (2020). Phenylpropanoid pathway engineering: An emerging approach towards plant defense. Pathogens 9, 312. doi: 10.3390/pathogens9040312
Yamada, Y., Kokabu, Y., Chaki, K., Yoshimoto, T., Ohgaki, M., Yoshida, S., et al. (2011). Isoquinoline alkaloid biosynthesis is regulated by a unique bHLH-type transcription factor in Coptis japonica. Plant Cell Physiol. 52, 1131–1141. doi: 10.1093/pcp/pcr062
Yang, Z. (2007). PAML 4: phylogenetic analysis by maximum likelihood. Mol. Biol. Evol. 24, 1586–1591. doi: 10.1093/molbev/msm088
Yu, G., Wang, L. G., Han, Y., He, Q. Y.. (2012). clusterProfiler: an r package for comparing biological themes among gene clusters. OMICS-A J. Integr. Biol. 16, 284–287. doi: 10.1089/omi.2011.0118
Zeng, L., Zhang, Q., Sun, R., Kong, H., Zhang, N., Ma, H., et al. (2014). Resolution of deep angiosperm phylogeny using conserved nuclear genes and estimates of early divergence times. Nat. Commun. 5, 1–12. doi: 10.1038/ncomms5956
Zhang, L., Chen, F., Zhang, X., Li, Z., Zhao, Y., Lohaus, R., et al. (2020b). The water lily genome and the early evolution of flowering plants. Nature 577, 79–84. doi: 10.1038/s41586-019-1852-5
Zhang, C., Scornavacca, C., Molloy, E. K., Mirarab, S.. (2020a). ASTRAL-pro: quartet-based species-tree inference despite paralogy. Mol. Biol. Evol. 37, 3292–3307. doi: 10.1093/molbev/msaa139
Zhao, T., Li, R., Yao, W., Wang, Y., Zhang, C., Li, Y., et al. (2021a). Genome-wide identification and characterisation of phenylalanine ammonia-lyase gene family in grapevine. J. Hortic. Sci. Biotechnol. 96, 456–468. doi: 10.1080/14620316.2021.1879685
Zhao, T., Schranz, M. E. (2019). Network-based microsynteny analysis identifies major differences and genomic outliers in mammalian and angiosperm genomes. Proc. Natl. Acad. Sci. United States America 116, 2165–2174. doi: 10.1073/pnas.1801757116
Zhao, T., Wu, T., Pei, T., Wang, Z., Yang, H., Jiang, J., et al. (2021b). Overexpression of SlGATA17 promotes drought tolerance in transgenic tomato plants by enhancing activation of the phenylpropanoid biosynthetic pathway. Front. Plant Sci. 12, 634888–634888. doi: 10.3389/fpls.2021.634888
Zhao, T., Zwaenepoel, A., Xue, J. Y., Kao, S. M., Li, Z., Schranz, M. E., et al. (2021c). Whole-genome microsynteny-based phylogeny of angiosperms. Nat. Commun. 12, 3498. doi: 10.1038/s41467-021-23665-0
Zhou, M., Memelink, J. (2016). Jasmonate-responsive transcription factors regulating plant secondary metabolism. Biotechnol. Adv. 34, 441–449. doi: 10.1016/j.biotechadv.2016.02.004
Keywords: Lauraceae, Lindera megaphylla, wood decay resistance (WDR), tandem and proximal duplications (TD/PD), gene microsynteny
Citation: Tian X-C, Guo J-F, Yan X-M, Shi T-L, Nie S, Zhao S-W, Bao Y-T, Li Z-C, Kong L, Su G-J, Mao J-F and Lin JX (2023) Unique gene duplications and conserved microsynteny potentially associated with resistance to wood decay in the Lauraceae. Front. Plant Sci. 14:1122549. doi: 10.3389/fpls.2023.1122549
Received: 13 December 2022; Accepted: 10 February 2023;
Published: 08 March 2023.
Edited by:
Guanjing Hu, Agricultural Genomics Institute at Shenzhen (CAAS), ChinaReviewed by:
Qi Wu, Chengdu University, ChinaXiwen Li, Institute of Chinese Materia Medica, China Academy of Chinese Medical Sciences, China
Copyright © 2023 Tian, Guo, Yan, Shi, Nie, Zhao, Bao, Li, Kong, Su, Mao and Lin. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Jinxing Lin, bGluanhAYmpmdS5lZHUuY24=; Jian-Feng Mao, amlhbmZlbmcubWFvQGJqZnUuZWR1LmNu
†These authors have contributed equally to this work