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.

Gait and behavior in an IL1β-mediated model of rat knee arthritis and effects of an IL1 antagonist.

Interleukin-1 beta (IL1β) is a proinflammatory cytokine that mediates arthritic pathologies. Our objectives were to evaluate pain and limb dysfunction resulting from IL1β over-expression in the rat knee and to investigate the ability of local IL1 receptor antagonist (IL1Ra) delivery to reverse-associated pathology. IL1β over-expression was induced in the right knees of 30 Wistar rats via intra-articular injection of rat fibroblasts retrovirally infected with human IL1β cDNA. A subset of animals received a 30 µl intra-articular injection of saline or human IL1Ra on day 1 after cell delivery (0.65 µg/µl hIL1Ra, n = 7 per group). Joint swelling, gait, and sensitivity were investigated over 1 week. On day 8, animals were sacrificed and joints were collected for histological evaluation. Joint inflammation and elevated levels of endogenous IL1β were observed in knees receiving IL1β-infected fibroblasts. Asymmetric gaits favoring the affected limb and heightened mechanical sensitivity (allodynia) reflected a unilateral pathology. Histopathology revealed cartilage loss on the femoral groove and condyle of affected joints. Intra-articular IL1Ra injection failed to restore gait and sensitivity to preoperative levels and did not reduce cartilage degeneration observed in histopathology. Joint swelling and degeneration subsequent to IL1β over-expression is associated limb hypersensitivity and gait compensation. Intra-articular IL1Ra delivery did not result in marked improvement for this model; this may be driven by rapid clearance of administered IL1Ra from the joint space. These results motivate work to further investigate the behavioral consequences of monoarticular arthritis and sustained release drug delivery strategies for the joint space.

Calcium dependent CAMTA1 in adult stem cell commitment to a myocardial lineage.

The phenotype of somatic cells has recently been found to be reversible. Direct reprogramming of one cell type into another has been achieved with transduction and over expression of exogenous defined transcription factors emphasizing their role in specifying cell fate. To discover early and novel endogenous transcription factors that may have a role in adult-derived stem cell acquisition of a cardiomyocyte phenotype, mesenchymal stem cells from human and mouse bone marrow and rat liver were co-cultured with neonatal cardiomyocytes as an in vitro cardiogenic microenvironment. Cell-cell communications develop between the two cell types as early as 24 hrs in co-culture and are required for elaboration of a myocardial phenotype in the stem cells 8-16 days later. These intercellular communications are associated with novel Ca(2+) oscillations in the stem cells that are synchronous with the Ca(2+) transients in adjacent cardiomyocytes and are detected in the stem cells as early as 24-48 hrs in co-culture. Early and significant up-regulation of Ca(2+)-dependent effectors, CAMTA1 and RCAN1 ensues before a myocardial program is activated. CAMTA1 loss-of-function minimizes the activation of the cardiac gene program in the stem cells. While the expression of RCAN1 suggests involvement of the well-characterized calcineurin-NFAT pathway as a response to a Ca(2+) signal, the CAMTA1 up-regulated expression as a response to such a signal in the stem cells was unknown. Cell-cell communications between the stem cells and adjacent cardiomyocytes induce Ca(2+) signals that activate a myocardial gene program in the stem cells via a novel and early Ca(2+)-dependent intermediate, up-regulation of CAMTA1.

Gene expression signatures of radiation response are specific, durable and accurate in mice and humans.

BACKGROUND: Previous work has demonstrated the potential for peripheral blood (PB) gene expression profiling for the detection of disease or environmental exposures. METHODS AND FINDINGS: We have sought to determine the impact of several variables on the PB gene expression profile of an environmental exposure, ionizing radiation, and to determine the specificity of the PB signature of radiation versus other genotoxic stresses. Neither genotype differences nor the time of PB sampling caused any lessening of the accuracy of PB signatures to predict radiation exposure, but sex difference did influence the accuracy of the prediction of radiation exposure at the lowest level (50 cGy). A PB signature of sepsis was also generated and both the PB signature of radiation and the PB signature of sepsis were found to be 100% specific at distinguishing irradiated from septic animals. We also identified human PB signatures of radiation exposure and chemotherapy treatment which distinguished irradiated patients and chemotherapy-treated individuals within a heterogeneous population with accuracies of 90% and 81%, respectively. CONCLUSIONS: We conclude that PB gene expression profiles can be identified in mice and humans that are accurate in predicting medical conditions, are specific to each condition and remain highly accurate over time.

