CRISPR/Cas9 screen in human iPSC-derived cortical neurons identifies NEK6 as a novel disease modifier of C9orf72 poly(PR) toxicity.

INTRODUCTION: The most common genetic cause of frontotemporal dementia (FTD) and amyotrophic lateral sclerosis (ALS) are hexanucleotide repeats in chromosome 9 open reading frame 72 (C9orf72). These repeats produce dipeptide repeat proteins with poly(PR) being the most toxic one. METHODS: We performed a kinome-wide CRISPR/Cas9 knock-out screen in human induced pluripotent stem cell (iPSC) -derived cortical neurons to identify modifiers of poly(PR) toxicity, and validated the role of candidate modifiers using in vitro, in vivo, and ex-vivo studies. RESULTS: Knock-down of NIMA-related kinase 6 (NEK6) prevented neuronal toxicity caused by poly(PR). Knock-down of nek6 also ameliorated the poly(PR)-induced axonopathy in zebrafish and NEK6 was aberrantly expressed in C9orf72 patients. Suppression of NEK6 expression and NEK6 activity inhibition rescued axonal transport defects in cortical neurons from C9orf72 patient iPSCs, at least partially by reversing p53-related DNA damage. DISCUSSION: We identified NEK6, which regulates poly(PR)-mediated p53-related DNA damage, as a novel therapeutic target for C9orf72 FTD/ALS.

Independence of repressive histone marks and chromatin compaction during senescent heterochromatic layer formation.

The expansion of repressive epigenetic marks has been implicated in heterochromatin formation during embryonic development, but the general applicability of this mechanism is unclear. Here we show that nuclear rearrangement of repressive histone marks H3K9me3 and H3K27me3 into nonoverlapping structural layers characterizes senescence-associated heterochromatic foci (SAHF) formation in human fibroblasts. However, the global landscape of these repressive marks remains unchanged upon SAHF formation, suggesting that in somatic cells, heterochromatin can be formed through the spatial repositioning of pre-existing repressively marked histones. This model is reinforced by the correlation of presenescent replication timing with both the subsequent layered structure of SAHFs and the global landscape of the repressive marks, allowing us to integrate microscopic and genomic information. Furthermore, modulation of SAHF structure does not affect the occupancy of these repressive marks, nor vice versa. These experiments reveal that high-order heterochromatin formation and epigenetic remodeling of the genome can be discrete events.

Prognostic relevance of molecular subtypes and master regulators in pancreatic ductal adenocarcinoma.

BACKGROUND: Pancreatic cancer is poorly characterized at genetic and non-genetic levels. The current study evaluates in a large cohort of patients the prognostic relevance of molecular subtypes and key transcription factors in pancreatic ductal adenocarcinoma (PDAC). METHODS: We performed gene expression analysis of whole-tumor tissue obtained from 118 surgically resected PDAC and 13 histologically normal pancreatic tissue samples. Cox regression models were used to study the effect on survival of molecular subtypes and 16 clinicopathological prognostic factors. In order to better understand the biology of PDAC we used iRegulon to identify transcription factors (TFs) as master regulators of PDAC and its subtypes. RESULTS: We confirmed the PDAssign gene signature as classifier of PDAC in molecular subtypes with prognostic relevance. We found molecular subtypes, but not clinicopathological factors, as independent predictors of survival. Regulatory network analysis predicted that HNF1A/B are among thousand TFs the top enriched master regulators of the genes expressed in the normal pancreatic tissue compared to the PDAC regulatory network. On immunohistochemistry staining of PDAC samples, we observed low expression of HNF1B in well differentiated towards no expression in poorly differentiated PDAC samples. We predicted IRF/STAT, AP-1, and ETS-family members as key transcription factors in gene signatures downstream of mutated KRAS. CONCLUSIONS: PDAC can be classified in molecular subtypes that independently predict survival. HNF1A/B seem to be good candidates as master regulators of pancreatic differentiation, which at the protein level loses its expression in malignant ductal cells of the pancreas, suggesting its putative role as tumor suppressor in pancreatic cancer. TRIAL REGISTRATION: The study was registered at ClinicalTrials.gov under the number NCT01116791 (May 3, 2010).

iRegulon: from a gene list to a gene regulatory network using large motif and track collections.

Identifying master regulators of biological processes and mapping their downstream gene networks are key challenges in systems biology. We developed a computational method, called iRegulon, to reverse-engineer the transcriptional regulatory network underlying a co-expressed gene set using cis-regulatory sequence analysis. iRegulon implements a genome-wide ranking-and-recovery approach to detect enriched transcription factor motifs and their optimal sets of direct targets. We increase the accuracy of network inference by using very large motif collections of up to ten thousand position weight matrices collected from various species, and linking these to candidate human TFs via a motif2TF procedure. We validate iRegulon on gene sets derived from ENCODE ChIP-seq data with increasing levels of noise, and we compare iRegulon with existing motif discovery methods. Next, we use iRegulon on more challenging types of gene lists, including microRNA target sets, protein-protein interaction networks, and genetic perturbation data. In particular, we over-activate p53 in breast cancer cells, followed by RNA-seq and ChIP-seq, and could identify an extensive up-regulated network controlled directly by p53. Similarly we map a repressive network with no indication of direct p53 regulation but rather an indirect effect via E2F and NFY. Finally, we generalize our computational framework to include regulatory tracks such as ChIP-seq data and show how motif and track discovery can be combined to map functional regulatory interactions among co-expressed genes. iRegulon is available as a Cytoscape plugin from http://iregulon.aertslab.org.

