No evidence for a role of MHC class II genotype in the chemical encoding of heterozygosity and relatedness in Antarctic fur seals.

Despite decades of research, surprisingly little is known about the mechanism(s) by which an individual's genotype is encoded in odour. Many studies have focused on the role of the major histocompatibility complex (MHC) owing to its importance for survival and mate choice. However, the salience of MHC-mediated odours compared to chemicals influenced by the rest of the genome remains unclear, especially in wild populations where it is challenging to quantify and control for the effects of the genomic background. We addressed this issue in Antarctic fur seals by analysing skin swabs together with full-length MHC DQB II exon 2 sequences and data from 41 genome-wide distributed microsatellites. We did not find any effects of MHC relatedness on chemical similarity and there was also no relationship between MHC heterozygosity and chemical diversity. However, multilocus heterozygosity showed a significant positive association with chemical diversity, even after controlling for MHC heterozygosity. Our results appear to rule out a dominant role of the MHC in the chemical encoding of genetic information in a wild vertebrate population and highlight the need for genome-wide approaches to elucidate the mechanism(s) and specific genes underlying genotype-odour associations.

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.].

Climate change selects for heterozygosity in a declining fur seal population.

Global environmental change is expected to alter selection pressures in many biological systems, but the long-term molecular and life history data required to quantify changes in selection are rare. An unusual opportunity is afforded by three decades of individual-based data collected from a declining population of Antarctic fur seals in the South Atlantic. Here, climate change has reduced prey availability and caused a significant decline in seal birth weight. However, the mean age and size of females recruiting into the breeding population are increasing. We show that such females have significantly higher heterozygosity (a measure of within-individual genetic variation) than their non-recruiting siblings and their own mothers. Thus, breeding female heterozygosity has increased by 8.5% per generation over the last two decades. Nonetheless, as heterozygosity is not inherited from mothers to daughters, substantial heterozygote advantage is not transmitted from one generation to the next and the decreasing viability of homozygous individuals causes the population to decline. Our results provide compelling evidence that selection due to climate change is intensifying, with far-reaching consequences for demography as well as phenotypic and genetic variation.

Experimental validation of in silico predicted RAD locus frequencies using genomic resources and short read data from a model marine mammal.

BACKGROUND: Restriction site-associated DNA sequencing (RADseq) has revolutionized the study of wild organisms by allowing cost-effective genotyping of thousands of loci. However, for species lacking reference genomes, it can be challenging to select the restriction enzyme that offers the best balance between the number of obtained RAD loci and depth of coverage, which is crucial for a successful outcome. To address this issue, PredRAD was recently developed, which uses probabilistic models to predict restriction site frequencies from a transcriptome assembly or other sequence resource based on either GC content or mono-, di- or trinucleotide composition. This program generates predictions that are broadly consistent with estimates of the true number of restriction sites obtained through in silico digestion of available reference genome assemblies. However, in practice the actual number of loci obtained could potentially differ as incomplete enzymatic digestion or patchy sequence coverage across the genome might lead to some loci not being represented in a RAD dataset, while erroneous assembly could potentially inflate the number of loci. To investigate this, we used genome and transcriptome assemblies together with RADseq data from the Antarctic fur seal (Arctocephalus gazella) to compare PredRAD predictions with empirical estimates of the number of loci obtained via in silico digestion and from de novo assemblies. RESULTS: PredRAD yielded consistently higher predicted numbers of restriction sites for the transcriptome assembly relative to the genome assembly. The trinucleotide and dinucleotide models also predicted higher frequencies than the mononucleotide or GC content models. Overall, the dinucleotide and trinucleotide models applied to the transcriptome and the genome assemblies respectively generated predictions that were closest to the number of restriction sites estimated by in silico digestion. Furthermore, the number of de novo assembled RAD loci mapping to restriction sites was similar to the expectation based on in silico digestion. CONCLUSIONS: Our study reveals generally high concordance between PredRAD predictions and empirical estimates of the number of RAD loci. This further supports the utility of PredRAD, while also suggesting that it may be feasible to sequence and assemble the majority of RAD loci present in an organism's genome.

Chemical fingerprints encode mother-offspring similarity, colony membership, relatedness, and genetic quality in fur seals.

