H3K36 dimethylation shapes the epigenetic interaction landscape by directing repressive chromatin modifications in embryonic stem cells.

Epigenetic modifications on the chromatin do not occur in isolation. Chromatin-associated proteins and their modification products form a highly interconnected network, and disturbing one component may rearrange the entire system. We see this increasingly clearly in epigenetically dysregulated cancers. It is important to understand the rules governing epigenetic interactions. Here, we use the mouse embryonic stem cell (mESC) model to describe in detail the relationships within the H3K27-H3K36-DNA methylation subnetwork. In particular, we focus on the major epigenetic reorganization caused by deletion of the histone 3 lysine 36 methyltransferase NSD1, which in mESCs deposits nearly all of the intergenic H3K36me2. Although disturbing the H3K27 and DNA methylation (DNAme) components also affects this network to a certain extent, the removal of H3K36me2 has the most drastic effect on the epigenetic landscape, resulting in full intergenic spread of H3K27me3 and a substantial decrease in DNAme. By profiling DNMT3A and CHH methylation (mCHH), we show that H3K36me2 loss upon Nsd1-KO leads to a massive redistribution of DNMT3A and mCHH away from intergenic regions and toward active gene bodies, suggesting that DNAme reduction is at least in part caused by redistribution of de novo methylation. Additionally, we show that pervasive acetylation of H3K27 is regulated by the interplay of H3K36 and H3K27 methylation. Our analysis highlights the importance of H3K36me2 as a major determinant of the developmental epigenome and provides a framework for further consolidating our knowledge of epigenetic networks.

The histone mark H3K36me2 recruits DNMT3A and shapes the intergenic DNA methylation landscape.

Enzymes that catalyse CpG methylation in DNA, including the DNA methyltransferases 1 (DNMT1), 3A (DNMT3A) and 3B (DNMT3B), are indispensable for mammalian tissue development and homeostasis1-4. They are also implicated in human developmental disorders and cancers5-8, supporting the critical role of DNA methylation in the specification and maintenance of cell fate. Previous studies have suggested that post-translational modifications of histones are involved in specifying patterns of DNA methyltransferase localization and DNA methylation at promoters and actively transcribed gene bodies9-11. However, the mechanisms that control the establishment and maintenance of intergenic DNA methylation remain poorly understood. Tatton-Brown-Rahman syndrome (TBRS) is a childhood overgrowth disorder that is defined by germline mutations in DNMT3A. TBRS shares clinical features with Sotos syndrome (which is caused by haploinsufficiency of NSD1, a histone methyltransferase that catalyses the dimethylation of histone H3 at K36 (H3K36me2)8,12,13), which suggests that there is a mechanistic link between these two diseases. Here we report that NSD1-mediated H3K36me2 is required for the recruitment of DNMT3A and maintenance of DNA methylation at intergenic regions. Genome-wide analysis shows that the binding and activity of DNMT3A colocalize with H3K36me2 at non-coding regions of euchromatin. Genetic ablation of Nsd1 and its paralogue Nsd2 in mouse cells results in a redistribution of DNMT3A to H3K36me3-modified gene bodies and a reduction in the methylation of intergenic DNA. Blood samples from patients with Sotos syndrome and NSD1-mutant tumours also exhibit hypomethylation of intergenic DNA. The PWWP domain of DNMT3A shows dual recognition of H3K36me2 and H3K36me3 in vitro, with a higher binding affinity towards H3K36me2 that is abrogated by TBRS-derived missense mutations. Together, our study reveals a trans-chromatin regulatory pathway that connects aberrant intergenic CpG methylation to human neoplastic and developmental overgrowth.

Sf3b4 regulates chromatin remodeler splicing and Hox expression.

