Targeted Apoptosis of Senescent Cells Restores Tissue Homeostasis in Response to Chemotoxicity and Aging.

The accumulation of irreparable cellular damage restricts healthspan after acute stress or natural aging. Senescent cells are thought to impair tissue function, and their genetic clearance can delay features of aging. Identifying how senescent cells avoid apoptosis allows for the prospective design of anti-senescence compounds to address whether homeostasis can also be restored. Here, we identify FOXO4 as a pivot in senescent cell viability. We designed a FOXO4 peptide that perturbs the FOXO4 interaction with p53. In senescent cells, this selectively causes p53 nuclear exclusion and cell-intrinsic apoptosis. Under conditions where it was well tolerated in vivo, this FOXO4 peptide neutralized doxorubicin-induced chemotoxicity. Moreover, it restored fitness, fur density, and renal function in both fast aging XpdTTD/TTD and naturally aged mice. Thus, therapeutic targeting of senescent cells is feasible under conditions where loss of health has already occurred, and in doing so tissue homeostasis can effectively be restored.

NOXA expression drives synthetic lethality to RUNX1 inhibition in pancreatic cancer.

Evasion from drug-induced apoptosis is a crucial mechanism of cancer treatment resistance. The proapoptotic protein NOXA marks an aggressive pancreatic ductal adenocarcinoma (PDAC) subtype. To identify drugs that unleash the death-inducing potential of NOXA, we performed an unbiased drug screening experiment. In NOXA-deficient isogenic cellular models, we identified an inhibitor of the transcription factor heterodimer CBFß/RUNX1. By genetic gain and loss of function experiments, we validated that the mode of action depends on RUNX1 and NOXA. Of note is that RUNX1 expression is significantly higher in PDACs compared to normal pancreas. We show that pharmacological RUNX1 inhibition significantly blocks tumor growth in vivo and in primary patient-derived PDAC organoids. Through genome-wide analysis, we detected that RUNX1-loss reshapes the epigenetic landscape, which gains H3K27ac enrichment at the NOXA promoter. Our study demonstrates a previously unknown mechanism of NOXA-dependent cell death, which can be triggered pharmaceutically. Therefore, our data show a way to target a therapy-resistant PDAC, an unmet clinical need.

Split Pool Ligation-based Single-cell Transcriptome sequencing (SPLiT-seq) data processing pipeline comparison.

BACKGROUND: Single-cell sequencing techniques are revolutionizing every field of biology by providing the ability to measure the abundance of biological molecules at a single-cell resolution. Although single-cell sequencing approaches have been developed for several molecular modalities, single-cell transcriptome sequencing is the most prevalent and widely applied technique. SPLiT-seq (split-pool ligation-based transcriptome sequencing) is one of these single-cell transcriptome techniques that applies a unique combinatorial-barcoding approach by splitting and pooling cells into multi-well plates containing barcodes. This unique approach required the development of dedicated computational tools to preprocess the data and extract the count matrices. Here we compare eight bioinformatic pipelines (alevin-fry splitp, LR-splitpipe, SCSit, splitpipe, splitpipeline, SPLiTseq-demultiplex, STARsolo and zUMI) that have been developed to process SPLiT-seq data. We provide an overview of the tools, their computational performance, functionality and impact on downstream processing of the single-cell data, which vary greatly depending on the tool used. RESULTS: We show that STARsolo, splitpipe and alevin-fry splitp can all handle large amount of data within reasonable time. In contrast, the other five pipelines are slow when handling large datasets. When using smaller dataset, cell barcode results are similar with the exception of SPLiTseq-demultiplex and splitpipeline. LR-splitpipe that is originally designed for processing long-read sequencing data is the slowest of all pipelines. Alevin-fry produced different down-stream results that are difficult to interpret. STARsolo functions nearly identical to splitpipe and produce results that are highly similar to each other. However, STARsolo lacks the function to collapse random hexamer reads for which some additional coding is required. CONCLUSION: Our comprehensive comparative analysis aids users in selecting the most suitable analysis tool for efficient SPLiT-seq data processing, while also detailing the specific prerequisites for each of these pipelines. From the available pipelines, we recommend splitpipe or STARSolo for SPLiT-seq data analysis.

