Differences in Mucosal Gene Expression in the Colon of Two Inbred Mouse Strains after Colonization with Commensal Gut Bacteria

The host genotype has been proposed to contribute to individually composed bacterial communities in the gut. To provide deeper insight into interactions between gut bacteria and host, we associated germ-free C3H and C57BL/10 mice with intestinal bacteria from a C57BL/10 donor mouse. Analysis of microbiota similarity between the animals with denaturing gradient gel electrophoresis revealed the development of a mouse strain-specific microbiota. Microarray-based gene expression analysis in the colonic mucosa identified 202 genes whose expression differed significantly by a factor of more than 2. Application of bioinformatics tools demonstrated that functional terms including signaling/secretion, lipid degradation/catabolism, guanine nucleotide/guanylate binding and immune response were significantly enriched in differentially expressed genes. We had a closer look at the 56 genes with expression differences of more than 4 and observed a higher expression in C57BL/10 mice of the genes coding for Tlr1 and Ang4 which are involved in the recognition and response to gut bacteria. A higher expression of Pla2g2a was detected in C3H mice. In addition, a number of interferon-inducible genes were higher expressed in C3H than in C57BL/10 mice including Gbp1, Mal, Oasl2, Ifi202b, Rtp4, Ly6g6c, Ifi27l2a, Usp18, Ifit1, Ifi44, and Ly6g indicating that interferons may play an essential role in microbiota regulation. However, genes coding for interferons, their receptors, factors involved in interferon expression regulation or signaling pathways were not differentially expressed between the two mouse strains. Taken together, our study confirms that the host genotype is involved in the establishment of host-specific bacterial communities in the gut. Based on expression differences after colonization with the same bacterial inoculum, we propose that Pla2g2a and interferon-dependent genes may contribute to this phenomenon.

DNA methyltransferase expression in the mouse germ line during periods of de novo methylation

DNA methyltransferase (DNMT) 3A and DNMT3B are both active de novo DNA methyltransferases required for development, whereas DNMT3L, which has no demonstrable methyltransferase activity, is required for methylation of imprinted genes in the oocyte. We show here that different mechanisms are used to restrict access by these proteins to their targets during germ cell development. Transcriptional control of the Dnmt3l promoter guarantees that message is low or absent except during periods of de novo activity. Use of an alternative promoter at the Dnmt3a locus produces the shorter Dnmt3a2 transcript in the germ line and postimplantation embryo only, whereas alternative splicing of the Dnmt3b transcript ensures that Dnmt3b1 is absent in the male prospermatogonia. Control of subcellular protein localization is a common theme for DNMT3A and DNMT3B, as proteins were seen in the nucleus only when methylation was occurring. These mechanisms converge to ensure that the only time that functional products from each locus are present in the germ cell nuclei is around embryonic day 17.5 in males and after birth in the growing oocytes in females.

Functional Consequences of Splicing of the Antisense Transcript COOLAIR on FLC Transcription

Antisense transcription is widespread in many genomes; however, how much is functional is hotly debated. We are investigating functionality of a set of long noncoding antisense transcripts, collectively called COOLAIR, produced at Arabidopsis FLOWERING LOCUS C (FLC). COOLAIR initiates just downstream of the major sense transcript poly(A) site and terminates either early or extends into the FLC promoter region. We now show that splicing of COOLAIR is functionally important. This was revealed through analysis of a hypomorphic mutation in the core spliceosome component PRP8. The prp8 mutation perturbs a cotranscriptional feedback mechanism linking COOLAIR processing to FLC gene body histone demethylation and reduced FLC transcription. The importance of COOLAIR splicing in this repression mechanism was confirmed by disrupting COOLAIR production and mutating the COOLAIR proximal splice acceptor site. Our findings suggest that altered splicing of a long noncoding transcript can quantitatively modulate gene expression through cotranscriptional coupling mechanisms.

A Burkholderia cenocepacia gene encoding a non-functional tyrosine phosphatase is required for the delayed maturation of the bacteria-containing vacuoles in macrophages

