Dataset Information

Phydms: software for phylogenetic analyses informed by deep mutational scanning.

ABSTRACT: It has recently become possible to experimentally measure the effects of all amino-acid point mutations to proteins using deep mutational scanning. These experimental measurements can inform site-specific phylogenetic substitution models of gene evolution in nature. Here we describe software that efficiently performs analyses with such substitution models. This software, phydms, can be used to compare the results of deep mutational scanning experiments to the selection on genes in nature. Given a phylogenetic tree topology inferred with another program, phydms enables rigorous comparison of how well different experiments on the same gene capture actual natural selection. It also enables re-scaling of deep mutational scanning data to account for differences in the stringency of selection in the lab and nature. Finally, phydms can identify sites that are evolving differently in nature than expected from experiments in the lab. As data from deep mutational scanning experiments become increasingly widespread, phydms will facilitate quantitative comparison of the experimental results to the actual selection pressures shaping evolution in nature.

SUBMITTER: Hilton SK

PROVIDER: S-EPMC5541924 | biostudies-literature | 2017

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

phydms: software for phylogenetic analyses informed by deep mutational scanning.

Hilton Sarah K SK Doud Michael B MB Bloom Jesse D JD

PeerJ 20170731

It has recently become possible to experimentally measure the effects of all amino-acid point mutations to proteins using deep mutational scanning. These experimental measurements can inform site-specific phylogenetic substitution models of gene evolution in nature. Here we describe software that efficiently performs analyses with such substitution models. This software, phydms, can be used to compare the results of deep mutational scanning experiments to the selection on genes in nature. Given ...[more]

PMID: 28785526

Similar Datasets

Project description:This dataset provides allele counts and raw fastqs for deep mutational scanning of the HIV-1 genes tat and rev when not-overlapped with one another (placed in the nef locus) as described in Fernandes et al. Functional segregation of overlapping genes in HIV Cell 2016 (in revision). Preselection (input) and post selection (replicate 1/2) files for every possible point mutant of these two HIV proteins from the NL4-3 background are given.Tab delimited files including codon counts across the amplicons are also included and are probably the most useful thing to most researchers. The data here was used to generate Figures 3 and 4 and 7 and might be of general use for people interested in deep mutational scanning, looking for signatures of epistasis in rev or tat, or reanalyzing and mining the data. FAQ: Why do the ends of each amplicon have such variation? In order to increase diversity across the flowcell, I pooled standard primers with N, NN, and NNN extensions to throw amplicons out of phase. When aligning you should trim the ends or ignore them. This means that the overlap between PE's can vary by 3 nt. Why are the filenames not easy to deal with? The filenames are tied to separate MiSeq runs. I hope to clean up the nomenclature and update this entry in the future while preserving the run information. You can get a sense of that as different residues will vary in Q-score, and that is mostly tied to the run they were pooled on and not any interesting biology. While this is makes it a little harder to follow, I think it's good to get a sense that doing this kind of analysis in high-throughput fashion leads to a reasonable amount of failure (i.e. RNA isolation, RT, fail) that led to repetition until we had good data for every position. Can you help me deal with this dataset? Yes. Please email me at jferna10@ucsc.edu, or contact me on twitter @jdf_ev. For reagents please contact Alan Frankel at frankel@cgl.ucsf.edu. Do you have the analysis scripts you used to process the data? Yes they are on github. https://github.com/nbstrauli/allele_frequency_trajectory_sim

Project description:Cytochrome P450s (CYPs) are a family of enzymes responsible for metabolizing nearly 80% of small molecule drugs. Variants in CYPs can substantially alter drug metabolism, which may result in improper dosing and severe adverse drug reactions. CYPs have low sequence conservation, making it difficult to anticipate whether variant effects measured in one CYP may extend to others based on sequence alone. Even closely related CYPs, like CYP2C9 and its closest homolog CYP2C19, have distinct phenotypic properties despite sharing 92% of their sequences. Thus, we used Variant Abundance by Massively Parallel sequencing (VAMP-seq) to measure the protein abundance of 7,660 missense variants in CYP2C19 expressed in cultured human cells. Our results confirmed positions and structural features critical for CYP function and revealed how variants at positions conserved across all eukaryotic CYPs influence abundance. We jointly analyzed 4,670 variants whose abundance was measured in both CYP2C19 and CYP2C9, finding that the homologs have different variant abundances in substrate recognition sites within the hydrophobic core, and that substitutions in some regions reduced abundance in CYP2C19 but not CYP2C9. We also measured the abundance of all single and some multiple WT amino acid exchanges between CYP2C19 and CYP2C9. While most exchanges had no effect, substitutions in substrate recognition site 4 (SRS4) reduced abundance in CYP2C19. When nearby amino acids were exchanged in double and triple mutants, we found distinct interactions between the sites in CYP2C19 and CYP2C9, revealing a region that is partially responsible for the difference in thermodynamic stability between the two homologs. Since these positions are also important for determining substrate specificity, there may be an evolutionary tradeoff between stability and altered enzymatic function. Finally, we used our data to analyze 368 previously unannotated human variants, finding that 43% had decreased abundance. Thus, by comparing variant effects between two closely related, important human genes, we have uncovered regions underlying their functional differences and paved the way for a more complete understanding of one of the most versatile families of enzymes.

Dataset Information

Phydms: software for phylogenetic analyses informed by deep mutational scanning.

Publications

phydms: software for phylogenetic analyses informed by deep mutational scanning.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets