Time-course RNAseq data of murine AB1 mesothelioma and Renca renal cancer following immune checkpoint therapy.

Time-critical transcriptional events in the immune microenvironment are important for response to immune checkpoint blockade (ICB), yet these events are difficult to characterise and remain incompletely understood. Here, we present whole tumor RNA sequencing data in the context of treatment with ICB in murine models of AB1 mesothelioma and Renca renal cell cancer. We sequenced 144 bulk RNAseq samples from these two cancer types across 4 time points prior and after treatment with ICB. We also performed single-cell sequencing on 12 samples of AB1 and Renca tumors an hour before ICB administration. Our samples were equally distributed between responders and non-responders to treatment. Additionally, we sequenced AB1-HA mesothelioma tumors treated with two sample dissociation protocols to assess the impact of these protocols on the quality transcriptional information in our samples. These datasets provide time-course information to transcriptionally characterize the ICB response and provide detailed information at the single-cell level of the early tumor microenvironment prior to ICB therapy.

Immune checkpoint therapy responders display early clonal expansion of tumor infiltrating lymphocytes.

Immune checkpoint therapy (ICT) causes durable tumour responses in a subgroup of patients, but it is not well known how T cell receptor beta (TCRß) repertoire dynamics contribute to the therapeutic response. Using murine models that exclude variation in host genetics, environmental factors and tumour mutation burden, limiting variation between animals to naturally diverse TCRß repertoires, we applied TCRseq, single cell RNAseq and flow cytometry to study TCRß repertoire dynamics in ICT responders and non-responders. Increased oligoclonal expansion of TCRß clonotypes was observed in responding tumours. Machine learning identified TCRß CDR3 signatures unique to each tumour model, and signatures associated with ICT response at various timepoints before or during ICT. Clonally expanded CD8+ T cells in responding tumours post ICT displayed effector T cell gene signatures and phenotype. An early burst of clonal expansion during ICT is associated with response, and we report unique dynamics in TCRß signatures associated with ICT response.

CD4+ T cells drive an inflammatory, TNF-α/IFN-rich tumor microenvironment responsive to chemotherapy.

While chemotherapy remains the first-line treatment for many cancers, it is still unclear what distinguishes responders from non-responders. Here, we characterize the chemotherapy-responsive tumor microenvironment in mice, using RNA sequencing on tumors before and after cyclophosphamide, and compare the gene expression profiles of responders with progressors. Responsive tumors have an inflammatory and highly immune infiltrated pre-treatment tumor microenvironment characterized by the enrichment of pathways associated with CD4+ T cells, interferons (IFNs), and tumor necrosis factor alpha (TNF-α). The same gene expression profile is associated with response to cyclophosphamide-based chemotherapy in patients with breast cancer. Finally, we demonstrate that tumors can be sensitized to cyclophosphamide and 5-FU chemotherapy by pre-treatment with recombinant TNF-α, IFNγ, and poly(I:C). Thus, a CD4+ T cell-inflamed pre-treatment tumor microenvironment is necessary for response to chemotherapy, and this state can be therapeutically attained by targeted immunotherapy.

Temporally restricted activation of IFNß signaling underlies response to immune checkpoint therapy in mice.

The biological determinants of the response to immune checkpoint blockade (ICB) in cancer remain incompletely understood. Little is known about dynamic biological events that underpin therapeutic efficacy due to the inability to frequently sample tumours in patients. Here, we map the transcriptional profiles of 144 responding and non-responding tumours within two mouse models at four time points during ICB. We find that responding tumours display on/fast-off kinetics of type-I-interferon (IFN) signaling. Phenocopying of this kinetics using time-dependent sequential dosing of recombinant IFNs and neutralizing antibodies markedly improves ICB efficacy, but only when IFNß is targeted, not IFNα. We identify Ly6C+/CD11b+ inflammatory monocytes as the primary source of IFNß and find that active type-I-IFN signaling in tumour-infiltrating inflammatory monocytes is associated with T cell expansion in patients treated with ICB. Together, our results suggest that on/fast-off modulation of IFNß signaling is critical to the therapeutic response to ICB, which can be exploited to drive clinical outcomes towards response.

