RADIS: analysis of RAD-seq data for interspecific phylogeny.
Ontology highlight
ABSTRACT: Unlabelled
In an attempt to make the processing of RAD-seq data easier and allow rapid and automated exploration of parameters/data for phylogenetic inference, we introduce the perl pipeline RADIS Users of RADIS can let their raw Illumina data be processed up to phylogenetic tree inference, or stop (and restart) the process at some point. Different values for key parameters can be explored in a single analysis (e.g. loci building, sample/loci selection), making possible a thorough exploration of data. RADIS relies on Stacks for demultiplexing of data, removing PCR duplicates and building individual and catalog loci. Scripts have been specifically written for trimming of reads and loci/sample selection. Finally, RAxML is used for phylogenetic inferences, though other software may be utilized.Availability and implementation
RADIS is written in perl, designed to run on Linux and Unix platforms. RADIS and its manual are freely available from http://www1.montpellier.inra.fr/CBGP/software/RADIS/.Contact
astrid.cruaud@supagro.inra.frSupplementary information
Supplementary data are available at Bioinformatics online.
SUBMITTER: Cruaud A
PROVIDER: S-EPMC5039923 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:Restriction site-associated DNA sequencing (RAD-seq) has become a powerful and widely used tool in molecular ecology studies as it allows to cost-effectively recover thousands of polymorphic sites across individuals of non-model organisms. However, its successful implementation in population genetics relies on correct data processing that would minimize potential loci-assembly biases and consequent genotyping error rates. RAD-seq data processing when no reference genome is available involves the assembly of hundreds of thousands high-throughput sequencing reads into orthologous loci, for which various key parameter values need to be selected by the researcher. Previous studies exploring the effect of these parameter values found or assumed that a larger number of recovered polymorphic loci is associated with a better assembly. Here, using three RAD-seq datasets from different species, we explore the effect of read filtering, loci assembly and polymorphic site selection on number of markers obtained and genetic differentiation inferred using the Stacks software. We find (i) that recovery of higher numbers of polymorphic loci is not necessarily associated with higher genetic differentiation, (ii) that the presence of PCR duplicates, selected loci assembly parameters and selected SNP filtering parameters affect the number of recovered polymorphic loci and degree of genetic differentiation, and (iii) that this effect is different in each dataset, meaning that defining a systematic universal protocol for RAD-seq data analysis may lead to missing relevant information about population differentiation.
| S-EPMC6549478 | biostudies-literature
Project description:Restriction-site associated DNA sequencing (RAD-seq) can identify and score thousands of genetic markers from a group of samples for population-genetics studies. One challenge of de novo RAD-seq analysis is to distinguish paralogous sequence variants (PSVs) from true single-nucleotide polymorphisms (SNPs) associated with orthologous loci. In the absence of a reference genome, it is difficult to differentiate true SNPs from PSVs, and their impact on downstream analysis remains unclear. Here, we introduce a network-based approach, PMERGE that connects fragments based on their DNA sequence similarity to identify probable PSVs. Applying our method to de novo RAD-seq data from 150 Atlantic salmon (Salmo salar) samples collected from 15 locations across the Southern Newfoundland coast allowed the identification of 87% of total PSVs identified through alignment to the Atlantic salmon genome. Removal of these paralogs altered the inferred population structure, highlighting the potential impact of filtering in RAD-seq analysis. PMERGE is also applied to a green crab (Carcinus maenas) data set consisting of 242 samples from 11 different locations and was successfully able to identify and remove the majority of paralogous loci (62%). The PMERGE software can be run as part of the widely used Stacks analysis package.
| S-EPMC6065343 | biostudies-literature
Project description:Hawaiian Melicope are one of the major adaptive radiations of the Hawaiian Islands comprising 54 endemic species. The lineage is monophyletic with an estimated crown age predating the rise of the current high islands. Phylogenetic inference based on Sanger sequencing has not been sufficient to resolve species or deeper level relationships. Here, we apply restriction site-associated DNA sequencing (RAD-seq) to the lineage to infer phylogenetic relationships. We employ Quartet Sampling to assess information content and statistical support, and to quantify discordance as well as partitioned ABBA-BABA tests to uncover evidence of introgression. Our new results drastically improved resolution of relationships within Hawaiian Melicope. The lineage is divided into five fully supported main clades, two of which correspond to morphologically circumscribed infrageneric groups. We provide evidence for both ancestral and current hybridization events. We confirm the necessity for a taxonomic revision of the Melicope section Pelea, as well as a re-evaluation of several species complexes by combining genomic and morphological data.
| S-EPMC6758601 | biostudies-literature
Project description:Previous phylogenetic studies in oaks (Quercus, Fagaceae) have failed to resolve the backbone topology of the genus with strong support. Here, we utilize next-generation sequencing of restriction-site associated DNA (RAD-Seq) to resolve a framework phylogeny of a predominantly American clade of oaks whose crown age is estimated at 23-33 million years old. Using a recently developed analytical pipeline for RAD-Seq phylogenetics, we created a concatenated matrix of 1.40 E06 aligned nucleotides, constituting 27,727 sequence clusters. RAD-Seq data were readily combined across runs, with no difference in phylogenetic placement between technical replicates, which overlapped by only 43-64% in locus coverage. 17% (4,715) of the loci we analyzed could be mapped with high confidence to one or more expressed sequence tags in NCBI Genbank. A concatenated matrix of the loci that BLAST to at least one EST sequence provides approximately half as many variable or parsimony-informative characters as equal-sized datasets from the non-EST loci. The EST-associated matrix is more complete (fewer missing loci) and has slightly lower homoplasy than non-EST subsampled matrices of the same size, but there is no difference in phylogenetic support or relative attribution of base substitutions to internal versus terminal branches of the phylogeny. We introduce a partitioned RAD visualization method (implemented in the R package RADami; http://cran.r-project.org/web/packages/RADami) to investigate the possibility that suboptimal topologies supported by large numbers of loci--due, for example, to reticulate evolution or lineage sorting--are masked by the globally optimal tree. We find no evidence for strongly-supported alternative topologies in our study, suggesting that the phylogeny we recover is a robust estimate of large-scale phylogenetic patterns in the American oak clade. Our study is one of the first to demonstrate the utility of RAD-Seq data for inferring phylogeny in a 23-33 million year-old clade.
| S-EPMC3976371 | biostudies-literature
Project description:BACKGROUND:New single-cell isolation technologies are facilitating studies on the transcriptomics of individual cells. Bio-Rad ddSEQ is a droplet-based microfluidic system that, when coupled with downstream Illumina library preparation and sequencing, enables the monitoring of thousands of genes per cell. Sequenced reads show unique features that do not permit the use of freely available tools to perform single cell demultiplexing. RESULTS:We present ddSeeker, a tool to perform initial processing and quality metrics of reads generated through Bio-Rad ddSEQ/Illumina experiments. Its application to the Illumina test dataset demonstrates that ddSeeker performs better than Illumina BaseSpace software, enabling a higher recovery of valid reads. We also show its utility in the analysis of an in-house dataset including two read sets characterized by low and high sequencing quality. ddSeeker and its source code are available at https://github.com/cgplab/ddSeeker . CONCLUSIONS:ddSeeker is a freely available tool to perform initial processing and quality metrics of reads generated through Bio-Rad ddSEQ/Illumina single cell transcriptomic experiments.
| S-EPMC6304778 | biostudies-other
Project description:BACKGROUND:In genotyping-by-sequencing (GBS) and restriction site-associated DNA sequencing (RAD-seq), read depth is important for assessing the quality of genotype calls and estimating allele dosage in polyploids. However, existing pipelines for GBS and RAD-seq do not provide read counts in formats that are both accurate and easy to access. Additionally, although existing pipelines allow previously-mined SNPs to be genotyped on new samples, they do not allow the user to manually specify a subset of loci to examine. Pipelines that do not use a reference genome assign arbitrary names to SNPs, making meta-analysis across projects difficult. RESULTS:We created the software TagDigger, which includes three programs for analyzing GBS and RAD-seq data. The first script, tagdigger_interactive.py, rapidly extracts read counts and genotypes from FASTQ files using user-supplied sets of barcodes and tags. Input and output is in CSV format so that it can be opened by spreadsheet software. Tag sequences can also be imported from the Stacks, TASSEL-GBSv2, TASSEL-UNEAK, or pyRAD pipelines, and a separate file can be imported listing the names of markers to retain. A second script, tag_manager.py, consolidates marker names and sequences across multiple projects. A third script, barcode_splitter.py, assists with preparing FASTQ data for deposit in a public archive by splitting FASTQ files by barcode and generating MD5 checksums for the resulting files. CONCLUSIONS:TagDigger is open-source and freely available software written in Python 3. It uses a scalable, rapid search algorithm that can process over 100 million FASTQ reads per hour. TagDigger will run on a laptop with any operating system, does not consume hard drive space with intermediate files, and does not require programming skill to use.
| S-EPMC4940913 | biostudies-other
Project description:Restriction-site associated DNA sequencing (RAD-seq) has recently become an important method to generate genome-wide molecular data for species delimitation, phylogeography, and population genetic studies. However, very few empirical studies have so far tested its applicability in phylogenetic reconstruction. The alpine-arctic genus Diapensia was selected to study the origin of the disjunction between the Arctic and the Himalayan-Hengduan Mountains (HHM). However, a previous phylogenetic analysis based on one nuclear and four plastid DNA regions failed to resolve the oldest divergences in Diapensia as well as the relationship between the two HHM species. Here we reconstruct a fully resolved phylogeny of Diapensia and address the conflict between the currently accepted taxonomy and the gene trees in the HHM species using RAD-seq. Based on a data set containing 2,650 loci selected to maximize the number of parsimony informative sites and allowing for a high level of missing data (51%), the phylogeny of Diapensia was fully resolved and each of the four species was reciprocally monophyletic. Whereas the arctic D. lapponica was inferred as sister to the HHM clade in the previous study, the RAD-seq data resolved the two arctic species as sisters to the HHM clade. Similar relationships were inferred from a differently filtered data set with far fewer loci (114) and less missing data (21%), but with lower support and with one of the two HHM species as non-monophyletic. Bayesian concordance analysis and Patterson's D-statistic tests suggested that admixture has occurred between the two HHM species.
| S-EPMC4598014 | biostudies-literature
Project description:Prunus is an economically important genus widely distributed in the temperate Northern Hemisphere. Previous studies on the genus using a variety of loci yielded conflicting phylogenetic hypotheses. Here, we generated nuclear reduced representation sequencing data and plastid genomes for 36 Prunus individuals and two outgroups. Both nuclear and plastome data recovered a well-resolved phylogeny. The species were divided into three main clades corresponding to their inflorescence types, - the racemose group, the solitary-flower group and the corymbose group - with the latter two sister to one another. Prunus was inferred to have diversified initially in the Late Cretaceous around 67.32 million years ago. The diversification of the three major clades began between the Paleocene and Miocene, suggesting that paleoclimatic events were an important driving force for Prunus diversification. Ancestral state reconstructions revealed that the most recent common ancestor of Prunus had racemose inflorescences, and the solitary-flower and corymb inflorescence types were derived by reduction of flower number and suppression of the rachis, respectively. We also tested the hybrid origin hypothesis of the racemose group proposed in previous studies. Prunus has undergone extensive hybridization events, although it is difficult to identify conclusively specific instances of hybridization when using SNP data, especially deep in the phylogeny. Our study provides well-resolved nuclear and plastid phylogenies of Prunus, reveals substantial cytonuclear discord at shallow scales, and sheds new light on inflorescence evolution in this economically important lineage.
| S-EPMC10435964 | biostudies-literature
Project description:Reduced representation libraries are being used as a preferred source of markers to address population genetic questions. However, libraries of RAD-Seq variants often suffer from significant percentage of missing data. In addition, algorithms used to mine SNPs from the raw data may also underscore biological variation. We investigate the effect of biological diversity in mining SNPs from the program STACKS and the effect of missing data on individual assignment implemented in STRUCTURE. We observed that changing diversity parameters in STACKS significantly alters the number of SNPs discovered and allowing for higher percentage of missing data retrieves more loci and possibly more power for individual assignment.
| S-EPMC4256836 | biostudies-literature
Project description:Hydrozoans are known for their complex life cycles, which can alternate between an asexually reproducing polyp stage and a sexually reproducing medusa stage. Most hydrozoan species, however, lack a free-living medusa stage and instead display a developmentally truncated form, called a medusoid or sporosac, which generally remains attached to the polyp. Although evolutionary transitions in medusa truncation and loss have been investigated phylogenetically, little is known about the genes involved in the development and loss of this life cycle stage. Here, we present a new workflow for evaluating differential expression (DE) between two species using short read Illumina RNA-seq data. Through interspecific DE analyses between two hydractiniid hydrozoans, Hydractinia symbiolongicarpus and Podocoryna carnea, we identified genes potentially involved in the developmental, functional, and morphological differences between the fully developed medusa of P. carnea and reduced sporosac of H. symbiolongicarpus. A total of 10,909 putative orthologs of H. symbiolongicarpus and P. carnea were identified from de novo assemblies of short read Illumina data. DE analysis revealed 938 of these are differentially expressed between P. carnea developing and adult medusa, when compared with H. symbiolongicarpus sporosacs, the majority of which have not been previously characterized in cnidarians. In addition, several genes with no corresponding ortholog in H. symbiolongicarpus were expressed in developing medusa of P. carnea. Results presented here show interspecific DE analyses of RNA-seq data to be a sensitive and reliable method for identifying genes and gene pathways potentially involved in morphological and life cycle differences between species.
| S-EPMC4558869 | biostudies-literature