Targeted resequencing of coding DNA sequences for SNP discovery in nonmodel species

Comparative assessment of genotyping-by-sequencing and whole-exome sequencing for estimating genetic diversity and geographic structure in small sample sizes: insights from wild jaguar populations

Biologists currently have an assortment of high-throughput sequencing techniques allowing the study of population dynamics in increasing detail. The utility of genetic estimates depends on their ability to recover meaningful approximations while filtering out noise produced by artifacts. In this study, we empirically compared the congruence of two reduced representation approaches (genotyping-by-sequencing, GBS, and whole-exome sequencing, WES) in estimating genetic diversity and population structure using SNP markers typed in a small number of wild jaguar (Panthera onca) samples from South America. Due to its targeted nature, WES allowed for a more straightforward reconstruction of loci compared to GBS, facilitating the identification of true polymorphisms across individuals. We therefore used WES-derived metrics as a benchmark against which GBS-derived indicators were compared, adjusting parameters for locus assembly and SNP filtering in the latter. We observed significant variation in SNP call rates across samples in GBS datasets, leading to a recurrent miscalling of heterozygous sites. This issue was further amplified by small sample sizes, ultimately impacting the consistency of summary statistics between genotyping methods. Recognizing that the genetic markers obtained from GBS and WES are intrinsically different due to varying evolutionary pressures, particularly selection, we consider that our empirical comparison offers valuable insights and highlights critical considerations for estimating population genetic attributes using reduced representation datasets. Our results emphasize the critical need for careful evaluation of missing data and stringent filtering to achieve reliable estimates of genetic diversity and differentiation in elusive wildlife species.

Coevolution with hosts underpins speciation in brood-parasitic cuckoos

Coevolution between interacting species is thought to increase biodiversity, but evidence linking microevolutionary processes to macroevolutionary patterns is scarce. We leveraged two decades of behavioral research coupled with historical DNA analysis to reveal that coevolution with hosts underpins speciation in brood-parasitic bronze-cuckoos. At a macroevolutionary scale, we show that highly virulent brood-parasitic taxa have higher speciation rates and are more likely to speciate in sympatry than less-virulent and nonparasitic relatives. We reveal the microevolutionary process underlying speciation: Hosts reject cuckoo nestlings, which selects for mimetic cuckoo nestling morphology. Where cuckoos exploit multiple hosts, selection for mimicry drives genetic and phenotypic divergence corresponding to host preference, even in sympatry. Our work elucidates perhaps the most common, but poorly characterized, evolutionary process driving biological diversification.

Whole-genome resequencing reveals the uniqueness of Subei yak

Subei yak is an essential local yak in the Gansu Province, which genetic resource has recently been discovered. It is a meat-milk dual-purpose variety with high fecundity and relatively stable population genetic structure. However, its population genetic structure and genetic diversity are yet to be reported. Therefore, this study aimed to identify molecular markers of Subei yak genome by whole-genome resequencing, and to analyze the population structure and genetic diversity of Subei yak. This study screened 12,079,496 single nucleotide polymorphism (SNP) molecular markers in the 20 Subei yaks genome using whole-genome resequencing technology. Of these SNPs, 32.09% were located in the intronic region of the genome. Principal component analysis, phylogenetic analysis, and population structure analysis revealed that the Subei yak belonged to an independent group in the domestic yak population. A selective clearance analysis was carried out on Subei yak and other domestic yaks, and the genes under positive selection were annotated. The functional enrichment analysis showed that Subei yak possessed prominent selection characteristics in terms of external environment perception, hypoxia adaptation, and muscle development. Furthermore, Subei yak showed excellent muscle fat deposition and meat quality traits. Thus, this study will serve as a reference for discovering population structure, genetic evolution, and other unique traits of Subei yak and for expanding the genetic variation catalog of yaks.

Comparative genomics of the waterfowl innate immune system

