Next Article in Journal
A Comparative Transcriptome Analysis of Volvariella volvacea Identified the Candidate Genes Involved in Fast Growth at the Mycelial Growth Stage
Previous Article in Journal
A Missense Variant Affecting the C-Terminal Tail of UNC93B1 in Dogs with Exfoliative Cutaneous Lupus Erythematosus (ECLE)
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Genome Assembly and Annotation of the Southern Elephant Seal Mirounga leonina

1
Unit of Research for Practical Application, Korea Polar Research Institute, Incheon 21990, Korea
2
Division of Biotechnology, College of Life Sciences and Biotechnology, Korea University, Seoul 02841; Korea
3
Division of Polar Life Science, Korea Polar Research Institute, Incheon 21990, Korea
4
Polar Sciences, University of Science & Technology, Yuseong-gu, Daejeon 34113, Korea
*
Author to whom correspondence should be addressed.
Equal contribution.
Genes 2020, 11(2), 160; https://doi.org/10.3390/genes11020160
Submission received: 13 January 2020 / Revised: 31 January 2020 / Accepted: 31 January 2020 / Published: 3 February 2020
(This article belongs to the Section Animal Genetics and Genomics)

Abstract

:
The southern elephant seal Mirounga leonina is the largest phocid seal and one of the two species of elephant seals. They are listed as ‘least concern’ by the International Union for Conservation of Nature (IUCN) Red List of Threatened Species 2015. Here, we have assembled the reference genome for M. leonina using the 10× chromium sequencing platform. The final genome assembly of M. leonina was 2.42 Gb long, with a contig N50 length of 54 Mb and a maximum length of 111.6 Mb. The M. leonina genome contained 20,457 predicted protein-coding genes and possessed 41.51% repeated sequences. The completeness of the M. leonina genome was evaluated using benchmarking universal single-copy orthologous genes (BUSCOs): the assembly was highly complete, containing 95.6% of the core set of mammalian genes. The high-quality genomic information on M. leonina will be essential for further understanding of adaptive metabolism upon repeated breath-hold dives and the exploration of molecular mechanisms contributing to its unique biochemical and physiological characteristics. The southern elephant seal genome project was deposited at NCBI (National Center for Biotechnology Information) under BioProject number PRJNA587380.

1. Introduction

