Positive epistasis between disease-causing missense mutations and silent polymorphism with effect on mRNA translation velocity

Rescue of Mutant CFTR Channel Activity by Investigational Co-Potentiator Therapy

Background: The potentiator VX-770 (ivacaftor) has been approved as a monotherapy for over 95 cystic fibrosis (CF)-causing variants associated with gating/conductance defects of the CF transmembrane conductance regulator (CFTR) channel. However, despite its therapeutic success, VX-770 only partially restores CFTR activity for many of these variants, indicating they may benefit from the combination of potentiators exhibiting distinct mechanisms of action (i.e., co-potentiators). We previously identified LSO-24, a hydroxy-1,2,3-triazole-based compound, as a modest potentiator of p.Arg334Trp-CFTR, a variant with a conductance defect for which no modulator therapy is currently approved. Objective/Methods: We synthesized a new set of LSO-24 structure-based compounds, screened their effects on p.Arg334Trp-CFTR activity, and assessed the additivity of hit compounds to VX-770, ABBV-974, ABBV-3067, and apigenin. After validation by electrophysiological assays, the most promising hits were also assessed in cells expressing other variants with defective gating/conductance, namely p.Pro205Ser, p.Ser549Arg, p.Gly551Asp, p.Ser945Leu, and p.Gly1349Asp. Results: We found that five compounds were able to increase p.Arg334Trp-CFTR activity with similar efficacy, but slightly greater potency promoted by LSO-150 and LSO-153 (EC50: 1.01 and 1.26 μM, respectively). These two compounds also displayed a higher rescue of p.Arg334Trp-CFTR activity in combination with VX-770, ABBV-974, and ABBV-3067, but not with apigenin. When tested in cells expressing other CFTR variants, LSO-24 and its derivative LSO-150 increased CFTR activity for the variants p.Ser549Arg, p.Gly551Asp, and p.Ser945Leu with a further effect in combination with VX-770 or ABBV-3067. No potentiator was able to rescue CFTR activity in p.Pro205Ser-expressing cells, while p.Gly1349Asp-CFTR responded to VX-770 and ABBV-3067 but not to LSO-24 or LSO-150. Conclusions: Our data suggest that these new potentiators might share a common mechanism with apigenin, which is conceivably distinct from that of VX-770 and ABBV-3067. The additive rescue of p.Arg334Trp-, p.Ser549Arg-, p.Gly551Asp-, and p.Ser945Leu-CFTR also indicates that these variants could benefit from the development of a co-potentiator therapy.

Distinct network patterns emerge from Cartesian and XOR epistasis models: a comparative network science analysis

BACKGROUND

Epistasis, the phenomenon where the effect of one gene (or variant) is masked or modified by one or more other genes, significantly contributes to the phenotypic variance of complex traits. Traditionally, epistasis has been modeled using the Cartesian epistatic model, a multiplicative approach based on standard statistical regression. However, a recent study investigating epistasis in obesity-related traits has identified potential limitations of the Cartesian epistatic model, revealing that it likely only detects a fraction of the genetic interactions occurring in natural systems. In contrast, the exclusive-or (XOR) epistatic model has shown promise in detecting a broader range of epistatic interactions and revealing more biologically relevant functions associated with interacting variants. To investigate whether the XOR epistatic model also forms distinct network structures compared to the Cartesian model, we applied network science to examine genetic interactions underlying body mass index (BMI) in rats (Rattus norvegicus).

RESULTS

Our comparative analysis of XOR and Cartesian epistatic models in rats reveals distinct topological characteristics. The XOR model exhibits enhanced sensitivity to epistatic interactions between the network communities found in the Cartesian epistatic network, facilitating the identification of novel trait-related biological functions via community-based enrichment analysis. Additionally, the XOR network features triangle network motifs, indicative of higher-order epistatic interactions. This research also evaluates the impact of linkage disequilibrium (LD)-based edge pruning on network-based epistasis analysis, finding that LD-based edge pruning may lead to increased network fragmentation, which may hinder the effectiveness of network analysis for the investigation of epistasis. We confirmed through network permutation analysis that most XOR and Cartesian epistatic networks derived from the data display distinct structural properties compared to randomly shuffled networks.

CONCLUSIONS

Collectively, these findings highlight the XOR model's ability to uncover meaningful biological associations and higher-order epistasis derived from lower-order network topologies. The introduction of community-based enrichment analysis and motif-based epistatic discovery emphasize network science as a critical approach for advancing epistasis research and understanding complex genetic architectures.

