An 85K SNP Array Uncovers Inbreeding and Cryptic Relatedness in an Antarctic Fur Seal Breeding Colony

Refinement of the Antarctic fur seal (Arctocephalus gazella) reference genome increases continuity and completeness

The Antarctic fur seal (Arctocephalus gazella) is an important top predator and indicator of the health of the Southern Ocean ecosystem. Although abundant, this species narrowly escaped extinction due to historical sealing and is currently declining as a consequence of climate change. Genomic tools are essential for understanding these anthropogenic impacts and for predicting long-term viability. However, the current reference genome ("arcGaz3") shows considerable room for improvement in terms of both completeness and contiguity. We therefore combined PacBio sequencing, haplotype-aware HiRise assembly, and scaffolding based on Hi-C information to generate a refined assembly of the Antarctic fur seal reference genome ("arcGaz4_h1"). The new assembly is 2.53 Gb long, has a scaffold N50 of 55.6 Mb and includes 18 chromosome-sized scaffolds, which correspond to the 18 chromosomes expected in otariids. Genome completeness is greatly improved, with 23,408 annotated genes and a Benchmarking Universal Single-Copy Orthologs score raised from 84.7% to 95.2%. We furthermore included the new genome in a reference-free alignment of the genomes of 11 pinniped species to characterize evolutionary conservation across the Pinnipedia using genome-wide Genomic Evolutionary Rate Profiling. We then implemented Gene Ontology enrichment analyses to identify biological processes associated with those genes showing the highest levels of either conservation or differentiation between the 2 major pinniped families, the Otariidae and Phocidae. We show that processes linked to neuronal development, the circulatory system, and osmoregulation are overrepresented both in conserved as well as in differentiated regions of the genome.

Genetic diversity analysis of Inner Mongolia cashmere goats (Erlangshan subtype) based on whole genome re-sequencing

BACKGROUND

Inner Mongolia cashmere goat (IMCG), renowned for its superior cashmere quality, is a Chinese indigenous goat breed that has been developed through natural and artificial selection over a long period. However, recently, the genetic resources of IMCGs have been significantly threatened by the introduction of cosmopolitan goat breeds and the absence of adequate breed protection systems.

RESULTS

In order to assess the conservation effectiveness of IMCGs and efficiently preserve and utilize the purebred germplasm resources, this study analyzed the genetic diversity, kinship, family structure, and inbreeding of IMCGs utilizing resequencing data from 225 randomly selected individuals analyzed using the Plink (v.1.90), GCTA (v.1.94.1), and R (v.4.2.1) software. A total of 12,700,178 high-quality SNPs were selected through quality control from 34,248,064 SNP sites obtained from 225 individuals. The average minor allele frequency (MAF), polymorphic information content (PIC), and Shannon information index (SHI) were 0.253, 0.284, and 0.530, respectively. The average observed heterozygosity (Ho) and the average expected heterozygosity (He) were 0.355 and 0.351, respectively. The analysis of the identity by state distance matrix and genomic relationship matrix has shown that most individuals' genetic distance and genetic relationship are far away, and the inbreeding coefficient is low. The family structure analysis identified 10 families among the 23 rams. A total of 14,109 runs of homozygosity (ROH) were identified in the 225 individuals, with an average ROH length of 1014.547 kb. The average inbreeding coefficient, calculated from ROH, was 0.026 for the overall population and 0.027 specifically among the 23 rams, indicating a low level of inbreeding within the conserved population.

CONCLUSIONS

The IMCGs exhibited moderate polymorphism and a low level of kinship with inbreeding occurring among a limited number of individuals. Simultaneously, it is necessary to prevent the loss of bloodline to guarantee the perpetuation of the IMCGs' germplasm resources.

Extant and extinct bilby genomes combined with Indigenous knowledge improve conservation of a unique Australian marsupial

Ninu (greater bilby, Macrotis lagotis) are desert-dwelling, culturally and ecologically important marsupials. In collaboration with Indigenous rangers and conservation managers, we generated the Ninu chromosome-level genome assembly (3.66 Gbp) and genome sequences for the extinct Yallara (lesser bilby, Macrotis leucura). We developed and tested a scat single-nucleotide polymorphism panel to inform current and future conservation actions, undertake ecological assessments and improve our understanding of Ninu genetic diversity in managed and wild populations. We also assessed the beneficial impact of translocations in the metapopulation (N = 363 Ninu). Resequenced genomes (temperate Ninu, 6; semi-arid Ninu, 6; and Yallara, 4) revealed two major population crashes during global cooling events for both species and differences in Ninu genes involved in anatomical and metabolic pathways. Despite their 45-year captive history, Ninu have fewer long runs of homozygosity than other larger mammals, which may be attributable to their boom-bust life history. Here we investigated the unique Ninu biology using 12 tissue transcriptomes revealing expression of all 115 conserved eutherian chorioallantoic placentation genes in the uterus, an XY1Y2 sex chromosome system and olfactory receptor gene expansions. Together, we demonstrate the holistic value of genomics in improving key conservation actions, understanding unique biological traits and developing tools for Indigenous rangers to monitor remote wild populations.

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.

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.

