Generation and analysis of a mouse intestinal metatranscriptome through Illumina based RNA-sequencing

With the advent of high through-put sequencing (HTS), the emerging science of metagenomics is transforming our understanding of the relationships of microbial communities with their environments. While metagenomics aims to catalogue the genes present in a sample through assessing which genes are actively expressed, metatranscriptomics can provide a mechanistic understanding of community inter-relationships. To achieve these goals, several challenges need to be addressed from sample preparation to sequence processing, statistical analysis and functional annotation. Here we use an inbred non-obese diabetic (NOD) mouse model in which germ-free animals were colonized with a defined mixture of eight commensal bacteria, to explore methods of RNA extraction and to develop a pipeline for the generation and analysis of metatranscriptomic data. Applying the Illumina HTS platform, we sequenced 12 NOD cecal samples prepared using multiple RNA-extraction protocols. The absence of a complete set of reference genomes necessitated a peptide-based search strategy. Up to 16% of sequence reads could be matched to a known bacterial gene. Phylogenetic analysis of the mapped ORFs revealed a distribution consistent with ribosomal RNA, the majority from Bacteroides or Clostridium species. To place these HTS data within a systems context, we mapped the relative abundance of corresponding Escherichia coli homologs onto metabolic and protein-protein interaction networks. These maps identified bacterial processes with components that were well-represented in the datasets. In summary this study highlights the potential of exploiting the economy of HTS platforms for metatranscriptomics.

RNA-Seq-based analysis of the physiologic cold shock-induced changes in Moraxella catarrhalis gene expression

BACKGROUND Moraxella catarrhalis, a major nasopharyngeal pathogen of the human respiratory tract, is exposed to rapid downshifts of environmental temperature when humans breathe cold air. The prevalence of pharyngeal colonization and respiratory tract infections caused by M. catarrhalis is greatest in winter. We investigated how M. catarrhalis uses the physiologic exposure to cold air to regulate pivotal survival systems that may contribute to M. catarrhalis virulence. RESULTS In this study we used the RNA-seq techniques to quantitatively catalogue the transcriptome of M. catarrhalis exposed to a 26 °C cold shock or to continuous growth at 37 °C. Validation of RNA-seq data using quantitative RT-PCR analysis demonstrated the RNA-seq results to be highly reliable. We observed that a 26 °C cold shock induces the expression of genes that in other bacteria have been related to virulence a strong induction was observed for genes involved in high affinity phosphate transport and iron acquisition, indicating that M. catarrhalis makes a better use of both phosphate and iron resources after exposure to cold shock. We detected the induction of genes involved in nitrogen metabolism, as well as several outer membrane proteins, including ompA, m35-like porin and multidrug efflux pump (acrAB) indicating that M. catarrhalis remodels its membrane components in response to downshift of temperature. Furthermore, we demonstrate that a 26 °C cold shock enhances the induction of genes encoding the type IV pili that are essential for natural transformation, and increases the genetic competence of M. catarrhalis, which may facilitate the rapid spread and acquisition of novel virulence-associated genes. CONCLUSION Cold shock at a physiologically relevant temperature of 26 °C induces in M. catarrhalis a complex of adaptive mechanisms that could convey novel pathogenic functions and may contribute to enhanced colonization and virulence.

Prevalence and predictors of kaposi sarcoma herpes virus seropositivity: a cross-sectional analysis of HIV-infected adults initiating ART in Johannesburg, South Africa

