Identification of the synonymous variant c.3141G > A in TNRC6B gene that altered RNA splicing by minigene assay.

BACKGROUND: Global developmental delay with speech and behavioral abnormalities (OMIM: 619243) is an autosomal dominant disease caused by variants in TNRC6B gene. METHOD: We reviewed and summarized clinical manifestations and genotypes in patients previously reported with TNRC6B gene variants. We used several prediction tools to predict pathogenicity and performed minigene assays to verify the function of the synonymous variant affecting RNA splicing. RESULT: The patient presented with convulsive seizures and developmental delay. WES combined with functional studies diagnosed a child with a synonymous variant in TNRC6B gene. Through minigene assay and Sanger sequencing, we demonstrated that c.3141G > A variant induced exon 7 skipping and the synonymous variant was pathogenic. CONCLUSION: Synonymous variants do not change the amino acids encoded by the codon, so we usually consider synonymous variants to be benign and ignore their pathogenicity. Minigene assay is a valuable tool to identify the effect of variation on RNA splicing and identify synonymous variants' benign or pathogenic. We showed that the synonymous variant was pathogenic by minigene assay. WES combined with minigene assay establishes a robust basis for genetic counseling and diagnosing diseases.

Evolutionary dependency of cancer mutations in gene pairs inferred by nonsynonymous-synonymous mutation ratios.

BACKGROUND: Determining the impact of somatic mutations requires understanding the functional relationship of genes acquiring mutations; however, it is largely unknown how mutations in functionally related genes influence each other. METHODS: We employed non-synonymous-to-synonymous or dNdS ratios to evaluate the evolutionary dependency (ED) of gene pairs, assuming a mutation in one gene of a gene pair can affect the evolutionary fitness of mutations in its partner genes as mutation context. We employed PanCancer- and tumor type-specific mutational profiles to infer the ED of gene pairs and evaluated their biological relevance with respect to gene dependency and drug sensitivity. RESULTS: We propose that dNdS ratios of gene pairs and their derived cdNS (context-dependent dNdS) scores as measure of ED distinguishing gene pairs either as synergistic (SYN) or antagonistic (ANT). Mutation contexts can induce substantial changes in the evolutionary fitness of mutations in the paired genes, e.g., IDH1 and IDH2 mutation contexts lead to substantial increase and decrease of dNdS ratios of ATRX indels and IDH1 missense mutations corresponding to SYN and ANT relationship with positive and negative cdNS scores, respectively. The impact of gene silencing or knock-outs on cell viability (genetic dependencies) often depends on ED, suggesting that ED can guide the selection of candidates for synthetic lethality such as TCF7L2-KRAS mutations. Using cell line-based drug sensitivity data, the effects of targeted agents on cell lines are often associated with mutations of genes exhibiting ED with the target genes, informing drug sensitizing or resistant mutations for targeted inhibitors, e.g., PRSS1 and CTCF mutations as resistant mutations to EGFR and BRAF inhibitors for lung adenocarcinomas and melanomas, respectively. CONCLUSIONS: We propose that the ED of gene pairs evaluated by dNdS ratios can advance our understanding of the functional relationship of genes with potential biological and clinical implications.

Synonymous variant at the terminal nucleotide in exon 3 of F7 causes abnormal splicing: A case report.

BACKGROUND: Synonymous variants are non-pathogenic due to non-substitution of amino acids. However, synonymous exonic terminal nucleotide substitutions may affect splicing. Splicing variants are easily analyzed at RNA level for genes expressed in blood cells. Minigene analysis provides another method for splicing variant analysis of genes that are poorly or not expressed in peripheral blood. METHODS: Whole exome sequencing was performed to screen for potential pathogenic mutations in the proband, which were validated within the family by Sanger sequencing. The pathogenicity of the synonymous mutation was analyzed using the minigene technology. RESULTS: The proband harbored the compound heterogeneous variants c. [291G >A; 572-50C >T] and c.681 + 1G >T in F7, of which the synonymous variant c.291G >A was located at the terminal position of exon 3. Minigene analysis revealed exon3 skipping due to this mutation, which may have subsequently affected protein sequence, structure, and function. CONCLUSION: Our finding confirmed the pathogenicity of c.291G >A, thus extending the pathogenic mutation spectrum of F7, and providing insights for effective reproductive counseling.

