A draft fur seal genome provides insights into factors affecting SNP validation and how to mitigate them

Multi-omics for studying and understanding polar life

Polar ecosystems are experiencing amongst the most rapid rates of regional warming on Earth. Here, we discuss 'omics' approaches to investigate polar biodiversity, including the current state of the art, future perspectives and recommendations. We propose a community road map to generate and more fully exploit multi-omics data from polar organisms. These data are needed for the comprehensive evaluation of polar biodiversity and to reveal how life evolved and adapted to permanently cold environments with extreme seasonality. We argue that concerted action is required to mitigate the impact of warming on polar ecosystems via conservation efforts, to sustainably manage these unique habitats and their ecosystem services, and for the sustainable bioprospecting of novel genes and compounds for societal gain.

Intronic primers reveal unexpectedly high major histocompatibility complex diversity in Antarctic fur seals

The major histocompatibility complex (MHC) is a group of genes comprising one of the most important components of the vertebrate immune system. Consequently, there has been much interest in characterising MHC variation and its relationship with fitness in a variety of species. Due to the exceptional polymorphism of MHC genes, careful PCR primer design is crucial for capturing all of the allelic variation present in a given species. We therefore developed intronic primers to amplify the full-length 267 bp protein-coding sequence of the MHC class II DQB exon 2 in the Antarctic fur seal. We then characterised patterns of MHC variation among mother-offspring pairs from two breeding colonies and detected 19 alleles among 771 clone sequences from 56 individuals. The distribution of alleles within and among individuals was consistent with a single-copy, classical DQB locus showing Mendelian inheritance. Amino acid similarity at the MHC was significantly associated with genome-wide relatedness, but no relationship was found between MHC heterozygosity and genome-wide heterozygosity. Finally, allelic diversity was several times higher than reported by a previous study based on partial exon sequences. This difference appears to be related to allele-specific amplification bias, implying that primer design can strongly impact the inference of MHC diversity.

Signatures of selection on mitonuclear integrated genes uncover hidden mitogenomic variation in fur seals

Nuclear copies of mitochondrial genes (numts) are commonplace in vertebrate genomes and have been characterized in many species. However, relatively little attention has been paid to understanding their evolutionary origins and to disentangling alternative sources of insertions. Numts containing genes with intact mitochondrial reading frames represent good candidates for this purpose. The sequences of the genes they contain can be compared to their mitochondrial homologs to characterize synonymous to non-synonymous substitution rates, which can shed light on the selection pressures these genes have been subjected to. Here, we characterise 25 numts in the Antarctic fur seal (Arctocephalus gazella) genome. Among those containing genes with intact mitochondrial reading frames, three carry multiple substitutions in comparison to their mitochondrial homologs. Our analyses reveal that one represents a historic insertion subjected to strong purifying selection since it colonized the Otarioidea in a genomic region enriched in retrotransposons. By contrast, the other two numts appear to be more recent and their large number of substitutions can be attributed to non-canonical insertions, either the integration of heteroplasmic mtDNA or hybridization. Our study sheds new light on the evolutionary history of pinniped numts and uncovers the presence of hidden sources of mitonuclear variation.

Phylogenomic Discordance in the Eared Seals is best explained by Incomplete Lineage Sorting following Explosive Radiation in the Southern Hemisphere

The phylogeny and systematics of fur seals and sea lions (Otariidae) have long been studied with diverse data types, including an increasing amount of molecular data. However, only a few phylogenetic relationships have reached acceptance because of strong gene-tree species tree discordance. Divergence times estimates in the group also vary largely between studies. These uncertainties impeded the understanding of the biogeographical history of the group, such as when and how trans-equatorial dispersal and subsequent speciation events occurred. Here, we used high-coverage genome-wide sequencing for 14 of the 15 species of Otariidae to elucidate the phylogeny of the family and its bearing on the taxonomy and biogeographical history. Despite extreme topological discordance among gene trees, we found a fully supported species tree that agrees with the few well-accepted relationships and establishes monophyly of the genus Arctocephalus. Our data support a relatively recent trans-hemispheric dispersal at the base of a southern clade, which rapidly diversified into six major lineages between 3 and 2.5 Ma. Otaria diverged first, followed by Phocarctos and then four major lineages within Arctocephalus. However, we found Zalophus to be nonmonophyletic, with California (Zalophus californianus) and Steller sea lions (Eumetopias jubatus) grouping closer than the Galapagos sea lion (Zalophus wollebaeki) with evidence for introgression between the two genera. Overall, the high degree of genealogical discordance was best explained by incomplete lineage sorting resulting from quasi-simultaneous speciation within the southern clade with introgresssion playing a subordinate role in explaining the incongruence among and within prior phylogenetic studies of the family. [Hybridization; ILS; phylogenomics; Pleistocene; Pliocene; monophyly.].

