A unifying model that explains the origins of human inverted copy number variants.

With the release of the telomere-to-telomere human genome sequence and the availability of both long-read sequencing and optical genome mapping techniques, the identification of copy number variants (CNVs) and other structural variants is providing new insights into human genetic disease. Different mechanisms have been proposed to account for the novel junctions in these complex architectures, including aberrant forms of DNA replication, non-allelic homologous recombination, and various pathways that repair DNA breaks. Here, we have focused on a set of structural variants that include an inverted segment and propose that they share a common initiating event: an inverted triplication with long, unstable palindromic junctions. The secondary rearrangement of these palindromes gives rise to the various forms of inverted structural variants. We postulate that this same mechanism (ODIRA: origin-dependent inverted-repeat amplification) that creates the inverted CNVs in inherited syndromes also generates the palindromes found in cancers.

Species-wide survey of the expressivity and complexity spectrum of traits in yeast.

Assessing the complexity and expressivity of traits at the species level is an essential first step to better dissect the genotype-phenotype relationship. As trait complexity behaves dynamically, the classic dichotomy between monogenic and complex traits is too simplistic. However, no systematic assessment of this complexity spectrum has been carried out on a population scale to date. In this context, we generated a large diallel hybrid panel composed of 190 unique hybrids coming from 20 natural isolates representative of the S. cerevisiae genetic diversity. For each of these hybrids, a large progeny of 160 individuals was obtained, leading to a total of 30,400 offspring individuals. Their mitotic growth was evaluated on 38 conditions inducing various cellular stresses. We developed a classification algorithm to analyze the phenotypic distributions of offspring and assess the trait complexity. We clearly found that traits are mainly complex at the population level. On average, we found that 91.2% of cross/trait combinations exhibit high complexity, while monogenic and oligogenic cases accounted for only 4.1% and 4.7%, respectively. However, the complexity spectrum is very dynamic, trait specific and tightly related to genetic backgrounds. Overall, our study provided greater insight into trait complexity as well as the underlying genetic basis of its spectrum in a natural population.

Template switching between the leading and lagging strands at replication forks generates inverted copy number variants through hairpin-capped extrachromosomal DNA.

Inherited and germ-line de novo copy number variants (CNVs) are increasingly found to be correlated with human developmental and cancerous phenotypes. Several models for template switching during replication have been proposed to explain the generation of these gross chromosomal rearrangements. We proposed a model of template switching (ODIRA-origin dependent inverted repeat amplification) in which simultaneous ligation of the leading and lagging strands at diverging replication forks could generate segmental inverted triplications through an extrachromosomal inverted circular intermediate. Here, we created a genetic assay using split-ura3 cassettes to trap the proposed inverted intermediate. However, instead of recovering circular inverted intermediates, we found inverted linear chromosomal fragments ending in native telomeres-suggesting that a template switch had occurred at the centromere-proximal fork of a replication bubble. As telomeric inverted hairpin fragments can also be created through double strand breaks we tested whether replication errors or repair of double stranded DNA breaks were the most likely initiating event. The results from CRISPR/Cas9 cleavage experiments and growth in the replication inhibitor hydroxyurea indicate that it is a replication error, not a double stranded break that creates the inverted junctions. Since inverted amplicons of the SUL1 gene occur during long-term growth in sulfate-limited chemostats, we sequenced evolved populations to look for evidence of linear intermediates formed by an error in replication. All of the data are compatible with a two-step version of the ODIRA model in which sequential template switching at short inverted repeats between the leading and lagging strands at a replication fork, followed by integration via homologous recombination, generates inverted interstitial triplications.

Functional interpretation, cataloging, and analysis of 1,341 glucose-6-phosphate dehydrogenase variants.

Glucose-6-phosphate dehydrogenase (G6PD) deficiency affects over 500 million individuals who can experience anemia in response to oxidative stressors such as certain foods and drugs. Recently, the World Health Organization (WHO) called for revisiting G6PD variant classification as a priority to implement genetic medicine in low- and middle-income countries. Toward this goal, we sought to collect reports of G6PD variants and provide interpretations. We identified 1,341 G6PD variants in population and clinical databases. Using the ACMG standards and guidelines for the interpretation of sequence variants, we provided interpretations for 268 variants, including 186 variants that were not reported or of uncertain significance in ClinVar, bringing the total number of variants with non-conflicting interpretations to 400. For 414 variants with functional or clinical data, we analyzed associations between activity, stability, and current classification systems, including the new 2022 WHO classification. We corroborated known challenges with classification systems, including phenotypic variation, emphasizing the importance of comparing variant effects across individuals and studies. Biobank data made available by All of Us illustrate the benefit of large-scale sequencing and phenotyping by adding additional support connecting variants to G6PD-deficient anemia. By leveraging available data and interpretation guidelines, we created a repository for information on G6PD variants and nearly doubled the number of variants with clinical interpretations. These tools enable better interpretation of G6PD variants for the implementation of genetic medicine.

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.

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.