Synonymous and non-synonymous codon substitutions can alleviate dependence on GroEL for folding.

The Escherichia coli GroEL/ES chaperonin system facilitates protein folding in an ATP-driven manner. There are <100 obligate clients of this system in E. coli although GroEL can interact and assist the folding of a multitude of proteins in vitro. It has remained unclear, however, which features distinguish obligate clients from all the other proteins in an E. coli cell. To address this question, we established a system for selecting mutations in mouse dihydrofolate reductase (mDHFR), a GroEL interactor, that diminish its dependence on GroEL for folding. Strikingly, both synonymous and non-synonymous codon substitutions were found to reduce mDHFR's dependence on GroEL. The non-synonymous substitutions increase the rate of spontaneous folding whereas computational analysis indicates that the synonymous substitutions appear to affect translation rates at specific sites.

Novel Synonymous Variant in IL7R Causes Preferential Expression of the Soluble Isoform.

PURPOSE: The interleukin-7 receptor (IL-7R) is primarily expressed on lymphoid cells and plays a crucial role in the development, proliferation, and survival of T cells. Autosomal recessive mutations that disrupt IL-7Rα chain expression give rise to a severe combined immunodeficiency (SCID), which is characterized by lymphopenia and a T-B+NK+ phenotype. The objective here was to diagnose two siblings displaying the T-B+NK+ SCID phenotype as initial clinical genetic testing did not detect any variants in known SCID genes. METHODS: Whole genome sequencing (WGS) was utilized to identify potential variants causing the SCID phenotype. Splicing prediction tools were employed to assess the deleterious impact of the mutation. Polymerase Chain Reaction (PCR), Sanger sequencing, flow cytometry, and ELISA were then used to validate the pathogenicity of the detected mutation. RESULTS: We discovered a novel homozygous synonymous mutation in the IL7R gene. Our functional studies indicate that this variant is pathogenic, causing exon 6, which encodes the transmembrane domain, to be preferentially spliced out. CONCLUSION: In this study, we identified a novel rare synonymous mutation causing a loss of IL-7Rα expression at the cellular membrane. This case demonstrates the value of reanalyzing genetic data based on the clinical phenotype and highlights the significance of functional studies in determining the pathogenicity of genetic variants.

Phenotypic Description of A Patient with ODLURO Syndrome and Functional Characterization of the Pathogenetic Role of A Synonymous Variant c.186G>A in KMT2E Gene.

O'Donnell-Luria-Rodan (ODLURO) syndrome is an autosomal dominant disorder caused by mutations in the KMT2E gene. The clinical phonotype of the affected individuals is typically characterized by global developmental delay, autism, epilepsy, hypotonia, macrocephaly, and very mild dysmorphic facial features. In this report, we describe the case of a 6-year-old boy with ODLURO syndrome who is a carrier of the synonymous mutation c.186G>A (p.Ala62=) in the KMT2E gene, predicted to alter splicing by in silico tools. Given the lack of functional studies on the c.186G>A variant, in order to assess its potential functional effect, we sequenced the patient's cDNA demonstrating its impact on the mechanism of splicing. To the best of our knowledge, our patient is the second to date reported carrying this synonymous mutation, but he is the first whose functional investigation has confirmed the deleterious consequence of the variant, resulting in exon 4 skipping. Additionally, we suggest a potential etiological mechanism that could be responsible for the aberrant splicing mechanism in KMT2E.

Re-examining Correlations Between Synonymous Codon Usage and Protein Bond Angles in Escherichia coli.

