Short Communication
Print
Short Communication
Complete mitochondrial genome of the Galápagos sea lion, Zalophus wollebaeki (Carnivora, Otariidae): paratype specimen confirms separate species status
expand article infoRita M. Austin§, Pia Merete Eriksen, Lutz Bachmann
‡ University of Oslo, Oslo, Norway
§ National Museum of Natural History, Smithsonian Institution, Washington D.C., United States of America
Open Access

Abstract

The endangered Galápagos sea lion (Zalophus wollebaeki) inhabits the Galápagos Islands off the coast of Ecuador. We present a complete mitochondrial genome (16 465 bp) of a female paratype from the collections of the Natural History Museum Oslo, Norway, assembled from next-generation sequencing reads. It contains all canonical protein-coding, rRNA, tRNA genes, and the D-loop region. Sequence similarity is 99.93% to a previously published conspecific mitogenome sequence and 99.37% to the mitogenome sequence of the sister species Z. californianus. Sequence similarity of the D-loop region of the Z. wollebaeki paratype mitogenome is >99%, while the sequence difference to the Z. californianus sequences exceeds 2.5%. The paratype mitogenome sequence supports the taxonomic status of Z. wollebaeki as a separate species.

Key words

integrative taxonomy, mitogenome, museomics, type sequencing, unique species

Introduction

Intrinsic to the fields of taxonomy and systematics is diagnosability. By convention, the name and description of any new species is unambiguously linked to the original name-bearing type specimen (types). Species identification of organisms is dependent on these specimens’ morphological and molecular attributes. While DNA sequencing and identification is increasingly used in modern biodiversity research, new challenges regarding taxonomic diagnosability have arisen because molecular data is infrequently available for type specimens. DNA sequencing of name-bearing type specimens is, thus, of particular importance for taxonomy as it enables explicit assignment of extant populations to known types (e.g., Renner 2016). The increasing number of reports clarifying evolutionary taxonomic and systematics issues exemplifies the importance and utility of type specimen sequencing (e.g., Kehlmaier et al. 2019; de Abreu-Jr et al. 2020; Zhang et al. 2020).

Here, we present a complete mitochondrial genome of the Galápagos sea lion (Zalophus wollebaeki Sivertsen, 1953) female paratype, collected in 1925 by Alf Wollebæk from Floreana Island in the Galápagos Archipelago off the coast of Ecuador. Primarily found in the Galápagos Islands archipelago, Z. wollebaeki often congregates in small groups on Isla de la Plata off mainland Ecuador (Melin et al. 2018). Considered Endangered by the International Union for Conservation of Nature (IUCN) (Trillmich 2015), the species has been negatively affected by anthropogenic pressures such as introduced diseases and the increasing frequency of El Niño events linked to climate change (Melin et al. 2018; Krüger et al. 2021; Páez-Rosas et al. 2021). In 2018, a limited population size of 17 000–24 000 individuals was reported after a decline of 23.8% in 2015, a year with an El Niño event (Páez-Rosas et al. 2021). As the number and intensity of El Niño events increases, Z. wollebaeki is also facing decreasing pup abundance, further impacting the species’ survivability and conservation status (Páez-Rosas et al. 2021).

The taxonomy of the genus Zalophus, including the three species Z. japonicus (Peters, 1866) Z. californianus (Lesson, 1828, as cited in Ellerman and Morrison-Scott 1966) and Z. wollebaeki, has been contentious. Zalophus wollebaeki was initially recognized as a new species (Sivertsen 1953), morphologically, but was later reconsidered to be a subspecies of Z. californianus (Scheffer 1958). More recently accumulated molecular evidence has suggested that Z. wollebaeki is a separate species (Wolf et al. 2007, 2008; Schramm et al. 2009; Berta and Churchill 2012; Krüger et al. 2021). When comparing mitochondrial D-loop and cytochrome B sequences, Z. californianus and Z. wollebaeki were found to be reciprocally monophyletic, and 25 microsatellite loci further revealed numerous private alleles (Wolf et al. 2007). Although, a recent study (Hassanin et al. 2021) reported an uncorrected pairwise distance between Z. californianus and Z. wollebaeki mitogenomes of only 0.5%, suggesting that they are, in fact, the same species; despite this low percent difference, they ultimately upheld the separate species taxonomy and instead suggested an adjusted divergence time between the two species to be 0.3–0.2 million years (Hassanin et al. 2021). Originally considered morphologically distinct (Sivertsen 1953), molecular assessment of a Z. wollebaeki type specimen can therefore provide a direct link to the initial taxonomic description, while also clarifying evolutionary relationships within Zalophus.

