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.

Latent factor analysis to discover pathway-associated putative segmental aneuploidies in human cancers.

Tumor microenvironmental stresses, such as hypoxia and lactic acidosis, play important roles in tumor progression. Although gene signatures reflecting the influence of these stresses are powerful approaches to link expression with phenotypes, they do not fully reflect the complexity of human cancers. Here, we describe the use of latent factor models to further dissect the stress gene signatures in a breast cancer expression dataset. The genes in these latent factors are coordinately expressed in tumors and depict distinct, interacting components of the biological processes. The genes in several latent factors are highly enriched in chromosomal locations. When these factors are analyzed in independent datasets with gene expression and array CGH data, the expression values of these factors are highly correlated with copy number alterations (CNAs) of the corresponding BAC clones in both the cell lines and tumors. Therefore, variation in the expression of these pathway-associated factors is at least partially caused by variation in gene dosage and CNAs among breast cancers. We have also found the expression of two latent factors without any chromosomal enrichment is highly associated with 12q CNA, likely an instance of "trans"-variations in which CNA leads to the variations in gene expression outside of the CNA region. In addition, we have found that factor 26 (1q CNA) is negatively correlated with HIF-1alpha protein and hypoxia pathways in breast tumors and cell lines. This agrees with, and for the first time links, known good prognosis associated with both a low hypoxia signature and the presence of CNA in this region. Taken together, these results suggest the possibility that tumor segmental aneuploidy makes significant contributions to variation in the lactic acidosis/hypoxia gene signatures in human cancers and demonstrate that latent factor analysis is a powerful means to uncover such a linkage.

MRP3: a molecular target for human glioblastoma multiforme immunotherapy.

BACKGROUND: Glioblastoma multiforme (GBM) is refractory to conventional therapies. To overcome the problem of heterogeneity, more brain tumor markers are required for prognosis and targeted therapy. We have identified and validated a promising molecular therapeutic target that is expressed by GBM: human multidrug-resistance protein 3 (MRP3). METHODS: We investigated MRP3 by genetic and immunohistochemical (IHC) analysis of human gliomas to determine the incidence, distribution, and localization of MRP3 antigens in GBM and their potential correlation with survival. To determine MRP3 mRNA transcript and protein expression levels, we performed quantitative RT-PCR, raising MRP3-specific antibodies, and IHC analysis with biopsies of newly diagnosed GBM patients. We used univariate and multivariate analyses to assess the correlation of RNA expression and IHC of MRP3 with patient survival, with and without adjustment for age, extent of resection, and KPS. RESULTS: Real-time PCR results from 67 GBM biopsies indicated that 59/67 (88%) samples highly expressed MRP3 mRNA transcripts, in contrast with minimal expression in normal brain samples. Rabbit polyvalent and murine monoclonal antibodies generated against an extracellular span of MRP3 protein demonstrated reactivity with defined MRP3-expressing cell lines and GBM patient biopsies by Western blotting and FACS analyses, the latter establishing cell surface MRP3 protein expression. IHC evaluation of 46 GBM biopsy samples with anti-MRP3 IgG revealed MRP3 in a primarily membranous and cytoplasmic pattern in 42 (91%) of the 46 samples. Relative RNA expression was a strong predictor of survival for newly diagnosed GBM patients. Hazard of death for GBM patients with high levels of MRP3 RNA expression was 2.71 (95% CI: 1.54-4.80) times that of patients with low/moderate levels (p = 0.002). CONCLUSIONS: Human GBMs overexpress MRP3 at both mRNA and protein levels, and elevated MRP3 mRNA levels in GBM biopsy samples correlated with a higher risk of death. These data suggest that the tumor-associated antigen MRP3 has potential use for prognosis and as a target for malignant glioma immunotherapy.

An integrated approach to the prediction of chemotherapeutic response in patients with breast cancer.

BACKGROUND: A major challenge in oncology is the selection of the most effective chemotherapeutic agents for individual patients, while the administration of ineffective chemotherapy increases mortality and decreases quality of life in cancer patients. This emphasizes the need to evaluate every patient's probability of responding to each chemotherapeutic agent and limiting the agents used to those most likely to be effective. METHODS AND RESULTS: Using gene expression data on the NCI-60 and corresponding drug sensitivity, mRNA and microRNA profiles were developed representing sensitivity to individual chemotherapeutic agents. The mRNA signatures were tested in an independent cohort of 133 breast cancer patients treated with the TFAC (paclitaxel, 5-fluorouracil, adriamycin, and cyclophosphamide) chemotherapy regimen. To further dissect the biology of resistance, we applied signatures of oncogenic pathway activation and performed hierarchical clustering. We then used mRNA signatures of chemotherapy sensitivity to identify alternative therapeutics for patients resistant to TFAC. Profiles from mRNA and microRNA expression data represent distinct biologic mechanisms of resistance to common cytotoxic agents. The individual mRNA signatures were validated in an independent dataset of breast tumors (P = 0.002, NPV = 82%). When the accuracy of the signatures was analyzed based on molecular variables, the predictive ability was found to be greater in basal-like than non basal-like patients (P = 0.03 and P = 0.06). Samples from patients with co-activated Myc and E2F represented the cohort with the lowest percentage (8%) of responders. Using mRNA signatures of sensitivity to other cytotoxic agents, we predict that TFAC non-responders are more likely to be sensitive to docetaxel (P = 0.04), representing a viable alternative therapy. CONCLUSIONS: Our results suggest that the optimal strategy for chemotherapy sensitivity prediction integrates molecular variables such as ER and HER2 status with corresponding microRNA and mRNA expression profiles. Importantly, we also present evidence to support the concept that analysis of molecular variables can present a rational strategy to identifying alternative therapeutic opportunities.

The effects of sleep hypoxia on coagulant factors and hepatic inflammation in emphysematous rats.

OBJECTIVES: To develop a sleep hypoxia (SH) in emphysema (SHE) rat model and to explore whether SHE results in more severe hepatic inflammation than emphysema alone and whether the inflammation changes levels of coagulant/anticoagulant factors synthesized in the liver. METHODS: Seventy-five rats were put into 5 groups: SH control (SHCtrl), treated with sham smoke exposure (16 weeks) and SH exposure (12.5% O(2), 3 h/d, latter 8 weeks); emphysema control (ECtrl), smoke exposure and sham SH exposure (21% O(2)); short SHE (SHEShort), smoke exposure and short SH exposure (1.5 h/d); mild SHE (SHEMild), smoke exposure and mild SH exposure (15% O(2)); standard SHE (SHEStand), smoke exposure and SH exposure. Therefore, ECtrl, SHEShort, SHEMild and SHEStand group were among emphysematous groups. Arterial blood gas (ABG) data was obtained during preliminary tests. After exposure, hepatic inflammation (interleukin -6 [IL-6] mRNA and protein, tumor necrosis factor α [TNFα] mRNA and protein) and liver coagulant/anticoagulant factors (antithrombin [AT], fibrinogen [FIB] and Factor VIII [F VIII]) were evaluated. SPSS 11.5 software was used for statistical analysis. RESULTS: Characteristics of emphysema were obvious in emphysematous groups and ABGs reached SH criteria on hypoxia exposure. Hepatic inflammation parameters and coagulant factors are the lowest in SHCtrl and the highest in SHEStand while AT is the highest in SHCtrl and the lowest in SHEStand. Inflammatory cytokines of liver correlate well with coagulant factors positively and with AT negatively. CONCLUSIONS: When SH is combined with emphysema, hepatic inflammation and coagulability enhance each other synergistically and produce a more significant liver-derivative inflammatory and prothrombotic status.

Genome-wide mRNA expression correlates of viral control in CD4+ T-cells from HIV-1-infected individuals.

There is great interindividual variability in HIV-1 viral setpoint after seroconversion, some of which is known to be due to genetic differences among infected individuals. Here, our focus is on determining, genome-wide, the contribution of variable gene expression to viral control, and to relate it to genomic DNA polymorphism. RNA was extracted from purified CD4+ T-cells from 137 HIV-1 seroconverters, 16 elite controllers, and 3 healthy blood donors. Expression levels of more than 48,000 mRNA transcripts were assessed by the Human-6 v3 Expression BeadChips (Illumina). Genome-wide SNP data was generated from genomic DNA using the HumanHap550 Genotyping BeadChip (Illumina). We observed two distinct profiles with 260 genes differentially expressed depending on HIV-1 viral load. There was significant upregulation of expression of interferon stimulated genes with increasing viral load, including genes of the intrinsic antiretroviral defense. Upon successful antiretroviral treatment, the transcriptome profile of previously viremic individuals reverted to a pattern comparable to that of elite controllers and of uninfected individuals. Genome-wide evaluation of cis-acting SNPs identified genetic variants modulating expression of 190 genes. Those were compared to the genes whose expression was found associated with viral load: expression of one interferon stimulated gene, OAS1, was found to be regulated by a SNP (rs3177979, p = 4.9E-12); however, we could not detect an independent association of the SNP with viral setpoint. Thus, this study represents an attempt to integrate genome-wide SNP signals with genome-wide expression profiles in the search for biological correlates of HIV-1 control. It underscores the paradox of the association between increasing levels of viral load and greater expression of antiviral defense pathways. It also shows that elite controllers do not have a fully distinctive mRNA expression pattern in CD4+ T cells. Overall, changes in global RNA expression reflect responses to viral replication rather than a mechanism that might explain viral control.

Cryptococcal cell morphology affects host cell interactions and pathogenicity.

Cryptococcus neoformans is a common life-threatening human fungal pathogen. The size of cryptococcal cells is typically 5 to 10 microm. Cell enlargement was observed in vivo, producing cells up to 100 microm. These morphological changes in cell size affected pathogenicity via reducing phagocytosis by host mononuclear cells, increasing resistance to oxidative and nitrosative stress, and correlated with reduced penetration of the central nervous system. Cell enlargement was stimulated by coinfection with strains of opposite mating type, and ste3aDelta pheromone receptor mutant strains had reduced cell enlargement. Finally, analysis of DNA content in this novel cell type revealed that these enlarged cells were polyploid, uninucleate, and produced daughter cells in vivo. These results describe a novel mechanism by which C. neoformans evades host phagocytosis to allow survival of a subset of the population at early stages of infection. Thus, morphological changes play unique and specialized roles during infection.

Two genes on A/J chromosome 18 are associated with susceptibility to Staphylococcus aureus infection by combined microarray and QTL analyses.

Although it has recently been shown that A/J mice are highly susceptible to Staphylococcus aureus sepsis as compared to C57BL/6J, the specific genes responsible for this differential phenotype are unknown. Using chromosome substitution strains (CSS), we found that loci on chromosomes 8, 11, and 18 influence susceptibility to S. aureus sepsis in A/J mice. We then used two candidate gene selection strategies to identify genes on these three chromosomes associated with S. aureus susceptibility, and targeted genes identified by both gene selection strategies. First, we used whole genome transcription profiling to identify 191 (56 on chr. 8, 100 on chr. 11, and 35 on chr. 18) genes on our three chromosomes of interest that are differentially expressed between S. aureus-infected A/J and C57BL/6J. Second, we identified two significant quantitative trait loci (QTL) for survival post-infection on chr. 18 using N(2) backcross mice (F(1) [C18A]xC57BL/6J). Ten genes on chr. 18 (March3, Cep120, Chmp1b, Dcp2, Dtwd2, Isoc1, Lman1, Spire1, Tnfaip8, and Seh1l) mapped to the two significant QTL regions and were also identified by the expression array selection strategy. Using real-time PCR, 6 of these 10 genes (Chmp1b, Dtwd2, Isoc1, Lman1, Tnfaip8, and Seh1l) showed significantly different expression levels between S. aureus-infected A/J and C57BL/6J. For two (Tnfaip8 and Seh1l) of these 6 genes, siRNA-mediated knockdown of gene expression in S. aureus-challenged RAW264.7 macrophages induced significant changes in the cytokine response (IL-1 beta and GM-CSF) compared to negative controls. These cytokine response changes were consistent with those seen in S. aureus-challenged peritoneal macrophages from CSS 18 mice (which contain A/J chromosome 18 but are otherwise C57BL/6J), but not C57BL/6J mice. These findings suggest that two genes, Tnfaip8 and Seh1l, may contribute to susceptibility to S. aureus in A/J mice, and represent promising candidates for human genetic susceptibility studies.

Hybrid transgenic mice reveal in vivo specificity of G protein-coupled receptor kinases in the heart.

