The study of changes in protein-DNA interactions measured by ChIP-seq on dynamic systems, such as cell differentiation, response to treatments or the comparison of healthy and diseased individuals, is still an open challenge. There are few computational methods comparing changes in ChIP-seq signals with replicates. Moreover, none of these previous approaches addresses ChIP-seq specific experimental artefacts arising from studies with biological replicates. We propose THOR, a Hidden Markov Model based approach, to detect differential peaks between pairs of biological conditions with replicates. THOR provides all pre- and post-processing steps required in ChIP-seq analyses. Moreover, we propose a novel normalization approach based on housekeeping genes to deal with cases where replicates have distinct signal-to-noise ratios. To evaluate differential peak calling methods, we delineate a methodology using both biological and simulated data. This includes an evaluation procedure that associates differential peaks with changes in gene expression as well as histone modifications close to these peaks. We evaluate THOR and seven competing methods on data sets with distinct characteristics from in vitro studies with technical replicates to clinical studies of cancer patients. Our evaluation analysis comprises of 13 comparisons between pairs of biological conditions. We show that THOR performs best in all scenarios.

Drug-target networks are receiving a lot of attention in late years, given its relevance for pharmaceutical innovation and drug lead discovery. Different in silico approaches have been proposed for the identification of new drug-target interactions, many of which are based on kernel methods. Despite technical advances in the latest years, these methods are not able to cope with large drug-target interaction spaces and to integrate multiple sources of biological information.

Cell fate decisions and pluripotency, but also malignancy depend on networks of key transcriptional regulators. The T-box transcription factor TBX3 has been implicated in the regulation of embryonic stem cell self-renewal and cardiogenesis. We have recently discovered that forced TBX3 expression in embryonic stem cells promotes mesendoderm specification directly by activating key lineage specification factors and indirectly by enhancing paracrine NODAL signalling. Interestingly, aberrant TBX3 expression is associated with breast cancer and melanoma formation. In other cancers, loss of TBX3 expression is associated with a more aggressive phenotype e.g. in gastric and cervical cancer. The precise function of TBX3 in pancreatic ductal adenocarcinoma remains to be determined. In the current study we provide conclusive evidence for TBX3 overexpression in pancreatic cancer samples as compared to healthy tissue. While proliferation remains unaltered, forced TBX3 expression strongly increases migration and invasion, but also angiogenesis in vitro and in vivo. Finally, we describe the TBX3-dependency of cancer stem cells that perpetuate themselves through an autocrine TBX3-ACTIVIN/NODAL signalling loop to sustain stemness. Thus, TBX3 is a new key player among pluripotency-related genes driving cancer formation.

Mammalian cells have developed intricate mechanisms to interpret, integrate, and respond to extracellular stimuli. For example, tumor necrosis factor (TNF) rapidly activates proinflammatory genes, but our understanding of how this occurs against the ongoing transcriptional program of the cell is far from complete. Here, we monitor the early phase of this cascade at high spatiotemporal resolution in TNF-stimulated human endothelial cells. NF-κB, the transcription factor complex driving the response, interferes with the regulatory machinery by binding active enhancers already in interaction with gene promoters. Notably, >50% of these enhancers do not encode canonical NF-κB binding motifs. Using a combination of genomics tools, we find that binding site selection plays a key role in NF-κΒ-mediated transcriptional activation and repression. We demonstrate the latter by describing the synergy between NF-κΒ and the corepressor JDP2. Finally, detailed analysis of a 2.8-Mbp locus using sub-kbp-resolution targeted chromatin conformation capture and genome editing uncovers how NF-κΒ that has just entered the nucleus exploits pre-existing chromatin looping to exert its multimodal role. This work highlights the involvement of topology in cis-regulatory element function during acute transcriptional responses, where primary DNA sequence and its higher-order structure constitute a regulatory context leading to either gene activation or repression.

