Skip to main content

Comparative analysis of mitochondrial genomes reveals marine adaptation in seagrasses



Seagrasses are higher marine flowering plants that evolved from terrestrial plants, but returned to the sea during the early evolution of monocotyledons through several separate lineages. Thus, they become a good model for studying the adaptation of plants to the marine environment. Sequencing of the mitochondrial (mt) genome of seagrasses is essential for understanding their evolutionary characteristics.


In this study, we sequenced the mt genome of two endangered seagrasses (Zostera japonica and Phyllospadix iwatensis). These data and data on previously sequenced mt genomes from monocotyledons provide new evolutionary evidence of genome size reduction, gene loss, and adaptive evolution in seagrasses. The mt genomes of Z. japonica and P. iwatensis are circular. The sizes of the three seagrasses (including Zostera marine) that have been sequenced to date are smaller than that of other monocotyledons. Additionally, we found a large number of repeat sequences in seagrasses. The most abundant long repeat sequences were 31–40 bp repeats. Our study also found that seagrass species lost extensive ribosomal protein genes during evolution. The rps7 gene and the rpl16 gene of P. iwatensis are exceptions to this trend. The phylogenetic analysis based on the mt genome strongly supports the previous results. Furthermore, we identified five positive selection genes (atp8, nad3, nad6, ccmFn, and matR) in seagrasses that may be associated with their adaptation to the marine environment.


In this study, we sequenced and annotated the mt genomes of Z. japonica and P. iwatensis and compared them with the genome of other monocotyledons. The results of this study will enhance our understanding of seagrass adaptation to the marine environment and can inform further investigations of the seagrass mt genome.

Peer Review reports


Seagrasses are higher marine flowering plants whose ancestors are thought to have been terrestrial monocotyledons that re-adapted to the marine submerged environment through multiple evolutions 70–100 million years ago [1, 2]. Currently, 74 species of seagrasses are known worldwide [3, 4]. Seagrass beds provide habitats for many marine organisms, drive energy flow and material transport of the entire marine ecosystem, and play an important role in regulating global environmental change [1, 5, 6]. Seagrasses have the characteristics and unique structure of both terrestrial and marine plants and can complete their entire life history in the seawater environment, mainly in estuaries, intertidal wetlands of bays, and shallow subtidal waters [7]. Since the environment of seagrasses is highly saline, relatively anoxic, and frequently disturbed by waves and tides, seagrasses have evolved unique morphological structures unknown in terrestrial monocotyledons or freshwater submerged monocotyledons [8, 9]. The emergence of seagrasses represents one of the most striking evolutionary transformations in the evolutionary history of angiosperms. Their unique evolutionary history is one of the most interesting scientific questions in the field of evolutionary biology.

As a member of the “second entry into the sea” club (Like cetaceans in the ocean), seagrass and its evolutionary biology have been a hot topic in recent years. Studies have attempted to reveal the adaptive evolutionary mechanisms of seagrass to the marine environment at the molecular and genomic levels. In 2016, Olsen et al. studied the genome of Z. marina for the first time and found that Z. marina has lost some genes commonly found in terrestrial angiosperms, such as genes regulating stomatal opening, infrared sensing and ultraviolet protection genes, and terpene and ethylene synthesis genes, and gained some genes regulating ion metabolism, gas exchange, nutrient uptake, and other life activities for adaptation to specific marine environments [10]. In the same year, Lee et al. also found that Zostera muelleri had lost or modified multiple genes (such as hormone biosynthesis and signalling genes, and cell wall catabolism genes) compared to terrestrial or floating aquatic plants in order to adapt to the sea [11]. Subsequently, Lee et al. explored the convergent evolutionary features of Halophila ovalis, Z. marina, and Z. muelleri to the marine environment. They found that all three seagrasses lost genes related to ethylene and terpenoid biosynthesis while retaining genes related to salt tolerance [12]. Our previous study showed that the coding sequences of the chloroplast genome of seagrasses have a higher evolutionary rate compared to terrestrial monocotyledons, which may be related to their habitat adaptation [13].

Mitochondria are unique organelles of eukaryotic cells with their own independent genetic system. Because they provide energy for cellular life activities, they are known as the “power house”. The plant mt genome encodes genes related to oxidative phosphorylation and respiratory metabolism, including cytochrome complex I-V subunit genes, cytochrome C synthesis-related genes, ribosomal protein-related genes, tRNAs, rRNAs, and some open reading frames of unknown function [14]. The functional gene sequences of plant mt genomes are very conservative, but the number, location, and order of functional genes vary among species [15]. Adaptive evolution is defined as the improved adaptability of a species to changing environmental conditions during the evolutionary process. Adaptive evolution refers to the ability of a species to adapt well to changing environmental conditions during the evolutionary process, and this process should occur as a distinct change in the sequence of this species. Given the conservation of functional genes in the mt genome, the study of these functional genes will enhance our understanding of the evolution of plant adaptations. To the best of our knowledge, no comparative analysis of the mt genomes of seagrasses adapted to the marine environment has been performed. In this study, we sequenced and annotated the mt genomes of Z. japonica and P. iwatensis and compared them with the genome of other monocotyledons (which contains all currently published species of Alismatales). Our results will provide additional information for a better understanding of seagrass genetics.


Genomic features of Zosteraceae species

The mt genomes of Z. japonica and P. iwatensis were circular, while the Z. marine was linear. Their lengths were 221,614 (Z. japonica), 178,929 (P. iwatensis), and 191,481 (Z. marine) (Fig. 1 and Table 1). The highest GC content was in P. iwatensis (48.69%), followed by Z. japonica (46.26%), and Z. marina (45.14%). There are 44, 50, and 38 genes annotated in the P. iwatensis, Z. japonica, and Z. marina mt genomes, respectively, including 26, 29, and 25 PCGs and 15, 18, and 10 tRNA genes, respectively; and three rRNA genes. The number of PCGs present in the three mt genomes of the Zosteraceae species was mostly single-copy PCGs, except for Z. japonica, which had three double-copy PCGs (matR, nad4 and nad4L) and one triple-copy PCG (atp6). Comparing the PCGs of the three seagrasses with three other Alismatales species and three terrestrial monocotyledons, we found that seagrass species have lost many rpl genes (except for rpl16 gene in P. iwatensis) and rps genes (except for the rps7 gene) during the evolutionary process (Fig. 2). Analysis of codon usage in the three species showed a strong preference for the ten codon families, Ala, Arg, Gly, Ile, Leu1, Pro, Ser2, Termination codons (TER), Thr, and Val (Fig. 3).

Fig. 1
figure 1

Gene maps of Z. japonica and P. iwatensis mt genomes

Table 1 Basic information on the mt genomes of the nine species in this study
Fig. 2
figure 2

Protein-coding genes present in the mt genomes of the nine species. Dark boxes indicate the presence of three genes, dark grey boxes indicate the presence of two genes, light grey boxes indicate the presence of one gene, and white boxes indicate the absence of the gene. Taxa in green show seagrass, and taxa in yellow show non-marine monocotyledons

Fig. 3
figure 3

The relative synonymous codon usage (RSCU) of three seagrasses. Each bar chart from left to right is Z. marina, Z. japonica, and P. iwatensis. Codon families are plotted on the X axis

Repeat sequences analysis

A total of 65 (Z. japonica and Z. marina) to 77 (P. iwatensis) SSRs were detected in the three mt genomes, including 38–49 mononucleotides (mono-), 6–9 dinucleotides (di-), 2–7 trinucleotides (tri-), 10–13 tetranucleotides (tetra-), 0–2 pentanucleotide (penta-), and 0–1 hexanucleotide (hexa-) (Fig. 4A). The highest percentage of SSRs was mononucleotides (63.28%), followed by tetranucleotides (16.91%), dinucleotides (11.11%), tetranucleotides (6.28%), pentanucleotides (1.93%), and hexanucleotides (0.48%) (Fig. 4B). Furthermore, we found that the A/T was the most common mononucleotide SSRs type in the three species, and that the repeat units of the other five SSRs were composed mainly of A or T (Fig. 4C). Tandem (T), forward (F), reverse (R), palindromic (P), and complement (C) repeat sequences in the three seagrass mt genomes were identified using a tandem repeat finder and REPuter. We identified three long repeat sequence types in the seagrass mt genome: tandem repeats, forward repeats, and palindromic repeats. We identified 3028 (P. iwatensis) to 5101 (Z. marina) long repeat elements, of which 19–91 were tandem repeats (Fig. 4D). The lengths of these long repeat sequences were highly variable, with 31–40 bp repeats being the most common among the three repeat types (Fig. 4E).

Fig. 4
figure 4

Analysis of repeated sequences of the three seagrass mt genomes compared in this study. A The number of SSRs in six SSR types; B distribution of SSR types; C the number of different SSR repeat unit types; D: the number of repeat elements in three long repeat types; and E: summary of long repeat types by length

Phylogenetic relationships

To understand the evolutionary status of seagrass mt genomes, a phylogenetic analysis was conducted on Z. japonica, P. iwatensis, and seven other monocots species, including four Alismatales (Z. marina, S. polyrhiza, S. aloides, and B. umbellatus), one Asparagales (A. officinalis), one Arecales (P. dactylifera), and one Poales (E. indica). A phylogenetic tree was obtained based on the alignment data matrix of 23 shared PCGs (Fig. 5). The phylogenetic tree shows a clear division of the Alismatales and other monocotyledons into two taxa. Regarding Alismatales, the phylogenetic tree shows that Z. japonica was a sister species to Z. marina, followed by P. iwatensis. These results strongly support the theory that marine seagrasses species evolved from the other freshwater species. These data are consistent with the results of a phylogenetic tree based on chloroplast genomes [13].

Fig. 5
figure 5

Phylogenetic tree constructed using the maximum likelihood (ML) methods, based on the 23 concatenate protein-coding sequences of whole mt genomes from nine monocotyledons. The green, blue, and purple colors of the arcs between the tree and the picture represent marine, freshwater, and terrestrial species, respectively

Adaptive evolution analysis

Analysis of the branch model for all shared PCGs showed that the ω values (ω = 0.41962) was significant (p <  0.01), indicating that monocotyledonous plants were strongly purified and selected for during their evolution at the whole-mitochondrion protein level (Table 2). Moreover, although the ω values of the seagrass group were lower than those of the terrestrial monocotyledonous group, the evolutionary rate of mt protein-encoding genes in seagrasses did not differ significantly from that of terrestrial monocotyledonous plants (p = 0.28). The branch-site model was used to examine the positive selection of individual codons in the seagrass lineage. We found that five genes (atp8, ccmFn, nad3, nad6, and matR) were positively selected in seagrasses (p-value < 0.05) and the BEB posterior probability for 20 amino acid sites larger than 0.90. Among these positively selected genes, the matR gene had the highest number of positively selected sites. To better understand the distribution of positive selection in seagrasses, we selected the PCGs of Z. japonica as the reference sequence and mapped the above inferred positive selection sites visually to the results of sequence alignment (Fig. 6).

Table 2 Selection pressure test results of all shared protein-coding genes in monocotyledons through a Branch model analysis
Fig. 6
figure 6

Selection pressure test results of shared protein-coding genes in seagrasses through a branch-site model analysis. The partial alignment of five genes suggesting sites with positive selection in the BEB test. The red blocks stand for the amino acids in seagrass with BEB posterior probability > 0.90. All positively selected sites used the gene location information of Z. japonica as the reference sequence


Mitogenome size evolution

Mitochondria are organelles with semi-autonomous genetic properties that perform respiratory metabolism to provide energy for life activities. Plant mitochondria possess more complex genomes than animals, with extensive size variations, sequence alignments, repeat content, and highly conserved coding sequences [16]. The circular mt genomes of the two seagrasses (Z. japonica and P. iwatensis) were sequenced for this study, and the previously published linear mt genome of Z. marina was roughly similar in size and gene composition. Their lengths were 221,614 (Z. japonica), 178,929 (P. iwatensis), and 191,481 (Z. marina). The size of the mt genomes of these three species was smaller than that of other monocotyledons. P. iwatensis has the smallest mt genome of all monocotyledons sequenced so far. Previous studies have shown that the relationship between the size of plant mt genomes and their gene numbers exhibits a C-value paradox, i.e., there is no positive correlation between increasing mt genome size and increasing gene numbers [17]. In the present study, we also found gene loss in seagrass species, but the reduced length of the lost genes was not the underlying cause of the reduced mt genome size. It has also been shown that the size of the mt genome is positively influenced by repeat sequences, which could explain why the mt genome of B. umbellatus, which had high repeat sequences (approximately 8.3%), was larger than the mt genome of S. aloides, which had low repeat sequences (approximately 0.6%) [18, 19]. Surprisingly, the three seagrasses with the smallest mt genomes had the largest number of repeat sequences, so there did not seem to be a clear correlation between the size of the mt genome and the number of repeat sequences. In general, the Alismatales mt genome is relatively small compared with other monocotyledons. Our study determined that the size of mt genomes of monocotyledons in different habitats differed significantly (Kruskal–Wallis, p = 0.027), with the mt genomes of marine seagrasses in Alismatales being smaller than those of other freshwater species and significantly smaller than those of other terrestrial monocotyledons (Nemeyi, p = 0.0199) (Fig. 7). These results suggest that different habitat types may have an important effect on the mt genome size of monocotyledons. We speculate that during the evolutionary process of “secondary entry into the sea,” seagrasses gradually discarded some “useless sequences” to streamline their mt genomes to better adapt to the marine environment. The GC content of the mt genomes of the three seagrasses ranged from 45.14 to 48.69%, similar to the range of mt genomes in other monocotyledons (Table 1). These results support the conclusion that GC content is conserved in higher plants.

Fig. 7
figure 7

Comparison of mt genome sizes in nine monocotyledonous species. A represents marine species including P. iwatensis, Z. japonica, and Z. marina; B represents freshwater species including S. polyrhiza, S. aloides, and B. umbellatus; and C represents terrestrial species including A. officinalis, P. dactylifera, and E. indica

Mitochondrial protein-coding genes in Seagrasses

Although the functional gene sequences of plant mt genomes are strongly conserved, the number, location, and order of functional genes vary among species [15]. In this study, we found that with the exception of the rps7 gene and the rpl16 gene of P. iwatensis, seagrass species lost extensive ribosomal protein genes during evolution. A previous survey of mt ribosomal protein genes in 280 angiosperms showed that Alismatales experienced significant gene loss [20]. Our study confirmed this and further found that the loss of ribosomal protein genes was more extensive in seagrasses. The above results seem to imply that ribosomal protein genes in seagrass mt genomes are dispensable. Studies have shown that there is DNA transfer between the organelle genome and the nuclear genome in plants. Where there is DNA sequence transfer between the nuclear genome and the mt genome, chloroplast DNA can be transferred to the mt genome [21, 22]. Whether the ribosomal protein genes missing from the seagrass mt genomes were completely lost or transferred to the nuclear genome was unknown until the complete nuclear genome sequence of Z. marina was published, providing evidence that at least four ribosomal protein-coding genes had been transferred to the nucleus [10, 19]. Interestingly, the rps7 gene is a ribosomal protein gene that is retained by all seagrass species. The reason for this phenomenon may be because the functional transfer to the nucleus is rare [23, 24]. ATP synthase in the plant mt endomembrane, which consists of two parts (F0 and F1), plays an important role in energy formation [25]. The atp9 gene codes for a structural protein component of complex V, which is part of the membrane-embedded F0 domain. The atp9 protein, along with atp8 and atp6, contributes to a ring-like structure, forming the channel through which protons diffuse down the mt membrane gradient. This allows ion transport across the membrane against the electrochemical potential gradient, which in turn drives the F1 part of the rotation motor to synthesize the energy molecule ATP [26,27,28]. Changes in the DNA sequence within the atp6 and atp9 genes and their upstream or downstream regions are known to be closely related to male sterility in plants [29, 30]. In our study, we found that the atp9 gene was lost in Z. japonica, whereas three copies of the atp6 gene existed. We speculate that the ATP6 gene may have replaced the ATP9 gene in the mt genome of Z. japonica and is performing some of its functions.

