Terrestrial Birth and Body Size Tune UCP1 Functionality in Seals.

The molecular evolution of the mammalian heater protein UCP1 is a powerful biomarker to understand thermoregulatory strategies during species radiation into extreme climates, such as aquatic life with high thermal conductivity. While fully aquatic mammals lost UCP1, most semiaquatic seals display intact UCP1 genes, apart from large elephant seals. Here, we show that UCP1 thermogenic activity of the small-bodied harbor seal is equally potent compared to terrestrial orthologs, emphasizing its importance for neonatal survival on land. In contrast, elephant seal UCP1 does not display thermogenic activity, not even when translating a repaired or a recently highlighted truncated version. Thus, the thermogenic benefits for neonatal survival during terrestrial birth in semiaquatic pinnipeds maintained evolutionary selection pressure on UCP1 function and were only outweighed by extreme body sizes among elephant seals, fully eliminating UCP1-dependent thermogenesis.

Variations in antimicrobial resistance genes present in the rectal faeces of seals in Scottish and Liverpool Bay coastal waters.

Antibiotic resistance genes originating from human activity are considered important environmental pollutants. Wildlife species can act as sentinels for coastal environmental contamination and in this study we used qPCR array technology to investigate the variety and abundance of antimicrobial resistance genes (ARGs), mobile genetic elements (MGEs) and integrons circulating within seal populations both near to and far from large human populations located around the Scottish and northwest English coast. Rectal swabs were taken from 50 live grey seals and nine live harbour seals. Nucleic acids were stabilised upon collection, enabling extraction of sufficient quality and quantity DNA for downstream analysis. 78 ARG targets, including genes of clinical significance, four MGE targets and three integron targets were used to monitor genes within 22 sample pools. 30 ARGs were detected, as well as the integrons intl1 and intl2 and tnpA transposase. Four ß-lactam, nine tetracycline, two phenicol, one trimethoprim, three aminoglycoside and ten multidrug resistance genes were detected as well as mcr-1 which confers resistance to colistin, an important drug of last resort. No sulphonamide, vancomycin, macrolide, lincosamide or streptogramin B (MLSB) resistance genes were detected. Resistance genes were detected in all sites but the highest number of ARGs (n = 29) was detected in samples derived from grey seals on the Isle of May, Scotland during the breeding season, and these genes also had the highest average abundance in relation to the 16S rRNA gene. This pilot study demonstrates the effectiveness of a culture-independent workflow for global analysis of ARGs within the microbiota of live, free-ranging, wild animals from habitats close to and remote from human habitation, and highlights seals as a valuable indicator species for monitoring the presence, abundance and land-sea transference of resistance genes within and between ecosystems.

Population genomics of harbour seal Phoca vitulina from northern British Columbia through California and comparison to the Atlantic subspecies.

The harbour seal Phoca vitulina is a ubiquitous pinniped species found throughout coastal waters of the Northern Hemisphere. Harbour seal impacts on ecosystem dynamics may be significant due to their high abundance and food web position. Two subspecies exist in North America, P. v. richardii in the Pacific Ocean and P. v. vitulina in the Atlantic. Strong natal philopatry of harbour seals can result in fine-scale genetic structure and isolation by distance. Management of harbour seals is expected to benefit from improved resolution of seal population structure and dynamics. Here, we use genotyping-by-sequencing to genotype 146 harbour seals from the eastern Pacific Ocean (i.e. British Columbia (BC), Oregon and California) and the western Atlantic Ocean (i.e. Québec, Newfoundland and Labrador). Using 12,742 identified variants, we confirm the recently identified elevated genetic diversity in the eastern Pacific relative to the western Atlantic and greatest differentiation between the subspecies. Further, we demonstrate that this is independent of reference genome bias or other potential technical artefacts. Coast-specific analyses with 8933 and 3828 variants in Pacific and Atlantic subspecies, respectively, identify divergence between BC and Oregon-California, and between Québec and Newfoundland-Labrador. Unexpected PCA outlier clusters were observed in two populations due to cryptic relatedness of individuals; subsequently, closely related samples were removed. Admixture analysis indicates an isolation-by-distance signature where Oregon seals contained some of the BC signature, whereas California did not. Additional sampling is needed in the central and north coast of BC to determine whether a discrete separation of populations exists within the region.