An 85K SNP Array Uncovers Inbreeding and Cryptic Relatedness in an Antarctic Fur Seal Breeding Colony

High density single nucleotide polymorphism (SNP) arrays allow large numbers of individuals to be rapidly and cost-effectively genotyped at large numbers of genetic markers. However, despite being widely used in studies of humans and domesticated plants and animals, SNP arrays are lacking for most wild organisms. We developed a custom 85K Affymetrix Axiom array for an intensively studied pinniped, the Antarctic fur seal (Arctocephalus gazella). SNPs were discovered from a combination of genomic and transcriptomic resources and filtered according to strict criteria. Out of a total of 85,359 SNPs tiled on the array, 75,601 (88.6%) successfully converted and were polymorphic in 270 animals from a breeding colony at Bird Island in South Georgia. Evidence was found for inbreeding, with three genomic inbreeding coefficients being strongly intercorrelated and the proportion of the genome in runs of homozygosity being non-zero in all individuals. Furthermore, analysis of genomic relatedness coefficients identified previously unknown first-degree relatives and multiple second-degree relatives among a sample of ostensibly unrelated individuals. Such “cryptic relatedness” within fur seal breeding colonies may increase the likelihood of consanguineous matings and could therefore have implications for understanding fitness variation and mate choice. Finally, we demonstrate the cross-amplification potential of the array in three related pinniped species. Overall, our SNP array will facilitate future studies of Antarctic fur seals and has the potential to serve as a more general resource for the wider pinniped research community.

Chromosomal-level genome assembly of the scimitar-horned oryx: Insights into diversity and demography of a species extinct in the wild

Captive populations provide a valuable insurance against extinctions in the wild. However, they are also vulnerable to the negative impacts of inbreeding, selection and drift. Genetic information is therefore considered a critical aspect of conservation management. Recent developments in sequencing technologies have the potential to improve the outcomes of management programmes; however, the transfer of these approaches to applied conservation has been slow. The scimitar-horned oryx (Oryx dammah) is a North African antelope that has been extinct in the wild since the early 1980s and is the focus of a large-scale and long-term reintroduction project. To enable the selection of suitable founder individuals, facilitate post-release monitoring and improve captive breeding management, comprehensive genomic resources are required. Here, we used 10X Chromium sequencing together with Hi-C contact mapping to develop a chromosomal-level genome assembly for the species. The resulting assembly contained 29 chromosomes with a scaffold N50 of 100.4 Mb, and displayed strong chromosomal synteny with the cattle genome. Using resequencing data from six additional individuals, we demonstrated relatively high genetic diversity in the scimitar-horned oryx compared to other mammals, despite it having experienced a strong founding event in captivity. Additionally, the level of diversity across populations varied according to management strategy. Finally, we uncovered a dynamic demographic history that coincided with periods of climate variation during the Pleistocene. Overall, our study provides a clear example of how genomic data can uncover valuable insights into captive populations and contributes important resources to guide future management decisions of an endangered species.

Complete mitochondrial genome of the Antarctic fur seal (Arctocephalus gazella)

The Antarctic fur seal (Arctocephalus gazella) is an abundant Antarctic otariid. Here, we present the complete mitochondrial DNA sequence of this species, which includes 13 protein-coding genes, 22 transfer RNA genes, 2 ribosomal RNA genes, and the control region for a total length of 16,156 bp. A phylogenetic analysis including all 25 publically available pinniped mitogenomes nested the Antarctic fur seal within the Otariid clade, which was clearly resolved from the Phocidae and Odobenidae.

Draft Genome Sequence of Phocine Herpesvirus 1 Isolated from the Brain of a Harbor Seal

