The Human Gene Mutation Database (HGMD®): optimizing its use in a clinical diagnostic or research setting

Genetic variant classification by predicted protein structure: A case study on IRF6

Next-generation genome sequencing has revolutionized genetic testing, identifying numerous rare disease-associated gene variants. However, to impute pathogenicity, computational approaches remain inadequate and functional testing of gene variant is required to provide the highest level of evidence. The emergence of AlphaFold2 has transformed the field of protein structure determination, and here we outline a strategy that leverages predicted protein structure to enhance genetic variant classification. We used the gene IRF6 as a case study due to its clinical relevance, its critical role in cleft lip/palate malformation, and the availability of experimental data on the pathogenicity of IRF6 gene variants through phenotype rescue experiments in irf6-/- zebrafish. We compared results from over 30 pathogenicity prediction tools on 37 IRF6 missense variants. IRF6 lacks an experimentally derived structure, so we used predicted structures to explore associations between mutational clustering and pathogenicity. We found that among these variants, 19 of 37 were unanimously predicted as deleterious by computational tools. Comparing in silico predictions with experimental findings, 12 variants predicted as pathogenic were experimentally determined as benign. Even with the recently published AlphaMissense model, 15/18 (83%) of the predicted pathogenic variants were experimentally determined as benign. In comparison, mapping variants to the protein revealed deleterious mutation clusters around the protein binding domain, whereas N-terminal variants tend to be benign, suggesting the importance of structural information in determining pathogenicity of mutations in this gene. In conclusion, incorporating gene-specific structural features of known pathogenic/benign mutations may provide meaningful insights into pathogenicity predictions in a gene-specific manner and facilitate the interpretation of variant pathogenicity.

Databases of ligand-binding pockets and protein-ligand interactions

Many research groups and institutions have created a variety of databases curating experimental and predicted data related to protein-ligand binding. The landscape of available databases is dynamic, with new databases emerging and established databases becoming defunct. Here, we review the current state of databases that contain binding pockets and protein-ligand binding interactions. We have compiled a list of such databases, fifty-three of which are currently available for use. We discuss variation in how binding pockets are defined and summarize pocket-finding methods. We organize the fifty-three databases into subgroups based on goals and contents, and describe standard use cases. We also illustrate that pockets within the same protein are characterized differently across different databases. Finally, we assess critical issues of sustainability, accessibility and redundancy.

Reclassification of an FBN1 variant emphasizes the importance of segregation analysis, information sharing, and multidisciplinary teamwork in understanding genetic variants in health and disease

Marfan syndrome (MFS) is a complex connective tissue disorder characterized by considerable clinical variability. The diagnosis of MFS is based on the Ghent criteria, which require the presence of both clinical and genetic features. MFS is primarily caused by pathogenic alterations in FBN1, which encodes the fibrillin-1 protein. Fibrillin-1 comprises multiple domains rich in cysteine residues, with disulfide bonds formed between these residues. It has long been recognized that variants that alter or introduce cysteine residues damage protein function, leading to the development of MFS. In this study, we report a cysteine-introducing variant: FBN1 variant, c.6724C>T (p.[Arg2242Cys]). We have observed this variant in several individuals without MFS, challenging our previous understanding of the underlying mechanism of MFS. This finding emphasizes the importance of revisiting and reevaluating our current knowledge in light of new and unexpected observations. Moreover, our study highlights the significance of incorporating local and national data on allele frequencies, as well as employing multidisciplinary phenotyping approaches, in the classification of genetic variants. By considering a wide range of information, we can enhance the accuracy and reliability of variant classification, ultimately improving the diagnosis and management of individuals with genetic disorders like MFS.

Comprehensive genetic analysis uncovers the mutational spectrum of MFRP and its genotype-phenotype correlation in a large cohort of Chinese microphthalmia patients

PURPOSE

Microphthalmia is a severe congenital ocular disease featured by abnormal ocular development. The aim of this study was to detail the genetic and clinical characteristics of a large cohort of Chinese patients with microphthalmia related to MFRP variants, focusing on uncovering genotype-phenotype correlations.

METHODS

Fifty microphthalmia patients from 44 unrelated Chinese families were recruited. Whole-exome sequencing (WES) was conducted to analyze the coding regions and adjacent intronic regions of MFRP. Axial lengths (AL) were measured for all probands and available family members. Protein structures of mutations with high frequency in our cohort were predicted. The genotype-phenotype correlations were explored by statistical analysis.

RESULTS

Sixteen MFRP variants were detected in 17 families, accounting for 38.64 % of all microphthalmia families. There were 9 novel mutations (c.427+1G>C, c.428-2A>C, c.561_575del:p.A188_E192del, c.836G>A:p.C279Y, c.1010_1021del:p.H337_E340del:p.Y479*, c.1516_1517del:p.S506Pfs*66, c.1561T>G:p.C521G, c.1616G>A:p.R539H, and c.1735C>T:p.P579S) and six previously reported variants in MFRP, with p.E496K and p.H337_E340del being highly frequent, found in eight (47.06 %) and two families (11.76 %), respectively. Seven variants (43.75 %) were located in the C-terminal cysteine-rich frizzled-related domain (CRD) (7/16, 43.75 %). Protein prediction implicated p.E496K and p.H337_E340del mutations might lead to a destabilization of the MFRP protein. The average AL of all 42 eyes was 16.02 ± 1.05 mm, and 78.36 % of eyes with AL < 16 mm harbored p.E496K variant. Twenty-six eyes with variant variant had shorter AL than that of the other 16 eyes without this variant (p = 0.006), highlighting a novel genotype-phenotype correlation.

CONCLUSIONS

In this largest cohort of Chinese patients with microphthalmia, the 9 novel variants, high frequency of p.E496W, and mutation hotspots in CRD reveals unique insights into the MFRP mutation spectrum among Chinese patients, indicating ethnic variability. A new genotype-phenotype correlation that p.E496K variant associated with a shorter AL is unveiled. Our findings enhance the current knowledge of MFRP-associated microphthalmia and provide valuable information for prenatal diagnosis as well as future therapy.

Molecular dynamics simulations involving different β-propeller mutations reported in Swiss and French patients correlate with their disease phenotypes

Integrin αIIbβ3 is the predominant receptor for fibrinogen which mediates platelet aggregation, an important step in hemostasis and thrombosis. Several mutations have been reported in the genes encoding αIIb and β3 subunits among patients with Glanzmann thrombasthenia, of which 177 are in the β-propeller domain. The two subunits form a heterodimer at the interface between β-propeller and β-I domains of αIIb and β3, respectively with their stability critical for intracellular trafficking, surface expression, and ligand binding. Our study was aimed at retrieving the β-propeller mutations from various databases and studying structural variations due to select mutations upon interaction with fibrinogen using molecular docking and molecular dynamics. Mutations were studied for their impact on phenotypic severity, structural stability, and evolutionary conservation. Molecular docking analysis and molecular dynamics simulations were carried out for αIIb-β3 complexes as well as αIIbβ3-fibrinogen complexes; in particular, E355K structure had more deviations, fluctuations, and other changes which compromised its structural stability and binding affinity when compared to both wild-type and G401C structures. Our comprehensive in silico analysis clearly reiterates that mutations in the β-propeller are not only responsible for structural changes in this domain but also have implications on the overall structure and function of integrin αIIbβ3.

Pathogenic cryptic variants detectable through exome data reanalysis significantly increase the diagnostic yield in Joubert syndrome

Joubert syndrome (JS) is a genetically heterogeneous neurodevelopmental ciliopathy. Despite exome sequencing (ES), several patients remain undiagnosed. This study aims to increase the diagnostic yield by uncovering cryptic variants through targeted ES reanalysis. We first focused on 26 patients in whom ES only disclosed heterozygous pathogenic coding variants in a JS gene. We reanalyzed raw ES data searching for copy number variants (CNVs) and intronic variants affecting splicing. We validated CNVs through real-time PCR or chromosomal microarray, and splicing variants through RT-PCR or minigenes. Cryptic variants were then searched in additional 44 ES-negative JS individuals. We identified cryptic "second hits" in 14 of 26 children (54%) and biallelic cryptic variants in 3 of 44 (7%), reaching a definite diagnosis in 17 of 70 (overall diagnostic gain 24%). We show that CNVs and intronic splicing variants are a common mutational mechanism in JS; more importantly, we demonstrate that a significant proportion of such variants can be disclosed simply through a focused reanalysis of available ES data, with a significantly increase of the diagnostic yield especially among patients previously found to carry heterozygous coding variants in the KIAA0586, CC2D2A and CPLANE1 genes.