SF3B proteins form a heptameric complex in the U2 small nuclear ribonucleoprotein, essential for pre-mRNA splicing. Heterozygous pathogenic variants in human SF3B4 are associated with head, face, limb, and vertebrae defects. Using the CRISPR/Cas9 system, we generated mice with constitutive heterozygous deletion of Sf3b4 and showed that mutant embryos have abnormal vertebral development. Vertebrae abnormalities were accompanied by changes in levels and expression pattern of Hox genes in the somites. RNA sequencing analysis of whole embryos and somites of Sf3b4 mutant and control litter mates revealed increased expression of other Sf3b4 genes. However, the mutants exhibited few differentially expressed genes and a large number of transcripts with differential splicing events (DSE), predominantly increased exon skipping and intron retention. Transcripts with increased DSE included several genes involved in chromatin remodeling that are known to regulate Hox expression. Our study confirms that Sf3b4 is required for normal vertebrae development and shows, for the first time, that like Sf3b1, Sf3b4 also regulates Hox expression. We propose that abnormal splicing of chromatin remodelers is primarily responsible for vertebral defects found in Sf3b4 heterozygous mutant embryos.

A report of two homozygous TERB1 protein-truncating variants in two unrelated women with primary infertility.

PURPOSE: To investigate the genetic etiology of patients with female infertility. METHODS: Whole Exome Sequencing was performed on genomic DNA extracted from the patient's blood. Exome data were filtered for damaging rare biallelic variants in genes with possible roles in reproduction. Sanger sequencing was used to validate the selected variants and segregate them in family members. RESULTS: A novel homozygous likely pathogenic variant, c.626G>A, p.Trp209*, was identified in the TERB1 gene of the patient. Additionally, we report a second homozygous pathogenic TERB1 variant, c.1703C>G, p.Ser568*, in an infertile woman whose azoospermic brother was previously described to be homozygous for her variant. CONCLUSIONS: Here, we report for the first time two homozygous likely pathogenic and pathogenic TERB1 variants, c.626G>A, p.Trp209* and c.1703C>G, p.Ser568*, respectively, in two unrelated women with primary infertility. TERB1 is known to play an essential role in homologous chromosome movement, synapsis, and recombination during the meiotic prophase I and has an established role in male infertility in humans. Our data add TERB1 to the shortlist of Meiosis I genes associated with human infertility in both sexes.

Mutation in Eftud2 causes craniofacial defects in mice via mis-splicing of Mdm2 and increased P53.

EFTUD2 is mutated in patients with mandibulofacial dysostosis with microcephaly (MFDM). We generated a mutant mouse line with conditional mutation in Eftud2 and used Wnt1-Cre2 to delete it in neural crest cells. Homozygous deletion of Eftud2 causes brain and craniofacial malformations, affecting the same precursors as in MFDM patients. RNAseq analysis of embryonic heads revealed a significant increase in exon skipping and increased levels of an alternatively spliced Mdm2 transcript lacking exon 3. Exon skipping in Mdm2 was also increased in O9-1 mouse neural crest cells after siRNA knock-down of Eftud2 and in MFDM patient cells. Moreover, we found increased nuclear P53, higher expression of P53-target genes and increased cell death. Finally, overactivation of the P53 pathway in Eftud2 knockdown cells was attenuated by overexpression of non-spliced Mdm2, and craniofacial development was improved when Eftud2-mutant embryos were treated with Pifithrin-α, an inhibitor of P53. Thus, our work indicates that the P53-pathway can be targeted to prevent craniofacial abnormalities and shows a previously unknown role for alternative splicing of Mdm2 in the etiology of MFDM.

Does the adenosine deaminase (ADA) gene confer risk of sleepwalking?

Sleepwalking is a common non-rapid eye movement (NREM) parasomnia and a significant cause of sleep-related injuries. While evidence suggest that the occurrence of this condition is partly determined by genetic factors, its pattern of inheritance remains unclear, and few molecular studies have been conducted. One promising candidate is the adenosine deaminase (ADA) gene. Adenosine and the ADA enzyme play an important role in the homeostatic regulation of NREM sleep. In a single sleepwalking family, genome-wide analysis identified a locus on chromosome 20, where ADA lies. In this study, we examined if variants in the ADA gene were associated with sleepwalking. In total, 251 sleepwalking patients were clinically assessed, and DNA samples were compared to those from 94 unaffected controls. Next-generation sequencing of the whole ADA gene was performed. Bio-informatic analysis enabled the identification of variants and assessed variants enrichment in our cohort compared to controls. We detected 25 different coding and non-coding variants, of which 22 were found among sleepwalkers. None were enriched in the sleepwalking population. However, many missense variants were predicted as likely pathogenic by at least two in silico prediction algorithms. This study involves the largest sleepwalking cohort in which the role of a susceptibility gene was investigated. Our results did not reveal an association between ADA gene and sleepwalking, thus ruling out the possibility of ADA as a major genetic factor for this condition. Future work is needed to identify susceptibility genes.

