A synonymous single nucleotide polymorphism in DeltaF508 CFTR alters the secondary structure of the mRNA and the expression of the mutant protein

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.

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.

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 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.

Exploring the influence of DNA methylation and single nucleotide polymorphisms of the Myostatin gene on growth traits in the hybrid grouper (Epinephelus fuscoguttatus (female) × Epinephelus polyphekadion (male))

Investigations into the correlation between growth characteristics and DNA methylation levels, along with genetic variations, can provide fundamental insights to enhance growth performance in groupers. The Myostatin (mstn) gene plays a vital role in regulating skeletal muscle development and growth. This study scrutinized the DNA methylation levels of the mstn gene across hybrid groupers (E. fuscoguttatus (♀) × E. polyphekadion (♂)) and their parental species, to evaluate its impact on growth attributes in grouper fish. The nucleotide sequence of the mstn gene was directly sequenced in the hybrid grouper, exhibiting different growth performance to identify the single nucleotide polymorphisms (SNPs) of the mstn gene and explore their correlation with growth characteristics. The findings revealed no significant differences in global DNA methylation levels within muscle tissue among the hybrid grouper and parents. However, significant differences in DNA methylation sites were discovered between the hybrid grouper and E. polyphekadion at sites 824 and 1521 (located at exon 2 and intron 2, respectively), and between E. fuscoguttatus and E. polyphekadion at site 1521. These variations could potentially influence the mRNA expression of the mstn gene. The study also identified that SNP g.1003 T > C in exon 2 of the mstn gene was significantly associated with various growth traits including body weight, total length, body length, head length, caudal peduncle height, and body height (p < 0.01). Specimens with the TT genotype at site 1003 demonstrated superior growth performance compared to those with the TC genotype. Furthermore, microstructural analyses of muscle tissue showed that the average area and diameter of muscle fibers in TT genotype individuals were significantly greater than those in TC genotype individuals. Therefore, this research provides robust evidence linking the DNA methylation level and polymorphisms of the mstn gene with growth traits, which could be beneficial for grouper breeding programs.

Interferon-stimulated gene 15 polymorphisms are associated with spontaneous preterm birth in Taiwanese women

PROBLEM

Immune and inflammatory responses are known to be major causes of preterm birth (PTB). The maternal genetic background plays an important role in the development of PTB. Interferon-stimulated gene 15 (ISG15) is an interferon-induced protein which can modulate immune cell activation and function. We aim to study if polymorphisms in the ISG15 gene are associated with spontaneous PTB (sPTB) risk in Taiwanese women.

METHOD OF STUDY

ISG15 rs4615788 C/G, rs1921 G/A, and rs8997 A/G polymorphisms were genotyped in a hospital-based study of 112 women with sPTB and 1120 term controls. The plasma concentrations of ISG15 were determined by enzyme-linked immunosorbent assay.

RESULTS

We found the ISG15 rs1921 G-rs8997 A haplotype was associated with decreased risk for PTB (χ2  = 6.26, p = .01, pc  = .04). The A/G genotype of ISG15 rs8997 polymorphism might have the potential to confer reduced risk of PTB women (χ2  = 4.09, p = .04, pc  = .08). Spontaneous PTB women displayed higher plasma ISG15 levels compared to term controls (p < .001). The plasma ISG15 levels among pregnant women with rs8997 A/G genotype were found significantly lower compared to G/G genotype (p = .03).

CONCLUSIONS

Women with the ISG15 rs1921 G-rs8997 A haplotype may associate with spontaneous PTB. These findings provide new insights into the etiology of preterm birth.

Polymorphism and expression level of the FADS3 gene and associated with the growth traits in Hu sheep

Growth traits are the economically important traits of sheep, and screening for genes related to growth and development is helpful for the genetic improvement of ovine growth traits. The fatty acid desaturase 3 (FADS3) is one of the important genes affecting the synthesis and accumulation of polyunsaturated fatty acids in animals. In this study, the expression levels of the FADS3 gene and polymorphism of the FADS3 gene associated with growth traits in Hu sheep were detected using quantitative real-time PCR (qRT-PCR), Sanger sequencing, and KAspar assay. The result showed that the expression levels of the FADS3 gene were widely expressed in all tissues, and the expression level of FADS3 in the lung was significantly higher than in other tissues (p < .05). Then, the polymorphism locus g. 2918 A > C was detected in intron 2 of the FADS3 gene, and associated analysis showed that the mutation in the FADS3 gene was associated significantly with growth traits (including body weight, body height, body length, and chest circumference, p < .05). Therefore, individuals with AA genotype showed significantly better growth traits than those with CC genotype, and FADS3 gene could be a candidate gene for improving growth traits in Hu sheep.

In silico methods for predicting functional synonymous variants

Single nucleotide variants (SNVs) contribute to human genomic diversity. Synonymous SNVs are previously considered to be "silent," but mounting evidence has revealed that these variants can cause RNA and protein changes and are implicated in over 85 human diseases and cancers. Recent improvements in computational platforms have led to the development of numerous machine-learning tools, which can be used to advance synonymous SNV research. In this review, we discuss tools that should be used to investigate synonymous variants. We provide supportive examples from seminal studies that demonstrate how these tools have driven new discoveries of functional synonymous SNVs.

Identification and in silico characterization of CSRP3 synonymous variants in dilated cardiomyopathy

BACKGROUND

Synonymous variations have always been ignored while studying the underlying genetic mechanisms for most of the human diseases. However, recent studies have suggested that these silent changes in the genome can alter the protein expression and folding.

METHODS AND RESULTS

CSRP3, which is a well-known candidate gene associated with dilated cardiomyopathy (DCM) and hypertrophic cardiomyopathy (HCM), was screened for 100 idiopathic DCM cases and 100 controls. Three synonymous variations were identified viz., c.96G > A, p.K32=; c.336G > A, p.A112=; c.354G > A, p.E118=. A comprehensive in silico analysis was performed using various web based widely accepted tools, Mfold, Codon Usage, HSF3.1 and RNA22. Mfold predicted structural changes in all the variants except c.96 G > A (p.K32=), however it predicted changes in the stability of mRNA due to all the synonymous variants. Codon bias was observed as evident by the Relative Synonymous Codon Usage and Log Ratio of Codon Usage Frequencies. The Human Splicing Finder also predicted remarkable changes in the regulatory elements in the variants c.336G > A and c.354 G > A. The miRNA target prediction using varied modes available in RNA22 revealed that 70.6% of the target sites of miRNAs in CSRP3 were altered due to variant c.336G > A while 29.41% sites were completely lost.

CONCLUSION

Findings of the present study suggest that synonymous variants revealed striking deviations in the structural conformation of mRNA, stability of mRNA, relative synonymous codon usage, splicing and miRNA binding sites from the wild type suggesting their possible role in the pathogenesis of DCM, either by destabilizing the mRNA structure, or codon usage bias or else altering the cis-acting regulatory elements during splicing.

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).

