An unusual diagnosis of alpha-mannosidosis with ocular anomalies: Behind the scenes of a hidden copy number variation.

Alpha-mannosidosis is a rare autosomal recessive lysosomal storage disorder caused by biallelic mutations in the MAN2B1 gene and characterized by a wide clinical heterogeneity. Diagnosis for this multisystemic disorder is confirmed by the presence of either a deficiency in the lysosomal enzyme acid alpha-mannosidase or biallelic mutations in the MAN2B1 gene. This diagnosis confirmation is crucial for both clinical management and genetic counseling purposes. Here we describe a late diagnosis of alpha-mannosidosis in a patient presenting with syndromic intellectual disability, and a rare retinopathy, where reverse phenotyping played a pivotal role in interpreting the exome sequencing result. While a first missense variant was classified as a variant of uncertain significance, the phenotype-guided analysis helped us detect and interpret an in-trans apparent alu-element insertion, which appeared to be a copy number variant (CNV) not identified by the CNV caller. A biochemical analysis showing abnormal excretion of urinary mannosyloligosaccharide and an enzyme assay permitted the re-classification of the missense variant to likely pathogenic, establishing the diagnosis of alpha-mannosidosis. This work emphasizes the importance of reverse phenotyping in the context of exome sequencing.

The p.E152K-STIM1 mutation deregulates Ca2+ signaling contributing to chronic pancreatitis.

Since deregulation of intracellular Ca2+ can lead to intracellular trypsin activation, and stromal interaction molecule-1 (STIM1) protein is the main regulator of Ca2+ homeostasis in pancreatic acinar cells, we explored the Ca2+ signaling in 37 STIM1 variants found in three pancreatitis patient cohorts. Extensive functional analysis of one particular variant, p.E152K, identified in three patients, provided a plausible link between dysregulated Ca2+ signaling within pancreatic acinar cells and chronic pancreatitis susceptibility. Specifically, p.E152K, located within the STIM1 EF-hand and sterile α-motif domain, increased the release of Ca2+ from the endoplasmic reticulum in patient-derived fibroblasts and transfected HEK293T cells. This event was mediated by altered STIM1-sarco/endoplasmic reticulum calcium transport ATPase (SERCA) conformational change and enhanced SERCA pump activity leading to increased store-operated Ca2+ entry (SOCE). In pancreatic AR42J cells expressing the p.E152K variant, Ca2+ signaling perturbations correlated with defects in trypsin activation and secretion, and increased cytotoxicity after cholecystokinin stimulation.This article has an associated First Person interview with the first author of the paper.

Coding and noncoding variants in EBF3 are involved in HADDS and simplex autism.

BACKGROUND: Previous research in autism and other neurodevelopmental disorders (NDDs) has indicated an important contribution of protein-coding (coding) de novo variants (DNVs) within specific genes. The role of de novo noncoding variation has been observable as a general increase in genetic burden but has yet to be resolved to individual functional elements. In this study, we assessed whole-genome sequencing data in 2671 families with autism (discovery cohort of 516 families, replication cohort of 2155 families). We focused on DNVs in enhancers with characterized in vivo activity in the brain and identified an excess of DNVs in an enhancer named hs737. RESULTS: We adapted the fitDNM statistical model to work in noncoding regions and tested enhancers for excess of DNVs in families with autism. We found only one enhancer (hs737) with nominal significance in the discovery (p = 0.0172), replication (p = 2.5 × 10-3), and combined dataset (p = 1.1 × 10-4). Each individual with a DNV in hs737 had shared phenotypes including being male, intact cognitive function, and hypotonia or motor delay. Our in vitro assessment of the DNVs showed they all reduce enhancer activity in a neuronal cell line. By epigenomic analyses, we found that hs737 is brain-specific and targets the transcription factor gene EBF3 in human fetal brain. EBF3 is genome-wide significant for coding DNVs in NDDs (missense p = 8.12 × 10-35, loss-of-function p = 2.26 × 10-13) and is widely expressed in the body. Through characterization of promoters bound by EBF3 in neuronal cells, we saw enrichment for binding to NDD genes (p = 7.43 × 10-6, OR = 1.87) involved in gene regulation. Individuals with coding DNVs have greater phenotypic severity (hypotonia, ataxia, and delayed development syndrome [HADDS]) in comparison to individuals with noncoding DNVs that have autism and hypotonia. CONCLUSIONS: In this study, we identify DNVs in the hs737 enhancer in individuals with autism. Through multiple approaches, we find hs737 targets the gene EBF3 that is genome-wide significant in NDDs. By assessment of noncoding variation and the genes they affect, we are beginning to understand their impact on gene regulatory networks in NDDs.

