A modular approach for integrative analysis of large-scale gene-expression and drug-response data

High-throughput technologies are now used to generate more than one type of data from the same biological samples. To properly integrate such data, we propose using co-modules, which describe coherent patterns across paired data sets, and conceive several modular methods for their identification. We first test these methods using in silico data, demonstrating that the integrative scheme of our Ping-Pong Algorithm uncovers drug-gene associations more accurately when considering noisy or complex data. Second, we provide an extensive comparative study using the gene-expression and drug-response data from the NCI-60 cell lines. Using information from the DrugBank and the Connectivity Map databases we show that the Ping-Pong Algorithm predicts drug-gene associations significantly better than other methods. Co-modules provide insights into possible mechanisms of action for a wide range of drugs and suggest new targets for therapy

An mRNA region of the canine distemper virus fusion protein gene lacking AUG codons can promote protein expression.

Canine distemper virus (CDV) produces a glycosylated type I fusion protein (F) with an internal hydrophobic signal sequence beginning around 115 residues downstream of the first AUG used for translation initiation. Cleavage of the signal sequence yields the F0 molecule, which is cleaved into the F1 and F2 subunits. Surprisingly, when all in-frame AUGs located in the first third of the F gene were mutated a protein of the same molecular size as the F0 molecule was still expressed from both the Onderstepoort (OP) and A75/17-CDV F genes. We designated this protein, which is initiated from a non-AUG codon protein Fx. Site-directed mutagenesis allowed to identify codon 85, a GCC codon coding for alanine, as the most likely position from which translation initiation of Fx occurs in OP-CDV. Deletion analysis demonstrated that at least 60 nucleotides upstream of the GCC codon are required for efficient Fx translation. This sequence is GC-rich, suggesting extensive folding. Secondary structure may therefore be important for translation initiation at codon 85.

Optimization of the piggyBac Transposon Using mRNA and Insulators: Toward a More Reliable Gene Delivery System.

Integrating and expressing stably a transgene into the cellular genome remain major challenges for gene-based therapies and for bioproduction purposes. While transposon vectors mediate efficient transgene integration, expression may be limited by epigenetic silencing, and persistent transposase expression may mediate multiple transposition cycles. Here, we evaluated the delivery of the piggyBac transposase messenger RNA combined with genetically insulated transposons to isolate the transgene from neighboring regulatory elements and stabilize expression. A comparison of piggyBac transposase expression from messenger RNA and DNA vectors was carried out in terms of expression levels, transposition efficiency, transgene expression and genotoxic effects, in order to calibrate and secure the transposition-based delivery system. Messenger RNA reduced the persistence of the transposase to a narrow window, thus decreasing side effects such as superfluous genomic DNA cleavage. Both the CTF/NF1 and the D4Z4 insulators were found to mediate more efficient expression from a few transposition events. We conclude that the use of engineered piggyBac transposase mRNA and insulated transposons offer promising ways of improving the quality of the integration process and sustaining the expression of transposon vectors.

Improvement of the HIV/AIDS vaccine candidate NYVAC-C by deletion of the viral type I IFN inhibitor and/or acquisition of replication competence

Background: Recombinant viruses based on the attenuated vaccinia virus strain NYVAC are promising HIV vaccine candidates as phase I/II clinical trials have shown good safety and immunogenicity profiles. However, this NYVAC strain is non-replicating in most human cell lines and encodes viral inhibitors of the immune system. Methods: With the aim to increase the immune potency of the current NYVAC-C vector (expressing the codon optimized clade C HIV-1 genes encoding gp120 and Gag-Pol-Nef polyprotein), we have generated and characterized three NYVAC-C-based vectors by, 1) deletion of the viral type I IFN inhibitor gene (NYVAC-CdeltaB19R), 2) restoration of virus replication competence in human cells by re-inserting K1L and C7L host range genes (NYVAC-C-KC) and, 3) combination of both strategies (NYVACC- KC-deltaB19R). Results: Insertion of the KC fragment restored the replication competence of the viruses in human cells (HeLa cells and primary dermal fibroblasts and keratinocytes), increased the expression of HIV antigens by more than 3-fold compared to the non-replicating homologs, inhibited apoptosis induced by the parental NYVAC-C and retained attenuation in a newborn mouse model. In adult mice, replication-competent viruses showed a limited capacity to replicate in tissues surrounding the inoculation site (ovaries and lymph nodes). After infection of keratinocytes, PBMCs and dendritic cells these viruses induced differential modulation in specific host cell signal transduction pathways, triggering genes important in immune modulation. Conclusion: We have developed improved NYVAC-C-based vectors with enhanced HIV-1 antigen expression, with the ability to replicate in cultured human cells and partially in some tissues, with an induced expression of cellular genes relevant to immune system activation, and which trigger IFN-dependent and independent signalling pathways, while maintaining a safety phenotype. These new vectors are promising new HIV vaccine candidates. These studies were performed within the Poxvirus Tcell Vaccine Discovery Consortium (PTVDC) which is part of the CAVD program.

