Microsatellite instability at U2AF-binding polypyrimidic tract sites perturbs alternative splicing during colorectal cancer initiation.

BACKGROUND: Microsatellite instability (MSI) due to mismatch repair deficiency (dMMR) is common in colorectal cancer (CRC). These cancers are associated with somatic coding events, but the noncoding pathophysiological impact of this genomic instability is yet poorly understood. Here, we perform an analysis of coding and noncoding MSI events at the different steps of colorectal tumorigenesis using whole exome sequencing and search for associated splicing events via RNA sequencing at the bulk-tumor and single-cell levels. RESULTS: Our results demonstrate that MSI leads to hundreds of noncoding DNA mutations, notably at polypyrimidine U2AF RNA-binding sites which are endowed with cis-activity in splicing, while higher frequency of exon skipping events are observed in the mRNAs of MSI compared to non-MSI CRC. At the DNA level, these noncoding MSI mutations occur very early prior to cell transformation in the dMMR colonic crypt, accounting for only a fraction of the exon skipping in MSI CRC. At the RNA level, the aberrant exon skipping signature is likely to impair colonic cell differentiation in MSI CRC affecting the expression of alternative exons encoding protein isoforms governing cell fate, while also targeting constitutive exons, making dMMR cells immunogenic in early stage before the onset of coding mutations. This signature is characterized by its similarity to the oncogenic U2AF1-S34F splicing mutation observed in several other non-MSI cancer. CONCLUSIONS: Overall, these findings provide evidence that a very early RNA splicing signature partly driven by MSI impairs cell differentiation and promotes MSI CRC initiation, far before coding mutations which accumulate later during MSI tumorigenesis.

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.

Microsatellite Instability and Aberrant Pre-mRNA Splicing: How Intimate Is It?

Cancers that belong to the microsatellite instability (MSI) class can account for up to 15% of all cancers of the digestive tract. These cancers are characterized by inactivation, through the mutation or epigenetic silencing of one or several genes from the DNA MisMatch Repair (MMR) machinery, including MLH1, MLH3, MSH2, MSH3, MSH6, PMS1, PMS2 and Exo1. The unrepaired DNA replication errors turn into mutations at several thousand sites that contain repetitive sequences, mainly mono- or dinucleotides, and some of them are related to Lynch syndrome, a predisposition condition linked to a germline mutation in one of these genes. In addition, some mutations shortening the microsatellite (MS) stretch could occur in the 3'-intronic regions, i.e., in the ATM (ATM serine/threonine kinase), MRE11 (MRE11 homolog) or the HSP110 (Heat shock protein family H) genes. In these three cases, aberrant pre-mRNA splicing was observed, and it was characterized by the occurrence of selective exon skipping in mature mRNAs. Because both the ATM and MRE11 genes, which as act as players in the MNR (MRE11/NBS1 (Nibrin)/RAD50 (RAD50 double strand break repair protein) DNA damage repair system, participate in double strand breaks (DSB) repair, their frequent splicing alterations in MSI cancers lead to impaired activity. This reveals the existence of a functional link between the MMR/DSB repair systems and the pre-mRNA splicing machinery, the diverted function of which is the consequence of mutations in the MS sequences.

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.

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.

Whole-Genome Sequencing Improves the Diagnosis of DFNB1 Monoallelic Patients.

Hearing loss is the most common sensory defect, due in most cases to a genetic origin. Variants in the GJB2 gene are responsible for up to 30% of non-syndromic hearing loss. Today, several deafness genotypes remain incomplete, confronting us with a diagnostic deadlock. In this study, whole-genome sequencing (WGS) was performed on 10 DFNB1 patients with incomplete genotypes. New variations on GJB2 were identified for four patients. Functional assays were realized to explore the function of one of them in the GJB2 promoter and confirm its impact on GJB2 expression. Thus, in this study WGS resolved patient genotypes, thus unlocking diagnosis. WGS afforded progress and bridged some gaps in our research.

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.

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.

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.

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.

The ratio of maximal handgrip force and maximal cycloergometry power as a diagnostic tool to screen for metabolic myopathies.

Metabolic myopathies comprise a diverse group of inborn errors of intermediary metabolism affecting skeletal muscle, and often present clinically as an inability to perform normal exercise. Our aim was to use the maximal mechanical performances achieved during two functional tests, isometric handgrip test and cycloergometer, to identify metabolic myopathies among patients consulting for exercise-induced myalgia. Eighty-three patients with exercise-induced myalgia and intolerance were evaluated, with twenty-three of them having a metabolic myopathy (McArdle, n = 9; complete myoadenylate deaminase deficiency, n = 10; respiratory chain deficiency, n = 4) and sixty patients with non-metabolic myalgia. In all patients, maximal power (MP) was determined during a progressive exercise test on a cycloergometer and maximal voluntary contraction force (MVC) was assessed using a handgrip dynamometer. The ratio between percent-predicted values for MVC and MP was calculated for each subject (MVC%pred:MP%pred ratio). In patients with metabolic myopathy, the MVC%pred:MP%pred ratio was significantly higher compared to non-metabolic myalgia (1.54 ± 0.62 vs. 0.92 ± 0.25; p < 0.0001). ROC analysis of MVC%pred:MP%pred ratio showed AUC of 0.843 (0.758-0.927, 95% CI) for differentiating metabolic myopathies against non-metabolic myalgia. The optimum cutoff was taken as 1.30 (se = 69.6%, sp = 96.7%), with a corresponding diagnostic odd ratio of 66.3 (12.5-350.7, 95% CI). For a pretest probability of 15% in our tertiary reference center, the posttest probability for metabolic myopathy is 78.6% when MVC%pred:MP%pred ratio is above 1.3. In conclusion, the MVC%pred:MP%pred ratio is appropriate as a screening test to distinguish metabolic myopathies from non-metabolic myalgia.

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.

Estimating the age of p.(Phe508del) with family studies of geographically distinct European populations and the early spread of cystic fibrosis.

The high incidence of cystic fibrosis (CF) is due to the frequency of the c.1521_1523delCTT variant in the cystic fibrosis transmembrane conductance regulator (CFTR), but its age and origin are uncertain. This gap limits attempts to shed light on the presumed heterozygote selective advantage that accounts for the variant's high prevalence among Caucasian Europeans and Europe-derived populations. In addition, explaining the nature of heterozygosity to screened individuals with one c.1521_1523delCTT variant is challenging when families raise questions about these issues. To address this gap, we obtained DNA samples from 190 patients bearing c.1521_1523delCTT and their parents residing in geographically distinct European populations plus a Germany-derived population in the USA. We identified microsatellites spanning CFTR and reconstructed haplotypes at 10 loci to estimate the time/age of the most recent common ancestor (tMRCA) with the Estiage program. We found that the age estimates differ between northwestern populations, where the mean tMRCA values vary between 4600 and 4725 years, and the southeastern populations where c.1521_1523delCTT seems to have been introduced only about 1000 years ago. The tMRCA values of Central Europeans were intermediate. Thus, our data resolve a controversy by establishing an early Bronze Age origin of the c.1521_1523delCTT allele and demonstrating its likely spread from northwest to southeast during ancient migrations. Moreover, taking the archeological record into account, our results introduce a novel concept by suggesting that Bell Beaker folk were the probable migrating population responsible for the early dissemination of c.1521_1523delCTT in prehistoric Europe.