Next Article in Journal
Genome Editing Tools in Plants
Previous Article in Journal
Comparative Transcriptome Analysis of Male and Female Conelets and Development of Microsatellite Markers in Pinus bungeana, an Endemic Conifer in China
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Genes
Volume 8
Issue 12
10.3390/genes8120395
genes-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Comparative Analysis of Four Calypogeia Species Revealed Unexpected Change in Evolutionarily-Stable Liverwort Mitogenomes

by
Monika Ślipiko
1,*,
Kamil Myszczyński
1,
Katarzyna Buczkowska-Chmielewska
2,
Alina Bączkiewicz
2,
Monika Szczecińska
1 and
Jakub Sawicki
1
1
Department of Botany and Nature Protection, University of Warmia and Mazury in Olsztyn, Plac Łódzki 1, 10-727 Olsztyn, Poland
2
Department of Biology, Institute of Experimental Biology, Adam Mickiewicz University in Poznań, Umultowska 89, 61-614 Poznań, Poland
*
Author to whom correspondence should be addressed.
Genes 2017, 8(12), 395; https://0-doi-org.brum.beds.ac.uk/10.3390/genes8120395
Submission received: 11 October 2017 / Revised: 8 December 2017 / Accepted: 12 December 2017 / Published: 19 December 2017
(This article belongs to the Section Plant Genetics and Genomics)
Browse Figures
Versions Notes

Abstract

:
Liverwort mitogenomes are considered to be evolutionarily stable. A comparative analysis of four Calypogeia species revealed differences compared to previously sequenced liverwort mitogenomes. Such differences involve unexpected structural changes in the two genes, cox1 and atp1, which have lost three and two introns, respectively. The group I introns in the cox1 gene are proposed to have been lost by two-step localized retroprocessing, whereas one-step retroprocessing could be responsible for the disappearance of the group II introns in the atp1 gene. These cases represent the first identified losses of introns in mitogenomes of leafy liverworts (Jungermanniopsida) contrasting the stability of mitochondrial gene order with certain changes in the gene content and intron set in liverworts.

1. Introduction

Group I and group II introns, next to spliceosomal and transfer RNA (tRNA) introns, belong to four main types of introns divided on the basis of splicing mechanism [1]. Although two transesterification reactions are used by group I and group II introns for their splicing, the reaction mechanisms are different. As a result, group I introns are removed in a linear form, and some of them can circularize, whereas group II introns are released as a lariat [2]. Both discussed intron groups are also known as mobile elements. Their mobility is possible thanks to internal encoded enzymes, but the movement mechanism in each of the groups is different. Group I introns can proliferate by a DNA-mediated homing mechanism, where intron-encoded endonucleases play a key role. In contrast, retrohoming (RNA-mediated mechanism) is used by group II introns for propagation [3]. Internal encoded enzymes of this intron class, maturase, reverse transcriptase and endonuclease, enable retrohoming [4].
In addition to the differences in splicing and mobility mechanisms, group I and group II introns also have a unique structure. Group I introns have a characteristic RNA fold consisting of 10 elements denoted from P1 to P10 [5], while group II introns typically have a secondary structure consisting of six double-helical domains [2,6].
In plant mitochondrial genomes (chondriomes), group II introns occur more commonly than group I introns [1]. Group II introns are also present in archaebacteria, bacteria and plastids, but are absent in a nuclear genome, unless a sequence with this type of introns was transferred from mitochondrion to nucleus [7,8]. On the other hand, group I introns are common in nuclear ribosomic RNA (rRNA)-encoding genes, frequent in fungal mitochondria and are also present in plant organellar genomes, bacteria and viruses [1,9]. Intron contents of the mitochondrial genome in angiosperms are reported to be rather stable [10,11], as well as in hornworts, mosses and in the earliest land plants: liverworts. On the other hand, intron contents among these four distinct lineages are reported to be significantly different [11,12]. The tendency towards a stable number of introns in the mitogenome among bryophytes is in accordance with the hypothesis that their mitochondrial genomes are slowly evolving [13].
Retroprocessing is the most frequently reported mechanism for removal of introns [14,15,16,17,18]. According to this model, also known as a reverse transcriptase-mediated model (RT-mediated model), spliced mRNA is reverse-transcribed, and then, the emergent intronless complementary DNA (cDNA) is integrated into the genome by homologous recombination [14,15,16,17]. A loss of introns is thus associated with a loss of editing sites in the genome [18]. Under the standard type of this mechanism, the introns located at 3’ ends of genes are more likely to be deleted, because reverse transcriptase polymerizes the RNA template from 3’ ends to 5’ends of genes and frequently dissociates from the template without completing its rewriting [14,15,19]. Although reports of cDNA production in vivo are very rare [15,20,21], reverse transcription has been widely discussed in gene evolution [18], and a 3’ bias of intron loss has been reported in some organisms [16], which may indicate an RT-mediated model of intron loss. Intron losses from quillwort Isoëtes engelmannii [22], and gymnosperms are the strongest [23] evidence of retroprocessing in the plant mitochondrial genome.
Other possible mechanisms of intron loss include genomic deletion, exonization, horizontal gene transfer (HGT) and gene conversion. In the case of genomic deletion, introns are excised imprecisely [24], resulting in the removal of adjacent exonic sequences or retaining small fragments of introns, which are exonized [10]. This exonization may also involve an entire intron, which is no longer cut out of the transcript, but preserved in the mature mRNA and translated [25]. Disorders in intron cutting may be caused by creating alternative splice sites, which are the result of genomic insertions or point mutations in a DNA sequence. In this way, the splicing system includes new sequences (here introns) as exons or elongated existing exons [26]. Intron loss by deletion has been reported only in the case of Petunia [27], whereas exonization has never been evidenced in plant mitochondrial genome [10]. Intron loss can be also explained by HGT of an intron-less gene and the following gene conversion with an intron-including gene. The above mechanism in the plant mitogenome has only been observed in Magnolia tripetala to date [10].
The genus Calypogeia Raddi belongs to leafy liverworts (Jungermanniopsida) and consists of about 90 described species [28]. The most characteristic taxonomical feature of this genus is the presence of oil bodies and their color, shape and pattern of distribution in the leaf and underleaf [29]. This genus is the most widespread (but hardly cosmopolitan) genus of the family Calypogeiaceae. Although it includes species with wide Holarctic distribution (e.g., Calypogeia suecica, Calypogeia sphagnicola, Calypogeia neesiana), it also contains neotropical species occurring in Central and South America (e.g., Calypogeia biapiculata, Calypogeia laxa, Calypogeia miquelii) [30,31]. The family Calypogeiaceae is easily recognized, but species identification, due to the presence of environmentally-induced modifications and atypical forms and the frequent absence of sporophyte, can often lead to misidentification [32,33]. Some morphologically similar (cryptic) species are possible to distinguish using molecular markers [34].
Genomic work on liverworts mainly concerns simple (metzgeriid) and complex (marchantiid) thalloid liverworts, and there are still only a few sequenced mitogenomes [12,35,36,37] and plastomes [38,39,40,41,42,43].
In this paper, we introduce for the first time the mitochondrial genome structure of the genus Calypogeia. The four presented Calypogeia mitogenomes are the first sequenced mitochondrial genomes of leafy liverworts. The main goal of our research was to examine the stability of the mitochondrial genome structure among liverworts mentioned in the literature. Unexpectedly, we discovered a lack of introns in the genes atp1 and cox1. It is the first case of intron loss within Jungermanniopsida leafy liverworts. We also discuss the possible mechanisms of intron disappearance.

