Detecting separate time scales in genetic expression data.

BACKGROUND: Biological processes occur on a vast range of time scales, and many of them occur concurrently. As a result, system-wide measurements of gene expression have the potential to capture many of these processes simultaneously. The challenge however, is to separate these processes and time scales in the data. In many cases the number of processes and their time scales is unknown. This issue is particularly relevant to developmental biologists, who are interested in processes such as growth, segmentation and differentiation, which can all take place simultaneously, but on different time scales. RESULTS: We introduce a flexible and statistically rigorous method for detecting different time scales in time-series gene expression data, by identifying expression patterns that are temporally shifted between replicate datasets. We apply our approach to a Saccharomyces cerevisiae cell-cycle dataset and an Arabidopsis thaliana root developmental dataset. In both datasets our method successfully detects processes operating on several different time scales. Furthermore we show that many of these time scales can be associated with particular biological functions. CONCLUSIONS: The spatiotemporal modules identified by our method suggest the presence of multiple biological processes, acting at distinct time scales in both the Arabidopsis root and yeast. Using similar large-scale expression datasets, the identification of biological processes acting at multiple time scales in many organisms is now possible.

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.

Impact of gene variants on sex-specific regulation of human Scavenger receptor class B type 1 (SR-BI) expression in liver and association with lipid levels in a population-based study.

BACKGROUND: Several studies have noted that genetic variants of SCARB1, a lipoprotein receptor involved in reverse cholesterol transport, are associated with serum lipid levels in a sex-dependent fashion. However, the mechanism underlying this gene by sex interaction has not been explored. METHODS: We utilized both epidemiological and molecular methods to study how estrogen and gene variants interact to influence SCARB1 expression and lipid levels. Interaction between 35 SCARB1 haplotype-tagged polymorphisms and endogenous estradiol levels was assessed in 498 postmenopausal Caucasian women from the population-based Rancho Bernardo Study. We further examined associated variants with overall and SCARB1 splice variant (SR-BI and SR-BII) expression in 91 human liver tissues using quantitative real-time PCR. RESULTS: Several variants on a haplotype block spanning intron 11 to intron 12 of SCARB1 showed significant gene by estradiol interaction affecting serum lipid levels, the strongest for rs838895 with HDL-cholesterol (p=9.2x10(-4)) and triglycerides (p=1.3x10(-3)) and the triglyceride:HDL cholesterol ratio (p=2.7x10(-4)). These same variants were associated with expression of the SR-BI isoform in a sex-specific fashion, with the strongest association found among liver tissue from 52 young women<45 years old (p=0.002). CONCLUSIONS: Estrogen and SCARB1 genotype may act synergistically to regulate expression of SCARB1 isoforms and impact serum levels of HDL cholesterol and triglycerides. This work highlights the importance of considering sex-dependent effects of gene variants on serum lipid levels.

Expression in aneuploid Drosophila S2 cells.

Extensive departures from balanced gene dose in aneuploids are highly deleterious. However, we know very little about the relationship between gene copy number and expression in aneuploid cells. We determined copy number and transcript abundance (expression) genome-wide in Drosophila S2 cells by DNA-Seq and RNA-Seq. We found that S2 cells are aneuploid for >43 Mb of the genome, primarily in the range of one to five copies, and show a male genotype ( approximately two X chromosomes and four sets of autosomes, or 2X;4A). Both X chromosomes and autosomes showed expression dosage compensation. X chromosome expression was elevated in a fixed-fold manner regardless of actual gene dose. In engineering terms, the system "anticipates" the perturbation caused by X dose, rather than responding to an error caused by the perturbation. This feed-forward regulation resulted in precise dosage compensation only when X dose was half of the autosome dose. Insufficient compensation occurred at lower X chromosome dose and excessive expression occurred at higher doses. RNAi knockdown of the Male Specific Lethal complex abolished feed-forward regulation. Both autosome and X chromosome genes show Male Specific Lethal-independent compensation that fits a first order dose-response curve. Our data indicate that expression dosage compensation dampens the effect of altered DNA copy number genome-wide. For the X chromosome, compensation includes fixed and dose-dependent components.

Gene expression disruptions of organism versus organ in Drosophila species hybrids.