Burkholderia cenocepacia infects patients with cystic fibrosis. We have previously shown that B. cenocepacia can survive in macrophages within membrane vacuoles (BcCVs) that preclude fusion with the lysosome. The bacterial factors involved in B. cenocepacia intracellular survival are not fully elucidated. We report here that deletion of BCAM0628, encoding a predicted low-molecular weight protein tyrosine phosphatase (LMW-PTP) that is restricted to B. cenocepacia strains of the transmissible ET-12 clone, accelerates the maturation of the BcCVs. Compared to parental strain and deletion mutants in other LMW-PTPs that are widely conserved in Burkholderia species, a greater proportion of BcCVs containing the BCAM0628 mutant were targeted to the lysosome. Accelerated BcCV maturation was not due to reduced intracellular viability since BCAM0628 survived and replicated in macrophages similarly to the parental strain. Therefore, BCAM0628 was referred to as dpm (delayed phagosome maturation). We provide evidence that the Dpm protein is secreted during growth in vitro and upon macrophage infection. Dpm secretion requires an N-terminal signal peptide. Heterologous expression of Dpm in B. multivorans confers to this bacterium a similar phagosomal maturation delay as found with B. cenocepacia. We demonstrate that Dpm is an inactive phosphatase, suggesting that its contribution to phagosomal maturation arrest must be unrelated to tyrosine phosphatase activity.

FLT1 Genetic Variation Predisposes to Neovascular AMD in Ethnically Diverse Populations and Alters Systemic FLT1 Expression

Purpose: Current understanding of the genetic risk factors for age-related macular degeneration (AMD) is not sufficiently predictive of the clinical course. The VEGF pathway is a key therapeutic target for treatment of neovascular AMD; however, risk attributable to genetic variation within pathway genes is unclear. We sought to identify single nucleotide polymorphisms (SNPs) associated with AMD within the VEGF pathway.
Methods: Using a tagSNP, direct sequencing and meta-analysis approach within four ethnically diverse cohorts, we identified genetic risk present in FLT1, though not within other VEGF pathway genes KDR, VEGFA, or VASH1. We used ChIP and ELISA in functional analysis.
Results: The FLT1 SNPs rs9943922, rs9508034, rs2281827, rs7324510, and rs9513115 were significantly associated with increased risk of neovascular AMD. Each association was more significant after meta-analysis than in any one of the four cohorts. All associations were novel, within noncoding regions of FLT1 that do not tag for coding variants in linkage disequilibrium. Analysis of soluble FLT1 demonstrated higher expression in unaffected individuals homozygous for the FLT1 risk alleles rs9943922 (P = 0.0086) and rs7324510 (P = 0.0057). In silico analysis suggests that these variants change predicted splice sites and RNA secondary structure, and have been identified in other neovascular pathologies. These data were supported further by murine chromatin immunoprecipitation demonstrating that FLT1 is a target of Nr2e3, a nuclear receptor gene implicated in regulating an AMD pathway.
Conclusions: Although exact variant functions are not known, these data demonstrate relevancy across ethnically diverse genetic backgrounds within our study and, therefore, hold potential for global efficacy.

Functional and genetic analysis of the colon cancer network

Cancer is a complex disease that has proven to be difficult to understand on the single-gene level. For this reason a functional elucidation needs to take interactions among genes on a systems-level into account. In this study, we infer a colon cancer network from a large-scale gene expression data set by using the method BC3Net. We provide a structural and a functional analysis of this network and also connect its molecular interaction structure with the chromosomal locations of the genes enabling the definition of cis- and trans-interactions. Furthermore, we investigate the interaction of genes that can be found in close neighborhoods on the chromosomes to gain insight into regulatory mechanisms. To our knowledge this is the first study analyzing the genome-scale colon cancer network.

Functional Genomic Screen Identifies Klebsiella pneumoniae Factors Implicated in Blocking Nuclear Factor κB (NF-κB) Signaling

Klebsiella pneumoniae is etiologic agent of community-acquired and nosocomial pneumonia. It has been shown that K. pneumoniae infections are characterized by reduced early inflammatory response. Recently our group have shown that K. pneumoniae dampens the activation of inflammatory responses by antagonizing the activation of the NF-κB canonical pathway. Our results revealed that K. pneumoniae capsule (CPS) was necessary but not sufficient to attenuate inflammation. To identify additional Klebsiella factors required to dampen inflammation, we standardized and applied a high-throughput gain-on-function screen to examine a Klebsiella transposon mutant library. We identified 114 mutants that triggered the activation of NF-κB. Two gene ontology categories accounted for half of the loci identified in the screening, that of metabolism and transport (32% of the mutants), and of enveloperelated genes (17%). Characterization of the mutants revealed that the lack of the enterobactin siderophore was linked to a reduced CPS expression which in turn underlined the NF- κB activation induced by the mutant. The lipopolysaccharide (LPS) O-polysaccharide and the pullulanase (PulA) type 2 secretion system (T2SS) are required for full effectiveness of immune evasion. Importantly, these factors do not play a redundant role. The fact that LPS Opolysaccharide and T2SS mutants-induced responses were dependent on TLR2-TLR4- MyD88 activation suggested that LPS Opolysaccharide and PulA perturbed TLRdependent recognition of K. pneumoniae. Finally, we demonstrate that LPS O-polysaccharide and pulA mutants are attenuated in the pneumonia mouse model. We propose that LPS Opolysaccharide and PulA T2SS could be new targets for designing new antimicrobials. Increasing TLR-governed defence responses might provide also selective alternatives for the management of K. pneumoniae pneumonia.

