Blimp-1 and c-Maf regulate Il10 and negatively regulate common and unique proinflammatory gene networks in IL-12 plus IL-27-driven T helper-1 cells.

Background: CD4 + Th1 cells producing IFN-γ are required to eradicate intracellular pathogens, however if uncontrolled these cells can cause immunopathology. The cytokine IL-10 is produced by multiple immune cells including Th1 cells during infection and regulates the immune response to minimise collateral host damage. In this study we aimed to elucidate the transcriptional network of genes controlling the expression of Il10 and proinflammatory cytokines, including Ifng in Th1 cells differentiated from mouse naive CD4 + T cells. Methods: We applied computational analysis of gene regulation derived from temporal profiling of gene expression clusters obtained from bulk RNA sequencing (RNA-seq) of flow cytometry sorted naïve CD4 + T cells from mouse spleens differentiated in vitro into Th1 effector cells with IL-12 and IL-27 to produce Ifng and Il10, compared to IL-27 alone which express Il10 only , or IL-12 alone which express Ifng and no Il10, or medium control driven-CD4 + T cells which do not express effector cytokines . Data were integrated with analysis of active genomic regions from these T cells using an assay for transposase-accessible chromatin with sequencing (ATAC)-seq, integrated with literature derived-Chromatin-immunoprecipitation (ChIP)-seq data and the RNA-seq data, to elucidate the transcriptional network of genes controlling expression of Il10 and pro-inflammatory effector genes in Th1 cells. The co-dominant role for the transcription factors, Prdm1 (encoding Blimp-1) and Maf (encoding c-Maf) , in cytokine gene regulation in Th1 cells, was confirmed using T cells obtained from mice with T-cell specific deletion of these transcription factors. Results: We show that the transcription factors Blimp-1 and c-Maf each have unique and common effects on cytokine gene regulation and not only co-operate to induce Il10 gene expression in IL-12 plus IL-27 differentiated mouse Th1 cells, but additionally directly negatively regulate key proinflammatory cytokines including Ifng, thus providing mechanisms for reinforcement of regulated Th1 cell responses. Conclusions: These data show that Blimp-1 and c-Maf positively and negatively regulate a network of both unique and common anti-inflammatory and pro-inflammatory genes to reinforce a Th1 response in mice that will eradicate pathogens with minimum immunopathology.

The predicted RNA-binding protein regulome of axonal mRNAs.

Neurons are morphologically complex cells that rely on the compartmentalization of protein expression to develop and maintain their cytoarchitecture. The targeting of RNA transcripts to axons is one of the mechanisms that allows rapid local translation of proteins in response to extracellular signals. 3' Untranslated regions (UTRs) of mRNA are noncoding sequences that play a critical role in determining transcript localization and translation by interacting with specific RNA-binding proteins (RBPs). However, how 3' UTRs contribute to mRNA metabolism and the nature of RBP complexes responsible for these functions remains elusive. We performed 3' end sequencing of RNA isolated from cell bodies and axons of sympathetic neurons exposed to either nerve growth factor (NGF) or neurotrophin 3 (NTF3, also known as NT-3). NGF and NTF3 are growth factors essential for sympathetic neuron development through distinct signaling mechanisms. Whereas NTF3 acts mostly locally, NGF signal is retrogradely transported from axons to cell bodies. We discovered that both NGF and NTF3 affect transcription and alternative polyadenylation in the nucleus and induce the localization of specific 3' UTR isoforms to axons, including short 3' UTR isoforms found exclusively in axons. The integration of our data with CLIP sequencing data supports a model whereby long 3' UTR isoforms associate with RBP complexes in the nucleus and, upon reaching the axons, are remodeled locally into shorter isoforms. Our findings shed new light into the complex relationship between nuclear polyadenylation, mRNA localization, and local 3' UTR remodeling in developing neurons.

Integrated transcriptome landscape of ALS identifies genome instability linked to TDP-43 pathology.

