Missense Changes Research Articles

A recurrent, de novo pathogenic variant in ARPC4 disrupts actin filament formation and causes microcephaly and speech delay

SummaryWe report seven affected individuals from six families with a recurrent, de novo variant in the ARPC4 gene (c.472C>T [p.Arg158Cys (GenBank: NM_005718.4)]). Core features in affected individuals include microcephaly, mild motor delays, and significant speech impairment. ARPC4 is a core subunit of the actin-related protein (ARP2/3) complex, which catalyzes the formation of F-actin networks. We show that the recurrent ARPC4 missense change is associated with a decreased amount of F-actin in cells from two affected individuals. Taken together, our results implicate heterozygous ARPC4 missense variants as a cause of neurodevelopmental disorders and microcephaly.

Spectrum and characterization of bi-allelic variants in MMAB causing cblB-type methylmalonic aciduria

Pathogenic variants in MMAB cause cblB-type methylmalonic aciduria, an autosomal-recessive disorder of propionate metabolism. MMAB encodes ATP:cobalamin adenosyltransferase, using ATP and cob(I)alamin to create 5’-deoxyadenosylcobalamin (AdoCbl), the cofactor of methylmalonyl-CoA mutase (MMUT). We identified bi-allelic disease-causing variants in MMAB in 97 individuals with cblB-type methylmalonic aciduria, including 33 different and 16 novel variants. Missense changes accounted for the most frequent pathogenic alleles (p.(Arg186Trp), N = 57; p.(Arg191Trp), N = 19); while c.700C > T (p.(Arg234*)) was the most frequently identified truncating variant (N = 14). In fibroblasts from 76 affected individuals, the ratio of propionate incorporation in the presence and absence of hydroxocobalamin (PI ratio) was associated to clinical cobalamin responsiveness and later disease onset. We found p.(Arg234*) to be associated with cobalamin responsiveness in vitro, and clinically with later onset; p.(Arg186Trp) and p.(Arg191Trp) showed no clear cobalamin responsiveness and early onset. Mapping these and novel variants onto the MMAB structure revealed their potential to affect ATP and AdoCbl binding. Follow-up biochemical characterization of recombinant MMAB identified its three active sites to be equivalent for ATP binding, determined by fluorescence spectroscopy (Kd = 21 µM) and isothermal calorimetry (Kd = 14 µM), but function as two non-equivalent AdoCbl binding sites (Kd1 = 0.55 μM; Kd2 = 8.4 μM). Ejection of AdoCbl was activated by ATP (Ka = 24 µM), which was sensitized by the presence of MMUT (Ka = 13 µM). This study expands the landscape of pathogenic MMAB variants, provides association of in vitro and clinical responsiveness, and facilitates insight into MMAB function, enabling better disease understanding.

Biallelic ZNFX1 variants are associated with a spectrum of immuno-hematological abnormalities.

Biallelic changes in the ZNFX1 gene have been recently reported to cause severe familial immunodeficiency. Through a search of our bio/databank with information from genetic testing of >55 000 individuals, we identified nine additional patients from seven families with six novel homozygous ZNFX1 variants. Consistent with the previously described phenotype, our patients suffered from monocytosis, thrombocytopenia, hepatosplenomegaly, recurrent infections, and lymphadenopathy. The two most severely affected probands also had renal involvement and clinical presentations compatible with hemophagocytic lymphohistiocytosis. The disease was less lethal among our patients than previously reported. We identified two missense changes, two variants predicted to result in complete protein loss through nonsense-mediated decay, and two frameshift changes that likely introduce a truncation. Our findings (i) independently confirm the role of ZNFX1 in primary genetic immunodeficiency, (ii) expand the genetic and clinical spectrum of ZNFX1-related disease, and (iii) illustrate the utility of large, well-curated, and continually updated genotype-phenotype databases in resolving molecular diagnoses of patients with initially negative genetic testing findings.

Assessing Prevalence and Carrier Frequency of Succinic Semialdehyde Dehydrogenase Deficiency.

Pathogenic variants in ALDH5A1 cause succinic semialdehyde dehydrogenase (SSADH) deficiency, with >180 cases reported worldwide. However, a nonspecific neurologic presentation and inconsistent variant nomenclature have limited diagnoses. In this study, pathogenic variants in ALDH5A1 were curated and variant prevalence assessed in the Genome Aggregation Database (gnomAD) to determine a minimum carrier frequency and to estimate disease prevalence. Stringent population variant analysis, including 98 reported disease-associated ALDH5A1 variants, indicates a pan-ethnic carrier frequency of ∼1/340, supporting a prevalence of SSADH deficiency of ∼1/460 000 worldwide, with highest carrier frequencies observed in East Asian and South Asian populations. Because heterozygous loss of function alleles are rare in gnomAD and >60% of reported disease-causing variants were missense changes that were not present in gnomAD, the pan-ethnic carrier frequency for SSADH deficiency is likely not fully represented in this study. Additional analyses to investigate the potential impact of more common ALDH5A1 variants with reduced but not deficient enzyme activity, including analysis in diverse populations, are needed to fully assess the prevalence of this ultra-rare disease.