Adaptive evolution

Previous studies have found that the mt genes in seed plants evolved more slowly than their chloroplast genes and nuclear genes [31]. In this study, we found ω values less than 1 in all monocotyledons, indicating that the main force affecting the evolution of mt genes was purifying selection. Furthermore, our study found no significant differences in the evolutionary rates of marine seagrasses and other monocotyledons, suggesting that mt PCGs may be conserved among different monocotyledons (p > 0.05). Seagrasses grow primarily in the transition zone between marine and terrestrial environments, and as a result, seagrass habitats experience large fluctuations in temperature, salinity, pH, and oxygen concentration. Tidal environments are periodically exposed to the air and susceptible to enhanced UV-B radiation. Previous studies have shown that environmental energy stimulates metabolism at many levels and that energy-rich habitats tend to have higher evolutionary rates [32, 33]. Solar radiation (especially UV radiation) plays a direct role in mutagenesis and may accelerate molecular evolution [34, 35]. Although no significant differences were seen in the evolutionary rates of marine seagrasses and other monocotyledons, the branch-site model analysis showed that 30 sites in the seagrasses lineage on atp8 (2/30), ccmFn (3/30), nad3 (3/30), nad6 (2/30), and matR (20/30) were under positive selection pressure. This suggested that there were inconsistent selection pressures across sites. Positive selection may provide important functional details related to adaptation to new environments [36]. Previous studies have shown that NADH dehydrogenase is the first as well as the largest complex enzyme in the respiratory chain, while ATP synthase is the last, directly producing ATP [37, 38]. Among them, NADH dehydrogenase acts as a proton pump, so mutations in this complex may affect the efficiency of the proton pump and should affect metabolic efficiency [39, 40]. ATP synthase couples with the electrochemical gradient of the inner mt membrane to synthesize ATP. Variation in these gene sequences may affect the production of ATP [37, 41]. Previous studies have also shown that adaptive evolution of the NADH dehydrogenase complex and ATP synthase play very important roles in the high altitude and intertidal environment adaptation of animals [42,43,44]. Likewise, seagrasses live in an oxygen-deficient marine environment, which implies an association between gene positive selection and adaptation to energy metabolism in marine environments. In this study, we propose that the positively selected genes nad3, nad6, and ATP8 play a role in adaptation to the marine environment in seagrasses. The two remaining positively selected genes were ccmFn and matR. The ccmFn gene primarily synthesizes the N terminus of the cytochrome C maturation protein subunit F and has been reported to increase transcript levels in wheat embryos and seedlings subjected to salt and osmotic stresses [45]. MatR is an mt maturase-related gene with the most positive selection sites and has been retained as a conserved ORF in the mtDNAs in nearly all angiosperms [20]. We speculate that these two positively selected genes may also play an important role in enabling seagrasses to adapt to the marine environment.


The mt genome of two seagrasses, Z. japonica and P. iwatensis, were assembled and annotated in this study, and extensive analyses were performed based on the DNA sequences and amino acid sequences of the annotated genes. The mt genomes of Z. japonica and P. iwatensis are circular, while the mt genome of Z. marine is linear. The mt genomes of the three seagrasses are smaller than those of other monocotyledons. Repeated sequence analysis showed that the seagrass mt genome contained a large number of repeated sequences. The 31–40 bp repeats were the most abundant long repeat sequences. We also found that seagrass species have lost extensive ribosomal protein genes during evolution, except for the rps7 gene and the rpl16 gene of P. iwatensis. Phylogenetic results showed that Z. japonica was a sister species to Z. marina, followed by P. iwatensis within Zosteraceae. Finally, we identified five positive selection genes in seagrasses that may be related to their adaptation to the marine environment. These findings will not only help to further investigate the seagrass mt genome, but also be valuable for future studies on the adaptation of seagrasses to the marine environment.

Materials and methods

Sample collection and sequencing assembly