Distinct IL-1α-responsive enhancers promote acute and coordinated changes in chromatin topology in a hierarchical manner.

How cytokine-driven changes in chromatin topology are converted into gene regulatory circuits during inflammation still remains unclear. Here, we show that interleukin (IL)-1α induces acute and widespread changes in chromatin accessibility via the TAK1 kinase and NF-κB at regions that are highly enriched for inflammatory disease-relevant SNPs. Two enhancers in the extended chemokine locus on human chromosome 4 regulate the IL-1α-inducible IL8 and CXCL1-3 genes. Both enhancers engage in dynamic spatial interactions with gene promoters in an IL-1α/TAK1-inducible manner. Microdeletions of p65-binding sites in either of the two enhancers impair NF-κB recruitment, suppress activation and biallelic transcription of the IL8/CXCL2 genes, and reshuffle higher-order chromatin interactions as judged by i4C interactome profiles. Notably, these findings support a dominant role of the IL8 "master" enhancer in the regulation of sustained IL-1α signaling, as well as for IL-8 and IL-6 secretion. CRISPR-guided transactivation of the IL8 locus or cross-TAD regulation by TNFα-responsive enhancers in a different model locus supports the existence of complex enhancer hierarchies in response to cytokine stimulation that prime and orchestrate proinflammatory chromatin responses downstream of NF-κB.

Epigenomic partitioning of a polygenic risk score for asthma reveals distinct genetically driven disease pathways.

BACKGROUND: Individual differences in susceptibility to develop asthma, a heterogeneous chronic inflammatory lung disease, are poorly understood. It remains debated whether genetics can predict asthma risk and how genetic variants modulate the complex pathophysiology of asthma. AIM: To build polygenic risk scores (PRSs) for asthma risk prediction and epigenomically link predictive genetic variants to pathophysiological mechanisms. METHODS: Restricted PRSs were constructed using single nucleotide variants derived from genome-wide association studies and validated using data generated in the Rotterdam Study, a Dutch prospective cohort of 14 926 individuals. Outcomes used were asthma, childhood-onset asthma (COA), adulthood-onset asthma (AOA), eosinophilic asthma, and asthma exacerbations. Genome-wide chromatin analysis data from 19 disease-relevant cell types were used for epigenomic PRS partitioning. RESULTS: PRSs obtained predicted asthma and related outcomes, with the strongest associations observed for COA (2.55 odds ratios per PRS standard deviation, area under the curve of 0.760). PRSs allowed for the classification of individuals into high and low-risk groups. PRS partitioning using epigenomic profiles identified 5 clusters of variants within putative gene regulatory regions linked to specific asthma-relevant cells, genes, and biological pathways. CONCLUSIONS: PRSs were associated with asthma(-related traits) in a Dutch prospective cohort, with substantially higher predictive power observed for COA than for AOA. Importantly, PRS variants could be epigenomically partitioned into clusters of regulatory variants with different pathophysiological association patterns and effect estimates, which likely represent distinct genetically driven disease pathways. Our findings have potential implications for personalized risk mitigation and treatment strategies.

LSD1/KDM1A and GFI1B repress endothelial fate and induce hematopoietic fate in induced pluripotent stem cell-derived hemogenic endotheliums.

Differentiation of induced pluripotent stem cells (iPSCs) into hematopoietic lineages offers great therapeutic potential. During embryogenesis, hemogenic endothelium (HE) gives rise to hematopoietic stem and progenitor cells through the endothelial-to-hematopoietic transition (EHT). Understanding this process using iPSCs is key to generating functional hematopoietic stem cells (HSCs), a currently unmet challenge. In this study, we examined the role of the transcriptional factor GFI1B and its co-factor LSD1/KDM1A in EHT. To this end, we employed patient-derived iPSC lines with a dominant negative dysfunctional GFI1BQ287* and irreversible pharmacological LSD1/KDM1A inhibition in healthy iPSC lines. The formation of HE remained unaffected; however, hematopoietic output was severely reduced in both conditions. Single-cell RNA sequencing (scRNAseq) performed on the CD144+/CD31+ population derived from healthy iPSCs revealed similar expression dynamics of genes associated with in vivo EHT. Interestingly, LSD1/KDM1A inhibition in healthy lines before EHT resulted in a complete absence of hematopoietic output. However, uncommitted HE cells did not display GFI1B expression, suggesting a timed transcriptional program. To test this hypothesis, we ectopically expressed GFI1B in uncommitted HE cells, leading to downregulation of endothelial genes and upregulation of hematopoietic genes, including GATA2, KIT, RUNX1, and SPI1. Thus, we demonstrate that LSD1/KDM1A and GFI1B can function at distinct temporal points in different cellular subsets during EHT. Although GFI1B is not detected in uncommitted HE cells, its ectopic expression allows for partial hematopoietic specification. These data indicate that precisely timed expression of specific transcriptional regulators during EHT is crucial to the eventual outcome of EHT.