Chemical communication underpins virtually all aspects of vertebrate social life, yet remains poorly understood because of its highly complex mechanistic basis. We therefore used chemical fingerprinting of skin swabs and genetic analysis to explore the chemical cues that may underlie mother-offspring recognition in colonially breeding Antarctic fur seals. By sampling mother-offspring pairs from two different colonies, using a variety of statistical approaches and genotyping a large panel of microsatellite loci, we show that colony membership, mother-offspring similarity, heterozygosity, and genetic relatedness are all chemically encoded. Moreover, chemical similarity between mothers and offspring reflects a combination of genetic and environmental influences, the former partly encoded by substances resembling known pheromones. Our findings reveal the diversity of information contained within chemical fingerprints and have implications for understanding mother-offspring communication, kin recognition, and mate choice.

Sex beyond species: the first genetically analyzed case of intergeneric fertile hybridization in pinnipeds.

A species, according to the biological concept, is a natural group of potentially interbreeding individuals isolated by diverse mechanisms. Hybridization is considered the production of offspring resulting from the interbreeding of two genetically distinct taxa. It has been documented in over 10% of wild animals, and at least in 34 cases for Artic marine mammals. In Otariids, intergeneric hybridization has been reported though neither confirming it through genetic analyses nor presenting evidence of fertile offspring. In this study, we report the finding of a hybrid adult female between a South American fur seal (Arctocephalus australis) and a South American sea lion (Otaria byronia), and its offspring, a male pup, in Uruguay. Further based on morphological constraints and breeding seasons, sex-biased hybridization between the two species is hypothesized. Morphological and genetic (nuclear and mitochondrial) results confirm de hybrid nature of the female-pup pair. Here we discuss a genetic dilution effect, considering other hybridization events must be occurring, and how isolation mechanisms could be circumvented. Moreover, the results obtained from stable isotope analysis suggest feeding habits may be a trait transmitted maternally, leading to consider broader issues regarding hybridization as an evolutionary innovation phenomenon.

Myth or relict: Does ancient DNA detect the enigmatic Upland seal?

The biological status of the so-called 'Upland seal' has remained contentious ever since historical records described a distinct seal from the uplands of New Zealand's (NZ) remote sub-Antarctic islands. Subsequent genetic surveys of the NZ fur seal (Arctocephalus forsteri) detected two highly-divergent mtDNA clades, hypothesized to represent a post-sealing hybrid swarm between 'mainland' (Australia-NZ; A. forsteri) and sub-Antarctic (putative 'Upland'; A. snaresensis) lineages. We present ancient-DNA analyses of prehistoric mainland NZ and sub-Antarctic fur seals, revealing that both of these genetic lineages were already widely distributed across the region at the time of human arrival. These findings indicate that anthropogenic factors did not contribute to the admixture of these lineages, and cast doubt on the validity of the Upland seal. Human-mediated impacts on Arctocephalus genetic diversity are instead highlighted by a dramatic temporal haplotype frequency-shift due to genetic drift in heavily bottlenecked populations following the cessation of industrial-scale harvesting. These extinction-recolonisation dynamics add to a growing picture of human-mediated change in NZ's coastal and marine ecosystems.

Low Spatial Genetic Differentiation Associated with Rapid Recolonization in the New Zealand Fur Seal Arctocephalus forsteri.

Population declines resulting from anthropogenic activities are of major consequence for the long-term survival of species because the resulting loss of genetic diversity can lead to extinction via the effects of inbreeding depression, fixation of deleterious mutations, and loss of adaptive potential. Otariid pinnipeds have been exploited commercially to near extinction with some species showing higher demographic resilience and recolonization potential than others. The New Zealand fur seal (NZFS) was heavily impacted by commercial sealing between the late 18th and early 19th centuries, but has recolonized its former range in southern Australia. The species has also recolonized its former range in New Zealand, yet little is known about the pattern of recolonization. Here, we first used 11 microsatellite markers (n = 383) to investigate the contemporary population structure and dispersal patterns in the NZFS (Arctocephalus forsteri). Secondly, we model postsealing recolonization with 1 additional mtDNA cytochrome b (n = 261) marker. Our data identified 3 genetic clusters: an Australian, a subantarctic, and a New Zealand one, with a weak and probably transient subdivision within the latter cluster. Demographic history scenarios supported a recolonization of the New Zealand coastline from remote west coast colonies, which is consistent with contemporary gene flow and with the species' high resilience. The present data suggest the management of distinct genetic units in the North and South of New Zealand along a genetic gradient. Assignment of individuals to their colony of origin was limited (32%) with the present data indicating the current microsatellite markers are unlikely sufficient to assign fisheries bycatch of NZFSs to colonies.