Detection of Brazilian hantavirus by reverse transcription polymerase chain reaction amplification of N gene in patients with hantavirus cardiopulmonary syndrome

We report a nested reverse transcription-polymerase chain reaction (RT-PCR) assay for hantavirus using primers selected to match high homology regions of hantavirus genomes detected from the whole blood of hantavirus cardiopulmonary syndrome (HCPS) patients from Brazil, also including the N gene nucleotide sequence of Araraquara virus. Hantavirus genomes were detected in eight out of nine blood samples from the HCPS patients by RT-PCR (88.9% positivity) and in all 9 blood samples (100% positivity) by nested-PCR. The eight amplicons obtained by RT-PCR (P1, P3-P9), including one obtained by nested-PCR (P-2) and not obtained by RT-PCR, were sequenced and showed high homology (94.8% to 99.1%) with the N gene of Araraquara hantavirus. Although the serologic method ELISA is the most appropriate test for HCPS diagnosis, the use of nested RT-PCR for hantavirus in Brazil would contribute to the diagnosis of acute hantavirus disease detecting viral genomes in patient specimens as well as initial genomic characterization of circulating hantaviruses.

Genetic variability of hepatitis A virus strain HAF-203 isolated in Brazil and expression of the VP1 gene in Escherichia coli

The hepatitis A virus (HAV) HAF-203 strain was isolated from an acute case of HAV infection. The primary isolation of HAF-203 in Brazil and its adaptation to the FRhK-4 cell lineage allowed the production of large amounts of viral particles enabling molecular characterization of the first HAV isolate in Brazil. The aim of our study was to determine the nucleotide sequence of the HAF-203 strain genome, compare it to other HAV genomes and highlight its genetic variability. The complete nucleotide sequence of the HAF-203 strain (7472 nucleotides) was compared to those obtained earlier by others for other HAV isolates. These analyses revealed 19 HAF-specific nucleotide sequence differences with 10 amino acid substitutions. Most of the non-conservative changes were located at VP1, 2C, and 3D genes, but the 3B region was the most variable. The availability of HAF-203 complementary DNA was useful for the production of the recombinant VP1 protein, which is a major determinant of viral infectivity. This recombinant protein was shown by enzyme-linked immunoassay and blotting, to be immunogenic and resemble the native protein, therefore suggesting its value as a reagent for incorporation into diagnostic tests.

Hitchhiking Trypanosoma cruzi minicircle DNA affects gene expression in human host cells via LINE-1 retrotransposon

The horizontal transfer of Trypanosoma cruzi mitochondrial minicircle DNA to the genomes of naturally infected humans may play an important role in the pathogenesis of Chagas disease. Minicircle integrations within LINE-1 elements create the potential for foreign DNA mobility within the host genome via the machinery associated with this retrotransposon. Here we document integration of minicircle DNA fragments in clonal human macrophage cell lines and their mobilization over time. The movement of an integration event in a clonal transfected cell line was tracked at three months and three years post-infection. The minicircle sequence integrated into a LINE-1 retrotransposon; one such foreign fragment subsequently relocated to another genomic location in association with associated LINE-1 elements. The p15 locus was altered at three years as a direct effect of minicircle/LINE-1 acquisition, resulting in elimination of p15 mRNA. Here we show for the first time a molecular pathology stemming from mobilization of a kDNA/LINE-1 mutation. These genomic changes and detected transcript variations are consistent with our hypothesis that minicircle integration is a causal component of parasite-independent, autoimmune-driven lesions seen in the heart and other target tissues associated with Chagas disease.

Re-mapping the molecular features of the human immunodeficiency virus type 1 and human T-cell lymphotropic virus type 1 Brazilian sequences using a bioinformatics unit established in Salvador, Bahia, Brazil, to give support to the viral epidemiology studies

The analysis of genetic data for human immunodeficiency virus type 1 (HIV-1) and human T-cell lymphotropic virus type 1 (HTLV-1) is essential to improve treatment and public health strategies as well as to select strains for vaccine programs. However, the analysis of large quantities of genetic data requires collaborative efforts in bioinformatics, computer biology, molecular biology, evolution, and medical science. The objective of this study was to review and improve the molecular epidemiology of HIV-1 and HTLV-1 viruses isolated in Brazil using bioinformatic tools available in the Laboratório Avançado de Sáude Pública (Lasp) bioinformatics unit. The analysis of HIV-1 isolates confirmed a heterogeneous distribution of the viral genotypes circulating in the country. The Brazilian HIV-1 epidemic is characterized by the presence of multiple subtypes (B, F1, C) and B/F1 recombinant virus while, on the other hand, most of the HTLV-1 sequences were classified as Transcontinental subgroup of the Cosmopolitan subtype. Despite the high variation among HIV-1 subtypes, protein glycosylation and phosphorylation domains were conserved in the pol, gag, and env genes of the Brazilian HIV-1 strains suggesting constraints in the HIV-1 evolution process. As expected, the functional protein sites were highly conservative in the HTLV-1 env gene sequences. Furthermore, the presence of these functional sites in HIV-1 and HTLV-1 strains could help in the development of vaccines that pre-empt the viral escape process.

The role of bile acid retention and a common polymorphism in the ABCB11 gene as host factors affecting antiviral treatment response in chronic hepatitis C.

Summary.  The outcome of hepatitis C virus (HCV) infection and the likelihood of a sustained virological response (SVR) to antiviral therapy depends on both viral and host characteristics. In vitro studies demonstrated that bile acids (BA) interfere with antiviral interferon effects. We investigate the influence of plasma BA concentrations and an ABCB11 polymorphism associated with lower transporter expression on viral load and SVR. Four hundred and fifty-one Caucasian HCV-patients treated with PEG-interferon and ribavirin were included in the study. ABCB11 1331T>C was genotyped, and plasma BA levels were determined. The 1331C allele was slightly overrepresented in HCV-patients compared to controls. In HCV-patients, a significant difference between patients achieving SVR vs non-SVR was observed for HCV-2/3 (5 vs 9 μm; P = 0.0001), while median BA levels in HCV-1 were marginally elevated. Normal BA levels <8 μm were significantly associated with SVR (58.3%vs 36.3%; OR 2.48; P = 0.0001). This difference was significant for HCV-2/3 (90.7%vs 67.6%; P = 0.002) but marginal in HCV-1 (38.7%vs 27.8%; P = 0.058). SVR rates were equivalent between ABCB11 genotypes for HCV-1, but increased for HCV-2/3 (TT 100%vs CC 78%; OR 2.01; P = 0.043). IL28B genotype had no influence on these associations. No correlation between BA levels and HCV RNA was detected for any HCV genotype. The higher allelic frequency of ABCB11 1331C in HCV-patients compared to controls may indirectly link increased BA to HCV chronicity. Our data support a role for BA as host factor affecting therapy response in HCV-2/3 patients, whereas a weaker association was found for HCV-1.

Characterization of polymorphisms in the mannose-binding lectin gene promoter among human immunodeficiency virus 1 infected subjects

The present study investigated the prevalence of mutations in the -550 (H/L) and -221 (X/Y) mannose-binding lectin (MBL) gene promoter regions and their impact on infection by human immunodeficiency virus 1 (HIV-1) in a population of 128 HIV-1 seropositive and 97 seronegative patients. The allele identification was performed through the sequence-specific primer polymerase chain reaction method, using primer sequences specific to each polymorphism. The evolution of the infection was evaluated through CD4+ T-lymphocyte counts and plasma viral load. The allele and haplotype frequencies among HIV-1-infected patients and seronegative healthy control patients did not show significant differences. CD4+ T-lymphocyte counts showed lower levels among seropositive patients carrying haplotypes LY, LX and HX, as compared to those carrying the HY haplotype. Mean plasma viral load was higher among seropositive patients with haplotypes LY, LX and HX than among those carrying the HY haplotype. When promoter and exon 1 mutations were matched, it was possible to identify a significantly higher viral load among HIV-1 infected individuals carrying haplotypes correlated to low serum levels of MBL. The current study shows that haplotypes related to medium and low MBL serum levels might directly influence the evolution of viral progression in patients. Therefore, it is suggested that the identification of haplotypes within the promoter region of the MBL gene among HIV-1 infected persons should be further evaluated as a prognostic tool for AIDS progression.

Methods for the ex vivo characterization of human CD8+ T subsets based on gene expression and replicative history analysis.