Unraveling networks of co-regulated genes on the sole basis of genome sequences.

With the growing number of available microbial genome sequences, regulatory signals can now be revealed as conserved motifs in promoters of orthologous genes (phylogenetic footprints). A next challenge is to unravel genome-scale regulatory networks. Using as sole input genome sequences, we predicted cis-regulatory elements for each gene of the yeast Saccharomyces cerevisiae by discovering over-represented motifs in the promoters of their orthologs in 19 Saccharomycetes species. We then linked all genes displaying similar motifs in their promoter regions and inferred a co-regulation network including 56,919 links between 3171 genes. Comparison with annotated regulons highlights the high predictive value of the method: a majority of the top-scoring predictions correspond to already known co-regulations. We also show that this inferred network is as accurate as a co-expression network built from hundreds of transcriptome microarray experiments. Furthermore, we experimentally validated 14 among 16 new functional links between orphan genes and known regulons. This approach can be readily applied to unravel gene regulatory networks from hundreds of microbial genomes for which no other information is available except the sequence. Long-term benefits can easily be perceived when considering the exponential increase of new genome sequences.

C. elegans ORFeome version 1.1: experimental verification of the genome annotation and resource for proteome-scale protein expression.

To verify the genome annotation and to create a resource to functionally characterize the proteome, we attempted to Gateway-clone all predicted protein-encoding open reading frames (ORFs), or the 'ORFeome,' of Caenorhabditis elegans. We successfully cloned approximately 12,000 ORFs (ORFeome 1.1), of which roughly 4,000 correspond to genes that are untouched by any cDNA or expressed-sequence tag (EST). More than 50% of predicted genes needed corrections in their intron-exon structures. Notably, approximately 11,000 C. elegans proteins can now be expressed under many conditions and characterized using various high-throughput strategies, including large-scale interactome mapping. We suggest that similar ORFeome projects will be valuable for other organisms, including humans.

FUS-ALS hiPSC-derived astrocytes impair human motor units through both gain-of-toxicity and loss-of-support mechanisms.

BACKGROUND: Astrocytes play a crucial, yet not fully elucidated role in the selective motor neuron pathology in amyotrophic lateral sclerosis (ALS). Among other responsibilities, astrocytes provide important neuronal homeostatic support, however this function is highly compromised in ALS. The establishment of fully human coculture systems can be used to further study the underlying mechanisms of the dysfunctional intercellular interplay, and has the potential to provide a platform for revealing novel therapeutic entry points. METHODS: In this study, we characterised human induced pluripotent stem cell (hiPSC)-derived astrocytes from FUS-ALS patients, and incorporated these cells into a human motor unit microfluidics model to investigate the astrocytic effect on hiPSC-derived motor neuron network and functional neuromuscular junctions (NMJs) using immunocytochemistry and live-cell recordings. FUS-ALS cocultures were systematically compared to their CRISPR-Cas9 gene-edited isogenic control systems. RESULTS: We observed a dysregulation of astrocyte homeostasis, which resulted in a FUS-ALS-mediated increase in reactivity and secretion of inflammatory cytokines. Upon coculture with motor neurons and myotubes, we detected a cytotoxic effect on motor neuron-neurite outgrowth, NMJ formation and functionality, which was improved or fully rescued by isogenic control astrocytes. We demonstrate that ALS astrocytes have both a gain-of-toxicity and loss-of-support function involving the WNT/ß-catenin pathway, ultimately contributing to the disruption of motor neuron homeostasis, intercellular networks and NMJs. CONCLUSIONS: Our findings shine light on a complex, yet highly important role of astrocytes in ALS, and provides further insight in to their pathological mechanisms.

Single-cell Transcriptomics Uncover a Novel Role of Myeloid Cells and T-lymphocytes in the Fibrotic Microenvironment in Peyronie's Disease.