Animal species differ considerably in their ability to fight off infections. Finding the genetic basis of these differences is not easy, as the immune response is comprised of a complex network of proteins that interact with one another to defend the body against infection. Here, we used population- and comparative genomics to study the evolutionary forces acting on the innate immune system in natural hosts of the avian influenza virus (AIV). For this purpose, we used a combination of hybrid capture, next- generation sequencing and published genomes to examine genetic diversity, divergence, and signatures of selection in 127 innate immune genes at a micro- and macroevolutionary time scale in 26 species of waterfowl. We show across multiple immune pathways (AIV-, toll-like-, and RIG-I -like receptors signalling pathways) that genes involved genes in pathogen detection (i.e., toll-like receptors) and direct pathogen inhibition (i.e., antimicrobial peptides and interferon-stimulated genes), as well as host proteins targeted by viral antagonist proteins (i.e., mitochondrial antiviral-signaling protein, [MAVS]) are more likely to be polymorphic, genetically divergent, and under positive selection than other innate immune genes. Our results demonstrate that selective forces vary across innate immune signaling signalling pathways in waterfowl, and we present candidate genes that may contribute to differences in susceptibility and resistance to infectious diseases in wild birds, and that may be manipulated by viruses. Our findings improve our understanding of the interplay between host genetics and pathogens, and offer the opportunity for new insights into pathogenesis and potential drug targets.

Coat Polymorphism in Eurasian Lynx: Adaptation to Environment or Phylogeographic Legacy?

We studied the relationship between the variability and contemporary distribution of pelage phenotypes in one of most widely distributed felid species and an array of environmental and demographic conditions. We collected 672 photographic georeferenced records of the Eurasian lynx throughout Eurasia. We assigned each lynx coat to one of five phenotypes. Then we fitted the coat patterns to different environmental and anthropogenic variables, as well as the effective geographic distances from inferred glacial refugia. A majority of lynx were either of the large spotted (41.5%) or unspotted (uniform, 36.2%) phenotype. The remaining patterns (rosettes, small spots and pseudo-rosettes) were represented in 11.0%, 7.4%, and 3.9% of samples, respectively. Although various environmental variables greatly affected lynx distribution and habitat suitability, it was the effect of least-cost distances from locations of the inferred refugia during the Last Glacial Maximum that explained the distribution of lynx coat patterns the best. Whereas the occurrence of lynx phenotypes with large spots was explained by the proximity to refugia located in the Caucasus/Middle East, the uniform phenotype was associated with refugia in the Far East and Central Asia. Despite the widely accepted hypothesis of adaptive functionality of coat patterns in mammals and exceptionally high phenotypic polymorphism in Eurasian lynx, we did not find well-defined signs of habitat matching in the coat pattern of this species. Instead, we showed how the global patterns of morphological variability in this large mammal and its environmental adaptations may have been shaped by past climatic change.

Capture by hybridization for full-length barcode-based eukaryotic and prokaryotic biodiversity inventories of deep sea ecosystems

Biodiversity inventory of marine systems remains limited due to unbalanced access to the three ocean dimensions. The use of environmental DNA (eDNA) for metabarcoding allows fast and effective biodiversity inventory and is forecast as a future biodiversity research and biomonitoring tool. However, in poorly understood ecosystems, eDNA results remain difficult to interpret due to large gaps in reference databases and PCR bias limiting the detection of some major phyla. Here, we aimed to circumvent these limitations by avoiding PCR and recollecting larger DNA fragments to improve assignment of detected taxa through phylogenetic reconstruction. We applied capture by hybridization (CBH) to enrich DNA from deep-sea sediment samples and compared the results with those obtained through an up-to-date metabarcoding PCR-based approach (MTB). Originally developed for bacterial communities and targeting 16S rDNA, the CBH approach was applied to 18S rDNA to improve the detection of species forming benthic communities of eukaryotes, with a particular focus on metazoans. The results confirmed the possibility of extending CBH to metazoans with two major advantages: (i) CBH revealed a broader spectrum of prokaryotic, eukaryotic, and particularly metazoan diversity, and (ii) CBH allowed much more robust phylogenetic reconstructions of full-length barcodes with up to 1900 base pairs. This is particularly important for taxa whose assignment is hampered by gaps in reference databases. This study provides a database and probes to apply 18S CBH to diverse marine systems, confirming this promising new tool to improve biodiversity assessments in data-poor ecosystems such as those in the deep sea.