2. Materials and Methods

2.1. Genome Sequencing, Assembly and Annotation

Total genomic liverwort DNA was extracted using ZR Plant/Seed DNA MiniPrepTM kit (Zymo Research Corp., Irvine, CA, USA) following the manufacturer’s recommendations. The DNA libraries of Calypogeia integristipula, Calypogeia fissa ssp. fissa and Calypogeia suecica were sequenced using the Illumina HiSeq 2000 platform (Illumina, San Diego, CA, USA) to generate 100 bp paired-end reads, and the DNA library of Calypogeia arguta was sequenced using HiSeqX (Illumina) to generate 150 bp paired-end reads (Table 1). After sequencing, the reads were cleaned by removing the adaptor sequences and low-quality reads with ambiguous sequences. Afterwards, the reads were mapped to the reference mitochondrial genome of Aneura pinguis using Geneious R8 software [44]. Then, gaps, in the obtained scaffold of the mitogenome sequence, were filled using an iterative fine-tuning step with minimum overlap and minimum overlap identity parameters set to 50 bp and 98%, respectively. The same genome assembly strategy was used for all Calypogeiaceae species analyzed in this study. The number of sequence reads and coverage depth for C. arguta deviates from other species since C. arguta has been deeply sequenced for nuclear genome research.
Genes were identified and annotated based on the closest known mitochondrial genomes of related species to the Calypogeiaceae, i.e., Aneura mirabilis, A. pinguis, Pellia endiviifolia, Ptilidium pulcherrimum and Marchantia polymorpha. Predictions were made using Geneious R8 software [44] and the BLAST tool [45]. Annotated sequences of the aforementioned Calypogeiaceae species were deposited in GenBank with the accession numbers specified in Table 1. A circular genome map, in order to visualize the structure of Calypogeiaceae mitogenome, was created on the basis of C. fissa ssp. fissa mitogenome using OGDraw software [46]. A MAUVE plot demonstrating the stability of the gene order in liverworts was made using Geneious R8 software [44] (Figure S1).

2.2. Prediction of RNA Editing Sites

Protein-coding sequences of cox1 and atp1 genes of four species (A. pinguis, M. polymorpha, Pleurozia purpurea and Treubia lacunosa) were downloaded from the NCBI (National Center for Biotechnology Information) database [47]. To predict editing sites and evaluate their editing rates within cox1 and atp1 genes, both the PREP-Mt [48] and PREPACT 2.0 [49] tools were used with a cutoff value of 0.8 with four downloaded and Calypogeia fissa ssp. fissa, Tritomaria quinquedentata cox1 and atp1 sequences.

3. Results and Discussion

3.1. Mitogenome Structure

The mitochondrial genome of Calypogeia is 159,061–163,057 bp in length (Table 1, Figure 1) and is slightly shorter than the closest related species with a known mitogenome structure: P. purpurea (168,526 bp) [36]. The length of the mitochondrial genome of Calypogeia is most similar to the mitogenome of A. pinguis (164,989 bp) [37]. The two other known mitochondrial genomes of the liverwort species, T. lacunosa [12] and M. polymorpha [38], differ more and are composed of 151,983 and 186,609 bp, respectively. The GC content in the studied genome (45%) is similar to other liverworts (42–45%) [12,36].
Seventy genes have been identified in the Calypogeia mitogenome: 42 protein-coding genes 25 tRNAs and three rRNAs (Table 2). Chondriomes of liverworts, similar to moss and hornwort mitogenomes, are reported to be rather static in gene content and order [12,13,50] and even with respect to pseudogene contents and retroposed pseudogene pieces [12]. The gene order in Calypogeia mitogenome is identical to the four aforementioned liverwort mitogenomes, although the gene content is slightly different between them. Mitogenomes of M. polymorpha, P. purpurea, A. pinguis and Calypogeia are very similar in gene composition. The differences occur mainly in the content of transfer RNA genes. M. polymorpha has only one copy of the trnRucu, but contains two more tRNA genes: trnRucg and trnTggu. One copy of the trnRucu has probably given rise to trnRucg [12]. The trnT gene is a part of the trnA-trnT-nad7 gene cluster, whose different forms were identified by Wahrmund et al. [51] in liverwort mitogenome evolution. In leafy (jungermanniid) liverworts and in simple thalloid (metzgeriid), trnT is lacking between trnA and nad7, whereas in Blasia pusilla, representing a sister lineage to all other complex thalloid (marchantiid) liverworts, this gene occurs in conserved Chara-like version. The trnT gene is also present in M. polymorpha and other complex thalloid (marchantiid) liverworts, but its sequence is inverted compared to Blasia. Furthermore, in the A. pinguis mitochondrial genome, there is only one copy of the trnYgua in contrast to other liverwort mitogenomes [36]. Another difference concerns the rtl gene. In all aforementioned mitogenomes, this gene is functional with nucleotide sequence similarity > 84%, whereas in P. purpurea, it may be a pseudogene because of the high level of sequence divergence and several indels in the open reading frame [36]. However, a big part of the reading frame in rtl is intact, so this gene in P. purpurea may be still functional. The other dissimilarities occur between mitogenomes of the above four species and T. lacunosa. In the mitochondrial genome of the latter, either some genes of the cytochrome c biogenesis (ccmC, ccmFN) are missing or some of them are pseudogenized (ccmB, ccmFC). Another conspicuous dissimilarity concerns the nad7 gene, which is only functional in T. lacunosa [12]. In most hornworts and liverworts, this gene is missing or occurs as a pseudogene with a degenerated structure [52,53,54], which is reflected in Calypogeia mitogenome and in the other sequenced mitochondrial genomes of liverworts. The only liverwort species with functional nad7 are Treubia and Haplomitrium [52] belonging to Haplomitriopsida, a sister clade to the rest of the liverworts: Marchantiopsida and Jungermanniopsida [55] with the inactive nad7 gene.
Ninety-four spacers and 24 introns (Table 3) have been found across the entire mitochondrial genome. The different length (1000–1300 bp) of the nad5-nad4 spacer, containing the inverted sequence of the second cob intron in M. polymorpha, was recognized in different liverwort groups. Almost the entire cob intron sequence is inserted in the nad5-nad4 spacer in marchantiid, while an internal region of this intron sequence copy is deleted among metzgeriid and jungermanniid taxa [56]. The above findings are also supported by the current study, because the nad5-nad4 spacer in Calypogeia has a structure typical of jungermanniid liverworts.
In the Calypogeia mitogenome, 22 introns (Table 3) are located in protein-coding genes, one in the rrn26 gene and one in the trnS gene. Seven genes (nad2, nad3, nad4, nad5, rpl2, rps14 and atp9) contain one intron. The genes cox2, cox3 and nad4L have two introns, whereas the coding sequence of the cob gene is divided into three introns. The largest number of introns is localized in the cox1 gene. However, surprisingly, only six introns occur in Calypogeia, whereas nine introns exist in the cox1 gene of thalloid liverworts sequenced to date [12,35,36,37]. The cox1 gene of Calypogeia lacks the cox1i395g1, cox1i624g1 and cox1i729g1 introns. The atp1 gene has also lost two introns (atp1i989g2 and atp1i1050g2) and become in Calypogeia intronless (Figure 2). The CDS (protein-coding sequence) structure of both genes has not been affected. The intron set among species within each of the three major lineages of bryophytes is reported to vary slightly [12]. The intron number in the previously sequenced liverwort mitogenomes ranges from 28 in T. lacunosa to 30 in M. polymorpha, not including introns in pseudogenes. In liverworts, three cases of changes in the intron number have been detected to date (marked with numbers (1)–(3)). The rrn18i1065gII intron present in M. polymorpha is lacking in P. purpurea [35,36], A. pinguis [37] and in the examined genus Calypogeia (1). The intron set of T. lacunosa differs most compared to the other liverworts. Apart from the absence of the mentioned rrn18i1065gII intron (1), it lacks one intron in the nad4L gene (2) and, as in Calypogeia, two introns of the atp1 gene (3) [12].

3.2. Possible Mechanisms of Intron Losses in Two Genes of Calypogeia

Among the currently sequenced liverwort mitogenomes, the Calypogeia mitochondrial genomes contain the fewest introns (24 introns). Three group I introns of the cox1 gene and two group II introns in atp1 are missing. Why have they disappeared? Deletion, as one of the intron loss mechanisms, can be ruled out. Introns in Calypogeia are precisely removed, and there are no intron fragments left or a deletion of small adjoining exon pieces. Exonization can also be rejected, because the introns physically disappeared from the cox1 gene and the exon structures in both atp1 and cox1 genes remain intact. Furthermore, because of the large size of these introns (over 1000 bp), it would be unlikely [10]. Horizontal gene transfer and gene conversion also seem improbable since no chimeric structure in any place of the genes was noted. However, all cox1 and atp1 exons among liverworts are so similar, that even if horizontal transport and gene conversion took place, it would be hard to notice it.
Taking into account the precision in the intron cut [10], the most probable mechanism of intron losses in Calypogeia genes is retroprocessing. However, this process in the cox1 gene must have occurred at least twice, because non-intronic fragments are located in two places of the gene and are separated by a group II intron (Figure 2). This is the so-called localized retroprocessing, affecting only a part of the native gene [10,18]. Retroprocessing usually involves the removal of the adjacent edited sites, together with introns [18,22,59,60,61]. Computational analysis revealed 10 RNA editing sites in the cox1 gene and five in the atp1 gene of Calypogeia (Figure 3). Despite the intron deletions, nearby nucleotide positions that require editing remain. While the remaining two editing sites in the fourth exon of the cox1 gene can be explained by embedding a cDNA fragment into a region between, but not including, the editing sites [18], the maintenance of edited sites in the fifth exon cannot be explained by this mechanism, because positions needing editing occur in the middle of this exon. Similarly, the loss of introns in the atp1 gene does not seem to be related to the disappearance of the editing sites, because editing sites occur very close to the previous intron-exon boundary. Perhaps post-transcriptional modification of RNA in Calypogeia takes place in two successive stages: splicing and editing. Perhaps immediately after splicing, but before editing, RNA is reverse transcribed, and then, a partially-processed cDNA fragment undergoes conversion with the native intron-bearing gene. As a result, introns are removed, but editing sites remain [10,18]. The loss of introns without concurrent loss of flanking editing sites has recently been reported in ferns [62]. Edited sites from the intron-missing genes of Calypogeia may also not be removed because they are crucial for the excision of the remaining introns [10]. On the other hand, the introns can also play an important role in the correct splicing of the other ones [63]. Perhaps this is why the cox1i511g2 intron in the cox1 gene of Calypogeia is preserved.
Apart from the precise intron excision, stronger evidence of retroprocessing is the loss of introns at the 3 ‘end of the gene. However, in the gene cox1 of Calypogeia, introns have been lost in the center of the gene. A similar pattern of intron disappearance was observed by Nielsen et al. [64] who suggested that this may be caused either by a mutational mechanism (e.g., reverse transcription primed internally) or selective pressure to maintain introns near the 5′ and 3′ ends of genes. In T. quinquedentata, whose cox1 gene also lacks the cox1i395gI, cox1i624gI and cox1i729gI introns, a loss of the cox1i1116g1 intron (3′-biased intron losses) was observed [65] (Figure 2), but the number and location of editing sites was very similar to those in cox1 of Calypogeia (Figure 3). Thus, despite a visible lack of intron loss bias towards the 3’ end of the cox1 gene in Calypogeia and the remaining editing sites, it is very likely that retroprocessing takes place in this gene.