Amyotrophic Lateral Sclerosis (ALS) causes motor neuron degeneration, with 97% of cases exhibiting TDP-43 proteinopathy. Elucidating pathomechanisms has been hampered by disease heterogeneity and difficulties accessing motor neurons. Human induced pluripotent stem cell-derived motor neurons (iPSMNs) offer a solution; however, studies have typically been limited to underpowered cohorts. Here, we present a comprehensive compendium of 429 iPSMNs from 15 datasets, and 271 post-mortem spinal cord samples. Using reproducible bioinformatic workflows, we identify robust upregulation of p53 signalling in ALS in both iPSMNs and post-mortem spinal cord. p53 activation is greatest with C9orf72 repeat expansions but is weakest with SOD1 and FUS mutations. TDP-43 depletion potentiates p53 activation in both post-mortem neuronal nuclei and cell culture, thereby functionally linking p53 activation with TDP-43 depletion. ALS iPSMNs and post-mortem tissue display enrichment of splicing alterations, somatic mutations, and gene fusions, possibly contributing to the DNA damage response.

Physiological intron retaining transcripts in the cytoplasm abound during human motor neurogenesis.

Intron retention (IR) is now recognized as a dominant splicing event during motor neuron (MN) development; however, the role and regulation of intron-retaining transcripts (IRTs) localized to the cytoplasm remain particularly understudied. Here we show that IR is a physiological process that is spatiotemporally regulated during MN lineage restriction and that IRTs in the cytoplasm are detected in as many as 13% (n = 2297) of the genes expressed during this process. We identify a major class of cytoplasmic IRTs that are not associated with reduced expression of their own genes but instead show a high capacity for RNA-binding protein and miRNA occupancy. Finally, we show that ALS-causing VCP mutations lead to a selective increase in cytoplasmic abundance of this particular class of IRTs, which in turn temporally coincides with an increase in the nuclear expression level of predicted miRNA target genes. Altogether, our study identifies a previously unrecognized class of cytoplasmic intronic sequences with potential regulatory function beyond gene expression.

Image-based deep learning reveals the responses of human motor neurons to stress and VCP-related ALS.

AIMS: Although morphological attributes of cells and their substructures are recognised readouts of physiological or pathophysiological states, these have been relatively understudied in amyotrophic lateral sclerosis (ALS) research. METHODS: In this study, we integrate multichannel fluorescence high-content microscopy data with deep learning imaging methods to reveal-directly from unsegmented images-novel neurite-associated morphological perturbations associated with (ALS-causing) VCP-mutant human motor neurons (MNs). RESULTS: Surprisingly, we reveal that previously unrecognised disease-relevant information is withheld in broadly used and often considered 'generic' biological markers of nuclei (DAPI) and neurons ( ß III-tubulin). Additionally, we identify changes within the information content of ALS-related RNA binding protein (RBP) immunofluorescence imaging that is captured in VCP-mutant MN cultures. Furthermore, by analysing MN cultures exposed to different extrinsic stressors, we show that heat stress recapitulates key aspects of ALS. CONCLUSIONS: Our study therefore reveals disease-relevant information contained in a range of both generic and more specific fluorescent markers and establishes the use of image-based deep learning methods for rapid, automated and unbiased identification of biological hypotheses.

Aberrant cytoplasmic intron retention is a blueprint for RNA binding protein mislocalization in VCP-related amyotrophic lateral sclerosis.