The novel c.634+4A>G splicing variant in RHCE results in weak C and e antigen expression in a pregnant woman originated from Japan.

BACKGROUND: In the RH blood group genes, molecular variants that alter antigen expression with potential clinical relevance are frequently identified and reported in the literature. STUDY DESIGN AND METHODS: A pregnant woman in her first pregnancy, who originates from Japan, was typed by routine serological testing. The RHCE gene was investigated to identify single nucleotide variants (SNVs) and/or structural variants by a commercial platform, Sanger sequencing, and quantitative multiplex PCR of short fluorescent fragments. The haplotypes were determined by sequencing PCR fragments generated from genomic DNA and subcloned into a plasmid vector. Effect on splicing was predicted by bioinformatics tools, including SpliceAI and the splicing module of Alamut. In parallel, functional analysis was carried out by a minigene splicing assay. RESULTS: A patient with no transfusion history was typed RH:1,2w,3,4,5w. An unreported single variant was identified in RHCE intron 4 at the heterozygous state: c.634+4A>G. Minigene splicing assay showed that this SNV decreases significantly the relative abundance of the full-length transcript, in accordance with the predictions made by the Alamut tools, but not SpliceAI, suggesting expression of a normal RhCE protein. CONCLUSION: Overall, the novel RHCE*02(c.634+4A>G) allele alters quantitatively, but not qualitatively, the expression of C and e in the RH blood group system, indicating that the patient is not at risk for alloimmunization and may safely receive C+e+ red blood cell units. This report illustrates the relevance of functional assays for the interpretation of rare variants and, specifically, how it may help guide transfusion management in patients.

A homozygous splice variant in ATP5PO, disrupts mitochondrial complex V function and causes Leigh syndrome in two unrelated families.

Mitochondrial complex V plays an important role in oxidative phosphorylation by catalyzing the generation of ATP. Most complex V subunits are nuclear encoded and not yet associated with recognized Mendelian disorders. Using exome sequencing, we identified a rare homozygous splice variant (c.87+3A>G) in ATP5PO, the complex V subunit which encodes the oligomycin sensitivity conferring protein, in three individuals from two unrelated families, with clinical suspicion of a mitochondrial disorder. These individuals had a similar, severe infantile and often lethal multi-systemic disorder that included hypotonia, developmental delay, hypertrophic cardiomyopathy, progressive epileptic encephalopathy, progressive cerebral atrophy, and white matter abnormalities on brain MRI consistent with Leigh syndrome. cDNA studies showed a predominant shortened transcript with skipping of exon 2 and low levels of the normal full-length transcript. Fibroblasts from the affected individuals demonstrated decreased ATP5PO protein, defective assembly of complex V with markedly reduced amounts of peripheral stalk proteins, and complex V hydrolytic activity. Further, expression of human ATP5PO cDNA without exon 2 (hATP5PO-∆ex2) in yeast cells deleted for yATP5 (ATP5PO homolog) was unable to rescue growth on media which requires oxidative phosphorylation when compared to the wild type construct (hATP5PO-WT), indicating that exon 2 deletion leads to a non-functional protein. Collectively, our findings support the pathogenicity of the ATP5PO c.87+3A>G variant, which significantly reduces but does not eliminate complex V activity. These data along with the recent report of an affected individual with ATP5PO variants, add to the evidence that rare biallelic variants in ATP5PO result in defective complex V assembly, function and are associated with Leigh syndrome.

Heterozygous HMGB1 loss-of-function variants are associated with developmental delay and microcephaly.

13q12.3 microdeletion syndrome is a rare cause of syndromic intellectual disability. Identification and genetic characterization of patients with 13q12.3 microdeletion syndrome continues to expand the phenotypic spectrum associated with it. Previous studies identified four genes within the approximately 300 Kb minimal critical region including two candidate protein coding genes: KATNAL1 and HMGB1. To date, no patients carrying a sequence-level variant or a single gene deletion in HMGB1 or KATNAL1 have been described. Here we report six patients with loss-of-function variants involving HMGB1 and who had phenotypic features similar to the previously described 13q12.3 microdeletion syndrome cases. Common features included developmental delay, language delay, microcephaly, obesity and dysmorphic features. In silico analyses suggest that HMGB1 is likely to be intolerant to loss-of-function, and previous in vitro data are in line with the role of HMGB1 in neurodevelopment. These results strongly suggest that haploinsufficiency of the HMGB1 gene may play a critical role in the pathogenesis of the 13q12.3 microdeletion syndrome.

Whole exome sequencing, a hypothesis-free approach to investigate recurrent early miscarriage.

RESEARCH QUESTION: Are there genetic determinants shared by unrelated women with unexplained recurrent early miscarriage (REM)? DESIGN: Thirty REM cases and 30 controls were selected with extreme phenotype among women from Eastern Brittany (France), previously enrolled in an incident case-control study on thrombophilic mutations. Cases and controls were selected based on the number of early miscarriages or live births, respectively. Peripheral blood was collected for DNA extraction at initial visit. The burden of low-frequency variants in the coding part of the genes was compared using whole exome sequencing (WES). RESULTS: Cases had 3 to 17 early miscarriages (20 cases: ≥5 previous losses). Controls had 1 to 4 live births (20 controls: ≥3 previous live births) and no miscarriages. WES data were available for 29 cases and 30 controls. A total of 209,387 variants were found (mean variant per patient: 59,073.05) with no difference between groups (P = 0.68). The top five most significantly associated genes were ABCA4, NFAM1, TCN2, AL078585.1 and EPS15. Previous studies suggest the involvement of vitamin B12 deficiency in REM. TCN2 encodes for vitamin B12 transporter into cells. Therefore, holotranscobalamin (active vitamin B12) was measured for both cases and controls (81.2 ± 32.1 versus 92.9 ± 34.3 pmol/l, respectively, P = 0.186). Five cases but no controls were below 50 pmol/l (P = 0.052). CONCLUSIONS: This study highlights four new genes of interest in REM, some of which belong to known networks of genes involved in embryonic development (clathrin-mediated endocytosis and ciliary pathway). The study also confirms the involvement of TCN2 (vitamin B12 pathway) in the early first trimester of pregnancy.

Pathogenic and likely pathogenic variants in at least five genes account for approximately 3% of mild isolated nonsyndromic thrombocytopenia.