Changes in midbrain pain receptor expression, gait and behavioral sensitivity in a rat model of radiculopathy.

Intervertebral disc herniation may contribute to inflammatory processes that associate with radicular pain and motor deficits. Molecular changes at the affected dorsal root ganglion (DRG), spinal cord, and even midbrain, have been documented in rat models of radiculopathy or nerve injury. The objective of this study was to evaluate gait and the expression of key pain receptors in the midbrain in a rodent model of radiculopathy. Radiculopathy was induced by harvesting tail nucleus pulposus (NP) and placing upon the right L5 DRG in rats (NP-treated, n=12). Tail NP was discarded in sham-operated animals (n=12). Mechanical allodynia, weight-bearing, and gait were evaluated in all animals over time. At 1 and 4 weeks after surgery, astrocyte and microglial activation was tested in DRG sections. Midbrain sections were similarly evaluated for immunoreactivity to serotonin (5HT(2B)), mu-opioid (µ-OR), and metabotropic glutamate (mGluR4 and 5) receptor antibodies. NP-treated animals placed less weight on the affected limb 1 week after surgery and experienced mechanical hypersensitivity over the duration of the study. Astroctye activation was observed at DRGs only at 4 weeks after surgery. Findings for pain receptors in the midbrain of NP-treated rats included an increased expression of 5HT(2B) at 1, but not 4 weeks; increased expression of µ-OR and mGluR5 at 1 and 4 weeks (periaqueductal gray region only); and no changes in expression of mGluR4 at any point in this study. These observations provide support for the hypothesis that the midbrain responds to DRG injury with a transient change in receptors regulating pain responses.

Analysis of the genome and transcriptome of Cryptococcus neoformans var. grubii reveals complex RNA expression and microevolution leading to virulence attenuation.

Cryptococcus neoformans is a pathogenic basidiomycetous yeast responsible for more than 600,000 deaths each year. It occurs as two serotypes (A and D) representing two varieties (i.e. grubii and neoformans, respectively). Here, we sequenced the genome and performed an RNA-Seq-based analysis of the C. neoformans var. grubii transcriptome structure. We determined the chromosomal locations, analyzed the sequence/structural features of the centromeres, and identified origins of replication. The genome was annotated based on automated and manual curation. More than 40,000 introns populating more than 99% of the expressed genes were identified. Although most of these introns are located in the coding DNA sequences (CDS), over 2,000 introns in the untranslated regions (UTRs) were also identified. Poly(A)-containing reads were employed to locate the polyadenylation sites of more than 80% of the genes. Examination of the sequences around these sites revealed a new poly(A)-site-associated motif (AUGHAH). In addition, 1,197 miscRNAs were identified. These miscRNAs can be spliced and/or polyadenylated, but do not appear to have obvious coding capacities. Finally, this genome sequence enabled a comparative analysis of strain H99 variants obtained after laboratory passage. The spectrum of mutations identified provides insights into the genetics underlying the micro-evolution of a laboratory strain, and identifies mutations involved in stress responses, mating efficiency, and virulence.

Radioactive in situ hybridization for detecting diverse gene expression patterns in tissue.

Knowing the timing, level, cellular localization, and cell type that a gene is expressed in contributes to our understanding of the function of the gene. Each of these features can be accomplished with in situ hybridization to mRNAs within cells. Here we present a radioactive in situ hybridization method modified from Clayton et al. (1988)(1) that has been working successfully in our lab for many years, especially for adult vertebrate brains(2-5). The long complementary RNA (cRNA) probes to the target sequence allows for detection of low abundance transcripts(6,7). Incorporation of radioactive nucleotides into the cRNA probes allows for further detection sensitivity of low abundance transcripts and quantitative analyses, either by light sensitive x-ray film or emulsion coated over the tissue. These detection methods provide a long-term record of target gene expression. Compared with non-radioactive probe methods, such as DIG-labeling, the radioactive probe hybridization method does not require multiple amplification steps using HRP-antibodies and/or TSA kit to detect low abundance transcripts. Therefore, this method provides a linear relation between signal intensity and targeted mRNA amounts for quantitative analysis. It allows processing 100-200 slides simultaneously. It works well for different developmental stages of embryos. Most developmental studies of gene expression use whole embryos and non-radioactive approaches(8,9), in part because embryonic tissue is more fragile than adult tissue, with less cohesion between cells, making it difficult to see boundaries between cell populations with tissue sections. In contrast, our radioactive approach, due to the larger range of sensitivity, is able to obtain higher contrast in resolution of gene expression between tissue regions, making it easier to see boundaries between populations. Using this method, researchers could reveal the possible significance of a newly identified gene, and further predict the function of the gene of interest.

