Detecting sequence dependent transcriptional pauses from RNA and protein number time series

Background: Evidence suggests that in prokaryotes sequence-dependent transcriptional pauses a?ect the dynamics of transcription and translation, as well as of small genetic circuits. So far, a few pause-prone sequences have been identi?ed from in vitro measurements of transcription elongation kinetics.

Results: Using a stochastic model of gene expression at the nucleotide and codon levels with realistic parameter values, we investigate three di?erent but related questions and present statistical methods for their analysis. First, we show that information from in vivo RNA and protein temporal numbers is su?cient to discriminate between models with and without a pause site in their coding sequence. Second, we demonstrate that it is possible to separate a large variety of models from each other with pauses of various durations and locations in the template by means of a hierarchical clustering and a random forest classi?er. Third, we introduce an approximate likelihood function that allows to estimate the location of a pause site.

Conclusions: This method can aid in detecting unknown pause-prone sequences from temporal measurements of RNA and protein numbers at a genome-wide scale and thus elucidate possible roles that these sequences play in the dynamics of genetic networks and phenotype.

Localization of erythropoietin gene expression in proximal renal tubular cells detected by digoxigenin-labelled oligonucleotide probes

Erythropoietin (EPO) is the main humoral stimulus of erythropoiesis. In adult mammals, the kidney releases EPO in response to hypoxic stress. Conflicting data have suggested either renal tubular or peritubular cell origins of EPO synthesis in vivo. In situ hybridization studies were performed to define further the kidney cell type(s) capable of increasing EPO gene expression during hypoxic stimulation. EPO gene expression was stimulated in mice exposed to acute hypobaric hypoxia. Kidneys from hypoxic and control normoxic mice were obtained. Six digoxigenin-labelled oligonucleotide probes complementary to EPO exon sequences were utilized for in situ hybridization for EPO messenger RNA. Positive hybridization signals were identified in some proximal tubular cells, confined to the inner third of the renal cortex of hypoxic mouse kidney.

RNA interference targeting cathepsin B of the carcinogenic liver fluke, Opisthorchis viverrini

Functional genomics have not been reported for Opisthorchis viverrini or the related fish-borne fluke, Clonorchis sinensis. Here we describe the introduction by square wave electroporation of Cy3-labeled small RNA into adult O. viverrini worms. Adult flukes were subjected to square wave electroporation employing a single pulse for 20 ms of 125V in the presence of 50 µg/ml of Cy3-siRNA. The parasites tolerated this manipulation and, at 24 and 48 h after electroporation, fluorescence from the Cy3-siRNA was evident throughout the parenchyma of the worms, with strong fluorescence evident in the guts and reproductive organs of the adult worms. Second, other worms were treated using the same electroporation settings with double stranded RNA targeting an endogenous papain-like cysteine protease, cathepsin B. This manipulation resulted in a significant reduction in specific mRNA levels encoding cathepsin B, and a significant reduction in cathepsin B activity against the diagnostic peptide, Z-Arg-Arg-AMC. This appears to be the first report of introduction of reporter genes into O. viverrini and the first report of experimental RNA interference (RNAi) in this fluke. The findings indicated the presence of an intact RNAi pathway in these parasites which, in turn, provides an opportunity to probe gene functions in this neglected tropical disease pathogen.

Isolation and characterization of a novel simazine-degrading bacterium from agricultural soil of central Chile, Pseudomonas sp MHP41

s-Triazine herbicides are used extensively in South America in agriculture and forestry. In this study, a bacterium designated as strain MHP41, capable of degrading simazine and atrazine, was isolated from agricultural soil in the Quillota valley, central Chile. Strain MHP41 is able to grow in minimal medium, using simazine as the sole nitrogen source. In this medium, the bacterium exhibited a growth rate of mu = 0.10 h(-1), yielding a high biomass of 4.2 x 10(8) CFU mL(-1). Resting cells of strain MHP41 degrade more than 80% of simazine within 60 min. The atzA, atzB, atzC, atzD, atzE and atzF genes encoding the enzymes of the simazine upper and lower pathways were detected in strain MHP41. The motile Gram-negative bacterium was identified as a Pseudomonas sp., based on the Biolog microplate system and comparative sequence analyses of the 16S rRNA gene. Amplified ribosomal DNA restriction analysis allowed the differentiation of strain MHP41 from Pseudomonas sp. ADP. The comparative 16S rRNA gene sequence analyses suggested that strain MHP41 is closely related to Pseudomonas nitroreducens and Pseudomonas multiresinovorans. This is the first s-triazine-degrading bacterium isolated in South America. Strain MHP41 is a potential biocatalyst for the remediation of s-triazine-contaminated environments.

