Comparison of amplification methods for transcriptomic analyses of low abundance prokaryotic RNA sources

Microarrays have established as instrumental for bacterial detection, identification, and genotyping as well as for transcriptomic studies. For gene expression analyses using limited numbers of bacteria (derived from in vivo or ex vivo origin, for example), RNA amplification is often required prior to labeling and hybridization onto microarrays. Evaluation of the fidelity of the amplification methods is crucial for the robustness and reproducibility of microarray results. We report here the first utilization of random primers and the highly processive Phi29 phage polymerase to amplify material for transcription profiling analyses. We compared two commercial amplification methods (GenomiPhi and MessageAmp kits) with direct reverse-transcription as the reference method, focusing on the robustness of mRNA quantification using either microarrays or quantitative RT-PCR. Both amplification methods using either poly-A tailing followed by in vitro transcription, or direct strand displacement polymerase, showed appreciable linearity. Strand displacement technique was particularly affordable compared to in vitro transcription-based (IVT) amplification methods and consisted in a single tube reaction leading to high amplification yields. Real-time measurements using low-, medium-, and highly expressed genes revealed that this simple method provided linear amplification with equivalent results in terms of relative messenger abundance as those obtained by conventional direct reverse-transcription.

Tissue metabolomics of hepatocellular carcinoma: Tumor energy metabolism and the role of transcriptomic classification

Hepatocellular carcinoma (HCC) is one of the commonest causes of death from cancer. A plethora of metabolomic investigations of HCC have yielded molecules in biofluids that are both up- and down-regulated but no real consensus has emerged regarding exploitable biomarkers for early detection of HCC. We report here a different approach, a combined transcriptomics and metabolomics study of energy metabolism in HCC. A panel of 31 pairs of HCC tumors and corresponding nontumor liver tissues from the same patients was investigated by gas chromatography-mass spectrometry (GCMS)-based metabolomics. HCC was characterized by ∼2-fold depletion of glucose, glycerol 3- and 2-phosphate, malate, alanine, myo-inositol, and linoleic acid. Data are consistent with a metabolic remodeling involving a 4-fold increase in glycolysis over mitochondrial oxidative phosphorylation. A second panel of 59 HCC that had been typed by transcriptomics and classified in G1 to G6 subgroups was also subjected to GCMS tissue metabolomics. No differences in glucose, lactate, alanine, glycerol 3-phosphate, malate, myo-inositol, or stearic acid tissue concentrations were found, suggesting that the Wnt/β-catenin pathway activated by CTNNB1 mutation in subgroups G5 and G6 did not exhibit specific metabolic remodeling. However, subgroup G1 had markedly reduced tissue concentrations of 1-stearoylglycerol, 1-palmitoylglycerol, and palmitic acid, suggesting that the high serum α-fetoprotein phenotype of G1, associated with the known overexpression of lipid catabolic enzymes, could be detected through metabolomics as increased lipid catabolism. Conclusion: Tissue metabolomics yielded precise biochemical information regarding HCC tumor metabolic remodeling from mitochondrial oxidation to aerobic glycolysis and the impact of molecular subtypes on this process.

Profiling invasiveness in head and neck cancer: recent contributions of genomic and transcriptomic approaches.

High-throughput molecular profiling approaches have emerged as precious research tools in the field of head and neck translational oncology. Such approaches have identified and/or confirmed the role of several genes or pathways in the acquisition/maintenance of an invasive phenotype and the execution of cellular programs related to cell invasion. Recently published new-generation sequencing studies in head and neck squamous cell carcinoma (HNSCC) have unveiled prominent roles in carcinogenesis and cell invasion of mutations involving NOTCH1 and PI3K-patwhay components. Gene-expression profiling studies combined with systems biology approaches have allowed identifying and gaining further mechanistic understanding into pathways commonly enriched in invasive HNSCC. These pathways include antigen-presenting and leucocyte adhesion molecules, as well as genes involved in cell-extracellular matrix interactions. Here we review the major insights into invasiveness in head and neck cancer provided by high-throughput molecular profiling approaches.

