Fast and cost-effective single nucleotide polymorphism (SNP) detection in the absence of a reference genome using semideep next-generation Random Amplicon Sequencing (RAMseq)

Population admixtures in medaka inferred by multiple arbitrary amplicon sequencing

Cost-effective genotyping can be achieved by sequencing PCR amplicons. Short 3-10 base primers can arbitrarily amplify thousands of loci using only a few primers. To improve the sequencing efficiency of the multiple arbitrary amplicon sequencing (MAAS) approach, we designed new primers and examined their efficiency in sequencing and genotyping. To demonstrate the effectiveness of our method, we applied it to examining the population structure of the small freshwater fish, medaka (Oryzias latipes). We obtained 2987 informative SNVs with no missing genotype calls for 67 individuals from 15 wild populations and three artificial strains. The estimated phylogenic and population genetic structures of the wild populations were consistent with previous studies, corroborating the accuracy of our genotyping method. We also attempted to reconstruct the genetic backgrounds of a commercial orange mutant strain, Himedaka, which has caused a genetic disturbance in wild populations. Our admixture analysis focusing on Himedaka showed that at least two wild populations had genetically been contributed to the nuclear genome of this mutant strain. Our genotyping methods and results will be useful in quantitative assessments of genetic disturbance by this commercially available strain.

Multi-year pair-bonding in Murray cod (Maccullochella peelii)

Mating strategies in fishes are known to include polygyny, polyandry and monogamy and provide valuable insights regarding powerful evolutionary forces such as sexual selection. Monogamy is a complex of mating systems that has been relatively neglected. Previous work on mating strategies in fishes has often been based on observation and focused on marine species rather than freshwater fishes. SNPs are increasingly being used as a molecular ecology tool in non-model organisms, and methods of probabilistic genetic analysis of such datasets are becoming available for use in the absence of parental genotypes. This approach can be used to infer mating strategies. The long-term pair bonding seen in mammals, reptiles and birds has not been recorded in freshwater fishes-in every other respect an extremely diverse group. This study shows that multi-year pair bonding occurs in an Australian Percichthyid fish that exhibits paternal care of eggs and larvae. Using SNPs, full sibling pairs of larvae were found over multiple years in a three-year study. Stable isotope signatures of the larvae support the genetic inference that full sibling pairs shared a common mother, the ultimate source of that isotopic signature during oogenesis. Spatial and temporal clustering also suggests that the full sibling larvae are unlikely to be false positive identifications of the probabilistic identification of siblings. For the first time, we show multi-year pair bonding in a wild freshwater fish. This will have important conservation and management implications for the species. This approach could provide insights into many behavioural, ecological and evolutionary questions, particularly if this is not a unique case. Our findings are likely to initiate interest in seeking more examples of monogamy and alternative mating strategies in freshwater fishes, particularly if others improve methods of analysis of SNP data for identification of siblings in the absence of parental genotypes.

Using DNA metabarcoding to investigate diet and niche partitioning in the native European otter (Lutra lutra) and invasive American mink (Neovison vison)

In the UK, the native European otter (Lutra lutra) and invasive American mink (Neovison vison) have experienced concurrent declines and expansions. Currently, the otter is recovering from persecution and waterway pollution, whereas the mink is in decline due to population control and probable interspecific interaction with the otter. We explored the potential of DNA metabarcoding for investigating diet and niche partitioning between these mustelids. Otter spraints (n = 171) and mink scats (n = 19) collected from three sites (Malham Tarn, River Hull and River Glaven) in northern and eastern England were screened for vertebrates using high-throughput sequencing. Otter diet mainly comprised aquatic fishes (81.0%) and amphibians (12.7%), whereas mink diet predominantly consisted of terrestrial birds (55.9%) and mammals (39.6%). The mink used a lower proportion (20%) of available prey (n = 40 taxa) than the otter and low niche overlap (0.267) was observed between these mustelids. Prey taxon richness of mink scats was lower than otter spraints and beta diversity of prey communities was driven by taxon turnover (i.e. the otter and mink consumed different prey taxa). Considering otter diet only, prey taxon richness was higher in spraints from the River Hull catchment and beta diversity of prey communities was driven by taxon turnover (i.e. the otter consumed different prey taxa at each site). Studies using morphological faecal analysis may misidentify the predator as well as prey items. Faecal DNA metabarcoding can resolve these issues and provide more accurate and detailed dietary information. When scaled up across multiple habitat types, DNA metabarcoding should greatly improve future understanding of resource use and niche overlap between the otter and mink.

