Pacybara: accurate long-read sequencing for barcoded mutagenized allelic libraries.

MOTIVATION: Long-read sequencing technologies, an attractive solution for many applications, often suffer from higher error rates. Alignment of multiple reads can improve base-calling accuracy, but some applications, e.g. sequencing mutagenized libraries where multiple distinct clones differ by one or few variants, require the use of barcodes or unique molecular identifiers. Unfortunately, sequencing errors can interfere with correct barcode identification, and a given barcode sequence may be linked to multiple independent clones within a given library. RESULTS: Here we focus on the target application of sequencing mutagenized libraries in the context of multiplexed assays of variant effects (MAVEs). MAVEs are increasingly used to create comprehensive genotype-phenotype maps that can aid clinical variant interpretation. Many MAVE methods use long-read sequencing of barcoded mutant libraries for accurate association of barcode with genotype. Existing long-read sequencing pipelines do not account for inaccurate sequencing or nonunique barcodes. Here, we describe Pacybara, which handles these issues by clustering long reads based on the similarities of (error-prone) barcodes while also detecting barcodes that have been associated with multiple genotypes. Pacybara also detects recombinant (chimeric) clones and reduces false positive indel calls. In three example applications, we show that Pacybara identifies and correctly resolves these issues. AVAILABILITY AND IMPLEMENTATION: Pacybara, freely available at https://github.com/rothlab/pacybara, is implemented using R, Python, and bash for Linux. It runs on GNU/Linux HPC clusters via Slurm, PBS, or GridEngine schedulers. A single-machine simplex version is also available.

Pacybara: Accurate long-read sequencing for barcoded mutagenized allelic libraries.

Long read sequencing technologies, an attractive solution for many applications, often suffer from higher error rates. Alignment of multiple reads can improve base-calling accuracy, but some applications, e.g. sequencing mutagenized libraries where multiple distinct clones differ by one or few variants, require the use of barcodes or unique molecular identifiers. Unfortunately, sequencing errors can interfere with correct barcode identification, and a given barcode sequence may be linked to multiple independent clones within a given library. Here we focus on the target application of sequencing mutagenized libraries in the context of multiplexed assays of variant effects (MAVEs). MAVEs are increasingly used to create comprehensive genotype-phenotype maps that can aid clinical variant interpretation. Many MAVE methods use long-read sequencing of barcoded mutant libraries for accurate association of barcode with genotype. Existing long-read sequencing pipelines do not account for inaccurate sequencing or non-unique barcodes. Here, we describe Pacybara, which handles these issues by clustering long reads based on the similarities of (error-prone) barcodes while also detecting barcodes that have been associated with multiple genotypes. Pacybara also detects recombinant (chimeric) clones and reduces false positive indel calls. In three example applications, we show that Pacybara identifies and correctly resolves these issues.

PacRAT: a program to improve barcode-variant mapping from PacBio long reads using multiple sequence alignment.