Implementing computational methods in tandem with synonymous gene recoding for therapeutic development

Synonymous gene recoding, the substitution of synonymous variants into the genetic sequence, has been used to overcome many production limitations in therapeutic development. However, the safety and efficacy of recoded therapeutics can be difficult to evaluate because synonymous codon substitutions can result in subtle, yet impactful changes in protein features and require sensitive methods for detection. Given that computational approaches have made significant leaps in recent years, we propose that machine-learning (ML) tools may be leveraged to assess gene-recoded therapeutics and foresee an opportunity to adapt codon contexts to enhance some powerful existing tools. Here, we examine how synonymous gene recoding has been used to address challenges in therapeutic development, explain the biological mechanisms underlying its effects, and explore the application of computational platforms to improve the surveillance of functional variants in therapeutic design.

Synonymous mutations in the phosphoglycerate kinase 1 gene induce an altered response to protein misfolding in Schizosaccharomyces pombe

Background

Proteostasis refers to the processes that regulate the biogenesis, folding, trafficking, and degradation of proteins. Any alteration in these processes can lead to cell malfunction. Protein synthesis, a key proteostatic process, is highly-regulated at multiple levels to ensure adequate adaptation to environmental and physiological challenges such as different stressors, proteotoxic conditions and aging, among other factors. Because alterations in protein translation can lead to protein misfolding, examining how protein translation is regulated may also help to elucidate in part how proteostasis is controlled. Codon usage bias has been implicated in the fine-tuning of translation rate, as more-frequent codons might be read faster than their less-frequent counterparts. Thus, alterations in codon usage due to synonymous mutations may alter translation kinetics and thereby affect the folding of the nascent polypeptide, without altering its primary structure. To date, it has been difficult to predict the effect of synonymous mutations on protein folding and cellular fitness due to a scarcity of relevant data. Thus, the purpose of this work was to assess the effect of synonymous mutations in discrete regions of the gene that encodes the highly-expressed enzyme 3-phosphoglycerate kinase 1 (pgk1) in the fission yeast Schizosaccharomyces pombe.

Results

By means of systematic replacement of synonymous codons along pgk1, we found slightly-altered protein folding and activity in a region-specific manner. However, alterations in protein aggregation, heat stress as well as changes in proteasome activity occurred independently of the mutated region. Concomitantly, reduced mRNA levels of the chaperones Hsp9 and Hsp16 were observed.

Conclusion

Taken together, these data suggest that codon usage bias of the gene encoding this highly-expressed protein is an important regulator of protein function and proteostasis.

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.

A minimal role for synonymous variation in human disease

Synonymous mutations change the DNA sequence of a gene without affecting the amino acid sequence of the encoded protein. Although some synonymous mutations can affect RNA splicing, translational efficiency, and mRNA stability, studies in human genetics, mutagenesis screens, and other experiments and evolutionary analyses have repeatedly shown that most synonymous variants are neutral or only weakly deleterious, with some notable exceptions. Based on a recent study in yeast, there have been claims that synonymous mutations could be as important as nonsynonymous mutations in causing disease, assuming the yeast findings hold up and translate to humans. Here, we argue that there is insufficient evidence to overturn the large, coherent body of knowledge establishing the predominant neutrality of synonymous variants in the human genome.

Ovine RAP1GAP and rBAT gene polymorphisms and their association with tail fat deposition in Hu sheep

Excessive fat deposition in the tail of sheep will affect its feed efficiency, which will increase the feeding cost. The purpose of this study was to identify the single nucleotide polymorphisms (SNPs) of RAP1GAP and rBAT genes by PCR amplification and Sanger sequencing, the SNPs were genotyped by KASP genotyping assays to evaluate their association with tail fat deposition traits. The results showed that two intronic mutations of g.13561 G > A and g.1460 T > C were found in RAP1GAP and rBAT, respectively. There were three genotypes of GG, AG, AA and CC, CT and TT at these two loci, respectively. Association analysis showed that g.13561 G > A of RAP1GAP was associated with tail width, tail fat weight and relative tail fat weight (P < 0.05). The g.1460 T > C of rBAT was associated with tail width and tail fat weight (P < 0.05). Different combinations of genotypes also differed significantly with tail fat deposition traits. In the tail fat tissue, the expression levels of RAP1GAP gene was significantly higher in small-tailed sheep than in big-tailed sheep, and the expression levels of rBAT gene was significantly higher in big-tailed sheep than in small-tailed sheep. In the liver, the expression levels of RAP1GAP and rBAT gene was significantly higher at 6 months than at 0 and 3 months. In conclusion, RAP1GAP and rBAT polymorphisms can be used as a candidate molecular marker to reduce tail fat deposition in sheep.

Structural, functional, and immunogenicity implications of F9 gene recoding

Hemophilia B is a blood clotting disorder caused by deficient activity of coagulation factor IX (FIX). Multiple recombinant FIX proteins are currently approved to treat hemophilia B, and several gene therapy products are currently being developed. Codon optimization is a frequently used technique in the pharmaceutical industry to improve recombinant protein expression by recoding a coding sequence using multiple synonymous codon substitutions. The underlying assumption of this gene recoding is that synonymous substitutions do not alter protein characteristics because the primary sequence of the protein remains unchanged. However, a critical body of evidence shows that synonymous variants can affect cotranslational folding and protein function. Gene recoding could potentially alter the structure, function, and in vivo immunogenicity of recoded therapeutic proteins. Here, we evaluated multiple recoded variants of F9 designed to further explore the effects of codon usage bias on protein properties. The detailed evaluation of these constructs showed altered conformations, and assessment of translation kinetics by ribosome profiling revealed differences in local translation kinetics. Assessment of wild-type and recoded constructs using a major histocompatibility complex (MHC)-associated peptide proteomics assay showed distinct presentation of FIX-derived peptides bound to MHC class II molecules, suggesting that despite identical amino acid sequence, recoded proteins could exhibit different immunogenicity risks. Posttranslational modification analysis indicated that overexpression from gene recoding results in suboptimal posttranslational processing. Overall, our results highlight potential functional and immunogenicity concerns associated with gene-recoded F9 products. These findings have general applicability and implications for other gene-recoded recombinant proteins.

Expression Profiles of the Ovine IL18 Gene and Association of Its Polymorphism With Hematologic Parameters in Hu Lambs

