Genome-wide cross-amplification of domestic sheep microsatellites in bighorn sheep and mountain goats

Contrasting whole-genome and reduced representation sequencing for population demographic and adaptive inference: an alpine mammal case study

Genomes capture the adaptive and demographic history of a species, but the choice of sequencing strategy and sample size can impact such inferences. We compared whole genome and reduced representation sequencing approaches to study the population demographic and adaptive signals of the North American mountain goat (Oreamnos americanus). We applied the restriction site-associated DNA sequencing (RADseq) approach to 254 individuals and whole genome resequencing (WGS) approach to 35 individuals across the species range at mid-level coverage (9X) and to 5 individuals at high coverage (30X). We used ANGSD to estimate the genotype likelihoods and estimated the effective population size (Ne), population structure, and explicitly modelled the demographic history with δaδi and MSMC2. The data sets were overall concordant in supporting a glacial induced vicariance and extremely low Ne in mountain goats. We evaluated a set of climatic variables and geographic location as predictors of genetic diversity using redundancy analysis. A moderate proportion of total variance (36% for WGS and 21% for RADseq data sets) was explained by geography and climate variables; both data sets support a large impact of drift and some degree of local adaptation. The empirical similarities of WGS and RADseq presented herein reassuringly suggest that both approaches will recover large demographic and adaptive signals in a population; however, WGS offers several advantages over RADseq, such as inferring adaptive processes and calculating runs-of-homozygosity estimates. Considering the predicted climate-induced changes in alpine environments and the genetically depauperate mountain goat, the long-term adaptive capabilities of this enigmatic species are questionable.

Genome-wide mining of microsatellites in king cobra (Ophiophagus hannah) and cross-species development of tetranucleotide SSR markers in Chinese cobra (Naja atra)

The complete genome sequence provides the opportunity for genome-wide and coding region analysis of SSRs in the king cobra and for cross-species identification of microsatellite markers in the Chinese cobra. In the Ophiophagus hannah genome, tetranucleotide repeats (38.03%) were the most abundant category, followed by dinucleotides (23.03%), pentanucleotides (13.07%), mononucleotides (11.78%), trinucleotides (11.49%) and hexanucleotides (2.6%). Twenty predominant motifs in the O. hannah genome were (A)n (C)n, (AC)n, (AG)n, (AT)n, (AGG)n, (AAT)n, (AAG)n, (AAC)n, (ATG)n, (ATAG)n, (AAGG)n, (ATCT)n, (CCTT)n, (ATTT)n, (AAAT)n, (AATAG)n, (ATTCT)n, (ATATGT)n, (AGATAT)n. In total, 4344 SSRs were found in coding sequences (CDSs). Tetranucleotides (52.79%) were the most abundant microsatellite type in CDS, followed by trinucleotides (28.50%), dinucleotides (11.02%), pentanucleotides (4.42%), mononucleotides (1.77%), and hexanucleotides (1.50%). A total of 984 CDSs containing microsatellites were assigned 11152 Gene Ontology (GO) functional terms. Gene Ontology (GO) analysis demonstrated that cellular process, cell and binding were the most frequent GO terms in biological process, cellular component and molecular function, respectively. Thirty-two novel highly polymorphic (PIC > 0.5) SSR markers for Naja atra were developed from cross-species amplification based on the tetranucleotide microsatellite sequences in the king cobra genome. The number of alleles (N_A) per locus had between 3 and 11 alleles with an average of 6.5, the polymorphism information content (PIC) value ranged from 0.521 to 0.858 (average = 0.707), the observed heterozygosity (Ho) of 32 microsatellite loci ranged from 0.292 to 0.875 (mean = 0.678), the expected heterozygosity (H_E) ranged from 0.561 to 0.889 (average = 0.761), and 3 microsatellite loci exhibited statistically significant departure from Hardy-Weinberg equilibrium (HWE) after Bonferroni correction (p < 0.003).

Building strong relationships between conservation genetics and primary industry leads to mutually beneficial genomic advances

Several reviews in the past decade have heralded the benefits of embracing high-throughput sequencing technologies to inform conservation policy and the management of threatened species, but few have offered practical advice on how to expedite the transition from conservation genetics to conservation genomics. Here, we argue that an effective and efficient way to navigate this transition is to capitalize on emerging synergies between conservation genetics and primary industry (e.g., agriculture, fisheries, forestry and horticulture). Here, we demonstrate how building strong relationships between conservation geneticists and primary industry scientists is leading to mutually-beneficial outcomes for both disciplines. Based on our collective experience as collaborative New Zealand-based scientists, we also provide insight for forging these cross-sector relationships.