The Clinical Genetics of Hemophilia B (Factor IX Deficiency).

Hemophilia B (HB) is a bleeding disorder caused by deficiency of or defect in blood coagulation factor IX (FIX) inherited in an X-linked manner. It results from one of over 1000 known pathogenic variants in the FIX gene, F9; missense and frameshift changes predominate. Although primarily males are affected with HB, heterozygous females may have excessive bleeding due to random or non-random X chromosome inactivation; in addition, homozygous, compound heterozygous, and hemizygous females have been reported. Somatic and germinal mosaicism for F9 variants has been observed. Development of antibodies to FIX treatment products (inhibitors) is rare and related to the type of causative variant present. Treatment is with products produced by recombinant DNA technology, and gene therapy is in clinical trials. Genetic counseling with up-to-date information is warranted for heterozygotes, potential heterozygotes, and men and women affected with HB.

An Updated PAH Mutational Spectrum of Phenylketonuria in Mexican Patients Attending a Single Center: Biochemical, Clinical-Genotyping Correlations.

Establishing the genotypes of patients with hyperphenylalaninemia (HPA)/phenylketonuria (PKU, MIM#261600) has been considered a cornerstone for rational medical management. However, knowledge of the phenylalanine hydroxylase gene (PAH) mutational spectrum in Latin American populations is still limited. Herein, we aim to update the mutational PAH spectrum in the largest cohort of HPA/PKU Mexican patients (N = 124) reported to date. The biallelic PAH genotype was investigated by Sanger automated sequencing, and genotypes were correlated with documented biochemical phenotypes and theoretical tetrahydrobiopterin (BH4) responsiveness. Patients were biochemically classified as having classic PKU (50%, 62/124), mild PKU (20.2%, 25/124) and mild HPA (29.8%, 37/124). Furthermore, 78.2% of the included patients (97/124) were identified by newborn screening. A total of 60 different pathogenic variants were identified, including three novel ones (c. 23del, c. 625_626insC and c. 1315 + 5_1315 + 6insGTGTAACAG), the main categories being missense changes (58%, 35/60) and those affecting the catalytic domain (56.6%, 34/60), and c. 60 + 5G > T was the most frequent variant (14.5%, 36/248) mainly restricted (69.2%) to patients from the central and western parts of Mexico. These 60 types of variants constituted 100 different biallelic PAH genotypes, with the predominance of compound-heterozygous ones (96/124, 77%). The expected BH4 responsiveness based on the PAH genotype was estimated in 52% of patients (65/124), mainly due to the p. (Val388Met) (rs62516101) allele. Instead, our study identified 27 null variants with an allelic phenotype value of zero, with a predominance of c. 60 + 5G > T, which predicts the absence of BH4 responsiveness. An identical genotype reported in BIOPKUdb was found in 92/124 (74%) of our patients, leading to a genotype–phenotype concordance in 80/92 (86.9%) of them. The high number of variants found confirms the heterogeneous and complex mutational landscape of HPA/PKU in Mexico.

Retinal Degeneration Associated With RPGRIP1: A Review of Natural History, Mutation Spectrum, and Genotype-Phenotype Correlation in 228 Patients.

Purpose: RPGRIP1 encodes a ciliary protein expressed in the photoreceptor connecting cilium. Mutations in this gene cause ∼5% of Leber congenital amaurosis (LCA) worldwide, but are also associated with cone–rod dystrophy (CRD) and retinitis pigmentosa (RP) phenotypes. Our purpose was to clinically characterize RPGRIP1 patients from our cohort, collect clinical data of additional RPGRIP1 patients reported previously in the literature, identify common clinical features, and seek genotype–phenotype correlations.Methods: Clinical data were collected from 16 patients of our cohort and 212 previously reported RPGRIP1 patients and included (when available) family history, best corrected visual acuity (BCVA), refraction, comprehensive ocular examination, optical coherence tomography (OCT) imaging, visual fields (VF), and full-field electroretinography (ffERG).Results: Out of 228 patients, the majority (197, 86%) were diagnosed with LCA, 18 (7%) with RP, and 13 (5%) with CRD. Age of onset was during early childhood (n = 133, average of 1.7 years). All patients but 6 had moderate hyperopia (n = 59, mean of 4.8D), and average BCVA was 0.06 Snellen (n = 124; only 10 patients had visual acuity [VA] > 0.10 Snellen). On funduscopy, narrowing of blood vessels was noted early in life. Most patients had mild bone spicule-like pigmentation starting in the midperiphery and later encroaching upon the posterior pole. OCT showed thinning of the outer nuclear layer (ONL), while cystoid changes and edema were relatively rare. VF were usually very constricted from early on. ffERG responses were non-detectable in the vast majority of cases. Most of the mutations are predicted to be null (363 alleles), and 93 alleles harbored missense mutations. Missense mutations were identified only in two regions: the RPGR-interacting domain and the C2 domains. Biallelic null mutations are mostly associated with a severe form of the disease, whereas biallelic missense mutations usually cause a milder disease (mostly CRD).Conclusion: Our results indicate that RPGRIP1 biallelic mutations usually cause severe retinal degeneration at an early age with a cone–rod pattern. However, most of the patients exhibit preservation of some (usually low) BCVA for a long period and can potentially benefit from gene therapy. Missense changes appear only in the conserved domains and are associated with a milder phenotype.

Biallelic TMEM260 variants cause truncus arteriosus, with or without renal defects.

Only two families have been reported with biallelic TMEM260 variants segregating with structural heart defects and renal anomalies syndrome (SHDRA). With a combination of genome, exome sequencing and RNA studies, we identified eight individuals from five families with biallelic TMEM260 variants. Variants included one multi-exon deletion, four nonsense/frameshifts, two splicing changes and one missense change. Together with the published cases, analysis of clinical data revealed ventricular septal defects (12/12), mostly secondary to truncus arteriosus (10/12), elevated creatinine levels (6/12), horse-shoe kidneys (1/12) and renal cysts (1/12) in patients. Three pregnancies were terminated on detection of severe congenital anomalies. Six patients died between the ages of 6 weeks and 5 years. Using a range of stringencies, carrier frequency for SHDRA was estimated at 0.0007-0.007 across ancestries. In conclusion, this study confirms the genetic basis of SHDRA, expands its known mutational spectrum and clarifies its clinical features. We demonstrate that SHDRA is a severe condition associated with substantial mortality in early childhood and characterised by congenital cardiac malformations with a variable renal phenotype.

Diagnostic interest of whole-body MRI in early- and late-onset LAMA2 muscular dystrophies: a large international cohort

LAMA2-related muscular dystrophy (LAMA2-RD) encompasses a group of recessive muscular dystrophies caused by mutations in the LAMA2 gene, which codes for the alpha-2 chain of laminin-211 (merosin). Diagnosis is straightforward in the classic congenital presentation with no ambulation and complete merosin deficiency in muscle biopsy, but is far more difficult in milder ambulant individuals with partial merosin deficiency. To investigate the diagnostic utility of muscle imaging in LAMA2-RD using whole-body magnetic resonance imaging (WBMRI). 27 patients (2-62years, 21-80% with acquisition of walking ability and 6 never ambulant) were included in an international collaborative study. All carried two pathogenic mutations, mostly private missense changes. An intronic variant (c.909 + 7A > G) was identified in all the Chilean cases. Three patients (two ambulant) showed intellectual disability, epilepsy, and brain structural abnormalities. WBMRI T1w sequences or T2 fat-saturated images (Dixon) revealed abnormal muscle fat replacement predominantly in subscapularis, lumbar paraspinals, gluteus minimus and medius, posterior thigh (adductor magnus, biceps femoris, hamstrings) and soleus. This involvement pattern was consistent for both ambulant and non-ambulant patients. The degree of replacement was predominantly correlated to the disease duration, rather than to the onset or the clinical severity. A "COL6-like sandwich sign" was observed in several muscles in ambulant adults, but different involvement of subscapularis, gluteus minimus, and medius changes allowed distinguishing LAMA2-RD from collagenopathies. The thigh muscles seem to be the best ones to assess disease progression. WBMRI in LAMA2-RD shows a homogeneous pattern of brain and muscle imaging, representing a supportive diagnostic tool.

Whole genome sequencing in the diagnosis of primary ciliary dyskinesia

BackgroundIt is estimated that 1–13% of cases of bronchiectasis in adults globally are attributable to primary ciliary dyskinesia (PCD) but many adult patients with bronchiectasis have not been investigated for PCD. PCD is a disorder caused by mutations in genes required for motile cilium structure or function, resulting in impaired mucociliary clearance. Symptoms appear in infancy but diagnosis is often late or missed, often due to the lack of a “gold standard” diagnostic tool and non-specific symptoms. Mutations in > 50 genes account for around 70% of cases, with additional genes, and non-coding, synonymous, missense changes or structural variants (SVs) in known genes presumed to account for the missing heritability.MethodsUK patients with no identified genetic confirmation for the cause of their PCD or bronchiectasis were eligible for whole genome sequencing (WGS) in the Genomics England Ltd 100,000 Genomes Project. 21 PCD probands and 52 non-cystic fibrosis (CF) bronchiectasis probands were recruited in Wessex Genome Medicine Centre (GMC). We carried out analysis of single nucleotide variants (SNVs) and SVs in all families recruited in Wessex GMC.Results16/21 probands in the PCD cohort received confirmed (n = 9), probable (n = 4) or possible (n = 3) diagnosis from WGS, although 13/16 of these could have been picked up by current standard of care gene panel testing. In the other cases, SVs were identified which were missed by panel testing. We identified variants in novel PCD candidate genes (IFT140 and PLK4) in 2 probands in the PCD cohort. 3/52 probands in the non-CF bronchiectasis cohort received a confirmed (n = 2) or possible (n = 1) diagnosis of PCD. We identified variants in novel PCD candidate genes (CFAP53 and CEP164) in 2 further probands in the non-CF bronchiectasis cohort.ConclusionsGenetic testing is an important component of diagnosing PCD, especially in cases of atypical disease history. WGS is effective in cases where prior gene panel testing has found no variants or only heterozygous variants. In these cases it may detect SVs and is a powerful tool for novel gene discovery.

NEW GENES, NEW TECHNIQUES IN NEUROMUSCULAR DISORDERS: O.10 Loss of function mutations in DNAJB4 cause a myopathy with early respiratory failure

Through next generation sequencing of patients with unexplained muscle disease, we identified four individuals from three unrelated families carrying homozygous variants in the DNAJB4 gene. Clinical presentation was similar with acute respiratory failure leading to life threatening hypercapnia and artificial ventilation. In addition, all patients presented some degree of spine rigidity. Disease onset was variable, between the 1st and 4th decade of life. MRI showed a similar selective pattern of muscle involvement. DNAJB4 is a heat shock protein belonging to the DNAJ/HSP40 family. DNAJ proteins are emerging as proteins critical for muscle protein homeostasis and another homologous protein, DNAJB6 – is causative for LGMD D1 (DNAJB6 related). LGMDD1 is characterized by proximal and/or distal muscle weakness, atrophy and protein aggregates, myofibrillar degeneration, and rimmed vacuoles on muscle biopsy. The DNAJB4 variants identified, a stop gain (c.856A>T; p.Lys286*) and two missense changes (c.74G>A; p.Arg25Gln and c.785T>C; p.Leu262Ser) are absent in the control population and are predicted to result in loss of protein, or protein function. Consistent with this, immunoblots from patient muscle tissue and fibroblasts confirmed the absence of DNAJB4 in the p.Lys286* and p.Leu262Ser patients. Functional studies confirm that the p.Lys286* and p.Leu262Ser variants are rapidly degraded in cells. In contrast the p.Arg25Gln is stable but behaves as a loss of function in cellular assays of DNAJB4 chaperone function. In muscle DNAJB4 localizes to the Z-disc and DNAJB4 knockout myotubes and mouse muscle accumulate desmin and other potential client proteins. These data demonstrate a novel chaperonopathy associated with DNAJB4 causing a myopathy with early respiratory failure. DNAJB4 loss of function mutations may lead to the accumulation of DNAJB4 client proteins resulting in muscle dysfunction and degeneration.

NEW GENES, NEW TECHNIQUES IN NEUROMUSCULAR DISORDERS: O.12 Muscle biopsy and RNAseq in the diagnosis of titin related diseases

The ACMG/AMP guidelines for variant interpretation underline that splicing defects are an important mechanism of pathogenicity. In a standard diagnostic pipeline, in silico tools are needed to assess the possible effect on the splice of any identified genetic change (including synonymous variants) and a direct analysis of RNA, cDNA and protein is crucial to support variant classification. Due the sheer size of the TTN transcripts, the interpretation of titin variants is challenging and many unsolved myopathy patients (up to 6%) carry a single heterozygous TTN causative variant. We re-evaluated by SpliceAI all the possible single nucleotide changes in the titin coding sequence, identifying 9,732 variants predicted to alter the splicing (threshold 0.5). Of these, 2,078 are synonymous or missense variants and 317 are nonsense variants. Few previously reported causative ‘missense’ variants are predicted to impact the splicing. At the same time, we analysed by RNA sequencing 14 samples from patients with a confirmed recessive titinopathy, 24 samples from patients with a single heterozygous disease-causing variant in TTN and six from patients with a clinical suspicion of titinopathy. Variants causing a premature stop codon, located in exons out of the M-band, result in a nonsense mediated decay of titin transcripts. Vice versa, the in-depth characterization of 15 splice variants demonstrates that most of them cause in-frame losses or gains, still resulting in a near–full length protein. Using current bioinformatic approaches, RNAseq fails to identify elusive variants in most unsolved patients with heterozygous TTNtv and/or TTN-compatible phenotype, suggesting three possible scenarios: a causative role of the identified missense changes, more complex mechanisms involving heterozygous titin truncating variants, or alternative molecular diagnoses. At the same time, our study demonstrates the clinical importance of TTN transcript and protein studies in the diagnostic pipeline of neuromuscular disorders.

Study to establish genetic association of cardiac conduction defect in Indian patients undergoing pacemaker implantation

Background: The aim was to study the genetic association of cardiac conduction defects (CCDs) by evaluating single nucleotide polymorphism (SNP) in genes of SCN1B and KCNJ2 and to evaluate baseline characteristics between cases and controls.Methods: Case group consisted of 81 individuals with diagnosis of conduction disturbances who underwent permanent pacemaker implantation. The control group consisted of 79 unrelated individuals above 18 years of age of the local population not having a present or past personal or family history (first degree relatives) of any cardiac ailment especially CCDs. Isolation of genomic deoxyribonucleic acid (DNA) from all samples was done, genomic DNA was checked to ensure the presence of intact DNA.Results: SCN1B: SNP rs55742440 had no bearing on the protein except in producing a splice variant. SNP rs67701503 does not lie in the splice-site region, thus not having any significance in the regulation of the gene as well. NetGene2 analysis of SNP rs67486287 negates its presence in the splice site. KCNJ2:SNP rs199473653 leads to a missense amino acid change, resulting in homozygous GG variant found in almost equal frequency in both groups. SNP rs199473653 gene has not been reported as a disease causing mutation.Conclusions: The alteration of nucleotide in SCN1B intron (SNP rs55742440, rs67701503, rs67486287) between cases and controls was found to have no odds of affecting the outcome of CCD. There was no variation or alteration in nucleotide bases of KCNJ2 (SNP rs786205813, rs199473653) between the groups.

Phenotypic effects of changes in the FTVTxK region of an Arabidopsis secondary wall cellulose synthase compared with results from analogous mutations in other isoforms.

Understanding protein structure and function relationships in cellulose synthase (CesA), including divergent isomers, is an important goal. Here, we report results from mutant complementation assays that tested the ability of sequence variants of AtCesA7, a secondary wall CesA of Arabidopsis thaliana, to rescue the collapsed vessels, short stems, and low cellulose content of the irx3-1 AtCesA7 null mutant. We tested a catalytic null mutation and seven missense or small domain changes in and near the AtCesA7 FTVTSK motif, which lies near the catalytic domain and may, analogously to bacterial CesA, exist within a substrate "gating loop." A low-to-high gradient of rescue occurred, and even inactive AtCesA7 had a small positive effect on stem cellulose content but not stem elongation. Overall, secondary wall cellulose content and stem length were moderately correlated, but the results were consistent with threshold amounts of cellulose supporting particular developmental processes. Vibrational sum frequency generation microscopy allowed tissue-specific analysis of cellulose content in stem xylem and interfascicular fibers, revealing subtle differences between selected genotypes that correlated with the extent of rescue of the collapsing xylem phenotype. Similar tests on PpCesA5 from the moss Physcomitrium (formerly Physcomitrella) patens helped us to synergize the AtCesA7 results with prior results on AtCesA1 and PpCesA5. The cumulative results show that the FTVTxK region is important for the function of an angiosperm secondary wall CesA as well as widely divergent primary wall CesAs, while differences in complementation results between isomers may reflect functional differences that can be explored in further work.

Genotype and phenotype distribution of 435 patients with Charcot–Marie–Tooth disease from central south China

The purpose was to provide an overview of genotype and phenotype distribution in a cohort of patients with Charcot-Marie-Tooth disease (CMT) and related disorders from central south China. In all, 435 patients were enrolled and detailed clinical data were collected. Multiplex ligation-dependent probe amplification for PMP22 duplication/deletion and CMT multi-gene panel sequencing were performed. Whole exome sequencing was further applied in the remaining patients who failed to achieve molecular diagnosis. Among the 435 patients, 216had CMT1, 14had hereditary neuropathy with pressure palsies (HNPP), 178had CMT2, 24had distal hereditary motor neuropathy (dHMN) and threehad hereditary sensory and autonomic neuropathy (HSAN). The overall molecular diagnosis rate was 70%: 75.7% in CMT1, 100% in HNPP, 64.6% in CMT2, 41.7% in dHMN and 33.3% in HSAN. The most common four genotypes accounted for 68.9% of molecular diagnosed patients. Relatively frequent causes were missense changes in PMP22 (4.6%) and SH3TC2 (2.3%) in CMT1; and GDAP1 (5.1%), IGHMBP2 (4.5%) and MORC2 (3.9%) in CMT2. Twenty of 160 detected pathogenic variants and the associated phenotypes have not been previously reported. Broad phenotype spectra were observed in six genes, amongst which the pathogenic variants in BAG3 and SPTLC1 were detected in two sporadic patients presenting with the CMT2 phenotype. Our results provided a unique genotypic and phenotypic landscape of patients with CMT and related disorders from central south China, including a relatively high proportion of CMT2 and lower occurrence of PMP22 duplication. The broad phenotype spectra in certain genes have advanced our understanding of CMT.

Dissection of pleiotropic effects of variants in and adjacent to F8 exon 19 and rescue of mRNA splicing and protein function

The pathogenic significance of nucleotide variants commonly relies on nucleotide position within the gene, with exonic changes generally attributed to quantitative or qualitative alteration of protein biosynthesis, secretion, activity, or clearance. However, these changes may exert pleiotropic effects on both protein biology and mRNA splicing due to the overlapping of the amino acid and splicing codes, thus shaping the disease phenotypes. Here, we focused on hemophilia A, in which the definition of F8 variants' causative role and association to bleeding phenotypes is crucial for proper classification, genetic counseling, and management of affected individuals. We extensively characterized a large panel of hemophilia A-causing variants (n = 30) within F8 exon 19 by combining and comparing in silico and recombinant expression analyses. We identified exonic variants with pleiotropic effects and dissected the altered protein features of all missense changes. Importantly, results from multiple prediction algorithms provided qualitative results, while recombinant assays allowed us to correctly infer the likely phenotype severity for 90% of variants. Molecular characterization of pathogenic variants was also instrumental for the development of tailored correction approaches to rescue splicing affecting variants or missense changes impairing protein folding. A single engineered U1snRNA rescued mRNA splicing of nine different variants and the use of a chaperone-like drug resulted in improved factor VIII protein secretion for four missense variants. Overall, dissection of the molecular mechanisms of a large panel of HA variants allowed precise classification of HA-affected individuals and favored the development of personalized therapeutic approaches.

RNA Editing By APOBEC3A Confers Dynamic Mutational Load in Monocytic Leukemias

Background: Although genomic mutations play an important role in the initiation or progression of myeloid leukemias, the role of mutations at the transcript level is less well understood. Myeloid neoplasms have recurrent somatic mutations in genes for pre-mRNA splicing, cohesin complex, histone modification, and DNA methylation. However, the number of somatic mutations identified in myeloid neoplasms (average 9-13) is usually lower than that of solid tumors. Recent analyses of baseline transcriptome sequences in various cancers have revealed abundant adenosine to inosine (A>I) RNA editing events, mainly at intronic and untranslated regions of genes. The role of condition-specific RNA editing in cancer, especially of the cytosine to uracil (C>U) type, however, remains essentially unknown. We recently discovered antiviral restriction enzyme APOBEC3A (A3A) cytidine deaminase as a novel cellular site-specific C>U RNA editing enzyme in monocytes. RNA editing by A3A is dynamically induced by interferons and/or severe hypoxia in monocytes, resulting in stop codon gains or missense changes in transcripts of scores of genes. The effect of hypoxia in inducing RNA editing by A3A can be mimicked by inhibition of mitochondrial respiration. Whether A3A-mediated RNA editing occurs in monocytic leukemias, including acute monocytic leukemia (AMoL) and chronic myelomonocytic leukemia (CMML), is unknown. We hypothesize that dynamic C>U RNA editing by A3Acontributes to mutational load of leukemic cells under stressful conditions of hypoxia and interferon exposure. In this study, we examined transcript sequences of monocytic leukemias for evidence of inducible RNA editing. We also examined the role of mitochondrial respiratory inhibition, which induces RNA editing by A3A, in growth of leukemic progenitors of CMML.Methods: Monocytic leukemia samples are obtained from the RPCI Leukemia Bank. Standard diagnostic information included classification (WHO 2008 and FAB), fraction of monocytic component (mature monocytes, promonocytes or monoblasts by flow cytometry and/or aspirate differential count), cytogenetic analyses (conventional karyotype and targeted FISH), and survival data. Leukemic BM mononuclear cells (BM-MNC) prepared by the Ficoll gradient method are thawed and cultured in hypoxia (1% O2) with interferon type 1 (IFN-1). DNA and RNA are isolated. AMoL samples are subjected to nextgen high throughput exome and RNA sequencing. Once variants from DNA and RNA of corresponding samples are obtained, the set difference of variants is extracted to keep those detected on RNA samples. We also examined 3 CMML bone marrow samples for evidence of hypoxia-inducible RNA editing by analyzing selected RNA editing sites. Colony forming unit (CFU) assays are also performed in CMML samples.Results: We find that the A3A gene is robustly expressed inboth in AMoL and CMML. Hypoxia (with or without IFN-1) induced A3A-mediated RNA editing in AMoL and CMML samples. We previously showed that several RNA edited genes by A3A in normal monocytes function in recurrently altered pathways in myeloid leukemias. Sanger sequencing confirmed RNA editing of such genes in monocytic leukemia samples as well (Table 1), suggesting that epitranscriptomic sequence changes may contribute to molecular defects driving leukomogenesis. RNA seq analysis of AMoL samples (n=3) showed that dozens of genes acquire nonsense or missense mutations at the transcript level under the conditions of severe hypoxia and IFN-1. CFU assays of CMML samples showed that inhibition of mitochondrial respiration, which triggers A3A-mediated RNA editing, inhibits the growth of progenitor cells (Figure 1).Conclusions: We show that A3A-mediated RNA editing is active in monocytic leukemias and targets genes functioning in pathways recurrently altered in myeloid leukemias. RNA editing by A3A may contribute to growth arrest and survival of leukemic progenitors under certain stress conditions. To our knowledge, these results provide the first experimental evidence that conditional RNA editing by A3A confers substantial mutational load in myeloid leukemias. Such RNA editing events will dynamically increase the total mutational load in monocytic leukemias, and may create neo-antigens for immune therapy by checkpoint inhibition. [Display omitted] DisclosuresNo relevant conflicts of interest to declare.

Role of MAX As a Tumor Suppressor Driver Gene in Multiple Myeloma

MYC associated factor X (MAX) is part of the basic-helix-loop-helix (bHLH) leucine zipper transcription factors. MAX forms homo- and hetero-dimers with other bHLH-like factors (e.g. MAD, MXL1 and MYC). The different dimers compete for CACGT DNA E-box consensus sequence (canonical binding site) and other binding sites (non-canonical binding sites) modulating transcriptional activation. The regulation of this process, among others, is based on the different affinity of the formed hetero-dimers to the canonical and non-canonical DNA binding sites and epigenetic modifications in these sequences. In the present study we aim to characterize and understand the role of MAX somatic mutations in the development of Multiple Myeloma (MM).Targeted amplicon sequencing was performed using the M3P (v2.0 or v3.0) gene selection with an average sequencing depth of 700X. We combined data from the different M3P cohorts for a total of 501 M3P patients, 337 newly diagnosed and 164 pre-treated. We also included data from the CoMMpas trial (804 newly diagnosed MM analyzed by WES with 100X average depth). In total 1.305 MM cases were analyzed. The MAX gene was mutated in 33 MM patients (2,5%). Incidence of mutations identified exceeded 5% variant read frequency (VR) and there was no difference in the incidence of the targeted M3P and the WES CoMMpass approach (M3P: 13 of 501 MM cases mutated [2.6%] vs. CoMMpass: 20 of 804 [2.5%]) or between newly diagnosed and pretreated cases (newly diagnosed M3P: 8 of 337 cases [2.4%] vs. pretreated M3P: 5 of 164; [3.0%]). In total, we identified 42 somatic mutations, with 7 patients showing more than one alteration. E.g. in a 73 y/o newly diagnosed MM patient we identified three missense mutations at known MAX hotspot location (p.E32V, p.R35L, p.R36K), all of them in subclonal range with variant reads of 26%, 21% and 10%, respectively. Each mutation exclusively occurred on a different sequencing read and no reads were shared between the mutations. This suggests parallel evolution of different subclones of the disease, and that mutant MAX was virtually present on all tumor cells (confirmed diploid karyotype by cytogenetics). In our cohort we did not observe incidence differences between newly diagnosed and pretreated patients, thus we believe MAX mutations are unlikely to be therapy-related. However, the existence of more than one variant in the gene in 21% of MAX mutated patients (7/33) may indicate a role in disease development. The loss of the protein activity of MAX may provide a survival advantage to the affected tumor clones, possibly due to the binding of other bHLH transcription factors (e.g. MYC) to the E-box, however underlying mechanisms need to be determined. Supporting this hypothesis, 16/42 MAX mutations (38%) affected the integrity of the protein (nonsense, stop gain/loss, tart loss and splice-site donor). The remaining were missense amino-acid changes located in the bHLH DNA binding domaine, including the hotspots R35 (8 mutations), R36 (6 mutations) and R60 (4 mutations), which, according to Wang et al. (Nucleic acid research, 2017) induce a loss of the DNA-binding capacity of the MAX protein. Thus, mutations in MAX in MM patients vastly induce the loss of function or the loss of the protein itself (at least 35/42 (85%)). This inactivation was related with disease development, confirming the role of MAX as a tumor suppressor driver gene in MM.Figure: A)Location of somatic missense mutations within MAX gene. All except E143G occurred in the DNA binding domaine. Mutations in a.a. G32, R35, R36 and R60 alter DNA binding properties (Wang d et al. Nucleic Acids Research, 2017). B) Three subclonal mutations in a single patient, located at the DNA binding site. The X-ray structure of MAX residues 22-107 bound to DNA (pdb: 5EYO) is shown. C)The mutations belong to different clones. Mutations were detected in independent sequencing reads, covered by the same amplicon, BAM file visualization in IGV sequencing reads aligned to hg19. [Display omitted] DisclosuresRaab:Heidelberg Pharma GMBH: Research Funding. Sonneveld:Celgene Corporation, Amgen, Janssen, Karyopharm, SkylineDx, PharmaMar: Consultancy; Celgene Corporation, Amgen, Janssen, Karyopharm, PharmaMar, SkylineDx: Honoraria; Celgene, Amgen, Janssen, Karyopharm, Takeda: Consultancy, Honoraria, Research Funding. Stewart:Roche: Consultancy; Amgen: Consultancy; Bristol-Myers Squibb: Consultancy; Celgene: Consultancy; Janssen: Consultancy.

