Exposing synonymous mutations

Toward a rational therapeutic for elastin related disease: Key considerations for elastin based regenerative medicine strategies

Elastin is a connective tissue protein, produced from the ELN gene, that provides elasticity and recoil to tissues that stretch, such as the large arteries of the body, lung parenchyma, skin, ligaments and elastic cartilages. It is produced as a soluble monomer, tropoelastin, that when cross-linked in the extracellular space generates a polymer that is extraordinarily stable, with a predicted half-life of >70 years. Although data suggest ongoing elastin transcription, it is rare to see new elastin deposited outside of its tight developmental window. Consequently, elastin-related disease comes about primarily in one of three scenarios: (1) inadequate elastin deposition, (2) production of poor-quality elastic fibers, or (3) increased destruction of previously deposited elastin. By understanding the pathways controlling elastin production and maintenance, we can design new therapeutics to thwart those abnormal processes. In this review, we will summarize the diseases arising from genetic and environmental alteration of elastin (Williams syndrome, supravalvar aortic stenosis, autosomal dominant cutis laxa, and ELN-related vascular and connective tissue dysfunction) and then describe the mechanisms controlling elastin production and maintenance that might be manipulated to generate novel therapeutics aimed at these conditions. We will end by summarizing existing therapeutic strategies targeting these disease mechanisms before outlining future approaches that may better solve the challenges associated with elastin based regenerative medicine.

Reanalysis of Exome Sequencing Data in the Indian Undiagnosed Diseases Program: Improving Diagnostic Yield and Ending Diagnostic Odyssey

In 2021, the Indian Undiagnosed Diseases Program was initiated for patients without a definite diagnosis despite extensive evaluation in four participating sites. Between February 2021 and March 2023, a total of 88 patients were recruited and underwent deep phenotyping. A uniform methodology for data re-analysis was implemented as the first step prior to conducting additional genomic testing. The largest cohort was of 38 patients with neurodevelopmental disorders (NDD). A genetic diagnosis was achieved in 24 of the 88 patients (27.2%), including 7 cases within the NDD cohort. Factors contributing to the increased diagnostic yield included: (a) identification of a novel disease association in DAAM2 gene, and (b) limitations of the standard analysis pipeline, particularly for synonymous variants in SELENOI and KIAA0753 genes, non-frameshift variant in GLRX5 gene, low-coverage variant in GJC2 gene, large deletions in PCNT and PHKG2 genes, and intronic variants in VPS33B and FBN1. Improved phenotyping led to a diagnosis in three cases, while genomic variants missed in the previous bioinformatics analysis were identified in 12 cases. The study also contributed to the development of enhanced bioinformatics scripts for variant prioritization and more refined literature search for novel disease associations. It highlights the importance of incorporating data reanalysis into clinical workflows before pursuing advanced diagnostic tests, particularly in resource-limited settings where healthcare expenses are often borne out of pocket.

Relationships between liver and rumen fluke infections, milk somatic cells and polymorphisms in the Toll-like receptor 5 gene and vitamin D metabolism-related genes in Holstein dairy cows

This study investigated polymorphisms in the genes CYP3A4, CYP2R1, and TLR5, and their associations with liver fluke (Fasciola hepatica) and rumen fluke (Calicophoron / Paramphistomum spp.) infections as well as with milk somatic cell count (SCC) as an indicator for mastitis in Holstein Friesian dairy cows. DNA sequencing of the genes' exons, 5'-, 3'-untranslated regions (UTR), introns, and flanking regions of 24 cows revealed 442 variants (388 SNPs and 54 InDels) including 116 previously unknown variants. We detected three novel non-synonymous variants leading to the derived amino acid exchanges, i.e. CYP3A4 p.Gly197Ser, CYP3A4 p.Ile388Val, and CYP2R1 p.Val11Ala. The newly identified SNP 25:36589922 T > C (ss11846100002) is positioned in the splice site of CYP3A4, but showed no impact on the binding score of the splice enzymes. The CYP2R1 and TLR5 genes presented 11 SNPs in the 5'- and 3'-UTR, partly influencing transcription factor binding or microRNA target sites. Associations between polymorphisms and constructed haplotypes with infection traits were analysed via (generalized) linear mixed models including further potential confounders. In total, 109 variants in CYP3A4, 37 variants in CYP2R1, and 18 variants in TLR5 were significantly associated with F. hepatica and rumen fluke infections, and with SCC. The CYP2R1 and TLR5 variants were mostly linked to SCC, indicating the genes' roles in immune responses to bacterial infections. Haplotype analysis revealed significant associations between specific CYP3A4 haplotypes and F. hepatica worm count and faecal egg counts. This study revealed significant insights into gene polymorphisms related to vitamin D metabolism and immune response, which seem to play a role in helminth and udder infections.

FDPSM: Feature-Driven Prediction Modeling of Pathogenic Synonymous Mutations

Synonymous mutations, once considered to be biologically neutral, are now recognized to affect protein expression and function by altering the RNA splicing, stability, or translation efficiency. These effects can contribute to disease, making the prediction of the pathogenicity a crucial task. Computational methods have been developed to analyze the sequence features and biological functions of synonymous mutations, but existing methods face limitations, including scarcity of labeled data, reliance on other prediction tools, and insufficient representation of feature interrelationships. Here, we present FDPSM, a novel prediction method specifically designed to predict pathogenic synonymous mutations. FDPSM was trained on a robust data set of 4251 positive and negative training samples to enhance predictive accuracy. The method leveraged a comprehensive set of features, including genomic context, conservation, splicing effects, functional effects, and epigenomics, without relying on prediction scores from other mutation pathogenicity tools. Recognizing that original features alone may not fully capture the distinctions between pathogenic and benign synonymous mutations, we enhanced the feature set by extracting effective information from the interactions and distribution of these features. The experimental results showed that FDPSM significantly outperformed existing methods in predicting the pathogenicity of synonymous mutations, offering a more accurate and reliable tool for this important task. FDPSM is available at https://github.com/xialab-ahu/FDPSM.

Novel non-synonymous and synonymous gene variants of SRD5A2 in patients with 46,XY-DSD and DSD-free subjects