Little evidence of inbreeding depression for birth mass, survival and growth in Antarctic fur seal pups.

Inbreeding depression, the loss of offspring fitness due to consanguineous mating, is generally detrimental for individual performance and population viability. We investigated inbreeding effects in a declining population of Antarctic fur seals (Arctocephalus gazella) at Bird Island, South Georgia. Here, localised warming has reduced the availability of the seal's staple diet, Antarctic krill, leading to a temporal increase in the strength of selection against inbred offspring, which are increasingly failing to recruit into the adult breeding population. However, it remains unclear whether selection operates before or after nutritional independence at weaning. We therefore used microsatellite data from 885 pups and their mothers, and SNP array data from 98 mother-offspring pairs, to quantify the effects of individual and maternal inbreeding on three important neonatal fitness traits: birth mass, survival and growth. We did not find any clear or consistent effects of offspring or maternal inbreeding on any of these traits. This suggests that selection filters inbred individuals out of the population as juveniles during the time window between weaning and recruitment. Our study brings into focus a poorly understood life-history stage and emphasises the importance of understanding the ecology and threats facing juvenile pinnipeds.

Transcriptome of the dead: characterisation of immune genes and marker development from necropsy samples in a free-ranging marine mammal.

BACKGROUND: Transcriptomes are powerful resources, providing a window on the expressed portion of the genome that can be generated rapidly and at low cost for virtually any organism. However, because many genes have tissue-specific expression patterns, developing a complete transcriptome usually requires a 'discovery pool' of individuals to be sacrificed in order to harvest mRNA from as many different types of tissue as possible. This hinders transcriptome development in large, charismatic and endangered species, many of which stand the most to gain from such approaches. To circumvent this problem in a model pinniped species, we 454 sequenced cDNA from testis, heart, spleen, intestine, kidney and lung tissues obtained from nine adult male Antarctic fur seals (Arctocephalus gazella) that died of natural causes at Bird Island, South Georgia. RESULTS: After applying stringent quality control criteria based on length and annotation, we obtained 12,397 contigs which, in combination with 454 data previously obtained from skin, gave a total of 23,096 unique contigs. Homology was found to 77.0% of dog (Canis lupus familiaris) transcripts, suggesting that the combined assembly represents a substantial proportion of this species' transcriptome. Moreover, only 0.5% of transcripts revealed sequence similarity to bacteria, implying minimal contamination, and the percentage of transcripts involved in cell death was low at 2.6%. Transcripts with immune-related annotations were almost five-fold enriched relative to skin and represented 13.2% of all spleen-specific contigs. By reference to the dog, we also identified transcripts revealing homology to five class I, ten class II and three class III genes of the Major Histocompatibility Complex and derived the putative genomic distribution of 17,121 contigs, 2,119 in silico mined microsatellites and 9,382 single nucleotide polymorphisms. CONCLUSIONS: Our findings suggest that transcriptome development based on samples collected post mortem may greatly facilitate genomic studies, not only of marine mammals but also more generally of species that are of conservation concern.

Female fur seals show active choice for males that are heterozygous and unrelated.

Much debate surrounds the exact rules that influence mating behaviour, and in particular the selective forces that explain the evolution of female preferences. A key example is the lek paradox, in which female choice is expected rapidly to become ineffective owing to loss of additive genetic variability for the preferred traits. Here we exploit a remarkable system in which female fur seals exert choice by moving across a crowded breeding colony to visit largely static males. We show that females move further to maximize the balance between male high multilocus heterozygosity and low relatedness. Such a system shows that female choice can be important even in a strongly polygynous species, and at the same time may help to resolve the lek paradox because heterozygosity has low heritability and inbreeding avoidance means there is no single 'best' male for all females.

