Abstract
Lindera aggregata, a member belongs to the genus Lindera of Lauraceae family. Its roots and leaves have been used as traditional Chinese medicine or functional food for thousands of years. However, its mitochondrial genome has not been explored. Our aim is to sequence and assemble the mitogenome of L. aggregata to elucidate the genetic mechanism and evolutionary pathway. The results had shown that the mitogenome was extremely complex and had a unique multi-branched conformation with total size of 912,473 bp. Comprehensive analysis of protein coding genes of 7 related species showed that there were 40 common genes in their mitogenome. Interestingly, positive selection had become an important factor in the evolution of ccmB, ccmFC, rps10, rps11 and rps7 genes. Furthermore, our data highlighted the repeated trend of homologous fragment migrations between chloroplast and mitochondrial organelles, and 38 homologous fragments were identified. Phylogenetic analysis identified a tree that was basically consistent with the phylogeny of Laurales species described in the APG IV system. To sum up, this study will be helpful to the study of population genetics and evolution of Lindera species.
1 Introduction
The genus Lindera is a member of the Lauraceae family, comprising approximately 100 species distributed in tropical, subtropical and temperate regions of Asia, as well as temperate regions of North America (Hung-Pin, 1987). There are more than 40 species of this genus in China, accounting for 46% of the whole genus (Tao et al., 2024). The plants of Lindera are widely used in traditional medicine and have high economic value, and Lindera aggregata (Sims) Kosterm is a representative one of them. It is also known as “Wu-Yao”, is a common folk medicinal plant (Lv et al., 2023).
The leaves of L. aggregata are edible and are often used as functional tea or dietary supplements, with health effects such as anti-liver injury and lipid-lowering (; Weng et al., 2022). The dried root of L. aggregata, as a traditional Chinese herbal medicine, was first recorded in the BenCaoShiYi in the Tang Dynasty, which had the function of regulating qi and relieving pain, warming the kidney and dispelling cold (). It mainly grows in eastern, central, southern and southwestern China, such as Zhejiang, Jiangxi, Anhui and other provinces. Among them, BenCaoGangMu recorded: “The Wu-Yao originated from Tiantai was the best”. Therefore, Zhejiang Tiantai has been a famous genuine producing area of L. aggregata since ancient times, which is called “Tiantai Wu-Yao”. The main active components of L. aggregata are sesquiterpenes, alkaloids, flavonoids and volatile oils, which have the effects of protecting liver, anti-inflammation, anti-virus, antibacterial, anti-tumor and anti-oxidation (Salleh, 2020). Due to the increasing use of authentic L. aggregata and the lack of centralized protection, the reserves of authentic L. aggregata are decreasing day by day, so intensive management is extremely urgent. According to the Chinese Pharmacopoeia, “the taproot, which is not spindle-shaped, cannot be used medicinally”. Moreover, the growing year of L. aggregata is very long, and the quality of traditional Chinese sold in the marker is often poor, which leads to great difficulties in germplasm selection. Therefore, it is essential to effectively identify authentic medicinal materials in traditional Chinese medicine. In recent years, DNA bar code technology has been widely used in the classification and identification of species, kinship and diversity analysis (Hebert et al., 2003). The primitive herbs used in traditional Chinese medicine are also the object of molecular research (Lu et al., 2016). For example, the nad7 gene of the mitogenome can be used to identify the Korean ginseng cultivar “Chunpoong” and other ginseng (Wang et al., 2009).
The genus Lindera is a large and complex genus with similar transitional characters among different species. In the past, the molecular phylogeny of Lindera was mainly focused on the chloroplast genome, which revealed that Lindera was a polyphyletic group (Zhao et al., 2018; Shi et al., 2024), and its systematic relationship with related species from Litsea and Laurus is still unclear. Mitochondrial genome analysis is very important for understanding the evolution and genome structure of various plants (Yang et al., 2023). Recent research on the mitogenome of the genus Cinnamomum has further discussed the evolutionary location of the Magnoliids species (; ), but information on the mitogenome of other genera from Magnoliids is still lacking. In particular, the mitochondrial genome of Lindera plants has not been reported so far. It is an urgent task to analyze the genomic characteristics of L. aggregata, compare and analyze the evolution of the mitogenome of Lindera and other Laurales plants.
Plant chloroplasts and mitochondria are two important semi-autonomous organelles except for the nucleus, which have their own independent genetic material and genetic system (). Mitochondria are important organelles responsible for energy conversion and the main place for aerobic respiration, which play an important role in plant growth and development and stress response (Millar et al., 2005; Wang et al., 2024). Compared with chloroplasts, the size and shape of the mitochondrial genome is very irregular, and its content is equivalent to 1%-2% of the nuclear genome in terms of size. Angiosperms are considered to be the largest group of mitochondrial genomes so far, and the sizes of mitochondrial genomes in different plants vary significantly, ranging from 66 kb to 11.3 Mb (Sloan et al., 2012; Skippington et al., 2015), but the coding genes are relatively conservative (Ogihara et al., 2005). Therefore, the differences in the length of mitogenome from different species mainly depend on the repeat sequences, especially the non-coding regions, including horizontal gene transfer (HGT) (). Due to the frequent recombination in the mitochondrial genome of higher plants, the structure of the mitogenome is also very different, resulting in the emergence of subring and isomers (). At the same time, the frequently recombined repeat sequences in most plant mitogenome are mainly non-coding DNA, which is closely related to genome recombination (Li et al., 2022). In addition, there are a large number of RNA editing phenomena in higher plant mitochondria, which are one of the necessary steps of plant mitochondrial gene expression (Kubo and Newton, 2008). At present, it has been found that the most types of RNA editing sites are C to U, which mainly occur in the protein coding region, and most editors change the type of amino acids to produce the desired function (; ), while the RNA editing sites in some chloroplasts are relatively conservative in evolution (; ).
In this study, the complete sequence of mitochondrial genome of L. aggregata was determined and assembled by the third generation PacBio HiFi sequencing technique. We obtained a multi-branched conformation map of mitogenome and dismantled it into a master circle and a linear fragment. Furthermore, seven mitogenome of Laurales were obtained from NCBI database and compared with the mitogenome of L. aggregata, including sequence size, GC content, haplotype network, Ka/Ks, codon usage bias, repeat sequences and RNA editing site prediction. The events of rearrangement in the evolution of mitogenome in 8 Laurales species were compared by using collinear maps. The fragment migrations between chloroplast and mitochondrial genome of L. aggregata were also analyzed. Finally, maximum likelihood and Bayesian inference phylogenetic trees were constructed using the conserved PCGs in the mitogenome of 30 species. These results provide valuable information for the mitogenome of Lindera and will be of great value to the classification, species identification and molecular breeding of this genus with high value.
2 Materials and methods
2.1 Sample collection, DNA extraction and sequencing
The sample materials were collected from the planting base of Lindera aggregata from Zhejiang Hongshiliang Group Tiantai Mountain Wu-Yao Co., Ltd., Zhejiang, China (120°48.23′ E, 29°12.71′ N). Fresh and tender leaves were collected and quick-frozen with liquid nitrogen. Genome DNA was extracted by improved CTAB method (). The quality and concentration of total DNA was detected by 1% agarose gel electrophoresis, NanoDrop spectrophotometer and Qubit v4 fluorometer. The qualified high-quality DNA samples (main band > 30 kb) were randomly broken into fragments (15-18 kb) by Covaris ultrasonic crusher; the segmented DNA were enriched and purified by magnetic beads, and the damage repair and end repair were carried out; the stem ring sequencing adapters were connected at both ends of the DNA fragment, and the failed fragments were removed by exonuclease. Finally, the constructed library was used for single molecule real-time (SMRT) sequencing on the PacBio Sequel II platform. The raw reads obtained subreads after removing the adapters and filtering, and the ccs software was used to generate high-precision HiFi reads (91.8 G) with parameter “min-passes = 3, min-rq = 0.99”.
2.2 Assembly and annotation of the mitochondrial genome
The mitochondrial genome was assembled by PMAT v1.5.3 software (), and the initial mitochondrial assembly map was obtained. Bandage v0.8.1 software (Wick et al., 2015) was used to display the mitochondrial genome and manually remove the extended segments of the chloroplast and nuclear genome. Then, the HiFi reads were compared with the mitochondrial genome by Minimap2 v2.24 tool (Li, 2021) to help analyze the repetitive regions in the mitochondrial genome, and finally produce a complete mitogenome of L. aggregata. Intelligent Plant Mitochondrial Genome Annotator (IPMGA, http://www.1kmpg.cn/ipmga/) was used to annotate the mitochondrial genome. tRNA and rRNA in mitochondrial genome were annotated by tRNAscan-SE v2.0 () and BLASTN software (), respectively. The mitochondrial genome was visualized by Plant Mitochondrial Genomes Map (PMGmap) online software (Zhang et al., 2024), and the cis- and trans-splicing gene maps were drawn. Finally, the mitogenome of L. aggregata was uploaded to the GenBank database to obtain the accession number.
2.3 Comparative analysis of 8 plants mitogenome from Laurales
We downloaded all the published mitogenome of Laurales (Cinnamomum camphora, C. chago, C. chekiangense, C. insularimontanum, Caryodaphnopsis henryi, Machilus pauhoi and Hernandia nymphaeifolia) from the NCBI database and re-annotation of them using the same method mentioned above for statistical and subsequent analysis. The protein coding gens (PCGs) of mitogenome from 7 Laurales plants and L. aggregata were extracted by PhyloSuit v1.2.3 software (Zhang et al., 2020). The multi-sequence alignment was performed among 40 common PCGs by MAFFT v7.149b software (Katoh and Standley, 2013), and the sequences were trimmed by trimAl v1.2 software (). The trimmed sequences were merged, and the maximum likelihood (ML) phylogenetic tree was constructed by IQ-TREE v2.2.0.3 software (Minh et al., 2020) with default parameters. Furthermore, through the Genepioneer Cloud (http://cloud.genepioneer.com:9929/#/home), the aligned sequences were removed from the sites with gap/missing and the sites with no different bases at the same position in all the sequences, while the remaining bases were haplotype sequences. Then, these haplotype sequences were analyzed using DnaSP v6.12.03 software (Rozas et al., 2003) to generate NEXUS haplotype data files. Haplotype networks of 40 conserved protein coding sequences were constructed via PopART v4.8.4 software (Leigh and Bryant, 2015). The ratio of non-synonymous to synonymous mutations (Ka/Ks) of 40 common PCGs in 8 species mitogenome from Laurales were calculated by KaKs_Calculator v2.0 software (Wang et al., 2010).
2.4 Analysis of codon usage bias
In order to analyze the codon usage bias of Laurales species, the relative synonymous codon usage (RSCU) values of PCGs and the number of codon corresponding to each amino acid in 8 species were counted by PhyloSuit v1.2.3 software (Zhang et al., 2020), and visualized by TBtools (). Moreover, in order to determine the GC content and effective number of codon (ENC) of the coding genes, we used the online tool ENC calculation (http://cloud.genepioneer.com:9929/#/tool) to obtain the parameters of 8 species and draw the ENC-plot. The ENC values indicated the extent to which the codon usage of the genome deviates from random selection (Wright, 1990). ENC-plot analysis can explore the relationship between ENC and GC3 distribution, which was an effective way to visualize the codon usage bias of genetic data. The standard curve (ENC = 2 + GC3 + 29/[GC32 + (1 - GC3)2]) (Romero et al., 2000) indicated that codon usage bias was completely determined by mutation when there was no selection pressure, that was, codon usage bias was completely determined by nucleic acid sequence composition. Each data point in the scatter graph represented a gene. The specific criterion was the distance between the scatters and the standard curve, and the close distance between the scatters and the standard curve indicated that the codon usage bias mainly determined by the base composition and was weakly affected by translation selection; the long distance indicated that the codon ENC value was low and had a strong significant correlation with the level of gene expression, and the codon usage bias strong (Sueoka, 1988).
2.5 Analysis of repeat fragments and prediction of RNA editing sites
Simple repetitive sequences (SSRs) were identified by the online tool MISA (https://webblast.ipk-gatersleben.de/misa/) (). The parameters were set as follows: the minimum number of repeats for the mononucleotide repeats, dinucleotide repeats, trinucleotide repeats, tetranucleotide repeats, pentanucleotide repeats and hexanucleotide repeats were set to 10, 5, 4, 3, 3 and 3, respectively. The tandem repeats were identified by Tandem Repeat Finder online tool (https://tandem.bu.edu/trf/home) () with default parameters. Dispersed repeats (including forward repeats, reverse repeats, palindromic repeats and complementary repeats) were analyzed by online tool REPuter (https://bibiserv.cebitec.uni-bielefeld.de/reputer) (Kurtz et al., 2001) with the parameters set as follows: minimum repeat length greater than 30 bp, maximum repeat length of 5000 bp and Hamming distance was 3. Moreover, in order to further observe the RNA editing events, the Deepred-mt tool () was used to predict C-to-U RNA editing sites in 40 PCGs of L. aggregata, and the prediction value was set to 0.9.
2.6 Identification and collinear analysis of homologous fragments
Based on our previously assembled and annotated chloroplast genome of L. aggregata (PP199190) (Shi et al., 2024), the homologous sequences were searched between chloroplast genome and mitochondrial genome using BLAST software with the parameter settings of word size 7 and E value = 1e−5. To further study the evolutionary relationship among species, we identified the homologous sequences among L. aggregata and other 7 species of Laurales by BLAST program with same parameters as above. The collinear maps within L. aggregata (mitochondria and chloroplast) and among different species were visualized by TBtools software (). Furthermore, we used Mauve v2.4.0 software () to detect the rearrangement of 7 Laurales mitogenome, and the L. aggregate mitogenome was used as the reference genome.
2.7 Phylogenetic analysis
We downloaded 29 complete mitogenome sequences of 9 different orders (Laurales, Magnoliales, Piperales, Chloranthales, Amborellales, Nymphaeales, Acorales, Poales, and Apiales) from NCBI database. Amborella trichopoda and Nymphaea colorata were used as outgroups to construct a phylogenetic tree with mitogenome of L. aggregata. Firstly, the conserved PCGs were extracted by PhyloSuit software. Then homologous sequence alignment was performed among conserved PCGs in 30 mitogenome by MAFFT software. The Gblocks software was used to trim the aligned sequences. Finally, the phylogenetic trees based on ML and Bayesian inference (BI) were constructed, respectively. The ML tree was built by IQ-TREE software with 5000 ultra-fast bootstraps were set, and the optimal model (GTR+F+I+G4) was obtained by ModelFinder program (Kalyaanamoorthy et al., 2017). The MrBayes v3.2.7a software (Ronquist et al., 2012) was used to construct BI phylogenetic tree, and the best model (GTR+I) was selected by jModelTest tool (). MCMC was run for 5,000,000 generations, with sampling every 1000 generations. The number of chains and runs was set to 4 and 3, respectively. In the total data, the top 25% of the trees aged and discarded, and the rest were used as a 50% consistent tree; convergence was determined by keeping the average standard deviation of the crossover frequencies < 0.01. ML and BI phylogenetic trees were visualized using Chiplot online tool (Xie et al., 2023). The accession numbers of 30 species were recorded in Supplementary Table S4.
3 Results
3.1 Assembly and annotation of L. aggregata mitochondrial genome
Using the PacBio Sequel II platform, we obtained a total of 5,327,958HiFi reads, with the maximum and average lengths of 59,940 bp and 17,227 bp, respectively. We used HiFi reads to assemble the complete mitochondrial genome of L. aggregata. The structure of mitogenome was very complex and had multi-branched conformation (Figure 1A), and we analyzed one type of conformation. The total length of the genome sequence was 912,473 bp, the CG content was 46.83%, and it was composed of a longer circular master circle and a shorter linear fragment (Figure 1B), with lengths of 868,093 bp (chromosome 1) and 44,380 bp (chromosome 2), respectively. The mitochondrial genome was further annotated, reveled a total of 40 different protein coding genes (PCGs), 28 tRNA genes and 3 rRNA genes (Figure 1C; Supplementary Table S1). Among the 40 PCGs, 23 were core genes and 17 were non-core genes. The 23 core genes include 5 ATP synthase genes (atp1, atp4, atp6, atp8, atp9), 4 cytochrome c biosynthesis genes (ccmB, ccmC, ccmFC, ccmFN), 1 ubichinol cytochrome c reductase gene (cob), 3 cytochrome c oxidase genes (cox1, cox2, cox3), 1 transport membrance protein (mttB) and 9 NADH dehydrogenase genes (nad1, nad2, nad3, nad4, nad4L, nad5, nad6, nad7, and nad9). However, it lacked a matures gene (matR) as same as the Cinnamomum species (; ). The non-core genes include 4 large subunits of the ribosomal proteins (rpl10, rpl16, rpl2, rpl5), 11 small subunits of the ribosomal proteins (rps1, rps10, rps11, rps12, rps13, rps14, rps19, rps2, rps3, rps4, rps7) and 2 succinate dehydrogenase genes (sdh3, sdh4). Furthermore, 8 cis-splicing genes (ccmFC, cox2, nad4, nad7, rpl2, rps10, rps3 and trnL-UAG) and 3 trans-splicing genes (nad1/2/5) were identified. The details of exon in each gene can be found in Supplementary Figure S1.
Figure 1
3.2 Comparative analysis of mitogenome among L. aggregata and other 7 Laurales species
In order to further explore the evolutionary characteristics of the mitochondrial genome, we compared it with seven other species of Laurales (Cinnamomum camphora, C. chago, C. chekiangense, C. insularimontanum, Caryodaphnopsis henryi, Machilus pauhoi and Hernandia nymphaeifolia). There were great differences in mitogenome size among 7 Laurales species, ranging from 535,805 bp (H. nymphaeifolia) to 1,168,029 bp (C. henry), and the GC content ranged from 45.73% (H. nymphaeifolia) to 47.12% (C. chago) (Supplementary Figure S2). In addition, the contents of GC, GC1, GC2 and GC3 of the PCGs were between 44.08% and 44.21%, 49.14% and 49.35%, 44.66% and 45.03%, 38.22% and 38.52%, respectively. For the specific values of each species were found in Supplementary Table S2.
It was worth noting that although there were wide differences in the total number of genes in the mitogenome from these eight Laurales plants, the number of PCGs remained relatively consistent (Supplementary Table S2). According to the comparative analysis, it can be found that H. nymphaeifolia had one more PCGs (matR) than other species, and had the largest number of genes (71), while C. camphora had the least number of genes (62). The total differences were mainly caused by tRNA, while the number of rRNA all were 3. In order to further analyze these mitogenome, we extracted the coding sequences. 40 shared PCGs were identified in all 8 species and used to construct ML phylogenetic tree (Supplementary Figure S3). Interestingly, our results showed that there was a close sister relationship between L. aggregata and four species of Cinnamomum.
We also analyzed 40 common PCGs mutation sites among these 8 Laurales species to construct haplotype networks for more in-depth study of genetic variation (Figure 2). It was worth noting that H. nymphaeifolia showed the most typical haplotype composition, followed by C. henryi, while other species showed more mixed composition. Although L. aggregata had 20 unique haplotypes, it also showed mixed haplotypes with M. pauhoi and C. chago, including 6 haplotypes shared with M. pauhoi, 4 haplotypes shared with C. chago, and 10 haplotypes shared with 2 of them. Furthermore, C. insularimontanum and C. chekiangense shared all genes.
Figure 2
In order to evaluate the effect of environmental stress on the evolution of mitochondrial genome. We calculated the Ka/Ks ratio values of 40 PCGs shared by 8 Laurales species to observe the evolutionary effect of mitogenome under environmental stress (Figure 3). The results showed that the Ka/Ks ratio values of most genes were less than 1, indicating that they under purifying selection, while the Ka/Ks ratio values of ccmB, ccmFC, rps10, rps11 and rps7 genes were greater than 1 in all species, indicating positive selection. Moreover, a small number of genes were selected in the opposite direction in different species.
Figure 3
3.3 Analysis of codon usage bias
Codon usage bias analysis of 8 mitochondrial genomes in Laurales showed that 61 codons were detected, encoding 20 amino acids (Figure 4). We determined that PCGs had a total of 91,645 codons in 8 mitogenome, distributed between 11,033 and 11,839. Among them, the most frequently used codon was GCU, and the trend was consistent in 8 species. In the 20 kinds of amino acids, serine, leucine and arginine had more codons and had 6 types of codons. Among them, serine had the largest codon number, with 9039 codons (9.86%), followed by leucine with 8780 codons (9.58%), and tryptophan had the least codon number, with only 1241 codons (1.35%) (Figure 4B). In addition, the relative synonymous codon usage (RSCU) values of 28 codons in 8 species was greater than 1, indicating that the frequency of these codons was higher, 28 codons were lower than expected (RSCU < 1), while 3 codons (AGU, ACC and ACA) had different usage bias in different species. Interestingly, methionine (AUG) and tryptophan (UGG) did not show codon usage bias (RSCU = 1) (Figure 4A). In terms of amino acids, with the exception of AUG (methionine) and UGG (tryptophan), most amino acids showed a bias in their codon usage patterns.
Figure 4
In order to further evaluate the codon usage bias of each gene in the mitogenome and the factors affecting the usage pattern, we extracted 40 PCGs to calculated the GC content of the first, second and third sites and effective number of codons (ENC) of these genes. The results showed that the GC1 values were ranged from 38.12% to 58.24%, the GC2 values were ranged from 34.71% to 56%, and the GC3 values were ranged from 27.35% to 46.88%. The average contents of GC in different locations (GC1, GC2 and GC3) were all less than 50%, indicating that there was a bias for A/T base and A/T ending codon in the mitogenome of L. aggregata. The values of ENC were between 39.21% and 60.77%, with an average of 52.98%, indicating that there was a weak codon usage bias in the mitogenome of L. aggregata (Figure 5A). Furthermore, through ENC-plot analysis, the results showed that most PCGs were far from the standard curve, and only a small number of genes were distributed beside the curve (Figure 5B), indicating that there was a strong significant correlation between most genes and the level of gene expression, and the codon had preference, which was consistent with the above results (Figure 4).
Figure 5
3.4 Analysis of repeat fragments and prediction of RNA editing sites
Simple sequence repeats (SSRs) were widely distributed in the mitochondrial genome. There were 169-401 SSRs detected across each of mitogenome from 8 Laurales species, for a total of 2259 SSRs. Although there were great differences in the number of SSRs among the 8 Laurales species, they all had 6 constituent types, and the occupation trend of each type was relatively consistent in each species (Figure 6A). Among them, the number of tetranucleotide repeats was the largest (55-144 loci), accounting for 32.54%-38.83% of the total SSRs. The number of hexanucleotide repeats was the least (5-10 loci), accounting for 1.91%-4.14% of the total SSRs. Furthermore, we also detected tandem and dispersed repeats in the mitogenome of 8 species. Among them, the distribution of tandem repeats was between 40 and 81, and there was a great difference in quantity. The number of dispersed repeats ranged from 675 to 2153, in which except H. nymphaeifolia only had forward (F) and palindromic (P) repeats, the other species contained four types repeats, and the number of F and P repeats were far more than reverse (R) and complementary (C) repeats (Figure 6B). The distribution trend of dispersed repeat sequence lengths among the 8 species was consistent, mainly between 30 and 60 bp, accounting for 72.13%-85.19% of the total (Figure 6D). Moreover, C. henryi and H. nymphaeifolia had the most and fewest repeat fragments respectively, which might be an important reason for the size of the mitochondrial genome.
Figure 6
In this study, we also predicted all possible C-to-U RNA editing sites in the L. aggregata mitogenome to reveal more deeply the gene expression in the mitochondrial genome. A total of 770 RNA editing sites were found, mainly involving the transition from nucleotide C to U (Figure 6C). It is worth noting that the nad4 gene showed the highest RNA editing frequency, predicting 68 editing sites, followed by the nad5 gene with 48 editing sites. However, most ribosomal proteins had a small number of editing sites, which was consistent with C. chekiangense and C. camphora (; ). In addition, we observed that RNA editing events mainly occurred at the first and second base positions, and the second base positions changed more frequently. These editing events lead to changes in amino acids, such as Proline (P) to Serine (S) or Leucine (L), Serine (S) to Phenylalanine (F) or Leucine (L), Arginine (R) to Cysteine (C), Tryptophan (Y) or Termination codon (*), Leucine (L) to Alanine (A), Threonine (T) to Isoleucine (I) or Methionine (M), Histidine (H) to Tyrosine (Y), Glutamine (Q) to Termination codon (*), Alanine (A) to Valine (V). Interestingly, most of these amino acid changes were related to the transformation of hydrophobic amino acids, which helped to improve the stability of proteins.
3.5 DNA fragment transfer analysis between mitochondrial and chloroplast genomes
In the evolution of higher plants, the transfer of genetic material between mitochondria and chloroplasts were a common phenomenon. However, it was worth noting that these sequence fragments originating from chloroplast organelles showed relatively low conservatism. Therefore, we searched the homologous fragments between the two organelles in order to explore the migration of sequences from chloroplasts to mitochondrial organelles (Figure 7A). Through sequence similarity analysis, a total of 38 homologous fragments were identified between chloroplast and mitochondrial genomes, with sequences length ranged from 32 bp to 1721 bp, with a total of 12,022 bp, accounting for 1.32% of the total mitogenome (Supplementary Table S3). The number of mismatches in these fragments varied from 0 to 119, and the number of gaps varied from 0 to 33. These fragments represented the migrations from chloroplasts to mitochondrial organelles, called MtPts. Of the 38 migration segments, MtPt-1 was the longest. We annotated these fragments and identified 15 complete genes, including 2 PCGs (petN (2)) and 13 tRNA genes (trnA-UGC, trnD-GUC, trnE-UUC, trnH-GUG, trnI-CAU, trnM-CAU, trnN-GUU (2), trnP-UGG, trnR-UCU, trnT-GGU, trnW-CCA, trnY-GUA). It was worth noting that the petN gene in the cpDNA contributed subunits of cytochrome b/f complex, while 13 tRNA genes might have been lost or experienced pseudogene changes in the cpDNA.
Figure 7
3.6 Collinearity analysis and structural rearrangement
In order to further study the homologous sequences of the mitogenome in 8 Laurales species, we identified the collinearity between L. aggregata and other Laurales species (Figure 7B and Supplementary Figure S4). The results showed that L. aggregata had a large number of collinear fragments with other 7 species in Laurales. Notably, the highest collinearity was detected between L. aggregate and C. camphora, which contained 1051 collinearity fragments totaling 598,706 bp, accounting for 65.61% of L. aggregate mitogenome. However, H. nymphaeifolia exhibited the weakest collinearity, with only 378 collinearity fragments spanning 234,099 bp, covering 25.65% of L. aggregate mitogenome. It was worth noting that although rearrangement events occurred frequently in intergenic regions, the protein-coding regions of the plant mitochondrial genome were still highly conserved. In the mitochondrial genomes of 8 species, most protein-coding regions overlap with colinear fragments. The most extensive collinearity was detected between L. aggregate and C. camphora, containing 99.15% of the protein-coding sequences from C. camphora mitogenome, while it was relatively lower with H. nymphaeifolia (95.93%). This showed that although the PCGs of these 8 Laurales mitogenome had high similarity in sequences, the overall mitogenome structure was unstable.
3.7 Phylogenetic analysis of different orders mitochondrial genomes
We constructed phylogenetic trees based on 24 conserved PCGs (atp1/4/6/8/9, ccmB/C/FN, cox1/2/3, cob, mttB, nad1/2/3/4/4L/5/6/7/9, rps12 and sdh4) of 30 species from 9 different orders (Amborellales, Nymphaeales, Apiales, Acorales, Poales, Chloranthales, Piperales, Magnoliales, and Laurales) in angiosperms. The topologies of maximum likelihood (ML) and Bayesian inference (BI) phylogenetic trees based on mitogenome were basically the same, which could closely reflect the taxonomic relationship between these species (Figures 8A, B). It was worth noting that the mitogenome of 8 Laurales species formed a separate branch and together with Magnoliales and Piperales to form Magnoliids. The overall structure of the phylogenetic trees based on mitogenome was completely consistent with the latest classification based on the Angiosperm Phylogeny Group (APG IV). Our results showed that L. aggregata belonged to Laurales and represented a relatively new differentiation. In addition, L. aggregata was closely related to M. pauhoi and C. chago, which was consistent with the results of the haplotype network mentioned above.
Figure 8
4 Discussion
Most mitochondrial genomes in higher plants have complex circular or linear structures due to frequent recombination mediated by repeat sequences, which can promote the formation of large/small subrings and isomers in plants (; Smith and Keeling, 2015). Recent studies have revealed great differences in the structure of mitogenome among different species. For example, the mitogenome of Angelica biserrata consists of six independent circular chromosomes, while the mitogenome of Mentha spicata consists of one linear chromosome and two circular chromosomes, each of which varies in size (Jiang et al., 2023; Wang et al., 2024). In this study, we sequenced, assembled and reported the complete mitochondrial genome of the important Chinese herbal medicine Wu-Yao for the first time. Its mitogenome consists of a master circle (868,093 bp) and a linear fragment (44,380 bp) with a total length of 912,473 bp. A total of 64 unique functional genes were identified in the mitogenome, including 40 PCGs, 21 tRNA genes and 3 rRNA genes. It was worth noting that the GC content of PCGs played an important role in determining amino acid composition in the evolution of terrestrial plants, and the GC content of the mitogenome in most higher plants was between 43% and 45% (Kubo and Mikami, 2007). The GC content of the PCGs in the mitogenome of L. aggregata was 44.10%, which was basically consistent with the GC content of the mitogenome from 7 related plants (44.08%-44.21%). The Ka/Ks analysis is an important evaluation method to understand species adaptability, selection factors, evaluate adaptive potential and study heredity and evolution. It is also an important evaluation index to understand how plant genes respond to environmental stress in the process of evolution (Xie et al., 2019). In our comparison, we found that most of the 40 common PCGs in 8 Laurales mitogenome were under negative selection, only a few genes (ccmB, ccmFC, rps10, rps11 and rps7) showed signs of positive selection in all species, while some genes (atp6, ccmFN, cob, nad3, nad4, nad5, rps2 and rps3) were affected by positive selection in individual species. The atp6, ccmB, ccmFC, ccmFN, cob, nad3, nad4 and nad5 genes may play an important role in redox reaction and other adaptation processes. The relatively high Ka/Ks ratio values of rps10 and ccmB genes indicated that they played an important role in species evolution and may be closely related to environmental stress.
Codons play an important role in identifying and transmitting genetic information of organisms and connecting proteins and DNA, and also play an important role in genetic and variation of organisms (; Xiao et al., 2024). Due to the influence of gene mutation and natural selection, some synonymous codons are often used frequently in the process of protein translation, which leads to codon usage preference. This preference for specific synonymous codons plays a vital role in shaping the genetic characteristics of these organisms (Hershberg and Petrov, 2008; Parvathy et al., 2022). In this study, we used 40 PCGs to evaluate the GC content and RSCU of mitogenome in 8 species from Laurales. The results showed that all the PCGs of these species had a preference for A/T base and A/T ending codon. This result is consistent with the previous studies on the mtNDA of Angelica and Mangifera species (Niu et al., 2022; Wang et al., 2024), indicating that there are some similarities in the codon usage bias of the mitogenome in different species. Moreover, in our previous study, we also found that the third position of the chloroplast genome had a similar trend to the mitogenome in codon usage pattern. This further highlights the common codon usage bias between mitogenome and cpDNA, and emphasizes the high similarity in codon usage preference between the two genetic components (Zhang et al., 2018; Wang et al., 2020; Tang et al., 2022). ENC-plot analysis and RSCU results showed that there was a preference for codon usage in L. aggregata mitogenome, and there was a strong correlation between codon usage and gene expression level. It is more affected by natural selection, which is consistent with previous findings in Chinese angelica and soybean plants (Tang et al., 2020; Wang et al., 2024).
The collinear analysis of homologous fragments among 8 species of Laurales showed that the mitogenome of L. aggregata and its related species had experienced a large number of rearrangement events. The discovery of these homologous collinear fragments has contributed to the evolution of the genome and the regulation of gene expression. In addition, the proportion and arrangement of repetitive fragments in plant mitogenome play an important role in promoting species evolution, and are also the dominant factors in the size of mitochondrial genome (Pfeifer et al., 2013; Li et al., 2022). A total of 303 SSRs,71 tandem repeats and 1598 dispersed repeats were detected in the mitogenome of L. aggregata. These repetitive fragments constitute the intricate branch structure of the mitochondrial genome of L. aggregata. It had been found that the gene transfer of homologous fragments in the mitochondria and plastid organelles of a single species and the horizontal transfer between the same organelles of different species highlighted the intercorrelation and flowing genetic structure of these organelles, which was a common phenomenon in higher plants (). For example, there is a similar trend of pseudogenes loss or alteration in the chloroplast genome of the Cistanche genus, affecting five genes related to photosynthesis and energy production (Miao et al., 2022). In our study, 38 fragments migrated from chloroplasts to mitochondria were successfully identified, and their lengths ranged from 32 bp to 1721 bp, accounting for 1.32% of the total mitochondrial genome length. Further annotation of these homologous fragments revealed 15 complete genes, including 2 PCGs (petN (2)) and 13 tRNA genes (trnA-UGC, trnD-GUC, trnE-UUC, trnH-GUG, trnI-CAU, trnM-CAU, trnN-GUU (2), trnP-UGG, trnR-UCU, trnT-GGU, trnW-CCA, trnY-GUA). It was worth noting that the petN gene in the cpDNA taken part in the formation of subunits within cytochrome b/f complex, while 13 tRNA genes might have been lost or experienced pseudogene changes in the cpDNA.
RNA editing is a common post-transcriptional phenomenon, which widely exists in the organelles of angiosperms, including insertion, deletion and substitution, leading to changes in genetic information (Sun et al., 2016). Studies have shown that RNA editing is mainly composed of base substitutions mediated by deaminase, in which cytosine substitutes for uracil (U-to-C) have the highest frequency, and a few of them are uracil substituted cytosine (C-to-U) and hypoxanthine substituted adenine (A-to-I) (Zhu et al., 2022). The identification of RNA editing is very important for the study of mitogenome in L. aggregata. Because RNA editing can change amino acids, it often increases the overall conservatism of amino acids and changes their physical and chemical properties, thus affecting the function of proteins. Previous studies have shown that the presence of hydrophilic amino acids can promote protein folding, which is related to the decrease of the overall stability of protein structure (; ). A total of 770 RNA editing sites were identified in the mitogenome of L. aggregata, of which the number of editing event of nad4 gene was the most (68 sites), while there were generally few RNA editing sites for genes related to ribosomal proteins. Furthermore, most of the amino acids in the mitogenome were transformed from hydrophilic amino acids to hydrophobic amino acids, resulting in the increase of protein hydrophobicity and overall structural stability. In order to further understand the close relationship between Laurales plants and other angiosperms, we used 24 shared PCGs of 30 species from 9 orders to construct phylogenetic trees. The results showed that it was consistent with the phylogenetic tree of APG IV and accurately showed the taxonomic relationship of angiosperms. It was worth noting that matR gene was missing from the mitogenome of Lauraceae species, while matR gene was found in its closest relatives, Hernandiaceae species. This shows that matR gene is lost after the differentiation of Lauraceae and Hernandiaceae, and can be used as an important marker to distinguish Lauraceae species from other species, and it may be great significance to the adaptive evolution of Lauraceae species.
5 Conclusion
In this study, we completed the assembly and annotation of the mitogenome of L. aggregata. The mitogenome showed a branch structure of 912,473 bp in length. Comparative analysis showed that the mitogenome of L. aggregata contained abundant repeated sequences and HGT events, which might be one of the important reasons for the difference in the size of the mitogenome of L. aggregata. Furthermore, the mitogenome of Laurales species had experienced frequent rearrangement events, and the non-coding regions showed poor collinearity, but the structure of the coding regions was very conservative. Phylogenetic analysis based on 30 mitochondrial PCGs showed that Magnoliids was sister group of the monocots and eudicots taxa, and the evolutionary position of Lindera species in Magnoliids. This study will lay a foundation for the application of population genetics and evolution of Lindera and other species from Magnoliids.
Statements
Data availability statement
The accession numbers of mitochondrial genome of L. aggregata in Gene Bank are PP848112 and PP848113. The raw HiFi sequencing data have been deposited in NCBI with accession number: SRR30136615.
Author contributions
YS: Writing – review & editing, Writing – original draft, Visualization, Validation, Supervision, Software, Resources, Project administration, Methodology, Investigation, Formal analysis, Data curation, Conceptualization. ZC: Data curation, Funding acquisition, Methodology, Writing – review & editing, Supervision, Software, Resources. JJ: Writing – review & editing, Visualization, Supervision, Software, Resources, Investigation, Formal analysis. WW: Writing – review & editing, Supervision, Software, Methodology, Formal analysis. WY: Writing – review & editing, Resources, Investigation, Funding acquisition. SZ: Formal analysis, Investigation, Writing – review & editing. WZ: Writing – review & editing, Visualization, Validation, Supervision, Software, Resources, Project administration, Methodology, Investigation, Funding acquisition, Formal analysis, Data curation, Conceptualization.
Funding
The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. This work was funded by the Basic Public Welfare Research Project of Zhejiang Province (LGN22C020001), Startup Funding of Taizhou University for the Biomass Polysaccharide Metabolism Institute (T20231801002) and Key Scientific and Technological Grant of Zhejiang for Breeding New Agricultural Varieties (2021C02074).
Conflict of interest
Author WY was employed by the company Zhejiang Hongshiliang Group Tiantai Mountain Wu-Yao Co., Ltd.
The remaining 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.2024.1439245/full#supplementary-material
References
1
AndargieM.CongyiZ. (2022). Genome-wide analysis of codon usage in sesame (Sesamum indicum L.). Heliyon8, e08687. doi: 10.1016/j.heliyon.2021.e08687
2
Arrieta-MontielM. P.MackenzieS. A. (2011). Plant mitochondrial genomes and recombination. Plant mitochondria1, 65–82. doi: 10.1007/978-0-387-89781-3
3
BeierS.ThielT.MünchT.ScholzU.MascherM. (2017). MISA-web: a web server for microsatellite prediction. Bioinformatics33, 2583–2585. doi: 10.1093/bioinformatics/btx198
4
BensonG. (1999). Tandem repeats finder: a program to analyze DNA sequences. Nucleic Acids Res.27, 573–580. doi: 10.1093/nar/27.2.573
5
BiC.PatersonA. H.WangX.XuY.WuD.QuY.et al. (2019). Analysis of the complete mitochondrial genome sequence of the diploid cotton Gossypium raimondii by comparative genomics approaches. BioMed. Res. Int.2019, 9691253. doi: 10.1155/2019/9691253
6
BiC.ShenF.HanF.QuY.HouJ.XuK.et al. (2024a). PMAT: an efficient plant mitogenome assembly toolkit using low-coverage HiFi sequencing data. Horticul Res.11, uhae023. doi: 10.1093/hr/uhae023
7
BiC.SunN.HanF.XuK.YangY.FergusonD. K. (2024b). The first mitogenome of Lauraceae (Cinnamomum chekiangense). Plant Diversity46, 144. doi: 10.1016/j.pld.2023.11.001
8
CaoY.XuanB.PengB.LiC.ChaiX.TuP. (2015). The genus Lindera: a source of structurally diverse molecules having pharmacological significance. Phytochem. Rev.15, 869–906. doi: 10.1007/s11101-015-9432-2
9
Capella-GutiérrezS.Silla-MartínezJ. M.GabaldónT. (2009). trimAl: a tool for automated alignment trimming in large-scale phylogenetic analyses. Bioinformatics25, 1972–1973. doi: 10.1093/bioinformatics/btp348
10
ChanP. p.LinB. y.MakA. j.LoweT. m. (2021). tRNAscan-SE 2.0: improved detection and functional classification of transfer RNA genes. Nucleic Acids Res.49, 9077–9096. doi: 10.1093/nar/gkab688
11
ChenC.ChenH.ZhangY.ThomasH. R.FrankM. H.HeY.et al. (2020). TBtools: An Integrative toolkit developed for interactive analyses of big biological data. Mol. Plant13, 1194–1202. doi: 10.1016/j.molp.2020.06.009
12
ChenY.YeW.ZhangY.XuY. (2015). High speed BLASTN: an accelerated MegaBLAST search tool. Nucleic Acids Res.43, 7762–7768. doi: 10.1093/nar/gkv784
13
CovelloP. S.GrayM. W. (1989). RNA editing in plant mitochondria. Nature341, 662–666. doi: 10.1038/341662a0
14
DarlingA. C.MauB.BlattnerF. R.PernaN. T. (2004). Mauve: multiple alignment of conserved genomic sequence with rearrangements. Genome Res.14, 1394–1403. doi: 10.1101/gr.2289704
15
DarribaD.TaboadaG. L.DoalloR.PosadaD. (2012). jModelTest 2: more models, new heuristics and high-performance computing. Nat. Methods9, 772. doi: 10.1038/nmeth.2109
16
DoyleJ. J.DoyleJ. L. (1987). A rapid DNA isolation procedure for small quantities of fresh leaf tissue. Phytochem. Bull., 19, 11–15.
17
DuanY.WangJ.KangF.LiJ. (2018). Research of molecular evolution characteristics of nrDNA ITS and cpDNA trnH-psbA sequences in Prunus mongolica. Genomics Appl. Biol.37, 3035–3041. doi: 10.13417/j.gab.037.003035
18
EderaA. A.SmallI.MiloneD. H.Sanchez-PuertaM. V. (2021). Deepred-Mt: Deep representation learning for predicting C-to-U RNA editing in plant mitochondria. Comput. Biol. Med.136, 104682. doi: 10.1016/j.compbiomed.2021.104682
19
FujiiS.KazamaT.YamadaM.ToriyamaK. (2010). Discovery of global genomic re-organization based on comparison of two newly sequenced rice mitochondrial genomes with cytoplasmic male sterility-related genes. BMC Genomics11, 1–15. doi: 10.1186/1471-2164-11-209
20
GandiniC. L.GarciaL. E.AbbonaC. C.Sanchez-PuertaM. V. (2019). The complete organelle genomes of Physochlaina orientalis: Insights into short sequence repeats across seed plant mitochondrial genomes. Mol. Phylogenet. Evol.137, 274–284. doi: 10.1016/j.ympev.2019.05.012
21
GrayM. W. (2009). RNA editing in plant mitochondria: 20 years later. IUBMB Life61, 1101–1104. doi: 10.1002/iub.272
22
GrayM. W.CovelloP. S. (1993). RNA editing in plant mitochondria and chloroplasts. FASEB J.7, 64–71. doi: 10.1096/fasebj.7.1.8422976
23
GreweF.HerresS.ViehöverP.PolsakiewiczM.WeisshaarB.KnoopV. (2011). A unique transcriptome: 1782 positions of RNA editing alter 1406 codon identities in mitochondrial mRNAs of the lycophyte Isoetes engelmannii. Nucleic Acids Res.39, 2890–2902. doi: 10.1093/nar/gkq1227
24
GualbertoJ. M.MileshinaD.WalletC.NiaziA. K.Weber-LotfiF.DietrichA. (2014). The plant mitochondrial genome: Dynamics and maintenance. Biochimie100, 107–120. doi: 10.1016/j.biochi.2013.09.016
25
HanF.BiC.ZhaoY.GaoM.WangY.ChenY. (2024). Unraveling the complex evolutionary features of the Cinnamomum camphora mitochondrial genome. Plant Cell Rep.43, 183. doi: 10.1007/s00299-024-03256-1
26
HanH.XuB.AminA.LiH.YuX.GongM.et al. (2019). Quercetin−3−O−α−L−rhamnopyranoside derived from the leaves of Lindera aggregata (Sims) Kosterm. evokes the autophagy−induced nuclear factor erythroid 2−related factor 2 antioxidant pathway in human umbilical vein endothelial cells. Int. J. Mol. Med.43, 461–474. doi: 10.3892/ijmm.2018.3976
27
HebertP. D. N.CywinskaA.BallS. L.DewaardJ. R. (2003). Biological identifications through DNA barcodes. Proc. R. Soc. London. Ser. B: Biol. Sci.270, 313–321. doi: 10.1098/rspb.2002.2218
28
HershbergR.PetrovD. A. (2008). Selection on codon bias. Annu. Rev. Genet.42, 287–299. doi: 10.1146/annurev.genet.42.110807.091442
29
Hung-PinT. (1987). A study on the system of lindera. J. Syst. Evol.25, 161–171.
30
JiangM.NiY.ZhangJ.LiJ.LiuC. (2023). Complete mitochondrial genome of Mentha spicata L. reveals multiple chromosomal configurations and RNA editing events. Int. J. Biol. Macromol251, 126257. doi: 10.1016/j.ijbiomac.2023.126257
31
KalyaanamoorthyS.MinhB. Q.WongT. K.Von HaeselerA.JermiinL. S. (2017). ModelFinder: fast model selection for accurate phylogenetic estimates. Nat. Methods 14, 587–589. doi: 10.1038/nmeth.4285
32
KatohK.StandleyD. M. (2013). MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Mol. Biol. Evol.30, 772–780. doi: 10.1093/molbev/mst010
33
KuboT.MikamiT. (2007). Organization and variation of angiosperm mitochondrial genome. Physiol Plantarum129, 6–13. doi: 10.1111/j.1399-3054.2006.00768.x
34
KuboT.NewtonK. J. (2008). Angiosperm mitochondrial genomes and mutations. Mitochondrion8, 5–14. doi: 10.1016/j.mito.2007.10.006
35
KurtzS.ChoudhuriJ. V.OhlebuschE.SchleiermacherC.StoyeJ.GiegerichR. (2001). REPuter: the manifold applications of repeat analysis on a genomic scale. Nucleic Acids Res.29, 4633–4642. doi: 10.1093/nar/29.22.4633
36
LeighJ. W.BryantD. (2015). popart: full-feature software for haplotype network construction. Methods Ecol. Evol.6, 1110–1116. doi: 10.1111/2041-210X.12410
37
LiH. (2021). New strategies to improve minimap2 alignment accuracy. Bioinformatics37, 4572–4574. doi: 10.1093/bioinformatics/btab705
38
LiJ.LiJ.MaY.KouL.WeiJ.WangW. (2022). The complete mitochondrial genome of okra (Abelmoschus esculentus): Using nanopore long reads to investigate gene transfer from chloroplast genomes and rearrangements of mitochondrial DNA molecules. BMC Genomics23, 481. doi: 10.1186/s12864-022-08706-2
39
LuJ.ZhaoZ.NiL.GaaweD. (2016). Identification of medicinal plants based on mitochondrial DNA sequences. Chin. Traditional Herbal Drugs47, 1791–1796. doi: 10.7501/j.issn.0253-2670.2016.10.027
40
LvY.ZouY.ZhangX.LiuB.PengX.ChuC. (2023). A review on the chemical constituents and pharmacological efficacies of Lindera aggregata (Sims) Kosterm. Front. Nutr.9, 1071276. doi: 10.3389/fnut.2022.1071276
41
MiaoY.ChenH.XuW.LiuC.HuangL. (2022). Cistanche Species mitogenomes suggest diversity and complexity in lamiales-order mitogenomes. Genes13, 1791. doi: 10.3390/genes13101791
42
MillarA. H.HeazlewoodJ. L.KristensenB. K.BraunH.-P.MøllerI. M. (2005). The plant mitochondrial proteome. Trends Plant Sci.10, 36–43. doi: 10.1016/j.tplants.2004.12.002
43
MinhB. Q.SchmidtH. A.ChernomorO.SchrempfD.WoodhamsM. D.Von HaeselerA.et al. (2020). IQ-TREE 2: new models and efficient methods for phylogenetic inference in the genomic era. Mol. Biol. Evol.37, 1530–1534. doi: 10.1093/molbev/msaa015
44
NiuY.GaoC.LiuJ. (2022). Complete mitochondrial genomes of three Mangifera species, their genomic structure and gene transfer from chloroplast genomes. BMC Genomics23, 147. doi: 10.1186/s12864-022-08383-1
45
OgiharaY.YamazakiY.MuraiK.KannoA.TerachiT.ShiinaT.et al. (2005). Structural dynamics of cereal mitochondrial genomes as revealed by complete nucleotide sequencing of the wheat mitochondrial genome. Nucleic Acids Res.33, 6235–6250. doi: 10.1093/nar/gki925
46
ParvathyS. T.UdayasuriyanV.BhadanaV. (2022). Codon usage bias. Mol. Biol. Rep.49, 539–565. doi: 10.1007/s11033-021-06749-4
47
PfeiferM.MartisM.AspT.MayerK. F.LübberstedtT.ByrneS.et al. (2013). The perennial ryegrass GenomeZipper: targeted use of genome resources for comparative grass genomics. Plant Physiol.161, 571–582. doi: 10.1104/pp.112.207282
48
RomeroH.ZavalaA.MustoH. (2000). Codon usage in Chlamydia trachomatis is the result of strand-specific mutational biases and a complex pattern of selective forces. Nucleic Acids Res.28, 2084–2090. doi: 10.1093/nar/28.10.2084
49
RonquistF.TeslenkoM.van der MarkP.AyresD. L.DarlingA.HöhnaS.et al. (2012). MrBayes 3.2: efficient Bayesian phylogenetic inference and model choice across a large model space. Systemat Biol.61, 539–542. doi: 10.1093/sysbio/sys029
50
RozasJ.Sánchez-DelbarrioJ. C.MesseguerX.RozasR. (2003). DnaSP, DNA polymorphism analyses by the coalescent and other methods. Bioinformatics19, 2496–2497. doi: 10.1093/bioinformatics/btg359
51
SallehW.M.N.H.W. (2020). Lindera aggregata (Sims) Kosterm: Review on phytochemistry and biological activities. Boletin Latinoamericano y del Caribe Plantas Medicinales y Aromaticas19, 527–541. doi: 10.37360/blacpma
52
ShiY.ChenZ.JiangJ.LiX.ZengW. (2024). Comparative analysis of chloroplast genomes of “Tiantai wu-yao”(Lindera aggregata) and taxa of the same genus and different genera. Genes15, 263. doi: 10.3390/genes15030263
53
SkippingtonE.BarkmanT. J.RiceD. W.PalmerJ. D. (2015). Miniaturized mitogenome of the parasitic plant Viscum scurruloideum is extremely divergent and dynamic and has lost all nad genes. Proc. Natl. Acad. Sci.112, E3515–E3524. doi: 10.1073/pnas.1504491112
54
SloanD. B.AlversonA. J.ChuckalovcakJ. P.WuM.MccauleyD. E.PalmerJ. D.et al. (2012). Rapid evolution of enormous, multichromosomal genomes in flowering plant mitochondria with exceptionally high mutation rates. PloS Biol.10, e1001241. doi: 10.1371/journal.pbio.1001241
55
SmithD. R.KeelingP. J. (2015). Mitochondrial and plastid genome architecture: reoccurring themes, but significant differences at the extremes. Proc. Natl. Acad. Sci.112, 10177–10184. doi: 10.1073/pnas.1422049112
56
SueokaN. (1988). Directional mutation pressure and neutral molecular evolution. Proc. Natl. Acad. Sci.85, 2653–2657. doi: 10.1073/pnas.85.8.2653
57
SunT.BentolilaS.HansonM. R. (2016). The unexpected diversity of plant organelle RNA editosomes. Trends Plant Sci.21, 962–973. doi: 10.1016/j.tplants.2016.07.005
58
TangD.WeiF.QuanC.HuangS.HuangY.WeiK.et al. (2022). Codon usage bias and evolution analysis in the mitochondrial genome of Mesona chinensis Benth. Acta Physiologiae Plantarum44, 118. doi: 10.1007/s11738-022-03453-z
59
TangX.-M.YangS.-Z.ChenH.-Z.SunZ.-D.LaiZ.-G.ZengW.-Y.et al. (2020). Comparative analysis on codon usage Bias inmitogenome of two species in genus Glycine. Guihaia.40, 926–934. doi: 10.11931/guihaia.gxzw201911033
60
TaoY.DengY.WangP. (2024). Traditional uses, phytochemistry, pharmacology, processing methods and quality control of Lindera aggregata (Sims) Kosterm: A critical review. J. Ethnopharmacology318, 116954. doi: 10.1016/j.jep.2023.116954
61
WangL.LiuX.XuY.ZhangZ.WeiY.HuY.et al. (2024). Assembly and comparative analysis of the first complete mitochondrial genome of a traditional Chinese medicine Angelica biserrata (Shan et Yuan) Yuan et Shan. Int. J. Biol. Macromol257, 128571. doi: 10.1016/j.ijbiomac.2023.128571
62
WangH.SunH.KwonW.-S.JinH.YangD.-C. (2009). Molecular identification of the Korean ginseng cultivar “Chunpoong” using the mitochondrial nad7 intron 4 region. Mitochondrial DNA20, 41–45. doi: 10.1080/19401730902856738
63
WangZ.XuB.LiB.ZhouQ.WangG.JiangX.et al. (2020). Comparative analysis of codon usage patterns in chloroplast genomes of six Euphorbiaceae species. PeerJ8, e8251. doi: 10.7717/peerj.8251
64
WangD.ZhangY.ZhangZ.ZhuJ.YuJ. (2010). KaKs_Calculator 2.0: a toolkit incorporating gamma-series methods and sliding window strategies. Genomics Proteomics Bioinf.8, 77–80. doi: 10.1016/S1672-0229(10)60008-3
65
WengM.YouS.LuoJ.LinZ.ChenT.PengX.et al. (2022). Antibacterial mechanism of polysaccharides from the leaves of Lindera aggregata (Sims) Kosterm. by metabolomics based on HPLC/MS. Int. J. Biol. Macromol221, 303–313. doi: 10.1016/j.ijbiomac.2022.09.009
66
WickR. R.SchultzM. B.ZobelJ.HoltK. E. (2015). Bandage: interactive visualization of de novo genome assemblies. Bioinformatics31, 3350–3352. doi: 10.1093/bioinformatics/btv383
67
WrightF. (1990). The ‘effective number of codons’ used in a gene. Gene87, 23–29. doi: 10.1016/0378-1119(90)90491-9
68
XiaoM.HuX.LiY.LiuQ.ShenS.JiangT.et al. (2024). Comparative analysis of codon usage patterns in the chloroplast genomes of nine forage legumes. Physiol. Mol. Biol. Plants30, 153–166. doi: 10.1007/s12298-024-01421-0
69
XieJ.ChenY.CaiG.CaiR.HuZ.WangH. (2023). Tree Visualization By One Table (tvBOT): a web application for visualizing, modifying and annotating phylogenetic trees. Nucleic Acids Res., 51, W587–W592. doi: 10.1093/nar/gkad359
70
XieD.-F.YuH.-X.PriceM.XieC.DengY.-Q.ChenJ.-P.et al. (2019). Phylogeny of Chinese Allium species in section Daghestanica and adaptive evolution of Allium (Amaryllidaceae, Allioideae) species revealed by the chloroplast complete genome. Front. Plant Sci.10, 460. doi: 10.3389/fpls.2019.00460
71
YangH.NiY.ZhangX.LiJ.ChenH.LiuC. (2023). The mitochondrial genomes of Panax notoginseng reveal recombination mediated by repeats associated with DNA replication. Int. J. Biol. Macromol252, 126359. doi: 10.1016/j.ijbiomac.2023.126359
72
ZhangX.ChenH.NiY.WuB.LiJ.BurzyńskiA.et al. (2024). Plant mitochondrial genome map (PMGmap): A software tool for the comprehensive visualization of coding, noncoding and genome features of plant mitochondrial genomes. Mol. Ecol. Resour.24, e13952. doi: 10.1111/1755-0998.13952
73
ZhangD.GaoF.JakovlićI.ZouH.ZhangJ.LiW. X.et al. (2020). PhyloSuite: An integrated and scalable desktop platform for streamlined molecular sequence data management and evolutionary phylogenetics studies. Mol. Ecol. Resour.20, 348–355. doi: 10.1111/1755-0998.13096
74
ZhangR.ZhangL.WangW.ZhangZ.DuH.QuZ.et al. (2018). Differences in codon usage bias between photosynthesis-related genes and genetic system-related genes of chloroplast genomes in cultivated and wild solanum species. Int. J. Mol. Sci.19, 3142. doi: 10.3390/ijms19103142
75
ZhaoM.-L.SongY.NiJ.YaoX.TanY.-H.XuZ.-F. (2018). Comparative chloroplast genomics and phylogenetics of nine Lindera species (Lauraceae). Sci. Rep.8, 8844. doi: 10.1038/s41598-018-27090-0
76
ZhuL.XianF.-J.ZhangQ.-N.HuJ. (2022). Research progress of RNA editing. Biotechnol. Bull.38, 1, 1–14. doi: 10.13560/j.cnki.biotech.bull.1985.2021-1326
Summary
Keywords
Lindera aggregata, mitochondrial genome, multi-branched conformation, homologous fragments, phylogeny
Citation
Shi Y, Chen Z, Jiang J, Wu W, Yu W, Zhang S and Zeng W (2024) The assembly and comparative analysis of the first complete mitogenome of Lindera aggregata. Front. Plant Sci. 15:1439245. doi: 10.3389/fpls.2024.1439245
Received
27 May 2024
Accepted
16 August 2024
Published
03 September 2024
Volume
15 - 2024
Edited by
Wenpan Dong, Beijing Forestry University, China
Reviewed by
Changwei Bi, Nanjing Forestry University, China
Fang He, Sichuan Agricultural University, China
Updates
Copyright
© 2024 Shi, Chen, Jiang, Wu, Yu, Zhang and Zeng.
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: Wei Zeng, zengw@tzc.edu.cn
Disclaimer
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.