We recently described aberrantly increased cytoplasmic SFPQ intron-retaining transcripts (IRTs) and concurrent SFPQ protein mislocalization as new hallmarks of amyotrophic lateral sclerosis (ALS). However, the generalizability and potential roles of cytoplasmic IRTs in health and disease remain unclear. Here, using time-resolved deep sequencing of nuclear and cytoplasmic fractions of human induced pluripotent stem cells undergoing motor neurogenesis, we reveal that ALS-causing VCP gene mutations lead to compartment-specific aberrant accumulation of IRTs. Specifically, we identify >100 IRTs with increased cytoplasmic abundance in ALS samples. Furthermore, these aberrant cytoplasmic IRTs possess sequence-specific attributes and differential predicted binding affinity to RNA binding proteins. Remarkably, TDP-43, SFPQ and FUS-RNA binding proteins known for nuclear-to-cytoplasmic mislocalization in ALS-abundantly and specifically bind to this aberrant cytoplasmic pool of IRTs. Our data are therefore consistent with a novel role for cytoplasmic IRTs in regulating compartment-specific protein abundance. This study provides new molecular insight into potential pathomechanisms underlying ALS and highlights aberrant cytoplasmic IRTs as potential therapeutic targets.

Cytoplasmic cleavage of IMPA1 3' UTR is necessary for maintaining axon integrity.

The 3' untranslated regions (3' UTRs) of messenger RNAs (mRNAs) are non-coding sequences involved in many aspects of mRNA metabolism, including intracellular localization and translation. Incorrect processing and delivery of mRNA cause severe developmental defects and have been implicated in many neurological disorders. Here, we use deep sequencing to show that in sympathetic neuron axons, the 3' UTRs of many transcripts undergo cleavage, generating isoforms that express the coding sequence with a short 3' UTR and stable 3' UTR-derived fragments of unknown function. Cleavage of the long 3' UTR of Inositol Monophosphatase 1 (IMPA1) mediated by a protein complex containing the endonuclease argonaute 2 (Ago2) generates a translatable isoform that is necessary for maintaining the integrity of sympathetic neuron axons. Thus, our study provides a mechanism of mRNA metabolism that simultaneously regulates local protein synthesis and generates an additional class of 3' UTR-derived RNAs.

Automated and unbiased discrimination of ALS from control tissue at single cell resolution.

Histopathological analysis of tissue sections is invaluable in neurodegeneration research. However, cell-to-cell variation in both the presence and severity of a given phenotype is a key limitation of this approach, reducing the signal to noise ratio and leaving unresolved the potential of single-cell scoring for a given disease attribute. Here, we tested different machine learning methods to analyse high-content microscopy measurements of hundreds of motor neurons (MNs) from amyotrophic lateral sclerosis (ALS) post-mortem tissue sections. Furthermore, we automated the identification of phenotypically distinct MN subpopulations in VCP- and SOD1-mutant transgenic mice, revealing common morphological cellular phenotypes. Additionally we established scoring metrics to rank cells and tissue samples for both disease probability and severity. By adapting this paradigm to human post-mortem tissue, we validated our core finding that morphological descriptors robustly discriminate ALS from control healthy tissue at single cell resolution. Determining disease presence, severity and unbiased phenotypes at single cell resolution might prove transformational in our understanding of ALS and neurodegeneration more broadly.

Mammary epithelial morphogenesis in 3D combinatorial microenvironments.

Human mammary epithelial cells can proliferate and reorganize into polarized multi-cellular constructs in-vitro, thereby functioning as an important model system in recapitulating key steps of in-vivo morphogenesis. Current approaches to constructing such three-dimensional mimics of the in-vivo microenvironment have involved the use of complex and ill-defined naturally derived matrices, whose properties are difficult to manipulate independently, and which have therefore limited our ability to understand the extrinsic regulation of morphogenesis. Here, we employ an automated, high-throughput approach to array modular building blocks of synthetic components, and develop a systematic approach to analyze colonies resulting from these varied microenvironmental combinations. This methodology allows us to systematically map the relationship between microenvironmental properties and ensuing morphogenetic phenotypes. Our analysis reveals that apico-basal polarity of mammary epithelial cells occurs within a narrow range of matrix stiffness, and that phenotypic homogeneity is favored in matrices which are insensitive to MMP-mediated degradation. Furthermore, combinations of extracellular proteins in the matrix finely tune the morphology of the mammary colonies, suggesting that subtle disregulations of the microenvironment may play a significant role in pathological disease states. This approach, which leverages the combinatorial possibilities of modular synthetic artificial extracellular matrices with an automated technology platform, demonstrates how morphogenesis can be assessed systematically in 3D, and provides new insights into mammary epithelial multicellularity.

