Opportunities and challenges of macrogenetic studies.

The rapidly emerging field of macrogenetics focuses on analysing publicly accessible genetic datasets from thousands of species to explore large-scale patterns and predictors of intraspecific genetic variation. Facilitated by advances in evolutionary biology, technology, data infrastructure, statistics and open science, macrogenetics addresses core evolutionary hypotheses (such as disentangling environmental and life-history effects on genetic variation) with a global focus. Yet, there are important, often overlooked, limitations to this approach and best practices need to be considered and adopted if macrogenetics is to continue its exciting trajectory and reach its full potential in fields such as biodiversity monitoring and conservation. Here, we review the history of this rapidly growing field, highlight knowledge gaps and future directions, and provide guidelines for further research.

The crucial role of genome-wide genetic variation in conservation.

The unprecedented rate of extinction calls for efficient use of genetics to help conserve biodiversity. Several recent genomic and simulation-based studies have argued that the field of conservation biology has placed too much focus on conserving genome-wide genetic variation, and that the field should instead focus on managing the subset of functional genetic variation that is thought to affect fitness. Here, we critically evaluate the feasibility and likely benefits of this approach in conservation. We find that population genetics theory and empirical results show that conserving genome-wide genetic variation is generally the best approach to prevent inbreeding depression and loss of adaptive potential from driving populations toward extinction. Focusing conservation efforts on presumably functional genetic variation will only be feasible occasionally, often misleading, and counterproductive when prioritized over genome-wide genetic variation. Given the increasing rate of habitat loss and other environmental changes, failure to recognize the detrimental effects of lost genome-wide genetic variation on long-term population viability will only worsen the biodiversity crisis.

Historical Introgression from Wild Relatives Enhanced Climatic Adaptation and Resistance to Pneumonia in Sheep.

How animals, particularly livestock, adapt to various climates and environments over short evolutionary time is of fundamental biological interest. Further, understanding the genetic mechanisms of adaptation in indigenous livestock populations is important for designing appropriate breeding programs to cope with the impacts of changing climate. Here, we conducted a comprehensive genomic analysis of diversity, interspecies introgression, and climate-mediated selective signatures in a global sample of sheep and their wild relatives. By examining 600K and 50K genome-wide single nucleotide polymorphism data from 3,447 samples representing 111 domestic sheep populations and 403 samples from all their seven wild relatives (argali, Asiatic mouflon, European mouflon, urial, snow sheep, bighorn, and thinhorn sheep), coupled with 88 whole-genome sequences, we detected clear signals of common introgression from wild relatives into sympatric domestic populations, thereby increasing their genomic diversities. The introgressions provided beneficial genetic variants in native populations, which were significantly associated with local climatic adaptation. We observed common introgression signals of alleles in olfactory-related genes (e.g., ADCY3 and TRPV1) and the PADI gene family including in particular PADI2, which is associated with antibacterial innate immunity. Further analyses of whole-genome sequences showed that the introgressed alleles in a specific region of PADI2 (chr2: 248,302,667-248,306,614) correlate with resistance to pneumonia. We conclude that wild introgression enhanced climatic adaptation and resistance to pneumonia in sheep. This has enabled them to adapt to varying climatic and environmental conditions after domestication.

Ewe are what ewe wear: bigger horns, better ewes and the potential consequence of trophy hunting on female fitness in bighorn sheep.

In polygynous species, secondary sexual traits such as weapons or elaborate ornaments have evolved through intrasexual competition for mates. In some species, these traits are present in both sexes but are underdeveloped in the sex facing lower intrasexual competition for mates. It is often assumed that these underdeveloped sexually selected traits are a vestige of strong sexual selection on the other sex. Here, we challenge this assumption and investigate whether the expression of secondary sexual traits is associated with fitness in female bighorn sheep. Analyses of 45 years of data revealed that female horn length at 2 years, while accounting for mass and environmental variables, is associated with younger age at primiparity, younger age of first offspring weaned, greater reproductive lifespan and higher lifetime reproductive success. There was no association between horn length and fecundity. These findings highlight a potential conservation issue. In this population, trophy hunting selects against males with fast-growing horns. Intersexual genetic correlations imply that intense selective hunting of large-horned males before they can reproduce can decrease female horn size. Therefore, intense trophy hunting of males based on horn size could reduce female reproductive performance through the associations identified here, and ultimately reduce population growth and viability.

Linking genetic, morphological, and behavioural divergence between inland island and mainland deer mice.

The island syndrome hypothesis (ISH) stipulates that, as a result of local selection pressures and restricted gene flow, individuals from island populations should differ from individuals within mainland populations. Specifically, island populations are predicted to contain individuals that are larger, less aggressive, more sociable, and that invest more in their offspring. To date, tests of the ISH have mainly compared oceanic islands to continental sites, and rarely smaller spatial scales such as inland watersheds. Here, using a novel set of genome-wide SNP markers in wild deer mice (Peromyscus maniculatus) we conducted a genomic assessment of predictions underlying the ISH in an inland riverine island system: analysing island-mainland population structure, and quantifying heritability of phenotypes thought to underlie the ISH. We found clear genomic differentiation between the island and mainland populations and moderate to high marker-based heritability estimates for overall variation in traits previously found to differ in line with the ISH between mainland and island locations. FST outlier analyses highlighted 12 loci associated with differentiation between mainland and island populations. Together these results suggest that the island populations examined are on independent evolutionary trajectories, the traits considered have a genetic basis (rather than phenotypic variation being solely due to phenotypic plasticity). Coupled with the previous results showing significant phenotypic differentiation between the island and mainland groups in this system, this study suggests that the ISH can hold even on a small spatial scale.

A new lineage of Galapagos giant tortoises identified from museum samples.

The Galapagos Archipelago is recognized as a natural laboratory for studying evolutionary processes. San Cristóbal was one of the first islands colonized by tortoises, which radiated from there across the archipelago to inhabit 10 islands. Here, we sequenced the mitochondrial control region from six historical giant tortoises from San Cristóbal (five long deceased individuals found in a cave and one found alive during an expedition in 1906) and discovered that the five from the cave are from a clade that is distinct among known Galapagos giant tortoises but closely related to the species from Española and Pinta Islands. The haplotype of the individual collected alive in 1906 is in the same clade as the haplotype in the contemporary population. To search for traces of a second lineage in the contemporary population on San Cristóbal, we closely examined the population by sequencing the mitochondrial control region for 129 individuals and genotyping 70 of these for both 21 microsatellite loci and >12,000 genome-wide single nucleotide polymorphisms [SNPs]. Only a single mitochondrial haplotype was found, with no evidence to suggest substructure based on the nuclear markers. Given the geographic and temporal proximity of the two deeply divergent mitochondrial lineages in the historical samples, they were likely sympatric, raising the possibility that the lineages coexisted. Without the museum samples, this important discovery of an additional lineage of Galapagos giant tortoise would not have been possible, underscoring the value of such collections and providing insights into the early evolution of this iconic radiation.

Macrogenetic studies must not ignore limitations of genetic markers and scale.

Millette et al. (Ecology Letters, 2020, 23:55-67) reported no consistent worldwide anthropogenic effects on animal genetic diversity using repurposed mitochondrial DNA sequences. We reexamine data from this study, describe genetic marker and scale limitations which might lead to misinterpretations with conservation implications, and provide advice to improve future macrogenetic studies.

Ancient hybridization patterns between bighorn and thinhorn sheep.

Whole-genome sequencing has advanced the study of species evolution, including the detection of genealogical discordant events such as ancient hybridization and incomplete lineage sorting (ILS). The evolutionary history of bighorn (Ovis canadensis) and thinhorn (Ovis dalli) sheep present an ideal system to investigate evolutionary discordance due to their recent and rapid radiation and putative secondary contact between bighorn and thinhorn sheep subspecies, specifically the dark pelage Stone sheep (O. dalli stonei) and predominately white Dall sheep (O. dalli dalli), during the last ice age. Here, we used multiple genomes of bighorn and thinhorn sheep, together with snow (O. nivicola) and the domestic sheep (O. aries) as outgroups, to assess their phylogenomic history, potential introgression patterns and their adaptive consequences. Among the Pachyceriforms (snow, bighorn and thinhorn sheep) a consistent monophyletic species tree was retrieved; however, many genealogical discordance patterns were observed. Alternative phylogenies frequently placed Stone and bighorn as sister clades. This relationship occurred more often and was less divergent than that between Dall and bighorn. We also observed many blocks containing introgression signal between Stone and bighorn genomes in which coat colour genes were present. Introgression signals observed between Dall and bighorn were more random and less frequent, and therefore probably due to ILS or intermediary secondary contact. These results strongly suggest that Stone sheep originated from a complex series of events, characterized by multiple, ancient periods of secondary contact with bighorn sheep.

