VarMap: a web tool for mapping genomic coordinates to protein sequence and structure and retrieving protein structural annotations

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.

Predicting human and viral protein variants affecting COVID-19 susceptibility and repurposing therapeutics

The COVID-19 disease is an ongoing global health concern. Although vaccination provides some protection, people are still susceptible to re-infection. Ostensibly, certain populations or clinical groups may be more vulnerable. Factors causing these differences are unclear and whilst socioeconomic and cultural differences are likely to be important, human genetic factors could influence susceptibility. Experimental studies indicate SARS-CoV-2 uses innate immune suppression as a strategy to speed-up entry and replication into the host cell. Therefore, it is necessary to understand the impact of variants in immunity-associated human proteins on susceptibility to COVID-19. In this work, we analysed missense coding variants in several SARS-CoV-2 proteins and their human protein interactors that could enhance binding affinity to SARS-CoV-2. We curated a dataset of 19 SARS-CoV-2: human protein 3D-complexes, from the experimentally determined structures in the Protein Data Bank and models built using AlphaFold2-multimer, and analysed the impact of missense variants occurring in the protein-protein interface region. We analysed 468 missense variants from human proteins and 212 variants from SARS-CoV-2 proteins and computationally predicted their impacts on binding affinities for the human viral protein complexes. We predicted a total of 26 affinity-enhancing variants from 13 human proteins implicated in increased binding affinity to SARS-CoV-2. These include key-immunity associated genes (TOMM70, ISG15, IFIH1, IFIT2, RPS3, PALS1, NUP98, AXL, ARF6, TRIMM, TRIM25) as well as important spike receptors (KREMEN1, AXL and ACE2). We report both common (e.g., Y13N in IFIH1) and rare variants in these proteins and discuss their likely structural and functional impact, using information on known and predicted functional sites. Potential mechanisms associated with immune suppression implicated by these variants are discussed. Occurrence of certain predicted affinity-enhancing variants should be monitored as they could lead to increased susceptibility and reduced immune response to SARS-CoV-2 infection in individuals/populations carrying them. Our analyses aid in understanding the potential impact of genetic variation in immunity-associated proteins on COVID-19 susceptibility and help guide drug-repurposing strategies.

Lethal Complications and Complex Genotypes in Shwachman Diamond Syndrome: Report of a Family with Recurrent Neonatal Deaths and a Case-Based Brief Review of the Literature

Shwachman Diamond Syndrome (SDS) is a multi-system disease characterized by exocrine pancreatic insufficiency with malabsorption, infantile neutropenia and aplastic anemia. Life-threatening complications include progression to acute myeloid leukemia (AML) or myelodysplastic syndrome (MDS), critical deep-tissue infections and asphyxiating thoracic dystrophy. In most patients, SDS results from biallelic pathogenic variants in the SBDS gene, different combinations of which contribute to heterogenous clinical presentations. Null variants are not well tolerated, supporting the theory that the loss of SBDS expression is likely lethal in both mice and humans. A novel complex genotype (SBDS:c.[242C>G;258+2T>C];[460-1G>A]/WFS1:c.[2327A>T];[1371G>T]) was detected in a family with recurrent neonatal deaths. A female neonate died three hours after birth with hemolytic anemia, and a male neonate with severe anemia, thrombocytopenia and neutropenia succumbed on day 40 after Staphylococcus epidermidis infection. A subsequent review of the literature focused on fatal complications, complex SBDS genotypes and/or unusual clinical presentations and disclosed rare cases, of which some had unexpected combinations of genetic and clinical findings. The impact of pathogenic variants and associated phenotypes is discussed in the context of data sharing towards expanding scientific expert networks, consolidating knowledge and advancing an understanding of novel underlying genotypes and complex phenotypes, facilitating informed clinical decisions and disease management.

SCN1A Channels a Wide Range of Epileptic Phenotypes: Report of Novel and Known Variants with Variable Presentations

SCN1A, the gene encoding for the Nav1.1 channel, exhibits dominant interneuron-specific expression, whereby variants disrupting the channel's function affect the initiation and propagation of action potentials and neuronal excitability causing various types of epilepsy. Dravet syndrome (DS), the first described clinical presentation of SCN1A channelopathy, is characterized by severe myoclonic epilepsy in infancy (SMEI). Variants' characteristics and other genetic or epigenetic factors lead to extreme clinical heterogeneity, ranging from non-epileptic conditions to developmental and epileptic encephalopathy (DEE). This current study reports on findings from 343 patients referred by physicians in hospitals and tertiary care centers in Greece between 2017 and 2023. Positive family history for specific neurologic disorders was disclosed in 89 cases and the one common clinical feature was the onset of seizures, at a mean age of 17 months (range from birth to 15 years old). Most patients were specifically referred for SCN1A investigation (Sanger Sequencing and MLPA) and only five for next generation sequencing. Twenty-six SCN1A variants were detected, including nine novel causative variants (c.4567A>Τ, c.5564C>A, c.2176+2T>C, c.3646G>C, c.4331C>A, c.1130_1131delGAinsAC, c.1574_1580delCTGAGGA, c.4620A>G and c.5462A>C), and are herein presented, along with subsequent genotype-phenotype associations. The identification of novel variants complements SCN1A databases extending our expertise on genetic counseling and patient and family management including gene-based personalized interventions.

Exploring the Micro-Mosaic Landscape of FGFR3 Mutations in the Ageing Male Germline and Their Potential Implications in Meiotic Differentiation