Temporal dynamics of linkage disequilibrium in two populations of bighorn sheep

Linkage disequilibrium (LD) is the nonrandom association of alleles at two markers. Patterns of LD have biological implications as well as practical ones when designing association studies or conservation programs aimed at identifying the genetic basis of fitness differences within and among populations. However, the temporal dynamics of LD in wild populations has received little empirical attention. In this study, we examined the overall extent of LD, the effect of sample size on the accuracy and precision of LD estimates, and the temporal dynamics of LD in two populations of bighorn sheep (Ovis canadensis) with different demographic histories. Using over 200 microsatellite loci, we assessed two metrics of multi-allelic LD, D', and χ ('2). We found that both populations exhibited high levels of LD, although the extent was much shorter in a native population than one that was founded via translocation, experienced a prolonged bottleneck post founding, followed by recent admixture. In addition, we observed significant variation in LD in relation to the sample size used, with small sample sizes leading to depressed estimates of the extent of LD but inflated estimates of background levels of LD. In contrast, there was not much variation in LD among yearly cross-sections within either population once sample size was accounted for. Lack of pronounced interannual variability suggests that researchers may not have to worry about interannual variation when estimating LD in a population and can instead focus on obtaining the largest sample size possible.

Harnessing cross-species alignment to discover SNPs and generate a draft genome sequence of a bighorn sheep (Ovis canadensis)

Background

Whole genome sequences (WGS) have proliferated as sequencing technology continues to improve and costs decline. While many WGS of model or domestic organisms have been produced, a growing number of non-model species are also being sequenced. In the absence of a reference, construction of a genome sequence necessitates de novo assembly which may be beyond the ability of many labs due to the large volumes of raw sequence data and extensive bioinformatics required. In contrast, the presence of a reference WGS allows for alignment which is more tractable than assembly. Recent work has highlighted that the reference need not come from the same species, potentially enabling a wide array of species WGS to be constructed using cross-species alignment. Here we report on the creation a draft WGS from a single bighorn sheep (Ovis canadensis) using alignment to the closely related domestic sheep (Ovis aries).

Results

Two sequencing libraries on SOLiD platforms yielded over 865 million reads, and combined alignment to the domestic sheep reference resulted in a nearly complete sequence (95% coverage of the reference) at an average of 12x read depth (104 SD). From this we discovered over 15 million variants and annotated them relative to the domestic sheep reference. We then conducted an enrichment analysis of those SNPs showing fixed differences between the reference and sequenced individual and found significant differences in a number of gene ontology (GO) terms, including those associated with reproduction, muscle properties, and bone deposition.

Conclusion

Our results demonstrate that cross-species alignment enables the creation of novel WGS for non-model organisms. The bighorn sheep WGS will provide a resource for future resequencing studies or comparative genomics.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1618-x) contains supplementary material, which is available to authorized users.

Beyond mice and men: environmental change, immunity and infections in wild ungulates

In the face of rapid environmental change, anticipating shifts in microparasite and macroparasite dynamics, including emergence events, is an enormous challenge. We argue that immunological studies in natural populations are pivotal to meeting this challenge: many components of environmental change--shifts in biotic assemblages, altered climate patterns and reduced environmental predictability--may affect host immunity. We suggest that wild ungulates can serve as model systems aiding the discovery of immunological mechanisms that link environmental change with parasite transmission dynamics. Our review of eco-immunological studies in wild ungulates reveals progress in understanding how co-infections affect immunity and parasite transmission and how environmental and genetic factors interact to shape immunity. Changes in bioavailability of micronutrients have been linked to immunity and health in wild ungulates. Although physiological stress in response to environmental change has been assessed, downstream effects on immunity have not been studied. Moreover, the taxonomic range of ungulates studied is limited to bovids (bighorn sheep, Soay sheep, chamois, musk oxen, bison, African buffalo) and a few cervids (red deer, black-tailed deer). We discuss areas where future studies in ungulates could lead to significant contributions in understanding the patterns of immunity and infection in natural populations and across species.

Ungulates as model systems for the study of disease processes in natural populations