The protein kinase D isoenzymes PKD1/2/3 are prominent downstream targets of PKCs (Protein Kinase Cs) and phospholipase D in various biological systems. Recently, we identified PKD isoforms as novel mediators of tumour cell-endothelial cell communication, tumour cell motility and metastasis. Although PKD isoforms have been implicated in physiological/tumour angiogenesis, a role of PKDs during embryonic development, vasculogenesis and angiogenesis still remains elusive. We investigated the role of PKDs in germ layer segregation and subsequent vasculogenesis and angiogenesis using mouse embryonic stem cells (ESCs). We show that mouse ESCs predominantly express PKD2 followed by PKD3 while PKD1 displays negligible levels. Furthermore, we demonstrate that PKD2 is specifically phosphorylated/activated at the time of germ layer segregation. Time-restricted PKD2-activation limits mesendoderm formation and subsequent cardiovasculogenesis during early differentiation while leading to branching angiogenesis during late differentiation. In line, PKD2 loss-of-function analyses showed induction of mesendodermal differentiation in expense of the neuroectodermal germ layer. Our in vivo findings demonstrate that embryoid bodies transplanted on chicken chorioallantoic membrane induced an angiogenic response indicating that timed overexpression of PKD2 from day 4 onwards leads to augmented angiogenesis in differentiating ESCs. Taken together, our results describe novel and time-dependent facets of PKD2 during early cell fate determination.

Dendritic cells (DC) are professional antigen presenting cells that develop from hematopoietic stem cells through successive steps of lineage commitment and differentiation. Multipotent progenitors (MPP) are committed to DC restricted common DC progenitors (CDP), which differentiate into specific DC subsets, classical DC (cDC) and plasmacytoid DC (pDC). To determine epigenetic states and regulatory circuitries during DC differentiation, we measured consecutive changes of genome-wide gene expression, histone modification and transcription factor occupancy during the sequel MPP-CDP-cDC/pDC. Specific histone marks in CDP reveal a DC-primed epigenetic signature, which is maintained and reinforced during DC differentiation. Epigenetic marks and transcription factor PU.1 occupancy increasingly coincide upon DC differentiation. By integrating PU.1 occupancy and gene expression we devised a transcription factor regulatory circuitry for DC commitment and subset specification. The circuitry provides the transcription factor hierarchy that drives the sequel MPP-CDP-cDC/pDC, including Irf4, Irf8, Tcf4, Spib and Stat factors. The circuitry also includes feedback loops inferred for individual or multiple factors, which stabilize distinct stages of DC development and DC subsets. In summary, here we describe the basic regulatory circuitry of transcription factors that drives DC development.

Human protein-coding genes are often accompanied by divergently transcribed non-coding RNAs whose functions, especially in cell fate decisions, are poorly understood. Using an hESC-based cardiac differentiation model, we define a class of divergent lncRNAs, termed yin yang lncRNAs (yylncRNAs), that mirror the cell-type-specific expression pattern of their protein-coding counterparts. yylncRNAs are preferentially encoded from the genomic loci of key developmental cell fate regulators. Most yylncRNAs are spliced polyadenylated transcripts showing comparable expression patterns in vivo in mouse and in human embryos. Signifying their developmental function, the key mesoderm specifier BRACHYURY (T) is accompanied by yylncT, which localizes to the active T locus during mesoderm commitment. yylncT binds the de novo DNA methyltransferase DNMT3B, and its transcript is required for activation of the T locus, with yylncT depletion specifically abolishing mesodermal commitment. Collectively, we report a lncRNA-mediated regulatory layer safeguarding embryonic cell fate transitions.

Processes like cellular senescence are characterized by complex events giving rise to heterogeneous cell populations. However, the early molecular events driving this cascade remain elusive. We hypothesized that senescence entry is triggered by an early disruption of the cells' three-dimensional (3D) genome organization. To test this, we combined Hi-C, single-cell and population transcriptomics, imaging, and in silico modeling of three distinct cells types entering senescence. Genes involved in DNA conformation maintenance are suppressed upon senescence entry across all cell types. We show that nuclear depletion of the abundant HMGB2 protein occurs early on the path to senescence and coincides with the dramatic spatial clustering of CTCF. Knocking down HMGB2 suffices for senescence-induced CTCF clustering and for loop reshuffling, while ectopically expressing HMGB2 rescues these effects. Our data suggest that HMGB2-mediated genomic reorganization constitutes a primer for the ensuing senescent program.

Fetuin-A / α2-Heremans-Schmid-glycoprotein (gene name Ahsg) is a systemic inhibitor of ectopic calcification. Due to its high affinity for calcium phosphate, fetuin-A is highly abundant in mineralized bone matrix. Foreshortened femora in fetuin-A-deficient Ahsg-/- mice indicated a role for fetuin-A in bone formation. We studied early postnatal bone development in fetuin-A-deficient mice and discovered that femora from Ahsg-/- mice exhibited severely displaced distal epiphyses and deformed growth plates, similar to the human disease slipped capital femoral epiphysis (SCFE). The growth plate slippage occurred in 70% of Ahsg-/- mice of both sexes around three weeks postnatal. At this time point, mice weaned and rapidly gained weight and mobility. Epiphysis slippage never occurred in wildtype and heterozygous Ahsg+/- mice. Homozygous fetuin-A-deficient Ahsg-/- mice and, to a lesser degree, heterozygous Ahsg+/- mice showed lesions separating the proliferative zone from the hypertrophic zone of the growth plate. The hypertrophic growth plate cartilage in long bones from Ahsg-/- mice was significantly elongated and V-shaped until three weeks of age and thus prior to the slippage. Genome-wide transcriptome analysis of laser-dissected distal femoral growth plates from 13-day-old Ahsg-/- mice revealed a JAK-STAT-mediated inflammatory response including a 550-fold induction of the chemokine Cxcl9. At this stage, vascularization of the elongated growth plates was impaired, which was visualized by immunofluorescence staining. Thus, fetuin-A-deficient mice may serve as a rodent model of growth plate pathologies including SCFE and inflammatory cartilage degradation.

Kidney failure is frequently observed during and after COVID-19, but it remains elusive whether this is a direct effect of the virus. Here, we report that SARS-CoV-2 directly infects kidney cells and is associated with increased tubule-interstitial kidney fibrosis in patient autopsy samples. To study direct effects of the virus on the kidney independent of systemic effects of COVID-19, we infected human-induced pluripotent stem-cell-derived kidney organoids with SARS-CoV-2. Single-cell RNA sequencing indicated injury and dedifferentiation of infected cells with activation of profibrotic signaling pathways. Importantly, SARS-CoV-2 infection also led to increased collagen 1 protein expression in organoids. A SARS-CoV-2 protease inhibitor was able to ameliorate the infection of kidney cells by SARS-CoV-2. Our results suggest that SARS-CoV-2 can directly infect kidney cells and induce cell injury with subsequent fibrosis. These data could explain both acute kidney injury in COVID-19 patients and the development of chronic kidney disease in long COVID.

Genomic and epigenomic studies revealed dysregulation of long non-coding RNAs in many cancer entities, including liver cancer. We identified an epigenetic mechanism leading to upregulation of the long intergenic non-coding RNA 152 (LINC00152) expression in human hepatocellular carcinoma (HCC). Here, we aimed to characterize a potential competing endogenous RNA (ceRNA) network, in which LINC00152 exerts oncogenic functions by sponging miRNAs, thereby affecting their target gene expression. Database and gene expression data of human HCC were integrated to develop a potential LINC00152-driven ceRNA in silico. RNA immunoprecipitation and luciferase assay were used to identify miRNA binding to LINC00152 in human HCC cells. Functionally active players in the ceRNA network were analyzed using gene editing, siRNA or miRNA mimic transfection, and expression vectors in vitro. RNA expression in human HCC in vivo was validated by RNA in situ hybridization. Let-7c-5p, miR-23a-3p, miR-125a-5p, miR-125b-5p, miR-143a-3p, miR-193-3p, and miR-195-5p were detected as new components of the potential LINC00152 ceRNA network in human HCC. LINC00152 was confirmed to sponge miR143a-3p in human HCC cell lines, thereby limiting its binding to their respective target genes, like KLC2. KLC2 was identified as a central mediator promoting pro-tumorigenic effects of LINC00152 overexpression in HCC cells. Furthermore, co-expression of LINC00152 and KLC2 was observed in human HCC cohorts and high KLC2 expression was associated with shorter patient survival. Functional assays demonstrated that KLC2 promoted cell proliferation, clonogenicity and migration in vitro. The LINC00152-miR-143a-3p-KLC2 axis may represent a therapeutic target in human HCC.

Long non-coding RNAs (lncRNAs) contribute to cardiac (patho)physiology. Aging is the major risk factor for cardiovascular disease with cardiomyocyte apoptosis as one underlying cause. Here, we report the identification of the aging-regulated lncRNA Sarrah (ENSMUST00000140003) that is anti-apoptotic in cardiomyocytes. Importantly, loss of SARRAH (OXCT1-AS1) in human engineered heart tissue results in impaired contractile force development. SARRAH directly binds to the promoters of genes downregulated after SARRAH silencing via RNA-DNA triple helix formation and cardiomyocytes lacking the triple helix forming domain of Sarrah show an increase in apoptosis. One of the direct SARRAH targets is NRF2, and restoration of NRF2 levels after SARRAH silencing partially rescues the reduction in cell viability. Overexpression of Sarrah in mice shows better recovery of cardiac contractile function after AMI compared to control mice. In summary, we identified the anti-apoptotic evolutionary conserved lncRNA Sarrah, which is downregulated by aging, as a regulator of cardiomyocyte survival.

The Eph receptor tyrosine kinases and their respective ephrin-ligands are an important family of membrane receptors, being involved in developmental processes such as proliferation, migration, and in the formation of brain cancer such as glioma. Intracellular signaling pathways, which are activated by Eph receptor signaling, are well characterized. In contrast, it is unknown so far whether ephrins modulate the expression of lncRNAs, which would enable the transduction of environmental stimuli into our genome through a great gene regulatory spectrum. Applying a combination of functional in vitro assays, RNA sequencing, and qPCR analysis, we found that the proliferation and migration promoting stimulation of mouse cerebellar granule cells (CB) with ephrinA5 diminishes the expression of the cancer-related lncRNA Snhg15. In a human medulloblastoma cell line (DAOY) ephrinA5 stimulation similarly reduced SNHG15 expression. Computational analysis identified triple-helix-mediated DNA-binding sites of Snhg15 in promoters of genes found up-regulated upon ephrinA5 stimulation and known to be involved in tumorigenic processes. Our findings propose a crucial role of Snhg15 downstream of ephrinA5-induced signaling in regulating gene transcription in the nucleus. These findings could be potentially relevant for the regulation of tumorigenic processes in the context of glioma.

Chronic inflammation in the liver provokes fibrosis and, on long-term, carcinogenesis. This sequence is prototypically recapitulated in mice with hepatocyte-specific knock-out of the NF-κB essential modulator (NEMO), termed NEMOLPC-KO mice, in which increased hepatocyte apoptosis and compensatory regeneration cause steatosis, inflammation and fibrosis. Natural killer T (NKT) cells carrying the chemokine receptor CXCR6 participate in liver inflammation and injury responses. Here, we investigated the role of CXCR6 in the NEMOLPC-KO mouse model. Unexpectedly, genetic deletion of CXCR6 enhanced hepatocyte death, inflammation and fibrosis in NEMOLPC-KO mice. Although CXCR6 expression is restricted to immune cells in the liver, the adoptive transfer of CXCR6+ cells did not protect NEMOLPC-KOCxcr6-/- mice from hepatic injury. Gene array analyses revealed up-regulated stress response and metabolism pathways in hepatocytes from NEMOLPC-KOCxcr6-/- mice, functionally corresponding to an increased susceptibility of these hepatocytes to TNFα-induced cell death in vitro. These data revealed a novel CXCR6-dependent mechanism of suppressing inflammatory hepatocytic responses to cellular stress.