Elephant seals have been highlighted as a crucial animal model for studying their unique behavior, physiology, population dynamics, and geographical distribution. Unlike other ocean-going mammals (e.g., whales, dugongs, manatees), elephant seals have to emerge from the water to rest, molt, mate, and rear pups [1]. Seals undergo a long fasting period, and the energy generated from the blubber-fat allows them stay on land for over a month without food and water. The northern elephant seal is known for losing up to 40% body mass during the prolonged fasting period [2]. Furthermore, seals have evolved to regulate osmotic challenge by producing endogenous water from lipid oxidation and concentration of urine [3].
The elephant seal genus Mirounga contains two species, the southern elephant seal Mirounga leonina and the northern elephant seal M. angustirostris. The southern elephant seal is the largest species in the clade Pinnipedia, as adult males can weigh up to four tons and grow up to 5.8 m in length. Four distinct populations were investigated in the circumpolar regions and sub-Antarctic islands [4,5]. Southern elephant seals have evolved to be able to deep- and long- dive for foraging. During predation, they spend most of their time underwater to prey on squid, mollusks, krill, and cephalopods, and they can undertake lengthy dives from 30~45 min to over 2 h and dive more than 2000 m in depth [6]. The enlarged eyes and sensitive vibrissae whiskers play a crucial role in predation in the deep sea [6].
M. leonina was almost exterminated by indiscriminate slaughter to harvest oil from the blubber in the early 19th century. Thus, they were listed as ‘vulnerable’ under the predecessor to the Environment Protection and Biodiversity Conservation Act 1999 (EPBC Act). Hunting was regulated by the International Convention for the Conservation of Antarctic seals and the Convention on Antarctic Marine Living Resources to conserve M. leonina individuals and their population, and population size has been recovered. They are currently listed as ‘least concern’ by the International Union for Conservation of Nature (IUCN) Red List of Threatened Species 2015 (e.T13583A45227247, http://www.iucnredlist.org, accessed on: 12. Dec. 2014). The southern elephant seal has encountered several threats, such as the depletion of prey stocks owing to intensive fishing in Antarctica, drastic fluctuations in habitat by climate change, and low genetic diversity [7]. Therefore, the consistent monitoring and highly regulated conservation of M. leonina populations are strongly required.
Genome information of M. leonina could help to understand the unique physiological and metabolic adaptations, such as the ability to spend large amounts of time in the sea as a mammal, survive in a frigid environment, perform repeated breath-hold dives, engage in deep-diving and long-ranging predation, tolerate routine hypoxia, and fast during the nursing period. In this study, the genome of the southern elephant seal was sequenced and analyzed for the first time in elephant seals. The genome of M. leonina was effectively sequenced from a small blood sample using the 10× chromium sequencing platform to diminish potential stress during sampling. This study will contribute to preserving the species by understanding its physiology and molecular mechanisms.

2. Materials and Methods

The genomic DNA of the southern elephant seal, Mirounga leonina, (Figure 1) was extracted from the blood specimen as described in our previous study [8]; approximately 1 mL of blood sample was collected from a single elephant seal on King George Island, South Shetland Islands, Antarctica. The blood sample (200 µL) was used to obtain high molecular weight gDNA using the QIAGEN MagAttract HMW DNA kit (QIAGEN, Germantown, MD, USA) according to the manufacturer’s protocol. The quality and quantity of the gDNA were analyzed using a 5400 Fragment analyzer (Agilent Technologies, Santa Clara, CA, USA) and Qubit 2.0 Fluorometer (Invitrogen, Life Technologies, Carlsbad, CA, USA). DNA libraries for the southern elephant seal individuals were generated using the 10× Genomics Chromium technology according to the manufacturer’s instructions. Gel bead-in-emulsions (GEMs) were created from a library of Genome Gel Beads combined with 1.5 ng of gDNA in a Master Mix and partitioning oil, using the 10× Genomics Chromium Controller instrument with a micro-fluidic Genome chip (PN-120257). The GEMs were then subjected to an isothermal incubation step. Bar-coded DNA fragments were extracted and underwent Illumina library construction, as detailed in the Chromium Genome Reagent Kits Version 2 User Guide (PN-120258). Library yield was measured through the Qubit dsDNA HS assay kit (Thermo Fisher Scientific, Waltham, MA, USA). Library fragment size and distribution were measured using an Agilent 2100 Bioanalyzer High Sensitivity DNA chip (Agilent, Santa Clara, CA, USA). The DNA was sequenced on a NovaSeq with a 2 × 250 bp read metric, and ∼1.6 billion paired-end (PE) reads were generated (Table 1).
The de novo genome assembly was performed using the paired-end sequence reads from the partitioned library as input for the Supernova assembler (RRID:SCR 016756) v2.1.1 (10× Genomics, San Francisco, CA, USA) [9] with default parameters.
We built a de novo repeat library using RepeatModeler v1.0.3, (RepeatModeler, RRID: SCR 015027) [10], including the RECON (RECON, RRID:SCR 006345) [10] and RepeatScout v1.0.5 (RepeatScout, RRID:SCR 014653) [11] software with default parameters. Tandem Repeats Finder [12] was used to predict consensus sequences, classification information for each repeat, and tandem repeats, including simple repeats, satellites, and low-complexity repeats.
We annotated the genome with MAKER (RRID:SCR_005309) pipeline [13], which implements both ab initio prediction and homology-based gene annotation. Ab initio gene prediction was performed using SNAP (SNAP, RRID:SCR 002127) [14] and Augustus (Augustus: Gene Prediction, RRID:SCR 008417). MAKER was initially run in est2genome mode, which was based on transcripts for M. leonina generated from the previous results [8]. As further evidence for the annotation, reference proteins from other sequenced mammal species (Arctocephalus gazelle, Leptonychotes weddellii, Odobenus rosmarus, Tursiops truncates, Orcinus orca, Trichechus manatus latirostris, Canis lupus familiaris, Felis catus, Equus caballus, and Homo sapiens) were used in the analysis. The Exonerate software package was used to polish the MAKER alignments, which provides integrated information for the SNAP software program (SNAP, RRID:SCR 007936). The Infernal software package ver. 1.1 (INFERNAL, RRID:SCR 011809) [15] and covariance models (CMs) from the Rfam database v12.1 (Rfam, RRID:SCR 007891) [16] were used to identify other non-coding RNAs. The putative tRNA genes were identified using tRNAscan-SE v1.3.1 (tRNAscan-SE, RRID:SCR 010835) [17], which uses a CM that scores candidates based on their sequence and predicted secondary structures. Functional annotations were conducted by aligning them to the NCBI (National Center for Biotechnology Information) non-redundant protein (nr), Swissprot (Swissprot, RRID:SCR_002380) [18], TrEMBL (TrEMBL, RRID:SCR_002380) [18] using BLAST v2.2.31 [19] with a maximal e-value of 1e-5, and to the Pfam database (Pfam, RRID:SCR_004726) [20] using HMMer V3.0 [21].
This study including sample collection and experimental research conducted on these animals was performed according to the law on activities and environmental protection to Antarctic approved by the Minister of Foreign Affairs and Trade of the Republic of Korea.

3. Results and Discussion

The final assembled southern elephant seal, Mirounga leonina, genome obtained a 2.42 Gb genome (GC content: 41.51%) estimated the heterozygosity to be 0.301% with an N50 scaffold length of 54.233 Mb and a maximum scaffold length of 111.625 Mb (Table 2 and Supplementary Figure S1). The number of scaffolds in the southern elephant seal genome assembly was 1.115, and 54 scaffolds were over 10 Mb long and occupied 90.6% of our assembly (Supplementary Table S1.). Genome completeness was estimated with Benchmarking Universal Single-copy Orthologs (BUSCO) v3.0. (RRID:SCR 015008, version 3.0) by using the vertebrata_odb9 database [22]. Of the 2586 total BUSCO groups searched, 2472 BUSCO core genes were completely, and 49 were partially identified, leading to a total of 98.1% of BUSCO genes being found in the M. leonina genome (Supplementary Table S2).
The analysis revealed that 41.51% of the assembled M. leonina genome consisted of repeat sequences, of which 33.28% were transposable elements (TE), including LINEs (22.12%), SINEs (6.06%), LTRs (3.06%), and DNA transposons (1.81%) (Table 3). The Kimura divergence (κ-values) estimates the age and transposition history of TEs [23]. The LINEs are the most abundant transposable elements. Kimura distances (κ-values) were calculated for all TE copies of each element in order to estimate the “relative age” and transposition history of TEs [24]. The calculated Kimura divergence for all the TE copies of M. leonina is shown (Figure 2) and strongly shaped by LINEs (κ-values ≤5), which means that transposable elements, and LINEs in particular, are a recent development in the southern elephant seal genome due to very similar copies (low κ-values) are indicative of rather recent activity.
A total of 20,457 protein-coding genes in the M. leonina genome were annotated based on the combination of ab initio gene prediction, homology search, and transcript mapping. The total length of exons occupied 30.5 Mb, with an average of 8.2 exons per gene (Table 4).
The functional classification of Gene Ontology (GO) categories (Gene Ontology, RRID:SCR_002811) was performed using the BLAST2GO v5.25 pipeline (Blast2GO, RRID:SCR_005828) [25]. A Kyoto Encyclopedia of Genes and Genomes (KEGG, RRID:SCR_001120) [26] pathway annotation analysis was performed using the KEGG Automatic Annotation Server (KAAS), and EuKaryotic Orthologous Groups [27] were annotated with the KOG database proteins using BLASTp v2.2.31 with a maximal e-value of 1e-5. As a result, totally, 19,439 genes were annotated by a minimum of one database (Table 4). A total of 3354, 15,354, and 11,501 genes were annotated by the gene ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG) databases, and Eukaryotic Orthogous Groups (KOG), respectively (Supplementary Figures S2–S4). The genome browser, local blast database and assembly data are accessible on https://antagen.kopri.re.kr/project/project.php, accessed on: 01.01.2020.
Synteny analyses of chromosomes between California sea lion (Zalophus californianus) reference genome and southern elephant seal were performed using SyMAP v3.4 [28]. Southern elephant seal assembly is highly contiguous compared to the 18 chromosomes of the California sea lion (Figure 3 and Supplementary Table S3). A total of 2.22 Gb of southern elephant seal assembled genome were mapped to California sea lion chromosomes. OrthoVenn2 [29] was used to identify paralogy and orthology in annotated proteins of the southern elephant seal and other marine mammals. The deduced southern elephant seal proteins have 13,473 orthologous groups (containing 14,335 genes of the southern elephant seal) within five Pinnipedia species, and 148 orthogroups (containing 353 genes) are only for southern elephant seal (Figure 4).
The southern elephant seal assembly shows that a high-quality reference genome from 10× linked reads sequence will be essential for the further understanding of adaptive metabolism upon repeated breath-hold dives and the exploration of molecular mechanisms contributing to its unique biochemical and physiological characteristics. Additionally, our genomic data will provide valuable genetic resources for evolutionary studies regarding the divergence of pinnipeds.

