Genome-wide analysis of mammalian promoter architecture and evolution

Mammalian promoters can be separated into two classes, conserved TATA box-enriched promoters, which initiate at a welldefined site, and more plastic, broad and evolvable CpG-rich promoters. We have sequenced tags corresponding to several hundred thousand transcription start sites (TSSs) in the mouse and human genomes, allowing precise analysis of the sequence architecture and evolution of distinct promoter classes. Different tissues and families of genes differentially use distinct types of promoters. Our tagging methods allow quantitative analysis of promoter usage in different tissues and show that differentially regulated alternative TSSs are a common feature in protein-coding genes and commonly generate alternative N termini. Among the TSSs, we identified new start sites associated with the majority of exons and with 3' UTRs. These data permit genome-scale identification of tissue-specific promoters and analysis of the cis-acting elements associated with them.

Transcriptional network dynamics in macrophage activation

Transcriptional regulatory networks govern cell differentiation and the cellular response to external stimuli. However, mammalian model systems have not yet been accessible for network analysis. Here, we present a genome-wide network analysis of the transcriptional regulation underlying the mouse macrophage response to bacterial lipopolysaccharide (LPS). Key to uncovering the network structure is our combination of time-series cap analysis of gene expression with in silico prediction of transcription factor binding sites. By integrating microarray and qPCR time-series expression data with a promoter analysis, we find dynamic subnetworks that describe how signaling pathways change dynamically during the progress of the macrophage LPS response, thus defining regulatory modules characteristic of the inflammatory response. In particular, our integrative analysis enabled us to suggest novel roles for the transcription factors ATF-3 and NRF-2 during the inflammatory response. We believe that our system approach presented here is applicable to understanding cellular differentiation in higher eukaryotes. (c) 2006 Elsevier Inc. All rights reserved.

Retrotransposon expression profiling: Unveiling the hidden SINE transcriptome through Next-Generation Sequencing data analysis

Of the ~1.7 million SINE elements in the human genome, only a tiny number are estimated to be active in transcription by RNA polymerase (Pol) III. Tracing the individual loci from which SINE transcripts originate is complicated by their highly repetitive nature. By exploiting RNA-Seq datasets and unique SINE DNA sequences, we devised a bioinformatic pipeline allowing us to identify Pol III-dependent transcripts of individual SINE elements. When applied to ENCODE transcriptomes of seven human cell lines, this search strategy identified ~1300 Alu loci and ~1100 MIR loci corresponding to detectable transcripts, with ~120 and ~60 respectively Alu and MIR loci expressed in at least three cell lines. In vitro transcription of selected SINEs did not reflect their in vivo expression properties, and required the native 5’-flanking region in addition to internal promoter. We also identified a cluster of expressed AluYa5-derived transcription units, juxtaposed to snaR genes on chromosome 19, formed by a promoter-containing left monomer fused to an Alu-unrelated downstream moiety. Autonomous Pol III transcription was also revealed for SINEs nested within Pol II-transcribed genes raising the possibility of an underlying mechanism for Pol II gene regulation by SINE transcriptional units. Moreover the application of our bioinformatic pipeline to both RNA-seq data of cells subjected to an in vitro pro-oncogenic stimulus and of in vivo matched tumor and non-tumor samples allowed us to detect increased Alu RNA expression as well as the source loci of such deregulation. The ability to investigate SINE transcriptomes at single-locus resolution will facilitate both the identification of novel biologically relevant SINE RNAs and the assessment of SINE expression alteration under pathological conditions.

Transcriptome analysis of a respiratory Saccharomyces cerevisiae strain suggests the expression of its phenotype is glucose insensitive and predominantly controlled by Hap4, Cat8 and Mig1

