Ayla Ergun

MicroRNAs as modulators of smoking-induced gene expression changes in human airway epithelium

We have shown that smoking impacts bronchial airway gene expression and that heterogeneity in this response associates with smoking-related disease risk. In this study, we sought to determine whether microRNAs (miRNAs) play a role in regulating the airway gene expression response to smoking. We examined whole-genome miRNA and mRNA expression in bronchial airway epithelium from current and never smokers (n = 20) and found 28 miRNAs to be differentially expressed (P < 0.05) with the majority being down-regulated in smokers. We further identified a number of mRNAs whose expression level is highly inversely correlated with miRNA expression in vivo. Many of these mRNAs contain potential binding sites for the differentially expressed miRNAs in their 3′-untranslated region (UTR) and are themselves affected by smoking. We found that either increasing or decreasing the levels of mir-218 (a miRNA that is strongly affected by smoking) in both primary bronchial epithelial cells and H1299 cells was sufficient to cause a corresponding decrease or increase in the expression of predicted mir-218 mRNA targets, respectively. Further, mir-218 expression is reduced in primary bronchial epithelium exposed to cigarette smoke condensate (CSC), and alteration of mir-218 levels in these cells diminishes the induction of the predicted mir-218 target MAFG in response to CSC. These data indicate that mir-218 levels modulate the airway epithelial gene expression response to cigarette smoke and support a role for miRNAs in regulating host response to environmental toxins.

A network biology approach to prostate cancer

There is a need to identify genetic mediators of solid‐tumor cancers, such as prostate cancer, where invasion and distant metastases determine the clinical outcome of the disease. Whole‐genome expression profiling offers promise in this regard, but can be complicated by the challenge of identifying the genes affected by a condition from the hundreds to thousands of genes that exhibit changes in expression. Here, we show that reverse‐engineered gene networks can be combined with expression profiles to compute the likelihood that genes and associated pathways are mediators of a disease. We apply our method to non‐recurrent primary and metastatic prostate cancer data, and identify the androgen receptor gene (AR) among the top genetic mediators and the AR pathway as a highly enriched pathway for metastatic prostate cancer. These results were not obtained on the basis of expression change alone. We further demonstrate that the AR gene, in the context of the network, can be used as a marker to detect the aggressiveness of primary prostate cancers. This work shows that a network biology approach can be used advantageously to identify the genetic mediators and mediating pathways associated with a disease.

microRNA Regulatory Network Inference Identifies miR-34a as a Novel Regulator of TGF-β Signaling in Glioblastoma

UNLABELLED Leveraging The Cancer Genome Atlas (TCGA) multidimensional data in glioblastoma, we inferred the putative regulatory network between microRNA and mRNA using the Context Likelihood of Relatedness modeling algorithm. Interrogation of the network in context of defined molecular subtypes identified 8 microRNAs with a strong discriminatory potential between proneural and mesenchymal subtypes. Integrative in silico analyses, a functional genetic screen, and experimental validation identified miR-34a as a tumor suppressor in proneural subtype glioblastoma. Mechanistically, in addition to its direct regulation of platelet-derived growth factor receptor-alpha (PDGFRA), promoter enrichment analysis of context likelihood of relatedness-inferred mRNA nodes established miR-34a as a novel regulator of a SMAD4 transcriptional network. Clinically, miR-34a expression level is shown to be prognostic, where miR-34a low-expressing glioblastomas exhibited better overall survival. This work illustrates the potential of comprehensive multidimensional cancer genomic data combined with computational and experimental models in enabling mechanistic exploration of relationships among different genetic elements across the genome space in cancer. SIGNIFICANCE We illustrate here that network modeling of complex multidimensional cancer genomic data can generate a framework in which to explore the biology of cancers, leading to discovery of new pathogenetic insights as well as potential prognostic biomarkers. Specifically in glioblastoma, within the context of the global network, promoter enrichment analysis of network edges uncovered a novel regulation of TGF-β signaling via a Smad4 transcriptomic network by miR-34a.

Interindividual variation in human T regulatory cells

