<?xml version="1.0" encoding="UTF-8"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing DTD v2.3 20070202//EN" "journalpublishing.dtd">
<article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" article-type="research-article" dtd-version="2.3" xml:lang="EN">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Front. Plant Sci.</journal-id>
<journal-title>Frontiers in Plant Science</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Plant Sci.</abbrev-journal-title>
<issn pub-type="epub">1664-462X</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="doi">10.3389/fpls.2024.1439245</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Plant Science</subject>
<subj-group>
<subject>Original Research</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>The assembly and comparative analysis of the first complete mitogenome of <italic>Lindera aggregata</italic>
</article-title>
</title-group>
<contrib-group>
<contrib contrib-type="author">
<name>
<surname>Shi</surname>
<given-names>Yujie</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/2826759"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-original-draft/"/>
<role content-type="https://credit.niso.org/contributor-roles/visualization/"/>
<role content-type="https://credit.niso.org/contributor-roles/validation/"/>
<role content-type="https://credit.niso.org/contributor-roles/supervision/"/>
<role content-type="https://credit.niso.org/contributor-roles/software/"/>
<role content-type="https://credit.niso.org/contributor-roles/resources/"/>
<role content-type="https://credit.niso.org/contributor-roles/project-administration/"/>
<role content-type="https://credit.niso.org/contributor-roles/methodology/"/>
<role content-type="https://credit.niso.org/contributor-roles/investigation/"/>
<role content-type="https://credit.niso.org/contributor-roles/formal-analysis/"/>
<role content-type="https://credit.niso.org/contributor-roles/data-curation/"/>
<role content-type="https://credit.niso.org/contributor-roles/conceptualization/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Chen</surname>
<given-names>Zhen</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/562619"/>
<role content-type="https://credit.niso.org/contributor-roles/data-curation/"/>
<role content-type="https://credit.niso.org/contributor-roles/funding-acquisition/"/>
<role content-type="https://credit.niso.org/contributor-roles/methodology/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
<role content-type="https://credit.niso.org/contributor-roles/supervision/"/>
<role content-type="https://credit.niso.org/contributor-roles/software/"/>
<role content-type="https://credit.niso.org/contributor-roles/resources/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Jiang</surname>
<given-names>Jingyong</given-names>
</name>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
<role content-type="https://credit.niso.org/contributor-roles/visualization/"/>
<role content-type="https://credit.niso.org/contributor-roles/supervision/"/>
<role content-type="https://credit.niso.org/contributor-roles/software/"/>
<role content-type="https://credit.niso.org/contributor-roles/resources/"/>
<role content-type="https://credit.niso.org/contributor-roles/investigation/"/>
<role content-type="https://credit.niso.org/contributor-roles/formal-analysis/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Wu</surname>
<given-names>Wenwu</given-names>
</name>
<xref ref-type="aff" rid="aff3">
<sup>3</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/1185151"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
<role content-type="https://credit.niso.org/contributor-roles/supervision/"/>
<role content-type="https://credit.niso.org/contributor-roles/software/"/>
<role content-type="https://credit.niso.org/contributor-roles/methodology/"/>
<role content-type="https://credit.niso.org/contributor-roles/formal-analysis/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Yu</surname>
<given-names>Weifu</given-names>
</name>
<xref ref-type="aff" rid="aff4">
<sup>4</sup>
</xref>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
<role content-type="https://credit.niso.org/contributor-roles/resources/"/>
<role content-type="https://credit.niso.org/contributor-roles/investigation/"/>
<role content-type="https://credit.niso.org/contributor-roles/funding-acquisition/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Zhang</surname>
<given-names>Shumeng</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/2752126"/>
<role content-type="https://credit.niso.org/contributor-roles/formal-analysis/"/>
<role content-type="https://credit.niso.org/contributor-roles/investigation/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
</contrib>
<contrib contrib-type="author" corresp="yes">
<name>
<surname>Zeng</surname>
<given-names>Wei</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<xref ref-type="author-notes" rid="fn001">
<sup>*</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/1330740"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
<role content-type="https://credit.niso.org/contributor-roles/visualization/"/>
<role content-type="https://credit.niso.org/contributor-roles/validation/"/>
<role content-type="https://credit.niso.org/contributor-roles/supervision/"/>
<role content-type="https://credit.niso.org/contributor-roles/software/"/>
<role content-type="https://credit.niso.org/contributor-roles/resources/"/>
<role content-type="https://credit.niso.org/contributor-roles/project-administration/"/>
<role content-type="https://credit.niso.org/contributor-roles/methodology/"/>
<role content-type="https://credit.niso.org/contributor-roles/investigation/"/>
<role content-type="https://credit.niso.org/contributor-roles/funding-acquisition/"/>
<role content-type="https://credit.niso.org/contributor-roles/formal-analysis/"/>
<role content-type="https://credit.niso.org/contributor-roles/data-curation/"/>
<role content-type="https://credit.niso.org/contributor-roles/conceptualization/"/>
</contrib>
</contrib-group>
<aff id="aff1">
<sup>1</sup>
<institution>Zhejiang Provincial Key Laboratory of Plant Evolutionary Ecology and Conservation, College of Life Sciences, Taizhou University</institution>, <addr-line>Taizhou</addr-line>, <country>China</country>
</aff>
<aff id="aff2">
<sup>2</sup>
<institution>Institute of Horticulture, Taizhou Academy of Agricultural Sciences</institution>, <addr-line>Linhai</addr-line>, <country>China</country>
</aff>
<aff id="aff3">
<sup>3</sup>
<institution>State Key Laboratory of Subtropical Silviculture, College of Forestry and Biotechnology, Zhejiang Agricultural and Forestry (A&amp;F) University</institution>, <addr-line>Hangzhou</addr-line>, <country>China</country>
</aff>
<aff id="aff4">
<sup>4</sup>
<institution>Zhejiang Hongshiliang Group Tiantai Mountain Wu-Yao Co., Ltd., RedRock Group</institution>, <addr-line>Taizhou</addr-line>, <country>China</country>
</aff>
<author-notes>
<fn fn-type="edited-by">
<p>Edited by: Wenpan Dong, Beijing Forestry University, China</p>
</fn>
<fn fn-type="edited-by">
<p>Reviewed by: Changwei Bi, Nanjing Forestry University, China</p>
<p>Fang He, Sichuan Agricultural University, China</p>
</fn>    <fn fn-type="corresp" id="fn001">
<p>*Correspondence: Wei Zeng, <email xlink:href="mailto:zengw@tzc.edu.cn">zengw@tzc.edu.cn</email>
</p>
</fn>
</author-notes>
<pub-date pub-type="epub">
<day>03</day>
<month>09</month>
<year>2024</year>
</pub-date>
<pub-date pub-type="collection">
<year>2024</year>
</pub-date>
<volume>15</volume>
<elocation-id>1439245</elocation-id>
<history>
<date date-type="received">
<day>27</day>
<month>05</month>
<year>2024</year>
</date>
<date date-type="accepted">
<day>16</day>
<month>08</month>
<year>2024</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#xa9; 2024 Shi, Chen, Jiang, Wu, Yu, Zhang and Zeng</copyright-statement>
<copyright-year>2024</copyright-year>
<copyright-holder>Shi, Chen, Jiang, Wu, Yu, Zhang and Zeng</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/">
<p>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.</p>
</license>
</permissions>
<abstract>
<p>
<italic>Lindera aggregata</italic>, a member belongs to the genus <italic>Lindera</italic> 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 <italic>L. aggregata</italic> 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 <italic>ccmB</italic>, <italic>ccmFC</italic>, <italic>rps10</italic>, <italic>rps11</italic> and <italic>rps7</italic> 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 <italic>Lindera</italic> species.</p>
</abstract>
<kwd-group>
<kwd>
<italic>Lindera aggregata</italic>
</kwd>
<kwd>mitochondrial genome</kwd>
<kwd>multi-branched conformation</kwd>
<kwd>homologous fragments</kwd>
<kwd>phylogeny</kwd>
</kwd-group>
<contract-num rid="cn001">LGN22C020001</contract-num>
<contract-sponsor id="cn001">Basic Public Welfare Research Program of Zhejiang Province<named-content content-type="fundref-id">10.13039/501100017577</named-content>
</contract-sponsor>
<counts>
<fig-count count="8"/>
<table-count count="0"/>
<equation-count count="0"/>
<ref-count count="76"/>
<page-count count="14"/>
<word-count count="6983"/>
</counts>
<custom-meta-wrap>
<custom-meta>
<meta-name>section-in-acceptance</meta-name>
<meta-value>Plant Systematics and Evolution</meta-value>
</custom-meta>
</custom-meta-wrap>
</article-meta>
</front>
<body>
<sec id="s1" sec-type="intro">
<label>1</label>
<title>Introduction</title>
<p>The genus <italic>Lindera</italic> 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 (<xref ref-type="bibr" rid="B29">Hung-Pin, 1987</xref>). There are more than 40 species of this genus in China, accounting for 46% of the whole genus (<xref ref-type="bibr" rid="B60">Tao et&#xa0;al., 2024</xref>). The plants of <italic>Lindera</italic> are widely used in traditional medicine and have high economic value, and <italic>Lindera aggregata</italic> (Sims) Kosterm is a representative one of them. It is also known as &#x201c;Wu-Yao&#x201d;, is a common folk medicinal plant (<xref ref-type="bibr" rid="B40">Lv et&#xa0;al., 2023</xref>).</p>
<p>The leaves of <italic>L. aggregata</italic> are edible and are often used as functional tea or dietary supplements, with health effects such as anti-liver injury and lipid-lowering (<xref ref-type="bibr" rid="B26">Han et&#xa0;al., 2019</xref>; <xref ref-type="bibr" rid="B65">Weng et&#xa0;al., 2022</xref>). The dried root of <italic>L. aggregata</italic>, 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 (<xref ref-type="bibr" rid="B8">Cao et&#xa0;al., 2015</xref>). It mainly grows in eastern, central, southern and southwestern China, such as Zhejiang, Jiangxi, Anhui and other provinces. Among them, BenCaoGangMu recorded: &#x201c;The Wu-Yao originated from Tiantai was the best&#x201d;. Therefore, Zhejiang Tiantai has been a famous genuine producing area of <italic>L. aggregata</italic> since ancient times, which is called &#x201c;Tiantai Wu-Yao&#x201d;. The main active components of <italic>L. aggregata</italic> are sesquiterpenes, alkaloids, flavonoids and volatile oils, which have the effects of protecting liver, anti-inflammation, anti-virus, antibacterial, anti-tumor and anti-oxidation (<xref ref-type="bibr" rid="B51">Salleh, 2020</xref>). Due to the increasing use of authentic <italic>L. aggregata</italic> and the lack of centralized protection, the reserves of authentic <italic>L. aggregata</italic> are decreasing day by day, so intensive management is extremely urgent. According to the Chinese Pharmacopoeia, &#x201c;the taproot, which is not spindle-shaped, cannot be used medicinally&#x201d;. Moreover, the growing year of <italic>L. aggregata</italic> 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 (<xref ref-type="bibr" rid="B27">Hebert et&#xa0;al., 2003</xref>). The primitive herbs used in traditional Chinese medicine are also the object of molecular research (<xref ref-type="bibr" rid="B39">Lu et&#xa0;al., 2016</xref>). For example, the <italic>nad7</italic> gene of the mitogenome can be used to identify the Korean ginseng cultivar &#x201c;Chunpoong&#x201d; and other ginseng (<xref ref-type="bibr" rid="B62">Wang et&#xa0;al., 2009</xref>).</p>
<p>The genus <italic>Lindera</italic> is a large and complex genus with similar transitional characters among different species. In the past, the molecular phylogeny of <italic>Lindera</italic> was mainly focused on the chloroplast genome, which revealed that <italic>Lindera</italic> was a polyphyletic group (<xref ref-type="bibr" rid="B75">Zhao et&#xa0;al., 2018</xref>; <xref ref-type="bibr" rid="B52">Shi et&#xa0;al., 2024</xref>), and its systematic relationship with related species from <italic>Litsea</italic> and <italic>Laurus</italic> is still unclear. Mitochondrial genome analysis is very important for understanding the evolution and genome structure of various plants (<xref ref-type="bibr" rid="B71">Yang et&#xa0;al., 2023</xref>). Recent research on the mitogenome of the genus <italic>Cinnamomum</italic> has further discussed the evolutionary location of the Magnoliids species (<xref ref-type="bibr" rid="B7">Bi et&#xa0;al., 2024b</xref>; <xref ref-type="bibr" rid="B25">Han et&#xa0;al., 2024</xref>), but information on the mitogenome of other genera from Magnoliids is still lacking. In particular, the mitochondrial genome of <italic>Lindera</italic> plants has not been reported so far. It is an urgent task to analyze the genomic characteristics of <italic>L. aggregata</italic>, compare and analyze the evolution of the mitogenome of <italic>Lindera</italic> and other Laurales plants.</p>
<p>Plant chloroplasts and mitochondria are two important semi-autonomous organelles except for the nucleus, which have their own independent genetic material and genetic system (<xref ref-type="bibr" rid="B24">Gualberto et&#xa0;al., 2014</xref>). 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 (<xref ref-type="bibr" rid="B42">Millar et&#xa0;al., 2005</xref>; <xref ref-type="bibr" rid="B61">Wang et&#xa0;al., 2024</xref>). 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 (<xref ref-type="bibr" rid="B54">Sloan et&#xa0;al., 2012</xref>; <xref ref-type="bibr" rid="B53">Skippington et&#xa0;al., 2015</xref>), but the coding genes are relatively conservative (<xref ref-type="bibr" rid="B45">Ogihara et&#xa0;al., 2005</xref>). 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) (<xref ref-type="bibr" rid="B20">Gandini et&#xa0;al., 2019</xref>). 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 (<xref ref-type="bibr" rid="B2">Arrieta-Montiel and Mackenzie, 2011</xref>). 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 (<xref ref-type="bibr" rid="B38">Li et&#xa0;al., 2022</xref>). 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 (<xref ref-type="bibr" rid="B34">Kubo and Newton, 2008</xref>). 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 (<xref ref-type="bibr" rid="B13">Covello and Gray, 1989</xref>; <xref ref-type="bibr" rid="B23">Grewe et&#xa0;al., 2011</xref>), while the RNA editing sites in some chloroplasts are relatively conservative in evolution (<xref ref-type="bibr" rid="B22">Gray and Covello, 1993</xref>; <xref ref-type="bibr" rid="B17">Duan et&#xa0;al., 2018</xref>).</p>
<p>In this study, the complete sequence of mitochondrial genome of <italic>L. aggregata</italic> 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 <italic>L. aggregata</italic>, 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 <italic>L. aggregata</italic> 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 <italic>Lindera</italic> and will be of great value to the classification, species identification and molecular breeding of this genus with high value.</p>
</sec>
<sec id="s2" sec-type="materials|methods">
<label>2</label>
<title>Materials and methods</title>
<sec id="s2_1">
<label>2.1</label>
<title>Sample collection, DNA extraction and sequencing</title>
<p>The sample materials were collected from the planting base of <italic>Lindera aggregata</italic> from Zhejiang Hongshiliang Group Tiantai Mountain Wu-Yao Co., Ltd., Zhejiang, China (120&#xb0;48.23&#x2032; E, 29&#xb0;12.71&#x2032; N). Fresh and tender leaves were collected and quick-frozen with liquid nitrogen. Genome DNA was extracted by improved CTAB method (<xref ref-type="bibr" rid="B16">Doyle and Doyle, 1987</xref>). 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 &gt; 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 &#x201c;min-passes = 3, min-rq = 0.99&#x201d;.</p>
</sec>
<sec id="s2_2">
<label>2.2</label>
<title>Assembly and annotation of the mitochondrial genome</title>
<p>The mitochondrial genome was assembled by PMAT v1.5.3 software (<xref ref-type="bibr" rid="B6">Bi et&#xa0;al., 2024a</xref>), and the initial mitochondrial assembly map was obtained. Bandage v0.8.1 software (<xref ref-type="bibr" rid="B66">Wick et&#xa0;al., 2015</xref>) 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 (<xref ref-type="bibr" rid="B37">Li, 2021</xref>) to help analyze the repetitive regions in the mitochondrial genome, and finally produce a complete mitogenome of <italic>L. aggregata</italic>. Intelligent Plant Mitochondrial Genome Annotator (IPMGA, <ext-link ext-link-type="uri" xlink:href="http://www.1kmpg.cn/ipmga/">http://www.1kmpg.cn/ipmga/</ext-link>) was used to annotate the mitochondrial genome. tRNA and rRNA in mitochondrial genome were annotated by tRNAscan-SE v2.0 (<xref ref-type="bibr" rid="B10">Chan et&#xa0;al., 2021</xref>) and BLASTN software (<xref ref-type="bibr" rid="B12">Chen et&#xa0;al., 2015</xref>), respectively. The mitochondrial genome was visualized by Plant Mitochondrial Genomes Map (PMGmap) online software (<xref ref-type="bibr" rid="B72">Zhang et&#xa0;al., 2024</xref>), and the <italic>cis-</italic> and <italic>trans-</italic>splicing gene maps were drawn. Finally, the mitogenome of <italic>L. aggregata</italic> was uploaded to the GenBank database to obtain the accession number.</p>
</sec>
<sec id="s2_3">
<label>2.3</label>
<title>Comparative analysis of 8 plants mitogenome from Laurales</title>
<p>We downloaded all the published mitogenome of Laurales (<italic>Cinnamomum camphora</italic>, <italic>C. chago</italic>, <italic>C. chekiangense</italic>, <italic>C. insularimontanum</italic>, <italic>Caryodaphnopsis henryi</italic>, <italic>Machilus pauhoi</italic> and <italic>Hernandia nymphaeifolia</italic>) 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 <italic>L. aggregata</italic> were extracted by PhyloSuit v1.2.3 software (<xref ref-type="bibr" rid="B73">Zhang et&#xa0;al., 2020</xref>). The multi-sequence alignment was performed among 40 common PCGs by MAFFT v7.149b software (<xref ref-type="bibr" rid="B32">Katoh and Standley, 2013</xref>), and the sequences were trimmed by trimAl v1.2 software (<xref ref-type="bibr" rid="B9">Capella-Guti&#xe9;rrez et&#xa0;al., 2009</xref>). The trimmed sequences were merged, and the maximum likelihood (ML) phylogenetic tree was constructed by IQ-TREE v2.2.0.3 software (<xref ref-type="bibr" rid="B43">Minh et&#xa0;al., 2020</xref>) with default parameters. Furthermore, through the Genepioneer Cloud (<ext-link ext-link-type="uri" xlink:href="http://cloud.genepioneer.com:9929/#/home">http://cloud.genepioneer.com:9929/#/home</ext-link>), 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 (<xref ref-type="bibr" rid="B50">Rozas et&#xa0;al., 2003</xref>) to generate NEXUS haplotype data files. Haplotype networks of 40 conserved protein coding sequences were constructed via PopART v4.8.4 software (<xref ref-type="bibr" rid="B36">Leigh and Bryant, 2015</xref>). 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 (<xref ref-type="bibr" rid="B64">Wang et&#xa0;al., 2010</xref>).</p>
</sec>
<sec id="s2_4">
<label>2.4</label>
<title>Analysis of codon usage bias</title>
<p>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 (<xref ref-type="bibr" rid="B73">Zhang et&#xa0;al., 2020</xref>), and visualized by TBtools (<xref ref-type="bibr" rid="B11">Chen et&#xa0;al., 2020</xref>). 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 (<ext-link ext-link-type="uri" xlink:href="http://cloud.genepioneer.com:9929/#/tool">http://cloud.genepioneer.com:9929/#/tool</ext-link>) 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 (<xref ref-type="bibr" rid="B67">Wright, 1990</xref>). ENC-plot analysis can explore the relationship between ENC and GC<sub>3</sub> distribution, which was an effective way to visualize the codon usage bias of genetic data. The standard curve (ENC = 2 + GC<sub>3</sub> + 29/[GC<sub>3</sub><sup>2</sup> + (1 - GC<sub>3</sub>)<sup>2</sup>]) (<xref ref-type="bibr" rid="B48">Romero et&#xa0;al., 2000</xref>) 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 (<xref ref-type="bibr" rid="B56">Sueoka, 1988</xref>).</p>
</sec>
<sec id="s2_5">
<label>2.5</label>
<title>Analysis of repeat fragments and prediction of RNA editing sites</title>
<p>Simple repetitive sequences (SSRs) were identified by the online tool MISA (<ext-link ext-link-type="uri" xlink:href="https://webblast.ipk-gatersleben.de/misa/">https://webblast.ipk-gatersleben.de/misa/</ext-link>) (<xref ref-type="bibr" rid="B3">Beier et&#xa0;al., 2017</xref>). 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 (<ext-link ext-link-type="uri" xlink:href="https://tandem.bu.edu/trf/home">https://tandem.bu.edu/trf/home</ext-link>) (<xref ref-type="bibr" rid="B4">Benson, 1999</xref>) with default parameters. Dispersed repeats (including forward repeats, reverse repeats, palindromic repeats and complementary repeats) were analyzed by online tool REPuter (<ext-link ext-link-type="uri" xlink:href="https://bibiserv.cebitec.uni-bielefeld.de/reputer">https://bibiserv.cebitec.uni-bielefeld.de/reputer</ext-link>) (<xref ref-type="bibr" rid="B35">Kurtz et&#xa0;al., 2001</xref>) 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 (<xref ref-type="bibr" rid="B18">Edera et&#xa0;al., 2021</xref>) was used to predict C-to-U RNA editing sites in 40 PCGs of <italic>L. aggregata</italic>, and the prediction value was set to 0.9.</p>
</sec>
<sec id="s2_6">
<label>2.6</label>
<title>Identification and collinear analysis of homologous fragments</title>
<p>Based on our previously assembled and annotated chloroplast genome of <italic>L. aggregata</italic> (PP199190) (<xref ref-type="bibr" rid="B52">Shi et&#xa0;al., 2024</xref>), 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<sup>&#x2212;5</sup>. To further study the evolutionary relationship among species, we identified the homologous sequences among <italic>L. aggregata</italic> and other 7 species of Laurales by BLAST program with same parameters as above. The collinear maps within <italic>L. aggregata</italic> (mitochondria and chloroplast) and among different species were visualized by TBtools software (<xref ref-type="bibr" rid="B11">Chen et&#xa0;al., 2020</xref>). Furthermore, we used Mauve v2.4.0 software (<xref ref-type="bibr" rid="B14">Darling et&#xa0;al., 2004</xref>) to detect the rearrangement of 7 Laurales mitogenome, and the <italic>L. aggregate</italic> mitogenome was used as the reference genome.</p>
</sec>
<sec id="s2_7">
<label>2.7</label>
<title>Phylogenetic analysis</title>
<p>We downloaded 29 complete mitogenome sequences of 9 different orders (Laurales, Magnoliales, Piperales, Chloranthales, Amborellales, Nymphaeales, Acorales, Poales, and Apiales) from NCBI database. <italic>Amborella trichopoda</italic> and <italic>Nymphaea colorata</italic> were used as outgroups to construct a phylogenetic tree with mitogenome of <italic>L. aggregata</italic>. 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 (<xref ref-type="bibr" rid="B31">Kalyaanamoorthy et&#xa0;al., 2017</xref>). The MrBayes v3.2.7a software (<xref ref-type="bibr" rid="B49">Ronquist et&#xa0;al., 2012</xref>) was used to construct BI phylogenetic tree, and the best model (GTR+I) was selected by jModelTest tool (<xref ref-type="bibr" rid="B15">Darriba et&#xa0;al., 2012</xref>). 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 &lt; 0.01. ML and BI phylogenetic trees were visualized using Chiplot online tool (<xref ref-type="bibr" rid="B69">Xie et&#xa0;al., 2023</xref>). The accession numbers of 30 species were recorded in <xref ref-type="supplementary-material" rid="ST1">
<bold>Supplementary Table S4</bold>
</xref>.</p>
</sec>
</sec>
<sec id="s3" sec-type="results">
<label>3</label>
<title>Results</title>
<sec id="s3_1">
<label>3.1</label>
<title>Assembly and annotation of <italic>L. aggregata</italic> mitochondrial genome</title>
<p>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 <italic>L. aggregata</italic>. The structure of mitogenome was very complex and had multi-branched conformation (<xref ref-type="fig" rid="f1">
<bold>Figure&#xa0;1A</bold>
</xref>), 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 (<xref ref-type="fig" rid="f1">
<bold>Figure&#xa0;1B</bold>
</xref>), 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 (<xref ref-type="fig" rid="f1">
<bold>Figure&#xa0;1C</bold>
</xref>; <xref ref-type="supplementary-material" rid="ST1">
<bold>Supplementary Table S1</bold>
</xref>). Among the 40 PCGs, 23 were core genes and 17 were non-core genes. The 23 core genes include 5 ATP synthase genes (<italic>atp1</italic>, <italic>atp4</italic>, <italic>atp6</italic>, <italic>atp8</italic>, <italic>atp9</italic>), 4 cytochrome c biosynthesis genes (<italic>ccmB</italic>, <italic>ccmC</italic>, <italic>ccmFC</italic>, <italic>ccmFN</italic>), 1 ubichinol cytochrome c reductase gene (<italic>cob</italic>), 3 cytochrome c oxidase genes (<italic>cox1</italic>, <italic>cox2</italic>, <italic>cox3</italic>), 1 transport membrance protein (<italic>mttB</italic>) and 9 NADH dehydrogenase genes (<italic>nad1</italic>, <italic>nad2</italic>, <italic>nad3</italic>, <italic>nad4</italic>, <italic>nad4L</italic>, <italic>nad5</italic>, <italic>nad6</italic>, <italic>nad7</italic>, and <italic>nad9</italic>). However, it lacked a matures gene (<italic>matR</italic>) as same as the <italic>Cinnamomum</italic> species (<xref ref-type="bibr" rid="B7">Bi et&#xa0;al., 2024b</xref>; <xref ref-type="bibr" rid="B25">Han et&#xa0;al., 2024</xref>). The non-core genes include 4 large subunits of the ribosomal proteins (<italic>rpl10</italic>, <italic>rpl16</italic>, <italic>rpl2</italic>, <italic>rpl5</italic>), 11 small subunits of the ribosomal proteins (<italic>rps1</italic>, <italic>rps10</italic>, <italic>rps11</italic>, <italic>rps12</italic>, <italic>rps13</italic>, <italic>rps14</italic>, <italic>rps19</italic>, <italic>rps2</italic>, <italic>rps3</italic>, <italic>rps4</italic>, <italic>rps7</italic>) and 2 succinate dehydrogenase genes (<italic>sdh3</italic>, <italic>sdh4</italic>). Furthermore, 8 <italic>cis</italic>-splicing genes (<italic>ccmFC</italic>, <italic>cox2</italic>, <italic>nad4</italic>, <italic>nad7</italic>, <italic>rpl2</italic>, <italic>rps10</italic>, <italic>rps3</italic> and <italic>trnL-UAG</italic>) and 3 <italic>trans</italic>-splicing genes (<italic>nad1/2/5</italic>) were identified. The details of exon in each gene can be found in <xref ref-type="supplementary-material" rid="ST1">
<bold>Supplementary Figure S1</bold>
</xref>.</p>
<fig id="f1" position="float">
<label>Figure&#xa0;1</label>
<caption>
<p>Map of mitochondrial genome assembly and annotation of <italic>L. aggregata.</italic> <bold>(A)</bold> The map of the preliminary assembly of the mitochondrial genome. <bold>(B)</bold> The master circle and linear fragment map of mitochondrial genome. Each contig is marked with different number, and the numbers order of connection is the order of unlocking circles. <bold>(C)</bold> The map of mitochondrial genome annotation. Genes with different functions are described in different colors. The colored parabola in the center circle represents the Dispersed Repeats.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fpls-15-1439245-g001.tif"/>
</fig>
</sec>
<sec id="s3_2">
<label>3.2</label>
<title>Comparative analysis of mitogenome among <italic>L. aggregata</italic> and other 7 Laurales species</title>
<p>In order to further explore the evolutionary characteristics of the mitochondrial genome, we compared it with seven other species of Laurales (<italic>Cinnamomum camphora</italic>, <italic>C. chago</italic>, <italic>C. chekiangense</italic>, <italic>C. insularimontanum</italic>, <italic>Caryodaphnopsis henryi</italic>, <italic>Machilus pauhoi</italic> and <italic>Hernandia nymphaeifolia</italic>). There were great differences in mitogenome size among 7 Laurales species, ranging from 535,805 bp (<italic>H. nymphaeifolia</italic>) to 1,168,029 bp (<italic>C. henry</italic>), and the GC content ranged from 45.73% (<italic>H. nymphaeifolia</italic>) to 47.12% (<italic>C. chago</italic>) (<xref ref-type="supplementary-material" rid="ST1">
<bold>Supplementary Figure S2</bold>
</xref>). 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 <xref ref-type="supplementary-material" rid="ST1">
<bold>Supplementary Table S2</bold>
</xref>.</p>
<p>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 (<xref ref-type="supplementary-material" rid="ST1">
<bold>Supplementary Table S2</bold>
</xref>). According to the comparative analysis, it can be found that <italic>H. nymphaeifolia</italic> had one more PCGs (<italic>matR</italic>) than other species, and had the largest number of genes (71), while <italic>C. camphora</italic> 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 (<xref ref-type="supplementary-material" rid="ST1">
<bold>Supplementary Figure S3</bold>
</xref>). Interestingly, our results showed that there was a close sister relationship between <italic>L. aggregata</italic> and four species of <italic>Cinnamomum</italic>.</p>
<p>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 (<xref ref-type="fig" rid="f2">
<bold>Figure&#xa0;2</bold>
</xref>). It was worth noting that <italic>H. nymphaeifolia</italic> showed the most typical haplotype composition, followed by <italic>C. henryi</italic>, while other species showed more mixed composition. Although <italic>L. aggregata</italic> had 20 unique haplotypes, it also showed mixed haplotypes with <italic>M. pauhoi</italic> and <italic>C. chago</italic>, including 6 haplotypes shared with <italic>M. pauhoi</italic>, 4 haplotypes shared with <italic>C. chago</italic>, and 10 haplotypes shared with 2 of them. Furthermore, <italic>C. insularimontanum</italic> and <italic>C. chekiangense</italic> shared all genes.</p>
<fig id="f2" position="float">
<label>Figure&#xa0;2</label>
<caption>
<p>Haplotype networks based on 40 shared PCGs of 8 Laurales species. Different species are given different colors.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fpls-15-1439245-g002.tif"/>
</fig>
<p>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 (<xref ref-type="fig" rid="f3">
<bold>Figure&#xa0;3</bold>
</xref>). 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 <italic>ccmB</italic>, <italic>ccmFC</italic>, <italic>rps10</italic>, <italic>rps11</italic> and <italic>rps7</italic> 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.</p>
<fig id="f3" position="float">
<label>Figure&#xa0;3</label>
<caption>
<p>The Ka/Ks ratio values of 40 shared PCGs in the mitochondrial genomes of 8 Laurales plants. Different species are given different colors.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fpls-15-1439245-g003.tif"/>
</fig>
</sec>
<sec id="s3_3">
<label>3.3</label>
<title>Analysis of codon usage bias</title>
<p>Codon usage bias analysis of 8 mitochondrial genomes in Laurales showed that 61 codons were detected, encoding 20 amino acids (<xref ref-type="fig" rid="f4">
<bold>Figure&#xa0;4</bold>
</xref>). 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%) (<xref ref-type="fig" rid="f4">
<bold>Figure&#xa0;4B</bold>
</xref>). 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 &lt; 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) (<xref ref-type="fig" rid="f4">
<bold>Figure&#xa0;4A</bold>
</xref>). 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.</p>
<fig id="f4" position="float">
<label>Figure&#xa0;4</label>
<caption>
<p>Codon usage bias analysis of 8 mitochondrial genomes in Laurales. <bold>(A)</bold> Analysis of relative synonymous codon usage (RSCU) from 8 species mitogenome. <bold>(B)</bold> The total codon usage of 40 shared PCGs in the mitochondrial genome from 8 Laurales plants.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fpls-15-1439245-g004.tif"/>
</fig>
<p>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 GC<sub>1</sub> values were ranged from 38.12% to 58.24%, the GC<sub>2</sub> values were ranged from 34.71% to 56%, and the GC<sub>3</sub> values were ranged from 27.35% to 46.88%. The average contents of GC in different locations (GC<sub>1</sub>, GC<sub>2</sub> and GC<sub>3</sub>) were all less than 50%, indicating that there was a bias for A/T base and A/T ending codon in the mitogenome of <italic>L. aggregata</italic>. 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 <italic>L. aggregata</italic> (<xref ref-type="fig" rid="f5">
<bold>Figure&#xa0;5A</bold>
</xref>). 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 (<xref ref-type="fig" rid="f5">
<bold>Figure&#xa0;5B</bold>
</xref>), 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 (<xref ref-type="fig" rid="f4">
<bold>Figure&#xa0;4</bold>
</xref>).</p>
<fig id="f5" position="float">
<label>Figure&#xa0;5</label>
<caption>
<p>Codon usage bias analysis of PCGs in mitochondrial genome of <italic>L. aggregata</italic>. <bold>(A)</bold> GC content of different sites from 40 PCGs. <bold>(B)</bold> ENC-plot of 40 PCGs. The black curve is the standard curve, the formula is ENC = 2 + GC<sub>3</sub> + 29/[GC<sub>3</sub>
<sup>2</sup> + (1 - GC<sub>3</sub>)<sup>2</sup>], and each red dot represents a gene.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fpls-15-1439245-g005.tif"/>
</fig>
</sec>
<sec id="s3_4">
<label>3.4</label>
<title>Analysis of repeat fragments and prediction of RNA editing sites</title>
<p>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 (<xref ref-type="fig" rid="f6">
<bold>Figure&#xa0;6A</bold>
</xref>). 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 <italic>H. nymphaeifolia</italic> 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 (<xref ref-type="fig" rid="f6">
<bold>Figure&#xa0;6B</bold>
</xref>). 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 (<xref ref-type="fig" rid="f6">
<bold>Figure&#xa0;6D</bold>
</xref>). Moreover, <italic>C. henryi</italic> and <italic>H. nymphaeifolia</italic> had the most and fewest repeat fragments respectively, which might be an important reason for the size of the mitochondrial genome.</p>
<fig id="f6" position="float">
<label>Figure&#xa0;6</label>
<caption>
<p>Repeat fragment analysis and prediction of RNA editing sites of mitochondrial genome. <bold>(A)</bold> Detection of mitogenome SSRs in 8 species of Laurales. Mono, Di, Tri, Tetra, Penta, and Hexa represent mononucleotide repeats, dinucleotide repeats, trinucleotide repeats, tetranucleotide repeats, pentanucleotide repeats and hexanucleotide repeats, respectively. <bold>(B)</bold> Prediction of dispersed and tandem repeats of mitogenome in 8 species of Laurales. F, R, P, C, TR represent forward repeats, reverse repeats, palindromic repeats, complementary repeats, and tandem repeats. <bold>(C)</bold> Prediction of RNA editing sites in the mitogenome of <italic>L. aggregata</italic>. <bold>(D)</bold> The size distribution of dispersed repeats in 8 Laurales mitogenome.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fpls-15-1439245-g006.tif"/>
</fig>
<p>In this study, we also predicted all possible C-to-U RNA editing sites in the <italic>L. aggregata</italic> 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 (<xref ref-type="fig" rid="f6">
<bold>Figure&#xa0;6C</bold>
</xref>). It is worth noting that the <italic>nad4</italic> gene showed the highest RNA editing frequency, predicting 68 editing sites, followed by the <italic>nad5</italic> gene with 48 editing sites. However, most ribosomal proteins had a small number of editing sites, which was consistent with <italic>C. chekiangense</italic> and <italic>C. camphora</italic> (<xref ref-type="bibr" rid="B7">Bi et&#xa0;al., 2024b</xref>; <xref ref-type="bibr" rid="B25">Han et&#xa0;al., 2024</xref>). 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.</p>
</sec>
<sec id="s3_5">
<label>3.5</label>
<title>DNA fragment transfer analysis between mitochondrial and chloroplast genomes</title>
<p>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 (<xref ref-type="fig" rid="f7">
<bold>Figure&#xa0;7A</bold>
</xref>). 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 (<xref ref-type="supplementary-material" rid="ST1">
<bold>Supplementary Table S3</bold>
</xref>). 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 (<italic>petN</italic> (2)) and 13 tRNA genes (<italic>trnA-UGC</italic>, <italic>trnD-GUC</italic>, <italic>trnE-UUC</italic>, <italic>trnH-GUG</italic>, <italic>trnI-CAU</italic>, <italic>trnM-CAU</italic>, <italic>trnN-GUU</italic> (2), <italic>trnP-UGG</italic>, <italic>trnR-UCU</italic>, <italic>trnT-GGU</italic>, <italic>trnW-CCA</italic>, <italic>trnY-GUA</italic>). It was worth noting that the <italic>petN</italic> 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.</p>
<fig id="f7" position="float">
<label>Figure&#xa0;7</label>
<caption>
<p>Collinearity analysis among different Laurales species. <bold>(A)</bold> Sequences transfer analysis of chloroplast genome (cpDNA) and mitochondrial genome in <italic>L. aggregata</italic>. The blue and green arcs represent mitogenome and cpDNA, respectively. Homologous fragments are represented by a yellow line between blue and green arcs. <bold>(B)</bold> Collinearity analysis between the mitogenome of <italic>L. aggregata</italic> and other 7 mitogenome of Laurales species. The arcs of different colors represent different mitochondrial genomes and are abbreviated to represent different species. <italic>Lindera aggregata</italic> (chr1 and chr2), <italic>Cinnamomum camphora</italic> (Cca), <italic>Cinnamomum chago</italic> (Ccha), <italic>Cinnamomum chekiangense</italic> (Cche), <italic>Cinnamomum insularimontanum</italic> (Cin), <italic>Caryodaphnopsis henryi</italic> (Che), <italic>Machilus pauhoi</italic> (Mpa_1 and Mpa_2) and <italic>Hernandia nymphaeifolia</italic> (Hny). The homologous fragments between 7 species and <italic>L. aggregata</italic> are represented by their corresponding colors, respectively.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fpls-15-1439245-g007.tif"/>
</fig>
</sec>
<sec id="s3_6">
<label>3.6</label>
<title>Collinearity analysis and structural rearrangement</title>
<p>In order to further study the homologous sequences of the mitogenome in 8 Laurales species, we identified the collinearity between <italic>L. aggregata</italic> and other Laurales species (<xref ref-type="fig" rid="f7">
<bold>Figure&#xa0;7B</bold>
</xref> and <xref ref-type="supplementary-material" rid="ST1">
<bold>Supplementary Figure S4</bold>
</xref>). The results showed that <italic>L. aggregata</italic> had a large number of collinear fragments with other 7 species in Laurales. Notably, the highest collinearity was detected between <italic>L. aggregate</italic> and <italic>C. camphora</italic>, which contained 1051 collinearity fragments totaling 598,706 bp, accounting for 65.61% of <italic>L. aggregate</italic> mitogenome. However, <italic>H. nymphaeifolia</italic> exhibited the weakest collinearity, with only 378 collinearity fragments spanning 234,099 bp, covering 25.65% of <italic>L. aggregate</italic> 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 <italic>L. aggregate</italic> and <italic>C. camphora</italic>, containing 99.15% of the protein-coding sequences from <italic>C. camphora</italic> mitogenome, while it was relatively lower with <italic>H. nymphaeifolia</italic> (95.93%). This showed that although the PCGs of these 8 Laurales mitogenome had high similarity in sequences, the overall mitogenome structure was unstable.</p>
</sec>
<sec id="s3_7">
<label>3.7</label>
<title>Phylogenetic analysis of different orders mitochondrial genomes</title>
<p>We constructed phylogenetic trees based on 24 conserved PCGs (<italic>atp1/4/6/8/9</italic>, <italic>ccmB/C/FN</italic>, <italic>cox1/2/3</italic>, <italic>cob</italic>, <italic>mttB</italic>, <italic>nad1/2/3/4/4L/5/6/7/9</italic>, <italic>rps12</italic> and <italic>sdh4</italic>) 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 (<xref ref-type="fig" rid="f8">
<bold>Figures&#xa0;8A, B</bold>
</xref>). 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 <italic>L. aggregata</italic> belonged to Laurales and represented a relatively new differentiation. In addition, <italic>L. aggregata</italic> was closely related to <italic>M. pauhoi</italic> and <italic>C. chago</italic>, which was consistent with the results of the haplotype network mentioned above.</p>
<fig id="f8" position="float">
<label>Figure&#xa0;8</label>
<caption>
<p>Phylogenetic trees of 30 plant mitochondrial genomes from 9 different Orders. <bold>(A)</bold> The maximum likelihood phylogenetic tree based on 24 common PCGs of 30 species. <bold>(B)</bold> The Bayesian inference phylogenetic tree based on 24 common PCGs of 30 species. The support rate values are shown on the branches on the trees.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fpls-15-1439245-g008.tif"/>
</fig>
</sec>
</sec>
<sec id="s4" sec-type="discussion">
<label>4</label>
<title>Discussion</title>
<p>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 (<xref ref-type="bibr" rid="B19">Fujii et&#xa0;al., 2010</xref>; <xref ref-type="bibr" rid="B55">Smith and Keeling, 2015</xref>). Recent studies have revealed great differences in the structure of mitogenome among different species. For example, the mitogenome of <italic>Angelica biserrata</italic> consists of six independent circular chromosomes, while the mitogenome of <italic>Mentha spicata</italic> consists of one linear chromosome and two circular chromosomes, each of which varies in size (<xref ref-type="bibr" rid="B30">Jiang et&#xa0;al., 2023</xref>; <xref ref-type="bibr" rid="B61">Wang et&#xa0;al., 2024</xref>). 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% (<xref ref-type="bibr" rid="B33">Kubo and Mikami, 2007</xref>). The GC content of the PCGs in the mitogenome of <italic>L. aggregata</italic> 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 (<xref ref-type="bibr" rid="B70">Xie et&#xa0;al., 2019</xref>). In our comparison, we found that most of the 40 common PCGs in 8 Laurales mitogenome were under negative selection, only a few genes (<italic>ccmB</italic>, <italic>ccmFC</italic>, <italic>rps10</italic>, <italic>rps11</italic> and <italic>rps7</italic>) showed signs of positive selection in all species, while some genes (<italic>atp6</italic>, <italic>ccmFN</italic>, <italic>cob</italic>, <italic>nad3</italic>, <italic>nad4</italic>, <italic>nad5</italic>, <italic>rps2</italic> and <italic>rps3</italic>) were affected by positive selection in individual species. The <italic>atp6</italic>, <italic>ccmB</italic>, <italic>ccmFC</italic>, <italic>ccmFN</italic>, <italic>cob</italic>, <italic>nad3</italic>, <italic>nad4</italic> and <italic>nad5</italic> genes may play an important role in redox reaction and other adaptation processes. The relatively high Ka/Ks ratio values of <italic>rps10</italic> and <italic>ccmB</italic> genes indicated that they played an important role in species evolution and may be closely related to environmental stress.</p>
<p>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 (<xref ref-type="bibr" rid="B1">Andargie and Congyi, 2022</xref>; <xref ref-type="bibr" rid="B68">Xiao et&#xa0;al., 2024</xref>). 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 (<xref ref-type="bibr" rid="B28">Hershberg and Petrov, 2008</xref>; <xref ref-type="bibr" rid="B46">Parvathy et&#xa0;al., 2022</xref>). 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 <italic>Angelica</italic> and <italic>Mangifera</italic> species (<xref ref-type="bibr" rid="B44">Niu et&#xa0;al., 2022</xref>; <xref ref-type="bibr" rid="B61">Wang et&#xa0;al., 2024</xref>), 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 (<xref ref-type="bibr" rid="B74">Zhang et&#xa0;al., 2018</xref>; <xref ref-type="bibr" rid="B63">Wang et&#xa0;al., 2020</xref>; <xref ref-type="bibr" rid="B58">Tang et&#xa0;al., 2022</xref>). ENC-plot analysis and RSCU results showed that there was a preference for codon usage in <italic>L. aggregata</italic> 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 (<xref ref-type="bibr" rid="B59">Tang et&#xa0;al., 2020</xref>; <xref ref-type="bibr" rid="B61">Wang et&#xa0;al., 2024</xref>).</p>
<p>The collinear analysis of homologous fragments among 8 species of Laurales showed that the mitogenome of <italic>L. aggregata</italic> 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 (<xref ref-type="bibr" rid="B47">Pfeifer et&#xa0;al., 2013</xref>; <xref ref-type="bibr" rid="B38">Li et&#xa0;al., 2022</xref>). A total of 303 SSRs,71 tandem repeats and 1598 dispersed repeats were detected in the mitogenome of <italic>L. aggregata</italic>. These repetitive fragments constitute the intricate branch structure of the mitochondrial genome of <italic>L. aggregata</italic>. 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 (<xref ref-type="bibr" rid="B24">Gualberto et&#xa0;al., 2014</xref>). For example, there is a similar trend of pseudogenes loss or alteration in the chloroplast genome of the <italic>Cistanche</italic> genus, affecting five genes related to photosynthesis and energy production (<xref ref-type="bibr" rid="B41">Miao et&#xa0;al., 2022</xref>). 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 (<italic>petN</italic> (2)) and 13 tRNA genes (<italic>trnA-UGC</italic>, <italic>trnD-GUC</italic>, <italic>trnE-UUC</italic>, <italic>trnH-GUG</italic>, <italic>trnI-CAU</italic>, <italic>trnM-CAU</italic>, <italic>trnN-GUU</italic> (2), <italic>trnP-UGG</italic>, <italic>trnR-UCU</italic>, <italic>trnT-GGU</italic>, <italic>trnW-CCA</italic>, <italic>trnY-GUA</italic>). It was worth noting that the <italic>petN</italic> 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.</p>
<p>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 (<xref ref-type="bibr" rid="B57">Sun et&#xa0;al., 2016</xref>). 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) (<xref ref-type="bibr" rid="B76">Zhu et&#xa0;al., 2022</xref>). The identification of RNA editing is very important for the study of mitogenome in <italic>L. aggregata</italic>. 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 (<xref ref-type="bibr" rid="B21">Gray, 2009</xref>; <xref ref-type="bibr" rid="B5">Bi et&#xa0;al., 2019</xref>). A total of 770 RNA editing sites were identified in the mitogenome of <italic>L. aggregata</italic>, of which the number of editing event of <italic>nad4</italic> 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 <italic>matR</italic> gene was missing from the mitogenome of Lauraceae species, while <italic>matR</italic> gene was found in its closest relatives, Hernandiaceae species. This shows that <italic>matR</italic> 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.</p>
</sec>
<sec id="s5" sec-type="conclusions">
<label>5</label>
<title>Conclusion</title>
<p>In this study, we completed the assembly and annotation of the mitogenome of <italic>L. aggregata</italic>. The mitogenome showed a branch structure of 912,473 bp in length. Comparative analysis showed that the mitogenome of <italic>L. aggregata</italic> 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 <italic>L. aggregata</italic>. 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 <italic>Lindera</italic> species in Magnoliids. This study will lay a foundation for the application of population genetics and evolution of <italic>Lindera</italic> and other species from Magnoliids.</p>
</sec>
</body>
<back>
<sec id="s6" sec-type="data-availability">
<title>Data availability statement</title>
<p>The accession numbers of mitochondrial genome of <italic>L. aggregata</italic> in Gene Bank are PP848112 and PP848113. The raw HiFi sequencing data have been deposited in NCBI with accession number: SRR30136615.</p>
</sec>
<sec id="s7" sec-type="author-contributions">
<title>Author contributions</title>
<p>YS: Writing &#x2013; review &amp; editing, Writing &#x2013; original draft, Visualization, Validation, Supervision, Software, Resources, Project administration, Methodology, Investigation, Formal analysis, Data curation, Conceptualization. ZC: Data curation, Funding acquisition, Methodology, Writing &#x2013; review &amp; editing, Supervision, Software, Resources. JJ: Writing &#x2013; review &amp; editing, Visualization, Supervision, Software, Resources, Investigation, Formal analysis. WW: Writing &#x2013; review &amp; editing, Supervision, Software, Methodology, Formal analysis. WY: Writing &#x2013; review &amp; editing, Resources, Investigation, Funding acquisition. SZ: Formal analysis, Investigation, Writing &#x2013; review &amp; editing. WZ: Writing &#x2013; review &amp; editing, Visualization, Validation, Supervision, Software, Resources, Project administration, Methodology, Investigation, Funding acquisition, Formal analysis, Data curation, Conceptualization.</p>
</sec>
<sec id="s8" sec-type="funding-information">
<title>Funding</title>
<p>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).</p>
</sec>
<sec id="s9" sec-type="COI-statement">
<title>Conflict of interest</title>
<p>Author WY was employed by the company Zhejiang Hongshiliang Group Tiantai Mountain Wu-Yao Co., Ltd.</p>
<p>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.</p>
</sec>
<sec id="s10" sec-type="disclaimer">
<title>Publisher&#x2019;s note</title>
<p>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.</p>
</sec>
<sec id="s11" sec-type="supplementary-material">
<title>Supplementary material</title>
<p>The Supplementary Material for this article can be found online at: <ext-link ext-link-type="uri" xlink:href="https://www.frontiersin.org/articles/10.3389/fpls.2024.1439245/full#supplementary-material">https://www.frontiersin.org/articles/10.3389/fpls.2024.1439245/full#supplementary-material</ext-link>
</p>
<supplementary-material xlink:href="Table1.docx" id="ST1" mimetype="application/vnd.openxmlformats-officedocument.wordprocessingml.document"/>
<supplementary-material xlink:href="Table2.xlsx" id="ST2" mimetype="application/vnd.openxmlformats-officedocument.spreadsheetml.sheet"/>
</sec>
<ref-list>
<title>References</title>
<ref id="B1">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Andargie</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Congyi</surname> <given-names>Z.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>Genome-wide analysis of codon usage in sesame (Sesamum indicum L.)</article-title>. <source>Heliyon</source> <volume>8</volume>, <elocation-id>e08687</elocation-id>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.heliyon.2021.e08687</pub-id>
</citation>
</ref>
<ref id="B2">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Arrieta-Montiel</surname> <given-names>M. P.</given-names>
</name>
<name>
<surname>Mackenzie</surname> <given-names>S. A.</given-names>
</name>
</person-group> (<year>2011</year>). <article-title>Plant mitochondrial genomes and recombination</article-title>. <source>Plant mitochondria</source> <volume>1</volume>, <fpage>65</fpage>&#x2013;<lpage>82</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1007/978-0-387-89781-3</pub-id>
</citation>
</ref>
<ref id="B3">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Beier</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Thiel</surname> <given-names>T.</given-names>
</name>
<name>
<surname>M&#xfc;nch</surname> <given-names>T.</given-names>
</name>
<name>
<surname>Scholz</surname> <given-names>U.</given-names>
</name>
<name>
<surname>Mascher</surname> <given-names>M.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>MISA-web: a web server for microsatellite prediction</article-title>. <source>Bioinformatics</source> <volume>33</volume>, <fpage>2583</fpage>&#x2013;<lpage>2585</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1093/bioinformatics/btx198</pub-id>
</citation>
</ref>
<ref id="B4">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Benson</surname> <given-names>G.</given-names>
</name>
</person-group> (<year>1999</year>). <article-title>Tandem repeats finder: a program to analyze DNA sequences</article-title>. <source>Nucleic Acids Res.</source> <volume>27</volume>, <fpage>573</fpage>&#x2013;<lpage>580</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1093/nar/27.2.573</pub-id>
</citation>
</ref>
<ref id="B5">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Bi</surname> <given-names>C.</given-names>
</name>
<name>
<surname>Paterson</surname> <given-names>A. H.</given-names>
</name>
<name>
<surname>Wang</surname> <given-names>X.</given-names>
</name>
<name>
<surname>Xu</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Wu</surname> <given-names>D.</given-names>
</name>
<name>
<surname>Qu</surname> <given-names>Y.</given-names>
</name>
<etal/>
</person-group>. (<year>2019</year>). <article-title>Analysis of the complete mitochondrial genome sequence of the diploid cotton Gossypium raimondii by comparative genomics approaches</article-title>. <source>BioMed. Res. Int.</source> <volume>2019</volume>, <fpage>9691253</fpage>. doi: <pub-id pub-id-type="doi">10.1155/2019/9691253</pub-id>
</citation>
</ref>
<ref id="B6">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Bi</surname> <given-names>C.</given-names>
</name>
<name>
<surname>Shen</surname> <given-names>F.</given-names>
</name>
<name>
<surname>Han</surname> <given-names>F.</given-names>
</name>
<name>
<surname>Qu</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Hou</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Xu</surname> <given-names>K.</given-names>
</name>
<etal/>
</person-group>. (<year>2024</year>a). <article-title>PMAT: an efficient plant mitogenome assembly toolkit using low-coverage HiFi sequencing data</article-title>. <source>Horticul Res.</source> <volume>11</volume>, <elocation-id>uhae023</elocation-id>. doi:&#xa0;<pub-id pub-id-type="doi">10.1093/hr/uhae023</pub-id>
</citation>
</ref>
<ref id="B7">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Bi</surname> <given-names>C.</given-names>
</name>
<name>
<surname>Sun</surname> <given-names>N.</given-names>
</name>
<name>
<surname>Han</surname> <given-names>F.</given-names>
</name>
<name>
<surname>Xu</surname> <given-names>K.</given-names>
</name>
<name>
<surname>Yang</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Ferguson</surname> <given-names>D. K.</given-names>
</name>
</person-group> (<year>2024</year>b). <article-title>The first mitogenome of Lauraceae (Cinnamomum chekiangense)</article-title>. <source>Plant Diversity</source> <volume>46</volume>, <fpage>144</fpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.pld.2023.11.001</pub-id>
</citation>
</ref>
<ref id="B8">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Cao</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Xuan</surname> <given-names>B.</given-names>
</name>
<name>
<surname>Peng</surname> <given-names>B.</given-names>
</name>
<name>
<surname>Li</surname> <given-names>C.</given-names>
</name>
<name>
<surname>Chai</surname> <given-names>X.</given-names>
</name>
<name>
<surname>Tu</surname> <given-names>P.</given-names>
</name>
</person-group> (<year>2015</year>). <article-title>The genus Lindera: a source of structurally diverse molecules having pharmacological significance</article-title>. <source>Phytochem. Rev.</source> <volume>15</volume>, <fpage>869</fpage>&#x2013;<lpage>906</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1007/s11101-015-9432-2</pub-id>
</citation>
</ref>
<ref id="B9">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Capella-Guti&#xe9;rrez</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Silla-Mart&#xed;nez</surname> <given-names>J. M.</given-names>
</name>
<name>
<surname>Gabald&#xf3;n</surname> <given-names>T.</given-names>
</name>
</person-group> (<year>2009</year>). <article-title>trimAl: a tool for automated alignment trimming in large-scale phylogenetic analyses</article-title>. <source>Bioinformatics</source> <volume>25</volume>, <fpage>1972</fpage>&#x2013;<lpage>1973</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1093/bioinformatics/btp348</pub-id>
</citation>
</ref>
<ref id="B10">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Chan</surname> <given-names>P. p.</given-names>
</name>
<name>
<surname>Lin</surname> <given-names>B. y.</given-names>
</name>
<name>
<surname>Mak</surname> <given-names>A. j.</given-names>
</name>
<name>
<surname>Lowe</surname> <given-names>T. m.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>tRNAscan-SE 2.0: improved detection and functional classification of transfer RNA genes</article-title>. <source>Nucleic Acids Res.</source> <volume>49</volume>, <fpage>9077</fpage>&#x2013;<lpage>9096</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1093/nar/gkab688</pub-id>
</citation>
</ref>
<ref id="B11">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Chen</surname> <given-names>C.</given-names>
</name>
<name>
<surname>Chen</surname> <given-names>H.</given-names>
</name>
<name>
<surname>Zhang</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Thomas</surname> <given-names>H. R.</given-names>
</name>
<name>
<surname>Frank</surname> <given-names>M. H.</given-names>
</name>
<name>
<surname>He</surname> <given-names>Y.</given-names>
</name>
<etal/>
</person-group>. (<year>2020</year>). <article-title>TBtools: An Integrative toolkit developed for interactive analyses of big biological data</article-title>. <source>Mol. Plant</source> <volume>13</volume>, <fpage>1194</fpage>&#x2013;<lpage>1202</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.molp.2020.06.009</pub-id>
</citation>
</ref>
<ref id="B12">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Chen</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Ye</surname> <given-names>W.</given-names>
</name>
<name>
<surname>Zhang</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Xu</surname> <given-names>Y.</given-names>
</name>
</person-group> (<year>2015</year>). <article-title>High speed BLASTN: an accelerated MegaBLAST search tool</article-title>. <source>Nucleic Acids Res.</source> <volume>43</volume>, <fpage>7762</fpage>&#x2013;<lpage>7768</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1093/nar/gkv784</pub-id>
</citation>
</ref>
<ref id="B13">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Covello</surname> <given-names>P. S.</given-names>
</name>
<name>
<surname>Gray</surname> <given-names>M. W.</given-names>
</name>
</person-group> (<year>1989</year>). <article-title>RNA editing in plant mitochondria</article-title>. <source>Nature</source> <volume>341</volume>, <fpage>662</fpage>&#x2013;<lpage>666</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1038/341662a0</pub-id>
</citation>
</ref>
<ref id="B14">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Darling</surname> <given-names>A. C.</given-names>
</name>
<name>
<surname>Mau</surname> <given-names>B.</given-names>
</name>
<name>
<surname>Blattner</surname> <given-names>F. R.</given-names>
</name>
<name>
<surname>Perna</surname> <given-names>N. T.</given-names>
</name>
</person-group> (<year>2004</year>). <article-title>Mauve: multiple alignment of conserved genomic sequence with rearrangements</article-title>. <source>Genome Res.</source> <volume>14</volume>, <fpage>1394</fpage>&#x2013;<lpage>1403</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1101/gr.2289704</pub-id>
</citation>
</ref>
<ref id="B15">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Darriba</surname> <given-names>D.</given-names>
</name>
<name>
<surname>Taboada</surname> <given-names>G. L.</given-names>
</name>
<name>
<surname>Doallo</surname> <given-names>R.</given-names>
</name>
<name>
<surname>Posada</surname> <given-names>D.</given-names>
</name>
</person-group> (<year>2012</year>). <article-title>jModelTest 2: more models, new heuristics and high-performance computing</article-title>. <source>Nat. Methods</source> <volume>9</volume>, <fpage>772</fpage>. doi: <pub-id pub-id-type="doi">10.1038/nmeth.2109</pub-id>
</citation>
</ref>
<ref id="B16">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Doyle</surname> <given-names>J. J.</given-names>
</name>
<name>
<surname>Doyle</surname> <given-names>J. L.</given-names>
</name>
</person-group> (<year>1987</year>). <article-title>A rapid DNA isolation procedure for small quantities of fresh leaf tissue</article-title>. <source>Phytochem. Bull.</source>, <volume>19</volume>, <fpage>11</fpage>&#x2013;<lpage>15</lpage>.</citation>
</ref>
<ref id="B17">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Duan</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Wang</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Kang</surname> <given-names>F.</given-names>
</name>
<name>
<surname>Li</surname> <given-names>J.</given-names>
</name>
</person-group> (<year>2018</year>). <article-title>Research of molecular evolution characteristics of nrDNA ITS and cpDNA trnH-psbA sequences in Prunus mongolica</article-title>. <source>Genomics Appl. Biol.</source> <volume>37</volume>, <fpage>3035</fpage>&#x2013;<lpage>3041</lpage>. doi: <pub-id pub-id-type="doi">10.13417/j.gab.037.003035</pub-id>
</citation>
</ref>
<ref id="B18">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Edera</surname> <given-names>A. A.</given-names>
</name>
<name>
<surname>Small</surname> <given-names>I.</given-names>
</name>
<name>
<surname>Milone</surname> <given-names>D. H.</given-names>
</name>
<name>
<surname>Sanchez-Puerta</surname> <given-names>M. V.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>Deepred-Mt: Deep representation learning for predicting C-to-U RNA editing in plant mitochondria</article-title>. <source>Comput. Biol. Med.</source> <volume>136</volume>, <fpage>104682</fpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.compbiomed.2021.104682</pub-id>
</citation>
</ref>
<ref id="B19">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Fujii</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Kazama</surname> <given-names>T.</given-names>
</name>
<name>
<surname>Yamada</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Toriyama</surname> <given-names>K.</given-names>
</name>
</person-group> (<year>2010</year>). <article-title>Discovery of global genomic re-organization based on comparison of two newly sequenced rice mitochondrial genomes with cytoplasmic male sterility-related genes</article-title>. <source>BMC Genomics</source> <volume>11</volume>, <fpage>1</fpage>&#x2013;<lpage>15</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1186/1471-2164-11-209</pub-id>
</citation>
</ref>
<ref id="B20">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Gandini</surname> <given-names>C. L.</given-names>
</name>
<name>
<surname>Garcia</surname> <given-names>L. E.</given-names>
</name>
<name>
<surname>Abbona</surname> <given-names>C. C.</given-names>
</name>
<name>
<surname>Sanchez-Puerta</surname> <given-names>M. V.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>The complete organelle genomes of Physochlaina orientalis: Insights into short sequence repeats across seed plant mitochondrial genomes</article-title>. <source>Mol. Phylogenet. Evol.</source> <volume>137</volume>, <fpage>274</fpage>&#x2013;<lpage>284</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.ympev.2019.05.012</pub-id>
</citation>
</ref>
<ref id="B21">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Gray</surname> <given-names>M. W.</given-names>
</name>
</person-group> (<year>2009</year>). <article-title>RNA editing in plant mitochondria: 20 years later</article-title>. <source>IUBMB Life</source> <volume>61</volume>, <fpage>1101</fpage>&#x2013;<lpage>1104</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1002/iub.272</pub-id>
</citation>
</ref>
<ref id="B22">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Gray</surname> <given-names>M. W.</given-names>
</name>
<name>
<surname>Covello</surname> <given-names>P. S.</given-names>
</name>
</person-group> (<year>1993</year>). <article-title>RNA editing in plant mitochondria and chloroplasts</article-title>. <source>FASEB J.</source> <volume>7</volume>, <fpage>64</fpage>&#x2013;<lpage>71</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1096/fasebj.7.1.8422976</pub-id>
</citation>
</ref>
<ref id="B23">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Grewe</surname> <given-names>F.</given-names>
</name>
<name>
<surname>Herres</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Vieh&#xf6;ver</surname> <given-names>P.</given-names>
</name>
<name>
<surname>Polsakiewicz</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Weisshaar</surname> <given-names>B.</given-names>
</name>
<name>
<surname>Knoop</surname> <given-names>V.</given-names>
</name>
</person-group> (<year>2011</year>). <article-title>A unique transcriptome: 1782 positions of RNA editing alter 1406 codon identities in mitochondrial mRNAs of the lycophyte Isoetes engelmannii</article-title>. <source>Nucleic Acids Res.</source> <volume>39</volume>, <fpage>2890</fpage>&#x2013;<lpage>2902</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1093/nar/gkq1227</pub-id>
</citation>
</ref>
<ref id="B24">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Gualberto</surname> <given-names>J. M.</given-names>
</name>
<name>
<surname>Mileshina</surname> <given-names>D.</given-names>
</name>
<name>
<surname>Wallet</surname> <given-names>C.</given-names>
</name>
<name>
<surname>Niazi</surname> <given-names>A. K.</given-names>
</name>
<name>
<surname>Weber-Lotfi</surname> <given-names>F.</given-names>
</name>
<name>
<surname>Dietrich</surname> <given-names>A.</given-names>
</name>
</person-group> (<year>2014</year>). <article-title>The plant mitochondrial genome: Dynamics and maintenance</article-title>. <source>Biochimie</source> <volume>100</volume>, <fpage>107</fpage>&#x2013;<lpage>120</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.biochi.2013.09.016</pub-id>
</citation>
</ref>
<ref id="B25">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Han</surname> <given-names>F.</given-names>
</name>
<name>
<surname>Bi</surname> <given-names>C.</given-names>
</name>
<name>
<surname>Zhao</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Gao</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Wang</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Chen</surname> <given-names>Y.</given-names>
</name>
</person-group> (<year>2024</year>). <article-title>Unraveling the complex evolutionary features of the Cinnamomum camphora mitochondrial genome</article-title>. <source>Plant Cell Rep.</source> <volume>43</volume>, <fpage>183</fpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1007/s00299-024-03256-1</pub-id>
</citation>
</ref>
<ref id="B26">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Han</surname> <given-names>H.</given-names>
</name>
<name>
<surname>Xu</surname> <given-names>B.</given-names>
</name>
<name>
<surname>Amin</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Li</surname> <given-names>H.</given-names>
</name>
<name>
<surname>Yu</surname> <given-names>X.</given-names>
</name>
<name>
<surname>Gong</surname> <given-names>M.</given-names>
</name>
<etal/>
</person-group>. (<year>2019</year>). <article-title>Quercetin&#x2212;3&#x2212;O&#x2212;&#x3b1;&#x2212;L&#x2212;rhamnopyranoside derived from the leaves of Lindera aggregata (Sims) Kosterm. evokes the autophagy&#x2212;induced nuclear factor erythroid 2&#x2212;related factor 2 antioxidant pathway in human umbilical vein endothelial cells</article-title>. <source>Int. J. Mol. Med.</source> <volume>43</volume>, <fpage>461</fpage>&#x2013;<lpage>474</lpage>. doi: <pub-id pub-id-type="doi">10.3892/ijmm.2018.3976</pub-id>
</citation>
</ref>
<ref id="B27">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Hebert</surname> <given-names>P. D. N.</given-names>
</name>
<name>
<surname>Cywinska</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Ball</surname> <given-names>S. L.</given-names>
</name>
<name>
<surname>Dewaard</surname> <given-names>J. R.</given-names>
</name>
</person-group> (<year>2003</year>). <article-title>Biological identifications through DNA barcodes</article-title>. <source>Proc. R. Soc. London. Ser. B: Biol. Sci.</source> <volume>270</volume>, <fpage>313</fpage>&#x2013;<lpage>321</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1098/rspb.2002.2218</pub-id>
</citation>
</ref>
<ref id="B28">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Hershberg</surname> <given-names>R.</given-names>
</name>
<name>
<surname>Petrov</surname> <given-names>D. A.</given-names>
</name>
</person-group> (<year>2008</year>). <article-title>Selection on codon bias</article-title>. <source>Annu. Rev. Genet.</source> <volume>42</volume>, <fpage>287</fpage>&#x2013;<lpage>299</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1146/annurev.genet.42.110807.091442</pub-id>
</citation>
</ref>
<ref id="B29">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Hung-Pin</surname> <given-names>T.</given-names>
</name>
</person-group> (<year>1987</year>). <article-title>A study on the system of lindera</article-title>. <source>J. Syst. Evol.</source> <volume>25</volume>, <fpage>161</fpage>&#x2013;<lpage>171</lpage>.</citation>
</ref>
<ref id="B30">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Jiang</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Ni</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Zhang</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Li</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Liu</surname> <given-names>C.</given-names>
</name>
</person-group> (<year>2023</year>). <article-title>Complete mitochondrial genome of Mentha spicata L. reveals multiple chromosomal configurations and RNA editing events</article-title>. <source>Int. J. Biol. Macromol</source> <volume>251</volume>, <fpage>126257</fpage>. doi: <pub-id pub-id-type="doi">10.1016/j.ijbiomac.2023.126257</pub-id>
</citation>
</ref>
<ref id="B31">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kalyaanamoorthy</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Minh</surname> <given-names>B. Q.</given-names>
</name>
<name>
<surname>Wong</surname> <given-names>T. K.</given-names>
</name>
<name>
<surname>Von Haeseler</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Jermiin</surname> <given-names>L. S.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>ModelFinder: fast model selection for accurate phylogenetic estimates</article-title>. <source>Nat. Methods 14</source>, <fpage>587</fpage>&#x2013;<lpage>589</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1038/nmeth.4285</pub-id>
</citation>
</ref>
<ref id="B32">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Katoh</surname> <given-names>K.</given-names>
</name>
<name>
<surname>Standley</surname> <given-names>D. M.</given-names>
</name>
</person-group> (<year>2013</year>). <article-title>MAFFT multiple sequence alignment software version 7: improvements in performance and usability</article-title>. <source>Mol. Biol. Evol.</source> <volume>30</volume>, <fpage>772</fpage>&#x2013;<lpage>780</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1093/molbev/mst010</pub-id>
</citation>
</ref>
<ref id="B33">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kubo</surname> <given-names>T.</given-names>
</name>
<name>
<surname>Mikami</surname> <given-names>T.</given-names>
</name>
</person-group> (<year>2007</year>). <article-title>Organization and variation of angiosperm mitochondrial genome</article-title>. <source>Physiol Plantarum</source> <volume>129</volume>, <fpage>6</fpage>&#x2013;<lpage>13</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1111/j.1399-3054.2006.00768.x</pub-id>
</citation>
</ref>
<ref id="B34">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kubo</surname> <given-names>T.</given-names>
</name>
<name>
<surname>Newton</surname> <given-names>K. J.</given-names>
</name>
</person-group> (<year>2008</year>). <article-title>Angiosperm mitochondrial genomes and mutations</article-title>. <source>Mitochondrion</source> <volume>8</volume>, <fpage>5</fpage>&#x2013;<lpage>14</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.mito.2007.10.006</pub-id>
</citation>
</ref>
<ref id="B35">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kurtz</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Choudhuri</surname> <given-names>J. V.</given-names>
</name>
<name>
<surname>Ohlebusch</surname> <given-names>E.</given-names>
</name>
<name>
<surname>Schleiermacher</surname> <given-names>C.</given-names>
</name>
<name>
<surname>Stoye</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Giegerich</surname> <given-names>R.</given-names>
</name>
</person-group> (<year>2001</year>). <article-title>REPuter: the manifold applications of repeat analysis on a genomic scale</article-title>. <source>Nucleic Acids Res.</source> <volume>29</volume>, <fpage>4633</fpage>&#x2013;<lpage>4642</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1093/nar/29.22.4633</pub-id>
</citation>
</ref>
<ref id="B36">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Leigh</surname> <given-names>J. W.</given-names>
</name>
<name>
<surname>Bryant</surname> <given-names>D.</given-names>
</name>
</person-group> (<year>2015</year>). <article-title>popart: full-feature software for haplotype network construction</article-title>. <source>Methods Ecol. Evol.</source> <volume>6</volume>, <fpage>1110</fpage>&#x2013;<lpage>1116</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1111/2041-210X.12410</pub-id>
</citation>
</ref>
<ref id="B37">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Li</surname> <given-names>H.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>New strategies to improve minimap2 alignment accuracy</article-title>. <source>Bioinformatics</source> <volume>37</volume>, <fpage>4572</fpage>&#x2013;<lpage>4574</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1093/bioinformatics/btab705</pub-id>
</citation>
</ref>
<ref id="B38">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Li</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Li</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Ma</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Kou</surname> <given-names>L.</given-names>
</name>
<name>
<surname>Wei</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Wang</surname> <given-names>W.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>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</article-title>. <source>BMC Genomics</source> <volume>23</volume>, <fpage>481</fpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1186/s12864-022-08706-2</pub-id>
</citation>
</ref>
<ref id="B39">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Lu</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Zhao</surname> <given-names>Z.</given-names>
</name>
<name>
<surname>Ni</surname> <given-names>L.</given-names>
</name>
<name>
<surname>Gaawe</surname> <given-names>D.</given-names>
</name>
</person-group> (<year>2016</year>). <article-title>Identification of medicinal plants based on mitochondrial DNA sequences</article-title>. <source>Chin. Traditional Herbal Drugs</source> <volume>47</volume>, <fpage>1791</fpage>&#x2013;<lpage>1796</lpage>. doi: <pub-id pub-id-type="doi">10.7501/j.issn.0253-2670.2016.10.027</pub-id>
</citation>
</ref>
<ref id="B40">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Lv</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Zou</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Zhang</surname> <given-names>X.</given-names>
</name>
<name>
<surname>Liu</surname> <given-names>B.</given-names>
</name>
<name>
<surname>Peng</surname> <given-names>X.</given-names>
</name>
<name>
<surname>Chu</surname> <given-names>C.</given-names>
</name>
</person-group> (<year>2023</year>). <article-title>A review on the chemical constituents and pharmacological efficacies of Lindera aggregata (Sims) Kosterm</article-title>. <source>Front. Nutr.</source> <volume>9</volume>, <elocation-id>1071276</elocation-id>. doi:&#xa0;<pub-id pub-id-type="doi">10.3389/fnut.2022.1071276</pub-id>
</citation>
</ref>
<ref id="B41">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Miao</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Chen</surname> <given-names>H.</given-names>
</name>
<name>
<surname>Xu</surname> <given-names>W.</given-names>
</name>
<name>
<surname>Liu</surname> <given-names>C.</given-names>
</name>
<name>
<surname>Huang</surname> <given-names>L.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>Cistanche Species mitogenomes suggest diversity and complexity in lamiales-order mitogenomes</article-title>. <source>Genes</source> <volume>13</volume>, <fpage>1791</fpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.3390/genes13101791</pub-id>
</citation>
</ref>
<ref id="B42">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Millar</surname> <given-names>A. H.</given-names>
</name>
<name>
<surname>Heazlewood</surname> <given-names>J. L.</given-names>
</name>
<name>
<surname>Kristensen</surname> <given-names>B. K.</given-names>
</name>
<name>
<surname>Braun</surname> <given-names>H.-P.</given-names>
</name>
<name>
<surname>M&#xf8;ller</surname> <given-names>I. M.</given-names>
</name>
</person-group> (<year>2005</year>). <article-title>The plant mitochondrial proteome</article-title>. <source>Trends Plant Sci.</source> <volume>10</volume>, <fpage>36</fpage>&#x2013;<lpage>43</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.tplants.2004.12.002</pub-id>
</citation>
</ref>
<ref id="B43">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Minh</surname> <given-names>B. Q.</given-names>
</name>
<name>
<surname>Schmidt</surname> <given-names>H. A.</given-names>
</name>
<name>
<surname>Chernomor</surname> <given-names>O.</given-names>
</name>
<name>
<surname>Schrempf</surname> <given-names>D.</given-names>
</name>
<name>
<surname>Woodhams</surname> <given-names>M. D.</given-names>
</name>
<name>
<surname>Von Haeseler</surname> <given-names>A.</given-names>
</name>
<etal/>
</person-group>. (<year>2020</year>). <article-title>IQ-TREE 2: new models and efficient methods for phylogenetic inference in the genomic era</article-title>. <source>Mol. Biol. Evol.</source> <volume>37</volume>, <fpage>1530</fpage>&#x2013;<lpage>1534</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1093/molbev/msaa015</pub-id>
</citation>
</ref>
<ref id="B44">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Niu</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Gao</surname> <given-names>C.</given-names>
</name>
<name>
<surname>Liu</surname> <given-names>J.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>Complete mitochondrial genomes of three Mangifera species, their genomic structure and gene transfer from chloroplast genomes</article-title>. <source>BMC Genomics</source> <volume>23</volume>, <fpage>147</fpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1186/s12864-022-08383-1</pub-id>
</citation>
</ref>
<ref id="B45">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ogihara</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Yamazaki</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Murai</surname> <given-names>K.</given-names>
</name>
<name>
<surname>Kanno</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Terachi</surname> <given-names>T.</given-names>
</name>
<name>
<surname>Shiina</surname> <given-names>T.</given-names>
</name>
<etal/>
</person-group>. (<year>2005</year>). <article-title>Structural dynamics of cereal mitochondrial genomes as revealed by complete nucleotide sequencing of the wheat mitochondrial genome</article-title>. <source>Nucleic Acids Res.</source> <volume>33</volume>, <fpage>6235</fpage>&#x2013;<lpage>6250</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1093/nar/gki925</pub-id>
</citation>
</ref>
<ref id="B46">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Parvathy</surname> <given-names>S. T.</given-names>
</name>
<name>
<surname>Udayasuriyan</surname> <given-names>V.</given-names>
</name>
<name>
<surname>Bhadana</surname> <given-names>V.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>Codon usage bias</article-title>. <source>Mol. Biol. Rep.</source> <volume>49</volume>, <fpage>539</fpage>&#x2013;<lpage>565</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1007/s11033-021-06749-4</pub-id>
</citation>
</ref>
<ref id="B47">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Pfeifer</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Martis</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Asp</surname> <given-names>T.</given-names>
</name>
<name>
<surname>Mayer</surname> <given-names>K. F.</given-names>
</name>
<name>
<surname>L&#xfc;bberstedt</surname> <given-names>T.</given-names>
</name>
<name>
<surname>Byrne</surname> <given-names>S.</given-names>
</name>
<etal/>
</person-group>. (<year>2013</year>). <article-title>The perennial ryegrass GenomeZipper: targeted use of genome resources for comparative grass genomics</article-title>. <source>Plant Physiol.</source> <volume>161</volume>, <fpage>571</fpage>&#x2013;<lpage>582</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1104/pp.112.207282</pub-id>
</citation>
</ref>
<ref id="B48">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Romero</surname> <given-names>H.</given-names>
</name>
<name>
<surname>Zavala</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Musto</surname> <given-names>H.</given-names>
</name>
</person-group> (<year>2000</year>). <article-title>Codon usage in Chlamydia trachomatis is the result of strand-specific mutational biases and a complex pattern of selective forces</article-title>. <source>Nucleic Acids Res.</source> <volume>28</volume>, <fpage>2084</fpage>&#x2013;<lpage>2090</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1093/nar/28.10.2084</pub-id>
</citation>
</ref>
<ref id="B49">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ronquist</surname> <given-names>F.</given-names>
</name>
<name>
<surname>Teslenko</surname> <given-names>M.</given-names>
</name>
<name>
<surname>van der Mark</surname> <given-names>P.</given-names>
</name>
<name>
<surname>Ayres</surname> <given-names>D. L.</given-names>
</name>
<name>
<surname>Darling</surname> <given-names>A.</given-names>
</name>
<name>
<surname>H&#xf6;hna</surname> <given-names>S.</given-names>
</name>
<etal/>
</person-group>. (<year>2012</year>). <article-title>MrBayes 3.2: efficient Bayesian phylogenetic inference and model choice across a large model space</article-title>. <source>Systemat Biol.</source> <volume>61</volume>, <fpage>539</fpage>&#x2013;<lpage>542</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1093/sysbio/sys029</pub-id>
</citation>
</ref>
<ref id="B50">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Rozas</surname> <given-names>J.</given-names>
</name>
<name>
<surname>S&#xe1;nchez-Delbarrio</surname> <given-names>J. C.</given-names>
</name>
<name>
<surname>Messeguer</surname> <given-names>X.</given-names>
</name>
<name>
<surname>Rozas</surname> <given-names>R.</given-names>
</name>
</person-group> (<year>2003</year>). <article-title>DnaSP, DNA polymorphism analyses by the coalescent and other methods</article-title>. <source>Bioinformatics</source> <volume>19</volume>, <fpage>2496</fpage>&#x2013;<lpage>2497</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1093/bioinformatics/btg359</pub-id>
</citation>
</ref>
<ref id="B51">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Salleh</surname> <given-names>W.M.N.H.W.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Lindera aggregata (Sims) Kosterm: Review on phytochemistry and biological activities</article-title>. <source>Boletin Latinoamericano y del Caribe Plantas Medicinales y Aromaticas</source> <volume>19</volume>, <fpage>527</fpage>&#x2013;<lpage>541</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.37360/blacpma</pub-id>
</citation>
</ref>
<ref id="B52">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Shi</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Chen</surname> <given-names>Z.</given-names>
</name>
<name>
<surname>Jiang</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Li</surname> <given-names>X.</given-names>
</name>
<name>
<surname>Zeng</surname> <given-names>W.</given-names>
</name>
</person-group> (<year>2024</year>). <article-title>Comparative analysis of chloroplast genomes of &#x201c;Tiantai wu-yao&#x201d;(Lindera aggregata) and taxa of the same genus and different genera</article-title>. <source>Genes</source> <volume>15</volume>, <fpage>263</fpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.3390/genes15030263</pub-id>
</citation>
</ref>
<ref id="B53">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Skippington</surname> <given-names>E.</given-names>
</name>
<name>
<surname>Barkman</surname> <given-names>T. J.</given-names>
</name>
<name>
<surname>Rice</surname> <given-names>D. W.</given-names>
</name>
<name>
<surname>Palmer</surname> <given-names>J. D.</given-names>
</name>
</person-group> (<year>2015</year>). <article-title>Miniaturized mitogenome of the parasitic plant Viscum scurruloideum is extremely divergent and dynamic and has lost all nad genes</article-title>. <source>Proc. Natl. Acad. Sci.</source> <volume>112</volume>, <fpage>E3515</fpage>&#x2013;<lpage>E3524</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1073/pnas.1504491112</pub-id>
</citation>
</ref>
<ref id="B54">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Sloan</surname> <given-names>D. B.</given-names>
</name>
<name>
<surname>Alverson</surname> <given-names>A. J.</given-names>
</name>
<name>
<surname>Chuckalovcak</surname> <given-names>J. P.</given-names>
</name>
<name>
<surname>Wu</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Mccauley</surname> <given-names>D. E.</given-names>
</name>
<name>
<surname>Palmer</surname> <given-names>J. D.</given-names>
</name>
<etal/>
</person-group>. (<year>2012</year>). <article-title>Rapid evolution of enormous, multichromosomal genomes in flowering plant mitochondria with exceptionally high mutation rates</article-title>. <source>PloS Biol.</source> <volume>10</volume>, <fpage>e1001241</fpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1371/journal.pbio.1001241</pub-id>
</citation>
</ref>
<ref id="B55">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Smith</surname> <given-names>D. R.</given-names>
</name>
<name>
<surname>Keeling</surname> <given-names>P. J.</given-names>
</name>
</person-group> (<year>2015</year>). <article-title>Mitochondrial and plastid genome architecture: reoccurring themes, but significant differences at the extremes</article-title>. <source>Proc. Natl. Acad. Sci.</source> <volume>112</volume>, <fpage>10177</fpage>&#x2013;<lpage>10184</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1073/pnas.1422049112</pub-id>
</citation>
</ref>
<ref id="B56">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Sueoka</surname> <given-names>N.</given-names>
</name>
</person-group> (<year>1988</year>). <article-title>Directional mutation pressure and neutral molecular evolution</article-title>. <source>Proc. Natl. Acad. Sci.</source> <volume>85</volume>, <fpage>2653</fpage>&#x2013;<lpage>2657</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1073/pnas.85.8.2653</pub-id>
</citation>
</ref>
<ref id="B57">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Sun</surname> <given-names>T.</given-names>
</name>
<name>
<surname>Bentolila</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Hanson</surname> <given-names>M. R.</given-names>
</name>
</person-group> (<year>2016</year>). <article-title>The unexpected diversity of plant organelle RNA editosomes</article-title>. <source>Trends Plant Sci.</source> <volume>21</volume>, <fpage>962</fpage>&#x2013;<lpage>973</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.tplants.2016.07.005</pub-id>
</citation>
</ref>
<ref id="B58">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Tang</surname> <given-names>D.</given-names>
</name>
<name>
<surname>Wei</surname> <given-names>F.</given-names>
</name>
<name>
<surname>Quan</surname> <given-names>C.</given-names>
</name>
<name>
<surname>Huang</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Huang</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Wei</surname> <given-names>K.</given-names>
</name>
<etal/>
</person-group>. (<year>2022</year>). <article-title>Codon usage bias and evolution analysis in the mitochondrial genome of Mesona chinensis Benth</article-title>. <source>Acta Physiologiae Plantarum</source> <volume>44</volume>, <fpage>118</fpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1007/s11738-022-03453-z</pub-id>
</citation>
</ref>
<ref id="B59">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Tang</surname> <given-names>X.-M.</given-names>
</name>
<name>
<surname>Yang</surname> <given-names>S.-Z.</given-names>
</name>
<name>
<surname>Chen</surname> <given-names>H.-Z.</given-names>
</name>
<name>
<surname>Sun</surname> <given-names>Z.-D.</given-names>
</name>
<name>
<surname>Lai</surname> <given-names>Z.-G.</given-names>
</name>
<name>
<surname>Zeng</surname> <given-names>W.-Y.</given-names>
</name>
<etal/>
</person-group>. (<year>2020</year>). <article-title>Comparative analysis on codon usage Bias inmitogenome of two species in genus Glycine</article-title>. <source>Guihaia.</source> <volume>40</volume>, <fpage>926</fpage>&#x2013;<lpage>934</lpage>. doi: 10.11931/guihaia.gxzw201911033</citation>
</ref>
<ref id="B60">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Tao</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Deng</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Wang</surname> <given-names>P.</given-names>
</name>
</person-group> (<year>2024</year>). <article-title>Traditional uses, phytochemistry, pharmacology, processing methods and quality control of Lindera aggregata (Sims) Kosterm: A critical review</article-title>. <source>J. Ethnopharmacology</source> <volume>318</volume>, <fpage>116954</fpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.jep.2023.116954</pub-id>
</citation>
</ref>
<ref id="B61">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wang</surname> <given-names>L.</given-names>
</name>
<name>
<surname>Liu</surname> <given-names>X.</given-names>
</name>
<name>
<surname>Xu</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Zhang</surname> <given-names>Z.</given-names>
</name>
<name>
<surname>Wei</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Hu</surname> <given-names>Y.</given-names>
</name>
<etal/>
</person-group>. (<year>2024</year>). <article-title>Assembly and comparative analysis of the first complete mitochondrial genome of a traditional Chinese medicine Angelica biserrata (Shan et Yuan) Yuan et Shan</article-title>. <source>Int. J. Biol. Macromol</source> <volume>257</volume>, <fpage>128571</fpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.ijbiomac.2023.128571</pub-id>
</citation>
</ref>
<ref id="B62">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wang</surname> <given-names>H.</given-names>
</name>
<name>
<surname>Sun</surname> <given-names>H.</given-names>
</name>
<name>
<surname>Kwon</surname> <given-names>W.-S.</given-names>
</name>
<name>
<surname>Jin</surname> <given-names>H.</given-names>
</name>
<name>
<surname>Yang</surname> <given-names>D.-C.</given-names>
</name>
</person-group> (<year>2009</year>). <article-title>Molecular identification of the Korean ginseng cultivar &#x201c;Chunpoong&#x201d; using the mitochondrial nad7 intron 4 region</article-title>. <source>Mitochondrial DNA</source> <volume>20</volume>, <fpage>41</fpage>&#x2013;<lpage>45</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1080/19401730902856738</pub-id>
</citation>
</ref>
<ref id="B63">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wang</surname> <given-names>Z.</given-names>
</name>
<name>
<surname>Xu</surname> <given-names>B.</given-names>
</name>
<name>
<surname>Li</surname> <given-names>B.</given-names>
</name>
<name>
<surname>Zhou</surname> <given-names>Q.</given-names>
</name>
<name>
<surname>Wang</surname> <given-names>G.</given-names>
</name>
<name>
<surname>Jiang</surname> <given-names>X.</given-names>
</name>
<etal/>
</person-group>. (<year>2020</year>). <article-title>Comparative analysis of codon usage patterns in chloroplast genomes of six Euphorbiaceae species</article-title>. <source>PeerJ</source> <volume>8</volume>, <elocation-id>e8251</elocation-id>. doi:&#xa0;<pub-id pub-id-type="doi">10.7717/peerj.8251</pub-id>
</citation>
</ref>
<ref id="B64">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wang</surname> <given-names>D.</given-names>
</name>
<name>
<surname>Zhang</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Zhang</surname> <given-names>Z.</given-names>
</name>
<name>
<surname>Zhu</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Yu</surname> <given-names>J.</given-names>
</name>
</person-group> (<year>2010</year>). <article-title>KaKs_Calculator 2.0: a toolkit incorporating gamma-series methods and sliding window strategies</article-title>. <source>Genomics Proteomics Bioinf.</source> <volume>8</volume>, <fpage>77</fpage>&#x2013;<lpage>80</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/S1672-0229(10)60008-3</pub-id>
</citation>
</ref>
<ref id="B65">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Weng</surname> <given-names>M.</given-names>
</name>
<name>
<surname>You</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Luo</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Lin</surname> <given-names>Z.</given-names>
</name>
<name>
<surname>Chen</surname> <given-names>T.</given-names>
</name>
<name>
<surname>Peng</surname> <given-names>X.</given-names>
</name>
<etal/>
</person-group>. (<year>2022</year>). <article-title>Antibacterial mechanism of polysaccharides from the leaves of Lindera aggregata (Sims) Kosterm. by metabolomics based on HPLC/MS</article-title>. <source>Int. J. Biol. Macromol</source> <volume>221</volume>, <fpage>303</fpage>&#x2013;<lpage>313</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.ijbiomac.2022.09.009</pub-id>
</citation>
</ref>
<ref id="B66">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wick</surname> <given-names>R. R.</given-names>
</name>
<name>
<surname>Schultz</surname> <given-names>M. B.</given-names>
</name>
<name>
<surname>Zobel</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Holt</surname> <given-names>K. E.</given-names>
</name>
</person-group> (<year>2015</year>). <article-title>Bandage: interactive visualization of <italic>de novo</italic> genome assemblies</article-title>. <source>Bioinformatics</source> <volume>31</volume>, <fpage>3350</fpage>&#x2013;<lpage>3352</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1093/bioinformatics/btv383</pub-id>
</citation>
</ref>
<ref id="B67">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wright</surname> <given-names>F.</given-names>
</name>
</person-group> (<year>1990</year>). <article-title>The &#x2018;effective number of codons&#x2019; used in a gene</article-title>. <source>Gene</source> <volume>87</volume>, <fpage>23</fpage>&#x2013;<lpage>29</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/0378-1119(90)90491-9</pub-id>
</citation>
</ref>
<ref id="B68">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Xiao</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Hu</surname> <given-names>X.</given-names>
</name>
<name>
<surname>Li</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Liu</surname> <given-names>Q.</given-names>
</name>
<name>
<surname>Shen</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Jiang</surname> <given-names>T.</given-names>
</name>
<etal/>
</person-group>. (<year>2024</year>). <article-title>Comparative analysis of codon usage patterns in the chloroplast genomes of nine forage legumes</article-title>. <source>Physiol. Mol. Biol. Plants</source> <volume>30</volume>, <fpage>153</fpage>&#x2013;<lpage>166</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1007/s12298-024-01421-0</pub-id>
</citation>
</ref>
<ref id="B69">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Xie</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Chen</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Cai</surname> <given-names>G.</given-names>
</name>
<name>
<surname>Cai</surname> <given-names>R.</given-names>
</name>
<name>
<surname>Hu</surname> <given-names>Z.</given-names>
</name>
<name>
<surname>Wang</surname> <given-names>H.</given-names>
</name>
</person-group> (<year>2023</year>). <article-title>Tree Visualization By One Table (tvBOT): a web application for visualizing, modifying and annotating phylogenetic trees</article-title>. <source>Nucleic Acids Res.</source>, <volume>51</volume>, <fpage>W587</fpage>&#x2013;<lpage>W592</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1093/nar/gkad359</pub-id>
</citation>
</ref>
<ref id="B70">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Xie</surname> <given-names>D.-F.</given-names>
</name>
<name>
<surname>Yu</surname> <given-names>H.-X.</given-names>
</name>
<name>
<surname>Price</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Xie</surname> <given-names>C.</given-names>
</name>
<name>
<surname>Deng</surname> <given-names>Y.-Q.</given-names>
</name>
<name>
<surname>Chen</surname> <given-names>J.-P.</given-names>
</name>
<etal/>
</person-group>. (<year>2019</year>). <article-title>Phylogeny of Chinese Allium species in section Daghestanica and adaptive evolution of Allium (Amaryllidaceae, Allioideae) species revealed by the chloroplast complete genome</article-title>. <source>Front. Plant Sci.</source> <volume>10</volume>, <elocation-id>460</elocation-id>. doi:&#xa0;<pub-id pub-id-type="doi">10.3389/fpls.2019.00460</pub-id>
</citation>
</ref>
<ref id="B71">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Yang</surname> <given-names>H.</given-names>
</name>
<name>
<surname>Ni</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Zhang</surname> <given-names>X.</given-names>
</name>
<name>
<surname>Li</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Chen</surname> <given-names>H.</given-names>
</name>
<name>
<surname>Liu</surname> <given-names>C.</given-names>
</name>
</person-group> (<year>2023</year>). <article-title>The mitochondrial genomes of Panax notoginseng reveal recombination mediated by repeats associated with DNA replication</article-title>. <source>Int. J. Biol. Macromol</source> <volume>252</volume>, <fpage>126359</fpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.ijbiomac.2023.126359</pub-id>
</citation>
</ref>
<ref id="B72">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhang</surname> <given-names>X.</given-names>
</name>
<name>
<surname>Chen</surname> <given-names>H.</given-names>
</name>
<name>
<surname>Ni</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Wu</surname> <given-names>B.</given-names>
</name>
<name>
<surname>Li</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Burzy&#x144;ski</surname> <given-names>A.</given-names>
</name>
<etal/>
</person-group>. (<year>2024</year>). <article-title>Plant mitochondrial genome map (PMGmap): A software tool for the comprehensive visualization of coding, noncoding and genome features of plant mitochondrial genomes</article-title>. <source>Mol. Ecol. Resour.</source> <volume>24</volume>, <elocation-id>e13952</elocation-id>. doi: <pub-id pub-id-type="doi">10.1111/1755-0998.13952</pub-id>
</citation>
</ref>
<ref id="B73">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhang</surname> <given-names>D.</given-names>
</name>
<name>
<surname>Gao</surname> <given-names>F.</given-names>
</name>
<name>
<surname>Jakovli&#x107;</surname> <given-names>I.</given-names>
</name>
<name>
<surname>Zou</surname> <given-names>H.</given-names>
</name>
<name>
<surname>Zhang</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Li</surname> <given-names>W. X.</given-names>
</name>
<etal/>
</person-group>. (<year>2020</year>). <article-title>PhyloSuite: An integrated and scalable desktop platform for streamlined molecular sequence data management and evolutionary phylogenetics studies</article-title>. <source>Mol. Ecol. Resour.</source> <volume>20</volume>, <fpage>348</fpage>&#x2013;<lpage>355</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1111/1755-0998.13096</pub-id>
</citation>
</ref>
<ref id="B74">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhang</surname> <given-names>R.</given-names>
</name>
<name>
<surname>Zhang</surname> <given-names>L.</given-names>
</name>
<name>
<surname>Wang</surname> <given-names>W.</given-names>
</name>
<name>
<surname>Zhang</surname> <given-names>Z.</given-names>
</name>
<name>
<surname>Du</surname> <given-names>H.</given-names>
</name>
<name>
<surname>Qu</surname> <given-names>Z.</given-names>
</name>
<etal/>
</person-group>. (<year>2018</year>). <article-title>Differences in codon usage bias between photosynthesis-related genes and genetic system-related genes of chloroplast genomes in cultivated and wild solanum species</article-title>. <source>Int. J. Mol. Sci.</source> <volume>19</volume>, <fpage>3142</fpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.3390/ijms19103142</pub-id>
</citation>
</ref>
<ref id="B75">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhao</surname> <given-names>M.-L.</given-names>
</name>
<name>
<surname>Song</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Ni</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Yao</surname> <given-names>X.</given-names>
</name>
<name>
<surname>Tan</surname> <given-names>Y.-H.</given-names>
</name>
<name>
<surname>Xu</surname> <given-names>Z.-F.</given-names>
</name>
</person-group> (<year>2018</year>). <article-title>Comparative chloroplast genomics and phylogenetics of nine Lindera species (Lauraceae)</article-title>. <source>Sci. Rep.</source> <volume>8</volume>, <fpage>8844</fpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1038/s41598-018-27090-0</pub-id>
</citation>
</ref>
<ref id="B76">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhu</surname> <given-names>L.</given-names>
</name>
<name>
<surname>Xian</surname> <given-names>F.-J.</given-names>
</name>
<name>
<surname>Zhang</surname> <given-names>Q.-N.</given-names>
</name>
<name>
<surname>Hu</surname> <given-names>J.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>Research progress of RNA editing</article-title>. <source>Biotechnol. Bull.</source> <volume>38</volume>, <fpage>1</fpage>, <fpage>1</fpage>&#x2013;<lpage>14</lpage>. doi: <pub-id pub-id-type="doi">10.13560/j.cnki.biotech.bull.1985.2021-1326</pub-id>
</citation>
</ref>
</ref-list>
</back>
</article>