TDP-43 stabilizes G3BP1 mRNA: relevance to amyotrophic lateral sclerosis/frontotemporal dementia.

TDP-43 nuclear depletion and concurrent cytoplasmic accumulation in vulnerable neurons is a hallmark feature of progressive neurodegenerative proteinopathies such as amyotrophic lateral sclerosis (ALS) and frontotemporal dementia (FTD). Cellular stress signalling and stress granule dynamics are now recognized to play a role in ALS/FTD pathogenesis. Defective stress granule assembly is associated with increased cellular vulnerability and death. Ras-GAP SH3-domain-binding protein 1 (G3BP1) is a critical stress granule assembly factor. Here, we define that TDP-43 stabilizes G3BP1 transcripts via direct binding of a highly conserved cis regulatory element within the 3' untranslated region. Moreover, we show in vitro and in vivo that nuclear TDP-43 depletion is sufficient to reduce G3BP1 protein levels. Finally, we establish that G3BP1 transcripts are reduced in ALS/FTD patient neurons bearing TDP-43 cytoplasmic inclusions/nuclear depletion. Thus, our data indicate that, in ALS/FTD, there is a compromised stress granule response in disease-affected neurons due to impaired G3BP1 mRNA stability caused by TDP-43 nuclear depletion. These data implicate TDP-43 and G3BP1 loss of function as contributors to disease.

Whole-exome sequencing reveals novel vacuolar ATPase genes' variants and variants in genes involved in lysosomal biology and autophagosomal formation in oral granular cell tumors.

BACKGROUND: Granular cell tumors (GCTs) are rare neuroectodermal soft tissue neoplasms that mainly affect the skin of the upper limbs and trunks and the oral cavity. GCTs are derived from Schwann cells and, ultrastructurally, their intracytoplasmic granules are considered autophagosomes or autophagolysosomes and are consistent with myelin accumulation. METHODS: In this study, a convenience set of 22 formalin-fixed, paraffin-embedded samples of oral GCTs, all but one sample located at the tongue, was screened for mutations by whole-exome (WES) or targeted sequencing. RESULTS: WES revealed two novel variants in genes of the vacuolar ATPase (V-ATPase) complex: ATP6AP1 frameshift c.746_749del, leading to p.P249Hfs*4, and ATP6V1A non-synonymous c.G868A, leading to p.D290N. Each of these mutations occurred in one case. With regard to the samples that were wild type for these V-ATPase variants, at least two samples presented variants in genes that are part of endosomal/lysosomal/autophagosomal networks including ABCA8, ABCC6, AGAP3, ATG9A, CTSB, DNAJC13, GALC, NPC1, SLC15A3, SLC31A2, and TMEM104. CONCLUSION: Although the mechanisms involved in oral GCT initiation and progression remain unclear, our results suggest that oral GCTs have V-ATPase variants similarly to GCTs from other tissues/organs, and additionally show variants in lysosomes/endosomes/autophagosomal genes.

Bmp signaling maintains a mesoderm progenitor cell state in the mouse tailbud.

Caudal somites are generated from a pool of progenitor cells located in the tailbud region. These progenitor cells form the presomitic mesoderm that gradually differentiates into somites under the action of the segmentation clock. The signals responsible for tailbud mesoderm progenitor pool maintenance during axial elongation are still elusive. Here, we show that Bmp signaling is sufficient to activate the entire mesoderm progenitor gene signature in primary cultures of caudal mesoderm cells. Bmp signaling acts through the key regulatory genes brachyury (T) and Nkx1-2 and contributes to the activation of several other regulators of the mesoderm progenitor gene network. In the absence of Bmp signaling, tailbud mesoderm progenitor cells acquire aberrant gene expression signatures of the heart, blood, muscle and skeletal embryonic lineages. Treatment of embryos with the Bmp inhibitor noggin confirmed the requirement for Bmp signaling for normal T expression and the prevention of abnormal lineage marker activation. Together, these results identify Bmp signaling as a non-cell-autonomous signal necessary for mesoderm progenitor cell homeostasis.