A fundamental and theoretical framework for mutation interactions and epistasis

Many pathological conditions are a result of intragenic epistasis; however, there are ambiguities in current epistasis models. Herein, the new Mutation Interaction Spectrum model defines a discrete outcome, named a Mutation Interaction, for each double point mutation in a gene and its component single mutations. The model is a universal genetic model of all types of mutation interactions and their functional outcomes and is derived from digital logic, commonly used in electrical engineering. Mutation interactions are normally classified as positive and negative epistasis. The model logics unifies common genetic relationships into one model, normalizing biological nomenclature, and disambiguates them with the 16 possible logic-based interactions. The model was tested by assaying transcriptional activity induced by HIV-1 Tat protein, for a random sampling of 3429 double mutations and all 1615 single mutations. All possible types of logic were observed for the Tat mutation interactions.

Ribosomal frameshifting selectively modulates the assembly, function, and pharmacological rescue of a misfolded CFTR variant

The cotranslational misfolding of the cystic fibrosis transmembrane conductance regulator chloride channel (CFTR) plays a central role in the molecular basis of CF. The misfolding of the most common CF variant (ΔF508) remodels both the translational regulation and quality control of CFTR. Nevertheless, it is unclear how the misassembly of the nascent polypeptide may directly influence the activity of the translation machinery. In this work, we identify a structural motif within the CFTR transcript that stimulates efficient -1 ribosomal frameshifting and triggers the premature termination of translation. Though this motif does not appear to impact the interactome of wild-type CFTR, silent mutations that disrupt this RNA structure alter the association of nascent ΔF508 CFTR with numerous translation and quality control proteins. Moreover, disrupting this RNA structure enhances the functional gating of the ΔF508 CFTR channel at the plasma membrane and its pharmacological rescue by the CFTR modulators contained in the CF drug Trikafta. The effects of the RNA structure on ΔF508 CFTR appear to be attenuated in the absence of the ER membrane protein complex, which was previously found to modulate ribosome collisions during "preemptive quality control" of a misfolded CFTR homolog. Together, our results reveal that ribosomal frameshifting selectively modulates the assembly, function, and pharmacological rescue of a misfolded CFTR variant. These findings suggest that interactions between the nascent chain, quality control machinery, and ribosome may dynamically modulate ribosomal frameshifting in order to tune the processivity of translation in response to cotranslational misfolding.

Ribosomal Frameshifting Selectively Modulates the Assembly, Function, and Pharmacological Rescue of a Misfolded CFTR Variant

The cotranslational misfolding of the cystic fibrosis transmembrane conductance regulator chloride channel (CFTR) plays a central role in the molecular basis of cystic fibrosis (CF). The misfolding of the most common CF variant (ΔF508) remodels both the translational regulation and quality control of CFTR. Nevertheless, it is unclear how the misassembly of the nascent polypeptide may directly influence the activity of the translation machinery. In this work, we identify a structural motif within the CFTR transcript that stimulates efficient -1 ribosomal frameshifting and triggers the premature termination of translation. Though this motif does not appear to impact the interactome of wild-type CFTR, silent mutations that disrupt this RNA structure alter the association of nascent ΔF508 CFTR with numerous translation and quality control proteins. Moreover, disrupting this RNA structure enhances the functional gating of the ΔF508 CFTR channel at the plasma membrane and its pharmacological rescue by the CFTR modulators contained in the CF drug Trikafta. The effects of the RNA structure on ΔF508 CFTR appear to be attenuated in the absence of the ER membrane protein complex (EMC), which was previously found to modulate ribosome collisions during "preemptive quality control" of a misfolded CFTR homolog. Together, our results reveal that ribosomal frameshifting selectively modulates the assembly, function, and pharmacological rescue of a misfolded CFTR variant. These findings suggest interactions between the nascent chain, quality control machinery, and ribosome may dynamically modulate ribosomal frameshifting in order to tune the processivity of translation in response to cotranslational misfolding.

The Folding Pathway of ABC Transporter CFTR: Effective and Robust