Identification of aneuploidy-tolerating mutations.

Aneuploidy causes a proliferative disadvantage in all normal cells analyzed to date, yet this condition is associated with a disease characterized by unabated proliferative potential, cancer. The mechanisms that allow cancer cells to tolerate the adverse effects of aneuploidy are not known. To probe this question, we identified aneuploid yeast strains with improved proliferative abilities. Their molecular characterization revealed strain-specific genetic alterations as well as mutations shared between different aneuploid strains. Among the latter, a loss-of-function mutation in the gene encoding the deubiquitinating enzyme Ubp6 improves growth rates in four different aneuploid yeast strains by attenuating the changes in intracellular protein composition caused by aneuploidy. Our results demonstrate the existence of aneuploidy-tolerating mutations that improve the fitness of multiple different aneuploidies and highlight the importance of ubiquitin-proteasomal degradation in suppressing the adverse effects of aneuploidy.

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.

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.

The environmental stress response causes ribosome loss in aneuploid yeast cells.

Aneuploidy, a condition characterized by whole chromosome gains and losses, is often associated with significant cellular stress and decreased fitness. However, how cells respond to the aneuploid state has remained controversial. In aneuploid budding yeast, two opposing gene-expression patterns have been reported: the "environmental stress response" (ESR) and the "common aneuploidy gene-expression" (CAGE) signature, in which many ESR genes are oppositely regulated. Here, we investigate this controversy. We show that the CAGE signature is not an aneuploidy-specific gene-expression signature but the result of normalizing the gene-expression profile of actively proliferating aneuploid cells to that of euploid cells grown into stationary phase. Because growth into stationary phase is among the strongest inducers of the ESR, the ESR in aneuploid cells was masked when stationary phase euploid cells were used for normalization in transcriptomic studies. When exponentially growing euploid cells are used in gene-expression comparisons with aneuploid cells, the CAGE signature is no longer evident in aneuploid cells. Instead, aneuploid cells exhibit the ESR. We further show that the ESR causes selective ribosome loss in aneuploid cells, providing an explanation for the decreased cellular density of aneuploid cells. We conclude that aneuploid budding yeast cells mount the ESR, rather than the CAGE signature, in response to aneuploidy-induced cellular stresses, resulting in selective ribosome loss. We propose that the ESR serves two purposes in aneuploid cells: protecting cells from aneuploidy-induced cellular stresses and preventing excessive cellular enlargement during slowed cell cycles by down-regulating translation capacity.

Fitness benefits of loss of heterozygosity in Saccharomyces hybrids.

With two genomes in the same organism, interspecific hybrids have unique fitness opportunities and costs. In both plants and yeasts, wild, pathogenic, and domesticated hybrids may eliminate portions of one parental genome, a phenomenon known as loss of heterozygosity (LOH). Laboratory evolution of hybrid yeast recapitulates these results, with LOH occurring in just a few hundred generations of propagation. In this study, we systematically looked for alleles that are beneficial when lost in order to determine how prevalent this mode of adaptation may be and to determine candidate loci that might underlie the benefits of larger-scale chromosome rearrangements. These aims were accomplished by mating Saccharomyces uvarum with the S. cerevisiae deletion collection to create hybrids such that each nonessential S. cerevisiae allele is deleted. Competitive fitness assays of these pooled, barcoded, hemizygous strains, and accompanying controls, revealed a large number of loci for which LOH is beneficial. We found that the fitness effects of hemizygosity are dependent on the species context, the selective environment, and the species origin of the deleted allele. Further, we found that hybrids have a wider distribution of fitness consequences versus matched S. cerevisiae hemizygous diploids. Our results suggest that LOH can be a successful strategy for adaptation of hybrids to new environments, and we identify candidate loci that drive the chromosomal rearrangements observed in evolution of yeast hybrids.

Transposon insertional mutagenesis in Saccharomyces uvarum reveals trans-acting effects influencing species-dependent essential genes.