SRD5A2 gene variants are associated with deficiency of steroid 5α-reductase type 2, which is an autosomal recessive disorder of sex development (DSD) present in 46,XY males with ambiguous genitalia. To determine the causality of the disorder, this study involved genetic screening of SRD5A2 in six unrelated patients with this condition. Polymerase chain reaction (PCR) assays excluded large duplications, insertions, or deletions, while bidirectional Sanger sequencing identified 15 single-nucleotide variants (SNVs), six patients with 46,XY-DSD carrying pathogenic non-synonymous SNVs (nsSNVs), and three subjects who were DSD-free with novel synonymous SNVs (sSNVs). Genomic outcomes showed that 9 non-synonymous coding SNVs are linked to patients with SRD5A2-associated steroid 5α-reductase type 2 deficiency (c.169G > C: p.E57Q; c.145G > A: p.A49T/c.686T > C: p.F229S; c.100G > A: p.G34R/c.344G > A: p.G115D; c.591G > T: p.E197D; c.92C > T: p.S31F/c.481A > C: p.I161L (a novel missense variant; Km,app =  1.19 ±  0.1 μM, Vmax,app =  688 ±  145.8 pmol/mg P/min); c.686T > C: p.F229S). This analysis also highlighted 2 non-disease-causing sSNVs in three DSD-free subjects (c.243G > T: p.T81 = ; c.594C > T: p.I198=). These silent mutations or sSNVs in the SRD5A2 gene have no functional consequences and might not be involved in steroid 5α-reductase 2 deficiency. The identification of these sSNVs in both healthy controls and patients might suggest natural genetic variability with a very low allele frequency in the Mexican population. Furthermore, these findings indicated that nsSNVs in the SRD5A2 gene altered normal development of external male genitalia, supporting their pathogenicity.

Assessing Splicing Variants in the PAX6 Gene: A Comprehensive Minigene Approach

Haploinsufficiency of the PAX6 gene causes aniridia, a congenital eye disorder characterised by the absence or malformation of the iris and foveal hypoplasia. Previous studies indicate that pathogenic splice variants account for up to 15% of all disease-causing PAX6 variants. However, this proportion may be significantly underestimated because the pathogenicity of splice variants can only be accurately established through experimental validation. In this study, we developed and validated a system of eight minigene constructions for the functional analysis of splicing variants in the PAX6 gene. This system covers all PAX6 coding exons and allows the analysis of any exon and most intronic variants of PAX6. Our comprehensive approach, employing fragment analysis and deep targeted sequencing, enabled us to accurately characterise 38 previously described PAX6 variants, including challenging cases with multiple splicing events. The application of our system revealed that the number of pathogenic splicing variants might be closer to 30% of all pathogenic PAX6 variants. This finding considerably reshapes our understanding of their significance in the genetic landscape of aniridia.

Role of silent mutations in KRAS -mutant tumors

ABSTRACT

Silent mutations within the RAS  gene have garnered increasing attention for their potential roles in tumorigenesis and therapeutic strategies. Kirsten-RAS ( KRAS ) mutations, predominantly oncogenic, are pivotal drivers in various cancers. While extensive research has elucidated the molecular mechanisms and biological consequences of active KRAS  mutations, the functional significance of silent mutations remains relatively understudied. This review synthesizes current knowledge on KRAS  silent mutations, highlighting their impact on cancer development. Silent mutations, which do not alter protein sequences but can affect RNA stability and translational efficiency, pose intriguing questions regarding their contribution to tumor biology. Understanding these mutations is crucial for comprehensively unraveling KRAS -driven oncogenesis and exploring novel therapeutic avenues. Moreover, investigations into the clinical implications of silent mutations in KRAS -mutant tumors suggest potential diagnostic and therapeutic strategies. Despite being in early stages, research on KRAS  silent mutations holds promise for uncovering novel insights that could inform personalized cancer treatments. In conclusion, this review underscores the evolving landscape of KRAS  silent mutations, advocating for further exploration to bridge fundamental biology with clinical applications in oncology.

An essential role for TASL in mouse autoimmune pathogenesis and Toll-like receptor signaling

TASL is an immune adaptor that binds to the solute carrier SLC15A4 and facilitates activation of the transcription factor IRF5 during Toll-like receptor (TLR) signaling. Similar to IRF5 and SLC15A4, single nucleotide polymorphisms (SNPs) within TASL have been implicated in increased susceptibility to systemic lupus erythematosus (SLE) in patients. However, the biological function of TASL in vivo and how SLE-associated SNPs increase disease risk is unknown. Here we report that mice deficient in Tasl lack responses to TLR7/9 stimulation and are protected from autoimmune symptoms induced by Aldara or pristane. Loss of Tasl reduces IRF5 phosphorylation and cytokine production in multiple immune cell types but has no effect on other aspects of TLR signaling. Conversely, an SLE-associated TASL risk variant increases TASL protein expression via codon usage, resulting in augmented cytokine production in human cells. Altogether, our study validates the essential function of TASL in TLR signaling and autoimmune pathogenesis.

The adaptive value of recombination in resolving intralocus sexual conflict by gene duplication

Recombination plays a key role in increasing the efficacy of selection. We investigate whether recombination can also play a role in resolving adaptive conflicts at loci coding for traits shared between the sexes. Errors during recombination events resulting in gene duplications may provide a long-term evolutionary advantage if those loci also experience sexually antagonistic (SA) selection since, after duplication, sex-specific expression profiles will be free to evolve, thereby reducing the load on population fitness and resolving the conflict. The potential advantage of gene duplication may be tempered by the short-term deleterious effects on gamete and zygote survival, which may be tolerable in a species with high reproductive output but not with low reproductive output. We used datasets of candidate SA loci from Drosophila melanogaster and humans to test these ideas. As in humans, sexually antagonistic alleles in flies with net positive effects across the two sexes occurred at higher frequencies than alleles with net negative effects. In flies, higher recombination rates were associated with more intense levels of sexual conflict and genes with paralogues occur in regions with higher recombination rates, indicating gene duplication events are associated with a history of SA selection. Genes experiencing higher levels of conflict also showed both a higher proportion with paralogues and higher numbers of paralogues. Together, our findings reveal multiple lines of evidence for a possible route towards the resolution of an adaptive conflict via gene duplication that is facilitated by higher recombination rates.

Association of Single-Nucleotide Variants in ACE2 with the Persistence of Positive qPCR Test for SARS-CoV-2 in Healthcare Professionals During the First Wave of the COVID-19 Pandemic

The persistence of qPCR positivity for SARS-CoV-2 in individuals who recovered from COVID-19 raised several questions regarding viral transmission, with a special interest in healthcare professionals who may pose a risk of transmitting SARS-CoV-2. This issue highlights the necessity for identifying the genetic risk factors associated with persistent SARS-CoV-2 infection. A promising target for achieving this goal is the angiotensin-converting enzyme 2 (ACE2) gene, which has been associated with clinical characteristics of COVID-19 infection, such as severity. The analysis of samples from the first wave of the COVID-19 pandemic represents the initial response of the immune human system against this new virus, without the effect of vaccination or the presence of multiple strains. The aim of this study was to analyze the association of genetic variants in ACE2 with persistent SARS-CoV-2 infection. We conducted a case-control study, including 151 healthcare workers who tested positive for SARS-CoV-2 by qPCR during the first wave of the COVID-19 pandemic, and who were followed up until their results were negative. ACE2 was sequenced through Sanger sequencing. The sequence was compared against a reference sequence and variants identified. Four ACE2 variants were associated with persistent SARS-CoV-2 qPCR positivity. Three of the variants with an effect on the resulting protein were associated with increased risk of persistent SARS-CoV-2 qPCR positivity, NG_012575.2:g.35481 C>T, NG_012575.2:g.35483 G>T and NG_012575.2:g.35498 G>T. On the other hand, the rs2285666 (NG_012575.2:g.14934 G>A) was associated with a higher risk for persistent SARS-CoV-2 qPCR positivity in women and rs4646150 (NG_012575.2:g.25701 G>A) in men. The NG_012575.2:g.35498 G>T variant represents an amino acid change with a possibly harmful effect on ACE2 function. Our results suggest that ACE2 variants might be useful for identifying the population at higher risk for developing persistent SARS-CoV-2-positive qPCR results. This knowledge can be helpful for designing health policies for protecting healthcare professionals and, in consequence, users of health services.