Tracking the transcriptional host response from the acute to the regenerative phase of experimental pneumococcal meningitis

Despite the availability of effective antibiotic therapies, pneumococcal meningitis (PM) has a case fatality rate of up to 30% and causes neurological sequelae in up to half of the surviving patients. The underlying brain damage includes apoptosis of neurons in the hippocampus and necrosis in the cortex. Therapeutic options to reduce acute injury and to improve outcome from PM are severely limited.With the aim to develop new therapies a number of pharmacologic interventions have been evaluated. However, the often unpredictable outcome of interventional studies suggests that the current concept of the pathophysiologic events during bacterial meningitis is fragmentary. The aim of this work is to describe the transcriptomic changes underlying the complex mechanisms of the host response to pneumococcal meningitis in a temporal and spatial context using a well characterized infant rat model.

A subset of metzincins and related genes constitutes a marker of human solid organ fibrosis

Metzincins and functionally related genes play important roles in extracellular matrix remodeling both in healthy and fibrotic conditions. We recently presented a transcriptomic classifier consisting of 19 metzincins and related genes (MARGS) discriminating biopsies from renal transplant patients with or without interstitial fibrosis/tubular atrophy (IF/TA) by virtue of gene expression measurement (Roedder et al., Am J Transplant 9:517-526, 2009). Here we demonstrate that the same algorithm has diagnostic value in non-transplant solid organ fibrosis. We used publically available microarray datasets of 325 human heart, liver, lung, kidney cortex, and pancreas microarray samples (265 with fibrosis, 60 healthy controls). Expression of nine commonly differentially expressed genes was confirmed by TaqMan low-density arrays (Applied Biosystems, USA) in 50 independent archival tissue specimens with matched histological diagnoses to microarray patients. In separate and in combined, integrated microarray data analyses of five datasets with 325 samples, the previously published MARGS classifier for renal post-transplant IF/TA had a mean AUC of 87% and 82%, respectively. These data demonstrate that the MARGS gene panel classifier not only discriminates IF/TA from normal renal transplant tissue, but also classifies solid organ fibrotic conditions of human pancreas, liver, heart, kidney, and lung tissue samples with high specificity and accuracy, suggesting that the MARGS classifier is a cross-platform, cross-organ classifier of fibrotic conditions of different etiologies when compared to normal tissue.

Life stage differences in mammary gland gene expression profile in non-human primates

Breast cancer (BC) is the most common malignancy of women in the developed world. To better understand its pathogenesis, knowledge of normal breast development is crucial, as BC is the result of disregulation of physiologic processes. The aim of this study was to investigate the impact of reproductive life stages on the transcriptional profile of the mammary gland in a primate model. Comparative transcriptomic analyses were carried out using breast tissues from 28 female cynomolgus macaques (Macaca fascicularis) at the following life stages: prepubertal (n = 5), adolescent (n = 4), adult luteal (n = 5), pregnant (n = 6), lactating (n = 3), and postmenopausal (n = 5). Mammary gland RNA was hybridized to Affymetrix GeneChip(®) Rhesus Macaque Genome Arrays. Differential gene expression was analyzed using ANOVA and cluster analysis. Hierarchical cluster analysis revealed distinct separation of life stage groups. More than 2,225 differentially expressed mRNAs were identified. Gene families or pathways that changed across life stages included those related to estrogen and androgen (ESR1, PGR, TFF1, GREB1, AR, 17HSDB2, 17HSDB7, STS, HSD11B1, AKR1C4), prolactin (PRLR, ELF5, STAT5, CSN1S1), insulin-like growth factor signaling (IGF1, IGFBP1, IGFBP5), extracellular matrix (POSTN, TGFB1, COL5A2, COL12A1, FOXC1, LAMC1, PDGFRA, TGFB2), and differentiation (CD24, CD29, CD44, CD61, ALDH1, BRCA1, FOXA1, POSTN, DICER1, LIG4, KLF4, NOTCH2, RIF1, BMPR1A, TGFB2). Pregnancy and lactation displayed distinct patterns of gene expression. ESR1 and IGF1 were significantly higher in the adolescent compared to the adult animals, whereas differentiation pathways were overrepresented in adult animals and pregnancy-associated life stages. Few individual genes were distinctly different in postmenopausal animals. Our data demonstrate characteristic patterns of gene expression during breast development. Several of the pathways activated during pubertal development have been implicated in cancer development and metastasis, supporting the idea that other developmental markers may have application as biomarkers for BC.

