Dataset Information

Experimental validation of in silico predicted RAD locus frequencies using genomic resources and short read data from a model marine mammal.

ABSTRACT: BACKGROUND:Restriction site-associated DNA sequencing (RADseq) has revolutionized the study of wild organisms by allowing cost-effective genotyping of thousands of loci. However, for species lacking reference genomes, it can be challenging to select the restriction enzyme that offers the best balance between the number of obtained RAD loci and depth of coverage, which is crucial for a successful outcome. To address this issue, PredRAD was recently developed, which uses probabilistic models to predict restriction site frequencies from a transcriptome assembly or other sequence resource based on either GC content or mono-, di- or trinucleotide composition. This program generates predictions that are broadly consistent with estimates of the true number of restriction sites obtained through in silico digestion of available reference genome assemblies. However, in practice the actual number of loci obtained could potentially differ as incomplete enzymatic digestion or patchy sequence coverage across the genome might lead to some loci not being represented in a RAD dataset, while erroneous assembly could potentially inflate the number of loci. To investigate this, we used genome and transcriptome assemblies together with RADseq data from the Antarctic fur seal (Arctocephalus gazella) to compare PredRAD predictions with empirical estimates of the number of loci obtained via in silico digestion and from de novo assemblies. RESULTS:PredRAD yielded consistently higher predicted numbers of restriction sites for the transcriptome assembly relative to the genome assembly. The trinucleotide and dinucleotide models also predicted higher frequencies than the mononucleotide or GC content models. Overall, the dinucleotide and trinucleotide models applied to the transcriptome and the genome assemblies respectively generated predictions that were closest to the number of restriction sites estimated by in silico digestion. Furthermore, the number of de novo assembled RAD loci mapping to restriction sites was similar to the expectation based on in silico digestion. CONCLUSIONS:Our study reveals generally high concordance between PredRAD predictions and empirical estimates of the number of RAD loci. This further supports the utility of PredRAD, while also suggesting that it may be feasible to sequence and assemble the majority of RAD loci present in an organism's genome.

SUBMITTER: Vendrami DLJ

PROVIDER: S-EPMC6341687 | biostudies-literature | 2019 Jan

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Experimental validation of in silico predicted RAD locus frequencies using genomic resources and short read data from a model marine mammal.

Vendrami David L J DLJ Forcada Jaume J Hoffman Joseph I JI

BMC genomics 20190122 1

<h4>Background</h4>Restriction site-associated DNA sequencing (RADseq) has revolutionized the study of wild organisms by allowing cost-effective genotyping of thousands of loci. However, for species lacking reference genomes, it can be challenging to select the restriction enzyme that offers the best balance between the number of obtained RAD loci and depth of coverage, which is crucial for a successful outcome. To address this issue, PredRAD was recently developed, which uses probabilistic mode ...[more]

PMID: 30669975

Similar Datasets

Project description:Multi-locus sequence typing (MLST) has become the gold standard for population analyses of bacterial pathogens. This method focuses on the sequences of a small number of loci (usually seven) to divide the population and is simple, robust and facilitates comparison of results between laboratories and over time. Over the last decade, researchers and population health specialists have invested substantial effort in building up public MLST databases for nearly 100 different bacterial species, and these databases contain a wealth of important information linked to MLST sequence types such as time and place of isolation, host or niche, serotype and even clinical or drug resistance profiles. Recent advances in sequencing technology mean it is increasingly feasible to perform bacterial population analysis at the whole genome level. This offers massive gains in resolving power and genetic profiling compared to MLST, and will eventually replace MLST for bacterial typing and population analysis. However given the wealth of data currently available in MLST databases, it is crucial to maintain backwards compatibility with MLST schemes so that new genome analyses can be understood in their proper historical context.We present a software tool, SRST, for quick and accurate retrieval of sequence types from short read sets, using inputs easily downloaded from public databases. SRST uses read mapping and an allele assignment score incorporating sequence coverage and variability, to determine the most likely allele at each MLST locus. Analysis of over 3,500 loci in more than 500 publicly accessible Illumina read sets showed SRST to be highly accurate at allele assignment. SRST output is compatible with common analysis tools such as eBURST, Clonal Frame or PhyloViz, allowing easy comparison between novel genome data and MLST data. Alignment, fastq and pileup files can also be generated for novel alleles.SRST is a novel software tool for accurate assignment of sequence types using short read data. Several uses for the tool are demonstrated, including quality control for high-throughput sequencing projects, plasmid MLST and analysis of genomic data during outbreak investigation. SRST is open-source, requires Python, BWA and SamTools, and is available from http://srst.sourceforge.net.

Project description:Marine mammals have greatly benefitted from a shift from resource exploitation towards conservation. Often lauded as symbols of conservation success, some marine mammal populations have shown remarkable recoveries after severe depletions. Others have remained at low abundance levels, continued to decline, or become extinct or extirpated. Here we provide a quantitative assessment of (1) publicly available population-level abundance data for marine mammals worldwide, (2) abundance trends and recovery status, and (3) historic population decline and recent recovery. We compiled 182 population abundance time series for 47 species and identified major data gaps. In order to compare across the largest possible set of time series with varying data quality, quantity and frequency, we considered an increase in population abundance as evidence of recovery. Using robust log-linear regression over three generations, we were able to classify abundance trends for 92 spatially non-overlapping populations as Significantly Increasing (42%), Significantly Decreasing (10%), Non-Significant Change (28%) and Unknown (20%). Our results were comparable to IUCN classifications for equivalent species. Among different groupings, pinnipeds and other marine mammals (sirenians, polar bears and otters) showed the highest proportion of recovering populations, likely benefiting from relatively fast life histories and nearshore habitats that provided visibility and protective management measures. Recovery was less frequent among cetaceans, but more common in coastal than offshore populations. For marine mammals with available historical abundance estimates (n = 47), larger historical population declines were associated with low or variable recent recoveries so far. Overall, our results show that many formerly depleted marine mammal populations are recovering. However, data-deficient populations and those with decreasing and non-significant trends require attention. In particular, increased study of populations with major data gaps, including offshore small cetaceans, cryptic species, and marine mammals in low latitudes and developing nations, is needed to better understand the status of marine mammal populations worldwide.

Dataset Information

Experimental validation of in silico predicted RAD locus frequencies using genomic resources and short read data from a model marine mammal.

Publications

Experimental validation of in silico predicted RAD locus frequencies using genomic resources and short read data from a model marine mammal.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets