GSVA: gene set variation analysis for microarray and RNA-seq data

Gene set enrichment (GSE) analysis is a popular framework for condensing information from gene expression profiles into a pathway or signature summary. The strengths of this approach over single gene analysis include noise and dimension reduction, as well as greater biological interpretability. As molecular profiling experiments move beyond simple case-control studies, robust and flexible GSE methodologies are needed that can model pathway activity within highly heterogeneous data sets. To address this challenge, we introduce Gene Set Variation Analysis (GSVA), a GSE method that estimates variation of pathway activity over a sample population in an unsupervised manner. We demonstrate the robustness of GSVA in a comparison with current state of the art sample-wise enrichment methods. Further, we provide examples of its utility in differential pathway activity and survival analysis. Lastly, we show how GSVA works analogously with data from both microarray and RNA-seq experiments. GSVA provides increased power to detect subtle pathway activity changes over a sample population in comparison to corresponding methods. While GSE methods are generally regarded as end points of a bioinformatic analysis, GSVA constitutes a starting point to build pathway-centric models of biology. Moreover, GSVA contributes to the current need of GSE methods for RNA-seq data. GSVA is an open source software package for R which forms part of the Bioconductor project and can be downloaded at http://www.bioconductor.org.

Exploring the effect of aminoglycoside guanidinylation on ligands for Tau exon 10 splicing regulatory element RNA

We describe the effect of guanidinylation of the aminoglycoside moiety on acridine-neamine-containing ligands for the stem-loop structure located at the exon 10-5′-intron junction of Tau pre-mRNA, an important regulatory element of tau gene alternative splicing. On the basis of dynamic combinatorial chemistry experiments, ligands that combine guanidinoneamine and two different acridines were synthesized and their RNA-binding properties were compared with those of their amino precursors. Fluorescence titration experiments and UV-monitored melting curves revealed that guanidinylation has a positive effect both on the binding affinity and specificity of the ligands for the stemloop RNA, as well as on the stabilization of all RNA sequences evaluated, particularly some mutated sequences associated with the development of FTDP-17 tauopathy. However, this correlation between binding affinity and stabilization due to guanidinylation was only found in ligands containing a longer spacer between the acridine and guanidinoneamine moieties, since a shorter spacer produced the opposite effect (e.g. lower binding affinity and lower stabilization). Furthermore, spectroscopic studies suggest that ligand binding does not significantly change the overall RNA structure upon binding (circular dichroism) and that the acridine moiety might intercalate near the bulged region of the stem->loop structure (UV-Vis and NMR spectroscopy).

Ip-10 Is Associated With Il28b Variation And Predicts The First Phase Decline Of Hcv Rna And Outcome Of Therapy In Chronic Hepatitis C

High systemic levels of IP-10 at onset of combination therapy for chronic hepatitis C mirror intrahepatic mRNA levels and predict a slower first phase decline in HCV RNA as well as poor outcome. Recently several genome wide association studies have revealed that single nucleotide polymorphisms (SNPs) on chromosome19 within proximity of IL28B predict spontaneous clearance of HCV infection and as therapeutic outcome among patients infected with HCV genotype 1, with three such SNPs being highly predictive: rs12979860, rs12980275, and rs8099917. In the present study, we correlated genetic variations in these SNPs from 253 Caucasian patients with pretreatment plasma levels of IP-10 and HCV RNA throughout therapy within a phase III treatment trial (HCV-DITTO). The favorable genetic variations in all three SNPs (CC, AA, and TT respectively) was significantly associated with lower baseline IP-10 (CC vs. CT/TT at rs12979860: median 189 vs. 258 pg/mL, P=0.02, AA vs. AG/GG at rs12980275: median 189 vs. 258 pg/mL, P=0.01, TT vs. TG/GG at rs8099917: median 224 vs. 288 pg/mL, P=0.04), were significantly less common among HCV genotype 1 infected patients than genotype 2/3 (P<0.0001, P<0.0001, and P=0.01 respectively) and had significantly higher baseline viral load than carriers of the SNP genotypes (6.3 vs. 5.9 log 10 IU/mL, P=0.0012, 6.3 vs. 6.0 log 10 IU/mL, P=0.026, and 6.3 vs. 5.8 log 10 IU/mL, P=0.0003 respectively). Among HCV genotype 1 infected homozygous or heterogeneous carriers of the favorable C, A, and T genotypes, lower baseline IP-10 was significantly associated with greater decline in HCV-RNA day 0-4, which translated into increased rates of achieving SVR among homozygous patients with baseline IP-10 below 150 pg/mL (85%, 75%, and 75% respectively). In a multivariate analysis among genotype 1 infected patients, both baseline IP-10 and the SNPs were significant independent predictors of SVR. Conclusion: Baseline plasma IP-10 is significantly associated with IL28B variations, and augments the predictiveness of the first phase decline in HCV RNA and final treatment outcome.

