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

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

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

Inbreeding depresses altruism in a cooperative society

In some animal species, individuals regularly breed with relatives, including siblings and parents. Given the high fitness costs of inbreeding, evolutionary biologists have found it challenging to understand the persistence of these inbred societies in nature. One appealing but untested explanation is that early life care may create a benign environment that offsets inbreeding depression, allowing inbred societies to evolve. We test this possibility using 21 years of data from a wild cooperatively breeding mammal, the banded mongoose, a species where almost one in ten young result from close inbreeding. We show that care provided by parents and alloparents mitigates inbreeding depression for early survival. However, as adults, inbred individuals provide less care, reducing the amount of help available to the next generation. Our results suggest that inbred cooperative societies are rare in nature partly because the protective care that enables elevated levels of inbreeding can be reduced by inbreeding depression

Unexpected genetic differentiation between recently recolonized populations of a long‐lived and highly vagile marine mammal

Many species have been heavily exploited by man leading to local extirpations, yet few studies have attempted to unravel subsequent recolonization histories. This has led to a significant gap in our knowledge of the long-term effects of exploitation on the amount and structure of contemporary genetic variation, with important implications for conservation. The Antarctic fur seal provides an interesting case in point, having been virtually exterminated in the nineteenth century but subsequently staged a dramatic recovery to recolonize much of its original range. Consequently, we evaluated the hypothesis that South Georgia (SG), where a few million seals currently breed, was the main source of immigrants to other locations including Livingston Island (LI), by genotyping 366 individuals from these two populations at 17 microsatellite loci and sequencing a 263 bp fragment of the mitochondrial hypervariable region 1. Contrary to expectations, we found highly significant genetic differences at both types of marker, with 51% of LI individuals carrying haplotypes that were not observed in 246 animals from SG. Moreover, the youngest of three sequentially founded colonies at LI showed greater similarity to SG at mitochondrial DNA than microsatellites, implying temporal and sex-specific variation in recolonization. Our findings emphasize the importance of relict populations and provide insights into the mechanisms by which severely depleted populations can recover while maintaining surprisingly high levels of genetic diversity

Genome-wide insights into introgression and its consequences for genome-wide heterozygosity in the Mytilus species complex across Europe.

The three mussel species comprising the Mytilus complex are widespread across Europe and readily hybridize when they occur in sympatry, resulting in a mosaic of populations with varying genomic backgrounds. Two of these species, M. edulis and M. galloprovincialis, are extensively cultivated across Europe, with annual production exceeding 230,000 tonnes. The third species, M. trossulus, is considered commercially damaging as hybridization with this species results in weaker shells and poor meat quality. We therefore used restriction site associated DNA sequencing to generate high-resolution insights into the structure of the Mytilus complex across Europe and to resolve patterns of introgression. Inferred species distributions were concordant with the results of previous studies based on smaller numbers of genetic markers, with M. edulis and M. galloprovincialis predominating in northern and southern Europe respectively, while introgression between these species was most pronounced in northern France and the Shetland Islands. We also detected traces of M. trossulus ancestry in several northern European populations, especially around the Baltic and in northern Scotland. Finally, genome-wide heterozygosity, whether quantified at the population or individual level, was lowest in M. edulis, intermediate in M. galloprovincialis, and highest in M. trossulus, while introgression was positively associated with heterozygosity in M. edulis but negatively associated with heterozygosity in M. galloprovincialis. Our study will help to inform mussel aquaculture by providing baseline information on the genomic backgrounds of different Mytilus populations across Europe and by elucidating the effects of introgression on genome-wide heterozygosity, which is known to influence commercially important traits such as growth, viability, and fecundity in mussels

A high-quality pedigree and genetic markers both reveal inbreeding depression for quality but not survival in a cooperative mammal