Supplementary Materials

The following are available online at https://www.mdpi.com/2073-4425/11/2/160/s1, Table S1: Lengths of the southern elephant seal genome assembly (over 1Mb). Table S2: Completeness of the southern elephant seal genome assembly evaluated with benchmarking universal single-copy orthologs (BUSCO). Table S3: Mapping of Southern elephant seal genome assembly to California sea lion chromosomes (assembly zalCal2.2). Figure S1: Graph of the k-mer distribution (K = 27) using GenomeScope. Figure S2: Gene ontology (GO). The horizontal axis indicates classes of the second-level GO-annotation, and the vertical axis indicates the number of genes in each class. Figure S3: Statistics of Kyoto Encyclopedia of Genes and Genomes (KEGG) classifications. Figure S4: Eukaryotic Orthologous Groups (KOG) classification of the predicted genes. Results are grouped into 24 functional classes according to their functions. The horizontal axis indicates each class, and the vertical axis indicates the number of genes in each class.

Author Contributions

H.P. and J.-H.K. conceived the study. B.-M.K., Y.J.L., S.K., E.J., S.J.L., J.H.L., J.-H.K., Y.M.C. and H.P. performed genome sequencing, assembly, and annotation. E.J., S.J.L. and J.-H.K. performed experiments. B.-M.K., Y.J.L., J.-H.K. and H.P. mainly wrote the paper. All authors contributed to writing and editing the manuscript as well as providing supplementary information and producing the figures. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the ‘Ecosystem Structure and Function of Marine Protected Area (MPA) in Antarctica’ project (PM18060), funded by the Ministry of Oceans and Fisheries (20170336), Post-Polar Genomics Project: Functional genomic study for securing of polar useful genes (PE20040), and Korea University Grant to H.P.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Le Boeuf, B.J.; Laws, R.M. Elephant Seals: Population Ecology, Behavior, and Physiology; University of California Press: Berkeley, CA, USA, 1994. [Google Scholar]
  2. Crocker, D.E.; Champagne, C.D.; Fowler, M.A.; Houser, D.S. Adiposity and fat metabolism in lactating and fasting northern elephant seals. Adv. Nutr. 2014, 5, 57–64. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Rossier, B.C. Osmoregulation during long-term fasting in lungfish and elephant seal: Old and new lessons for the nephrologist. Nephron 2016, 134, 5–9. [Google Scholar] [CrossRef]
  4. Laws, R.M. Antarctic Seals: Research Methods and Techniques; Cambridge University Press: New York, NY, USA, 1993. [Google Scholar]
  5. Hofmeyr, G.J.G. Mirounga leonina. In The IUCN Red List of Threatened Species; Species Survival Commission: Gland, Switzerland, 2015. [Google Scholar]
  6. Biuw, M.; Boehme, L.; Guinet, C.; Hindell, M.; Costa, D.; Charrassin, J.B.; Roquet, F.; Bailleul, F.; Meredith, M.; Thorpe, S.; et al. Variations in behavior and condition of a Southern Ocean top predator in relation to in situ oceanographic conditions. Proc. Natl. Acad. Sci. USA 2007, 104, 13705. [Google Scholar] [CrossRef] [Green Version]
  7. Kovacs, K.M.; Aguilar, A.; Aurioles, D.; Burkanov, V.; Campagna, C.; Gales, N.; Gelatt, T.; Goldsworthy, S.D.; Goodman, S.J.; Hofmeyr, G.J.; et al. Global threats to pinnipeds. Mar. Mammal Sci. 2012, 28, 414–436. [Google Scholar] [CrossRef]
  8. Kim, B.-M.; Ahn, D.-H.; Kang, S.; Jeong, J.; Jo, E.; Kim, J.-H.; Rhee, J.-S.; Park, H. De novo Assembly and Annotation of the Blood Transcriptome of the Southern Elephant Seal Mirounga leonina from the South Shetland Islands, Antarctica. Ocean Sci. J. 2019, 54, 307–315. [Google Scholar] [CrossRef]
  9. Weisenfeld, N.I.; Kumar, V.; Shah, P.; Church, D.M.; Jaffe, D.B. Direct determination of diploid genome sequences. Genome Res. 2017, 27, 757–767. [Google Scholar] [CrossRef] [Green Version]
  10. Bao, Z.; Eddy, S.R. Automated de novo identification of repeat sequence families in sequenced genomes. Genome Res. 2002, 12, 1269–1276. [Google Scholar] [CrossRef] [Green Version]
  11. Price, A.L.; Jones, N.C.; Pevzner, P.A. De novo identification of repeat families in large genomes. Bioinformatics 2005, 21, i351–i358. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Benson, G. Tandem repeats finder: A program to analyze DNA sequences. Nucleic Acids Res. 1999, 27, 573–580. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Holt, C.; Yandell, M. MAKER2: An annotation pipeline and genome-database management tool for second-generation genome projects. BMC Bioinform. 2011, 12, 491. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Korf, I. Gene finding in novel genomes. BMC Bioinform. 2004, 5, 59. [Google Scholar] [CrossRef] [Green Version]
  15. Nawrocki, E.P.; Eddy, S.R. Infernal 1.1: 100-fold faster RNA homology searches. Bioinformatics 2013, 29, 2933–2935. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Griffiths-Jones, S.; Moxon, S.; Marshall, M.; Khanna, A.; Eddy, S.R.; Bateman, A. Rfam: Annotating non-coding RNAs in complete genomes. Nucleic Acids Res. 2005, 33, D121–D124. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Lowe, T.M.; Eddy, S.R. tRNAscan-SE: A program for improved detection of transfer RNA genes in genomic sequence. Nucleic Acids Res. 1997, 25, 955–964. [Google Scholar] [CrossRef] [PubMed]
  18. Boeckmann, B.; Bairoch, A.; Apweiler, R.; Blatter, M.-C.; Estreicher, A.; Gasteiger, E.; Martin, M.J.; Michoud, K.; O’Donovan, C.; Phan, I. The SWISS-PROT protein knowledgebase and its supplement TrEMBL in 2003. Nucleic Acids Res. 2003, 31, 365–370. [Google Scholar] [CrossRef] [PubMed]
  19. Altschul, S.F.; Gish, W.; Miller, W.; Myers, E.W.; Lipman, D.J. Basic local alignment search tool. J. Mol. Biol. 1990, 215, 403–410. [Google Scholar] [CrossRef]
  20. El-Gebali, S.; Mistry, J.; Bateman, A.; Eddy, S.R.; Luciani, A.; Potter, S.C.; Qureshi, M.; Richardson, L.J.; Salazar, G.A.; Smart, A. The Pfam protein families database in 2019. Nucleic Acids Res. 2018, 47, D427–D432. [Google Scholar] [CrossRef]
  21. Eddy, S.R.; Mitchison, G.; Durbin, R. Maximum discrimination hidden Markov models of sequence consensus. J. Comput. Biol. 1995, 2, 9–23. [Google Scholar] [CrossRef]
  22. Simão, F.A.; Waterhouse, R.M.; Ioannidis, P.; Kriventseva, E.V.; Zdobnov, E.M. BUSCO: Assessing genome assembly and annotation completeness with single-copy orthologs. Bioinformatics 2015, 31, 3210–3212. [Google Scholar] [CrossRef] [Green Version]
  23. Kimura, M. A simple method for estimating evolutionary rates of base substitutions through comparative studies of nucleotide sequences. J. Mol. Evol. 1980, 16, 111–120. [Google Scholar] [CrossRef]
  24. Chalopin, D.; Naville, M.; Plard, F.; Galiana, D.; Volff, J.-N. Comparative analysis of transposable elements highlights mobilome diversity and evolution in vertebrates. Genome Biol. Evol. 2015, 7, 567–580. [Google Scholar] [CrossRef] [PubMed]
  25. Conesa, A.; Götz, S.; García-Gómez, J.M.; Terol, J.; Talón, M.; Robles, M. Blast2GO: A universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics 2005, 21, 3674–3676. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Kanehisa, M.; Goto, S. KEGG: Kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000, 28, 27–30. [Google Scholar] [CrossRef] [PubMed]
  27. Tatusov, R.L.; Natale, D.A.; Garkavtsev, I.V.; Tatusova, T.A.; Shankavaram, U.T.; Rao, B.S.; Kiryutin, B.; Galperin, M.Y.; Fedorova, N.D.; Koonin, E.V. The COG database: New developments in phylogenetic classification of proteins from complete genomes. Nucleic Acids Res. 2001, 29, 22–28. [Google Scholar] [CrossRef] [PubMed]
  28. Soderlund, C.; Bomhoff, M.; Nelson, W.M. SyMAP v3.4: A turnkey synteny system with application to plant genomes. Nucleic Acids Res. 2011, 39, e68. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Xu, L.; Dong, Z.; Fang, L.; Luo, Y.; Wei, Z.; Guo, H.; Zhang, G.; Gu, Y.Q.; Coleman-Derr, D.; Xia, Q. OrthoVenn2: A web server for whole-genome comparison and annotation of orthologous clusters across multiple species. Nucleic Acids Res. 2019, 47, W52–W58. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Example of Mirounga leonina, southern elephant seal. Photo by Hyun Park.