4. Conclusions

In summary, the structure of the Calypogeia mitogenome is in line with reports of the stability of mitochondrial genomes in bryophyte lineages. The gene order is identical to other liverworts, while gene content fits the patterns that emerged in liverwort evolution. Dissimilarities in gene content among liverwort mitogenomes occur with regard to tRNA genes and protein-coding genes such as rtl, nad7, ccmB and ccmFC that are either functional or pseudogenized. Furthermore, Calypogeia species, like other sequenced liverworts, have functional ccmC and ccmFN genes, that are absent in T. lacunosa. The most unexpected difference occurs in the intron set. The cox1 gene lacks three introns and atp1 two. This is the first notification of intron losses in Jungermanniopsida leafy liverworts, although the disappearance of introns in atp1 has already been reported in T. lacunosa. The mechanism responsible for intron disappearance seems to be retroprocessing: one-step retroprocessing in atp1 and two-step localized retroprocessing in the cox1 gene. The above findings indicate that the gene order in liverwort mitogenomes is stable, but the gene and intron contents may vary.

Supplementary Materials:

The following are available online at www.mdpi.com/2073-4425/8/12/395/s1. Figure S1: Mauve alignment of six mitogenomes of liverworts.

Acknowledgments

The sequencing of liverwort mitogenomes was financially supported by The National Science Center Kraków, Poland: Calypogeia mitogenome, Grant No. 2015/19/B/NZ8/03970, Aneura pinguis mitogenome, Grant No. 2016/21/B/NZ8/03325, Tritomaria quinquedentata mitogenome, Grant No. 2017/01/X/NZ8/01094.

Author Contributions

