PHF3 regulates neuronal gene expression through the Pol II CTD reader domain SPOC

Uncovering the dynamics and consequences of RNA isoform changes during neuronal differentiation

Static gene expression programs have been extensively characterized in stem cells and mature human cells. However, the dynamics of RNA isoform changes upon cell-state-transitions during cell differentiation, the determinants and functional consequences have largely remained unclear. Here, we established an improved model for human neurogenesis in vitro that is amenable for systems-wide analyses of gene expression. Our multi-omics analysis reveals that the pronounced alterations in cell morphology correlate strongly with widespread changes in RNA isoform expression. Our approach identifies thousands of new RNA isoforms that are expressed at distinct differentiation stages. RNA isoforms mainly arise from exon skipping and the alternative usage of transcription start and polyadenylation sites during human neurogenesis. The transcript isoform changes can remodel the identity and functions of protein isoforms. Finally, our study identifies a set of RNA binding proteins as a potential determinant of differentiation stage-specific global isoform changes. This work supports the view of regulated isoform changes that underlie state-transitions during neurogenesis.

Exploring the Early Molecular Pathogenesis of Osteoarthritis Using Differential Network Analysis of Human Synovial Fluid

The molecular mechanisms that drive the onset and development of osteoarthritis (OA) remain largely unknown. In this exploratory study, we used a proteomic platform (SOMAscan assay) to measure the relative abundance of more than 6000 proteins in synovial fluid (SF) from knees of human donors with healthy or mildly degenerated tissues, and knees with late-stage OA from patients undergoing knee replacement surgery. Using a linear mixed effects model, we estimated the differential abundance of 6251 proteins between the three groups. We found 583 proteins upregulated in the late-stage OA, including MMP1, collagenase 3 and interleukin-6. Further, we selected 760 proteins (800 aptamers) based on absolute fold changes between the healthy and mild degeneration groups. To those, we applied Gaussian Graphical Models (GGMs) to analyze the conditional dependence of proteins and to identify key proteins and subnetworks involved in early OA pathogenesis. After regularization and stability selection, we identified 102 proteins involved in GGM networks. Notably, network complexity was lost in the protein graph for mild degeneration when compared to controls, suggesting a disruption in the regular protein interplay. Furthermore, among our main findings were several downregulated (in mild degeneration versus healthy) proteins with unique interactions in the healthy group, one of which, SLCO5A1, has not previously been associated with OA. Our results suggest that this protein is important for healthy joint function. Further, our data suggests that SF proteomics, combined with GGMs, can reveal novel insights into the molecular pathogenesis and identification of biomarker candidates for early-stage OA.

Genomic analysis of 116 autism families strengthens known risk genes and highlights promising candidates

Autism spectrum disorder (ASD) is a complex neurodevelopmental condition with a strong genetic component in which rare variants contribute significantly to risk. We performed whole genome and/or exome sequencing (WGS and WES) and SNP-array analysis to identify both rare sequence and copy number variants (SNVs and CNVs) in 435 individuals from 116 ASD families. We identified 37 rare potentially damaging de novo SNVs (pdSNVs) in the cases (n = 144). Interestingly, two of them (one stop-gain and one missense variant) occurred in the same gene, BRSK2. Moreover, the identification of 8 severe de novo pdSNVs in genes not previously implicated in ASD (AGPAT3, IRX5, MGAT5B, RAB8B, RAP1A, RASAL2, SLC9A1, YME1L1) highlighted promising candidates. Potentially damaging CNVs (pdCNVs) provided support to the involvement of inherited variants in PHF3, NEGR1, TIAM1 and HOMER1 in neurodevelopmental disorders (NDD), although mostly acting as susceptibility factors with incomplete penetrance. Interpretation of identified pdSNVs/pdCNVs according to the ACMG guidelines led to a molecular diagnosis in 19/144 cases, although this figure represents a lower limit and is expected to increase thanks to further clarification of the role of likely pathogenic variants in ASD/NDD candidate genes not yet established. In conclusion, our study highlights promising ASD candidate genes and contributes to characterize the allelic diversity, mode of inheritance and phenotypic impact of de novo and inherited risk variants in ASD/NDD genes.

C19ORF84 connects piRNA and DNA methylation machineries to defend the mammalian germ line