Multiple species-specific molecular markers using nanofluidic array as a tool to detect prey DNA from carnivore scats

Large carnivore feeding ecology plays a crucial role for management and conservation for predators and their prey. One of the keys to this kind of research is to identify the species composition in the predator diet, for example, prey determination from scat content. DNA-based methods applied to detect prey in predators' scats are viable alternatives to traditional macroscopic approaches, showing an increased reliability and higher prey detection rate. Here, we developed a molecular method for prey species identification in wolf (Canis lupus) scats using multiple species-specific marker loci on the cytochrome b gene for 18 target species. The final panel consisted of 80 assays, with a minimum of four markers per target species, and that amplified specifically when using a high-throughput Nanofluidic array technology (Fluidigm Inc.). As a practical example, we applied the method to identify target prey species DNA in 80 wolf scats collected in Sweden. Depending on the number of amplifying markers required to obtain a positive species call in a scat, the success in determining at least one prey species from the scats ranged from 44% to 92%. Although we highlight the need to evaluate the optimal number of markers for sensitive target species detection, the developed method is a fast and cost-efficient tool for prey identification in wolf scats and it also has the potential to be further developed and applied to other areas and large carnivores as well.

Mitogenome Phylogeny Including Data from Additional Subspecies Provides New Insights into the Historical Biogeography of the Eurasian lynx Lynx lynx

Previous molecular studies of the wide-ranging Eurasian lynx Lynx lynx focused mainly on its northern Palearctic populations, with the consequence that the reconstruction of this species’ evolutionary history did not include genetic variation present in its southern Palearctic distribution. We sampled a previously not considered Asian subspecies (L. l. dinniki), added published data from another Asian subspecies (L. l. isabellinus), and reassessed the Eurasian lynx mtDNA phylogeny along with previously published data from northern Palearctic populations. Our mitogenome-based analyses revealed the existence of three major clades (A: Central Asia, B: SE Europe/SW Asia, C: Europe and Northern Asia) and at least five lineages, with diversification in Lynx lynx commencing at least 28kyr earlier than hitherto estimated. The subspecies L. l. isabellinus harbors the most basal matriline, consistent with the origin of Lynx lynx in this subspecies’ current range. L. l. dinniki harbors the second most basal matriline, which is related to, and may be the source of, the mtDNA diversity of the critically endangered Balkan lynx L. l. balcanicus. Our results suggest that the Anatolian peninsula was a glacial refugium for Eurasian lynx, with previously unconsidered implications for the colonization of Europe by this species.

Sample identification and pedigree reconstruction in Wolverine (Gulo gulo) using SNP genotyping of non-invasive samples

For conservation genetic studies using non-invasively collected samples, genome-wide data may be hard to acquire. Until now, such studies have instead mostly relied on analyses of traditional genetic markers such as microsatellites (SSRs). Recently, high throughput genotyping of single nucleotide polymorphisms (SNPs) has become available, expanding the use of genomic methods to include non-model species of conservation concern. We have developed a 96-marker SNP array for use in applied conservation monitoring of the Scandinavian wolverine (Gulo gulo) population. By genotyping more than a thousand non-invasively collected samples, we were able to obtain precise estimates of different types of genotyping errors and sample dropout rates. The SNP panel significantly outperforms the SSR markers (and DBY intron markers for sexing) both in terms of precision in genotyping, sex assignment and individual identification, as well as in the proportion of samples successfully genotyped. Furthermore, SNP genotyping offers a simplified laboratory and analysis pipeline with fewer samples needed to be repeatedly genotyped in order to obtain reliable consensus data. In addition, we utilised a unique opportunity to successfully demonstrate the application of SNP genotype data for reconstructing pedigrees in wild populations, by validating the method with samples from wild individuals with known relatedness. By offering a simplified workflow with improved performance, we anticipate this methodology will facilitate the use of non-invasive samples to improve genetic management of many different types of populations that have previously been challenging to survey.

Development of a Panel of Genotyping-in-Thousands by Sequencing in Capsicum