Widespread FUS mislocalization is a molecular hallmark of amyotrophic lateral sclerosis.

Mutations causing amyotrophic lateral sclerosis (ALS) clearly implicate ubiquitously expressed and predominantly nuclear RNA binding proteins, which form pathological cytoplasmic inclusions in this context. However, the possibility that wild-type RNA binding proteins mislocalize without necessarily becoming constituents of cytoplasmic inclusions themselves remains relatively unexplored. We hypothesized that nuclear-to-cytoplasmic mislocalization of the RNA binding protein fused in sarcoma (FUS), in an unaggregated state, may occur more widely in ALS than previously recognized. To address this hypothesis, we analysed motor neurons from a human ALS induced-pluripotent stem cell model caused by the VCP mutation. Additionally, we examined mouse transgenic models and post-mortem tissue from human sporadic ALS cases. We report nuclear-to-cytoplasmic mislocalization of FUS in both VCP-mutation related ALS and, crucially, in sporadic ALS spinal cord tissue from multiple cases. Furthermore, we provide evidence that FUS protein binds to an aberrantly retained intron within the SFPQ transcript, which is exported from the nucleus into the cytoplasm. Collectively, these data support a model for ALS pathogenesis whereby aberrant intron retention in SFPQ transcripts contributes to FUS mislocalization through their direct interaction and nuclear export. In summary, we report widespread mislocalization of the FUS protein in ALS and propose a putative underlying mechanism for this process.

Publisher Correction: c-Maf controls immune responses by regulating disease-specific gene networks and repressing IL-2 in CD4+ T cells.

In the version of this article initially published, the Supplementary Data file was an incorrect version. The correct version is now provided. The error has been corrected in the HTML and PDF version of the article.

Intron retention and nuclear loss of SFPQ are molecular hallmarks of ALS.

Mutations causing amyotrophic lateral sclerosis (ALS) strongly implicate ubiquitously expressed regulators of RNA processing. To understand the molecular impact of ALS-causing mutations on neuronal development and disease, we analysed transcriptomes during in vitro differentiation of motor neurons (MNs) from human control and patient-specific VCP mutant induced-pluripotent stem cells (iPSCs). We identify increased intron retention (IR) as a dominant feature of the splicing programme during early neural differentiation. Importantly, IR occurs prematurely in VCP mutant cultures compared with control counterparts. These aberrant IR events are also seen in independent RNAseq data sets from SOD1- and FUS-mutant MNs. The most significant IR is seen in the SFPQ transcript. The SFPQ protein binds extensively to its retained intron, exhibits lower nuclear abundance in VCP mutant cultures and is lost from nuclei of MNs in mouse models and human sporadic ALS. Collectively, we demonstrate SFPQ IR and nuclear loss as molecular hallmarks of familial and sporadic ALS.

c-Maf controls immune responses by regulating disease-specific gene networks and repressing IL-2 in CD4+ T cells.

The transcription factor c-Maf induces the anti-inflammatory cytokine IL-10 in CD4+ T cells in vitro. However, the global effects of c-Maf on diverse immune responses in vivo are unknown. Here we found that c-Maf regulated IL-10 production in CD4+ T cells in disease models involving the TH1 subset of helper T cells (malaria), TH2 cells (allergy) and TH17 cells (autoimmunity) in vivo. Although mice with c-Maf deficiency targeted to T cells showed greater pathology in TH1 and TH2 responses, TH17 cell-mediated pathology was reduced in this context, with an accompanying decrease in TH17 cells and increase in Foxp3+ regulatory T cells. Bivariate genomic footprinting elucidated the c-Maf transcription-factor network, including enhanced activity of NFAT; this led to the identification and validation of c-Maf as a negative regulator of IL-2. The decreased expression of the gene encoding the transcription factor RORγt (Rorc) that resulted from c-Maf deficiency was dependent on IL-2, which explained the in vivo observations. Thus, c-Maf is a positive and negative regulator of the expression of cytokine-encoding genes, with context-specific effects that allow each immune response to occur in a controlled yet effective manner.