Characteristics of TCR Repertoire Associated With Successful Immune Checkpoint Therapy Responses.

Immunotherapies have revolutionized cancer treatment. In particular, immune checkpoint therapy (ICT) leads to durable responses in some patients with some cancers. However, the majority of treated patients do not respond. Understanding immune mechanisms that underlie responsiveness to ICT will help identify predictive biomarkers of response and develop treatments to convert non-responding patients to responding ones. ICT primarily acts at the level of adaptive immunity. The specificity of adaptive immune cells, such as T and B cells, is determined by antigen-specific receptors. T cell repertoires can be comprehensively profiled by high-throughput sequencing at the bulk and single-cell level. T cell receptor (TCR) sequencing allows for sensitive tracking of dynamic changes in antigen-specific T cells at the clonal level, giving unprecedented insight into the mechanisms by which ICT alters T cell responses. Here, we review how the repertoire influences response to ICT and conversely how ICT affects repertoire diversity. We will also explore how changes to the repertoire in different anatomical locations can better correlate and perhaps predict treatment outcome. We discuss the advantages and limitations of current metrics used to characterize and represent TCR repertoire diversity. Discovery of predictive biomarkers could lie in novel analysis approaches, such as network analysis of amino acids similarities between TCR sequences. Single-cell sequencing is a breakthrough technology that can link phenotype with specificity, identifying T cell clones that are crucial for successful ICT. The field of immuno-sequencing is rapidly developing and cross-disciplinary efforts are required to maximize the analysis, application, and validation of sequencing data. Unravelling the dynamic behavior of the TCR repertoire during ICT will be highly valuable for tracking and understanding anti-tumor immunity, biomarker discovery, and ultimately for the development of novel strategies to improve patient outcomes.

The bone marrow microenvironment of pre-B acute lymphoblastic leukemia at single-cell resolution.

The bone marrow microenvironment (BMM) plays a key role in leukemia progression, but its molecular complexity in pre-B cell acute lymphoblastic leukemia (B-ALL), the most common cancer in children, remains poorly understood. To gain further insight, we used single-cell RNA sequencing to characterize the kinetics of the murine BMM during B-ALL progression. Normal pro- and pre-B cells were found to be the most affected at the earliest stages of disease and this was associated with changes in expression of genes regulated by the AP1-transcription factor complex and regulatory factors NELFE, MYC and BCL11A. Granulocyte-macrophage progenitors show reduced expression of the tumor suppressor long non-coding RNA Neat1 and disruptions in the rate of transcription. Intercellular communication networks revealed monocyte-dendritic precursors to be consistently active during B-ALL progression, with enriched processes including cytokine-mediated signaling pathway, neutrophil-mediated immunity and regulation of cell migration and proliferation. In addition, we confirmed that the hematopoietic stem and progenitor cell compartment was perturbed during leukemogenesis. These findings extend our understanding of the complexity of changes and molecular interactions among the normal cells of the BMM during B-ALL progression.

Rhinovirus Infection Drives Complex Host Airway Molecular Responses in Children With Cystic Fibrosis.

Early-life viral infections are responsible for pulmonary exacerbations that can contribute to disease progression in young children with cystic fibrosis (CF). The most common respiratory viruses detected in the CF airway are human rhinoviruses (RV), and augmented airway inflammation in CF has been attributed to dysregulated airway epithelial responses although evidence has been conflicting. Here, we exposed airway epithelial cells from children with and without CF to RV in vitro. Using RNA-Seq, we profiled the transcriptomic differences of CF and non-CF airway epithelial cells at baseline and in response to RV. There were only modest differences between CF and non-CF cells at baseline. In response to RV, there were 1,442 and 896 differentially expressed genes in CF and non-CF airway epithelial cells, respectively. The core antiviral responses in CF and non-CF airway epithelial cells were mediated through interferon signaling although type 1 and 3 interferon signaling, when measured, were reduced in CF airway epithelial cells following viral challenge consistent with previous reports. The transcriptional responses in CF airway epithelial cells were more complex than in non-CF airway epithelial cells with diverse over-represented biological pathways, such as cytokine signaling and metabolic and biosynthetic pathways. Network analysis highlighted that the differentially expressed genes of CF airway epithelial cells' transcriptional responses were highly interconnected and formed a more complex network than observed in non-CF airway epithelial cells. We corroborate observations in fully differentiated air-liquid interface (ALI) cultures, identifying genes involved in IL-1 signaling and mucin glycosylation that are only dysregulated in the CF airway epithelial response to RV infection. These data provide novel insights into the CF airway epithelial cells' responses to RV infection and highlight potential pathways that could be targeted to improve antiviral and anti-inflammatory responses in CF.

Bilateral murine tumor models for characterizing the response to immune checkpoint blockade.

The therapeutic response to immune checkpoint blockade (ICB) is highly variable, not only between different cancers but also between patients with the same cancer type. The biological mechanisms underlying these differences in response are incompletely understood. Identifying correlates in patient tumor samples is challenging because of genetic and environmental variability. Murine studies usually compare different tumor models or treatments, introducing potential confounding variables. This protocol describes bilateral murine tumor models, derived from syngeneic cancer cell lines, that display a symmetrical yet dichotomous response to ICB. These models enable detailed analysis of whole tumors in a highly homogeneous background, combined with knowledge of the therapeutic outcome within a few weeks, and could potentially be used for mechanistic studies using other (immuno-)therapies. We discuss key considerations and describe how to use two cell lines as fully optimized models. We discuss experimental details, including proper inoculation technique to achieve symmetry and one-sided surgical tumor removal, which takes only 5 min per mouse. Furthermore, we outline the preparation of bulk tissue or single-cell suspensions for downstream analyses such as bulk RNA-seq, immunohistochemistry, single-cell RNA-seq and flow cytometry.

Expression Levels of Therapeutic Targets as Indicators of Sensitivity to Targeted Therapeutics.

Cancer precision medicine aims to predict the drug likely to yield the best response for a patient. Genomic sequencing of tumors is currently being used to better inform treatment options; however, this approach has had a limited clinical impact due to the paucity of actionable mutations. An alternative to mutation status is the use of gene expression signatures to predict response. Using data from two large-scale studies, The Genomics of Drug Sensitivity of Cancer (GDSC) and The Cancer Therapeutics Response Portal (CTRP), we investigated the relationship between the sensitivity of hundreds of cell lines to hundreds of drugs, and the relative expression levels of the targets these drugs are directed against. For approximately one third of the drugs considered (73/222 in GDSC and 131/360 in CTRP), sensitivity was significantly correlated with the expression of at least one of the known targets. Surprisingly, for 8% of the annotated targets, there was a significant anticorrelation between target expression and sensitivity. For several cases, this corresponded to drugs targeting multiple genes in the same family, with the expression of one target significantly correlated with sensitivity and another significantly anticorrelated suggesting a possible role in resistance. Furthermore, we identified nontarget genes that are significantly correlated or anticorrelated with drug sensitivity, and find literature linking several to sensitization and resistance. Our analyses provide novel and important insights into both potential mechanisms of resistance and relative efficacy of drugs against the same target.

Sensitization to immune checkpoint blockade through activation of a STAT1/NK axis in the tumor microenvironment.