Parasites and pathogens are a fundamental driving force in the ecology and evolution of mammalian populations, and understanding disease processes in natural populations is an urgent priority in the face of increased rates of infectious disease emergence. In this review, we argue that mammalogists are uniquely placed to contribute to addressing these challenges because in-depth knowledge of mammal species is fundamental to the development of wild model systems that could accelerate discovery in disease ecology. The use of animal models-species for which a broad range of diagnostic, molecular, and genetic tools have been developed-in tightly controlled laboratory environments has been instrumental in driving progress in the biomedical sciences. However, in natural populations, disease processes operate in the context of enormous genetic, phenotypic, and environmental variability. Understanding diseases in animal populations (including humans) thus requires investment in "wild animal models" that explicitly include individual variation and relevant environmental gradients. Wild mammal groups such as primates and rodents have already been identified as potentially useful models of infectious diseases in the wild. Here, we discuss the enormous potential that ungulates hold as candidates for wild model systems. The diversity, broad geographic distribution, and often high abundance of species in this group make them a highly accessible target for disease research. Moreover, a depth of background knowledge, close relationships to domesticated animals, and ongoing management of many wild ungulate species provide context, tools, and opportunity for cutting-edge research at the interface of ecological and biomedical sciences. Studies of wild ungulates are already helping to unravel some key challenges in infectious disease research, including the role of parasites in trophic cascades, the consequences of climate change for disease dynamics, and the systems biology of host-parasite interactions. Other areas where ungulate studies may provide new insight include research on the sources and drivers of emerging infectious diseases.

Genetic structure of Chinese indigenous goats and the special geographical structure in the Southwest China as a geographic barrier driving the fragmentation of a large population

Background

China has numerous native domestic goat breeds, however, extensive studies are focused on the genetic diversity within the fewer breeds and limited regions, the population demograogic history and origin of Chinese goats are still unclear. The roles of geographical structure have not been analyzed in Chinese goat domestic process. In this study, the genetic relationships of Chinese indigenous goat populations were evaluated using 30 microsatellite markers.

Methodology/Principal Findings

Forty Chinese indigenous populations containing 2078 goats were sampled from different geographic regions of China. Moderate genetic diversity at the population level (H_S of 0.644) and high population diversity at the species level (H_T value of 0.737) were estimated. Significant moderate population differentiation was detected (F_ST value of 0.129). Significant excess homozygosity (F_IS of 0.105) and recent population bottlenecks were detected in thirty-six populations. Neighbour-joining tree, principal components analysis and Bayesian clusters all revealed that Chinese goat populations could be subdivided into at least four genetic clusters: Southwest China, South China, Northwest China and East China. It was observed that the genetic diversity of Northern China goats was highest among these clusters. The results here suggested that the goat populations in Southwest China might be the earliest domestic goats in China.

Conclusions/Significance

Our results suggested that the current genetic structure of Chinese goats were resulted from the special geographical structure, especially in the Western China, and the Western goat populations had been separated by the geographic structure (Hengduan Mountains and Qinling Mountains-Huaihe River Line) into two clusters: the Southwest and Northwest. It also indicated that the current genetic structure was caused by the geographical origin mainly, in close accordance with the human’s migration history throughout China. This study provides a fundamental genetic profile for the conservation of these populations and better to understand the domestication process and origin of Chinese goats.

Estimating genome-wide heterozygosity: effects of demographic history and marker type

Heterozygosity-fitness correlations (HFCs) are often used to link individual genetic variation to differences in fitness. However, most studies examining HFCs find weak or no correlations. Here, we derive broad theoretical predictions about how many loci are needed to adequately measure genomic heterozygosity assuming different levels of identity disequilibrium (ID), a proxy for inbreeding. We then evaluate the expected ability to detect HFCs using an empirical data set of 200 microsatellites and 412 single nucleotide polymorphisms (SNPs) genotyped in two populations of bighorn sheep (Ovis canadensis), with different demographic histories. In both populations, heterozygosity was significantly correlated across marker types, although the strength of the correlation was weaker in a native population compared with one founded via translocation and later supplemented with additional individuals. Despite being bi-allelic, SNPs had similar correlations to genome-wide heterozygosity as microsatellites in both populations. For both marker types, this association became stronger and less variable as more markers were considered. Both populations had significant levels of ID; however, estimates were an order of magnitude lower in the native population. As with heterozygosity, SNPs performed similarly to microsatellites, and precision and accuracy of the estimates of ID increased as more loci were considered. Although dependent on the demographic history of the population considered, these results illustrate that genome-wide heterozygosity, and therefore HFCs, are best measured by a large number of markers, a feat now more realistically accomplished with SNPs than microsatellites.

A quantitative trait locus analysis of personality in wild bighorn sheep