Sequence Diversity and Differences at the Highly Duplicated MHC-I Gene Reflect Viral Susceptibility in Sympatric Pinniped Species.

Differences in disease susceptibility among species can result from rapid host-pathogen coevolution and differences in host species ecology that affect the strength and direction of natural selection. Among 2 sympatric pinniped species that differ in sociality and putative disease exposure, we investigate observed differences in susceptibility through an analysis of a highly variable, duplicated gene family involved in the vertebrate immune response. Using high-throughput amplicon sequencing, we characterize diversity at the 2 exons that encode the peptide binding region of the major histocompatibility complex class I (MHC-I) gene in harbor (N = 60) and gray (N = 90) seal populations from the Northwest Atlantic. Across species, we identified 106 full-length exon 2 and 103 exon 3 sequence variants and a minimum of 11 duplicated MHC-I loci. The sequence variants clustered in 15 supertypes defined by the physiochemical properties of the peptide binding region, including a putatively novel Northwest Atlantic MHC-I diversity sublineage. Trans-species polymorphisms, dN/dS ratios, and evidence of gene conversion among supertypes are consistent with balancing selection acting on this gene. High functional redundancy suggests particularly strong selection among gray seals at the novel Northwest Atlantic MHC-I diversity sublineage. At exon 2, harbor seals had a significantly greater number of variants per individual than gray seals, but fewer supertypes. Supertype richness and private supertypes are hypothesized to contribute to observed differences in disease resistance between species, as consistently, across the North Atlantic and many disease outbreaks, gray seals appear to be more resistant to respiratory viruses than harbor seals.

Origin and expansion of the world's most widespread pinniped: Range-wide population genomics of the harbour seal (Phoca vitulina).

The harbour seal (Phoca vitulina) is the most widely distributed pinniped, occupying a wide variety of habitats and climatic zones across the Northern Hemisphere. Intriguingly, the harbour seal is also one of the most philopatric seals, raising questions as to how it colonized its current range. To shed light on the origin, remarkable range expansion, population structure and genetic diversity of this species, we used genotyping-by-sequencing to analyse ~13,500 biallelic single nucleotide polymorphisms from 286 individuals sampled from 22 localities across the species' range. Our results point to a Northeast Pacific origin of the harbour seal, colonization of the North Atlantic via the Canadian Arctic, and subsequent stepping-stone range expansions across the North Atlantic from North America to Europe, accompanied by a successive loss of genetic diversity. Our analyses further revealed a deep divergence between modern North Pacific and North Atlantic harbour seals, with finer-scale genetic structure at regional and local scales consistent with strong philopatry. The study provides new insights into the harbour seal's remarkable ability to colonize and adapt to a wide range of habitats. Furthermore, it has implications for current harbour seal subspecies delineations and highlights the need for international and national red lists and management plans to ensure the protection of genetically and demographically isolated populations.

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.

Mitochondrial DNA reveals secondary contact in Japanese harbour seals, the southernmost population in the western Pacific.

In this study, we used relatively large number of samples (n = 178) and control region of mtDNA (454bp) to clearify the divergence history of Japanese harbour seals (Phoca vitulina stejnegeri) and phylogenetic relationship between the seals in Japan and other countries. Our results suggested that Japanese harbour seals possibly consisted of more than two lineages and secondary contact of populations after a long isolation. Furthermore, one of the lineage was made only by Japanese harbour seals (Group P1). The proportion of Group P1 was the highest at the South West and gradually decreased towards the North East of Hokkaido, Japan. On the other hand, the haplotypes do not belonged to Group P1 showed close relationship to the seals in the North Pacific. Based on the fossil record of harbour seal in Japan and the range of sea ice during the Last Glacial Maximum (LGM), Group P1 might have entered Japan before the LGM and became isolated due to the geographical boundary, and gradually extended its range from the South West towards the North East of Hokkaido after the disappearance of the sea ice, while the seals which are not in Group P1 immigrated into Japan from the North Pacific.