The generation of an antigen-specific T-lymphocyte response is a complex multi-step process. Upon T-cell receptor-mediated recognition of antigen presented by activated dendritic cells, naive T-lymphocytes enter a program of proliferation and differentiation, during the course of which they acquire effector functions and may ultimately become memory T-cells. A major goal of modern immunology is to precisely identify and characterize effector and memory T-cell subpopulations that may be most efficient in disease protection. Sensitive methods are required to address these questions in exceedingly low numbers of antigen-specific lymphocytes recovered from clinical samples, and not manipulated in vitro. We have developed new techniques to dissect immune responses against viral or tumor antigens. These allow the isolation of various subsets of antigen-specific T-cells (with major histocompatibility complex [MHC]-peptide multimers and five-color FACS sorting) and the monitoring of gene expression in individual cells (by five-cell reverse transcription-polymerase chain reaction [RT-PCR]). We can also follow their proliferative life history by flow-fluorescence in situ hybridization (FISH) analysis of average telomere length. Recently, using these tools, we have identified subpopulations of CD8+ T-lymphocytes with distinct proliferative history and partial effector-like properties. Our data suggest that these subsets descend from recently activated T-cells and are committed to become differentiated effector T-lymphocytes.

Sequence analysis of the 2009 pandemic influenza A H1N1 virus haemagglutinin gene from 2009-2010 Brazilian clinical samples

In this paper, we analysed the haemagglutinin (HA) gene identified by polymerase chain reaction from 90 influenza A H1N1 virus strains that circulated in Brazil from April 2009-June 2010. A World Health Organization sequencing protocol allowed us to identify amino acid mutations in the HA protein at positions S220T (71%), D239G/N/S (20%), Y247H (4.5%), E252K (3.3%), M274V (2.2%), Q310H (26.7%) and E391K (12%). A fatal outcome was associated with the D239G mutation (p < 0.0001). Brazilian HA genetic diversity, in comparison to a reference strain from California, highlights the role of influenza virus surveillance for study of viral evolution, in addition to monitoring the spread of the virus worldwide.

Strategies for the successful implementation of viral laboratory automation.

It has been estimated that more than 70% of all medical activity is directly related to information providing analytical data. Substantial technological advances have taken place recently, which have allowed a previously unimagined number of analytical samples to be processed while offering high quality results. Concurrently, yet more new diagnostic determinations have been introduced - all of which has led to a significant increase in the prescription of analytical parameters. This increased workload has placed great pressure on the laboratory with respect to health costs. The present manager of the Clinical Laboratory (CL) has had to examine cost control as well as rationing - meaning that the CL's focus has not been strictly metrological, as if it were purely a system producing results, but instead has had to concentrate on its efficiency and efficacy. By applying re-engineering criteria, an emphasis has had to be placed on improved organisation and operating practice within the CL, focussing on the current criteria of the Integrated Management Areas where the technical and human resources are brought together. This re-engineering has been based on the concepts of consolidating and integrating the analytical platforms, while differentiating the production areas (CORE Laboratory) from the information areas. With these present concepts in mind, automation and virological treatment, along with serology in general, follow the same criteria as the rest of the operating methodology in the Clinical Laboratory.

Equine herpesvirus-2 E10 gene product, but not its cellular homologue, activates NF-kappaB transcription factor and c-Jun N-terminal kinase.

We have previously reported on the death effector domain containing E8 gene product from equine herpesvirus-2, designated FLICE inhibitory protein (v-FLIP), and on its cellular homologue, c-FLIP, which inhibit the activation of caspase-8 by death receptors. Here we report on the structure and function of the E10 gene product of equine herpesvirus-2, designated v-CARMEN, and on its cellular homologue, c-CARMEN, which contain a caspase-recruiting domain (CARD) motif. c-CARMEN is highly homologous to the viral protein in its N-terminal CARD motif but differs in its C-terminal extension. v-CARMEN and c-CARMEN interact directly in a CARD-dependent manner yet reveal different binding specificities toward members of the tumor necrosis factor receptor-associated factor (TRAF) family. v-CARMEN binds to TRAF6 and weakly to TRAF3 and, upon overexpression, potently induces the c-Jun N-terminal kinase (JNK), p38, and nuclear factor (NF)-kappaB transcriptional pathways. c-CARMEN or truncated versions thereof do not appear to induce JNK and NF-kappaB activation by themselves, nor do they affect the JNK and NF-kappaB activating potential of v-CARMEN. Thus, in contrast to the cellular homologue, v-CARMEN may have additional properties in its unique C terminus that allow for an autonomous activator effect on NF-kappaB and JNK. Through activation of NF-kappaB, v-CARMEN may regulate the expression of the cellular and viral genes important for viral replication.

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.