Anupam Sinha

Combining transcription factor binding affinities with open-chromatin data for accurate gene expression prediction

The binding and contribution of transcription factors (TF) to cell specific gene expression is often deduced from open-chromatin measurements to avoid costly TF ChIP-seq assays. Thus, it is important to develop computational methods for accurate TF binding prediction in open-chromatin regions (OCRs). Here, we report a novel segmentation-based method, TEPIC, to predict TF binding by combining sets of OCRs with position weight matrices. TEPIC can be applied to various open-chromatin data, e.g. DNaseI-seq and NOMe-seq. Additionally, Histone-Marks (HMs) can be used to identify candidate TF binding sites. TEPIC computes TF affinities and uses open-chromatin/HM signal intensity as quantitative measures of TF binding strength. Using machine learning, we find low affinity binding sites to improve our ability to explain gene expression variability compared to the standard presence/absence classification of binding sites. Further, we show that both footprints and peaks capture essential TF binding events and lead to a good prediction performance. In our application, gene-based scores computed by TEPIC with one open-chromatin assay nearly reach the quality of several TF ChIP-seq data sets. Finally, these scores correctly predict known transcriptional regulators as illustrated by the application to novel DNaseI-seq and NOMe-seq data for primary human hepatocytes and CD4+ T-cells, respectively.

Uncoupling of mucosal gene regulation, mRNA splicing and adherent microbiota signatures in inflammatory bowel disease

Objective An inadequate host response to the intestinal microbiota likely contributes to the manifestation and progression of human inflammatory bowel disease (IBD). However, molecular approaches to unravelling the nature of the defective crosstalk and its consequences for intestinal metabolic and immunological networks are lacking. We assessed the mucosal transcript levels, splicing architecture and mucosa-attached microbial communities of patients with IBD to obtain a comprehensive view of the underlying, hitherto poorly characterised interactions, and how these are altered in IBD. Design Mucosal biopsies from Crohns disease and patients with UC, disease controls and healthy individuals (n=63) were subjected to microbiome, transcriptome and splicing analysis, employing next-generation sequencing. The three data levels were integrated by different bioinformatic approaches, including systems biology-inspired network and pathway analysis. Results Microbiota, host transcript levels and host splicing patterns were influenced most strongly by tissue differences, followed by the effect of inflammation. Both factors point towards a substantial disease-related alteration of metabolic processes. We also observed a strong enrichment of splicing events in inflamed tissues, accompanied by an alteration of the mucosa-attached bacterial taxa. Finally, we noted a striking uncoupling of the three molecular entities when moving from healthy individuals via disease controls to patients with IBD. Conclusions Our results provide strong evidence that the interplay between microbiome and host transcriptome, which normally characterises a state of intestinal homeostasis, is drastically perturbed in Crohns disease and UC. Consequently, integrating multiple OMICs levels appears to be a promising approach to further disentangle the complexity of IBD.

DNA Methylation and Transcription Patterns in Intestinal Epithelial Cells From Pediatric Patients With Inflammatory Bowel Diseases Differentiate Disease Subtypes and Associate With Outcome.

Background & Aims We analyzed DNA methylation patterns and transcriptomes of primary intestinal epithelial cells (IEC) of children newly diagnosed with inflammatory bowel diseases (IBD) to learn more about pathogenesis. Methods We obtained mucosal biopsies (N = 236) collected from terminal ileum and ascending and sigmoid colons of children (median age 13 years) newly diagnosed with IBD (43 with Crohn’s disease [CD], 23 with ulcerative colitis [UC]), and 30 children without IBD (controls). Patients were recruited and managed at a hospital in the United Kingdom from 2013 through 2016. We also obtained biopsies collected at later stages from a subset of patients. IECs were purified and analyzed for genome-wide DNA methylation patterns and gene expression profiles. Adjacent microbiota were isolated from biopsies and analyzed by 16S gene sequencing. We generated intestinal organoid cultures from a subset of samples and genome-wide DNA methylation analysis was performed. Results We found gut segment-specific differences in DNA methylation and transcription profiles of IECs from children with IBD vs controls; some were independent of mucosal inflammation. Changes in gut microbiota between IBD and control groups were not as large and were difficult to assess because of large amounts of intra-individual variation. Only IECs from patients with CD had changes in DNA methylation and transcription patterns in terminal ileum epithelium, compared with controls. Colon epithelium from patients with CD and from patients with ulcerative colitis had distinct changes in DNA methylation and transcription patterns, compared with controls. In IECs from patients with IBD, changes in DNA methylation, compared with controls, were stable over time and were partially retained in ex-vivo organoid cultures. Statistical analyses of epithelial cell profiles allowed us to distinguish children with CD or UC from controls; profiles correlated with disease outcome parameters, such as the requirement for treatment with biologic agents. Conclusions We identified specific changes in DNA methylation and transcriptome patterns in IECs from pediatric patients with IBD compared with controls. These data indicate that IECs undergo changes during IBD development and could be involved in pathogenesis. Further analyses of primary IECs from patients with IBD could improve our understanding of the large variations in disease progression and outcomes.