Figure 1. Example of Mirounga leonina, southern elephant seal. Photo by Hyun Park.
Genes 11 00160 g001
Figure 2. Kimura distance-based copy divergence analysis of transposable elements in the M. leonina genome. Graphs represent genome coverage (Y-axis) for each type of TE (DNA transposons, SINE, LINE, and LTR retrotransposons).
Figure 2. Kimura distance-based copy divergence analysis of transposable elements in the M. leonina genome. Graphs represent genome coverage (Y-axis) for each type of TE (DNA transposons, SINE, LINE, and LTR retrotransposons).
Genes 11 00160 g002
Figure 3. Synteny blocks between Southern elephant seal scaffolds and 19 chromosomes of California sea lion (assembly zalCal2.2). Southern elephant seal scaffolds over 10 Mb in length were selected. Alignment was accomplished with SyMAP v3.4. The color blocks represent California sea lion chromosomes and empty blocks represent Southern elephant seal scaffolds. Green and blue bars represent genes of California sea lion and Southern elephant seal, respectively. Connections within the circle represent alignment between the two assemblies.
Figure 3. Synteny blocks between Southern elephant seal scaffolds and 19 chromosomes of California sea lion (assembly zalCal2.2). Southern elephant seal scaffolds over 10 Mb in length were selected. Alignment was accomplished with SyMAP v3.4. The color blocks represent California sea lion chromosomes and empty blocks represent Southern elephant seal scaffolds. Green and blue bars represent genes of California sea lion and Southern elephant seal, respectively. Connections within the circle represent alignment between the two assemblies.
Genes 11 00160 g003
Figure 4. Venn diagram representing paralogous and orthologous groups between Southern elephant seal (Mirounga leonina), Pacific walrus (Odobenus rosmarus), Weddell seal (Leptonychotes weddellii), Steller sea lion (Eumetopias jubatus) and California sea lion (Zalophus californianus) obtained with OrthoVenn2.
Figure 4. Venn diagram representing paralogous and orthologous groups between Southern elephant seal (Mirounga leonina), Pacific walrus (Odobenus rosmarus), Weddell seal (Leptonychotes weddellii), Steller sea lion (Eumetopias jubatus) and California sea lion (Zalophus californianus) obtained with OrthoVenn2.
Genes 11 00160 g004
Table 1. Sequencing data generated for Mirounga leonina genome assembly and annotation.
Table 1. Sequencing data generated for Mirounga leonina genome assembly and annotation.
Number of ReadsTotal Read Bases
Genome1,606,582,076242,593,893,476
Transcriptome80,958,40019,430,016,000
Table 2. Statistics for Mirounga leonina genome assembly.
Table 2. Statistics for Mirounga leonina genome assembly.
AssemblySupernova 2.0
Number of scaffolds1115
Total size of scaffolds2,417,339,903
Longest scaffolds111,625,095
N50 scaffolds length54,232,831
Number of scaffolds >10M54
Gap (%)0.65
Table 3. Statistics for annotated Mirounga leonina transposable elements.
Table 3. Statistics for annotated Mirounga leonina transposable elements.
Number
of Elements
Length (bp)Percentage
of Sequence (%)
SINEs:888,301148,028,3856.06
MIRs254,250 35,222,085 1.44
LINEs:1,624,363540,626,27822.12
LINE11,473,069 508,962,965 20.83
LINE2150,067 31,135,376 1.27
L3/CR1634 81,165 0.00
LTR elements:268,71374,825,3843.06
ERVL81,307 29,432,080 1.20
ERVL-MaLRs134,246 29,421,479 1.20
ERV_classI52,451 15,768,764 0.65
ERV_classII370 158,839 0.01
DNA elements:261,76944,140,5971.81
hAT-Charlie146,299 24,083,632 0.99
TcMar-Tigger43,576 10,220,361 0.42
Unclassified:20,6305,747,8080.24
Total interspersed repeats: 813,368,45233.28
Small RNA:637,161113,111,8974.63
Satellites:186785,9130.00
Simple repeats:618,43726,612,5571.09
Low complexity:99,3115,218,0080.21
Bold are total of each classes of transposable elements.
Table 4. Mirounga leonina genome annotation statistics.
Table 4. Mirounga leonina genome annotation statistics.
Annotation Database Annotated Number Percentage (%)
No. of genes20,457
nr annotation 19,90197.3
GO annotation 310915.2
KEGG annotation 11,50156.2
KOG annotation 15,49875.8
Pfam annotation 16,73381.8
Swissprot annotation 19,65396.1
TrEMBL annotation 17,60086.0
CountLength Sum (bp)
Exon 168,37530,497,777
CDS 167,91928,718,340

