Reaction dynamics under confinement: an exact path integral treatment of a two-stage model of stochastic gene expression

Gene expression in living systems is inherently stochastic, and tends to produce varying numbers of proteins over repeated cycles of transcription and translation. In this paper, an expression is derived for the steady-state protein number distribution starting from a two-stage kinetic model of the gene expression process involving p proteins and r mRNAs. The derivation is based on an exact path integral evaluation of the joint distribution, P(p, r, t), of p and r at time t, which can be expressed in terms of the coupled Langevin equations for p and r that represent the two-stage model in continuum form. The steady-state distribution of p alone, P(p), is obtained from P(p, r, t) (a bivariate Gaussian) by integrating out the r degrees of freedom and taking the limit t -> infinity. P(p) is found to be proportional to the product of a Gaussian and a complementary error function. It provides a generally satisfactory fit to simulation data on the same two-stage process when the translational efficiency (a measure of intrinsic noise levels in the system) is relatively low; it is less successful as a model of the data when the translational efficiency (and noise levels) are high.

The influence of reduced oxygen availability on gene expression in laboratory (H37Rv) and clinical strains (S7 and S10) of Mycobacterium tuberculosis

Mycobacterium tuberculosis has the ability to persist within the host in a dormant stage. One important condition believed to contribute to dormancy is reduced access to oxygen known as hypoxia. However, the response of M. tuberculosis to such hypoxia condition is not fully characterized. Virtually all dormant models against tuberculosis tested in animals used laboratory strain H37Rv or Erdman strain. But major outbreaks of tuberculosis (TB) occur with the strains that have widely different genotypes and phenotypes compared to H37Rv. In this study, we used a custom oligonucleotide microarray to determine the overall transcriptional response of laboratory strain (H37Rv) and most prevalent clinical strains (S7 and S10) of M. tuberculosis from South India to hypoxia. Analysis of microarray results revealed that a total of 1161 genes were differentially regulated (>= 1.5 fold change) in H37Rv, among them 659 genes upregulated and 502 genes down regulated. Microarray data of clinical isolates showed that a total of 790 genes were differentially regulated in S7 among which 453 genes were upregulated and 337 down regulated. Interestingly, numerous genes were also differentially regulated in S10 (total 2805 genes) of which 1463 genes upregulated and 1342 genes down regulated during reduced oxygen condition (Wayne's model). One hundred and thirty-four genes were found common and upregulated among all three strains (H37Rv, S7, and S10) and can be targeted for drug/vaccine development against TB. (C) 2015 Published by Elsevier B.V.

ICGE: an R package for detecting relevant clusters and atypical units in gene expression

Background: Gene expression technologies have opened up new ways to diagnose and treat cancer and other diseases. Clustering algorithms are a useful approach with which to analyze genome expression data. They attempt to partition the genes into groups exhibiting similar patterns of variation in expression level. An important problem associated with gene classification is to discern whether the clustering process can find a relevant partition as well as the identification of new genes classes. There are two key aspects to classification: the estimation of the number of clusters, and the decision as to whether a new unit (gene, tumor sample ... ) belongs to one of these previously identified clusters or to a new group. Results: ICGE is a user-friendly R package which provides many functions related to this problem: identify the number of clusters using mixed variables, usually found by applied biomedical researchers; detect whether the data have a cluster structure; identify whether a new unit belongs to one of the pre-identified clusters or to a novel group, and classify new units into the corresponding cluster. The functions in the ICGE package are accompanied by help files and easy examples to facilitate its use. Conclusions: We demonstrate the utility of ICGE by analyzing simulated and real data sets. The results show that ICGE could be very useful to a broad research community.

Gene expression profiles in liver of zebrafish treated with microcystin-LR

Microcystin-LR (MC-LR) is the most frequently studied cyclic heptatoxin produced by cyanobacteria, which has tremendous negative impacts on fish, while its molecular mechanism behind remained unclear at present. Here, Affymetrix Zebrafish GeneChip was used to identify alterations in gene expression of zebrafish (Danio rerio) after MC-LR exposure. Among the 14,900 transcripts in the microarray, 273 genes were differentially expressed, in which 243 genes were elevated and 30 were decreased. According to GOstat analysis, MC-LR mainly influenced the cell cycle and mitogen-activated protein kinases (MAPK) signaling pathways. In addition, many immune-related genes were also influenced. These data suggest that MC-LR could promote tumorigenesis and cause immunotoxicity in fish. (C) 2008 Elsevier B.V. All rights reserved.

Comparison of Gene Expression Profiles of Fenneropenaeus chinensis Challenged with WSSV and Vibrio

Microarray technique was used to analyze the gene expression profiles of shrimp when they were challenged by WSSV and heat-inactivated Vibrio anguillarum, respectively. At 6 h post challenge (HPC), a total of 806 clones showed differential expression profile in WSSV-challenged samples, but not in Vibrio-challenged samples. The genes coding energy metabolism enzyme and structure protein were the most downregulated elements in 6 h post WSSV-challenged (HPC-WSSV) tissues. However, a total of 155 clones showed differential expression in the Vibrio-challenged samples, but not in WSSV-challenged samples. Serine-type endopeptidase and lysosome-related genes were the most upregulated elements in tissues 6 h post Vibrio challenge (HPC-Vibrio). Totally, 188 clones showed differential expression in both 6 and 12 HPC-WSSV and HPC-Vibrio samples. Most of the differentially expressed genes (185/188) were downregulated in the samples of 12 HPC-WSSV, whereas upregulated in the samples at 6 and 12 HPC-Vibrio and 6 HPC-WSSV. The expression profiles of three differentially expressed genes identified in microarray hybridization were analyzed in hemocytes, lymphoid organ, and hepatopancreas of shrimp challenged by WSSV or Vibrio through real-time PCR. The results further confirmed the microarray hybridization results. The data will provide great help for us in understanding the immune mechanism of shrimp responding to WSSV or Vibrio.

Rubisco small subunit, chlorophyll a/b-binding protein and sucrose : fructan-6-fructosyl transferase gene expression and sugar status in single barley leaf cells in situ. Cell type specificity and induction by light

Chungui Lu, Olga A. Koroleva, John F. Farrar, Joe Gallagher, Chris J. Pollock, and A. Deri Tomos (2002). Rubisco small subunit, chlorophyll a/b-binding protein and sucrose : fructan-6-fructosyl transferase gene expression and sugar status in single barley leaf cells in situ. Cell type specificity and induction by light. Plant Physiology, 130 (3) pp.1335-1348 Sponsorship: BBSRC RAE2008

SplicerAV: a tool for mining microarray expression data for changes in RNA processing.

BACKGROUND: Over the past two decades more than fifty thousand unique clinical and biological samples have been assayed using the Affymetrix HG-U133 and HG-U95 GeneChip microarray platforms. This substantial repository has been used extensively to characterize changes in gene expression between biological samples, but has not been previously mined en masse for changes in mRNA processing. We explored the possibility of using HG-U133 microarray data to identify changes in alternative mRNA processing in several available archival datasets. RESULTS: Data from these and other gene expression microarrays can now be mined for changes in transcript isoform abundance using a program described here, SplicerAV. Using in vivo and in vitro breast cancer microarray datasets, SplicerAV was able to perform both gene and isoform specific expression profiling within the same microarray dataset. Our reanalysis of Affymetrix U133 plus 2.0 data generated by in vitro over-expression of HRAS, E2F3, beta-catenin (CTNNB1), SRC, and MYC identified several hundred oncogene-induced mRNA isoform changes, one of which recognized a previously unknown mechanism of EGFR family activation. Using clinical data, SplicerAV predicted 241 isoform changes between low and high grade breast tumors; with changes enriched among genes coding for guanyl-nucleotide exchange factors, metalloprotease inhibitors, and mRNA processing factors. Isoform changes in 15 genes were associated with aggressive cancer across the three breast cancer datasets. CONCLUSIONS: Using SplicerAV, we identified several hundred previously uncharacterized isoform changes induced by in vitro oncogene over-expression and revealed a previously unknown mechanism of EGFR activation in human mammary epithelial cells. We analyzed Affymetrix GeneChip data from over 400 human breast tumors in three independent studies, making this the largest clinical dataset analyzed for en masse changes in alternative mRNA processing. The capacity to detect RNA isoform changes in archival microarray data using SplicerAV allowed us to carry out the first analysis of isoform specific mRNA changes directly associated with cancer survival.

Altered gene expression and DNA damage in peripheral blood cells from Friedreich's ataxia patients: cellular model of pathology.

The neurodegenerative disease Friedreich's ataxia (FRDA) is the most common autosomal-recessively inherited ataxia and is caused by a GAA triplet repeat expansion in the first intron of the frataxin gene. In this disease, transcription of frataxin, a mitochondrial protein involved in iron homeostasis, is impaired, resulting in a significant reduction in mRNA and protein levels. Global gene expression analysis was performed in peripheral blood samples from FRDA patients as compared to controls, which suggested altered expression patterns pertaining to genotoxic stress. We then confirmed the presence of genotoxic DNA damage by using a gene-specific quantitative PCR assay and discovered an increase in both mitochondrial and nuclear DNA damage in the blood of these patients (p<0.0001, respectively). Additionally, frataxin mRNA levels correlated with age of onset of disease and displayed unique sets of gene alterations involved in immune response, oxidative phosphorylation, and protein synthesis. Many of the key pathways observed by transcription profiling were downregulated, and we believe these data suggest that patients with prolonged frataxin deficiency undergo a systemic survival response to chronic genotoxic stress and consequent DNA damage detectable in blood. In conclusion, our results yield insight into the nature and progression of FRDA, as well as possible therapeutic approaches. Furthermore, the identification of potential biomarkers, including the DNA damage found in peripheral blood, may have predictive value in future clinical trials.

Intergenic and genic sequence lengths have opposite relationships with respect to gene expression.

Eukaryotic genomes are mostly composed of noncoding DNA whose role is still poorly understood. Studies in several organisms have shown correlations between the length of the intergenic and genic sequences of a gene and the expression of its corresponding mRNA transcript. Some studies have found a positive relationship between intergenic sequence length and expression diversity between tissues, and concluded that genes under greater regulatory control require more regulatory information in their intergenic sequences. Other reports found a negative relationship between expression level and gene length and the interpretation was that there is selection pressure for highly expressed genes to remain small. However, a correlation between gene sequence length and expression diversity, opposite to that observed for intergenic sequences, has also been reported, and to date there is no testable explanation for this observation. To shed light on these varied and sometimes conflicting results, we performed a thorough study of the relationships between sequence length and gene expression using cell-type (tissue) specific microarray data in Arabidopsis thaliana. We measured median gene expression across tissues (expression level), expression variability between tissues (expression pattern uniformity), and expression variability between replicates (expression noise). We found that intergenic (upstream and downstream) and genic (coding and noncoding) sequences have generally opposite relationships with respect to expression, whether it is tissue variability, median, or expression noise. To explain these results we propose a model, in which the lengths of the intergenic and genic sequences have opposite effects on the ability of the transcribed region of the gene to be epigenetically regulated for differential expression. These findings could shed light on the role and influence of noncoding sequences on gene expression.

Diagnosis of partial body radiation exposure in mice using peripheral blood gene expression profiles.

In the event of a terrorist-mediated attack in the United States using radiological or improvised nuclear weapons, it is expected that hundreds of thousands of people could be exposed to life-threatening levels of ionizing radiation. We have recently shown that genome-wide expression analysis of the peripheral blood (PB) can generate gene expression profiles that can predict radiation exposure and distinguish the dose level of exposure following total body irradiation (TBI). However, in the event a radiation-mass casualty scenario, many victims will have heterogeneous exposure due to partial shielding and it is unknown whether PB gene expression profiles would be useful in predicting the status of partially irradiated individuals. Here, we identified gene expression profiles in the PB that were characteristic of anterior hemibody-, posterior hemibody- and single limb-irradiation at 0.5 Gy, 2 Gy and 10 Gy in C57Bl6 mice. These PB signatures predicted the radiation status of partially irradiated mice with a high level of accuracy (range 79-100%) compared to non-irradiated mice. Interestingly, PB signatures of partial body irradiation were poorly predictive of radiation status by site of injury (range 16-43%), suggesting that the PB molecular response to partial body irradiation was anatomic site specific. Importantly, PB gene signatures generated from TBI-treated mice failed completely to predict the radiation status of partially irradiated animals or non-irradiated controls. These data demonstrate that partial body irradiation, even to a single limb, generates a characteristic PB signature of radiation injury and thus may necessitate the use of multiple signatures, both partial body and total body, to accurately assess the status of an individual exposed to radiation.

A noisy linear map underlies oscillations in cell size and gene expression in bacteria.

During bacterial growth, a cell approximately doubles in size before division, after which it splits into two daughter cells. This process is subjected to the inherent perturbations of cellular noise and thus requires regulation for cell-size homeostasis. The mechanisms underlying the control and dynamics of cell size remain poorly understood owing to the difficulty in sizing individual bacteria over long periods of time in a high-throughput manner. Here we measure and analyse long-term, single-cell growth and division across different Escherichia coli strains and growth conditions. We show that a subset of cells in a population exhibit transient oscillations in cell size with periods that stretch across several (more than ten) generations. Our analysis reveals that a simple law governing cell-size control-a noisy linear map-explains the origins of these cell-size oscillations across all strains. This noisy linear map implements a negative feedback on cell-size control: a cell with a larger initial size tends to divide earlier, whereas one with a smaller initial size tends to divide later. Combining simulations of cell growth and division with experimental data, we demonstrate that this noisy linear map generates transient oscillations, not just in cell size, but also in constitutive gene expression. Our work provides new insights into the dynamics of bacterial cell-size regulation with implications for the physiological processes involved.

Characterization of autonomous thyroid adenoma: metabolism, gene expression, and pathology.

Fifty-one in vivo characterized autonomous single adenomas have been studied for functional parameters in vitro, for gene and protein expression and for pathology, and have been systematically compared to the corresponding extratumoral quiescent tissue. The adenomas were characterized by a high level of iodide trapping that corresponds to a high level of Na+ /iodide symporter gene expression, a high thyroperoxidase mRNA and protein content, and a low H2O2 generation. This explains the iodide metabolism characteristics demonstrated before, ie, the main cause of the "hot" character of the adenomas is their increased iodide transport. The adenomas spontaneously secreted higher amounts of thyroid hormone than the quiescent tissue and in agreement with previous in vivo data, this secretion could be further enhanced by thyrotropin (TSH). Inositol uptake was also increased but there was no spontaneous increase of the generation of inositol phosphates and this metabolism could be further activated by TSH. These positive responses to TSH are in agreement with the properties of TSH-stimulated thyroid cells in vitro and in vivo. They are compatible with the characteristics of mutated TSH receptors whose constitutive activation accounts for the majority of autonomous thyroid adenomas in Europe. The number of cycling cells, as evaluated by MIB-1 immunolabeling was low but increased in comparison with the corresponding quiescent tissue or normal tissue. The cycling cells are observed mainly at the periphery; there was very little apoptosis. Both findings account for the slow growth of these established adenomas. On the other hand, by thyroperoxidase immunohistochemistry, the whole lesion appeared hyperfunctional, which demonstrates a dissociation of mitogenic and functional stimulations. Thyroglobulin, TSH receptor, and E-cadherin mRNA accumulations were not modified in a consistent way, which confirms the near-constitutive expression of the corresponding genes in normal differentiated tissue. On the contrary, early immediate genes expressions (c-myc, NGF1B, egr 1, genes of the fos and jun families) were decreased. This may be explained by the proliferative heterogeneity of the lesion and the previously described short, biphasic expression of these genes when induced by mitogenic agents. All the characteristics of the autonomous adenomas can therefore be explained by the effect of the known activating mutations of genes coding for proteins of the TSH cyclic adenosine monophosphate (cAMP) cascade, all cells being functionally activated while only those at the periphery multiply. The reason of this heterogeneity is unknown.

sscMap: An extensible Java application for connecting small-molecule drugs using gene-expression signatures

Background
Connectivity mapping is a process to recognize novel pharmacological and toxicological properties in small molecules by comparing their gene expression signatures with others in a database. A simple and robust method for connectivity mapping with increased specificity and sensitivity was recently developed, and its utility demonstrated using experimentally derived gene signatures.

Results
This paper introduces sscMap (statistically significant connections' map), a Java application designed to undertake connectivity mapping tasks using the recently published method. The software is bundled with a default collection of reference gene-expression profiles based on the publicly available dataset from the Broad Institute Connectivity Map 02, which includes data from over 7000 Affymetrix microarrays, for over 1000 small-molecule compounds, and 6100 treatment instances in 5 human cell lines. In addition, the application allows users to add their custom collections of reference profiles and is applicable to a wide range of other 'omics technologies.

Conclusion
The utility of sscMap is two fold. First, it serves to make statistically significant connections between a user-supplied gene signature and the 6100 core reference profiles based on the Broad Institute expanded dataset. Second, it allows users to apply the same improved method to custom-built reference profiles which can be added to the database for future referencing. The software can be freely downloaded from http://purl.oclc.org/NET/sscMap

A statistical framework for the design of microarray experiments and effective detection of differential gene expression

Motivation: Microarray experiments generate a high data volume. However, often due to financial or experimental considerations, e.g. lack of sample, there is little or no replication of the experiments or hybridizations. These factors combined with the intrinsic variability associated with the measurement of gene expression can result in an unsatisfactory detection rate of differential gene expression (DGE). Our motivation was to provide an easy to use measure of the success rate of DGE detection that could find routine use in the design of microarray experiments or in post-experiment assessment.

Novel breast cancer biomarkers identified by integrative proteomic and gene expression mapping

Proteomic and transcriptomic platforms both play important roles in cancer research, with differing strengths and limitations. Here, we describe a proteo-transcriptomic integrative strategy for discovering novel cancer biomarkers, combining the direct visualization of differentially expressed proteins with the high-throughput scale of gene expression profiling. Using breast cancer as a case example, we generated comprehensive two-dimensional electrophoresis (2DE)/mass spectrometry (MS) proteomic maps of cancer (MCF-7 and HCC-38) and control (CCD-1059Sk) cell lines, identifying 1724 expressed protein spots representing 484 different protein species. The differentially expressed cell-line proteins were then mapped to mRNA transcript databases of cancer cell lines and primary breast tumors to identify candidate biomarkers that were concordantly expressed at the gene expression level. Of the top nine selected biomarker candidates, we reidentified ANX1, a protein previously reported to be differentially expressed in breast cancers and normal tissues, and validated three other novel candidates, CRAB, 6PGL, and CAZ2, as differentially expressed proteins by immunohistochemistry on breast tissue microarrays. In total, close to half (4/9) of our protein biomarker candidates were successfully validated. Our study thus illustrates how the systematic integration of proteomic and transcriptomic data from both cell line and primary tissue samples can prove advantageous for accelerating cancer biomarker discovery.