The genetic basis of animal behavioural diversity in natural populations.

Individual differences in animal behaviour influence ecological and evolutionary processes. Much behavioural variation has a heritable component, suggesting that genetics may play a role in its development. Yet, the study of the mechanistic description linking genes to behaviour in nature remains in its infancy, and such research is considered a challenge in contemporary biology. Here, we performed a literature review and meta-analysis to assess trends in analytical approaches used to investigate the relationship between genes and behaviour in natural systems, specifically candidate gene approaches, quantitative trait locus (QTL) mapping, and genome-wide association studies (GWAS). We aimed to determine the efficacy and success of each approach, while also describing which behaviours and species were examined by researchers most often. We found that the majority of QTL mapping and GWAS results revealed a significant or suggestive effect (Zr = 0.3 [95% CI: 0.25:0.35] and Zr = 0.39 [0.33:0.46], respectively) between the trait of interest and genetic marker(s) tested, while over half of candidate gene accounts (Zr = 0.16 [0.11:0.21]) did not find a significant association. Approximately a third of all study estimates investigated animal personality traits; though, reproductive and migratory behaviours were also well-represented. Our findings show that despite widespread accessibility of molecular approaches given current sequencing technologies, efforts to elucidate the genetic basis of behaviour in free-ranging systems has been limited to relatively few species. We discuss challenges encountered by researchers, and recommend integration of novel genomic methods with longitudinal studies to usher in the next wave of behavioural genomic research.

Confidently identifying the correct K value using the ΔK method: When does K = 2?

Populations delineated based on genetic data are commonly used for wildlife conservation and management. Many studies use the program structure combined with the ΔK method to identify the most probable number of populations (K). We recently found K = 2 was identified more often when studies used ΔK compared to studies that did not. We suggested two reasons for this: hierarchical population structure leads to underestimation, or the ΔK method does not evaluate K = 1 causing an overestimation. The present contribution aims to develop a better understanding of the limits of the method using one, two and three population simulations across migration scenarios. From these simulations we identified the "best K" using model likelihood and ΔK. Our findings show that mean probability plots and ΔK are unable to resolve the correct number of populations once migration rate exceeds 0.005. We also found a strong bias towards selecting K = 2 using the ΔK method. We used these data to identify the range of values where the ΔK statistic identifies a value of K that is not well supported. Finally, using the simulations and a review of empirical data, we found that the magnitude of ΔK corresponds to the level of divergence between populations. Based on our findings, we suggest researchers should use the ΔK method cautiously; they need to report all relevant data, including the magnitude of ΔK, and an estimate of connectivity for the research community to assess whether meaningful genetic structure exists within the context of management and conservation.

The influence of a priori grouping on inference of genetic clusters: simulation study and literature review of the DAPC method.

Inference of genetic clusters is a key aim of population genetics, sparking development of numerous analytical methods. Within these, there is a conceptual divide between finding de novo structure versus assessment of a priori groups. Recently developed, Discriminant Analysis of Principal Components (DAPC), combines discriminant analysis (DA) with principal component (PC) analysis. When applying DAPC, the groups used in the DA (specified a priori or described de novo) need to be carefully assessed. While DAPC has rapidly become a core technique, the sensitivity of the method to misspecification of groups and how it is being empirically applied, are unknown. To address this, we conducted a simulation study examining the influence of a priori versus de novo group designations, and a literature review of how DAPC is being applied. We found that with a priori groupings, distance between genetic clusters reflected underlying FST. However, when migration rates were high and groups were described de novo there was considerable inaccuracy, both in terms of the number of genetic clusters suggested and placement of individuals into those clusters. Nearly all (90.1%) of 224 studies surveyed used DAPC to find de novo clusters, and for the majority (62.5%) the stated goal matched the results. However, most studies (52.3%) omit key run parameters, preventing repeatability and transparency. Therefore, we present recommendations for standard reporting of parameters used in DAPC analyses. The influence of groupings in genetic clustering is not unique to DAPC, and researchers need to consider their goal and which methods will be most appropriate.