From STRs to SNPs via ddRAD-seq: Geographic assignment of confiscated tortoises at reduced costs

Assigning individuals to their source populations is crucial for conservation research, especially for endangered species threatened by illegal trade and translocations. Genetic assignment can be achieved with different types of molecular markers, but technical advantages and cost saving are recently promoting the shift from short tandem repeats (STRs) to single nucleotide polymorphisms (SNPs). Here, we designed, developed, and tested a small panel of SNPs for cost-effective geographic assignment of individuals with unknown origin of the endangered Mediterranean tortoise Testudo hermanni. We started by performing a ddRAD-seq experiment on 70 wild individuals of T. hermanni from 38 locations. Results obtained using 3182 SNPs are comparable to those previously obtained using STR markers in terms of genetic structure and power to identify the macro-area of origin. However, our SNPs revealed further insights into the substructure in Western populations, especially in Southern Italy. A small panel of highly informative SNPs was then selected and tested by genotyping 190 individuals using the KASP genotyping chemistry. All the samples from wild populations of known geographic origin were genetically re-assigned with high accuracy to the original population. This reduced SNPs panel represents an efficient molecular tool that enables individuals to be genotyped at low cost (less than €15 per sample) for geographical assignment and identification of hybrids. This information is crucial for the management in-situ of confiscated animals and their possible re-allocation in the wild. Our methodological pipeline can easily be extended to other species.

Dynamics of reduced genetic diversity in increasingly fragmented populations of Florida scrub jays,Aphelocoma coerulescens

Understanding the genomic consequences of population decline is important for predicting species' vulnerability to intensifying global change. Empirical information about genomic changes in populations in the early stages of decline, especially for those still experiencing immigration, remains scarce. We used 7834 autosomal SNPs and demographic data for 288 Florida scrub jays (Aphelocoma coerulescens; FSJ) sampled in 2000 and 2008 to compare levels of genetic diversity, inbreeding, relatedness, and lengths of runs of homozygosity (ROH) between two subpopulations within dispersal distance of one another but have experienced contrasting demographic trajectories. At Archbold Biological Station (ABS), the FSJ population has been stable because of consistent habitat protection and management, while at nearby Placid Lakes Estates (PLE), the population declined precipitously due to suburban development. By the onset of our sampling in 2000, birds in PLE were already less heterozygous, more inbred, and on average more related than birds in ABS. No significant changes occurred in heterozygosity or inbreeding across the 8-year sampling interval, but average relatedness among individuals decreased in PLE, thus by 2008 average relatedness did not differ between sites. PLE harbored a similar proportion of short ROH but a greater proportion of long ROH than ABS, suggesting one continuous population of shared demographic history in the past, which is now experiencing more recent inbreeding. These results broadly uphold the predictions of simple population genetic models based on inferred effective population sizes and rates of immigration. Our study highlights how, in just a few generations, formerly continuous populations can diverge in heterozygosity and levels of inbreeding with severe local population decline despite ongoing gene flow.

Low heritability and high phenotypic plasticity of salivary cortisol in response to environmental heterogeneity in a wild pinniped

Individuals are unique in how they interact with and respond to their environment. Correspondingly, unpredictable challenges or environmental stressors often produce an individualized response of the hypothalamic-pituitary-adrenal (HPA) axis and its downstream effector cortisol. We used a fully crossed, repeated measures design to investigate the factors shaping individual variation in baseline cortisol in Antarctic fur seal pups and their mothers. Saliva samples were collected from focal individuals at two breeding colonies, one with low and the other with high density, during two consecutive years of contrasting food availability. Mothers and pups were sampled concurrently at birth and shortly before weaning, while pups were additionally sampled every 20 days. We found that heritability was low for baseline cortisol, while within-individual repeatability and among-individual variability were high. A substantial proportion of the variation in baseline cortisol could be explained in pups and mothers by a combination of intrinsic and extrinsic factors including sex, weight, day, season, and colony of birth. Our findings provide detailed insights into the individualization of endocrine phenotypes and their genetic and environmental drivers in a wild pinniped. Furthermore, the strong associations between cortisol and life history traits that we report in fur seals could have important implications for understanding the population dynamics of species impacted by environmental change.

The design and application of a 50 K SNP chip for a threatened Aotearoa New Zealand passerine, the hihi

