Dataset Information

Is RAD-seq suitable for phylogenetic inference? An in silico assessment and optimization.

ABSTRACT: INFERRING PHYLOGENETIC RELATIONSHIPS BETWEEN CLOSELY RELATED TAXA CAN BE HINDERED BY THREE FACTORS: (1) the lack of informative molecular variation at short evolutionary timescale; (2) the lack of established markers in poorly studied taxa; and (3) the potential phylogenetic conflicts among different genomic regions due to incomplete lineage sorting or introgression. In this context, Restriction site Associated DNA sequencing (RAD-seq) seems promising as this technique can generate sequence data from numerous DNA fragments scattered throughout the genome, from a large number of samples, and without preliminary knowledge on the taxa under study. However, divergence beyond the within-species level will necessarily reduce the number of conserved and non-duplicated restriction sites, and therefore the number of loci usable for phylogenetic inference. Here, we assess the suitability of RAD-seq for phylogeny using a simulated experiment on the 12 Drosophila genomes, with divergence times ranging from 5 to 63 million years. These simulations show that RAD-seq allows the recovery of the known Drosophila phylogeny with strong statistical support, even for relatively ancient nodes. Notably, this conclusion is robust to the potentially confounding effects of sequencing errors, heterozygosity, and low coverage. We further show that clustering RAD-seq data using the BLASTN and SiLiX programs significantly improves the recovery of orthologous RAD loci compared with previously proposed approaches, especially for distantly related species. This study therefore validates the view that RAD sequencing is a powerful tool for phylogenetic inference.

SUBMITTER: Cariou M

PROVIDER: S-EPMC3631399 | biostudies-literature | 2013 Apr

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Is RAD-seq suitable for phylogenetic inference? An in silico assessment and optimization.

Cariou Marie M Duret Laurent L Charlat Sylvain S

Ecology and evolution 20130227 4

INFERRING PHYLOGENETIC RELATIONSHIPS BETWEEN CLOSELY RELATED TAXA CAN BE HINDERED BY THREE FACTORS: (1) the lack of informative molecular variation at short evolutionary timescale; (2) the lack of established markers in poorly studied taxa; and (3) the potential phylogenetic conflicts among different genomic regions due to incomplete lineage sorting or introgression. In this context, Restriction site Associated DNA sequencing (RAD-seq) seems promising as this technique can generate sequence data ...[more]

PMID: 23610629

Similar Datasets

Project description:Accurate genome-wide identification of orthologs is a central problem in comparative genomics, a fact reflected by the numerous orthology identification projects developed in recent years. However, only a few reports have compared their accuracy, and indeed, several recent efforts have not yet been systematically evaluated. Furthermore, orthology is typically only assessed in terms of function conservation, despite the phylogeny-based original definition of Fitch. We collected and mapped the results of nine leading orthology projects and methods (COG, KOG, Inparanoid, OrthoMCL, Ensembl Compara, Homologene, RoundUp, EggNOG, and OMA) and two standard methods (bidirectional best-hit and reciprocal smallest distance). We systematically compared their predictions with respect to both phylogeny and function, using six different tests. This required the mapping of millions of sequences, the handling of hundreds of millions of predicted pairs of orthologs, and the computation of tens of thousands of trees. In phylogenetic analysis or in functional analysis where high specificity is required, we find that OMA and Homologene perform best. At lower functional specificity but higher coverage level, OrthoMCL outperforms Ensembl Compara, and to a lesser extent Inparanoid. Lastly, the large coverage of the recent EggNOG can be of interest to build broad functional grouping, but the method is not specific enough for phylogenetic or detailed function analyses. In terms of general methodology, we observe that the more sophisticated tree reconstruction/reconciliation approach of Ensembl Compara was at times outperformed by pairwise comparison approaches, even in phylogenetic tests. Furthermore, we show that standard bidirectional best-hit often outperforms projects with more complex algorithms. First, the present study provides guidance for the broad community of orthology data users as to which database best suits their needs. Second, it introduces new methodology to verify orthology. And third, it sets performance standards for current and future approaches.