wisepair: a computer program for individual matching in genetic tracking studies.

Individual-based data sets tracking organisms over space and time are fundamental to answering broad questions in ecology and evolution. A 'permanent' genetic tag circumvents a need to invasively mark or tag animals, especially if there are little phenotypic differences among individuals. However, genetic tracking of individuals does not come without its limits; correctly matching genotypes and error rates associated with laboratory work can make it difficult to parse out matched individuals. In addition, defining a sampling design that effectively matches individuals in the wild can be a challenge for researchers. Here, we combine the two objectives of defining sampling design and reducing genotyping error through an efficient Python-based computer-modelling program, wisepair. We describe the methods used to develop the computer program and assess its effectiveness through three empirical data sets, with and without reference genotypes. Our results show that wisepair outperformed similar genotype matching programs using previously published from reference genotype data of diurnal poison frogs (Allobates femoralis) and without-reference (faecal) genotype sample data sets of harbour seals (Phoca vitulina) and Eurasian otters (Lutra lutra). In addition, due to limited sampling effort in the harbour seal data, we present optimal sampling designs for future projects. wisepair allows for minimal sacrifice in the available methods as it incorporates sample rerun error data, allelic pairwise comparisons and probabilistic simulations to determine matching thresholds. Our program is the lone tool available to researchers to define parameters a priori for genetic tracking studies.

Structural characterisation of Toll-like receptor 1 (TLR1) and Toll-like receptor 6 (TLR6) in elephant and harbor seals.

Pinnipeds are a diverse clade of semi-aquatic mammals, which act as key indicators of ecosystem health. Their transition from land to marine environments provides a complex microbial milieu, making them vulnerable to both aquatic and terrestrial pathogens, thereby contributing to pinniped population decline. Indeed, viral pathogens such as influenza A virus and phocine distemper virus (PDV) have been identified as the cause of several of these mass mortality events. Furthermore, bacterial infection with mammalian Brucella sp. and methicillin-resistant Staphylococcus aureus strains have also been observed in marine mammals, posing further risk to both co-habiting endangered species and public health. During these disease outbreaks, mortality rates have varied amongst different pinniped species. Analyses of innate immune receptors at the host-pathogen interface have previously identified variants which may drive these species-specific responses. Through a combination of both sequence- and structure-based methods, this study characterises members of the Toll-like receptor (TLR) 1 superfamily from both harbour and elephant seals, identifying variations which will help us to understand these species-specific innate immune responses, potentially aiding the development of specific vaccine-adjuvants for these species.

Brain Meta-Transcriptomics from Harbor Seals to Infer the Role of the Microbiome and Virome in a Stranding Event.

Marine diseases are becoming more frequent, and tools for identifying pathogens and disease reservoirs are needed to help prevent and mitigate epizootics. Meta-transcriptomics provides insights into disease etiology by cataloguing and comparing sequences from suspected pathogens. This method is a powerful approach to simultaneously evaluate both the viral and bacterial communities, but few studies have applied this technique in marine systems. In 2009 seven harbor seals, Phoca vitulina, stranded along the California coast from a similar brain disease of unknown cause of death (UCD). We evaluated the differences between the virome and microbiome of UCDs and harbor seals with known causes of death. Here we determined that UCD stranded animals had no viruses in their brain tissue. However, in the bacterial community, we identified Burkholderia and Coxiella burnetii as important pathogens associated with this stranding event. Burkholderia were 100% prevalent and ~2.8 log2 fold more abundant in the UCD animals. Further, while C. burnetii was found in only 35.7% of all samples, it was highly abundant (~94% of the total microbial community) in a single individual. In this harbor seal, C. burnetii showed high transcription rates of invading and translation genes, implicating it in the pathogenesis of this animal. Based on these data we propose that Burkholderia taxa and C. burnetii are potentially important opportunistic neurotropic pathogens in UCD stranded harbor seals.

Uridine Diphosphate-Glucuronosyltransferase (UGT) Xenobiotic Metabolizing Activity and Genetic Evolution in Pinniped Species.