Genome-wide analysis toward the epigenetic aetiology of myelodysplastic syndrome disease progression and pharmacoepigenomic basis of hypomethylating agents drug treatment response.

Myelodysplastic syndromes (MDS) consist of a group of hematological malignancies characterized by ineffective hematopoiesis, cytogenetic abnormalities, and often a high risk of transformation to acute myeloid leukemia (AML). So far, there have been only a very limited number of studies assessing the epigenetics component contributing to the pathophysiology of these disorders, but not a single study assessing this at a genome-wide level. Here, we implemented a generic high throughput epigenomics approach, using methylated DNA sequencing (MeD-seq) of LpnPI digested fragments to identify potential epigenomic targets associated with MDS subtypes. Our results highlighted that PCDHG and ZNF gene families harbor potential epigenomic targets, which have been shown to be differentially methylated in a variety of comparisons between different MDS subtypes. Specifically, CpG islands, transcription start sites and post-transcriptional start sites within ZNF124, ZNF497 and PCDHG family are differentially methylated with fold change above 3,5. Overall, these findings highlight important aspects of the epigenomic component of MDS syndromes pathogenesis and the pharmacoepigenomic basis to the hypomethylating agents drug treatment response, while this generic high throughput whole epigenome sequencing approach could be readily implemented to other genetic diseases with a strong epigenetic component.

Redundant and specific roles of cohesin STAG subunits in chromatin looping and transcriptional control.

Cohesin is a ring-shaped multiprotein complex that is crucial for 3D genome organization and transcriptional regulation during differentiation and development. It also confers sister chromatid cohesion and facilitates DNA damage repair. Besides its core subunits SMC3, SMC1A, and RAD21, cohesin in somatic cells contains one of two orthologous STAG subunits, STAG1 or STAG2. How these variable subunits affect the function of the cohesin complex is still unclear. STAG1- and STAG2-cohesin were initially proposed to organize cohesion at telomeres and centromeres, respectively. Here, we uncover redundant and specific roles of STAG1 and STAG2 in gene regulation and chromatin looping using HCT116 cells with an auxin-inducible degron (AID) tag fused to either STAG1 or STAG2. Following rapid depletion of either subunit, we perform high-resolution Hi-C, gene expression, and sequential ChIP studies to show that STAG1 and STAG2 do not co-occupy individual binding sites and have distinct ways by which they affect looping and gene expression. These findings are further supported by single-molecule localizations via direct stochastic optical reconstruction microscopy (dSTORM) super-resolution imaging. Since somatic and congenital mutations of the STAG subunits are associated with cancer (STAG2) and intellectual disability syndromes with congenital abnormalities (STAG1 and STAG2), we verified STAG1-/STAG2-dependencies using human neural stem cells, hence highlighting their importance in particular disease contexts.

Multifaceted actions of Zeb2 in postnatal neurogenesis from the ventricular-subventricular zone to the olfactory bulb.

The transcription factor Zeb2 controls fate specification and subsequent differentiation and maturation of multiple cell types in various embryonic tissues. It binds many protein partners, including activated Smad proteins and the NuRD co-repressor complex. How Zeb2 subdomains support cell differentiation in various contexts has remained elusive. Here, we studied the role of Zeb2 and its domains in neurogenesis and neural differentiation in the young postnatal ventricular-subventricular zone (V-SVZ), in which neural stem cells generate olfactory bulb-destined interneurons. Conditional Zeb2 knockouts and separate acute loss- and gain-of-function approaches indicated that Zeb2 is essential for controlling apoptosis and neuronal differentiation of V-SVZ progenitors before and after birth, and we identified Sox6 as a potential downstream target gene of Zeb2. Zeb2 genetic inactivation impaired the differentiation potential of the V-SVZ niche in a cell-autonomous fashion. We also provide evidence that its normal function in the V-SVZ also involves non-autonomous mechanisms. Additionally, we demonstrate distinct roles for Zeb2 protein-binding domains, suggesting that Zeb2 partners co-determine neuronal output from the mouse V-SVZ in both quantitative and qualitative ways in early postnatal life.