Rosenberg AA, Marx A, Bronstein AM (Codon-specific Ramachandran plots show amino acid backbone conformation depends on identity of the translated codon. Nat Commun. 2022:13:2815) recently found a surprising correlation between synonymous codon usage and the dihedral bond angles of the resulting amino acid. However, their analysis did not account for the strongest known correlate of codon usage: gene expression. We re-examined the relationship between bond angles and codon usage by applying the approach of Rosenberg et al. to simulated protein-coding sequences that (i) have random codon usage, (ii) codon usage determined by mutation biases, and (iii) maintain the general relationship between codon usage and gene expression via the assumption of selection-mutation-drift equilibrium. We observed correlations between dihedral bond angle and codon usage when codon usage is entirely random, indicating possible conflation of noise with differences in bond angle distributions between synonymous codons. More relevant to the general analysis of codon usage patterns, we found surprisingly good agreement between the analysis of the real sequences and the analysis of sequences simulated assuming selection-mutation-drift equilibrium, with 91% of significant synonymous codon pairs detected in the former were also detected in the latter. We believe the correlation between codon usage and dihedral bond angles resulted from the variation in codon usage across genes due to the interplay between mutation bias, natural selection for translation efficiency, and gene expression, further underscoring these factors must be controlled for when looking for novel patterns related to codon usage.

Fitness difference between two synonymous mutations of Phytophthora infestans ATP6 gene.

BACKGROUND: Sequence variation produced by mutation provides the ultimate source of natural selection for species adaptation. Unlike nonsynonymous mutation, synonymous mutations are generally considered to be selectively neutral but accumulating evidence suggests they also contribute to species adaptation by regulating the flow of genetic information and the development of functional traits. In this study, we analysed sequence characteristics of ATP6, a housekeeping gene from 139 Phytophthora infestans isolates, and compared the fitness components including metabolic rate, temperature sensitivity, aggressiveness, and fungicide tolerance among synonymous mutations. RESULTS: We found that the housekeeping gene exhibited low genetic variation and was represented by two major synonymous mutants at similar frequency (0.496 and 0.468, respectively). The two synonymous mutants were generated by a single nucleotide substitution but differed significantly in fitness as well as temperature-mediated spatial distribution and expression. The synonymous mutant ending in AT was more common in cold regions and was more expressed at lower experimental temperature than the synonymous mutant ending in GC and vice versa. CONCLUSION: Our results are consistent with the argument that synonymous mutations can modulate the adaptive evolution of species including pathogens and have important implications for sustainable disease management, especially under climate change.

Synonymous Mutations Can Alter Protein Dimerization Through Localized Interface Misfolding Involving Self-entanglements.

Synonymous mutations in messenger RNAs (mRNAs) can reduce protein-protein binding substantially without changing the protein's amino acid sequence. Here, we use coarse-grain simulations of protein synthesis, post-translational dynamics, and dimerization to understand how synonymous mutations can influence the dimerization of two E. coli homodimers, oligoribonuclease and ribonuclease T. We synthesize each protein from its wildtype, fastest- and slowest-translating synonymous mRNAs in silico and calculate the ensemble-averaged interaction energy between the resulting dimers. We find synonymous mutations alter oligoribonuclease's dimer properties. Relative to wildtype, the dimer interaction energy becomes 4% and 10% stronger, respectively, when translated from its fastest- and slowest-translating mRNAs. Ribonuclease T dimerization, however, is insensitive to synonymous mutations. The structural and kinetic origin of these changes are misfolded states containing non-covalent lasso-entanglements, many of which structurally perturb the dimer interface, and whose probability of occurrence depends on translation speed. These entangled states are kinetic traps that persist for long time scales. Entanglements cause altered dimerization energies for oligoribonuclease, as there is a large association (odds ratio: 52) between the co-occurrence of non-native self-entanglements and weak-binding dimer conformations. Simulated at all-atom resolution, these entangled structures persist for long timescales, indicating the conclusions are independent of model resolution. Finally, we show that regions of the protein we predict to have changes in entanglement are also structurally perturbed during refolding, as detected by limited-proteolysis mass spectrometry. Thus, non-native changes in entanglement at dimer interfaces is a mechanism through which oligomer structure and stability can be altered.

Quantifying negative selection on synonymous variants.

