A resource of variant effect predictions of single nucleotide variants in model organisms

ProtVar: mapping and contextualizing human missense variation

Genomic variation can impact normal biological function in complex ways and so understanding variant effects requires a broad range of data to be coherently assimilated. Whilst the volume of human variant data and relevant annotations has increased, the corresponding increase in the breadth of participating fields, standards and versioning mean that moving between genomic, coding, protein and structure positions is increasingly complex. In turn this makes investigating variants in diverse formats and assimilating annotations from different resources challenging. ProtVar addresses these issues to facilitate the contextualization and interpretation of human missense variation with unparalleled flexibility and ease of accessibility for use by the broadest range of researchers. By precalculating all possible variants in the human proteome it offers near instantaneous mapping between all relevant data types. It also combines data and analyses from a plethora of resources to bring together genomic, protein sequence and function annotations as well as structural insights and predictions to better understand the likely effect of missense variation in humans. It is offered as an intuitive web server https://www.ebi.ac.uk/protvar where data can be explored and downloaded, and can be accessed programmatically via an API.

Guidelines for releasing a variant effect predictor

Computational methods for assessing the likely impacts of mutations, known as variant effect predictors (VEPs), are widely used in the assessment and interpretation of human genetic variation, as well as in other applications like protein engineering. Many different VEPs have been released to date, and there is tremendous variability in their underlying algorithms and outputs, and in the ways in which the methodologies and predictions are shared. This leads to considerable challenges for end users in knowing which VEPs to use and how to use them. Here, to address these issues, we provide guidelines and recommendations for the release of novel VEPs. Emphasising open-source availability, transparent methodologies, clear variant effect score interpretations, standardised scales, accessible predictions, and rigorous training data disclosure, we aim to improve the usability and interpretability of VEPs, and promote their integration into analysis and evaluation pipelines. We also provide a large, categorised list of currently available VEPs, aiming to facilitate the discovery and encourage the usage of novel methods within the scientific community.

BRCA1, BRCA2, and TP53 germline and somatic variants and clinicopathological characteristics of Brazilian patients with epithelial ovarian cancer

BACKGROUND

Approximately 3/4 of ovarian cancers are diagnosed in advanced stages, with the high-grade epithelial ovarian carcinoma (EOC) accounting for 90% of the cases. EOC present high genomic instability and somatic loss-of-function variants in genes associated with homologous recombination mutational repair pathway (HR), such as BRCA1 and BRCA2, and in TP53. The identification of germline variants in HR genes in EOC is relevant for treatment of platinum resistant tumors and relapsed tumors with therapies based in synthetic lethality such as PARP inhibitors. Patients with somatic variants in HR genes may also benefit from these therapies. In this work was analyzed the frequency of somatic variants in BRCA1, BRCA2, and TP53 in an EOC cohort of Brazilian patients, estimating the proportion of variants in tumoral tissue and their association with progression-free survival and overall survival.

METHODS

The study was conducted with paired blood/tumor samples from 56 patients. Germline and tumoral sequences of BRCA1, BRCA2, and TP53 were obtained by massive parallel sequencing. The HaplotypeCaller method was used for calling germline variants, and somatic variants were called with Mutect2.

RESULTS

A total of 26 germline variants were found, and seven patients presented germline pathogenic or likely pathogenic variants in BRCA1 or BRCA2. The analysis of tumoral tissue identified 52 somatic variants in 41 patients, being 43 somatic variants affecting or likely affecting protein functionality. Survival analyses showed that tumor staging was associated with overall survival (OS), while the presence of somatic mutation in TP53 was not associated with OS or progression-free survival.

CONCLUSION

Frequency of pathogenic or likely pathogenic germline variants in BRCA1 and BRCA2 (12.5%) was lower in comparison with other studies. TP53 was the most altered gene in tumors, with 62.5% presenting likely non-functional or non-functional somatic variants, while eight 14.2% presented likely non-functional or non-functional somatic variants in BRCA1 or BRCA2.

PanEffect: a pan-genome visualization tool for variant effects in maize

SUMMARY

Understanding the effects of genetic variants is crucial for accurately predicting traits and functional outcomes. Recent approaches have utilized artificial intelligence and protein language models to score all possible missense variant effects at the proteome level for a single genome, but a reliable tool is needed to explore these effects at the pan-genome level. To address this gap, we introduce a new tool called PanEffect. We implemented PanEffect at MaizeGDB to enable a comprehensive examination of the potential effects of coding variants across 50 maize genomes. The tool allows users to visualize over 550 million possible amino acid substitutions in the B73 maize reference genome and to observe the effects of the 2.3 million natural variations in the maize pan-genome. Each variant effect score, calculated from the Evolutionary Scale Modeling (ESM) protein language model, shows the log-likelihood ratio difference between B73 and all variants in the pan-genome. These scores are shown using heatmaps spanning benign outcomes to potential functional consequences. In addition, PanEffect displays secondary structures and functional domains along with the variant effects, offering additional functional and structural context. Using PanEffect, researchers now have a platform to explore protein variants and identify genetic targets for crop enhancement.

AVAILABILITY AND IMPLEMENTATION