Hybrid dysfunctions, such as sterility, may result in part from disruptions in the regulation of gene expression. Studies of hybrids within the Drosophila simulans clade have reported genes expressed above or below the expression observed in their parent species, and such misexpression is associated with male sterility in multigenerational backcross hybrids. However, these studies often examined whole bodies rather than testes or had limited replication using less-sensitive but global techniques. Here, we use a new RNA isolation technique to re-examine hybrid gene expression disruptions in both testes and whole bodies from single Drosophila males by real-time quantitative RT-PCR. We find two early-spermatogenesis transcripts are underexpressed in hybrid whole-bodies but not in assays of testes alone, while two late-spermatogenesis transcripts seem to be underexpressed in both whole-bodies and testes alone. Although the number of transcripts surveyed is limited, these results provide some support for a previous hypothesis that the spermatogenesis pathway in these sterile hybrids may be disrupted sometime after the expression of the early meiotic arrest genes.

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.

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.

Overexpression of the cardiac beta(2)-adrenergic receptor and expression of a beta-adrenergic receptor kinase-1 (betaARK1) inhibitor both increase myocardial contractility but have differential effects on susceptibility to ischemic injury.

Cardiac beta(2)-adrenergic receptor (beta(2)AR) overexpression is a potential contractile therapy for heart failure. Cardiac contractility was elevated in mice overexpressing beta(2)ARs (TG4s) with no adverse effects under normal conditions. To assess the consequences of beta(2)AR overexpression during ischemia, perfused hearts from TG4 and wild-type mice were subjected to 20-minute ischemia and 40-minute reperfusion. During ischemia, ATP and pH fell lower in TG4 hearts than wild type. Ischemic injury was greater in TG4 hearts, as indicated by lower postischemic recoveries of contractile function, ATP, and phosphocreatine. Because beta(2)ARs, unlike beta(1)ARs, couple to G(i) as well as G(s), we pretreated mice with the G(i) inhibitor pertussis toxin (PTX). PTX treatment increased basal contractility in TG4 hearts and abolished the contractile resistance to isoproterenol. During ischemia, ATP fell lower in TG4+PTX than in TG4 hearts. Recoveries of contractile function and ATP were lower in TG4+PTX than in TG4 hearts. We also studied mice that overexpressed either betaARK1 (TGbetaARK1) or a betaARK1 inhibitor (TGbetaARKct). Recoveries of function, ATP, and phosphocreatine were higher in TGbetaARK1 hearts than in wild-type hearts. Despite basal contractility being elevated in TGbetaARKct hearts to the same level as that of TG4s, ischemic injury was not increased. In summary, beta(2)AR overexpression increased ischemic injury, whereas betaARK1 overexpression was protective. Ischemic injury in the beta(2)AR overexpressors was exacerbated by PTX treatment, implying that it was G(s) not G(i) activity that enhanced injury. Unlike beta(2)AR overexpression, basal contractility was increased by betaARK1 inhibitor expression without increasing ischemic injury, thus implicating a safer potential therapy for heart failure.

Expression of a beta-adrenergic receptor kinase 1 inhibitor prevents the development of myocardial failure in gene-targeted mice.

Heart failure is accompanied by severely impaired beta-adrenergic receptor (betaAR) function, which includes loss of betaAR density and functional uncoupling of remaining receptors. An important mechanism for the rapid desensitization of betaAR function is agonist-stimulated receptor phosphorylation by the betaAR kinase (betaARK1), an enzyme known to be elevated in failing human heart tissue. To investigate whether alterations in betaAR function contribute to the development of myocardial failure, transgenic mice with cardiac-restricted overexpression of either a peptide inhibitor of betaARK1 or the beta2AR were mated into a genetic model of murine heart failure (MLP-/-). In vivo cardiac function was assessed by echocardiography and cardiac catheterization. Both MLP-/- and MLP-/-/beta2AR mice had enlarged left ventricular (LV) chambers with significantly reduced fractional shortening and mean velocity of circumferential fiber shortening. In contrast, MLP-/-/betaARKct mice had normal LV chamber size and function. Basal LV contractility in the MLP-/-/betaARKct mice, as measured by LV dP/dtmax, was increased significantly compared with the MLP-/- mice but less than controls. Importantly, heightened betaAR desensitization in the MLP-/- mice, measured in vivo (responsiveness to isoproterenol) and in vitro (isoproterenol-stimulated membrane adenylyl cyclase activity), was completely reversed with overexpression of the betaARK1 inhibitor. We report here the striking finding that overexpression of this inhibitor prevents the development of cardiomyopathy in this murine model of heart failure. These findings implicate abnormal betaAR-G protein coupling in the pathogenesis of the failing heart and point the way toward development of agents to inhibit betaARK1 as a novel mode of therapy.