Advanced paternal age increases the risk of transmitting de novo germline mutations, particularly missense mutations activating the receptor tyrosine kinase (RTK) signalling pathway, as exemplified by the FGFR3 mutation, which is linked to achondroplasia (ACH). This risk is attributed to the expansion of spermatogonial stem cells carrying the mutation, forming sub-clonal clusters in the ageing testis, thereby increasing the frequency of mutant sperm and the number of affected offspring from older fathers. While prior studies proposed a correlation between sub-clonal cluster expansion in the testis and elevated mutant sperm production in older donors, limited data exist on the universality of this phenomenon. Our study addresses this gap by examining the testis-expansion patterns, as well as the increases in mutations in sperm for two FGFR3 variants-c.1138G>A (p.G380R) and c.1948A>G (p.K650E)-which are associated with ACH or thanatophoric dysplasia (TDII), respectively. Unlike the ACH mutation, which showed sub-clonal expansion events in an aged testis and a significant increase in mutant sperm with the donor's age, as also reported in other studies, the TDII mutation showed focal mutation pockets in the testis but exhibited reduced transmission into sperm and no significant age-related increase. The mechanism behind this divergence remains unclear, suggesting potential pleiotropic effects of aberrant RTK signalling in the male germline, possibly hindering differentiation requiring meiosis. This study provides further insights into the transmission risks of micro-mosaics associated with advanced paternal age in the male germline.

Neurodevelopmental disorders and cancer networks share pathways, but differ in mechanisms, signaling strength, and outcome

Epidemiological studies suggest that individuals with neurodevelopmental disorders (NDDs) are more prone to develop certain types of cancer. Notably, however, the case statistics can be impacted by late discovery of cancer in individuals afflicted with NDDs, such as intellectual disorders, autism, and schizophrenia, which may bias the numbers. As to NDD-associated mutations, in most cases, they are germline while cancer mutations are sporadic, emerging during life. However, somatic mosaicism can spur NDDs, and cancer-related mutations can be germline. NDDs and cancer share proteins, pathways, and mutations. Here we ask (i) exactly which features they share, and (ii) how, despite their commonalities, they differ in clinical outcomes. To tackle these questions, we employed a statistical framework followed by network analysis. Our thorough exploration of the mutations, reconstructed disease-specific networks, pathways, and transcriptome levels and profiles of autism spectrum disorder (ASD) and cancers, point to signaling strength as the key factor: strong signaling promotes cell proliferation in cancer, and weaker (moderate) signaling impacts differentiation in ASD. Thus, we suggest that signaling strength, not activating mutations, can decide clinical outcome.

Unraveling Vulnerabilities in Endocrine Therapy-Resistant HER2+/ER+ Breast Cancer

Breast tumors overexpressing human epidermal growth factor receptor (HER2) confer intrinsic resistance to endocrine therapy (ET), and patients with HER2/estrogen receptor-positive (HER2+/ER+) breast cancer (BCa) are less responsive to ET than HER2-/ER+. However, real-world evidence reveals that a large subset of patients with HER2+/ER+ receive ET as monotherapy, positioning this treatment pattern as a clinical challenge. In the present study, we developed and characterized 2 in vitro models of ET-resistant (ETR) HER2+/ER+ BCa to identify possible therapeutic vulnerabilities. To mimic ETR to aromatase inhibitors (AIs), we developed 2 long-term estrogen deprivation (LTED) cell lines from BT-474 (BT474) and MDA-MB-361 (MM361). Growth assays, PAM50 subtyping, and genomic and transcriptomic analyses, followed by validation and functional studies, were used to identify targetable differences between ET-responsive parental and ETR-LTED HER2+/ER+ cells. Compared to their parental cells, MM361 LTEDs grew faster, lost ER, and increased HER2 expression, whereas BT474 LTEDs grew slower and maintained ER and HER2 expression. Both LTED variants had reduced responsiveness to fulvestrant. Whole-genome sequencing of aggressive MM361 LTEDs identified mutations in genes encoding transcription factors and chromatin modifiers. Single-cell RNA sequencing demonstrated a shift towards non-luminal phenotypes, and revealed metabolic remodeling of MM361 LTEDs, with upregulated lipid metabolism and ferroptosis-associated antioxidant genes, including GPX4. Combining a GPX4 inhibitor with anti-HER2 agents induced significant cell death in both MM361 and BT474 LTEDs. The BT474 and MM361 AI-resistant models capture distinct phenotypes of HER2+/ER+ BCa and identify altered lipid metabolism and ferroptosis remodeling as vulnerabilities of this type of ETR BCa.

Unraveling Vulnerabilities in Endocrine Therapy-Resistant HER2+/ER+ Breast Cancer

Background

Breast tumors overexpressing human epidermal growth factor receptor (HER2) confer intrinsic resistance to endocrine therapy (ET), and patients with HER2/ estrogen receptor-positive (HER2+/HR+) breast cancer (BCa) are less responsive to ET than HER2-/ER+. However, real-world evidence reveals that a large subset of HER2+/ER+ patients receive ET as monotherapy, positioning this treatment pattern as a clinical challenge. In the present study, we developed and characterized two distinct in vitro models of ET-resistant (ETR) HER2+/ER+ BCa to identify possible therapeutic vulnerabilities.

Methods

To mimic ETR to aromatase inhibitors (AI), we developed two long-term estrogen-deprived (LTED) cell lines from BT-474 (BT474) and MDA-MB-361 (MM361). Growth assays, PAM50 molecular subtyping, genomic and transcriptomic analyses, followed by validation and functional studies, were used to identify targetable differences between ET-responsive parental and ETR-LTED HER2+/ER+ cells.

Results

Compared to their parental cells, MM361 LTEDs grew faster, lost ER, and increased HER2 expression, whereas BT474 LTEDs grew slower and maintained ER and HER2 expression. Both LTED variants had reduced responsiveness to fulvestrant. Whole-genome sequencing of the more aggressive MM361 LTED model system identified exonic mutations in genes encoding transcription factors and chromatin modifiers. Single-cell RNA sequencing demonstrated a shift towards non-luminal phenotypes, and revealed metabolic remodeling of MM361 LTEDs, with upregulated lipid metabolism and antioxidant genes associated with ferroptosis, including GPX4. Combining the GPX4 inhibitor RSL3 with anti-HER2 agents induced significant cell death in both the MM361 and BT474 LTEDs.