De novo protein folding into a native three-dimensional structure is indispensable for biological function, is instructed by its amino acid sequence, and occurs along a vectorial trajectory. The human proteome contains thousands of membrane-spanning proteins, whose biosynthesis begins on endoplasmic reticulum-associated ribosomes. Nearly half of all membrane proteins traverse the membrane more than once, including therapeutically important protein families such as solute carriers, G-protein-coupled receptors, and ABC transporters. These mediate a variety of functions like signal transduction and solute transport and are often of vital importance for cell function and tissue homeostasis. Missense mutations in multispan membrane proteins can lead to misfolding and cause disease; an example is the ABC transporter Cystic Fibrosis Transmembrane Conductance Regulator (CFTR). Even though our understanding of multispan membrane-protein folding still is rather rudimental, the cumulative knowledge of 20 years of basic research on CFTR folding has led to development of drugs that modulate the misfolded protein. This has provided the prospect of a life without CF to the vast majority of patients. In this review we describe our understanding of the folding pathway of CFTR in cells, which is modular and tolerates many defects, making it effective and robust. We address how modulator drugs affect folding and function of CFTR, and distinguish protein stability from its folding process. Since the domain architecture of (mammalian) ABC transporters are highly conserved, we anticipate that the insights we discuss here for folding of CFTR may lay the groundwork for understanding the general rules of ABC-transporter folding.

The Effects of Codon Usage on Protein Structure and Folding

The rate of protein synthesis is slower than many folding reactions and varies depending on the synonymous codons encoding the protein sequence. Synonymous codon substitutions thus have the potential to regulate cotranslational protein folding mechanisms, and a growing number of proteins have been identified with folding mechanisms sensitive to codon usage. Typically, these proteins have complex folding pathways and kinetically stable native structures. Kinetically stable proteins may fold only once over their lifetime, and thus, codon-mediated regulation of the pioneer round of protein folding can have a lasting impact. Supporting an important role for codon usage in folding, conserved patterns of codon usage appear in homologous gene families, hinting at selection. Despite these exciting developments, there remains few experimental methods capable of quantifying translation elongation rates and cotranslational folding mechanisms in the cell, which challenges the development of a predictive understanding of how biology uses codons to regulate protein folding.

Navigating the complexities of multi-domain protein folding

Proteome complexity has expanded tremendously over evolutionary time, enabling biological diversification. Much of this complexity is achieved by combining a limited set of structural units into long polypeptides. This widely used evolutionary strategy poses challenges for folding of the resulting multi-domain proteins. As a consequence, their folding differs from that of small single-domain proteins, which generally fold quickly and reversibly. Co-translational processes and chaperone interactions are important aspects of multi-domain protein folding. In this review, we discuss some of the recent experimental progress toward understanding these processes.

Racial inequities and rare CFTR variants: Impact on cystic fibrosis diagnosis and treatment

Cystic fibrosis (CF) has been traditionally viewed as a disease that affects White individuals. However, CF occurs among all races, ethnicities, and geographic ancestries. The disorder results from mutations in the CF transmembrane conductance regulator (CFTR). Varying incidence of CF is reported among Black, Indigenous, and People of Color (BIPOC), who typically exhibit worse clinical outcomes. These populations are more likely to carry rare CFTR variants omitted from newborn screening panels, leading to disparities in care such as delayed diagnosis and treatment. In this study, we present a case-in-point describing an individual of Gambian descent identified with CF. Patient genotype includes a premature termination codon (PTC) (c.2353C>T) and previously undescribed single nucleotide deletion (c.1970delG), arguing against effectiveness of currently available CFTR modulator-based interventions. Strategies for overcoming these two variants will likely include combinations of PTC suppressors, nonsense mediated decay inhibitors, and/or alternative approaches (e.g. gene therapy). Investigations such as the present study establish a foundation from which therapeutic treatments may be developed. Importantly, c.2353C>T and c.1970delG were not detected in the patient by traditional CFTR screening panels, which include an implicit racial and ethnic diagnostic bias as these tests are comprised of mutations largely observed in people of European ancestry. We suggest that next-generation sequencing of CFTR should be utilized to confirm or exclude a CF diagnosis, in order to equitably serve BIPOC individuals. Additional epidemiologic data, basic science investigations, and translational work are imperative for improving understanding of disease prevalence and progression, CFTR variant frequency, genotype-phenotype correlation, pharmacologic responsiveness, and personalized medicine approaches for patients with African ancestry and other historically understudied geographic lineages.

Distinct Network Patterns Emerge from Cartesian and XOR Epistasis Models: A Comparative Network Science Analysis

Background