Myocardial expression of a constitutively active alpha 1B-adrenergic receptor in transgenic mice induces cardiac hypertrophy.

Transgenic mice were generated by using the alpha-myosin heavy chain promoter coupled to the coding sequence of a constitutively active mutant alpha 1B-adrenergic receptor (AR). These transgenic animals demonstrated cardiac-specific expression of this alpha 1-AR with resultant activation of phospholipase C as shown by increased myocardial diacylglycerol content. A phenotype consistent with cardiac hypertrophy developed in adult transgenic mice with increased heart/body weight ratios, myocyte cross-sectional areas, and ventricular atrial natriuretic factor mRNA levels relative to nontransgenic controls. These transgenic animals may provide insight into the biochemical triggers that induce hypertrophy in cardiac disease and serve as a convenient experimental model for studies of this condition.

Cloning and expression of a human kidney cDNA for an alpha 2-adrenergic receptor subtype.

An alpha 2-adrenergic receptor subtype has been cloned from a human kidney cDNA library using the gene for the human platelet alpha 2-adrenergic receptor as a probe. The deduced amino acid sequence resembles the human platelet alpha 2-adrenergic receptor and is consistent with the structure of other members of the family of guanine nucleotide-binding protein-coupled receptors. The cDNA was expressed in a mammalian cell line (COS-7), and the alpha 2-adrenergic ligand [3H]rauwolscine was bound. Competition curve analysis with a variety of adrenergic ligands suggests that this cDNA clone represents the alpha 2B-adrenergic receptor. The gene for this receptor is on human chromosome 4, whereas the gene for the human platelet alpha 2-adrenergic receptor (alpha 2A) lies on chromosome 10. This ability to express the receptor in mammalian cells, free of other adrenergic receptor subtypes, should help in developing more selective alpha-adrenergic ligands.

Genetic variants in the TEP1 gene are associated with prostate cancer risk and recurrence.

BACKGROUND: Telomere-related genes play an important role in carcinogenesis and progression of prostate cancer (PCa). It is not fully understood whether genetic variations in telomere-related genes are associated with development and progression in PCa patients. METHODS: Six potentially functional single-nucleotide polymorphisms (SNPs) of three key telomere-related genes were evaluated in 1015 PCa cases and 1052 cancer-free controls, to test their associations with risk of PCa. Among 426 PCa patients who underwent radical prostatectomy (RP), the prognostic significance of the studied SNPs on biochemical recurrence (BCR) was also assessed using the Kaplan-Meier analysis and Cox proportional hazards regression model. The relative telomere lengths (RTLs) were measured in peripheral blood leukocytes using real-time PCR in the RP patients. RESULTS: TEP1 rs1760904 AG/AA genotypes were significantly associated with a decreased risk of PCa (odds ratio (OR): 0.77, 95% confidence interval (CI): 0.64-0.93, P=0.005) compared with the GG genotype. By using median RTL as a cutoff level, RP patients with TEP1 rs1760904 AG/AA genotypes tended to have a longer RTL than those with the GG genotype (OR: 1.55, 95% CI: 1.04-2.30, P=0.031). A significant interaction between TEP1 rs1713418 and age in modifying PCa risk was observed (P=0.005). After adjustment for clinicopathologic risk factors, the presence of heterozygotes or rare homozygotes of TEP1 rs1760904 and TNKS2 rs1539042 were associated with BCR in the RP cohorts (hazard ratio: 0.53, 95% CI: 0.36-0.79, P=0.002 and hazard ratio: 1.67, 95% CI: 1.07-2.48, P=0.017, respectively). CONCLUSIONS: These data suggest that genetic variations in the TEP1 gene may be biomarkers for risk of PCa and BCR after RP.