Material and methods

A female paratype, collected in 1925 by Alf Wollebæk from Floreana Island (Natural History Museum of Oslo, Norway, voucher number NHMO-30317) was used for this study. Interestingly, the given type series collection location differs among the original catalog record (“Chatham/San Cristóbol”), Sivertsen's (1953) holotype description (“Floreana/Sancta Maria”), and the Wiig and Bachmann (2013) publication (“Isla San Cristóbol”) for this specimen. Here, we specify Floreana/Santa Maria Island in corroboration with other specimen collection locations and dates from the same expedition (Wollebæk 1934) that specify “Santa Maria”.

Total genomic DNA was extracted from a left front flipper skin biopsy (257 mg) using the QIAamp DNA Micro Kit (Qiagen, Germany) according to the manufacturer’s instructions. Remaining tissue and DNA extract are stored in the scientific collections of the Natural History Museum, University of Oslo, Norway (voucher number NHMO-DMA-30317/6-D). Extracted DNA (2.4 μg) was submitted for custom sequencing (Illumina NovaSeq 2×150 bp) at the Norwegian Sequencing Centre (https://www.sequencing.uio.no).

The obtained 26 164 466 raw reads (SRA number PRJNA805083) were adapter-trimmed and quality filtered using AdapterRemoval2 (Schubert et al. 2016), with a total of 23 297 925 trimmed and merged reads de novo assembled with SPAdes v. 3.13.1 (Bankevich et al. 2012). An initial mitochondrial sequence assembly using MITObim v. 1.9.1 (Hahn et al. 2013), with a Z. wollebaeki D-loop sequence (GenBank accession number AM422173.1) as a bait, yielded a 6042 bp contig. The final mitogenome sequence was obtained by blasting the MITObim sequence against the 257 673 SPAdes-generated scaffolds >500 bp in length. The final mitochondrial genome sequence (GenBank accession number OM636180) was annotated using MITOS2 (Donath et al. 2019) alongside other published Otariidae mitogenomes.

Excluding the D-loop region, the Z. wollebaeki paratype mitogenome was aligned to 13 other otariid species [including the previously sequenced Z. wollebaeki specimen (Hassanin et al. 2021)] and three phocid outgroup species using MUSCLE on the EMBL-EBI server (Kanz et al. 2005) (https://www.ebi.ac.uk/Tools/msa/muscle/).

The maximum likelihood analysis was run on the ATGC Montpellier Bioinformatics platform (http://www.atgc-montpellier.fr) using PhyML (Guindon et al. 2010) under the GTR model and the Akaike information criterion.

Results and discussion

The mitogenome of the Z. wollebaeki paratype was assembled with an average coverage of 36.2X. It is 16 465 bp long and includes all canonical protein-coding sequences, rRNAs, tRNAs, and the D-loop region. Sequence similarity to a previously reported Z. wollebaeki (SRR4431565) mitogenome (Hassanin et al. 2021) was 99.93% (12 nucleotide substitutions) and 99.37% to the sister species Z. californianus. Accordingly, a maximum likelihood analysis placed the Z. wollebaeki paratype sequence in a clade consisting of Z. wollebaeki, Z. californianus, and Z. japonicus (Fig. 1), confirming the authenticity of the paratype sequence.

Figure 1. 

Maximum likelihood tree of several otariid species. Relationships are depicted among the Z. wollebaeki paratype mitogenome (bolded and blue), 13 otariid species (incl. the previously described Z. wollebaeki mitogenome sequence, GenBank accession number AM422173.1), and three phocid outgroup species. The scale indicates genetic distance. Nodes with 100% bootstrap support (1000 replicates) are depicted in green, with the remaining nodes’ respective bootstrap values indicated to the left.

Comparison with previously published D-loop sequences of Z. wollebaeki (Wolf et al. 2007, 2008) revealed a new haplotype that is differentiated from the other 36 haplotypes with one C→T transition (position 15 610). Authenticity of the substitution was confirmed with all 37 reads covering the region sharing the T nucleotide, indicating that erroneous incorporation due to DNA damage can be excluded. The new haplotype is most similar to haplotype Zw_H_10, one of the most common haplotypes reported for Z. wollebaeki (Wolf et al. 2008). The sequence similarity of the D-loop region of the Z. wollebaeki paratype mitogenome is >99%, whereas sequence differences exceed 2.5% when compared to the 12 Z. californianus sequences >600 bp that are listed in GenBank (Wolf et al. 2007). This indicates the Z. wollebaeki paratype mitogenome sequence matches well with all nucleotide sequence data reported for more recently collected samples of the species.

The final assembly of the Z. wollebaeki paratype mitogenome sequence included 5152 raw reads, which is 0.022% of the adapter-trimmed and quality filtered Illumina readpool. This seemingly low proportion is in the same order of magnitude as observed in other studies. For example, for total DNA extracted from mouse embryonic fibroblasts Quispe-Tintaya et al. (2013) reported a proportion of 0.1% mitochondrial DNA reads, while Anmarkrud and Lifjeld (2017) found 0.03–0.18% mapping mitochondrial DNA reads for historical DNA extracted from museum specimens of extinct birds. The estimated proportion of 0.022% mitochondrial reads is much less than the frequently used rule of thumb stating that mitochondrial DNA represents ~1% of a total DNA extract of mammals. However, the proportion of mitochondrial reads in a next-generation sequencing readpool depends on many parameters such as the extraction protocol, the sequencing methodology, and DNA quality or contamination from other sources. Museum samples are frequently contaminated with external DNA of human origin. Therefore, a control mapping of the adapter-trimmed and quality filtered Illumina readpool to a human mitochondrial genome was conducted. The very low number of 1031 mapping reads indicated that contamination with human DNA is a minor issue, however, this does not mean that the overall sample contamination from external sources is low.

Overall, the complete mitochondrial genome sequence and newly identified haplotype represent valuable genetic references in support of a species distinction between Z. californianus and Z. wollebaeki using a museum paratype specimen, which may also be constructive for conservation efforts geared toward this charismatic and unique species, and its habitat. With evolutionary relationships within Zalophus clarified, genetic assignment of extant populations can now be made more accurately and readily. Furthermore, having an unambiguous connection between the species’ genetic information and the original taxonomic description of Z. wollebaeki and paratype specimen, fulfills recent recommendations for incorporating DNA-based species identifications and diagnoses into the various Codes of Nomenclature.

Acknowledgements

We would like to thank the Natural History Museum of Oslo’s collections manager Lars Erik Johannessen and Mammals curator Kjetil Lysne Voje for granting sample access. We are also thankful to the referees and editor who contributed their time and valuable insight to improving this manuscript.

Additional information

Conflict of interest

No conflict of interest was declared.

Ethical statement

No approval of research ethics committees was required to accomplish the goals of this study as no live animals were considered in this study using, instead, historic museum collection material.

Funding

Funding for this project was provided by the Natural History Museum, University of Oslo. Access to Saga, the high-performance computer at the Norwegian Research Infrastructure Services (NRIS), was granted through the Norwegian Metacenter for Computational Science (NOTUR; project nn9201k).

Author contributions

RMA and LB formulated the original research goals and concept; RMA, PME, and LB contributed equally to the methodologies; RMA gathered samples and conducted all laboratory work; RMA and PME conducted formal analyses, visualization, and synthesis. All authors contributed to writing the original draft, editing, and approving the manuscript.

Author ORCIDs

Rita M. Austin https://orcid.org/0000-0002-9315-517X

Pia Merete Eriksen https://orcid.org/0000-0003-1832-4794

Lutz Bachmann https://orcid.org/0000-0001-7451-2074

Data availability

The complete mitochondrial genome and sequencing reads that support this study are openly available on GenBank of NCBI under the accession no. OM636180). The associated BioProject, SRA, and Bio-Sample numbers are PRJNA805083, SRP359162, and SAMN25827133, respectively. The complete mitogenome sequence of Z. wollebaeki reported by Hassanin et al. (2021) was obtained from the authors.

References

  • Anmarkrud JA, Lifjeld JT (2017) Complete mitochondrial genomes of eleven extinct or possibly extinct bird species. Molecular Ecology Resources 17(2): 334–341. https://doi.org/10.1111/1755-0998.12600
  • Bankevich A, Nurk S, Antipov D, Gurevich AA, Dvorkin M, Kulikov AS, Lesin VM, Nikolenko SI, Prjibelski AD, Pyshkin AV, Sirotkin A, Vyahhi N, Tesler G, Alekseyev MA, Pevzner PA (2012) SPAdes: A new genome assembly algorithm and its applications to single-cell sequencing. Journal of Computational Biology 19(5): 455–477. https://doi.org/10.1089/cmb.2012.0021
  • de Abreu-Jr EF, Pavan SE, Tsuchiya MN, Wilson DE, Percequillo AR, Maldonado JE (2020) Museomics of tree squirrels: A dense taxon sampling of mitogenomes reveals hidden diversity, phenotypic convergence, and the need of a taxonomic overhaul. BMC Evolutionary Biology 20(1): 77. https://doi.org/10.1186/s12862-020-01639-y
  • Donath A, Jühling F, Al-Arab M, Bernhart SH, Reinhardt F, Stadler PF, Middendorf M, Bernt M (2019) Improved annotation of protein-coding genes boundaries in metazoan mitochondrial genomes. Nucleic Acids Research 47(20): 10543–10552. https://doi.org/10.1093/nar/gkz833
  • Guindon S, Dufayard JF, Lefort V, Anisimova M, Hordijk W, Gascuel O (2010) New algorithms and methods to estimate maximum-likelihood phylogenies: Assessing the performance of PhyML 3.0. Systematic Biology 59(3): 307–321. https://doi.org/10.1093/sysbio/syq010
  • Hahn C, Bachmann L, Chevreux B (2013) Reconstructing mitochondrial genomes directly from genomic next-generation sequencing reads—A baiting and iterative mapping approach. Nucleic Acids Research 41(13): e129–e129. https://doi.org/10.1093/nar/gkt371
  • Hassanin A, Veron G, Ropiquet A, Jansen van Vuuren B, Lécu A, Goodman SM, Haider J, Nguyen TT (2021) Evolutionary history of Carnivora (Mammalia, Laurasiatheria) inferred from mitochondrial genomes. PLoS ONE 16(2): e0240770. https://doi.org/10.1371/journal.pone.0240770
  • Kanz C, Aldebert P, Althorpe N, Baker W, Baldwin A, Bates K, Browne P, van den Broek A, Castro M, Cochrane G, Duggan K, Eberhardt R, Faruque N, Gamble J, Garcia Diez F, Harte N, Kulikova T, Lin Q, Lombard V, Lopez R, Mancuso R, McHale M, Nardone F, Silventoinen V, Sobhany S, Stoehr P, Tuli MA, Tzouvara K, Vaughan R, Wu D, Zhu W, Apweiler R (2005) The EMBL nucleotide sequence database. Nucleic Acids Research 33(Database issue): D29–D33. https://doi.org/10.1093/nar/gki098
  • Kehlmaier C, Zhang X, Georges A, Campbell PD, Thomson S, Fritz U (2019) Mitogenomics of historical type specimens of Australasian turtles: Clarification of taxonomic confusion and old mitochondrial introgression. Scientific Reports 9(1): 5841. https://doi.org/10.1038/s41598-019-42310-x
  • Krüger O, Kalberer S, Meise K, Schwarz J, DeRango E, Piedrahita P, Páez-Rosas D, Trillmich F (2021) Galápagos Sea Lion (Zalophus wollebaeki, Sivertsen 1953). In: Heckel G, Schramm Y (Eds) Ecology and Conservation of Pinnipeds in Latin America. Springer Nature, 145–163. https://doi.org/10.1007/978-3-030-63177-2
  • Melin SR, Trillmich F, Aurioles-Gamboa D (2018) California, Galapagos, and Japanese Sea Lions: Zalophus californianus, Z. wollebaeki, and Z. japonicus. In: Wursig B, Thewissen JGM, Kovacs K (Eds) Encyclopedia of Marine Mammals. Elsevier, 153–157. https://doi.org/10.1016/B978-0-12-804327-1.00003-0
  • Páez-Rosas D, Torre J, Espinoza E, Marchetti A, Seim H, Riofrio-Lazo M (2021) Declines and recovery in endangered Galapagos pinnipeds during the El Niño event. Scientific Reports 11(1): 8785. https://doi.org/10.1038/s41598-021-88350-0
  • Peters W (1866) Einen Nachtrag zu Seiner abhandlung uber die ohrenrobben, otariae (Monatsbericht 1866 pp. 261). Monatsberichte der koniglich preussischen akademie der wissenschaften zu berlin 1866: 665–672.
  • Quispe-Tintaya W, White RR, Popov VN, Vijg J, Maslov AY (2013) Fast mitochondrial DNA isolation from mammalian cells for next-generation sequencing. BioTechniques 55(3): 133–136. https://doi.org/10.2144/000114077
  • Renner SS (2016) A return to Linnaeus’s focus on diagnosis, not description: The use of DNA characters in the formal naming of species. Systematic Biology 65(6): 1085–1095. https://doi.org/10.1093/sysbio/syw032
  • Scheffer VB (1958) Seals, Sea Lions and Walruses: A Review of the Pinnipedia. Stanford University Press.
  • Schramm Y, Mesnick S, De la Rosa J, Palacios D, Lowry M, Aurioles-Gamboa D, Snell H, Escorza-Trevino S (2009) Phylogeography of California and Galápagos sea lions and population structure within the California sea lion. Marine Biology 156(7): 1375–1387. https://doi.org/10.1007/s00227-009-1178-1
  • Sivertsen E (1953) A new species of sea lion, Zalophus wollebaeki, from the Galapagos Islands. Det Kongelige Norske Videnskabers Selskabs Forhandlinger 26: 1–3.
  • Wolf JBW, Tautz D, Trillmich F (2007) Galápagos and Californian sea lions are separate species: Genetic analysis of the genus Zalophus and its implications for conservation management. Frontiers in Zoology 4(1): 1–13. https://doi.org/10.1186/1742-9994-4-20
  • Wolf JBW, Brunner S, Salazar S, Trillmich F, Tautz D (2008) Tracing early stages of species differentiation: Ecological, morphological and genetic divergence of Galápagos sea lion populations. BMC Evolutionary Biology 8(1): 150. https://doi.org/10.1186/1471-2148-8-150
  • Wollebæk A (1934) De forheksede øer. J. W. Cappelens Forlag, Oslo.
  • Zhang J, Lees DC, Shen J, Cong Q, Huertas B, Martin G, Grishin NV (2020) The mitogenome of a Malagasy butterfly Malaza fastuosus (Mabille, 1884) recovered from the holotype collected over 140 years ago adds support for a new subfamily of Hesperiidae (Lepidoptera). Genome 63(4): 195–202. https://doi.org/10.1139/gen-2019-0189
login to comment