Connectivity Mapping for Candidate Therapeutics Identification Using Next Generation Sequencing RNA-Seq Data

The advent of next generation sequencing technologies (NGS) has expanded the area of genomic research, offering high coverage and increased sensitivity over older microarray platforms. Although the current cost of next generation sequencing is still exceeding that of microarray approaches, the rapid advances in NGS will likely make it the platform of choice for future research in differential gene expression. Connectivity mapping is a procedure for examining the connections among diseases, genes and drugs by differential gene expression initially based on microarray technology, with which a large collection of compound-induced reference gene expression profiles have been accumulated. In this work, we aim to test the feasibility of incorporating NGS RNA-Seq data into the current connectivity mapping framework by utilizing the microarray based reference profiles and the construction of a differentially expressed gene signature from a NGS dataset. This would allow for the establishment of connections between the NGS gene signature and those microarray reference profiles, alleviating the associated incurring cost of re-creating drug profiles with NGS technology. We examined the connectivity mapping approach on a publicly available NGS dataset with androgen stimulation of LNCaP cells in order to extract candidate compounds that could inhibit the proliferative phenotype of LNCaP cells and to elucidate their potential in a laboratory setting. In addition, we also analyzed an independent microarray dataset of similar experimental settings. We found a high level of concordance between the top compounds identified using the gene signatures from the two datasets. The nicotine derivative cotinine was returned as the top candidate among the overlapping compounds with potential to suppress this proliferative phenotype. Subsequent lab experiments validated this connectivity mapping hit, showing that cotinine inhibits cell proliferation in an androgen dependent manner. Thus the results in this study suggest a promising prospect of integrating NGS data with connectivity mapping. © 2013 McArt et al.

Cytosolic phospholipase A2 gene expression in rat mesangial cells is regulated post-transcriptionally

Cytosolic phospholipase A2 (cPLA2) is thought to be the rate-limiting enzyme in the arachidonic acid/eicosanoid cascade. The ability of various agonists to increase steady-state cPLA2 mRNA levels has previously been reported. The current study delineates the contributions of transcriptional and post-transcriptional processes to the regulation of cPLA2 gene expression in response to a variety of agonists in cultured rat glomerular mesangial cells. Epidermal growth factor, platelet-derived growth factor, serum and phorbol myristate acetate all increase the half-life of cPLA2 mRNA transcripts, indicating a role for post-transcriptional modulation of gene expression. The presence of three ATTTA motifs in the 3' untranslated region (3'UTR) of the rat cPLA2 cDNA is ascertained. Heterologous expression of chimeric constructs with different 3'UTRs ligated into the 3' end of the luciferase coding region reveals that the presence of the cPLA2 3'UTR results in reduced luciferase activity compared with constructs without the cPLA2 3'UTR. Furthermore, the luciferase activity in the constructs with the cPLA2 3'UTR is increased in response to the same agonists which stabilize endogenous cPLA2 mRNA. A negligible effect of these agonists on transcriptional control of cPLA2 is evident using promoter-reporter constructs expressed in transient and stable transfectants. Taken together, these results indicate predominant post-transcriptional regulation of cPLA2 mRNA levels.

RNAi dynamics in Juvenile Fasciola spp. Liver flukes reveals the persistence of gene silencing in vitro

Background: Fasciola spp. liver fluke cause pernicious disease in humans and animals. Whilst current control is unsustainable due to anthelmintic resistance, gene silencing (RNA interference, RNAi) has the potential to contribute to functional validation of new therapeutic targets. The susceptibility of juvenile Fasciola hepatica to double stranded (ds)RNA-induced RNAi has been reported. To exploit this we probe RNAi dynamics, penetrance and persistence with the aim of building a robust platform for reverse genetics in liver fluke. We describe development of standardised RNAi protocols for a commercially-available liver fluke strain (the US Pacific North West Wild Strain), validated via robust transcriptional silencing of seven virulence genes, with in-depth experimental optimisation of three: cathepsin L (FheCatL) and B (FheCatB) cysteine proteases, and a σ-class glutathione transferase (FheσGST).

