Structural Biology Helps Interpret Variants of Uncertain Significance in Genes Causing Endocrine and Metabolic Disorders

Co-inheritance of Rare Variants of Maturity-Onset Diabetes of the Young (MODY): A Case Report and Review of the Literature

Maturity onset diabetes of the young (MODY) is a rare, monogenic autosomal dominant form of diabetes that is characterized by early-onset, non-insulin-dependent hyperglycemia, strong family history, and is often misdiagnosed as type 1 or type 2 diabetes. Co-inheritance of multiple MODY genes, however, is rare. We describe here a case of MODY involving co-inherited rare variants in the ABCC8 and B-lymphocyte kinase (BLK) genes. A 55-year-old non-obese man with a past medical history of dyslipidemia and premature ischemic heart disease was initially misdiagnosed with type 2 diabetes for more than 18 years. He is a smoker with a strong family history of diabetes affecting both of his parents and most of his siblings. Despite treatment with different oral antihyperglycemics, his diabetes remained uncontrolled with glycated hemoglobin (HBA1c) between 8 and 10% until the addition of gliclazide, which improved his HBA1c to 5.7%. Based on all the previous information, MODY was suspected, and genetic testing was done, which showed rare variants in the BLK and ABCC8 genes and suggested a co-occurrence of MODY11 and MODY12. This case highlights the importance of accurate genetic testing, which is crucial for proper MODY subtyping, enabling tailored treatment strategies and potentially improving patient outcomes. Moreover, the consistent presence of the BLK gene variant in limited cases of co-inheritance raises questions about its causative role in MODY, suggesting a need for further investigation into its clinical significance.

Safety and efficacy of canakinumab treatment for undifferentiated autoinflammatory diseases: the data of a retrospective cohort two-centered study

Introduction

The blockade of interleukine-1 (anakinra and canakinumab) is a well-known highly effective tool for monogenic autoinflammatory diseases (AIDs), such as familial Mediterranean fever, tumor necrosis factor receptor-associated periodic syndrome, hyperimmunoglobulinaemia D syndrome, and cryopyrin-associated periodic syndrome, but this treatment has not been assessed for patients with undifferentiated AIDs (uAIDs). Our study aimed to assess the safety and efficacy of canakinumab for patients with uAIDs.

Methods

Information on 32 patients with uAIDs was retrospectively collected and analyzed. Next-generation sequencing and Federici criteria were used for the exclusion of the known monogenic AID.

Results

The median age of the first episode was 2.5 years (IQR: 1.3; 5.5), that of the disease diagnosis was 5.7 years (IQR: 2.5;12.7), and that of diagnostic delay was 1.1 years (IQR: 0.4; 6.1). Patients had variations in the following genes: IL10, NLRP12, STAT2, C8B, LPIN2, NLRC4, PSMB8, PRF1, CARD14, IFIH1, LYST, NFAT5, PLCG2, COPA, IL23R, STXBP2, IL36RN, JAK1, DDX58, LACC1, LRBA, TNFRSF11A, PTHR1, STAT4, TNFRSF1B, TNFAIP3, TREX1, and SLC7A7. The main clinical features were fever (100%), rash (91%; maculopapular predominantly), joint involvement (72%), splenomegaly (66%), hepatomegaly (59%), lymphadenopathy (50%), myalgia (28%), heart involvement (31%), intestinal involvement (19%); eye involvement (9%), pleuritis (16%), ascites (6%), deafness, hydrocephalia (3%), and failure to thrive (25%). Initial treatment before canakinumab consisted of non-biologic therapies: non-steroidal anti-inflammatory drugs (NSAID) (91%), corticosteroids (88%), methotrexate (38%), intravenous immunoglobulin (IVIG) (34%), cyclosporine A (25%), colchicine (6%) cyclophosphamide (6%), sulfasalazine (3%), mycophenolate mofetil (3%), hydroxychloroquine (3%), and biologic drugs: tocilizumab (62%), sarilumab, etanercept, adalimumab, rituximab, and infliximab (all 3%). Canakinumab induced complete remission in 27 patients (84%) and partial remission in one patient (3%). Two patients (6%) were primary non-responders, and two patients (6%) further developed secondary inefficacy. All patients with partial efficacy or inefficacy were switched to tocilizumab (n = 4) and sarilumab (n = 1). The total duration of canakinumab treatment was 3.6 (0.1; 8.7) years. During the study, there were no reported Serious Adverse Events (SAEs). The patients experienced non-frequent mild respiratory infections at a rate that is similar as before canakinumab is administered. Additionally, one patient developed leucopenia, but it was not necessary to stop canakinumab for this patient.

Conclusion

The treatment of patients with uAIDs using canakinumab was safe and effective. Further randomized clinical trials are required to confirm the efficacy and safety.

Novel Genomic Variants, Atypical Phenotypes and Evidence of a Digenic/Oligogenic Contribution to Disorders/Differences of Sex Development in a Large North African Cohort