Reduction of choroidal neovascularization in mice by adeno-associated virus-delivered anti-vascular endothelial growth factor short hairpin RNA.

BACKGROUND: Strategies leading to the long-term suppression of inappropriate ocular angiogenesis are required to avoid the need for repetitive monthly injections for treatment of diseases of the eye, such as age-related macular degeneration (AMD). The present study aimed to develop a strategy for the sustained repression of vascular endothelial growth factor (VEGF), which is identified as the key player in exudative AMD. METHODS: We have employed short hairpin (sh)RNAs combined with adeno-associated virus (AAV) delivery to obtain the targeted expression of potent gene-regulatory molecules. Anti-VEGF shRNAs were analyzed in human retinal pigment epithelial (RPE) cells using Renilla luciferase screening. For in vivo delivery of the most potent shRNA, self-complementary AAV vectors were packaged in serotype 8 capsids (scAAV2/8-hU6-sh9). In vivo efficacy was evaluated either by injection of scAAV2/8-hU6-sh9 into murine hind limb muscles or in a laser-induced murine model of choroidal neovascularization (CNV) following scAAV2/8-hU6-sh9 subretinal delivery. RESULTS: Plasmids encoding anti-VEGF shRNAs showed efficient knockdown of human VEGF in RPEs. Intramuscular administration led to localized expression and 91% knockdown of endogenous murine (m)VEGF. Subsequently, the ability of AAV2/8-encoded shRNAs to impair vessel formation was evaluated in the murine model of CNV. In this model, the sizes of the CNV were significantly reduced (up to 48%) following scAAV2/8-hU6-sh9 subretinal delivery. CONCLUSIONS: Using anti-VEGF vectors, we have demonstrated efficient silencing of endogenous mVEGF and showed that subretinal administration of scAAV2/8-hU6-sh9 has the ability to impair vessel formation in an AMD animal model. Thus, AAV-encoded shRNA can be used for the inhibition of neovascularization, leading to the development of sustained anti-VEGF therapy. Copyright © 2012 John Wiley & Sons, Ltd.

Structure-function analysis of the human TFIIB-related factor II protein reveals an essential role for the C-terminal domain in RNA polymerase III transcription.

The transcription factors TFIIB, Brf1, and Brf2 share related N-terminal zinc ribbon and core domains. TFIIB bridges RNA polymerase II (Pol II) with the promoter-bound preinitiation complex, whereas Brf1 and Brf2 are involved, as part of activities also containing TBP and Bdp1 and referred to here as Brf1-TFIIIB and Brf2-TFIIIB, in the recruitment of Pol III. Brf1-TFIIIB recruits Pol III to type 1 and 2 promoters and Brf2-TFIIIB to type 3 promoters such as the human U6 promoter. Brf1 and Brf2 both have a C-terminal extension absent in TFIIB, but their C-terminal extensions are unrelated. In yeast Brf1, the C-terminal extension interacts with the TBP/TATA box complex and contributes to the recruitment of Bdp1. Here we have tested truncated Brf2, as well as Brf2/TFIIB chimeric proteins for U6 transcription and for assembly of U6 preinitiation complexes. Our results characterize functions of various human Brf2 domains and reveal that the C-terminal domain is required for efficient association of the protein with U6 promoter-bound TBP and SNAP(c), a type 3 promoter-specific transcription factor, and for efficient recruitment of Bdp1. This in turn suggests that the C-terminal extensions in Brf1 and Brf2 are crucial to specific recruitment of Pol III over Pol II.

