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.

The oldest known fur seal.

The poorly known fossil record of fur seals and sea lions (Otariidae) does not reflect their current diversity and widespread abundance. This limited fossil record contrasts with the more complete fossil records of other pinnipeds such as walruses (Odobenidae). The oldest known otariids appear 5-6 Ma after the earliest odobenids, and the remarkably derived craniodental morphology of otariids offers few clues to their early evolutionary history and phylogenetic affinities among pinnipeds. We report a new otariid, Eotaria crypta, from the lower middle Miocene 'Topanga' Formation (15-17.1 Ma) of southern California, represented by a partial mandible with well-preserved dentition. Eotaria crypta is geochronologically intermediate between 'enaliarctine' stem pinnipedimorphs (16.6-27 Ma) and previously described otariid fossils (7.3-12.5 Ma), as well as morphologically intermediate by retaining an M2 and a reduced M1 metaconid cusp and lacking P2-4 metaconid cusps. Eotaria crypta eliminates the otariid ghost lineage and confirms that otariids evolved from an 'enaliarctine'-like ancestor.

The genetic legacy of extreme exploitation in a polar vertebrate.

Understanding the effects of human exploitation on the genetic composition of wild populations is important for predicting species persistence and adaptive potential. We therefore investigated the genetic legacy of large-scale commercial harvesting by reconstructing, on a global scale, the recent demographic history of the Antarctic fur seal (Arctocephalus gazella), a species that was hunted to the brink of extinction by 18th and 19th century sealers. Molecular genetic data from over 2,000 individuals sampled from all eight major breeding locations across the species' circumpolar geographic distribution, show that at least four relict populations around Antarctica survived commercial hunting. Coalescent simulations suggest that all of these populations experienced severe bottlenecks down to effective population sizes of around 150-200. Nevertheless, comparably high levels of neutral genetic variability were retained as these declines are unlikely to have been strong enough to deplete allelic richness by more than around 15%. These findings suggest that even dramatic short-term declines need not necessarily result in major losses of diversity, and explain the apparent contradiction between the high genetic diversity of this species and its extreme exploitation history.

GCalignR: An R package for aligning gas-chromatography data for ecological and evolutionary studies.

Chemical cues are arguably the most fundamental means of animal communication and play an important role in mate choice and kin recognition. Consequently, there is growing interest in the use of gas chromatography (GC) to investigate the chemical basis of eco-evolutionary interactions. Both GC-MS (mass spectrometry) and FID (flame ionization detection) are commonly used to characterise the chemical composition of biological samples such as skin swabs. The resulting chromatograms comprise peaks that are separated according to their retention times and which represent different substances. Across chromatograms of different samples, homologous substances are expected to elute at similar retention times. However, random and often unavoidable experimental variation introduces noise, making the alignment of homologous peaks challenging, particularly with GC-FID data where mass spectral data are lacking. Here we present GCalignR, a user-friendly R package for aligning GC-FID data based on retention times. The package was developed specifically for ecological and evolutionary studies that seek to investigate similarity patterns across multiple and often highly variable biological samples, for example representing different sexes, age classes or reproductive stages. The package also implements dynamic visualisations to facilitate inspection and fine-tuning of the resulting alignments and can be integrated within a broader workflow in R to facilitate downstream multivariate analyses. We demonstrate an example workflow using empirical data from Antarctic fur seals and explore the impact of user-defined parameter values by calculating alignment error rates for multiple datasets. The resulting alignments had low error rates for most of the explored parameter space and we could also show that GCalignR performed equally well or better than other available software. We hope that GCalignR will help to simplify the processing of chemical datasets and improve the standardization and reproducibility of chemical analyses in studies of animal chemical communication and related fields.

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.

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

Custom genotyping arrays provide a flexible and accurate means of genotyping single nucleotide polymorphisms (SNPs) in a large number of individuals of essentially any organism. However, validation rates, defined as the proportion of putative SNPs that are verified to be polymorphic in a population, are often very low. A number of potential causes of assay failure have been identified, but none have been explored systematically. In particular, as SNPs are often developed from transcriptomes, parameters relating to the genomic context are rarely taken into account. Here, we assembled a draft Antarctic fur seal (Arctocephalus gazella) genome (assembly size: 2.41 Gb; scaffold/contig N50 : 3.1 Mb/27.5 kb). We then used this resource to map the probe sequences of 144 putative SNPs genotyped in 480 individuals. The number of probe-to-genome mappings and alignment length together explained almost a third of the variation in validation success, indicating that sequence uniqueness and proximity to intron-exon boundaries play an important role. The same pattern was found after mapping the probe sequences to the Walrus and Weddell seal genomes, suggesting that the genomes of species divergent by as much as 23 million years can hold information relevant to SNP validation outcomes. Additionally, reanalysis of genotyping data from seven previous studies found the same two variables to be significantly associated with SNP validation success across a variety of taxa. Finally, our study reveals considerable scope for validation rates to be improved, either by simply filtering for SNPs whose flanking sequences align uniquely and completely to a reference genome, or through predictive modelling.

Complete mitochondrial genome of the stoat (Mustela erminea) and New Zealand fur seal (Arctocephalus forsteri) and their significance for mammalian phylogeny.

The complete mitochondrial genome of three mustelid species, stoats (Mustela erminea), weasels (Mustela nivalis) and ferrets (Mustela furo), and the New Zealand fur seal (Arctocephalus forsteri) were sequenced using direct mitochondrial DNA extraction and overlapping long PCRs. The usual 37 mammalian mitochondrial genes (13 protein coding genes, 22 t-RNA and 2 r-RNA) were identified in all four mitogenomes. The divergence of stoats from other members of the sub-family Mustelinae was dated 4.5 million years ago. The mitogenomic data were consistent with a bear-like origin of seals.

Rates of assay success and genotyping error when single nucleotide polymorphism genotyping in non-model organisms: a case study in the Antarctic fur seal.

Although single nucleotide polymorphisms (SNPs) are increasingly being recognized as powerful molecular markers, their application to non-model organisms can bring significant challenges. Among these are imperfect conversion rates of assays designed from in silico resources and the enhanced potential for genotyping error relative to pre-validated, highly optimized human SNPs. To explore these issues, we used Illumina's GoldenGate assay to genotype 480 Antarctic fur seal (Arctocephalus gazella) individuals at 144 putative SNPs derived from a 454 transcriptome assembly. One hundred and thirty-five polymorphic SNPs (93.8%) were automatically validated by the program GenomeStudio, and the initial genotyping error rate, estimated from nine replicate samples, was 0.004 per reaction. However, an almost tenfold further reduction in the error rate was achieved by excluding 31 loci (21.5%) that exhibited unclear clustering patterns, manually editing clusters to allow rescoring of ambiguous or incorrect genotypes, and excluding 18 samples (3.8%) with unreliable genotypes. After stringent quality filtering, we also found a counter-intuitive negative relationship between in silico minor allele frequency and the conversion rate, suggesting that some of our assays may have been designed from paralogous loci. Nevertheless, we obtained over 45 000 individual SNP genotypes with a final error rate of 0.0005, indicating that the GoldenGate assay is eminently capable of generating large, high-quality data sets for non-model organisms. This has positive implications for future studies of the evolutionary, behavioural and conservation genetics of natural populations.

Differentiating between Steller sea lion (Eumetopias jubatus) and northern fur seal (Callorhinus ursinus) scats through analysis of faecal DNA.

We describe a method to determine the species of pinniped from faeces collected from sympatric Steller sea lion (Eumetopias jubatus) and northern fur seal (Callorhinus ursinus) rookeries using newly developed species-specific primers that amplify a 667-669-base pair segment from the mitochondrial DNA (mtDNA) cytochrome B (cytB) gene region. The primers yielded the correct species in 100% of tissue samples from 10 known animals and 100% of faecal samples from 13 known animals. Species could be identified unequivocally for 87.7% of faecal samples from 122 unknown individuals. The ability to differentiate between scats of sympatrically breeding Steller sea lions and northern fur seals will contribute to the range-wide knowledge of the foraging strategies of both species as well as allow researchers to examine the niche partitioning and potential resource competition between the two predators.