M.Ś. initiated the study, analyzed the data and wrote the paper. K.M. assembled the mitogenomes, performed the prediction of editing sites and wrote the paper. K.B.-C. and A.B. collected and delimited specimens and performed DNA extraction. M.S. prepared the figures and tables. J.S. supervised the project and revised the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References and Note

  1. Lang, B.F.; Laforest, M.-J.; Burger, G. Mitochondrial introns: A critical view. Trends Genet. 2007, 23, 119–125. [Google Scholar] [CrossRef] [PubMed]
  2. Saldanha, R.; Mohr, G.; Belfort, M.; Lambowitz, A.M. Group I and group II introns. FASEB J. 1993, 7, 15–24. [Google Scholar] [PubMed]
  3. Eickbush, T.H. Mobile introns: Retrohoming by complete reverse splicing. Curr. Biol. 1999, 9, R11–R14. [Google Scholar] [CrossRef]
  4. Simon, D.M.; Clarke, N.A.; McNeil, B.A.; Johnson, I.; Pantuso, D.; Dai, L.; Chai, D.; Zimmerly, S. Group II introns in eubacteria and archea: ORF-less introns and new varieties. RNA 2008, 14, 1704–1713. [Google Scholar] [CrossRef] [PubMed]
  5. Haugen, P.; Reeb, V.; Lutzoni, F.; Bhattacharya, D. The evolution of homing endonuclease genes and group I introns in nuclear rDNA. Mol. Biol. Evol. 2004, 21, 129–140. [Google Scholar] [CrossRef] [PubMed]
  6. Michel, F.; Ferat, J.L. Structure and activities of group II introns. Annu. Rev. Biochem. 1995, 64, 435–461. [Google Scholar] [CrossRef] [PubMed]
  7. Knoop, V.; Brennicke, A. Promiscuous mitochondrial group II intron sequences in plant nuclear genomes. J. Mol. Evol. 1997, 39, 144–150. [Google Scholar]
  8. Lin, X.; Kaul, S.; Rounsley, S.; Shea, T.P.; Benito, M.I.; Town, C.D.; Fujii, C.Y.; Mason, T.; Bowman, C.L.; Barnstead, M.; et al. Sequence and analysis of chromosome 2 of the plant Arabidopsis thaliana. Nature 1999, 402, 761–768. [Google Scholar] [CrossRef] [PubMed]
  9. Haugen, P.; Simon, D.M.; Bhattacharya, D. The natural history of group I introns. Trends Genet. 2005, 21, 111–119. [Google Scholar] [CrossRef] [PubMed]
  10. Hepburn, N.J.; Schmidt, D.W.; Mower, J.P. Loss of two introns from the Magnolia tripetala mitochondrial cox2 gene implicates horizontal gene transfer and gene conversion as a novel mechanism of intron loss. Mol. Biol. Evol. 2012, 29, 3111–3120. [Google Scholar] [CrossRef] [PubMed]
  11. Pruchner, D.; Nassal, B.; Schindler, M.; Knoop, V. Mosses share mitochondrial group II introns with flowering plants, not with liverworts. Mol. Genet. Genom. 2001, 266, 608–613. [Google Scholar] [CrossRef] [PubMed]
  12. Liu, Y.; Xue, J.-Y.; Wang, B.; Li, L.; Qiu, Y.-L. The mitochondrial genomes of the early land plants Treubia lacunosa and Anomodon rugelii: Dynamic and conservative evolution. PLoS ONE 2011, 6, e25836. [Google Scholar] [CrossRef] [PubMed]
  13. Liu, Y.; Wang, B.; Li, L.; Qiu, Y.-L.; Xue, J.-Y. Conservative and dynamic evolution of mitochondrial genomes in early land plants. In Genomics of Chloroplasts and Mitochondria; Bock, R., Knoop, V., Eds.; Springer: Dordrecht, The Netherlands, 2012; pp. 159–174. [Google Scholar]
  14. Fink, G.R. Pseudogenes in yeast? Cell 1987, 49, 5–6. [Google Scholar] [CrossRef]
  15. Derr, L.K.; Strathern, J.N. A role for reverse transcripts in gene conversion. Nature 1993, 361, 170–173. [Google Scholar] [CrossRef] [PubMed]
  16. Mourier, T.; Jeffares, D.C. Eucaryotic intron loss. Science 2003, 300, 1393. [Google Scholar] [CrossRef] [PubMed]
  17. Cohen, N.E.; Shen, R.; Carmel, L. The role of reverse transcriptase in intron gain and loss machanisms. Mol. Biol. Evol. 2013, 29, 179–186. [Google Scholar] [CrossRef] [PubMed]
  18. Cuenca, A.; Ross, G.; Graham, S.W.; Barrett, C.F.; Davis, J.I.; Seberg, O.; Petersen, G. Localized retroprocessing as a model of intron loss in the plant mitochondrial genome. Genome Biol. Evol. 2016, 8, 2176–2189. [Google Scholar] [CrossRef] [PubMed]
  19. Zhang, L.-Y.; Yang, Y.-F.; Niu, D.-K. Evaluation of models of the mechanisms underlying intron loss and gain in Aspergillus fungi. J. Mol. Evol. 2010, 71, 364–373. [Google Scholar] [CrossRef] [PubMed]
  20. Dhellin, O.; Maestre, J.; Heidmann, T. Functional differences between the human LINE retrotrasposon and retroviral reverse transcriptases for in vivo mRNA reverse transcription. EMBO J. 1997, 16, 6590–6602. [Google Scholar] [CrossRef] [PubMed]
  21. Odom, O.W.; Herrin, D.L. Reverse transcription of spliced psbA mRNA in Chlamydomonas spp. and its possible role in evolutionary intron loss. Mol. Biol. Evol. 2013, 30, 2666–2675. [Google Scholar] [CrossRef] [PubMed]
  22. Grewe, F.; Herres, S.; Viehöver, P.; Polsakiewicz, M.; Weisshaar, B.; Knoop, V. A unique transcriptome: 1782 positions of RNA editing alter 1406 codon indetities in mitochondrial mRNAs of the lycophyte Isoëtes engelmannii. Nucleic Acids Res. 2011, 39, 2890–2902. [Google Scholar] [CrossRef] [PubMed]
  23. Ran, J.H.; Gao, H.; Wang, X.Q. Fast evolution of the retroprocessed mitochondrial rps3 gene in Conifer II and further evidence for the phylogeny of gymnosperms. Mol. Phylogenet. Evol. 2010, 54, 136–149. [Google Scholar] [CrossRef] [PubMed]
  24. Coulombe-Huntington, J.; Majewski, J. Intron loss and gain in Drosophila. Mol. Biol. Evol. 2007, 24, 2842–2850. [Google Scholar] [CrossRef] [PubMed]
  25. Wang, W.; Zheng, H.; Yang, S.; Yu, S.; Li, J.; Jiang, H.; Su, J.; Yang, L.; Zhang, J.; McDermott, J.; et al. Origin and evolution of new exons in rodents. Genome Res. 2005, 15, 1258–1264. [Google Scholar] [CrossRef] [PubMed]
  26. Schmitz, J.; Brosius, J. Exonization of transposed elements: A challenge and opportunity for evolution. Biochimie 2011, 93, 1928–1934. [Google Scholar] [CrossRef] [PubMed]
  27. Pruitt, K.D.; Hanson, M.R. Cytochrome oxidase subunit II sequences in Petunia mitochondria: Two intron-containing genes and an intron-less pseudogene associated with cytoplasmic male sterility. Curr. Genet. 1989, 16, 281–291. [Google Scholar] [CrossRef] [PubMed]
  28. Schuster, R.M. The Hepaticae and Anthocerotae of North America East of the Hundredth Meridian; Columbia Univeristy Press: New York, NY, USA, 1969; Volume 2, p. 1062. ISBN 0231089821. [Google Scholar]
  29. Buch, H. Vorarbeiten zu einer Lebermoosflora Fenno-Scandias III. Die Gattung Calypogeia Raddi. Mem. Soc. F. Fl. Fenn. 1935, 11, 197–214. [Google Scholar]
  30. Bischler, H. The genus Calypogeia Raddi in Central and South America. I-III. Candollea 1963, 18, 19–128. [Google Scholar]
  31. Dauphin, G. Catalogue of Costa Rican Hepaticae and Anthocerotae. Trop. Bryol. 2005, 26, 141–218. [Google Scholar]
  32. Szweykowski, J. Species problems and taxonomic methods in Bryophytes. In The New Manual of Bryology; Schuster, R.M., Ed.; Hattori Botanical Laboratory: Nichinan, Japan, 1984; Volume 2, pp. 1130–1171. [Google Scholar]
  33. Buczkowska, K.; Odrzykoski, I.J.; Chudzińska, E. Delimitation of some European species of Calypogeia Raddi (Hepaticae, Jungermanniales) based on cytological characters of oil modies and multienzyme phenotype. Nova Hedwigia 2004, 78, 147–163. [Google Scholar] [CrossRef]
  34. Buczkowska, K.; Dabert, M. The development of species-specific SCAR markers for delimitation of Calypogeia species. J. Bryol. 2011, 33, 291–299. [Google Scholar] [CrossRef]
  35. Oda, K.; Yamato, K.; Ohta, E.; Nakamura, Y.; Takemura, M.; Nozato, N.; Akashi, K.; Kanegae, T.; Ogura, Y.; Kohchi, T.; et al. Gene organization deduced from the complete sequence of liverwort Marchantia polymorpha mitochondrial DNA—A primitive form of plant mitochondrial genome. J. Mol. Biol. 1992, 223, 1–7. [Google Scholar] [CrossRef]
  36. Wang, B.; Xue, J.; Li, L.; Liu, Y.; Qiu, Y.-L. The complete mitochondrial genome sequence of the liverwort Pleurozia purpurea reveals extremely conservative mitochondrial genome evolution in liverworts. Curr. Genet. 2009, 55, 601. [Google Scholar] [CrossRef] [PubMed]
  37. Myszczyński, K.; Bączkiewicz, A.; Szczecińska, M.; Buczkowska, K.; Kulik, T.; Sawicki, J. The complete mitochondrial genome of the cryptic species C of Aneura pinguis. Mitochondrial DNA Part A 2017, 28, 112–113. [Google Scholar] [CrossRef] [PubMed]
  38. Ohyama, K.; Fukuzawa, H.; Kohchi, T.; Shirai, H.; Sano, T.; Sano, S.; Umesono, K.; Shiki, Y.; Takeuchi, M.; Chang, Z.; et al. Chloroplast gene organization deduced from complete sequence of liverwort Marchantia polymorpha chloroplast DNA. Nature 1984, 322, 572–574. [Google Scholar] [CrossRef]
  39. Ohyama, K.; Fukuzawa, H.; Kohchi, T.; Sano, T.; Sano, S.; Shirai, H.; Umesono, K.; Shiki, Y.; Takeuchi, M.; Chang, Z.; et al. Structure and organization of Marchantia polymorpha chloroplast genome. I. Cloning and gene identification. J. Mol. Biol. 1988, 203, 281–298. [Google Scholar] [CrossRef]
  40. Wickett, N.J.; Zhang, Y.; Hansen, S.K.; Roper, J.M.; Kuehl, J.V.; Plock, S.A.; Wolf, P.G.; dePamphilis, C.W.; Boore, J.L.; Goffinet, B. Functional gene losses occur with minimal size reduction in the plastid genome of the parasitic liverwort Aneura mirabilis. Mol. Biol. Evol. 2008, 25, 393–401. [Google Scholar] [CrossRef] [PubMed]
  41. Forrest, L.L.; Wickett, N.J.; Cox, C.J.; Goffinet, B. Deep sequencing of Ptilidium (Ptilidiaceae) suggests evolutionary stasis in liverwort plastid genome structure. Plant Ecol. Evol. 2011, 144, 29–43. [Google Scholar] [CrossRef]
  42. Grosche, C.; Funk, H.T.; Maier, U.G.; Zauner, S. The chloroplast genome of Pellia endiviifolia: Gene content, RNA-editing pattern, and the origin of chloroplast editing. Genome Biol. Evol. 2012, 4, 1349–1357. [Google Scholar] [CrossRef] [PubMed]
  43. Myszczyński, K.; Bączkiewicz, A.; Buczkowska, K.; Ślipiko, M.; Szczecińska, M.; Sawicki, J. The extraordinary variation of the organellar genomes of the Aneura pinguis revealed advanced cryptic speciation of the early land plants. Sci. Rep. 2017, 7, 9804. [Google Scholar] [CrossRef] [PubMed]
  44. Kearse, M.; Moir, R.; Wilson, A.; Stones-Havas, S.; Cheung, M.; Sturrock, S.; Buxton, S.; Cooper, A.; Markowitz, S.; Duran, C.; et al. Geneious Basic: An integrated and extendable desktop software platform for the organization and analysis of sequence data. Bioinformatics 2012, 28, 1647–1649. [Google Scholar] [CrossRef] [PubMed]
  45. Johnson, M.; Zaretskaya, I.; Raytselis, Y.; Merezhuk, Y.; McGinnis, S.; Madden, T.L. NCBI BLAST: A better web inerface. Nucleic Acids Res. 2008, 36, W5–W9. [Google Scholar] [CrossRef] [PubMed]
  46. Lohse, M.; Drechsel, O.; Kahlau, S.; Bock, R. OrganellarGenomeDRAW—A suite of tools for generating physical maps of plastid and mitochondrial genomes and visualizing expression data sets. Nucleic Acids Res. 2013, 41, W575–W581. [Google Scholar] [CrossRef] [PubMed]
  47. National Center for Biotechnology Information (NCBI) [Internet]. Bethesda (MD): National Library of Medicine (US), National Center for Biotechnology Information. 1988. Available online: www.ncbi.nlm.nih.gov (accessed on 10 March 2017).
  48. Mower, J.P. The PREP Suite: Predictive RNA editors for plant mitochondrial genes, chloroplast genes and user-defined alignments. Nucleic Acids Res. 2009, 37, W253–W259. [Google Scholar] [CrossRef] [PubMed]
  49. Lenz, H.; Knoop, V. PERPACT 2.0: Predicting C-to-U and U-to-C RNA editing in organelle genome sequences with multiple references and curated RNA editing annotation. Bioinform. Biol. Insights 2013, 7, 1–19. [Google Scholar] [CrossRef] [PubMed]
  50. Sawicki, J.; Szczecińska, M.; Kulik, T.; Gomolińska, A.M.; Plášek, V. The complete mitochondrial genome of the epiphytic moss Orthotrichum speciosum. Mitochondrial DNA Part A. 2016, 27, 1709–1710. [Google Scholar]
  51. Wahrmund, U.; Groth-Malonek, M.; Knoop, V. Tracing plant mitochondrial DNA evolution: Rearrangements of the ancient mitochondrial gene cluster trnA-trnT-nad7 in liverwort phylogeny. J. Mol. Evol. 2008, 66, 621–629. [Google Scholar] [CrossRef] [PubMed]
  52. Groth-Malonek, M.; Wahrmund, U.; Polsakiewicz, M.; Knoop, V. Evolution of a pseudogene: Exclusive survival of a functional mitochondria nad7 gene supports Haplomitrium as the earliest liverwort lineage and proposes a secondary loss of RNA editing in Marchantiidae. Miol. Biol. Evol. 2007, 24, 1068–1074. [Google Scholar] [CrossRef] [PubMed]
  53. Li, L.; Wang, B.; Liu, Y.; Qiu, Y.-L. The complete mitochondrial genome sequence of the hornwort Megaceros aenigmaticus shows a mixed mode of conservative yet dynamic evolution in early land plant mitochondrial genomes. J. Mol. Evol. 2009, 68, 665–678. [Google Scholar] [CrossRef] [PubMed]
  54. Xue, J.Y.; Liu, Y.; Li, L.; Wang, B.; Qiu, Y.L. The complete mitochondria genome sequence of the hornwort Phaeoceres laevis: Retention of many ancient pseudogenes and conservative evolution of mitochondrial genomes in hornworts. Curr. Genet. 2010, 56, 53–62. [Google Scholar] [CrossRef] [PubMed]
  55. Crandall-Stotler, B.; Stotler, R.E.; Long, D.G. Phylogeny and classification of the Marchantiophyta. Edinb. J. Bot. 2009, 66, 155–198. [Google Scholar] [CrossRef]
  56. Groth-Malonek, M.; Rein, T.; Wilson, R.; Groth, H.; Heinrichs, J.; Knoop, V. Different fates of two mitochondrial gene spacers in early land plant evolution. Int. J. Plant Sci. 2007, 168, 709–717. [Google Scholar] [CrossRef]
  57. Dombrovska, O.; Qiu, Y.L. Distribution of introns in the mitochondrial gene nad1 in land plants: Phylogenetic and molecular evolutionary implications. Mol. Phylogenet. Evol. 2004, 32, 246–263. [Google Scholar] [CrossRef] [PubMed]
  58. Knoop, V. The mitochondria DNA of land plants: Peculiarities in phylogenetic perspective. Curr. Genet. 2004, 46, 123–139. [Google Scholar] [CrossRef] [PubMed]
  59. Geiss, K.T.; Abbas, G.M.; Makaroff, C.A. Intron loss from the NADH dehydrogenase subunit 4 gene of lettuce mitochondrial DNA: Evidence for homologous recombination of a cDNA intermediate. Mol. Gen. Genet. 1994, 243, 97–105. [Google Scholar] [CrossRef] [PubMed]
  60. Itchoda, N.; Nishizawa, S.; Nagano, H.; Kubo, T.; Mikami, T. The sugar beet mitochondrial nad4 gene: An intron loss and its phylogenetic implication in the Caryophyllales. Theor. Appl. Genet. 2002, 104, 209–213. [Google Scholar] [CrossRef] [PubMed]
  61. Sloan, D.B.; MacQueen, A.H.; Alverson, A.J.; Palmer, J.D.; Taylor, D.R. Extensive loss of RNA editing sites in rapidly evolving Silene mitochondrial genomes: Selection vs. retroprocessing as the driving force. Genetics 2010, 185, 1369–1380. [Google Scholar] [CrossRef] [PubMed]
  62. Zumkeller, S.M.; Knoop, V.; Knie, N. Convergent evolution of fern-specific mitochondrial group II intron atp1i361g2 and its ancient source paralogue rps3i249g2 and independent losses of intron and RNA editing among Pteridaceae. Genome Biol. Evol. 2016, 8, 2505–2519. [Google Scholar] [CrossRef] [PubMed]
  63. Petersen, K.; Schöttler, M.A.; Karcher, D.; Thiele, W.; Bock, R. Elimination of a group II intron from a plastid gene causes a mutant phenotype. Nucleic Acids Res. 2011, 39, 5181–5192. [Google Scholar] [CrossRef] [PubMed]
  64. Nielsen, C.B.; Friedman, B.; Birren, B.; Burge, C.B.; Galagan, J.E. Patterns of intron gain and loss in fungi. PLoS Biol. 2004, 2, e422. [Google Scholar] [CrossRef] [PubMed]
  65. Mitogenome of Tritomaria quinquedentata (GenBank, accession number MG640570).