Early infantile epileptic encephalopathy due to biallelic pathogenic variants in PIGQ: Report of seven new subjects and review of the literature.

We investigated seven children from six families to expand the phenotypic spectrum associated with an early infantile epileptic encephalopathy caused by biallelic pathogenic variants in the phosphatidylinositol glycan anchor biosynthesis class Q (PIGQ) gene. The affected children were all identified by clinical or research exome sequencing. Clinical data, including EEGs and MRIs, was comprehensively reviewed and flow cytometry and transfection experiments were performed to investigate PIGQ function. Pathogenic biallelic PIGQ variants were associated with increased mortality. Epileptic seizures, axial hypotonia, developmental delay and multiple congenital anomalies were consistently observed. Seizure onset occurred between 2.5 months and 7 months of age and varied from treatable seizures to recurrent episodes of status epilepticus. Gastrointestinal issues were common and severe, two affected individuals had midgut volvulus requiring surgical correction. Cardiac anomalies including arrythmias were observed. Flow cytometry using granulocytes and fibroblasts from affected individuals showed reduced expression of glycosylphosphatidylinositol (GPI)-anchored proteins. Transfection of wildtype PIGQ cDNA into patient fibroblasts rescued this phenotype. We expand the phenotypic spectrum of PIGQ-related disease and provide the first functional evidence in human cells of defective GPI-anchoring due to pathogenic variants in PIGQ.

Recessive mutation in CD2AP causes focal segmental glomerulosclerosis in humans and mice.

Although sequence variants in CD2-associated protein (CD2AP) have been identified in patients with focal segmental glomerulosclerosis (FSGS), definitive proof of causality in human disease is meager. By whole-exome sequencing, we identified a homozygous frame-shift mutation in CD2AP (p.S198fs) in three siblings born of consanguineous parents who developed childhood-onset FSGS and end stage renal disease. When the same frameshift mutation was introduced in mice by gene editing, the mice developed FSGS and kidney failure. These results provide conclusive evidence that homozygous mutation of CD2AP causes FSGS in humans.

SHORT syndrome due to a novel de novo mutation in PRKCE (Protein Kinase Cɛ) impairing TORC2-dependent AKT activation.

SHORT syndrome is a rare, recognizable syndrome resulting from heterozygous mutations in PIK3R1 encoding a regulatory subunit of phosphoinositide-3-kinase (PI3K). The condition is characterized by short stature, intrauterine growth restriction, lipoatrophy and a facial gestalt involving a triangular face, deep set eyes, low hanging columella and small chin. PIK3R1 mutations in SHORT syndrome result in reduced signaling through the PI3K-AKT-mTOR pathway. We performed whole exome sequencing for an individual with clinical features of SHORT syndrome but negative for PIK3R1 mutation and her parents. A rare de novo variant in PRKCE was identified. The gene encodes PKCε and, as such, the AKT-mTOR pathway function was assessed using phospho-specific antibodies with patient lymphoblasts and following ectopic expression of the mutant in HEK293 cells. Kinase analysis showed that the variant resulted in a partial loss-of-function. Whilst interaction with PDK1 and the mTORC2 complex component SIN1 was preserved in the mutant PKCε, it bound to SIN1 with a higher affinity than wild-type PKCε and the dynamics of mTORC2-dependent priming of mutant PKCε was altered. Further, mutant PKCε caused impaired mTORC2-dependent pAKT-S473 following rapamycin treatment. Reduced pFOXO1-S256 and pS6-S240/244 levels were also observed in the patient LCLs. To date, mutations in PIK3R1 causing impaired PI3K-dependent AKT activation are the only known cause of SHORT syndrome. We identify a SHORT syndrome child with a novel partial loss-of-function defect in PKCε. This variant causes impaired AKT activation via compromised mTORC2 complex function.