Cancer immunotherapy using antibodies that target immune checkpoints has delivered outstanding results. However, responses only occur in a subset of patients, and it is not fully understood what biological processes determine an effective outcome. This lack of understanding hinders the development of rational combination treatments. We set out to define the pretreatment microenvironment associated with an effective outcome by using the fact that inbred mouse strains bearing monoclonal cancer cell line-derived tumors respond in a dichotomous manner to immune checkpoint blockade (ICB). We compared the cellular composition and gene expression profiles of responsive and nonresponsive tumors from mice before ICB and validated the findings in cohorts of patients with cancer treated with ICB antibodies. We found that responsive tumors were characterized by an inflammatory gene expression signature consistent with up-regulation of signal transducer and activator of transcription 1 (STAT1) and Toll-like receptor 3 (TLR3) signaling and down-regulation of interleukin-10 (IL-10) signaling. In addition, responsive tumors had more infiltrating-activated natural killer (NK) cells, which were necessary for response. Pretreatment of mice with large established tumors using the STAT1-activating cytokine interferon-γ (IFNγ), the TLR3 ligand poly(I:C), and an anti-IL-10 antibody sensitized tumors to ICB by attracting IFNγ-producing NK cells into the tumor, resulting in increased cure rates. Our results identify a pretreatment tumor microenvironment that predicts response to ICB, which can be therapeutically attained. These data suggest a biomarker-driven approach to patient management to establish whether a patient would benefit from treatment with sensitizing therapeutics before ICB.

C1 CAGE detects transcription start sites and enhancer activity at single-cell resolution.

Single-cell transcriptomic profiling is a powerful tool to explore cellular heterogeneity. However, most of these methods focus on the 3'-end of polyadenylated transcripts and provide only a partial view of the transcriptome. We introduce C1 CAGE, a method for the detection of transcript 5'-ends with an original sample multiplexing strategy in the C1TM microfluidic system. We first quantifiy the performance of C1 CAGE and find it as accurate and sensitive as other methods in the C1 system. We then use it to profile promoter and enhancer activities in the cellular response to TGF-ß of lung cancer cells and discover subpopulations of cells differing in their response. We also describe enhancer RNA dynamics revealing transcriptional bursts in subsets of cells with transcripts arising from either strand in a mutually exclusive manner, validated using single molecule fluorescence in situ hybridization.

Antisense Transcription in Loci Associated to Hereditary Neurodegenerative Diseases.

Natural antisense transcripts are common features of mammalian genes providing additional regulatory layers of gene expression. A comprehensive description of antisense transcription in loci associated to familial neurodegenerative diseases may identify key players in gene regulation and provide tools for manipulating gene expression. We take advantage of the FANTOM5 sequencing datasets that represent the largest collection to date of genome-wide promoter usage in almost 2000 human samples. Transcription start sites (TSSs) are mapped at high resolution by the use of a modified protocol of cap analysis of gene expression (CAGE) for high-throughput single molecule next-generation sequencing with Helicos (hCAGE). Here we present the analysis of antisense transcription at 17 loci associated to hereditary Alzheimer's disease, Frontotemporal Dementia, Parkinson's disease, Amyotrophic Lateral Sclerosis, and Huntington's disease. We focused our analysis on libraries derived from brain tissues and primary cells. We also screened libraries from total blood and blood cell populations in the quest for peripheral biomarkers of neurodegenerative diseases. We identified 63 robust promoters in antisense orientation to genes associated to familial neurodegeneration. When applying a less stringent cutoff, this number increases to over 400. A subset of these promoters represents alternative TSSs for 24 FANTOM5 annotated long noncoding RNA (lncRNA) genes, in antisense orientation to 13 of the loci analyzed here, while the remaining contribute to the expression of additional transcript variants. Intersection with GWAS studies, sample ontology, and dynamic expression reveals association to specific genetic traits as well as cell and tissue types, not limited to neurodegenerative diseases. Antisense transcription was validated for a subset of genes, including those encoding for Microtubule-Associated Protein Tau, α-synuclein, Parkinsonism-associated deglycase DJ-1, and Leucin-Rich Repeat Kinase 2. This work provides evidence for the existence of additional regulatory mechanisms of the expression of neurodegenerative disease-causing genes by previously not-annotated and/or not-validated antisense long noncoding RNAs.