Correlation between male social status, testosterone levels, and parasitism in a dimorphic polygynous mammal.

Life history trade-offs have often been assumed to be the consequence of restrictions in the availability of critical resources such as energy and nutrients, which necessitate the differential allocation of resources to costly traits. Here, we examined endocrine (testosterone) and health (parasite burdens) parameters in territorial and non-territorial New Zealand fur seal males. We documented intra-sexual differences in sexual behaviours, testosterone levels, and parasitism that suggest a trade-off exists between reproductive success and physical health, particularly susceptibility to helminths and acanthocephalans, in males displaying different mating tactics (i.e., territorial and non-territorial tactics). Levels of testosterone were higher in territorial males and correlated positively with reproductive effort (i.e., intra- and inter-sexual interactions). However, these territorial males also exhibited high levels of parasitic infection, which may impair survival in the long-term. Our study, while limited in sample size, provides preliminary evidence for a link between male mating tactics, testosterone levels and parasite loads, and potential effects on reproductive success and life history that should be explored further.

The monophyletic origin of sea lions and fur seals (Carnivora; Otariidae) in the Southern Hemisphere.

The pinniped family Otariidae (sea lions and fur seals) is composed of 7 extant genera with 14 species. They are mainly distributed in the Southern Hemisphere, but the fossil record is only known from the Northern Hemisphere until Pliocene. To clarify the biological and zoogeographical events during their evolution, it is necessary to reconstruct a robust phylogenetic tree. However, phylogenetic relationships among otariids continue to be controversial, except for the basal position of the northern fur seal among the extant otariids. We reconstructed phylogenetic trees of otariids based on mitochondrial genomes and multiple nuclear genes (IRBP and type I STS markers). The monophyly of the otariids including both sea lions and fur seals in the Southern Hemisphere was strongly supported by both the mitochondrial and nuclear evidence. We propose a novel evolutionary and dispersal scenario of otariids based on this phylogenetic hypothesis, estimated divergence times, and fossil records. According to our results, the center of origin of the southern otariids is hypothesized to be the eastern South Pacific along the west coast of South America.

Clareamento da melanina na epiderme do lobo-marinho-sul-americano e sua aplicação na imunohistoquímica enzimática / Bleaching of melanin in the epidermis of south American fur seal and its application on enzyme immunohistochemistry

O Lobo-marinho-sul-americano (Arctocephalus australis) é um mamífero marinho anfíbio distribuído ao longo da Costa do Atlântico e do Pacífico da América do Sul. Esta espécie está bem adaptada a diferentes habitats devido à morfologia dos membros em forma de nadadeira e de seu sistema tegumentar. Estudos imuno-histoquímicos são importantes para avaliar os mecanismos de adaptação da pele devido a diferencial expressão dos antígenos presentes no tecido dependendo da região da superfície corporal. Entretanto, sua epiderme altamente pigmentada (melanina) impede a visualização dos marcadores cromógenos utilizados na imunohistoquímica. Neste trabalho foi desenvolvido um método de clarear a melanina para permitir a visualização dos cromógenos sem alterar a afinidade antígeno-anticorpo para a imuno-histoquímica. A análise do índice do PCNA (proliferating cell nuclear antigen) na epiderme de A. australis, com diaminobenzidina (DAB) como cromógeno foi usada para testar o método. O clareamento da melanina permitiu obter o índice de proliferação celular na epiderme e evitar resultados falso-positivos sem afetar os resultados imuno-histoquímicos.