Share and Cite

MDPI and ACS Style

Kim, B.-M.; Lee, Y.J.; Kim, J.-H.; Kim, J.-H.; Kang, S.; Jo, E.; Lee, S.J.; Lee, J.H.; Chi, Y.M.; Park, H. The Genome Assembly and Annotation of the Southern Elephant Seal Mirounga leonina. Genes 2020, 11, 160. https://doi.org/10.3390/genes11020160

AMA Style

Kim B-M, Lee YJ, Kim J-H, Kim J-H, Kang S, Jo E, Lee SJ, Lee JH, Chi YM, Park H. The Genome Assembly and Annotation of the Southern Elephant Seal Mirounga leonina. Genes. 2020; 11(2):160. https://doi.org/10.3390/genes11020160

Chicago/Turabian Style

Kim, Bo-Mi, Yoon Jin Lee, Jeong-Hoon Kim, Jin-Hyoung Kim, Seunghyun Kang, Euna Jo, Seung Jae Lee, Jun Hyuck Lee, Young Min Chi, and Hyun Park. 2020. "The Genome Assembly and Annotation of the Southern Elephant Seal Mirounga leonina" Genes 11, no. 2: 160. https://doi.org/10.3390/genes11020160

APA Style

Kim, B. -M., Lee, Y. J., Kim, J. -H., Kim, J. -H., Kang, S., Jo, E., Lee, S. J., Lee, J. H., Chi, Y. M., & Park, H. (2020). The Genome Assembly and Annotation of the Southern Elephant Seal Mirounga leonina. Genes, 11(2), 160. https://doi.org/10.3390/genes11020160

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