The PanEffect code is freely available on GitHub (https://github.com/Maize-Genetics-and-Genomics-Database/PanEffect). A maize implementation of PanEffect and underlying datasets are available at MaizeGDB (https://www.maizegdb.org/effect/maize/).

Saturation genome editing of DDX3X clarifies pathogenicity of germline and somatic variation

Loss-of-function of DDX3X is a leading cause of neurodevelopmental disorders (NDD) in females. DDX3X is also a somatically mutated cancer driver gene proposed to have tumour promoting and suppressing effects. We perform saturation genome editing of DDX3X, testing in vitro the functional impact of 12,776 nucleotide variants. We identify 3432 functionally abnormal variants, in three distinct classes. We train a machine learning classifier to identify functionally abnormal variants of NDD-relevance. This classifier has at least 97% sensitivity and 99% specificity to detect variants pathogenic for NDD, substantially out-performing in silico predictors, and resolving up to 93% of variants of uncertain significance. Moreover, functionally-abnormal variants can account for almost all of the excess nonsynonymous DDX3X somatic mutations seen in DDX3X-driven cancers. Systematic maps of variant effects generated in experimentally tractable cell types have the potential to transform clinical interpretation of both germline and somatic disease-associated variation.

The regulatory landscape of the yeast phosphoproteome

The cellular ability to react to environmental fluctuations depends on signaling networks that are controlled by the dynamic activities of kinases and phosphatases. Here, to gain insight into these stress-responsive phosphorylation networks, we generated a quantitative mass spectrometry-based atlas of early phosphoproteomic responses in Saccharomyces cerevisiae exposed to 101 environmental and chemical perturbations. We report phosphosites on 59% of the yeast proteome, with 18% of the proteome harboring a phosphosite that is regulated within 5 min of stress exposure. We identify shared and perturbation-specific stress response programs, uncover loss of phosphorylation as an integral early event, and dissect the interconnected regulatory landscape of kinase-substrate networks, as we exemplify with target of rapamycin signaling. We further reveal functional organization principles of the stress-responsive phosphoproteome based on phosphorylation site motifs, kinase activities, subcellular localizations, shared functions and pathway intersections. This information-rich map of 25,000 regulated phosphosites advances our understanding of signaling networks.

The fitness cost of spurious phosphorylation

The fidelity of signal transduction requires the binding of regulatory molecules to their cognate targets. However, the crowded cell interior risks off-target interactions between proteins that are functionally unrelated. How such off-target interactions impact fitness is not generally known, but quantifying this is required to understand the constraints faced by cell systems as they evolve. Here, we use the model organism S. cerevisiae to inducibly express tyrosine kinases. Because yeast lacks bona fide tyrosine kinases, most of the resulting tyrosine phosphorylation is spurious. This provides a suitable system to measure the impact of artificial protein interactions on fitness. We engineered 44 yeast strains each expressing a tyrosine kinase, and quantitatively analysed their phosphoproteomes. This analysis resulted in ~30,000 phosphosites mapping to ~3,500 proteins. Examination of the fitness costs in each strain revealed a strong correlation between the number of spurious pY sites and decreased growth. Moreover, the analysis of pY effects on protein structure and on protein function revealed over 1000 pY events that we predict to be deleterious. However, we also find that a large number of the spurious pY sites have a negligible effect on fitness, possibly because of their low stoichiometry. This result is consistent with our evolutionary analyses demonstrating a lack of phosphotyrosine counter-selection in species with bona fide tyrosine kinases. Taken together, our results suggest that, alongside the risk for toxicity, the cell can tolerate a large degree of non-functional crosstalk as interaction networks evolve.

Cross-feeding affects the target of resistance evolution to an antifungal drug

Pathogenic fungi are a cause of growing concern. Developing an efficient and safe antifungal is challenging because of the similar biological properties of fungal and host cells. Consequently, there is an urgent need to better understand the mechanisms underlying antifungal resistance to prolong the efficacy of current molecules. A major step in this direction would be to be able to predict or even prevent the acquisition of resistance. We leverage the power of experimental evolution to quantify the diversity of paths to resistance to the antifungal 5-fluorocytosine (5-FC), commercially known as flucytosine. We generated hundreds of independent 5-FC resistant mutants derived from two genetic backgrounds from wild isolates of Saccharomyces cerevisiae. Through automated pin-spotting, whole-genome and amplicon sequencing, we identified the most likely causes of resistance for most strains. Approximately a third of all resistant mutants evolved resistance through a pleiotropic drug response, a potentially novel mechanism in response to 5-FC, marked by cross-resistance to fluconazole. These cross-resistant mutants are characterized by a loss of respiration and a strong tradeoff in drug-free media. For the majority of the remaining two thirds, resistance was acquired through loss-of-function mutations in FUR1, which encodes an important enzyme in the metabolism of 5-FC. We describe conditions in which mutations affecting this particular step of the metabolic pathway are favored over known resistance mutations affecting a step upstream, such as the well-known target cytosine deaminase encoded by FCY1. This observation suggests that ecological interactions may dictate the identity of resistance hotspots.

Discovering functionally important sites in proteins

Proteins play important roles in biology, biotechnology and pharmacology, and missense variants are a common cause of disease. Discovering functionally important sites in proteins is a central but difficult problem because of the lack of large, systematic data sets. Sequence conservation can highlight residues that are functionally important but is often convoluted with a signal for preserving structural stability. We here present a machine learning method to predict functional sites by combining statistical models for protein sequences with biophysical models of stability. We train the model using multiplexed experimental data on variant effects and validate it broadly. We show how the model can be used to discover active sites, as well as regulatory and binding sites. We illustrate the utility of the model by prospective prediction and subsequent experimental validation on the functional consequences of missense variants in HPRT1 which may cause Lesch-Nyhan syndrome, and pinpoint the molecular mechanisms by which they cause disease.

High-throughput deep learning variant effect prediction with Sequence UNET

Understanding coding mutations is important for many applications in biology and medicine but the vast mutation space makes comprehensive experimental characterisation impossible. Current predictors are often computationally intensive and difficult to scale, including recent deep learning models. We introduce Sequence UNET, a highly scalable deep learning architecture that classifies and predicts variant frequency from sequence alone using multi-scale representations from a fully convolutional compression/expansion architecture. It achieves comparable pathogenicity prediction to recent methods. We demonstrate scalability by analysing 8.3B variants in 904,134 proteins detected through large-scale proteomics. Sequence UNET runs on modest hardware with a simple Python package.

Rare variant aggregation in 148,508 exomes identifies genes associated with proxy dementia

Proxy phenotypes allow for the utilization of genetic data from large population cohorts to analyze late-onset diseases by using parental diagnoses as a proxy for genetic disease risk. Proxy phenotypes based on parental diagnosis status have been used in previous studies to identify common variants associated with Alzheimer's disease. As of yet, proxy phenotypes have not been used to identify genes associated with Alzheimer's disease through rare variants. Here we show that a proxy Alzheimer's disease/dementia phenotype can capture known Alzheimer's disease risk genes through rare variant aggregation. We generated a proxy Alzheimer's disease/dementia phenotype for 148,508 unrelated individuals of European ancestry in the UK biobank in order to perform exome-wide rare variant aggregation analyses to identify genes associated with proxy Alzheimer's disease/dementia. We identified four genes significantly associated with the proxy phenotype, three of which were significantly associated with proxy Alzheimer's disease/dementia in an independent replication cohort consisting of 197,506 unrelated individuals of European ancestry in the UK biobank. All three of the replicated genes have been previously associated with clinically diagnosed Alzheimer's disease (SORL1, TREM2, and TOMM40/APOE). We show that proxy Alzheimer's disease/dementia can be used to identify genes associated with Alzheimer's disease through rare variant aggregation.

Understanding Functional Redundancy and Promiscuity of Multidrug Transporters in E. coli under Lipophilic Cation Stress

Multidrug transporters (MDTs) are major contributors to microbial drug resistance and are further utilized for improving host phenotypes in biotechnological applications. Therefore, the identification of these MDTs and the understanding of their mechanisms of action in vivo are of great importance. However, their promiscuity and functional redundancy represent a major challenge towards their identification. Here, a multistep tolerance adaptive laboratory evolution (TALE) approach was leveraged to achieve this goal. Specifically, a wild-type E. coli K-12-MG1655 and its cognate knockout individual mutants ΔemrE, ΔtolC, and ΔacrB were evolved separately under increasing concentrations of two lipophilic cations, tetraphenylphosphonium (TPP⁺), and methyltriphenylphosphonium (MTPP⁺). The evolved strains showed a significant increase in MIC values of both cations and an apparent cross-cation resistance. Sequencing of all evolved mutants highlighted diverse mutational mechanisms that affect the activity of nine MDTs including acrB, mdtK, mdfA, acrE, emrD, tolC, acrA, mdtL, and mdtP. Besides regulatory mutations, several structural mutations were recognized in the proximal binding domain of acrB and the permeation pathways of both mdtK and mdfA. These details can aid in the rational design of MDT inhibitors to efficiently combat efflux-based drug resistance. Additionally, the TALE approach can be scaled to different microbes and molecules of medical and biotechnological relevance.

d-StructMAn: Containerized structural annotation on the scale from genetic variants to whole proteomes

Background

Structural annotation of genetic variants in the context of intermolecular interactions and protein stability can shed light onto mechanisms of disease-related phenotypes. Three-dimensional structures of related proteins in complexes with other proteins, nucleic acids, or ligands enrich such functional interpretation, since intermolecular interactions are well conserved in evolution.

Results

We present d-StructMAn, a novel computational method that enables structural annotation of local genetic variants, such as single-nucleotide variants and in-frame indels, and implements it in a highly efficient and user-friendly tool provided as a Docker container. Using d-StructMAn, we annotated several very large sets of human genetic variants, including all variants from ClinVar and all amino acid positions in the human proteome. We were able to provide annotation for more than 46% of positions in the human proteome representing over 60% proteins.

Conclusions

d-StructMAn is the first of its kind and a highly efficient tool for structural annotation of protein-coding genetic variation in the context of observed and potential intermolecular interactions. d-StructMAn is readily applicable to proteome-scale datasets and can be an instrumental building machine-learning tool for predicting genotype-to-phenotype relationships.

Discovering single-cell eQTLs from scRNA-seq data only

eQTL studies are essential for understanding genomic regulation. The effects of genetic variations on gene regulation are cell-type-specific and cellular-context-related, so studying eQTLs at a single-cell level is crucial. The ideal solution is to use both mutation and expression data from the same cells. However, the current technology of such paired data in single cells is still immature. We present a new method, eQTLsingle, to discover eQTLs only with single-cell RNA-seq (scRNA-seq) data, without genomic data. It detects mutations from scRNA-seq data and models gene expression of different genotypes with the zero-inflated negative binomial (ZINB) model to find associations between genotypes and phenotypes at the single-cell level. On a glioblastoma and gliomasphere scRNA-seq dataset, eQTLsingle discovered hundreds of cell-type-specific tumor-related eQTLs, most of which cannot be found in bulk eQTL studies. Detailed analyses on examples of the discovered eQTLs revealed important underlying regulatory mechanisms. eQTLsingle is a uniquely powerful tool for utilizing the vast scRNA-seq resources for single-cell eQTL studies, and it is available for free academic use at https://github.com/horsedayday/eQTLsingle.

Laboratory evolution of synthetic electron transport system variants reveals a larger metabolic respiratory system and its plasticity

The bacterial respiratory electron transport system (ETS) is branched to allow condition-specific modulation of energy metabolism. There is a detailed understanding of the structural and biochemical features of respiratory enzymes; however, a holistic examination of the system and its plasticity is lacking. Here we generate four strains of Escherichia coli harboring unbranched ETS that pump 1, 2, 3, or 4 proton(s) per electron and characterized them using a combination of synergistic methods (adaptive laboratory evolution, multi-omic analyses, and computation of proteome allocation). We report that: (a) all four ETS variants evolve to a similar optimized growth rate, and (b) the laboratory evolutions generate specific rewiring of major energy-generating pathways, coupled to the ETS, to optimize ATP production capability. We thus define an Aero-Type System (ATS), which is a generalization of the aerobic bioenergetics and is a metabolic systems biology description of respiration and its inherent plasticity.

The bacterial respiratory electron transport system (ETS) is branched to allow condition-specific modulation of energy metabolism. Here the authors examine the systems level properties of aerobic electron transport system using adaptive laboratory evolution and multi-omics analyses.

Evolutionary action of mutations reveals antimicrobial resistance genes in Escherichia coli

Since antibiotic development lags, we search for potential drug targets through directed evolution experiments. A challenge is that many resistance genes hide in a noisy mutational background as mutator clones emerge in the adaptive population. Here, to overcome this noise, we quantify the impact of mutations through evolutionary action (EA). After sequencing ciprofloxacin or colistin resistance strains grown under different mutational regimes, we find that an elevated sum of the evolutionary action of mutations in a gene identifies known resistance drivers. This EA integration approach also suggests new antibiotic resistance genes which are then shown to provide a fitness advantage in competition experiments. Moreover, EA integration analysis of clinical and environmental isolates of antibiotic resistant of E. coli identifies gene drivers of resistance where a standard approach fails. Together these results inform the genetic basis of de novo colistin resistance and support the robust discovery of phenotype-driving genes via the evolutionary action of genetic perturbations in fitness landscapes.

Predicting and interpreting large-scale mutagenesis data using analyses of protein stability and conservation

Understanding and predicting the functional consequences of single amino acid changes is central in many areas of protein science. Here, we collect and analyze experimental measurements of effects of >150,000 variants in 29 proteins. We use biophysical calculations to predict changes in stability for each variant and assess them in light of sequence conservation. We find that the sequence analyses give more accurate prediction of variant effects than predictions of stability and that about half of the variants that show loss of function do so due to stability effects. We construct a machine learning model to predict variant effects from protein structure and sequence alignments and show how the two sources of information support one another and enable mechanistic interpretations. Together, our results show how one can leverage large-scale experimental assessments of variant effects to gain deeper and general insights into the mechanisms that cause loss of function.

Multiple epistatic DNA variants in a single gene affect gene expression in trans

DNA variants that alter gene expression in trans are important sources of phenotypic variation. Nevertheless, the identity of trans-acting variants remains poorly understood. Single causal variants in several genes have been reported to affect the expression of numerous distant genes in trans. Whether these simple molecular architectures are representative of trans-acting variation is unknown. Here, we studied the large RAS signaling regulator gene IRA2, which contains variants with extensive trans-acting effects on gene expression in the yeast Saccharomyces cerevisiae. We used systematic CRISPR-based genome engineering and a sensitive phenotyping strategy to dissect causal variants to the nucleotide level. In contrast to the simple molecular architectures known so far, IRA2 contained at least seven causal nonsynonymous variants. The effects of these variants were modulated by nonadditive, epistatic interactions. Two variants at the 5'-end affected gene expression and growth only when combined with a third variant that also had no effect in isolation. Our findings indicate that the molecular basis of trans-acting genetic variation may be considerably more complex than previously appreciated.

OUP accepted manuscript

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.

Escherichia coli Data-Driven Strain Design Using Aggregated Adaptive Laboratory Evolution Mutational Data

Microbes are being engineered for an increasingly large and diverse set of applications. However, the designing of microbial genomes remains challenging due to the general complexity of biological systems. Adaptive Laboratory Evolution (ALE) leverages nature’s problem-solving processes to generate optimized genotypes currently inaccessible to rational methods. The large amount of public ALE data now represents a new opportunity for data-driven strain design. This study describes how novel strain designs, or genome sequences not yet observed in ALE experiments or published designs, can be extracted from aggregated ALE data and demonstrates this by designing, building, and testing three novel Escherichia coli strains with fitnesses comparable to ALE mutants. These designs were achieved through a meta-analysis of aggregated ALE mutations data (63 Escherichia coli K-12 MG1655 based ALE experiments, described by 93 unique environmental conditions, 357 independent evolutions, and 13 957 observed mutations), which additionally revealed global ALE mutation trends that inform on ALE-derived strain design principles. Such informative trends anticipate ALE-derived strain designs as largely gene-centric, as opposed to noncoding, and composed of a relatively small number of beneficial variants (approximately 6). These results demonstrate how strain design efforts can be enhanced by the meta-analysis of aggregated ALE data.

Aborting meiosis allows recombination in sterile diploid yeast hybrids

Hybrids between diverged lineages contain novel genetic combinations but an impaired meiosis often makes them evolutionary dead ends. Here, we explore to what extent an aborted meiosis followed by a return-to-growth (RTG) promotes recombination across a panel of 20 Saccharomyces cerevisiae and S. paradoxus diploid hybrids with different genomic structures and levels of sterility. Genome analyses of 275 clones reveal that RTG promotes recombination and generates extensive regions of loss-of-heterozygosity in sterile hybrids with either a defective meiosis or a heavily rearranged karyotype, whereas RTG recombination is reduced by high sequence divergence between parental subgenomes. The RTG recombination preferentially arises in regions with low local heterozygosity and near meiotic recombination hotspots. The loss-of-heterozygosity has a profound impact on sexual and asexual fitness, and enables genetic mapping of phenotypic differences in sterile lineages where linkage analysis would fail. We propose that RTG gives sterile yeast hybrids access to a natural route for genome recombination and adaptation.

Restoring fertility in yeast hybrids: Breeding and quantitative genetics of beneficial traits

Hybrids between species can harbor a combination of beneficial traits from each parent and may exhibit hybrid vigor, more readily adapting to new harsher environments. Interspecies hybrids are also sterile and therefore an evolutionary dead end unless fertility is restored, usually via auto-polyploidisation events. In the Saccharomyces genus, hybrids are readily found in nature and in industrial settings, where they have adapted to severe fermentative conditions. Due to their hybrid sterility, the development of new commercial yeast strains has so far been primarily conducted via selection methods rather than via further breeding. In this study, we overcame infertility by creating tetraploid intermediates of Saccharomyces interspecies hybrids to allow continuous multigenerational breeding. We incorporated nuclear and mitochondrial genetic diversity within each parental species, allowing for quantitative genetic analysis of traits exhibited by the hybrids and for nuclear-mitochondrial interactions to be assessed. Using pooled F12 generation segregants of different hybrids with extreme phenotype distributions, we identified quantitative trait loci (QTLs) for tolerance to high and low temperatures, high sugar concentration, high ethanol concentration, and acetic acid levels. We identified QTLs that are species specific, that are shared between species, as well as hybrid specific, in which the variants do not exhibit phenotypic differences in the original parental species. Moreover, we could distinguish between mitochondria-type-dependent and -independent traits. This study tackles the complexity of the genetic interactions and traits in hybrid species, bringing hybrids into the realm of full genetic analysis of diploid species, and paves the road for the biotechnological exploitation of yeast biodiversity.

A global map of associations between types of protein posttranslational modifications and human genetic diseases

There are >200 types of protein posttranslational modifications (PTMs) described in eukaryotes, each with unique proteome coverage and functions. We hypothesized that some genetic diseases may be caused by the removal of a specific type of PTMs by genomic variants and the consequent deregulation of particular functions. We collected >320,000 human PTMs representing 59 types and crossed them with >4M nonsynonymous DNA variants annotated with predicted pathogenicity and disease associations. We report >1.74M PTM-variant co-occurrences that an enrichment analysis distributed into 215 pairwise associations between 18 PTM types and 148 genetic diseases. Of them, 42% were not previously described. Removal of lysine acetylation exerts the most pronounced effect, and less studied PTM types such as S-glutathionylation or S-nitrosylation show relevance. Using pathogenicity predictions, we identified PTM sites that may produce particular diseases if prevented. Our results provide evidence of a substantial impact of PTM-specific removal on the pathogenesis of genetic diseases and phenotypes.

•

There is an enrichment of disease-associated nsSNVs preventing certain types of PTMs

•

We report 215 pairwise associations between 18 PTM types and 148 genetic diseases

•

The removal of lysine acetylation exerts the most pronounced effect

•

We report a set of PTM sites that may produce particular diseases if prevented

Exploring absent protein function in yeast: assaying post translational modification and human genetic variation

Yeast is a valuable eukaryotic model organism that has evolved many processes conserved up to humans, yet many protein functions, including certain DNA and protein modifications, are absent. It is this absence of protein function that is fundamental to approaches using yeast as an in vivo test system to investigate human proteins. Functionality of the heterologous expressed proteins is connected to a quantitative, selectable phenotype, enabling the systematic analyses of mechanisms and specificity of DNA modification, post-translational protein modifications as well as the impact of annotated cancer mutations and coding variation on protein activity and interaction. Through continuous improvements of yeast screening systems, this is increasingly carried out on a global scale using deep mutational scanning approaches. Here we discuss the applicability of yeast systems to investigate absent human protein function with a specific focus on the impact of protein variation on protein-protein interaction modulation.

MTR3D: identifying regions within protein tertiary structures under purifying selection

The identification of disease-causal variants is non-trivial. By mapping population variation from over 448,000 exome and genome sequences to over 81,000 experimental structures and homology models of the human proteome, we have calculated both regional intolerance to missense variation (Missense Tolerance Ratio, MTR), using a sliding window of 21–41 codons, and introduce a new 3D spatial intolerance to missense variation score (3D Missense Tolerance Ratio, MTR3D), using spheres of 5–8 Å. We show that the MTR3D is less biased by regions with limited data and more accurately identifies regions under purifying selection than estimates relying on the sequence alone. Intolerant regions were highly enriched for both ClinVar pathogenic and COSMIC somatic missense variants (Mann–Whitney U test P < 2.2 × 10⁻¹⁶). Further, we combine sequence- and spatial-based scores to generate a consensus score, MTRX, which distinguishes pathogenic from benign variants more accurately than either score separately (AUC = 0.85). The MTR3D server enables easy visualisation of population variation, MTR, MTR3D and MTRX scores across the entire gene and protein structure for >17,000 human genes and >42,000 alternative alternate transcripts, including both Ensembl and RefSeq transcripts. MTR3D is freely available by user-friendly web-interface and API at http://biosig.unimelb.edu.au/mtr3d/.

Integrated Multi-Omics Analysis of Mechanisms Underlying Yeast Ethanol Tolerance

For yeast cells, tolerance to high levels of ethanol is vital both in their natural environment and in industrially relevant conditions. We recently genotyped experimentally evolved yeast strains adapted to high levels of ethanol and identified mutations linked to ethanol tolerance. In this study, by integrating genomic sequencing data with quantitative proteomics profiles from six evolved strains (data set identifier PXD006631) and construction of protein interaction networks, we elucidate exactly how the genotype and phenotype are related at the molecular level. Our multi-omics approach points to the rewiring of numerous metabolic pathways affected by genomic and proteomic level changes, from energy-producing and lipid pathways to differential regulation of transposons and proteins involved in cell cycle progression. One of the key differences is found in the energy-producing metabolism, where the ancestral yeast strain responds to ethanol by switching to respiration and employing the mitochondrial electron transport chain. In contrast, the ethanol-adapted strains appear to have returned back to energy production mainly via glycolysis and ethanol fermentation, as supported by genomic and proteomic level changes. This work is relevant for synthetic biology where systems need to function under stressful conditions, as well as for industry and in cancer biology, where it is important to understand how the genotype relates to the phenotype.

Stress tolerance enhancement via SPT15 base editing in Saccharomyces cerevisiae

BACKGROUND

Saccharomyces cerevisiae is widely used in traditional brewing and modern fermentation industries to produce biofuels, chemicals and other bioproducts, but challenged by various harsh industrial conditions, such as hyperosmotic, thermal and ethanol stresses. Thus, its stress tolerance enhancement has been attracting broad interests. Recently, CRISPR/Cas-based genome editing technology offers unprecedented tools to explore genetic modifications and performance improvement of S. cerevisiae.

RESULTS

Here, we presented that the Target-AID (activation-induced cytidine deaminase) base editor of enabling C-to-T substitutions could be harnessed to generate in situ nucleotide changes on the S. cerevisiae genome, thereby introducing protein point mutations in cells. The general transcription factor gene SPT15 was targeted, and total 36 mutants with diversified stress tolerances were obtained. Among them, the 18 tolerant mutants against hyperosmotic, thermal and ethanol stresses showed more than 1.5-fold increases of fermentation capacities. These mutations were mainly enriched at the N-terminal region and the convex surface of the saddle-shaped structure of Spt15. Comparative transcriptome analysis of three most stress-tolerant (A140G, P169A and R238K) and two most stress-sensitive (S118L and L214V) mutants revealed common and distinctive impacted global transcription reprogramming and transcriptional regulatory hubs in response to stresses, and these five amino acid changes had different effects on the interactions of Spt15 with DNA and other proteins in the RNA Polymerase II transcription machinery according to protein structure alignment analysis.

CONCLUSIONS

Taken together, our results demonstrated that the Target-AID base editor provided a powerful tool for targeted in situ mutagenesis in S. cerevisiae and more potential targets of Spt15 residues for enhancing yeast stress tolerance.

How Bacteria Change after Exposure to Silver Nanoformulations: Analysis of the Genome and Outer Membrane Proteome

OBJECTIVE

the main purpose of this work was to compare the genetic and phenotypic changes of E. coli treated with silver nanoformulations (E. coli BW25113 wt, E. coli BW25113 AgR, E. coli J53, E. coli ATCC 11229 wt, E. coli ATCC 11229 var. S2 and E. coli ATCC 11229 var. S7). Silver, as the metal with promising antibacterial properties, is currently widely used in medicine and the biomedical industry, in both ionic and nanoparticles forms. Silver nanoformulations are usually considered as one type of antibacterial agent, but their physical and chemical properties determine the way of interactions with the bacterial cell, the mode of action, and the bacterial cell response to silver.

METHODS

the changes in the bacterial genome, resulting from the treatment of bacteria with various silver nanoformulations, were verified by analyzing of genes (selected with mutfunc) and their conservative and non-conservative mutations selected with BLOSUM62. The phenotype was verified using an outer membrane proteome analysis (OMP isolation, 2-DE electrophoresis, and MS protein identification).

RESULTS

the variety of genetic and phenotypic changes in E. coli strains depends on the type of silver used for bacteria treatment. The most changes were identified in E. coli ATCC 11229 treated with silver nanoformulation signed as S2 (E. coli ATCC 11229 var. S2). We pinpointed 39 genes encoding proteins located in the outer membrane, 40 genes of their regulators, and 22 genes related to other outer membrane structures, such as flagellum, fimbria, lipopolysaccharide (LPS), or exopolysaccharide in this strain. Optical density of OmpC protein in E. coli electropherograms decreased after exposure to silver nanoformulation S7 (noticed in E. coli ATCC 11229 var. S7), and increased after treatment with the other silver nanoformulations (SNF) marked as S2 (noticed in E. coli ATCC 11229 var. S2). Increase of FliC protein optical density was identified in turn after Ag+ treatment (noticed in E.coli AgR).

CONCLUSION

the results show that silver nanoformulations (SNF) exerts a selective pressure on bacteria causing both conservative and non-conservative mutations. The proteomic approach revealed that the levels of some proteins have changed after treatment with appropriate SNF.

Mapping the degradation pathway of a disease-linked aspartoacylase variant

Canavan disease is a severe progressive neurodegenerative disorder that is characterized by swelling and spongy degeneration of brain white matter. The disease is genetically linked to polymorphisms in the aspartoacylase (ASPA) gene, including the substitution C152W. ASPA C152W is associated with greatly reduced protein levels in cells, yet biophysical experiments suggest a wild-type like thermal stability. Here, we use ASPA C152W as a model to investigate the degradation pathway of a disease-causing protein variant. When we expressed ASPA C152W in Saccharomyces cerevisiae, we found a decreased steady state compared to wild-type ASPA as a result of increased proteasomal degradation. However, molecular dynamics simulations of ASPA C152W did not substantially deviate from wild-type ASPA, indicating that the native state is structurally preserved. Instead, we suggest that the C152W substitution interferes with the de novo folding pathway resulting in increased proteasomal degradation before reaching its stable conformation. Systematic mapping of the protein quality control components acting on misfolded and aggregation-prone species of C152W, revealed that the degradation is highly dependent on the molecular chaperone Hsp70, its co-chaperone Hsp110 as well as several quality control E3 ubiquitin-protein ligases, including Ubr1. In addition, the disaggregase Hsp104 facilitated refolding of aggregated ASPA C152W, while Cdc48 mediated degradation of insoluble ASPA protein. In human cells, ASPA C152W displayed increased proteasomal turnover that was similarly dependent on Hsp70 and Hsp110. Our findings underscore the use of yeast to determine the protein quality control components involved in the degradation of human pathogenic variants in order to identify potential therapeutic targets.

Understanding the Origins of Loss of Protein Function by Analyzing the Effects of Thousands of Variants on Activity and Abundance

Understanding and predicting how amino acid substitutions affect proteins are keys to our basic understanding of protein function and evolution. Amino acid changes may affect protein function in a number of ways including direct perturbations of activity or indirect effects on protein folding and stability. We have analyzed 6,749 experimentally determined variant effects from multiplexed assays on abundance and activity in two proteins (NUDT15 and PTEN) to quantify these effects and find that a third of the variants cause loss of function, and about half of loss-of-function variants also have low cellular abundance. We analyze the structural and mechanistic origins of loss of function and use the experimental data to find residues important for enzymatic activity. We performed computational analyses of protein stability and evolutionary conservation and show how we may predict positions where variants cause loss of activity or abundance. In this way, our results link thermodynamic stability and evolutionary conservation to experimental studies of different properties of protein fitness landscapes.

My personal mutanome: a computational genomic medicine platform for searching network perturbing alleles linking genotype to phenotype

Massive genome sequencing data have inspired new challenges in personalized treatments and facilitated oncological drug discovery. We present a comprehensive database, My Personal Mutanome (MPM), for accelerating the development of precision cancer medicine protocols. MPM contains 490,245 mutations from over 10,800 tumor exomes across 33 cancer types in The Cancer Genome Atlas mapped to 94,563 structure-resolved/predicted protein-protein interaction interfaces ("edgetic") and 311,022 functional sites ("nodetic"), including ligand-protein binding sites and 8 types of protein posttranslational modifications. In total, 8884 survival results and 1,271,132 drug responses are obtained for these mapped interactions. MPM is available at https://mutanome.lerner.ccf.org .

Folliculin variants linked to Birt-Hogg-Dubé syndrome are targeted for proteasomal degradation

Germline mutations in the folliculin (FLCN) tumor suppressor gene are linked to Birt-Hogg-Dubé (BHD) syndrome, a dominantly inherited genetic disease characterized by predisposition to fibrofolliculomas, lung cysts, and renal cancer. Most BHD-linked FLCN variants include large deletions and splice site aberrations predicted to cause loss of function. The mechanisms by which missense variants and short in-frame deletions in FLCN trigger disease are unknown. Here, we present an integrated computational and experimental study that reveals that the majority of such disease-causing FLCN variants cause loss of function due to proteasomal degradation of the encoded FLCN protein, rather than directly ablating FLCN function. Accordingly, several different single-site FLCN variants are present at strongly reduced levels in cells. In line with our finding that FLCN variants are protein quality control targets, several are also highly insoluble and fail to associate with the FLCN-binding partners FNIP1 and FNIP2. The lack of FLCN binding leads to rapid proteasomal degradation of FNIP1 and FNIP2. Half of the tested FLCN variants are mislocalized in cells, and one variant (ΔE510) forms perinuclear protein aggregates. A yeast-based stability screen revealed that the deubiquitylating enzyme Ubp15/USP7 and molecular chaperones regulate the turnover of the FLCN variants. Lowering the temperature led to a stabilization of two FLCN missense proteins, and for one (R362C), function was re-established at low temperature. In conclusion, we propose that most BHD-linked FLCN missense variants and small in-frame deletions operate by causing misfolding and degradation of the FLCN protein, and that stabilization and resulting restoration of function may hold therapeutic potential of certain disease-linked variants. Our computational saturation scan encompassing both missense variants and single site deletions in FLCN may allow classification of rare FLCN variants of uncertain clinical significance.

MISCAST: MIssense variant to protein StruCture Analysis web SuiTe

Human genome sequencing efforts have greatly expanded, and a plethora of missense variants identified both in patients and in the general population is now publicly accessible. Interpretation of the molecular-level effect of missense variants, however, remains challenging and requires a particular investigation of amino acid substitutions in the context of protein structure and function. Answers to questions like 'Is a variant perturbing a site involved in key macromolecular interactions and/or cellular signaling?', or 'Is a variant changing an amino acid located at the protein core or part of a cluster of known pathogenic mutations in 3D?' are crucial. Motivated by these needs, we developed MISCAST (missense variant to protein structure analysis web suite; http://miscast.broadinstitute.org/). MISCAST is an interactive and user-friendly web server to visualize and analyze missense variants in protein sequence and structure space. Additionally, a comprehensive set of protein structural and functional features have been aggregated in MISCAST from multiple databases, and displayed on structures alongside the variants to provide users with the biological context of the variant location in an integrated platform. We further made the annotated data and protein structures readily downloadable from MISCAST to foster advanced offline analysis of missense variants by a wide biological community.

Population-scale predictions of DPD and TPMT phenotypes using a quantitative pharmacogene-specific ensemble classifier

BACKGROUND

Inter-individual differences in dihydropyrimidine dehydrogenase (DPYD encoding DPD) and thiopurine S-methyltransferase (TPMT) activity are important predictors for fluoropyrimidine and thiopurine toxicity. While several variants in these genes are known to decrease enzyme activities, many additional genetic variations with unclear functional consequences have been identified, complicating informed clinical decision-making in the respective carriers.

METHODS

We used a novel pharmacogenetically trained ensemble classifier to analyse DPYD and TPMT genetic variability based on sequencing data from 138,842 individuals across eight populations.

RESULTS

The algorithm accurately predicted in vivo consequences of DPYD and TPMT variants (accuracy 91.4% compared to 95.3% in vitro). Further analysis showed high genetic complexity of DPD deficiency, advocating for sequencing-based DPYD profiling, whereas genotyping of four variants in TPMT was sufficient to explain >95% of phenotypic TPMT variability. Lastly, we provided population-scale profiles of ethnogeographic variability in DPD and TPMT phenotypes, and revealed striking interethnic differences in frequency and genetic constitution of DPD and TPMT deficiency.

CONCLUSION

These results provide the most comprehensive data set of DPYD and TPMT variability published to date with important implications for population-adjusted genetic profiling strategies of fluoropyrimidine and thiopurine risk factors and precision public health.

Population-scale proteome variation in human induced pluripotent stem cells

Human disease phenotypes are driven primarily by alterations in protein expression and/or function. To date, relatively little is known about the variability of the human proteome in populations and how this relates to variability in mRNA expression and to disease loci. Here, we present the first comprehensive proteomic analysis of human induced pluripotent stem cells (iPSC), a key cell type for disease modelling, analysing 202 iPSC lines derived from 151 donors, with integrated transcriptome and genomic sequence data from the same lines. We characterised the major genetic and non-genetic determinants of proteome variation across iPSC lines and assessed key regulatory mechanisms affecting variation in protein abundance. We identified 654 protein quantitative trait loci (pQTLs) in iPSCs, including disease-linked variants in protein-coding sequences and variants with trans regulatory effects. These include pQTL linked to GWAS variants that cannot be detected at the mRNA level, highlighting the utility of dissecting pQTL at peptide level resolution.

Causal mutations from adaptive laboratory evolution are outlined by multiple scales of genome annotations and condition-specificity

BACKGROUND

Adaptive Laboratory Evolution (ALE) has emerged as an experimental approach to discover mutations that confer phenotypic functions of interest. However, the task of finding and understanding all beneficial mutations of an ALE experiment remains an open challenge for the field. To provide for better results than traditional methods of ALE mutation analysis, this work applied enrichment methods to mutations described by a multiscale annotation framework and a consolidated set of ALE experiment conditions. A total of 25,321 unique genome annotations from various sources were leveraged to describe multiple scales of mutated features in a set of 35 Escherichia coli based ALE experiments. These experiments totalled 208 independent evolutions and 2641 mutations. Additionally, mutated features were statistically associated across a total of 43 unique experimental conditions to aid in deconvoluting mutation selection pressures.

RESULTS

Identifying potentially beneficial, or key, mutations was enhanced by seeking coding and non-coding genome features significantly enriched by mutations across multiple ALE replicates and scales of genome annotations. The median proportion of ALE experiment key mutations increased from 62%, with only small coding and non-coding features, to 71% with larger aggregate features. Understanding key mutations was enhanced by considering the functions of broader annotation types and the significantly associated conditions for key mutated features. The approaches developed here were used to find and characterize novel key mutations in two ALE experiments: one previously unpublished with Escherichia coli grown on glycerol as a carbon source and one previously published with Escherichia coli tolerized to high concentrations of L-serine.

CONCLUSIONS

The emergent adaptive strategies represented by sets of ALE mutations became more clear upon observing the aggregation of mutated features across small to large scale genome annotations. The clarification of mutation selection pressures among the many experimental conditions also helped bring these strategies to light. This work demonstrates how multiscale genome annotation frameworks and data-driven methods can help better characterize ALE mutations, and thus help elucidate the genotype-to-phenotype relationship of the studied organism.

VarMap: a web tool for mapping genomic coordinates to protein sequence and structure and retrieving protein structural annotations

Motivation

Understanding the protein structural context and patterning on proteins of genomic variants can help to separate benign from pathogenic variants and reveal molecular consequences. However, mapping genomic coordinates to protein structures is non-trivial, complicated by alternative splicing and transcript evidence.

Results

Here we present VarMap, a web tool for mapping a list of chromosome coordinates to canonical UniProt sequences and associated protein 3D structures, including validation checks, and annotating them with structural information.

Availability and implementation

https://www.ebi.ac.uk/thornton-srv/databases/VarMap.

Supplementary information

Supplementary data are available at Bioinformatics online.

Effects of a high-cultivation temperature on the physiology of three different Yarrowia lipolytica strains

Despite the increasing relevance, ranging from academic research to industrial applications, only a limited number of non-conventional, oleaginous Yarrowia lipolytica strains are characterized in detail. Therefore, we analyzed three strains in regard to their metabolic and physiological properties, especially with respect to important characteristics of a production strain. By investigating different cultivation conditions and media compositions, similarities and differences between the distinct strain backgrounds could be derived. Especially sugar alcohol production, as well as an agglomeration of cells were found to be connected with growth at high temperatures. In addition, sugar alcohol production was independent of high substrate concentrations under these conditions. To investigate the genotypic basis of particular traits, including growth characteristics and metabolite concentrations, genomic analysis were performed. We found sequence variations for one third of the annotated proteins but no obvious link to all phenotypic features.

Precision medicine - networks to the rescue

Genetic variants are often not predictive of the phenotypic outcome. Individuals carrying the same pathogenic variant, associated with Mendelian or complex disease, can manifest to different extents, from severe-to-mild to no disease. Improving the accuracy of predicted clinical manifestations of genetic variants has emerged as one of the biggest challenges in precision medicine, which can only be addressed by understanding the mechanisms underlying genotype-phenotype relationships. Efforts to understand the molecular basis of these relationships have identified complex systems of interacting biomolecules that underlie cellular function. Here, we review recent advances in how modeling cellular systems as networks of interacting proteins has fueled identification of disease-associated processes, delineation of underlying molecular mechanisms, and prediction of the pathogenicity of variants. This review is intended to be inspiring for clinicians, geneticists, and network biologists alike who aim to jointly advance our understanding of human disease and accelerate progress toward precision medicine.

Perturbing proteomes at single residue resolution using base editing

Base editors derived from CRISPR-Cas9 systems and DNA editing enzymes offer an unprecedented opportunity for the precise modification of genes, but have yet to be used at a genome-scale throughput. Here, we test the ability of the Target-AID base editor to systematically modify genes genome-wide by targeting yeast essential genes. We mutate around 17,000 individual sites in parallel across more than 1500 genes. We identify over 700 sites at which mutations have a significant impact on fitness. Using previously determined and preferred Target-AID mutational outcomes, we find that gRNAs with significant effects on fitness are enriched in variants predicted to be deleterious based on residue conservation and predicted protein destabilization. We identify key features influencing effective gRNAs in the context of base editing. Our results show that base editing is a powerful tool to identify key amino acid residues at the scale of proteomes.

VarSite: Disease variants and protein structure

VarSite is a web server mapping known disease-associated variants from UniProt and ClinVar, together with natural variants from gnomAD, onto protein 3D structures in the Protein Data Bank. The analyses are primarily image-based and provide both an overview for each human protein, as well as a report for any specific variant of interest. The information can be useful in assessing whether a given variant might be pathogenic or benign. The structural annotations for each position in the protein include protein secondary structure, interactions with ligand, metal, DNA/RNA, or other protein, and various measures of a given variant's possible impact on the protein's function. The 3D locations of the disease-associated variants can be viewed interactively via the 3dmol.js JavaScript viewer, as well as in RasMol and PyMOL. Users can search for specific variants, or sets of variants, by providing the DNA coordinates of the base change(s) of interest. Additionally, various agglomerative analyses are given, such as the mapping of disease and natural variants onto specific Pfam or CATH domains. The server is freely accessible to all at: https://www.ebi.ac.uk/thornton-srv/databases/VarSite.

The impact of the genetic background on gene deletion phenotypes in Saccharomyces cerevisiae

Loss-of-function (LoF) mutations associated with disease do not manifest equally in different individuals. The impact of the genetic background on the consequences of LoF mutations remains poorly characterized. Here, we systematically assessed the changes in gene deletion phenotypes for 3,786 gene knockouts in four Saccharomyces cerevisiae strains and 38 conditions. We observed 18.5% of deletion phenotypes changing between pairs of strains on average with a small fraction conserved in all four strains. Conditions causing higher wild-type growth differences and the deletion of pleiotropic genes showed above-average changes in phenotypes. In addition, we performed a genome-wide association study (GWAS) for growth under the same conditions for a panel of 925 yeast isolates. Gene-condition associations derived from GWAS were not enriched for genes with deletion phenotypes under the same conditions. However, cases where the results were congruent indicate the most likely mechanism underlying the GWAS signal. Overall, these results show a high degree of genetic background dependencies for LoF phenotypes.

DNA variants affecting the expression of numerous genes in trans have diverse mechanisms of action and evolutionary histories

DNA variants that alter gene expression contribute to variation in many phenotypic traits. In particular, trans-acting variants, which are often located on different chromosomes from the genes they affect, are an important source of heritable gene expression variation. However, our knowledge about the identity and mechanism of causal trans-acting variants remains limited. Here, we developed a fine-mapping strategy called CRISPR-Swap and dissected three expression quantitative trait locus (eQTL) hotspots known to alter the expression of numerous genes in trans in the yeast Saccharomyces cerevisiae. Causal variants were identified by engineering recombinant alleles and quantifying the effects of these alleles on the expression of a green fluorescent protein-tagged gene affected by the given locus in trans. We validated the effect of each variant on the expression of multiple genes by RNA-sequencing. The three variants differed in their molecular mechanism, the type of genes they reside in, and their distribution in natural populations. While a missense leucine-to-serine variant at position 63 in the transcription factor Oaf1 (L63S) was almost exclusively present in the reference laboratory strain, the two other variants were frequent among S. cerevisiae isolates. A causal missense variant in the glucose receptor Rgt2 (V539I) occurred at a poorly conserved amino acid residue and its effect was strongly dependent on the concentration of glucose in the culture medium. A noncoding variant in the conserved fatty acid regulated (FAR) element of the OLE1 promoter influenced the expression of the fatty acid desaturase Ole1 in cis and, by modulating the level of this essential enzyme, other genes in trans. The OAF1 and OLE1 variants showed a non-additive genetic interaction, and affected cellular lipid metabolism. These results demonstrate that the molecular basis of trans-regulatory variation is diverse, highlighting the challenges in predicting which natural genetic variants affect gene expression.

Differences in the DNA sequence of individual genomes contribute to differences in many traits, such as appearance, physiology, and the risk for common diseases. An important group of these DNA variants influences how individual genes across the genome are turned on or off. In this paper, we describe a strategy for identifying such “trans-acting” variants in different strains of baker’s yeast. We used this strategy to reveal three single DNA base changes that each influences the expression of dozens of genes. These three DNA variants were very different from each other. Two of them changed the protein sequence, one in a transcription factor and the other in a sugar sensor. The third changed the expression of an enzyme, a change that in turn caused other genes to alter their expression. One variant existed in only a few yeast isolates, while the other two existed in many isolates collected from around the world. This diversity of DNA variants that influence the expression of many other genes illustrates how difficult it is to predict which DNA variants in an individual’s genome will have effects on the organism.

Computational and cellular studies reveal structural destabilization and degradation of MLH1 variants in Lynch syndrome

Defective mismatch repair leads to increased mutation rates, and germline loss-of-function variants in the repair component MLH1 cause the hereditary cancer predisposition disorder known as Lynch syndrome. Early diagnosis is important, but complicated by many variants being of unknown significance. Here we show that a majority of the disease-linked MLH1 variants we studied are present at reduced cellular levels. We show that destabilized MLH1 variants are targeted for chaperone-assisted proteasomal degradation, resulting also in degradation of co-factors PMS1 and PMS2. In silico saturation mutagenesis and computational predictions of thermodynamic stability of MLH1 missense variants revealed a correlation between structural destabilization, reduced steady-state levels and loss-of-function. Thus, we suggest that loss of stability and cellular degradation is an important mechanism underlying many MLH1 variants in Lynch syndrome. Combined with analyses of conservation, the thermodynamic stability predictions separate disease-linked from benign MLH1 variants, and therefore hold potential for Lynch syndrome diagnostics.

Future Preventive Gene Therapy of Polygenic Diseases from a Population Genetics Perspective

With the accumulation of scientific knowledge of the genetic causes of common diseases and continuous advancement of gene-editing technologies, gene therapies to prevent polygenic diseases may soon become possible. This study endeavored to assess population genetics consequences of such therapies. Computer simulations were used to evaluate the heterogeneity in causal alleles for polygenic diseases that could exist among geographically distinct populations. The results show that although heterogeneity would not be easily detectable by epidemiological studies following population admixture, even significant heterogeneity would not impede the outcomes of preventive gene therapies. Preventive gene therapies designed to correct causal alleles to a naturally-occurring neutral state of nucleotides would lower the prevalence of polygenic early- to middle-age-onset diseases in proportion to the decreased population relative risk attributable to the edited alleles. The outcome would manifest differently for late-onset diseases, for which the therapies would result in a delayed disease onset and decreased lifetime risk; however, the lifetime risk would increase again with prolonging population life expectancy, which is a likely consequence of such therapies. If the preventive heritable gene therapies were to be applied on a large scale, the decreasing frequency of risk alleles in populations would reduce the disease risk or delay the age of onset, even with a fraction of the population receiving such therapies. With ongoing population admixture, all groups would benefit over generations.