Epistasis, the phenomenon where the effect of one gene (or variant) is masked or modified by one or more other genes, can significantly contribute to the observed phenotypic variance of complex traits. To date, it has been generally assumed that genetic interactions can be detected using a Cartesian, or multiplicative, interaction model commonly utilized in standard regression approaches. However, a recent study investigating epistasis in obesity-related traits in rats and mice has identified potential limitations of the Cartesian model, revealing that it only detects some of the genetic interactions occurring in these systems. By applying an alternative approach, the exclusive-or (XOR) model, the researchers detected a greater number of epistatic interactions and identified more biologically relevant ontological terms associated with the interacting loci. This suggests that the XOR model may provide a more comprehensive understanding of epistasis in these species and phenotypes. To further explore these findings and determine if different interaction models also make up distinct epistatic networks, we leverage network science to provide a more comprehensive view into the genetic interactions underlying BMI in this system.

Results

Our comparative analysis of networks derived from Cartesian and XOR interaction models in rats (Rattus norvegicus) uncovers distinct topological characteristics for each model-derived network. Notably, we discover that networks based on the XOR model exhibit an enhanced sensitivity to epistatic interactions. This sensitivity enables the identification of network communities, revealing novel trait-related biological functions through enrichment analysis. Furthermore, we identify triangle network motifs in the XOR epistatic network, suggestive of higher-order epistasis, based on the topology of lower-order epistasis.

Conclusions

These findings highlight the XOR model's ability to uncover meaningful biological associations as well as higher-order epistasis from lower-order epistatic networks. Additionally, our results demonstrate that network approaches not only enhance epistasis detection capabilities but also provide more nuanced understandings of genetic architectures underlying complex traits. The identification of community structures and motifs within these distinct networks, especially in XOR, points to the potential for network science to aid in the discovery of novel genetic pathways and regulatory networks. Such insights are important for advancing our understanding of phenotype-genotype relationships.

Evolving precision: rRNA expansion segment 7S modulates translation velocity and accuracy in eukaryal ribosomes

Ribosome-enhanced translational miscoding of the genetic code causes protein dysfunction and loss of cellular fitness. During evolution, open reading frame length increased, necessitating mechanisms for enhanced translation fidelity. Indeed, eukaryal ribosomes are more accurate than bacterial counterparts, despite their virtually identical, conserved active centers. During the evolution of eukaryotic organisms ribosome expansions at the rRNA and protein level occurred, which potentially increases the options for translation regulation and cotranslational events. Here we tested the hypothesis that ribosomal RNA expansions can modulate the core function of the ribosome, faithful protein synthesis. We demonstrate that a short expansion segment present in all eukaryotes' small subunit, ES7S, is crucial for accurate protein synthesis as its presence adjusts codon-specific velocities and guarantees high levels of cognate tRNA selection. Deletion of ES7S in yeast enhances mistranslation and causes protein destabilization and aggregation, dramatically reducing cellular fitness. Removal of ES7S did not alter ribosome architecture but altered the structural dynamics of inter-subunit bridges thus affecting A-tRNA selection. Exchanging the yeast ES7S sequence with the human ES7S increases accuracy whereas shortening causes the opposite effect. Our study demonstrates that ES7S provided eukaryal ribosomes with higher accuracy without perturbing the structurally conserved decoding center.

Translation velocity determines the efficacy of engineered suppressor tRNAs on pathogenic nonsense mutations

Nonsense mutations - the underlying cause of approximately 11% of all genetic diseases - prematurely terminate protein synthesis by mutating a sense codon to a premature stop or termination codon (PTC). An emerging therapeutic strategy to suppress nonsense defects is to engineer sense-codon decoding tRNAs to readthrough and restore translation at PTCs. However, the readthrough efficiency of the engineered suppressor tRNAs (sup-tRNAs) largely varies in a tissue- and sequence context-dependent manner and has not yet yielded optimal clinical efficacy for many nonsense mutations. Here, we systematically analyze the suppression efficacy at various pathogenic nonsense mutations. We discover that the translation velocity of the sequence upstream of PTCs modulates the sup-tRNA readthrough efficacy. The PTCs most refractory to suppression are embedded in a sequence context translated with an abrupt reversal of the translation speed leading to ribosomal collisions. Moreover, modeling translation velocity using Ribo-seq data can accurately predict the suppression efficacy at PTCs. These results reveal previously unknown molecular signatures contributing to genotype-phenotype relationships and treatment-response heterogeneity, and provide the framework for the development of personalized tRNA-based gene therapies.

Pan-cellular organelles and suborganelles-from common functions to cellular diversity?