Long non-coding RNAs (lncRNAs) can act as scaffolds that promote the interaction of proteins, RNA, and DNA. There is increasing evidence of sequence-specific interactions of lncRNAs with DNA via triple-helix (triplex) formation. This process allows lncRNAs to recruit protein complexes to specific genomic regions and regulate gene expression. Here we propose a computational method called Triplex Domain Finder (TDF) to detect triplexes and characterize DNA-binding domains and DNA targets statistically. Case studies showed that this approach can detect the known domains of lncRNAs Fendrr, HOTAIR and MEG3. Moreover, we validated a novel DNA-binding domain in MEG3 by a genome-wide sequencing method. We used TDF to perform a systematic analysis of the triplex-forming potential of lncRNAs relevant to human cardiac differentiation. We demonstrated that the lncRNA with the highest triplex-forming potential, GATA6-AS, forms triple helices in the promoter of genes relevant to cardiac development. Moreover, down-regulation of GATA6-AS impairs GATA6 expression and cardiac development. These data indicate the unique ability of our computational tool to identify novel triplex-forming lncRNAs and their target genes.

Transposase-Accessible Chromatin followed by sequencing (ATAC-seq) is a simple protocol for detection of open chromatin. Computational footprinting, the search for regions with depletion of cleavage events due to transcription factor binding, is poorly understood for ATAC-seq. We propose the first footprinting method considering ATAC-seq protocol artifacts. HINT-ATAC uses a position dependency model to learn the cleavage preferences of the transposase. We observe strand-specific cleavage patterns around transcription factor binding sites, which are determined by local nucleosome architecture. By incorporating all these biases, HINT-ATAC is able to significantly outperform competing methods in the prediction of transcription factor binding sites with footprints.

Primary cells enter replicative senescence after a limited number of cell divisions. This process needs to be considered in cell culture experiments, and it is particularly important for regenerative medicine. Replicative senescence is associated with reproducible changes in DNA methylation (DNAm) at specific sites in the genome. The mechanism that drives senescence-associated DNAm changes remains unknown - it may involve stochastic DNAm drift due to imperfect maintenance of epigenetic marks or it is directly regulated at specific sites in the genome.

The relevance of topographic cues for commitment of induced pluripotent stem cells (iPSCs) is largely unknown. In this study, we demonstrate that groove-ridge structures with a periodicity in the submicrometer range induce elongation of iPSC colonies, guide the orientation of apical actin fibers, and direct the polarity of cell division. Elongation of iPSC colonies impacts also on their intrinsic molecular patterning, which seems to be orchestrated from the rim of the colonies. BMP4-induced differentiation is enhanced in elongated colonies, and the submicron grooves impact on the spatial modulation of YAP activity upon induction with this morphogen. Interestingly, TAZ, a YAP paralog, shows distinct cytoskeletal localization in iPSCs. These findings demonstrate that topography can guide orientation and organization of iPSC colonies, which may affect the interaction between mechanosensors and mechanotransducers in iPSCs.

Nephrotic syndrome (NS) is characterized by severe proteinuria as a consequence of kidney glomerular injury due to podocyte damage. In vitro models mimicking in vivo podocyte characteristics are a prerequisite to resolve NS pathogenesis. The detailed characterization of organoid podocytes resulting from a hybrid culture protocol showed a podocyte population that resembles adult podocytes and was superior compared with 2D counterparts, based on single-cell RNA sequencing, super-resolution imaging and electron microscopy. In this study, these next-generation podocytes in kidney organoids enabled personalized idiopathic nephrotic syndrome modeling, as shown by activated slit diaphragm signaling and podocyte injury following protamine sulfate, puromycin aminonucleoside treatment and exposure to NS plasma containing pathogenic permeability factors. Organoids cultured from cells of a patient with heterozygous NPHS2 mutations showed poor NPHS2 expression and aberrant NPHS1 localization, which was reversible after genetic correction. Repaired organoids displayed increased VEGFA pathway activity and transcription factor activity known to be essential for podocyte physiology, as shown by RNA sequencing. This study shows that organoids are the preferred model of choice to study idiopathic and congenital podocytopathies.

The advent of multi-modal single-cell sequencing techniques have shed new light on molecular mechanisms by simultaneously inspecting transcriptomes, epigenomes and proteomes of the same cell. However, to date, the existing computational approaches for integration of multimodal single-cell data are either computationally expensive, require the delineation of parameters or can only be applied to particular modalities.