Inbreeding depression, the reduced fitness of offspring of closely related parents, is commonplace in both captive and wild populations and has important consequences for conservation and mating system evolution. However, because of the difficulty of collecting pedigree and life history data from wild populations, relatively few studies have been able to compare inbreeding depression for traits at different points in the life cycle. Moreover, pedigrees give the expected proportion of the genome that is identical by descent (IBDg) whereas in theory with enough molecular markers realised IBDg can be quantified directly. We therefore investigated inbreeding depression for multiple life-history traits in a wild population of banded mongooses using pedigree-based inbreeding coefficients (fped) and standardised multilocus heterozygosity (sMLH) measured at 35-43 microsatellites. Within an information theoretic framework, we evaluated support for either fped or sMLH as inbreeding terms and used sequential regression to determine whether the residuals of sMLH on fped explain fitness variation above and beyond fped. We found no evidence of inbreeding depression for survival, either before or after nutritional independence. By contrast, inbreeding was negatively associated with two quality related traits, yearling body mass and annual male reproductive success. Yearling body mass was associated with fped but not sMLH, while male annual reproductive success was best explained by both fped and residual sMLH. Thus, our study not only uncovers variation in the extent to which different traits show inbreeding depression, but also reveals trait-specific differences in the ability of pedigrees and molecular markers to explain fitness variation and suggests that for certain traits genetic markers may capture variation in realised IBDg above and beyond the pedigree expectation. This article is protected by copyright. All rights reserved

Cooperatively breeding banded mongooses do not avoid inbreeding through familiarity-based kin recognition

In species that live in family groups, such as cooperative breeders, inbreeding is usually avoided through the recognition of familiar kin. For example, individuals may avoid mating with conspecifics encountered regularly in infancy, as these likely include parents, siblings, and closely related alloparents. Other mechanisms have also been reported, albeit rarely; for example, individuals may compare their own phenotype to that of others, with close matches representing likely relatives (“phenotype matching”). However, determinants of the primary inbreeding avoidance mechanisms used by a given species remain poorly understood. We use 24 years of life history and genetic data to investigate inbreeding avoidance in wild cooperatively breeding banded mongooses (Mungos mungo). We find that inbreeding avoidance occurs within social groups but is far from maximised (mean pedigree relatedness between 351 breeding pairs = 0.144). Unusually for a group-living vertebrate, we find no evidence that females avoid breeding with males with which they are familiar in early life. This is probably explained by communal breeding; females give birth in tight synchrony and pups are cared for communally, thus reducing the reliability of familiarity-based proxies of relatedness. We also found little evidence that inbreeding is avoided by preferentially breeding with males of specific age classes. Instead, females may exploit as-yet unknown proxies of relatedness, for example, through phenotype matching, or may employ postcopulatory inbreeding avoidance mechanisms. Investigation of species with unusual breeding systems helps to identify constraints against inbreeding avoidance and contributes to our understanding of the distribution of inbreeding across species. Significance statement: Choosing the right mate is never easy, but it may be particularly difficult for banded mongooses. In most social animals, individuals avoid mating with those that were familiar to them as infants, as these are likely to be relatives. However, we show that this rule does not work in banded mongooses. Here, the offspring of several mothers are raised in large communal litters by their social group, and parents seem unable to identify or direct care towards their own pups. This may make it difficult to recognise relatives based on their level of familiarity and is likely to explain why banded mongooses frequently inbreed. Nevertheless, inbreeding is lower than expected if mates are chosen at random, suggesting that alternative pre- or post-copulatory inbreeding avoidance mechanisms are used

Heterogeneous Genomic Divergence Landscape in Two Commercially Important European Scallop Species

Two commercially important scallop species of the genus Pecten are found in Europe: the north Atlantic Pecten maximus and the Mediterranean Pecten jacobaeus whose distributions abut at the Almeria–Orán front. Whilst previous studies have quantified genetic divergence between these species, the pattern of differentiation along the Pecten genome is unknown. Here, we mapped RADseq data from 235 P. maximus and 27 P. jacobaeus to a chromosome-level reference genome, finding a heterogeneous landscape of genomic differentiation. Highly divergent genomic regions were identified across 14 chromosomes, while the remaining five showed little differentiation. Demographic and comparative genomics analyses suggest that this pattern resulted from an initial extended period of isolation, which promoted divergence, followed by differential gene flow across the genome during secondary contact. Single nucleotide polymorphisms present within highly divergent genomic regions were located in areas of low recombination and contrasting patterns of LD decay were found between the two species, hinting at the presence of chromosomal inversions in P. jacobaeus. Functional annotations revealed that highly differentiated regions were enriched for immune-related processes and mRNA modification. While future work is necessary to characterize structural differences, this study provides new insights into the speciation genomics of P. maximus and P. jacobaeus