712
Views
0
CrossRef citations to date
0
Altmetric
Mitogenome Announcement

The complete mitogenome of the inarticulate brachiopod Glottidia pyramidata reveals insights into gene order variation, deviant ATP8 and mtORFans in the Brachiopoda

, &
Pages 2701-2703 | Received 14 Jun 2021, Accepted 04 Aug 2021, Published online: 18 Aug 2021

Abstract

Brachiopods are a clade of marine organisms with a tremendously diverse and abundant fossil record but with fewer than 500 species extant today. Even if a better understanding of their biology and genetics could help to test hypotheses about their impressive decline, knowledge of genetics and evolutionary genomics in extant brachiopods is very poor. Here, we present the complete mitochondrial genome sequence of the inarticulate Glottidia pyramidata, an eastern North American extant representative of the phylum Brachiopoda. Besides the general characteristics of the sequenced mitogenome, we present its unusual features such as deviant ATP8 protein sequence and supernumerary ORFs, and also unique gene order, considering the available genome sequences of other brachiopod species.

Also known as lamp shells or the ‘other’ bivalves, brachiopods are a phylum of marine invertebrates with a rich paleontological history, with more than 30,000 extinct species known to science (Thayer Citation1985; Carlson Citation2016). Even though they resemble bivalved mollusks, brachiopods are in fact bivalved lophophorates. Unlike their diverse and abundant ancestral counterparts, fewer than 500 species are alive today and they live largely hidden, at great depths and locations with low diversity and abundance (Carlson Citation2016). Competition from ‘mussel-like bivalves’ has been proposed as a plausible cause of the spectacular decline of this animal group in the end-Permian mass extinction (Thayer Citation1985). A better understanding of the biology and genetics of the few extant species may help to test hypotheses about this impressive decline, but knowledge of genetics and evolutionary genomics in extant brachiopods is still relatively poor (Carlson Citation2016), including at the mitochondrial level.

Interestingly, some of the few brachiopod mitochondrial DNAs (mtDNAs) sequenced to date resemble those of bivalves with DUI (i.e. with doubly uniparental inheritance of mtDNA characterized by the presence of two sex-associated mtDNA lineages inherited through males and females, respectively; Breton et al. Citation2007; Stewart et al. Citation2020). As for DUI bivalves (e.g. Breton et al. Citation2010, Citation2011a, Citation2011b, Citation2014; Milani et al. Citation2013; Lubośny et al. Citation2018; Guerra et al. Citation2019; Capt et al. Citation2020), brachiopods present important variation in mtDNA size (17 to >28 kb) and gene arrangement, with longer and deviant protein-coding sequences compared to their homologues in other animals (e.g. longer cox2 gene and deviant ATP8 protein sequences) as well as supernumerary open reading frames (ORFs) or mtORFans (ORFs without recognizable homologies to other known genes) (Helfenbein et al. Citation2001; Endo et al. Citation2005; Luo et al. Citation2015; Karagozlu et al. Citation2017a). There are seven complete brachiopod mitogenomes available to date, three of them are from Lingula anatina (order Lingulida), and the four others are from species from the order Terebratulida, which contains >75% of all extant brachiopod species (Emig et al. Citation2013). Three other orders exist in this phylum: Craniida, Rhynchonellida, Thecideida (Emig et al. Citation2013). In the traditional two-class system, Terebratulida + Rhynchonellida + Thecideida constitute the class Articulata (brachiopods with two valves connected by a tooth and socket hinge) and Lingulida + Craniida the class Inarticulata (unhinged valves connected by muscles alone) (Carlson Citation2016). It is mainly the inarticulate L. anatina that presents unusual mitochondrial features, such as intraspecific variation in size and arrangement, and mtORFans (Endo et al. Citation2005; Luo et al. Citation2015; Karagozlu et al. Citation2017a). This raises the question if some of these features are also shared by other inarticulate species, thus highlighting the need for more complete mtDNAs to expand our knowledge on brachiopod mitogenomics.