Widespread adoption of DNA sequencing has resulted in large numbers of genetic variants, whose contribution to disease is not easily determined. Although many types of variation are known to disrupt cellular processes in predictable ways, for some categories of variants, the effects may not be directly detectable. A particular example is synonymous variants, that is, those single-nucleotide variants that create a codon substitution, such that the produced amino acid sequence is unaffected. Contrary to the original theory suggesting that synonymous variants are benign, there is a growing volume of research showing that, despite their "silent" mechanism of action, some synonymous variation may be deleterious. Here, we studied the extent of the negative selective pressure acting on different classes of synonymous variants by analyzing the relative enrichment of synonymous singleton variants in the human exomes provided by gnomAD. Using a modification of the mutability-adjusted proportion of singletons (MAPS) metric as a measure of purifying selection, we found that some classes of synonymous variants are subject to stronger negative selection than others. For instance, variants that reduce codon optimality undergo stronger selection than optimality-increasing variants. Besides, selection affects synonymous variants implicated in splice-site-loss or splice-site-gain events. To understand what drives this negative selection, we tested a number of predictors in the aim to explain the variability in the selection scores. Our findings provide insights into the effects of synonymous variants at the population level, highlighting the specifics of the role that these variants play in health and disease.

Selection on synonymous sites: the unwanted transcript hypothesis.

Although translational selection to favour codons that match the most abundant tRNAs is not readily observed in humans, there is nonetheless selection in humans on synonymous mutations. We hypothesize that much of this synonymous site selection can be explained in terms of protection against unwanted RNAs - spurious transcripts, mis-spliced forms or RNAs derived from transposable elements or viruses. We propose not only that selection on synonymous sites functions to reduce the rate of creation of unwanted transcripts (for example, through selection on exonic splice enhancers and cryptic splice sites) but also that high-GC content (but low-CpG content), together with intron presence and position, is both particular to functional native mRNAs and used to recognize transcripts as native. In support of this hypothesis, transcription, nuclear export, liquid phase condensation and RNA degradation have all recently been shown to promote GC-rich transcripts and suppress AU/CpG-rich ones. With such 'traps' being set against AU/CpG-rich transcripts, the codon usage of native genes has, in turn, evolved to avoid such suppression. That parallel filters against AU/CpG-rich transcripts also affect the endosomal import of RNAs further supports the unwanted transcript hypothesis of synonymous site selection and explains the similar design rules that have enabled the successful use of transgenes and RNA vaccines.

Molecular Mechanisms and the Significance of Synonymous Mutations.

Synonymous mutations result from the degeneracy of the genetic code. Most amino acids are encoded by two or more codons, and mutations that change a codon to another synonymous codon do not change the amino acid in the gene product. Historically, such mutations have been considered silent because they were assumed to have no to very little impact. However, research in the last few decades has produced several examples where synonymous mutations play important roles. These include optimizing expression by enhancing translation initiation and accelerating or decelerating translation elongation via codon usage and mRNA secondary structures, stabilizing mRNA molecules and preventing their breakdown before translation, and faulty protein folding or increased degradation due to enhanced ubiquitination and suboptimal secretion of proteins into the appropriate cell compartments. Some consequences of synonymous mutations, such as mRNA stability, can lead to different outcomes in prokaryotes and eukaryotes. Despite these examples, the significance of synonymous mutations in evolution and in causing disease in comparison to nonsynonymous mutations that do change amino acid residues in proteins remains controversial. Whether the molecular mechanisms described by which synonymous mutations affect organisms can be generalized remains poorly understood and warrants future research in this area.

When a synonymous mutation breaks the silence in a thalassaemia patient.