Personality, the presence of persistent behav105ioral differences among individuals over time or contexts, potentially has important ecological and evolutionary consequences. However, a lack of knowledge about its genetic architecture limits our ability to understand its origin, evolution, and maintenance. Here, we report on a genome-wide quantitative trait locus (QTL) analysis for two personality traits, docility and boldness, in free-living female bighorn sheep from Ram Mountain, Alberta, Canada. Our variance component linkage analysis based on 238 microsatellite loci genotyped in 310 pedigreed individuals identified suggestive docility and boldness QTL on sheep chromosome 2 and 6, respectively. A lack of QTL overlap indicated that genetic covariance between traits was not modulated by pleiotropic effects at a major locus and may instead result from linkage disequilibrium or pleiotropic effects at QTL of small effects. To our knowledge, this study represents the first attempt to dissect the genetic architecture of personality in a free-living wildlife population, an important step toward understanding the link between molecular genetic variation in personality and fitness and the evolutionary processes maintaining this variation.

Genomic consequences of genetic rescue in an insular population of bighorn sheep (Ovis canadensis)

Genetic rescue is a management intervention whereby a small population is supplemented with individuals from other populations in an attempt to reverse the effects of inbreeding and increased genetic load. One such rescue was recently documented in the population of bighorn sheep (Ovis canadensis) within the National Bison Range wildlife refuge (Montana, USA). Here, we examine the locus-specific effects of rescue in this population using a newly developed genome-wide set of 195 microsatellite loci and first-generation linkage map. We found that the rate of introgression varied among loci and that 111 loci, 57% of those examined, deviated from patterns of neutral inheritance. The most common deviation was an excess of homozygous genotypes relative to neutral expectations, indicative of directional selection. As in previous study of this rescue, individuals with more introduced alleles had higher reproductive success and longevity. In addition, we found 30 loci, distributed throughout the genome, which seem to have individual effects on these life history traits. Although the potential for outbreeding depression is a major concern when translocating individuals between populations, we found no evidence of such effects in this population.

QTL mapping for sexually dimorphic fitness-related traits in wild bighorn sheep

Dissecting the genetic architecture of fitness-related traits in wild populations is key to understanding evolution and the mechanisms maintaining adaptive genetic variation. We took advantage of a recently developed genetic linkage map and phenotypic information from wild pedigreed individuals from Ram Mountain, Alberta, Canada, to study the genetic architecture of ecologically important traits (horn volume, length, base circumference and body mass) in bighorn sheep. In addition to estimating sex-specific and cross-sex quantitative genetic parameters, we tested for the presence of quantitative trait loci (QTLs), colocalization of QTLs between bighorn sheep and domestic sheep, and sex × QTL interactions. All traits showed significant additive genetic variance and genetic correlations tended to be positive. Linkage analysis based on 241 microsatellite loci typed in 310 pedigreed animals resulted in no significant and five suggestive QTLs (four for horn dimension on chromosomes 1, 18 and 23, and one for body mass on chromosome 26) using genome-wide significance thresholds (Logarithm of odds (LOD) >3.31 and >1.88, respectively). We also confirmed the presence of a horn dimension QTL in bighorn sheep at the only position known to contain a similar QTL in domestic sheep (on chromosome 10 near the horns locus; nominal P<0.01) and highlighted a number of regions potentially containing weight-related QTLs in both species. As expected for sexually dimorphic traits involved in male-male combat, loci with sex-specific effects were detected. This study lays the foundation for future work on adaptive genetic variation and the evolutionary dynamics of sexually dimorphic traits in bighorn sheep.

Temporal dynamics of genetic variability in a mountain goat (Oreamnos americanus) population

The association between population dynamics and genetic variability is of fundamental importance for both evolutionary and conservation biology. We combined long-term population monitoring and molecular genetic data from 123 offspring and their parents at 28 microsatellite loci to investigate changes in genetic diversity over 14 cohorts in a small and relatively isolated population of mountain goats (Oreamnos americanus) during a period of demographic increase. Offspring heterozygosity decreased while parental genetic similarity and inbreeding coefficients (F(IS) ) increased over the study period (1995-2008). Immigrants introduced three novel alleles into the population and matings between residents and immigrants produced more heterozygous offspring than local crosses, suggesting that immigration can increase population genetic variability. The population experienced genetic drift over the study period, reflected by a reduced allelic richness over time and an 'isolation-by-time' pattern of genetic structure. The temporal decline of individual genetic diversity despite increasing population size probably resulted from a combination of genetic drift due to small effective population size, inbreeding and insufficient counterbalancing by immigration. This study highlights the importance of long-term genetic monitoring to understand how demographic processes influence temporal changes of genetic diversity in long-lived organisms.