G protein-coupled receptor kinases (GRKs) phosphorylate activated G protein-coupled receptors, including alpha(1B)-adrenergic receptors (ARs), resulting in desensitization. In vivo analysis of GRK substrate selectivity has been limited. Therefore, we generated hybrid transgenic mice with myocardium-targeted overexpression of 1 of 3 GRKs expressed in the heart (GRK2 [commonly known as the beta-AR kinase 1], GRK3, or GRK5) with concomitant cardiac expression of a constitutively activated mutant (CAM) or wild-type alpha(1B)AR. Transgenic mice with cardiac CAMalpha(1B)AR overexpression had enhanced myocardial alpha(1)AR signaling and elevated heart-to-body weight ratios with ventricular atrial natriuretic factor expression denoting myocardial hypertrophy. Transgenic mouse hearts overexpressing only GRK2, GRK3, or GRK5 had no hypertrophy. In hybrid transgenic mice, enhanced in vivo signaling through CAMalpha(1B)ARs, as measured by myocardial diacylglycerol content, was attenuated by concomitant overexpression of GRK3 but not GRK2 or GRK5. CAMalpha(1B)AR-induced hypertrophy and ventricular atrial natriuretic factor expression were significantly attenuated with either concurrent GRK3 or GRK5 overexpression. Similar GRK selectivity was seen in hybrid transgenic mice with wild-type alpha(1B)AR overexpression concurrently with a GRK. GRK2 overexpression was without effect on any in vivo CAM or wild-type alpha(1B)AR cardiac phenotype, which is in contrast to previously reported in vitro findings. Furthermore, endogenous myocardial alpha(1)AR mitogen-activated protein kinase signaling in single-GRK transgenic mice also exhibited selectivity, as GRK3 and GRK5 desensitized in vivo alpha(1)AR mitogen-activated protein kinase responses that were unaffected by GRK2 overexpression. Thus, these results demonstrate that GRKs differentially interact with alpha(1B)ARs in vivo such that GRK3 desensitizes all alpha(1B)AR signaling, whereas GRK5 has partial effects and, most interestingly, GRK2 has no effect on in vivo alpha(1B)AR signaling in the heart.

WNT3 is a biomarker capable of predicting the definitive endoderm differentiation potential of hESCs.

Generation of functional cells from human pluripotent stem cells (PSCs) through in vitro differentiation is a promising approach for drug screening and cell therapy. However, the observed large and unavoidable variation in the differentiation potential of different human embryonic stem cell (hESC)/induced PSC (iPSC) lines makes the selection of an appropriate cell line for the differentiation of a particular cell lineage difficult. Here, we report identification of WNT3 as a biomarker capable of predicting definitive endoderm (DE) differentiation potential of hESCs. We show that the mRNA level of WNT3 in hESCs correlates with their DE differentiation efficiency. In addition, manipulations of hESCs through WNT3 knockdown or overexpression can respectively inhibit or promote DE differentiation in a WNT3 level-dependent manner. Finally, analysis of several hESC lines based on their WNT3 expression levels allowed accurate prediction of their DE differentiation potential. Collectively, our study supports the notion that WNT3 can serve as a biomarker for predicting DE differentiation potential of hESCs.

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.

Sensitive and precise quantification of insulin-like mRNA expression in Caenorhabditis elegans.

Insulin-like signaling regulates developmental arrest, stress resistance and lifespan in the nematode Caenorhabditis elegans. However, the genome encodes 40 insulin-like peptides, and the regulation and function of individual peptides is largely uncharacterized. We used the nCounter platform to measure mRNA expression of all 40 insulin-like peptides as well as the insulin-like receptor daf-2, its transcriptional effector daf-16, and the daf-16 target gene sod-3. We validated the platform using 53 RNA samples previously characterized by high density oligonucleotide microarray analysis. For this set of genes and the standard nCounter protocol, sensitivity and precision were comparable between the two platforms. We optimized conditions of the nCounter assay by varying the mass of total RNA used for hybridization, thereby increasing sensitivity up to 50-fold and reducing the median coefficient of variation as much as 4-fold. We used deletion mutants to demonstrate specificity of the assay, and we used optimized conditions to assay insulin-like gene expression throughout the C. elegans life cycle. We detected expression for nearly all insulin-like genes and find that they are expressed in a variety of distinct patterns suggesting complexity of regulation and specificity of function. We identified insulin-like genes that are specifically expressed during developmental arrest, larval development, adulthood and embryogenesis. These results demonstrate that the nCounter platform provides a powerful approach to analyzing insulin-like gene expression dynamics, and they suggest hypotheses about the function of individual insulin-like genes.

Polymorphisms in the kinesin-like factor 1 B gene and risk of epithelial ovarian cancer in Eastern Chinese women.