Further Evidence That Defects in Main Thyroid Dysgenesis-Related Genes Are an Uncommon Etiology for Primary Congenital Hypothyroidism in Mexican Patients: Report of Rare Variants in FOXE1, NKX2-5 and TSHR.

Mexico shows a high birth prevalence of congenital hypothyroidism (CH) due to thyroid dysgenesis (TD). PAX8 defects underlie only 1% of these cases and NKX2-1 does not seem to be involved. Here, we analyzed other TD-related genes in 128 non-related Mexican patients (females 77.3%; 6 months to 16.6 years) with non-syndromic CH-TD diagnosis established by clinical evaluation, thyroid hormone serum profiling, and scintigraphy (74%) or ultrasonography (26%). We performed Sanger sequencing of FOXE1, NKX2-5, and TSHR and evaluated copy number variations (CNVs) in TSHR, FOXE1, PAX8, and NKX2-1 by multiplex ligation-dependent probe amplification. Odds ratios for TD risk were explored for FOXE1 polyalanine stretches [polyAla-rs71369530] in cases and controls (N = 116). Five rare missense changes cataloged as benign (NKX2-5:p.(Ala119Ser)-rs137852684), of unknown significance (FOXE1:p.(Ala335Gly)-rs543372757; TSHR:p.(Asp118Asn)-rs1414102266), and likely pathogenic (FOXE1:p.(Gly124Arg)-rs774035532; TSHR:p.(Trp422Arg)-rs746029360) accounted for 1.5% (N = 2/128) of clinically relevant genotypes (supported in part by protein modeling) in CH-TD. No CNVs were identified, nor did polyAla > 14 alanines in FOXE1 significantly protect against TD. The present and previously published data collectively show that small clinically relevant germline variants in PAX8, FOXE1, and TSHR are found in only a very small proportion (2.5%) of isolated CH-TD Mexican patients.