Methodology/Principal Findings: Robust transcriptional silencing of targets in both F. hepatica and Fasciola gigantica juveniles is achievable following exposure to long (200–320 nt) dsRNAs or 27 nt short interfering (si)RNAs. Although juveniles are highly RNAi-susceptible, they display slower transcript and protein knockdown dynamics than those reported previously. Knockdown was detectable following as little as 4h exposure to trigger (target-dependent) and in all cases silencing persisted for ≥25 days following long dsRNA exposure. Combinatorial silencing of three targets by mixing multiple long dsRNAs was similarly efficient. Despite profound transcriptional suppression, we found a significant time-lag before the occurrence of protein suppression; FheσGST and FheCatL protein suppression were only detectable after 9 and 21 days, respectively.

Conclusions/Significance: In spite of marked variation in knockdown dynamics, we find that a transient exposure to long dsRNA or siRNA triggers robust RNAi penetrance and persistence in liver fluke NEJs supporting the development of multiple-throughput phenotypic screens for control target validation. RNAi persistence in fluke encourages in vivo studies on gene function using worms exposed to RNAi-triggers prior to infection.

S100A2 is a BRCA1/p63 coregulated tumour suppressor gene with roles in the regulation of mutant p53 stability

Here, we show for the first time that the familial breast/ovarian cancer susceptibility gene, BRCA1, along with interacting ΔNp63 proteins, transcriptionally upregulate the putative tumour suppressor protein, S100A2. Both BRCA1 and ΔNp63 proteins are required for S100A2 expression. BRCA1 requires ΔNp63 proteins for recruitment to the S100A2 proximal promoter region, while exogenous expression of individual ΔNp63 proteins cannot activate S100A2 transcription in the absence of a functional BRCA1. Consequently, mutation of the ΔNp63/p53 response element within the S100A2 promoter completely abrogates the ability of BRCA1 to upregulate S100A2. S100A2 shows growth control features in a range of cell models. Transient or stable exogenous S100A2 expression inhibits the growth of BRCA1 mutant and basal-like breast cancer cell lines, while short interfering RNA (siRNA) knockdown of S100A2 in non-tumorigenic cells results in enhanced proliferation. S100A2 modulates binding of mutant p53 to HSP90, which is required for efficient folding of mutant p53 proteins, by competing for binding to HSP70/HSP90 organising protein (HOP). HOP is a cochaperone that is required for the efficient transfer of proteins from HSP70 to HSP90. Loss of S100A2 leads to an HSP90-dependent stabilisation of mutant p53 with a concomitant loss of p63. Accordingly, S100A2-deficient cells are more sensitive to the HSP-90 inhibitor, 17-N-allylamino-17-demethoxygeldanamycin, potentially representing a novel therapeutic strategy for S100A2- and BRCA1-deficient cancers. Taken together, these data demonstrate the importance of S100A2 downstream of the BRCA1/ΔNp63 signalling axis in modulating transcriptional responses and enforcing growth control mechanisms through destabilisation of mutant p53.

Gene expression in normal Urothelium depends on location within the bladder:A possible link to bladder carcinogenesis

Objectives: Clinical studies have shown that more than 70% of primary bladder tumours arise in the area around the ureteric orifice. In this study a genomic approach was taken to explore the molecular mechanisms that may influence this phenomenon.

Methods: RNA was isolated from each individual normal ureteric orifice and the dome biopsy from 33 male patients. Equal amounts of the pooled ureteric orifice and dome mRNAs were labelled with Cy3 and Cy5, respectively before hybridising to the gene chip (UniGEM 2.0, Incyte Genomics Inc., Wilmington, Delaware, USA). Results: Significant changes (more than a twofold difference) in gene expression were observed in 3.1% (312) of the 10,176 gene array: 211 genes upregulated and 101 downregulated. Analysis of Cdc25B, TK1, PKM, and PDGFra with RT-PCR supported the reliability of the microarray result. Seladin-1 was the most upregulated gene in the ureteric orifice: 8.3-fold on the microarray and 11.4-fold by real time PCR.

Conclusions: Overall, this study suggests significant altered gene expression between these two anatomically distinct areas of the normal human bladder. Of particular note is Seladin-1, whose significance in cancer is yet to be clarified. Further studies of the genes discovered by this work will help clarify which of these differences influence primary bladder carcinogenesis. (c) 2006 European Association of Urology. Published by Elsevier B.V. All rights reserved.

Three motAB stator gene products in Bdellovibrio bacteriovorus contribute to motility of a single flagellum during predatory and prey-independent growth

The predatory bacterium Bdellovibrio bacteriovorus uses flagellar motility to locate regions rich in Gram-negative prey bacteria, colliding and attaching to prey and then ceasing flagellar motility. Prey are then invaded to form a "bdelloplast" in a type IV pilus-dependent process, and prey contents are digested, allowing Bdellovibrio growth and septation. After septation, Bdellovibrio flagellar motility resumes inside the prey bdelloplast prior to its lysis and escape of Bdellovibrio progeny. Bdellovibrio can also grow slowly outside prey as long flagellate host-independent (HI) cells, cultured on peptone-rich media. The B. bacteriovorus HD100 genome encodes three pairs of MotAB flagellar motor proteins, each of which could potentially form an inner membrane ion channel, interact with the FliG flagellar rotor ring, and produce flagellar rotation. In 2004, Flannagan and coworkers (R. S. Flannagan, M. A. Valvano, and S. F. Koval, Microbiology 150:649-656, 2004) used antisense RNA and green fluorescent protein (GFP) expression to downregulate a single Bdellovibrio motA gene and reported slowed release from the bdelloplast and altered motility of the progeny. Here we inactivated each pair of motAB genes and found that each pair contributes to motility, both predatorily, inside the bdelloplast and during HI growth; however, each pair was dispensable, and deletion of no pair abolished motility totally. Driving-ion studies with phenamil, carbonyl cyanide m-chlorophenylhydrazone (CCCP), and different pH and sodium conditions indicated that all Mot pairs are proton driven, although the sequence similarities of each Mot pair suggests that some may originate from halophilic species. Thus, Bdellovibrio is a "dedicated motorist," retaining and expressing three pairs of mot genes.

RNA interference in adult Ascaris suum - an opportunity for the development of a functional genomics platform that supports organism-, tissue- and cell-based biology in a nematode parasite

The sustainable control of animal parasitic nematodes requires the development of efficient functional genomics platforms to facilitate target validation and enhance anthelmintic discovery. Unfortunately, the utility of RNA interference (RNAi) for the validation of novel drug targets in nematode parasites remains problematic. Ascaris suum is an important veterinary parasite and a zoonotic pathogen. Here we show that adult A. suum is RNAi competent, and highlight the induction, spread and consistency of RNAi across multiple tissue types. This platform provides a new opportunity to undertake whole organism-, tissue- and cell-level gene function studies to enhance target validation processes for nematode parasites of veterinary/medical significance.

Inactivating mutations in an SH2 domain-encoding gene in X-linked lymphoproliferative syndrome

X-linked lymphoproliferative syndrome (XLP) is an inherited immunodeficiency characterized by increased susceptibility to Epstein-Barr virus (EBV). In affected males, primary EBV infection leads to the uncontrolled proliferation of virus-containing B cells and reactive cytotoxic T cells, often culminating in the development of high-grade lymphoma. The XLP gene has been mapped to chromosome band Xq25 through linkage analysis and the discovery of patients harboring large constitutional genomic deletions. We describe here the presence of small deletions and intragenic mutations that specifically disrupt a gene named DSHP in 6 of 10 unrelated patients with XLP. This gene encodes a predicted protein of 128 amino acids composing a single SH2 domain with extensive homology to the SH2 domain of SHIP, an inositol polyphosphate 5-phosphatase that functions as a negative regulator of lymphocyte activation. DSHP is expressed in transformed T cell lines and is induced following in vitro activation of peripheral blood T lymphocytes. Expression of DSHP is restricted in vivo to lymphoid tissues, and RNA in situ hybridization demonstrates DSHP expression in activated T and B cell regions of reactive lymph nodes and in both T and B cell neoplasms. These observations confirm the identity of DSHP as the gene responsible for XLP, and suggest a role in the regulation of lymphocyte activation and proliferation. Induction of DSHP may sustain the immune response by interfering with SHIP-mediated inhibition of lymphocyte activation, while its inactivation in XLP patients results in a selective immunodeficiency to EBV.

Molecular Medicine and Molecular Pathology for Haematologists:I. Gene Speak made Easy: An Overview of Genes and Gene Expression

