Histone methylation antagonism drives tumor immune evasion in squamous cell carcinomas.

How cancer-associated chromatin abnormalities shape tumor-immune interaction remains incompletely understood. Recent studies have linked DNA hypomethylation and de-repression of retrotransposons to anti-tumor immunity through the induction of interferon response. Here, we report that inactivation of the histone H3K36 methyltransferase NSD1, which is frequently found in squamous cell carcinomas (SCCs) and induces DNA hypomethylation, unexpectedly results in diminished tumor immune infiltration. In syngeneic and genetically engineered mouse models of head and neck SCCs, NSD1-deficient tumors exhibit immune exclusion and reduced interferon response despite high retrotransposon expression. Mechanistically, NSD1 loss results in silencing of innate immunity genes, including the type III interferon receptor IFNLR1, through depletion of H3K36 di-methylation (H3K36me2) and gain of H3K27 tri-methylation (H3K27me3). Inhibition of EZH2 restores immune infiltration and impairs the growth of Nsd1-mutant tumors. Thus, our work uncovers a druggable chromatin cross talk that regulates the viral mimicry response and enables immune evasion of DNA hypomethylated tumors.

USP15 regulates type I interferon response and is required for pathogenesis of neuroinflammation.

Genes and pathways in which inactivation dampens tissue inflammation present new opportunities for understanding the pathogenesis of common human inflammatory diseases, including inflammatory bowel disease, rheumatoid arthritis and multiple sclerosis. We identified a mutation in the gene encoding the deubiquitination enzyme USP15 (Usp15L749R) that protected mice against both experimental cerebral malaria (ECM) induced by Plasmodium berghei and experimental autoimmune encephalomyelitis (EAE). Combining immunophenotyping and RNA sequencing in brain (ECM) and spinal cord (EAE) revealed that Usp15L749R-associated resistance to neuroinflammation was linked to dampened type I interferon responses in situ. In hematopoietic cells and in resident brain cells, USP15 was coexpressed with, and functionally acted together with the E3 ubiquitin ligase TRIM25 to positively regulate type I interferon responses and to promote pathogenesis during neuroinflammation. The USP15-TRIM25 dyad might be a potential target for intervention in acute or chronic states of neuroinflammation.

Nucleome programming is required for the foundation of totipotency in mammalian germline development.

Germ cells are unique in engendering totipotency, yet the mechanisms underlying this capacity remain elusive. Here, we perform comprehensive and in-depth nucleome analysis of mouse germ-cell development in vitro, encompassing pluripotent precursors, primordial germ cells (PGCs) before and after epigenetic reprogramming, and spermatogonia/spermatogonial stem cells (SSCs). Although epigenetic reprogramming, including genome-wide DNA de-methylation, creates broadly open chromatin with abundant enhancer-like signatures, the augmented chromatin insulation safeguards transcriptional fidelity. These insulatory constraints are then erased en masse for spermatogonial development. Notably, despite distinguishing epigenetic programming, including global DNA re-methylation, the PGCs-to-spermatogonia/SSCs development entails further euchromatization. This accompanies substantial erasure of lamina-associated domains, generating spermatogonia/SSCs with a minimal peripheral attachment of chromatin except for pericentromeres-an architecture conserved in primates. Accordingly, faulty nucleome maturation, including persistent insulation and improper euchromatization, leads to impaired spermatogenic potential. Given that PGCs after epigenetic reprogramming serve as oogenic progenitors as well, our findings elucidate a principle for the nucleome programming that creates gametogenic progenitors in both sexes, defining a basis for nuclear totipotency.

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.

Depletion of H3K36me2 recapitulates epigenomic and phenotypic changes induced by the H3.3K36M oncohistone mutation.