BACKGROUND: Thrombocytopenia has a variety of different etiologies, both acquired and hereditary. Inherited thrombocytopenia may be associated with other symptoms (syndromic forms) or may be strictly isolated. To date, only about half of all the familial forms of thrombocytopenia have been accounted for in terms of well-defined genetic abnormalities. However, data are limited on the nature and frequency of the underlying causative genetic variants in individuals with mild isolated nonsyndromic thrombocytopenia. STUDY DESIGN AND METHODS: Thirteen known or candidate genes for isolated thrombocytopenia were included in a gene panel analysis in which targeted next-generation sequencing was performed on 448 French blood donors with mild isolated nonsyndromic thrombocytopenia. RESULTS: A total of 68 rare variants, including missense, splice site, frameshift, nonsense, and in-frame variants (all heterozygous) were identified in 11 of the 13 genes screened. Twenty-nine percent (N = 20) of the variants detected were absent from both the French Exome Project and gnomAD exome databases. Using stringent criteria and an unbiased approach, we classified seven predicted loss-of-function variants (three in ITGA2B and four in TUBB1) and four missense variants (one in GP1BA, two in ITGB3 and one in ACTN1) as being pathogenic or likely pathogenic. Altogether, they were found in 13 members (approx. 3%) of our studied cohort. CONCLUSION: We present the results of gene panel sequencing of known and candidate thrombocytopenia genes in mild isolated nonsyndromic thrombocytopenia. Pathogenic and likely pathogenic variants in five known thrombocytopenia genes were identified, accounting for approximately 3% of individuals with the condition.

Implementation of a molecular tumor board at a regional level to improve access to targeted therapy.

BACKGROUND: With the development of precision oncology, Molecular Tumor Boards (MTB) are developing in many institutions. However, the implementation of MTB in routine clinical practice has still not been thoroughly studied. MATERIAL AND METHODS: Since the first drugs approved for targeted therapies, patient tumor samples were centralized to genomic testing platforms. In our institution, all tumor samples have been analyzed since 2014 by Next Generation Sequencing (NGS). In 2015, we established a regional MTB to discuss patient cases with 1 or more alterations identified by NGS, in genes different from those related to drug approval. We conducted a retrospective comparative analysis to study whether our MTB increased the prescriptions of Molecular Targeted Therapies (MTT) and the inclusions of patients in clinical trials with MTT, in comparison with patients with available NGS data but no MTB discussion. RESULTS: In 2014, 86 patients had UGA, but the results were not available to clinicians and not discussed in MTB. During the years 2015 and 2016, 113 patients with an UGA (unreferenced genomic alteration) were discussed in MTB. No patients with an UGA were included in 2014 in a clinical trial, versus 2 (2%) in 2015-2016. 13 patients with an UGA (12%) were treated in 2015-2016 with a MTT whereas in 2014, no patient (p = 0.001). CONCLUSIONS: In this retrospective analysis, we showed that the association of large-scale genomic testing and MTB was feasible, and could increase the prescription of MTT. However, in routine clinical practice, the majority of patients with UGA still do not have access to MTT.

Characterization of GJB2 cis-regulatory elements in the DFNB1 locus.

Although most disease-causing variants are within coding region of genes, it is now well established that cis-acting regulatory sequences, depending on 3D-chromatin organization, are required for temporal and spatial control of gene expression. Disruptions of such regulatory elements and/or chromatin conformation are likely to play a critical role in human genetic disease. Hence, recurrent monoallelic cases, who present the most common hereditary type of nonsyndromic hearing loss (i.e., DFNB1), carry only one identified pathogenic allele. This strongly suggests the presence of uncharacterized distal cis-acting elements in the missing allele. Here within, we study the spatial organization of a large DFNB1 locus encompassing the gap junction protein beta 2 (GJB2) gene, the most frequently mutated gene in this inherited hearing loss phenotype, with the chromosome conformation capture carbon copy technology (5C). By combining this approach with functional activity reporter assays and mapping of CCCTC-binding factor (CTCF) along the DFNB1 locus, we identify a novel set of cooperating GJB2 cis-acting elements and suggest a DFNB1 three-dimensional looping regulation model.

Molecular basis of weak D expression in the Indian population and report of a novel, predominant variant RHD allele.