Genes 08 00395 g001 550
Figure 1. Gene map of the mitogenome of Calypogeia fissa ssp. fissa. Genes inside and outside the outer circle are transcribed in counterclockwise and clockwise directions, respectively. The genes are color-coded based on their function. The inner circle visualizes the G/C content. Pseudogenes have been marked with an asterisk.
Figure 1. Gene map of the mitogenome of Calypogeia fissa ssp. fissa. Genes inside and outside the outer circle are transcribed in counterclockwise and clockwise directions, respectively. The genes are color-coded based on their function. The inner circle visualizes the G/C content. Pseudogenes have been marked with an asterisk.
Genes 08 00395 g001
Genes 08 00395 g002 550
Figure 2. Comparison of the atp1 and cox1 gene structures of different liverwort species. Lost introns are marked with lines.
Figure 2. Comparison of the atp1 and cox1 gene structures of different liverwort species. Lost introns are marked with lines.
Genes 08 00395 g002
Genes 08 00395 g003 550
Figure 3. Occurrence of editing sites in atp1 and cox1 genes of different liverwort species (marked in red).
Figure 3. Occurrence of editing sites in atp1 and cox1 genes of different liverwort species (marked in red).
Genes 08 00395 g003
Table
Table 1. Sample species, voucher specimens and sequencing results of species used in this study. Mean coverage calculated as follows: number of reads that mapped to the mitogenome x average read length (bp)/mitogenome length (bp).
Table 1. Sample species, voucher specimens and sequencing results of species used in this study. Mean coverage calculated as follows: number of reads that mapped to the mitogenome x average read length (bp)/mitogenome length (bp).
SpeciesVoucherSequencing Results (Total Number of Reads)Reads Mapped to the Mitogenome (Number of Reads)Average Read Length (bp)Mitogenome Length (bp)Mitogenome Mean CoverageGenBank Accession Number
Calypogeia argutaUnited Kingdom, collection. DC1420383,522,8542,665,062147.9159,0612478.0MF401630
Calypogeia integristipulaSE Poland, Bieszczady Mts, W slope of Mt Rozsypaniec Wołosacki, 1214 m, collection POZW 4192819,604,00069,73999.6163,05742.6MF401629
Calypogeia fissa ssp. fissaW Poland, Lubuskie Province, Biecz forestry, collection POZW 4230622,601,000123,62099.8162,17576.1MF401632
Calypogeia suecicaSE Poland, Beskid Sądecki Mts, Potok Czarny stream, 717 m, collection POZW 4236612,862,00025,96399.8161,96016.0MF401631
Table
Table 2. Gene contents in liverwort mitochondrial genomes. The black circle or “ψ” indicate the presence of a functional gene or a pseudogene, respectively. The white square indicates a gene lacking. Two black circles indicate the presence of a duplicated copy. a The gene rtl may be pseudogenized [36].
Table 2. Gene contents in liverwort mitochondrial genomes. The black circle or “ψ” indicate the presence of a functional gene or a pseudogene, respectively. The white square indicates a gene lacking. Two black circles indicate the presence of a duplicated copy. a The gene rtl may be pseudogenized [36].
Gene/SpeciesTreubia lacunosaMarchantia polymorphaPleurozia purpureaAneura pinguisCalypogeia fissa ssp. fissaTritomaria quinquedentata
atp1●●●●●●
atp4●●●●●●
atp6●●●●●●
atp8●●●●●●
atp9●●●●●●
ccmBψ●●●●●
ccmC□●●●●●
ccmFCψ●●●●●
ccmFN□●●●●●
cob●●●●●●
cox1●●●●●●
cox2●●●●●●
cox3●●●●●●
nad1●●●●●●
nad2●●●●●●
nad3●●●●●●
nad4●●●●●●
nad4L●●●●●●
nad5●●●●●●
nad6●●●●●●
nad7●ψψψψψ
nad9●●●●●●
rpl2●●●●●●
rpl5●●●●●●
rpl6●●●●●●
rpl10●●●●●●
rpl16●●●●●●
rps1●●●●●●
rps2●●●●●●
rps3●●●●●●
rps4●●●●●●
rps7●●●●●●
rps8●●●●●●
rps10●●●●●●
rps11●●●●●●
rps12●●●●●●
rps13●●●●●●
rps14●●●●●●
rps19●●●●●●
rrn5●●●●●●
rrn18●●●●●●
rrn26●●●●●●
rtl●●● a●●●
sdh3●●●●●●
sdh4●●●●●●
tatC●●●●●●
trnAugc●●●●●●
trnCgca●●●●●●
trnDguc●●●●●●
trnEuuc●●●●●●
trnFgaa●●●●●●
trnGgcc●●●●●●
trnGucc●●●●●●
trnHgug●●●●●●
trnIcau●●●●●●
trnKuuu●●●●●●
trnLcaa●●●●●●
trnLuaa●●●●●●
trnLuag●●●●●●
trnMcau●●●●●●
trnMfcau●●●●●●●●●
trnNguu●●●●●●
trnPugg●●●●●●
trnQuug●●●●●●
trnRacg●●●●●●
trnRucg□●□□□□
trnRucu●●●●●●●●●●●
trnSgcu●●●●●●
trnSuga●●●●●●
trnTggu□●□□□□
trnVuac●●●●●●
trnWcca●●●●●●
trnYgua●●●●●●●●●●●●
ATP synthase subunits: atp1, atp4, atp6, atp8, atp9; Cytochrome c biogenesis subunits: ccmB, ccmC, ccmFC, ccmFN; Cytochrome c reductase subunits: cob; Cytochrome c oxidase subunits: cox1, cox2, cox3; NADH dehydrogenase subunits: nad1, nad2, nad3, nad4, nad4L, nad5, nad6, nad7, nad9; Large ribosomal protein units: rpl2, rpl5, rpl6, rpl10, rpl16; Small ribosomal protein units: rps1, rps2, rps3, rps4, rps7, rps8, rps10, rps11, rps12, rps13, rps14, rps19; Ribosomal RNAs: rrn5, rrn18, rrn26; Reverse transcriptase: rtl; Succinate dehydrogenase subunits: sdh3, sdh4; Twin arginine subunit c: tatC; Transfer RNAs: trnAugc, trnCgca, trnDguc, trnEuuc, trnFgaa, trnGgcc, trnGucc, trnHgug, trnIcau, trnKuuu, trnLcaa, trnLuaa, trnLuag, trnMcau, trnMfcau, trnNguu, trnPugg, trnQuug, trnRacg, trnRucg, trnRucu, trnSgcu, trnSuga, trnTggu, trnVuac, trnWcca, trnYgua.
Table
Table 3. Intron contents in liverwort mitochondrial genomes. The black circle indicates the presence of an intron. The white square indicates an intron lacking. Intron nomenclature follows Dombrovska and Qiu [57] and Knoop [58].
Table 3. Intron contents in liverwort mitochondrial genomes. The black circle indicates the presence of an intron. The white square indicates an intron lacking. Intron nomenclature follows Dombrovska and Qiu [57] and Knoop [58].
Intron/SpeciesTreubia lacunosaMarchantia polymorphaPleurozia purpureaAneura pinguisCalypogeia fissa ssp. fissaTritomaria quinquedentata
atp1i989g2□●●●□□
atp1i1050g2□●●●□□
atp9i87g2●●●●●●
cobi372g2●●●●●●
cobi783g2●●●●●●
cobi824g2●●●●●●
cox1i44g2●●●●●●
cox1i178g2●●●●●●
cox1i375g1●●●●●●
cox1i395g1●●●●□□
cox1i511g2●●●●●●
cox1i624g1●●●●□□
cox1i729g1●●●●□□
cox1i1116g1●●●●●□
cox1i1305g1●●●●●●
cox2i97g2●●●●●●
cox2i250g2●●●●●●
cox3i171g2●●●●●●
cox3i625g2●●●●●●
nad2i709g2●●●●●●
nad3i140g2●●●●●●
nad4i548g2●●●●●●
nad4Li100g2●●●●●●
nad4Li283g2□●●●●●
nad5i753g1●●●●●●
nad7i336g2●●●●●●
nad7i1113g2●●●●●●
rpl2i28g2●●●●●●
rps14i114g2●●●●●●
rrn18i1065g2□●□□□□
rrn26i827g2●●●●●●
trnSgcui43g2●●●●●●