Cell diversification is at the base of increasing multicellular organism complexity in phylogeny achieved during ontogeny. However, there are also functions common to all cells, such as cell division, cell migration, translation, endocytosis, exocytosis, etc. Here we revisit the organelles involved in such common functions, reviewing recent evidence of unexpected differences of proteins at these organelles. For instance, centrosomes or mitochondria differ significantly in their protein composition in different, sometimes closely related, cell types. This has relevance for development and disease. Particularly striking is the high amount and diversity of RNA-binding proteins at these and other organelles, which brings us to review the evidence for RNA at different organelles and suborganelles. We include a discussion about (sub)organelles involved in translation, such as the nucleolus and ribosomes, for which unexpected cell type-specific diversity has also been reported. We propose here that the heterogeneity of these organelles and compartments represents a novel mechanism for regulating cell diversity. One reason is that protein functions can be multiplied by their different contributions in distinct organelles, as also exemplified by proteins with moonlighting function. The specialized organelles still perform pan-cellular functions but in a cell type-specific mode, as discussed here for centrosomes, mitochondria, vesicles, and other organelles. These can serve as regulatory hubs for the storage and transport of specific and functionally important regulators. In this way, they can control cell differentiation, plasticity, and survival. We further include examples highlighting the relevance for disease and propose to examine organelles in many more cell types for their possible differences with functional relevance.

Interaction models matter: an efficient, flexible computational framework for model-specific investigation of epistasis

PURPOSE

Epistasis, the interaction between two or more genes, is integral to the study of genetics and is present throughout nature. Yet, it is seldom fully explored as most approaches primarily focus on single-locus effects, partly because analyzing all pairwise and higher-order interactions requires significant computational resources. Furthermore, existing methods for epistasis detection only consider a Cartesian (multiplicative) model for interaction terms. This is likely limiting as epistatic interactions can evolve to produce varied relationships between genetic loci, some complex and not linearly separable.

METHODS

We present new algorithms for the interaction coefficients for standard regression models for epistasis that permit many varied models for the interaction terms for loci and efficient memory usage. The algorithms are given for two-way and three-way epistasis and may be generalized to higher order epistasis. Statistical tests for the interaction coefficients are also provided. We also present an efficient matrix based algorithm for permutation testing for two-way epistasis. We offer a proof and experimental evidence that methods that look for epistasis only at loci that have main effects may not be justified. Given the computational efficiency of the algorithm, we applied the method to a rat data set and mouse data set, with at least 10,000 loci and 1,000 samples each, using the standard Cartesian model and the XOR model to explore body mass index.

RESULTS

This study reveals that although many of the loci found to exhibit significant statistical epistasis overlap between models in rats, the pairs are mostly distinct. Further, the XOR model found greater evidence for statistical epistasis in many more pairs of loci in both data sets with almost all significant epistasis in mice identified using XOR. In the rat data set, loci involved in epistasis under the XOR model are enriched for biologically relevant pathways.

CONCLUSION

Our results in both species show that many biologically relevant epistatic relationships would have been undetected if only one interaction model was applied, providing evidence that varied interaction models should be implemented to explore epistatic interactions that occur in living systems.

Ribosome inactivation regulates translation elongation in neurons

Cellular plasticity is crucial for adapting to ever-changing stimuli. As a result, cells consistently reshape their translatome, and, consequently, their proteome. The control of translational activity has been thoroughly examined at the stage of translation initiation. However, the regulation of ribosome speed in cells is widely unknown. In this study, we utilized a timed ribosome runoff approach, along with proteomics and transmission electron microscopy, to investigate global translation kinetics in cells. We found that ribosome speeds vary among various cell types, such as astrocytes, induced pluripotent human stem cells, human neural stem cells, and human and rat neurons. Of all cell types studied, mature cortical neurons exhibit the highest rate of translation. This finding is particularly remarkable because mature cortical neurons express the eukaryotic elongation factor 2 (eEF2) at lower levels than other cell types. Neurons solve this conundrum by inactivating a fraction of their ribosomes. As a result, the increase in eEF2 levels leads to a reduction of inactive ribosomes and an enhancement of active ones. Processes that alter the demand for active ribosomes, like neuronal excitation, cause increased inactivation of redundant ribosomes in an eEF2-dependent manner. Our data suggest a novel regulatory mechanism in which neurons dynamically inactivate ribosomes to facilitate translational remodeling. These findings have important implications for developmental brain disorders characterized by, among other things, aberrant translation.

Effect of mRNA/tRNA mutations on translation speed: Implications for human diseases