Hematological traits are important indexes to evaluate health status and immunological conditions in human and livestock. In this study, we measured the hematologic indexes of 819 male Hu lambs and performed the descriptive statistical analysis. The results showed the coefficients of variation of partial indexes were >10%, and the heritability for mean erythrocyte volume (MCV), white blood cell count (WBC), hemoglobin concentration (HGB), hematocrit (HCT), and red blood cell (RBC) distribution-standard deviation (RDW_SD) were medium to high, ranging from 0.17 to 0.43. In addition, Interleukin 18 (IL18), as an important regulator of both innate and acquired immune responses, was selected as candidate gene and subjected to the expression profile analysis, single nucleotide polymorphism (SNP) scanning and association analysis by using quantitative real-time PCR (qRT-PCR), PCR amplification, Sanger sequencing, and KASP genotyping. The results of qRT-PCR indicated that IL18 is predominantly expressed in lymph and lung compared with that in the other tested tissues. In addition, three novel polymorphisms (g. 24991544 A > G, g. 24991651 A > G, and g. 24991749 C > T) were identified in IL18, and the three SNPs were in a strong linkage state. Therefore, only a SNP was genotyped and performed association analysis in the enlarged experimental population, the result of association analysis demonstrated that the polymorphism g. 24991651 A > G was significantly associated with RBC, MCV, MCHC, and RDW_CV. These results will provide the reference values and the novel genetic markers of hematological parameters in sheep.

Polymorphisms in common antihypertensive targets: Pharmacogenomic implications for the treatment of cardiovascular disease

The idea of personalized medicine came to fruition with sequencing the human genome; however, aside from a few cases, the genetic revolution has yet to materialize. Cardiovascular diseases are the leading cause of death globally, and hypertension is a common prelude to nearly all cardiovascular diseases. Thus, hypertension is an ideal candidate disease to apply tenants of personalized medicine to lessen cardiovascular disease. Herein is a survey that visually depicts the polymorphisms in the top eight antihypertensive targets. Although there are numerous genome-wide association studies regarding cardiovascular disease, few studies look at the effects of receptor polymorphisms on drug treatment. With 17,000+ polymorphisms in the combined target proteins examined, it is expected that some of the clinical variability in the treatment of hypertension is due to polymorphisms in the drug targets. Recent advances in techniques and technology, such as high throughput examination of single mutations, structure prediction, computational power for modeling, and CRISPR models of point mutations, allow for a relatively rapid and comprehensive examination of the effects of known and future polymorphisms on drug affinity and effects. As hypertension is easy to measure and has a plethora of clinically viable ligands, hypertension makes an excellent disease to study pharmacogenomics in the lab and the clinic. If the promises of personalized medicine are to materialize, a concerted effort to examine the effects polymorphisms have on drugs is required. A clinician with the knowledge of a patient's genotype can then prescribe drugs that are optimal for treating that specific patient.

Expression features of the ovine FTO gene and association between FTO polymorphism and tail fat deposition related-traits in Hu sheep

Fat mass and obesity associated (FTO) gene is a famous dominant predictor of obesity, and plays pivotal roles in the regulating fatty acid transport and fat metabolism. In the current study, the quantitative real-time PCR (qRT-PCR), DNA-pooled sequencing and KAspar assay were performed to detect the expression features of FTO and the polymorphisms of FTO associated with the tail fat weight related-traits of Hu sheep. The results indicated that the expression of FTO gene is widely expressed in the tissues tested, and the expression level of FTO in tail fat tissue was evidently higher compared with that in the other tissues. In addition, FTO showed the highest expression level in tail fat tissue at three months. The expression of FTO mRNA was lower in the large-tail fat group compared with that in the small-tail fat group in 6 months old. Subsequently, the polymorphism loci g. 23704451C > A detected in the FTO gene was confirmed to be significantly related to the tail length and the weight of tail fat. These results suggested that the polymorphism might be regarded as novel molecular marker for the breeding of small tail sheep.

Codon-specific Ramachandran plots show amino acid backbone conformation depends on identity of the translated codon

Synonymous codons translate into chemically identical amino acids. Once considered inconsequential to the formation of the protein product, there is evidence to suggest that codon usage affects co-translational protein folding and the final structure of the expressed protein. Here we develop a method for computing and comparing codon-specific Ramachandran plots and demonstrate that the backbone dihedral angle distributions of some synonymous codons are distinguishable with statistical significance for some secondary structures. This shows that there exists a dependence between codon identity and backbone torsion of the translated amino acid. Although these findings cannot pinpoint the causal direction of this dependence, we discuss the vast biological implications should coding be shown to directly shape protein conformation and demonstrate the usefulness of this method as a tool for probing associations between codon usage and protein structure. Finally, we urge for the inclusion of exact genetic information into structural databases.

Nanomechanics combined with HDX reveals allosteric drug binding sites of CFTR NBD1

Cystic fibrosis (CF) is a frequent genetic disease in Caucasians that is caused by the deletion of F508 (ΔF508) in the nucleotide binding domain 1 (NBD1) of the CF transmembrane conductance regulator (CFTR). The ΔF508 compromises the folding energetics of the NBD1, as well as the folding of three other CFTR domains. Combination of FDA approved corrector molecules can efficiently but incompletely rescue the ΔF508-CFTR folding and stability defect. Thus, new pharmacophores that would reinstate the wild-type-like conformational stability of the ΔF508-NBD1 would be highly beneficial. The most prominent molecule, 5-bromoindole-3-acetic acid (BIA) that can thermally stabilize the NBD1 has low potency and efficacy. To gain insights into the NBD1 (un)folding dynamics and BIA binding site localization, we combined molecular dynamics (MD) simulations, atomic force spectroscopy (AFM) and hydrogen-deuterium exchange (HDX) experiments. We found that the NBD1 α-subdomain with three adjacent strands from the β-subdomain plays an important role in early folding steps, when crucial non-native interactions are formed via residue F508. Our AFM and HDX experiments showed that BIA associates with this α-core region and increases the resistance of the ΔF508-NBD1 against mechanical unfolding, a phenomenon that could be exploited in future developments of folding correctors.

Experimental demonstration and pan-structurome prediction of climate-associated riboSNitches in Arabidopsis

Background

Genome-wide association studies (GWAS) aim to correlate phenotypic changes with genotypic variation. Upon transcription, single nucleotide variants (SNVs) may alter mRNA structure, with potential impacts on transcript stability, macromolecular interactions, and translation. However, plant genomes have not been assessed for the presence of these structure-altering polymorphisms or “riboSNitches.”

Results

We experimentally demonstrate the presence of riboSNitches in transcripts of two Arabidopsis genes, ZINC RIBBON 3 (ZR3) and COTTON GOLGI-RELATED 3 (CGR3), which are associated with continentality and temperature variation in the natural environment. These riboSNitches are also associated with differences in the abundance of their respective transcripts, implying a role in regulating the gene's expression in adaptation to local climate conditions. We then computationally predict riboSNitches transcriptome-wide in mRNAs of 879 naturally inbred Arabidopsis accessions. We characterize correlations between SNPs/riboSNitches in these accessions and 434 climate descriptors of their local environments, suggesting a role of these variants in local adaptation. We integrate this information in CLIMtools V2.0 and provide a new web resource, T-CLIM, that reveals associations between transcript abundance variation and local environmental variation.

Conclusion