Novel Biallelic Synonymous Exonic Variant in VPS13A Affecting mRNA Splicing: Case Report

Objectives

Chorea-acanthocytosis is an autosomal recessively inherited condition caused by loss-of-function pathogenic variants in VPS13A. We identified a novel synonymous exonic variant leading to abnormal mRNA splicing in a patient with chorea-acanthocytosis.

Methods

A patient with focal epilepsy developed generalized chorea with orolingual dystonia, cognitive decline, and peripheral neuropathy, consistent with chorea-acanthocytosis. Her parents were first cousins, but there was otherwise no family history. Targeted gene sequencing for variants in VPS13A, mRNA splicing analysis, and Western blot for chorein were performed.

Results

A homozygous synonymous variant in exon 41 of VPS13A (NM_033305.3): c.5157C>T; p.Gly1719 = was identified; this was previously classified as a variant of uncertain significance. SpliceAI predicted a splice donor gain with a score of 0.75 2 base pairs upstream of the reported variant. RNA splicing analysis revealed the creation of a type III splice variant, resulting in a frameshift and a premature termination codon. Western blot showed absent chorein/VPS13A protein.

Discussion

The variant is reclassified as likely pathogenic based on the American College of Medical Genetics criteria. This is the first reported case of ChAc caused by a synonymous variant in VPS13A proven to affect splicing. Our report further expands the spectrum of variants known to cause ChAc.

Low nucleotide diversity of the Plasmodium falciparum AP2-EXP2 gene among clinical samples from Ghana

BACKGROUND

PfAP2-EXP2 is located within chromosome 6 of Plasmodium falciparum recently identified to be undergoing an extensive selective sweep in West African isolates. The gene encoding this transcription factor, PfAP2-EXP2, is essential and thus likely subject to purifying selection that limits variants in the parasite population despite its genomic location.

METHODS

72 Plasmodium falciparum field samples and 801 clinical sequences from the Pf6 MalariaGEN dataset of Ghanaian origin, were integrated and analysed.

RESULTS

A total of 14 single nucleotide variants of which 5 were missense variants, were identified after quality checks and filtering. Except for one, all identified variants were rare among the clinical samples obtained in this study (Minor allelic frequency < 0.01). Further results revealed a considerably low dN/dS value (0.208) suggesting the presence of purifying selection. Further, all the mutant amino acids were wildtype residues in AP2-EXP2 orthologous proteins-tentatively suggesting a genus-level conservation of amino acid residues. Computational analysis and predictions corroborated these findings.

CONCLUSIONS

Despite the recent extensive selective sweep within chromosome 6 of West African isolates, PfAP2-EXP2 of Ghanaian origin exhibits low nucleotide diversity and very low dN/dS consistent with purifying selection acting to maintain the function of an essential gene. The conservation of AP2-EXP2 is an important factor that makes it a potential drug target.

The protein domains of vertebrate species in which selection is more effective have greater intrinsic structural disorder

The nearly neutral theory of molecular evolution posits variation among species in the effectiveness of selection. In an idealized model, the census population size determines both this minimum magnitude of the selection coefficient required for deleterious variants to be reliably purged, and the amount of neutral diversity. Empirically, an 'effective population size' is often estimated from the amount of putatively neutral genetic diversity and is assumed to also capture a species' effectiveness of selection. A potentially more direct measure of the effectiveness of selection is the degree to which selection maintains preferred codons. However, past metrics that compare codon bias across species are confounded by among-species variation in %GC content and/or amino acid composition. Here, we propose a new Codon Adaptation Index of Species (CAIS), based on Kullback-Leibler divergence, that corrects for both confounders. We demonstrate the use of CAIS correlations, as well as the Effective Number of Codons, to show that the protein domains of more highly adapted vertebrate species evolve higher intrinsic structural disorder.

Computational exploration of protein structure dynamics and RNA structural consequences of PKD1 missense variants: implications in ADPKD pathogenesis

We analyzed the impact of nine previously identified missense PKD1 variants from our studies, including c.6928G > A p.G2310R, c.8809G > A p.E2937K, c.2899 T > C p.W967R, c.6284A > G p.D2095G, c.6644G > A p.R2215Q, c.7810G > A p.D2604N, c.11249G > C p.R3750P, c.1001C > T p.T334M, and c.3101A > G p.N1034S on RNA structures and PC1 protein structure dynamics utilizing computational tools. RNA structure analysis was done using short RNA snippets of 41 nucleotides with the variant position at the 21st nucleotide, ensuring 20 bases on both sides. The secondary structures of these RNA snippets were predicted using RNAstructure. Structural changes of the mutants compared to the wild type were analyzed using the MutaRNA webserver. Molecular dynamics (MD) simulation of PC1 wild-type and mutant protein regions were performed using GROMACS 2018 (GROMOS96 54a7 force field). Findings revealed that five variants including c.8809G > A (p.E2937K), c.11249G > C (p.R3750P), c.3101A > G (p.N1034S), c.6928G > A (p.G2310R), c.6644G > A (p.R2215Q) exhibited major alterations in RNA structures and thereby their interactions with other proteins or RNAs affecting protein structure dynamics. While certain variants have minimal impact on RNA conformations, their observed alterations in MD simulations indicate impact on protein structure dynamics highlighting the importance of evaluating the functional consequences of genetic variants by considering both RNA and protein levels. The study also emphasizes that each missense variant exerts a unique impact on RNA stability, and protein structure dynamics, potentially contributing to the heterogeneous clinical manifestations and progression observed in Autosomal Dominant Polycystic Kidney Disease (ADPKD) patients offering a novel perspective in this direction. Thus, the utility of studying the structure dynamics through computational tools can help in prioritizing the variants for their functional implications, understanding the molecular mechanisms underlying variability in ADPKD presentation and developing targeted therapeutic interventions.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-024-04057-9.

Functional characterization of 2,832 JAG1 variants supports reclassification for Alagille syndrome and improves guidance for clinical variant interpretation