Fabry disease in female monozygotic twins with complex intronic haplotype variants: a case report

BACKGROUND

Fabry disease is an X-linked lysosomal storage disease caused by the impairment of α-galactosidase A. The complex intronic haplotype (CIH) variants, located in promoter and intronic regulatory lesions, has been found in patients with classical forms of Fabry disease. We present a case of Fabry disease in female monozygotic twins exhibiting the CIH mutation and classical manifestations.

CASE PRESENTATION

A 61-year-old woman with a history of stroke, carotid artery occlusion, hypertrophic cardiomyopathy, and chronic kidney disease was referred to the nephrology clinic for management of her chronic kidney disease. Her monozygotic twin sister also presented with hypertrophic cardiomyopathy, atrial flutter, carotid stenosis, and proteinuria. Clinical symptoms and a comprehensive family history strongly suggested the presence of Fabry disease. Genetic analysis revealed the presence of 5 variants within a complex intronic haplotype (CIH): c.-10 C > T, c.369 + 990 C > A, c.370 - 81_370-77delCAGCC, c.640-16 A > G, and c.1000-22 C > T. We conducted a review of the patient's previous kidney biopsy findings, which demonstrated the presence of lamellated inclusion bodies in electron microscopy. Remarkably, both the monozygotic twin sister and her son exhibited the same genetic mutation. Enzyme replacement therapy was initiated for the patient. Her kidney function decreased throughout a thorough 2-year follow-up period, while there was a slight decrease in the left ventricular mass index.

CONCLUSIONS

This is the first reported case of female monozygotic twins with the CIH variants representing cardiac, cerebrovascular, and renal manifestations suggestive of Fabry disease.

Thirteen Indians with camptodactyly-arthropathy-coxa vara-pericarditis syndrome

Camptodactyly-arthropathy-coxa vara-pericarditis (CACP) syndrome (MIM# 208250) is a rare monogenic disorder, characterized by early onset of camptodactyly, progressive coxa vara, bilateral arthropathy and constrictive pericarditis. The syndrome is caused by biallelic loss-of-function variants in PRG4 . Deficiency of PRG4 results in progressive worsening of joint deformity with age. Thirteen individuals with CACP syndrome from eight consanguineous Indian families were evaluated. We used exome sequencing to elucidate disease-causing variants in all the probands. These variants were further validated and segregated by Sanger sequencing, confirming the diagnosis of CACP syndrome in them. Seven females and six males aged 2-23 years were studied. Camptodactyly (13/13), coxa vara (11/13), short femoral neck (11/13) and arthritis in large joints (12/13) [wrists (11/13), ankle (11/13), elbow (10/13) and knee (10/13)] were observed commonly. Five novel disease-causing variants (c.3636G>T, c.1935del, c.1134dup, c.1699del and c.962T>A) and two previously reported variants (c.1910_1911del and c.2816_2817del) were identified in homozygous state in PRG4 . We describe the phenotype and mutations in one of the large cohorts of patients with CACP syndrome, from India.

Global Presence and Penetrance of CSF1R-Related Disorder

Objectives

To highlight the worldwide presence of CSF1R-related disorder (CSF1R-RD), discuss its penetrance, and provide the first haplotype analysis.

Methods

Data on patients worldwide were collected, including demographics, genotype, family history, and clinical status. For haplotype analysis, polymorphisms of short tandem repeats in 3 distinct families with CSF1R p.Ile794Thr variant were examined.

Results

Nineteen new patients were included, at a mean age of 38.7 years (ranging from 11 to 74 years), from 14 families from the Americas, Asia, Australia, and Europe, including the first from Mexico, North Macedonia, and Ukraine. Fifteen CSF1R variants were found, including 8 novel. Three patients were compound heterozygotes with disease onset at 1, 4, and 22 years. Patients with heterozygous CSF1R variants developed symptoms at a mean of 39.0 years (range 8-71 years). Four patients died at a mean of 3.3 years from onset (range 2-5 years). Negative family history was noted in 7 patients. In haplotype analysis, 2 families exhibited shared haplotype encompassing ∼6-Mb region downstream of the CSF1R while the third family displayed a different haplotype.

Discussion

CSF1R-RD has a global prevalence. The reasons for negative family history include de novo variants (as shown by the haplotype analysis), mosaicism, and incomplete penetrance, which are possibly modulated by environmental and genetic factors.

Genomics 2 Proteins portal: a resource and discovery tool for linking genetic screening outputs to protein sequences and structures

Recent advances in AI-based methods have revolutionized the field of structural biology. Concomitantly, high-throughput sequencing and functional genomics have generated genetic variants at an unprecedented scale. However, efficient tools and resources are needed to link disparate data types-to 'map' variants onto protein structures, to better understand how the variation causes disease, and thereby design therapeutics. Here we present the Genomics 2 Proteins portal ( https://g2p.broadinstitute.org/ ): a human proteome-wide resource that maps 20,076,998 genetic variants onto 42,413 protein sequences and 77,923 structures, with a comprehensive set of structural and functional features. Additionally, the Genomics 2 Proteins portal allows users to interactively upload protein residue-wise annotations (for example, variants and scores) as well as the protein structure beyond databases to establish the connection between genomics to proteins. The portal serves as an easy-to-use discovery tool for researchers and scientists to hypothesize the structure-function relationship between natural or synthetic variations and their molecular phenotypes.

Mapping the Pharmacogenetic Landscape in a Ugandan Population: Implications for Personalized Medicine in an Underrepresented Population

Africans are extremely underrepresented in global genomic research. African populations face high burdens of communicable and non-communicable diseases and experience widespread polypharmacy. As population-specific genetic studies are crucial to understanding unique genetic profiles and optimizing treatments to reduce medication-related complications in this diverse population, the present study aims to characterize the pharmacogenomics profile of a rural Ugandan population. We analyzed low-pass whole genome sequencing data from 1998 Ugandans to investigate 18 clinically actionable pharmacogenes in this population. We utilized PyPGx to identify star alleles (haplotype patterns) and compared allele frequencies across populations using the Pharmacogenomics Knowledgebase PharmGKB. Clinical interpretations of the identified alleles were conducted following established dosing guidelines. Over 99% of participants displayed actionable phenotypes across the 18 pharmacogenes, averaging 3.5 actionable genotypes per individual. Several variant alleles known to affect drug metabolism (i.e., CYP3A5*1, CYP2B6*9, CYP3A5*6, CYP2D6*17, CYP2D6*29, and TMPT*3C)-which are generally more prevalent in African individuals-were notably enriched in the Ugandan cohort, beyond reported frequencies in other African peoples. More than half of the cohort exhibited a predicted impaired drug response associated with CFTR, IFNL3, CYP2B6, and CYP2C19, and approximately 31% predicted altered CYP2D6 metabolism. Potentially impaired CYP2C9, SLCO1B1, TPMT, and DPYD metabolic phenotypes were also enriched in Ugandans compared with other African populations. Ugandans exhibit distinct allele profiles that could impact drug efficacy and safety. Our findings have important implications for pharmacogenomics in Uganda, particularly with respect to the treatment of prevalent communicable and non-communicable diseases, and they emphasize the potential of pharmacogenomics-guided therapies to optimize healthcare outcomes and precision medicine in Uganda.

Comprehensive Proteomic Analysis of Dysferlinopathy Unveiling Molecular Mechanisms and Biomarkers Linked to Pathological Progression

AIMS

Previous proteomics studies in dysferlinopathy muscle have been limited in scope, often utilizing 2D-electrophoresis and yielding only a small number of differential expression calls. To address this gap, this study aimed to employ high-resolution proteomics to explore the proteomic landscapes of dysferlinopathy and analyze the correlation between muscle pathological changes and alterations in protein expression in muscle biopsies.

METHODS

We conducted a comprehensive approach to investigate the proteomic profile and disease-associated changes in the muscle tissue proteome from 15 patients with dysferlinopathy, exhibiting varying degrees of dystrophic pathology, alongside age-matched controls. Our methodology encompasses tandem mass tag (TMT)-labeled liquid chromatography-mass spectrometry (LC-MS/MS)-based proteomics, protein-protein interaction (PPI) network analysis, weighted gene co-expression network analysis, and differential expression analysis. Subsequently, we examined the correlation between the expression of key proteins and the clinical characteristics of the patients to identify pathogenic targets associated with DYSF mutations in dysferlinopathy.

RESULTS

A total of 1600 differentially expressed proteins were identified, with 1321 showing high expression levels and 279 expressed at lower levels. Our investigation yields a molecular profile delineating the altered protein networks in dysferlinopathy-afflicted skeletal muscle, uncovering dysregulation across numerous cellular pathways and molecular processes, including mRNA metabolic processes, regulated exocytosis, immune response, muscle system processes, energy metabolic processes, and calcium transmembrane transport. Moreover, we observe significant associations between the protein expression of ANXA1, ANXA2, ANXA4, ANXA5, LMNA, PYGM, and the extent of histopathologic changes in muscle biopsies from patients with dysferlinopathy, validated through immunoblotting and immunofluorescence assays.

CONCLUSIONS

Through the aggregation of expression data from dysferlinopathy-impacted muscles exhibiting a range of pathological alterations, we identified multiple key proteins associated with the dystrophic pathology of patients with dysferlinopathy. These findings provide novel insights into the pathogenesis of dysferlinopathy and propose promising targets for future therapeutic endeavors.

Preclinical alternative drug discovery programs for monogenic rare diseases. Should small molecules or gene therapy be used? The case of hereditary spastic paraplegias

Patients diagnosed with rare diseases and their and families search desperately to organize drug discovery campaigns. Alternative models that differ from default paradigms offer real opportunities. There are, however, no clear guidelines for the development of such models, which reduces success rates and raises costs. We address the main challenges in making the discovery of new preclinical treatments more accessible, using rare hereditary paraplegia as a paradigmatic case. First, we discuss the necessary expertise, and the patients' clinical and genetic data. Then, we revisit gene therapy, de novo drug development, and drug repurposing, discussing their applicability. Moreover, we explore a pool of recommended in silico tools for pathogenic variant and protein structure prediction, virtual screening, and experimental validation methods, discussing their strengths and weaknesses. Finally, we focus on successful case applications.

Investigating therapeutic nonsense suppression in a neurofibromatosis mouse model

Neurofibromatosis type 1 (NF1) is a human genetic disorder caused by variants in the NF1 gene. Plexiform neurofibromas, one of many NF1 manifestations, are benign peripheral nerve sheath tumors occurring in up to 50% of NF1 patients. A substantial fraction of NF1 pathogenetic variants are nonsense mutations, which result in the synthesis of truncated non-functional NF1 protein (neurofibromin). To date, no therapeutics have restored neurofibromin expression or addressed the consequences of this protein's absence in NF1 nonsense mutation patients, but nonsense suppression is a potential approach to the problem. Ataluren is a small molecule drug that has been shown to stimulate functional nonsense codon readthrough in several models of nonsense mutation diseases, as well as in Duchenne muscular dystrophy patients. To test ataluren's potential applicability in nonsense mutation NF1 patients, we evaluated its therapeutic effects using three treatment regimens in a previously established NF1 patient-derived (c.2041C > T; p.Arg681X) nonsense mutation mouse model. Collectively, our experiments indicate that: i) ataluren appeared to slow the growth of neurofibromas and alleviate some paralysis phenotypes, ii) female Nf1-nonsense mutation mice manifested more severe paralysis and neurofibroma phenotypes than male mice, iii) ataluren doses with apparent effectiveness were lower in female mice than in male mice, and iv) age factors also influenced ataluren's effectiveness.

The co-occurrence of genetic variants in the TYR and OCA2 genes confers susceptibility to albinism

Although rare genetic conditions are mostly caused by DNA sequence alterations that functionally disrupt individual genes, large-scale studies using genome sequencing have started to unmask additional complexity. Understanding how combinations of variants in different genes shape human phenotypes is expected to provide important insights into the clinical and genetic heterogeneity of rare disorders. Here, we use albinism, an archetypal rare condition associated with hypopigmentation, as an exemplar for the study of genetic interactions. We analyse data from the Genomics England 100,000 Genomes Project alongside a cohort of 1120 individuals with albinism, and investigate the effect of dual heterozygosity for the combination of two established albinism-related variants: TYR:c.1205 G > A (p.Arg402Gln) [rs1126809] and OCA2:c.1327 G > A (p.Val443Ile) [rs74653330]. As each of these changes alone is insufficient to cause disease when present in the heterozygous state, we sought evidence of synergistic effects. We show that, when both variants are present, the probability of receiving a diagnosis of albinism is significantly increased (odds ratio 12.8; 95% confidence interval 6.0 - 24.7; p-value 2.1 ×10-8). Further analyses in an independent cohort, the UK Biobank, support this finding and highlight that heterozygosity for the TYR:c.1205 G > A and OCA2:c.1327 G > A variant combination is associated with statistically significant alterations in visual acuity and central retinal thickness (traits that are considered albinism endophenotypes). The approach discussed in this report opens up new avenues for the investigation of oligogenic patterns in apparently Mendelian disorders.

An intronic variant in LAMB3 contributes to junctional epidermolysis bullosa and enamel hypoplasia via translational attenuation

OBJECTIVES

This study aimed to investigate the genetic etiology of a family affected by junctional epidermolysis bullosa (JEB) and generalized enamel hypoplasia, and to explore how an intronic variant influenced the 5' untranslated region (5'UTR), thereby affecting LAMB3 expression and contributing to the pathogenesis of the disease.

DESIGN

Whole-exome and whole-genome sequencing were used to screen for genetic defects in the patient. Mutational consequences were characterized through luciferase assays, splice assay, in silico analyses, and verification using the patient's gingival sample.

RESULTS

A nonsense variant (c.2983 C>T; p.Gln995*) and an intronic variant (c.-38+2 T>C) of LAMB3 were identified. In vitro assays demonstrated that the intronic variant activated a cryptic splice site, resulting in a 120 bp intronic inclusion. This splicing alteration significantly reduced the translation efficiency of the downstream coding sequence, while overall mRNA expression remained unaffected. Bioinformatic analysis unveiled the creation of three upstream AUG codons, leading to the presence of two upstream open reading frames (uORFs) and one overlapping ORF. The longer uORF's AUG exhibited a moderate Kozak strength similar to that of the main ORF's AUG. Structural analysis of the mutant 5'UTR sequence revealed a more complex secondary structure, characterized by a large branch loop and a stem-loop preceding the coding sequence's start codon.

CONCLUSION

This study suggests that variants affecting the 5'UTR may contribute to the genetic etiology of JEB. These findings could help enhance the diagnostic accuracy and efficiency in JEB patients.

Decoding Missense Variants by Incorporating Phase Separation via Machine Learning

Computational models have made significant progress in predicting the effect of protein variants. However, deciphering numerous variants of uncertain significance (VUS) located within intrinsically disordered regions (IDRs) remains challenging. To address this issue, we introduce phase separation, which is tightly linked to IDRs, into the investigation of missense variants. Phase separation is vital for multiple physiological processes. By leveraging missense variants that alter phase separation propensity, we develop a machine learning approach named PSMutPred to predict the impact of missense mutations on phase separation. PSMutPred demonstrates robust performance in predicting missense variants that affect natural phase separation. In vitro experiments further underscore its validity. By applying PSMutPred on over 522,000 ClinVar missense variants, it significantly contributes to decoding the pathogenesis of disease variants, especially those in IDRs. Our work provides insights into the understanding of a vast number of VUSs in IDRs, expediting clinical interpretation and diagnosis.

Novel variants and genotype-phenotype correlation in a multicentre cohort of GNE myopathy in China

BACKGROUND

GlcNAc2-epimerase (GNE) myopathy is a rare autosomal recessive disorder caused by pathogenic variants in the GNE gene, which is essential for the sialic acid biosynthesis pathway.

OBJECTIVE

This multi-centre study aimed to delineate the clinical phenotype and GNE variant spectrum in Chinese patients, enhancing our understanding of the genetic diversity and clinical manifestation across different populations.

METHODS

We retrospectively analysed GNE variants from 113 patients, integrating these data with external GNE variants from online databases for a global perspective, examining their consequences, distribution, ethnicity and severity.

RESULTS

This study revealed 97 distinct GNE variants, including 35 (36.08%) novel variants. Two more patients with deep intronic variant c.862+870C>T were identified, while whole genome sequencing (WGS) uncovered another two novel intronic variants: c.52-8924G>T and c.1505-12G>A. Nanopore long reads sequencing (LRS) and further PCR analysis verified a 639 bp insertion at chr9:36249241. Missense variants predominantly located in the epimerase/kinase domain coding region, indicating the impairment of catalytic function as a key pathogenic consequence. Comparative studies with Japanese, Korean and Jewish, our cohorts showed later onset ages by 2 years. The high allele frequency of the non-catalytic GNE variant, c.620A>T, might underlie the milder phenotype of Chinese patients.

CONCLUSIONS

Comprehensive techniques such as WGS and Nanopore LRS warrants the identifying of GNE variants. Patients with the non-catalytic GNE variant, c.620A>T, had a milder disease progression and later wheelchair use.

Generative language models on nucleotide sequences of human genes

Language models, especially transformer-based ones, have achieved colossal success in natural language processing. To be precise, studies like BERT for natural language understanding and works like GPT-3 for natural language generation are very important. If we consider DNA sequences as a text written with an alphabet of four letters representing the nucleotides, they are similar in structure to natural languages. This similarity has led to the development of discriminative language models such as DNABERT in the field of DNA-related bioinformatics. To our knowledge, however, the generative side of the coin is still largely unexplored. Therefore, we have focused on the development of an autoregressive generative language model such as GPT-3 for DNA sequences. Since working with whole DNA sequences is challenging without extensive computational resources, we decided to conduct our study on a smaller scale and focus on nucleotide sequences of human genes, i.e. unique parts of DNA with specific functions, rather than the whole DNA. This decision has not significantly changed the structure of the problem, as both DNA and genes can be considered as 1D sequences consisting of four different nucleotides without losing much information and without oversimplification. First of all, we systematically studied an almost entirely unexplored problem and observed that recurrent neural networks (RNNs) perform best, while simple techniques such as N-grams are also promising. Another beneficial point was learning how to work with generative models on languages we do not understand, unlike natural languages. The importance of using real-world tasks beyond classical metrics such as perplexity was noted. In addition, we examined whether the data-hungry nature of these models can be altered by selecting a language with minimal vocabulary size, four due to four different types of nucleotides. The reason for reviewing this was that choosing such a language might make the problem easier. However, in this study, we found that this did not change the amount of data required very much.

Pervasive mislocalization of pathogenic coding variants underlying human disorders

Widespread sequencing has yielded thousands of missense variants predicted or confirmed as disease causing. This creates a new bottleneck: determining the functional impact of each variant-typically a painstaking, customized process undertaken one or a few genes and variants at a time. Here, we established a high-throughput imaging platform to assay the impact of coding variation on protein localization, evaluating 3,448 missense variants of over 1,000 genes and phenotypes. We discovered that mislocalization is a common consequence of coding variation, affecting about one-sixth of all pathogenic missense variants, all cellular compartments, and recessive and dominant disorders alike. Mislocalization is primarily driven by effects on protein stability and membrane insertion rather than disruptions of trafficking signals or specific interactions. Furthermore, mislocalization patterns help explain pleiotropy and disease severity and provide insights on variants of uncertain significance. Our publicly available resource extends our understanding of coding variation in human diseases.

NeuroBooster Array: A Genome-Wide Genotyping Platform to Study Neurological Disorders Across Diverse Populations

BACKGROUND

Commercial genome-wide genotyping arrays have historically neglected coverage of genetic variation across populations.

OBJECTIVE

We aimed to create a multi-ancestry genome-wide array that would include a wide range of neuro-specific genetic content to facilitate genetic research in neurological disorders across multiple ancestral groups, fostering diversity and inclusivity in research studies.

METHODS

We developed the Illumina NeuroBooster Array (NBA), a custom high-throughput and cost-effective platform on a backbone of 1,914,934 variants from the Infinium Global Diversity Array and added custom content comprising 95,273 variants associated with more than 70 neurological conditions or traits, and we further tested its performance on more than 2000 patient samples. This novel platform includes approximately 10,000 tagging variants to facilitate imputation and analyses of neurodegenerative disease-related genome-wide association study loci across diverse populations.

RESULTS

In this article, we describe NBA's potential as an efficient means for researchers to assess known and novel disease genetic associations in a multi-ancestry framework. The NBA can identify rare genetic variants and accurately impute more than 15 million common variants across populations. Apart from enabling sample prioritization for further whole-genome sequencing studies, we envisage that NBA will play a pivotal role in recruitment for interventional studies in the precision medicine space.

CONCLUSIONS

From a broader perspective, the NBA serves as a promising means to foster collaborative research endeavors in the field of neurological disorders worldwide. Ultimately, this carefully designed tool is poised to make a substantial contribution to uncovering the genetic etiology underlying these debilitating conditions. © 2024 The Author(s). Movement Disorders published by Wiley Periodicals LLC on behalf of International Parkinson and Movement Disorder Society. This article has been contributed to by U.S. Government employees and their work is in the public domain in the USA.

Genetic polymorphism in untranslated regions of PRKCZ influences mRNA structure, stability and binding sites

BACKGROUND

Variations in untranslated regions (UTR) alter regulatory pathways impacting phenotype, disease onset, and course of disease. Protein kinase C Zeta (PRKCZ), a serine-threonine kinase, is implicated in cardiovascular, neurological and oncological disorders. Due to limited research on PRKCZ, this study aimed to investigate the impact of UTR genetic variants' on binding sites for transcription factors and miRNA. RNA secondary structure, eQTLs, and variation tolerance analysis were also part of the study.

METHODS

The data related to PRKCZ gene variants was downloaded from the Ensembl genome browser, COSMIC and gnomAD. The RegulomeDB database was used to assess the functional impact of 5' UTR and 3'UTR variants. The analysis of the transcription binding sites (TFBS) was done through the Alibaba tool, and the Kyoto Encyclopaedia of Genes and Genomes (KEGG) was employed to identify pathways associated with PRKCZ. To predict the effect of variants on microRNA binding sites, PolymiRTS was utilized for 3' UTR variants, and the SNPinfo tool was used for 5' UTR variants.

RESULTS

The results obtained indicated that a total of 24 variants present in the 3' UTR and 25 variants present in the 5' UTR were most detrimental. TFBS analysis revealed that 5' UTR variants added YY1, repressor, and Oct1, whereas 3' UTR variants added AP-2alpha, AhR, Da, GR, and USF binding sites. The study predicted TFs that influenced PRKCZ expression. RNA secondary structure analysis showed that eight 5' UTR and six 3' UTR altered the RNA structure by either removal or addition of the stem-loop. The microRNA binding site analysis highlighted that seven 3' UTR and one 5' UTR variant altered the conserved site and also created new binding sites. eQTLs analysis showed that one variant was associated with PRKCZ expression in the lung and thyroid. The variation tolerance analysis revealed that PRKCZ was an intolerant gene.

CONCLUSION

This study laid the groundwork for future studies aimed at targeting PRKCZ as a therapeutic target.

Mitochondrial bioenergetics and cardiolipin remodeling abnormalities in mitochondrial trifunctional protein deficiency

Mitochondrial trifunctional protein (TFP) deficiency is an inherited metabolic disorder leading to a block in long-chain fatty acid β-oxidation. Mutations in HADHA and HADHB, which encode the TFP α and β subunits, respectively, usually result in combined TFP deficiency. A single common mutation, HADHA c.1528G>C (p.E510Q), leads to isolated 3-hydroxyacyl-CoA dehydrogenase deficiency. TFP also catalyzes a step in the remodeling of cardiolipin (CL), a phospholipid critical to mitochondrial membrane stability and function. We explored the effect of mutations in TFP subunits on CL and other phospholipid content and composition and the consequences of these changes on mitochondrial bioenergetics in patient-derived fibroblasts. Abnormalities in these parameters varied extensively among different fibroblasts, and some cells were able to maintain basal oxygen consumption rates similar to controls. Although CL reduction was universally identified, a simultaneous increase in monolysocardiolipins was discrepant among cells. A similar profile was seen in liver mitochondria isolates from a TFP-deficient mouse model. Response to new potential drugs targeting CL metabolism might be dependent on patient genotype.

Genome-scale quantification and prediction of pathogenic stop codon readthrough by small molecules

Premature termination codons (PTCs) cause ~10-20% of inherited diseases and are a major mechanism of tumor suppressor gene inactivation in cancer. A general strategy to alleviate the effects of PTCs would be to promote translational readthrough. Nonsense suppression by small molecules has proven effective in diverse disease models, but translation into the clinic is hampered by ineffective readthrough of many PTCs. Here we directly tackle the challenge of defining drug efficacy by quantifying the readthrough of ~5,800 human pathogenic stop codons by eight drugs. We find that different drugs promote the readthrough of complementary subsets of PTCs defined by local sequence context. This allows us to build interpretable models that accurately predict drug-induced readthrough genome-wide, and we validate these models by quantifying endogenous stop codon readthrough. Accurate readthrough quantification and prediction will empower clinical trial design and the development of personalized nonsense suppression therapies.

Congenital anterior segment ocular disorders: Genotype-phenotype correlations and emerging novel mechanisms

Development of the anterior segment of the eye requires reciprocal sequential interactions between the arising tissues, facilitated by numerous genetic factors. Disruption of any of these processes results in congenital anomalies in the affected tissue(s) leading to anterior segment disorders (ASD) including aniridia, Axenfeld-Rieger anomaly, congenital corneal opacities (Peters anomaly, cornea plana, congenital primary aphakia), and primary congenital glaucoma. Current understanding of the genetic factors involved in ASD remains incomplete, with approximately 50% overall receiving a genetic diagnosis. While some genes are strongly associated with a specific clinical diagnosis, the majority of known factors are linked with highly variable phenotypic presentations, with pathogenic variants in FOXC1, CYP1B1, and PITX2 associated with the broadest spectrum of ASD conditions. This review discusses typical clinical presentations including associated systemic features of various forms of ASD; the latest functional data and genotype-phenotype correlations related to 25 ASD factors including newly identified genes; promising novel candidates; and current and emerging treatments for these complex conditions. Recent developments of interest in the genetics of ASD include identification of phenotypic expansions for several factors, discovery of multiple modes of inheritance for some genes, and novel mechanisms including a growing number of non-coding variants and alleles affecting specific domains/residues and requiring further studies.

Diverse biophysical mechanisms in voltage-gated sodium channel Nav1.4 variants associated with myotonia

Mutations in SCN4A gene encoding Nav1.4 channel α-subunit, are known to cause neuromuscular disorders such as myotonia or paralysis. Here, we study the effect of two amino acid replacements, K1302Q and G1306E, in the DIII-IV loop of the channel, corresponding to mutations found in patients with myotonia. We combine clinical, electrophysiological, and molecular modeling data to provide a holistic picture of the molecular mechanisms operating in mutant channels and eventually leading to pathology. We analyze the existing clinical data for patients with the K1302Q substitution, which was reported for adults with or without myotonia phenotypes, and report two new unrelated patients with the G1306E substitution, who presented with severe neonatal episodic laryngospasm and childhood-onset myotonia. We provide a functional analysis of the mutant channels by expressing Nav1.4 α-subunit in Xenopus oocytes in combination with β1 subunit and recording sodium currents using two-electrode voltage clamp. The K1302Q variant exhibits abnormal voltage dependence of steady-state fast inactivation, being the likely cause of pathology. K1302Q does not lead to decelerated fast inactivation, unlike several other myotonic mutations such as G1306E. For both mutants, we observe increased window currents corresponding to a larger population of channels available for activation. To elaborate the structural rationale for our experimental data, we explore the contacts involving K/Q1302 and E1306 in the AlphaFold2 model of wild-type Nav1.4 and Monte Carlo-minimized models of mutant channels. Our data provide the missing evidence to support the classification of K1302Q variant as likely pathogenic and may be used by clinicians.

Variant Impact Predictor database (VIPdb), version 2: trends from three decades of genetic variant impact predictors

BACKGROUND

Variant interpretation is essential for identifying patients' disease-causing genetic variants amongst the millions detected in their genomes. Hundreds of Variant Impact Predictors (VIPs), also known as Variant Effect Predictors (VEPs), have been developed for this purpose, with a variety of methodologies and goals. To facilitate the exploration of available VIP options, we have created the Variant Impact Predictor database (VIPdb).

RESULTS

The Variant Impact Predictor database (VIPdb) version 2 presents a collection of VIPs developed over the past three decades, summarizing their characteristics, ClinGen calibrated scores, CAGI assessment results, publication details, access information, and citation patterns. We previously summarized 217 VIPs and their features in VIPdb in 2019. Building upon this foundation, we identified and categorized an additional 190 VIPs, resulting in a total of 407 VIPs in VIPdb version 2. The majority of the VIPs have the capacity to predict the impacts of single nucleotide variants and nonsynonymous variants. More VIPs tailored to predict the impacts of insertions and deletions have been developed since the 2010s. In contrast, relatively few VIPs are dedicated to the prediction of splicing, structural, synonymous, and regulatory variants. The increasing rate of citations to VIPs reflects the ongoing growth in their use, and the evolving trends in citations reveal development in the field and individual methods.

CONCLUSIONS

VIPdb version 2 summarizes 407 VIPs and their features, potentially facilitating VIP exploration for various variant interpretation applications. VIPdb is available at  https://genomeinterpretation.org/vipdb.

Genotypic spectrum of ABCA4-associated retinal degenerations in 211 unrelated Mexican patients: identification of 22 novel disease-causing variants

The purpose of this study was to analyze and molecularly describe the largest group of patients with ABCA4-associated retinal degeneration in Latin America. Pathogenic variants in ABCA4, a member of the ATP Binding Cassette (ABC) transporters superfamily, is one of the most common causes of inherited visual deficiency in humans. Retinal phenotypes associated with genetic defects in ABCA4 are collectively known as ABCA4-associated retinal degenerations (ABCA4R), a group of recessively inherited disorders associated with a high allelic heterogeneity. While large groups of Caucasian and Asiatic individuals suffering from ABCA4R have been well characterized, molecular information from certain ethnic groups is limited or unavailable, precluding a more realistic knowledge of ABCA4-related mutational profile worldwide. In this study, we describe the molecular findings of a large group of 211 ABCA4R index cases from Mexico. Genotyping was performed using either next generation sequencing (NGS) of a retinal dystrophy genes panel or exome. ABCA4 targeted mutation testing was applied to a subgroup of subjects in whom founder mutations were suspected. A total of 128 different ABCA4 pathogenic variants were identified, including 22 previously unpublished variants. The most common type of genetic variation was single nucleotide substitutions which occurred in 92.7% (408/440 alleles). According to the predicted protein effect, the most frequent variant type was missense, occurring in 83.5% of disease-causing alleles (368/440). Mutations such as p.Ala1773Val are fully demonstrated as a founder effect in native inhabitants of certain regions of Mexico. This study also gives us certain indications of other founder effects that need to be further studied in the near future. This is the largest molecularly characterized ABCA4R Latin American cohort, and our results supports the value of conducting genetic screening in underrepresented populations for a better knowledge of the mutational profile leading to monogenic diseases.

Fast and deep phosphoproteome analysis with the Orbitrap Astral mass spectrometer

Owing to its roles in cellular signal transduction, protein phosphorylation plays critical roles in myriad cell processes. That said, detecting and quantifying protein phosphorylation has remained a challenge. We describe the use of a novel mass spectrometer (Orbitrap Astral) coupled with data-independent acquisition (DIA) to achieve rapid and deep analysis of human and mouse phosphoproteomes. With this method, we map approximately 30,000 unique human phosphorylation sites within a half-hour of data collection. The technology is benchmarked to other state-of-the-art MS platforms using both synthetic peptide standards and with EGF-stimulated HeLa cells. We apply this approach to generate a phosphoproteome multi-tissue atlas of the mouse. Altogether, we detect 81,120 unique phosphorylation sites within 12 hours of measurement. With this unique dataset, we examine the sequence, structural, and kinase specificity context of protein phosphorylation. Finally, we highlight the discovery potential of this resource with multiple examples of phosphorylation events relevant to mitochondrial and brain biology.

Major Causes of Conflicting Interpretations of Variant Pathogenicity in Rare Disease: A Systematic Analysis

The identification of the genetic causes of inherited disorders from next-generation sequencing (NGS) data remains a complicated process, in particular due to challenges in interpretation of the vast amount of generated data and hundreds of candidate variants identified. Inconsistencies in variant classification, where genetic centers classify the same variant differently, can hinder accurate diagnoses for rare diseases. Publicly available databases that collect data on human genetic variations and their association with diseases provide ample opportunities to discover conflicts in variant interpretation worldwide. In this study, we explored patterns of variant classification discrepancies using data from ClinVar, a public archive of variant interpretations. We found that 5.7% of variants have conflicting interpretations (COIs) reported, and the vast majority of interpretation conflicts arise for variants of uncertain significance (VUS). As many as 78% of clinically relevant genes harbor variants with COIs, and genes with high COI rates tended to have more exons and longer transcripts, with a greater proportion of genes linked to several distinct conditions. The enrichment analysis of COI-enriched genes revealed that the products of these genes are involved in cardiac disorders, muscle development, and function. To improve diagnoses, we believe that specific variant interpretation rules could be developed for such genes. Additionally, our findings underscore the need for the publication of variant pathogenicity evidence and the importance of considering every variant as VUS unless proven otherwise.

Case report: First diagnosis of Fabry disease in North Macedonia in a patient presenting with kidney failure on hemodialysis

Introduction

Fabry disease is a rare X-linked lysosomal storage disorder caused by α-galactosidase A (α-Gal A) deficiency. Reduced or absent enzyme activity causes progressive lysosomal accumulation of globotriaosylceramide (Lyso-Gb3) in various cells throughout the body to trigger inflammation and fibrosis.

Case description

We present the first familial case of Fabry Disease in North Macedonia identified based on clinical manifestations and confirmed through enzyme, biomarker, and genetic tests. The index case in the family was a 45-year-old male undergoing hemodialysis therapy. He has had chronic burning uncontrolled limb pain since childhood, intermittent abdominal cramps, anhidrosis, and hypertension. The constellation of clinical presentations accompanied by similar symptoms in close family members prompted the enzyme, biomarker, and genetic analyses for Fabry disease. Genetic testing identified a known pathogenic GLA missense variant c.443G>A or p.(Ser148Asn) in the hemizygous state. Subsequent family studies allowed identification of another hemizygous male and five heterozygous female carriers affected by this X-linked disorder.

Conclusion

We report identification of the first familial case of Fabry disease in North Macedonia and describe the phenotype associated with the Ser148Asn GLA variant. Greater awareness of this rare disease linked to continuous medical education is crucial for timely diagnosis and treatment.

Principles of clinical genetics for rheumatologists: clinical indications and interpretation of broad-based genetic testing

Advances in DNA sequencing technologies, especially next-generation sequencing (NGS), which is the basis for whole-exome sequencing (WES) and whole-genome sequencing (WGS), have profoundly transformed immune-mediated rheumatic disease diagnosis. Recently, substantial cost reductions have facilitated access to these diagnostic tools, expanded the capacity of molecular diagnostics and enabled the pursuit of precision medicine in rheumatology. Understanding the fundamental principles of genetics and diversity in genetic variant classification is a crucial milestone in rheumatology. However, despite the growing availability of DNA sequencing platforms, a significant number of autoinflammatory diseases (AIDs), neuromuscular disorders, hereditary collagen diseases, and monogenic bone diseases remain unsolved, and variants of uncertain significance (VUS) pose a formidable challenge to addressing these unmet needs in the coming decades. This article aims to provide an overview of the clinical indications and interpretation of comprehensive genetic testing in the medical field, addressing the related complexities and implications.

Potentiation by Protein Synthesis Inducers of Translational Readthrough of Pathogenic Premature Termination Codons in PTEN Isoforms

The PTEN tumor suppressor is frequently targeted in tumors and patients with PTEN hamartoma tumor syndrome (PHTS) through nonsense mutations generating premature termination codons (PTC) that may cause the translation of truncated non-functional PTEN proteins. We have previously described a global analysis of the readthrough reconstitution of the protein translation and function of the human canonical PTEN isoform by aminoglycosides. Here, we report the efficient functional readthrough reconstitution of the PTEN translational isoform PTEN-L, which displays a minimal number of PTC in its specific N-terminal extension in association with disease. We illustrate the importance of the specific PTC and its nucleotide proximal sequence for optimal readthrough and show that the more frequent human PTEN PTC variants and their mouse PTEN PTC equivalents display similar patterns of readthrough efficiency. The heterogeneous readthrough response of the different PTEN PTC variants was independent of the length of the PTEN protein being reconstituted, and we found a correlation between the amount of PTEN protein being synthesized and the PTEN readthrough efficiency. Furthermore, combination of aminoglycosides and protein synthesis inducers increased the readthrough response of specific PTEN PTC. Our results provide insights with which to improve the functional reconstitution of human-disease-related PTC pathogenic variants from PTEN isoforms by increasing protein synthesis coupled to translational readthrough.

RNA variant assessment using transactivation and transdifferentiation

Understanding the impact of splicing and nonsense variants on RNA is crucial for the resolution of variant classification as well as their suitability for precision medicine interventions. This is primarily enabled through RNA studies involving transcriptomics followed by targeted assays using RNA isolated from clinically accessible tissues (CATs) such as blood or skin of affected individuals. Insufficient disease gene expression in CATs does however pose a major barrier to RNA based investigations, which we show is relevant to 1,436 Mendelian disease genes. We term these "silent" Mendelian genes (SMGs), the largest portion (36%) of which are associated with neurological disorders. We developed two approaches to induce SMG expression in human dermal fibroblasts (HDFs) to overcome this limitation, including CRISPR-activation-based gene transactivation and fibroblast-to-neuron transdifferentiation. Initial transactivation screens involving 40 SMGs stimulated our development of a highly multiplexed transactivation system culminating in the 6- to 90,000-fold induction of expression of 20/20 (100%) SMGs tested in HDFs. Transdifferentiation of HDFs directly to neurons led to expression of 193/516 (37.4%) of SMGs implicated in neurological disease. The magnitude and isoform diversity of SMG expression following either transactivation or transdifferentiation was comparable to clinically relevant tissues. We apply transdifferentiation and/or gene transactivation combined with short- and long-read RNA sequencing to investigate the impact that variants in USH2A, SCN1A, DMD, and PAK3 have on RNA using HDFs derived from affected individuals. Transactivation and transdifferentiation represent rapid, scalable functional genomic solutions to investigate variants impacting SMGs in the patient cell and genomic context.

Using computational approaches to enhance the interpretation of missense variants in the PAX6 gene

The PAX6 gene encodes a highly-conserved transcription factor involved in eye development. Heterozygous loss-of-function variants in PAX6 can cause a range of ophthalmic disorders including aniridia. A key molecular diagnostic challenge is that many PAX6 missense changes are presently classified as variants of uncertain significance. While computational tools can be used to assess the effect of genetic alterations, the accuracy of their predictions varies. Here, we evaluated and optimised the performance of computational prediction tools in relation to PAX6 missense variants. Through inspection of publicly available resources (including HGMD, ClinVar, LOVD and gnomAD), we identified 241 PAX6 missense variants that were used for model training and evaluation. The performance of ten commonly used computational tools was assessed and a threshold optimization approach was utilized to determine optimal cut-off values. Validation studies were subsequently undertaken using PAX6 variants from a local database. AlphaMissense, SIFT4G and REVEL emerged as the best-performing predictors; the optimized thresholds of these tools were 0.967, 0.025, and 0.772, respectively. Combining the prediction from these top-three tools resulted in lower performance compared to using AlphaMissense alone. Tailoring the use of computational tools by employing optimized thresholds specific to PAX6 can enhance algorithmic performance. Our findings have implications for PAX6 variant interpretation in clinical settings.

Store-operated calcium entry dysfunction in CRAC channelopathy: Insights from a novel STIM1 mutation

Store-operated calcium entry (SOCE) plays a crucial role in maintaining cellular calcium homeostasis. This mechanism involves proteins, such as stromal interaction molecule 1 (STIM1) and ORAI1. Mutations in the genes encoding these proteins, especially STIM1, can lead to various diseases, including CRAC channelopathies associated with severe combined immunodeficiency. Herein, we describe a novel homozygous mutation, NM_003156 c.792-3C > G, in STIM1 in a patient with a clinical profile of CRAC channelopathy, including immune system deficiencies and muscle weakness. Functional analyses revealed three distinct spliced forms in the patient cells: wild-type, exon 7 skipping, and intronic retention. Calcium influx analysis revealed impaired SOCE in the patient cells, indicating a loss of STIM1 function. We developed an antisense oligonucleotide treatment that improves STIM1 splicing and highlighted its potential as a therapeutic approach. Our findings provide insights into the complex effects of STIM1 mutations and shed light on the multifaceted clinical presentation of the patient.

A novel LTBP2 gene variant in a Turkish family with juvenile-onset open-angle glaucoma

BACKGROUND

Juvenile-onset open-angle glaucoma (JOAG) is a rare form of primary open-angle glaucoma (POAG) with an early age of onset before 40 years. Latent transforming growth factor-beta binding protein 2 (LTBP-2) is an extracellular matrix protein with a multi-domain structure and homology to fibrillins. LTBP2 gene variants have been associated with JOAG in a small number of patients. Herein, we report a novel missense variant in the LTBP2 gene in a Turkish family with JOAG.

MATERIALS AND METHODS

Blood samples were obtained from three siblings (a 20-year-old woman with JOAG, 26-year-old man with JOAG, and 15-year-old girl with posterior embryotoxon) for genetic analysis. Their father had moderate-severe POAG and the 24-year-old brother had JOAG. The mother and 32-year-old sister were healthy. Although the parents reported no consanguinity, they come from the same village.

RESULTS

Clinical exome sequencing analysis of the two siblings with JOAG revealed a novel c.607C>T p.(R203C) (rs777450651) homozygous LTBP2 variant, while the variant was heterozygous in their 15-year-old sister. There were no mutations in the MYOC, CYP1B1, or FBN1 genes.

CONCLUSION

We documented a novel missense mutation in the LTBP2 gene leading to a severe form of JOAG with refractory IOP and progressive optic nerve damage, which seems to show autosomal recessive inheritance.

Elucidating the functional impact of G137V and G144R variants in Maroteaux Lamy's Syndrome by Molecular Dynamics Simulation

Mucopolysaccharidoses VI (Maroteaux Lamy syndrome) is a metabolic disorder due to the loss of enzyme activity of N-acetyl galactosamine-4-sulphatase arising from mutations in the ARSB gene. The mutated ARSB is the origin for the accumulation of GAGs within the lysosome leading to severe growth deformities, causing lysosomal storage disease. The main focus of this study is to identify the deleterious variants by applying bioinformatics tools to predict the conservation, pathogenicity, stability, and effect of the ARSB variants. We examined 170 missense variants, of which G137V and G144R were the resultant variants predicted detrimental to the progression of the disease. The native along with G137V and G144R structures were fixed as the receptors and subjected to Molecular docking with the small molecule Odiparcil to analyze the binding efficiency and the varied interactions of the receptors towards the drug. The interaction resulted in similar docking scores of - 7.3 kcal/mol indicating effective binding and consistent interactions of the drug with residues CYS117, GLN118, THR182, and GLN517 for native, along with G137V and G144R structures. Molecular Dynamics were conducted to validate the stability and flexibility of the native and variant structures on ligand binding. The overall study indicates that the drug has similar therapeutic towards the native and variant based on the higher binding affinity and also the complexes show stability with an average of 0.2 nm RMS value. This can aid in the future development therapeutics for the Maroteaux Lamy syndrome.

The Systemic Genotype-Phenotype Characterization of PAX6-Related Eye Disease in 164 Chinese Families

Purpose

This study aims to evaluate the genetic and phenotypic characteristics and elucidate the genotype-phenotype correlations of a large Chinese cohort with PAX6-related disorders.

Methods

Variants detected with exome sequencing were filtered through multistep bioinformatic and co-segregation analyses, and validated by Sanger sequencing. The related clinical data were collected, and cluster analysis and statistical analysis of the PAX6-related phenotypes across different variant groups were carried out. Parental mosaicism was investigated using cloning analysis and Droplet digital PCR.

Results

A total of 119 pathogenic or likely pathogenic PAX6 variants, including 74 truncation, 31 missense, and 14 others, were identified in 228 patients from 164 unrelated families. The most common phenotypes were foveal hypoplasia (97.8%), nystagmus (92.6%), aniridia (76.7%), cataract (36.8%), and iris hypoplasia (22.4%). Mosaicism ranging from 13.9% to 18.8% was identified in 3 unrelated patients' parents with relatively mild phenotypes. Missense variants in the linker region of the paired domain were associated with high myopia, whereas truncation variants in the homeodomain and proline-serine-threonine-rich domain were associated with hyperopia. Similarly, the degree of iris defects, visual acuity, and associated ocular comorbidity varied among the different types and locations of PAX6 variants.

Conclusions

Our data indicate that foveal hypoplasia but not aniridia is the most common sign of PAX6-related disorders, contributing to subtle iris changes that might easily be overlooked in clinical practice. Recognition of mosaicism in atypical cases or parents with very mild phenotypes is important in genetic counseling as their offspring are at increased risk of typical aniridia. Recognition of the genotype-phenotype relationship emphasizes involvement of PAX6 regulation in shaping complex ocular phenotypes.

Undetected Neuromuscular Disease in Patients after Heart Transplantation

(1) Heart transplantation (HTX) improves the overall survival and functional status of end-stage heart failure patients with cardiomyopathies (CMPs). The majority of CMPs have genetic causes, and the overlap between CMPs and inherited myopathies is well documented. However, the long-term outcome in skeletal muscle function and possibility of an undiagnosed underlying genetic cause of both a cardiac and skeletal pathology remain unknown. (2) Thirty-nine patients were assessed using open and standardized interviews on muscle function, a quality-of-life (EuroQol EQ-5D-3L) questionnaire, and a physical examination (Medical Research Council Muscle scale). Whole-exome sequencing was completed in three stages for those with skeletal muscle weakness. (3) Seven patients (17.9%) reported new-onset muscle weakness and motor limitations. Objective muscle weakness in the upper and lower extremities was seen in four patients. In three of them, exome sequencing revealed pathogenic/likely pathogenic variants in the genes encoding nexilin, myosin heavy chain, titin, and SPG7. (4) Our findings support a positive long-term outcome of skeletal muscle function in HTX patients. However, 10% of patients showed clinical signs of myopathy due to a possible genetic cause. The integration of genetic testing and standardized neurological assessment of motor function during the peri-HTX period should be considered.

A de novo ARIH2 gene mutation was detected in a patient with autism spectrum disorders and intellectual disability

E3 ubiquitin protein ligase encoded by ARIH2 gene catalyses the ubiquitination of target proteins and plays a crucial role in posttranslational modifications across various cellular processes. As prior documented, mutations in genes involved in the ubiquitination process are often associated with autism spectrum disorder (ASD) and/or intellectual disability (ID). In the current study, a de novo heterozygous mutation was identified in the splicing intronic region adjacent to the last exon of the ARIH2 gene using whole exome sequencing (WES). We hypothesize that this mutation, found in an ASD/ID patient, disrupts the protein Ariadne domain which is involved in the autoinhibition of ARIH2 enzyme. Predictive analyses elucidated the implications of the novel mutation in the splicing process and confirmed its autosomal dominant inheritance model. Nevertheless, we cannot exclude the possibility that other genetic factors, undetectable by WES, such as mutations in non-coding regions and polygenic risk in inter-allelic complementation, may contribute to the patient's phenotype. This work aims to suggest potential relationship between the detected mutation in ARIH2 gene and both ASD and ID, even though functional studies combined with new sequencing approaches will be necessary to validate this hypothesis.

Elucidating the Role of Human ALAS2 C-terminal Mutations Resulting in Loss of Function and Disease

The conserved enzyme aminolevulinic acid synthase (ALAS) initiates heme biosynthesis in certain bacteria and eukaryotes by catalyzing the condensation of glycine and succinyl-CoA to yield aminolevulinic acid. In humans, the ALAS isoform responsible for heme production during red blood cell development is the erythroid-specific ALAS2 isoform. Owing to its essential role in erythropoiesis, changes in human ALAS2 (hALAS2) function can lead to two different blood disorders. X-linked sideroblastic anemia results from loss of ALAS2 function, while X-linked protoporphyria results from gain of ALAS2 function. Interestingly, mutations in the ALAS2 C-terminal extension can be implicated in both diseases. Here, we investigate the molecular basis for enzyme dysfunction mediated by two previously reported C-terminal loss-of-function variants, hALAS2 V562A and M567I. We show that the mutations do not result in gross structural perturbations, but the enzyme stability for V562A is decreased. Additionally, we show that enzyme stability moderately increases with the addition of the pyridoxal 5'-phosphate (PLP) cofactor for both variants. The variants display differential binding to PLP and the individual substrates compared to wild-type hALAS2. Although hALAS2 V562A is a more active enzyme in vitro, it is less efficient concerning succinyl-CoA binding. In contrast, the M567I mutation significantly alters the cooperativity of substrate binding. In combination with previously reported cell-based studies, our work reveals the molecular basis by which hALAS2 C-terminal mutations negatively affect ALA production necessary for proper heme biosynthesis.

Evaluating predictors of kinase activity of STK11 variants identified in primary human non-small cell lung cancers

Critical evaluation of computational tools for predicting variant effects is important considering their increased use in disease diagnosis and driving molecular discoveries. In the sixth edition of the Critical Assessment of Genome Interpretation (CAGI) challenge, a dataset of 28 STK11 rare variants (27 missense, 1 single amino acid deletion), identified in primary non-small cell lung cancer biopsies, was experimentally assayed to characterize computational methods from four participating teams and five publicly available tools. Predictors demonstrated a high level of performance on key evaluation metrics, measuring correlation with the assay outputs and separating loss-of-function (LoF) variants from wildtype-like (WT-like) variants. The best participant model, 3Cnet, performed competitively with well-known tools. Unique to this challenge was that the functional data was generated with both biological and technical replicates, thus allowing the assessors to realistically establish maximum predictive performance based on experimental variability. Three out of the five publicly available tools and 3Cnet approached the performance of the assay replicates in separating LoF variants from WT-like variants. Surprisingly, REVEL, an often-used model, achieved a comparable correlation with the real-valued assay output as that seen for the experimental replicates. Performing variant interpretation by combining the new functional evidence with computational and population data evidence led to 16 new variants receiving a clinically actionable classification of likely pathogenic (LP) or likely benign (LB). Overall, the STK11 challenge highlights the utility of variant effect predictors in biomedical sciences and provides encouraging results for driving research in the field of computational genome interpretation.

Molecular Heterogeneity of Osteopetrosis in India: Report of 17 Novel Variants

Osteopetrosis is a clinically and genetically heterogeneous group of inherited bone disorders that is caused by defects in osteoclast formation or function. Treatment options vary with the disease severity and an accurate molecular diagnosis helps in prognostication and treatment decisions. We investigated the genetic causes of osteopetrosis in 31 unrelated patients of Indian origin. Screening for the genetic variants was done by Sanger sequencing or next generation sequencing in 48 samples that included 31 samples from index patients, 16 from parents' and 1 chorionic villus sample. A total of 30 variants, including 29 unique variants, were identified in 26 of the 31 patients in the study. TCIRG1 was the most involved gene (n = 14) followed by TNFRSF11A (n = 4) and CLCN7 (n = 3). A total of 17 novel variants were identified. Prenatal diagnosis was done in one family and the foetus showed homozygous c.807 + 2T > G variant in TCIRG1. Molecular diagnosis of osteopetrosis aids in therapeutic decisions including the need for a stem cell transplantation and gives a possible option of performing prenatal diagnosis in affected families. Further studies would help in understanding the genetic etiology in patients where no variants were identified.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12288-023-01732-4.

Implementation, Evolution, and Laboratory Performance of Methods-Based Proficiency Testing for Next-Generation Sequencing Detection of Germline Sequence Variants

CONTEXT.—

Next-generation sequencing (NGS)-based assays are used for diagnosis of diverse inherited disorders. Limited data are available pertaining to interlaboratory analytical performance of these assays.

OBJECTIVE.—

To report on the College of American Pathologists (CAP) NGS Germline Program, which is methods based, and explore the evolution in laboratory testing practices.

DESIGN.—

Results from the NGS Germline Program from 2016-2020 were analyzed for interlaboratory analytical performance. Self-reported laboratory testing practices were also evaluated.

RESULTS.—

From 2016-2020, a total of 297 laboratories participated in at least 1 program mailing. Of the 289 laboratories that provided information on tests offered, 138 (47.8%) offered only panel testing throughout their enrollment, while 35 (12.1%) offered panels and exome testing, 30 (10.4%) offered only exomes, 9 (3.1%) offered only genomes, and 15 (5.2%) offered panels, exomes, and genomes. The remainder (62 laboratories, 21.4%) changed their test offerings during the 2016-2020 timeframe. Considering each genomic position/interval, the median detection percentage at variant positions across the 2016-2020 mailings ranged from 94.3% to 100%, while at reference positions (no variant detected), the median correct response percentage was 100% across all mailings. When considering performance of individual laboratories, 89.5% (136 of 152) to 98.0% (149 of 152) of laboratories successfully met the detection threshold (≥90% of the variants present), while 94.6% (87 of 92) to 100% (163 of 163) of laboratories met the 95% specificity threshold across mailings.

CONCLUSIONS.—

Since the inception of this program, laboratories have consistently performed well. The median sensitivity and specificity of detection of sequence variants included in this program (eg, single nucleotide variants, insertions, and deletions) were 100.0%.

Improving transparency of computational tools for variant effect prediction

Efforts to integrate computational tools for variant effect prediction into the process of clinical decision-making are in progress. However, for such efforts to succeed and help to provide more informed clinical decisions, it is necessary to enhance transparency and address the current limitations of computational predictors.

Uncovering genetic mimics in multiple sclerosis: A single-center clinical exome sequencing study

Background

Multiple sclerosis (MS) shares clinical/radiological features with several monogenic diseases that can mimic MS.

Objective

We aimed to determine if exome sequencing can identify monogenic diseases in patients diagnosed with MS according to the McDonald criteria thus uncovering them as being misdiagnosed.

Methods

We performed whole exome sequencing in a cohort of 278 patients with MS, clinically or radiologically isolated syndrome without cerebrospinal fluid-specific oligoclonal bands (CSF-OCBs) (n = 228), a positive family history of MS (n = 44), or both (n = 6), thereby focusing on individuals potentially more likely to have underlying monogenic conditions mimicking MS. We prioritized 495 genes associated with monogenic diseases sharing features with MS.

Results

A disease-causing variant in NOTCH3 was identified in one patient without CSF-OCBs, no spinal lesions, with non-response to immunotherapy, and a family history of dementia, thereby converting the diagnosis to cerebral autosomal-dominant arteriopathy with subcortical infarcts and leukoencephalopathy (CADASIL). Moreover, 18 patients (6.5% of total) carried variants of unclear significance.

Conclusion

Monogenic diseases being misdiagnosed as MS seem rare in patients diagnosed with MS according to the McDonald criteria, even in CSF-OCB negative cases. The detected pathogenic NOTCH3 variant emphasizes CADASIL as a rare differential diagnosis and highlights the relevance of genetic testing in selected MS cases with atypical presentations.

Variant Effect Prediction in the Age of Machine Learning

Over the years, many computational methods have been created for the analysis of the impact of single amino acid substitutions resulting from single-nucleotide variants in genome coding regions. Historically, all methods have been supervised and thus limited by the inadequate sizes of experimentally curated data sets and by the lack of a standardized definition of variant effect. The emergence of unsupervised, deep learning (DL)-based methods raised an important question: Can machines learn the language of life from the unannotated protein sequence data well enough to identify significant errors in the protein "sentences"? Our analysis suggests that some unsupervised methods perform as well or better than existing supervised methods. Unsupervised methods are also faster and can, thus, be useful in large-scale variant evaluations. For all other methods, however, their performance varies by both evaluation metrics and by the type of variant effect being predicted. We also note that the evaluation of method performance is still lacking on less-studied, nonhuman proteins where unsupervised methods hold the most promise.