We functionally validate two plant riboSNitches and, for the first time, demonstrate riboSNitch conditionality dependent on temperature, coining the term “conditional riboSNitch.” We provide the first pan-genome-wide prediction of riboSNitches in plants. We expand our previous CLIMtools web resource with riboSNitch information and with 1868 additional Arabidopsis genomes and 269 additional climate conditions, which will greatly facilitate in silico studies of natural genetic variation, its phenotypic consequences, and its role in local adaptation.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13059-022-02656-4.

PPARγ Gene Polymorphisms, Metabolic Disorders, and Coronary Artery Disease

Being activated by endogenous and exogenous ligands, nuclear receptor peroxisome proliferator-activated receptor gamma (PPARγ) enhances insulin sensitivity, promotes adipocyte differentiation, stimulates adipogenesis, and has the properties of anti-atherosclerosis, anti-inflammation, and anti-oxidation. The Human PPARγ gene (PPARG) contains thousands of polymorphic loci, among them two polymorphisms (rs10865710 and rs7649970) in the promoter region and two polymorphisms (rs1801282 and rs3856806) in the exonic region were widely reported to be significantly associated with coronary artery disease (CAD). Mechanistically, PPARG polymorphisms lead to abnormal expression of PPARG gene and/or dysfunction of PPARγ protein, causing metabolic disorders such as hypercholesterolemia and hypertriglyceridemia, and thereby increasing susceptibility to CAD.

Tetra-primers ARMS-PCR Based Association Analyses of Synonymous and Intronic Variants in the ADAM12 Gene with Susceptibility to Knee Osteoarthritis: A Case-Control Study

Genetic variations in a disintegrin and metalloprotease 12 (ADAM12) gene may contribute to develop Osteoarthritis (OA) that is characterized by cartilage matrix degradation and osteophytes formation. Therefore, the aim of present study was to analyze the association between the ADAM12 gene variants and knee OA predisposition. Tetra-primers ARMS-PCR was employed, to genotype the ADAM12 gene polymorphisms (rs1044122 and rs1871054) in 400 knee OA patients and equal number of age-matched controls. The association between ADAM12 gene variants and OA susceptibility was estimated using the Chi-square, logistic regression, haplotypes and linkage analyses. A significant association of rs1044122 (genotype: χ² = 18.94; P < 0.001, allele: χ² = 19.10; P < 0.001) and rs1871054 (genotype: χ² = 10.04; P = 0.007, allele: χ² = 10.57; P = 0.001) was observed with increased OA susceptibility. The variant genotype of rs1044122 increased OA risk more than twice [odds ratio (OR) 2.20; P = 0.001] and the risk was higher in females (OR 2.43; P = 0.001). The variant genotype of rs1871054 was perceived to almost double the risk in females (OR 1.97; P = 0.003). Moreover, a significant association of rs1044122 and rs1871054 under the additive genetic model (P < 0.001 and P = 0.002, respectively) was observed. The targeted ADAM12 gene polymorphisms, showed significant association with knee OA susceptibility. Females harboring the polymorphisms might be at risk. Besides, the haplotype CC of rs1044122 and rs1871054 in the ADAM12 gene may double knee OA risk. These findings may help in determining the etiology of OA and recognizing the people at risk of developing knee OA.

Clinical and molecular characterization of a patient with MBTPS1 related spondyloepiphyseal dysplasia: Evidence of pathogenicity for a synonymous variant

BACKGROUND

A novel autosomal recessive skeletal dysplasia resulting from pathogenic variants in membrane-bound transcription factor peptidase, site 1 (MBTPS1) has been recently delineated. To date, only three patients have been reported.

METHODS

In this study, we reported the clinical and molecular features of a Chinese boy who was diagnosed with spondyloepiphyseal dysplasia. The effects of variants on mRNA splicing were analyzed through transcript analysis in vivo and minigene splice assay in vitro.

RESULTS

The proband mainly showed short stature, special facial features, cataract, hernias, and serious sleep apnea syndrome. Growth hormone stimulation tests suggested the boy had growth hormone deficiency. Imaging examinations suggested abnormal thoracolumbar vertebrae and severely decreased bone mineral density. Genetic analysis of MBTPS1 gene revealed two novel heterozygous variants, a nonsense mutation c.2656C > T (p.Q886*, 167) in exon 20 and a synonymous variant c.774C > T (p.A258=) in exon 6. The transcript analysis in vivo exhibited that the synonymous variant c.774C > T caused exon 6 skipping. The minigene splice assay in vitro confirmed the alteration of MBTPS1 mRNA splicing and the exon skipping was partially restored by an antisense oligonucleotide (ASO) treatment.

CONCLUSION

Notably, we report a Chinese rare case of spondyloepiphyseal dysplasia and validate its pathogenic synonymous variant in the MBTPS1 gene.

A synonymous mutation in IGF-1 impacts the transcription and translation process of gene expression

Insulin-like growth factor 1 (IGF-1) is considered to be a crucial gene in the animal development of bone and body size. In this study, a unique synonymous mutation (c.258 A > G) of the IGF-1 gene was modified with an adenine base editor to observe the growth and developmental situation of mutant mice. Significant expression differences and molecular mechanisms among vectors with different alanine synonymous codons were explored. Although modification of a single synonymous codon rarely interferes with animal phenotypes, we observed that the expression and secretion of IGF-1 were different between 8-week-old homozygous (Ho) and wild-type (WT) mice. In addition, the IGF-1 with optimal codon combinations showed a higher expression content than other codon combination modes at both transcription and translation levels and performed proliferation promotion. The gene stability and translation initiation efficiency also changed significantly. Our findings illustrated that the synonymous mutation altered the IGF-1 gene expression in individual mice and suggested that the synonymous mutation affected the IGF-1 expression and biological function through the transcription and translation processes.

Functional analysis of novel genetic variants of NKX2-5 associated with nonsyndromic congenital heart disease

NKX2-5, a master cardiac regulatory transcription factor was the first known genetic cause of congenital heart diseases (CHDs). To further investigate its role in CHD pathogenesis, we performed mutational screening of 285 CHD probands and 200 healthy controls. Five coding sequence variants were identified in six CHD cases (2.1%), including three in the N-terminal region (p.A61G, p.R95L, and p.E131K) and one each in homeodomain (HD) (p.A148E) and tyrosine-rich domain (p.P247A). Variant-p.A148E showed tertiary structure changes and differential DNA binding affinity of mutant compared to wild type. Two N-terminal variants-p.A61G and p.E131K along with HD variant p.A148E demonstrated significantly reduced transcriptional activity of Nppa and Actc1 promoters in dual luciferase promoter assay supported by their reduced expression in qRT-PCR. Nonetheless, variant p.R95L affected the synergy of NKX2-5 with serum response factor and TBX5 leading to significantly decreased Actc1 promoter activity depicting a distinctive role of this region. The aberrant expression of other target genes-Irx4, Mef2c, Bmp10, Myh6, Myh7, and Myocd is also observed in response to NKX2-5 variants, possibly due to the defective gene regulatory network. Severely impaired downstream promoter activities and abnormal expression of target genes due to N-terminal variants supports the emerging role of this region during cardiac-developmental pathways.