Disruption of SF3B1 results in deregulated expression and splicing of key genes and pathways in myelodysplastic syndrome hematopoietic stem and progenitor cells

The splicing factor SF3B1 is the most commonly mutated gene in the myelodysplastic syndrome (MDS), particularly in patients with refractory anemia with ring sideroblasts (RARS). We investigated the functional effects of SF3B1 disruption in myeloid cell lines: SF3B1 knockdown resulted in growth inhibition, cell cycle arrest and impaired erythroid differentiation and deregulation of many genes and pathways, including cell cycle regulation and RNA processing. MDS is a disorder of the hematopoietic stem cell and we thus studied the transcriptome of CD34 + cells from MDS patients with SF3B1 mutations using RNA sequencing. Genes significantly differentially expressed at the transcript andor exon level in SF3B1 mutant compared with wild-type cases include genes that are involved in MDS pathogenesis (ASXL1 and CBL), iron homeostasis and mitochondrial metabolism (ALAS2, ABCB7 and SLC25A37) and RNA splicingprocessing (PRPF8 and HNRNPD). Many genes regulated by a DNA damage-induced BRCA1-BCLAF1-SF3B1 protein complex showed differential expressionsplicing in SF3B1 mutant cases. This is the first study to determine the target genes of SF3B1 mutation in MDS CD34 + cells. Our data indicate that SF3B1 has a critical role in MDS by affecting the expression and splicing of genes involved in specific cellular processespathways, many of which are relevant to the known RARS pathophysiology, suggesting a causal link.

Urothelial cancer gene regulatory networks inferred from large-scale RNAseq, Bead and Oligo gene expression data

BACKGROUND: Urothelial pathogenesis is a complex process driven by an underlying network of interconnected genes. The identification of novel genomic target regions and gene targets that drive urothelial carcinogenesis is crucial in order to improve our current limited understanding of urothelial cancer (UC) on the molecular level. The inference of genome-wide gene regulatory networks (GRN) from large-scale gene expression data provides a promising approach for a detailed investigation of the underlying network structure associated to urothelial carcinogenesis.

METHODS: In our study we inferred and compared three GRNs by the application of the BC3Net inference algorithm to large-scale transitional cell carcinoma gene expression data sets from Illumina RNAseq (179 samples), Illumina Bead arrays (165 samples) and Affymetrix Oligo microarrays (188 samples). We investigated the structural and functional properties of GRNs for the identification of molecular targets associated to urothelial cancer.

RESULTS: We found that the urothelial cancer (UC) GRNs show a significant enrichment of subnetworks that are associated with known cancer hallmarks including cell cycle, immune response, signaling, differentiation and translation. Interestingly, the most prominent subnetworks of co-located genes were found on chromosome regions 5q31.3 (RNAseq), 8q24.3 (Oligo) and 1q23.3 (Bead), which all represent known genomic regions frequently deregulated or aberated in urothelial cancer and other cancer types. Furthermore, the identified hub genes of the individual GRNs, e.g., HID1/DMC1 (tumor development), RNF17/TDRD4 (cancer antigen) and CYP4A11 (angiogenesis/ metastasis) are known cancer associated markers. The GRNs were highly dataset specific on the interaction level between individual genes, but showed large similarities on the biological function level represented by subnetworks. Remarkably, the RNAseq UC GRN showed twice the proportion of significant functional subnetworks. Based on our analysis of inferential and experimental networks the Bead UC GRN showed the lowest performance compared to the RNAseq and Oligo UC GRNs.

CONCLUSION: To our knowledge, this is the first study investigating genome-scale UC GRNs. RNAseq based gene expression data is the data platform of choice for a GRN inference. Our study offers new avenues for the identification of novel putative diagnostic targets for subsequent studies in bladder tumors.

Identification of lactic acid bacteria strains modulating incretin hormone secretion and gene expression in enteroendocrine cells