Nutritional control of mRNA isoform expression during developmental arrest and recovery in C. elegans.

Nutrient availability profoundly influences gene expression. Many animal genes encode multiple transcript isoforms, yet the effect of nutrient availability on transcript isoform expression has not been studied in genome-wide fashion. When Caenorhabditis elegans larvae hatch without food, they arrest development in the first larval stage (L1 arrest). Starved larvae can survive L1 arrest for weeks, but growth and post-embryonic development are rapidly initiated in response to feeding. We used RNA-seq to characterize the transcriptome during L1 arrest and over time after feeding. Twenty-seven percent of detectable protein-coding genes were differentially expressed during recovery from L1 arrest, with the majority of changes initiating within the first hour, demonstrating widespread, acute effects of nutrient availability on gene expression. We used two independent approaches to track expression of individual exons and mRNA isoforms, and we connected changes in expression to functional consequences by mining a variety of databases. These two approaches identified an overlapping set of genes with alternative isoform expression, and they converged on common functional patterns. Genes affecting mRNA splicing and translation are regulated by alternative isoform expression, revealing post-transcriptional consequences of nutrient availability on gene regulation. We also found that phosphorylation sites are often alternatively expressed, revealing a common mode by which alternative isoform expression modifies protein function and signal transduction. Our results detail rich changes in C. elegans gene expression as larvae initiate growth and post-embryonic development, and they provide an excellent resource for ongoing investigation of transcriptional regulation and developmental physiology.

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.

The Role of TET Proteins in the Epigenetic Regulation of Neural Gene Expression and Behavior

Understanding how genes affect behavior is critical to develop precise therapies for human behavioral disorders. The ability to investigate the relationship between genes and behavior has been greatly advanced over the last few decades due to progress in gene-targeting technology. Recently, the Tet gene family was discovered and implicated in epigenetic modification of DNA methylation by converting 5-methylcytosine to 5-hydroxymethylcytosine (5hmC). 5hmC and its catalysts, the TET proteins, are highly abundant in the postnatal brain but with unclear functions. To investigate their neural functions, we generated new lines of Tet1 and Tet3 mutant mice using a gene targeting approach. We designed both mutations to cause a frameshift by deleting the largest coding exon of Tet1 (Tet1Δe4) and the catalytic domain of Tet3 (Tet3Δe7-9). As Tet1 is also highly expressed in embryonic stem cells (ESCs), we generated Tet1 homozygous deleted ESCs through sequential targeting to compare the function of Tet1 in the brain to its role in ESCs. To test our hypothesis that TET proteins epigenetically regulate transcription of key neural genes important for normal brain function, we examined transcriptional and epigenetic differences in the Tet1Δe4 mouse brain. The oxytocin receptor (OXTR), a neural gene implicated in social behaviors, is suggested to be epigenetically regulated by an unknown mechanism. Interestingly, several human studies have found associations between OXTR DNA hypermethylation and a wide spectrum of behavioral traits and neuropsychiatric disorders including autism spectrum disorders. Here we report the first evidence for an epigenetic mechanism of Oxtr transcription as expression of Oxtr is reduced in the brains of Tet1Δe4-/- mice. Likewise, the CpG island overlapping the promoter of Oxtr is hypermethylated during early embryonic development and persists into adulthood. We also discovered altered histone modifications at the hypermethylated regions, indicating the loss of TET1 has broad effects on the chromatin structure at Oxtr. Unexpectedly, we discovered an array of novel mRNA isoforms of Oxtr that are selectively reduced in Tet1Δe4-/- mice. Additionally, Tet1Δe4-/- mice display increased agonistic behaviors and impaired maternal care and short-term memory. Our findings support a novel role for TET1 in regulating Oxtr expression by preventing DNA hypermethylation and implicate TET1 in social behaviors, offering novel insight into Oxtr epigenetic regulation and its role in neuropsychiatric disorders.