CDKL5 deficiency disorder in males: Five new variants and review of the literature.

The X-linked Cyclin-Dependent Kinase-Like 5 (CDKL5) gene encodes a serine-threonine kinase highly expressed in the developing brain. Loss of function of CDKL5 is pointed out to underlie the CDKL5 Deficiency Disorder (CDD), an X-linked dominant disease characterized by early-onset epileptic encephalopathy and developmental delay, usually affecting females more than males. To the best to our knowledge, only 45 males with CDD have been reported so far. Type and position of CDKL5 variants with different impact on the protein are reported to influence the clinical presentation. X-chromosome inactivation occurring in females and post-zygotic mosaicism in males are also believed to contribute to this variability. Based on these issues, genotype-phenotype correlations are still challenging. Here, we describe clinical features of five additional affected males with unreported CDKL5 variants, expanding the molecular spectrum of the disorder. We also reviewed the clinical profile of the previously reported 45 males with molecularly confirmed CDD. Severe developmental delay, cortical visual impairment, and early-onset refractory epilepsy characterize the CDD picture in males. By assessing the molecular spectrum, we confirm that germ-line truncating CDKL5 variants, equally distributed across the coding sequence, are the most recurrent mutations in CDD, and cause the worsen phenotype. While recurrence and relevance of missense substitutions within C-terminal remain still debated, disease-causing missense changes affecting the N-terminal catalytic domain correlate to a severe clinical phenotype. Finally, our data provide evidence that post-zygotic CDKL5 mosaicism may result in milder phenotypes and, at least in a subset of subjects, in variable response to antiepileptic treatments.