Elucidating the Interaction of CF Airway Epithelial Cells and Rhinovirus: Using the Host-Pathogen Relationship to Identify Future Therapeutic Strategies.

Chronic lung disease remains the primary cause of mortality in cystic fibrosis (CF). Growing evidence suggests respiratory viral infections are often more severe in CF compared to healthy peers and contributes to pulmonary exacerbations (PEx) and deterioration of lung function. Rhinovirus is the most prevalent respiratory virus detected, particularly during exacerbations in children with CF <5 years old. However, even though rhinoviral infections are likely to be one of the factors initiating the onset of CF lung disease, there is no effective targeted treatment. A better understanding of the innate immune responses by CF airway epithelial cells, the primary site of infection for viruses, is needed to identify why viral infections are more severe in CF. The aim of this review is to present the clinical impact of virus infection in both young children and adults with CF, focusing on rhinovirus infection. Previous in vitro and in vivo investigations looking at the mechanisms behind virus infection will also be summarized. The review will finish on the potential of transcriptomics to elucidate the host-pathogen responses by CF airway cells to viral infection and identify novel therapeutic targets.

Conserved temporal ordering of promoter activation implicates common mechanisms governing the immediate early response across cell types and stimuli.

The promoters of immediate early genes (IEGs) are rapidly activated in response to an external stimulus. These genes, also known as primary response genes, have been identified in a range of cell types, under diverse extracellular signals and using varying experimental protocols. Whereas genomic dissection on a case-by-case basis has not resulted in a comprehensive catalogue of IEGs, a rigorous meta-analysis of eight genome-wide FANTOM5 CAGE (cap analysis of gene expression) time course datasets reveals successive waves of promoter activation in IEGs, recapitulating known relationships between cell types and stimuli: we obtain a set of 57 (42 protein-coding) candidate IEGs possessing promoters that consistently drive a rapid but transient increase in expression over time. These genes show significant enrichment for known IEGs reported previously, pathways associated with the immediate early response, and include a number of non-coding RNAs with roles in proliferation and differentiation. Surprisingly, we also find strong conservation of the ordering of activation for these genes, such that 77 pairwise promoter activation orderings are conserved. Using the leverage of comprehensive CAGE time series data across cell types, we also document the extensive alternative promoter usage by such genes, which is likely to have been a barrier to their discovery until now. The common activation ordering of the core set of early-responding genes we identify may indicate conserved underlying regulatory mechanisms. By contrast, the considerably larger number of transiently activated genes that are specific to each cell type and stimulus illustrates the breadth of the primary response.

Transcriptional landscape of Mycobacterium tuberculosis infection in macrophages.

Mycobacterium tuberculosis (Mtb) infection reveals complex and dynamic host-pathogen interactions, leading to host protection or pathogenesis. Using a unique transcriptome technology (CAGE), we investigated the promoter-based transcriptional landscape of IFNγ (M1) or IL-4/IL-13 (M2) stimulated macrophages during Mtb infection in a time-kinetic manner. Mtb infection widely and drastically altered macrophage-specific gene expression, which is far larger than that of M1 or M2 activations. Gene Ontology enrichment analysis for Mtb-induced differentially expressed genes revealed various terms, related to host-protection and inflammation, enriched in up-regulated genes. On the other hand, terms related to dis-regulation of cellular functions were enriched in down-regulated genes. Differential expression analysis revealed known as well as novel transcription factor genes in Mtb infection, many of them significantly down-regulated. IFNγ or IL-4/IL-13 pre-stimulation induce additional differentially expressed genes in Mtb-infected macrophages. Cluster analysis uncovered significant numbers, prolonging their expressional changes. Furthermore, Mtb infection augmented cytokine-mediated M1 and M2 pre-activations. In addition, we identified unique transcriptional features of Mtb-mediated differentially expressed lncRNAs. In summary we provide a comprehensive in depth gene expression/regulation profile in Mtb-infected macrophages, an important step forward for a better understanding of host-pathogen interaction dynamics in Mtb infection.