To understand how complex genetic networks perform and regulate diverse cellular processes, the function of each individual component must be defined. Comprehensive phenotypic studies of mutant alleles have been successful in model organisms in determining what processes depend on the normal function of a gene. These results are often ported to newly sequenced genomes by using sequence homology. However, sequence similarity does not always mean identical function or phenotype, suggesting that new methods are required to functionally annotate newly sequenced species. We have implemented comparative analysis by high-throughput experimental testing of gene dispensability in Saccharomyces uvarum, a sister species of Saccharomyces cerevisiae. We created haploid and heterozygous diploid Tn7 insertional mutagenesis libraries in S. uvarum to identify species-dependent essential genes, with the goal of detecting genes with divergent functions and/or different genetic interactions. Comprehensive gene dispensability comparisons with S. cerevisiae predicted diverged dispensability at 12% of conserved orthologs, and validation experiments confirmed 22 differentially essential genes. Despite their differences in essentiality, these genes were capable of cross-species complementation, demonstrating that trans-acting factors that are background-dependent contribute to differential gene essentiality. This study shows that direct experimental testing of gene disruption phenotypes across species can inform comparative genomic analyses and improve gene annotations. Our method can be widely applied in microorganisms to further our understanding of genome evolution.

Directed evolution as an approach to increase fructose utilization in synthetic grape juice by wine yeast AWRI 796.

In winemaking, slow or stuck alcoholic fermentation can impact processing efficiency and wine quality. Residual fructose in the later stages of fermentation can leave the wine 'out of specification' unless removed, which requires reinoculation or use of a more fructophilic yeast. As such, robust, fermentation efficient strains are still highly desirable to reduce this risk. We report on a combined EMS mutagenesis and Directed Evolution (DE) approach as a 'proof of concept' to improve fructose utilization and decrease fermentation duration. One evolved isolate, Tee 9, was evaluated against the parent, AWRI 796 in defined medium (CDGJM) and Semillon juice. Interestingly, Tee 9 exhibited improved fermentation in CDGJM at several nitrogen contents, but not in juice. Genomic comparison between AWRI 796 and Tee 9 identified 371 mutations, but no chromosomal copy number variation. A total of 95 noncoding and 276 coding mutations were identified in 297 genes (180 of which encode proteins with one or more substitutions). Whilst introduction of two of these, Gid7 (E726K) or Fba1 (G135S), into AWRI 796 did not lead to the fermentation improvement seen in Tee 9, similar allelic swaps with the other mutations are needed to understand Tee 9's adaption to CDGJM. Furthermore, the 378 isolates, potentially mutagenized but with the same genetic background, are likely a useful resource for future phenotyping and genome-wide association studies.

A polyploid admixed origin of beer yeasts derived from European and Asian wine populations.

Strains of Saccharomyces cerevisiae used to make beer, bread, and wine are genetically and phenotypically distinct from wild populations associated with trees. The origins of these domesticated populations are not always clear; human-associated migration and admixture with wild populations have had a strong impact on S. cerevisiae population structure. We examined the population genetic history of beer strains and found that ale strains and the S. cerevisiae portion of allotetraploid lager strains were derived from admixture between populations closely related to European grape wine strains and Asian rice wine strains. Similar to both lager and baking strains, ale strains are polyploid, providing them with a passive means of remaining isolated from other populations and providing us with a living relic of their ancestral hybridization. To reconstruct their polyploid origin, we phased the genomes of two ale strains and found ale haplotypes to both be recombinants between European and Asian alleles and to also contain novel alleles derived from extinct or as yet uncharacterized populations. We conclude that modern beer strains are the product of a historical melting pot of fermentation technology.

Temperature preference can bias parental genome retention during hybrid evolution.

Interspecific hybridization can introduce genetic variation that aids in adaptation to new or changing environments. Here, we investigate how hybrid adaptation to temperature and nutrient limitation may alter parental genome representation over time. We evolved Saccharomyces cerevisiae x Saccharomyces uvarum hybrids in nutrient-limited continuous culture at 15°C for 200 generations. In comparison to previous evolution experiments at 30°C, we identified a number of responses only observed in the colder temperature regime, including the loss of the S. cerevisiae allele in favor of the cryotolerant S. uvarum allele for several portions of the hybrid genome. In particular, we discovered a genotype by environment interaction in the form of a loss of heterozygosity event on chromosome XIII; which species' haplotype is lost or maintained is dependent on the parental species' temperature preference and the temperature at which the hybrid was evolved. We show that a large contribution to this directionality is due to a temperature dependent fitness benefit at a single locus, the high affinity phosphate transporter gene PHO84. This work helps shape our understanding of what forces impact genome evolution after hybridization, and how environmental conditions may promote or disfavor the persistence of hybrids over time.

Genome-wide consequences of deleting any single gene.

Loss or duplication of chromosome segments can lead to further genomic changes associated with cancer. However, it is not known whether only a select subset of genes is responsible for driving further changes. To determine whether perturbation of any given gene in a genome suffices to drive subsequent genetic changes, we analyzed the yeast knockout collection for secondary mutations of functional consequence. Unlike wild-type, most gene knockout strains were found to have one additional mutant gene affecting nutrient responses and/or heat-stress-induced cell death. Moreover, independent knockouts of the same gene often evolved mutations in the same secondary gene. Genome sequencing identified acquired mutations in several human tumor suppressor homologs. Thus, mutation of any single gene may cause a genomic imbalance, with consequences sufficient to drive adaptive genetic changes. This complicates genetic analyses but is a logical consequence of losing a functional unit originally acquired under pressure during evolution.