Kidney Organoids Are Capable of Forming Tumors, but Not Teratomas.

Induced pluripotent stem cell (iPSC)-derived kidney organoids are a potential tool for the regeneration of kidney tissue. They represent an early stage of nephrogenesis and have been shown to successfsully vascularize and mature further in vivo. However, there are concerns regarding the long-term safety and stability of iPSC derivatives. Specifically, the potential for tumorigenesis may impede the road to clinical application. To study safety and stability of kidney organoids, we analyzed their potential for malignant transformation in a teratoma assay and following long-term subcutaneous implantation in an immune-deficient mouse model. We did not detect fully functional residual iPSCs in the kidney organoids as analyzed by gene expression analysis, single-cell sequencing and immunohistochemistry. Accordingly, kidney organoids failed to form teratoma. Upon long-term subcutaneous implantation of whole organoids in immunodeficient IL2Ry-/-RAG2-/- mice, we observed tumor formation in 5 out of 103 implanted kidney organoids. These tumors were composed of WT1+CD56+ immature blastemal cells and showed histological resemblance with Wilms tumor. No genetic changes were identified that contributed to the occurrence of tumorigenic cells within the kidney organoids. However, assessment of epigenetic changes revealed a unique cluster of differentially methylated genes that were also present in undifferentiated iPSCs. We discovered that kidney organoids have the capacity to form tumors upon long-term implantation. The presence of epigenetic modifications combined with the lack of environmental cues may have caused an arrest in terminal differentiation. Our results indicate that the safe implementation of kidney organoids should exclude the presence of pro-tumorigenic methylation in kidney organoids.

Mutation Analysis of Pancreatic Juice and Plasma for the Detection of Pancreatic Cancer.

Molecular profiling may enable earlier detection of pancreatic cancer (PC) in high-risk individuals undergoing surveillance and allow for personalization of treatment. We hypothesized that the detection rate of DNA mutations is higher in pancreatic juice (PJ) than in plasma due to its closer contact with the pancreatic ductal system, from which pancreatic cancer cells originate, and higher overall cell-free DNA (cfDNA) concentrations. In this study, we included patients with pathology-proven PC or intraductal papillary mucinous neoplasm (IPMN) with high-grade dysplasia (HGD) from two prospective clinical trials (KRASPanc and PACYFIC) for whom both PJ and plasma were available. We performed next-generation sequencing on PJ, plasma, and tissue samples and described the presence (and concordance) of mutations in these biomaterials. This study included 26 patients (25 PC and 1 IPMN with HGD), of which 7 were women (27%), with a median age of 71 years (IQR 12) and a median BMI of 23 kg/m2 (IQR 4). Ten patients with PC (40%) were (borderline) resectable at baseline. Tissue was available from six patients (resection n = 5, biopsy n = 1). A median volume of 2.9 mL plasma (IQR 1.0 mL) and 0.7 mL PJ (IQR 0.1 mL, p < 0.001) was used for DNA isolation. PJ had a higher median cfDNA concentration (2.6 ng/µL (IQR 4.2)) than plasma (0.29 ng/µL (IQR 0.40)). A total of 41 unique somatic mutations were detected: 24 mutations in plasma (2 KRAS, 15 TP53, 2 SMAD4, 3 CDKN2A 1 CTNNB1, and 1 PIK3CA), 19 in PJ (3 KRAS, 15 TP53, and 1 SMAD4), and 8 in tissue (2 KRAS, 2 CDKN2A, and 4 TP53). The mutation detection rate (and the concordance with tissue) did not differ between plasma and PJ. In conclusion, while the concentration of cfDNA was indeed higher in PJ than in plasma, the mutation detection rate was not different. A few cancer-associated genetic variants were detected in both biomaterials. Further research is needed to increase the detection rate and assess the performance and suitability of plasma and PJ for PC (early) detection.