Significance Control of immunologic tolerance and homeostasis rely on regulatory T lymphocytes that express the transcription factor FOXP3. To characterize the interindividual variation of human Treg cells, we performed a genome-wide expression and genotypic analysis of 168 human donors, healthy or affected by type-1 or type-2 diabetes (T1D, T2D). We identify cis-acting genetic variants that condition Treg effector but not specification genes, and gene clusters that suggest Treg-specific regulatory pathways for some key signature genes (CTLA4, DUSP4). We also identify factors that may control FOXP3 mRNA or protein expression, the specification of the Treg signature, and Treg suppressive efficacy. Although no single transcript correlates with diabetes, overall expression of the Treg signature is perturbed in T1D, but not T2D, patients. FOXP3+ regulatory T (Treg) cells enforce immune self-tolerance and homeostasis, and variation in some aspects of Treg function may contribute to human autoimmune diseases. Here, we analyzed population-level Treg variability by performing genome-wide expression profiling of CD4+ Treg and conventional CD4+ T (Tconv) cells from 168 donors, healthy or with established type-1 diabetes (T1D) or type-2 diabetes (T2D), in relation to genetic and immunologic screening. There was a range of variability in Treg signature transcripts, some almost invariant, others more variable, with more extensive variability for genes that control effector function (ENTPD1, FCRL1) than for lineage-specification factors like FOXP3 or IKZF2. Network analysis of Treg signature genes identified coregulated clusters that respond similarly to genetic and environmental variation in Treg and Tconv cells, denoting qualitative differences in otherwise shared regulatory circuits whereas other clusters are coregulated in Treg, but not Tconv, cells, suggesting Treg-specific regulation of genes like CTLA4 or DUSP4. Dense genotyping identified 110 local genetic variants (cis-expression quantitative trait loci), some of which are specifically active in Treg, but not Tconv, cells. The Treg signature became sharper with age and with increasing body-mass index, suggesting a tuning of Treg function with repertoire selection and/or chronic inflammation. Some Treg signature transcripts correlated with FOXP3 mRNA and/or protein, suggesting transcriptional or posttranslational regulatory relationships. Although no single transcript showed significant association to diabetes, overall expression of the Treg signature was subtly perturbed in T1D, but not T2D, patients.

Chromatin topology is coupled to Polycomb group protein subnuclear organization

The genomes of metazoa are organized at multiple scales. Many proteins that regulate genome architecture, including Polycomb group (PcG) proteins, form subnuclear structures. Deciphering mechanistic links between protein organization and chromatin architecture requires precise description and mechanistic perturbations of both. Using super-resolution microscopy, here we show that PcG proteins are organized into hundreds of nanoscale protein clusters. We manipulated PcG clusters by disrupting the polymerization activity of the sterile alpha motif (SAM) of the PcG protein Polyhomeotic (Ph) or by increasing Ph levels. Ph with mutant SAM disrupts clustering of endogenous PcG complexes and chromatin interactions while elevating Ph level increases cluster number and chromatin interactions. These effects can be captured by molecular simulations based on a previously described chromatin polymer model. Both perturbations also alter gene expression. Organization of PcG proteins into small, abundant clusters on chromatin through Ph SAM polymerization activity may shape genome architecture through chromatin interactions.

A cluster of coregulated genes determines TGF-β–induced regulatory T-cell (Treg) dysfunction in NOD mice

Foxp3+ regulatory T cells (Tregs) originate in the thymus, but the Treg phenotype can also be induced in peripheral lymphoid organs or in vitro by stimulation of conventional CD4+ T cells with IL-2 and TGF-β. There have been divergent reports on the suppressive capacity of these TGF-Treg cells. We find that TGF-Tregs derived from diabetes-prone NOD mice, although expressing normal Foxp3 levels, are uniquely defective in suppressive activity, whereas TGF-Tregs from control strains (B6g7) or ex vivo Tregs from NOD mice all function normally. Most Treg-typical transcripts were shared by NOD or B6g7 TGF-Tregs, except for a small group of differentially expressed genes, including genes relevant for suppressive activity (Lrrc32, Ctla4, and Cd73). Many of these transcripts form a coregulated cluster in a broader analysis of T-cell differentiation. The defect does not map to idd3 or idd5 regions. Whereas Treg cells from NOD mice are normal in spleen and lymph nodes, the NOD defect is observed in locations that have been tied to pathogenesis of diabetes (small intestine lamina propria and pancreatic lymph node). Thus, a genetic defect uniquely affects a specific Treg subpopulation in NOD mice, in a manner consistent with a role in determining diabetes susceptibility.

Differential splicing across immune system lineages

Alternative splicing (AS) allows increased diversity and orthogonal regulation of the transcriptional products of mammalian genomes. To assess the distribution and variation of alternative splicing across cell lineages of the immune system, we comprehensively analyzed RNA sequencing and microarray data generated by the Immunological Genome Project Consortium. AS is pervasive: 60% of genes showed frequent AS isoforms in T or B lymphocytes, with 7,599 previously unreported isoforms. Distinct cell specificity was observed, with differential exon skipping in 5% of genes otherwise coexpressed in both B and T cells. The distribution of isoforms was mostly all or none, suggesting on/off switching as a frequent mode of AS regulation in lymphocytes. From the identification of differential exon use in the microarray data, clustering of exon inclusion/exclusion patterns across all Immunological Genome Project cell types showed that ∼70% of AS exons are distributed along a common pattern linked to lineage differentiation and cell cycling. Other AS events distinguished myeloid from lymphoid cells or affected only a small set of exons without clear lineage specificity (e.g., Ptprc). Computational analysis predicted specific associations between AS exons and splicing regulators, which were verified by detection of the hnRPLL/Ptprc connection.