Synonymous mutations were considered to lack functional roles in human diseases; however, distinguishing deleterious synonymous mutations from benign ones is still a challenge. In this article, we identified a deleterious synonymous mutation ß-codon 16 (C>T). HBB: c.51C>T, in compound heterozygous form with known ß-thalassaemia mutation patients who clinically presented as non-transfusion-dependent thalassaemia (NTDT). A total of 9 families with 11 compound heterozygous index cases were reported. In the heterozygous state, codon 16 (C>T) mutation results in borderline HbA2 (3.18 ± 0.36%) and slightly reduced RBC indices (RBCs: 4.73 ± 0.75 × 106 /µL, Hb: 12.26 ± 2.60 g/dL, MCV: 79.48 ± 8.40 fL, MCH: 25.95 ± 4.15 pg). The compound heterozygous patients showed elevated HbA2 (5.98 ± 1.17%) and HbF (12.75 ± 7.51%) and presented clinically as NTDT with a mean Hb of 6.95 ± 1.29 g/dL. Many of them were dependent on few transfusions and had mild splenomegaly. Of the 11 patients, 5 (45.4%) were treated with hydroxyurea. This study highlights the clinical significance of synonymous mutation, when inherited with other ß-thalassaemia mutations leading to the phenotype of NTDT. Thus, the study would help to improve screening protocols for ß-thalassaemia carriers, which will ultimately improve the prevention programme.

Effect Predictor of Driver Synonymous Mutations Based on Multi-Feature Fusion and Iterative Feature Representation Learning.

Accurate identification of driver mutations is crucial in genetic studies of human cancers. While numerous cancer driver missense mutations have been identified, research into potential cancer drivers for synonymous mutations has shown limited success to date. Here, we developed a novel machine learning framework, epSMic, for predicting cancer driver synonymous mutations. epSMic employs an iterative feature representation scheme that facilitates the learning of discriminative features from various sequential models in a supervised iterative mode. We constructed the benchmark datasets and encoded the embedding sequence, physicochemical property, and basic information such as conservation and splicing feature. The evaluation results on benchmark test datasets demonstrate that epSMic outperforms existing methods, making it a valuable tool for researchers in identifying functional synonymous mutations in cancer. We hope epSMic can enable researchers to concentrate on synonymous mutations that have a functional impact on cancer.

A synonymous variant is unmasked in thalassaemia.

The genetic underpinnings of beta-thalassaemia encompass a myriad of molecular mechanisms. The ability of synonymous mutations, an often-overlooked category of variants, to influence ß-globin expression and phenotypic disease is highlighted by this report by Gorivale et al. Commentary on: Gorivale et al. When a synonymous mutation breaks the silence in a thalassaemia patient. Br J Haematol 2024;204:677-682.

Molecular bases for strong phenotypic effects of single synonymous codon substitutions in the E. coli ccdB toxin gene.

BACKGROUND: Single synonymous codon mutations typically have only minor or no effects on gene function. Here, we estimate the effects on cell growth of ~ 200 single synonymous codon mutations in an operonic context by mutating almost all positions of ccdB, the 101-residue long cytotoxin of the ccdAB Toxin-Antitoxin (TA) operon to most degenerate codons. Phenotypes were assayed by transforming the mutant library into CcdB sensitive and resistant E. coli strains, isolating plasmid pools, and subjecting them to deep sequencing. Since autoregulation is a hallmark of TA operons, phenotypes obtained for ccdB synonymous mutants after transformation in a RelE toxin reporter strain followed by deep sequencing provided information on the amount of CcdAB complex formed. RESULTS: Synonymous mutations in the N-terminal region involved in translation initiation showed the strongest non-neutral phenotypic effects. We observe an interplay of numerous factors, namely, location of the codon, codon usage, t-RNA abundance, formation of anti-Shine Dalgarno sequences, predicted transcript secondary structure, and evolutionary conservation in determining phenotypic effects of ccdB synonymous mutations. Incorporation of an N-terminal, hyperactive synonymous mutation, in the background of the single synonymous codon mutant library sufficiently increased translation initiation, such that mutational effects on either folding or termination of translation became more apparent. Introduction of putative pause sites not only affects the translational rate, but might also alter the folding kinetics of the protein in vivo. CONCLUSION: In summary, the study provides novel insights into diverse mechanisms by which synonymous mutations modulate gene function. This information is useful in optimizing heterologous gene expression in E. coli and understanding the molecular bases for alteration in gene expression that arise due to synonymous mutations.

Incorporating mutational heterogeneity to identify genes that are enriched for synonymous mutations in cancer.