A global view of Staphylococcus aureus whole genome expression upon internalization in human epithelial cells

BACKGROUND: Staphylococcus aureus, a leading cause of chronic or acute infections, is traditionally considered an extracellular pathogen despite repeated reports of S. aureus internalization by a variety of non-myeloid cells in vitro. This property potentially contributes to bacterial persistence, protection from antibiotics and evasion of immune defenses. Mechanisms contributing to internalization have been partly elucidated, but bacterial processes triggered intracellularly are largely unknown. RESULTS: We have developed an in vitro model using human lung epithelial cells that shows intracellular bacterial persistence for up to 2 weeks. Using an original approach we successfully collected and amplified low amounts of bacterial RNA recovered from infected eukaryotic cells. Transcriptomic analysis using an oligoarray covering the whole S. aureus genome was performed at two post-internalization times and compared to gene expression of non-internalized bacteria. No signs of cellular death were observed after prolonged internalization of Staphylococcus aureus 6850 in epithelial cells. Following internalization, extensive alterations of bacterial gene expression were observed. Whereas major metabolic pathways including cell division, nutrient transport and regulatory processes were drastically down-regulated, numerous genes involved in iron scavenging and virulence were up-regulated. This initial adaptation was followed by a transcriptional increase in several metabolic functions. However, expression of several toxin genes known to affect host cell integrity appeared strictly limited. CONCLUSION: These molecular insights correlated with phenotypic observations and demonstrated that S. aureus modulates gene expression at early times post infection to promote survival. Staphylococcus aureus appears adapted to intracellular survival in non-phagocytic cells.

UPLC-ESI-TOFMS-based metabolomics and gene expression dynamics inspector self-organizing metabolomic maps as tools for understanding the cellular response to ionizing radiation

Global transcriptomic and proteomic profiling platforms have yielded important insights into the complex response to ionizing radiation (IR). Nonetheless, little is known about the ways in which small cellular metabolite concentrations change in response to IR. Here, a metabolomics approach using ultraperformance liquid chromatography coupled with electrospray time-of-flight mass spectrometry was used to profile, over time, the hydrophilic metabolome of TK6 cells exposed to IR doses ranging from 0.5 to 8.0 Gy. Multivariate data analysis of the positive ions revealed dose- and time-dependent clustering of the irradiated cells and identified certain constituents of the water-soluble metabolome as being significantly depleted as early as 1 h after IR. Tandem mass spectrometry was used to confirm metabolite identity. Many of the depleted metabolites are associated with oxidative stress and DNA repair pathways. Included are reduced glutathione, adenosine monophosphate, nicotinamide adenine dinucleotide, and spermine. Similar measurements were performed with a transformed fibroblast cell line, BJ, and it was found that a subset of the identified TK6 metabolites were effective in IR dose discrimination. The GEDI (Gene Expression Dynamics Inspector) algorithm, which is based on self-organizing maps, was used to visualize dynamic global changes in the TK6 metabolome that resulted from IR. It revealed dose-dependent clustering of ions sharing the same trends in concentration change across radiation doses. "Radiation metabolomics," the application of metabolomic analysis to the field of radiobiology, promises to increase our understanding of cellular responses to stressors such as radiation.

Proteomic analysis of the Theileria annulata schizont

