Insights into the prevalence and underlying causes of clonal variation through transcriptomic analysis in Pichia pastoris.

Clonal variation, wherein a range of specific productivities of secreted proteins are observed from supposedly identical transformants, is an accepted aspect of working with Pichia pastoris. It means that a significant number of transformants need to be tested to obtain a representative sample, and in commercial protein production, companies regularly screen thousands of transformants to select for the highest secretor. Here, we have undertaken a detailed investigation of this phenomenon by characterising clones transformed with the human serum albumin gene. The titers of nine clones, each containing a single copy of the human serum albumin gene (identified by qPCR), were measured and the clones grouped into three categories, namely, high-, mid- and low-level secretors. Transcriptomic analysis, using microarrays, showed that no regulatory patterns consistently correlated with titer, suggesting that the causes of clonal variation are varied. However, a number of physiological changes appeared to underlie the differences in titer, suggesting there is more than one biochemical signature for a high-secreting strain. An anomalous low-secreting strain displaying high transcript levels that appeared to be nutritionally starved further emphasises the complicated nature of clonal variation.

Core binding factor (CBF) is required for Epstein-Barr virus EBNA3 proteins to regulate target gene expression.

ChIP-seq performed on lymphoblastoid cell lines (LCLs), expressing epitope-tagged EBNA3A, EBNA3B or EBNA3C from EBV-recombinants, revealed important principles of EBNA3 binding to chromatin. When combined with global chromatin looping data, EBNA3-bound loci were found to have a singular character, each directly associating with either EBNA3-repressed or EBNA3-activated genes, but not with both. EBNA3A and EBNA3C showed significant association with repressed and activated genes. Significant direct association for EBNA3B loci could only be shown with EBNA3B-repressed genes. A comparison of EBNA3 binding sites with known transcription factor binding sites in LCL GM12878 revealed substantial co-localization of EBNA3s with RUNX3-a protein induced by EBV during B cell transformation. The beta-subunit of core binding factor (CBFß), that heterodimerizes with RUNX3, could co-immunoprecipitate robustly EBNA3B and EBNA3C, but only weakly EBNA3A. Depletion of either RUNX3 or CBFß with lentivirus-delivered shRNA impaired epitope-tagged EBNA3B and EBNA3C binding at multiple regulated gene loci, indicating a requirement for CBF heterodimers in EBNA3 recruitment during target-gene regulation. ShRNA-mediated depletion of CBFß in an EBNA3C-conditional LCL confirmed the role of CBF in the regulation of EBNA3C-induced and -repressed genes. These results reveal an important role for RUNX3/CBF during B cell transformation and EBV latency that was hitherto unexplored.

Deciphering the metabolic response of Mycobacterium tuberculosis to nitrogen stress.

A key component to the success of Mycobacterium tuberculosis as a pathogen is the ability to sense and adapt metabolically to the diverse range of conditions encountered in vivo, such as oxygen tension, environmental pH and nutrient availability. Although nitrogen is an essential nutrient for every organism, little is known about the genes and pathways responsible for nitrogen assimilation in M. tuberculosis. In this study we have used transcriptomics and chromatin immunoprecipitation and high-throughput sequencing to address this. In response to nitrogen starvation, a total of 185 genes were significantly differentially expressed (96 up-regulated and 89 down regulated; 5% genome) highlighting several significant areas of metabolic change during nitrogen limitation such as nitrate/nitrite metabolism, aspartate metabolism and changes in cell wall biosynthesis. We identify GlnR as a regulator involved in the nitrogen response, controlling the expression of at least 33 genes in response to nitrogen limitation. We identify a consensus GlnR binding site and relate its location to known transcriptional start sites. We also show that the GlnR response regulator plays a very different role in M. tuberculosis to that in non-pathogenic mycobacteria, controlling genes involved in nitric oxide detoxification and intracellular survival instead of genes involved in nitrogen scavenging.

Translational arrest due to cytoplasmic redox stress delays adaptation to growth on methanol and heterologous protein expression in a typical fed-batch culture of Pichia pastoris.

RESULTS: We have followed a typical fed-batch induction regime for heterologous protein production under the control of the AOX1 promoter using both microarray and metabolomic analysis. The genetic constructs involved 1 and 3 copies of the TRY1 gene, encoding human trypsinogen. In small-scale laboratory cultures, expression of the 3 copy-number construct induced the unfolded protein response (UPR) sufficiently that titres of extracellular trypsinogen were lower in the 3-copy construct than with the 1-copy construct. In the fed-batch-culture, a similar pattern was observed, with higher expression from the 1-copy construct, but in this case there was no significant induction of UPR with the 3-copy strain. Analysis of the microarray and metabolomic information indicates that the 3-copy strain was undergoing cytoplasmic redox stress at the point of induction with methanol. In this Crabtree-negative yeast, this redox stress appeared to delay the adaptation to growth on methanol and supressed heterologous protein production, probably due to a block in translation. CONCLUSION: Although redox imbalance as a result of artificially imposed hypoxia has previously been described, this is the first time that it has been characterised as a result of a transient metabolic imbalance and shown to involve a stress response which can lead to translational arrest. Without detailed analysis of the underlying processes it could easily have been mis-interpreted as secretion stress, transmitted through the UPR.

Nitrogen stress response and stringent response are coupled in Escherichia coli.

Assimilation of nitrogen is an essential process in bacteria. The nitrogen regulation stress response is an adaptive mechanism used by nitrogen-starved Escherichia coli to scavenge for alternative nitrogen sources and requires the global transcriptional regulator NtrC. In addition, nitrogen-starved E. coli cells synthesize a signal molecule, guanosine tetraphosphate (ppGpp), which serves as an effector molecule of many processes including transcription to initiate global physiological changes, collectively termed the stringent response. The regulatory mechanisms leading to elevated ppGpp levels during nutritional stresses remain elusive. Here, we show that transcription of relA, a key gene responsible for the synthesis of ppGpp, is activated by NtrC during nitrogen starvation. The results reveal that NtrC couples these two major bacterial stress responses to manage conditions of nitrogen limitation, and provide novel mechanistic insights into how a specific nutritional stress leads to elevating ppGpp levels in bacteria.

Combinatorial stress responses: direct coupling of two major stress responses in Escherichia coli.

Nitrogen is an essential element for all life, and this is no different for the bacterial cell. Numerous cellular macromolecules contain nitrogen, including proteins, nucleic acids and cell wall components. In Escherichia coli and related bacteria, the nitrogen stress (Ntr) response allows cells to rapidly sense and adapt to nitrogen limitation by scavenging for alternative nitrogen sources through the transcriptional activation of transport systems and catabolic and biosynthetic operons by the global transcriptional regulator NtrC. Nitrogen-starved bacterial cells also synthesize the (p)ppGpp effector molecules of a second global bacterial stress response - the stringent response. Recently, we showed that the transcription of relA, the gene which encodes the major (p)ppGpp synthetase in E. coli, is activated by NtrC during nitrogen starvation. Our results revealed that in E. coli and related bacteria, NtrC functions in combinatorial stress and serves to couple two major stress responses, the Ntr response and stringent response.

Deciphering the response of Mycobacterium smegmatis to nitrogen stress using bipartite active modules.

BACKGROUND: The ability to adapt to environments with fluctuating nutrient availability is vital for bacterial survival. Although essential for growth, few nitrogen metabolism genes have been identified or fully characterised in mycobacteria and nitrogen stress survival mechanisms are unknown. RESULTS: A global transcriptional analysis of the mycobacterial response to nitrogen stress, showed a significant change in the differential expression of 16% of the Mycobacterium smegmatis genome. Gene expression changes were mapped onto the metabolic network using Active Modules for Bipartite Networks (AMBIENT) to identify metabolic pathways showing coordinated transcriptional responses to the stress. AMBIENT revealed several key features of the metabolic response not identified by KEGG enrichment alone. Down regulated reactions were associated with the general reduction in cellular metabolism as a consequence of reduced growth rate. Up-regulated modules highlighted metabolic changes in nitrogen assimilation and scavenging, as well as reactions involved in hydrogen peroxide metabolism, carbon scavenging and energy generation. CONCLUSIONS: Application of an Active Modules algorithm to transcriptomic data identified key metabolic reactions and pathways altered in response to nitrogen stress, which are central to survival under nitrogen limiting environments.