Hotspot histone H3 mutations have emerged as drivers of oncogenesis in cancers of multiple lineages. Specifically, H3 lysine 36 to methionine (H3K36M) mutations are recurrently identified in chondroblastomas, undifferentiated sarcomas, and head and neck cancers. While the mutation reduces global levels of both H3K36 dimethylation (H3K36me2) and trimethylation (H3K36me3) by dominantly inhibiting their respective specific methyltransferases, the relative contribution of these methylation states to the chromatin and phenotypic changes associated with H3K36M remains unclear. Here, we specifically deplete H3K36me2 or H3K36me3 in mesenchymal cells, using CRISPR-Cas9 to separately knock out the corresponding methyltransferases NSD1/2 or SETD2. By profiling and comparing the epigenomic and transcriptomic landscapes of these cells with cells expressing the H3.3K36M oncohistone, we find that the loss of H3K36me2 could largely recapitulate H3.3K36M's effect on redistribution of H3K27 trimethylation (H3K27me3) and gene expression. Consistently, knockout of Nsd1/2, but not Setd2, phenocopies the differentiation blockade and hypersensitivity to the DNA-hypomethylating agent induced by H3K36M. Together, our results support a functional divergence between H3K36me2 and H3K36me3 and their nonredundant roles in H3K36M-driven oncogenesis.

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.

Infantile onset encephalomyopathy, retinopathy, optic atrophy, and mitochondrial DNA depletion associated with a novel pathogenic DHX16 variant.

We studied a patient with mitochondrial DNA depletion in skeletal muscle and a multiorgan phenotype, including fatal encephalomyopathy, retinopathy, optic atrophy, and sensorineural hearing loss. Instead of pathogenic variants in the mitochondrial maintenance genes, we identified previously unpublished variant in DHX16 gene, a de novo heterozygous c.1360C>T (p. Arg454Trp). Variants in DHX16 encoding for DEAH-box RNA helicase have previously been reported only in five patients with a phenotype called as neuromuscular oculoauditory syndrome including developmental delay, neuromuscular symptoms, and ocular or auditory defects with or without seizures. We performed functional studies on patient-derived fibroblasts and skeletal muscle revealing, that the DHX16 expression was decreased. Clinical features together with functional data suggest, that our patient's disease is associated with a novel pathogenic DHX16 variant, and mtDNA depletion could be a secondary manifestation of the disease.

Detection and genomic analysis of BRAF fusions in Juvenile Pilocytic Astrocytoma through the combination and integration of multi-omic data.

BACKGROUND: Juvenile Pilocytic Astrocytomas (JPAs) are one of the most common pediatric brain tumors, and they are driven by aberrant activation of the mitogen-activated protein kinase (MAPK) signaling pathway. RAF-fusions are the most common genetic alterations identified in JPAs, with the prototypical KIAA1549-BRAF fusion leading to loss of BRAF's auto-inhibitory domain and subsequent constitutive kinase activation. JPAs are highly vascular and show pervasive immune infiltration, which can lead to low tumor cell purity in clinical samples. This can result in gene fusions that are difficult to detect with conventional omics approaches including RNA-Seq. METHODS: To this effect, we applied RNA-Seq as well as linked-read whole-genome sequencing and in situ Hi-C as new approaches to detect and characterize low-frequency gene fusions at the genomic, transcriptomic and spatial level. RESULTS: Integration of these datasets allowed the identification and detailed characterization of two novel BRAF fusion partners, PTPRZ1 and TOP2B, in addition to the canonical fusion with partner KIAA1549. Additionally, our Hi-C datasets enabled investigations of 3D genome architecture in JPAs which showed a high level of correlation in 3D compartment annotations between JPAs compared to other pediatric tumors, and high similarity to normal adult astrocytes. We detected interactions between BRAF and its fusion partners exclusively in tumor samples containing BRAF fusions. CONCLUSIONS: We demonstrate the power of integrating multi-omic datasets to identify low frequency fusions and characterize the JPA genome at high resolution. We suggest that linked-reads and Hi-C could be used in clinic for the detection and characterization of JPAs.