BACKGROUND: Peyronie's disease (PD) is an acquired fibrotic disease affecting the penile tunica albuginea that can lead to curvature and deformities, shortening, and erectile dysfunction. Although immunological mechanisms have been suggested for the pathophysiology of PD, these have not been investigated using single-cell transcriptomics. OBJECTIVE: To investigate the immunological signature of plaques from PD patients using immunohistochemistry (IHC) and single-cell RNA sequencing (scRNA-Seq). DESIGN, SETTING, AND PARTICIPANTS: Tunica albuginea biopsy was performed in patients undergoing penile surgery for either PD (n = 12) or plication or penile cancer (control, n = 6). The inclusion criteria for PD patients were stable chronic disease (≥12 mo in duration) and no previous penile surgery or intralesional injection therapy. OUTCOME MEASUREMENTS AND STATISTICAL ANALYSIS: IHC was performed on surgical samples from ten patients with PD and five control subjects. An additional two PD and one control sample were used for scRNA-Seq (droplet-based; 10X Genomics). Cell clusters were visualised using heatmaps and t-distributed stochastic neighbour embedding plots (BioTuring v2.7.5). RESULTS AND LIMITATIONS: IHC revealed the presence of myeloid dendritic cells (DCs; CD68+, TLR4+, CD206+), cytotoxic T lymphocytes (CTLs; CD3+, CD8+), and B lymphocytes (CD20+) in PD plaques, which were absent in controls. scRNA-Seq yielded results for 3312 PD and 5658 control cells. Cell clusters contained fibroblasts (COL1A2+), myofibroblasts (COL1A2+, ACTA2+), smooth muscle cells (ACTA2+, DES+), endothelial cells (VWF+), myeloid cells (CD14+), T lymphocytes (CD3D+), and neutrophils (ALPL+). Myeloid cell subclustering showed infiltration of monocyte-derived cells; control tissue contained classical DCs and resident macrophages. Lymphocyte subclustering revealed mucosal-associated invariant T (MAIT) cells and CTLs in PD. Differential gene expression suggests an increase in inflammatory and immune responses in chronic PD. The study is limited by the small scRNA-seq sample size (n = 3) for IHC, mitigated by a larger cohort of historic paraffin-embedded samples (n = 15), which showed largely parallel findings. Owing to tissue stiffness and extracellular matrix adhesion, our single-cell yield was lower for PD than for the control sample. CONCLUSIONS: Our data suggest that even in the chronic PD stage (painless and stable curvature) there is a sustained inflammatory reaction. While vascularisation and collagen production are elevated, the inflammation is driven by specialised monocyte-derived CTL and MAIT cells. These findings could uncover new avenues for medical treatment of PD. PATIENT SUMMARY: We looked at the role of the immune system in patients suffering from Peyronie's disease, a condition causing shortening and curvature of the penis. We found that even in a stable, chronic stage of the disease, there is activation of the immune system. Our results suggest that there is potential for novel treatments for this condition.

Unraveling the Molecular Basis of the Dystrophic Process in Limb-Girdle Muscular Dystrophy LGMD-R12 by Differential Gene Expression Profiles in Diseased and Healthy Muscles.

Limb-girdle muscular dystrophy R12 (LGMD-R12) is caused by two mutations in anoctamin-5 (ANO5). Our aim was to identify genes and pathways that underlie LGMD-R12 and explain differences in the molecular predisposition and susceptibility between three thigh muscles that are severely (semimembranosus), moderately (vastus lateralis) or mildly (rectus femoris) affected in this disease. We performed transcriptomics on these three muscles in 16 male LGMD-R12 patients and 15 age-matched male controls. Our results showed that LGMD-R12 dystrophic muscle is associated with the expression of genes indicative of fibroblast and adipocyte replacement, such as fibroadipogenic progenitors and immune cell infiltration, while muscle protein synthesis and metabolism were downregulated. Muscle degeneration was associated with an increase in genes involved in muscle injury and inflammation, and muscle repair/regeneration. Baseline differences between muscles in healthy individuals indicated that muscles that are the most affected by LGMD-R12 have the lowest expression of transcription factor networks involved in muscle (re)generation and satellite stem cell activation. Instead, they show relative high levels of fetal/embryonic myosins, all together indicating that muscles differ in their baseline regenerative potential. To conclude, we profiled the gene expression landscape in LGMD-R12, identified baseline differences in expression levels between differently affected muscles and characterized disease-associated changes.

The transcriptional program controlled by the stem cell leukemia gene Scl/Tal1 during early embryonic hematopoietic development.

The basic helix-loop-helix transcription factor Scl/Tal1 controls the development and subsequent differentiation of hematopoietic stem cells (HSCs). However, because few Scl target genes have been validated to date, the underlying mechanisms have remained largely unknown. In this study, we have used ChIP-Seq technology (coupling chromatin immunoprecipitation with deep sequencing) to generate a genome-wide catalog of Scl-binding events in a stem/progenitor cell line, followed by validation using primary fetal liver cells and comprehensive transgenic mouse assays. Transgenic analysis provided in vivo validation of multiple new direct Scl target genes and allowed us to reconstruct an in vivo validated network consisting of 17 factors and their respective regulatory elements. By coupling ChIP-Seq in model cell lines with in vivo transgenic validation and sophisticated bioinformatic analysis, we have identified a widely applicable strategy for the reconstruction of stem cell regulatory networks in which biologic material is otherwise limiting. Moreover, in addition to revealing multiple previously unrecognized links to known HSC regulators, as well as novel links to genes not previously implicated in HSC function, comprehensive transgenic analysis of regulatory elements provided substantial new insights into the transcriptional control of several important hematopoietic regulators, including Cbfa2t3h/Eto2, Cebpe, Nfe2, Zfpm1/Fog1, Erg, Mafk, Gfi1b, and Myb.