Genetically informed captive breeding of hybrids of an extinct species of Galapagos tortoise.

Hybridization poses a major challenge for species conservation because it threatens both genetic integrity and adaptive potential. Yet, hybridization can occasionally offer unprecedented opportunity for species recovery if the genome of an extinct taxon is present among living hybrids such that selective breeding could recapture it. We explored the design elements for establishing a captive-breeding program for Galapagos tortoises (Chelonoidis spp.) built around individuals with admixed ancestry involving an extinct species. The target individuals were hybrids between the extinct species from Floreana Island, C. niger, and an extant species, C. becki, which were recently found in the endemic range of C. becki, from Wolf Volcano on Isabela Island. We combined genotypic data from 35 tortoises with high ancestry from C. niger with forward-in-time simulations to explore captive breeding strategies that maximized overall genetic diversity and ancestry from C. niger while accommodating resource constraints, species biology, and the urgency to return tortoises to Floreana Island for facilitating ecosystem restoration. Overall genetic diversity was maximized when in the simulation tortoises were organized in relatively small breeding groups. Substantial amounts of the C. niger genome were captured despite limited resources available for selectively breeding tortoises in captivity. Genetic diversity was maximized when captive-bred offspring were released to the wild rather than being used as additional breeders. Our results provide genetic-based and practical guidance on the inclusion of hybrids with genomic representation from extinct taxa into species restoration programs and informs the ongoing debate on the value of hybrids in biodiversity conservation.

Reproducción en Cautiverio Informada Genéticamente de Híbridos de una Especie Extinta de Tortuga de las Galápagos Resumen La hibridación representa un obstáculo importante para la conservación de especies ya que amenaza tanto a la integridad genética como al potencial adaptativo. Aun así, la hibridación ocasionalmente puede ofrecer una oportunidad sin precedentes para la recuperación de una especie si el genoma de un taxón extinto está presente entre los híbridos vivientes de tal manera que la reproducción selectiva pudiera recuperarlo. Exploramos los elementos de diseño para el establecimiento de un programa de reproducción en cautiverio de la tortuga de las Galápagos (Chelonoidis spp.) construido en torno a los individuos con linajes mixtos que incluyeran una especie extinta. Los individuos fueron los híbridos de la especie extinta en la Isla Floreana, C. niger, y la especie viviente C. becki, encontrados recientemente en la distribución geográfica endémica de la segunda especie en el Volcán Wolf (Isla Isabela). Combinamos los datos genotípicos de 35 tortugas con un linaje cargado de C. niger usando simulaciones futuras de la descendencia generada por el programa para explorar las estrategias de reproducción en cautiverio que maximizaran en general la diversidad genética y el linaje de C. niger a la vez que se ajustaba a las restricciones de recursos, la biología de la especie y la urgencia por regresar las tortugas a la Isla Floreana para facilitar la restauración del ecosistema. En general, la diversidad genética se maximizó cuando en la simulación las tortugas estuvieron organizadas en grupos de reproducción relativamente pequeños y cuando cantidades sustanciales del genoma de C. niger fueron capturados con base en los recursos disponibles para reproducir selectivamente a las tortugas en cautiverio. La diversidad genética se vio especialmente maximizada cuando las crías reproducidas en cautiverio fueron liberadas en lugar de ser utilizadas como reproductoras adicionales. Nuestros resultados proporcionan una guía práctica y basada en la genética para la inclusión de híbridos con representación genómica de un taxón extinto en los programas de restauración de especies. Cuando incorporamos a los híbridos con diversidad genética que previamente se creía perdida en los programas con el propósito de la reintroducción de especies, nuestro estudio informa al debate continuo sobre el valor de los híbridos para la conservación de la biodiversidad.

From transects to transcripts: Teasing apart the architecture of reproductive isolation.