Progressive Motor Neuron Pathology and the Role of Astrocytes in a Human Stem Cell Model of VCP-Related ALS.

Motor neurons (MNs) and astrocytes (ACs) are implicated in the pathogenesis of amyotrophic lateral sclerosis (ALS), but their interaction and the sequence of molecular events leading to MN death remain unresolved. Here, we optimized directed differentiation of induced pluripotent stem cells (iPSCs) into highly enriched (> 85%) functional populations of spinal cord MNs and ACs. We identify significantly increased cytoplasmic TDP-43 and ER stress as primary pathogenic events in patient-specific valosin-containing protein (VCP)-mutant MNs, with secondary mitochondrial dysfunction and oxidative stress. Cumulatively, these cellular stresses result in synaptic pathology and cell death in VCP-mutant MNs. We additionally identify a cell-autonomous VCP-mutant AC survival phenotype, which is not attributable to the same molecular pathology occurring in VCP-mutant MNs. Finally, through iterative co-culture experiments, we uncover non-cell-autonomous effects of VCP-mutant ACs on both control and mutant MNs. This work elucidates molecular events and cellular interplay that could guide future therapeutic strategies in ALS.

Phenobarbital induces cell cycle transcriptional responses in mouse liver humanized for constitutive androstane and pregnane x receptors.

The constitutive androstane receptor (CAR) and the pregnane X receptor (PXR) are closely related nuclear receptors involved in drug metabolism and play important roles in the mechanism of phenobarbital (PB)-induced rodent nongenotoxic hepatocarcinogenesis. Here, we have used a humanized CAR/PXR mouse model to examine potential species differences in receptor-dependent mechanisms underlying liver tissue molecular responses to PB. Early and late transcriptomic responses to sustained PB exposure were investigated in liver tissue from double knock-out CAR and PXR (CAR(KO)-PXR(KO)), double humanized CAR and PXR (CAR(h)-PXR(h)), and wild-type C57BL/6 mice. Wild-type and CAR(h)-PXR(h) mouse livers exhibited temporally and quantitatively similar transcriptional responses during 91 days of PB exposure including the sustained induction of the xenobiotic response gene Cyp2b10, the Wnt signaling inhibitor Wisp1, and noncoding RNA biomarkers from the Dlk1-Dio3 locus. Transient induction of DNA replication (Hells, Mcm6, and Esco2) and mitotic genes (Ccnb2, Cdc20, and Cdk1) and the proliferation-related nuclear antigen Mki67 were observed with peak expression occurring between 1 and 7 days PB exposure. All these transcriptional responses were absent in CAR(KO)-PXR(KO) mouse livers and largely reversible in wild-type and CAR(h)-PXR(h) mouse livers following 91 days of PB exposure and a subsequent 4-week recovery period. Furthermore, PB-mediated upregulation of the noncoding RNA Meg3, which has recently been associated with cellular pluripotency, exhibited a similar dose response and perivenous hepatocyte-specific localization in both wild-type and CAR(h)-PXR(h) mice. Thus, mouse livers coexpressing human CAR and PXR support both the xenobiotic metabolizing and the proliferative transcriptional responses following exposure to PB.

Ha-ras and ß-catenin oncoproteins orchestrate metabolic programs in mouse liver tumors.

