Assembly and comparative analysis of the first complete mitochondrial genome of Salix psammophila, a good windbreak and sand fixation shrub

Abstract

Salix psammophila, commonly known as the sandlive willow, is a vital shrub species within the Salicaceae family, particularly significant for its ecological role in regions susceptible to desertification and sandy soils. In this study, we assembled the complete S. psammophila mitochondrial genome using Pacbio HiFi third-generation sequencing data. The genome was found to be a typical single circular structure, with a total length of 715,555 bp and a GC content of 44.89%. We annotated 33 unique protein-coding genes (PCGs), which included 24 core mitochondrial genes and 9 variable genes, as well as 18 tRNA genes (5 of which were multicopy genes) and 3 rRNA genes. Comparative analysis of the PCGs from the mitochondrial genomes of S. psammophila, Populus deltoides, Populus simonii, Salix wilsonii, and Salix suchowensis revealed that these genes are relatively conserved within the Salicaceae family, with variability primarily occurring in the ribosomal protein genes. The absence of the rps14, which encodes a ribosomal protein, may have played a role in the evolution of stress tolerance in Salicaceae plants. Additionally, we identified 232 SSRs, 19 tandem repeat sequences, and 236 dispersed repeat sequences in the S. psammophila mitochondrial genome, with palindromic and forward repeats being the most abundant. The longest palindromic repeat measured 260 bp, while the longest forward repeat was 86,068 bp. Furthermore, 324 potential RNA editing sites were discovered, all involving C-to-U edits, with the nad4 having the highest number of edits. These findings provide valuable insights into the phylogenetic and genetic research of Salicaceae plants.

Introduction

Salix psammophila is a perennial fast-growing and multi-resistance shrub in the genus Salix of the Salicaceae family, naturally distributed in arid and semi-arid desert areas, produced in the Maowusu Sandland and the Kubuqi Desert, and is the first choice of species for the prevention and control of desertification areas (Bao and Zhang, 2012; Jia et al., 2019; Lu et al., 2019; Yang et al., 2022). In the S. psammophila asexual line, the 17–38 asexual line is an excellent asexual variety, and in the morphology of the performance of the upright, short branchlets, the landscape effect is good. Currently, research on S. psammophila mainly focuses on strain type and sand barriers, while the genetic resources and phylogenetic status of S. psammophila mitochondria are still unclear. In order to further develop and utilize S. psammophila, the 17–38 asexual line of S. psammophila was used as a germplasm resource to assemble and annotate the S. psammophila mitochondrial genome and to analyze the structural characteristics and phylogenetic status of the S. psammophila mitochondrial genome.

The mitochondria are organelles with semi-autonomous genetic characteristics, relatively independent of the nucleus, with a relatively independent genetic system (Gray, 2012; Chen et al., 2017). It is an important organelle in eukaryotic cells that produce ATP by oxidative phosphorylation and is a source of energy for a variety of biochemical processes, playing a crucial role in plant development and reproduction (Wu et al., 2019). Plant mitochondrial genomes have a variety of structures due to the possession of a large number of repetitive sequences, and assembled mitochondrial genomes may be unicyclic, polycyclic, linear, or multibranching in structure (Ye et al., 2017; Cheng et al., 2021). For example, the mitochondrial genome of Populus deltoides and P. simonii has been assembled into three circular chromosomes. In contrast, the mitochondrial genome of P. tremula and P. tremula × P. alba is each organized into a single circular chromosome (Kersten et al., 2016; Bi et al., 2022; Qu et al., 2022). In addition, the size of plant mitochondrial genomes varies greatly, for instance, the small 66-kb Viscum scurruloideum (mistletoe) and the large 11.3-Mb Silene conica (Sloan et al., 2012; Skippington et al., 2015). It has been shown that the main cause of mitochondrial genome amplification in plants is the insertion of foreign sequences and the presence of a large number of repetitive sequences (Richardson and Palmer, 2006; Li, 2009; Smith, 2011). The number of genes in plant mitogenomes also varied, from 19 in the Viscum album to 221 in the pepper (Capsicum annuum) (Jo et al., 2011; Petersen et al., 2015).

Although the size, structure, and number of plant mitogenomes vary greatly, the encoded genes are relatively conserved, such as ATP synthase genes, NADH dehydrogenase genes, cytochrome C reductase genes, and cytochrome C oxidase genes (Bi et al., 2020). In addition, there are protein-coding genes (PCGs) that are lost during the evolution of the plant mitochondrial genome, such as the succinate dehydrogenase genes and the ribosomal protein (RP) genes (Bi et al., 2016). The genes rps12, sdh3, and sdh4 appear to be absent in the monocotyledonous plants, whereas the rps2 gene is notably missing in dicotyledonous plants (Zhang et al., 2012). RPs have important roles in metabolism, cell division, and growth processes (Wang et al., 2013a). RP genes are among the most highly expressed genes in most cells, and their expressions are differentially regulated by environmental and signaling molecules (Vemanna et al., 2019; Petibon et al., 2020; Fakih et al., 2023). Wang et al. explored the changes of ribosomes in phosphate- and iron-deficient Arabidopsis roots and found that most RP genes were related to translation and ribosome assembly, and some were related to low temperature, UV-B, and salt stress (Wang et al., 2013a).

It is well known that plant mitochondrial genomes vary greatly in size, structure, incorporation of exogenous DNA, and mutation rates, and there are a large number of repetitive sequences in the genomes, making it difficult to sequence plant mitochondrial genomes (Petersen et al., 2020). However, with the rapid development of third-generation sequencing technology, especially the emergence of PacBio and Oxford Nanopore, more and more plant mitochondrial genomes have been assembled and submitted to the NCBI database. In this study, we obtained, assembled, and annotated the complete S. psammophila mitochondrial genome for the first time based on Pacbio HiFi third-generation sequencing data. Characteristics of the S. psammophila mitochondrial genome were analyzed, including codon preference, RNA editing events, repetitive sequence analysis, phylogenetic analysis, and sequence migration analysis. These results will further lay the foundation for understanding and studying the evolution and inheritance of mitochondria in Salicaceae plants.

Materials and methods

Plant materials, DNA extraction, and sequencing

The experimental material used in this study was mature leaves of the 3-year-old 17–38 asexual clone of S. psammophila from the Germplasm Resources gene bank of S. psammophila in Ordos Dalad, the Inner Mongolia Autonomous Region of China (E 110°38′59.1″, N 40°14′15.5″). The leaves were collected, snap frozen in liquid nitrogen, and stored in a −80°C refrigerator. The total genomic DNA was extracted using a plant DNA extraction kit and sequenced with Pacbio HiFi sequencing technology following the manufacturer’s kits and protocols (Zhang et al., 2012, 2024).