RNA interference: a new and powerful tool for functional genomic analysis.

In many species, the introduction of double-stranded RNA induces potent and specific gene silencing, referred to as RNA interference. This phenomenon, which is based on targeted degradation of mRNAs and occurs in almost any eukaryote, from trypanosomes to mice including plants and fungi, has sparked general interest from both applied and fundamental standpoints. RNA interference, which is currently used to investigate gene function in a variety of systems, is linked to natural resistance to viruses and transposon silencing, as if it were a primitive immune system involved in genome surveillance. Here, we review the mechanism of RNA interference in post-transcriptional gene silencing, its function in nature, its value for functional genomic analysis, and the modifications and improvements that may make it more efficient and inheritable. We also discuss the future directions of this versatile technique in both fundamental and applied science.

RNA interference mitigates motor and neuropathological deficits in a cerebellar mouse model of machado-joseph disease.

Machado-Joseph disease or Spinocerebellar ataxia type 3 is a progressive fatal neurodegenerative disorder caused by the polyglutamine-expanded protein ataxin-3. Recent studies demonstrate that RNA interference is a promising approach for the treatment of Machado-Joseph disease. However, whether gene silencing at an early time-point is able to prevent the appearance of motor behavior deficits typical of the disease when initiated before onset of the disease had not been explored. Here, using a lentiviral-mediated allele-specific silencing of mutant ataxin-3 in an early pre-symptomatic cerebellar mouse model of Machado-Joseph disease we show that this strategy hampers the development of the motor and neuropathological phenotypic characteristics of the disease. At the histological level, the RNA-specific silencing of mutant ataxin-3 decreased formation of mutant ataxin-3 aggregates, preserved Purkinje cell morphology and expression of neuronal markers while reducing cell death. Importantly, gene silencing prevented the development of impairments in balance, motor coordination, gait and hyperactivity observed in control mice. These data support the therapeutic potential of RNA interference for Machado-Joseph disease and constitute a proof of principle of the beneficial effects of early allele-specific silencing for therapy of this disease.

Agents anti-infectieux et fonction réna!e: vers des posologies sur mesure [Antimicrobial agents and renal elimination: towards individual dosage adjustment?].

The efficacy and safety of anti-infective treatments are associated with the drug blood concentration profile, which is directly correlated with a dosing adjustment to the individual patient's condition. Dosing adjustments to the renal function recommended in reference books are often imprecise and infrequently applied in clinical practice. The recent generalisation of the KDOQI (Kidney Disease Outcome Quality Initiative) staging of chronically impaired renal function represents an opportunity to review and refine the dosing recommendations in patients with renal insufficiency. The literature has been reviewed and compared to a predictive model of the fraction of drug cleared by the kidney based on the Dettli's principle. Revised drug dosing recommendations integrating these predictive parameters are proposed.

Combining RT-PCR-seq and RNA-seq to catalog all genic elements encoded in the human genome.

Within the ENCODE Consortium, GENCODE aimed to accurately annotate all protein-coding genes, pseudogenes, and noncoding transcribed loci in the human genome through manual curation and computational methods. Annotated transcript structures were assessed, and less well-supported loci were systematically, experimentally validated. Predicted exon-exon junctions were evaluated by RT-PCR amplification followed by highly multiplexed sequencing readout, a method we called RT-PCR-seq. Seventy-nine percent of all assessed junctions are confirmed by this evaluation procedure, demonstrating the high quality of the GENCODE gene set. RT-PCR-seq was also efficient to screen gene models predicted using the Human Body Map (HBM) RNA-seq data. We validated 73% of these predictions, thus confirming 1168 novel genes, mostly noncoding, which will further complement the GENCODE annotation. Our novel experimental validation pipeline is extremely sensitive, far more than unbiased transcriptome profiling through RNA sequencing, which is becoming the norm. For example, exon-exon junctions unique to GENCODE annotated transcripts are five times more likely to be corroborated with our targeted approach than with extensive large human transcriptome profiling. Data sets such as the HBM and ENCODE RNA-seq data fail sampling of low-expressed transcripts. Our RT-PCR-seq targeted approach also has the advantage of identifying novel exons of known genes, as we discovered unannotated exons in ~11% of assessed introns. We thus estimate that at least 18% of known loci have yet-unannotated exons. Our work demonstrates that the cataloging of all of the genic elements encoded in the human genome will necessitate a coordinated effort between unbiased and targeted approaches, like RNA-seq and RT-PCR-seq.

TLR3 and Rig-like receptor on myeloid dendritic cells and Rig-like receptor on human NK cells are both mandatory for production of IFN-gamma in response to double-stranded RNA.

Cross-talk between NK cells and dendritic cells (DCs) is critical for the potent therapeutic response to dsRNA, but the receptors involved remained controversial. We show in this paper that two dsRNAs, polyadenylic-polyuridylic acid and polyinosinic-polycytidylic acid [poly(I:C)], similarly engaged human TLR3, whereas only poly(I:C) triggered human RIG-I and MDA5. Both dsRNA enhanced NK cell activation within PBMCs but only poly(I:C) induced IFN-gamma. Although myeloid DCs (mDCs) were required for NK cell activation, induction of cytolytic potential and IFN-gamma production did not require contact with mDCs but was dependent on type I IFN and IL-12, respectively. Poly(I:C) but not polyadenylic-polyuridylic acid synergized with mDC-derived IL-12 for IFN-gamma production by acting directly on NK cells. Finally, the requirement of both TLR3 and Rig-like receptor (RLR) on mDCs and RLRs but not TLR3 on NK cells for IFN-gamma production was demonstrated using TLR3- and Cardif-deficient mice and human RIG-I-specific activator. Thus, we report the requirement of cotriggering TLR3 and RLR on mDCs and RLRs on NK cells for a pathogen product to induce potent innate cell activation.

Detection of Leishmania RNA virus in Leishmania parasites.

BACKGROUND: Patients suffering from cutaneous leishmaniasis (CL) caused by New World Leishmania (Viannia) species are at high risk of developing mucosal (ML) or disseminated cutaneous leishmaniasis (DCL). After the formation of a primary skin lesion at the site of the bite by a Leishmania-infected sand fly, the infection can disseminate to form secondary lesions. This metastatic phenotype causes significant morbidity and is often associated with a hyper-inflammatory immune response leading to the destruction of nasopharyngeal tissues in ML, and appearance of nodules or numerous ulcerated skin lesions in DCL. Recently, we connected this aggressive phenotype to the presence of Leishmania RNA virus (LRV) in strains of L. guyanensis, showing that LRV is responsible for elevated parasitaemia, destructive hyper-inflammation and an overall exacerbation of the disease. Further studies of this relationship and the distribution of LRVs in other Leishmania strains and species would benefit from improved methods of viral detection and quantitation, especially ones not dependent on prior knowledge of the viral sequence as LRVs show significant evolutionary divergence. METHODOLOGY/PRINCIPAL FINDINGS: This study reports various techniques, among which, the use of an anti-dsRNA monoclonal antibody (J2) stands out for its specific and quantitative recognition of dsRNA in a sequence-independent fashion. Applications of J2 include immunofluorescence, ELISA and dot blot: techniques complementing an arsenal of other detection tools, such as nucleic acid purification and quantitative real-time-PCR. We evaluate each method as well as demonstrate a successful LRV detection by the J2 antibody in several parasite strains, a freshly isolated patient sample and lesion biopsies of infected mice. CONCLUSIONS/SIGNIFICANCE: We propose that refinements of these methods could be transferred to the field for use as a diagnostic tool in detecting the presence of LRV, and potentially assessing the LRV-related risk of complications in cutaneous leishmaniasis.

Heavy metal ions are potent inhibitors of protein folding.

Environmental and occupational exposure to heavy metals such as cadmium, mercury and lead results in severe health hazards including prenatal and developmental defects. The deleterious effects of heavy metal ions have hitherto been attributed to their interactions with specific, particularly susceptible native proteins. Here, we report an as yet undescribed mode of heavy metal toxicity. Cd2+, Hg2+ and Pb2+ proved to inhibit very efficiently the spontaneous refolding of chemically denatured proteins by forming high-affinity multidentate complexes with thiol and other functional groups (IC(50) in the nanomolar range). With similar efficacy, the heavy metal ions inhibited the chaperone-assisted refolding of chemically denatured and heat-denatured proteins. Thus, the toxic effects of heavy metal ions may result as well from their interaction with the more readily accessible functional groups of proteins in nascent and other non-native form. The toxic scope of heavy metals seems to be substantially larger than assumed so far.

Hepatitis C virus RNA replication requires a conserved structural motif within the transmembrane domain of the NS5B RNA-dependent RNA polymerase.

Hepatitis C virus (HCV) nonstructural protein 5B (NS5B), the viral RNA-dependent RNA polymerase (RdRp), is a tail-anchored protein with a highly conserved C-terminal transmembrane domain (TMD) that is required for the assembly of a functional replication complex. Here, we report that the TMD of the HCV RdRp can be functionally replaced by a newly identified analogous membrane anchor of the GB virus B (GBV-B) NS5B RdRp. Replicons with a chimeric RdRp consisting of the HCV catalytic domain and the GBV-B membrane anchor replicated with reduced efficiency. Compensatory amino acid changes at defined positions within the TMD improved the replication efficiency of these chimeras. These observations highlight a conserved structural motif within the TMD of the HCV NS5B RdRp that is required for RNA replication.

Aminoterminal amphipathic α-helix AH1 of hepatitis C virus nonstructural protein 4B possesses a dual role in RNA replication and virus production.

Nonstructural protein 4B (NS4B) is a key organizer of hepatitis C virus (HCV) replication complex formation. In concert with other nonstructural proteins, it induces a specific membrane rearrangement, designated as membranous web, which serves as a scaffold for the HCV replicase. The N-terminal part of NS4B comprises a predicted and a structurally resolved amphipathic α-helix, designated as AH1 and AH2, respectively. Here, we report a detailed structure-function analysis of NS4B AH1. Circular dichroism and nuclear magnetic resonance structural analyses revealed that AH1 folds into an amphipathic α-helix extending from NS4B amino acid 4 to 32, with positively charged residues flanking the helix. These residues are conserved among hepaciviruses. Mutagenesis and selection of pseudorevertants revealed an important role of these residues in RNA replication by affecting the biogenesis of double-membrane vesicles making up the membranous web. Moreover, alanine substitution of conserved acidic residues on the hydrophilic side of the helix reduced infectivity without significantly affecting RNA replication, indicating that AH1 is also involved in virus production. Selective membrane permeabilization and immunofluorescence microscopy analyses of a functional replicon harboring an epitope tag between NS4B AH1 and AH2 revealed a dual membrane topology of the N-terminal part of NS4B during HCV RNA replication. Luminal translocation was unaffected by the mutations introduced into AH1, but was abrogated by mutations introduced into AH2. In conclusion, our study reports the three-dimensional structure of AH1 from HCV NS4B, and highlights the importance of positively charged amino acid residues flanking this amphipathic α-helix in membranous web formation and RNA replication. In addition, we demonstrate that AH1 possesses a dual role in RNA replication and virus production, potentially governed by different topologies of the N-terminal part of NS4B.

The complementary strand of the human T-cell leukemia virus type 1 RNA genome encodes a bZIP transcription factor that down-regulates viral transcription.

The RNA genome of the human T-cell leukemia virus type 1 (HTLV-1) codes for proteins involved in infectivity, replication, and transformation. We report in this study the characterization of a novel viral protein encoded by the complementary strand of the HTLV-1 RNA genome. This protein, designated HBZ (for HTLV-1 bZIP factor), contains a N-terminal transcriptional activation domain and a leucine zipper motif in its C terminus. We show here that HBZ is able to interact with the bZIP transcription factor CREB-2 (also called ATF-4), known to activate the HTLV-1 transcription by recruiting the viral trans-activator Tax on the Tax-responsive elements (TxREs). However, we demonstrate that the HBZ/CREB-2 heterodimers are no more able to bind to the TxRE and cyclic AMP response element sites. Taking these findings together, the functional inactivation of CREB-2 by HBZ is suggested to contribute to regulation of the HTLV-1 transcription. Moreover, the characterization of a minus-strand gene protein encoded by HTLV-1 has never been reported until now.