There are various interspecies differences in xenobiotic-metabolizing enzymes. It is known that cats show slow glucuronidation of drugs such as acetaminophen and strong side effects due to the UGT1A6 pseudogene. Recently, the UGT1A6 pseudogene was found in the Northern elephant seal and Otariidae was suggested to be UGT1A6-deficient. From the results of measurements of uridine diphosphate-glucuronosyltransferase (UGT) activity using liver microsomes, the Steller sea lion, Northern fur seal, and Caspian seal showed UGT activity toward 1-hydroxypyrene and acetaminophen as low as in cats, which was significantly lower than in rat and dog. Furthermore, UGT1A6 pseudogenes were found in Steller sea lion and Northern fur seal, and all Otariidae species were suggested to have the UGT1A6 pseudogene. The UGT1 family genes appear to have undergone birth-and-death evolution based on a phylogenetic and synteny analysis of the UGT1 family in mammals including Carnivora. UGT1A2-1A5 and UGT1A7-1A10 are paralogous genes to UGT1A1 and UGTA6, respectively, and their numbers were lower in cat, ferret and Pacific walrus than in human, rat, and dog. Felidae and Pinnipedia, which are less exposed to natural xenobiotics such as plant-derived toxins due to their carnivorous diet, have experienced fewer gene duplications of xenobiotic-metabolizing UGT genes, and even possess UGT1A6 pseudogenes. Artificial environmental pollutants and drugs conjugated by UGT are increasing dramatically, and their elimination to the environment can be of great consequence to cat and Pinnipedia species, whose low xenobiotic glucuronidation capacity makes them highly sensitive to these compounds.

Sequence variation and gene duplication at the MHC DRB loci of the spotted seal Phoca largha.

The major histocompatibility complex (MHC) is one of the most important genetic systems associated with resistance to infectious diseases in vertebrates. The spotted seal (Phoca largha) is one of the most endangered species in China. In this study, we present the first step in the molecular characterization of a DRB-like locus in the spotted seal by analyzing the nucleotide sequence of the polymorphic exon 2 segments, a 288-nucleotide sequence. By examining the segment from a group of 41 individuals, 28 alleles were identified. No deletion, insertion, or exceptional stop codon was detected, suggesting that these alleles could be functional in vivo. The nucleotide and amino acid sequences of the segment both showed a relatively high level of similarity (nucleotides 97%; amino acids 98%) to those of Meles meles and Zalophus californianus. The high level of spotted seal MHC-DRB polymorphism revealed in the present study has not been reported for the Phocidae and could be a consequence of the small spotted seal population adapting to the Bohai Sea, which probably has a relatively high level of pathogens.

High-throughput sequencing reveals inbreeding depression in a natural population.

Proxy measures of genome-wide heterozygosity based on approximately 10 microsatellites have been used to uncover heterozygosity fitness correlations (HFCs) for a wealth of important fitness traits in natural populations. However, effect sizes are typically very small and the underlying mechanisms remain contentious, as a handful of markers usually provides little power to detect inbreeding. We therefore used restriction site associated DNA (RAD) sequencing to accurately estimate genome-wide heterozygosity, an approach transferrable to any organism. As a proof of concept, we first RAD sequenced oldfield mice (Peromyscus polionotus) from a known pedigree, finding strong concordance between the inbreeding coefficient and heterozygosity measured at 13,198 single-nucleotide polymorphisms (SNPs). When applied to a natural population of harbor seals (Phoca vitulina), a weak HFC for parasite infection based on 27 microsatellites strengthened considerably with 14,585 SNPs, the deviance explained by heterozygosity increasing almost fivefold to a remarkable 49%. These findings arguably provide the strongest evidence to date of an HFC being due to inbreeding depression in a natural population lacking a pedigree. They also suggest that under some circumstances heterozygosity may explain far more variation in fitness than previously envisaged.

Integrating genetic data and population viability analyses for the identification of harbour seal (Phoca vitulina) populations and management units.