Next-generation sequencing has transformed the fields of ecological and evolutionary genetics by allowing for cost-effective identification of genome-wide variation. Single nucleotide polymorphism (SNP) arrays, or "SNP chips", enable very large numbers of individuals to be consistently genotyped at a selected set of these identified markers, and also offer the advantage of being able to analyse samples of variable DNA quality. We used reduced representation restriction-aided digest sequencing (RAD-seq) of 31 birds of the threatened hihi (Notiomystis cincta; stitchbird) and low-coverage whole genome sequencing (WGS) of 10 of these birds to develop an Affymetrix 50 K SNP chip. We overcame the limitations of having no hihi reference genome and a low quantity of sequence data by separate and pooled de novo assembly of each of the 10 WGS birds. Reads from all individuals were mapped back to these de novo assemblies to identify SNPs. A subset of RAD-seq and WGS SNPs were selected for inclusion on the chip, prioritising SNPs with the highest quality scores whose flanking sequence uniquely aligned to the zebra finch (Taeniopygia guttata) genome. Of the 58,466 SNPs manufactured on the chip, 72% passed filtering metrics and were polymorphic. By genotyping 1,536 hihi on the array, we found that SNPs detected in multiple assemblies were more likely to successfully genotype, representing a cost-effective approach to identify SNPs for genotyping. Here, we demonstrate the utility of the SNP chip by describing the high rates of linkage disequilibrium in the hihi genome, reflecting the history of population bottlenecks in the species.

Genomic data of different resolutions reveal consistent inbreeding estimates but contrasting homozygosity landscapes for the threatened Aotearoa New Zealand hihi

Inbreeding can lead to a loss of heterozygosity in a population and when combined with genetic drift may reduce the adaptive potential of a species. However, there is uncertainty about whether resequencing data can provide accurate and consistent inbreeding estimates. Here, we performed an in-depth inbreeding analysis for hihi (Notiomystis cincta), an endemic and nationally vulnerable passerine bird of Aotearoa New Zealand. We first focused on subsampling variants from a reference genome male, and found that low-density data sets tend to miss runs of homozygosity (ROH) in some places and overestimate ROH length in others, resulting in contrasting homozygosity landscapes. Low-coverage resequencing and 50 K SNP array densities can yield comparable inbreeding results to high-coverage resequencing approaches, but the results for all data sets are highly dependent on the software settings employed. Second, we extended our analysis to 10 hihi where low-coverage whole genome resequencing, RAD-seq and SNP array genotypes are available. We inferred ROH and individual inbreeding to evaluate the relative effects of sequencing depth versus SNP density on estimating inbreeding coefficients and found that high rates of missingness downwardly bias both the number and length of ROH. In summary, when using genomic data to evaluate inbreeding, studies must consider that ROH estimates are heavily dependent on analysis parameters, data set density and individual sequencing depth.

Chemical patterns of colony membership and mother-offspring similarity in Antarctic fur seals are reproducible

Replication studies are essential for evaluating the validity of previous research findings. However, it has proven challenging to reproduce the results of ecological and evolutionary studies, partly because of the complexity and lability of many of the phenomena being investigated, but also due to small sample sizes, low statistical power and publication bias. Additionally, replication is often considered too difficult in field settings where many factors are beyond the investigator’s control and where spatial and temporal dependencies may be strong. We investigated the feasibility of reproducing original research findings in the field of chemical ecology by performing an exact replication of a previous study of Antarctic fur seals (Arctocephalus gazella). In the original study, skin swabs from 41 mother-offspring pairs from two adjacent breeding colonies on Bird Island, South Georgia, were analyzed using gas chromatography-mass spectrometry. Seals from the two colonies differed significantly in their chemical fingerprints, suggesting that colony membership may be chemically encoded, and mothers were also chemically similar to their pups, hinting at the possible involvement of phenotype matching in mother-offspring recognition. In the current study, we generated and analyzed chemical data from a non-overlapping sample of 50 mother-offspring pairs from the same two colonies 5 years later. The original results were corroborated in both hypothesis testing and estimation contexts, with p-values remaining highly significant and effect sizes, standardized between studies by bootstrapping the chemical data over individuals, being of comparable magnitude. However, exact replication studies are only capable of showing whether a given effect can be replicated in a specific setting. We therefore investigated whether chemical signatures are colony-specific in general by expanding the geographic coverage of our study to include pups from a total of six colonies around Bird Island. We detected significant chemical differences in all but a handful of pairwise comparisons between colonies. This finding adds weight to our original conclusion that colony membership is chemically encoded, and suggests that chemical patterns of colony membership not only persist over time but can also be generalized over space. Our study systematically confirms and extends our previous findings, while also implying more broadly that spatial and temporal heterogeneity need not necessarily negate the reproduction and generalization of ecological research findings.