Glucagon-like peptide-1 (GLP-1) and glucose-dependent insulinotropic polypeptide (GIP) are incretin hormones released from intestinal enteroendocrine (EE) cells and have well-established glucose-lowering actions. Lactic acid bacteria (LAB) colonise the human intestine, but it is unknown whether LAB and EE cells interact. Acute co-culture of LAB with EE cells showed that certain LAB strains elicit GLP-1 and GIP secretion (13-194-fold) and upregulate their gene expression. LAB-induced incretin hormone secretion did not appear to involve nutrient mechanisms, nor was there any evidence of cytolysis. Instead PCR array studies implicated signalling agents of the toll-like receptor system, e.g. adaptor protein MyD88 was decreased 23-fold and cell surface antigen CD14 was increased 17-fold. Mechanistic studies found that blockade of MyD88 triggered significant GLP-1 secretion. Furthermore, blocking of CD14 completely attenuated LAB-induced secretion. A recent clinical trial clearly shows that LAB have potential for alleviating type 2 diabetes, and further characterisation of this bioactivity is warranted.

Prokineticin ligands and receptors are expressed in the human fetal ovary and regulate germ cell expression of COX2

CONTEXT: Fetal ovarian development and primordial follicle formation underpin future female fertility. Prokineticin (PROK) ligands regulate cell survival, proliferation and angiogenesis in adult reproductive tissues including the ovary. However, their expression and function during fetal ovarian development remains unclear.

OBJECTIVE: To investigate expression and localization of the PROK ligands, receptors and their downstream transcriptional targets in the human fetal ovary.

SETTING: This study was conducted at the University of Edinburgh.

PARTICIPANTS: Ovaries were collected from 37 morphologically normal human fetuses.

DESIGN AND MAIN OUTCOME MEASURES: mRNA and protein expression of PROK ligands and receptors was determined in human fetal ovaries using qRT-PCR, immunoblotting and immunohistochemistry. Functional studies were performed using a human germ tumour cell line (TCam-2) stably transfected with PROKR1.

RESULTS: Expression of PROK1 and PROKR1 was significantly higher in mid-gestation ovaries (17-20 weeks) than at earlier gestations (8-11 and 14-16 weeks). PROK2 significantly increased across the gestations examined. PROKR2 expression remained unchanged. PROK ligand and receptor proteins were predominantly localised to germ cells (including oocytes within primordial follicles) and endothelial cells, indicating these cell types to be the targets of PROK signalling in the human fetal ovary. PROK1 treatment of a germ cell line stably-expressing PROKR1 resulted in ERK phosphorylation, and elevated COX2 expression.

CONCLUSIONS: Developmental changes in expression and regulation of COX2 and pERK by PROK1 suggest that PROK ligands may be novel regulators of germ cell development in the human fetal ovary, interacting within a network of growth and survival factors prior to primordial follicle formation.

Convergent genetic and expression data implicate immunity in Alzheimer's disease

Background: Late-onset Alzheimer's disease (AD) is heritable with 20 genes showing genome-wide association in the International Genomics of Alzheimer's Project (IGAP). To identify the biology underlying the disease, we extended these genetic data in a pathway analysis.

Methods: The ALIGATOR and GSEA algorithms were used in the IGAP data to identify associated functional pathways and correlated gene expression networks in human brain.

Results: ALIGATOR identified an excess of curated biological pathways showing enrichment of association. Enriched areas of biology included the immune response (P = 3.27 X 10(-12) after multiple testing correction for pathways), regulation of endocytosis (P = 1.31 X 10(-11)), cholesterol transport (P = 2.96 X 10(-9)), and proteasome-ubiquitin activity (P = 1.34 X 10(-6)). Correlated gene expression analysis identified four significant network modules, all related to the immune response (corrected P = .002-.05).

Conclusions: The immime response, regulation of endocytosis, cholesterol transport, and protein ubiquitination represent prime targets for AD therapeutics. (C) 2015 Published by Elsevier Inc. on behalf of The Alzheimer's Association.

Uncovering BRCA1-regulated signalling pathways by microarray-based expression profiling

The introduction of microarray technology to the scientific and medical communities has dramatically changed the way in which we now address basic biomedical questions. Expression profiling using microarrays facilitates an experimental approach where alterations in the transcript level of entire transcriptomes can be simultaneously assayed in response to defined stimuli. We have used microarray analysis to identify downstream transcriptional targets of the BRCA1 (Breast Cancer 1) tumour-suppressor gene as a means of defining its function. BRCA1 has been implicated in the predisposition to early onset breast and ovarian cancer and while its exact function remains to be defined, roles in DNA repair, cell-cycle control and transcriptional regulation have been implied. In the current study we have generated cell lines with tetracycline-regulated, inducible expression of BRCA1 as a tool to identify genes, which might represent important effectors of BRCA1 function. Oligonucleotide array-based expression profiling identified a number of genes that were upregulated at various times following inducible expression of BRCA1 including the DNA damage-responsive gene GADD45 (Growth Arrest after DNA Damage). Identified targets were confirmed by Northern blot analysis and their functional significance as BRCA1 targets examined.