Identification of populations and management units is an essential step in the study of natural systems. Still, there is limited consensus regarding how to define populations and management units, and whether genetic methods allow for inference at the relevant spatial and temporal scale. Here, we present a novel approach, integrating genetic, life history and demographic data to identify populations and management units in southern Scandinavian harbour seals. First, 15 microsatellite markers and model- and distance-based genetic clustering methods were used to determine the population genetic structure in harbour seals. Second, we used harbour seal demographic and life history data to conduct population viability analyses (PVAs) in the vortex simulation model in order to determine whether the inferred genetic units could be classified as management units according to Lowe and Allendorf's (Molecular Ecology, 19, 2010, 3038) 'population viability criterion' for demographic independence. The genetic analyses revealed fine-scale population structuring in southern Scandinavian harbour seals and pointed to the existence of several genetic units. The PVAs indicated that the census population size of each of these genetic units was sufficiently large for long-term population viability, and hence that the units could be classified as demographically independent management units. Our study suggests that population genetic inference can offer the same degree of temporal and spatial resolution as 'nongenetic' methods and that the combined use of genetic data and PVAs constitutes a promising approach for delineating populations and management units.

De novo assembly and characterization of spotted seal Phoca largha transcriptome using Illumina paired-end sequencing.

Spotted seal (Phoca largha) is categorized as a critically endangered species in China. The aim of this study was to investigate spotted seal transcriptome by the approach of Illumina paired-end sequencing technology. We obtained a total of 52,146,394 reads for the mixed tissues of liver and spleen from the spotted seal. The de novo assemblies yielded 354,014 contigs and 178,466 unigenes. In the transcriptome, 193 unigenes were assigned to defense mechanisms. Three unigenes encoded MHC class I and 17 unigenes encoded MHC class II. In addition, bioinformatics analysis revealed a total of 4425 simple sequence repeats (SSRs). Fifty SSRs were randomly selected to validate amplification and determine the degree of polymorphism in the genomic DNA pools. Thirty-five primer pairs successfully amplified the expected DNA fragments and detected significant polymorphism among 28 spotted seal individuals. These results would contribute to the understanding of the genetic makeup of spotted seal transcriptome and provide useful information for functional genomic research in this species.

Characterization of the spotted seal Phoca largha transcriptome using Illumina paired-end sequencing and development of SSR markers.

Next-generation sequencing provides a powerful new approach for developing functional genomic tools for non-model species. This study aims to investigate the spotted seal (Phoca largha) transcriptome by the approach of Illumina paired-end sequencing technology. We obtained a total of 52,146,394 reads for the mixed tissues of liver and spleen from spotted seal. The de novo assemblies yielded 354,014 contigs and 178,466 unigenes. In addition, bioinformatics analysis revealed a total of 4425 simple sequence repeats (SSRs). Fifty SSRs were randomly selected to validate amplification and determine the degree of polymorphism in the genomic DNA pools. Thirty-five primer pairs successfully amplified expected DNA fragments and detected significant polymorphism among 28 spotted seal individuals. These results contribute to the understanding of the genetic makeup of spotted seal transcriptome and provide useful information for functional genomic research in this species.

Variation in European harbour seal immune response genes and susceptibility to phocine distemper virus (PDV).

Phocine distemper virus (PDV) has caused two mass mortalities of European harbour seals (Phoca vitulina) in recent decades. Levels of mortality varied considerably among European populations in both the 1988 and 2002 epidemics, with higher mortality in continental European populations in comparison to UK populations. High levels of genetic differentiation at neutral makers among seal populations allow for the possibility that there could be potential genetic differences at functional loci that may account for some of the variation in mortality. Recent genome sequencing of carnivore species and development of genomic tools have now made it possible to explore the possible contribution of variation in candidate genes from harbour seals in relation to the differential mortality patterns. We assessed variation in eight genes (CD46, IFNG, IL4, IL8, IL10, RARa, SLAM and TLR2) encoding key proteins involved in host cellular interactions with Morbilliviruses and the relationship of variants to disease status. This work constitutes the first genetic association study for Morbillivirus disease susceptibility in a non-model organism, and for a natural mortality event. We found no variation in harbour seals from across Europe in the protein coding domains of the viral receptors SLAM and CD46, but SNPs were present in SLAM intron 2. SNPs were also present in IL8 p2 and RARa exon 1. There was no significant association of SLAM or RARa polymorphisms with disease status implying no role of these genes in determining resistance to PDV induced mortality, that could be detected with the available samples and the small number of polymorphisms indentified. However there was significant differentiation of allele frequencies among populations. PDV and other morbilliviruses are important models for wildlife epidemiology, host switches and viral evolution. Despite a negative result in this case, full sequencing of pinniped and other 'non-model' carnivore genomes will help in refining understanding the role of host genetics in disease susceptibility for these viruses.