Project description:Drug-induced Torsade-de-Pointes (TdP) has been responsible for the withdrawal of many drugs from the market and is therefore of major concern to global regulatory agencies and the pharmaceutical industry. The Comprehensive in vitro Proarrhythmia Assay (CiPA) was proposed to improve prediction of TdP risk, using in silico models and in vitro multi-channel pharmacology data as integral parts of this initiative. Previously, we reported that combining dynamic interactions between drugs and the rapid delayed rectifier potassium current (IKr) with multi-channel pharmacology is important for TdP risk classification, and we modified the original O'Hara Rudy ventricular cell mathematical model to include a Markov model of IKr to represent dynamic drug-IKr interactions (IKr-dynamic ORd model). We also developed a novel metric that could separate drugs with different TdP liabilities at high concentrations based on total electronic charge carried by the major inward ionic currents during the action potential. In this study, we further optimized the IKr-dynamic ORd model by refining model parameters using published human cardiomyocyte experimental data under control and drug block conditions. Using this optimized model and manual patch clamp data, we developed an updated version of the metric that quantifies the net electronic charge carried by major inward and outward ionic currents during the steady state action potential, which could classify the level of drug-induced TdP risk across a wide range of concentrations and pacing rates. We also established a framework to quantitatively evaluate a system's robustness against the induction of early afterdepolarizations (EADs), and demonstrated that the new metric is correlated with the cell's robustness to the pro-EAD perturbation of IKr conductance reduction. In summary, in this work we present an optimized model that is more consistent with experimental data, an improved metric that can classify drugs at concentrations both near and higher than clinical exposure, and a physiological framework to check the relationship between a metric and EAD. These findings provide a solid foundation for using in silico models for the regulatory assessment of TdP risk under the CiPA paradigm.

Project description:BACKGROUND:Oseltamivir is a first-line antiviral drug, especially in primary hospitals. During the ongoing outbreak of coronavirus disease 2019 (COVID-19), most patients with COVID-19 who are symptomatic have used oseltamivir. Considering its popular and important role as an antiviral drug, it is necessary to evaluate oseltamivir in the treatment of COVID-19. OBJECTIVE:To evaluate the effect of oseltamivir against COVID-19. METHODS:Swiss-model was used to construct the structure of the N-terminal RNA-binding domain (NRBD) of the nucleoprotein (NC), papain-like protease (PLpro), and RNA-directed RNA polymerase (RdRp) of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). TM-align program was performed to compare the structure of the viral proteins with the structure of the neuraminidase of influenza A. Molecular docking was used to analyze the theoretical possibility of effective binding of oseltamivir with the active centers of the viral proteins. In vitro study was used to evaluate the antiviral efficiency of oseltamivir against SARS-CoV-2. By clinical case analysis, we statistically evaluated whether the history of oseltamivir use influenced the progression of the disease. RESULTS:The structures of NRBD, PLpro, and RdRp were built successfully. The results from TM-align suggested that the S protein, NRBD, 3C-like protease (3CLpro), PLPrO, and RdRp were structurally similar to the influenza A neuraminidase, with TM-scores of 0.30077, 0.19254, 0.28766, 0.30666, and 0.34047, respectively. Interestingly, the active center of 3CL pro was found to be similar to the active center from the neuraminidase of influenza A. Through an analysis of molecular docking, we discovered that oseltamivir carboxylic acid was more favorable to bind to the active site of 3CLpro effectively, but its inhibitory effect was not strong compared with the positive group. Finally, we used in vitro study and retrospective case analysis to verify our speculations. We found that oseltamivir is ineffective against SARS-CoV-2 in vitro study and the clinical use of oseltamivir did not improve the patients' symptoms and signs and did not slow the disease progression. CONCLUSIONS:We consider that oseltamivir isn't suitable for the treatment of COVID-19. During the outbreak of novel coronavirus, when oseltamivir is not effective for the patients after they take it, health workers should be highly vigilant about the possibility of COVID-19.

Dataset Information

Is RAD-seq suitable for phylogenetic inference? An in silico assessment and optimization.

Publications

Is RAD-seq suitable for phylogenetic inference? An in silico assessment and optimization.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets