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.

Yeast screens identify the RNA polymerase II CTD and SPT5 as relevant targets of BRCA1 interaction.

BRCA1 has been implicated in numerous DNA repair pathways that maintain genome integrity, however the function responsible for its tumor suppressor activity in breast cancer remains obscure. To identify the most highly conserved of the many BRCA1 functions, we screened the evolutionarily distant eukaryote Saccharomyces cerevisiae for mutants that suppressed the G1 checkpoint arrest and lethality induced following heterologous BRCA1 expression. A genome-wide screen in the diploid deletion collection combined with a screen of ionizing radiation sensitive gene deletions identified mutants that permit growth in the presence of BRCA1. These genes delineate a metabolic mRNA pathway that temporally links transcription elongation (SPT4, SPT5, CTK1, DEF1) to nucleopore-mediated mRNA export (ASM4, MLP1, MLP2, NUP2, NUP53, NUP120, NUP133, NUP170, NUP188, POM34) and cytoplasmic mRNA decay at P-bodies (CCR4, DHH1). Strikingly, BRCA1 interacted with the phosphorylated RNA polymerase II (RNAPII) carboxy terminal domain (P-CTD), phosphorylated in the pattern specified by the CTDK-I kinase, to induce DEF1-dependent cleavage and accumulation of a RNAPII fragment containing the P-CTD. Significantly, breast cancer associated BRCT domain defects in BRCA1 that suppressed P-CTD cleavage and lethality in yeast also suppressed the physical interaction of BRCA1 with human SPT5 in breast epithelial cells, thus confirming SPT5 as a relevant target of BRCA1 interaction. Furthermore, enhanced P-CTD cleavage was observed in both yeast and human breast cells following UV-irradiation indicating a conserved eukaryotic damage response. Moreover, P-CTD cleavage in breast epithelial cells was BRCA1-dependent since damage-induced P-CTD cleavage was only observed in the mutant BRCA1 cell line HCC1937 following ectopic expression of wild type BRCA1. Finally, BRCA1, SPT5 and hyperphosphorylated RPB1 form a complex that was rapidly degraded following MMS treatment in wild type but not BRCA1 mutant breast cells. These results extend the mechanistic links between BRCA1 and transcriptional consequences in response to DNA damage and suggest an important role for RNAPII P-CTD cleavage in BRCA1-mediated cancer suppression.

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.

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.

A switch in the control of growth of the wing imaginal disks of Manduca sexta.

BACKGROUND: Insulin and ecdysone are the key extrinsic regulators of growth for the wing imaginal disks of insects. In vitro tissue culture studies have shown that these two growth regulators act synergistically: either factor alone stimulates only limited growth, but together they stimulate disks to grow at a rate identical to that observed in situ. It is generally thought that insulin signaling links growth to nutrition, and that starvation stops growth because it inhibits insulin secretion. At the end of larval life feeding stops but the disks continue to grow, so at that time disk growth has become uncoupled from nutrition. We sought to determine at exactly what point in development this uncoupling occurs. METHODOLOGY: Growth and cell proliferation in the wing imaginal disks and hemolymph carbohydrate concentrations were measured at various stages in the last larval instar under experimental conditions of starvation, ligation, rescue, and hormone treatment. PRINCIPAL FINDINGS: Here we show that in the last larval instar of M. sexta, the uncoupling of nutrition and growth occurs as the larva passes the critical weight. Before this time, starvation causes a decline in hemolymph glucose and trehalose and a cessation of wing imaginal disks growth, which can be rescued by injections of trehalose. After the critical weight the trehalose response to starvation disappears, and the expression of insulin becomes decoupled from nutrition. After the critical weight the wing disks loose their sensitivity to repression by juvenile hormone, and factors from the abdomen, but not the brain, are required to drive continued growth. CONCLUSIONS: During the last larval instar imaginal disk growth becomes decoupled from somatic growth at the time that the endocrine events of metamorphosis are initiated. These regulatory changes ensure that disk growth continues uninterrupted when the nutritive and endocrine signals undergo the drastic changes associated with metamorphosis.

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.

Use of 16S ribosomal RNA gene analyses to characterize the bacterial signature associated with poor oral health in West Virginia.

BACKGROUND: West Virginia has the worst oral health in the United States, but the reasons for this are unclear. This pilot study explored the etiology of this disparity using culture-independent analyses to identify bacterial species associated with oral disease. METHODS: Bacteria in subgingival plaque samples from twelve participants in two independent West Virginia dental-related studies were characterized using 16S rRNA gene sequencing and Human Oral Microbe Identification Microarray (HOMIM) analysis. Unifrac analysis was used to characterize phylogenetic differences between bacterial communities obtained from plaque of participants with low or high oral disease, which was further evaluated using clustering and Principal Coordinate Analysis. RESULTS: Statistically different bacterial signatures (P<0.001) were identified in subgingival plaque of individuals with low or high oral disease in West Virginia based on 16S rRNA gene sequencing. Low disease contained a high frequency of Veillonella and Streptococcus, with a moderate number of Capnocytophaga. High disease exhibited substantially increased bacterial diversity and included a large proportion of Clostridiales cluster bacteria (Selenomonas, Eubacterium, Dialister). Phylogenetic trees constructed using 16S rRNA gene sequencing revealed that Clostridiales were repeated colonizers in plaque associated with high oral disease, providing evidence that the oral environment is somehow influencing the bacterial signature linked to disease. CONCLUSIONS: Culture-independent analyses identified an atypical bacterial signature associated with high oral disease in West Virginians and provided evidence that the oral environment influenced this signature. Both findings provide insight into the etiology of the oral disparity in West Virginia.

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.

Desensitization, internalization, and signaling functions of beta-arrestins demonstrated by RNA interference.

Beta-arrestins bind to activated G protein-coupled receptor kinase-phosphorylated receptors, which leads to their desensitization with respect to G proteins, internalization via clathrin-coated pits, and signaling via a growing list of "scaffolded" pathways. To facilitate the discovery of novel adaptor and signaling roles of beta-arrestins, we have developed and validated a generally applicable interfering RNA approach for selectively suppressing beta-arrestins 1 or 2 expression by up to 95%. Beta-arrestin depletion in HEK293 cells leads to enhanced cAMP generation in response to beta(2)-adrenergic receptor stimulation, markedly reduced beta(2)-adrenergic receptor and angiotensin II receptor internalization and impaired activation of the MAP kinases ERK 1 and 2 by angiotensin II. This approach should allow discovery of novel signaling and regulatory roles for the beta-arrestins in many seven-membrane-spanning receptor systems.

Receptor-specific in vivo desensitization by the G protein-coupled receptor kinase-5 in transgenic mice.

Transgenic mice were generated with cardiac-specific overexpression of the G protein-coupled receptor kinase-5 (GRK5), a serine/threonine kinase most abundantly expressed in the heart compared with other tissues. Animals overexpressing GRK5 showed marked beta-adrenergic receptor desensitization in both the anesthetized and conscious state compared with nontransgenic control mice, while the contractile response to angiotensin II receptor stimulation was unchanged. In contrast, the angiotensin II-induced rise in contractility was significantly attenuated in transgenic mice overexpressing the beta-adrenergic receptor kinase-1, another member of the GRK family. These data suggest that myocardial overexpression of GRK5 results in selective uncoupling of G protein-coupled receptors and demonstrate that receptor specificity of the GRKs may be important in determining the physiological phenotype.

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.

The receptor kinase family: primary structure of rhodopsin kinase reveals similarities to the beta-adrenergic receptor kinase.

Light-dependent deactivation of rhodopsin as well as homologous desensitization of beta-adrenergic receptors involves receptor phosphorylation that is mediated by the highly specific protein kinases rhodopsin kinase (RK) and beta-adrenergic receptor kinase (beta ARK), respectively. We report here the cloning of a complementary DNA for RK. The deduced amino acid sequence shows a high degree of homology to beta ARK. In a phylogenetic tree constructed by comparing the catalytic domains of several protein kinases, RK and beta ARK are located on a branch close to, but separate from the cyclic nucleotide-dependent protein kinase and protein kinase C subfamilies. From the common structural features we conclude that both RK and beta ARK are members of a newly delineated gene family of guanine nucleotide-binding protein (G protein)-coupled receptor kinases that may function in diverse pathways to regulate the function of such receptors.