In the male mouse germ line, PIWI-interacting RNAs (piRNAs), bound by the PIWI protein MIWI2 (PIWIL4), guide DNA methylation of young active transposons through SPOCD1. However, the underlying mechanisms of SPOCD1-mediated piRNA-directed transposon methylation and whether this pathway functions to protect the human germ line remain unknown. We identified loss-of-function variants in human SPOCD1 that cause defective transposon silencing and male infertility. Through the analysis of these pathogenic alleles, we discovered that the uncharacterized protein C19ORF84 interacts with SPOCD1. DNMT3C, the DNA methyltransferase responsible for transposon methylation, associates with SPOCD1 and C19ORF84 in fetal gonocytes. Furthermore, C19ORF84 is essential for piRNA-directed DNA methylation and male mouse fertility. Finally, C19ORF84 mediates the in vivo association of SPOCD1 with the de novo methylation machinery. In summary, we have discovered a conserved role for the human piRNA pathway in transposon silencing and C19ORF84, an uncharacterized protein essential for orchestrating piRNA-directed DNA methylation.

The phosphorylated trimeric SOSS1 complex and RNA polymerase II trigger liquid-liquid phase separation at double-strand breaks

Double-strand breaks (DSBs) are the most severe type of DNA damage. Previously, we demonstrated that RNA polymerase II (RNAPII) phosphorylated at the tyrosine 1 (Y1P) residue of its C-terminal domain (CTD) generates RNAs at DSBs. However, the regulation of transcription at DSBs remains enigmatic. Here, we show that the damage-activated tyrosine kinase c-Abl phosphorylates hSSB1, enabling its interaction with Y1P RNAPII at DSBs. Furthermore, the trimeric SOSS1 complex, consisting of hSSB1, INTS3, and c9orf80, binds to Y1P RNAPII in response to DNA damage in an R-loop-dependent manner. Specifically, hSSB1, as a part of the trimeric SOSS1 complex, exhibits a strong affinity for R-loops, even in the presence of replication protein A (RPA). Our in vitro and in vivo data reveal that the SOSS1 complex and RNAPII form dynamic liquid-like repair compartments at DSBs. Depletion of the SOSS1 complex impairs DNA repair, underscoring its biological role in the R-loop-dependent DNA damage response.

The SPOC proteins DIDO3 and PHF3 co-regulate gene expression and neuronal differentiation

Transcription is regulated by a multitude of activators and repressors, which bind to the RNA polymerase II (Pol II) machinery and modulate its progression. Death-inducer obliterator 3 (DIDO3) and PHD finger protein 3 (PHF3) are paralogue proteins that regulate transcription elongation by docking onto phosphorylated serine-2 in the C-terminal domain (CTD) of Pol II through their SPOC domains. Here, we show that DIDO3 and PHF3 form a complex that bridges the Pol II elongation machinery with chromatin and RNA processing factors and tethers Pol II in a phase-separated microenvironment. Their SPOC domains and C-terminal intrinsically disordered regions are critical for transcription regulation. PHF3 and DIDO exert cooperative and antagonistic effects on the expression of neuronal genes and are both essential for neuronal differentiation. In the absence of PHF3, DIDO3 is upregulated as a compensatory mechanism. In addition to shared gene targets, DIDO specifically regulates genes required for lipid metabolism. Collectively, our work reveals multiple layers of gene expression regulation by the DIDO3 and PHF3 paralogues, which have specific, co-regulatory and redundant functions in transcription.

Whole genome analysis of rare deleterious variants adds further evidence to BRSK2 and other risk genes in Autism Spectrum Disorder