Background Kaposi sarcoma (KS) is the most common AIDS-defining tumour in HIV-infected individuals in Africa. Kaposi sarcoma herpes virus (KSHV) infection precedes development of KS. KSHV co-infection may be associated with worse outcomes in HIV disease and elevated KSHV viral load may be an early marker for advanced HIV disease among untreated patients. We examined the prevalence of KSHV among adults initiating antiretroviral therapy (ART) and compared immunological, demographic and clinical factors between patients seropositive and seronegative for KSHV. Results We analyzed cross-sectional data collected from 404 HIV-infected treatment-naïve adults initiating ART at the Themba Lethu Clinic, Johannesburg, South Africa between November 2008 and March 2009. Subjects were screened at ART initiation for antibodies to KSHV lytic K8.1 and latent Orf73 antigens. Seropositivity to KSHV was defined as positive to either lytic KSHV K8.1 or latent KSHV Orf73 antibodies. KSHV viremia was determined by quantitative PCR and CD3, 4 and 8 lymphocyte counts were determined with flow cytometry. Of the 404 participants, 193 (48%) tested positive for KSHV at ART initiation; with 76 (39%) reactive to lytic K8.1, 35 (18%) to latent Orf73 and 82 (42%) to both. One individual presented with clinical KS at ART initiation. The KSHV infected group was similar to those without KSHV in terms of age, race, gender, ethnicity, smoking and alcohol use. KSHV infected individuals presented with slightly higher median CD3 (817 vs. 726 cells/mm3) and CD4 (90 vs. 80 cells/mm3) counts than KSHV negative subjects. We found no associations between KSHV seropositivity and body mass index, tuberculosis status, WHO stage, HIV RNA levels, full blood count or liver function tests at initiation. Those with detectable KSHV viremia (n = 19), however, appeared to present with signs of more advanced HIV disease including anemia and WHO stage 3 or 4 defining conditions compared to those in whom the virus was undetectable. Conclusions We demonstrate a high prevalence of KSHV among HIV-infected adults initiating ART in a large urban public-sector HIV clinic. KSHV viremia but not KSHV seropositivity may be associated with markers of advanced HIV disease.

A sensitized RNA interference screen identifies a novel role for the PI3K p110γ isoform in medulloblastoma cell proliferation and chemoresistance

Medulloblastoma is the most common malignant brain tumor in children and is associated with a poor outcome. We were interested in gaining further insight into the potential of targeting the human kinome as a novel approach to sensitize medulloblastoma to chemotherapeutic agents. A library of small interfering RNA (siRNA) was used to downregulate the known human protein and lipid kinases in medulloblastoma cell lines. The analysis of cell proliferation, in the presence or absence of a low dose of cisplatin after siRNA transfection, identified new protein and lipid kinases involved in medulloblastoma chemoresistance. PLK1 (polo-like kinase 1) was identified as a kinase involved in proliferation in medulloblastoma cell lines. Moreover, a set of 6 genes comprising ATR, LYK5, MPP2, PIK3CG, PIK4CA, and WNK4 were identified as contributing to both cell proliferation and resistance to cisplatin treatment in medulloblastoma cells. An analysis of the expression of the 6 target genes in primary medulloblastoma tumor samples and cell lines revealed overexpression of LYK5 and PIK3CG. The results of the siRNA screen were validated by target inhibition with specific pharmacological inhibitors. A pharmacological inhibitor of p110γ (encoded by PIK3CG) impaired cell proliferation in medulloblastoma cell lines and sensitized the cells to cisplatin treatment. Together, our data show that the p110γ phosphoinositide 3-kinase isoform is a novel target for combinatorial therapies in medulloblastoma.

RNA interference screening identifies a novel role for autocrine fibroblast growth factor signaling in neuroblastoma chemoresistance

Chemotherapeutic drug resistance is one of the major causes for treatment failure in high-risk neuroblastoma (NB), the most common extra cranial solid tumor in children. Poor prognosis is typically associated with MYCN amplification. Here, we utilized a loss-of-function kinome-wide RNA interference screen to identify genes that cause cisplatin sensitization. We identified fibroblast growth factor receptor 2 (FGFR2) as an important determinant of cisplatin resistance. Pharmacological inhibition of FGFR2 confirmed the importance of this kinase in NB chemoresistance. Silencing of FGFR2 sensitized NB cells to cisplatin-induced apoptosis, which was regulated by the downregulation of the anti-apoptotic proteins BCL2 and BCLX(L). Mechanistically, FGFR2 was shown to activate protein kinase C-δ to induce BCL2 expression. FGFR2, as well as the ligand fibroblast growth factor-2, were consistently expressed in primary NB and NB cell lines, indicating the presence of an autocrine loop. Expression analysis revealed that FGFR2 correlates with MYCN amplification and with advanced stage disease, demonstrating the clinical relevance of FGFR2 in NB. These findings suggest a novel role for FGFR2 in chemoresistance and provide a rational to combine pharmacological inhibitors against FGFR2 with chemotherapeutic agents for the treatment of NB.Oncogene advance online publication, 1 October 2012; doi:10.1038/onc.2012.416.

Comparative human-horse sequence analysis of the CYP3A subfamily gene cluster