Are Nonsynonymous Transversions Generally More Deleterious than Nonsynonymous Transitions?

It has been suggested that, due to the structure of the genetic code, nonsynonymous transitions are less likely than transversions to cause radical changes in amino acid physicochemical properties so are on average less deleterious. This view was supported by some but not all mutagenesis experiments. Because laboratory measures of fitness effects have limited sensitivities and relative frequencies of different mutations in mutagenesis studies may not match those in nature, we here revisit this issue using comparative genomics. We extend the standard codon model of sequence evolution by adding the parameter η that quantifies the ratio of the fixation probability of transitional nonsynonymous mutations to that of transversional nonsynonymous mutations. We then estimate η from the concatenated alignment of all protein-coding DNA sequences of two closely related genomes. Surprisingly, η ranges from 0.13 to 2.0 across 90 species pairs sampled from the tree of life, with 51 incidences of η < 1 and 30 incidences of η >1 that are statistically significant. Hence, whether nonsynonymous transversions are overall more deleterious than nonsynonymous transitions is species-dependent. Because the corresponding groups of amino acid replacements differ between nonsynonymous transitions and transversions, η is influenced by the relative exchangeabilities of amino acid pairs. Indeed, an extensive search reveals that the large variation in η is primarily explainable by the recently reported among-species disparity in amino acid exchangeabilities. These findings demonstrate that genome-wide nucleotide substitution patterns in coding sequences have species-specific features and are more variable among evolutionary lineages than are currently thought.

A synonymous variant in MYO15A enriched in the Ashkenazi Jewish population causes autosomal recessive hearing loss due to abnormal splicing

Nonsyndromic hearing loss is genetically heterogeneous. Despite comprehensive genetic testing, many cases remain unsolved because the clinical significance of identified variants is uncertain or because biallelic pathogenic variants are not identified for presumed autosomal recessive cases. Common synonymous variants are often disregarded. Determining the pathogenicity of synonymous variants may improve genetic diagnosis. We report a synonymous variant c.9861 C > T/p.(Gly3287=) in MYO15A in homozygosity or compound heterozygosity with another pathogenic or likely pathogenic MYO15A variant in 10 unrelated families with nonsyndromic sensorineural hearing loss. Biallelic variants in MYO15A were identified in 21 affected and were absent in 22 unaffected siblings. A mini-gene assay confirms that the synonymous variant leads to abnormal splicing. The variant is enriched in the Ashkenazi Jewish population. Individuals carrying biallelic variants involving c.9861 C > T often exhibit progressive post-lingual hearing loss distinct from the congenital profound deafness typically associated with biallelic loss-of-function MYO15A variants. This study establishes the pathogenicity of the c.9861 C > T variant in MYO15A and expands the phenotypic spectrum of MYO15A-related hearing loss. Our work also highlights the importance of multicenter collaboration and data sharing to establish the pathogenicity of a relatively common synonymous variant for improved diagnosis and management of hearing loss.

Role of Synonymous Mutations in the Evolution of TEM β-Lactamase Genes

Nonsynonymous mutations are well documented in TEM β-lactamases. The resulting amino acid changes often alter the conferred phenotype from broad spectrum (2b) conferred by TEM-1 to extended spectrum (2be), inhibitor resistant (2br), or both extended spectrum and inhibitor resistant (2ber). The encoding bla _TEM genes also deviate in numerous synonymous mutations, which are not well understood. bla _TEM-3 (2be), bla _TEM-33 (2br), and bla _TEM-109 (2ber) were studied in comparison to bla _TEM-1 bla _TEM-33 was chosen for more detailed studies because it deviates from bla _TEM-1 by a single nonsynonymous mutation and three additional synonymous mutations. Genes encoding the enzymes with only nonsynonymous or all (including synonymous) mutations plus all permutations between bla _TEM-1 and bla _TEM-33 were expressed in Escherichia coli cells. In disc diffusion assays, genes encoding TEM-3, TEM-33, and TEM-109 with all synonymous mutations resulted in higher resistance levels than genes without synonymous mutations. Disc diffusion assays with the 16 genes carrying all possible nucleotide change combinations between bla _TEM-1 and bla _TEM-33 indicated different susceptibilities for different variants. Nucleotide BLAST searches did not identify genes without synonymous mutations but did identify some without nonsynonymous mutations. Energies of possible secondary mRNA structures calculated with mfold are generally higher with synonymous mutations, suggesting that their role could be to destabilize the mRNA and facilitate its unfolding for efficient translation. In summary, our data indicate that transition from bla _TEM-1 to other variant genes by simply acquiring the nonsynonymous mutations is not favored. Instead, synonymous mutations seem to support the transition to other variant genes with nonsynonymous mutations leading to different phenotypes.

Synonymous variants that disrupt messenger RNA structure are significantly constrained in the human population

Background

The role of synonymous single-nucleotide variants in human health and disease is poorly understood, yet evidence suggests that this class of “silent” genetic variation plays multiple regulatory roles in both transcription and translation. One mechanism by which synonymous codons direct and modulate the translational process is through alteration of the elaborate structure formed by single-stranded mRNA molecules. While tools to computationally predict the effect of non-synonymous variants on protein structure are plentiful, analogous tools to systematically assess how synonymous variants might disrupt mRNA structure are lacking.

Results

We developed novel software using a parallel processing framework for large-scale generation of secondary RNA structures and folding statistics for the transcriptome of any species. Focusing our analysis on the human transcriptome, we calculated 5 billion RNA-folding statistics for 469 million single-nucleotide variants in 45,800 transcripts. By considering the impact of all possible synonymous variants globally, we discover that synonymous variants predicted to disrupt mRNA structure have significantly lower rates of incidence in the human population.

Conclusions

These findings support the hypothesis that synonymous variants may play a role in genetic disorders due to their effects on mRNA structure. To evaluate the potential pathogenic impact of synonymous variants, we provide RNA stability, edge distance, and diversity metrics for every nucleotide in the human transcriptome and introduce a “Structural Predictivity Index” (SPI) to quantify structural constraint operating on any synonymous variant. Because no single RNA-folding metric can capture the diversity of mechanisms by which a variant could alter secondary mRNA structure, we generated a SUmmarized RNA Folding (SURF) metric to provide a single measurement to predict the impact of secondary structure altering variants in human genetic studies.

Co-Translational Folding of the First Transmembrane Domain of ABC-Transporter CFTR is Supported by Assembly with the First Cytosolic Domain