Conclusions

The BT474 and MM361 AI-resistant models capture distinct phenotypes of HER2+/ER+ BCa and identify altered lipid metabolism and ferroptosis remodeling as vulnerabilities of this type of ETR BCa.

The clinical and genetic spectrum of autosomal-recessive TOR1A-related disorders

In the field of rare diseases, progress in molecular diagnostics led to the recognition that variants linked to autosomal-dominant neurodegenerative diseases of later onset can, in the context of biallelic inheritance, cause devastating neurodevelopmental disorders and infantile or childhood-onset neurodegeneration. TOR1A-associated arthrogryposis multiplex congenita 5 (AMC5) is a rare neurodevelopmental disorder arising from biallelic variants in TOR1A, a gene that in the heterozygous state is associated with torsion dystonia-1 (DYT1 or DYT-TOR1A), an early-onset dystonia with reduced penetrance. While 15 individuals with AMC5-TOR1A have been reported (less than 10 in detail), a systematic investigation of the full disease-associated spectrum has not been conducted. Here, we assess the clinical, radiological and molecular characteristics of 57 individuals from 40 families with biallelic variants in TOR1A. Median age at last follow-up was 3 years (0-24 years). Most individuals presented with severe congenital flexion contractures (95%) and variable developmental delay (79%). Motor symptoms were reported in 79% and included lower limb spasticity and pyramidal signs, as well as gait disturbances. Facial dysmorphism was an integral part of the phenotype, with key features being a broad/full nasal tip, narrowing of the forehead and full cheeks. Analysis of disease-associated manifestations delineated a phenotypic spectrum ranging from normal cognition and mild gait disturbance to congenital arthrogryposis, global developmental delay, intellectual disability, absent speech and inability to walk. In a subset, the presentation was consistent with foetal akinesia deformation sequence with severe intrauterine abnormalities. Survival was 71%, with higher mortality in males. Death occurred at a median age of 1.2 months (1 week-9 years), due to respiratory failure, cardiac arrest or sepsis. Analysis of brain MRI studies identified non-specific neuroimaging features, including a hypoplastic corpus callosum (72%), foci of signal abnormality in the subcortical and periventricular white matter (55%), diffuse white matter volume loss (45%), mega cisterna magna (36%) and arachnoid cysts (27%). The molecular spectrum included 22 distinct variants, defining a mutational hotspot in the C-terminal domain of the Torsin-1A protein. Genotype-phenotype analysis revealed an association of missense variants in the 3-helix bundle domain to an attenuated phenotype, while missense variants near the Walker A/B motif as well as biallelic truncating variants were linked to early death. In summary, this systematic cross-sectional analysis of a large cohort of individuals with biallelic TOR1A variants across a wide age-range delineates the clinical and genetic spectrum of TOR1A-related autosomal-recessive disease and highlights potential predictors for disease severity and survival.

Rare variant contribution to cholestatic liver disease in a South Asian population in the United Kingdom

This study assessed the contribution of five genes previously known to be involved in cholestatic liver disease in British Bangladeshi and Pakistani people. Five genes (ABCB4, ABCB11, ATP8B1, NR1H4, TJP2) were interrogated by exome sequencing data of 5236 volunteers. Included were non-synonymous or loss of function (LoF) variants with a minor allele frequency < 5%. Variants were filtered, and annotated to perform rare variant burden analysis, protein structure, and modelling analysis in-silico. Out of 314 non-synonymous variants, 180 fulfilled the inclusion criteria and were mostly heterozygous unless specified. 90 were novel and of those variants, 22 were considered likely pathogenic and 9 pathogenic. We identified variants in volunteers with gallstone disease (n = 31), intrahepatic cholestasis of pregnancy (ICP, n = 16), cholangiocarcinoma and cirrhosis (n = 2). Fourteen novel LoF variants were identified: 7 frameshift, 5 introduction of premature stop codon and 2 splice acceptor variants. The rare variant burden was significantly increased in ABCB11. Protein modelling demonstrated variants that appeared to likely cause significant structural alterations. This study highlights the significant genetic burden contributing to cholestatic liver disease. Novel likely pathogenic and pathogenic variants were identified addressing the underrepresentation of diverse ancestry groups in genomic research.

The clinical and molecular spectrum of ZFYVE26-associated hereditary spastic paraplegia: SPG15