Genetic interplay between human longevity and metabolic pathways: a large-scale eQTL study

Human longevity is a complex phenotype influenced by genetic and environmental components. Unraveling the contribution of genetic vs. nongenetic factors to longevity is a challenging task. Here, we conducted a large‐scale RNA‐sequencing‐based expression quantitative trait loci study (eQTL) with subsequent heritability analysis. The investigation was performed on blood samples from 244 individuals from Germany and Denmark, representing various age groups including long‐lived subjects up to the age of 104 years. Our eQTL‐based approach revealed for the first time that human longevity is associated with a depletion of metabolic pathways in a genotype‐dependent and independent manner. Further analyses indicated that 20% of the differentially expressed genes are influenced by genetic variants in cis. The subsequent study of twins showed that the transcriptional activity of a third of the differentially regulated genes is heritable. These findings suggest that longevity‐associated biological processes such as altered metabolism are, to a certain extent, also the driving force of longevity rather than just a consequence of old age.

ADAM17 is required for EGF-R-induced intestinal tumors via IL-6 trans-signaling.

Colorectal cancer is treated with antibodies blocking epidermal growth factor receptor (EGF-R), but therapeutic success is limited. EGF-R is stimulated by soluble ligands, which are derived from transmembrane precursors by ADAM17-mediated proteolytic cleavage. In mouse intestinal cancer models in the absence of ADAM17, tumorigenesis was almost completely inhibited, and the few remaining tumors were of low-grade dysplasia. RNA sequencing analysis demonstrated down-regulation of STAT3 and Wnt pathway components. Because EGF-R on myeloid cells, but not on intestinal epithelial cells, is required for intestinal cancer and because IL-6 is induced via EGF-R stimulation, we analyzed the role of IL-6 signaling. Tumor formation was equally impaired in IL-6−/− mice and sgp130Fc transgenic mice, in which only trans-signaling via soluble IL-6R is abrogated. ADAM17 is needed for EGF-R–mediated induction of IL-6 synthesis, which via IL-6 trans-signaling induces &bgr;-catenin–dependent tumorigenesis. Our data reveal the possibility of a novel strategy for treatment of colorectal cancer that could circumvent intrinsic and acquired resistance to EGF-R blockade.

Vedolizumab is associated with changes in innate rather than adaptive immunity in patients with inflammatory bowel disease

Objective Vedolizumab, a monoclonal antibody directed against the integrin heterodimer α4β7, is approved for the treatment of Crohn’s disease and ulcerative colitis. The efficacy of vedolizumab has been suggested to result from inhibition of intestinal T cell trafficking although human data to support this conclusion are scarce. We therefore performed a comprehensive analysis of vedolizumab-induced alterations in mucosal and systemic immunity in patients with inflammatory bowel disease (IBD), using anti-inflammatory therapy with the TNFα antibody infliximab as control. Design Immunophenotyping, immunohistochemistry, T cell receptor profiling and RNA sequencing were performed using blood and colonic biopsies from patients with IBD before and during treatment with vedolizumab (n=18) or, as control, the anti-TNFα antibody infliximab (n=20). Leucocyte trafficking in vivo was assessed using single photon emission computed tomography and endomicroscopy. Results Vedolizumab was not associated with alterations in the abundance or phenotype of lamina propria T cells and did not affect the mucosal T cell repertoire or leucocyte trafficking in vivo. Surprisingly, however, α4β7 antibody treatment was associated with substantial effects on innate immunity including changes in macrophage populations and pronounced alterations in the expression of molecules involved in microbial sensing, chemoattraction and regulation of the innate effector response. These effects were specific to vedolizumab, not observed in response to the TNFα antibody infliximab, and associated with inhibition of intestinal inflammation. Conclusion Our findings suggest that modulation of innate immunity contributes to the therapeutic efficacy of vedolizumab in IBD. Trial registration number NCT02694588

Integrative analysis of single cell expression data reveals distinct regulatory states in bidirectional promoters

Background Bidirectional promoters (BPs) are prevalent in eukaryotic genomes. However, it is poorly understood how the cell integrates different epigenomic information, such as transcription factor (TF) binding and chromatin marks, to drive gene expression at BPs. Single cell sequencing technologies are revolutionizing the field of genome biology. Therefore, this study focuses on the integration of single cell RNA-seq data with bulk ChIP-seq and other epigenetics data, for which single cell technologies are not yet established, in the context of BPs. Results We performed integrative analyses of novel human single cell RNA-seq (scRNA-seq) data with bulk ChIP-seq and other epigenetics data. scRNA-seq data revealed distinct transcription states of BPs that were previously not recognized. We find associations between these transcription states to distinct patterns in structural gene features, DNA accessibility, histone modification, DNA methylation and TF binding profiles. Conclusions Our results suggest that a complex interplay of all of these elements is required to achieve BP-specific transcriptional output in this specialized promoter configuration. Further, our study implies that novel statistical methods can be developed to deconvolute masked subpopulations of cells measured with different bulk epigenomic assays using scRNA-seq data.