ABC transporters transport a wealth of molecules across membranes and consist of transmembrane and cytosolic domains. Their activity cycle involves a tightly regulated and concerted domain choreography. Regulation is driven by the cytosolic domains and function by the transmembrane domains. Folding of these polytopic multidomain proteins to their functional state is a challenge for cells, which is mitigated by co-translational and sequential events. We here reveal the first stages of co-translational domain folding and assembly of CFTR, the ABC transporter defective in the most abundant rare inherited disease cystic fibrosis. We have combined biosynthetic radiolabeling with protease-susceptibility assays and domain-specific antibodies. The most N-terminal domain, TMD1 (transmembrane domain 1), folds both its hydrophobic and soluble helices during translation: the transmembrane helices pack tightly and the cytosolic N- and C-termini assemble with the first cytosolic helical loop ICL1, leaving only ICL2 exposed. This N-C-ICL1 assembly is strengthened by two independent events: (i) assembly of ICL1 with the N-terminal subdomain of the next domain, cytosolic NBD1 (nucleotide-binding domain 1); and (ii) in the presence of corrector drug VX-809, which rescues cell-surface expression of a range of disease-causing CFTR mutants. Both lead to increased shielding of the CFTR N-terminus, and their additivity implies different modes of action. Early assembly of NBD1 and TMD1 is essential for CFTR folding and positions both domains for the required assembly with TMD2. Altogether, we have gained insights into this first, nucleating, VX-809-enhanced domain-assembly event during and immediately after CFTR translation, involving structures conserved in type-I ABC exporters.

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

Epistasis refers to the dependence of a mutation on other mutation(s) and the genetic context in general. In the context of human disorders, epistasis complicates the spectrum of disease symptoms and has been proposed as a major contributor to variations in disease outcome. The nonadditive relationship between mutations and the lack of complete understanding of the underlying physiological effects limit our ability to predict phenotypic outcome. Here, we report positive epistasis between intragenic mutations in the cystic fibrosis transmembrane conductance regulator (CFTR)-the gene responsible for cystic fibrosis (CF) pathology. We identified a synonymous single-nucleotide polymorphism (sSNP) that is invariant for the CFTR amino acid sequence but inverts translation speed at the affected codon. This sSNP in cis exhibits positive epistatic effects on some CF disease-causing missense mutations. Individually, both mutations alter CFTR structure and function, yet when combined, they lead to enhanced protein expression and activity. The most robust effect was observed when the sSNP was present in combination with missense mutations that, along with the primary amino acid change, also alter the speed of translation at the affected codon. Functional studies revealed that synergistic alteration in ribosomal velocity is the underlying mechanism; alteration of translation speed likely increases the time window for establishing crucial domain-domain interactions that are otherwise perturbed by each individual mutation.

Synonymous variants in holoprosencephaly alter codon usage and impact the Sonic Hedgehog protein

Synonymous single nucleotide variants (sSNVs) have been implicated in various genetic disorders through alterations of pre-mRNA splicing, mRNA structure and miRNA regulation. However, their impact on synonymous codon usage and protein translation remains to be elucidated in clinical context. Here, we explore the functional impact of sSNVs in the Sonic Hedgehog (SHH) gene, identified in patients affected by holoprosencephaly, a congenital brain defect resulting from incomplete forebrain cleavage. We identified eight sSNVs in SHH, selectively enriched in holoprosencephaly patients as compared to healthy individuals, and systematically assessed their effect at both transcriptional and translational levels using a series of in silico and in vitro approaches. Although no evidence of impact of these sSNVs on splicing, mRNA structure or miRNA regulation was found, five sSNVs introduced significant changes in codon usage and were predicted to impact protein translation. Cell assays demonstrated that these five sSNVs are associated with a significantly reduced amount of the resulting protein, ranging from 5% to 23%. Inhibition of the proteasome rescued the protein levels for four out of five sSNVs, confirming their impact on protein stability and folding. Remarkably, we found a significant correlation between experimental values of protein reduction and computational measures of codon usage, indicating the relevance of in silico models in predicting the impact of sSNVs on translation. Considering the critical role of SHH in brain development, our findings highlight the clinical relevance of sSNVs in holoprosencephaly and underline the importance of investigating their impact on translation in human pathologies.

Is CFTR an exchanger?: Regulation of HCO3 -Transport and extracellular pH by CFTR

The disease, cystic fibrosis, is caused by the malfunction of the cystic fibrosis transmembrane conductance regulator. Expression of functional CFTR may normally regulate extracellular pH via control of bicarbonate efflux. Reports also suggest that the CFTR may be a Cl-/HCO3- exchanger. If true, this could be very important for treatment of CF given the airway host defense system is quite sensitive to pH, and acidic pH been found to increase mucus viscosity. We compared evidentiary support of four possible models of CFTR's role in the transport of bicarbonate: 1) CFTR as a Cl-channel that permits bicarbonate conductance, 2) CFTR as an anion Cl-/HCO3- exchanger (AE), 3.) CFTR as both a Cl-channel and an AE, and 4.) CFTR as a Cl-channel that allows for transport of bicarbonate and regulates an independent AE. The effect of stimulators and inhibitors of CFTR and AEs were evaluated via iodide efflux and studies of extracellular pH. This data, as well as that published by others, suggest that while CFTR may support and regulate bicarbonate flux it is unlikely it directly performs Cl-/HCO3- anion exchange.

•

Four possible models of CFTR's role in the transport of HC03- were evaluated.

•

Cells exposed to modulators of CFTR/AE activity were studied by efflux and pHo.

•

Efflux: While CFTR inhibition reduced anion efflux, AE inhibition did not.

•

pHo: In Cl-free buffer conditions CFTR activation reduced pHo acidification.

•

CFTR may regulate but unlikely it directly performs Cl-/HCO3- exchange.

Dysregulated PDGFR alpha expression and novel somatic mutations in colorectal cancer: association to RAS wild type status and tumor size

Background

Platelet derived growth factor receptor alpha (PDGFRα) has been considered as a relevant factor in tumor proliferation, angiogenesis and metastatic dissemination. It was a target of tyrosine kinase (TK) inhibitors emerged in the therapy of diverse cancers. In colorectal cancer, the commonly used therapy is anti-epithelial growth factor receptor (EGFR). However, both RAS mutated and a subgroup of RAS wild type patients resist to such therapy. The aim of this study is to investigate PDGFRα protein expression and mutational status in colorectal adenocarcinoma and their association with clinicopathological features and molecular RAS status to provide useful information for the identification of an effective biomarker that might be implicated in prognosis and treatment prediction.

Methods

Our study enrolled 103 formalin fixed paraffin-embedded (FFPE) colorectal adenocarcinoma. PDGFRα expression was investigated by immunohistochemistry (IHC). Hotspot exon 18 of PDGFRA was studied by PCR followed by Sanger sequencing and RAS status was determined by real time quantitative PCR. Thirteen normal colon tissues were used as negative controls.

Results

PDGFRα staining was detected in the cytoplasm of all tissues. Low expression was observed in all normal colon mucosa. In adenocarcinoma, 45% (45/100) of cases showed PDGFRα overexpression. This overexpression was significantly associated with mutations in exon 18 (P = 0.024), RAS wild type status (P < 10^–3), tumor diameter (P = 0.048), whereas there was no association with tumor side (P = 0.13) and other clinicopathological features.