Understanding the processes underlying speciation has long been a challenge to evolutionary biologists. This spurs from difficulties teasing apart the various mechanisms that contribute to the evolution of barriers to reproduction. The study by Rafati et al. () in this issue of Molecular Ecology combines spatially explicit whole-genome resequencing with evaluation of differential gene expression across individuals with mixed ancestry to associate the genomic architecture of reproductive barriers with expression of reproductive incompatibilities. In a natural hybrid zone between rabbit subspecies, Oryctolagus cuniculus cuniculus and O. c. algirus (Figure ), Rafati et al. () use landscape-level patterns of allele frequency variation to identify potential candidate regions of the genome associated with reproductive isolation. These candidate regions are used to test predictions associated with the genomic architecture of reproductive barriers, including the role of structural rearrangements, enrichment of functional categories associated with incompatibilities, and the contribution of protein-coding versus regulatory changes. A lack of structural rearrangements and limited protein-coding changes in candidate regions point towards the importance of regulatory variation as major contributors to genetic incompatibilities, while functional enrichments indicate overrepresentation of genes associated with male infertility. To quantify phenotypic expression of proposed incompatibilities, the authors assess gene expression of experimental crosses. Extensive misregulation of gene expression within the testes of backcross hybrids relative to F1 and parental individuals provides an important link between genotype and phenotype, validating hypotheses developed from assessment of genomic architectures. Together, this work shows how pairing natural hybrid zones with experimental crosses can be used to link observations in nature to mechanistic underpinnings that may be tested experimentally.

Genetic Pedigree Analysis of the Pilot Breeding Program for the Rediscovered Galapagos Giant Tortoise from Floreana Island.

An aim of many captive breeding programs is to increase population sizes for reintroduction and establishment of self-sustaining wild populations. Genetic analyses play a critical role in these programs: monitoring genetic variation, identifying the origin of individuals, and assigning parentage to track family sizes. Here, we use genetic pedigree analyses to examine 3 seasons of a pilot breeding program for the Floreana island Galapagos giant tortoise, C. niger, that had been declared extinct for ~150 years until individuals with mixed ancestry were recently discovered. We determined that 8 of 9 founding individuals were assigned parentage to at least 1 of 130 offspring produced, though there was considerable reproductive skew. In addition, we observed that genetic diversity of the progeny was lower than that of the founders. Despite the observed reproductive skew, we did not see evidence for assortative mating based on relatedness, but there was a trend toward reduced fitness when more related individuals bred. Finally, we found that the majority of progeny had ancestry assigned to the Floreana species (mean ± SE = 0.51 ± 0.02), though individual estimates varied. The success of these pilot seasons bodes well for a larger breeding program to help restore the previously extinct tortoise from Floreana island. Future efforts should continue to monitor for reproductive skew and assortative mating to maintain allelic diversity. We would also recommend forming smaller breeding groups and rotating individuals among them to prevent long-term reproductive skew among pairs.

Temporal Mitogenomics of the Galapagos Giant Tortoise from Pinzón Reveals Potential Biases in Population Genetic Inference.

Empirical population genetic studies generally rely on sampling subsets of the population(s) of interest and of the nuclear or organellar genome targeted, assuming each is representative of the whole. Violations of these assumptions may impact population-level parameter estimation and lead to spurious inferences. Here, we used targeted capture to sequence the full mitochondrial genome from 123 individuals of the Galapagos giant tortoise endemic to Pinzón Island (Chelonoidis duncanensis) sampled at 2 time points pre- and postbottleneck (circa 1906 and 2014) to explicitly assess differences in diversity estimates and demographic reconstructions based on subsets of the mitochondrial genome versus the full sequences and to evaluate potential biases associated with diversity estimates and demographic reconstructions from postbottlenecked samples alone. Haplotypic diversities were equal between the temporal samples based on the full mitochondrial genome, but single gene estimates suggested either decreases or increases in diversity depending upon the region. Demographic reconstructions based on the full sequence were more similar between the temporal samples than those based on the control region alone, or a subset of 3 regions, where the trends in population size changes shifted in magnitude and direction between the temporal samples. In all cases, the estimated coalescent point was more distant for the historical than contemporary sample. In summary, our results empirically demonstrate the influence of sampling bias when interpreting population genetic patterns and punctuate the need for careful consideration of potentially conflicting evolutionary signal across the mitochondrial genome.

Genome-Wide Assessment of Diversity and Divergence Among Extant Galapagos Giant Tortoise Species.