A region-based gene association study combined with a leave-one-out sensitivity analysis identifies SMG1 as a pancreatic cancer susceptibility gene.

Pancreatic adenocarcinoma (PC) is a lethal malignancy that is familial or associated with genetic syndromes in 10% of cases. Gene-based surveillance strategies for at-risk individuals may improve clinical outcomes. However, familial PC (FPC) is plagued by genetic heterogeneity and the genetic basis for the majority of FPC remains elusive, hampering the development of gene-based surveillance programs. The study was powered to identify genes with a cumulative pathogenic variant prevalence of at least 3%, which includes the most prevalent PC susceptibility gene, BRCA2. Since the majority of known PC susceptibility genes are involved in DNA repair, we focused on genes implicated in these pathways. We performed a region-based association study using the Mixed-Effects Score Test, followed by leave-one-out characterization of PC-associated gene regions and variants to identify the genes and variants driving risk associations. We evaluated 398 cases from two case series and 987 controls without a personal history of cancer. The first case series consisted of 109 patients with either FPC (n = 101) or PC at ≤50 years of age (n = 8). The second case series was composed of 289 unselected PC cases. We validated this discovery strategy by identifying known pathogenic BRCA2 variants, and also identified SMG1, encoding a serine/threonine protein kinase, to be significantly associated with PC following correction for multiple testing (p = 3.22x10-7). The SMG1 association was validated in a second independent series of 532 FPC cases and 753 controls (p<0.0062, OR = 1.88, 95%CI 1.17-3.03). We showed segregation of the c.4249A>G SMG1 variant in 3 affected relatives in a FPC kindred, and we found c.103G>A to be a recurrent SMG1 variant associating with PC in both the discovery and validation series. These results suggest that SMG1 is a novel PC susceptibility gene, and we identified specific SMG1 gene variants associated with PC risk.

ClinPred: Prediction Tool to Identify Disease-Relevant Nonsynonymous Single-Nucleotide Variants.

Advances in high-throughput DNA sequencing have revolutionized the discovery of variants in the human genome; however, interpreting the phenotypic effects of those variants is still a challenge. While several computational approaches to predict variant impact are available, their accuracy is limited and further improvement is needed. Here, we introduce ClinPred, an efficient tool for identifying disease-relevant nonsynonymous variants. Our predictor incorporates two machine learning algorithms that use existing pathogenicity scores and, notably, benefits from inclusion of normal population allele frequency from the gnomAD database as an input feature. Another major strength of our approach is the use of ClinVar-a rapidly growing database that allows selection of confidently annotated disease-causing variants-as a training set. Compared to other methods, ClinPred showed superior accuracy for predicting pathogenicity, achieving the highest area under the curve (AUC) score and increasing both the specificity and sensitivity in different test datasets. It also obtained the best performance according to various other metrics. Moreover, ClinPred performance remained robust with respect to disease type (cancer or rare disease) and mechanism (gain or loss of function). Importantly, we observed that adding allele frequency as a predictive feature-as opposed to setting fixed allele frequency cutoffs-boosts the performance of prediction. We provide pre-computed ClinPred scores for all possible human missense variants in the exome to facilitate its use by the community.

Causative Mutations and Mechanism of Androgenetic Hydatidiform Moles.

Androgenetic complete hydatidiform moles are human pregnancies with no embryos and affect 1 in every 1,400 pregnancies. They have mostly androgenetic monospermic genomes with all the chromosomes originating from a haploid sperm and no maternal chromosomes. Androgenetic complete hydatidiform moles were described in 1977, but how they occur has remained an open question. We identified bi-allelic deleterious mutations in MEI1, TOP6BL/C11orf80, and REC114, with roles in meiotic double-strand breaks formation in women with recurrent androgenetic complete hydatidiform moles. We investigated the occurrence of androgenesis in Mei1-deficient female mice and discovered that 8% of their oocytes lose all their chromosomes by extruding them with the spindles into the first polar body. We demonstrate that Mei1-/- oocytes are capable of fertilization and 5% produce androgenetic zygotes. Thus, we uncover a meiotic abnormality in mammals and a mechanism for the genesis of androgenetic zygotes that is the extrusion of all maternal chromosomes and their spindles into the first polar body.