All-in-one approaches for rapid and highly specific quantifcation of single nucleotide polymorphisms based on ligase detection reaction using molecular beacons as turn-on probes

Rapid, simple, specific and sensitive approaches for single nucleotide polymorphisms (SNPs) detection are essential for clinical diagnosis. In this study, all-in-one approaches, consisting of the whole detection process including ligase detection reaction (LDR) and real time quantitative polymerase chain reaction performed in one PCR tube by a one-step operation on a real-time PCR system using molecular beacon (MB) as turn-on probe, were developed for rapid, simple, specific and sensitive quantifcation of SNPs. High specificity of the all-in-one approach was achieved by using the LDR, which employs a thermostable and single-base discerning Hifi Taq DNA ligase to ligate adjacently hybridized LDR-specific probes. In addition, a highly specific probe, MB, was used to detect the products of all-in-one approach, which doubly enhances the specificity of the all-in-one approach. The linear dynamic range and high sensitivity of mutant DNA (MutDNA) and wild-type DNA (WtDNA) all-in-one approaches for the detection of MutDNA and WtDNA were studied in vitro, with a broad linear dynamic range of 0.1 fM to 1 pM and detection limits of 65.3 aM and 31.2 aM, respectively. In addition, the MutDNA and WtDNA all-in-one approaches were able to accurately detect allele frequency changes as low as 0.1%. In particular, the epidermal growth factor receptor T790M MutDNA frequency in the tissue of five patients with non-small cell lung cancer detected by all-in-one approaches were in agreement with clinical detection results, indicating the excellent practicability of the developed approaches for the quantification of SNPs in real samples. In summary, the developed all-in-one approaches exhibited promising potential for further applications in clinical diagnosis.

Developing and quality testing of microsatellite loci for four species of Glossina

Microsatellite loci still represent valuable resources for the study of the population biology of non-model organisms. Discovering or adapting new suitable microsatellite markers in species of interest still represents a useful task, especially so for non-model organisms as tsetse flies (genus Glossina), which remain a serious threat to the health of humans and animals in sub-Saharan Africa. In this paper, we present the development of new microsatellite loci for four species of Glossina: two from the Morsitans group, G. morsitans morsitans (Gmm) from Zimbabwe, G. pallidipes (Gpalli) from Tanzania; and the other two from the Palpalis group, G. fuscipes fuscipes (Gff) from Chad, and G. palpalis gambiensis (Gpg) from Guinea. We found frequent short allele dominance and null alleles. Stuttering could also be found and amended when possible. Cryptic species seemed to occur frequently in all taxa but Gff. This explains why it may be difficult finding ecumenical primers, which thus need adaptation according to each taxonomic and geographic context. Amplification problems occurred more often in published old markers, and Gmm and Gpg were the most affected (stronger heterozygote deficits). Trinucleotide markers displayed selection signature in some instances (Gmm). Combining old and new loci, for Gmm, eight loci can be safely used (with correction for null alleles); and five seem particularly promising; for Gpalli, only five to three loci worked well, depending on the clade, which means that the use of loci from other species (four morsitans loci seemed to work well), or other new primers will need to be used; for Gff, 14 loci behaved well, but with null alleles, seven of which worked very well; and for G. palpalis sl, only four loci, needing null allele and stuttering corrections seem to work well, and other loci from the literature are thus needed, including X-linked markers, five of which seem to work rather well (in females only), but new markers will probably be needed. Finally, the high proportion of X-linked markers (around 30%) was explained by the non-Y DNA quantity and chromosome structure of tsetse flies studied so far.

