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

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

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