Does reduced heterozygosity influence dispersal? A test using spatially structured populations in an alpine ungulate

Despite having a profound effect on population dynamics, the reasons that animals disperse are poorly understood. Evolutionary explanations have focused on inbreeding and competition, where the potential cost of philopatry is negated through dispersal. Such scenarios lead to the prediction that less successful individuals preferentially disperse, termed 'fitness-associated dispersal'. Since heterozygosity is associated with fitness, we assessed whether dispersed animals had less observed heterozygosity (H(O)) than residents. We tested this prediction using both genetic and population-monitoring data of mountain goats (Oreamnos americanus). Individuals classified as dispersers through cross-assignment had the lowest mean H(O), followed by residents, and then admixed individuals. Dispersed individuals had 6.3 per cent less H(O) than their subpopulation of origin. In the long-term study of the mountain goat herd at Caw Ridge, Alberta, immigrants had the lowest H(O); however, the opposite pattern was seen in emigrants, which may be related to density dependence. This study is the first to provide empirical evidence that heterozygosity is associated with dispersal.

Genetic linkage map of a wild genome: genomic structure, recombination and sexual dimorphism in bighorn sheep

Background

The construction of genetic linkage maps in free-living populations is a promising tool for the study of evolution. However, such maps are rare because it is difficult to develop both wild pedigrees and corresponding sets of molecular markers that are sufficiently large. We took advantage of two long-term field studies of pedigreed individuals and genomic resources originally developed for domestic sheep (Ovis aries) to construct a linkage map for bighorn sheep, Ovis canadensis. We then assessed variability in genomic structure and recombination rates between bighorn sheep populations and sheep species.

Results

Bighorn sheep population-specific maps differed slightly in contiguity but were otherwise very similar in terms of genomic structure and recombination rates. The joint analysis of the two pedigrees resulted in a highly contiguous map composed of 247 microsatellite markers distributed along all 26 autosomes and the X chromosome. The map is estimated to cover about 84% of the bighorn sheep genome and contains 240 unique positions spanning a sex-averaged distance of 3051 cM with an average inter-marker distance of 14.3 cM. Marker synteny, order, sex-averaged interval lengths and sex-averaged total map lengths were all very similar between sheep species. However, in contrast to domestic sheep, but consistent with the usual pattern for a placental mammal, recombination rates in bighorn sheep were significantly greater in females than in males (~12% difference), resulting in an autosomal female map of 3166 cM and an autosomal male map of 2831 cM. Despite differing genome-wide patterns of heterochiasmy between the sheep species, sexual dimorphism in recombination rates was correlated between orthologous intervals.

Conclusions

We have developed a first-generation bighorn sheep linkage map that will facilitate future studies of the genetic architecture of trait variation in this species. While domestication has been hypothesized to be responsible for the elevated mean recombination rate observed in domestic sheep, our results suggest that it is a characteristic of Ovis species. However, domestication may have played a role in altering patterns of heterochiasmy. Finally, we found that interval-specific patterns of sexual dimorphism were preserved among closely related Ovis species, possibly due to the conserved position of these intervals relative to the centromeres and telomeres. This study exemplifies how transferring genomic resources from domesticated species to close wild relative can benefit evolutionary ecologists while providing insights into the evolution of genomic structure and recombination rates of domesticated species.

A genome-wide set of SNPs detects population substructure and long range linkage disequilibrium in wild sheep

The development of genomic resources for wild species is still in its infancy. However, cross-species utilization of technologies developed for their domestic counterparts has the potential to unlock the genomes of organisms that currently lack genomic resources. Here, we apply the OvineSNP50 BeadChip, developed for domestic sheep, to two related wild ungulate species: the bighorn sheep (Ovis canadensis) and the thinhorn sheep (Ovis dalli). Over 95% of the domestic sheep markers were successfully genotyped in a sample of fifty-two bighorn sheep while over 90% were genotyped in two thinhorn sheep. Pooling the results from both species identified 868 single-nucleotide polymorphisms (SNPs), 570 were detected in bighorn sheep, while 330 SNPs were identified in thinhorn sheep. The total panel of SNPs was able to discriminate between the two species, assign population of origin for bighorn sheep and detect known relationship classes within one population of bighorn sheep. Using an informative subset of these SNPs (n=308), we examined the extent of genome-wide linkage disequilibrium (LD) within one population of bighorn sheep and found that high levels of LD persist over 4 Mb.