Applying genomic data in wildlife monitoring: Development guidelines for genotyping degraded samples with reduced single nucleotide polymorphism panels

The genomic era has led to an unprecedented increase in the availability of genome-wide data for a broad range of taxa. Wildlife management strives to make use of these vast resources to enable refined genetic assessments that enhance biodiversity conservation. However, as new genomic platforms emerge, problems remain in adapting the usually complex approaches for genotyping of noninvasively collected wildlife samples. Here, we provide practical guidelines for the standardized development of reduced single nucleotide polymorphism (SNP) panels applicable for microfluidic genotyping of degraded DNA samples, such as faeces or hairs. We demonstrate how microfluidic SNP panels can be optimized to efficiently monitor European wildcat (Felis silvestris S.) populations. We show how panels can be set up in a modular fashion to accommodate informative markers for relevant population genetics questions, such as individual identification, hybridization assessment and the detection of population structure. We discuss various aspects regarding the implementation of reduced SNP panels and provide a framework that will allow both molecular ecologists and practitioners to help bridge the gap between genomics and applied wildlife conservation.

A Fast and Scalable Workflow for SNPs Detection in Genome Sequences Using Hadoop Map-Reduce

Next generation sequencing (NGS) technologies produce a huge amount of biological data, which poses various issues such as requirements of high processing time and large memory. This research focuses on the detection of single nucleotide polymorphism (SNP) in genome sequences. Currently, SNPs detection algorithms face several issues, e.g., computational overhead cost, accuracy, and memory requirements. In this research, we propose a fast and scalable workflow that integrates Bowtie aligner with Hadoop based Heap SNP caller to improve the SNPs detection in genome sequences. The proposed workflow is validated through benchmark datasets obtained from publicly available web-portals, e.g., NCBI and DDBJ DRA. Extensive experiments have been performed and the results obtained are compared with Bowtie and BWA aligner in the alignment phase, while compared with GATK, FaSD, SparkGA, Halvade, and Heap in SNP calling phase. Experimental results analysis shows that the proposed workflow outperforms existing frameworks e.g., GATK, FaSD, Heap integrated with BWA and Bowtie aligners, SparkGA, and Halvade. The proposed framework achieved 22.46% more efficient F-score and 99.80% consistent accuracy on average. More, comparatively 0.21% mean higher accuracy is achieved. Moreover, SNP mining has also been performed to identify specific regions in genome sequences. All the frameworks are implemented with the default configuration of memory management. The observations show that all workflows have approximately same memory requirement. In the future, it is intended to graphically show the mined SNPs for user-friendly interaction, analyze and optimize the memory requirements as well.

Applying genomic data in wildlife monitoring: Development guidelines for genotyping degraded samples with reduced single nucleotide polymorphism panels

The genomic era has led to an unprecedented increase in the availability of genome-wide data for a broad range of taxa. Wildlife management strives to make use of these vast resources to enable refined genetic assessments that enhance biodiversity conservation. However, as new genomic platforms emerge, problems remain in adapting the usually complex approaches for genotyping of noninvasively collected wildlife samples. Here, we provide practical guidelines for the standardized development of reduced single nucleotide polymorphism (SNP) panels applicable for microfluidic genotyping of degraded DNA samples, such as faeces or hairs. We demonstrate how microfluidic SNP panels can be optimized to efficiently monitor European wildcat (Felis silvestris S.) populations. We show how panels can be set up in a modular fashion to accommodate informative markers for relevant population genetics questions, such as individual identification, hybridization assessment and the detection of population structure. We discuss various aspects regarding the implementation of reduced SNP panels and provide a framework that will allow both molecular ecologists and practitioners to help bridge the gap between genomics and applied wildlife conservation.

Genotyping-in-Thousands by sequencing (GT-seq) panel development and application to minimally invasive DNA samples to support studies in molecular ecology