Here, we present the complete mitogenome of the inarticulate brachiopod Glottidia pyramidata Stimpson 1860 (Brachiopoda: Inarticulata). Four adult specimens were collected by Gulf Specimen Marine Laboratories and Aquarium in June 2016 in the Dickerson Bay (USA; near N30.023633 W84.385280) and sent alive to the Université de Montréal. Each individual was dissected and sexed by inspecting the gonads under a light microscope (×100) for the presence of eggs or sperm, three females and one male were unambiguously sexed and total DNA was extracted from one male gonad, sequenced, assembled, and annotated following Guerra et al. (Citation2018). The specimen and DNA sample were deposited in the Department of Biological Sciences Tissue & DNA Collection at the Université de Montréal (https://bio.umontreal.ca/english/home/; specimen code B1408 and DNA sample code B1408g; person in charge Breton S, [email protected]). Since brachiopods present characteristics that resemble those of DUI bivalves, we first looked for the presence of DUI in G. pyramidata, but our results revealed only one set of mitochondrial contigs (consistent with the existence of only one mitochondrial lineage), suggesting that it might possess a strictly maternal mitochondrial transmission like other metazoans (Breton and Stewart Citation2015). The mitogenome sequence has been deposited in GenBank (accession number MW732171).

Figure 1. Maximum-likelihood phylogeny of the inarticulate Glottidia pyramidata in the phylum Brachiopoda. The tree was obtained with MEGAX and is based on concatenated sequences of 12 protein-coding genes (atp8 was excluded) from eight brachiopod species and two outgroups chosen from the phylum Bryozoa. GenBank accession numbers are presented after species names. The number below or above the branches indicate bootstrap support values above 95%. All known gene orders of Brachiopoda were mapped on the obtained phylogeny.

Figure 1. Maximum-likelihood phylogeny of the inarticulate Glottidia pyramidata in the phylum Brachiopoda. The tree was obtained with MEGAX and is based on concatenated sequences of 12 protein-coding genes (atp8 was excluded) from eight brachiopod species and two outgroups chosen from the phylum Bryozoa. GenBank accession numbers are presented after species names. The number below or above the branches indicate bootstrap support values above 95%. All known gene orders of Brachiopoda were mapped on the obtained phylogeny.

The complete mtDNA of G. pyramidata is 19,709 bp long, which is larger than the mtDNAs of the articulate brachiopods Terebratulina retusa (15,451 bp; Stechmann and Schlegel Citation1999), Terebratalia transversa (14,291 bp; Helfenbein et al. Citation2001), Laqueus rubellus (14,017 bp; Noguchi et al. Citation2000), and Laqueus japonicus (14,267 bp; Karagozlu et al. Citation2017b), but in the range of the sizes observed for the inarticulate L. anatina (17,970 bp, Luo et al. Citation2015; 25,790 bp, Karagozlu et al. Citation2017a, and 28,818 bp, Endo et al. Citation2005). It contains 36 genes, which are all located on the same strand as in other brachiopod species. The annotation of atp8 is unsure, otherwise all other typical mitochondrial genes, i.e. 12 coding for proteins, 22 for tRNAs, and two for rRNAs, were retrieved and annotated without any ambiguities. The atp8 gene has potentially been localized in an unassigned region preceding atp6 (as for the inarticulate L. anatina), and it is shorter (43aa) than what has been found in other brachiopods (51aa to 59aa). It does not possess the typical MPQL amino acid signature present at the N-terminus of metazoan ATP8 (only L. anatina also presents this feature), but it is characterized by a hydrophobic N-terminus domain like typical ATP8 proteins (Breton et al. Citation2010). Moreover, as L. anatina, the mitogenome of G. pyramidata contains mtORFans, which are found between the genes nad4L and nad6 (i.e. three mtORFans of 200, 99, and 64 amino acids).

For our phylogenetic analysis, all brachiopod mitogenomes publicly available, as well as two outgroup mitogenomes (the bryozoans Bugula neritina [AY690838] and Celleporella hyalina [JQ839275]) were retrieved from GenBank (01-07-2020) and their 12 protein-coding genes (atp8 was excluded) were translated with the invertebrate mitochondrial genetic code and individually aligned using MEGAX (Kumar et al. Citation2018) and the MUSCLE algorithm (Edgar Citation2004). These alignments were afterwards concatenated with MEGAX and used to construct a maximum-likelihood (ML) tree using MEGAX and the JTT model for amino acid substitution. We mapped all known brachiopod gene orders on the phylogeny. Our results show that the inarticulate G. pyramidata clustered together with the inarticulate L. anatina with high support value, and that the articulate brachiopods cluster together (). The gene order in G. pyramidata is unique among brachiopod species sequenced to date, adding to the gene order variability previously reported in this group (Luo et al. Citation2015).

Disclosure statement

The authors report no conflict of interest.

Data availability statement

The genome sequence data that support the findings of this study are openly available in GenBank of NCBI at https://www.ncbi.nlm.nih.gov under the accession no. MW732171. The associated BioProject, SRA, and Bio-Sample numbers are PRJNA737150, SRP323932, and SAMN19679991, respectively.

Additional information

Funding

This work was supported by the National Sciences and Engineering Research Council (NSERC) under the Grant Number RGPIN-2019-04076 to S.B. S.B. also holds the Canada Research Chair in Evolutionary Mitochondrial Biology.

Notes on contributors

Thierry Niaison

Thierry Niaison is currently a BSc candidate at the Université de Montréal. His fields of interest include evolutionary mitochondrial biology and physiology.

Davide Guerra

Davide Guerra was a postdoctoral fellow in S Breton’s lab from 2015 until 2018. He is currently working for a governmental agency in Italy.

Sophie Breton

Sophie Breton’s research program is integrative, using approaches including molecular biology, physiology, genomics, transcriptomics, proteomics, and evolutionary biology to better understand the role of mitochondria and their genomes in speciation. In this perspective, her research aims at better understand the function, structure, expression, evolution, and inheritance of mitochondrial genomes, the co-adaptation of mitochondrial and nuclear genomes and the physiological mechanisms of animal adaptation to environmental conditions.

References

  • Breton S, Beaupré HD, Stewart DT, Hoeh WR, Blier PU. 2007. The unusual system of doubly uniparental inheritance of mtDNA: isn't one enough? Trends Genet. 23(9):465–474.
  • Breton S, Ghiselli F, Passamonti M, Milani L, Stewart DT, Hoeh WR. 2011a. Evidence for a fourteenth mtDNA-encoded protein in the female-transmitted mtDNA of marine mussels (Bivalvia: Mytilidae). PLOS One. 6(4):e19365.
  • Breton S, Milani L, Ghiselli F, Guerra D, Stewart DT, Passamonti M. 2014. A resourceful genome: updating the functional repertoire and evolutionary role of animal mitochondrial DNAs. Trends Genet. 30(12):555–564.
  • Breton S, Stewart DT. 2015. Atypical mitochondrial inheritance patterns in eukaryotes. Genome. 58(10):423–431.
  • Breton S, Stewart DT, Hoeh WR. 2010. Characterization of a mitochondrial ORF from the gender-associated mtDNAs of Mytilus spp. (Bivalvia: Mytilidae): identification of the “missing” ATPase 8 gene. Mar Genomics. 3(1):11–18.
  • Breton S, Stewart DT, Shepardson S, Trdan RJ, Bogan AE, Chapman EG, Ruminas AJ, Piontkivska H, Hoeh WR. 2011b. Novel protein genes in animal mtDNA: a new sex determination system in freshwater mussels (Bivalvia: Unionoida)? Mol Biol Evol. 28(5):1645–1659.
  • Capt C, Bouvet K, Guerra D, Robicheau BM, Stewart DT, Pante E, Breton S. 2020. Unorthodox features in two venerid bivalves with doubly uniparental inheritance of mitochondria. Sci Rep. 10(1):1087.
  • Carlson SJ. 2016. The evolution of Brachiopoda. Annu Rev Earth Planet Sci. 44(1):409–438.
  • Edgar RC. 2004. MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 32(5):1792–1797.
  • Emig CC, Bitner MA, Álvarez F. 2013. Phylum Brachiopoda. In: Zhang Z-Q, editor. Animal biodiversity: an outline of higher-level classification and survey of taxonomic richness. Auckland: Magnolia Press; p. 75–78.
  • Endo K, Noguchi Y, Ueshima R, Jacobs HT. 2005. Novel repetitive structures, deviant protein-encoding sequences and unidentified ORFs in the mitochondrial genome of the brachiopod Lingula anatina. J Mol Evol. 61(1):36–53.
  • Guerra D, Bouvet K, Breton S. 2018. Mitochondrial gene order evolution in Mollusca: inference of the ancestral state from the mtDNA of Chaetopleura apiculata (Polyplacophora, Chaetopleuridae). Mol Phylogenet Evol. 120:233–239.
  • Guerra D, Lopes-Lima M, Froufe E, Ming Gan H, Paz O, Amaro R, Klunzinger MW, Callil C, Prié V, Bogan AE, et al. 2019. Variability of mitochondrial ORFans hints at possible differences in the system of doubly uniparental inheritance of mitochondria among families of freshwater mussels (Bivalvia: Unionida). BMC Evol Biol. 19(1):229.
  • Helfenbein KG, Brown WM, Boore JL. 2001. The complete mitochondrial genome of the articulate brachiopod Terebratalia transversa. Mol Biol Evol. 18(9):1734–1744.
  • Karagozlu MZ, Kim S-G, Thinh DD, Kim C-B. 2017a. Complete mitochondrial genome analysis of Lingula anatina from Korea (Brachiopoda, Lingulida, Lingulidae). Mitochondrial DNA B Resour. 2(2):829–830.
  • Karagozlu MZ, Kim S-G, Thinh DD, Kim C-B. 2017b. Complete mitochondrial genome of Laqueus japonicus (Brachiopoda, Terebratulida, Laqueidae). Mitochondrial DNA B Resour. 2(2):883–884.
  • Kumar S, Stecher G, Li M, Knyaz C, Tamura K. 2018. MEGA X: molecular evolutionary genetics analysis across computing platforms. Mol Biol Evol. 35(6):1547–1549.
  • Lubośny M, Przyłucka A, Śmietanka B, Breton S, Burzyński A. 2018. Actively transcribed and expressed atp8 gene in Mytilus edulis mussels. PeerJ. 6:e4897.
  • Luo Y-J, Satoh M, Endo K. 2015. Mitochondrial gene order variation in the brachiopod Lingula anatina and its implications for mitochondrial evolution in lophotrochozoans. Mar Genomics. 24 Pt 1:31–40.
  • Milani L, Ghiselli F, Guerra D, Breton S, Passamonti M. 2013. A comparative analysis of mitochondrial ORFans: new clues on their origin and role in species with doubly uniparental inheritance of mitochondria. Genome Biol Evol. 5(7):1408–1434.
  • Noguchi Y, Endo K, Tajima F, Ueshima R. 2000. The mitochondrial genome of the brachiopod Laqueus rubellus. Genetics. 155(1):245–259.
  • Stechmann A, Schlegel M. 1999. Analysis of the complete mitochondrial DNA sequence of the brachiopod Terebratulina retusa places Brachiopoda within the protostomes. Proc Biol Sci. 266(1433):2043–2052.
  • Stewart DT, Breton S, Chase EE, Robicheau BM, Bettinazzi S, Pante E, Youssef N, Garrido-Ramos MA. 2020. An unusual evolutionary strategy: the origins, genetic repertoire, and implications of doubly uniparental inheritance of mitochondrial DNA in bivalves. In: Pontarotti P, editor. Evolutionary biology from concept to application. Berlin: Springer-Verlag; p. 301–323.
  • Thayer CW. 1985. Brachiopods vs. mussels: competition, predation, and palatability. Science. 228(4707):1527–1528.