Shared activity patterns arising at genetic susceptibility loci reveal underlying genomic and cellular architecture of human disease.

Genetic variants underlying complex traits, including disease susceptibility, are enriched within the transcriptional regulatory elements, promoters and enhancers. There is emerging evidence that regulatory elements associated with particular traits or diseases share similar patterns of transcriptional activity. Accordingly, shared transcriptional activity (coexpression) may help prioritise loci associated with a given trait, and help to identify underlying biological processes. Using cap analysis of gene expression (CAGE) profiles of promoter- and enhancer-derived RNAs across 1824 human samples, we have analysed coexpression of RNAs originating from trait-associated regulatory regions using a novel quantitative method (network density analysis; NDA). For most traits studied, phenotype-associated variants in regulatory regions were linked to tightly-coexpressed networks that are likely to share important functional characteristics. Coexpression provides a new signal, independent of phenotype association, to enable fine mapping of causative variants. The NDA coexpression approach identifies new genetic variants associated with specific traits, including an association between the regulation of the OCT1 cation transporter and genetic variants underlying circulating cholesterol levels. NDA strongly implicates particular cell types and tissues in disease pathogenesis. For example, distinct groupings of disease-associated regulatory regions implicate two distinct biological processes in the pathogenesis of ulcerative colitis; a further two separate processes are implicated in Crohn's disease. Thus, our functional analysis of genetic predisposition to disease defines new distinct disease endotypes. We predict that patients with a preponderance of susceptibility variants in each group are likely to respond differently to pharmacological therapy. Together, these findings enable a deeper biological understanding of the causal basis of complex traits.

Transcriptome Analysis Uncovers a Growth-Promoting Activity of Orosomucoid-1 on Hepatocytes.

The acute phase protein orosomucoid-1 (Orm1) is mainly expressed by hepatocytes (HPCs) under stress conditions. However, its specific function is not fully understood. Here, we report a role of Orm1 as an executer of HPC proliferation. Increases in serum levels of Orm1 were observed in patients after surgical resection for liver cancer and in mice undergone partial hepatectomy (PH). Transcriptome study showed that Orm1 became the most abundant in HPCs isolated from regenerating mouse liver tissues after PH. Both in vitro and in vivo siRNA-induced knockdown of Orm1 suppressed proliferation of mouse regenerating HPCs and human hepatic cells. Microarray analysis in regenerating mouse livers revealed that the signaling pathways controlling chromatin replication, especially the minichromosome maintenance protein complex genes were uniformly down-regulated following Orm1 knockdown. These data suggest that Orm1 is induced in response to hepatic injury and executes liver regeneration by activating cell cycle progression in HPCs.

Integrative CAGE and DNA Methylation Profiling Identify Epigenetically Regulated Genes in NSCLC.

Lung cancer is the leading cause of cancer-related deaths worldwide. The majority of cancer driver mutations have been identified; however, relevant epigenetic regulation involved in tumorigenesis has only been fragmentarily analyzed. Epigenetically regulated genes have a great theranostic potential, especially in tumors with no apparent driver mutations. Here, epigenetically regulated genes were identified in lung cancer by an integrative analysis of promoter-level expression profiles from Cap Analysis of Gene Expression (CAGE) of 16 non-small cell lung cancer (NSCLC) cell lines and 16 normal lung primary cell specimens with DNA methylation data of 69 NSCLC cell lines and 6 normal lung epithelial cells. A core set of 49 coding genes and 10 long noncoding RNAs (lncRNA), which are upregulated in NSCLC cell lines due to promoter hypomethylation, was uncovered. Twenty-two epigenetically regulated genes were validated (upregulated genes with hypomethylated promoters) in the adenocarcinoma and squamous cell cancer subtypes of lung cancer using The Cancer Genome Atlas data. Furthermore, it was demonstrated that multiple copies of the REP522 DNA repeat family are prominently upregulated due to hypomethylation in NSCLC cell lines, which leads to cancer-specific expression of lncRNAs, such as RP1-90G24.10, AL022344.4, and PCAT7. Finally, Myeloma Overexpressed (MYEOV) was identified as the most promising candidate. Functional studies demonstrated that MYEOV promotes cell proliferation, survival, and invasion. Moreover, high MYEOV expression levels were associated with poor prognosis.Implications: This report identifies a robust list of 22 candidate driver genes that are epigenetically regulated in lung cancer; such genes may complement the known mutational drivers.Visual Overview: http://mcr.aacrjournals.org/content/molcanres/15/10/1354/F1.large.jpg Mol Cancer Res; 15(10); 1354-65. ©2017 AACR.