Secreted Metabolites of Bifidobacterium infantis and Lactobacillus acidophilus Protect Immature Human Enterocytes from IL-1β-Induced Inflammation: A Transcription Profiling Analysis

Combination regimens of Bifidobacterium infantis and Lactobacillus acidophilus have been demonstrated to prevent necrotizing enterocolitis (NEC) in clinical trials. However, the molecular mechanisms responsible for this protective effect are not well understood. Additionally, conditioned media from individual cultures of these two probiotics show strain specific modulation of inflammation using in vitro human intestinal NEC models. Here we report a transcription profiling analysis of gene expression in immature human fetal intestinal epithelial cells (H4 cells) pretreated with conditioned media from B. infantis (BCM) or L. acidophilus (LCM) prior to IL-1β stimulation. Compared with control media, the two probiotic-conditioned media (PCM) treatments altered the expression of hundreds of genes involved in the immune response, apoptosis and cell survival, cell adhesion, the cell cycle, development and angiogenesis. In IL-1β-stimulated cells, PCM treatment decreased the upregulation of genes in the NF-κB activation pathway and downregulated genes associated with extracellular matrix (ECM) remodeling. Compared with LCM, BCM showed more significant modulatory effects on ECM remodeling, reflected by a lower p value. IL-6 and IL-8 production was significantly reduced in IL-1β-stimulated cells pretreated with PCM (p<0.05), which was consistent with their altered gene expression. Western blot analysis showed that compared with IL-1β stimulation alone, PCM treatment attenuated the decrease of cytoplasmic IκBα and NF-κB p65 levels as well as the increase of nuclear NF-κB p65 levels in the stimulated cells (p<0.05). In conclusion, PCM treatment exerted anti-inflammatory effects in immature human fetal enterocytes primarily by modulating genes in the NF-κB signaling and ECM remodeling pathways. Additionally, some components of these signaling pathways, particularly the ECM remodeling pathway, were more profoundly affected by BCM than LCM.

CAT7 and cat7l long non-coding RNAs Tune Polycomb Repressive Complex 1 Function During Human and Zebrafish Development.

The essential functions of polycomb repressive complex 1 (PRC1) in development and gene silencing are thought to involve long non-coding RNAs (lncRNAs), but few specific lncRNAs that guide PRC1 activity are known. We screened for lncRNAs, which co-precipitate with PRC1 from chromatin and found candidates that impact polycomb group protein (PcG)-regulated gene expression in vivo. A novel lncRNA from this screen, CAT7, regulates expression and polycomb group binding at the MNX1 locus during early neuronal differentiation. CAT7 contains a unique tandem repeat domain that shares high sequence similarity to a non-syntenic zebrafish analog, cat7l. Defects caused by interference of cat7l RNA during zebrafish embryogenesis were rescued by human CAT7 RNA, enhanced by interference of a PRC1 component, and suppressed by interference of a known PRC1 target gene, demonstrating cat7l genetically interacts with a PRC1. We propose a model whereby PRC1 acts in concert with specific lncRNAs and that CAT7/cat7l represents convergent lncRNAs that independently evolved to tune PRC1 repression at individual loci.

Combined administration of testosterone plus an ornithine decarboxylase inhibitor as a selective prostate-sparing anabolic therapy

Because of its anabolic effects on muscle, testosterone is being explored as a function‐promoting anabolic therapy for functional limitations associated with aging; however, concerns about testosterones adverse effects on prostate have inspired efforts to develop strategies that selectively increase muscle mass while sparing the prostate. Testosterones promyogenic effects are mediated through upregulation of follistatin. We show here that the administration of recombinant follistatin (rFst) increased muscle mass in mice, but had no effect on prostate mass. Consistent with the results of rFst administration, follistatin transgenic mice with constitutively elevated follistatin levels displayed greater muscle mass than controls, but had similar prostate weights. To elucidate signaling pathways regulated differentially by testosterone and rFst in prostate and muscle, we performed microarray analysis of mRNAs from prostate and levator ani of castrated male mice treated with vehicle, testosterone, or rFst. Testosterone and rFst shared the regulation of many transcripts in levator ani; however, in prostate, 593 transcripts in several growth‐promoting pathways were differentially expressed after testosterone treatment, while rFst showed a negligible effect with only 9 transcripts differentially expressed. Among pathways that were differentially responsive to testosterone in prostate, we identified ornithine decarboxylase (Odc1), an enzyme in polyamine biosynthesis, as a testosterone‐responsive gene that is unresponsive to rFst. Accordingly, we administered testosterone with and without α‐difluoromethylornithine (DFMO), an Odc1 inhibitor, to castrated mice. DFMO selectively blocked testosterones effects on prostate, but did not affect testosterones anabolic effects on muscle. Co‐administration of testosterone and Odc1 inhibitor presents a novel therapeutic strategy for prostate‐sparing anabolic therapy.