Autism spectrum disorder (ASD) is a complex neurodevelopmental condition with a strong genetic component in which rare variants contribute significantly to risk. We have performed whole genome and/or exome sequencing (WGS and WES) and SNP-array analysis to identify both rare sequence and copy number variants (SNVs and CNVs) in 435 individuals from 116 ASD families. We identified 37 rare potentially damaging de novo SNVs (pdSNVs) in cases (n = 144). Interestingly, two of them (one stop-gain and one missense variant) occurred in the same gene, BRSK2. Moreover, the identification of 9 severe de novo pdSNVs in genes not previously implicated in ASD (RASAL2, RAP1A, IRX5, SLC9A1, AGPAT3, MGAT3, RAB8B, MGAT5B, YME1L1), highlighted novel candidates. Potentially damaging CNVs (pdCNVs) provided support to the involvement of inherited variants in PHF3, NEGR1, TIAM1 and HOMER1 in neurodevelopmental disorders (NDD), although mostly acting as susceptibility factors with incomplete penetrance. Interpretation of identified pdSNVs/pdCNVs according to the ACMG guidelines led to a molecular diagnosis in 19/144 cases, but this figure represents a lower limit and is expected to increase thanks to further clarification of the role of likely pathogenic variants in new ASD/NDD candidates. In conclusion, our study strengthens the role of BRSK2 and other neurodevelopmental genes in ASD risk, highlights novel candidates and contributes to characterize the allelic diversity, mode of inheritance and phenotypic impact of de novo and inherited risk variants in ASD/NDD genes.

Distinctive interactomes of RNA polymerase II phosphorylation during different stages of transcription

During eukaryotic transcription, RNA polymerase II undergoes dynamic post-translational modifications on the C-terminal domain (CTD) of the largest subunit, generating an information-rich PTM landscape that transcriptional regulators bind. The phosphorylation of Ser5 and Ser2 of CTD heptad occurs spatiotemporally with the transcriptional stages, recruiting different transcriptional regulators to Pol II. To delineate the protein interactomes at different transcriptional stages, we reconstructed phosphorylation patterns of the CTD at Ser5 and Ser2 in vitro. Our results showed that distinct protein interactomes are recruited to RNA polymerase II at different stages of transcription by the phosphorylation of Ser2 and Ser5 of the CTD heptads. In particular, we characterized calcium homeostasis endoplasmic reticulum protein (CHERP) as a regulator bound by phospho-Ser2 heptad. Pol II association with CHERP recruits an accessory splicing complex whose loss results in broad changes in alternative splicing events. Our results shed light on the PTM-coded recruitment process that coordinates transcription.

RNA polymerase II depletion from the inactive X chromosome territory is not mediated by physical compartmentalization

Subnuclear compartmentalization has been proposed to play an important role in gene regulation by segregating active and inactive parts of the genome in distinct physical and biochemical environments. During X chromosome inactivation (XCI), the noncoding Xist RNA coats the X chromosome, triggers gene silencing and forms a dense body of heterochromatin from which the transcription machinery appears to be excluded. Phase separation has been proposed to be involved in XCI, and might explain the exclusion of the transcription machinery by preventing its diffusion into the Xist-coated territory. Here, using quantitative fluorescence microscopy and single-particle tracking, we show that RNA polymerase II (RNAPII) freely accesses the Xist territory during the initiation of XCI. Instead, the apparent depletion of RNAPII is due to the loss of its chromatin stably bound fraction. These findings indicate that initial exclusion of RNAPII from the inactive X reflects the absence of actively transcribing RNAPII, rather than a consequence of putative physical compartmentalization of the inactive X heterochromatin domain.

Loss of circIGF1R Suppresses Cardiomyocytes Proliferation by Sponging miR-362-5p

Circular RNAs (circRNAs) are generally formed by the back-splicing of precursor mRNA. Increasing evidence implicates the important role of circRNAs in cardiovascular diseases. However, the role of circ-insulin-like growth factor 1 receptor (circIGF1R) in cardiomyocyte (CM) proliferation remains unclear. Here, we investigated the potential role of the circIGF1R in the proliferation of CMs. We found that circIGF1R expression in heart tissues and primary CMs from adult mice was significantly lower than that in neonatal mice at postnatal 1 day (p1). Increased circIGF1R expression was detected in the injured neonatal heart at 0.5 and 1 days post-resection. circIGF1R knockdown significantly decreased the proliferation of primary CMs. Combined prediction software, luciferase reporter gene analysis, and quantitative real time-PCR (qPCR) revealed that circIGF1R interacted with miR-362-5p. A significant increase in miR-362-5p expression was detected in the adult heart compared with that in the neonatal heart. Further, heart injury significantly decreased the expression of miR-362-5p in neonatal mice. Treatment with miR-362-5p mimics significantly suppressed the proliferation of primary CMs, whereas knockdown of miR-362-5p promoted the CMs proliferation. Meanwhile, miR-362-5p silencing can rescue the proliferation inhibition of CMs induced by circIGF1R knockdown. Target prediction and qPCR validation revealed that miR-362-5p significantly inhibited the expression of Phf3 in primary CMs. In addition, decreased Phf3 expression was detected in adult hearts compared with neonatal hearts. Consistently, increased Phf3 expression was detected in injured neonatal hearts compared with that in sham hearts. Knockdown of Phf3 markedly repressed CMs proliferation. Taken together, these findings suggest that circIGF1R might contribute to cardiomyocyte proliferation by promoting Pfh3 expression by sponging miR-362-5p and provide an important experimental basis for the regulation of heart regeneration.