In a majority of individuals with disorders/differences of sex development (DSD) a genetic etiology is often elusive. However, new genes causing DSD are routinely reported and using the unbiased genomic approaches, such as whole exome sequencing (WES) should result in an increased diagnostic yield. Here, we performed WES on a large cohort of 125 individuals all of Algerian origin, who presented with a wide range of DSD phenotypes. The study excluded individuals with congenital adrenal hypoplasia (CAH) or chromosomal DSD. Parental consanguinity was reported in 36% of individuals. The genetic etiology was established in 49.6% (62/125) individuals of the total cohort, which includes 42.2% (35/83) of 46, XY non-syndromic DSD and 69.2% (27/39) of 46, XY syndromic DSD. No pathogenic variants were identified in the 46, XX DSD cases (0/3). Variants in the AR, HSD17B3, NR5A1 and SRD5A2 genes were the most common causes of DSD. Other variants were identified in genes associated with congenital hypogonadotropic hypogonadism (CHH), including the CHD7 and PROKR2. Previously unreported pathogenic/likely pathogenic variants (n = 30) involving 25 different genes were identified in 22.4% of the cohort. Remarkably 11.5% of the 46, XY DSD group carried variants classified as pathogenic/likely pathogenic variant in more than one gene known to cause DSD. The data indicates that variants in PLXNA3, a candidate CHH gene, is unlikely to be involved in CHH. The data also suggest that NR2F2 variants may cause 46, XY DSD.

Computational biology insights into genotype-clinical phenotype-protein phenotype relationships between novel SLC26A2 variants identified in inherited skeletal dysplasias

BACKGROUND

Pathogenic variants in the transmembrane sulfate transporter protein SLC26A2 are associated with different phenotypes of inherited chondrodysplasias. As limited data is published from India, in this study we sought to elucidate the molecular basis of inherited chondrodysplasias in an Indian cohort.

METHODS

Molecular screening of 32 fetuses with antenatally diagnosed lethal skeletal dysplasia was performed by next generation sequencing and Sanger sequencing. The genotype-protein phenotype characterization was done using computational biology techniques like homology modelling, stability and pathogenicity predictions.

RESULTS

We identified five rare autosomal recessive SLC26A2 [NM_000112.4] variants, including three homozygous c.796dupA(p.Thr266Asnfs*12), c.1724delA(p.Lys575Serfs*10), and c.1375_1377dup(p.Val459dup) and two heterozygous variants (c.532C > T(p.Arg178*)) and (c.1382C > T(p.Ala461Val)) in compound heterozygous form in a total of four foetuses. Genotype-protein phenotype annotations highlighted that the clinically severe achondrogenesis 1B causative c.796dupA(p.Thr266Asnfs*12) and c.1724delA(p.Lys575Serfs*10)variants impact SLC26A2 protein structure by deletion of the protein core and transmembrane STAS domains, respectively. In clinically moderate atelosteogenesis type 2 phenotype, the c.1382C > T(p.Ala461Val) variant is predicted to distort alpha helix conformation and alter the bonding properties and free energy dynamics of transmembrane domains and the c.532C > T(p.Arg178*) variant results in loss of both core transmembrane and STAS domains of the SLC26A2 protein. The c.1375_1377dup(p.Val459dup) variant identified in clinically milder atelosteogenesis type II-diastrophic dysplasia spectrum lethal phenotype is predicted to decrease the Qualitative Model Energy Analysis (QMean), which affects major geometrical aspects of the SLC26A2 protein structure.

CONCLUSION

We expand the spectrum of SLC26A2 related lethal chondrodysplasia and report three novel variants correlating clinical severity and protein phenotype within the lethal spectrum of this rare dysplasia. We demonstrate the relevance of structural characterization to aid novel variant reclassification to provide better prenatal management and reproductive options to families with lethal antenatal skeletal disorder.

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.

Missense3D-DB web catalogue: an atom-based analysis and repository of 4M human protein-coding genetic variants

The interpretation of human genetic variation is one of the greatest challenges of modern genetics. New approaches are urgently needed to prioritize variants, especially those that are rare or lack a definitive clinical interpretation. We examined 10,136,597 human missense genetic variants from GnomAD, ClinVar and UniProt. We were able to perform large-scale atom-based mapping and phenotype interpretation of 3,960,015 of these variants onto 18,874 experimental and 84,818 in house predicted three-dimensional coordinates of the human proteome. We demonstrate that 14% of amino acid substitutions from the GnomAD database that could be structurally analysed are predicted to affect protein structure (n = 568,548, of which 566,439 rare or extremely rare) and may, therefore, have a yet unknown disease-causing effect. The same is true for 19.0% (n = 6266) of variants of unknown clinical significance or conflicting interpretation reported in the ClinVar database. The results of the structural analysis are available in the dedicated web catalogue Missense3D-DB ( http://missense3d.bc.ic.ac.uk/ ). For each of the 4 M variants, the results of the structural analysis are presented in a friendly concise format that can be included in clinical genetic reports. A detailed report of the structural analysis is also available for the non-experts in structural biology. Population frequency and predictions from SIFT and PolyPhen are included for a more comprehensive variant interpretation. This is the first large-scale atom-based structural interpretation of human genetic variation and offers geneticists and the biomedical community a new approach to genetic variant interpretation.

Can Predicted Protein 3D Structures Provide Reliable Insights into whether Missense Variants Are Disease Associated?

Knowledge of protein structure can be used to predict the phenotypic consequence of a missense variant. Since structural coverage of the human proteome can be roughly tripled to over 50% of the residues if homology-predicted structures are included in addition to experimentally determined coordinates, it is important to assess the reliability of using predicted models when analyzing missense variants. Accordingly, we assess whether a missense variant is structurally damaging by using experimental and predicted structures. We considered 606 experimental structures and show that 40% of the 1965 disease-associated missense variants analyzed have a structurally damaging change in the mutant structure. Only 11% of the 2134 neutral variants are structurally damaging. Importantly, similar results are obtained when 1052 structures predicted using Phyre2 algorithm were used, even when the model shares low (<40%) sequence identity to the template. Thus, structure-based analysis of the effects of missense variants can be effectively applied to homology models. Our in-house pipeline, Missense3D, for structurally assessing missense variants was made available at http://www.sbg.bio.ic.ac.uk/~missense3d.