BACKGROUND: Synonymous mutations, which change the DNA sequence but not the encoded protein sequence, can affect protein structure and function, mRNA maturation, and mRNA half-lives. The possibility that synonymous mutations might be enriched in cancer has been explored in several recent studies. However, none of these studies control for all three types of mutational heterogeneity (patient, histology, and gene) that are known to affect the accurate identification of non-synonymous cancer-associated genes. Our goal is to adopt the current standard for non-synonymous mutations in an investigation of synonymous mutations. RESULTS: Here, we create an algorithm, MutSigCVsyn, an adaptation of MutSigCV, to identify cancer-associated genes that are enriched for synonymous mutations based on a non-coding background model that takes into account the mutational heterogeneity across these levels. Using MutSigCVsyn, we first analyzed 2572 cancer whole-genome samples from the Pan-cancer Analysis of Whole Genomes (PCAWG) to identify non-synonymous cancer drivers as a quality control. Indicative of the algorithm accuracy we find that 58.6% of these candidate genes were also found in Cancer Census Gene (CGC) list, and 66.2% were found within the PCAWG cancer driver list. We then applied it to identify 30 putative cancer-associated genes that are enriched for synonymous mutations within the same samples. One of the promising gene candidates is the B cell lymphoma 2 (BCL-2) gene. BCL-2 regulates apoptosis by antagonizing the action of proapoptotic BCL-2 family member proteins. The synonymous mutations in BCL2 are enriched in its anti-apoptotic domain and likely play a role in cancer cell proliferation. CONCLUSION: Our study introduces MutSigCVsyn, an algorithm that accounts for mutational heterogeneity at patient, histology, and gene levels, to identify cancer-associated genes that are enriched for synonymous mutations using whole genome sequencing data. We identified 30 putative candidate genes that will benefit from future experimental studies on the role of synonymous mutations in cancer biology.

Investigation of a synonymous mutation in Btk in a patient with agammaglobulinemia: A case report.

BACKGROUND: X-linked agammaglobulinemia (XLA) is the most common form of agammaglobulinemia and is caused by mutations in Btk, which encodes Bruton tyrosine kinase (BTK). CASE DESCRIPTION: We describe a 36-year-old male who presented as an infant with hypogammaglobulinemia and sinopulmonary infections and was initially diagnosed with common variable immunodeficiency. Genetic testing showed he was hemizygous for Btk c.240G > A. This synonymous variant affecting the last nucleotide of exon 3 leads to aberrant splicing of most but not all mRNA transcripts. CONCLUSION: We demonstrated reduced BTK protein expression confirming the pathogenicity of the variant and related our findings to genotype-phenotype relationship studies ina XLA caused by synonymous mutations.

Gene recoding by synonymous mutations creates promiscuous intragenic transcription initiation in mycobacteria.

IMPORTANCE: Mycobacterium tuberculosis (Mtb) is the causative agent of tuberculosis, one of the deadliest infectious diseases worldwide. Previous studies have established that synonymous recoding to introduce rare codon pairings can attenuate viral pathogens. We hypothesized that non-optimal codon pairing could be an effective strategy for attenuating gene expression to create a live vaccine for Mtb. We instead discovered that these synonymous changes enabled the transcription of functional mRNA that initiated in the middle of the open reading frame and from which many smaller protein products were expressed. To our knowledge, this is one of the first reports that synonymous recoding of a gene in any organism can create or induce intragenic transcription start sites.

Synonymous Variants of Uncertain Silence.

Synonymous variants, traditionally regarded as silent mutations due to their lack of impact on protein sequence, structure and function, have been the subject of increasing scrutiny. This commentary explores the emerging evidence challenging the notion of synonymous variants as functionally inert. Analysis of the activity of 70 synonymous variants in the HIV Tat transcription factor revealed that 50% of the variants exhibited significant deviations from wild-type activity. Our analysis supports previous work and raises important questions about the broader impact of non-silent synonymous variants in human genes. Considering the potential functional implications, the authors propose classifying such variants as "synonymous variants of uncertain silence" (sVUS), highlighting the need for cautious interpretation and further investigations in clinical and genetic testing settings.