Genotyping by sequencing (GBS) enables genotyping of multiple loci at low cost. However, the single nucleotide polymorphisms (SNPs) revealed by GBS tend to be randomly distributed between individuals, limiting their direct comparisons without applying the various filter options to obtain a comparable dataset of SNPs. Here, we developed a panel of a multiplex targeted sequencing method, genotyping-in-thousands by sequencing (GT-seq), to genotype SNPs in Capsicum spp. Previously developed Fluidigm^® SNP markers were converted to GT-seq markers and combined with new GT-seq markers developed using SNP information obtained through GBS. We then optimized multiplex PCR conditions: we obtained the highest genotyping rate when the first PCR consisted of 25 cycles. In addition, we determined that 101 primer pairs performed best when amplifying target sequences of 79 bp. We minimized interference of multiplex PCR by primer dimer formation using the PrimerPooler program. Using our GT-seq pipeline on Illumina Miseq and Nextseq platforms, we genotyped up to 1,500 (Miseq) and 1,300 (Nextseq) samples for the optimum panel size of 100 loci. To allow the genotyping of Capsicum species, we designed 332 informative GT-seq markers from Fluidigm SNP markers and GBS-derived SNPs. This study illustrates the first application of GT-seq in crop plants. The GT-seq marker set developed here will be a useful tool for molecular breeding of peppers in the future.

Population genomics for wildlife conservation and management

Biodiversity is under threat worldwide. Over the past decade, the field of population genomics has developed across nonmodel organisms, and the results of this research have begun to be applied in conservation and management of wildlife species. Genomics tools can provide precise estimates of basic features of wildlife populations, such as effective population size, inbreeding, demographic history and population structure, that are critical for conservation efforts. Moreover, population genomics studies can identify particular genetic loci and variants responsible for inbreeding depression or adaptation to changing environments, allowing for conservation efforts to estimate the capacity of populations to evolve and adapt in response to environmental change and to manage for adaptive variation. While connections from basic research to applied wildlife conservation have been slow to develop, these connections are increasingly strengthening. Here we review the primary areas in which population genomics approaches can be applied to wildlife conservation and management, highlight examples of how they have been used, and provide recommendations for building on the progress that has been made in this field.

Ancestral mitogenome capture of the Southeast Asian banded linsang

Utilising a reconstructed ancestral mitochondrial genome of a clade to design hybridisation capture baits can provide the opportunity for recovering mitochondrial sequences from all its descendent and even sister lineages. This approach is useful for taxa with no extant close relatives, as is often the case for rare or extinct species, and is a viable approach for the analysis of historical museum specimens. Asiatic linsangs (genus Prionodon) exemplify this situation, being rare Southeast Asian carnivores for which little molecular data is available. Using ancestral capture we recover partial mitochondrial genome sequences for seven banded linsangs (P. linsang) from historical specimens, representing the first intraspecific genetic dataset for this species. We additionally assemble a high quality mitogenome for the banded linsang using shotgun sequencing for time-calibrated phylogenetic analysis. This reveals a deep divergence between the two Asiatic linsang species (P. linsang, P. pardicolor), with an estimated divergence of ~12 million years (Ma). Although our sample size precludes any robust interpretation of the population structure of the banded linsang, we recover two distinct matrilines with an estimated tMRCA of ~1 Ma. Our results can be used as a basis for further investigation of the Asiatic linsangs, and further demonstrate the utility of ancestral capture for studying divergent taxa without close relatives.

Genomic patterns in the widespread Eurasian lynx shaped by Late Quaternary climatic fluctuations and anthropogenic impacts