BACKGROUND: We previously described the first respiratory Saccharomyces cerevisiae strain, KOY.TM6*P, by integrating the gene encoding a chimeric hexose transporter, Tm6*, into the genome of an hxt null yeast. Subsequently we transferred this respiratory phenotype in the presence of up to 50 g/L glucose to a yeast strain, V5 hxt1-7Delta, in which only HXT1-7 had been deleted. In this study, we compared the transcriptome of the resultant strain, V5.TM6*P, with that of its wild-type parent, V5, at different glucose concentrations. RESULTS: cDNA array analyses revealed that alterations in gene expression that occur when transitioning from a respiro-fermentative (V5) to a respiratory (V5.TM6*P) strain, are very similar to those in cells undergoing a diauxic shift. We also undertook an analysis of transcription factor binding sites in our dataset by examining previously-published biological data for Hap4 (in complex with Hap2, 3, 5), Cat8 and Mig1, and used this in combination with verified binding consensus sequences to identify genes likely to be regulated by one or more of these. Of the induced genes in our dataset, 77% had binding sites for the Hap complex, with 72% having at least two. In addition, 13% were found to have a binding site for Cat8 and 21% had a binding site for Mig1. Unexpectedly, both the up- and down-regulation of many of the genes in our dataset had a clear glucose dependence in the parent V5 strain that was not present in V5.TM6*P. This indicates that the relief of glucose repression is already operable at much higher glucose concentrations than is widely accepted and suggests that glucose sensing might occur inside the cell. CONCLUSION: Our dataset gives a remarkably complete view of the involvement of genes in the TCA cycle, glyoxylate cycle and respiratory chain in the expression of the phenotype of V5.TM6*P. Furthermore, 88% of the transcriptional response of the induced genes in our dataset can be related to the potential activities of just three proteins: Hap4, Cat8 and Mig1. Overall, our data support genetic remodelling in V5.TM6*P consistent with a respiratory metabolism which is insensitive to external glucose concentrations.

Monitoring of transcriptional dynamics of HIF and NFκB activities

Genetic experiments over the last few decades have identified many regulatory proteins critical for DNA transcription. The dynamics of their transcriptional activities shape the differential expression of the genes they control. Here we describe a simple method, based on the secreted luciferase, to measure the activities of two transcription factors NF?B and HIF. This technique can effectively monitor dynamics of transcriptional events in a population of cells and be up-scaled for high-throughput screening and promoter analysis, making it ideal for data-demanding applications such as mathematical modelling.

Pichia pastoris regulates its gene-specific response to different carbon sources at the transcriptional, rather than the translational, level

Background: The methylotrophic, Crabtree-negative yeast Pichia pastoris is widely used as a heterologous protein production host. Strong inducible promoters derived from methanol utilization genes or constitutive glycolytic promoters are typically used to drive gene expression. Notably, genes involved in methanol utilization are not only repressed by the presence of glucose, but also by glycerol. This unusual regulatory behavior prompted us to study the regulation of carbon substrate utilization in different bioprocess conditions on a genome wide scale. Results: We performed microarray analysis on the total mRNA population as well as mRNA that had been fractionated according to ribosome occupancy. Translationally quiescent mRNAs were defined as being associated with single ribosomes (monosomes) and highly-translated mRNAs with multiple ribosomes (polysomes). We found that despite their lower growth rates, global translation was most active in methanol-grown P. pastoris cells, followed by excess glycerol- or glucose-grown cells. Transcript-specific translational responses were found to be minimal, while extensive transcriptional regulation was observed for cells grown on different carbon sources. Due to their respiratory metabolism, cells grown in excess glucose or glycerol had very similar expression profiles. Genes subject to glucose repression were mainly involved in the metabolism of alternative carbon sources including the control of glycerol uptake and metabolism. Peroxisomal and methanol utilization genes were confirmed to be subject to carbon substrate repression in excess glucose or glycerol, but were found to be strongly de-repressed in limiting glucose-conditions (as are often applied in fed batch cultivations) in addition to induction by methanol. Conclusions: P. pastoris cells grown in excess glycerol or glucose have similar transcript profiles in contrast to S. cerevisiae cells, in which the transcriptional response to these carbon sources is very different. The main response to different growth conditions in P. pastoris is transcriptional; translational regulation was not transcript-specific. The high proportion of mRNAs associated with polysomes in methanol-grown cells is a major finding of this study; it reveals that high productivity during methanol induction is directly linked to the growth condition and not only to promoter strength.

Pseudomonas aeruginosa AMPR transcriptional regulatory network

In Enterobacteriaceae, the transcriptional regulator AmpR, a member of the LysR family, regulates the expression of a chromosomal β-lactamase AmpC. The regulatory repertoire of AmpR is broader in Pseudomonas aeruginosa, an opportunistic pathogen responsible for numerous acute and chronic infections including cystic fibrosis. Previous studies showed that in addition to regulating ampC, P. aeruginosa AmpR regulates the sigma factor AlgT/U and production of some quorum sensing (QS)-regulated virulence factors. In order to better understand the ampR regulon, the transcriptional profiles generated using DNA microarrays and RNA-Seq of the prototypic P. aeruginosa PAO1 strain with its isogenic ampR deletion mutant, PAOΔampR were analyzed. Transcriptome analysis demonstrates that the AmpR regulon is much more extensive than previously thought influencing the differential expression of over 500 genes. In addition to regulating resistance to β-lactam antibiotics via AmpC, AmpR also regulates non-β-lactam antibiotic resistance by modulating the MexEF-OprN efflux pump. Virulence mechanisms including biofilm formation, QS-regulated acute virulence, and diverse physiological processes such as oxidative stress response, heat-shock response and iron uptake are AmpR-regulated. Real-time PCR and phenotypic assays confirmed the transcriptome data. Further, Caenorhabditis elegans model demonstrates that a functional AmpR is required for full pathogenicity of P. aeruginosa. AmpR, a member of the core genome, also regulates genes in the regions of genome plasticity that are acquired by horizontal gene transfer. The extensive AmpR regulon included other transcriptional regulators and sigma factors, accounting for the extensive AmpR regulon. Gene expression studies demonstrate AmpR-dependent expression of the QS master regulator LasR that controls expression of many virulence factors. Using a chromosomally tagged AmpR, ChIP-Seq studies show direct AmpR binding to the lasR promoter. The data demonstrates that AmpR functions as a global regulator in P. aeruginosa and is a positive regulator of acute virulence while negatively regulating chronic infection phenotypes. In summary, my dissertation sheds light on the complex regulatory circuit in P. aeruginosa to provide a better understanding of the bacterial response to antibiotics and how the organism coordinately regulates a myriad of virulence factors.

Single-Cell Transcriptome Analysis of Olfactory Sensory Neurons

Olfactory sensory neurons (OSNs), which detect a myriad of odorants, are known to express one allele of one olfactory receptor (OR) gene (Olfr) from the largest gene family in the mammalian genome. The OSNs expressing the same OR project their axons to the main olfactory bulb where they converge to form glomeruli. This “One neuron-one receptor rule” makes the olfactory epithelium (OE), which consists of a vast number of OSNs expressing unique ORs, one of the most heterogeneous cell populations. However, the mechanism of how the single OR allele is chosen remains unclear along with the question of whether one OSN only expresses a single OR gene, a hypothesis that has not been rigorously verified while we performed the experiments. Moreover, failure of axonal targeting to single glomerulus was observed in MeCP2 deficient OSNs where delayed development was proposed as an explanation for the phenotype. How Mecp2 mutation caused this aberrant targeting is not entirely understood.

In this dissertation, we explored the transcriptomes of single and mature OSNs by single-cell RNA-Seq to reveal their heterogeneity and further studied the OR gene expression from these isolated OSNs. The singularity of sequenced OSNs was ensured by the observation of monoallelic expression of X-linked genes from the hybrid samples from crosses between mice of different strains where strain-specific polymorphisms could be used to track the allelic origins of SNP-containing reads. The clustering of expression profiles from triplicates that originated from the same cell assured that the transcriptomic identities of OSNs were maintained through the experimental process. The average gene expression profiles of sequenced OSNs correlated well to the conventional transcriptome data of FACS-sorted Omp-positive cells, and the top-ranked expression of OR was conceded in the single-OSN transcriptomes. While exploring cellular diversity, in addition to OR genes, we revealed nearly 200 differentially expressed genes among the sequenced OSNs in this study. Among the 36 sequenced OSNs, eight cells (22.2%) showed multiple OR gene expression and the presences of additional ORs were not restricted to the neighbor loci that shared the transcriptional effect of the primary OR expression, suggesting that the “One neuron-one receptor rule” might not be strictly true at the transcription level. All of the inferable ORs, including additional co-expressed ORs, were shown to be monoallelic. Our sequencing of 21 Mecp2308 mutant OSNs, of which 62% expressed more than one OR genes, and the expression levels of the additional ORs were significantly higher than those in the wild-type, suggested that MeCP2 plays a role in the regulation of singular OR gene expression. Dual label in situ hybridization along with the sequence data revealed that dorsal and ventral ORs were co-expressed in the same Mecp2 mutant OSN, further implying that MeCP2 might be involved in regulation of OR territories in the OE. Our results suggested a new role of MeCP2 in OR gene choice and ratified that this multiple-OR expression caused by Mecp2 mutation did not accompany delayed OSN development that has been observed in the previous studies on the Mecp2 mutants.

Transcriptomic Analysis of the Host Response and Innate Resilience to Enterotoxigenic Escherichia coli Infection in Humans.

BACKGROUND: Enterotoxigenic Escherichia coli (ETEC) is a globally prevalent cause of diarrhea. Though usually self-limited, it can be severe and debilitating. Little is known about the host transcriptional response to infection. We report the first gene expression analysis of the human host response to experimental challenge with ETEC. METHODS: We challenged 30 healthy adults with an unattenuated ETEC strain, and collected serial blood samples shortly after inoculation and daily for 8 days. We performed gene expression analysis on whole peripheral blood RNA samples from subjects in whom severe symptoms developed (n = 6) and a subset of those who remained asymptomatic (n = 6) despite shedding. RESULTS: Compared with baseline, symptomatic subjects demonstrated significantly different expression of 406 genes highlighting increased immune response and decreased protein synthesis. Compared with asymptomatic subjects, symptomatic subjects differentially expressed 254 genes primarily associated with immune response. This comparison also revealed 29 genes differentially expressed between groups at baseline, suggesting innate resilience to infection. Drug repositioning analysis identified several drug classes with potential utility in augmenting immune response or mitigating symptoms. CONCLUSIONS: There are statistically significant and biologically plausible differences in host gene expression induced by ETEC infection. Differential baseline expression of some genes may indicate resilience to infection.

Analysis of expression quantitative trait loci in D. melanogaster

Expression Quantitative Trait Loci (eQTL) analysis allows for the identification of genetic variation associated with variation in gene expression. It is often unclear however, which of the associated variants are causal, and by what mechanism. Integrating functional genomic data with eQTL data can provide insight into the impact of natural variation in the population, and the nature of the transcriptional machinery itself. In this thesis, I integrate functional genomic data with eQTL data derived from both 5’ CAGE and 3’ TagXseq expression assays, in developing embryos. I first use both datasets to analyse the transcription landscape in embryonic D., melanogaster, and then carry out an analysis of sequence motifs associated with transcription factor binding sites, promoters, and 3’ polyadenylation sites. Finally, I integrate functional genomic data, including these novel sequence motifs, to shed light on the mechanisms of gene expression variation in D.,melanogaster. I am able to demonstrate that some variants effecting gene regulation in Drosophila are found within haplotypes which buffer their effects.

An in vitro investigation on the cytotoxic and nuclear receptor transcriptional activity of the mycotoxins fumonisin B1 and beauvericin

Fumonisin B1 (FB1) and beauvericin (BEA) are secondary metabolites of filamentous fungi, which under appropriate temperature and humidity conditions may develop on various foods and feeds. To date few studies have been performed to evaluate the toxicological and endocrine disrupting effects of FB1 and BEA. The present study makes use of various in vitro bioassays including; oestrogen, androgen, progestagen and glucocorticoid reporter gene assays (RGAs) for the study of nuclear receptor transcriptional activity, the thiazolyl blue tetrazolium bromide (MTT) assay to monitor cytotoxicity and high content analysis (HCA) for the detection of pre-lethal toxicity in the RGA and Caco-2 human colon adenocarcinoma cells. At the receptor level, 0.001-10μM BEA or FB1 did not induce any agonist responses in the RGAs. However at non-cytotoxic concentrations, an antagonistic effect was exhibited by FB1 on the androgen nuclear receptor transcriptional activity at 10μM and BEA on the progestagen and glucocorticoid receptors at 1μM. MTT analysis showed no decrease in cell viability at any concentration of FB1, whereas BEA showed a significant decrease in viability at 10μM. HCA analysis confirmed that the reduction in the progestagen receptor transcriptional activity at 1μM BEA was not due to pre-lethal toxicity. In addition, BEA (10μM) induced significant toxicity in both the TM-Luc (progestagen responsive) and Caco-2 cells.

Transcriptional profiling of chronic clinical hepatic schistosomiasis japonica indicates reduced metabolism and immune responses