The apicomplexan parasite, Theileria annulata, is the causative agent of tropical theileriosis, a devastating lymphoproliferative disease of cattle. The schizont stage transforms bovine leukocytes and provides an intriguing model to study host/pathogen interactions. The genome of T. annulata has been sequenced and transcriptomic data are rapidly accumulating. In contrast, little is known about the proteome of the schizont, the pathogenic, transforming life cycle stage of the parasite. Using one-dimensional (1-D) gel LC-MS/MS, a proteomic analysis of purified T. annulata schizonts was carried out. In whole parasite lysates, 645 proteins were identified. Proteins with transmembrane domains (TMDs) were under-represented and no proteins with more than four TMDs could be detected. To tackle this problem, Triton X-114 treatment was applied, which facilitates the extraction of membrane proteins, followed by 1-D gel LC-MS/MS. This resulted in the identification of an additional 153 proteins. Half of those had one or more TMD and 30 proteins with more than four TMDs were identified. This demonstrates that Triton X-114 treatment can provide a valuable additional tool for the identification of new membrane proteins in proteomic studies. With two exceptions, all proteins involved in glycolysis and the citric acid cycle were identified. For at least 29% of identified proteins, the corresponding transcripts were not present in the existing expressed sequence tag databases. The proteomics data were integrated into the publicly accessible database resource at EuPathDB (www.eupathdb.org) so that mass spectrometry-based protein expression evidence for T. annulata can be queried alongside transcriptional and other genomics data available for these parasites.

Early gene expression analysis in 9L orthotopic tumor-bearing rats identifies immune modulation in molecular response to synchrotron microbeam radiation therapy

Synchrotron Microbeam Radiation Therapy (MRT) relies on the spatial fractionation of the synchrotron photon beam into parallel micro-beams applying several hundred of grays in their paths. Several works have reported the therapeutic interest of the radiotherapy modality at preclinical level, but biological mechanisms responsible for the described efficacy are not fully understood to date. The aim of this study was to identify the early transcriptomic responses of normal brain and glioma tissue in rats after MRT irradiation (400Gy). The transcriptomic analysis of similarly irradiated normal brain and tumor tissues was performed 6 hours after irradiation of 9 L orthotopically tumor-bearing rats. Pangenomic analysis revealed 1012 overexpressed and 497 repressed genes in the irradiated contralateral normal tissue and 344 induced and 210 repressed genes in tumor tissue. These genes were grouped in a total of 135 canonical pathways. More than half were common to both tissues with a predominance for immunity or inflammation (64 and 67% of genes for normal and tumor tissues, respectively). Several pathways involving HMGB1, toll-like receptors, C-type lectins and CD36 may serve as a link between biochemical changes triggered by irradiation and inflammation and immunological challenge. Most immune cell populations were involved: macrophages, dendritic cells, natural killer, T and B lymphocytes. Among them, our results highlighted the involvement of Th17 cell population, recently described in tumor. The immune response was regulated by a large network of mediators comprising growth factors, cytokines, lymphokines. In conclusion, early response to MRT is mainly based on inflammation and immunity which appear therefore as major contributors to MRT efficacy.

A review on bovine besnoitiosis: a disease with economic impact in herd health management, caused by Besnoitia besnoiti (Franco and Borges, 1916).

Bovine besnoitiosis is caused by the largely unexplored apicomplexan parasite Besnoitia besnoiti. In cows, infection during pregnancy often results in abortion, and chronically infected bulls become infertile. Similar to other apicomplexans B. besnoiti has acquired a largely intracellular lifestyle, but its complete life cycle is still unknown, modes of transmission have not been entirely resolved and the definitive host has not been identified. Outbreaks of bovine besnoitiosis in cattle were described in the 1990s in Portugal and Spain, and later several cases were also detected in France. More cases have been reported recently in hitherto unaffected countries, including Italy, Germany, Switzerland, Hungary and Croatia. To date, there is still no effective pharmaceutical compound available for the treatment of besnoitiosis in cattle, and progress in the identification of novel targets for intervention through pharmacological or immunological means is hampered by the lack of molecular data on the genomic and transcriptomic level. In addition, the lack of an appropriate small animal laboratory model, and wide gaps in our knowledge on the host-parasite interplay during the life cycle of this parasite, renders vaccine and drug development a cost- and labour-intensive undertaking.