The process of hepatocarcinogenesis in the diethylnitrosamine (DEN) initiation/phenobarbital (PB) promotion mouse model involves the selective clonal outgrowth of cells harboring oncogene mutations in Ctnnb1, while spontaneous or DEN-only-induced tumors are often Ha-ras- or B-raf-mutated. The molecular mechanisms and pathways underlying these different tumor sub-types are not well characterized. Their identification may help identify markers for xenobiotic promoted versus spontaneously occurring liver tumors. Here, we have characterized mouse liver tumors harboring either Ctnnb1 or Ha-ras mutations via integrated molecular profiling at the transcriptional, translational and post-translational levels. In addition, metabolites of the intermediary metabolism were quantified by high resolution (1)H magic angle nuclear magnetic resonance. We have identified tumor genotype-specific differences in mRNA and miRNA expression, protein levels, post-translational modifications, and metabolite levels that facilitate the molecular and biochemical stratification of tumor phenotypes. Bioinformatic integration of these data at the pathway level led to novel insights into tumor genotype-specific aberrant cell signaling and in particular to a better understanding of alterations in pathways of the cell intermediary metabolism, which are driven by the constitutive activation of the ß-Catenin and Ha-ras oncoproteins in tumors of the two genotypes.

Computational modeling identifies key gene regulatory interactions underlying phenobarbital-mediated tumor promotion.

Gene regulatory interactions underlying the early stages of non-genotoxic carcinogenesis are poorly understood. Here, we have identified key candidate regulators of phenobarbital (PB)-mediated mouse liver tumorigenesis, a well-characterized model of non-genotoxic carcinogenesis, by applying a new computational modeling approach to a comprehensive collection of in vivo gene expression studies. We have combined our previously developed motif activity response analysis (MARA), which models gene expression patterns in terms of computationally predicted transcription factor binding sites with singular value decomposition (SVD) of the inferred motif activities, to disentangle the roles that different transcriptional regulators play in specific biological pathways of tumor promotion. Furthermore, transgenic mouse models enabled us to identify which of these regulatory activities was downstream of constitutive androstane receptor and ß-catenin signaling, both crucial components of PB-mediated liver tumorigenesis. We propose novel roles for E2F and ZFP161 in PB-mediated hepatocyte proliferation and suggest that PB-mediated suppression of ESR1 activity contributes to the development of a tumor-prone environment. Our study shows that combining MARA with SVD allows for automated identification of independent transcription regulatory programs within a complex in vivo tissue environment and provides novel mechanistic insights into PB-mediated hepatocarcinogenesis.

Identification of Dlk1-Dio3 imprinted gene cluster noncoding RNAs as novel candidate biomarkers for liver tumor promotion.

The molecular events during nongenotoxic carcinogenesis and their temporal order are poorly understood but thought to include long-lasting perturbations of gene expression. Here, we have investigated the temporal sequence of molecular and pathological perturbations at early stages of phenobarbital (PB) mediated liver tumor promotion in vivo. Molecular profiling (mRNA, microRNA [miRNA], DNA methylation, and proteins) of mouse liver during 13 weeks of PB treatment revealed progressive increases in hepatic expression of long noncoding RNAs and miRNAs originating from the Dlk1-Dio3 imprinted gene cluster, a locus that has recently been associated with stem cell pluripotency in mice and various neoplasms in humans. PB induction of the Dlk1-Dio3 cluster noncoding RNA (ncRNA) Meg3 was localized to glutamine synthetase-positive hypertrophic perivenous hepatocytes, suggesting a role for ß-catenin signaling in the dysregulation of Dlk1-Dio3 ncRNAs. The carcinogenic relevance of Dlk1-Dio3 locus ncRNA induction was further supported by in vivo genetic dependence on constitutive androstane receptor and ß-catenin pathways. Our data identify Dlk1-Dio3 ncRNAs as novel candidate early biomarkers for mouse liver tumor promotion and provide new opportunities for assessing the carcinogenic potential of novel compounds.