BACKGROUND: The Rh blood group system is the most polymorphic system and is implicated in hemolytic transfusion reaction and hemolytic disease of the fetus and newborn. Molecular genetics of the RH genes have been extensively studied in Caucasians, Africans, and East Asians and the variant alleles giving rise to weak and partial D phenotypes have been reported. However, limited genetic studies have been carried out in the large Indian population, even though the variability of Rh expression has been documented. STUDY DESIGN AND METHODS: In this study we sought to characterize the molecular bases of weak D expression in Indians. RHD gene in samples presenting with a weak D phenotype by serologic analyses (n = 223) was genotyped by conventional molecular approaches. RESULTS: In addition to referenced and novel single-nucleotide variations, a novel approximately 12-kb duplication event, including Exon 3, was identified predominantly in variant D samples (130/223, 58.3%) and characterized at the nucleotide sequence level. Functional analyses suggested that this genetic variation quantitatively affects the expression of the normal transcript and then subsequently the expression of the normal RhD protein. CONCLUSION: We describe a major novel, variant RHD allele in Indians that can be easily identified routinely by implementing a simple genotyping assay. Although we may consider this variation as a weak partial D variant, further studies and observations are needed to confirm the same. These findings may contribute to improve significantly Rh blood group diagnostics in more than one billion Indians.

Sequential analysis of 18 genes in polycythemia vera and essential thrombocythemia reveals an association between mutational status and clinical outcome.

Philadelphia-negative classical myeloproliferative neoplasms (MPN) are clonal diseases characterized by driver mutations of JAK2, MPL, or CALR. Additional mutations may occur in epigenetic regulators, signaling, or splicing genes that may be useful in the prognostic assessment of MPN patients. In primary myelofibrosis, molecular-based prognostic scoring systems have been recently proposed, but few data are available to date for polycythemia vera (PV) and essential thrombocythemia (ET). In this study, we used a next generation sequencing-based 18-gene panel in 50 JAK2V617F positive PV and JAK2V617F positive ET patients from an institutional cohort investigated at diagnosis and at 3-year follow-up (3y). Disease progression at 3y was defined by a composite criterion. Patients (28 PV and 22 ET) were included according to their clinical status, with or without disease progression. At diagnosis, we found 28 additional mutations in 21 of the 50 patients. Patients with disease progression were more likely to have at least one additional mutation. There was no difference between PV and ET. All patients with two or more additional mutations exhibited disease progression at 3y. No novel mutations appeared at 3y. The allele burden increase by at least one mutation at 3y was more frequent in patients with disease progression. Our data suggest that screening for additional mutations in PV and ET could identify patients at a higher risk of disease progression. © 2017 Wiley Periodicals, Inc.

Extensive functional analyses of RHD splice site variants: Insights into the potential role of splicing in the physiology of Rh.

BACKGROUND: Among more than 300 mutated alleles identified so far within the RHD gene, almost 40 are assumed to alter cellular splicing and therefore may have a direct effect on Rh phenotype both at the quantitative and at the qualitative levels. Functional data are, however, mostly unavailable to assess the direct involvement of splicing defect in the underlying physiology. STUDY DESIGN AND METHODS: We generated plasmid constructs to carry out an exhaustive investigation of 38 RHD variants located within or in the vicinity of exon-intron junctions by a minigene splicing assay, further characterized the transcript structures by sequencing, and identified cryptic sites activated by the genetic defect. Bioinformatics predictions were carried out in parallel and compared with the functional data. RESULTS: For the first time we demonstrate that a product including the full-length Exon 9 is transcribed in the presence of the c.1227G>A substitution frequently carried by Asians with DEL phenotype and confirmed that splicing is altered in the RHD*weak D Type 2 allele, a rare variant most commonly found in Caucasians. CONCLUSION: Overall we 1) show significant correlation between functional analyses, bioinformatics predictions, and phenotypes, when available, especially for variants in close proximity of the consensus splice sites; 2) classify the variations as splicing or nonsplicing variants; and 3) provide functional data to further improve bioinformatics splicing tools. Conversely assessment of seven silent exonic variants was mainly inconclusive.

Insights into RHCE Molecular Analysis in Samples with Partial D Variants: the Experience of Western France.