Assembly and annotation of the mitochondrial genome

We used Pacbio HiFi third-generation sequencing data to assemble the S. psammophila mitochondrial genome. The long-reads data were assembled using the default parameters of the Flye software (v2.9.4) (Kolmogorov et al., 2019) to obtain a graphical mitochondrial genome in GFA format. Use makeblastdb to build a library of all fasta-formatted contigs assembled. Then, the contig fragments containing the mitochondrial genome were identified using the BLASTn program (Chen et al., 2015) and using the plant mitochondrial genome from Arabidopsis thaliana as a reference genome with the parameters “-evalue 1e-5 -outfmt 6 -max_hsps 10 -word_size 7 -task blastn-short”. In addition, we also included two closely related species, S. wilsonii and S. suchowensis, as references for re-assembly. This mitochondrial genome was visualized by Bandage software (v0.8.1) (Wick et al., 2015) and screened for mitochondrial contigs based on BLASTn results to obtain a sketch of the mitochondrial genome of S. psammophila.

A. thaliana (NC_037304), Liriodendron tulipifera (NC_021152.1), P. simonii (MZ905370), and S. wilsonii (NC_064688.1) were used as reference genomes, and the mitochondrial genome of S. psammophila was annotated using Geseq software (v2.03) (Tillich et al., 2017). The tRNA and rRNA of the genome were annotated using Tr4Nascan-SE software (v.2.0.11) (Lowe and Eddy, 1997) and BLASTn software (v2.13.0) (Chen et al., 2015). Finally, each mitochondrial genome annotation was manually corrected using Apollo software (v1.11.8) (Lewis et al., 2002).

Structural analysis and codon preference analysis of the genome

To resolve the obtained repetitive regions in the graphical mitochondrial genome, the long reads were compared to the repetitive sequences to determine whether any long reads crossed the repetitive region, especially the ctg3 repetitive region. Selection of suitable pathways is based on long reads to hypothesize the most likely mitochondrial genome structure formed in S. psammophila.

Genomic PCGs were extracted by the Phylosuite software (v1.1.16) (Zhang et al., 2020) and then analyzed for codon preference by the Mega software (v7.0) (Kumar et al., 2016), and RSCU values were calculated.

Sequence repeat analysis