The kinesin-like factor 1 B (KIF1B) gene plays an important role in the process of apoptosis and the transformation and progression of malignant cells. Genetic variations in KIF1B may contribute to risk of epithelial ovarian cancer (EOC). In this study of 1,324 EOC patients and 1,386 cancer-free female controls, we investigated associations between two potentially functional single nucleotide polymorphisms in KIF1B and EOC risk by the conditional logistic regression analysis. General linear regression model was used to evaluate the correlation between the number of variant alleles and KIF1B mRNA expression levels. We found that the rs17401966 variant AG/GG genotypes were significantly associated with a decreased risk of EOC (adjusted odds ratio (OR) = 0.81, 95 % confidence interval (CI) = 0.68-0.97), compared with the AA genotype, but no associations were observed for rs1002076. Women who carried both rs17401966 AG/GG and rs1002076 AG/AA genotypes of KIF1B had a 0.82-fold decreased risk (adjusted 95 % CI = 0.69-0.97), compared with others. Additionally, there was no evidence of possible interactions between about-mentioned co-variants. Further genotype-phenotype correlation analysis indicated that the number of rs17401966 variant G allele was significantly associated with KIF1B mRNA expression levels (P for GLM = 0.003 and 0.001 in all and Chinese subjects, respectively), with GG carriers having the lowest level of KIF1B mRNA expression. Taken together, the rs17401966 polymorphism likely regulates KIF1B mRNA expression and thus may be associated with EOC risk in Eastern Chinese women. Larger, independent studies are warranted to validate our findings.

Pax3 expression enhances PDGF-B-induced brainstem gliomagenesis and characterizes a subset of brainstem glioma.

High-grade Brainstem Glioma (BSG), also known as Diffuse Intrinsic Pontine Glioma (DIPG), is an incurable pediatric brain cancer. Increasing evidence supports the existence of regional differences in gliomagenesis such that BSG is considered a distinct disease from glioma of the cerebral cortex (CG). In an effort to elucidate unique characteristics of BSG, we conducted expression analysis of mouse PDGF-B-driven BSG and CG initiated in Nestin progenitor cells and identified a short list of expression changes specific to the brainstem gliomagenesis process, including abnormal upregulation of paired box 3 (Pax3). In the neonatal mouse brain, Pax3 expression marks a subset of brainstem progenitor cells, while it is absent from the cerebral cortex, mirroring its regional expression in glioma. Ectopic expression of Pax3 in normal brainstem progenitors in vitro shows that Pax3 inhibits apoptosis. Pax3-induced inhibition of apoptosis is p53-dependent, however, and in the absence of p53, Pax3 promotes proliferation of brainstem progenitors. In vivo, Pax3 enhances PDGF-B-driven gliomagenesis by shortening tumor latency and increasing tumor penetrance and grade, in a region-specific manner, while loss of Pax3 function extends survival of PDGF-B-driven;p53-deficient BSG-bearing mice by 33%. Importantly, Pax3 is regionally expressed in human glioma as well, with high PAX3 mRNA characterizing 40% of human BSG, revealing a subset of tumors that significantly associates with PDGFRA alterations, amplifications of cell cycle regulatory genes, and is exclusive of ACVR1 mutations. Collectively, these data suggest that regional Pax3 expression not only marks a novel subset of BSG but also contributes to PDGF-B-induced brainstem gliomagenesis.

Traumatic brain injury exacerbates neurodegenerative pathology: improvement with an apolipoprotein E-based therapeutic.

Cognitive impairment is common following traumatic brain injury (TBI), and neuroinflammatory mechanisms may predispose to the development of neurodegenerative disease. Apolipoprotein E (apoE) polymorphisms modify neuroinflammatory responses, and influence both outcome from acute brain injury and the risk of developing neurodegenerative disease. We demonstrate that TBI accelerates neurodegenerative pathology in double-transgenic animals expressing the common human apoE alleles and mutated amyloid precursor protein, and that pathology is exacerbated in the presence of the apoE4 allele. The administration of an apoE-mimetic peptide markedly reduced the development of neurodegenerative pathology in mice homozygous for apoE3 as well as apoE3/E4 heterozygotes. These results demonstrate that TBI accelerates the cardinal neuropathological features of neurodegenerative disease, and establishes the potential for apoE mimetic therapies in reducing pathology associated with neurodegeneration.