Molecular Medicine and Molecular Pathology are integral parts of Haematology as we enter the new millennium. Their origins can be linked to fundamental developments in the basic sciences, particularly genetics, chemistry and biochemistry. The structure of DNA and the genetic code that it encrypts are the critical starting points to our understanding of these new disciplines. The genetic alphabet is a simple one, consisting of just 4 letters, buts its influence is crucial to human development and differentiation. The concept of a gene is not a new one but the Human Genome Project (a joint world-wide effort to characterise our entire genetic make-up) is providing an invaluable understanding of how genes function in normal cellular processes and pinpointing how disruption of these processes can lead to disease. Transcription and translation are the key events by which our genotype is converted to our phenotype (via a messenger RNA intermediate), producing the myriad proteins and enzymes which populate the cellular factory of our body. Unlike the bacterial or prokaryotic genome, the human genome contains a large amount of non coding DNA (less than 1% of our genome codes for proteins), and our genes are interrupted, with the coding regions or exons separated by non coding introns. Precise removal of the intronic material after transcription (though a process called splicing) is critical for efficient translation to occur. Incorrect splicing can lead to the generation of mutant proteins, which can have a dilaterious effect on the phenotype of the individual. Thus the 100,000-200,000 genes which are present in each cell in our body have a defined control mechanism permitting efficient and appropriate expression of proteins and enzymes and yet a single base change in just one of those genes can lead to diseases such as haemophilia or fanconis anaemia.

Characterization of a reference material for BCR-ABL (M-BCR) mRNA quantitation by real-time amplification assays:towards new standards for gene expression measurements

Monitoring of BCR-ABL transcripts has become established practice in the management of chronic myeloid leukemia. However, nucleic acid amplification techniques are prone to variations which limit the reliability of real-time quantitative PCR (RQ-PCR) for clinical decision making, highlighting the need for standardization of assays and reporting of minimal residual disease (MRD) data. We evaluated a lyophilized preparation of a leukemic cell line (K562) as a potential quality control reagent. This was found to be relatively stable, yielding comparable respective levels of ABL, GUS and BCR-ABL transcripts as determined by RQ-PCR before and after accelerated degradation experiments as well as following 5 years storage at -20 degrees C. Vials of freeze-dried cells were sent at ambient temperature to 22 laboratories on four continents, with RQ-PCR analyses detecting BCR-ABL transcripts at levels comparable to those observed in primary patient samples. Our results suggest that freeze-dried cells can be used as quality control reagents with a range of analytical instrumentations and could enable the development of urgently needed international standards simulating clinically relevant levels of MRD.

Meta-analysis of prostate cancer gene expression data identifies a novel discriminatory signature enriched for glycosylating enzymes

BACKGROUND: Tumorigenesis is characterised by changes in transcriptional control. Extensive transcript expression data have been acquired over the last decade and used to classify prostate cancers. Prostate cancer is, however, a heterogeneous multifocal cancer and this poses challenges in identifying robust transcript biomarkers.

METHODS: In this study, we have undertaken a meta-analysis of publicly available transcriptomic data spanning datasets and technologies from the last decade and encompassing laser capture microdissected and macrodissected sample sets.

RESULTS: We identified a 33 gene signature that can discriminate between benign tissue controls and localised prostate cancers irrespective of detection platform or dissection status. These genes were significantly overexpressed in localised prostate cancer versus benign tissue in at least three datasets within the Oncomine Compendium of Expression Array Data. In addition, they were also overexpressed in a recent exon-array dataset as well a prostate cancer RNA-seq dataset generated as part of the The Cancer Genomics Atlas (TCGA) initiative. Biologically, glycosylation was the single enriched process associated with this 33 gene signature, encompassing four glycosylating enzymes. We went on to evaluate the performance of this signature against three individual markers of prostate cancer, v-ets avian erythroblastosis virus E26 oncogene homolog (ERG) expression, prostate specific antigen (PSA) expression and androgen receptor (AR) expression in an additional independent dataset. Our signature had greater discriminatory power than these markers both for localised cancer and metastatic disease relative to benign tissue, or in the case of metastasis, also localised prostate cancer.

CONCLUSION: In conclusion, robust transcript biomarkers are present within datasets assembled over many years and cohorts and our study provides both examples and a strategy for refining and comparing datasets to obtain additional markers as more data are generated.