Minimally invasive sampling (MIS) is widespread in wildlife studies; however, its utility for massively parallel DNA sequencing (MPS) is limited. Poor sample quality and contamination by exogenous DNA can make MIS challenging to use with modern genotyping-by-sequencing approaches, which have been traditionally developed for high-quality DNA sources. Given that MIS is often more appropriate in many contexts, there is a need to make such samples practical for harnessing MPS. Here, we test the ability for Genotyping-in-Thousands by sequencing (GT-seq), a multiplex amplicon sequencing approach, to effectively genotype minimally invasive cloacal DNA samples collected from the Western Rattlesnake (Crotalus oreganus), a threatened species in British Columbia, Canada. As there was no previous genetic information for this species, an optimized panel of 362 SNPs was selected for use with GT-seq from a de novo restriction site-associated DNA sequencing (RADseq) assembly. Comparisons of genotypes generated within and among RADseq and GT-seq for the same individuals found low rates of genotyping error (GT-seq: 0.50%; RADseq: 0.80%) and discordance (2.57%), the latter likely due to the different genotype calling models employed. GT-seq mean genotype discordance between blood and cloacal swab samples collected from the same individuals was also minimal (1.37%). Estimates of population diversity parameters were similar across GT-seq and RADseq data sets, as were inferred patterns of population structure. Overall, GT-seq can be effectively applied to low-quality DNA samples, minimizing the inefficiencies presented by exogenous DNA typically found in minimally invasive samples and continuing the expansion of molecular ecology and conservation genetics in the genomics era.

Genetic and genomic monitoring with minimally invasive sampling methods

The decreasing cost and increasing scope and power of emerging genomic technologies are reshaping the field of molecular ecology. However, many modern genomic approaches (e.g., RAD-seq) require large amounts of high-quality template DNA. This poses a problem for an active branch of conservation biology: genetic monitoring using minimally invasive sampling (MIS) methods. Without handling or even observing an animal, MIS methods (e.g., collection of hair, skin, faeces) can provide genetic information on individuals or populations. Such samples typically yield low-quality and/or quantities of DNA, restricting the type of molecular methods that can be used. Despite this limitation, genetic monitoring using MIS is an effective tool for estimating population demographic parameters and monitoring genetic diversity in natural populations. Genetic monitoring is likely to become more important in the future as many natural populations are undergoing anthropogenically driven declines, which are unlikely to abate without intensive adaptive management efforts that often include MIS approaches. Here, we profile the expanding suite of genomic methods and platforms compatible with producing genotypes from MIS, considering factors such as development costs and error rates. We evaluate how powerful new approaches will enhance our ability to investigate questions typically answered using genetic monitoring, such as estimating abundance, genetic structure and relatedness. As the field is in a period of unusually rapid transition, we also highlight the importance of legacy data sets and recommend how to address the challenges of moving between traditional and next-generation genetic monitoring platforms. Finally, we consider how genetic monitoring could move beyond genotypes in the future. For example, assessing microbiomes or epigenetic markers could provide a greater understanding of the relationship between individuals and their environment.

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

Targeted capture coupled with high-throughput sequencing can be used to gain information about nuclear sequence variation at hundreds to thousands of loci. Divergent reference capture makes use of molecular data of one species to enrich target loci in other (related) species. This is particularly valuable for nonmodel organisms, for which often no a priori knowledge exists regarding these loci. Here, we have used targeted capture to obtain data for 809 nuclear coding DNA sequences (CDS) in a nonmodel organism, the Eurasian lynx Lynx lynx, using baits designed with the help of the published genome of a related model organism (the domestic cat Felis catus). Using this approach, we were able to survey intraspecific variation at hundreds of nuclear loci in L. lynx across the species' European range. A large set of biallelic candidate SNPs was then evaluated using a high-throughput SNP genotyping platform (Fluidigm), which we then reduced to a final 96 SNP-panel based on assay performance and reliability; validation was carried out with 100 additional Eurasian lynx samples not included in the SNP discovery phase. The 96 SNP-panel developed from CDS performed very successfully in the identification of individuals and in population genetic structure inference (including the assignment of individuals to their source population). In keeping with recent studies, our results show that genic SNPs can be valuable for genetic monitoring of wildlife species.