Share and Cite

MDPI and ACS Style

Ślipiko, M.; Myszczyński, K.; Buczkowska-Chmielewska, K.; Bączkiewicz, A.; Szczecińska, M.; Sawicki, J. Comparative Analysis of Four Calypogeia Species Revealed Unexpected Change in Evolutionarily-Stable Liverwort Mitogenomes. Genes 2017, 8, 395. https://0-doi-org.brum.beds.ac.uk/10.3390/genes8120395

AMA Style

Ślipiko M, Myszczyński K, Buczkowska-Chmielewska K, Bączkiewicz A, Szczecińska M, Sawicki J. Comparative Analysis of Four Calypogeia Species Revealed Unexpected Change in Evolutionarily-Stable Liverwort Mitogenomes. Genes. 2017; 8(12):395. https://0-doi-org.brum.beds.ac.uk/10.3390/genes8120395

Chicago/Turabian Style

Ślipiko, Monika, Kamil Myszczyński, Katarzyna Buczkowska-Chmielewska, Alina Bączkiewicz, Monika Szczecińska, and Jakub Sawicki. 2017. "Comparative Analysis of Four Calypogeia Species Revealed Unexpected Change in Evolutionarily-Stable Liverwort Mitogenomes" Genes 8, no. 12: 395. https://0-doi-org.brum.beds.ac.uk/10.3390/genes8120395

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop