Unite and conquer: univariate and multivariate approaches for finding differentially expressed gene sets

Motivation: Recently, many univariate and several multivariate approaches have been suggested for testing differential expression of gene sets between different phenotypes. However, despite a wealth of literature studying their performance on simulated and real biological data, still there is a need to quantify their relative performance when they are testing different null hypotheses.

Results: In this article, we compare the performance of univariate and multivariate tests on both simulated and biological data. In the simulation study we demonstrate that high correlations equally affect the power of both, univariate as well as multivariate tests. In addition, for most of them the power is similarly affected by the dimensionality of the gene set and by the percentage of genes in the set, for which expression is changing between two phenotypes. The application of different test statistics to biological data reveals that three statistics (sum of squared t-tests, Hotelling's T2, N-statistic), testing different null hypotheses, find some common but also some complementing differentially expressed gene sets under specific settings. This demonstrates that due to complementing null hypotheses each test projects on different aspects of the data and for the analysis of biological data it is beneficial to use all three tests simultaneously instead of focusing exclusively on just one.

The 67-kd laminin receptor is preferentially expressed by proliferating retinal vessels in a murine model of ischemic retinopathy.

Endothelial cell association with vascular basement membranes is complex and plays a critical role in regulation of cell adhesion and proliferation. The interaction between the membrane-associated 67-kd receptor (67LR) and the basement membrane protein laminin has been studied in several cell systems where it was shown to be crucial for adhesion and attachment during angiogenesis. As angiogenesis in the pathological setting of proliferative retinopathy is a major cause of blindness in the Western world we examined the expression of 67LR in a murine model of hyperoxia-induced retinopathy that exhibits retinal neovascularization. Mice exposed to hyperoxia for 5 days starting at postnatal day 7 (P7) and returned to room air (at P12) showed closure of the central retinal vasculature. In response to the ensuing retinal ischemia, there was consistent preretinal neovascularization starting around P17, which persisted until P21, after which the new vessels regressed. Immunohistochemistry was performed on these retinas using an antibody specific for 67LR. At P12, immunoreactivity for 67LR was absent in the retina, but by P17 it was observed in preretinal proliferating vessels and also within the adjacent intraretinal vasculature. Intraretinal 67LR immunoreactivity diminished beyond P17 until by P21 immunoreactivity was almost completely absent, although it persisted in the preretinal vasculature. Control P17 mice (not exposed to hyperoxia) failed to demonstrate any 67LR immunoreactivity in their retinas. Parallel in situ hybridization studies demonstrated 67LR gene expression in the retinal ganglion cells of control and hyperoxia-exposed mice. In addition, the neovascular intra- and preretinal vessels of hyperoxia-treated P17 and P21 mice labeled strongly for 67LR mRNA. This study has characterized 67LR immunolocalization and gene expression in a murine model of ischemic retinopathy. Results suggest that, although the 67LR gene is expressed at high levels in the retinal ganglion cells, the mature receptor protein is preferentially localized to the proliferating retinal vasculature and is almost completely absent from quiescent vessels. The differential expression of 67LR between proliferating and quiescent retinal vessels suggests that this laminin receptor is an important and novel target for future chemotherapeutic intervention during proliferative vasculopathies.

Electrophysiological, pharmacological and molecular profile of the transient outward rectifying conductance in neonatal rat symapthetic preganglionic neurons in vitro.

Transient outward rectifying conductances or A-like conductances in sympathetic preganglionic neurons (SPN) are prolonged, lasting for hundreds of milliseconds to seconds and are thought to play a key role in the regulation of SPN firing frequency. Here, a multidisciplinary electrophysiological, pharmacological and molecular single-cell rt-PCR approach was used to investigate the kinetics, pharmacological profile and putative K + channel subunits underlying the transient outward rectifying conductance expressed in SPN. SPN expressed a 4-aminopyridine (4-AP) sensitive transient outward rectification with significantly longer decay kinetics than reported for many other central neurons. The conductance and corresponding current in voltage-clamp conditions was also sensitive to the Kv4.2 and Kv4.3 blocker phrixotoxin-2 (1-10 µM) and the blocker of rapidly inactivating Kv channels, pandinotoxin-Ka (50 nM). The conductance and corresponding current was only weakly sensitive to the Kv1 channel blocker tityustoxin-Ka and insensitive to dendrotoxin I (200 nM) and the Kv3.4 channel blocker BDS-II (1 µM). Single-cell RT-PCR revealed mRNA expression for the a-subunits Kv4.1 and Kv4.3 in the majority and Kv1.5 in less than half of SPN. mRNA for accessory ß-subunits was detected for Kvß2 in all SPN with differential expression of mRNA for KChIP1, Kvß1 and Kvß3 and the peptidase homologue DPP6. These data together suggest that the transient outwardly rectifying conductance in SPN is mediated by members of the Kv4 subfamily (Kv4.1 and Kv4.3) in association with the ß-subunit Kvß2. Differential expression of the accessory ß subunits, which may act to modulate channel density and kinetics in SPN, may underlie the prolonged and variable time-course of this conductance in these neurons. © 2011 IBRO.

Tigecycline challenge triggers sRNA production in Salmonella enterica serovar Typhimurium

Background: Bacteria employ complex transcriptional networks involving multiple genes in response to stress, which is not limited to gene and protein networks but now includes small RNAs (sRNAs). These regulatory RNA molecules are increasingly shown to be able to initiate regulatory cascades and modulate the expression of multiple genes that are involved in or required for survival under environmental challenge. Despite mounting evidence for the importance of sRNAs in stress response, their role upon antibiotic exposure remains unknown. In this study, we sought to determine firstly, whether differential expression of sRNAs occurs upon antibiotic exposure and secondly, whether these sRNAs could be attributed to microbial tolerance to antibiotics.

Results: A small scale sRNA cloning strategy of Salmonella enterica serovar Typhimurium SL1344 challenged with half the minimal inhibitory concentration of tigecycline identified four sRNAs (sYJ5, sYJ20, sYJ75 and sYJ118) which were reproducibly upregulated in the presence of either tigecycline or tetracycline. The coding sequences of the four sRNAs were found to be conserved across a number of species. Genome analysis found that sYJ5 and sYJ118 mapped between the 16S and 23S rRNA encoding genes. sYJ20 (also known as SroA) is encoded upstream of the tbpAyabKyabJ operon and is classed as a riboswitch, whilst its role in antibiotic stress-response appears independent of its riboswitch function. sYJ75 is encoded between genes that are involved in enterobactin transport and metabolism. Additionally we find that the genetic deletion of sYJ20 rendered a reduced viability phenotype in the presence of tigecycline, which was recovered when complemented. The upregulation of some of these sRNAs were also observed when S. Typhimurium was challenged by ampicillin (sYJ5, 75 and 118); or when Klebsiella pneumoniae was challenged by tigecycline (sYJ20 and 118).

Conclusions: Small RNAs are overexpressed as a result of antibiotic exposure in S. Typhimurium where the same molecules are upregulated in a related species or after exposure to different antibiotics. sYJ20, a riboswitch, appears to possess a trans-regulatory sRNA role in antibiotic tolerance. These findings imply that the sRNA mediated response is a component of the bacterial response to antibiotic challenge.

New susceptibility Loci associated with kidney disease in type 1 diabetes

Diabetic kidney disease, or diabetic nephropathy (DN), is a major complication of diabetes and the leading cause of end-stage renal disease (ESRD) that requires dialysis treatment or kidney transplantation. In addition to the decrease in the quality of life, DN accounts for a large proportion of the excess mortality associated with type 1 diabetes (T1D). Whereas the degree of glycemia plays a pivotal role in DN, a subset of individuals with poorly controlled T1D do not develop DN. Furthermore, strong familial aggregation supports genetic susceptibility to DN. However, the genes and the molecular mechanisms behind the disease remain poorly understood, and current therapeutic strategies rarely result in reversal of DN. In the GEnetics of Nephropathy: an International Effort (GENIE) consortium, we have undertaken a meta-analysis of genome-wide association studies (GWAS) of T1D DN comprising ~2.4 million single nucleotide polymorphisms (SNPs) imputed in 6,691 individuals. After additional genotyping of 41 top ranked SNPs representing 24 independent signals in 5,873 individuals, combined meta-analysis revealed association of two SNPs with ESRD: rs7583877 in the AFF3 gene (P?=?1.2×10(-8)) and an intergenic SNP on chromosome 15q26 between the genes RGMA and MCTP2, rs12437854 (P?=?2.0×10(-9)). Functional data suggest that AFF3 influences renal tubule fibrosis via the transforming growth factor-beta (TGF-ß1) pathway. The strongest association with DN as a primary phenotype was seen for an intronic SNP in the ERBB4 gene (rs7588550, P?=?2.1×10(-7)), a gene with type 2 diabetes DN differential expression and in the same intron as a variant with cis-eQTL expression of ERBB4. All these detected associations represent new signals in the pathogenesis of DN.

Rice-arsenate interactions in hydroponics:whole genome transcriptional analysis

Rice (Oryza sativa) varieties that are arsenate-tolerant (Bala) and -sensitive (Azucena) were used to conduct a transcriptome analysis of the response of rice seedlings to sodium arsenate (AsV) in hydroponic solution. RNA extracted from the roots of three replicate experiments of plants grown for 1 week in phosphate-free nutrient with or without 13.3 muM AsV was used to challenge the Affymetrix (52K) GeneChip Rice Genome array. A total of 576 probe sets were significantly up-regulated at least 2-fold in both varieties, whereas 622 were down-regulated. Ontological classification is presented. As expected, a large number of transcription factors, stress proteins, and transporters demonstrated differential expression. Striking is the lack of response of classic oxidative stress-responsive genes or phytochelatin synthases/synthatases. However, the large number of responses from genes involved in glutathione synthesis, metabolism, and transport suggests that glutathione conjugation and arsenate methylation may be important biochemical responses to arsenate challenge. In this report, no attempt is made to dissect differences in the response of the tolerant and sensitive variety, but analysis in a companion article will link gene expression to the known tolerance loci available in the BalaxAzucena mapping population.

Elucidating the Regulon of the Multidrug Resistance Regulator RarA in Klebsiella pneumoniae

RarA is an AraC-type regulator in Klebsiella pneumoniae, which, when overexpressed, confers a low-level multidrug-resistant (MDR) phenotype linked to the upregulation of both the acrAB and oqxAB efflux genes. Increased rarA expression has also been shown to be integral in the development of tigecycline resistance in the absence of ramA in K. pneumoniae. Given its phenotypic role in MDR, microarray analyses were performed to determine the RarA regulon. Transcriptome analysis was undertaken using strains Ecl8?rarA/pACrarA-2 (rarA-expressing construct) and Ecl8?rarA/pACYC184 (vector-only control) using bespoke microarray slides consisting of probes derived from the genomic sequences of K. pneumoniae MGH 78578 (NC_009648.1) and Kp342 (NC_011283.1). Our results show that rarA overexpression resulted in the differential expression of 66 genes (42 upregulated and 24 downregulated). Under the COG (clusters of orthologous groups) functional classification, the majority of affected genes belonged to the category of cell envelope biogenesis and posttranslational modification, along with genes encoding the previously uncharacterized transport proteins (e.g., KPN_03141, sdaCB, and leuE) and the porin OmpF. However, genes associated with energy production and conversion and amino acid transport/metabolism (e.g., nuoA, narJ, and proWX) were found to be downregulated. Biolog phenotype analyses demonstrated that rarA overexpression confers enhanced growth of the overexpresser in the presence of several antibiotic classes (i.e., beta-lactams and fluoroquinolones), the antifungal/antiprotozoal compound clioquinol, disinfectants (8-hydroxyquinoline), protein synthesis inhibitors (i.e., minocycline and puromycin), membrane biogenesis agents (polymyxin B and amitriptyline), DNA synthesis (furaltadone), and the cytokinesis inhibitor (sanguinarine). Both our transcriptome and phenotypic microarray data support and extend the role of RarA in the MDR phenotype of K. pneumoniae.

Proteomic analysis of colorectal cancer metastasis:stathmin-1 revealed as a player in cancer cell migration and prognostic marker

Metastasis accounts largely for the high mortality rate of colorectal cancer (CRC) patients. In this study, we performed comparative proteome analysis of primary CRC cell lines HCT-116 and its metastatic derivative E1 using 2-D DIGE. We identified 74 differentially expressed proteins, many of which function in transcription, translation, angiogenesis signal transduction, or cytoskeletal remodeling pathways, which are indispensable cellular processes involved in the metastatic cascade. Among these proteins, stathmin-1 (STMN1) was found to be highly up-regulated in E1 as compared to HCT-116 and was thus selected for further functional studies. Our results showed that perturbations in STMN1 levels resulted in significant changes in cell migration, invasion, adhesion, and colony formation. We further showed that the differential expression of STMN1 correlated with the cells' metastatic potential in other paradigms of CRC models. Using immunohistochemistry, we also showed that STMN1 was highly expressed in colorectal primary tumors and metastatic tissues as compared to the adjacent normal colorectal tissues. Furthermore, we also showed via tissue microarray analyses of 324 CRC tissues and Kaplan-Meier survival plot that CRC patients with higher expression of STMN1 have poorer prognosis.

Re-analysis of microarray data reveals insights into altered transcriptional activity of TH17 and TReg signalling in psoriasis

Identifying differential expression of genes in psoriatic and healthy skin by microarray data analysis is a key approach to understand the pathogenesis of psoriasis. Analysis of more than one dataset to identify genes commonly upregulated reduces the likelihood of false positives and narrows down the possible signature genes. Genes controlling the critical balance between T helper 17 and regulatory T cells are of special interest in psoriasis. Our objectives were to identify genes that are consistently upregulated in lesional skin from three published microarray datasets. We carried out a reanalysis of gene expression data extracted from three experiments on samples from psoriatic and nonlesional skin using the same stringency threshold and software and further compared the expression levels of 92 genes related to the T helper 17 and regulatory T cell signaling pathways. We found 73 probe sets representing 57 genes commonly upregulated in lesional skin from all datasets. These included 26 probe sets representing 20 genes that have no previous link to the etiopathogenesis of psoriasis. These genes may represent novel therapeutic targets and surely need more rigorous experimental testing to be validated. Our analysis also identified 12 of 92 genes known to be related to the T helper 17 and regulatory T cell signaling pathways, and these were found to be differentially expressed in the lesional skin samples.

Multiple Regulatory Pathways Associated with High-Level Ciprofloxacin and Multidrug Resistance in Salmonella enterica Serovar Enteritidis: Involvement of ramA and Other Global Regulators

Mechanisms of antibiotic resistance were examined in nalidixic acid-resistant Salmonella enterica serovar Enteritidis field isolates displaying decreased susceptibility to ciprofloxacin and in in vitro-derived ciprofloxacin-resistant mutants (104-cip and 5408-cip). All field isolates harbored a single gyrA mutation (D87Y). Deletion of acrB and complementation with wild-type gyrA increased quinolone susceptibility. Selection for ciprofloxacin resistance was associated with the development of an additional gyrA (S83F) mutation in 104-cip, novel gyrB (E466D) and parE (V461G) mutations in 5408-cip, overexpression of acrB and decreased susceptibility to nonquinolone antibiotics in both mutants, and decreased OmpF production and altered lipopoly- saccharide in 104-cip. Complementation of mutated gyrA and gyrB with wild-type alleles restored susceptibility to quinolones in 104-cip and significantly decreased the ciprofloxacin MIC in 5408-cip. Complementation of parE had no effect on quinolone MICs. Deletion of acrB restored susceptibility to ciprofloxacin and other antibiotics tested. Both soxS and marA were overexpressed in 104-cip, and ramA was overexpressed in 5408-cip. Inactivation of each of these global regulators lowered ciprofloxacin MICs, decreased expression of acrB, and restored susceptibility to other antibiotics. Mutations were found in soxR (R20H) and in soxS (E52K) in 104-cip and in ramR (G25A) in 5408-cip. In conclusion, both efflux activity and a single gyrA mutation contribute to nalidixic acid resistance and reduced ciprofloxacin sensitivity. Ciprofloxacin resistance and decreased susceptibility to multiple antibiotics can result from different genetic events leading to development of target gene mutations, increased efflux activity resulting from differential expression of global regulators associated with mutations in their regulatory genes, and possible altered membrane permeability.

The Fasciola hepatica genome: gene duplication and polymorphism reveals adaptation to the host environment and the capacity for rapid evolution

BACKGROUND: The liver fluke Fasciola hepatica is a major pathogen of livestock worldwide, causing huge economic losses to agriculture, as well as 2.4 million human infections annually.

RESULTS: Here we provide a draft genome for F. hepatica, which we find to be among the largest known pathogen genomes at 1.3 Gb. This size cannot be explained by genome duplication or expansion of a single repeat element, and remains a paradox given the burden it may impose on egg production necessary to transmit infection. Despite the potential for inbreeding by facultative self-fertilisation, substantial levels of polymorphism were found, which highlights the evolutionary potential for rapid adaptation to changes in host availability, climate change or to drug or vaccine interventions. Non-synonymous polymorphisms were elevated in genes shared with parasitic taxa, which may be particularly relevant for the ability of the parasite to adapt to a broad range of definitive mammalian and intermediate molluscan hosts. Large-scale transcriptional changes, particularly within expanded protease and tubulin families, were found as the parasite migrated from the gut, across the peritoneum and through the liver to mature in the bile ducts. We identify novel members of anti-oxidant and detoxification pathways and defined their differential expression through infection, which may explain the stage-specific efficacy of different anthelmintic drugs.

CONCLUSIONS: The genome analysis described here provides new insights into the evolution of this important pathogen, its adaptation to the host environment and external selection pressures. This analysis also provides a platform for research into novel drugs and vaccines.

Molecular subtyping of primary prostate cancer reveals specific and shared target genes of different ETS rearrangements

This work aimed to evaluate whether ETS transcription factors frequently involved in rearrangements in prostate carcinomas (PCa), namely ERG and ETV1, regulate specific or shared target genes. We performed differential expression analysis on nine normal prostate tissues and 50 PCa enriched for different ETS rearrangements using exon-level expression microarrays, followed by in vitro validation using cell line models. We found specific deregulation of 57 genes in ERG-positive PCa and 15 genes in ETV1-positive PCa, whereas deregulation of 27 genes was shared in both tumor subtypes. We further showed that the expression of seven tumor-associated ERG target genes (PLA1A, CACNA1D, ATP8A2, HLA-DMB, PDE3B, TDRD1, and TMBIM1) and two tumor-associated ETV1 target genes (FKBP10 and GLYATL2) was significantly affected by specific ETS silencing in VCaP and LNCaP cell line models, respectively, whereas the expression of three candidate ERG and ETV1 shared targets (GRPR, KCNH8, and TMEM45B) was significantly affected by silencing of either ETS. Interestingly, we demonstrate that the expression of TDRD1, the topmost overexpressed gene of our list of ERG-specific candidate targets, is inversely correlated with the methylation levels of a CpG island found at -66 bp of the transcription start site in PCa and that TDRD1 expression is regulated by direct binding of ERG to the CpG island in VCaP cells. We conclude that ETS transcription factors regulate specific and shared target genes and that TDRD1, FKBP10, and GRPR are promising therapeutic targets and can serve as diagnostic markers for molecular subtypes of PCa harboring specific fusion gene rearrangements.

The androgen receptor fuels prostate cancer by regulating central metabolism and biosynthesis

The androgen receptor (AR) is a key regulator of prostate growth and the principal drug target for the treatment of prostate cancer. Previous studies have mapped AR targets and identified some candidates which may contribute to cancer progression, but did not characterize AR biology in an integrated manner. In this study, we took an interdisciplinary approach, integrating detailed genomic studies with metabolomic profiling and identify an anabolic transcriptional network involving AR as the core regulator. Restricting flux through anabolic pathways is an attractive approach to deprive tumours of the building blocks needed to sustain tumour growth. Therefore, we searched for targets of the AR that may contribute to these anabolic processes and could be amenable to therapeutic intervention by virtue of differential expression in prostate tumours. This highlighted calcium/calmodulin-dependent protein kinase kinase 2, which we show is overexpressed in prostate cancer and regulates cancer cell growth via its unexpected role as a hormone-dependent modulator of anabolic metabolism. In conclusion, it is possible to progress from transcriptional studies to a promising therapeutic target by taking an unbiased interdisciplinary approach.

AURKA overexpression accompanies dysregulation of DNA-damage response genes in invasive urothelial cell carcinoma

Invasive urothelial cell carcinoma (UCC) is characterized by increased chromosomal instability and follows an aggressive clinical course in contrast to non-invasive disease. To identify molecular processes that confer and maintain an aggressive malignant phenotype, we used a high-throughput genome-wide approach to interrogate a cohort of high and low clinical risk UCC tumors. Differential expression analyses highlighted cohesive dysregulation of critical genes involved in the G(2)/M checkpoint in aggressive UCC. Hierarchical clustering based on DNA Damage Response (DDR) genes separated tumors according to a pre-defined clinical risk phenotype. Using array-comparative genomic hybridization, we confirmed that the DDR was disrupted in tumors displaying high genomic instability. We identified DNA copy number gains at 20q13.2-q13.3 (AURKA locus) and determined that overexpression of AURKA accompanied dysregulation of DDR genes in high risk tumors. We postulated that DDR-deficient UCC tumors are advantaged by a selective pressure for AURKA associated override of M phase barriers and confirmed this in an independent tissue microarray series. This mechanism that enables cancer cells to maintain an aggressive phenotype forms a rationale for targeting AURKA as a therapeutic strategy in advanced stage UCC.

Genome-wide association study for sight-threatening diabetic retinopathy reveals association with genetic variation near the GRB2 gene

AIMS/HYPOTHESIS: Diabetic retinopathy is a serious complication of diabetes mellitus and can lead to blindness. A genetic component, in addition to traditional risk factors, has been well described although strong genetic factors have not yet been identified. Here, we aimed to identify novel genetic risk factors for sight-threatening diabetic retinopathy using a genome-wide association study.

METHODS: Retinopathy was assessed in white Australians with type 2 diabetes mellitus. Genome-wide association analysis was conducted for comparison of cases of sight-threatening diabetic retinopathy (n = 336) with diabetic controls with no retinopathy (n = 508). Top ranking single nucleotide polymorphisms were typed in a type 2 diabetes replication cohort, a type 1 diabetes cohort and an Indian type 2 cohort. A mouse model of proliferative retinopathy was used to assess differential expression of the nearby candidate gene GRB2 by immunohistochemistry and quantitative western blot.

RESULTS: The top ranked variant was rs3805931 with p = 2.66 × 10(-7), but no association was found in the replication cohort. Only rs9896052 (p = 6.55 × 10(-5)) was associated with sight-threatening diabetic retinopathy in both the type 2 (p = 0.035) and the type 1 (p = 0.041) replication cohorts, as well as in the Indian cohort (p = 0.016). The study-wide meta-analysis reached genome-wide significance (p = 4.15 × 10(-8)). The GRB2 gene is located downstream of this variant and a mouse model of retinopathy showed increased GRB2 expression in the retina.

CONCLUSIONS/INTERPRETATION: Genetic variation near GRB2 on chromosome 17q25.1 is associated with sight-threatening diabetic retinopathy. Several genes in this region are promising candidates and in particular GRB2 is upregulated during retinal stress and neovascularisation.