Recent discoveries establish tRNAs as central regulators of mRNA translation dynamics, and therefore cotranslational folding and function of the encoded protein. The tRNA pool, whose composition and abundance change in a cell- and tissue-dependent manner, is the main factor which determines mRNA translation velocity. In this review, we discuss a group of pathogenic mutations, in the coding sequences of either protein-coding genes or in tRNA genes, that alter mRNA translation dynamics. We also summarize advances in tRNA biology that have uncovered how variations in tRNA levels on account of genetic mutations affect protein folding and function, and thereby contribute to phenotypic diversity in clinical manifestations.

Gene Interactions in Human Disease Studies-Evidence Is Mounting

Despite monumental advances in molecular technology to generate genome sequence data at scale, there is still a considerable proportion of heritability in most complex diseases that remains unexplained. Because many of the discoveries have been single-nucleotide variants with small to moderate effects on disease, the functional implication of many of the variants is still unknown and, thus, we have limited new drug targets and therapeutics. We, and many others, posit that one primary factor that has limited our ability to identify novel drug targets from genome-wide association studies may be due to gene interactions (epistasis), gene-environment interactions, network/pathway effects, or multiomic relationships. We propose that many of these complex models explain much of the underlying genetic architecture of complex disease. In this review, we discuss the evidence from multiple research avenues, ranging from pairs of alleles to multiomic integration studies and pharmacogenomics, that supports the need for further investigation of gene interactions (or epistasis) in genetic and genomic studies of human disease. Our goal is to catalog the mounting evidence for epistasis in genetic studies and the connections between genetic interactions and human health and disease that could enable precision medicine of the future.

Multi-animal-model study reveals mutations in neural plasticity and nociception genes linked to excessive alcohol drinking

BACKGROUND

The basis for familial alcohol use disorder (AUD) remains an enigma due to various biological and societal confounds. The present study used three of the most adopted and documented rat models, combining the alcohol-preferring/non-alcohol-preferring (P/NP) lines and high alcohol-drinking/low alcohol-drinking (HAD/LAD) replicated lines, of AUD as examined through the lens of whole genomic analyses.

METHODS

We used complete genome sequencing of the P/NP lines and previously published sequences of the HAD/LAD replicates to enhance the discovery of variants associated with AUD and to remove confounding with genetic background and random genetic drift. Specifically, we used high-order statistical methods to search for genetic variants whose frequency changes in whole sets of gene ontologies corresponded with phenotypic changes in the direction of selection, that is, ethanol-drinking preference.

RESULTS

Our first finding was that in addition to variants causing translational changes, the principal genetic changes associated with drinking predisposition were silent mutations and mutations in the 3' untranslated regions (3'UTR) of genes. Neither of these types of mutations alters the amino acid sequence of the translated protein but they influence both the rate and conformation of gene transcription, including its stability and posttranslational events that alter gene efficacy. This finding argues for refocusing human genomic studies on changes in gene efficacy. Among the key ontologies identified were the central genes associated with the Na+ voltage-gated channels of neurons and glia (including the Scn1a, Scn2a, Scn2b, Scn3a, Scn7a, and Scn9a subtypes) and excitatory glutamatergic secretion (including Grm2 and Myo6), both of which are essential in neuroplasticity. In addition, we identified "Nociception or Sensory Perception of Pain," which contained variants in nociception (Arrb1, Ccl3, Ephb1) and enlist sodium (Scn1a, Scn2a, Scn2b, Scn3a, Scn7a), pain activation (Scn9a), and potassium channel (Kcna1) genes.

CONCLUSION

The multi-model analyses used herein reduced the confounding effects of random drift and the "founders" genetic background. The most differentiated bidirectionally selected genes across all three animal models were Scn9a, Scn1a, and Kcna, all of which are annotated in the nociception ontology. The complexity of neuroplasticity and nociception adds strength to the hypothesis that neuroplasticity and pain (physical or psychological) are prominent phenotypes genetically linked to the development of AUD.

Monosomes buffer translational stress to allow for active ribosome elongation