A complex intronic enhancer regulates expression of the CFTR gene by direct interaction with the promoter.

Genes can maintain spatiotemporal expression patterns by long-range interactions between cis-acting elements. The cystic fibrosis transmembrane conductance regulator gene (CFTR) is expressed primarily in epithelial cells. An element located within a DNase I-hypersensitive site (DHS) 10 kb into the first intron was previously shown to augment CFTR promoter activity in a tissue-specific manner. Here, we reveal the mechanism by which this element influences CFTR transcription. We employed a high-resolution method of mapping DHS using tiled microarrays to accurately locate the intron 1 DHS. Transfection of promoter-reporter constructs demonstrated that the element displays classical tissue-specific enhancer properties and can independently recruit factors necessary for transcription initiation. In vitro DNase I footprinting analysis identified a protected region that corresponds to a conserved, predicted binding site for hepatocyte nuclear factor 1 (HNF1). We demonstrate by electromobility shift assays (EMSA) and chromatin immunoprecipitation (ChIP) that HNF1 binds to this element both in vitro and in vivo. Moreover, using chromosome conformation capture (3C) analysis, we show that this element interacts with the CFTR promoter in CFTR-expressing cells. These data provide the first insight into the three- dimensional (3D) structure of the CFTR locus and confirm the contribution of intronic cis-acting elements to the regulation of CFTR gene expression.

Inactivation of the von Hippel-Lindau tumor suppressor leads to selective expression of a human endogenous retrovirus in kidney cancer.

A human endogenous retrovirus type E (HERV-E) was recently found to be selectively expressed in most renal cell carcinomas (RCCs). Importantly, antigens derived from this provirus are immunogenic, stimulating cytotoxic T cells that kill RCC cells in vitro and in vivo. Here, we show HERV-E expression is restricted to the clear cell subtype of RCC (ccRCC) characterized by an inactivation of the von Hippel-Lindau (VHL) tumor-suppressor gene with subsequent stabilization of hypoxia-inducible transcription factors (HIFs)-1α and -2α. HERV-E expression in ccRCC linearly correlated with HIF-2α levels and could be silenced in tumor cells by either transfection of normal VHL or small interfering RNA inhibition of HIF-2α. Using chromatin immunoprecipitation, we demonstrated that HIF-2α can serve as transcriptional factor for HERV-E by binding with HIF response element (HRE) localized in the proviral 5' long terminal repeat (LTR). Remarkably, the LTR was found to be hypomethylated only in HERV-E-expressing ccRCC while other tumors and normal tissues possessed a hypermethylated LTR preventing proviral expression. Taken altogether, these findings provide the first evidence that inactivation of a tumor suppressor gene can result in aberrant proviral expression in a human tumor and give insights needed for translational research aimed at boosting human immunity against antigenic components of this HERV-E.

The upstream enhancer elements of the G6PC promoter are critical for optimal G6PC expression in murine glycogen storage disease type Ia.

Glycogen storage disease type-Ia (GSD-Ia) patients deficient in glucose-6-phosphatase-α (G6Pase-α or G6PC) manifest impaired glucose homeostasis characterized by fasting hypoglycemia, growth retardation, hepatomegaly, nephromegaly, hyperlipidemia, hyperuricemia, and lactic acidemia. Two efficacious recombinant adeno-associated virus pseudotype 2/8 (rAAV8) vectors expressing human G6Pase-α have been independently developed. One is a single-stranded vector containing a 2864-bp of the G6PC promoter/enhancer (rAAV8-GPE) and the other is a double-stranded vector containing a shorter 382-bp minimal G6PC promoter/enhancer (rAAV8-miGPE). To identify the best construct, a direct comparison of the rAAV8-GPE and the rAAV8-miGPE vectors was initiated to determine the best vector to take forward into clinical trials. We show that the rAAV8-GPE vector directed significantly higher levels of hepatic G6Pase-α expression, achieved greater reduction in hepatic glycogen accumulation, and led to a better toleration of fasting in GSD-Ia mice than the rAAV8-miGPE vector. Our results indicated that additional control elements in the rAAV8-GPE vector outweigh the gains from the double-stranded rAAV8-miGPE transduction efficiency, and that the rAAV8-GPE vector is the current choice for clinical translation in human GSD-Ia.