Conclusion

Overexpression of PDGFRα in adenocarcinoma suggests its potential role in tumor cells growth and invasion. The association between PDGFRα overexpression in both tumor and stromal adenocarcinoma cells with RAS wild type status suggests its potential role in anti-EGFR therapy resistance and the relevance of using it as specific or adjuvant therapeutic target.

The Effects of Single Nucleotide Polymorphisms in Cancer RNAi Therapies

Tremendous progress in RNAi delivery methods and design has allowed for the effective development of siRNA-based therapeutics that are currently under clinical investigation for various cancer treatments. This approach has the potential to revolutionize cancer therapy by providing the ability to specifically downregulate or upregulate the mRNA of any protein of interest. This exquisite specificity, unfortunately, also has a downside. Genetic variations in the human population are common because of the presence of single nucleotide polymorphisms (SNPs). SNPs lead to synonymous and non-synonymous changes and they occur once in every 300 base pairs in both coding and non-coding regions in the human genome. Much less common are the somatic mosaicism variations associated with genetically distinct populations of cells within an individual that is derived from postzygotic mutations. These heterogeneities in the population can affect the RNAi's efficacy or more problematically, which can lead to unpredictable and sometimes adverse side effects. From a more positive viewpoint, both SNPs and somatic mosaicisms have also been implicated in human diseases, including cancer, and these specific changes could offer the ability to effectively and, more importantly, selectively target the cancer cells. In this review, we discuss how SNPs in the human population can influence the development and success of novel anticancer RNAi therapies and the importance of why SNPs should be carefully considered.

A synonymous UPF3B variant causing a speech disorder implicates NMD as a regulator of neurodevelopmental disorder gene networks

Loss-of-function mutations of the X-chromosome gene UPF3B cause male neurodevelopmental disorders (NDDs) via largely unknown mechanisms. We investigated initially by interrogating a novel synonymous UPF3B variant in a male with absent speech. In silico and functional studies using cell lines derived from this individual show altered UPF3B RNA splicing. The resulting mRNA species encodes a frame-shifted protein with a premature termination codon (PTC) predicted to elicit degradation via nonsense-mediated mRNA decay (NMD). UPF3B mRNA was reduced in the cell line, and no UPF3B protein was produced, confirming a loss-of-function allele. UPF3B is itself involved in the NMD mechanism which degrades both PTC-bearing mutant transcripts and also many physiological transcripts. RNAseq analysis showed that ~1.6% of mRNAs exhibited altered expression. These mRNA changes overlapped and correlated with those we identified in additional cell lines obtained from individuals harbouring other UPF3B mutations, permitting us to interrogate pathogenic mechanisms of UPF3B-associated NDDs. We identified 102 genes consistently deregulated across all UPF3B mutant cell lines. Of the 51 upregulated genes, 75% contained an NMD-targeting feature, thus identifying high-confidence direct NMD targets. Intriguingly, 22 of the dysregulated genes encoded known NDD genes, suggesting UPF3B-dependent NMD regulates gene networks critical for cognition and behaviour. Indeed, we show that 78.5% of all NDD genes encode a transcript predicted to be targeted by NMD. These data describe the first synonymous UPF3B mutation in a patient with prominent speech and language disabilities and identify plausible mechanisms of pathology downstream of UPF3B mutations involving the deregulation of NDD-gene networks.

Identification and analysis of RNA structural disruptions induced by single nucleotide variants using Riprap and RiboSNitchDB

RNA conformational alteration has significant impacts on cellular processes and phenotypic variations. An emerging genetic factor of RNA conformational alteration is a new class of single nucleotide variant (SNV) named riboSNitch. RiboSNitches have been demonstrated to be involved in many genetic diseases. However, identifying riboSNitches is notably difficult as the signals of RNA structural disruption are often subtle. Here, we introduce a novel computational framework–RIboSNitch Predictor based on Robust Analysis of Pairing probabilities (Riprap). Riprap identifies structurally disrupted regions around any given SNVs based on robust analysis of local structural configurations between wild-type and mutant RNA sequences. Compared to previous approaches, Riprap shows higher accuracy when assessed on hundreds of known riboSNitches captured by various experimental RNA structure probing methods including the parallel analysis of RNA structure (PARS) and the selective 2′-hydroxyl acylation analyzed by primer extension (SHAPE). Further, Riprap detects the experimentally validated riboSNitch that regulates human catechol-O-methyltransferase haplotypes and outputs structurally disrupted regions precisely at base resolution. Riprap provides a new approach to interpreting disease-related genetic variants. In addition, we construct a database (RiboSNitchDB) that includes the annotation and visualization of all presented riboSNitches in this study as well as 24 629 predicted riboSNitches from human expression quantitative trait loci.

EST-SNP Study of Olea europaea L. Uncovers Functional Polymorphisms between Cultivated and Wild Olives

BACKGROUND

The species Olea europaea includes cultivated varieties (subsp. europaea var. europaea), wild plants (subsp. europaea var. sylvestris), and five other subspecies spread over almost all continents. Single nucleotide polymorphisms in the expressed sequence tag able to underline intra-species differentiation are not yet identified, beyond a few plastidial markers.

METHODS

In the present work, more than 1000 transcript-specific SNP markers obtained by the genotyping of 260 individuals were studied. These genotypes included cultivated, oleasters, and samples of subspecies guanchica, and were analyzed in silico, in order to identify polymorphisms on key genes distinguishing different Olea europaea forms.

RESULTS

Phylogeny inference and principal coordinate analysis allowed to detect two distinct clusters, clearly separating wilds and guanchica samples from cultivated olives, meanwhile the structure analysis made possible to differentiate these three groups. Sequences carrying the polymorphisms that distinguished wild and cultivated olives were analyzed and annotated, allowing to identify 124 candidate genes that have a functional role in flower development, stress response, or involvement in important metabolic pathways. Signatures of selection that occurred during olive domestication, were detected and reported.

CONCLUSION

This deep EST-SNP analysis provided important information on the genetic and genomic diversity of the olive complex, opening new opportunities to detect gene polymorphisms with potential functional and evolutionary roles, and to apply them in genomics-assisted breeding, highlighting the importance of olive germplasm conservation.

A Synonymous Exonic Splice Silencer Variant in IRF6 as a Novel and Cryptic Cause of Non-Syndromic Cleft Lip and Palate