Introduction: The synthesis of proteins is a fundamental process in the life-span of all cells. The activation of ribosomes on transcripts is the starting signal for elongation and, in turn, the translation of an mRNA. Thereby, most mRNAs circulate between single (monosomes) and multi ribosomal particles (polysomes), a process that defines their translational activity. The interplay between monosomes and polysomes is thought to crucially impact translation rate. How monosomes and polysomes are balanced during stress remains, however, elusive. Methods: Here, we set out to investigate the monosome and polysome levels as well as their kinetics under different translational stress conditions including mTOR inhibition, downregulation of the eukaryotic elongation factor 2 (eEF2) and amino acid depletion. Results: By using a timed ribosome runoff approach in combination with polysome profiling, we found that the used translational stressors show very distinct effects on translation. However, they all had in common that the activity of monosomes was preferentially affected. This adaptation seems to be needed for sufficient translation elongation. Even under harsh conditions such as amino acid starvation, we detected active polysomes while monosomes were mostly inactive. Hence, it is plausible that cells compensate the reduced availability of essential factors during stress by adapting the levels of active monosomes to favor sufficient elongation. Discussion: These results suggest that monosome and polysome levels are balanced under stress conditions. Together, our data argue for the existence of translational plasticity that ensure sufficient protein synthesis under stress conditions, a process that is necessary for cell survival and recovery.

Synonymous and Nonsynonymous Substitutions in Dictyostelium discoideum Ammonium Transporter amtA Are Necessary for Functional Complementation in Saccharomyces cerevisiae

Ammonium transporters are present in all three domains of life. They have undergone extensive horizontal gene transfer (HGT), gene duplication, and functional diversification and therefore offer an excellent paradigm to study protein evolution. We attempted to complement a mep1Δmep2Δmep3Δ strain of Saccharomyces cerevisiae (triple-deletion strain), which otherwise cannot grow on ammonium as a sole nitrogen source at concentrations of <3 mM, with amtA of Dictyostelium discoideum, an orthologue of S. cerevisiae MEP2. We observed that amtA did not complement the triple-deletion strain of S. cerevisiae for growth on low-ammonium medium. We isolated two mutant derivatives of amtA (amtA M1 and amtA M2) from a PCR-generated mutant plasmid library that complemented the triple-deletion strain of S. cerevisiae. amtA M1 bears three nonsynonymous and two synonymous substitutions, which are necessary for its functionality. amtA M2 bears two nonsynonymous substitutions and one synonymous substitution, all of which are necessary for functionality. Interestingly, AmtA M1 transports ammonium but does not confer methylamine toxicity, while AmtA M2 transports ammonium and confers methylamine toxicity, demonstrating functional diversification. Preliminary biochemical analyses indicated that the mutants differ in their conformations as well as their mechanisms of ammonium transport. These intriguing results clearly point out that protein evolution cannot be fathomed by studying nonsynonymous and synonymous substitutions in isolation. The above-described observations have significant implications for various facets of biological processes and are discussed in detail. IMPORTANCE Functional diversification following gene duplication is one of the major driving forces of protein evolution. While the role of nonsynonymous substitutions in the functional diversification of proteins is well recognized, knowledge of the role of synonymous substitutions in protein evolution is in its infancy. Using functional complementation, we isolated two functional alleles of the D. discoideum ammonium transporter gene (amtA), which otherwise does not function in S. cerevisiae as an ammonium transporters. One of them is an ammonium transporter, while the other is an ammonium transporter that also confers methylammonium (ammonium analogue) toxicity, suggesting functional diversification. Surprisingly, both alleles require a combination of synonymous and nonsynonymous substitutions for their functionality. These results bring out a hitherto-unknown pathway of protein evolution and pave the way for not only understanding protein evolution but also interpreting single nucleotide polymorphisms (SNPs).

Features of CFTR mRNA and implications for therapeutics development

Cystic fibrosis (CF) is an autosomal recessive disease impacting ∼100,000 people worldwide. This lethal disorder is caused by mutation of the CF transmembrane conductance regulator (CFTR) gene, which encodes an ATP-binding cassette-class C protein. More than 2,100 variants have been identified throughout the length of CFTR. These defects confer differing levels of severity in mRNA and/or protein synthesis, folding, gating, and turnover. Drug discovery efforts have resulted in recent development of modulator therapies that improve clinical outcomes for people living with CF. However, a significant portion of the CF population has demonstrated either no response and/or adverse reactions to small molecules. Additional therapeutic options are needed to restore underlying genetic defects for all patients, particularly individuals carrying rare or refractory CFTR variants. Concerted focus has been placed on rescuing variants that encode truncated CFTR protein, which also harbor abnormalities in mRNA synthesis and stability. The current mini-review provides an overview of CFTR mRNA features known to elicit functional consequences on final protein conformation and function, including considerations for RNA-directed therapies under investigation. Alternative exon usage in the 5'-untranslated region, polypyrimidine tracts, and other sequence elements that influence splicing are discussed. Additionally, we describe mechanisms of CFTR mRNA decay and post-transcriptional regulation mediated through interactions with the 3'-untranslated region (e.g. poly-uracil sequences, microRNAs). Contributions of synonymous single nucleotide polymorphisms to CFTR transcript utilization are also examined. Comprehensive understanding of CFTR RNA biology will be imperative for optimizing future therapeutic endeavors intended to address presently untreatable forms of CF.