Genome-wide assessments allow for fuller characterization of genetic diversity, finer-scale population delineation, and better detection of demographically significant units to guide conservation compared with those based on "traditional" markers. Galapagos giant tortoises (Chelonoidis spp.) have long provided a case study for how evolutionary genetics may be applied to advance species conservation. Ongoing efforts to bolster tortoise populations, which have declined by 90%, have been informed by analyses of mitochondrial DNA sequence and microsatellite genotypic data, but could benefit from genome-wide markers. Taking this next step, we used double-digest restriction-site associated DNA sequencing to collect genotypic data at >26000 single nucleotide polymorphisms (SNPs) for 117 individuals representing all recognized extant Galapagos giant tortoise species. We then quantified genetic diversity, population structure, and compared results to estimates from mitochondrial DNA and microsatellite loci. Our analyses detected 12 genetic lineages concordant with the 11 named species as well as previously described structure within one species, C. becki. Furthermore, the SNPs provided increased resolution, detecting admixture in 4 individuals. SNP-based estimates of diversity and differentiation were significantly correlated with those derived from nuclear microsatellite loci and mitochondrial DNA sequences. The SNP toolkit presented here will serve as a resource for advancing efforts to understand tortoise evolution, species radiations, and aid conservation of the Galapagos tortoise species complex.

The K = 2 conundrum.

Assessments of population genetic structure have become an increasing focus as they can provide valuable insight into patterns of migration and gene flow. structure, the most highly cited of several clustering-based methods, was developed to provide robust estimates without the need for populations to be determined a priori. structure introduces the problem of selecting the optimal number of clusters, and as a result, the ΔK method was proposed to assist in the identification of the "true" number of clusters. In our review of 1,264 studies using structure to explore population subdivision, studies that used ΔK were more likely to identify K = 2 (54%, 443/822) than studies that did not use ΔK (21%, 82/386). A troubling finding was that very few studies performed the hierarchical analysis recommended by the authors of both ΔK and structure to fully explore population subdivision. Furthermore, extensions of earlier simulations indicate that, with a representative number of markers, ΔK frequently identifies K = 2 as the top level of hierarchical structure, even when more subpopulations are present. This review suggests that many studies may have been over- or underestimating population genetic structure; both scenarios have serious consequences, particularly with respect to conservation and management. We recommend publication standards for population structure results so that readers can assess the implications of the results given their own understanding of the species biology.

Adaptive introgression as a resource for management and genetic conservation in a changing climate.

Current rates of climate change require organisms to respond through migration, phenotypic plasticity, or genetic changes via adaptation. We focused on questions regarding species' and populations' ability to respond to climate change through adaptation. Specifically, the role adaptive introgression, movement of genetic material from the genome of 1 species into the genome of another through repeated interbreeding, may play in increasing species' ability to respond to a changing climate. Such interspecific gene flow may mediate extinction risk or consequences of limited adaptive potential that result from standing genetic variation and mutation alone, enabling a quicker demographic recovery in response to changing environments. Despite the near dismissal of the potential benefits of hybridization by conservation practitioners, we examined a number of case studies across different taxa that suggest gene flow between sympatric or parapatric sister species or within species that exhibit strong ecotypic differentiation may represent an underutilized management option to conserve evolutionary potential in a changing environment. This will be particularly true where advanced-generation hybrids exhibit adaptive traits outside the parental phenotypic range, a phenomenon known as transgressive segregation. The ideas presented in this essay are meant to provoke discussion regarding how we maintain evolutionary potential, the conservation value of natural hybrid zones, and consideration of their important role in adaptation to climate.

Cross-Species Application of SNP Chips is Not Suitable for Identifying Runs of Homozygosity.

Cross-species application of single-nucleotide polymorphism (SNP) chips is a valid, relatively cost-effective alternative to the high-throughput sequencing methods generally required to obtain a genome-wide sampling of polymorphisms. Kharzinova et al. (2015) examined the applicability of SNP chips developed in domestic bovids (cattle and sheep) to a semi-wild cervid (reindeer). The ancestors of bovids and cervids diverged between 20 and 30 million years ago (Hassanin and Douzery 2003; Bibi et al. 2013). Empirical work has shown that for a SNP chip developed in a bovid and applied to a cervid species, approximately 50% genotype success with 1% of the loci being polymorphic is expected (Miller et al. 2012). The genotyping of Kharzinova et al. (2015) follows this pattern; however, these data are not appropriate for identifying runs of homozygosity (ROH) and can be problematic for estimating linkage disequilibrium (LD) and we caution readers in this regard.

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.