Variant Interpretation: Functional Assays to the Rescue.

Classical genetic approaches for interpreting variants, such as case-control or co-segregation studies, require finding many individuals with each variant. Because the overwhelming majority of variants are present in only a few living humans, this strategy has clear limits. Fully realizing the clinical potential of genetics requires that we accurately infer pathogenicity even for rare or private variation. Many computational approaches to predicting variant effects have been developed, but they can identify only a small fraction of pathogenic variants with the high confidence that is required in the clinic. Experimentally measuring a variant's functional consequences can provide clearer guidance, but individual assays performed only after the discovery of the variant are both time and resource intensive. Here, we discuss how multiplex assays of variant effect (MAVEs) can be used to measure the functional consequences of all possible variants in disease-relevant loci for a variety of molecular and cellular phenotypes. The resulting large-scale functional data can be combined with machine learning and clinical knowledge for the development of "lookup tables" of accurate pathogenicity predictions. A coordinated effort to produce, analyze, and disseminate large-scale functional data generated by multiplex assays could be essential to addressing the variant-interpretation crisis.

Effect of the Ala234Asp replacement in mitochondrial branched-chain amino acid aminotransferase on the production of BCAAs and fusel alcohols in yeast.

In the yeast Saccharomyces cerevisiae, the mitochondrial branched-chain amino acid (BCAA) aminotransferase Bat1 plays an important role in the synthesis of BCAAs (valine, leucine, and isoleucine). Our upcoming study (Large et al. bioRχiv. 10.1101/2020.06.26.166157, Large et al. 2020) will show that the heterozygous tetraploid beer yeast strain, Wyeast 1056, which natively has a variant causing one amino acid substitution of Ala234Asp in Bat1 on one of the four chromosomes, produced higher levels of BCAA-derived fusel alcohols in the brewer's wort medium than a derived strain lacking this mutation. Here, we investigated the physiological role of the A234D variant Bat1 in S. cerevisiae. Both bat1∆ and bat1A234D cells exhibited the same phenotypes relative to the wild-type Bat1 strain-namely, a repressive growth rate in the logarithmic phase; decreases in intracellular valine and leucine content in the logarithmic and stationary growth phases, respectively; an increase in fusel alcohol content in culture medium; and a decrease in the carbon dioxide productivity. These results indicate that amino acid change from Ala to Asp at position 234 led to a functional impairment of Bat1, although homology modeling suggests that Asp234 in the variant Bat1 did not inhibit enzymatic activity directly. KEY POINTS: • Yeast cells expressing Bat1A234D exhibited a slower growth phenotype. • The Val and Leu levels were decreased in yeast cells expressing Bat1A234D. • The A234D substitution causes a loss-of-function in Bat1. • The A234D substitution in Bat1 increased fusel alcohol production in yeast cells.

Differential paralog divergence modulates genome evolution across yeast species.

Evolutionary outcomes depend not only on the selective forces acting upon a species, but also on the genetic background. However, large timescales and uncertain historical selection pressures can make it difficult to discern such important background differences between species. Experimental evolution is one tool to compare evolutionary potential of known genotypes in a controlled environment. Here we utilized a highly reproducible evolutionary adaptation in Saccharomyces cerevisiae to investigate whether experimental evolution of other yeast species would select for similar adaptive mutations. We evolved populations of S. cerevisiae, S. paradoxus, S. mikatae, S. uvarum, and interspecific hybrids between S. uvarum and S. cerevisiae for ~200-500 generations in sulfate-limited continuous culture. Wild-type S. cerevisiae cultures invariably amplify the high affinity sulfate transporter gene, SUL1. However, while amplification of the SUL1 locus was detected in S. paradoxus and S. mikatae populations, S. uvarum cultures instead selected for amplification of the paralog, SUL2. We measured the relative fitness of strains bearing deletions and amplifications of both SUL genes from different species, confirming that, converse to S. cerevisiae, S. uvarum SUL2 contributes more to fitness in sulfate limitation than S. uvarum SUL1. By measuring the fitness and gene expression of chimeric promoter-ORF constructs, we were able to delineate the cause of this differential fitness effect primarily to the promoter of S. uvarum SUL1. Our data show evidence of differential sub-functionalization among the sulfate transporters across Saccharomyces species through recent changes in noncoding sequence. Furthermore, these results show a clear example of how such background differences due to paralog divergence can drive changes in genome evolution.