Pathogenic variants in the JAG1 gene are a primary cause of the multi-system disorder Alagille syndrome. Although variant detection rates are high for this disease, there is uncertainty associated with the classification of missense variants that leads to reduced diagnostic yield. Consequently, up to 85% of reported JAG1 missense variants have uncertain or conflicting classifications. We generated a library of 2,832 JAG1 nucleotide variants within exons 1-7, a region with a high number of reported missense variants, and designed a high-throughput assay to measure JAG1 membrane expression, a requirement for normal function. After calibration using a set of 175 known or predicted pathogenic and benign variants included within the variant library, 486 variants were characterized as functionally abnormal (n = 277 abnormal and n = 209 likely abnormal), of which 439 (90.3%) were missense. We identified divergent membrane expression occurring at specific residues, indicating that loss of the wild-type residue itself does not drive pathogenicity, a finding supported by structural modeling data and with broad implications for clinical variant classification both for Alagille syndrome and globally across other disease genes. Of 144 uncertain variants reported in patients undergoing clinical or research testing, 27 had functionally abnormal membrane expression, and inclusion of our data resulted in the reclassification of 26 to likely pathogenic. Functional evidence augments the classification of genomic variants, reducing uncertainty and improving diagnostics. Inclusion of this repository of functional evidence during JAG1 variant reclassification will significantly affect resolution of variant pathogenicity, making a critical impact on the molecular diagnosis of Alagille syndrome.

Population Diversity at the Single-Cell Level

Population-scale single-cell genomics is a transformative approach for unraveling the intricate links between genetic and cellular variation. This approach is facilitated by cutting-edge experimental methodologies, including the development of high-throughput single-cell multiomics and advances in multiplexed environmental and genetic perturbations. Examining the effects of natural or synthetic genetic variants across cellular contexts provides insights into the mutual influence of genetics and the environment in shaping cellular heterogeneity. The development of computational methodologies further enables detailed quantitative analysis of molecular variation, offering an opportunity to examine the respective roles of stochastic, intercellular, and interindividual variation. Future opportunities lie in leveraging long-read sequencing, refining disease-relevant cellular models, and embracing predictive and generative machine learning models. These advancements hold the potential for a deeper understanding of the genetic architecture of human molecular traits, which in turn has important implications for understanding the genetic causes of human disease.

Severity of effect considerations regarding the use of mutation as a toxicological endpoint for risk assessment: A report from the 8th International Workshop on Genotoxicity Testing (IWGT)

Exposure levels without appreciable human health risk may be determined by dividing a point of departure on a dose-response curve (e.g., benchmark dose) by a composite adjustment factor (AF). An "effect severity" AF (ESAF) is employed in some regulatory contexts. An ESAF of 10 may be incorporated in the derivation of a health-based guidance value (HBGV) when a "severe" toxicological endpoint, such as teratogenicity, irreversible reproductive effects, neurotoxicity, or cancer was observed in the reference study. Although mutation data have been used historically for hazard identification, this endpoint is suitable for quantitative dose-response modeling and risk assessment. As part of the 8th International Workshops on Genotoxicity Testing, a sub-group of the Quantitative Analysis Work Group (WG) explored how the concept of effect severity could be applied to mutation. To approach this question, the WG reviewed the prevailing regulatory guidance on how an ESAF is incorporated into risk assessments, evaluated current knowledge of associations between germline or somatic mutation and severe disease risk, and mined available data on the fraction of human germline mutations expected to cause severe disease. Based on this review and given that mutations are irreversible and some cause severe human disease, in regulatory settings where an ESAF is used, a majority of the WG recommends applying an ESAF value between 2 and 10 when deriving a HBGV from mutation data. This recommendation may need to be revisited in the future if direct measurement of disease-causing mutations by error-corrected next generation sequencing clarifies selection of ESAF values.

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.

Homozygous EPRS1 missense variant causing hypomyelinating leukodystrophy-15 alters variant-distal mRNA m6A site accessibility

Hypomyelinating leukodystrophy (HLD) is an autosomal recessive disorder characterized by defective central nervous system myelination. Exome sequencing of two siblings with severe cognitive and motor impairment and progressive hypomyelination characteristic of HLD revealed homozygosity for a missense single-nucleotide variant (SNV) in EPRS1 (c.4444 C > A; p.Pro1482Thr), encoding glutamyl-prolyl-tRNA synthetase, consistent with HLD15. Patient lymphoblastoid cell lines express markedly reduced EPRS1 protein due to dual defects in nuclear export and cytoplasmic translation of variant EPRS1 mRNA. Variant mRNA exhibits reduced METTL3 methyltransferase-mediated writing of N6-methyladenosine (m6A) and reduced reading by YTHDC1 and YTHDF1/3 required for efficient mRNA nuclear export and translation, respectively. In contrast to current models, the variant does not alter the sequence of m6A target sites, but instead reduces their accessibility for modification. The defect was rescued by antisense morpholinos predicted to expose m6A sites on target EPRS1 mRNA, or by m6A modification of the mRNA by METTL3-dCas13b, a targeted RNA methylation editor. Our bioinformatic analysis predicts widespread occurrence of SNVs associated with human health and disease that similarly alter accessibility of distal mRNA m6A sites. These results reveal a new RNA-dependent etiologic mechanism by which SNVs can influence gene expression and disease, consequently generating opportunities for personalized, RNA-based therapeutics targeting these disorders.

Population-specific variations in KCNH2 predispose patients to delayed ventricular repolarization upon dihydroartemisinin-piperaquine therapy

Dihydroartemisinin-piperaquine is efficacious for the treatment of uncomplicated malaria and its use is increasing globally. Despite the positive results in fighting malaria, inhibition of the Kv11.1 channel (hERG; encoded by the KCNH2 gene) by piperaquine has raised concerns about cardiac safety. Whether genetic factors could modulate the risk of piperaquine-mediated QT prolongations remained unclear. Here, we first profiled the genetic landscape of KCNH2 variability using data from 141,614 individuals. Overall, we found 1,007 exonic variants distributed over the entire gene body, 555 of which were missense. By optimizing the gene-specific parametrization of 16 partly orthogonal computational algorithms, we developed a KCNH2-specific ensemble classifier that identified a total of 116 putatively deleterious missense variations. To evaluate the clinical relevance of KCNH2 variability, we then sequenced 293 Malian patients with uncomplicated malaria and identified 13 variations within the voltage sensing and pore domains of Kv11.1 that directly interact with channel blockers. Cross-referencing of genetic and electrocardiographic data before and after piperaquine exposure revealed that carriers of two common variants, rs1805121 and rs41314375, experienced significantly higher QT prolongations (ΔQTc of 41.8 ms and 61 ms, respectively, vs 14.4 ms in controls) with more than 50% of carriers having increases in QTc >30 ms. Furthermore, we identified three carriers of rare population-specific variations who experienced clinically relevant delayed ventricular repolarization. Combined, our results map population-scale genetic variability of KCNH2 and identify genetic biomarkers for piperaquine-induced QT prolongation that could help to flag at-risk patients and optimize efficacy and adherence to antimalarial therapy.