Analysis of the human monocyte-derived macrophage transcriptome and response to lipopolysaccharide provides new insights into genetic aetiology of inflammatory bowel disease.

The FANTOM5 consortium utilised cap analysis of gene expression (CAGE) to provide an unprecedented insight into transcriptional regulation in human cells and tissues. In the current study, we have used CAGE-based transcriptional profiling on an extended dense time course of the response of human monocyte-derived macrophages grown in macrophage colony-stimulating factor (CSF1) to bacterial lipopolysaccharide (LPS). We propose that this system provides a model for the differentiation and adaptation of monocytes entering the intestinal lamina propria. The response to LPS is shown to be a cascade of successive waves of transient gene expression extending over at least 48 hours, with hundreds of positive and negative regulatory loops. Promoter analysis using motif activity response analysis (MARA) identified some of the transcription factors likely to be responsible for the temporal profile of transcriptional activation. Each LPS-inducible locus was associated with multiple inducible enhancers, and in each case, transient eRNA transcription at multiple sites detected by CAGE preceded the appearance of promoter-associated transcripts. LPS-inducible long non-coding RNAs were commonly associated with clusters of inducible enhancers. We used these data to re-examine the hundreds of loci associated with susceptibility to inflammatory bowel disease (IBD) in genome-wide association studies. Loci associated with IBD were strongly and specifically (relative to rheumatoid arthritis and unrelated traits) enriched for promoters that were regulated in monocyte differentiation or activation. Amongst previously-identified IBD susceptibility loci, the vast majority contained at least one promoter that was regulated in CSF1-dependent monocyte-macrophage transitions and/or in response to LPS. On this basis, we concluded that IBD loci are strongly-enriched for monocyte-specific genes, and identified at least 134 additional candidate genes associated with IBD susceptibility from reanalysis of published GWA studies. We propose that dysregulation of monocyte adaptation to the environment of the gastrointestinal mucosa is the key process leading to inflammatory bowel disease.

Transcriptome analysis of periodontitis-associated fibroblasts by CAGE sequencing identified DLX5 and RUNX2 long variant as novel regulators involved in periodontitis.

Periodontitis is affecting over half of the adult population, and represents a major public health problem. Previously, we isolated a subset of gingival fibroblasts (GFs) from periodontitis patients, designated as periodontitis-associated fibroblasts (PAFs), which were highly capable of collagen degradation. To elucidate their molecular profiles, GFs isolated form healthy and periodontitis-affected gingival tissues were analyzed by CAGE-seq and integrated with the FANTOM5 atlas. GFs from healthy gingival tissues displayed distinctive patterns of CAGE profiles as compared to fibroblasts from other organ sites and characterized by specific expression of developmentally important transcription factors such as BARX1, PAX9, LHX8, and DLX5. In addition, a novel long non-coding RNA associated with LHX8 was described. Furthermore, we identified DLX5 regulating expression of the long variant of RUNX2 transcript, which was specifically active in GFs but not in their periodontitis-affected counterparts. Knockdown of these factors in GFs resulted in altered expression of extracellular matrix (ECM) components. These results indicate activation of DLX5 and RUNX2 via its distal promoter represents a unique feature of GFs, and is important for ECM regulation. Down-regulation of these transcription factors in PAFs could be associated with their property to degrade collagen, which may impact on the process of periodontitis.