BACKGROUND: Although systematic blood group genotyping of patients/donors is virtually possible, serological studies remain the gold standard to identify samples of clinical interest that may be further genotyped. In this context, we sought to identify variant D alleles that are likely to be clinically relevant in terms of other Rh antigens in a subset of population genotyped in Western France. METHODS: Samples presenting with the RHD*weak D type 4.2.2 allele (n = 47) were selected for the study. RHCE exons 1-7 were directly sequenced, and expression of Rh antigens was predicted on the basis of the molecular data. RESULTS: Of the 47 samples tested, 19 (40.4%) were predicted to be of potential clinical interest. Moreover, we could show that selecting the samples to be genotyped by the nature of their variant D allele (i.e., RHD*weak D type 4.2.2 allele) rather than by their Duffy-null status appears to increase significantly the likelihood of identifying clinically relevant individuals for Rh status. CONCLUSION: On the basis of our findings we suggest that all individuals genotyped as weak D type 4.2.2 should be systematically screened for RHCE variants by molecular analysis on a routine basis.

Hereditary pancreatitis caused by triplication of the trypsinogen locus.

Hereditary pancreatitis has been reported to be caused by 'gain-of-function' missense mutations in the cationic trypsinogen gene (PRSS1). Here we report the triplication of a approximately 605-kb segment containing the PRSS1 gene on chromosome 7 in five families with hereditary pancreatitis. This triplication, which seems to result in a gain of trypsin through a gene dosage effect, represents a previously unknown molecular mechanism causing hereditary pancreatitis.

A degradation-sensitive anionic trypsinogen (PRSS2) variant protects against chronic pancreatitis.

Chronic pancreatitis is a common inflammatory disease of the pancreas. Mutations in the genes encoding cationic trypsinogen (PRSS1) and the pancreatic secretory trypsin inhibitor (SPINK1) are associated with chronic pancreatitis. Because increased proteolytic activity owing to mutated PRSS1 enhances the risk for chronic pancreatitis, mutations in the gene encoding anionic trypsinogen (PRSS2) may also predispose to disease. Here we analyzed PRSS2 in individuals with chronic pancreatitis and controls and found, to our surprise, that a variant of codon 191 (G191R) is overrepresented in control subjects: G191R was present in 220/6,459 (3.4%) controls but in only 32/2,466 (1.3%) affected individuals (odds ratio 0.37; P = 1.1 x 10(-8)). Upon activation by enterokinase or trypsin, purified recombinant G191R protein showed a complete loss of trypsin activity owing to the introduction of a new tryptic cleavage site that renders the enzyme hypersensitive to autocatalytic proteolysis. In conclusion, the G191R variant of PRSS2 mitigates intrapancreatic trypsin activity and thereby protects against chronic pancreatitis.

Next-generation sequencing is a credible strategy for blood group genotyping.

Although several medium/high-throughput tools have been engineered for molecular analysis of blood group genes, they usually rely on the targeting of single nucleotide polymorphisms, while other variants remain unidentified. To circumvent this limitation a strategy for genotyping blood group genes by next-generation sequencing (NGS) was set up. Libraries consisting of exons, flanking introns and untranslated regions of 18 genes involved in 15 blood systems were generated by the Ion AmpliSeq(™) Library Kit 2.0 and by fragmenting polymerase chain reaction products, normalized by two different approaches, mixed and sequenced by the Ion Torrent Personal Genome Machine (PGM(™) ) Sequencer. In our conditions, defined to limit both intra- and inter-sample variability, sequences from mixed libraries were read in a single run for a total coverage of 86·03% of the coding DNA sequences, including all loci defining the most clinically relevant antigens in all genes, except ABO. Importantly, the challenging attempt to generate gene-specific data for the homologous genes was successful. This work, which combines two complementary approaches to generate libraries, defines technical conditions for genotyping blood group genes, illustrates that NGS is suitable for such an application and suggests that, after automation, this novel tool could be used for molecular typing at the laboratory level.