The protein domains of vertebrate species in which selection is more effective have greater intrinsic structural disorder

The nearly neutral theory of molecular evolution posits variation among species in the effectiveness of selection. In an idealized model, the census population size determines both this minimum magnitude of the selection coefficient required for deleterious variants to be reliably purged, and the amount of neutral diversity. Empirically, an "effective population size" is often estimated from the amount of putatively neutral genetic diversity and is assumed to also capture a species' effectiveness of selection. A potentially more direct measure of the effectiveness of selection is the degree to which selection maintains preferred codons. However, past metrics that compare codon bias across species are confounded by among-species variation in %GC content and/or amino acid composition. Here we propose a new Codon Adaptation Index of Species (CAIS), based on Kullback-Leibler divergence, that corrects for both confounders. We demonstrate the use of CAIS correlations, as well as the Effective Number of Codons, to show that the protein domains of more highly adapted vertebrate species evolve higher intrinsic structural disorder.

Identification of Synonymous Pathogenic Variants in Monogenic Disorders by Integrating Exome with Transcriptome Sequencing

Exome sequencing is becoming a first-tier clinical diagnostic test for Mendelian diseases, drastically reducing the time and cost of diagnostic odyssey and improving the diagnosis rate. Despite its success, exome sequencing faces practical challenges in assessing the pathogenicity of numerous intronic and synonymous variants, leaving a significant proportion of patients undiagnosed. In this study, a whole-blood transcriptome database was constructed that showed the expression profile of 2981 Online Mendelian Inheritance in Man disease genes in blood samples. Meanwhile, a workflow integrating exome sequencing, blood transcriptome sequencing, and in silico prediction tools to identify and validate splicing-altering intronic or synonymous variants was proposed. Following this pipeline, seven synonymous variants in eight patients were discovered. Of these, the functional evidence of c.981G>A (PIGN), c.1161A>G (ALPL), c.858G>A (ATP6AP2), and c.1011G>T (MTHFR) have not been reported previously. RNA sequencing validation confirmed that these variants induced aberrant splicing, expanding the disease-causing variant spectrum of these genes. Overall, this study shows the feasibility of combining multi-omics data to identify splicing-altering variants, especially the power of RNA sequencing. It also reveals that synonymous variants, which often are overlooked in standard diagnostic approaches, comprise an important portion of unresolved genetic diseases.

Functional genomics in inborn errors of immunity

Inborn errors of immunity (IEI) comprise a diverse spectrum of 485 disorders as recognized by the International Union of Immunological Societies Committee on Inborn Error of Immunity in 2022. While IEI are monogenic by definition, they illuminate various pathways involved in the pathogenesis of polygenic immune dysregulation as in autoimmune or autoinflammatory syndromes, or in more common infectious diseases that may not have a significant genetic basis. Rapid improvement in genomic technologies has been the main driver of the accelerated rate of discovery of IEI and has led to the development of innovative treatment strategies. In this review, we will explore various facets of IEI, delving into the distinctions between PIDD and PIRD. We will examine how Mendelian inheritance patterns contribute to these disorders and discuss advancements in functional genomics that aid in characterizing new IEI. Additionally, we will explore how emerging genomic tools help to characterize new IEI as well as how they are paving the way for innovative treatment approaches for managing and potentially curing these complex immune conditions.

Identification of a novel and plant height-independent QTL for coleoptile length in barley and validation of its effect using near isogenic lines

KEY MESSAGE

This study reported the identification and validation of novel QTL conferring coleoptile length in barley and predicted candidate genes underlying the largest effect QTL based on orthologous analysis and comparison of the whole genome assemblies for both parental genotypes of the mapping population. Coleoptile length (CL) is one of the most important agronomic traits in cereal crops due to its direct influence on the optimal depth for seed sowing which facilitates better seedling establishment. Varieties with longer coleoptiles are preferred in drought-prone areas where less moisture maintains at the top layer of the soil. Compared to wheat, genetic study on coleoptile length is limited in barley. Here, we reported a study on detecting the genomic regions associated with CL in barley by assessing a population consisting of 201 recombinant inbred lines. Four putative QTL conferring CL were consistently identified on chromosomes 1H, 5H, 6H, and 7H in each of the trials conducted. Of these QTL, the two located on chromosomes 5H and 6H (designated as Qcl.caf-5H and Qcl.caf-6H) are likely novel and Qcl.caf-5H showed the most significant effect explaining up to 30.9% of phenotypic variance with a LOD value of 15.1. To further validate the effect of this putative QTL, five pairs of near isogenic lines (NILs) were then developed and assessed. Analysis of the NILs showed an average difference of 21.0% in CL between the two isolines. Notably, none of the other assessed morphological characteristics showed consistent differences between the two isolines for each pair of the NILs. Candidate genes underlying the Qcl.caf-5H locus were also predicted by employing orthologous analysis and comparing the genome assemblies for both parental genotypes of the mapping population in the present study. Taken together, these findings expand our understanding on genetic basis of CL and will be indicative for further gene cloning and functional analysis underly this locus in barley.

Broken silence: 22,841 predicted deleterious synonymous variants identified in the human exome through computational analysis

Synonymous single nucleotide variants (sSNVs) do not alter the primary structure of a protein, thus it was previously accepted that they were neutral. Recently, several studies demonstrated their significance to a range of diseases. Still, variant prioritization strategies lack focus on sSNVs. Here, we identified 22,841 deleterious synonymous variants in 125,748 human exomes using two in silico predictors (SilVA and CADD). While 98.2% of synonymous variants are classified as neutral, 1.8% are predicted to be deleterious, yielding an average of 9.82 neutral and 0.18 deleterious sSNVs per exome. Further investigation of prediction features via Heterogeneous Ensemble Feature Selection revealed that impact on amino acid sequence and conservation carry the most weight for a deleterious prediction. Thirty nine detrimental sSNVs are not rare and are located on disease associated genes. Ten distinct putatively non-deleterious sSNVs are likely to be under positive selection in the North-Western European and East Asian populations. Taken together our analysis gives voice to the so-called silent mutations as we propose a robust framework for evaluating the deleteriousness of sSNVs in variant prioritization studies.

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.

Inherited Retinal Degeneration Caused by Dehydrodolichyl Diphosphate Synthase Mutation-Effect of an ALG6 Modifier Variant

Modern advances in disease genetics have uncovered numerous modifier genes that play a role in the severity of disease expression. One such class of genetic conditions is known as inherited retinal degenerations (IRDs), a collection of retinal degenerative disorders caused by mutations in over 300 genes. A single missense mutation (K42E) in the gene encoding the enzyme dehydrodolichyl diphosphate synthase (DHDDS), which is required for protein N-glycosylation in all cells and tissues, causes DHDDS-IRD (retinitis pigmentosa type 59 (RP59; OMIM #613861)). Apart from a retinal phenotype, however, DHDDS-IRD is surprisingly non-syndromic (i.e., without any systemic manifestations). To explore disease pathology, we selected five glycosylation-related genes for analysis that are suggested to have disease modifier variants. These genes encode glycosyltransferases (ALG6, ALG8), an ER resident protein (DDOST), a high-mannose oligosaccharyl transferase (MPDU1), and a protein N-glycosylation regulatory protein (TNKS). DNA samples from 11 confirmed DHDDS (K42E)-IRD patients were sequenced at the site of each candidate genetic modifier. Quantitative measures of retinal structure and function were performed across five decades of life by evaluating foveal photoreceptor thickness, visual acuity, foveal sensitivity, macular and extramacular rod sensitivity, and kinetic visual field extent. The ALG6 variant, (F304S), was correlated with greater macular cone disease severity and less peripheral rod disease severity. Thus, modifier gene polymorphisms may account for a significant portion of phenotypic variation observed in human genetic disease. However, the consequences of the polymorphisms may be counterintuitively complex in terms of rod and cone populations affected in different regions of the retina.

HECT-type ubiquitin ligase KAKTUS mediates the proteasome-dependent degradation of cyclin-dependent kinase inhibitor KRP2 during trichome morphogenesis in Arabidopsis

Trichome development is a fascinating model to elaborate the plant cell differentiation and growth processes. A wealth of information has pointed to the contributions of the components associated with cell cycle control and ubiquitin/26S proteasome system (UPS) to trichome morphogenesis, but how these two pathways are connected remains obscure. Here, we report that HECT‐type ubiquitin ligase KAKTUS (KAK) targets the cyclin‐dependent kinase (CDK) inhibitor KRP2 (for kip‐related protein 2) for proteasome‐dependent degradation during trichome branching in Arabidopsis. We show that over‐expression of KRP2 promotes trichome branching and endoreduplication which is similar to kak loss of function mutants. KAK directly interacts with KRP2 and mediates KRP2 degradation. Mutation of KAK results in the accumulation of steady‐state KRP2. Consistently, in kak pKRP2:KRP2‐GFP plants, the trichome branching is further induced compared with the single mutant. Taken together, our studies bridge the cell cycle control and UPS pathways during trichome development and underscore the importance of post‐translational control in epidermal differentiation.

Screening and verifying the mutations in the LDLR and APOB genes in a Chinese family with familial hypercholesterolemia

BACKGROUND

Familial hypercholesterolemia (FH) is an autosomal dominant genetic disorder. The primary objective of this study was to identify the major pathogenic mutations in a Chinese family with FH.

METHODS

Whole-genome sequencing (WGS) was used to identify variants of FH-related genes, including low-density lipoprotein receptor (LDLR), apolipoprotein B (APOB), and proprotein convertase subtilisin/kexin 9 (PCSK9). Bioinformatics software was used to predict signal peptides, transmembrane structures, and spatial construction information of the mutated sequences. Western blotting was performed on the mutant protein to determine the presence of the major structural domains of the LDLR. The PCSK9 and APOB genes were screened and analyzed. Moreover, the proband and his brother were treated with a PCSK9 inhibitor for 1 year, and the effect of the treatment on lipid levels was assessed.

RESULTS

WGS revealed two potentially pathogenic mutations in the LDLR gene. One was a novel mutation, c.497delinsGGATCCCCCAGCTGCATCCCCCAG (p. Ala166fs), and the other was a known pathogenic mutation, c.2054C>T (p. Pro685Leu). Bioinformatics prediction and in vitro experiments revealed that the novel mutation could not be expressed on the cell membrane. Numerous gene variants were identified in the APOB gene that may have a significant impact on the family members with FH. Thus, it is suggested that the severe manifestation of FH in the proband primarily resulted from the cumulative genetic effects of variants in both LDLR and APOB. However, a subsequent study indicated that treatment with a PCSK9 inhibitor (Evolocumab) did not significantly reduce the blood lipid levels in the proband or his brother.

CONCLUSIONS

The cumulative effect of LDLR and APOB variants was the primary cause of elevated blood lipid levels in this family. However, PCSK9 inhibitor therapy did not appear to be beneficial for the proband. This study emphasizes the importance of genetic testing in determining the most suitable treatment options for patients with FH.

Revisiting the Role of Genetic Variation in Adaptation

AbstractEvolutionary biologists have thought about the role of genetic variation during adaptation for a very long time-before we understood the organization of the genetic code, the provenance of genetic variation, and how such variation influenced the phenotypes on which natural selection acts. Half a century after the discovery of the structure of DNA and the unraveling of the genetic code, we have a rich understanding of these problems and the means to both delve deeper and widen our perspective across organisms and natural populations. The 2022 Vice Presidential Symposium of the American Society of Naturalists highlighted examples of recent insights into the role of genetic variation in adaptive processes, which are compiled in this special section. The work was conducted in different parts of the world, included theoretical and empirical studies with diverse organisms, and addressed distinct aspects of how genetic variation influences adaptation. In our introductory article to the special section, we discuss some important recent insights about the generation and maintenance of genetic variation, its impacts on phenotype and fitness, its fate in natural populations, and its role in driving adaptation. By placing the special section articles in the broader context of recent developments, we hope that this overview will also serve as a useful introduction to the field.

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

BACKGROUND

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

METHODS

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

RESULTS

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

CONCLUSION

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

Highly Abundant Proteins Are Highly Thermostable

Highly abundant proteins tend to evolve slowly (a trend called E-R anticorrelation), and a number of hypotheses have been proposed to explain this phenomenon. The misfolding avoidance hypothesis attributes the E-R anticorrelation to the abundance-dependent toxic effects of protein misfolding. To avoid these toxic effects, protein sequences (particularly those of highly expressed proteins) would be under selection to fold properly. One prediction of the misfolding avoidance hypothesis is that highly abundant proteins should exhibit high thermostability (i.e., a highly negative free energy of folding, ΔG). Thus far, only a handful of analyses have tested for a relationship between protein abundance and thermostability, producing contradictory results. These analyses have been limited by 1) the scarcity of ΔG data, 2) the fact that these data have been obtained by different laboratories and under different experimental conditions, 3) the problems associated with using proteins' melting energy (Tm) as a proxy for ΔG, and 4) the difficulty of controlling for potentially confounding variables. Here, we use computational methods to compare the free energy of folding of pairs of human-mouse orthologous proteins with different expression levels. Even though the effect size is limited, the most highly expressed ortholog is often the one with a more negative ΔG of folding, indicating that highly expressed proteins are often more thermostable.

Delayed-onset adenosine deaminase deficiency with a novel synonymous mutation and a case series from China

BACKGROUND

Adenosine deaminase (ADA) is a key enzyme in the purine salvage pathway. Genetic defects of the ADA gene can cause a subtype of severe combined immunodeficiency. To date, few Chinese cases have been reported.

METHODS

We retrospectively reviewed the medical records of patients diagnosed with ADA deficiency in Beijing Children's Hospital and summarized the previously published ADA deficiency cases from China in the literature.

RESULTS

Nine patients were identified with two novel mutations (W272X and Q202 =). Early-onset infection, thymic abnormalities and failure to thrive were the most common manifestations of Chinese ADA-deficient patients. The ADA genotype has a major effect on the clinical phenotype. Notably, a novel synonymous mutation (c.606G>A, p.Q202=) was identified in a delayed-onset patient, which affected pre-mRNA splicing leading to a frameshift and premature truncation of the protein. Furthermore, the patient showed γδT cells expansion with an increased effect or phenotype, which may be associated with the delayed onset of disease. In addition, we reported cerebral aneurysm and intracranial artery stenosis for the first time in ADA deficiency. Five patients died with a median age of four months, while two patients received stem cell transplantation and are alive.

CONCLUSIONS

This study described the first case series of Chinese ADA-deficient patients. Early-onset infection, thymic abnormalities and failure to thrive were the most common manifestations in our patients. We identified a synonymous mutation that affected pre-mRNA splicing in the ADA gene, which had never been reported in ADA deficiency. Furthermore, we reported cerebral aneurysm in a delayed-onset patient for the first time. Further study is warranted to investigate the underlying mechanisms.

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.

How soluble misfolded proteins bypass chaperones at the molecular level

Subpopulations of soluble, misfolded proteins can bypass chaperones within cells. The extent of this phenomenon and how it happens at the molecular level are unknown. Through a meta-analysis of the experimental literature we find that in all quantitative protein refolding studies there is always a subpopulation of soluble but misfolded protein that does not fold in the presence of one or more chaperones, and can take days or longer to do so. Thus, some misfolded subpopulations commonly bypass chaperones. Using multi-scale simulation models we observe that the misfolded structures that bypass various chaperones can do so because their structures are highly native like, leading to a situation where chaperones do not distinguish between the folded and near-native-misfolded states. More broadly, these results provide a mechanism by which long-time scale changes in protein structure and function can persist in cells because some misfolded states can bypass components of the proteostasis machinery.

Case report: Autosomal recessive type 3 Stickler syndrome caused by compound heterozygous mutations in COL11A2

Background: Stickler syndrome (SS) is a group of hereditary collagenopathies caused by a variety of collagen and non-collagen genes. Affected patients have characteristic manifestations involving ophthalmic, articular, craniofacial and auditory disorders. SS is classified into several subtypes according to clinical and molecular features. Type 3 SS is an ultra-rare disease, known as non-ocular SS or otospondylomegaepiphyseal dysplasia (OSMED) with only a few pathogenic COL11A2 variants reported to date. Case presentation: A 29-year-old Chinese male was referred to our hospital for hearing loss and multiple joint pain. He presented a phenotype highly suggestive of OSMED, including progressive sensorineural deafness, spondyloepiphyseal dysplasia with large epiphyses, platyspondyly, degenerative osteoarthritis, and sunken nasal bridge. We detected compound heterozygous mutations in COL11A2, both of which were predicted to be splicing mutations. One is synonymous mutation c.3774C>T (p.Gly1258Gly) supposed to be a splice site mutation, the other is a novel intron mutation c.4750 + 5 G>A, which is a highly conservative site across several species. We also present a review of the current known pathogenic mutation spectrum of COL11A2 in patients with type 3 SS. Conclusion: Both synonymous extonic and intronic variants are easily overlooked by whole-exome sequencing. For patients with clinical manifestations suspected of SS syndrome, next-generation whole-genome sequencing is necessary for precision diagnosis and genetic counseling.

Fine mapping of a Fusarium crown rot resistant locus on chromosome arm 6HL in barley by exploiting near isogenic lines, transcriptome profiling, and a large near isogenic line-derived population

KEY MESSAGE

This study reported validation and fine mapping of a Fusarium crown rot resistant locus on chromosome arm 6HL in barley using near isogenic lines, transcriptome sequences, and a large near isogenic line-derived population. Fusarium crown rot (FCR), caused by Fusarium pseudograminearum, is a chronic and serious disease affecting cereal production in semi-arid regions globally. The increasing prevalence of this disease in recent years is attributed to the widespread adoption of minimum tillage and stubble retention practices. In the study reported here, we generated eight pairs of near isogenic lines (NILs) targeting a putative QTL (Qcrs.caf-6H) conferring FCR resistance in barley. Assessing the NILs confirmed the large effect of this locus. Aimed to develop markers that can be reliably used in incorporating this resistant allele into breeding programs and identify candidate genes, transcriptomic analyses were conducted against three of the NIL pairs and a large NIL-derived population consisting of 1085 F7 recombinant inbred lines generated. By analyzing the transcriptomic data and the fine mapping population, Qcrs.caf-6H was delineated into an interval of 0.9 cM covering a physical distance of ~ 547 kb. Six markers co-segregating with this locus were developed. Based on differential gene expression and SNP variations between the two isolines among the three NIL pairs, candidate genes underlying the resistance at this locus were detected. These results would improve the efficiency of incorporating the targeted locus into barley breeding programs and facilitate the cloning of causal gene(s) responsible for the resistance.

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.

Most synonymous allelic variants in HIV tat are not silent

The genetic code has degenerate codons that produce no change in the translated protein sequence and are generally thought to be silent. However, some synonymous variants are clearly not silent. Herein, we questioned the frequency of non-silent synonymous variants. We tested how random synonymous variants in the HIV Tat transcription factor effect transcription of an LTR-GFP reporter. Our model system has the advantage of directly measuring the function of the gene in human cells. Approximately, 67% of synonymous variants in Tat were non-silent, either having reduced activity or were full loss-of-function alleles. Eight mutant codons had higher codon usage than wild type, accompanied by reduced transcriptional activity. These were clustered on a loop in the Tat structure. We conclude that most synonymous Tat variants are not silent in human cells, and 25% are associated with changes in codon usage, likely effecting protein folding.

Aberrant Cyclin D1 splicing in cancer: from molecular mechanism to therapeutic modulation

Cyclin D1 (CCND1), a crucial mediator of cell cycle progression, possesses many mutation types with different mutation frequencies in human cancers. The G870A mutation is the most common mutation in CCND1, which produces two isoforms: full-length CCND1a and divergent C-terminal CCND1b. The dysregulation of the CCND1 isoforms is associated with multiple human cancers. Exploring the molecular mechanism of CCND1 isoforms has offer new insight for cancer treatment. On this basis, the alterations of CCND1 gene are described, including amplification, overexpression, and mutation, especially the G870A mutation. Subsequently, we review the characteristics of CCND1 isoforms caused by G870A mutation. Additionally, we summarize cis-regulatory elements, trans-acting factors, and the splice mutation involved in splicing regulation of CCND1. Furthermore, we highlight the function of CCND1 isoforms in cell cycle, invasion, and metastasis in cancers. Importantly, the clinical role of CCND1 isoforms is also discussed, particularly concerning prognosis, chemotherapy, and radiotherapy. Last, emphasis is given to the corrective strategies that modulate the cancerous CCND1 isoforms. Thus, it is highlighting significance of aberrant isoforms of CCND1 as targets for cancer therapy.

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

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.

Association analysis between the TLR9 gene polymorphism rs352140 and type 1 diabetes

Background

To a great extent, genetic factors contribute to the susceptibility to type 1 diabetes (T1D) development, and by triggering immune imbalance, Toll-like receptor (TLR) 9 is involved in the development of T1D. However, there is a lack of evidence supporting a genetic association between polymorphisms in the TLR9 gene and T1D.

Methods

In total, 1513 individuals, including T1D patients (n=738) and healthy control individuals (n=775), from the Han Chinese population were recruited for an association analysis of the rs352140 polymorphism of the TLR9 gene and T1D. rs352140 was genotyped by MassARRAY. The allele and genotype distributions of rs352140 in the T1D and healthy groups and those in different T1D subgroups were analyzed by the chi-squared test and binary logistic regression model. The chi-square test and Kruskal-Wallis H test were performed to explore the association between genotype and phenotype in T1D patients.

Results

The allele and genotype distributions of rs352140 were significantly different in T1D patients and healthy control individuals (p=0.019, p=0.035). Specifically, the T allele and TT genotype of rs352140 conferred a higher risk of T1D (OR=1.194, 95% CI=1.029-1.385, p=0.019, OR=1.535, 95% CI=1.108-2.126, p=0.010). The allele and genotype distributions of rs352140 were not significantly different between childhood-onset and adult-onset T1D and between T1D with a single islet autoantibody and T1D with multiple islet autoantibodies (p=0.603, p=0.743). rs352140 was associated with T1D susceptibility according to the recessive and additive models (p=0.015, p=0.019) but was not associated with T1D susceptibility in the dominant and overdominant models (p=0.117, p=0.928). Moreover, genotype-phenotype association analysis showed that the TT genotype of rs352140 was associated with higher fasting C-peptide levels (p=0.017).

Conclusion

In the Han Chinese population, the TLR9 polymorphism rs352140 is associated with T1D and is a risk factor for susceptibility to T1D.

Selection of Diagnostically Significant Regions of the SLC26A4 Gene Involved in Hearing Loss

Screening pathogenic variants in the SLC26A4 gene is an important part of molecular genetic testing for hearing loss (HL) since they are one of the common causes of hereditary HL in many populations. However, a large size of the SLC26A4 gene (20 coding exons) predetermines the difficulties of its complete mutational analysis, especially in large samples of patients. In addition, the regional or ethno-specific prevalence of SLC26A4 pathogenic variants has not yet been fully elucidated, except variants c.919-2A>G and c.2168A>G (p.His723Arg), which have been proven to be most common in Asian populations. We explored the distribution of currently known pathogenic and likely pathogenic (PLP) variants across the SLC26A4 gene sequence presented in the Deafness Variation Database for the selection of potential diagnostically important parts of this gene. As a result of this bioinformatic analysis, we found that molecular testing ten SLC26A4 exons (4, 6, 10, 11, 13−17 and 19) with flanking intronic regions can provide a diagnostic rate of 61.9% for all PLP variants in the SLC26A4 gene. The primary sequencing of these SLC26A4 regions may be applied as an initial effective diagnostic testing in samples of patients of unknown ethnicity or as a subsequent step after the targeted testing of already-known ethno- or region-specific pathogenic SLC26A4 variants.

The Enigmatic Etiology of Oculo-Auriculo-Vertebral Spectrum (OAVS): An Exploratory Gene Variant Interaction Approach in Candidate Genes

The clinical diagnosis of oculo-auriculo-vertebral spectrum (OAVS) is established when microtia is present in association with hemifacial hypoplasia (HH) and/or ocular, vertebral, and/or renal malformations. Genetic and non-genetic factors have been associated with microtia/OAVS. Although the etiology remains unknown in most patients, some cases may have an autosomal dominant, autosomal recessive, or multifactorial inheritance. Among the possible genetic factors, gene−gene interactions may play important roles in the etiology of complex diseases, but the literature lacks related reports in OAVS patients. Therefore, we performed a gene−variant interaction analysis within five microtia/OAVS candidate genes (HOXA2, TCOF1, SALL1, EYA1 and TBX1) in 49 unrelated OAVS Mexican patients (25 familial and 24 sporadic cases). A statistically significant intergenic interaction (p-value < 0.001) was identified between variants p.(Pro1099Arg) TCOF1 (rs1136103) and p.(Leu858=) SALL1 (rs1965024). This intergenic interaction may suggest that the products of these genes could participate in pathways related to craniofacial alterations, such as the retinoic acid (RA) pathway. The absence of clearly pathogenic variants in any of the analyzed genes does not support a monogenic etiology for microtia/OAVS involving these genes in our patients. Our findings could suggest that in addition to high-throughput genomic approaches, future gene−gene interaction analyses could contribute to improving our understanding of the etiology of microtia/OAVS.

Signals of adaptation to agricultural stress in the genomes of two European bumblebees

Human-induced environmental impacts on wildlife are widespread, causing major biodiversity losses. One major threat is agricultural intensification, typically characterised by large areas of monoculture, mechanical tillage, and the use of agrochemicals. Intensification leads to the fragmentation and loss of natural habitats, native vegetation, and nesting and breeding sites. Understanding the adaptability of insects to these changing environmental conditions is critical to predicting their survival. Bumblebees, key pollinators of wild and cultivated plants, are used as model species to assess insect adaptation to anthropogenic stressors. We investigated the effects of agricultural pressures on two common European bumblebees, Bombus pascuorum and B. lapidarius. Restriction-site Associated DNA Sequencing was used to identify loci under selective pressure across agricultural-natural gradients over 97 locations in Europe. 191 unique loci in B. pascuorum and 260 in B. lapidarius were identified as under selective pressure, and associated with agricultural stressors. Further investigation suggested several candidate proteins including several neurodevelopment, muscle, and detoxification proteins, but these have yet to be validated. These results provide insights into agriculture as a stressor for bumblebees, and signal for conservation action in light of ongoing anthropogenic changes.

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.