Phocine herpesvirus 1 (PhHV-1) is a viral pathogen with high prevalence, morbidity, and mortality in harbor seals. In this study, we used a metagenomic approach to assemble the PhHV-1 genome from the brain tissue of a harbor seal.

Phocine herpesvirus 1 (PhHV-1) is a viral pathogen with high prevalence, morbidity, and mortality in harbor seals. In this study, we used a metagenomic approach to assemble the PhHV-1 genome from the brain tissue of a harbor seal. Here, we present a 119-kb draft genome of PhHV-1 comprising 76 open reading frames.

Deciphering the evolutionary signatures of pinnipeds using novel genome sequences: The first genomes of Phoca largha, Callorhinus ursinus, and Eumetopias jubatus

The pinnipeds, which comprise seals, sea lions, and walruses, are a remarkable group of marine animals with unique adaptations to semi-aquatic life. However, their genomes are poorly characterized. In this study, we sequenced and characterized the genomes of three pinnipeds (Phoca largha, Callorhinus ursinus, and Eumetopias jubatus), focusing on site-wise sequence changes. We detected rapidly evolving genes in pinniped lineages and substitutions unique to pinnipeds associated with amphibious sound perception. Phenotypic convergence-related sequence convergences are not common in marine mammals. For example, FASN, KCNA5, and IL17RA contain substitutions specific to pinnipeds, yet are potential candidates of phenotypic convergence (blubber, response to hypoxia, and immunity to pathogens) in all marine mammals. The outcomes of this study will provide insight into targets for future studies of convergent evolution or gene function.

A global cline in a colour polymorphism suggests a limited contribution of gene flow towards the recovery of a heavily exploited marine mammal

Evaluating how populations are connected by migration is important for understanding species resilience because gene flow can facilitate recovery from demographic declines. We therefore investigated the extent to which migration may have contributed to the global recovery of the Antarctic fur seal (Arctocephalus gazella), a circumpolar distributed marine mammal that was brought to the brink of extinction by the sealing industry in the eighteenth and nineteenth centuries. It is widely believed that animals emigrating from South Georgia, where a relict population escaped sealing, contributed to the re-establishment of formerly occupied breeding colonies across the geographical range of the species. To investigate this, we interrogated a genetic polymorphism (S291F) in the melanocortin 1 receptor gene, which is responsible for a cream-coloured phenotype that is relatively abundant at South Georgia and which appears to have recently spread to localities as far afield as Marion Island in the sub-Antarctic Indian Ocean. By sequencing a short region of this gene in 1492 pups from eight breeding colonies, we showed that S291F frequency rapidly declines with increasing geographical distance from South Georgia, consistent with locally restricted gene flow from South Georgia mainly to the South Shetland Islands and Bouvetøya. The S291F allele was not detected farther afield, suggesting that although emigrants from South Georgia may have been locally important, they are unlikely to have played a major role in the recovery of geographically more distant populations.

Ancient convergent losses of Paraoxonase 1 yield potential risks for modern marine mammals

Mammals diversified by colonizing drastically different environments, with each transition yielding numerous molecular changes, including losses of protein function. Though not initially deleterious, these losses could subsequently carry deleterious pleiotropic consequences. We have used phylogenetic methods to identify convergent functional losses across independent marine mammal lineages. In one extreme case, Paraoxonase 1 (PON1) accrued lesions in all marine lineages, while remaining intact in all terrestrial mammals. These lesions coincide with PON1 enzymatic activity loss in marine species' blood plasma. This convergent loss is likely explained by parallel shifts in marine ancestors' lipid metabolism and/or bloodstream oxidative environment affecting PON1's role in fatty acid oxidation. PON1 loss also eliminates marine mammals' main defense against neurotoxicity from specific man-made organophosphorus compounds, implying potential risks in modern environments.

Analysis of ZP1 gene reveals differences in zona pellucida composition in carnivores