Missense, nonsense, splice site and regulatory region variants in interferon regulatory factor 6 (IRF6) have been shown to contribute to both syndromic and non-syndromic forms of cleft lip and/or palate (CL/P). We report the diagnostic evaluation of a complex multigeneration family of Honduran ancestry with a pedigree structure consistent with autosomal-dominant inheritance with both incomplete penetrance and variable expressivity. The proband's grandmother bore children with two partners and CL/P segregates on both sides of each lineage. Through whole-exome sequencing of five members of the family, we identified a single shared synonymous variant, located in the middle of exon 7 of IRF6 (p.Ser307Ser; g.209963979 G>A; c.921C>T). The variant was shown to segregate in the seven affected individuals and through three unaffected obligate carriers, spanning both sides of this pedigree. This variant is very rare, only being found in three (all of Latino ancestry) of 251,352 alleles in the gnomAD database. While the variant did not create a splice acceptor/donor site, in silico analysis predicted it to impact an exonic splice silencer element and the binding of major splice regulatory factors. In vitro splice assays supported this by revealing multiple abnormal splicing events, estimated to impact >60% of allelic transcripts. Sequencing of the alternate splice products demonstrated the unmasking of a cryptic splice site six nucleotides 5' of the variant, as well as variable utilization of cryptic splice sites in intron 6. The ectopic expression of different splice regulatory proteins altered the proportion of abnormal splicing events seen in the splice assay, although the alteration was dependent on the splice factor. Importantly, each alternatively spliced mRNA is predicted to result in a frame shift and prematurely truncated IRF6 protein. This is the first study to identify a synonymous variant as a likely cause of NS-CL/P and highlights the care that should be taken by laboratories when considering and interpreting variants.

Slowing ribosome velocity restores folding and function of mutant CFTR

Cystic fibrosis (CF) is caused by mutations in the CF transmembrane conductance regulator (CFTR), with approximately 90% of patients harboring at least one copy of the disease-associated variant F508del. We utilized a yeast phenomic system to identify genetic modifiers of F508del-CFTR biogenesis, from which ribosomal protein L12 (RPL12/uL11) emerged as a molecular target. In the present study, we investigated mechanism(s) by which suppression of RPL12 rescues F508del protein synthesis and activity. Using ribosome profiling, we found that rates of translation initiation and elongation were markedly slowed by RPL12 silencing. However, proteolytic stability and patch-clamp assays revealed RPL12 depletion significantly increased F508del-CFTR steady-state expression, interdomain assembly, and baseline open-channel probability. We next evaluated whether Rpl12-corrected F508del-CFTR could be further enhanced with concomitant pharmacologic repair (e.g., using clinically approved modulators lumacaftor and tezacaftor) and demonstrated additivity of these treatments. Rpl12 knockdown also partially restored maturation of specific CFTR variants in addition to F508del, and WT Cftr biogenesis was enhanced in the pancreas, colon, and ileum of Rpl12 haplosufficient mice. Modulation of ribosome velocity therefore represents a robust method for understanding both CF pathogenesis and therapeutic response.

In silico analysis on the functional and structural impact of Rad50 mutations involved in DNA strand break repair

BACKGROUND

DNA double strand break repair is important to preserve the fidelity of our genetic makeup after DNA damage. Rad50 is one of the components in MRN complex important for DNA repair mechanism. Rad50 mutations can lead to microcephaly, mental retardation and growth retardation in human. However, Rad50 mutations in human and other organisms have never been gathered and heuristically compared for their deleterious effects. It is important to assess the conserved region in Rad50 and its homolog to identify vital mutations that can affect functions of the protein.

METHOD

In this study, Rad50 mutations were retrieved from SNPeffect 4.0 database and literature. Each of the mutations was analyzed using various bioinformatic analyses such as PredictSNP, MutPred, SNPeffect 4.0, I-Mutant and MuPro to identify its impact on molecular mechanism, biological function and protein stability, respectively.

RESULTS

We identified 103 mostly occurred mutations in the Rad50 protein domains and motifs, which only 42 mutations were classified as most deleterious. These mutations are mainly situated at the specific motifs such as Walker A, Q-loop, Walker B, D-loop and signature motif of the Rad50 protein. Some of these mutations were predicted to negatively affect several important functional sites that play important roles in DNA repair mechanism and cell cycle signaling pathway, highlighting Rad50 crucial role in this process. Interestingly, mutations located at non-conserved regions were predicted to have neutral/non-damaging effects, in contrast with previous experimental studies that showed deleterious effects. This suggests that software used in this study may have limitations in predicting mutations in non-conserved regions, implying further improvement in their algorithm is needed. In conclusion, this study reveals the priority of acid substitution associated with the genetic disorders. This finding highlights the vital roles of certain residues such as K42E, C681A/S, CC684R/S, S1202R, E1232Q and D1238N/A located in Rad50 conserved regions, which can be considered for a more targeted future studies.

A review of databases predicting the effects of SNPs in miRNA genes or miRNA-binding sites

Modern precision medicine comprises the knowledge and understanding of individual differences in the genomic sequence of patients to provide tailor-made treatments. Regularly, such variants are considered in coding regions only, and their effects are predicted based on their impact on the amino acid sequence of expressed proteins. However, assessing the effects of variants in noncoding elements, in particular microRNAs (miRNAs) and their binding sites, is important as well, as a single miRNA can influence the expression patterns of many genes at the same time. To analyze the effects of variants in miRNAs and their target sites, several databases storing variant impact predictions have been published. In this review, we will compare the core functionalities and features of these databases and discuss the importance of up-to-date data resources in the context of web applications. Finally, we will outline some recommendations for future developments in the field.

TERT-rs33963617 and CLPTM1L-rs77518573 reduce the risk of non-small cell lung cancer in Chinese population

BACKGROUND

Genome-wide association studies (GWAS) have identified 5p15.33 as a susceptible locus for lung cancer. However, for non-small cell lung cancer (NSCLC), low-frequency risk variants in this region have not been systematically studied. We intended to explore the associations between low-frequency variants on 5p15.33 and NSCLC using a next-generation sequencing based approach in this study.

METHODS

We have acquisited the variation spectrum of 400 NSCLC patients on 5p15.33 by sequencing the targeted region before. Candidate variants were primarily selected by restricting the minor allele frequency (MAF 1-5%) and then by comparing their frequency in 400 NSCLC patients with 1008 East Asians from The genome Aggregation Database (gnomAD). The associations between candidate variants and NSCLC were discovered and replicated in two case-control sets: discovery stage with 960 cases and 916 controls, and replication stage with 1596 cases and 1614 controls in total.

RESULTS

Five low-frequency variants were selected as our candidates and subsequent association analyses showed that 2 polymorphisms were significantly associated with risk of NSCLC, including rs33963617 (OR = 0.63, 95% CI: 0.53-0.76, P = 3.80 × 10^-7) in TERT and rs77518573 (OR = 0.73, 95% CI: 0.63-0.84, P = 2.00 × 10^-5) in upstream of CLPTM1L. When stratified by histologic subtype, a significant association was only investigated in adenocarcinoma for rs77518573. We also observed an obvious cumulative effect of the two significant variants.

CONCLUSIONS

We newly identified two NSCLC related variants on chromosome 5p15.33. Both TERT-rs33963617 and CLPTM1L-rs77518573 conferred reduced risk for NSCLC in Chinese Han population.