The South American fur seal (Arctocephalus australis) is an amphibious marine mammal distributed along the Atlantic and Pacific coasts of South America. The species is well adjusted to different habitats due to the morphology of its fin-like members and due to some adaptations in their integumentary system. Immunohistochemical studies are very important to evaluate the mechanisms of skin adaptation due the differential expression of the antigens present in the tissue depending of the region of the body surface. However, its strongly pigmented (melanin) epidermis prevents the visualization of the immunohistochemical chromogens markers. In this study a melanin bleaching method was developed aimed to allow the visualization of the chromogens without interfering in the antigen-antibody affinity for immunohistochemistry. The analysis of PCNA (proliferating cell nuclear antigen) index in the epidermis of A. australis by immunohistochemistry with diaminobenzidine (DAB) as chromogen was used to test the method. The bleaching of the melanin allowed to obtain the cell proliferation index in epidermis and to avoid false positive results without affecting the immunohistochemical results.

Ménage à trois on Macquarie Island: hybridization among three species of fur seal (Arctocephalus spp.) following historical population extinction.

Human-induced changes to natural systems can cause major disturbances to fundamental ecological and population processes and result in local extinctions and secondary contacts between formerly isolated species. Extensive fur seal harvesting during the nineteenth century on Macquarie Island (subantarctic) resulted in extinction of the original population. Recolonization by three species has been slow and complex, characterized by the establishment of breeding groups of Antarctic and subantarctic fur seals (Arctocephalus gazella and Arctocephalus tropicalis) and presumed nonbreeding (itinerant) male New Zealand fur seals (Arctocephalus forsteri). One thousand and seven pups from eight annual cohorts (1992-2003) were analysed using mitochondrial control region data (RFLP) and 10 microsatellite loci to estimate species composition and hybridization. Antarctic fur seals predominated, but hybridization occurred between all three species (17-30% of all pups). Involvement of New Zealand fur seals was unexpected as females are absent and males are not observed to hold territories during the breeding season. The proportion of hybrids in the population has fallen over time, apparently owing to substantial influxes of pure Antarctic and subantarctic individuals and non-random mating. Over 50% of New Zealand hybrids and 43% of Antarctic-subantarctic hybrids were not F(1), which indicates some degree of hybrid reproductive success, and this may be underestimated: simulations showed that hybrids become virtually undetectable by the third generation of backcrossing. While human impacts seem to have driven novel hybridization in this population, the present 'time slices' analysis suggests some biological resistance to complete homogenization.

Genetic variability in Guadalupe fur seals.

The Guadalupe fur seal (Arctocephalus townsendi) population underwent one or two severe bottlenecks due to commercial sealing in the late 19th century. Since then the protected population has been growing steadily around their only rookery, Isla de Guadalupe, Mexico. We probed both nuclear and mitochondrial genomes using multilocus nuclear DNA profiling and mitochondrial DNA sequencing to estimate the level of genetic variability of the present population. Unlike other pinniped populations that have experienced similar historical bottlenecks, such as Hawaiian monk seals and northern elephant seals, high levels of genetic variability were found.

Phylogenetic relationships within the eared seals (Otariidae: Carnivora): implications for the historical biogeography of the family.

Phylogenetic relationships within the family Otariidae were investigated using two regions of the mitochondrial genome. A 360-bp region of the cytochrome b gene was employed for the primary phylogenetic analysis, while a 356-bp segment of the control region was used to enhance resolution of the terminal nodes. Traditional classification of the family into the subfamilies Arctocephalinae (fur seals) and Otariinae (sea lions) is not supported, with the fur seal Callorhinus ursinus having a basal relationship relative to the rest of the family. This is consistent with the fossil record which suggests that this genus diverged from the line leading to the remaining fur seals and sea lions about 6 million years ago (mya). There is also little evidence to support or refute the monophyly of sea lions. Four sea lion clades and five fur seal clades were observed, but relationships among these clades are unclear. Similar genetic divergences between the sea lion clades (D(a) = 0.054-0.078), as well as between the major Arctocephalus fur seal clades (D(a) = 0.040-0.069) suggest that these groups underwent periods of rapid radiation at about the time they diverged from each other. Rapid radiations of this type make the resolution of relationships between the resulting species difficult and indicate the requirement for additional molecular data from both nuclear and mitochondrial genes. The phylogenetic relationships within the family and the genetic distances among some taxa highlight inconsistencies in the current taxonomic classification of the family.