Neu-Laxova syndrome presenting prenatally with increased nuchal translucency and cystic hygroma: The utility of exome sequencing in deciphering the diagnosis.

Neu-Laxova syndrome (NLS) is a lethal autosomal recessive microcephaly syndrome associated with intrauterine growth restriction (IUGR) and multiple congenital anomalies. Clinical features include central nervous system malformations, joint contractures, ichthyosis, edema, and dysmorphic facial features. Biallelic pathogenic variants in either the PHGDH or PSAT1 genes have been shown to cause NLS. Using exome sequencing, we aimed to identify the underlying genetic diagnosis in three fetuses (from one family) with prenatal skin edema, severe IUGR, micrognathia, renal anomalies, and arthrogryposis and identified a homozygous c.1A>C (p.Met1?, NM_006623.3) variant in the PHGDH gene. Loss of the translation start codon is a novel genetic mechanism for the development of NLS. Prenatal diagnosis of NLS is challenging and few reports describe the fetal pathology. Fetal neuropathologic examination revealed: delayed brain development, congenital agenesis of the corticospinal tracts, and hypoplasia of the hippocampus, cerebellum and brainstem. Each pregnancy also showed increased nuchal translucency (NT) or cystic hygroma. While NLS is rare, it may be a cause of recurrent increased NT/cystic hygroma. This finding provides further support that cystic hygroma has many different genetic causes and that exome sequencing may shed light on the underlying genetic diagnoses in this group of prenatal patients.

Recessive mutations in the kinase ZAK cause a congenital myopathy with fibre type disproportion.

Congenital myopathies define a heterogeneous group of neuromuscular diseases with neonatal or childhood hypotonia and muscle weakness. The genetic cause is still unknown in many patients, precluding genetic counselling and better understanding of the physiopathology. To identify novel genetic causes of congenital myopathies, exome sequencing was performed in three consanguineous families. We identified two homozygous frameshift mutations and a homozygous nonsense mutation in the mitogen-activated protein triple kinase ZAK. In total, six affected patients carry these mutations. Reverse transcription polymerase chain reaction and transcriptome analyses suggested nonsense mRNA decay as a main impact of mutations. The patients demonstrated a generalized slowly progressive muscle weakness accompanied by decreased vital capacities. A combination of proximal contractures with distal joint hyperlaxity is a distinct feature in one family. The low endurance and compound muscle action potential amplitude were strongly ameliorated on treatment with anticholinesterase inhibitor in another patient. Common histopathological features encompassed fibre size variation, predominance of type 1 fibre and centralized nuclei. A peculiar subsarcolemmal accumulation of mitochondria pointing towards the centre of the fibre was a novel histological hallmark in one family. These findings will improve the molecular diagnosis of congenital myopathies and implicate the mitogen-activated protein kinase (MAPK) signalling as a novel pathway altered in these rare myopathies.

Autosomal dominant cutis laxa with progeroid features due to a novel, de novo mutation in ALDH18A1.

De novo dominant mutations in the aldehyde dehydrogenase 18 family member A1 (ALDH18A1) gene have recently been shown to cause autosomal dominant cutis laxa with progeroid features (MIM 616603). To date, all de novo dominant mutations have been found in a single highly conserved amino acid residue at position p.Arg138. We report an 8-year-old male with a clinical diagnosis of autosomal dominant cutis laxa (ADCL) with progeroid features and a novel de novo missense mutation in ALDH18A1 (NM_002860.3: c.377G>A (p.Arg126His)). This is the first report of an individual with ALDH18A1-ADCL due to a substitution at a residue other than p.Arg138. Knowledge of the complete spectrum of dominant-acting mutations that cause this rare syndrome will have implications for molecular diagnosis and genetic counselling of these families.

H3.1 K36M mutation in a congenital-onset soft tissue neoplasm.