The zona pellucida (ZP) is an extracellular envelope that surrounds mammalian oocytes. This coat participates in the interaction between gametes, induction of the acrosome reaction, block of polyspermy and protection of the oviductal embryo. Previous studies suggested that carnivore ZP was formed by three glycoproteins (ZP2, ZP3 and ZP4), with ZP1 being a pseudogene. However, a recent study in the cat found that all four proteins were expressed. In the present study, in silico and molecular analyses were performed in several carnivores to clarify the ZP composition in this order of mammals. The in silico analysis demonstrated the presence of the ZP1 gene in five carnivores: cheetah, panda, polar bear, tiger and walrus, whereas in the Antarctic fur seal and the Weddell seal there was evidence of pseudogenisation. Molecular analysis showed the presence of four ZP transcripts in ferret ovaries (ZP1, ZP2, ZP3 and ZP4) and three in fox ovaries (ZP2, ZP3 and ZP4). Analysis of the fox ZP1 gene showed the presence of a stop codon. The results strongly suggest that all four ZP genes are expressed in most carnivores, whereas ZP1 pseudogenisation seems to have independently affected three families (Canidae, Otariidae and Phocidae) of the carnivore tree.

RAD Sequencing and a Hybrid Antarctic Fur Seal Genome Assembly Reveal Rapidly Decaying Linkage Disequilibrium, Global Population Structure and Evidence for Inbreeding

Recent advances in high throughput sequencing have transformed the study of wild organisms by facilitating the generation of high quality genome assemblies and dense genetic marker datasets. These resources have the potential to significantly advance our understanding of diverse phenomena at the level of species, populations and individuals, ranging from patterns of synteny through rates of linkage disequilibrium (LD) decay and population structure to individual inbreeding. Consequently, we used PacBio sequencing to refine an existing Antarctic fur seal (Arctocephalus gazella) genome assembly and genotyped 83 individuals from six populations using restriction site associated DNA (RAD) sequencing. The resulting hybrid genome comprised 6,169 scaffolds with an N50 of 6.21 Mb and provided clear evidence for the conservation of large chromosomal segments between the fur seal and dog (Canis lupus familiaris). Focusing on the most extensively sampled population of South Georgia, we found that LD decayed rapidly, reaching the background level by around 400 kb, consistent with other vertebrates but at odds with the notion that fur seals experienced a strong historical bottleneck. We also found evidence for population structuring, with four main Antarctic island groups being resolved. Finally, appreciable variance in individual inbreeding could be detected, reflecting the strong polygyny and site fidelity of the species. Overall, our study contributes important resources for future genomic studies of fur seals and other pinnipeds while also providing a clear example of how high throughput sequencing can generate diverse biological insights at multiple levels of organization.

Genetic sex assignment in wild populations using genotyping-by-sequencing data: A statistical threshold approach

Establishing the sex of individuals in wild systems can be challenging and often requires genetic testing. Genotyping-by-sequencing (GBS) and other reduced-representation DNA sequencing (RRS) protocols (e.g., RADseq, ddRAD) have enabled the analysis of genetic data on an unprecedented scale. Here, we present a novel approach for the discovery and statistical validation of sex-specific loci in GBS data sets. We used GBS to genotype 166 New Zealand fur seals (NZFS, Arctocephalus forsteri) of known sex. We retained monomorphic loci as potential sex-specific markers in the locus discovery phase. We then used (i) a sex-specific locus threshold (SSLT) to identify significantly male-specific loci within our data set; and (ii) a significant sex-assignment threshold (SSAT) to confidently assign sex in silico the presence or absence of significantly male-specific loci to individuals in our data set treated as unknowns (98.9% accuracy for females; 95.8% for males, estimated via cross-validation). Furthermore, we assigned sex to 86 individuals of true unknown sex using our SSAT and assessed the effect of SSLT adjustments on these assignments. From 90 verified sex-specific loci, we developed a panel of three sex-specific PCR primers that we used to ascertain sex independently of our GBS data, which we show amplify reliably in at least two other pinniped species. Using monomorphic loci normally discarded from large SNP data sets is an effective way to identify robust sex-linked markers for nonmodel species. Our novel pipeline can be used to identify and statistically validate monomorphic and polymorphic sex-specific markers across a range of species and RRS data sets.

Transcriptomic SNP discovery for custom genotyping arrays: impacts of sequence data, SNP calling method and genotyping technology on the probability of validation success

BACKGROUND