Novel pathogenic variants in NLRP7, NLRP5, and PADI6 in patients with recurrent hydatidiform moles and reproductive failure.

Recurrent hydatidiform moles (RHMs) are human pregnancies with abnormal embryonic development and hyperproliferating trophoblast. Biallelic mutations in NLRP7 and KHDC3L, members of the subcortical maternal complex (SCMC), explain the etiology of RHMs in only 60% of patients. Here we report the identification of seven functional variants in a recessive state in three SCMC members, five in NLRP7, one in NLRP5, and one in PADI6. In NLRP5, we report the first patient with RHMs and biallelic mutations. In PADI6, the patient had four molar pregnancies, two of which had fetuses with various abnormalities including placental mesenchymal dysplasia and intra-uterine growth restriction, which are features of Beckwith-Wiedemann syndrome and Silver Russell syndrome, respectively. Our findings corroborate recent studies and highlight the common oocyte origin of all these conditions and the continuous spectrum of abnormalities associated with deficiencies in the SCMC genes.

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.

Implementation of the Configuration Structure of an Integrated Computational Core of a Pulsed NQR Sensor Based on FPGA.

This paper presents a method for implementing the configuration structure of an integrated computational core of a pulsed nuclear quadrupole resonance (NQR) sensor based on a field-programmable gate array (FPGA), which comprises the following modules: a three-channel direct digital synthesizer (DDS), a pulse sequence shaper and a software-defined radio. Experimental studies carried out using the in-circuit analyzer SignalTap Logic Analyzer have confirmed the reliability of the correct and stable operation of the functional modules of the configuration structure at all stages of signal transformations, starting from the formation of the envelope of the excitation pulses and ending with the obtainment of low-frequency quadrature signals at the outlet of the compensating filters. The time and frequency dependences of the amplitude of the output signals generated using the DDS based on a 48 bit phase accumulator are investigated. This development can be used when creating pulsed coherent NQR sensors in the frequency range of 1 MHz-50 MHz.

ZBTB7B (ThPOK) Is Required for Pathogenesis of Cerebral Malaria and Protection against Pulmonary Tuberculosis.

We used a genome-wide screen in N-ethyl-N-nitrosourea (ENU)-mutagenized mice to identify genes in which recessive loss-of-function mutations protect against pathological neuroinflammation. We identified an R367Q mutation in the ZBTB7B (ThPOK) protein in which homozygosity causes protection against experimental cerebral malaria (ECM) caused by infection with Plasmodium berghei ANKA. Zbtb7bR367Q homozygous mice show a defect in the lymphoid compartment expressed as severe reduction in the number of single-positive CD4 T cells in the thymus and in the periphery, reduced brain infiltration of proinflammatory leukocytes in P. berghei ANKA-infected mice, and reduced production of proinflammatory cytokines by primary T cells ex vivo and in vivo Dampening of proinflammatory immune responses in Zbtb7bR367Q mice is concomitant to increased susceptibility to infection with avirulent (Mycobacterium bovis BCG) and virulent (Mycobacterium tuberculosis H37Rv) mycobacteria. The R367Q mutation maps to the first DNA-binding zinc finger domain of ThPOK and causes loss of base contact by R367 in the major groove of the DNA, which is predicted to impair DNA binding. Global immunoprecipitation of ThPOK-containing chromatin complexes coupled to DNA sequencing (ChIP-seq) identified transcriptional networks and candidate genes likely to play key roles in CD4+ CD8+ T cell development and in the expression of lineage-specific functions of these cells. This study highlights ThPOK as a global regulator of immune function in which alterations may affect normal responses to infectious and inflammatory stimuli.

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.