The transcriptome of equine peripheral blood mononuclear cells.

Complete transcriptomic data at high resolution are available only for a few model organisms with medical importance. The gene structures of non-model organisms are mostly computationally predicted based on comparative genomics with other species. As a result, more than half of the horse gene models are known only by projection. Experimental data supporting these gene models are scarce. Moreover, most of the annotated equine genes are single-transcript genes. Utilizing RNA sequencing (RNA-seq) the experimental validation of predicted transcriptomes has become accessible at reasonable costs. To improve the horse genome annotation we performed RNA-seq on 561 samples of peripheral blood mononuclear cells (PBMCs) derived from 85 Warmblood horses. The mapped sequencing reads were used to build a new transcriptome assembly. The new assembly revealed many alternative isoforms associated to known genes or to those predicted by the Ensembl and/or Gnomon pipelines. We also identified 7,531 transcripts not associated with any horse gene annotated in public databases. Of these, 3,280 transcripts did not have a homologous match to any sequence deposited in the NCBI EST database suggesting horse specificity. The unknown transcripts were categorized as coding and noncoding based on predicted coding potential scores. Among them 230 transcripts had high coding potential score, at least 2 exons, and an open reading frame of at least 300 nt. We experimentally validated 9 new equine coding transcripts using RT-PCR and Sanger sequencing. Our results provide valuable detailed information on many transcripts yet to be annotated in the horse genome.

TrypanoCyc: a community-led biochemical pathways database for Trypanosoma brucei.

The metabolic network of a cell represents the catabolic and anabolic reactions that interconvert small molecules (metabolites) through the activity of enzymes, transporters and non-catalyzed chemical reactions. Our understanding of individual metabolic networks is increasing as we learn more about the enzymes that are active in particular cells under particular conditions and as technologies advance to allow detailed measurements of the cellular metabolome. Metabolic network databases are of increasing importance in allowing us to contextualise data sets emerging from transcriptomic, proteomic and metabolomic experiments. Here we present a dynamic database, TrypanoCyc (http://www.metexplore.fr/trypanocyc/), which describes the generic and condition-specific metabolic network of Trypanosoma brucei, a parasitic protozoan responsible for human and animal African trypanosomiasis. In addition to enabling navigation through the BioCyc-based TrypanoCyc interface, we have also implemented a network-based representation of the information through MetExplore, yielding a novel environment in which to visualise the metabolism of this important parasite.

Integrated mRNA and miRNA expression profiling in blood reveals candidate biomarkers associated with endurance exercise in the horse.

The adaptive response to extreme endurance exercise might involve transcriptional and translational regulation by microRNAs (miRNAs). Therefore, the objective of the present study was to perform an integrated analysis of the blood transcriptome and miRNome (using microarrays) in the horse before and after a 160 km endurance competition. A total of 2,453 differentially expressed genes and 167 differentially expressed microRNAs were identified when comparing pre- and post-ride samples. We used a hypergeometric test and its generalization to gain a better understanding of the biological functions regulated by the differentially expressed microRNA. In particular, 44 differentially expressed microRNAs putatively regulated a total of 351 depleted differentially expressed genes involved variously in glucose metabolism, fatty acid oxidation, mitochondrion biogenesis, and immune response pathways. In an independent validation set of animals, graphical Gaussian models confirmed that miR-21-5p, miR-181b-5p and miR-505-5p are candidate regulatory molecules for the adaptation to endurance exercise in the horse. To the best of our knowledge, the present study is the first to provide a comprehensive, integrated overview of the microRNA-mRNA co-regulation networks that may have a key role in controlling post-transcriptomic regulation during endurance exercise.