Epithelial RNase H2 Maintains Genome Integrity and Prevents Intestinal Tumorigenesis in Mice

Background & Aims RNase H2 is a holoenzyme, composed of 3 subunits (ribonuclease H2 subunits A, B, and C), that cleaves RNA:DNA hybrids and removes mis-incorporated ribonucleotides from genomic DNA through ribonucleotide excision repair. Ribonucleotide incorporation by eukaryotic DNA polymerases occurs during every round of genome duplication and produces the most frequent type of naturally occurring DNA lesion. We investigated whether intestinal epithelial proliferation requires RNase H2 function and whether RNase H2 activity is disrupted during intestinal carcinogenesis. Methods We generated mice with epithelial-specific deletion of ribonuclease H2 subunit B (H2bΔIEC) and mice that also had deletion of tumor-suppressor protein p53 (H2b/p53ΔIEC); we compared phenotypes with those of littermate H2bfl/fl or H2b/p53fl/fl (control) mice at young and old ages. Intestinal tissues were collected and analyzed by histology. We isolated epithelial cells, generated intestinal organoids, and performed RNA sequence analyses. Mutation signatures of spontaneous tumors from H2b/p53ΔIEC mice were characterized by exome sequencing. We collected colorectal tumor specimens from 467 patients, measured levels of ribonuclease H2 subunit B, and associated these with patient survival times and transcriptome data. Results The H2bΔIEC mice had DNA damage to intestinal epithelial cells and proliferative exhaustion of the intestinal stem cell compartment compared with controls and H2b/p53ΔIEC mice. However, H2b/p53ΔIEC mice spontaneously developed small intestine and colon carcinomas. DNA from these tumors contained T>G base substitutions at GTG trinucleotides. Analyses of transcriptomes of human colorectal tumors associated lower levels of RNase H2 with shorter survival times. Conclusions In analyses of mice with disruption of the ribonuclease H2 subunit B gene and colorectal tumors from patients, we provide evidence that RNase H2 functions as a colorectal tumor suppressor. H2b/p53ΔIEC mice can be used to study the roles of RNase H2 in tissue-specific carcinogenesis.

A Drosophila model of cigarette smoke induced COPD identifies Nrf2 signaling as an expedient target for intervention

Chronic obstructive pulmonary disease (COPD) is among the most important causes of death. Signaling systems that are relevant for tissue repair and detoxification of reactive oxygen species or xenobiotics are thought to be impaired in lungs of patients suffering from this disease. Here, we developed a simple cigarette smoke induced Drosophila model of COPD based on chronic cigarette smoke exposure that recapitulates major pathological hallmarks of the disease and thus can be used to investigate new therapeutic strategies. Chronic cigarette smoke exposure led to premature death of the animals and induced a set of phenotypes reminiscent of those seen in COPD patients, including reduced physical activity, reduced body fat, increased metabolic rate and a substantial reduction of the respiratory surface. A detailed transcriptomic analysis revealed that especially the TGF-β, Nrf2 and the JAK/STAT signaling pathways are altered by chronic cigarette smoke exposure. Based on these results, we focused on Nrf2 signaling. A pharmacological intervention study performed with oltipraz, an activator of Nrf2 signaling, increased survival of cigarette smoke exposed animals significantly. Thus, the Drosophila COPD model recapitulates many major hallmarks of COPD and it is highly useful to evaluate the potential of alternative therapeutic strategies.

Regions of common inter-individual DNA methylation differences in human monocytes: genetic basis and potential function

BackgroundThere is increasing evidence for inter-individual methylation differences at CpG dinucleotides in the human genome, but the regional extent and function of these differences have not yet been studied in detail. For identifying regions of common methylation differences, we used whole genome bisulfite sequencing data of monocytes from five donors and a novel bioinformatic strategy.ResultsWe identified 157 differentially methylated regions (DMRs) with four or more CpGs, almost none of which has been described before. The DMRs fall into different chromatin states, where methylation is inversely correlated with active, but not repressive histone marks. However, methylation is not correlated with the expression of associated genes. High-resolution single nucleotide polymorphism (SNP) genotyping of the five donors revealed evidence for a role of cis-acting genetic variation in establishing methylation patterns. To validate this finding in a larger cohort, we performed genome-wide association studies (GWAS) using SNP genotypes and 450k array methylation data from blood samples of 1128 individuals. Only 30/157 (19%) DMRs include at least one 450k CpG, which shows that these arrays miss a large proportion of DNA methylation variation. In most cases, the GWAS peak overlapped the CpG position, and these regions are enriched for CREB group, NF-1, Sp100 and CTCF binding motifs. In two cases, there was tentative evidence for a trans-effect by KRAB zinc finger proteins.ConclusionsAllele-specific DNA methylation occurs in discrete chromosomal regions and is driven by genetic variation in cis and trans, but in general has little effect on gene expression.