Genomic evidence for homoploid hybrid speciation in a marine mammal apex predator.

Hybridization is widespread and constitutes an important source of genetic variability and evolution. In animals, its role in generating novel and independent lineages (hybrid speciation) has been strongly debated, with only a few cases supported by genomic data. The South American fur seal (SAfs) Arctocephalus australis is a marine apex predator of Pacific and Atlantic waters, with a disjunct set of populations in Peru and Northern Chile [Peruvian fur seal (Pfs)] with controversial taxonomic status. We demonstrate, using complete genome and reduced representation sequencing, that the Pfs is a genetically distinct species with an admixed genome that originated from hybridization between the SAfs and the Galapagos fur seal (Arctocephalus galapagoensis) ~400,000 years ago. Our results strongly support the origin of Pfs by homoploid hybrid speciation over alternative introgression scenarios. This study highlights the role of hybridization in promoting species-level biodiversity in large vertebrates.

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 with their mitochondrial homologs to characterize synonymous to nonsynonymous substitution rates, which can shed light on the selection pressures these genes have been subjected to. Here, we characterize 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 noncanonical 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.

Demographic Reconstruction of Antarctic Fur Seals Supports the Krill Surplus Hypothesis.

Much debate surrounds the importance of top-down and bottom-up effects in the Southern Ocean, where the harvesting of over two million whales in the mid twentieth century is thought to have produced a massive surplus of Antarctic krill. This excess of krill may have allowed populations of other predators, such as seals and penguins, to increase, a top-down hypothesis known as the 'krill surplus hypothesis'. However, a lack of pre-whaling population baselines has made it challenging to investigate historical changes in the abundance of the major krill predators in relation to whaling. Therefore, we used reduced representation sequencing and a coalescent-based maximum composite likelihood approach to reconstruct the recent demographic history of the Antarctic fur seal, a pinniped that was hunted to the brink of extinction by 18th and 19th century sealers. In line with the known history of this species, we found support for a demographic model that included a substantial reduction in population size around the time period of sealing. Furthermore, maximum likelihood estimates from this model suggest that the recovered, post-sealing population at South Georgia may have been around two times larger than the pre-sealing population. Our findings lend support to the krill surplus hypothesis and illustrate the potential of genomic approaches to shed light on long-standing questions in population biology.

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.

Bayesian inference of a historical bottleneck in a heavily exploited marine mammal.

Emerging Bayesian analytical approaches offer increasingly sophisticated means of reconstructing historical population dynamics from genetic data, but have been little applied to scenarios involving demographic bottlenecks. Consequently, we analysed a large mitochondrial and microsatellite dataset from the Antarctic fur seal Arctocephalus gazella, a species subjected to one of the most extreme examples of uncontrolled exploitation in history when it was reduced to the brink of extinction by the sealing industry during the late eighteenth and nineteenth centuries. Classical bottleneck tests, which exploit the fact that rare alleles are rapidly lost during demographic reduction, yielded ambiguous results. In contrast, a strong signal of recent demographic decline was detected using both Bayesian skyline plots and Approximate Bayesian Computation, the latter also allowing derivation of posterior parameter estimates that were remarkably consistent with historical observations. This was achieved using only contemporary samples, further emphasizing the potential of Bayesian approaches to address important problems in conservation and evolutionary biology.

Selection for increased body length in Subantarctic fur seals on Amsterdam Island.

The traditional distinction between ecological and evolutionary times is eroding, calling for tighter links between ecology and evolution. An example of such a brigde between the two disciplines is the so-called 'animal model', a methodology initially developed by animal breeders, which has become very popular among ecologists studying contemporary microevolution. Using a Bayesian multi-trait 'animal model', we investigated the quantitative genetics of body size, a fitness-related trait, in Subantarctic fur seals (Arctocephalus tropicalis) breeding on Amsterdam Island, Southern Ocean. Our approach jointly modelled the growth and selection processes at work in this population. Body length is heritable for both sexes, and females are under selection for increased body length in this population. We strongly suspect the peculiar ecological context of impoverished, suitable prey availability exacerbated by density-dependence phenomena to be an important selective agent on females breeding on Amsterdam Island.