Microsatellite repeats, tandem repeats, and dispersed repeats were identified using MISA (v2.1) (https://webblast.ipk-gatersleben.de/misa/) (Beier et al., 2017), TRF (v4.09) (https://tandem.bu.edu/trf/trf.unix.help.html) (Benson, 1999), and the REPuter web server (https://bibiserv.cebitec.uni-bielefeld.de/reputer/) (Kurtz, 2001), and the results were analyzed using Excel software and the Circos package (v0.69.9) (Zhang et al., 2013) for visualization.

Sequence transfer analysis and synteny analysis

The chloroplast genome was assembled and annotated using the GetOrganelle software (Jin et al., 2020) and CPGAVAS2 software (Shi et al., 2019). The annotation results were subsequently corrected using the CPGView software (Liu et al., 2023). Homologous fragments of the mitochondrial and chloroplast genomes were analyzed using the BLASTn software (v2.13.0) (Chen et al., 2015), and the results were visualized using the Circos package (v0.69.9) (Zhang et al., 2013).

Six tree species, including S. psammophila, were selected for analysis, and based on the BLASTn results of two–two comparisons between individual mitochondrial genomes, homologous sequences of more than 500 bp were retained as conserved covariate blocks, and Multiple Synteny Plot was plotted.

Prediction of RNA editing sites

We used the sequences of all PCGs encoded by the S. psammophila mitochondrial genome as the input file and used Deepred-mt (Edera et al., 2021) to predict RNA editing sites from C to U of the mitochondrial PCGs, and retained all results with probability values greater than 0.9.

Phylogenetic analysis

Thirty-four mitochondrial genomes from three families, namely, Salicaceae, Euphorbiaceae, and Rhizophoraceae, were selected for phylogenetic analysis based on their affinities, with two Rhizophoraceae mitochondrial genomes set as the outgroup. Common genes were extracted using the PhyloSuite software (v1.1.16) (Zhang et al., 2020). Multiple sequence alignment was performed using MAFFT software (v7.505) (Katoh and Standley, 2013). Then, phylogenetic trees were constructed using the IQ-TREE software (v1.6.12) (Nguyen et al., 2015) based on the maximum likelihood (ML) method with the parameters “–alrt 1000 -B 1000” and the analysis results were visualized using ITOL software (v6) (Letunic and Bork, 2019).

Results

Assembly and annotation of the mitochondrial genome

A total of 265,714 reads were utilized in assembling the S. psammophila mitochondrial genome, with the longest read length reaching 48,394 bp, an average read length of 18,699 bp, an N50 read length of 18,657 bp, and total read bases amounting to 4,968,733,600 bp. After exclusion of duplicated regions, one molecular circular was obtained with a total length of 715,555 bp and 44.89% GC content (Figure 1). We annotated the mitochondrial genome of the S. psammophila and the annotation results are shown in Table 1. A total of 33 unique PCGs were annotated, including 24 unique core genes and 9 variable genes, 18 tRNA genes (5 of which were multicopy genes), and 3 rRNA genes. The core genes include five ATP synthase genes (atp1, atp4, atp6, atp8, and atp9), nine NADH dehydrogenase genes (nad1, nad2, nad3, nad4, nad4L, nad5, nad6, nad7, and nad9), four cytochrome C biogenesis genes (ccmB, ccmC, ccmFC, and ccmFN), three cytochrome C oxidase genes (cox1, cox2, and cox3), one protein transport subunit gene (mttB), one maturases gene (matR), and one cytochrome b gene (cob). Variable genes include three RP large subunit genes (rpl2, rpl10, and rpl16), five RP small subunit genes (rps1, rps3, rps4, rps7, and rps12), and one succinate dehydrogenase (sdh4).

Figure 1

The structure of the mitochondrial genome

A sketch of the assembled S. psammophila mitochondrial genome was visualized using Bandage software (v0.8.1), and the results are shown in Figure 2A. The sketch contains three nodes (see Supplementary Table S2A for details). Each node represents a sequence (contig) obtained by assembly, and two nodes are connected by a black line, which represents the overlapping region between the two sequences. Together, all sequences form a complex multibranched closed genome structure that represents the complete mitochondrial genome sequence of S. psammophila. We use long reads to address several key nodes where branching exists, derive the relevant sequences that are at branching nodes, and map them to long reads. Two sequences connected by a black line appear on the same long reads and are connected at the beginning and end, which means that the long reads support the interconnection of these two sequences. A number of different connections exist on different branch nodes, but we prioritize connections that are supported by more long reads when assembling. The resolution of branching nodes due to repeated sequences (ctg3) based on long reads yields a circular contig sequence (Figure 2B), the resolution of which is summarized in Supplementary Table S2B. In addition, ctg3 can be formed into a circular structure with ctg1 and ctg2, respectively, to form a potential secondary conformation of the two circular (Figure 2C).

Figure 2

Codon preference of the mitochondrial genome

Codon preference analysis was performed on 33 unique PCGs of S. psammophila mitochondria, and the codon usage of individual amino acids is shown in Supplementary Table S3. All amino acids except Met (AUG) and Trp (UGG) consist of multiple codons. Codons with relative synonymous codon use (RSCU) greater than 1 are considered to be used preferentially by amino acids. As shown in Supplementary Figure S1, codon usage preferences are also prevalent in mitochondrial PCGs, except for the start codons AUG and Trp (UGG), which have an RSCU value of 1. For example, the termination codon has a high codon usage preference for UAA with an RSCU value of 1.69, which is the highest among mitochondrial PCGs. This is followed by alanine (Ala), which has a high codon usage preference for GCU with an RSCU value of 1.64.

Sequence repeat analysis

Repeat sequence analysis showed that a total of 232 SSRs were identified in the mitochondrial genome of the S. psammophila, with 86, 36, 24, 74, 9, and 3 monomers, dimers, trimers, tetramers, pentamers, and hexamers, respectively (Figure 3A). Among them, monomeric and dimeric SSRs accounted for 52.59% of the total SSRs. Thymine (T) monomeric repeat sequences accounted for 55.81% (48) of the 86 monomeric SSRs. Repeat sequences are categorized into tandem repeat sequences and dispersed repeat sequences. Tandem repeat sequences, also known as satellite DNA, are core repetitive units of approximately 7 to 200 bases that are repeated many times in tandem and are widely found in prokaryotic and eukaryotic genomes. There were 19 tandem repeats in the mitochondrial genes with greater than 85% match and lengths ranging from 12 to 28 bp, and 362 pairs of dispersed repeats with lengths greater than or equal to 30 bp. The dispersed repeats contained 180 pairs of palindromic repeats, 179 pairs of forward repeats and 2 pairs of complementary repeats (Figure 3B). Of these, the longest palindromic repeat is 260 bp and the longest forward repeat is 86,068 bp.

Figure 3

Sequence transfer analysis

During mitochondrial evolution, some fragments migrate from chloroplasts to mitochondrial genes, and the length and sequence similarity of the migrated fragments change from species to species. Based on the sequence similarity, a total of 38 fragments from the chloroplast genome were found in the mitochondrial genome of S. psammophila, with a total length of 16,536 bp, which accounted for 2.31% of the total length of the mitochondrial genome (Figure 4). Among these homologous fragments, MTPT 16 was the longest, with a length of 2,579 bp, and MTPT 8 and MTPT 9 were the shortest, both at 33 bp. Annotating these homologous sequences, we identified eight complete genes, namely, one PCG (atpE) and seven tRNA genes (trnD-GUC, trnM-CAU, trnN-GUU, trnP-UGG, trnS-GGA, trnV-UAC, and trnW-CCA).

Figure 4

Prediction of RNA editing sites

RNA editing events were predicted for 33 unique PCGs of the mitochondrial genome by setting the following criterion: cutoff value = 0.9. A total of 324 potential RNA editing sites were identified, all of which were C-to-U editing and basically appeared in the second or third position, with a few occurring in the third position (Supplementary Table S8). As shown in Figure 5, among these mitochondrial genes, the most edited gene was the nad4 gene, with 38 RNA editing sites identified. The least edited genes were atp1, cox2, rpl10, and rps4, with 1 RNA editing site identified.

Figure 5

Phylogenetic analysis

To further explore the evolutionary relationships of the mitochondria of S. psammophila, phylogenetic trees were constructed for 34 species of plants from three families, namely, Salicaceae, Euphorbiaceae, and Rhizophoraceae, based on the DNA sequences of 24 conserved mitochondrial PCGs, in which Bruguiera sexangula and Kandelia obovata of the Rhizophoraceae family were set up as the outgroups (Figure 6). The PCGs shared by these species are atp1, atp4, atp6, atp8, atp9, ccmB, ccmC, ccmFC, ccmFN, cob, cox1, cox2, cox3, matR, mttB, nad1, nad2, nad3, nad4, nad4L, nad5, nad6, nad7, and nad9. The topology of the mitochondrial DNA phylogeny coincides with the latest classification of the Angiosperm Phylogeny Group (APG). S. psammophila belongs to the family Salicaceae in the order Salicales and is most closely related to S. polaris.

Figure 6

Synteny analysis

S. koriyanagi, S. viminalis, S. suchowensis, S. integra, S. psammophila, and S. polaris were selected for the synteny analysis. As shown in Figure 7, there are a large number of homologous collinear blocks between the mitochondrial genomes of S. psammophila and the other five willow species. In addition, there are some gap regions where sequences are species-specific and have no homology with the rest of the species. The results showed that the mitochondrial genome of S. psammophila underwent genomic rearrangement with S. integra but not with S. polaris.

Figure 7

Discussion

Complexity and variability mitochondrial genome size and structure in plants

The plant mitochondrial genome is an evolutionary dynamic entity with great plant diversity in terms of size, structure, gene content, etc (Petersen et al., 2020). S. psammophila mitochondrial genome is a circular structure with a total length of 715,555 bp and a GC content of 44.89%. The mitochondrial genome structure and GC content were similar to those of the known P. deltoides, P. simonii, S. wilsonii, and S. suchowensis, all of which were cyclic in structure, and the GC content was approximately 44.80%. However, there are also some differences in the genome structure and size of the Salicaceae plants. The mitochondrial genomes of both P. deltoides and P. simonii are three circulars, the P. tremula, P. tremula × P. alba, P. alba, and P. davidiana are the typical single circular sequences, as are the S. suchowensis, S. wilsonii, and S. purpurea (Kersten et al., 2016; Choi et al., 2017; Ye et al., 2017; Brenner et al., 2019; Chen et al., 2020; Bi et al., 2022; Han et al., 2022; Qu et al., 2022). In the genus Populus, the mitochondrial genome size is very conserved, and poplars such as P. simonii (781,478 bp), P. deltoides (802,637 bp), P. tremula (783,442 bp), and P. alba (838,420 bp) are essentially similar in mitochondrial genome size (Kersten et al., 2016; Brenner et al., 2019; Bi et al., 2022; Qu et al., 2022). Numerous studies have shown that a large number of repetitive units, the incorporation of foreign sequences, and the acquisition or loss of large genetic fragments of their own are the main causes of large differences in plant mitochondrial genomes (Wu et al., 2022). In the mitochondrial genome of the genus Populus, there are no repetitive sequences >350 bp, which may be the reason for the small differences in size and structure of the genus Populus mitochondrial genome (Bi et al., 2022). In contrast, the mitochondrial genome size of the genus Salix is highly variable, ranging from 599 kb (S. purpurea) to 735 kb (S. cardiophylla) (Wei et al., 2016; Chen et al., 2020). The size of plant mitochondrial genomes varies dramatically between families and even between different species of the same genus. For example, the mitochondrial genome size difference between woody cotton (Gossypium arboreum) and upland cotton (Gossypium hirsutum) is 65.6 kb (Ala et al., 2023). Furthermore, there was no positive correlation between genome size and gene number; P. deltoides mitochondrial genomes had the largest size but not the largest number of genes, with the largest number of genes being found in S. brachista, and the difference in the number of genes between their mitochondrial genomes was small (Supplementary Table S11). The mitochondrial genome of Welwitschia mirabilis is more than twice the size of Cycas taitungensis, but had fewer genes than Cycas taitungensis (Guo et al., 2016), Meanwhile, the sizes of the S. psammophila and S. wilsonii mitochondrial genomes are similar, and there are significant differences in the types and copy numbers of tRNA, suggesting that different species within the genus Salix possess their unique biological functions. Previous studies have demonstrated that the phylogeny of the mitogenome in Populus can effectively distinguish between different groups of poplar trees (Qu et al., 2022). Similarly, species within the genus Salix can also be categorized using phylogenetic trees. For instance, S. wilsonii and S. dunnii are classified under Sect. Wilsonianae, while S. brachista and S. paraflabellaris are grouped under Sect. Lindleyanae.

Repeat sequences in the mitochondrial genome of Salicaceae plants

A large number of repetitive sequences are distributed in the plant mitochondrial genome, and their size and arrangement are closely related to the differences and rearrangements among the plant mitochondrial genomes, which play an important role in the evolution of plant mitochondria (Tanaka et al., 2014; Wynn and Christensen, 2019; Yang et al., 2021). Repetitive sequences in plant mitochondrial genomes can be categorized into two groups, dispersed repeats and tandem repeats (Wu et al., 2019). Simple sequence repeat (SSR) is a highly variable tandem repeat sequence of 1–6 bp, commonly used as a molecular marker for species identification and genetic diversity studies (Ma et al., 2017; Li and Ye, 2020). A total of 232 SSRs were detected in the S. psammophila mitochondrial genome, and the analysis revealed that 95% of the mononucleotide SSRs were A/T, which was similar to the mitochondrial genomes of other genus Salix, with the highest proportion of A/T repeats among the mononucleotide repeats (Ye et al., 2017; Han et al., 2022). In addition to tandem repeat sequences, plant mitochondrial genomes are characterized by the presence of a large number of dispersed repeat sequences, including palindromic repeat sequences, forward repeat sequences, reverse repeat sequences, and complementary repeat sequences. There are 362 dispersed repeat sequences in the S. psammophila mitochondrial genome, with palindromic and forward repeats being the most abundant. As shown in Supplementary Table S6, a total of 236 (65.19%) dispersed repeats range from 30 to 39 bp, 108 (29.83%) repeats are between 40 and 99 bp in length, 4.7% of the repeats are within the range of 100 to 300 bp in length, and one of the repeats is 86,068 bp, which is a forward repeat. Larger repetitive sequences are critical in genomic structural changes, and paired forward and reverse repeat sequences (>1 kb) can produce subgenomic or isomeric conformations (Bi et al., 2020). In addition to S. psammophila, S. cardiophylla, S. paraflabellaris, and S. suchowensis mitochondrial genomes also contain repetitive sequences greater than 1 kb in length, and their impact on the structure and size expansion of mitochondrial genomes is of interest (Chang et al., 2013; Dong et al., 2018; Chen et al., 2020).

The rps14 function during the evolution of the genome of Salicaceae plants

The PCG analysis based on the mitochondrial genome of S. psammophila and S. suchowensis, S. wilsonii, P. simonii, and P. deltoides revealed that the mitochondrial genome of Salicaceae plants, like other plants, is relatively conserved. It has been shown that the number of PCGs gradually decreases during the evolution of species, and that the genes lost in most cases are RP genes and succinate dehydrogenase genes (Adams and Palmer, 2003). Through comparative analysis of the mitochondrial genomes of several Salicaceae plants, it was found that the rps14 gene was lost in S. psammophila, S. suchowensis, S. wilsonii, and P. simonii (Ye et al., 2017; Bi et al., 2022; Han et al., 2022; Qu et al., 2022). RPs are proteins in ribosomes that play a crucial role in ribosome biogenesis, protein synthesis, cell growth, development, and apoptosis (Ramakrishnan and White, 1998; Naora and Naora, 1999; Maguire and Zimmermann, 2001; Wang et al., 2013b). S. psammophila and P. simonii play important roles in resisting wind damage and fixing sand dunes and have excellent characteristics such as drought resistance, tolerance to barrenness, and strong adaptability. They are the main tree species in the Northwest and North China regions (Jia et al., 2019; Wu et al., 2020; Yang et al., 2022). S. suchowensis is mainly distributed in central and eastern China, and has the function of windbreak and sand fixation (Jia et al., 2020; Wang et al., 2024). The original ecological area of P. deltoides is from the southeastern United States to southern Canada and plays an important role in riparian ecosystems (Fahrenkrog et al., 2017). Based on their different geographical distribution and physiological characteristics, we speculate that the loss of rps14 in the evolution of Salicaceae plants may lead to their sensitivity to adversities such as drought (Supplementary Table S12). RPs are important for the translation of diverse proteins and are involved in overall cellular adaptation under stress conditions (Vemanna et al., 2019). It has been shown that 46 RPs are downregulated for expression under drought conditions (Narsai et al., 2013). Bi et al. found that rps14 and the ccmB, ccmFC, rps1, and rps10 genes may play an important role in the evolutionary process by analyzing selection pressures on PCGs of the Phaseolus vulgaris mitochondrial genome (Bi et al., 2020).

Conservation of RNA editing events in plant mitochondrial genomes

The oxidative phosphorylation system in mitochondria consists of five complexes (I–V), the normal assembly of which requires standard processing of pre-mRNAs, involving RNA editing and intron shearing, and is important for the maintenance of mitochondrial functions (Dudkina et al., 2006; Li et al., 2014; Sun et al., 2018). RNA editing is a phenomenon in which genetic information is altered through the insertion of nucleotides, or nucleotide deletions and transitions, into the messenger RNA of a functional gene (Keller et al., 1999). RNA editing can occur in a variety of forms, including C-to-U, U-to-C, and A-to-I conversions; U insertions and deletions; and G insertions (Small et al., 2020). In plants, RNA editing events occur mainly in plant mitochondria and chloroplasts, and take place mainly in the form of a C-to-U transition (Hao et al., 2021). We analyzed the 33 PCGs within the mitochondrial genome of S. psammophila and identified a total of 324 RNA editing sites. All of these sites exhibited C to U transitions, which is consistent with findings in S. suchowensis and P. deltoides (Ye et al., 2017; Qu et al., 2022). In addition, we found the highest number of edits to nad4 in S. psammophila and S. suchowensis mitochondria. It has also been found in plants such as Corydalis saxicola, with the most editing of the nad4 in the mitochondrial genome (Maldonado et al., 2022; Li et al., 2024). The nad4 is one of the components of complex I, which is the largest complex related to respiration in the mitochondria of land plants and is important for maintaining plant growth and development (Møller et al., 2021). Thus, extensive editing of the nad4 may be important for the functioning of complex I in plant mitochondria.

Conclusions

In this study, we successfully assembled and annotated the mitochondrial genome of S. psammophila for the first time, identifying it as a single circular structure with a total length of 715,555 bp and a GC content of 44.89%. Comparative analysis with mitochondrial genomes of S. wilsonii, S. suchowensis, P. simonii, and P. deltoides revealed that PCGs within the Salicaceae family are highly conserved, with significant variations primarily observed in RPs. Notably, rps14 is absent in S. psammophila, a species well-adapted to arid environments. This absence might indicate a role in enhancing tolerance to environmental stresses such as drought. These findings contribute valuable insights into the phylogenetic and genetic research of Salicaceae plants.

References

• 1

  AdamsK.PalmerD. (2003). Evolution of mitochondrial gene content: gene loss and transfer to the nucleus. Mol. Phylogenet. Evol.29, 380–395. doi: 10.1016/S1055-7903(03)00194-5

  • CrossRef
  • Google Scholar

• 2

  AlaK. G.ZhaoZ.NiL.WangZ. (2023). Comparative analysis of mitochondrial genomes of two alpine medicinal plants of Gentiana (Gentianaceae). PloS One18, e0281134. doi: 10.1371/journal.pone.0281134

  • CrossRef
  • Google Scholar

• 3

  BaoY.ZhangG. (2012). Study of adsorption characteristics of methylene blue onto activated carbon made by Salix Psammophila. Energy Proc.16, 1141–1146. doi: 10.1016/j.egypro.2012.01.182

  • CrossRef
  • Google Scholar

• 4

  BeierS.ThielT.MünchT.ScholzU.MascherM. (2017). MISA-web: a web server for microsatellite prediction. Bioinformatics33, 2583–2585. doi: 10.1093/bioinformatics/btx198

  • CrossRef
  • Google Scholar

• 5

  BensonG. (1999). Tandem repeats finder: a program to analyze DNA sequences. Nucleic Acids Res.27, 573–580. doi: 10.1093/nar/27.2.573

  • CrossRef
  • Google Scholar

• 6

  BiC.LuN.XuY.HeC.LuZ. (2020). Characterization and analysis of the mitochondrial genome of common bean (Phaseolus vulgaris) by Comparative Genomic Approaches. Int. J. Mol. Sci.21, 3778. doi: 10.3390/ijms21113778

  • CrossRef
  • Google Scholar

• 7

  BiC.PatersonA. H.WangX.XuY.WuD.QuY.et al. (2016). Analysis of the complete mitochondrial genome sequence of the diploid cotton Gossypium raimondii by comparative genomics approaches. BioMed. Res. Int.2016, 1–18. doi: 10.1155/2016/5040598

  • CrossRef
  • Google Scholar

• 8

  BiC.QuY.HouJ.WuK.YeN.YinT. (2022). Deciphering the multi-chromosomal mitochondrial genome of Populus simonii. Front. Plant Sci.13. doi: 10.3389/fpls.2022.914635

  • CrossRef
  • Google Scholar

• 9

  BrennerW. G.MaderM.MüllerN. A.HoenickaH.SchroederH.ZornI.et al. (2019). High level of conservation of mitochondrial RNA editing sites among Four Populus Species. G3 (Bethesda)9, 709–717. doi: 10.1534/g3.118.200763

  • CrossRef
  • Google Scholar

• 10

  ChangS.WangY.LuJ.GaiJ.LiJ.ChuP.et al. (2013). The mitochondrial genome of soybean reveals complex genome structures and gene evolution at intercellular and phylogenetic levels. PloS One8, e56502. doi: 10.1371/journal.pone.0056502

  • CrossRef
  • Google Scholar

• 11

  ChenY.YeW.ZhangY.XuY. (2015). High speed BLASTN: an accelerated MegaBLAST search tool. Nucleic Acids Res.43, 7762–7768. doi: 10.1093/nar/gkv784

  • CrossRef
  • Google Scholar

• 12

  ChenX.ZhangL.HuangY.ZhaoF. (2020). Mitochondrial genome of Salix cardiophylla and its implications for infrageneric division of the genus of Salix. Mitochondrial DNA B Resour.5, 3485–3486. doi: 10.1080/23802359.2020.1827065

  • CrossRef
  • Google Scholar

• 13

  ChenZ.ZhaoN.LiS.GroverC. E.NieH.WendelJ. F.et al. (2017). Plant mitochondrial genome evolution and cytoplasmic male sterility. CRC Crit. Rev. Plant Sci.36, 55–69. doi: 10.1080/07352689.2017.1327762

  • CrossRef
  • Google Scholar

• 14

  ChengY.HeX.PriyadarshaniS. V. G. N.WangY.YeL.ShiC.et al. (2021). Assembly and comparative analysis of the complete mitochondrial genome of Suaeda glauca. BMC Genomics22, 167. doi: 10.1186/s12864-021-07490-9

  • CrossRef
  • Google Scholar

• 15

  ChoiM. N.HanM.LeeH.ParkH.-S.KimM.-Y.KimJ.-S.et al. (2017). The complete mitochondrial genome sequence of Populus davidiana dode. Mitochondrial DNA B Resour.2, 113–114. doi: 10.1080/23802359.2017.1289346

  • CrossRef
  • Google Scholar

• 16

  DongS.ZhaoC.ChenF.LiuY.ZhangS.WuH.et al. (2018). The complete mitochondrial genome of the early flowering plant Nymphaea colorata is highly repetitive with low recombination. BMC Genomics19, 614. doi: 10.1186/s12864-018-4991-4

  • CrossRef
  • Google Scholar

• 17

  DudkinaN. V.HeinemeyerJ.SunderhausS.BoekemaE. J.BraunH.-P. (2006). Respiratory chain supercomplexes in the plant mitochondrial membrane. Trends Plant Sci.11, 232–240. doi: 10.1016/j.tplants.2006.03.007

  • CrossRef
  • Google Scholar

• 18

  EderaA. A.SmallI.MiloneD. H.Sanchez-PuertaM. V. (2021). Deepred-Mt: deep representation learning for predicting C-to-U RNA editing in plant mitochondria. Comput. Biol. Med.136, 104682. doi: 10.1016/j.compbiomed.2021.104682

  • CrossRef
  • Google Scholar

• 19

  FahrenkrogA. M.NevesL. G.ResendeM. F. R.DervinisC.DavenportR.BarbazukW. B.et al. (2017). Population genomics of the eastern cottonwood ( Populus deltoides ). Ecol. Evol.7, 9426–9440. doi: 10.1002/ece3.3466

  • CrossRef
  • Google Scholar

• 20

  FakihZ.PlourdeM. B.GermainH. (2023). Differential participation of plant ribosomal proteins from the small ribosomal subunit in protein translation under stress. Biomolecules13, 1160. doi: 10.3390/biom13071160

  • CrossRef
  • Google Scholar

• 21

  GrayM. W. (2012). Mitochondrial evolution. Cold Spring Harb. Perspect. Biol.4, a011403–a011403. doi: 10.1101/cshperspect.a011403

  • CrossRef
  • Google Scholar

• 22

  GuoW.GreweF.FanW.YoungG. J.KnoopV.PalmerJ. D.et al. (2016). Ginkgo and Welwitschia mitogenomes reveal extreme contrasts in gymnosperm mitochondrial evolution. Mol. Biol. Evol.33, 1448–1460. doi: 10.1093/molbev/msw024

  • CrossRef
  • Google Scholar

• 23

  HanF.QuY.ChenY.XuL.BiC. (2022). Assembly and comparative analysis of the complete mitochondrial genome of Salix wilsonii using PacBio HiFi sequencing. Front. Plant Sci.13. doi: 10.3389/fpls.2022.1031769

  • CrossRef
  • Google Scholar

• 24

  HaoW.LiuG.WangW.ShenW.ZhaoY.SunJ.et al. (2021). RNA editing and its roles in plant organelles. Front. Genet.12. doi: 10.3389/fgene.2021.757109

  • CrossRef
  • Google Scholar

• 25

  JiaH.LiJ.ZhangJ.SunP.LuM.HuJ. (2019). The Salix psammophila SpRLCK1 involved in drought and salt tolerance. Plant Physiol. Biochem.144, 222–233. doi: 10.1016/j.plaphy.2019.09.042

  • CrossRef
  • Google Scholar

• 26

  JiaH.WangL.LiJ.SunP.LuM.HuJ. (2020). Physiological and metabolic responses of Salix sinopurpurea and Salix suchowensis to drought stress. Trees34, 563–577. doi: 10.1007/s00468-019-01937-z

  • CrossRef
  • Google Scholar

• 27

  JinJ.-J.YuW.-B.YangJ.-B.SongY.dePamphilisC. W.YiT.-S.et al. (2020). GetOrganelle: a fast and versatile toolkit for accurate de novo assembly of organelle genomes. Genome Biol.21, 241. doi: 10.1186/s13059-020-02154-5

  • CrossRef
  • Google Scholar

• 28

  JoY. D.ParkJ.KimJ.SongW.HurC.-G.LeeY.-H.et al. (2011). Complete sequencing and comparative analyses of the pepper (Capsicum annuum L.) plastome revealed high frequency of tandem repeats and large insertion/deletions on pepper plastome. Plant Cell Rep.30, 217–229. doi: 10.1007/s00299-010-0929-2

  • CrossRef
  • Google Scholar

• 29

  KatohK.StandleyD. M. (2013). MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Mol. Biol. Evol.30, 772–780. doi: 10.1093/molbev/mst010

  • CrossRef
  • Google Scholar

• 30

  KellerW.WolfJ.GerberA. (1999). Editing of messenger RNA precursors and of tRNAs by adenosine to inosine conversion. FEBS Lett.452, 71–76. doi: 10.1016/S0014-5793(99)00590-6

  • CrossRef
  • Google Scholar

• 31

  KerstenB.Faivre RampantP.MaderM.Le PaslierM.-C.BounonR.BerardA.et al. (2016). Genome sequences of Populus tremula chloroplast and mitochondrion: implications for holistic poplar breeding. PloS One11, e0147209. doi: 10.1371/journal.pone.0147209

  • CrossRef
  • Google Scholar

• 32

  KolmogorovM.YuanJ.LinY.PevznerP. A. (2019). Assembly of long, error-prone reads using repeat graphs. Nat. Biotechnol.37, 540–546. doi: 10.1038/s41587-019-0072-8

  • CrossRef
  • Google Scholar

• 33

  KumarS.StecherG.TamuraK. (2016). MEGA7: Molecular evolutionary genetics analysis version 7.0 for bigger datasets. Mol. Biol. Evol.33, 1870–1874. doi: 10.1093/molbev/msw054

  • CrossRef
  • Google Scholar

• 34

  KurtzS. (2001). REPuter: the manifold applications of repeat analysis on a genomic scale. Nucleic Acids Res.29, 4633–4642. doi: 10.1093/nar/29.22.4633

  • CrossRef
  • Google Scholar

• 35

  LetunicI.BorkP. (2019). Interactive Tree Of Life (iTOL) v4: recent updates and new developments. Nucleic Acids Res.47, W256–W259. doi: 10.1093/nar/gkz239

  • CrossRef
  • Google Scholar

• 36

  LewisS.SearleS.HarrisN.GibsonM.IyerV.RichterJ.et al. (2002). Apollo: a sequence annotation editor. Genome Biol.3, research0082.1. doi: 10.1186/gb-2002-3-12-research0082

  • CrossRef
  • Google Scholar

• 37

  LiZ.-J. (2009). Methods for the identification of horizontal gene transfer (HGT) events and progress in related fields: methods for the identification of horizontal gene transfer (HGT) events and progress in related fields. Hereditas (Beijing)30, 1108–1114. doi: 10.3724/SP.J.1005.2008.01108

  • CrossRef
  • Google Scholar

• 38

  LiC.LiuH.QinM.TanY.OuX.ChenX.et al. (2024). RNA editing events and expression profiles of mitochondrial protein-coding genes in the endemic and endangered medicinal plant, Corydalis saxicola. Front. Plant Sci.15. doi: 10.3389/fpls.2024.1332460

  • CrossRef
  • Google Scholar

• 39

  LiJ.YeC. (2020). Genome-wide analysis of microsatellite and sex-linked marker identification in Gleditsia sinensis. BMC Plant Biol.20, 338. doi: 10.1186/s12870-020-02551-9

  • CrossRef
  • Google Scholar

• 40

  LiX.ZhangY.HouM.SunF.ShenY.XiuZ.et al. (2014). Small kernel 1 encodes a pentatricopeptide repeat protein required for mitochondrial nad7 transcript editing and seed development in maize (Zea mays) and rice (Oryza sativa). Plant J.79, 797–809. doi: 10.1111/tpj.12584

  • CrossRef
  • Google Scholar

• 41

  LiuS.NiY.LiJ.ZhangX.YangH.ChenH.et al. (2023). CPGView: a package for visualizing detailed chloroplast genome structures. Mol. Ecol. Resour.23, 694–704. doi: 10.1111/1755-0998.13729

  • CrossRef
  • Google Scholar

• 42

  LoweT. M.EddyS. R. (1997). tRNAscan-SE: A program for improved detection of transfer RNA genes in genomic sequence. Nucleic Acids Res.25, 955–964. doi: 10.1093/nar/25.5.955

  • CrossRef
  • Google Scholar

• 43

  LuD.HaoL.HuangH.ZhangG. (2019). The complete chloroplast genome of Salix psamaphila, a desert shrub in northwest China. Mitochondrial DNA B Resour.4, 3432–3433. doi: 10.1080/23802359.2019.1675485

  • CrossRef
  • Google Scholar

• 44

  MaQ.LiS.BiC.HaoZ.SunC.YeN. (2017). Complete chloroplast genome sequence of a major economic species, Ziziphus jujuba (Rhamnaceae). Curr. Genet.63, 117–129. doi: 10.1007/s00294-016-0612-4

  • CrossRef
  • Google Scholar

• 45

  MaguireB. A.ZimmermannR. A. (2001). The ribosome in focus. Cell104, 813–816. doi: 10.1016/S0092-8674(01)00278-1

  • CrossRef
  • Google Scholar

• 46

  MaldonadoM.AbeK. M.LettsJ. A. (2022). A structural perspective on the RNA editing of plant respiratory complexes. Int. J. Mol. Sci.23, 684. doi: 10.3390/ijms23020684

  • CrossRef
  • Google Scholar

• 47

  MøllerI. M.RasmussonA. G.Van AkenO. (2021). Plant mitochondria – past, present and future. Plant J.108, 912–959. doi: 10.1111/tpj.15495

  • CrossRef
  • Google Scholar

• 48

  NaoraH.NaoraH. (1999). Involvement of ribosomal proteins in regulating cell growth and apoptosis: translational modulation or recruitment for extraribosomal activity? Immunol. Cell Biol.77, 197–205. doi: 10.1046/j.1440-1711.1999.00816.x

  • CrossRef
  • Google Scholar

• 49

  NarsaiR.WangC.ChenJ.WuJ.ShouH.WhelanJ. (2013). Antagonistic, overlapping and distinct responses to biotic stress in rice (Oryza sativa) and interactions with abiotic stress. BMC Genomics14, 93. doi: 10.1186/1471-2164-14-93

  • CrossRef
  • Google Scholar

• 50

  NguyenL.-T.SchmidtH. A.Von HaeselerA.MinhB. Q. (2015). IQ-TREE: a fast and effective stochastic algorithm for estimating maximum-likelihood phylogenies. Mol. Biol. Evol.32, 268–274. doi: 10.1093/molbev/msu300

  • CrossRef
  • Google Scholar

• 51

  PetersenG.AndersonB.BraunH.-P.MeyerE. H.MøllerI. M. (2020). Mitochondria in parasitic plants. Mitochondrion52, 173–182. doi: 10.1016/j.mito.2020.03.008

  • CrossRef
  • Google Scholar

• 52

  PetersenG.CuencaA.MøllerI. M.SebergO. (2015). Massive gene loss in mistletoe (Viscum, Viscaceae) mitochondria. Sci. Rep.5, 17588. doi: 10.1038/srep17588

  • CrossRef
  • Google Scholar

• 53

  PetibonC.Malik GhulamM.CatalaM.Abou ElelaS. (2020). Regulation of ribosomal protein genes: an ordered anarchy. Wiley Interdiscip. Rev. RNA12, e1632. doi: 10.1002/wrna.1632

  • CrossRef
  • Google Scholar

• 54

  QuY.ZhouP.TongC.BiC.XuL. (2022). Assembly and analysis of the Populus deltoides mitochondrial genome: the first report of a multicircular mitochondrial conformation for the genus Populus. J. For Res. (Harbin)34, 717–733. doi: 10.1007/s11676-022-01511-3

  • CrossRef
  • Google Scholar

• 55

  RamakrishnanV.WhiteS. W. (1998). Ribosomal protein structures: insights into the architecture, machinery and evolution of the ribosome. Trends Biochem. Sci.23, 208–212. doi: 10.1016/S0968-0004(98)01214-6

  • CrossRef
  • Google Scholar

• 56

  RichardsonA. O.PalmerJ. D. (2006). Horizontal gene transfer in plants. J. Exp. Bot.58, 1–9. doi: 10.1093/jxb/erl148

  • CrossRef
  • Google Scholar

• 57

  ShiL.ChenH.JiangM.WangL.WuX.HuangL.et al. (2019). CPGAVAS2, an integrated plastome sequence annotator and analyzer. Nucleic Acids Res.47, W65–W73. doi: 10.1093/nar/gkz345

  • CrossRef
  • Google Scholar

• 58

  SkippingtonE.BarkmanT. J.RiceD. W.PalmerJ. D. (2015). Miniaturized mitogenome of the parasitic plant Viscum scurruloideum is extremely divergent and dynamic and has lost all nad genes. Proc. Natl. Acad. Sci. U.S.A.112, E3515–E 3524. doi: 10.1073/pnas.1504491112

  • CrossRef
  • Google Scholar

• 59

  SloanD. B.AlversonA. J.ChuckalovcakJ. P.WuM.McCauleyD. E.PalmerJ. D.et al. (2012). Rapid evolution of enormous, multichromosomal genomes in flowering plant mitochondria with exceptionally high mutation rates. PloS Biol.10, e1001241. doi: 10.1371/journal.pbio.1001241

  • CrossRef
  • Google Scholar

• 60

  SmallI. D.Schallenberg-RüdingerM.TakenakaM.MireauH.Ostersetzer-BiranO. (2020). Plant organellar RNA editing: what 30 years of research has revealed. Plant J.101, 1040–1056. doi: 10.1111/tpj.14578

  • CrossRef
  • Google Scholar

• 61

  SmithD. R. (2011). Extending the limited transfer window hypothesis to inter-organelle DNA migration. Genome Biol. Evol.3, 743–748. doi: 10.1093/gbe/evr068

  • CrossRef
  • Google Scholar

• 62

  SunF.ZhangX.ShenY.WangH.LiuR.WangX.et al. (2018). The pentatricopeptide repeat protein EMPTY PERICARP 8 is required for the splicing of three mitochondrial introns and seed development in maize. Plant J.95, 919–932. doi: 10.1111/tpj.14030

  • CrossRef
  • Google Scholar

• 63

  TanakaY.TsudaM.YasumotoK.TerachiT.YamagishiH. (2014). The complete mitochondrial genome sequence of Brassica oleracea and analysis of coexisting mitotypes. Curr. Genet.60, 277–284. doi: 10.1007/s00294-014-0433-2

  • CrossRef
  • Google Scholar

• 64

  TillichM.LehwarkP.PellizzerT.Ulbricht-JonesE. S.FischerA.BockR.et al. (2017). GeSeq - versatile and accurate annotation of organelle genomes. Nucleic Acids Res.45, W6–W11. doi: 10.1093/nar/gkx391

  • CrossRef
  • Google Scholar

• 65

  VemannaR. S.BakadeR.BhartiP.KumarM. K. P.SreemanS. M.Senthil-KumarM.et al. (2019). Cross-talk signaling in rice during combined drought and bacterial blight stress. Front. Plant Sci.10. doi: 10.3389/fpls.2019.00193

  • CrossRef
  • Google Scholar

• 66

  WangZ.HouJ.LuL.QiZ.SunJ.GaoW.et al. (2013b). Small ribosomal protein subunit S7 suppresses ovarian tumorigenesis through regulation of the PI3K/AKT and MAPK pathways. PloS One8, e79117. doi: 10.1371/journal.pone.0079117

  • CrossRef
  • Google Scholar

• 67

  WangJ.LanP.GaoH.ZhengL.LiW.SchmidtW. (2013a). Expression changes of ribosomal proteins in phosphate- and iron-deficient Arabidopsis roots predict stress-specific alterations in ribosome composition. BMC Genomics14, 783. doi: 10.1186/1471-2164-14-783

  • CrossRef
  • Google Scholar

• 68

  WangY.WangH.YuC.YanX.ChuJ.JiangB.et al. (2024). Comprehensive bioinformation analysis of homeodomain-leucine zipper gene family and expression pattern of HD-Zip I under abiotic stress in Salix suchowensis. BMC Genomics25, 182. doi: 10.1186/s12864-024-10067-x

  • CrossRef
  • Google Scholar

• 69

  WeiS.WangX.BiC.XuY.WuD.YeN. (2016). Assembly and analysis of the complete Salix purpurea L. (Salicaceae) mitochondrial genome sequence. Springerplus5, 1894. doi: 10.1186/s40064-016-3521-6

  • CrossRef
  • Google Scholar

• 70

  WickR. R.SchultzM. B.ZobelJ.HoltK. E. (2015). Bandage: interactive visualization of de novo genome assemblies. Bioinformatics31, 3350–3352. doi: 10.1093/bioinformatics/btv383

  • CrossRef
  • Google Scholar

• 71

  WuZ.HuK.YanM.SongL.WenJ.MaC.et al. (2019). Mitochondrial genome and transcriptome analysis of five alloplasmic male-sterile lines in Brassica juncea. BMC Genomics20, 348. doi: 10.1186/s12864-019-5721-2

  • CrossRef
  • Google Scholar

• 72

  WuZ.LiaoX.ZhangX.TembrockL. R.BrozA. (2022). Genomic architectural variation of plant mitochondria-a review of multichromosomal structuring. J. Syst. Evol.60, 160–168. doi: 10.1111/jse.12655

  • CrossRef
  • Google Scholar

• 73

  WuH.YaoD.ChenY.YangW.ZhaoW.GaoH.et al. (2020). De Novo genome assembly of Populus simonii further supports that Populus simonii and Populus trichocarpa belong to different sections. G3 (Bethesda)10, 455–466. doi: 10.1534/g3.119.400913

  • CrossRef
  • Google Scholar

• 74

  WynnE. L.ChristensenA. C. (2019). Repeats of unusual size in plant mitochondrial genomes: identification, incidence and evolution. G3 (Bethesda)9, 549–559. doi: 10.1534/g3.118.200948

  • CrossRef
  • Google Scholar

• 75

  YangH.FanL.YuX.ZhangX.HaoP.WeiD.et al. (2022). Analysis of the NAC gene family in Salix and the identification of SpsNAC005 gene contributing to salt and drought tolerance. Forests13, 971. doi: 10.3390/f13070971

  • CrossRef
  • Google Scholar

• 76

  YangH.LiW.YuX.ZhangX.ZhangZ.LiuY.et al. (2021). Insights into molecular structure, genome evolution and phylogenetic implication through mitochondrial genome sequence of Gleditsia sinensis. Sci. Rep.11, 14850. doi: 10.1038/s41598-021-93480-6

  • CrossRef
  • Google Scholar

• 77

  YeN.WangX.LiJ.BiC.XuY.WuD.et al. (2017). Assembly and comparative analysis of complete mitochondrial genome sequence of an economic plant Salix suchowensis. PeerJ5, e3148. doi: 10.7717/peerj.3148

  • CrossRef
  • Google Scholar

• 78

  ZhangT.FangY.WangX.DengX.ZhangX.HuS.et al. (2012). The complete chloroplast and mitochondrial genome sequences of Boea hygrometrica: insights into the evolution of plant organellar genomes. PloS One7, e30531. doi: 10.1371/journal.pone.0030531

  • CrossRef
  • Google Scholar

• 79

  ZhangD.GaoF.JakovlićI.ZouH.ZhangJ.LiW. X.et al. (2020). PhyloSuite: an integrated and scalable desktop platform for streamlined molecular sequence data management and evolutionary phylogenetics studies. Mol. Ecol. Resour.20, 348–355. doi: 10.1111/1755-0998.13096

  • CrossRef
  • Google Scholar

• 80

  ZhangJ.LiuG.WeiJ. (2024). Assembly and comparative analysis of the first complete mitochondrial genome of Setaria italica. Planta260, 23. doi: 10.1007/s00425-024-04386-2

  • CrossRef
  • Google Scholar

• 81

  ZhangH.MeltzerP.DavisS. (2013). RCircos: an R package for Circos 2D track plots. BMC Bioinf.14, 244. doi: 10.1186/1471-2105-14-244

  • CrossRef
  • Google Scholar