Schistosomiasis is a significant cause of human morbidity and mortality. We performed a genome-wide transcriptional survey of liver biopsies obtained from Chinese patients with chronic schistosomiasis only, or chronic schistosomiasis with a current or past history of viral hepatitis B. Both disease groups were compared with patients with no prior history or indicators of any liver disease. Analysis showed in the main, downregulation in gene expression, particularly those involved in signal transduction via EIF2 signalling and mTOR signalling, as were genes associated with cellular remodelling. Focusing on immune associated pathways, genes were generally downregulated. However, a set of three genes associated with granulocytes, MMP7, CLDN7, CXCL6 were upregulated. Differential gene profiles unique to schistosomiasis included the gene Granulin which was decreased despite being generally considered a marker for liver disease, and IGBP2 which is associated with increased liver size, and was the most upregulated gene in schistosomiasis only patients, all of which presented with hepatomegaly. The unique features of gene expression, in conjunction with previous reports in the murine model of the cellular composition of granulomas, granuloma formation and recovery, provide an increased understanding of the molecular immunopathology and general physiological processes underlying hepatic schistosomiasis.

Transcriptional profiling of the oesophageal gland region of male worms of Schistosoma mansoni

The intestinal tract of schistosomes opens at the mouth and leads into the foregut or oesophageal region that is lined with syncytium continuous with the apical cytoplasm of the tegument. The oesophagus is surrounded by a specialised gland, the oesophageal gland. This gland releases materials into the lumen of the oesophagus and the region is thought to initiate the lysis of erythrocytes and neutralisation of immune effectors of the host. The oesophageal region is present in the early invasive schistosomulum, a stage potentially targetable by anti-schistosome vaccines. We used a 44k oligonucleotide microarray to identify highly up-regulated genes in microdissected frozen sections of the oesophageal gland of male worms of S. mansoni. We show that 122 genes were up-regulated 2-fold or higher in the oesophageal gland compared with a whole male worm tissue control. The enriched genes included several associated with lipid metabolism and transmembrane transport as well as some micro-exon genes. Since the oesophageal gland is important in the initiation of digestion and the fact that it develops early after invasion of the mammalian host, further study of selected highly up-regulated functionally important genes in this tissue may reveal new anti-schistosome intervention targets for schistosomiasis control.

Functional analysis of the 5′ flanking region of the human G6PC3 gene: regulation of promoter activity by glucose, pyruvate, AMP kinase and the pentose phosphate pathway

G6PC3 is a widely expressed isoform of glucose-6-phosphatase, found in many foetal and adult tissues. Mutations in this gene cause developmental abnormalities and severe neutropenia due to abolition of glucose recycling between the cytoplasm and endoplasmic reticulum. Low G6PC3 expression as a result of promoter polymorphisms or dysregulation could produce similar outcomes. Here we investigated the regulation of human G6PC3 promoter activity. HeLa and H4IIE cells were transiently transfected with G6PC3 promoter coupled to the firefly luciferase gene, and promoter activity was measured by dual luciferase assay. Activity was highest in a 453 bp segment of the G6PC3 promoter, from − 455 to − 3 relative to the transcriptional start site. This promoter was unresponsive to glucostatic hormones. Its activity increased significantly between 1 and 5.5 mM glucose, and was not elevated further by glucose concentrations up to 25 mM. Pyruvate increased its activity, but β-hydroxybutyrate and sodium acetate did not. Promoter activity was reduced by inhibitors of hexokinase, glyceraldehyde phosphate dehydrogenase and the oxidative branch of the pentose phosphate pathway, but not by a transketolase inhibitor. Deletion of two adjacent Enhancer-boxes (− 274 to − 279 and − 299 to − 304) reduced promoter activity and abolished the glucose effect, suggesting they could function as a glucose response element. Deletion of an additional downstream 140 bp (− 140 to − 306) restored activity, but not the glucose response, suggesting the presence of repressor elements in this region. 5-Aminoimidazole-4-carboxamide 1-β-d-ribofuranoside (AICAR) reduced promoter activity, showing dependence on AMP-kinase. Regulation of the G6PC3 promoter is thus radically different to that of the hepatic isoform, G6PC. It is sensitive to carbohydrate, but not to fatty acid metabolites, and at much lower physiological concentrations. Based on these findings, we speculate that reduced G6PC3 expression could occur during hypoglycemic episodes in vivo, which are common in utero and in the postnatal period. If such episodes lower G6PC3 expression they could place the foetus or infant at risk of impaired immune function and development, and this possibility requires further examination both in vitro and in vivo.