Dysregulated phosphate metabolism in autism spectrum disorder: associations and insights for future research

Studies of autism spectrum disorder (ASD) related to exposure to toxic levels of dietary phosphate are lacking. Phosphate toxicity from dysregulated phosphate metabolism can negatively impact almost every major organ system of the body, including the central nervous system. The present paper used a grounded theory-literature review method to synthesise associations of dysregulated phosphate metabolism with the aetiology of ASD. Cell signalling in autism has been linked to an altered balance between phosphoinositide kinases, which phosphorylate proteins, and the counteracting effect of phosphatases in neuronal membranes. Glial cell overgrowth in the developing ASD brain can lead to disturbances in neuro-circuitry, neuroinflammation and immune responses which are potentially related to excessive inorganic phosphate. The rise in ASD prevalence has been suggested to originate in changes to the gut microbiome from increasing consumption of additives in processed food, including phosphate additives. Ketogenic diets and dietary patterns that eliminate casein also reduce phosphate intake, which may account for many of the suggested benefits of these diets in children with ASD. Dysregulated phosphate metabolism is causatively linked to comorbid conditions associated with ASD such as cancer, tuberous sclerosis, mitochondrial dysfunction, diabetes, epilepsy, obesity, chronic kidney disease, tauopathy, cardiovascular disease and bone mineral disorders. Associations and proposals presented in this paper offer novel insights and directions for future research linking the aetiology of ASD with dysregulated phosphate metabolism and phosphate toxicity from excessive dietary phosphorus intake.

Late-Stage Metastatic Melanoma Emerges through a Diversity of Evolutionary Pathways

Understanding the evolutionary pathways to metastasis and resistance to immune-checkpoint inhibitors (ICI) in melanoma is critical for improving outcomes. Here, we present the most comprehensive intrapatient metastatic melanoma dataset assembled to date as part of the Posthumous Evaluation of Advanced Cancer Environment (PEACE) research autopsy program, including 222 exome sequencing, 493 panel-sequenced, 161 RNA sequencing, and 22 single-cell whole-genome sequencing samples from 14 ICI-treated patients. We observed frequent whole-genome doubling and widespread loss of heterozygosity, often involving antigen-presentation machinery. We found KIT extrachromosomal DNA may have contributed to the lack of response to KIT inhibitors of a KIT-driven melanoma. At the lesion-level, MYC amplifications were enriched in ICI nonresponders. Single-cell sequencing revealed polyclonal seeding of metastases originating from clones with different ploidy in one patient. Finally, we observed that brain metastases that diverged early in molecular evolution emerge late in disease. Overall, our study illustrates the diverse evolutionary landscape of advanced melanoma.

SIGNIFICANCE

Despite treatment advances, melanoma remains a deadly disease at stage IV. Through research autopsy and dense sampling of metastases combined with extensive multiomic profiling, our study elucidates the many mechanisms that melanomas use to evade treatment and the immune system, whether through mutations, widespread copy-number alterations, or extrachromosomal DNA. See related commentary by Shain, p. 1294. This article is highlighted in the In This Issue feature, p. 1275.

Genetic predisposition and evolutionary traces of pediatric cancer risk: a prospective 5-year population-based genome sequencing study of children with CNS tumors

BACKGROUND

The etiology of central nervous system (CNS) tumors in children is largely unknown and population-based studies of genetic predisposition are lacking.

METHODS

In this prospective, population-based study, we performed germline whole-genome sequencing in 128 children with CNS tumors, supplemented by a systematic pedigree analysis covering 3543 close relatives.

RESULTS