Cytochrome P450 enzymes (CYP450s) represent a superfamily of haem-thiolate proteins. CYP450s are most abundant in the liver, a major site of drug metabolism, and play key roles in the metabolism of a variety of substrates, including drugs and environmental contaminants. Interaction of two or more different drugs with the same enzyme can account for adverse effects and failure of therapy. Human CYP3A4 metabolizes about 50% of all known drugs, but little is known about the orthologous CYP450s in horses. We report here the genomic organization of the equine CYP3A gene cluster as well as a comparative analysis with the human CYP3A gene cluster. The equine CYP450 genes of the 3A family are located on ECA 13 between 6.97-7.53 Mb, in a region syntenic to HSA 7 99.05-99.35 Mb. Seven potential, closely linked equine CYP3A genes were found, in contrast to only four genes in the human genome. RNA was isolated from an equine liver sample, and the approximately 1.5-kb coding sequence of six CYP3A genes could be amplified by RT-PCR. Sequencing of the RT-PCR products revealed numerous hitherto unknown single nucleotide polymorphisms (SNPs) in these six CYP3A genes, and one 6-bp deletion compared to the reference sequence (EquCab2.0). The presence of the variants was confirmed in a sample of genomic DNA from the same horse. In conclusion, orthologous genes for the CYP3A family exist in horses, but their number differs from those of the human CYP3A gene family. CYP450 genes of the same family show high homology within and between mammalian species, but can be highly polymorphic.

Differences between RNA and DNA due to RNA editing in temporal lobe epilepsy

To investigate whether alterations in RNA editing (an enzymatic base-specific change to the RNA sequence during primary transcript formation from DNA) of neurotransmitter receptor genes and of transmembrane ion channel genes play a role in human temporal lobe epilepsy (TLE), this exploratory study analyzed 14 known cerebral editing sites in RNA extracted from the brain tissue of 41 patients who underwent surgery for mesial TLE, 23 with hippocampal sclerosis (MTLE+HS). Because intraoperatively sampled RNA cannot be obtained from healthy controls and the best feasible control is identically sampled RNA from patients with a clinically shorter history of epilepsy, the primary aim of the study was to assess the correlation between epilepsy duration and RNA editing in the homogenous group of MTLE+HS. At the functionally relevant I/V site of the voltage-gated potassium channel Kv1.1, an inverse correlation of RNA editing was found with epilepsy duration (r=-0.52, p=0.01) but not with patient age at surgery, suggesting a specific association with either the epileptic process itself or its antiepileptic medication history. No significant correlations were found between RNA editing and clinical parameters at other sites within glutamate receptor or serotonin 2C receptor gene transcripts. An "all-or-none" (≥95% or ≤5%) editing pattern at most or all sites was discovered in 2 patients. As a secondary part of the study, RNA editing was also analyzed as in the previous literature where up to now, few single editing sites were compared with differently obtained RNA from inhomogenous patient groups and autopsies, and by measuring editing changes in our mouse model. The present screening study is first to identify an editing site correlating with a clinical parameter, and to also provide an estimate of the possible effect size at other sites, which is a prerequisite for power analysis needed in planning future studies.

RNA degradation differentially affects quantitative mRNA measurements of endogenous reference genes in human placenta

It has been highlighted that RNA quality and appropriate reference gene selection is crucial for the interpretation of RT-qPCR results in human placental samples. In this context we investigated the effect of RNA degradation on the mRNA abundance of seven frequently used reference genes in 119 human placental samples. Combining RNA integrity measurements, RT-qPCR analysis and mathematical modeling we found major differences regarding the effect of RNA degradation on the measured expression levels between the different reference genes. Furthermore, we demonstrated that a modified RNA extraction method significantly improved RNA quality and consequently increased transcript levels of all reference genes.

HIV treatment response and prognosis in Europe and North America in the first decade of highly active antiretroviral therapy: a collaborative analysis