The seagrasses Z. japonica (120.6801 E, 37.93384 N) and P. iwatensis (120.6476 E, 37.97221 N) were collected from collected from Yantai, China, in 2020 and identified by Prof. Xuexi Tang. The Z. japonica (voucher number: OUC-S118) and P. iwatensis (voucher number: OUC-S119) are currently stored at the Marine Ecology Laboratory of Ocean University of China. The data for Zostera marina were downloaded from NCBI ( Total DNA was extracted from fresh leaves using TRIzol® Reagent (Invitrogen). Total DNA was sequenced using the PacBio Sequel II and an Illumina NovaSeq 6000 platform (150 bp*2) (Shanghai BIOZERON Co., Ltd., Shanghai, China) after library construction. Illumina sequences were assembled using GetOrganelle v1.7.1 [46]. They were compared to PacBio sequences using BWA v0.7.17 [47] in order to extract the trigeneration data of the target samples. The extracted PacBio sequences were mixed with the Illumina sequences to assemble the mt genome using SPAdes v3.14.1 [48]. The assembled mt genomes were polished with clean reads to correct assembled errors using Pilon v1.18 [49]. Zostera marina was used as a reference to determine its scaffold start position and orientation to obtain the final mt genome sequences.

Genome annotation

We compared the protein sequence of Zostera marina to the sample genome sequence, to filter out the bad comparison results and remove redundancy. We then ran GeneWise v2.1 [50] for an accurate comparison. Next, we used AUGUSTUS v3.2.2 software [51] to perform De novo gene prediction on the mt genome. Finally, we corrected and integrated the gene set to obtain the sample genome coding genes. The noncoding RNAs (ncRNAs) contained in the genome were predicted using RNAmmer and tRNAscan-SE software [52, 53]. A circular map of the two seagrass species was obtained using Organellar Genome DRAW [54]. The relative synonymous codon usage (RSCU) for all mtDNA of three seagrasses was estimated based on protein-coding genes using CodonW (

Repeat sequences analysis

Simple sequence repeat (SSR) markers were identified in the three seagrasses sequences using MISA [55] ( The repeats of 1, 2, 3, 4, 5, and 6 bases with 8, 5, 4, 3, 3, and 3 repeat numbers, respectively, were identified in this analysis. Forward (F), reverse (R), palindromic (P), and complement (C) repeats were searched using the online program REPuter [56] ( with Hamming distance set to 3 and a minimum repeat size of 30 bp. The tandem repeats with > 6 bp repeat unit were detected using Tandem Repeats Finder v4.09 software [57] ( with default parameters.

Phylogenetic construction

In order to determine the phylogenetic position of seagrasses in the Alismatales, we analyzed a total of nine species in this study. Among them, the mt genome sequences of two species (Z. japonica and P. iwatensis) were newly sequenced, and those of the other Monocotyledons species were downloaded from NCBI database. These sampled Monocotyledons species included all four current species of Alismatales (Z. marina, NC_035345; Spirodela polyrhiza, NC_017840; Butomus umbellatus, NC_021399; and Stratiotes aloides, NC_035317), one representative species of Poales (Eleusine indica, NC_040989), one representative species of Arecales (Phoenix dactylifera, NC_016740), and one representative species of Asparagales (Asparagus officinalis, NC_053642). We counted all protein-coding gene sequences (PCGs) of the nine mt genomes. All shared single-copy PCGs were concatenated into a super-matrix and aligned using MUSCLE [58]. The phylogenetic tree was constructed by using the maximum likelihood (ML) method based on the GTRGAMMA model with 1000 bootstrap replicates using RAxML [59]. Finally, iTOL ( was used to visualize and annotate the tree.

Adaptive evolution analysis

The non-synonymous (dN)/synonymous (dS) substitution rate (ω = dN/dS) between single-copy PCGs are usually used to determine the positive Darwinian selection pressure acting on genes. Ratios of ω > 1, = 1 and < 1 represent positive selection, neutral selection and negative selection, respectively. The CODEML program implemented in PAML 4.7 package was used to determine positive selection sites or branches [60]. A likelihood ratio test (LRT) with χ2 distribution was performed to determine and compare nested models with a p <  0.05 significance threshold. Bayes empirical Bayes (BEB) analysis was used to identify sites under positive selection with Bayesian posterior probability > 0.90 [61].

The branches with seagrasses were used as foreground branches, and the other monocotyledons were the background branches. To determine whether the shared PCGs of seagrasses faced different evolutionary forces in different habitats, we used a null one-ratio model (model = 0, NSsites = 0). This model estimated the same ω for all branches against the two-ratio model (model = 2, NSsites = 0), which estimated a variable ω in a specific branch using LRT. A P-value < 0.05 was selected to reject the null one-ratio model and evaluate the significance of the alternative hypothesis. We also used a branch-site model in PAML to detect positive selection genes in seagrasses. The LRT was used to contrast a null model Ma0 (model = 2, NSsites = 2, fix-omega = 1, omega = 1) and a model Ma (model = 2, NSsites = 2, fix-omega = 0, omega = 1) of positive selection pressures on the foreground branch. These analyses were based on the ML tree in our study (Fig. 5).

Availability of data and materials

The datasets generated for this study can be found in National Center for Biotechnology Information (NCBI) under the accession numbers: OP441720 ( and OP441721 (





Relative synonymous codon usage


Termination codons


Noncoding RNAs


Simple sequence repeat


Protein-coding gene sequences


Maximum likelihood

dN :


dS :



Likelihood ratio test


Bayes empirical Bayes


  1. Orth RJ, Carruthers TJ, Dennison WC, Duarte CM, Fourqurean JW, Heck KL, et al. A global crisis for seagrass ecosystems. Bioscience. 2006;56(12):987–96.

    Article  Google Scholar 

  2. Waycott M, Procaccini G, Les DH, Reusch T. Seagrass evolution, ecology and conservation: a genetic perspective. Netherlands: Springer; 2007.

    Google Scholar 

  3. Short FT, Polidoro B, Livingstone SR, Carpenter KE, Bandeira S, Bujang JS, et al. Extinction risk assessment of the world’s seagrass species. Biol Conserv. 2011;144(7):1961–71.

    Article  Google Scholar 

  4. Huang X-P, Jiang Z-J, Fan H-Q, Chen X-Y, Zhou Y, Zhang J-P, et al. The nomenclature of the “algae” name of seagrasses in China. Oceanol Limnol Sin. 2016;47(1):290–4.

    Google Scholar 

  5. Duarte CM, Middelburg JJ, Caraco N. Major role of marine vegetation on the oceanic carbon cycle. Biogeosciences. 2005;2(1):1–8.

    Article  CAS  Google Scholar 

  6. Hemminga MA, Duarte CM. Seagrass ecology. Cambridge: Cambridge University Press; 2000.

    Book  Google Scholar 

  7. Phillips RC, Menez EG. Seagrasses. Smithson Contrib Mar Sci. Vol. 34. Smithsonian institution press; 1988.

  8. Kuo J, Cambridge M, Kirkman H. Anatomy and structure of Australian seagrasses. In: Seagrasses of Australia. Cham: Springer; 2018. p. 93–125.

    Chapter  Google Scholar 

  9. Kuo J, Den Hartog C. Seagrass morphology, anatomy, and ultrastructure. In: Seagrasses: biology, ecologyand conservation. Dordrecht: Springer; 2007. p. 51–87.

    Google Scholar 

  10. Olsen JL, Rouzé P, Verhelst B, Lin Y-C, Bayer T, Collen J, et al. The genome of the seagrass Zostera marina reveals angiosperm adaptation to the sea. Nature. 2016;530(7590):331–5.

    Article  CAS  PubMed  Google Scholar 

  11. Lee H, Golicz AA, Bayer PE, Jiao Y, Tang H, Paterson AH, et al. The genome of a southern hemisphere seagrass species (Zostera muelleri). Plant Physiol. 2016;172(1):272–83.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Lee H, Golicz AA, Bayer PE, Severn-Ellis AA, Chan C-KK, Batley J, et al. Genomic comparison of two independent seagrass lineages reveals habitat-driven convergent evolution. J Exp Bot. 2018;69(15):3689–702.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. Chen J, Zang Y, Shang S, Liang S, Zhu M, Wang Y, et al. Comparative chloroplast genomes of Zosteraceae species provide adaptive evolution insights into Seagrass. Front Plant Sci. 2021;12:741152.

    Article  PubMed  PubMed Central  Google Scholar 

  14. Millar AH, Whelan J, Soole KL, Day DA. Organization and regulation of mitochondrial respiration in plants. Annu Rev Plant Biol. 2011;62:79–104.

    Article  CAS  PubMed  Google Scholar 

  15. Palmer JD, Herbon LA. Plant mitochondrial DNA evolved rapidly in structure, but slowly in sequence. J Mol Evol. 1988;28(1):87–97.

    Article  CAS  PubMed  Google Scholar 

  16. Kozik A, Rowan BA, Lavelle D, Berke L, Schranz ME, Michelmore RW, et al. The alternative reality of plant mitochondrial DNA: one ring does not rule them all. PLoS Genet. 2019;15(8):e1008373.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Palmer JD, Adams KL, Cho Y, Parkinson CL, Qiu Y-L, Song K. Dynamic evolution of plant mitochondrial genomes: mobile genes and introns and highly variable mutation rates. Proc Natl Acad Sci. 2000;97(13):6960–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Allen JO, Fauron CM, Minx P, Roark L, Oddiraju S, Lin GN, et al. Comparisons among two fertile and three male-sterile mitochondrial genomes of maize. Genetics. 2007;177(2):1173–92.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Petersen G, Cuenca A, Zervas A, Ross GT, Graham SW, Barrett CF, et al. Mitochondrial genome evolution in Alismatales: size reduction and extensive loss of ribosomal protein genes. PLoS One. 2017;12(5):e0177606.

    Article  PubMed  PubMed Central  Google Scholar 

  20. Adams KL, Qiu Y-L, Stoutemyer M, Palmer JD. Punctuated evolution of mitochondrial gene content: high and variable rates of mitochondrial gene loss and transfer to the nucleus during angiosperm evolution. Proc Natl Acad Sci. 2002;99(15):9905–12.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Tsunewaki K. Interorganellar DNA transfer in wheat: dynamics and phylogenetic origin. Proc Jpn Acad Ser B Phys Biol Sci. 2011;87(8):529–49.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Hao W, Palmer JD. Fine-scale mergers of chloroplast and mitochondrial genes create functional, transcompartmentally chimeric mitochondrial genes. Proc Natl Acad Sci. 2009;106(39):16728–33.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Park S, Grewe F, Zhu A, Ruhlman TA, Sabir J, Mower JP, et al. Dynamic evolution of Geranium mitochondrial genomes through multiple horizontal and intracellular gene transfers. New Phytol. 2015;208(2):570–83.

    Article  CAS  PubMed  Google Scholar 

  24. Liu S-L, Zhuang Y, Zhang P, Adams KL. Comparative analysis of structural diversity and sequence evolution in plant mitochondrial genes transferred to the nucleus. Mol Biol Evol. 2009;26(4):875–91.

    Article  PubMed  Google Scholar 

  25. Boyer PD. The binding change mechanism for ATP synthase—some probabilities and possibilities. Biochim Biophys Acta. 1993;1140(3):215–50.

    Article  CAS  PubMed  Google Scholar 

  26. Bietenhader M, Martos A, Tetaud E, Aiyar RS, Sellem CH, Kucharczyk R, et al. Experimental relocation of the mitochondrial ATP9 gene to the nucleus reveals forces underlying mitochondrial genome evolution. PLoS Genet. 2012;1817(8):S6.

    Google Scholar 

  27. Jonckheere AI, Smeitink JA, Rodenburg RJ. Mitochondrial ATP synthase: architecture, function and pathology. J Inherit Metab Dis. 2012;35(2):211–25.

    Article  CAS  PubMed  Google Scholar 

  28. Rak M, Zeng X, Brière J-J, Tzagoloff A. Assembly of F0 in Saccharomyces cerevisiae. Biochim Biophys Acta. 2009;1793(1):108–16.

    Article  CAS  PubMed  Google Scholar 

  29. Kim DH, Kim B-D. The organization of mitochondrial atp6 gene region in male fertile and CMS lines of pepper (Capsicum annuum L.). Curr Genet. 2006;49(1):59–67.

    Article  CAS  PubMed  Google Scholar 

  30. Dieterich J-H, Braun H-P, Schmitz U. Alloplasmic male sterility in Brassica napus (CMS'Tournefortii-Stiewe') is associated with a special gene arrangement around a novel atp9 gene. Mol Gen Genomics. 2003;269(6):723–31.

    Article  CAS  Google Scholar 

  31. Drouin G, Daoud H, Xia J. Relative rates of synonymous substitutions in the mitochondrial, chloroplast and nuclear genomes of seed plants. Mol Phylogenet Evol. 2008;49(3):827–31.

    Article  CAS  PubMed  Google Scholar 

  32. Davies TJ, Savolainen V, Chase MW, Moat J, Barraclough TG. Environmental energy and evolutionary rates in flowering plants. Proc R Soc Lond Ser B Biol Sci. 2004;271(1553):2195–200.

    Article  Google Scholar 

  33. Clarke A, Gaston KJ. Climate, energy and diversity. Proc R Soc B Biol Sci. 2006;273(1599):2257–66.

    Article  Google Scholar 

  34. Rothschild LJ. The influence of UV radiation on protistan evolution. J Eukaryot Microbiol. 1999;46(5):548–55.

    Article  CAS  PubMed  Google Scholar 

  35. Willis K, Bennett K, Birks H. Variability in thermal and UV-B energy fluxes through time and their influence on plant diversity and speciation. J Biogeogr. 2009;36(9):1630–44.

    Article  Google Scholar 

  36. Nielsen R. Molecular signatures of natural selection. Annu Rev Genet. 2005;39:197–218.

    Article  CAS  PubMed  Google Scholar 

  37. Mishmar D, Ruiz-Pesini E, Golik P, Macaulay V, Clark AG, Hosseini S, et al. Natural selection shaped regional mtDNA variation in humans. Proc Natl Acad Sci. 2003;100(1):171–6.

    Article  CAS  PubMed  Google Scholar 

  38. Weiss H, Friedrich T, Hofhaus G, Preis D. The respiratory-chain NADH dehydrogenase (complex I) of mitochondria. Eur J Biochem. 1991;1991:55–68.

    Google Scholar 

  39. Da Fonseca RR, Johnson WE, O'Brien SJ, Ramos MJ, Antunes A. The adaptive evolution of the mammalian mitochondrial genome. BMC Genomics. 2008;9(1):1–22.

    Article  Google Scholar 

  40. Hassanin A, Ropiquet A, Couloux A, Cruaud C. Evolution of the mitochondrial genome in mammals living at high altitude: new insights from a study of the tribe Caprini (Bovidae, Antilopinae). J Mol Evol. 2009;68(4):293–310.

    Article  CAS  PubMed  Google Scholar 

  41. Wallace DC. Why do we still have a maternally inherited mitochondrial DNA? Insights from evolutionary medicine. Annu Rev Biochem. 2007;76:781–821.

    Article  CAS  PubMed  Google Scholar 

  42. Zhou T, Shen X, Irwin DM, Shen Y, Zhang Y. Mitogenomic analyses propose positive selection in mitochondrial genes for high-altitude adaptation in galliform birds. Mitochondrion. 2014;18:70–5.

    Article  CAS  PubMed  Google Scholar 

  43. Li F, Lv Y, Wen Z, Bian C, Zhang X, Guo S, et al. The complete mitochondrial genome of the intertidal spider (Desis jiaxiangi) provides novel insights into the adaptive evolution of the mitogenome and the evolution of spiders. BMC Ecol Evol. 2021;21(1):1–15.

    Google Scholar 

  44. Wang J, Shi Y, Elzo MA, Dang S, Jia X, Lai S. Genetic diversity of ATP8 and ATP6 genes is associated with high-altitude adaptation in yak. Mitochondrial DNA A DNA Mapp Seq Anal. 2018;29(3):385–93.

    CAS  PubMed  Google Scholar 

  45. Naydenov NG, Khanam S, Siniauskaya M, Nakamura C. Profiling of mitochondrial transcriptome in germinating wheat embryos and seedlings subjected to cold, salinity and osmotic stresses. Genes Genet Syst. 2010;85(1):31–42.

    Article  CAS  PubMed  Google Scholar 

  46. Jin J-J, Yu W-B, Yang J-B, Song Y, DePamphilis CW, Yi T-S, et al. GetOrganelle: a fast and versatile toolkit for accurate de novo assembly of organelle genomes. Genome Biol. 2020;21(1):1–31.

    Article  Google Scholar 

  47. Li H, Durbin R. Fast and accurate short read alignment with burrows-wheeler transform. Bioinformatics. 2009;25(14):1754–60.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  48. Bankevich A, Nurk S, Antipov D, Gurevich AA, Dvorkin M, Kulikov AS, et al. SPAdes: a new genome assembly algorithm and its applications to single-cell sequencing. J Comput Biol. 2012;19(5):455–77.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  49. Walker BJ, Abeel T, Shea T, Priest M, Abouelliel A, Sakthikumar S, et al. Pilon: an integrated tool for comprehensive microbial variant detection and genome assembly improvement. PLoS One. 2014;9(11):e112963.

    Article  PubMed  PubMed Central  Google Scholar 

  50. Birney E, Clamp M, Durbin R. GeneWise and genomewise. Genome Res. 2004;14(5):988–95.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  51. Stanke M, Diekhans M, Baertsch R, Haussler D. Using native and syntenically mapped cDNA alignments to improve de novo gene finding. Bioinformatics. 2008;24(5):637–44.

    Article  CAS  PubMed  Google Scholar 

  52. Lagesen K, Hallin P, Rødland EA, Stærfeldt H-H, Rognes T, Ussery DW. RNAmmer: consistent and rapid annotation of ribosomal RNA genes. Nucleic Acids Res. 2007;35(9):3100–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  53. Lowe TM, Eddy SR. tRNAscan-SE: a program for improved detection of transfer RNA genes in genomic sequence. Nucleic Acids Res. 1997;25(5):955–64.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  54. Greiner S, Lehwark P, Bock R. OrganellarGenomeDRAW (OGDRAW) version 1.3. 1: expanded toolkit for the graphical visualization of organellar genomes. Nucleic Acids Res. 2019;47(W1):W59–64.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  55. Beier S, Thiel T, Münch T, Scholz U, Mascher M. MISA-web: a web server for microsatellite prediction. Bioinformatics. 2017;33(16):2583–5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  56. Kurtz S, Choudhuri JV, Ohlebusch E, Schleiermacher C, Stoye J, Giegerich R. REPuter: the manifold applications of repeat analysis on a genomic scale. Nucleic Acids Res. 2001;29(22):4633–42.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  57. Benson G. Tandem repeats finder: a program to analyze DNA sequences. Nucleic Acids Res. 1999;27(2):573–80.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  58. Edgar RC. MUSCLE: multiple sequence alignment with improved accuracy and speed. In: Computational Systems Bioinformatics Conference; 2004.

    Google Scholar 

  59. Stamatakis A. RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics. 2014;30(9):1312–3.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  60. Yang Z. PAML 4: phylogenetic analysis by maximum likelihood. Mol Biol Evol. 2007;24(8):1586–91.

    Article  CAS  PubMed  Google Scholar 

  61. Yang Z, Wong WS, Nielsen R. Bayes empirical Bayes inference of amino acid sites under positive selection. Mol Biol Evol. 2005;22(4):1107–18.

    Article  CAS  PubMed  Google Scholar 

Download references


We are very grateful to the editors and reviewers for critically evaluating the manuscript and providing constructive comments for its improvement.


This research was funded by the National Natural Science Foundation of China (42206117, 42176154, 42006144), the NSFC-Shandong Joint Fund (U1806213), and the National Key R&D Program of China (2019YFD0901204).

Author information

Authors and Affiliations



Jun Chen, Yu Zang and Xuexi Tang conducted the project. Song Xue, Liang Shuo and Meiling Zhu supervised the study. Jun Chen, Shuai Shang performed analysis of the data. Jun Chen and Ying Wang wrote the manuscript. The author(s) read and approved the final manuscript.

Corresponding authors

Correspondence to Ying Wang or Xuexi Tang.

Ethics declarations

Ethics approval and consent to participate

All the experiments on plants were conducted in accordance with the guidelines of Ocean University of China.

Consent for publication

Not applicable.

Competing interests

All authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Chen, J., Zang, Y., Liang, S. et al. Comparative analysis of mitochondrial genomes reveals marine adaptation in seagrasses. BMC Genomics 23, 800 (2022).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • Seagrass
  • Mitochondrial genome
  • Comparative genomics
  • Adaptive evolution