Single-Cell RNA Sequencing of Donor-Reactive T Cells Reveals Role of Apoptosis in Donor-Specific Hyporesponsiveness of Kidney Transplant Recipients.

After kidney transplantation (KT), donor-specific hyporesponsiveness (DSH) of recipient T cells develops over time. Recently, apoptosis was identified as a possible underlying mechanism. In this study, both transcriptomic profiles and complete V(D)J variable regions of TR transcripts from individual alloreactive T cells of kidney transplant recipients were determined with single-cell RNA sequencing. Alloreactive T cells were identified by CD137 expression after stimulation of peripheral blood mononuclear cells obtained from KT recipients (N = 7) prior to and 3-5 years after transplantation with cells of their donor or a third party control. The alloreactive T cells were sorted, sequenced and the transcriptome and T cell receptor profiles were analyzed using unsupervised clustering. Alloreactive T cells retain a highly polyclonal T Cell Receptor Alpha/Beta repertoire over time. Post transplantation, donor-reactive CD4+ T cells had a specific downregulation of genes involved in T cell cytokine-mediated pathways and apoptosis. The CD8+ donor-reactive T cell profile did not change significantly over time. Single-cell expression profiling shows that activated and pro-apoptotic donor-reactive CD4+ T cell clones are preferentially lost after transplantation in stable kidney transplant recipients.

Targeted chromatin conformation analysis identifies novel distal neural enhancers of ZEB2 in pluripotent stem cell differentiation.

The transcription factor zinc finger E-box binding protein 2 (ZEB2) controls embryonic and adult cell fate decisions and cellular maturation in many stem/progenitor cell types. Defects in these processes in specific cell types underlie several aspects of Mowat-Wilson syndrome (MOWS), which is caused by ZEB2 haplo-insufficiency. Human ZEB2, like mouse Zeb2, is located on chromosome 2 downstream of a ±3.5 Mb-long gene-desert, lacking any protein-coding gene. Using temporal targeted chromatin capture (T2C), we show major chromatin structural changes based on mapping in-cis proximities between the ZEB2 promoter and this gene desert during neural differentiation of human-induced pluripotent stem cells, including at early neuroprogenitor cell (NPC)/rosette state, where ZEB2 mRNA levels increase significantly. Combining T2C with histone-3 acetylation mapping, we identified three novel candidate enhancers about 500 kb upstream of the ZEB2 transcription start site. Functional luciferase-based assays in heterologous cells and NPCs reveal co-operation between these three enhancers. This study is the first to document in-cis Regulatory Elements located in ZEB2's gene desert. The results further show the usability of T2C for future studies of ZEB2 REs in differentiation and maturation of multiple cell types and the molecular characterization of newly identified MOWS patients that lack mutations in ZEB2 protein-coding exons.

Rapid in vitro generation of bona fide exhausted CD8+ T cells is accompanied by Tcf7 promotor methylation.

Exhaustion is a dysfunctional state of cytotoxic CD8+ T cells (CTL) observed in chronic infection and cancer. Current in vivo models of CTL exhaustion using chronic viral infections or cancer yield very few exhausted CTL, limiting the analysis that can be done on these cells. Establishing an in vitro system that rapidly induces CTL exhaustion would therefore greatly facilitate the study of this phenotype, identify the truly exhaustion-associated changes and allow the testing of novel approaches to reverse or prevent exhaustion. Here we show that repeat stimulation of purified TCR transgenic OT-I CTL with their specific peptide induces all the functional (reduced cytokine production and polyfunctionality, decreased in vivo expansion capacity) and phenotypic (increased inhibitory receptors expression and transcription factor changes) characteristics of exhaustion. Importantly, in vitro exhausted cells shared the transcriptomic characteristics of the gold standard of exhaustion, CTL from LCMV cl13 infections. Gene expression of both in vitro and in vivo exhausted CTL was distinct from T cell anergy. Using this system, we show that Tcf7 promoter DNA methylation contributes to TCF1 downregulation in exhausted CTL. Thus this novel in vitro system can be used to identify genes and signaling pathways involved in exhaustion and will facilitate the screening of reagents that prevent/reverse CTL exhaustion.

Dominant-acting CSF1R variants cause microglial depletion and altered astrocytic phenotype in zebrafish and adult-onset leukodystrophy.

Tissue-resident macrophages of the brain, including microglia, are implicated in the pathogenesis of various CNS disorders and are possible therapeutic targets by their chemical depletion or replenishment by hematopoietic stem cell therapy. Nevertheless, a comprehensive understanding of microglial function and the consequences of microglial depletion in the human brain is lacking. In human disease, heterozygous variants in CSF1R, encoding the Colony-stimulating factor 1 receptor, can lead to adult-onset leukoencephalopathy with axonal spheroids and pigmented glia (ALSP) possibly caused by microglial depletion. Here, we investigate the effects of ALSP-causing CSF1R variants on microglia and explore the consequences of microglial depletion in the brain. In intermediate- and late-stage ALSP post-mortem brain, we establish that there is an overall loss of homeostatic microglia and that this is predominantly seen in the white matter. By introducing ALSP-causing missense variants into the zebrafish genomic csf1ra locus, we show that these variants act dominant negatively on the number of microglia in vertebrate brain development. Transcriptomics and proteomics on relatively spared ALSP brain tissue validated a downregulation of microglia-associated genes and revealed elevated astrocytic proteins, possibly suggesting involvement of astrocytes in early pathogenesis. Indeed, neuropathological analysis and in vivo imaging of csf1r zebrafish models showed an astrocytic phenotype associated with enhanced, possibly compensatory, endocytosis. Together, our findings indicate that microglial depletion in zebrafish and human disease, likely as a consequence of dominant-acting pathogenic CSF1R variants, correlates with altered astrocytes. These findings underscore the unique opportunity CSF1R variants provide to gain insight into the roles of microglia in the human brain, and the need to further investigate how microglia, astrocytes, and their interactions contribute to white matter homeostasis.

Genome-wide DNA methylation profiling using the methylation-dependent restriction enzyme LpnPI.

DNA methylation is a well-known epigenetic modification that plays a crucial role in gene regulation, but genome-wide analysis of DNA methylation remains technically challenging and costly. DNA methylation-dependent restriction enzymes can be used to restrict CpG methylation analysis to methylated regions of the genome only, which significantly reduces the required sequencing depth and simplifies subsequent bioinformatics analysis. Unfortunately, this approach has been hampered by complete digestion of DNA in CpG methylation-dense regions, resulting in fragments that are too small for accurate mapping. Here, we show that the activity of DNA methylation-dependent enzyme, LpnPI, is blocked by a fragment size smaller than 32 bp. This unique property prevents complete digestion of methylation-dense DNA and allows accurate genome-wide analysis of CpG methylation at single-nucleotide resolution. Methylated DNA sequencing (MeD-seq) of LpnPI digested fragments revealed highly reproducible genome-wide CpG methylation profiles for >50% of all potentially methylated CpGs, at a sequencing depth less than one-tenth required for whole-genome bisulfite sequencing (WGBS). MeD-seq identified a high number of patient and tissue-specific differential methylated regions (DMRs) and revealed that patient-specific DMRs observed in both blood and buccal samples predict DNA methylation in other tissues and organs. We also observed highly variable DNA methylation at gene promoters on the inactive X Chromosome, indicating tissue-specific and interpatient-specific escape of X Chromosome inactivation. These findings highlight the potential of MeD-seq for high-throughput epigenetic profiling.

The core spliceosome as target and effector of non-canonical ATM signalling.

In response to DNA damage, tissue homoeostasis is ensured by protein networks promoting DNA repair, cell cycle arrest or apoptosis. DNA damage response signalling pathways coordinate these processes, partly by propagating gene-expression-modulating signals. DNA damage influences not only the abundance of messenger RNAs, but also their coding information through alternative splicing. Here we show that transcription-blocking DNA lesions promote chromatin displacement of late-stage spliceosomes and initiate a positive feedback loop centred on the signalling kinase ATM. We propose that initial spliceosome displacement and subsequent R-loop formation is triggered by pausing of RNA polymerase at DNA lesions. In turn, R-loops activate ATM, which signals to impede spliceosome organization further and augment ultraviolet-irradiation-triggered alternative splicing at the genome-wide level. Our findings define R-loop-dependent ATM activation by transcription-blocking lesions as an important event in the DNA damage response of non-replicating cells, and highlight a key role for spliceosome displacement in this process.