Thirteen children (10%) harbored pathogenic variants in known cancer genes. These children were more likely to have medulloblastoma (OR 5.9, CI 1.6-21.2) and develop metasynchronous CNS tumors (P = 0.01). Similar carrier frequencies were seen among children with low-grade glioma (12.8%) and high-grade tumors (12.2%). Next, considering the high mortality of childhood CNS tumors throughout most of human evolution, we explored known pediatric-onset cancer genes, showing that they are more evolutionarily constrained than genes associated with risk of adult-onset malignancies (P = 5e-4) and all other genes (P = 5e-17). Based on this observation, we expanded our analysis to 2986 genes exhibiting high evolutionary constraint in 141,456 humans. This analysis identified eight directly causative loss-of-functions variants, and showed a dose-response association between degree of constraint and likelihood of pathogenicity-raising the question of the role of other highly constrained gene alterations detected.

CONCLUSIONS

Approximately 10% of pediatric CNS tumors can be attributed to rare variants in known cancer genes. Genes associated with high risk of childhood cancer show evolutionary evidence of constraint.

SPOC domain proteins in health and disease

Since it was first described >20 yr ago, the SPOC domain (Spen paralog and ortholog C-terminal domain) has been identified in many proteins all across eukaryotic species. SPOC-containing proteins regulate gene expression on various levels ranging from transcription to RNA processing, modification, export, and stability, as well as X-chromosome inactivation. Their manifold roles in controlling transcriptional output implicate them in a plethora of developmental processes, and their misregulation is often associated with cancer. Here, we provide an overview of the biophysical properties of the SPOC domain and its interaction with phosphorylated binding partners, the phylogenetic origin of SPOC domain proteins, the diverse functions of mammalian SPOC proteins and their homologs, the mechanisms by which they regulate differentiation and development, and their roles in cancer.

The SPOC domain is a phosphoserine binding module that bridges transcription machinery with co- and post-transcriptional regulators

The heptad repeats of the C-terminal domain (CTD) of RNA polymerase II (Pol II) are extensively modified throughout the transcription cycle. The CTD coordinates RNA synthesis and processing by recruiting transcription regulators as well as RNA capping, splicing and 3'end processing factors. The SPOC domain of PHF3 was recently identified as a CTD reader domain specifically binding to phosphorylated serine-2 residues in adjacent CTD repeats. Here, we establish the SPOC domains of the human proteins DIDO, SHARP (also known as SPEN) and RBM15 as phosphoserine binding modules that can act as CTD readers but also recognize other phosphorylated binding partners. We report the crystal structure of SHARP SPOC in complex with CTD and identify the molecular determinants for its specific binding to phosphorylated serine-5. PHF3 and DIDO SPOC domains preferentially interact with the Pol II elongation complex, while RBM15 and SHARP SPOC domains engage with writers and readers of m6A, the most abundant RNA modification. RBM15 positively regulates m6A levels and mRNA stability in a SPOC-dependent manner, while SHARP SPOC is essential for its localization to inactive X-chromosomes. Our findings suggest that the SPOC domain is a major interface between the transcription machinery and regulators of transcription and co-transcriptional processes.

Epigenetic regulator genes direct lineage switching in MLL/AF4 leukemia

The fusion gene MLL/AF4 defines a high-risk subtype of pro-B acute lymphoblastic leukemia. Relapse can be associated with a lineage switch from acute lymphoblastic to acute myeloid leukemia, resulting in poor clinical outcomes caused by resistance to chemotherapies and immunotherapies. In this study, the myeloid relapses shared oncogene fusion breakpoints with their matched lymphoid presentations and originated from various differentiation stages from immature progenitors through to committed B-cell precursors. Lineage switching is linked to substantial changes in chromatin accessibility and rewiring of transcriptional programs, including alternative splicing. These findings indicate that the execution and maintenance of lymphoid lineage differentiation is impaired. The relapsed myeloid phenotype is recurrently associated with the altered expression, splicing, or mutation of chromatin modifiers, including CHD4 coding for the ATPase/helicase of the nucleosome remodelling and deacetylation complex. Perturbation of CHD4 alone or in combination with other mutated epigenetic modifiers induces myeloid gene expression in MLL/AF4+ cell models, indicating that lineage switching in MLL/AF4 leukemia is driven and maintained by disrupted epigenetic regulation.