Human Gene Mutation Database (HGMD): 2003 update

An expansion of phenotype: novel homozygous variant in the MED17 identified in patients with progressive microcephaly and global developmental delay

Global developmental delay (GDD) is a lifelong disability that affects 1-3% of the population around the globe. It is phenotypically variable and highly heterogeneous in terms of the underlying genetics. Patients with GDD are intellectually disabled (ID) manifesting cognitive impairment and deficient adaptive behavior. Here, we investigated a two-looped consanguineous family segregating severe ID, seizure, and progressive microcephaly. Magnetic resonance imaging (MRI) of the brain showed mild brain atrophy and myelination defect. Whole exome sequencing (WES) was performed on the DNA samples of two patients and a novel homozygous missense variant (Chr11:g0.93528085; NM_004268.5_c.871T > C; p. Trp291Gly) was identified in the MED17 gene. Sanger sequencing revealed that the identified variant is heterozygous in both parents and healthy siblings. This variant is conserved among different species, causes a non-conserved amino acid change, and is predicted deleterious by various in silico tools. The variant is not reported in population variant databases. MED17 (OMIM: 613668) encodes for the mediator of RNA polymerase II transcription complex subunit 17. Structure modeling of MED17 protein revealed that Trp291 is involved in different inter-helical interactions, providing structural stability. Replacement of Trp291Gly, a less hydrophobic amino acid loses the inter-helical interaction leading to a perturb variant of MED17 protein.

GWAS meta-analysis of intrahepatic cholestasis of pregnancy implicates multiple hepatic genes and regulatory elements

Intrahepatic cholestasis of pregnancy (ICP) is a pregnancy-specific liver disorder affecting 0.5-2% of pregnancies. The majority of cases present in the third trimester with pruritus, elevated serum bile acids and abnormal serum liver tests. ICP is associated with an increased risk of adverse outcomes, including spontaneous preterm birth and stillbirth. Whilst rare mutations affecting hepatobiliary transporters contribute to the aetiology of ICP, the role of common genetic variation in ICP has not been systematically characterised to date. Here, we perform genome-wide association studies (GWAS) and meta-analyses for ICP across three studies including 1138 cases and 153,642 controls. Eleven loci achieve genome-wide significance and have been further investigated and fine-mapped using functional genomics approaches. Our results pinpoint common sequence variation in liver-enriched genes and liver-specific cis-regulatory elements as contributing mechanisms to ICP susceptibility.

A homozygous variant in the GPIHBP1 gene in a child with severe hypertriglyceridemia and a systematic literature review

Background: Due to nonspecific symptoms, rare dyslipidaemias are frequently misdiagnosed, overlooked, and undertreated, leading to increased risk for severe cardiovascular disease, pancreatitis and/or multiple organ failures before diagnosis. Better guidelines for the recognition and early diagnosis of rare dyslipidaemias are urgently required.

Methods: Genomic DNA was isolated from blood samples of a Pakistani paediatric patient with hypertriglyceridemia, and from his parents and siblings. Next-generation sequencing (NGS) was performed, and an expanded dyslipidaemia panel was employed for genetic analysis.

Results: The NGS revealed the presence of a homozygous missense pathogenic variant c.230G>A (NM_178172.6) in exon 3 of the GPIHBP1 (glycosylphosphatidylinositol-anchored high-density lipoprotein-binding protein 1) gene resulting in amino acid change p.Cys77Tyr (NP_835466.2). The patient was 5.5 years old at the time of genetic diagnosis. The maximal total cholesterol and triglyceride levels were measured at the age of 10 months (850.7 mg/dl, 22.0 mmol/L and 5,137 mg/dl, 58.0 mmol/L, respectively). The patient had cholesterol deposits at the hard palate, eruptive xanthomas, lethargy, poor appetite, and mild splenomegaly. Both parents and sister were heterozygous for the familial variant in the GPIHBP1 gene. Moreover, in the systematic review, we present 62 patients with pathogenic variants in the GPIHBP1 gene and clinical findings, associated with hyperlipoproteinemia.

Conclusion: In a child with severe hypertriglyceridemia, we identified a pathogenic variant in the GPIHBP1 gene causing hyperlipoproteinemia (type 1D). In cases of severe elevations of plasma cholesterol and/or triglycerides genetic testing for rare dyslipidaemias should be performed as soon as possible for optimal therapy and patient management.

Translation termination codons in protein synthesis and disease

Fidelity of protein synthesis, a process shaped by several mechanisms involving specialized ribosome regions and external factors, ensures the precise reading of sense as well as stop codons (UGA, UAG, UAA), which are usually localized at the 3' of mRNA and drive the release of the polypeptide chain. However, either natural (NTCs) or premature (PTCs) termination codons, the latter arising from nucleotide changes, can undergo a recoding process named ribosome or translational readthrough, which insert specific amino acids (NTCs) or subset(s) depending on the stop codon type (PTCs). This process is particularly relevant for nonsense mutations, a relatively frequent cause of genetic disorders, which impair gene expression at different levels by potentially leading to mRNA degradation and/or synthesis of truncated proteins. As a matter of fact, many efforts have been made to develop efficient and safe readthrough-inducing compounds, which have been challenged in several models of human disease to provide with a therapy. In this view, the dissection of the molecular determinants shaping the outcome of readthrough, namely nucleotide and protein contexts as well as their interplay and impact on protein structure/function, is crucial to identify responsive nonsense mutations resulting in functional full-length proteins. The interpretation of experimental and mechanistic findings is also important to define a possibly clear picture of potential readthrough-favorable features useful to achieve rescue profiles compatible with therapeutic thresholds typical of each targeted disorder, which is of primary importance for the potential translatability of readthrough into a personalized and mutation-specific, and thus patient-oriented, therapeutic strategy.

A Comprehensive Evaluation of the Performance of Prediction Algorithms on Clinically Relevant Missense Variants

The rapid integration of genomic technologies in clinical diagnostics has resulted in the detection of a multitude of missense variants whose clinical significance is often unknown. As a result, a plethora of computational tools have been developed to facilitate variant interpretation. However, choosing an appropriate software from such a broad range of tools can be challenging; therefore, systematic benchmarking with high-quality, independent datasets is critical. Using three independent benchmarking datasets compiled from the ClinVar database, we evaluated the performance of ten widely used prediction algorithms with missense variants from 21 clinically relevant genes, including BRCA1 and BRCA2. A fourth dataset consisting of 1053 missense variants was also used to investigate the impact of type 1 circularity on their performance. The performance of the prediction algorithms varied widely across datasets. Based on Matthews Correlation Coefficient and Area Under the Curve, SNPs&GO and PMut consistently displayed an overall above-average performance across the datasets. Most of the tools demonstrated greater sensitivity and negative predictive values at the expense of lower specificity and positive predictive values. We also demonstrated that type 1 circularity significantly impacts the performance of these tools and, if not accounted for, may confound the selection of the best performing algorithms.

Molecular characterization of triple-negative myeloproliferative neoplasms by next-generation sequencing

The role of next-generation sequencing (NGS) in identifying mutations in the driver, epigenetic regulator, RNA splicing, and signaling pathway genes in myeloproliferative neoplasms (MPNs) has contributed substantially to our understanding of the disease pathogenesis as well as disease evolution. NGS aids in determining the clonal nature of the disease in a subset of these disorders where mutations in the driver genes are not detected. There is a paucity of real-world data on the utility of this test in the characterization of triple-negative myeloproliferative neoplasms (TN-MPN). In this study, 46 samples of TN-MPN (essential thrombocythemia (ET) = 17; primary myelofibrosis (PMF) = 23; & myeloproliferative neoplasm unclassified (MPN-u) = 6) were screened for markers of clonality using targeted NGS. Among these, 25 (54.3%) patients had mutations that would help determine the clonal nature of the disease. Eight of the 17 TN-ET (47%) and 13 of the 23 TN-PMF (56.5%) patients had noncanonical mutations in the driver genes and mutations in the genes involved in epigenetic regulation. Identification of mutations categorized as high molecular markers (HMR) in 2 patients helped classify them as PMF with high risk according to the MIPSS 70 scoring system. A novel mutation in the MPIG6B (C6orf25) gene associated with childhood myelofibrosis was detected in a 14-year-old girl. The presence of clonal hematopoiesis could be confirmed in four of the six MPN-u patients in this cohort. This study demonstrates the utility of NGS in improving the characterization of TN-MPN by establishing clonality and detecting noncanonical mutations in driver genes, thereby aiding in clinical decision-making.

An Extended C-Terminus, the Possible Culprit for Differential Regulation of 5-Aminolevulinate Synthase Isoforms

5-Aminolevulinate synthase (ALAS; E.C. 2.3.1.37) is a pyridoxal 5′-phosphate (PLP)-dependent enzyme that catalyzes the key regulatory step of porphyrin biosynthesis in metazoa, fungi, and α-proteobacteria. ALAS is evolutionarily related to transaminases and is therefore classified as a fold type I PLP-dependent enzyme. As an enzyme controlling the key committed and rate-determining step of a crucial biochemical pathway ALAS is ideally positioned to be subject to allosteric feedback inhibition. Extensive kinetic and mutational studies demonstrated that the overall enzyme reaction is limited by subtle conformational changes of a hairpin loop gating the active site. These findings, coupled with structural information, facilitated early prediction of allosteric regulation of activity via an extended C-terminal tail unique to eukaryotic forms of the enzyme. This prediction was subsequently supported by the discoveries that mutations in the extended C-terminus of the erythroid ALAS isoform (ALAS2) cause a metabolic disorder known as X-linked protoporphyria not by diminishing activity, but by enhancing it. Furthermore, kinetic, structural, and molecular modeling studies demonstrated that the extended C-terminal tail controls the catalytic rate by modulating conformational flexibility of the active site loop. However, the precise identity of any such molecule remains to be defined. Here we discuss the most plausible allosteric regulators of ALAS activity based on divergences in AlphaFold-predicted ALAS structures and suggest how the mystery of the mechanism whereby the extended C-terminus of mammalian ALASs allosterically controls the rate of porphyrin biosynthesis might be unraveled.

Somatic RET Indels in Sporadic Medullary Thyroid Cancer: Prevalence and Response to Selpercatinib

CONTEXT

Although the majority of RET alterations are single nucleotide variants (SNV), small deletions and/or insertions have been reported at variable prevalence. No information about the efficacy of RET-specific inhibitors in patients harboring RET indels has been provided.

OBJECTIVE

We present an update on the prevalence of RET indels in medullary thyroid cancer (MTC) and describe the efficacy of selpercatinib in patients with advanced MTC with RET indels.

METHODS

The MTC tissues of 287 patients were analyzed using an Ion S5 targeted sequencing. The functional role of the reported indels have been evaluated by MutationTaster. Clinical and pathological data of MTC patients harboring a RET indel were collected and analyzed. Two patients with a RET indel were treated with selpercatinib.

RESULTS

Among 178 RET-positive cases, 147 (82.6%) harbored a SNV and 31 (17.4%) a RET in-frame indel. Nine indels were not previously reported and were found to be disease causing by MutationTaster. Patients harboring an indel were found to have an aggressive disease and 2 of them were treated with selpercatinib, experiencing a good response to the treatment.

CONCLUSION

These data show that RET indels are not infrequent and correlate with an aggressive disease. Two RET indel-positive patients showed a partial response to the treatment with a highly selective RET inhibitor; thus, these RET indels can be considered actionable mutations. In order to not miss these alterations, the analysis of the full gene is recommended.

Template-independent genome editing in the Pcdh15av-3j mouse, a model of human DFNB23 nonsyndromic deafness

Although frameshift mutations lead to 22% of inherited Mendelian disorders in humans, there is no efficient in vivo gene therapy strategy available to date, particularly in nondividing cells. Here, we show that nonhomologous end-joining (NHEJ)-mediated nonrandom editing profiles compensate the frameshift mutation in the Pcdh15 gene and restore the lost mechanotransduction function in postmitotic hair cells of Pcdh15^av-3J mice, an animal model of human nonsyndromic deafness DFNB23. Identified by an ex vivo evaluation system in cultured cochlear explants, the selected guide RNA restores reading frame in approximately 50% of indel products and recovers mechanotransduction in more than 70% of targeted hair cells. In vivo treatment shows that half of the animals gain improvements in auditory responses, and balance function is restored in the majority of injected mutant mice. These results demonstrate that NHEJ-mediated reading-frame restoration is a simple and efficient strategy in postmitotic systems.

Loss-of-function, gain-of-function and dominant-negative mutations have profoundly different effects on protein structure

Most known pathogenic mutations occur in protein-coding regions of DNA and change the way proteins are made. Taking protein structure into account has therefore provided great insight into the molecular mechanisms underlying human genetic disease. While there has been much focus on how mutations can disrupt protein structure and thus cause a loss of function (LOF), alternative mechanisms, specifically dominant-negative (DN) and gain-of-function (GOF) effects, are less understood. Here, we investigate the protein-level effects of pathogenic missense mutations associated with different molecular mechanisms. We observe striking differences between recessive vs dominant, and LOF vs non-LOF mutations, with dominant, non-LOF disease mutations having much milder effects on protein structure, and DN mutations being highly enriched at protein interfaces. We also find that nearly all computational variant effect predictors, even those based solely on sequence conservation, underperform on non-LOF mutations. However, we do show that non-LOF mutations could potentially be identified by their tendency to cluster in three-dimensional space. Overall, our work suggests that many pathogenic mutations that act via DN and GOF mechanisms are likely being missed by current variant prioritisation strategies, but that there is considerable scope to improve computational predictions through consideration of molecular disease mechanisms.

Most known pathogenic mutations occur in protein-coding regions of DNA and change the way proteins are made. Here the authors analyse the locations of thousands of human disease mutations and their predicted effects on protein structure and show that,while loss-of-function mutations tend to be highly disruptive, non-loss-of-function mutations are in general much milder at a protein structural level.

Glutaric aciduria type 1: Diagnosis, clinical features and long‐term outcome in a large cohort of 34 Irish patients

Glutaric aciduria type 1 (GA1) is a rare neurometabolic disorder that can lead to encephalopathic crises and severe dystonic movement disorders. Adherence to strict dietary restriction, in particular a diet low in lysine, carnitine supplementation and emergency treatment in pre‐symptomatic patients diagnosed by high‐risk screen (HRS) or newborn screen (NBS) leads to a favourable outcome. We present biochemical and clinical characteristics and long‐term outcome data of 34 Irish patients with GA1 aged 1–40 years. Sixteen patients were diagnosed clinically, and 17 patients by HRS, prior to introduction of NBS for GA1 in the Republic of Ireland in 2018. One patient was diagnosed by NBS. Clinical diagnosis was at a median of 1 year (range 1 month to 8 years) and by HRS was at a median of 4 days (range 3 days to 11 years). 14/18 (77.8%) diagnosed by HRS or NBS had neither clinical manifestations nor radiological features of GA1, or had radiological features only, compared to 0/16 (0%) diagnosed clinically (p < 0.001). Patients diagnosed clinically who survived to school‐age were more likely to have significant cerebral palsy and dystonia (7/11; 63.6% vs. 0/13; 0%, p < 0.001). They were less likely to be in mainstream school versus the HRS group (5/10; 50% vs. 12/13; 92.3%; p = 0.012). Clinical events occurring after 6 years of age were unusual, but included spastic diplegia, thalamic haemorrhage, Chiari malformation, pituitary hormone deficiency and epilepsy. The exact aetiology of these events is unclear.

A RAS-Independent Biomarker Panel to Reliably Predict Response to MEK Inhibition in Colorectal Cancer

Simple Summary

Today, clinical management for the majority of cancer patients is still based on a “one-size-fits-all” approach. To improve the outcomes in the era of personalized medicine, it is essential to stratify patients based on established and novel biomarkers. In the present study, we investigated a SMAD4 loss-of-function mutation, which is associated with chemoresistance and decreased overall survival in colorectal cancer (CRC). To investigate the molecular mechanism behind the impact on drug response, we used CRISPR technology on patient-derived organoid models (PDOs) of CRC. We showed that PDOs with loss-of-function SMAD4 mutations are sensitive to MEK-inhibitors. Using a novel four-gene signature reliably predicts sensitivity towards MEK-inhibitors, regardless of the RAS and BRAF status. The present study is a significant step towards personalized cancer therapy by identifying a new biomarker.

Abstract

Background: In colorectal cancer (CRC), mutations of genes associated with the TGF-β/BMP signaling pathway, particularly affecting SMAD4, are known to correlate with decreased overall survival and it is assumed that this signaling axis plays a key role in chemoresistance. Methods: Using CRISPR technology on syngeneic patient-derived organoids (PDOs), we investigated the role of a loss-of-function of SMAD4 in sensitivity to MEK-inhibitors. CRISPR-engineered SMAD4^R361H PDOs were subjected to drug screening, RNA-Sequencing, and multiplex protein profiling (DigiWest^®). Initial observations were validated on an additional set of 62 PDOs with known mutational status. Results: We show that loss-of-function of SMAD4 renders PDOs sensitive to MEK-inhibitors. Multiomics analyses indicate that disruption of the BMP branch within the TGF-β/BMP pathway is the pivotal mechanism of increased drug sensitivity. Further investigation led to the identification of the SFAB-signature (SMAD4, FBXW7, ARID1A, or BMPR2), coherently predicting sensitivity towards MEK-inhibitors, independent of both RAS and BRAF status. Conclusion: We identified a novel mutational signature that reliably predicts sensitivity towards MEK-inhibitors, regardless of the RAS and BRAF status. This finding poses a significant step towards better-tailored cancer therapies guided by the use of molecular biomarkers.

Congenital insensitivity to pain: a novel mutation affecting a U12-type intron causes multiple aberrant splicing of SCN9A

ABSTRACT

Mutations in the alpha subunit of voltage-gated sodium channel 1.7 (NaV1.7), encoded by SCN9A gene, play an important role in the regulation of nociception and can lead to a wide range of clinical outcomes, ranging from extreme pain syndromes to congenital inability to experience pain. To expand the phenotypic and genotypic spectrum of SCN9A-related channelopathies, we describe the proband, a daughter born from consanguineous parents, who had pain insensitivity, diminished temperature sensation, foot burns, and severe loss of nociceptive nerve fibers in the epidermis. Next-generation sequencing of SCN9A (NM_002977.3) revealed a novel homozygous substitution (c.377+7T>G) in the donor splice site of intron 3. As the RNA functional testing is challenging, the in silico analysis is the first approach to predict possible alterations. In this case, the computational analysis was unable to identify the splicing consensus and could not provide any prediction for splicing defects. The affected intron indeed belongs to the U12 type, a family of introns characterised by noncanonical consensus at the splice sites, accounting only for 0.35% of all human introns, and is not included in most of the training sets for splicing prediction. A functional study on proband RNA showed different aberrant transcripts, where exon 3 was missing and an intron fragment was included. A quantification study using real-time polymerase chain reaction showed a significant reduction of the NaV1.7 canonical transcript. Collectively, these data widen the spectrum of SCN9A-related insensitivity to pain by describing a mutation causing NaV1.7 deficiency, underlying the nociceptor dysfunction, and highlight the importance of molecular investigation of U12 introns' mutations despite the silent prediction.

Polarized trafficking and copper transport activity of ATP7B: a mutational approach to establish genotype-phenotype correlation in Wilson disease

Mutation in ATP7B gene causes Wilson disease (WD) that is characterized by severe hepatic and neurological symptoms. ATP7B localizes at the trans-Golgi Network (TGN) transporting copper to copper-dependent enzymes and traffics in apically targeted vesicles upon intracellular copper elevation. To decode the cellular underpinnings of WD manifestation we investigated copper-responsive polarized trafficking and copper transport activity of fifteen WD causing point mutations in ATP7B. Amino-terminal mutations Gly85Val, Leu168Pro and Gly591Asp displayed TGN and sub-apical localization whereas, Leu492Ser mislocalized at the basolateral region. The actuator domain mutation Gly875Arg shows retention in the endoplasmic reticulum (ER), Ala874Val and Leu795Phe show partial targeting to TGN and post-Golgi vesicles. The Nucleotide-Binding Domain mutations His1069Gln and Leu1083Phe also display impaired targeting. The C-terminal mutations Leu1373Pro/Arg is arrested at ER but Ser1423Asn shows TGN localization. Transmembrane mutant Arg778Leu resides in ER and TGN while Arg969Gln is exclusively ER localized. Cellular Cu level does not alter the targeting of any of the studied mutations. Mutants that traffic to TGN exhibits biosynthetic function. Finally, we correlated cellular phenotypes with the clinical manifestation of the two most prevalent mutations; the early onset and more aggressive WD caused by Arg778Leu and the milder form of WD caused by mutation His1069Gln. This article is protected by copyright. All rights reserved.

Retrospective Natural History Study of RPGR-Related Cone- and Cone-Rod Dystrophies While Expanding the Mutation Spectrum of the Disease

Variants in the X-linked retinitis pigmentosa GTPase regulator gene (RPGR) and, specifically, in its retinal opening reading frame-15 isoform (RPGRORF15) may cause rod-cone (RCD), cone, and cone-rod dystrophies (CDs and CRDs). While RPGR-related RCDs have been frequently evaluated, the characteristics and progression of RPGR-related CD/CRDs are largely unknown. Therefore, the goal of our work was to perform genotype-phenotype correlations specifically in RPGRORF15-related CD/CRDs. This retrospective longitudinal study included 34 index patients and two affected relatives with a molecular diagnosis of RPGR-related CD/CRDs. Patients were recruited at the "Quinze-Vingts" Hospital, Paris, France and screened for mutations in RPGRORF15 at the Institut de la Vision, Paris, France. We identified 29 distinct variants, of which 27 were truncating. All were located in the 3' half of the RPGRORF15 transcript. Twenty of them were novel. Fifteen subjects were affected by CD, the remaining had CRD. When analyzing the longitudinal data, a progressive decline in visual acuity (VA) was noted, with more than 60% of the patients reaching VA ≥ 1 LogMar in the best eye after the fifth decade of life. To our knowledge, this is the largest described study of a cohort of CD/CRD patients affected by RPGRORF15 variants. Longitudinal data showed a rapidly progressive disease, possibly locating an optimal window of intervention for future therapies in younger ages.

Whole-Exome Sequencing Identified Genes Responsible for Thoracic Aortic Aneurysms and Dissections in three Chinese Families

Thoracic aortic aneurysms and dissections are precarious conditions that often cannot be diagnosed with fatal outcomes. Over the last few years, pathogenic variants in numerous genes have been identified that predispose to heritable presentations of TAAD. An evidence-based strategy for the selection of genes to test in familial TAAD helps inform family screening and intervention to prevent life-threatening events. Using whole-exome sequencing, four members of three unrelated families clinically diagnosed with TAAD were used to identify the genetic origin of the disorder. Variant evaluation was carried out to detect the pathogenic mutation. Our studies suggest that mutations of COL3A1 and ACTA2 are responsible for familial TAAD. In addition, we highlight FBLN5, FBN1, SLC2A10, FBN2, and NOTCH1 as candidate genes. Future studies of crosstalk among the pathways may provide us a step toward understanding the pathogenic mechanism. This finding indicates the necessity of obtaining family medical history and screening of extended relatives of patients with TAAD for the early identification and treatment of TAAD.

Quantitative prediction of variant effects on alternative splicing in MAPT using endogenous pre-messenger RNA structure probing

Splicing is highly regulated and is modulated by numerous factors. Quantitative predictions for how a mutation will affect precursor mRNA (pre-mRNA) structure and downstream function are particularly challenging. Here, we use a novel chemical probing strategy to visualize endogenous precursor and mature MAPT mRNA structures in cells. We used these data to estimate Boltzmann suboptimal structural ensembles, which were then analyzed to predict consequences of mutations on pre-mRNA structure. Further analysis of recent cryo-EM structures of the spliceosome at different stages of the splicing cycle revealed that the footprint of the B^act complex with pre-mRNA best predicted alternative splicing outcomes for exon 10 inclusion of the alternatively spliced MAPT gene, achieving 74% accuracy. We further developed a β-regression weighting framework that incorporates splice site strength, RNA structure, and exonic/intronic splicing regulatory elements capable of predicting, with 90% accuracy, the effects of 47 known and 6 newly discovered mutations on inclusion of exon 10 of MAPT. This combined experimental and computational framework represents a path forward for accurate prediction of splicing-related disease-causing variants.

Rare and population-specific functional variation across pig lines

BACKGROUND

It is expected that functional, mainly missense and loss-of-function (LOF), and regulatory variants are responsible for most phenotypic differences between breeds and genetic lines of livestock species that have undergone diverse selection histories. However, there is still limited knowledge about the existing missense and LOF variation in commercial livestock populations, in particular regarding population-specific variation and how it can affect applications such as across-breed genomic prediction.

METHODS

We re-sequenced the whole genome of 7848 individuals from nine commercial pig lines (average sequencing coverage: 4.1×) and imputed whole-genome genotypes for 440,610 pedigree-related individuals. The called variants were categorized according to predicted functional annotation (from LOF to intergenic) and prevalence level (number of lines in which the variant segregated; from private to widespread). Variants in each category were examined in terms of their distribution along the genome, alternative allele frequency, per-site Wright's fixation index (F_ST), individual load, and association to production traits.

RESULTS

Of the 46 million called variants, 28% were private (called in only one line) and 21% were widespread (called in all nine lines). Genomic regions with a low recombination rate were enriched with private variants. Low-prevalence variants (called in one or a few lines only) were enriched for lower allele frequencies, lower F_ST, and putatively functional and regulatory roles (including LOF and deleterious missense variants). On average, individuals carried fewer private deleterious missense alleles than expected compared to alleles with other predicted consequences. Only a small subset of the low-prevalence variants had intermediate allele frequencies and explained small fractions of phenotypic variance (up to 3.2%) of production traits. The significant low-prevalence variants had higher per-site F_ST than the non-significant ones. These associated low-prevalence variants were tagged by other more widespread variants in high linkage disequilibrium, including intergenic variants.

CONCLUSIONS

Most low-prevalence variants have low minor allele frequencies and only a small subset of low-prevalence variants contributed detectable fractions of phenotypic variance of production traits. Accounting for low-prevalence variants is therefore unlikely to noticeably benefit across-breed analyses, such as the prediction of genomic breeding values in a population using reference populations of a different genetic background.

Whole-exome sequencing identifies rare genetic variants associated with human plasma metabolites

Metabolite levels measured in the human population are endophenotypes for biological processes. We combined sequencing data for 3,924 (whole-exome sequencing, WES, discovery) and 2,805 (whole-genome sequencing, WGS, replication) donors from a prospective cohort of blood donors in England. We used multiple approaches to select and aggregate rare genetic variants (minor allele frequency [MAF] < 0.1%) in protein-coding regions and tested their associations with 995 metabolites measured in plasma by using ultra-high-performance liquid chromatography-tandem mass spectrometry. We identified 40 novel associations implicating rare coding variants (27 genes and 38 metabolites), of which 28 (15 genes and 28 metabolites) were replicated. We developed algorithms to prioritize putative driver variants at each locus and used mediation and Mendelian randomization analyses to test directionality at associations of metabolite and protein levels at the ACY1 locus. Overall, 66% of reported associations implicate gene targets of approved drugs or bioactive drug-like compounds, contributing to drug targets' validating efforts.

A novel missense variant in the LMNB2 gene causes progressive myoclonus epilepsy

Progressive myoclonus epilepsies (PMEs) are a group of disorders embracing myoclonus, seizures, and neurological dysfunctions. Because of the genetic and clinical heterogeneity, a large proportion of PMEs cases have remained molecularly undiagnosed. The present study aimed to determine the underlying genetic factors that contribute to the PME phenotype in an Iranian female patient. We describe a consanguineous Iranian family with autosomal recessive PME that had remained undiagnosed despite extensive genetic and pathological tests. After performing neuroimaging and clinical examinations, due to heterogeneity of PMEs, the proband was subjected to paired-end whole-exome sequencing and the candidate variant was confirmed by Sanger sequencing. Various in-silico tools were also used to predict the pathogenicity of the variant. In this study, we identified a novel homozygous missense variant (NM_032737.4:c.472C > T; p.(Arg158Trp)) in the LMNB2 gene (OMIM: 150341) as the most likely disease-causing variant. Neuroimaging revealed a progressive significant generalized atrophy in the cerebral and cerebellum without significant white matter signal changes. Video-electroencephalography monitoring showed a generalized pattern of high-voltage sharp waves in addition to multifocal spikes and waves compatible with mixed type seizures and epileptic encephalopathic pattern. Herein, we introduce the second case of PME caused by a novel variant in the LMNB2 gene. This study also underscores the potentiality of next-generation sequencing in the genetic diagnosis of patients with neurologic diseases with an unknown cause.

Interpreting protein variant effects with computational predictors and deep mutational scanning

Computational predictors of genetic variant effect have advanced rapidly in recent years. These programs provide clinical and research laboratories with a rapid and scalable method to assess the likely impacts of novel variants. However, it can be difficult to know to what extent we can trust their results. To benchmark their performance, predictors are often tested against large datasets of known pathogenic and benign variants. These benchmarking data may overlap with the data used to train some supervised predictors, which leads to data re-use or circularity, resulting in inflated performance estimates for those predictors. Furthermore, new predictors are usually found by their authors to be superior to all previous predictors, which suggests some degree of computational bias in their benchmarking. Large-scale functional assays known as deep mutational scans provide one possible solution to this problem, providing independent datasets of variant effect measurements. In this Review, we discuss some of the key advances in predictor methodology, current benchmarking strategies and how data derived from deep mutational scans can be used to overcome the issue of data circularity. We also discuss the ability of such functional assays to directly predict clinical impacts of mutations and how this might affect the future need for variant effect predictors.

Mutation spectrum of ATP7B gene in pediatric patients with Wilson disease in Vietnam

Background

Wilson disease (WD) is caused by mutations in the copper-transporting P-type adenosine triphosphatase encoded by the ATP7B gene. In this study, we screened and identified the ATP7B mutations among unrelated Vietnamese pediatric patients.

Methods

One-hundred-thirteen pediatric patients with clinically diagnosed WD were recruited. DNA samples were extracted from peripheral blood. Mutations in the ATP7B gene were identified by Sanger sequencing.

Results

Approximately 98% of the clinically diagnosed WD patients carried ATP7B mutations. A total of 35 different ATP7B variants were detected, including five novel mutations (L658P, L792P, T977K, IVS4 + 1G > A and IVS20 + 4A > G). Remarkably, this study revealed that S105* was the most prevalent variant (32.27%), followed by L1371P (9.09%), I1148T (7.27%), R778L (6.36%), T850I (5.45%), V176Sfs*28 and IVS14-2A > G (4.55%). Most ATP7B mutations were located in the exon 2 (37.73%), exon 16 (10.00%), exon 8 (9.55%), exon 20 (9.09%), exon 10 and exon 18 (5.45%), exon 14 (5.00%), exon 13 and intron 14 (4.55%). We developed a streamlined procedure to quickly characterize mutations in the ATP7B gene in the Vietnamese children, starting with sequencing exon 2 and subsequently to exons 8,10,13-16,18, and 20 to allow quick diagnosis of clinically suspected patients.

Conclusion

The mutational spectrum and hotspots of ATP7B gene in the Vietnamese population were fairly different from other East Asian populations. A streamlined procedure was developed to screen exon 2 in ATP7B gene among suspected WD patients to reduce genetically diagnostic cost, to facilitate early detection and intervention in countries with limited resources.

Clinical and genetic findings in TRPM1-related congenital stationary night blindness

PURPOSE

Congenital stationary night blindness (CSNB) is a heterogeneous group of Mendelian retinal disorders that present in childhood. Biallelic variants altering the protein-coding region of the TRPM1 gene are one of the commonest causes of CSNB. Here, we report the clinical and genetic findings in 10 unrelated individuals with TRPM1-retinopathy.

METHODS

Study subjects were recruited through a tertiary clinical ophthalmic genetic service at Manchester, UK. All participants underwent visual electrodiagnostic testing and panel-based genetic analysis.

RESULTS

Study subjects had a median age of 8 years (range: 3-20 years). All probands were myopic and had electroretinographic findings in keeping with complete CSNB. Notably, three probands reported no night vision problems. Fourteen different disease-associated TRPM1 variants were detected. One individual was homozygous for the NM_001252024.2 (TRPM1):c.965 + 29G>A variant and a mini-gene assay highlighted that this change results in mis-splicing and premature protein termination. Additionally, two unrelated probands who had CSNB and mild neurodevelopmental abnormalities were found to carry a 15q13.3 microdeletion. This copy number variant encompasses seven genes, including TRPM1, and was encountered in the heterozygous state and in trans with a missense TRPM1 variant in each case.

CONCLUSION

Our findings highlight the importance of comprehensive genomic analysis, beyond the exons and protein-coding regions of genes, for individuals with CSNB. When this characteristic retinal phenotype is accompanied by extraocular findings (including learning and/or behavioural difficulties), a 15q13.3 microdeletion should be suspected. Focused analysis (e.g. microarray testing) is recommended to look for large-scale deletions encompassing TRPM1 in patients with CSNB and neurodevelopmental abnormalities.

Phenotype expansion of variants affecting p38 MAPK signaling in hypospadias patients

Background

Hypospadias is a congenital anomaly of the male urogenital system. Genetics factors play an important role in its pathogenesis. To search for potential causal genes/variants for hypospadias, we performed exome sequencing in a pedigree with three patients across two generations and a cohort of 49 sporadic patients with hypospadias.

Results

A novel BRAF variant (NM_004333.6: c.362C > A) was found to co-segregate with the hypospadias phenotype in the disease pedigree. In cells overexpressing the BRAF mutant, the phosphorylation level of p38 MAPK was significantly increased as compared with the cells overexpressing the wild-type BRAF or RASopathy-related BRAF mutant. This variant further led to a reduced transcription level of the SRY gene, which is essential for the normal development of the male reproductive system. In the cohort of sporadic patients, we identified two additional variants in p38 MAPK signaling-related genes (TRIM67 and DAB2IP) potentially associated with hypospadias.

Conclusion

Our study expands the phenotypic spectrum of variants affecting p38 MAPK signaling toward the involvement of hypospadias.

MiRLog and dbmiR: prioritization and functional annotation tools to study human microRNA sequence variants

The recent identification of noncoding variants with pathogenic effects suggests that these variations could underlie a significant number of undiagnosed cases. Several computational methods have been developed to predict the functional impact of noncoding variants, but they exhibit only partial concordance and are not integrated with functional annotation resources, making the interpretation of these variants still challenging. MicroRNAs (miRNAs) are small noncoding RNA molecules that act as fine regulators of gene expression and play crucial functions in several biological processes, such as cell proliferation and differentiation. An increasing number of studies demonstrate a significant impact of miRNA single nucleotide variants (SNVs) both in Mendelian diseases and complex traits. To predict the functional effect of miRNA SNVs, we implemented a new meta‐predictor, MiRLog, and we integrated it into a comprehensive database, dbmiR, which includes a precompiled list of all possible miRNA allelic SNVs, providing their biological annotations at nucleotide and miRNA levels. MiRLog and dbmiR were used to explore the genetic variability of miRNAs in 15,708 human genomes included in the gnomAD project, finding several ultra‐rare SNVs with a potentially deleterious effect on miRNA biogenesis and function representing putative contributors to human phenotypes.

Mutation analysis reveals novel and known mutations in SAG gene in first two Egyptian families with Oguchi disease

Background

Oguchi disease is a rare type of congenital stationary night blindness associated with an abnormal fundus appearance. It is inherited in an autosomal recessive manner where two types exist according to the gene affected; type 1 associated with S-antigen (SAG) gene mutations and type 2 associated with rhodopsin kinase (GRK1) gene mutations.

Purpose

The aim of this work was to describe the clinical and genetic findings of the first two reported families of Oguchi disease in Egypt and African region.

Methods

Four members of two consanguineous Egyptian families with history of night blindness since childhood underwent complete ophthalmological examination, standard automated static perimetry, fundus color photography, fundus autofluorescence (FAF), fundus fluorescein angiography (FFA) in light-adapted state and spectral-domain optical coherence tomography (SD-OCT) of both the macula and the optic nerve head as well as central corneal thickness with repeated fundus photography following prolonged dark adaptation. Mutation screening of 7 coding exons of GRK1 gene and 15 coding exons of SAG gene as well as some flanking regions were performed using Sanger sequencing technique. The variants were tested for pathogenicity using different in silico functional analysis tools.

Results

The clinical examination and investigations confirmed Oguchi disease phenotype. One patient showed p.R193* (c.577C > T) which is a previously reported SAG gene mutation in a homozygous form. The other three patients from a different family showed (c.649–1 G > C), a novel canonical splice site SAG gene mutation in a homozygous form.

Conclusion

The identification of the novel canonical splice site SAG gene variant in three members of the same family with clinically confirmed Oguchi disease reinforces its pathogenicity. A fourth patient from another family carried a previously reported mutation in the same gene. SAG gene variants may be the underlying genetic cause for Oguchi disease in Egypt. Our findings have expanded the spectrum of Oguchi disease-associated mutations in SAG gene and may serve as a basis for genetic diagnosis for Oguchi disease.

Integrating convolution and self-attention improves language model of human genome for interpreting non-coding regions at base-resolution

Interpretation of non-coding genome remains an unsolved challenge in human genetics due to impracticality of exhaustively annotating biochemically active elements in all conditions. Deep learning based computational approaches emerge recently to help interpret non-coding regions. Here, we present LOGO (Language of Genome), a self-attention based contextualized pre-trained language model containing only two self-attention layers with 1 million parameters as a substantially light architecture that applies self-supervision techniques to learn bidirectional representations of the unlabelled human reference genome. LOGO is then fine-tuned for sequence labelling task, and further extended to variant prioritization task via a special input encoding scheme of alternative alleles followed by adding a convolutional module. Experiments show that LOGO achieves 15% absolute improvement for promoter identification and up to 4.5% absolute improvement for enhancer-promoter interaction prediction. LOGO exhibits state-of-the-art multi-task predictive power on thousands of chromatin features with only 3% parameterization benchmarking against the fully supervised model, DeepSEA and 1% parameterization against a recent BERT-based DNA language model. For allelic-effect prediction, locality introduced by one dimensional convolution shows improved sensitivity and specificity for prioritizing non-coding variants associated with human diseases. In addition, we apply LOGO to interpret type 2 diabetes (T2D) GWAS signals and infer underlying regulatory mechanisms. We make a conceptual analogy between natural language and human genome and demonstrate LOGO is an accurate, fast, scalable, and robust framework to interpret non-coding regions for global sequence labeling as well as for variant prioritization at base-resolution.

Evaluation of Evidence for Pathogenicity Demonstrates That BLK, KLF11, and PAX4 Should Not Be Included in Diagnostic Testing for MODY

Maturity-onset diabetes of the young (MODY) is an autosomal dominant form of monogenic diabetes, reported to be caused by variants in 16 genes. Concern has been raised about whether variants in BLK (MODY11), KLF11 (MODY7), and PAX4 (MODY9) cause MODY. We examined variant-level genetic evidence (cosegregation with diabetes and frequency in population) for published putative pathogenic variants in these genes and used burden testing to test gene-level evidence in a MODY cohort (n = 1,227) compared with a control population (UK Biobank [n = 185,898]). For comparison we analyzed well-established causes of MODY, HNF1A, and HNF4A. The published variants in BLK, KLF11, and PAX4 showed poor cosegregation with diabetes (combined logarithm of the odds [LOD] scores ≤1.2), compared with HNF1A and HNF4A (LOD scores >9), and are all too common to cause MODY (minor allele frequency >4.95 × 10-5). Ultra-rare missense and protein-truncating variants (PTV) were not enriched in a MODY cohort compared with the UK Biobank population (PTV P > 0.05, missense P > 0.1 for all three genes) while HNF1A and HNF4A were enriched (P < 10-6). Findings of sensitivity analyses with different population cohorts supported our results. Variant and gene-level genetic evidence does not support BLK, KLF11, or PAX4 as a cause of MODY. They should not be included in MODY diagnostic genetic testing.

Clinical and genetic characterization of three Russian patients with pycnodysostosis due to pathogenic variants in the CTSK gene

BACKGROUND

Pycnodysostosis (PD, OMIM # 265800) is a rare variant of skeletal dysplasia with an autosomal recessive type of inheritance, characterized by a combination of specific features such as disproportionate nanism, generalized osteosclerosis, and distinct craniofacial dysmorphism. Radiographic features include acro-osteolysis of the distal phalanges in association with sclerosing bone lesions with multiple fractures. The polymorphism of the clinical manifestations of pycnodysostosis and low prevalence of the disorder lead to the difficulties with early.

METHODS

The following tests were used for diagnostics: genealogical analysis, clinical examination, neurological examination according to the standard method with an assessment of the psychoemotional sphere, radiological analysis, searching for pathogenic variants in the CTSK gene by the automated Sanger sequencing.

RESULTS

We describe first clinical and genetic characteristics of three Russian patients with pycnodysostosis from unrelated families. Two patients have a novel homozygous nucleotide substitution c.746T>A (p. Ile249Asn), and one has a previously described homozygous pathogenic variant c.746T>C (p.Ile249Thr) in the CTSK gene. In all three cases, a transition or transversion was found at nucleotide position 746 in exon 6 of the CTSK gene, leading to two different amino acid substitutions in the polypeptide chain. The obtained results may indicate the presence of a major pathogenic variant in the CTSK gene, leading to the typical manifestation of the disease.

CONCLUSION

The data presented in the study enlarge the clinical, radiological, and mutational spectrum of pycnodysostosis. Typical clinical manifestations and the small size of the CTSK gene make the automated Sanger sequencing the optimal method for diagnosis of pycnodysostosis.

Computational Saturation Mutagenesis to Investigate the Effects of Neurexin-1 Mutations on AlphaFold Structure

Neurexin-1 (NRXN1) is a membrane protein essential in synapse formation and cell signaling as a cell-adhesion molecule and cell-surface receptor. NRXN1 and its binding partner neuroligin have been associated with deficits in cognition. Recent genetics research has linked NRXN1 missense mutations to increased risk for brain disorders, including schizophrenia (SCZ) and autism spectrum disorder (ASD). Investigation of the structure-function relationship in NRXN1 has proven difficult due to a lack of the experimental full-length membrane protein structure. AlphaFold, a deep learning-based predictor, succeeds in high-quality protein structure prediction and offers a solution for membrane protein model construction. In the study, we applied a computational saturation mutagenesis method to analyze the systemic effects of missense mutations on protein functions in a human NRXN1 structure predicted from AlphaFold and an experimental Bos taurus structure. The folding energy changes were calculated to estimate the effects of the 29,540 mutations of AlphaFold model on protein stability. The comparative study on the experimental and computationally predicted structures shows that these energy changes are highly correlated, demonstrating the reliability of the AlphaFold structure for the downstream bioinformatics analysis. The energy calculation revealed that some target mutations associated with SCZ and ASD could make the protein unstable. The study can provide helpful information for characterizing the disease-causing mutations and elucidating the molecular mechanisms by which the variations cause SCZ and ASD. This methodology could provide the bioinformatics protocol to investigate the effects of target mutations on multiple AlphaFold structures.

The Contribution of Rare Copy Number Variants in FAS Towards Pathogenesis of Autoimmune Lymphoproliferative Syndrome

Rare copy number variants in FAS can cause ALPS by method of haploinsufficiency.

CNV analysis can provide a molecular diagnosis for patients with ALPS when whole exome sequencing or panel-based testing are inconclusive.

Autoimmune lymphoproliferative syndrome (ALPS) is characterized by chronic nonmalignant lymphadenopathy, splenomegaly, cytopenias, and other autoimmune manifestations. ALPS is caused by lymphocyte accumulation from defects in FAS-mediated apoptosis. Heterozygous germline or somatic pathogenic single nucleotide variants in FAS are the most common molecular etiology of ALPS. Through the Centralized Sequencing Program at the National Institute of Allergy and Infectious Diseases, we performed exome sequencing on subjects with a clinical diagnosis of ALPS, with a subset receiving copy number variant (CNV) analysis. In this cohort, we identified 3 subjects from unrelated families with CNVs at the FAS locus. One subject had a de novo ∼0.828 Mb copy number loss encompassing all of FAS. The second subject had a maternally inherited ∼1.004 Mb copy number loss encompassing all of FAS. The third subject had a paternally inherited ∼0.044 Mb copy number loss encompassing exons 7 through 9 of FAS. Subjects with deletions in FAS had clinical presentations and biomarker profiles similar to those with ALPS and with germline and somatic FAS variants. We demonstrate that CNV analysis should be pursued if there is clinical and biomarker evidence of ALPS because it can lead to a molecular diagnosis and appropriate treatment when FAS sequencing is inconclusive.

Ferroptosis and Apoptosis Are Involved in the Formation of L-Selenomethionine-Induced Ocular Defects in Zebrafish Embryos

Selenium is an essential trace element for humans and other vertebrates, playing an important role in antioxidant defense, neurobiology and reproduction. However, the toxicity of excessive selenium has not been thoroughly evaluated, especially for the visual system of vertebrates. In this study, fertilized zebrafish embryos were treated with 0.5 µM L-selenomethionine to investigate how excessive selenium alters zebrafish eye development. Selenium-stressed zebrafish embryos showed microphthalmia and altered expression of genes required for retinal neurogenesis. Moreover, ectopic proliferation, disrupted mitochondrial morphology, elevated ROS-induced oxidative stress, apoptosis and ferroptosis were observed in selenium-stressed embryos. Two antioxidants—reduced glutathione (GSH) and N-acetylcysteine (NAC)—and the ferroptosis inhibitor ferrostatin (Fer-1) were unable to rescue selenium-induced eye defects, but the ferroptosis and apoptosis activator cisplatin (CDDP) was able to improve microphthalmia and the expression of retina-specific genes in selenium-stressed embryos. In summary, our results reveal that ferroptosis and apoptosis might play a key role in selenium-induced defects of embryonic eye development. The findings not only provide new insights into selenium-induced cellular damage and death, but also important implications for studying the association between excessive selenium and ocular diseases in the future.

Structural pharmacology of PTH and PTHrP

Parathyroid hormone (PTH) and PTH-related peptide (PTHrP) regulate extracellular phosphate and calcium homeostasis as well as bone remodeling. PTH is a classic endocrine peptide hormone whose synthesis and negative feedback by multiple factors control release from the parathyroid glands. PTHrP is ubiquitously expressed (pre- and postnatally) and acts in an autocrine/paracrine manner. This review considers the structural pharmacology and actions of PTH and PTHrP, biological consequences of inherited mutations, engineered analogs that illuminate similarities and differences in physiologic actions, and targeted therapeutic opportunities.

In Silico Analysis of the L-2-Hydroxyglutarate Dehydrogenase Gene Mutations and Their Biological Impact on Disease Etiology

The L-2-hydroxyglutarate dehydrogenase (L2HGDH) gene encodes an important mitochondrial enzyme. However, its altered activity results in excessive levels of L-2-hydroxyglutarate, which results in diverse psychiatric features of intellectual disability. In the current study, we executed an in-silico analysis of all reported L2HGDH missense and nonsense variants in order to investigate their biological significance. Among the superimposed 3D models, the highest similarity index for a wild-type structure was shown by the mutant Glu336Lys (87.26%), while the lowest similarity index value was shown by Arg70* (10.00%). Three large active site pockets were determined using protein active site prediction, in which the 2nd largest pocket was shown to encompass the substrate L-2-hydroxyglutarate (L2HG) binding residues, i.e., 89Gln, 195Tyr, 402Ala, 403Gly and 404Val. Moreover, interactions of wild-type and mutant L2HGDH variants with the close functional interactor D2HGDH protein resulted in alterations in the position, number and nature of networking residues. We observed that the binding of L2HG with the L2HGDH enzyme is affected by the nature of the amino acid substitution, as well as the number and nature of bonds between the substrate and protein molecule, which are able to affect its biological activity.

Clinical Sphingolipids Pathway in Parkinson’s Disease: From GCase to Integrated-Biomarker Discovery

Alterations in the sphingolipid metabolism of Parkinson’s Disease (PD) could be a potential diagnostic feature. Only around 10–15% of PD cases can be diagnosed through genetic alterations, while the remaining population, idiopathic PD (iPD), manifest without validated and specific biomarkers either before or after motor symptoms appear. Therefore, clinical diagnosis is reliant on the skills of the clinician, which can lead to misdiagnosis. IPD cases present with a spectrum of non-specific symptoms (e.g., constipation and loss of the sense of smell) that can occur up to 20 years before motor function loss (prodromal stage) and formal clinical diagnosis. Prodromal alterations in metabolites and proteins from the pathways underlying these symptoms could act as biomarkers if they could be differentiated from the broad values seen in a healthy age-matched control population. Additionally, these shifts in metabolites could be integrated with other emerging biomarkers/diagnostic tests to give a PD-specific signature. Here we provide an up-to-date review of the diagnostic value of the alterations in sphingolipids pathway in PD by focusing on the changes in definitive PD (postmortem confirmed brain data) and their representation in “probable PD” cerebrospinal fluid (CSF) and blood. We conclude that the trend of holistic changes in the sphingolipid pathway in the PD brain seems partly consistent in CSF and blood, and could be one of the most promising pathways in differentiating PD cases from healthy controls, with the potential to improve early-stage iPD diagnosis and distinguish iPD from other Parkinsonism when combined with other pathological markers.

Predominant Founder Effect among Recurrent Pathogenic Variants for an X-Linked Disorder

For disorders with X-linked inheritance, variants may be transmitted through multiple generations of carrier females before an affected male is ascertained. Pathogenic RS1 variants exclusively cause X-linked retinoschisis (XLRS). While RS1 is constrained to variation, recurrent variants are frequently observed in unrelated probands. Here, we investigate recurrent pathogenic variants to determine the relative burden of mutational hotspot and founder allele events to this phenomenon. A cohort RS1 variant analysis and standardized classification, including variant enrichment in the XLRS cohort and in RS1 functional domains, were performed on 332 unrelated XLRS probands. A total of 108 unique RS1 variants were identified. A subset of 19 recurrently observed RS1 variants were evaluated in 190 probands by a haplotype analysis, using microsatellite and single nucleotide polymorphisms. Fourteen variants had at least two probands with common variant-specific haplotypes over ~1.95 centimorgans (cM) flanking RS1. Overall, 99/190 of reportedly unrelated probands had 25 distinct shared haplotypes. Examination of this XLRS cohort for common RS1 haplotypes indicates that the founder effect plays a significant role in this disorder, including variants in mutational hotspots. This improves the accuracy of clinical variant classification and may be generalizable to other X-linked disorders.

Identification, function, and biological relevance of POGLUT2 and POGLUT3

O-glycosylation of Epidermal Growth Factor-like (EGF) repeats plays crucial roles in protein folding, trafficking and function. The Notch extracellular domain has been used as a model to study these mechanisms due to its many O-glycosylated EGF repeats. Three enzymes were previously known to O-glycosylate Notch EGF repeats: Protein O-Glucosyltransferase 1 (POGLUT1), Protein O-Fucosyltransferase 1 (POFUT1), and EGF Domain Specific O-Linked N-Acetylglucosamine Transferase (EOGT). All of these modifications affect Notch activity. Recently, POGLUT2 and POGLUT3 were identified as two novel O-glucosyltransferases that modify a few Notch EGF repeats at sites distinct from those modified by POGLUT1. Comparison of these modification sites revealed a putative consensus sequence which predicted modification of many extracellular matrix proteins including fibrillins (FBNs) and Latent TGFβ-binding proteins (LTBPs). Glycoproteomic analysis revealed that approximately half of the 47 EGF repeats in FBN1 and FBN2, and half of the 18 EGF repeats in LTBP1, are modified by POGLUT2 and/or POGLUT3. Cellular assays showed that loss of modifications by POGLUT2 and/or POGLUT3 significantly reduces FBN1 secretion. There is precedent for EGF modifications to affect protein-protein interactions, as has been demonstrated by research of POGLUT1 and POFUT1 modifications on Notch. Here we discuss the identification and characterization of POGLUT2 and POGLUT3 and the ongoing research that continues to elucidate the biological significance of these novel enzymes.

Structural basis of activation of the tumor suppressor protein neurofibromin

Mutations in the NF1 gene cause the familial genetic disease neurofibromatosis type I, as well as predisposition to cancer. The NF1 gene product, neurofibromin, is a GTPase-activating protein and acts as a tumor suppressor by negatively regulating the small GTPase, Ras. However, structural insights into neurofibromin activation remain incompletely defined. Here, we provide cryoelectron microscopy (cryo-EM) structures that reveal an extended neurofibromin homodimer in two functional states: an auto-inhibited state with occluded Ras-binding site and an asymmetric open state with an exposed Ras-binding site. Mechanistically, the transition to the active conformation is stimulated by nucleotide binding, which releases a lock that tethers the catalytic domain to an extended helical repeat scaffold in the occluded state. Structure-guided mutational analysis supports functional relevance of allosteric control. Disease-causing mutations are mapped and primarily impact neurofibromin stability. Our findings suggest a role for nucleotides in neurofibromin regulation and may lead to therapeutic modulation of Ras signaling.

Neurodevelopmental disorder with microcephaly, hypotonia, and variable brain anomalies in a consanguineous Iranian family is associated with a homozygous start loss variant in the PRUNE1 gene

Background

Homozygous or compound heterozygous PRUNE1 mutations cause a neurodevelopmental disorder with microcephaly, hypotonia, and variable brain malformations (NMIHBA) (OMIM #617481). The PRUNE1 gene encodes a member of the phosphoesterase (DHH) protein superfamily that is involved in the regulation of cell migration. To date, most of the described mutations in the PRUNE1 gene are clustered in DHH domain.

Methods

We subjected 4 members (two affected and two healthy) of a consanguineous Iranian family in the study. The proband underwent whole-exome sequencing and a start loss identified variant was confirmed by Sanger sequencing. Co-segregation of the detected variant with the disease in family was confirmed.

Results

By whole-exome sequencing, we identified the a start loss variant, NM_021222.3:c.3G>A; p.(Met1?), in the PRUNE1 in two patients of a consanguineous Iranian family with spastic quadriplegic cerebral palsy (CP), hypotonia, developmental regression, and cerebellar atrophy. Sanger sequencing confirmed the segregation of the variant with the disease in the family. Protein structure analysis also revealed that the variant probably leads to the deletion of DHH (Asp-His-His) domain, the active site of the protein, and loss of PRUNE1 function.

Conclusion

We identified a start loss variant, NM_021222.3:c.3G>A; p.(Met1?) in the PRUNE1 gene in two affected members as a possible cause of NMIHBA in an Iranian family. We believe that the study adds a new pathogenic variant in spectrum of mutations in the PRUNE1 gene as a cause of PRUNE1-related syndrome.

Dynamic coupling of residues within proteins as a mechanistic foundation of many enigmatic pathogenic missense variants

Many pathogenic missense mutations are found in protein positions that are neither well-conserved nor fall in any known functional domains. Consequently, we lack any mechanistic underpinning of dysfunction caused by such mutations. We explored the disruption of allosteric dynamic coupling between these positions and the known functional sites as a possible mechanism for pathogenesis. In this study, we present an analysis of 591 pathogenic missense variants in 144 human enzymes that suggests that allosteric dynamic coupling of mutated positions with known active sites is a plausible biophysical mechanism and evidence of their functional importance. We illustrate this mechanism in a case study of β-Glucocerebrosidase (GCase) in which a vast majority of 94 sites harboring Gaucher disease-associated missense variants are located some distance away from the active site. An analysis of the conformational dynamics of GCase suggests that mutations on these distal sites cause changes in the flexibility of active site residues despite their distance, indicating a dynamic communication network throughout the protein. The disruption of the long-distance dynamic coupling caused by missense mutations may provide a plausible general mechanistic explanation for biological dysfunction and disease.

Variants in CDH23 cause a broad spectrum of hearing loss: from non-syndromic to syndromic hearing loss as well as from congenital to age-related hearing loss

Variants in the CDH23 gene are known to be responsible for both syndromic hearing loss (Usher syndrome type ID: USH1D) and non-syndromic hearing loss (DFNB12). Our series of studies demonstrated that CDH23 variants cause a broad range of phenotypes of non-syndromic hearing loss (DFNB12); from congenital profound hearing loss to late-onset high-frequency-involved progressive hearing loss. In this study, based on the genetic and clinical data from more than 10,000 patients, the mutational spectrum, clinical characteristics and genotype/phenotype correlations were evaluated. The present results reconfirmed that the variants in CDH23 are an important cause of non-syndromic sensorineural hearing loss. In addition, we showed that the mutational spectrum in the Japanese population, which is probably representative of the East Asian population in general, as well as frequent CDH23 variants that might be due to some founder effects. The present study demonstrated CDH23 variants cause a broad range of phenotypes, from non-syndromic to syndromic hearing loss as well as from congenital to age-related hearing loss. Genotype (variant combinations) and phenotype (association with retinal pigmentosa, onset age) are shown to be well correlated and are thought to be related to the residual function defined by the CDH23 variants.

Postnatal Progressive Craniosynostosis in Syndromic Conditions: Two Patients With Saethre-Chotzen Due to TWIST1 Gene Deletions and Review of the Literature

Saethre-Chotzen syndrome (SCS) is a known craniosynostosis syndrome with a variable presentation of craniofacial and somatic involvement. Congenital coronal craniosynostosis is most commonly observed in SCS; however, progressive postnatal craniosynostosis of other sutures has been reported. The authors present 2 infants with progressive postnatal craniosynostosis and SCS caused by chromosome 7p deletions including the TWIST1 gene. The evolution of their clinical features and a literature review of patients with syndromic, postnatal progressive craniosynostosis illustrate the importance of longitudinal observation and management of these patients.

Human iPSC-Derived Retinal Organoids and Retinal Pigment Epithelium for Novel Intronic RPGR Variant Assessment for Therapy Suitability

The RPGR gene encodes Retinitis Pigmentosa GTPase Regulator, a known interactor with ciliary proteins, which is involved in maintaining healthy photoreceptor cells. Variants in RPGR are the main contributor to X-linked rod-cone dystrophy (RCD), and RPGR gene therapy approaches are in clinical trials. Hence, elucidation of the pathogenicity of novel RPGR variants is important for a patient therapy opportunity. Here, we describe a novel intronic RPGR variant, c.1415 − 9A>G, in a patient with RCD, which was classified as a variant of uncertain significance according to current clinical diagnostic criteria. The variant lay several base pairs intronic to the canonical splice acceptor site, raising suspicion of an RPGR RNA splicing abnormality and consequent protein dysfunction. To investigate disease causation in an appropriate disease model, induced pluripotent stem cells were generated from patient fibroblasts and differentiated to retinal pigment epithelium (iPSC-RPE) and retinal organoids (iPSC-RO). Abnormal RNA splicing of RPGR was demonstrated in patient fibroblasts, iPSC-RPE and iPSC-ROs, leading to a predicted frameshift and premature stop codon. Decreased RPGR expression was demonstrated in these cell types, with a striking loss of RPGR localization at the ciliary transitional zone, critically in the photoreceptor cilium of the patient iPSC-ROs. Mislocalisation of rhodopsin staining was present in the patient’s iPSC-RO rod photoreceptor cells, along with an abnormality of L/M opsin staining affecting cone photoreceptor cells and increased photoreceptor apoptosis. Additionally, patient iPSC-ROs displayed an increase in F-actin expression that was consistent with an abnormal actin regulation phenotype. Collectively, these studies indicate that the splicing abnormality caused by the c.1415 − 9A>G variant has an impact on RPGR function. This work has enabled the reclassification of this variant to pathogenic, allowing the consideration of patients with this variant having access to gene therapy clinical trials. In addition, we have identified biomarkers of disease suitable for the interrogation of other RPGR variants of uncertain significance.

Screening of OTULIN gene mutation with targeted next generation sequencing in Turkish populations and in silico analysis of these mutations

BACKGROUND

OTULIN-related autoinflammatory syndrome (ORAS) is an autosomal recessive disease characterized by systemic inflammation, recurrent fever. Due to limited knowledge about the OTULIN DNA variants that cause ORAS, the diagnosis and treatment of this disease is difficult. In this study, we aim to identify OTULIN DNA variants responsible for the genetic pathology of ORAS and observe the effects of these variants on the OTULIN protein structure and the function with different bioinformatics approaches.

METHODS

The present study included 3230 individuals with the suspicion of an autoinflammatory disease who were referred to Ege University Children's Hospital Molecular Medicine Laboratory. OTULIN variants were detected using a panel consisting of 37 different autoinflammatory diseases (AID) genes via targeted Next-Generation Sequencing.

RESULTS

As a result of the study, DNA variants associated with various AID were detected in 65% of the individuals to whom the panel was applied. Among these variants, only three different OTULIN variants (p.Val82Ile, p.Gln115His and p.Leu131_Arg132insLeuCysThrGlu) were detected. The pathogenic effects of the variants detected in the OTULIN gene were determined by using Polyphen2 as "Probably Pathogenic" for the p.Val82Ile and "benign" for the p.Gln115His. At the same time, the effects of these variants on the structure and function of the OTULIN protein were investigated by in silico approaches. Both variants reduce protein stability and binding affinity.

CONCLUSION

The results of the current study suggest that the evaluation of OTULIN variants with in silico approaches will contribute to the development of personalized treatments by diagnosing the disease specific to the variant.

LIM domain-wide comprehensive virtual mutagenesis provides structural rationale for cardiomyopathy mutations in CSRP3

Cardiomyopathies are a severe and chronic cardiovascular burden worldwide, affecting a large cohort in the general population. Cysteine and glycine-rich protein 3 (CSRP3) is one of key proteins implicated in dominant dilated cardiomyopathy (DCM) and hypertrophic cardiomyopathy (HCM). In this study, we device a rapid in silico screening protocol that creates a mutational landscape map for all possible allowed and disallowed substitutions in the protein of interest. This map provides the structural and functional insights on the stability of LIM domains of CSRP3. Further, the sequence analysis delineates the eukaryotic CSRP3 protein orthologs which complements the mutational map, but provide limited information of amino acid exchanges. Next, we also evaluated the effect of HCM/DCM mutations on these domains. One of highly destabilising mutations—L44P (also disease causing) and a neutral mutation—L44M were further subjected to molecular dynamics (MD) simulations. The results establish that L44P substitution affects the LIM domain structure by altering secondary structure and due to loss of hydrophobic interaction with Phenylananine 35. The present study provides a useful perspective to our understanding of the role of mutations in the CSRP3 LIM domains and their evolution. This study provides a novel computational screening method for quick identification of key mutation sites for specific protein structures that can reduce the burden on experimental research.

De novo PBX1 variant in a patient with glaucoma, kidney anomalies, and developmental delay: An expansion of the CAKUTHED phenotype

An infant was referred for evaluation of congenital glaucoma and corneal clouding. In addition, he had a pelvic kidney, hypotonia, patent ductus arteriosus, abnormal pinnae, and developmental delay. Exome sequencing identified a previously unpublished de novo single nucleotide insertion in PBX1 c.400dupG (NM_002585.3), predicted to cause a frameshift resulting in a truncated protein with loss of function (p.Ala134Glyfs*65). Identification of this loss of function variant supports the diagnosis of congenital anomalies of the kidney and urinary tract syndrome with or without hearing loss, abnormal ears, or developmental delay (CAKUTHED). Here, we propose glaucoma as an extra-renal manifestation associated with PBX1-related disease due to the relationship of PBX1 with MEIS1, MEIS2, and FOXC1 transcription factors associated with eye development.

Molecular Diagnostic Outcomes from 700 Cases: What Can We Learn from a Retrospective Analysis of Clinical Exome Sequencing?

Clinical exome sequencing (CES) aids in the diagnosis of rare genetic disorders. Herein, we report the molecular diagnostic yield and spectrum of genetic alterations contributing to disease in 700 pediatric cases analyzed at the Children's Hospital of Philadelphia. The overall diagnostic yield was 23%, with three cases having more than one molecular diagnosis and 2.6% having secondary/additional findings. A candidate gene finding was reported in another 8.4% of cases. The clinical indications with the highest diagnostic yield were neurodevelopmental disorders (including seizures), whereas immune- and oncology-related indications were negatively associated with molecular diagnosis. The rapid expansion of knowledge regarding the genome's role in human disease necessitates reanalysis of CES samples. To capture these new discoveries, a subset of cases (n = 240) underwent reanalysis, with an increase in diagnostic yield. We describe our experience reporting CES results in a pediatric setting, including reporting of secondary findings, reporting newly discovered genetic conditions, and revisiting negative test results. Finally, we highlight the challenges associated with implementing critical updates to the CES workflow. Although these updates are necessary, they demand an investment of time and resources from the laboratory. In summary, these data demonstrate the clinical utility of exome sequencing and reanalysis, while highlighting the critical considerations for continuous improvement of a CES test in a clinical laboratory.

Association study of CLDN14 variations in patients with kidney stones

Claudin-14 protein plays an essential role in regulating calcium ions in the kidney and ear. Two phenotypes, hearing loss and kidney stones, were reportedly associated with variations in the CLDN14 gene. This study aimed to understand CLDN14 mutations’ contribution to hearing loss and renal stone formation in a Pakistani cohort. We analyzed CLDN14 sequence variations in 100 patients, along with healthy individuals, to assess whether specific polymorphisms were associated with the disease. Also, we performed an in silico analysis using a mutation database and protein annotation. The rs219779’s genotype CT (p = 0.0020) and rs219780’s genotype AG (p = 0.0012) were significantly associated with kidney stones. We also found that a novel haplotype, “TA” associated with kidney stone formation, has moderate linkage disequilibrium. The TA haplotype was significantly correlated with a kidney stone risk formation of 3.76-fold (OR (CI 95%) = 3.76 (1.83–7.72)) and p = 0.0016 compared to other haplotypes. In silico analysis revealed that mutations associated with hearing loss were not correlated with renal stone formation but affected claudin-14 protein stability. We structurally mapped a novel TA haplotype of CLDN14 that, based on our analysis, likely contributes to the pathogenesis of renal stones.