SUMMARY: Use of PacBio sequencing for characterizing barcoded libraries of genetic variants is on the rise. However, current approaches in resolving PacBio sequencing artifacts can result in a high number of incorrectly identified or unusable reads. Here, we developed a PacBio Read Alignment Tool (PacRAT) that improves the accuracy of barcode-variant mapping through several steps of read alignment and consensus calling. To quantify the performance of our approach, we simulated PacBio reads from eight variant libraries of various lengths and showed that PacRAT improves the accuracy in pairing barcodes and variants across these libraries. Analysis of real (non-simulated) libraries also showed an increase in the number of reads that can be used for downstream analyses when using PacRAT. AVAILABILITY AND IMPLEMENTATION: PacRAT is written in Python and is freely available (https://github.com/dunhamlab/PacRAT). SUPPLEMENTARY INFORMATION: Supplemental data are available at Bioinformatics online.

Measuring Pharmacogene Variant Function at Scale Using Multiplexed Assays.

As costs of next-generation sequencing decrease, identification of genetic variants has far outpaced our ability to understand their functional consequences. This lack of understanding is a central challenge to a key promise of pharmacogenomics: using genetic information to guide drug selection and dosing. Recently developed multiplexed assays of variant effect enable experimental measurement of the function of thousands of variants simultaneously. Here, we describe multiplexed assays that have been performed on nearly 25,000 variants in eight key pharmacogenes (ADRB2, CYP2C9, CYP2C19, NUDT15, SLCO1B1, TMPT, VKORC1, and the LDLR promoter), discuss advances in experimental design, and explore key challenges that must be overcome to maximize the utility of multiplexed functional data.

Massively parallel characterization of CYP2C9 variant enzyme activity and abundance.

CYP2C9 encodes a cytochrome P450 enzyme responsible for metabolizing up to 15% of small molecule drugs, and CYP2C9 variants can alter the safety and efficacy of these therapeutics. In particular, the anti-coagulant warfarin is prescribed to over 15 million people annually and polymorphisms in CYP2C9 can affect individual drug response and lead to an increased risk of hemorrhage. We developed click-seq, a pooled yeast-based activity assay, to test thousands of variants. Using click-seq, we measured the activity of 6,142 missense variants in yeast. We also measured the steady-state cellular abundance of 6,370 missense variants in a human cell line by using variant abundance by massively parallel sequencing (VAMP-seq). These data revealed that almost two-thirds of CYP2C9 variants showed decreased activity and that protein abundance accounted for half of the variation in CYP2C9 function. We also measured activity scores for 319 previously unannotated human variants, many of which may have clinical relevance.

Expression and Functional Characterization of Breast Cancer-Associated Cytochrome P450 4Z1 in Saccharomyces cerevisiae.

CYP4Z1 is an "orphan" cytochrome P450 (P450) enzyme that has provoked interest because of its hypothesized role in breast cancer through formation of the signaling molecule 20-hydroxyeicosatetraenoic acid (20-HETE). We expressed human CYP4Z1 in Saccharomyces cerevisiae and evaluated its catalytic capabilities toward arachidonic and lauric acids (AA and LA). Specific and sensitive mass spectrometry assays enabled discrimination of the regioselectivity of hydroxylation of these two fatty acids. CYP4Z1 generated 7-, 8-, 9-, 10-, and 11-hydroxy LA, whereas the 12-hydroxy metabolite was not detected. HET0016, the prototypic CYP4 inhibitor, only weakly inhibited laurate metabolite formation (IC50 ∼15 µM). CYP4Z1 preferentially oxidized AA to the 14(S),15(R)-epoxide with high regioselectivity and stereoselectivity, a reaction that was also insensitive to HET0016, but neither 20-HETE nor 20-carboxy-AA were detectable metabolites. Docking of LA and AA into a CYP4Z1 homology model was consistent with this preference for internal fatty acid oxidation. Thus, human CYP4Z1 has an inhibitor profile and product regioselectivity distinct from most other CYP4 enzymes, consistent with CYP4Z1's lack of a covalently linked heme. These data suggest that, if CYP4Z1 modulates breast cancer progression, it does so by a mechanism other than direct production of 20-HETE.

Experimental Evolution Reveals Favored Adaptive Routes to Cell Aggregation in Yeast.

Yeast flocculation is a community-building cell aggregation trait that is an important mechanism of stress resistance and a useful phenotype for brewers; however, it is also a nuisance in many industrial processes, in clinical settings, and in the laboratory. Chemostat-based evolution experiments are impaired by inadvertent selection for aggregation, which we observe in 35% of populations. These populations provide a testing ground for understanding the breadth of genetic mechanisms Saccharomyces cerevisiae uses to flocculate, and which of those mechanisms provide the biggest adaptive advantages. In this study, we employed experimental evolution as a tool to ask whether one or many routes to flocculation are favored, and to engineer a strain with reduced flocculation potential. Using a combination of whole genome sequencing and bulk segregant analysis, we identified causal mutations in 23 independent clones that had evolved cell aggregation during hundreds of generations of chemostat growth. In 12 of those clones, we identified a transposable element insertion in the promoter region of known flocculation gene FLO1, and, in an additional five clones, we recovered loss-of-function mutations in transcriptional repressor TUP1, which regulates FLO1 and other related genes. Other causal mutations were found in genes that have not been previously connected to flocculation. Evolving a flo1 deletion strain revealed that this single deletion reduces flocculation occurrences to 3%, and demonstrated the efficacy of using experimental evolution as a tool to identify and eliminate the primary adaptive routes for undesirable traits.