In the field of hereditary spastic paraplegia (HSP), progress in molecular diagnostics needs to be translated into robust phenotyping studies to understand genetic and phenotypic heterogeneity and to support interventional trials. ZFYVE26-associated hereditary spastic paraplegia (HSP-ZFYVE26, SPG15) is a rare, early-onset complex HSP, characterized by progressive spasticity and a variety of other neurological symptoms. While prior reports, often in populations with high rates of consanguinity, have established a general phenotype, there is a lack of systematic investigations and a limited understanding of age-dependent manifestation of symptoms. Here we delineate the clinical, neuroimaging and molecular features of 44 individuals from 36 families, the largest cohort assembled to date. Median age at last follow-up was 23.8 years covering a wide age range (11-61 years). While symptom onset often occurred in early childhood [median: 24 months, interquartile range (IQR) = 24], a molecular diagnosis was reached at a median age of 18.8 years (IQR = 8), indicating significant diagnostic delay. We demonstrate that most patients present with motor and/or speech delay or learning disabilities. Importantly, these developmental symptoms preceded the onset of motor symptoms by several years. Progressive spasticity in the lower extremities, the hallmark feature of HSP-ZFYVE26, typically presents in adolescence and involves the distal lower limbs before progressing proximally. Spasticity in the upper extremities was seen in 64%. We found a high prevalence of extrapyramidal movement disorders including cerebellar ataxia (64%) and dystonia (11%). Parkinsonism (16%) was present in a subset and showed no sustained response to levodopa. Cognitive decline and neurogenic bladder dysfunction progressed over time in most patients. A systematic analysis of brain MRI features revealed a common diagnostic signature consisting of thinning of the anterior corpus callosum, signal changes of the anterior forceps and non-specific cortical and cerebellar atrophy. The molecular spectrum included 45 distinct variants, distributed across the protein structure without mutational hotspots. Spastic Paraplegia Rating Scale scores, SPATAX Disability Scores and the Four Stage Functional Mobility Score showed moderate strength in representing the proportion of variation between disease duration and motor dysfunction. Plasma neurofilament light chain levels were significantly elevated in all patients (Mann-Whitney U-test, P < 0.0001) and were correlated inversely with age (Spearman's rank correlation coefficient r = -0.65, P = 0.01). In summary, our systematic cross-sectional analysis of HSP-ZFYVE26 patients across a wide age-range, delineates core clinical, neuroimaging and molecular features and identifies markers of disease severity. These results raise awareness to this rare disease, facilitate an early diagnosis and create clinical trial readiness.

Delineation of functionally essential protein regions for 242 neurodevelopmental genes

Neurodevelopmental disorders (NDDs), including severe paediatric epilepsy, autism and intellectual disabilities are heterogeneous conditions in which clinical genetic testing can often identify a pathogenic variant. For many of them, genetic therapies will be tested in this or the coming years in clinical trials. In contrast to first-generation symptomatic treatments, the new disease-modifying precision medicines require a genetic test-informed diagnosis before a patient can be enrolled in a clinical trial. However, even in 2022, most identified genetic variants in NDD genes are 'variants of uncertain significance'. To safely enrol patients in precision medicine clinical trials, it is important to increase our knowledge about which regions in NDD-associated proteins can 'tolerate' missense variants and which ones are 'essential' and will cause a NDD when mutated. In addition, knowledge about functionally indispensable regions in the 3D structure context of proteins can also provide insights into the molecular mechanisms of disease variants. We developed a novel consensus approach that overlays evolutionary, and population based genomic scores to identify 3D essential sites (Essential3D) on protein structures. After extensive benchmarking of AlphaFold predicted and experimentally solved protein structures, we generated the currently largest expert curated protein structure set for 242 NDDs and identified 14 377 Essential3D sites across 189 gene disorders associated proteins. We demonstrate that the consensus annotation of Essential3D sites improves prioritization of disease mutations over single annotations. The identified Essential3D sites were enriched for functional features such as intermembrane regions or active sites and discovered key inter-molecule interactions in protein complexes that were otherwise not annotated. Using the currently largest autism, developmental disorders, and epilepsies exome sequencing studies including >360 000 NDD patients and population controls, we found that missense variants at Essential3D sites are 8-fold enriched in patients. In summary, we developed a comprehensive protein structure set for 242 NDDs and identified 14 377 Essential3D sites in these. All data are available at https://es-ndd.broadinstitute.org for interactive visual inspection to enhance variant interpretation and development of mechanistic hypotheses for 242 NDDs genes. The provided resources will enhance clinical variant interpretation and in silico drug target development for NDD-associated genes and encoded proteins.

Effect of exogenous IL-37 on immune cells from a patient carrying a potential IL37 loss-of-function variant: A case study

INTRODUCTION

Chronic inflammatory or autoimmune diseases are commonly treated with immunosuppressive medication such as NSAIDs, corticosteroids, or antibodies against specific cytokines (TNF, IL-1 IL-17, IL-23, etc.) or signalling cascades (e.g. JAK-STAT inhibitors). Using sequencing data to locate genetic mutations in relevant genes allows the identification of alternative targets in a patient-tailored therapy setting. Interleukin (IL)-37 is an anti-inflammatory cytokine with broad effects on innate and adaptive immune cell function. Dysfunctional IL-37 expression or signalling is linked to various autoinflammatory disorders. The administration of recombinant IL-37 to hyperinflammatory patients that are non-responsive to standard treatment bears the potential to alleviate symptoms.

METHODS

In this case study, the (hyper)responsiveness of immune cell subsets was investigated in a single patient with a seronegative autoimmune disorder who carries a heterozygous stop-gain variant in IL37 (IL37 Chr2(GRCh37):g.113670640G > A NM_014439.3:c.51G > A p.(Trp17*)). As the patient has been non-responsive to blockage of TNF or IL-1 by Etanercept or Anakinra, respectively, additional in-vitro experiments were set out to elucidate whether treatment with recombinant IL-37 could normalise observed immune cell functions.

FINDINGS

Characterisation of immune cell function showed no elevated overall production of acute-phase pro-inflammatory cytokines by patient PBMCs and neutrophils at baseline or upon stimulation. T-cell responses were elevated, as was the metabolic activity and IL-1Ra production of PBMCs at baseline. The identified stop-gain variant in IL37 does not result in the absence of the protein in circulation. In line with this, treatment with recombinant IL-37 did overall not dampen immune responses with the exception of the complete suppression of IL-17.

CONCLUSION

The heterozygous stop-gain variant in IL37 (IL37 NM_014439.3:c.51G > A p.(Trp17*)) is not of functional relevance as we observed no clear pro-inflammatory phenotype in immune cells of a patient carrying this variant.

Genetic and phenotypic continuum of HOXA genes: A case with double HOXA9/HOXA13 mutations

The HOXA genes cluster plays a key role in embryologic development. Mutations in HOXA genes have been linked to different human phenotypes, including developmental delay, limb anomalies, and urogenital malformations. The present study reported a clinical and genetic investigation of a female patient with polymalformative syndrome including left arm agenesis, bicornuate uterus and bicuspid aortic valve. Using whole exome sequencing, two heterozygous missense variants were identified. Of these, one was a novel variant in the HOXA13 gene [p.(Tyr290Ser)] and the second a heterozygous variant in the HOXA9 gene [p.(Ala102Pro)]. To the best of our knowledge, this is the first association of HOXA9/HOXA13 point mutations linked to a syndromic case. In conclusion, the present study suggested that the phenotypic spectrum of vertebral anomalies, anal atresia, cardiac defects, tracheo‑esophageal fistula, renal anomalies and limb abnormalities/hand‑foot‑genital syndrome may be attributable to the combination of different HOXA variants, particularly in patients with a severe clinical presentation. The current report contributed as well to the molecular understanding of HOXA genes‑related phenotypes via the identification of novel variant and genes associations.

Muscle spasms as presenting feature of Nivelon-Nivelon-Mabile syndrome

Hedgehog acyltransferase gene (HHAT)-associated Nivelon-Nivelon-Mabile syndrome (NNMS) is a rare genetic disorder of multiple system involvement with microcephaly, central nervous system malformations, skeletal dysplasia, and 46,XY sex reversal. Other variable and inconsistent features reported in this disorder are muscle spasms, facial dysmorphism, prenatal onset growth restriction, microphthalmia, and holoprosencephaly. This is the sixth postnatal reported patient with biallelic variants in HHAT gene, who presented with microcephaly, short stature, muscle hypertrophy, muscle spasms, and facial dysmorphism. The most prominent and presenting finding in this patient were muscle hypertrophy and muscle spasms which had a clinical response to phenytoin and acetazolamide treatment. Our report emphasizes the phenotypic variability of NNMS and further reiterates muscle spasms as an important clinical manifestation of this extremely rare condition.

A novel pathogenic ATP6V1B2 variant: Widening the genotypic spectrum of the epileptic neurodevelopmental phenotype

ATP6V1B2 pathogenic variants are linked with variable phenotypes, such as dominant deafness-onychodystrophy syndrome (DDOD), autosomal dominant Zimmermann-Laband syndrome type 2 (ZLS2), and some cases of DOORS (deafness, onychodystrophy, osteodystrophy, intellectual disability [ID], and seizures). Epilepsy was first linked to ATP6V1B2, when the p.(Glu374Gln) missense variant was detected in a patient with ID and seizures, but without characteristic features of DDOD or ZLS2 syndromes. We herein report a novel pathogenic ATP6V1B2:p.Glu374Gly variant detected in an adult patient with ID and myoclonic-atonic seizures. The (re)occurrence of different variants affecting the same highly conserved hydrophilic glutamic acid on position 374 of the V-proton ATPase subunit B, indicates a potential novel pathogenic hotspot and a critical role for the specific residue in the development of epilepsy. ATP6V1B2 gene defects should be considered when analyzing patients with epilepsy, even in the absence of most cardinal features of DDOD, DOORS, or ZLS such as deafness, onychodystrophy, and osteodystrophy.

VIVID: a web application for variant interpretation and visualisation in multidimensional analyses

Large-scale comparative genomics- and population genetic studies generate enormous amounts of polymorphism data in the form of DNA variants. Ultimately, the goal of many of these studies is to associate genetic variants to phenotypes or fitness. We introduce VIVID, an interactive, user-friendly web application that integrates a wide range of approaches for encoding genotypic to phenotypic information in any organism or disease, from an individual or population, in three-dimensional (3D) space. It allows mutation mapping and annotation, calculation of interactions and conservation scores, prediction of harmful effects, analysis of diversity and selection, and 3D visualization of genotypic information encoded in Variant Call Format on AlphaFold2 protein models. VIVID enables the rapid assessment of genes of interest in the study of adaptive evolution and the genetic load, and it helps prioritizing targets for experimental validation. We demonstrate the utility of VIVID by exploring the evolutionary genetics of the parasitic protist Plasmodium falciparum, revealing geographic variation in the signature of balancing selection in potential targets of functional antibodies.

Genetic association and Mendelian randomization for hypothyroidism highlight immune molecular mechanisms

We carried out a genome-wide association analysis including 51,194 cases of hypothyroidism and 443,383 controls. In total, 139 risk loci were associated to hypothyroidism with genes involved in lymphocyte function. Candidate genes associated with hypothyroidism were identified by using molecular quantitative trait loci, colocalization, and enhancer-promoter chromatin looping. Mendelian randomization (MR) identified 42 blood expressed genes and circulating proteins as candidate causal molecules in hypothyroidism. Drug-gene interaction analysis provided evidence that immune checkpoint and tyrosine kinase inhibitors used in cancer therapy increase the risk of hypothyroidism. Hence, integrative mapping and MR support that expression of genes and proteins enriched in lymphocyte function are associated with the risk of hypothyroidism and provide genetic evidence for drug-induced hypothyroidism and identify actionable potential drug targets.

•

GWAS for hypothyroidism identified 139 risk loci including 76 novel associations

•

GWAS was enriched in pathways related to lymphocyte function

•

In total, 28 potentially deleterious missense variants were identified

•

Mendelian randomization and colocalization identified 61 blood causal candidate genes

Functional validation of novel variants in B4GALNT1 associated with early-onset complex hereditary spastic paraplegia with impaired ganglioside synthesis

Childhood-onset forms of hereditary spastic paraplegia are ultra-rare diseases and often present with complex features. Next-generation-sequencing allows for an accurate diagnosis in many cases but the interpretation of novel variants remains challenging, particularly for missense mutations. Where sufficient knowledge of the protein function and/or downstream pathways exists, functional studies in patient-derived cells can aid the interpretation of molecular findings. We here illustrate the case of a 13-year-old female who presented with global developmental delay and later mild intellectual disability, progressive spastic diplegia, spastic-ataxic gait, dysarthria, urinary urgency, and loss of deep tendon reflexes of the lower extremities. Exome sequencing showed a novel splice-site variant in trans with a novel missense variant in B4GALNT1 [NM_001478.5: c.532-1G>C/c.1556G>C (p.Arg519Pro)]. Functional studies in patient-derived fibroblasts and cell models of GM2 synthase deficiency confirmed a loss of B4GALNT1 function with no synthesis of GM2 and other downstream gangliosides. Collectively these results established the diagnosis of B4GALNT1-associated HSP (SPG26). Our approach illustrates the importance of careful phenotyping and functional characterization of novel gene variants, particularly in the setting of ultra-rare diseases, and expands the clinical and molecular spectrum of SPG26, a disorder of complex ganglioside biosynthesis.

Assessing the clinical utility of protein structural analysis in genomic variant classification: experiences from a diagnostic laboratory

BACKGROUND

The widespread clinical application of genome-wide sequencing has resulted in many new diagnoses for rare genetic conditions, but testing regularly identifies variants of uncertain significance (VUS). The remarkable rise in the amount of genomic data has been paralleled by a rise in the number of protein structures that are now publicly available, which may have clinical utility for the interpretation of missense and in-frame insertions or deletions.

METHODS

Within a UK National Health Service genomic medicine diagnostic laboratory, we investigated the number of VUS over a 5-year period that were evaluated using protein structural analysis and how often this analysis aided variant classification.

RESULTS

We found 99 novel missense and in-frame variants across 67 genes that were initially classified as VUS by our diagnostic laboratory using standard variant classification guidelines and for which further analysis of protein structure was requested. Evidence from protein structural analysis was used in the re-assessment of 64 variants, of which 47 were subsequently reclassified as pathogenic or likely pathogenic and 17 remained as VUS. We identified several case studies where protein structural analysis aided variant interpretation by predicting disease mechanisms that were consistent with the observed phenotypes, including loss-of-function through thermodynamic destabilisation or disruption of ligand binding, and gain-of-function through de-repression or escape from proteasomal degradation.

CONCLUSIONS

We have shown that using in silico protein structural analysis can aid classification of VUS and give insights into the mechanisms of pathogenicity. Based on our experience, we propose a generic evidence-based workflow for incorporating protein structural information into diagnostic practice to facilitate variant classification.

Analysis of Novel Variants Associated with Three Human Ovarian Cancer Cell Lines

Background:

Identification of mutations is of great significance in cancer research, as it can contribute to development of therapeutic strategies and prevention of cancer formation. Ovarian cancer is one of the leading cancer-related causes of death in Taiwan. Accumulation of genetic mutations can lead to cancer.

Objective:

We utilized whole-exome sequencing to explore cancer-associated missense variants in three human ovarian cancer cell lines derived from Taiwanese patients.

Methods:

We use (i) cell line whole-exome sequencing data, (ii) 188 patients’ whole-exome sequencing data, and (iii) use of in vitro experiments to verify predicted variant results. We establish an effective analysis workflow for discovery of novel ovarian cancer variants, comprising three steps: (i) use of public databases and in-house hospital data to select novel variants (ii) investigation of protein structural stability caused by genetic mutations, and (iii) use of in vitro experiments to verify predictions.

Results:

Our study enumerated 296 novel variants by imposing specific criteria and using sophisticated bioinformatics tools for further analysis. Eleven and 54 missense novel variants associated with cancerous and non-cancerous genes, respectively, were identified. We show that 13 missense mutations affect the stability of protein 3D structure, while 11 disease-causing novel variants were confirmed by PCR sequencing. Among these, ten variants were predicted to be pathogenic, while the pathogenicity of one was uncertain.

Conclusion:

We confirm that novel variant genes play a crucial role in ovarian cancer patients, with 11 novel variants that may promote progression and development of ovarian cancer.

A Novel Osteochondrodysplasia With Empty Sella Associates With a TBX2 Variant

Skeletal dysplasias comprise a heterogenous group of developmental disorders of skeletal and cartilaginous tissues. Several different forms have been described and the full spectrum of their clinical manifestations and underlying genetic causes are still incompletely understood. We report a three-generation Finnish family with an unusual, autosomal dominant form of osteochondrodysplasia and an empty sella. Affected individuals (age range 24-44 years) exhibit unusual codfish-shaped vertebrae, severe early-onset and debilitating osteoarthritis and an empty sella without endocrine abnormalities. Clinical characteristics also include mild dysmorphic features, reduced sitting height ratio, and obesity. Whole-exome sequencing excluded known skeletal dysplasias and identified a novel heterozygous missense mutation c.899C>T (p.Thr300Met) in TBX2, confirmed by Sanger sequencing. TBX2 is important for development of the skeleton and the brain and three prior reports have described variations in TBX2 in patients portraying a complex phenotype with vertebral anomalies, craniofacial dysmorphism and endocrine dysfunctions. Our mutation lies near a previously reported disease-causing variant and is predicted pathogenic with deleterious effects on protein function. Our findings expand the current spectrum of skeletal dysplasias, support the association of TBX2 mutations with skeletal dysplasia and suggest a role for TBX2 in development of the spinal and craniofacial structures and the pituitary gland.

Novel de novo pathogenic variant in the GNAI1 gene as a cause of severe disorders of intellectual development

Pathogenic sequence variant in the GNAI1 gene were recently introduced as a cause of novel syndrome with a manifestation of variable developmental delay and autistic features. In our study, we report a case of monozygotic twins with severe intellectual disability and motor delay and developmental dysphasia. Both probands and their parents were examined using multi-step molecular diagnostic algorithm including whole-exome sequencing (WES), resulting in the identification of a novel, de novo pathogenic sequence variant in the GNAI1 gene, NM_002069.6:c.815 A>G, p.(Asp272Gly) in probands. Using WES we also verified the microarray findings of a familial 8q24.23q24.3 duplication and heterozygous 5q13.2 deletion, not associated with clinical symptoms in probands. Our results confirmed the role of the GNAI1 gene in the pathogenesis of syndromic neurodevelopmental disorders. They support trio- or quatro-based WES as a suitable molecular diagnostics method for the simultaneous detection of clinically relevant sequence variants and CNVs in individuals with neurodevelopmental disorders and rare diseases.

Variant PNLDC1, Defective piRNA Processing, and Azoospermia

BACKGROUND

P-element-induced wimpy testis (PIWI)-interacting RNAs (piRNAs) are short (21 to 35 nucleotides in length) and noncoding and are found almost exclusively in germ cells, where they regulate aberrant expression of transposable elements and postmeiotic gene expression. Critical to the processing of piRNAs is the protein poly(A)-specific RNase-like domain containing 1 (PNLDC1), which trims their 3' ends and, when disrupted in mice, causes azoospermia and male infertility.

METHODS

We performed exome sequencing on DNA samples from 924 men who had received a diagnosis of nonobstructive azoospermia. Testicular-biopsy samples were analyzed by means of histologic and immunohistochemical tests, in situ hybridization, reverse-transcriptase-quantitative-polymerase-chain-reaction assay, and small-RNA sequencing.

RESULTS

Four unrelated men of Middle Eastern descent who had nonobstructive azoospermia were found to carry mutations in PNLDC1: the first patient had a biallelic stop-gain mutation, p.R452Ter (rs200629089; minor allele frequency, 0.00004); the second, a novel biallelic missense variant, p.P84S; the third, two compound heterozygous mutations consisting of p.M259T (rs141903829; minor allele frequency, 0.0007) and p.L35PfsTer3 (rs754159168; minor allele frequency, 0.00004); and the fourth, a novel biallelic canonical splice acceptor site variant, c.607-2A→T. Testicular histologic findings consistently showed error-prone meiosis and spermatogenic arrest with round spermatids of type Sa as the most advanced population of germ cells. Gene and protein expression of PNLDC1, as well as the piRNA-processing proteins PIWIL1, PIWIL4, MYBL1, and TDRKH, were greatly diminished in cells of the testes. Furthermore, the length distribution of piRNAs and the number of pachytene piRNAs was significantly altered in men carrying PNLDC1 mutations.

CONCLUSIONS

Our results suggest a direct mechanistic effect of faulty piRNA processing on meiosis and spermatogenesis in men, ultimately leading to male infertility. (Funded by Innovation Fund Denmark and others.).

Identification of Novel Genomic-Variant Patterns of OR56A5, OR52L1, and CTSD in Retinitis Pigmentosa Patients by Whole-Exome Sequencing

Inherited retinal dystrophies (IRDs) are rare but highly heterogeneous genetic disorders that affect individuals and families worldwide. However, given its wide variability, its analysis of the driver genes for over 50% of the cases remains unexplored. The present study aims to identify novel driver genes, disease-causing variants, and retinitis pigmentosa (RP)-associated pathways. Using family-based whole-exome sequencing (WES) to identify putative RP-causing rare variants, we identified a total of five potentially pathogenic variants located in genes OR56A5, OR52L1, CTSD, PRF1, KBTBD13, and ATP2B4. Of the variants present in all affected individuals, genes OR56A5, OR52L1, CTSD, KBTBD13, and ATP2B4 present as missense mutations, while PRF1 and CTSD present as frameshift variants. Sanger sequencing confirmed the presence of the novel pathogenic variant PRF1 (c.124_128del) that has not been reported previously. More causal-effect or evidence-based studies will be required to elucidate the precise roles of these SNPs in the RP pathogenesis. Taken together, our findings may allow us to explore the risk variants based on the sequencing data and upgrade the existing variant annotation database in Taiwan. It may help detect specific eye diseases such as retinitis pigmentosa in East Asia.

Common genetic associations between age-related diseases

Age is a common risk factor in many diseases, but the molecular basis for this relationship is elusive. In this study we identified 4 disease clusters from 116 diseases in the UK Biobank data, defined by their age-of-onset profiles, and found that diseases with the same onset profile are genetically more similar, suggesting a common etiology. This similarity was not explained by disease categories, co-occurrences or disease cause-effect relationships. Two of the four disease clusters had an increased risk of occurrence from age 20 and 40 years respectively. They both showed an association with known aging-related genes, yet differed in functional enrichment and evolutionary profiles. Moreover, they both had age-related expression and methylation changes. We also tested mutation accumulation and antagonistic pleiotropy theories of aging and found support for both.

ADDRESS: A Database of Disease-associated Human Variants Incorporating Protein Structure and Folding Stabilities

Numerous human diseases are caused by mutations in genomic sequences. Since amino acid changes affect protein function through mechanisms often predictable from protein structure, the integration of structural and sequence data enables us to estimate with greater accuracy whether and how a given mutation will lead to disease. Publicly available annotated databases enable hypothesis assessment and benchmarking of prediction tools. However, the results are often presented as summary statistics or black box predictors, without providing full descriptive information. We developed a new semi-manually curated human variant database presenting information on the protein contact-map, sequence-to-structure mapping, amino acid identity change, and stability prediction for the popular UniProt database. We found that the profiles of pathogenic and benign missense polymorphisms can be effectively deduced using decision trees and comparative analyses based on the presented dataset. The database is made publicly available through https://zhanglab.ccmb.med.umich.edu/ADDRESS.

From chemoproteomic-detected amino acids to genomic coordinates: insights into precise multi-omic data integration

The integration of proteomic, transcriptomic, and genetic variant annotation data will improve our understanding of genotype-phenotype associations. Due, in part, to challenges associated with accurate inter-database mapping, such multi-omic studies have not extended to chemoproteomics, a method that measures the intrinsic reactivity and potential "druggability" of nucleophilic amino acid side chains. Here, we evaluated mapping approaches to match chemoproteomic-detected cysteine and lysine residues with their genetic coordinates. Our analysis revealed that database update cycles and reliance on stable identifiers can lead to pervasive misidentification of labeled residues. Enabled by this examination of mapping strategies, we then integrated our chemoproteomics data with computational methods for predicting genetic variant pathogenicity, which revealed that codons of highly reactive cysteines are enriched for genetic variants that are predicted to be more deleterious and allowed us to identify and functionally characterize a new damaging residue in the cysteine protease caspase-8. Our study provides a roadmap for more precise inter-database mapping and points to untapped opportunities to improve the predictive power of pathogenicity scores and to advance prioritization of putative druggable sites.

MISCAST: MIssense variant to protein StruCture Analysis web SuiTe

Human genome sequencing efforts have greatly expanded, and a plethora of missense variants identified both in patients and in the general population is now publicly accessible. Interpretation of the molecular-level effect of missense variants, however, remains challenging and requires a particular investigation of amino acid substitutions in the context of protein structure and function. Answers to questions like 'Is a variant perturbing a site involved in key macromolecular interactions and/or cellular signaling?', or 'Is a variant changing an amino acid located at the protein core or part of a cluster of known pathogenic mutations in 3D?' are crucial. Motivated by these needs, we developed MISCAST (missense variant to protein structure analysis web suite; http://miscast.broadinstitute.org/). MISCAST is an interactive and user-friendly web server to visualize and analyze missense variants in protein sequence and structure space. Additionally, a comprehensive set of protein structural and functional features have been aggregated in MISCAST from multiple databases, and displayed on structures alongside the variants to provide users with the biological context of the variant location in an integrated platform. We further made the annotated data and protein structures readily downloadable from MISCAST to foster advanced offline analysis of missense variants by a wide biological community.

Splicing impact of deep exonic missense variants in CAPN3 explored systematically by minigene functional assay

Improving the accuracy of variant interpretation during diagnostic sequencing is a major goal for genomic medicine. To explore an often-overlooked splicing effect of missense variants, we developed the functional assay ("minigene") for the majority of exons of CAPN3, the gene responsible for limb girdle muscular dystrophy. By systematically screening 21 missense variants distributed along the gene, we found that eight clinically relevant missense variants located at a certain distance from the exon-intron borders (deep exonic missense variants) disrupted normal splicing of CAPN3 exons. Several recent machine learning-based computational tools failed to predict splicing impact for the majority of these deep exonic missense variants, highlighting the importance of including variants of this type in the training sets during the future algorithm development. Overall, 24 variants in CAPN3 gene were explored, leading to the change in the American College of Medical Genetics and Genomics classification of seven of them when results of the "minigene" functional assay were considered. Our findings reveal previously unknown splicing impact of several clinically important variants in CAPN3 and draw attention to the existence of deep exonic variants with a disruptive effect on gene splicing that could be overlooked by the current approaches in clinical genetics.

Structural analysis of pathogenic missense mutations in GABRA2 and identification of a novel de novo variant in the desensitization gate

Background

Cys‐loop receptors control neuronal excitability in the brain and their dysfunction results in numerous neurological disorders. Recently, six missense variants in GABRA2, a member of this family, have been associated with early infantile epileptic encephalopathy (EIEE). We identified a novel de novo missense variant in GABRA2 in a patient with EIEE and performed protein structural analysis of the seven variants.

Methods

The novel variant was identified by trio whole‐genome sequencing. We performed protein structural analysis of the seven variants, and compared them to previously reported pathogenic mutations at equivalent positions in other Cys‐loop receptors. Additionally, we studied the distribution of disease‐associated variants in the transmembrane helices of these proteins.

Results

The seven variants are in the transmembrane domain, either close to the desensitization gate, the activation gate, or in inter‐subunit interfaces. Six of them have pathogenic mutations at equivalent positions in other Cys‐loop receptors, emphasizing the importance of these residues. Also, pathogenic mutations are more common in the pore‐lining helix, consistent with this region being highly constrained for variation in control populations.

Conclusion

Our study reports a novel pathogenic variant in GABRA2, characterizes the regions where pathogenic mutations are in the transmembrane helices, and underscores the value of considering sequence, evolutionary, and structural information as a strategy for variant interpretation of novel missense mutations.

The Clinical Genome and Ancestry Report: An interactive web application for prioritizing clinically implicated variants from genome sequencing data with ancestry composition

Genome sequencing is positioned as a routine clinical work-up for diverse clinical conditions. A commonly used approach to highlight candidate variants with potential clinical implication is to search over locus- and gene-centric knowledge databases. Most web-based applications allow a federated query across diverse databases for a single variant; however, sifting through a large number of genomic variants with combination of filtering criteria is a substantial challenge. Here we describe the Clinical Genome and Ancestry Report (CGAR), an interactive web application developed to follow clinical interpretation workflows by organizing variants into seven categories: (1) reported disease-associated variants, (2) rare- and high-impact variants in putative disease-associated genes, (3) secondary findings that the American College of Medical Genetics and Genomics recommends reporting back to patients, (4) actionable pharmacogenomic variants, (5) focused reports for candidate genes, (6) de novo variant candidates for trio analysis, and (7) germline and somatic variants implicated in cancer risk, diagnosis, treatment and prognosis. For each variant, a comprehensive list of external links to variant-centric and phenotype databases are provided. Furthermore, genotype-derived ancestral composition is used to highlight allele frequencies from a matched population since some disease-associated variants show a wide variation between populations. CGAR is an open-source software and is available at https://tom.tch.harvard.edu/apps/cgar/.

VarSite: Disease variants and protein structure

VarSite is a web server mapping known disease-associated variants from UniProt and ClinVar, together with natural variants from gnomAD, onto protein 3D structures in the Protein Data Bank. The analyses are primarily image-based and provide both an overview for each human protein, as well as a report for any specific variant of interest. The information can be useful in assessing whether a given variant might be pathogenic or benign. The structural annotations for each position in the protein include protein secondary structure, interactions with ligand, metal, DNA/RNA, or other protein, and various measures of a given variant's possible impact on the protein's function. The 3D locations of the disease-associated variants can be viewed interactively via the 3dmol.js JavaScript viewer, as well as in RasMol and PyMOL. Users can search for specific variants, or sets of variants, by providing the DNA coordinates of the base change(s) of interest. Additionally, various agglomerative analyses are given, such as the mapping of disease and natural variants onto specific Pfam or CATH domains. The server is freely accessible to all at: https://www.ebi.ac.uk/thornton-srv/databases/VarSite.