In vivo live imaging of RNA polymerase II transcription factories in primary cells.

Transcription steps are marked by different modifications of the C-terminal domain of RNA polymerase II (RNAPII). Phosphorylation of Ser5 and Ser7 by cyclin-dependent kinase 7 (CDK7) as part of TFIIH marks initiation, whereas phosphorylation of Ser2 by CDK9 marks elongation. These processes are thought to take place in localized transcription foci in the nucleus, known as "transcription factories," but it has been argued that the observed clusters/foci are mere fixation or labeling artifacts. We show that transcription factories exist in living cells as distinct foci by live-imaging fluorescently labeled CDK9, a kinase known to associate with active RNAPII. These foci were observed in different cell types derived from CDK9-mCherry knock-in mice. We show that these foci are very stable while highly dynamic in exchanging CDK9. Chromatin immunoprecipitation (ChIP) coupled with deep sequencing (ChIP-seq) data show that the genome-wide binding sites of CDK9 and initiating RNAPII overlap on transcribed genes. Immunostaining shows that CDK9-mCherry foci colocalize with RNAPII-Ser5P, much less with RNAPII-Ser2P, and not with CDK12 (a kinase reported to be involved in the Ser2 phosphorylation) or with splicing factor SC35. In conclusion, transcription factories exist in living cells, and initiation and elongation of transcripts takes place in different nuclear compartments.

The Somatic Mutation Paradigm in Congenital Malformations: Hirschsprung Disease as a Model.

Patients with Hirschsprung disease (HSCR) do not always receive a genetic diagnosis after routine screening in clinical practice. One of the reasons for this could be that the causal mutation is not present in the cell types that are usually tested-whole blood, dermal fibroblasts or saliva-but is only in the affected tissue. Such mutations are called somatic, and can occur in a given cell at any stage of development after conception. They will then be present in all subsequent daughter cells. Here, we investigated the presence of somatic mutations in HSCR patients. For this, whole-exome sequencing and copy number analysis were performed in DNA isolated from purified enteric neural crest cells (ENCCs) and blood or fibroblasts of the same patient. Variants identified were subsequently validated by Sanger sequencing. Several somatic variants were identified in all patients, but causative mutations for HSCR were not specifically identified in the ENCCs of these patients. Larger copy number variants were also not found to be specific to ENCCs. Therefore, we believe that somatic mutations are unlikely to be identified, if causative for HSCR. Here, we postulate various modes of development following the occurrence of a somatic mutation, to describe the challenges in detecting such mutations, and hypothesize how somatic mutations may contribute to 'missing heritability' in developmental defects.

Multi-Omics Profiling in Marfan Syndrome: Further Insights into the Molecular Mechanisms Involved in Aortic Disease.

Thoracic aortic aneurysm is a potentially life-threatening disease with a strong genetic contribution. Despite identification of multiple genes involved in aneurysm formation, little is known about the specific underlying mechanisms that drive the pathological changes in the aortic wall. The aim of our study was to unravel the molecular mechanisms underlying aneurysm formation in Marfan syndrome (MFS). We collected aortic wall samples from FBN1 variant-positive MFS patients (n = 6) and healthy donor hearts (n = 5). Messenger RNA (mRNA) expression levels were measured by RNA sequencing and compared between MFS patients and controls, and between haploinsufficient (HI) and dominant negative (DN) FBN1 variants. Immunohistochemical staining, proteomics and cellular respiration experiments were used to confirm our findings. FBN1 mRNA expression levels were highly variable in MFS patients and did not significantly differ from controls. Moreover, we did not identify a distinctive TGF-ß gene expression signature in MFS patients. On the contrary, differential gene and protein expression analysis, as well as vascular smooth muscle cell respiration measurements, pointed toward inflammation and mitochondrial dysfunction. Our findings confirm that inflammatory and mitochondrial pathways play important roles in the pathophysiological processes underlying MFS-related aortic disease, providing new therapeutic options.