Development and characterization of a glycine biosensor system for fine-tuned metabolic regulation in Escherichia coli

Background

In vivo biosensors have a wide range of applications, ranging from the detection of metabolites to the regulation of metabolic networks, providing versatile tools for synthetic biology and metabolic engineering. However, in view of the vast array of metabolite molecules, the existing number and performance of biosensors is far from sufficient, limiting their potential applications in metabolic engineering. Therefore, we developed the synthetic glycine-ON and -OFF riboswitches for metabolic regulation and directed evolution of enzyme in Escherichia coli.

Results

The results showed that a synthetic glycine-OFF riboswitch (glyOFF6) and an increased-detection-range synthetic glycine-ON riboswitch (glyON14) were successfully screened from a library based on the Bacillus subtilis glycine riboswitch using fluorescence-activated cell sorting (FACS) and tetA-based dual genetic selection. The two synthetic glycine riboswitches were successfully used in tunable regulation of lactate synthesis, dynamic regulation of serine synthesis and directed evolution of alanine-glyoxylate aminotransferase in Escherichia coli, respectively. Mutants AGXT22 and AGXT26 of alanine-glyoxylate aminotransferase with an increase of 58% and 73% enzyme activity were obtained by using a high-throughput screening platform based on the synthetic glycine-OFF riboswitch, and successfully used to increase the 5-aminolevulinic acid yield of engineered Escherichia coli.

Conclusions

A synthetic glycine-OFF riboswitch and an increased-detection-range synthetic glycine-ON riboswitch were successfully designed and screened. The developed riboswitches showed broad application in tunable regulation, dynamic regulation and directed evolution of enzyme in E. coli.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12934-022-01779-4.

Synonymous mutations that regulate translation speed might play a non-negligible role in liver cancer development

Background

Synonymous mutations do not change the protein sequences. Automatically, they have been regarded as neutral events and are ignored in the mutation-based cancer studies. However, synonymous mutations will change the codon optimality, resulting in altered translational velocity.

Methods

We fully utilized the transcriptome and translatome of liver cancer and normal tissue from ten patients. We profiled the mutation spectrum and examined the effect of synonymous mutations on translational velocity.

Results

Synonymous mutations that increase the codon optimality significantly enhanced the translational velocity, and were enriched in oncogenes. Meanwhile, synonymous mutations decreasing codon optimality slowed down translation, and were enriched in tumor suppressor genes. These synonymous mutations significantly contributed to the translational changes in tumor samples compared to normal samples.

Conclusions

Synonymous mutations might play a role in liver cancer development by altering codon optimality and translational velocity. Synonymous mutations should no longer be ignored in the genome-wide studies.

The Effect of Synonymous Single-Nucleotide Polymorphisms on an Atypical Cystic Fibrosis Clinical Presentation

Synonymous single nucleotide polymorphisms (sSNPs), which change a nucleotide, but not the encoded amino acid, are perceived as neutral to protein function and thus, classified as benign. We report a patient who was diagnosed with cystic fibrosis (CF) at an advanced age and presented very mild CF symptoms. The sequencing of the whole cystic fibrosis transmembrane conductance regulator (CFTR) gene locus revealed that the patient lacks known CF-causing mutations. We found a homozygous sSNP (c.1584G>A) at the end of exon 11 in the CFTR gene. Using sensitive molecular methods, we report that the c.1584G>A sSNP causes cognate exon skipping and retention of a sequence from the downstream intron, both of which, however, occur at a relatively low frequency. In addition, we found two other sSNPs (c.2562T>G (p.Thr854=) and c.4389G>A (p.Gln1463=)), for which the patient is also homozygous. These two sSNPs stabilize the CFTR protein expression, compensating, at least in part, for the c.1584G>A-triggered inefficient splicing. Our data highlight the importance of considering sSNPs when assessing the effect(s) of complex CFTR alleles. sSNPs may epistatically modulate mRNA and protein expression levels and consequently influence disease phenotype and progression.