Disentangling the contribution of long-term evolutionary processes and recent anthropogenic impacts to current genetic patterns of wildlife species is key to assessing genetic risks and designing conservation strategies. Here, we used 80 whole nuclear genomes and 96 mitogenomes from populations of the Eurasian lynx covering a range of conservation statuses, climatic zones and subspecies across Eurasia to infer the demographic history, reconstruct genetic patterns, and discuss the influence of long-term isolation and/or more recent human-driven changes. Our results show that Eurasian lynx populations shared a common history until 100,000 years ago, when Asian and European populations started to diverge and both entered a period of continuous and widespread decline, with western populations, except Kirov, maintaining lower effective sizes than eastern populations. Population declines and increased isolation in more recent times probably drove the genetic differentiation between geographically and ecologically close westernmost European populations. By contrast, and despite the wide range of habitats covered, populations are quite homogeneous genetically across the Asian range, showing a pattern of isolation by distance and providing little genetic support for the several proposed subspecies. Mitogenomic and nuclear divergences and population declines starting during the Late Pleistocene can be mostly attributed to climatic fluctuations and early human influence, but the widespread and sustained decline since the Holocene is more probably the consequence of anthropogenic impacts which intensified in recent centuries, especially in western Europe. Genetic erosion in isolated European populations and lack of evidence for long-term isolation argue for the restoration of lost population connectivity.

Novel genomic resources for shelled pteropods: a draft genome and target capture probes for Limacina bulimoides, tested for cross-species relevance

BACKGROUND

Pteropods are planktonic gastropods that are considered as bio-indicators to monitor impacts of ocean acidification on marine ecosystems. In order to gain insight into their adaptive potential to future environmental changes, it is critical to use adequate molecular tools to delimit species and population boundaries and to assess their genetic connectivity. We developed a set of target capture probes to investigate genetic variation across their large-sized genome using a population genomics approach. Target capture is less limited by DNA amount and quality than other genome-reduced representation protocols, and has the potential for application on closely related species based on probes designed from one species.

RESULTS

We generated the first draft genome of a pteropod, Limacina bulimoides, resulting in a fragmented assembly of 2.9 Gbp. Using this assembly and a transcriptome as a reference, we designed a set of 2899 genome-wide target capture probes for L. bulimoides. The set of probes includes 2812 single copy nuclear targets, the 28S rDNA sequence, ten mitochondrial genes, 35 candidate biomineralisation genes, and 41 non-coding regions. The capture reaction performed with these probes was highly efficient with 97% of the targets recovered on the focal species. A total of 137,938 single nucleotide polymorphism markers were obtained from the captured sequences across a test panel of nine individuals. The probes set was also tested on four related species: L. trochiformis, L. lesueurii, L. helicina, and Heliconoides inflatus, showing an exponential decrease in capture efficiency with increased genetic distance from the focal species. Sixty-two targets were sufficiently conserved to be recovered consistently across all five species.

CONCLUSION

The target capture protocol used in this study was effective in capturing genome-wide variation in the focal species L. bulimoides, suitable for population genomic analyses, while providing insights into conserved genomic regions in related species. The present study provides new genomic resources for pteropods and supports the use of target capture-based protocols to efficiently characterise genomic variation in small non-model organisms with large genomes.

Terrestrial mammal surveillance using hybridization capture of environmental DNA from African waterholes

Determining species distributions can be extremely challenging but is crucial to ecological and conservation research. Environmental DNA (eDNA) approaches have shown particular promise in aquatic systems for several vertebrate and invertebrate species. For terrestrial animals, however, eDNA-based surveys are considerably more difficult due to the lack of or difficulty in obtaining appropriate sampling substrate. In water-limited ecosystem where terrestrial mammals are often forced to congregate at waterholes, water and sediment from shared water sources may be a suitable substrate for noninvasive eDNA approaches. We characterized mitochondrial DNA sequences from a broad range of terrestrial mammal species in two different African ecosystems (in Namibia and Tanzania) using eDNA isolated from native water, sediment and water filtered through glass fibre filters. A hybridization capture enrichment with RNA probes targeting the mitochondrial genomes of 38 mammal species representing the genera/families expected at the respective ecosystems was employed, and 16 species were identified, with a maximum mitogenome coverage of 99.8%. Conventional genus-specific PCRs were tested on environmental samples for two genera producing fewer positive results than hybridization capture enrichment. An experiment with mock samples using DNA from non-African mammals showed that baits covering 30% of nontarget mitogenomes produced 91% mitogenome coverage after capture. In the mock samples, over-representation of DNA of one species still allowed for the detection of DNA of other species that was at a 100-fold lower concentration. Hybridization capture enrichment of eDNA is therefore an effective method for monitoring terrestrial mammal species from shared water sources.