BACKGROUND: Highly active antiretroviral therapy (HAART) for the treatment of HIV infection was introduced a decade ago. We aimed to examine trends in the characteristics of patients starting HAART in Europe and North America, and their treatment response and short-term prognosis. METHODS: We analysed data from 22,217 treatment-naive HIV-1-infected adults who had started HAART and were followed up in one of 12 cohort studies. The probability of reaching 500 or less HIV-1 RNA copies per mL by 6 months, and the change in CD4 cell counts, were analysed for patients starting HAART in 1995-96, 1997, 1998, 1999, 2000, 2001, and 2002-03. The primary endpoints were the hazard ratios for AIDS and for death from all causes in the first year of HAART, which were estimated using Cox regression. RESULTS: The proportion of heterosexually infected patients increased from 20% in 1995-96 to 47% in 2002-03, and the proportion of women from 16% to 32%. The median CD4 cell count when starting HAART increased from 170 cells per muL in 1995-96 to 269 cells per muL in 1998 but then decreased to around 200 cells per muL. In 1995-96, 58% achieved HIV-1 RNA of 500 copies per mL or less by 6 months compared with 83% in 2002-03. Compared with 1998, adjusted hazard ratios for AIDS were 1.07 (95% CI 0.84-1.36) in 1995-96 and 1.35 (1.06-1.71) in 2002-03. Corresponding figures for death were 0.87 (0.56-1.36) and 0.96 (0.61-1.51). INTERPRETATION: Virological response after starting HAART improved over calendar years, but such improvement has not translated into a decrease in mortality.

Prognosis of HIV-1-infected patients up to 5 years after initiation of HAART: collaborative analysis of prospective studies

OBJECTIVE: To estimate the prognosis over 5 years of HIV-1-infected, treatment-naive patients starting HAART, taking into account the immunological and virological response to therapy. DESIGN: A collaborative analysis of data from 12 cohorts in Europe and North America on 20,379 adults who started HAART between 1995 and 2003. METHODS: Parametric survival models were used to predict the cumulative incidence at 5 years of a new AIDS-defining event or death, and death alone, first from the start of HAART and second from 6 months after the start of HAART. Data were analysed by intention-to-continue-treatment, ignoring treatment changes and interruptions. RESULTS: During 61 798 person-years of follow-up, 1005 patients died and an additional 1303 developed AIDS. A total of 10 046 (49%) patients started HAART either with a CD4 cell count of less than 200 cells/microl or with a diagnosis of AIDS. The 5-year risk of AIDS or death (death alone) from the start of HAART ranged from 5.6 to 77% (1.8-65%), depending on age, CD4 cell count, HIV-1-RNA level, clinical stage, and history of injection drug use. From 6 months the corresponding figures were 4.1-99% for AIDS or death and 1.3-96% for death alone. CONCLUSION: On the basis of data collected routinely in HIV care, prognostic models with high discriminatory power over 5 years were developed for patients starting HAART in industrialized countries. A risk calculator that produces estimates for progression rates at years 1 to 5 after starting HAART is available from www.art-cohort-collaboration.org.

The double-stranded RNA-activated protein kinase mediates radiation resistance in mouse embryo fibroblasts through nuclear factor kappaB and Akt activation

PURPOSE: Activation of the double-stranded RNA-activated protein kinase (PKR) leads to the induction of various pathways including the down-regulation of translation through phosphorylation of the eukaryotic translation initiation factor 2alpha (eIF-2alpha). There have been no reports to date about the role of PKR in radiation sensitivity. EXPERIMENTAL DESIGN: A clonogenic survival assay was used to investigate the sensitivity of PKR mouse embryo fibroblasts (MEF) to radiation therapy. 2-Aminopurine (2-AP), a chemical inhibitor of PKR, was used to inhibit PKR activation. Nuclear factor-kappaB (NF-kappaB) activation was assessed by electrophoretic mobility shift assay (EMSA). Expression of PKR and downstream targets was examined by Western blot analysis and immunofluorescence. RESULTS: Ionizing radiation leads to dose- and time-dependent increases in PKR expression and function that contributes to increased cellular radiation resistance as shown by clonogenic survival and terminal nucleotidyl transferase-mediated nick end labeling (TUNEL) apoptosis assays. Specific inhibition of PKR with the chemical inhibitor 2-AP restores radiation sensitivity. Plasmid transfection of the PKR wild-type (wt) gene into PKR(-/-) MEFs leads to increased radiation resistance. The protective effect of PKR to radiation may be mediated in part through NF-kappaB and Akt because both NF-kappaB and Akt are activated after ionizing radiation in PKR+/+ but not PKR-/- cells. CONCLUSIONS: We suggest a novel role for PKR as a mediator of radiation resistance modulated in part through the protective effects of NF-kappaB and Akt activation. The modification of PKR activity may be a novel strategy in the future to overcome radiation resistance.