We describe a patient who presented with a congenital soft tissue lesion initially diagnosed as infantile fibromatosis at 15 days of age. Unusually, the mass demonstrated malignant progression leading to death at 20 months of age. Biological progression to malignancy is not known to occur in fibromatosis, and fibrosarcoma is not known to progress from a benign lesion. Whole-exome sequencing of the tumor identified a driver mutation in histone H3.1 at lysine (K)36. Our findings support the link between oncohistones and infantile soft tissue tumors and provide additional evidence for the oncogenic effects of p.K36M in H3 variants.

Germline and somatic FGFR1 abnormalities in dysembryoplastic neuroepithelial tumors.

Dysembryoplastic neuroepithelial tumor (DNET) is a benign brain tumor associated with intractable drug-resistant epilepsy. In order to identify underlying genetic alterations and molecular mechanisms, we examined three family members affected by multinodular DNETs as well as 100 sporadic tumors from 96 patients, which had been referred to us as DNETs. We performed whole-exome sequencing on 46 tumors and targeted sequencing for hotspot FGFR1 mutations and BRAF p.V600E was used on the remaining samples. FISH, copy number variation assays and Sanger sequencing were used to validate the findings. By whole-exome sequencing of the familial cases, we identified a novel germline FGFR1 mutation, p.R661P. Somatic activating FGFR1 mutations (p.N546K or p.K656E) were observed in the tumor samples and further evidence for functional relevance was obtained by in silico modeling. The FGFR1 p.K656E mutation was confirmed to be in cis with the germline p.R661P variant. In 43 sporadic cases, in which the diagnosis of DNET could be confirmed on central blinded neuropathology review, FGFR1 alterations were also frequent and mainly comprised intragenic tyrosine kinase FGFR1 duplication and multiple mutants in cis (25/43; 58.1 %) while BRAF p.V600E alterations were absent (0/43). In contrast, in 53 cases, in which the diagnosis of DNET was not confirmed, FGFR1 alterations were less common (10/53; 19 %; p < 0.0001) and hotspot BRAF p.V600E (12/53; 22.6 %) (p < 0.001) prevailed. We observed overexpression of phospho-ERK in FGFR1 p.R661P and p.N546K mutant expressing HEK293 cells as well as FGFR1 mutated tumor samples, supporting enhanced MAP kinase pathway activation under these conditions. In conclusion, constitutional and somatic FGFR1 alterations and MAP kinase pathway activation are key events in the pathogenesis of DNET. These findings point the way towards existing targeted therapies.

Next-generation sequencing for diagnosis of rare diseases in the neonatal intensive care unit.

BACKGROUND: Rare diseases often present in the first days and weeks of life and may require complex management in the setting of a neonatal intensive care unit (NICU). Exhaustive consultations and traditional genetic or metabolic investigations are costly and often fail to arrive at a final diagnosis when no recognizable syndrome is suspected. For this pilot project, we assessed the feasibility of next-generation sequencing as a tool to improve the diagnosis of rare diseases in newborns in the NICU. METHODS: We retrospectively identified and prospectively recruited newborns and infants admitted to the NICU of the Children's Hospital of Eastern Ontario and the Ottawa Hospital, General Campus, who had been referred to the medical genetics or metabolics inpatient consult service and had features suggesting an underlying genetic or metabolic condition. DNA from the newborns and parents was enriched for a panel of clinically relevant genes and sequenced on a MiSeq sequencing platform (Illumina Inc.). The data were interpreted with a standard informatics pipeline and reported to care providers, who assessed the importance of genotype-phenotype correlations. RESULTS: Of 20 newborns studied, 8 received a diagnosis on the basis of next-generation sequencing (diagnostic rate 40%). The diagnoses were renal tubular dysgenesis, SCN1A-related encephalopathy syndrome, myotubular myopathy, FTO deficiency syndrome, cranioectodermal dysplasia, congenital myasthenic syndrome, autosomal dominant intellectual disability syndrome type 7 and Denys-Drash syndrome. INTERPRETATION: This pilot study highlighted the potential of next-generation sequencing to deliver molecular diagnoses rapidly with a high success rate. With broader use, this approach has the potential to alter health care delivery in the NICU.