Genome wide analysis of the complete GlnR nitrogen-response regulon in Mycobacterium smegmatis.

BACKGROUND: Nitrogen is an essential element for bacterial growth and an important component of biological macromolecules. Consequently, responding to nitrogen limitation is critical for bacterial survival and involves the interplay of signalling pathways and transcriptional regulation of nitrogen assimilation and scavenging genes. In the soil dwelling saprophyte Mycobacterium smegmatis the OmpR-type response regulator GlnR is thought to mediate the transcriptomic response to nitrogen limitation. However, to date only ten genes have been shown to be in the GlnR regulon, a vastly reduced number compared to other organisms. RESULTS: We investigated the role of GlnR in the nitrogen limitation response and determined the entire GlnR regulon, by combining expression profiling of M. smegmatis wild type and glnR deletion mutant, with GlnR-specific chromatin immunoprecipitation and high throughput sequencing. We identify 53 GlnR binding sites during nitrogen limitation that control the expression of over 100 genes, demonstrating that GlnR is the regulator controlling the assimilation and utilisation of nitrogen. We also determine a consensus GlnR binding motif and identify key residues within the motif that are required for specific GlnR binding. CONCLUSIONS: We have demonstrated that GlnR is the global nitrogen response regulator in M. smegmatis, directly regulating the expression of more than 100 genes. GlnR controls key nitrogen stress survival processes including primary nitrogen metabolism pathways, the ability to utilise nitrate and urea as alternative nitrogen sources, and the potential to use cellular components to provide a source of ammonium. These studies further our understanding of how mycobacteria survive nutrient limiting conditions.

A non-bacterial transcription factor inhibits bacterial transcription by a multipronged mechanism.

The process of transcription initiation is the major target for regulation of gene expression in bacteria and is performed by a multi-subunit RNA polymerase enzyme (RNAp). A complex network of regulatory elements controls the activity of the RNAp to fine-tune transcriptional output. Thus, RNAp is a nexus for controlling bacterial gene expression at the transcription level. Many bacteriophages, viruses that infect bacteria, encode transcription factors that specifically target and modulate the activity of the host RNAp and, thereby, facilitate the acquisition of the host bacteria by the phage. Here, we describe the modus operandi of a T7 bacteriophage-encoded small protein called Gp2 and define Gp2 as a non-bacterial regulator of bacterial transcription.

Adenylylation of mycobacterial Glnk (PII) protein is induced by nitrogen limitation.

PII proteins are pivotal regulators of nitrogen metabolism in most prokaryotes, controlling the activities of many targets, including nitrogen assimilation enzymes, two component regulatory systems and ammonium transport proteins. Escherichia coli contains two PII-like proteins, PII (product of glnB) and GlnK, both of which are uridylylated under nitrogen limitation at a conserved Tyrosine-51 residue by GlnD (a uridylyl transferase). PII-uridylylation in E. coli controls glutamine synthetase (GS) adenylylation by GlnE and mediates the NtrB/C transcriptomic response. Mycobacteria contain only one PII protein (GlnK) which in environmental Actinomycetales is adenylylated by GlnD under nitrogen limitation. However in mycobacteria, neither the type of GlnK (PII) covalent modification nor its precise role under nitrogen limitation is known. In this study, we used LC-Tandem MS to analyse the modification state of mycobacterial GlnK (PII), and demonstrate that during nitrogen limitation GlnK from both non-pathogenic Mycobacterium smegmatis and pathogenic Mycobacterium tuberculosis is adenylylated at the Tyrosine-51 residue; we also show that GlnD is the adenylyl transferase enzyme responsible. Further analysis shows that in contrast to E. coli, GlnK (PII) adenylylation in M. tuberculosis does not regulate GS adenylylation, nor does it mediate the transcriptomic response to nitrogen limitation.

XperimentR: painless annotation of a biological experiment for the laboratory scientist.

BACKGROUND: Today's biological experiments often involve the collaboration of multidisciplinary researchers utilising several high throughput 'omics platforms. There is a requirement for the details of the experiment to be adequately described using standardised ontologies to enable data preservation, the analysis of the data and to facilitate the export of the data to public repositories. However there are a bewildering number of ontologies, controlled vocabularies, and minimum standards available for use to describe experiments. There is a need for user-friendly software tools to aid laboratory scientists in capturing the experimental information. RESULTS: A web application called XperimentR has been developed for use by laboratory scientists, consisting of a browser-based interface and server-side components which provide an intuitive platform for capturing and sharing experimental metadata. Information recorded includes details about the biological samples, procedures, protocols, and experimental technologies, all of which can be easily annotated using the appropriate ontologies. Files and raw data can be imported and associated with the biological samples via the interface, from either users' computers, or commonly used open-source data repositories. Experiments can be shared with other users, and experiments can be exported in the standard ISA-Tab format for deposition in public databases. XperimentR is freely available and can be installed natively or by using a provided pre-configured Virtual Machine. A guest system is also available for trial purposes. CONCLUSION: We present a web based software application to aid the laboratory scientist to capture, describe and share details about their experiments.

Structure and evolution of barley powdery mildew effector candidates.

BACKGROUND: Protein effectors of pathogenicity are instrumental in modulating host immunity and disease resistance. The powdery mildew pathogen of grasses Blumeria graminis causes one of the most important diseases of cereal crops. B. graminis is an obligate biotrophic pathogen and as such has an absolute requirement to suppress or avoid host immunity if it is to survive and cause disease. RESULTS: Here we characterise a superfamily predicted to be the full complement of Candidates for Secreted Effector Proteins (CSEPs) in the fungal barley powdery mildew parasite B. graminis f.sp. hordei. The 491 genes encoding these proteins constitute over 7% of this pathogen's annotated genes and most were grouped into 72 families of up to 59 members. They were predominantly expressed in the intracellular feeding structures called haustoria, and proteins specifically associated with the haustoria were identified by large-scale mass spectrometry-based proteomics. There are two major types of effector families: one comprises shorter proteins (100-150 amino acids), with a high relative expression level in the haustoria and evidence of extensive diversifying selection between paralogs; the second type consists of longer proteins (300-400 amino acids), with lower levels of differential expression and evidence of purifying selection between paralogs. An analysis of the predicted protein structures underscores their overall similarity to known fungal effectors, but also highlights unexpected structural affinities to ribonucleases throughout the entire effector super-family. Candidate effector genes belonging to the same family are loosely clustered in the genome and are associated with repetitive DNA derived from retro-transposons. CONCLUSIONS: We employed the full complement of genomic, transcriptomic and proteomic analyses as well as structural prediction methods to identify and characterize the members of the CSEPs superfamily in B. graminis f.sp. hordei. Based on relative intron position and the distribution of CSEPs with a ribonuclease-like domain in the phylogenetic tree we hypothesize that the associated genes originated from an ancestral gene, encoding a secreted ribonuclease, duplicated successively by repetitive DNA-driven processes and diversified during the evolution of the grass and cereal powdery mildew lineage.

CombFunc: predicting protein function using heterogeneous data sources.

Only a small fraction of known proteins have been functionally characterized, making protein function prediction essential to propose annotations for uncharacterized proteins. In recent years many function prediction methods have been developed using various sources of biological data from protein sequence and structure to gene expression data. Here we present the CombFunc web server, which makes Gene Ontology (GO)-based protein function predictions. CombFunc incorporates ConFunc, our existing function prediction method, with other approaches for function prediction that use protein sequence, gene expression and protein-protein interaction data. In benchmarking on a set of 1686 proteins CombFunc obtains precision and recall of 0.71 and 0.64 respectively for gene ontology molecular function terms. For biological process GO terms precision of 0.74 and recall of 0.41 is obtained. CombFunc is available at http://www.sbg.bio.ic.ac.uk/combfunc.

Interactions between inflammatory signals and the progesterone receptor in regulating gene expression in pregnant human uterine myocytes.

The absence of a fall in circulating progesterone levels has led to the concept that human labour is associated with 'functional progesterone withdrawal' caused through changes in the expression or function of progesterone receptor (PR). At the time of labour, the human uterus is heavily infiltrated with inflammatory cells, which release cytokines to create a 'myometrial inflammation' via NF-κB activation. The negative interaction between NF-κB and PR, may represent a mechanism to account for 'functional progesterone withdrawal' at term. Conversely, PR may act to inhibit NF-κB function and so play a role in inhibition of myometrial inflammation during pregnancy. To model this inter-relationship, we have used small interfering (si) RNA-mediated knock-down of PR in human pregnant myocytes and whole genome microarray analysis to identify genes regulated through PR. We then activated myometrial inflammation using IL-1ß stimulation to determine the role of PR in myometrial inflammation regulation. Through PR-knock-down, we found that PR regulates gene networks involved in myometrial quiescence and extracellular matrix integrity. Activation of myometrial inflammation was found to antagonize PR-induced gene expression, of genes normally upregulated via PR. We found that PR does not play a role in repression of pro-inflammatory gene networks induced by IL-1ß and that only MMP10 was significantly regulated in opposite directions by IL-1ß and PR. We conclude that progesterone acting through PR does not generally inhibit myometrial inflammation. Activation of myometrial inflammation does cause 'functional progesterone withdrawal' but only in the context of genes normally upregulated via PR.

Gene expression profiling in male genital lichen sclerosus.

Male genital lichen sclerosus (MGLSc) has a bimodal distribution in boys and men. It is associated with squamous cell carcinoma (SCC). The pathogenesis of MGLSc is unknown. HPV and autoimmune mechanisms have been mooted. Anti extracellular matrix protein (ECM)1 antibodies have been identified in women with GLSc. The gene expression pattern of LSc is unknown. Using DNA microarrays we studied differences in gene expression in healthy and diseased prepuces obtained at circumcision in adult males with MGLSc (n = 4), paediatric LSc (n = 2) and normal healthy paediatric foreskin (n = 4). In adult samples 51 genes with significantly increased expression and 87 genes with significantly reduced expression were identified; paediatric samples revealed 190 genes with significantly increased expression and 148 genes with significantly reduced expression. Concordance of expression profiles between adult and paediatric samples indicates the same disease process. Functional analysis revealed increased expression in the adult and child MGSLc samples in the immune response/cellular defence gene ontology (GO) category and reduced expression in other categories including genes related to squamous cancer. No specific HPV, autoimmune or squamous carcinogenesis-associated gene expression patterns were found. ECM1 and CABLES1 expression were significantly reduced in paediatric and adult samples respectively.

Genome expansion and gene loss in powdery mildew fungi reveal tradeoffs in extreme parasitism.

Powdery mildews are phytopathogens whose growth and reproduction are entirely dependent on living plant cells. The molecular basis of this life-style, obligate biotrophy, remains unknown. We present the genome analysis of barley powdery mildew, Blumeria graminis f.sp. hordei (Blumeria), as well as a comparison with the analysis of two powdery mildews pathogenic on dicotyledonous plants. These genomes display massive retrotransposon proliferation, genome-size expansion, and gene losses. The missing genes encode enzymes of primary and secondary metabolism, carbohydrate-active enzymes, and transporters, probably reflecting their redundancy in an exclusively biotrophic life-style. Among the 248 candidate effectors of pathogenesis identified in the Blumeria genome, very few (less than 10) define a core set conserved in all three mildews, suggesting that most effectors represent species-specific adaptations.

Identification of genes selectively disallowed in the pancreatic islet.

We have previously identified two genes, encoding lactate dehydrogenase (Ldha) and the monocarboxylate carrier, MCT1 (Slc16a1) whose expression is remarkably low in pancreatic ß-cells and islets. We sought here to determine whether these may be part of a larger family of genes selectively repressed ("disallowed") in the pancreatic islet. Using new and publicly available microarray data, we undertook a bioinformatic analysis of gene expression in islets and a range of other murine tissues. We compared data sets from three sources of mouse pancreatic islets with a total of 30 datasets from nine tissues, to identify genes with at least five-fold down-regulation in islets. 39 genes were revealed as being specifically repressed in islets. These included Ldha and Slc16a1 as expected but also genes involved in several other metabolic pathways which could affect glucose stimulated insulin secretion. Of these, adenylate kinase 3 (AK3) is a mitochondrial enzyme which acts on GTP, and ornithine aminotransferase (OAT) lies on the pathway converting glutamate to ornithine. The removal of an enzyme which could dissipate mitochondrial GTP levels in beta cells provides support for the theory that mitochondrial GTP may be an important for regulating insulin secretion, whilst blocking an alternative metabolic fate for glutamate is consistent with a signalling role for glutamate. The identification of these genes should inform efforts to generate fully functional ß-cells from stem cell sources, and may provide new targets in type 2 diabetes.

The estrogen receptor-alpha-induced microRNA signature regulates itself and its transcriptional response.

Following estrogenic activation, the estrogen receptor-alpha (ERalpha) directly regulates the transcription of target genes via DNA binding. MicroRNAs (miRNAs) modulated by ERalpha have the potential to fine tune these regulatory systems and also provide an alternate mechanism that could impact on estrogen-dependent developmental and pathological systems. Through a microarray approach, we identify the subset of microRNAs (miRNAs) modulated by ERalpha, which include upregulation of miRNAs derived from the processing of the paralogous primary transcripts (pri-) mir-17-92 and mir-106a-363. Characterization of the mir-17-92 locus confirms that the ERalpha target protein c-MYC binds its promoter in an estrogen-dependent manner. We observe that levels of pri-mir-17-92 increase earlier than the mature miRNAs derived from it, implicating precursor cleavage modulation after transcription. Pri-mir-17-92 is immediately cleaved by DROSHA to pre-miR-18a, indicating that its regulation occurs during the formation of the mature molecule from the precursor. The clinical implications of this novel regulatory system were confirmed by demonstrating that pre-miR-18a was significantly upregulated in ERalpha-positive compared to ERalpha-negative breast cancers. Mechanistically, miRNAs derived from these paralogous pri-miRNAs (miR-18a, miR-19b, and miR-20b) target and downregulate ERalpha, while a subset of pri-miRNA-derived miRNAs inhibit protein translation of the ERalpha transcriptional p160 coactivator, AIB1. Therefore, different subsets of miRNAs identified act as part of a negative autoregulatory feedback loop. We propose that ERalpha, c-MYC, and miRNA transcriptional programs invoke a sophisticated network of interactions able to provide the wide range of coordinated cellular responses to estrogen.

MiMiR--an integrated platform for microarray data sharing, mining and analysis.

BACKGROUND: Despite considerable efforts within the microarray community for standardising data format, content and description, microarray technologies present major challenges in managing, sharing, analysing and re-using the large amount of data generated locally or internationally. Additionally, it is recognised that inconsistent and low quality experimental annotation in public data repositories significantly compromises the re-use of microarray data for meta-analysis. MiMiR, the Microarray data Mining Resource was designed to tackle some of these limitations and challenges. Here we present new software components and enhancements to the original infrastructure that increase accessibility, utility and opportunities for large scale mining of experimental and clinical data. RESULTS: A user friendly Online Annotation Tool allows researchers to submit detailed experimental information via the web at the time of data generation rather than at the time of publication. This ensures the easy access and high accuracy of meta-data collected. Experiments are programmatically built in the MiMiR database from the submitted information and details are systematically curated and further annotated by a team of trained annotators using a new Curation and Annotation Tool. Clinical information can be annotated and coded with a clinical Data Mapping Tool within an appropriate ethical framework. Users can visualise experimental annotation, assess data quality, download and share data via a web-based experiment browser called MiMiR Online. All requests to access data in MiMiR are routed through a sophisticated middleware security layer thereby allowing secure data access and sharing amongst MiMiR registered users prior to publication. Data in MiMiR can be mined and analysed using the integrated EMAAS open source analysis web portal or via export of data and meta-data into Rosetta Resolver data analysis package. CONCLUSION: The new MiMiR suite of software enables systematic and effective capture of extensive experimental and clinical information with the highest MIAME score, and secure data sharing prior to publication. MiMiR currently contains more than 150 experiments corresponding to over 3000 hybridisations and supports the Microarray Centre's large microarray user community and two international consortia. The MiMiR flexible and scalable hardware and software architecture enables secure warehousing of thousands of datasets, including clinical studies, from microarray and potentially other -omics technologies.