Single nucleotide polymorphism (SNP) discovery is an important goal of many studies. However, the number of 'putative' SNPs discovered from a sequence resource may not provide a reliable indication of the number that will successfully validate with a given genotyping technology. For this it may be necessary to account for factors such as the method used for SNP discovery and the type of sequence data from which it originates, suitability of the SNP flanking sequences for probe design, and genomic context. To explore the relative importance of these and other factors, we used Illumina sequencing to augment an existing Roche 454 transcriptome assembly for the Antarctic fur seal (Arctocephalus gazella). We then mapped the raw Illumina reads to the new hybrid transcriptome using BWA and BOWTIE2 before calling SNPs with GATK. The resulting markers were pooled with two existing sets of SNPs called from the original 454 assembly using NEWBLER and SWAP454. Finally, we explored the extent to which SNPs discovered using these four methods overlapped and predicted the corresponding validation outcomes for both Illumina Infinium iSelect HD and Affymetrix Axiom arrays.

RESULTS

Collating markers across all discovery methods resulted in a global list of 34,718 SNPs. However, concordance between the methods was surprisingly poor, with only 51.0 % of SNPs being discovered by more than one method and 13.5 % being called from both the 454 and Illumina datasets. Using a predictive modeling approach, we could also show that SNPs called from the Illumina data were on average more likely to successfully validate, as were SNPs called by more than one method. Above and beyond this pattern, predicted validation outcomes were also consistently better for Affymetrix Axiom arrays.

CONCLUSIONS

Our results suggest that focusing on SNPs called by more than one method could potentially improve validation outcomes. They also highlight possible differences between alternative genotyping technologies that could be explored in future studies of non-model organisms.

Genomic Methods Take the Plunge: Recent Advances in High-Throughput Sequencing of Marine Mammals

The dramatic increase in the application of genomic techniques to non-model organisms (NMOs) over the past decade has yielded numerous valuable contributions to evolutionary biology and ecology, many of which would not have been possible with traditional genetic markers. We review this recent progression with a particular focus on genomic studies of marine mammals, a group of taxa that represent key macroevolutionary transitions from terrestrial to marine environments and for which available genomic resources have recently undergone notable rapid growth. Genomic studies of NMOs utilize an expanding range of approaches, including whole genome sequencing, restriction site-associated DNA sequencing, array-based sequencing of single nucleotide polymorphisms and target sequence probes (e.g., exomes), and transcriptome sequencing. These approaches generate different types and quantities of data, and many can be applied with limited or no prior genomic resources, thus overcoming one traditional limitation of research on NMOs. Within marine mammals, such studies have thus far yielded significant contributions to the fields of phylogenomics and comparative genomics, as well as enabled investigations of fitness, demography, and population structure. Here we review the primary options for generating genomic data, introduce several emerging techniques, and discuss the suitability of each approach for different applications in the study of NMOs.

Born blonde: a recessive loss-of-function mutation in the melanocortin 1 receptor is associated with cream coat coloration in Antarctic fur seals

Although the genetic basis of color variation has been extensively studied in humans and domestic animals, the genetic polymorphisms responsible for different color morphs remain to be elucidated in many wild vertebrate species. For example, hypopigmentation has been observed in numerous marine mammal species but the underlying mutations have not been identified. A particularly compelling candidate gene for explaining color polymorphism is the melanocortin 1 receptor (MC1R), which plays a key role in the regulation of pigment production. We therefore used Antarctic fur seals (Arctocephalus gazella) as a highly tractable marine mammal system with which to test for an association between nucleotide variation at the MC1R and melanin‐based coat color phenotypes. By sequencing 70 wild‐type individuals with dark‐colored coats and 26 hypopigmented individuals with cream‐colored coats, we identified a nonsynonymous mutation that results in the substitution of serine with phenylalanine at an evolutionarily highly conserved structural domain. All of the hypopigmented individuals were homozygous for the allele coding for phenylalanine, consistent with a recessive loss‐of‐function allele. In order to test for cryptic population structure, which can generate artefactual associations, and to evaluate whether homozygosity at the MC1R could be indicative of low genome‐wide heterozygosity, we also genotyped all of the individuals at 50 polymorphic microsatellite loci. We were unable to detect any population structure and also found that wild‐type and hypopigmented individuals did not differ significantly in their standardized multilocus heterozygosity. Such a lack of association implies that hypopigmented individuals are unlikely to suffer disproportionately from inbreeding depression, and hence, we have no reason to believe that they are at a selective disadvantage in the wider population.