Identification and phylogenetic analysis of novel cytochrome P450 1A genes from ungulate species.

As part of an ongoing effort to understand the biological response of wild and domestic ungulates to different environmental pollutants such as dioxin-like compounds, cDNAs encoding for CYP1A1 and CYP1A2 were cloned and characterized. Four novel CYP1A cDNA fragments from the livers of four wild ungulates (elephant, hippopotamus, tapir and deer) were identified. Three fragments from hippopotamus, tapir and deer were classified as CYP1A2, and the other fragment from elephant was designated as CYP1A1/2. The deduced amino acid sequences of these fragment CYP1As showed identities ranging from 76 to 97% with other animal CYP1As. The phylogenetic analysis of these fragments showed that both elephant and hippopotamus CYP1As made separate branches, while tapir and deer CYP1As were located beside that of horse and cattle respectively in the phylogenetic tree. Analysis of dN/dS ratio among the identified CYP1As indicated that odd toed ungulate CYP1A2s were exposed to different selection pressure.

A computational model of teeth and the developmental origins of morphological variation.

The relationship between the genotype and the phenotype, or the genotype-phenotype map, is generally approached with the tools of multivariate quantitative genetics and morphometrics. Whereas studies of development and mathematical models of development may offer new insights into the genotype-phenotype map, the challenge is to make them useful at the level of microevolution. Here we report a computational model of mammalian tooth development that combines parameters of genetic and cellular interactions to produce a three-dimensional tooth from a simple tooth primordia. We systematically tinkered with each of the model parameters to generate phenotypic variation and used geometric morphometric analyses to identify, or developmentally ordinate, parameters best explaining population-level variation of real teeth. To model the full range of developmentally possible morphologies, we used a population sample of ringed seals (Phoca hispida ladogensis). Seal dentitions show a high degree of variation, typically linked to the lack of exact occlusion. Our model suggests that despite the complexity of development and teeth, there may be a simple basis for dental variation. Changes in single parameters regulating signalling during cusp development may explain shape variation among individuals, whereas a parameter regulating epithelial growth may explain serial, tooth-to-tooth variation along the jaw. Our study provides a step towards integrating the genotype, development and the phenotype.

Genetic variation in the harbor seal (Phoca vitulina) and spotted seal (Phoca largha) around Hokkaido, Japan, based on mitochondrial cytochrome b sequences.

The harbor seal (Phoca vitulina) and spotted seal (Phoca largha) are the main seal species around Hokkaido, Japan. While some investigations have been conducted on the ecology and morphology of these two species, there is a lack of genetic information. We studied variation in mitochondrial cytochrome b sequences in the two species. Fifteen haplotypes were observed in 39 harbor seals from Erimo, Akkeshi, and Nosappu, and 23 were observed in 31 spotted seals from Erimo, Akkeshi, Nosappu, Rausu, Yagishiri Island, and Hamamasu. Phylogenetic trees showed two harbor seal lineages: Group I contained primarily haplotypes from Erimo, and Group II contained haplotypes from Akkeshi and Nosappu. Because the Erimo population had fewer haplotypes and less nucleotide diversity than the Akkeshi and Nosappu populations, we considered it to be Isolated from the others. In contrast, genetic variance within populations of spotted seals (97.3%) was far higher than that among populations (2.7%), determined by analysis of molecular variance. There were no significant difference among the spotted seal populations, indicating the absence of distinct lineages around Hokkaido. The differences in the genetic population structure between the two species could have been generated by their ecological differences. This study provides basic genetic information on these seal species and will contribute to the conservation and management of fisheries and seals throughout Hokkaido.