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

Metatranscriptomics of N2-fixing cyanobacteria in the Amazon River plume

Biological N2 fixation is an important nitrogen source for surface ocean microbial communities. However, nearly all information on the diversity and gene expression of organisms responsible for oceanic N2 fixation in the environment has come from targeted approaches that assay only a small number of genes and organisms. Using genomes of diazotrophic cyanobacteria to extract reads from extensive meta-genomic and -transcriptomic libraries, we examined diazotroph diversity and gene expression from the Amazon River plume, an area characterized by salinity and nutrient gradients. Diazotroph genome and transcript sequences were most abundant in the transitional waters compared with lower salinity or oceanic water masses. We were able to distinguish two genetically divergent phylotypes within the Hemiaulus-associated Richelia sequences, which were the most abundant diazotroph sequences in the data set. Photosystem (PS)-II transcripts in Richelia populations were much less abundant than those in Trichodesmium, and transcripts from several Richelia PS-II genes were absent, indicating a prominent role for cyclic electron transport in Richelia. In addition, there were several abundant regulatory transcripts, including one that targets a gene involved in PS-I cyclic electron transport in Richelia. High sequence coverage of the Richelia transcripts, as well as those from Trichodesmium populations, allowed us to identify expressed regions of the genomes that had been overlooked by genome annotations. High-coverage genomic and transcription analysis enabled the characterization of distinct phylotypes within diazotrophic populations, revealed a distinction in a core process between dominant populations and provided evidence for a prominent role for noncoding RNAs in microbial communities.

Microarray analysis of microbiota of gingival lesions in noma patients

Noma (cancrum oris) is a gangrenous disease of unknown etiology affecting the maxillo-facial region of young children in extremely limited resource countries. In an attempt to better understand the microbiological events occurring during this disease, we used phylogenetic and low-density microarrays targeting the 16S rRNA gene to characterize the gingival flora of acute noma and acute necrotizing gingivitis (ANG) lesions, and compared them to healthy control subjects of the same geographical and social background. Our observations raise doubts about Fusobacterium necrophorum, a previously suspected causative agent of noma, as this species was not associated with noma lesions. Various oral pathogens were more abundant in noma lesions, notably Atopobium spp., Prevotella intermedia, Peptostreptococcus spp., Streptococcus pyogenes and Streptococcus anginosus. On the other hand, pathogens associated with periodontal diseases such as Aggregatibacter actinomycetemcomitans, Capnocytophaga spp., Porphyromonas spp. and Fusobacteriales were more abundant in healthy controls. Importantly, the overall loss of bacterial diversity observed in noma samples as well as its homology to that of ANG microbiota supports the hypothesis that ANG might be the immediate step preceding noma.

Extraction of bacterial RNA from soil: challenges and solutions

Detection of bacterial gene expression in soil emerged in the early 1990s and provided information on bacterial responses in their original soil environments. As a key procedure in the detection, extraction of bacterial RNA from soil has attracted much interest, and many methods of soil RNA extraction have been reported in the past 20 years. In addition to various RT-PCR-based technologies, new technologies for gene expression analysis, such as microarrays and high-throughput sequencing technologies, have recently been applied to examine bacterial gene expression in soil. These technologies are driving improvements in RNA extraction protocols. In this mini-review, progress in the extraction of bacterial RNA from soil is summarized with emphasis on the major difficulties in the development of methodologies and corresponding strategies to overcome them.

Changes in the Staphylococcus aureus transcriptome during early adaptation to the lung

Staphylococcus aureus is a common inhabitant of the human nasopharynx. It is also a cause of life-threatening illness, producing a potent array of virulence factors that enable survival in normally sterile sites. The transformation of S. aureus from commensal to pathogen is poorly understood. We analyzed S. aureus gene expression during adaptation to the lung using a mouse model of S. aureus pneumonia. Bacteria were isolated by bronchoalveolar lavage after residence in vivo for up to 6 hours. S. aureus in vivo RNA transcription was compared by microarray to that of shake flask grown stationary phase and early exponential phase cells. Compared to in vitro conditions, the in vivo transcriptome was dramatically altered within 30 minutes. Expression of central metabolic pathways changed significantly in response to the lung environment. Gluconeogenesis (fbs, pckA) was down regulated, as was TCA cycle and fermentation pathway gene expression. Genes associated with amino acid synthesis, RNA translation and nitrate respiration were upregulated, indicative of a highly active metabolic state during the first 6 hours in the lung. Virulence factors regulated by agr were down regulated in vivo and in early exponential phase compared to stationary phase cells. Over time in vivo, expression of ahpCF, involved in H₂O₂ scavenging, and uspA, which encodes a universal stress regulator, increased. Transcription of leukotoxic α and β-type phenol-soluble modulins psmα1-4 and psmβ1-2 increased 13 and 8-fold respectively; hld mRNA, encoding δ-hemolysin, was increased 9-fold. These were the only toxins to be significantly upregulated in vivo. These data provide the first complete survey of the S. aureus transcriptome response to the mammalian airway. The results present intriguing contrasts with previous work in other in vitro and in vivo models and provide novel insights into the adaptive and temporal response of S. aureus early in the pathogenesis of pneumonia.

Influence of RNA labeling on expression profiling of microRNAs

Although a number of technical parameters are now being examined to optimize microRNA profiling experiments, it is unknown whether reagent or component changes to the labeling step affect starting RNA requirements or microarray performance. Human brain/lung samples were each labeled in duplicate, at 1.0, 0.5, 0.2, and 0.1 μg of total RNA, by means of two kits that use the same labeling procedure but differ in the reagent composition used to label microRNAs. Statistical measures of reliability and validity were used to evaluate microarray data. Cross-platform confirmation was accomplished using TaqMan microRNA assays. Synthetic microRNA spike-in experiments were also performed to establish the microarray signal dynamic range using the ligation-modified kit. Technical replicate correlations of signal intensity values were high using both kits, but improved with the ligation-modified assay. The drop in detection call sensitivity and miRNA gene list correlations, when using reduced amounts of standard-labeled RNA, was considerably improved with the ligation-modified kit. Microarray signal dynamic range was found to be linear across three orders of magnitude from 4.88 to 5000 attomoles. Thus, optimization of the microRNA labeling reagent can result in at least a 10-fold decrease in microarray total RNA requirements with little compromise to data quality. Clinical investigations bottlenecked by the amount of starting material may use a ligation mix modification strategy to reduce total RNA requirements.

In-vivo expression profiling of Pseudomonas aeruginosa infections reveals niche-specific and strain-independent transcriptional programs

Pseudomonas aeruginosa is a threatening, opportunistic pathogen causing disease in immunocompromised individuals. The hallmark of P. aeruginosa virulence is its multi-factorial and combinatorial nature. It renders such bacteria infectious for many organisms and it is often resistant to antibiotics. To gain insights into the physiology of P. aeruginosa during infection, we assessed the transcriptional programs of three different P. aeruginosa strains directly after isolation from burn wounds of humans. We compared the programs to those of the same strains using two infection models: a plant model, which consisted of the infection of the midrib of lettuce leaves, and a murine tumor model, which was obtained by infection of mice with an induced tumor in the abdomen. All control conditions of P. aeruginosa cells growing in suspension and as a biofilm were added to the analysis. We found that these different P. aeruginosa strains express a pool of distinct genetic traits that are activated under particular infection conditions regardless of their genetic variability. The knowledge herein generated will advance our understanding of P. aeruginosa virulence and provide valuable cues for the definition of prospective targets to develop novel intervention strategies.

Transcriptome analysis of Neisseria meningitidis in human whole blood and mutagenesis studies identify virulence factors involved in blood survival

During infection Neisseria meningitidis (Nm) encounters multiple environments within the host, which makes rapid adaptation a crucial factor for meningococcal survival. Despite the importance of invasion into the bloodstream in the meningococcal disease process, little is known about how Nm adapts to permit survival and growth in blood. To address this, we performed a time-course transcriptome analysis using an ex vivo model of human whole blood infection. We observed that Nm alters the expression of ≈30% of ORFs of the genome and major dynamic changes were observed in the expression of transcriptional regulators, transport and binding proteins, energy metabolism, and surface-exposed virulence factors. In particular, we found that the gene encoding the regulator Fur, as well as all genes encoding iron uptake systems, were significantly up-regulated. Analysis of regulated genes encoding for surface-exposed proteins involved in Nm pathogenesis allowed us to better understand mechanisms used to circumvent host defenses. During blood infection, Nm activates genes encoding for the factor H binding proteins, fHbp and NspA, genes encoding for detoxifying enzymes such as SodC, Kat and AniA, as well as several less characterized surface-exposed proteins that might have a role in blood survival. Through mutagenesis studies of a subset of up-regulated genes we were able to identify new proteins important for survival in human blood and also to identify additional roles of previously known virulence factors in aiding survival in blood. Nm mutant strains lacking the genes encoding the hypothetical protein NMB1483 and the surface-exposed proteins NalP, Mip and NspA, the Fur regulator, the transferrin binding protein TbpB, and the L-lactate permease LctP were sensitive to killing by human blood. This increased knowledge of how Nm responds to adaptation in blood could also be helpful to develop diagnostic and therapeutic strategies to control the devastating disease cause by this microorganism.

Integrative and comparative genomics analysis of early hepatocellular carcinoma differentiated from liver regeneration in young and old

BACKGROUND

Hepatocellular carcinoma (HCC) is the third-leading cause of cancer-related deaths worldwide. It is often diagnosed at an advanced stage, and hence typically has a poor prognosis. To identify distinct molecular mechanisms for early HCC we developed a rat model of liver regeneration post-hepatectomy, as well as liver cells undergoing malignant transformation and compared them to normal liver using a microarray approach. Subsequently, we performed cross-species comparative analysis coupled with copy number alterations (CNA) of independent early human HCC microarray studies to facilitate the identification of critical regulatory modules conserved across species.

RESULTS

We identified 35 signature genes conserved across species, and shared among different types of early human HCCs. Over 70% of signature genes were cancer-related, and more than 50% of the conserved genes were mapped to human genomic CNA regions. Functional annotation revealed genes already implicated in HCC, as well as novel genes which were not previously reported in liver tumors. A subset of differentially expressed genes was validated using quantitative RT-PCR. Concordance was also confirmed for a significant number of genes and pathways in five independent validation microarray datasets. Our results indicated alterations in a number of cancer related pathways, including p53, p38 MAPK, ERK/MAPK, PI3K/AKT, and TGF-beta signaling pathways, and potential critical regulatory role of MYC, ERBB2, HNF4A, and SMAD3 for early HCC transformation.

CONCLUSIONS

The integrative analysis of transcriptional deregulation, genomic CNA and comparative cross species analysis brings new insights into the molecular profile of early hepatoma formation. This approach may lead to robust biomarkers for the detection of early human HCC.

An optimized RNA amplification method for prokaryotic expression profiling analysis

DNA microarray technology has been extensively used for gene expression analysis of both eukaryotic and prokaryotic organisms. For eukaryotic gene expression profiling, the poly(A)-based reverse transcription of messenger RNA (mRNA) followed by T7 RNA polymerase-based in vitro transcription is generally required to produce enough RNA targets for hybridization with the microarray chips. However, the same method cannot be directly applied to prokaryotic mRNAs due to the lack of poly(A) sequences at the 3' ends. Conventional methods usually require large amounts of starting RNAs and lead to high background noise. Recently developed amplification methods enable smaller amounts of prokaryotic RNA to be used from samples with species-specific primers, oligo(dT) primers, or random primers. In this study, three target preparation methods, including the direct labeling, polyadenylation-involved oligo-dT priming, and random priming amplification (respectively referred to as DL, PAOD, and RPA hereafter) were evaluated through expression profiling of a heat shock model of Escherichia coli. The PAOD method was found to be more sensitive and more specific in differential gene expression measurements than either DL and RPA, even when the E. coli RNA was only a small proportion of the simulated eukaryotic host RNA. The results suggest that PAOD is the preferred target preparation method for prokaryotic transcriptome.

Review of microarray studies for host-intracellular pathogen interactions

Obligate intracellular bacteria are privileged soldiers on the battlefield that represent host-pathogen interactions. Microarrays are a powerful technology that can increase our knowledge about how bacteria respond to and interact with their hosts. This review summarizes the limitations inherent to host-pathogen interaction studies and essential strategies to improve microarray investigations of intracellular bacteria. We have compiled the comparative genomic and gene expression analyses of obligate intracellular bacteria currently available from microarrays. In this review we explore ways in which microarrays can be used to identify polymorphisms in different obligate intracellular bacteria such as Coxiella burnetii, Chlamydia trachomatis, Ehrlichia chaffeensis, Rickettsia prowazekii and Tropheryma whipplei. These microarray studies reveal that, while genomic content is highly conserved in obligate intracellular bacteria, genetic polymorphisms can potentially occur to increase bacterial pathogenesis. Additionally, changes in the gene expression of C. trachomatis throughout its life cycle, as well as changes in the gene expression profile of the pathogens R. prowazekii, Rickettsia rickettsii, Rickettsia typhi, T. whipplei and C. trachomatis following environmental changes, are discussed. Finally, an in vivo model of Rickettsia conorii within the skin is discussed. The gene expression analyses highlight the capacity of obligate intracellular bacteria to adapt to environmental changes and potentially to thwart the host response.

Left ventricular global transcriptional profiling in human end-stage dilated cardiomyopathy

We employed ABI high-density oligonucleotide microarrays containing 31,700 sixty-mer probes (representing 27,868 annotated human genes) to determine differential gene expression in idiopathic dilated cardiomyopathy (DCM). We identified 626 up-regulated and 636 down-regulated genes in DCM compared to controls. Most significant changes occurred in the tricarboxylic acid cycle, angiogenesis, and apoptotic signaling pathways, among which 32 apoptosis- and 13 MAPK activity-related genes were altered. Inorganic cation transporter, catalytic activities, energy metabolism and electron transport-related processes were among the most critically influenced pathways. Among the up-regulated genes were HTRA1 (6.9-fold), PDCD8(AIFM1) (5.2) and PRDX2 (4.4) and the down-regulated genes were NR4A2 (4.8), MX1 (4.3), LGALS9 (4), IFNA13 (4), UNC5D (3.6) and HDAC2 (3) (p<0.05), all of which have no clearly defined cardiac-related function yet. Gene ontology and enrichment analysis also revealed significant alterations in mitochondrial oxidative phosphorylation, metabolism and Alzheimer's disease pathways. Concordance was also confirmed for a significant number of genes and pathways in an independent validation microarray dataset. Furthermore, verification by real-time RT-PCR showed a high degree of consistency with the microarray results. Our data demonstrate an association of DCM with alterations in various cellular events and multiple yet undeciphered genes that may contribute to heart muscle disease pathways.

Rickettsia conorii transcriptional response within inoculation eschar

Background

Rickettsia conorii, the causative agent of the Mediterranean spotted fever, is transmitted to humans by the bite of infected ticks Rhipicephalus sanguineus. The skin thus constitutes an important barrier for the entry and propagation of R. conorii. Given this, analysis of the survival strategies used by the bacterium within infected skin is critical for our understanding of rickettsiosis.

Methodology/Principal Findings

Here, we report the first genome-wide analysis of R. conorii gene expression from infected human skin biopsies. Our data showed that R. conorii exhibited a striking transcript signature that is remarkably conserved across patients, regardless of genotype. The expression profiles obtained using custom Agilent microarrays were validated by quantitative RT-PCR. Within eschars, the amount of detected R. conorii transcripts was of 55%, this value being of 74% for bacteria grown in Vero cells. In such infected host tissues, approximately 15% (n = 211) of the total predicted R. conorii ORFs appeared differentially expressed compared to bacteria grown in standard laboratory conditions. These genes are mostly down-regulated and encode proteins essential for bacterial replication. Some of the strategies displayed by rickettsiae to overcome the host defense barriers, thus avoiding killing, were also pointed out. The observed up-regulation of rickettsial genes associated with DNA repair is likely to correspond to a DNA-damaging agent enriched environment generated by the host cells to eradicate the pathogens. Survival of R. conorii within eschars also involves adaptation to osmotic stress, changes in cell surface proteins and up-regulation of some virulence factors. Interestingly, in contrast to down-regulated transcripts, we noticed that up-regulated ones rather exhibit a small nucleotide size, most of them being exclusive for the spotted fever group rickettsiae.

Conclusion/Significance

Because eschar is a site for rickettsial introduction, the pattern of rickettsial gene expression observed here may define how rickettsiae counteract the host defense.

Detection of large numbers of novel sequences in the metatranscriptomes of complex marine microbial communities

Background

Sequencing the expressed genetic information of an ecosystem (metatranscriptome) can provide information about the response of organisms to varying environmental conditions. Until recently, metatranscriptomics has been limited to microarray technology and random cloning methodologies. The application of high-throughput sequencing technology is now enabling access to both known and previously unknown transcripts in natural communities.

Methodology/Principal Findings

We present a study of a complex marine metatranscriptome obtained from random whole-community mRNA using the GS-FLX Pyrosequencing technology. Eight samples, four DNA and four mRNA, were processed from two time points in a controlled coastal ocean mesocosm study (Bergen, Norway) involving an induced phytoplankton bloom producing a total of 323,161,989 base pairs. Our study confirms the finding of the first published metatranscriptomic studies of marine and soil environments that metatranscriptomics targets highly expressed sequences which are frequently novel. Our alternative methodology increases the range of experimental options available for conducting such studies and is characterized by an exceptional enrichment of mRNA (99.92%) versus ribosomal RNA. Analysis of corresponding metagenomes confirms much higher levels of assembly in the metatranscriptomic samples and a far higher yield of large gene families with >100 members, ∼91% of which were novel.

Conclusions/Significance

This study provides further evidence that metatranscriptomic studies of natural microbial communities are not only feasible, but when paired with metagenomic data sets, offer an unprecedented opportunity to explore both structure and function of microbial communities – if we can overcome the challenges of elucidating the functions of so many never-seen-before gene families.

Comparison of detection and signal amplification methods for DNA microarrays

One of the factors limiting the use of DNA microarray technology for the detection of pathogenic organisms from clinical and environmental matrices has been inadequate assay sensitivity. To assess the effectiveness of post-hybridization secondary detection steps to enhance the sensitivity of DNA microarray-based pathogen detection, we evaluated a panel of 11 commercial and novel hybridization detection and signal amplification methods (direct labeling, indirect aminoallyl labeling, antibody, DNA dendrimers, viral particles, internally fluorescent nanoparticles, tyramide signal amplification, resonance light scattering nanoparticles and quantum dots) using a multiplex PCR and spotted long oligonucleotide microarray for Vibrio cholerae. Quantitative parameters such as sensitivity, signal intensity, background, assay complexity, time and cost were assessed and provide comparative criteria to be considered for DNA microarray experimental design. While the most important parameter is likely to vary based on the assay, when weighted equally, the findings suggest that recognition element- and dye-functionalized viral particles provide the most attractive option for microarray detection and signal amplification.

Regulation of whole bacterial pathogen transcription within infected hosts

DNA microarrays are a powerful and promising approach to gain a detailed understanding of the bacterial response and the molecular cross-talk that can occur as a consequence of host-pathogen interactions. However, published studies mainly describe the host response to infection. Analysis of bacterial gene regulation in the course of infection has confronted many challenges. This review summarizes the different strategies used over the last few years to investigate, at the genomic scale, and using microarrays, the alterations in the bacterial transcriptome in response to interactions with host cells. Thirty-seven studies involving 19 different bacterial pathogens were compiled and analyzed. Our in silico comparison of the transcription profiles of bacteria grown in broth or in contact with eukaryotic cells revealed some features commonly observed when bacteria interact with host cells, including stringent response and cell surface remodeling.

Development of a method for recovering rickettsial RNA from infected cells to analyze gene expression profiling of obligate intracellular bacteria

The Rickettsia genus is composed of Gram-negative bacteria responsible for Typhus and spotted fevers. Because of the limitations imposed by their obligate intracellular location, the molecular mechanisms responsible for their pathogenicity remain poorly understood. Several rickettsial genomes are now available, thus providing the foundation for a new era of post-genomic research. Here, using Rickettsia conorii as model, we developed a suitable method for microarray-based transcriptome analysis of rickettsiae. Total RNA was extracted from infected Vero cells using a protocol preserving its integrity, as observed by Bioanalyzer (Agilent) profiles. By a subtractive hybridization method, the samples were subsequently depleted of eukaryotic RNA that represents up to 90% of the whole extract and that hampers fluorochrome labeling of rickettsial nucleic acids. To obtain the amount of material required for microarray hybridization, the bacterial RNA was then amplified using random primers. Hybridizations were carried out on microarrays specific for R. conorii but containing a limited number of selected targets. Our results show that this method yielded reproducible signals. Transcriptional changes observed following exposure of R. conorii to a nutrient stress were verified by real-time quantitative PCR and by quantitative reverse transcription PCR starting from amplified cDNA and total RNA as templates, respectively. We conclude that this approach has great potential for the study of mechanisms behind the virulence and intracellular survival of members of the genus Rickettsia.

Genetics-squared: combining host and pathogen genetics in the analysis of innate immunity and bacterial virulence

The interaction of bacterial pathogens with their hosts' innate immune systems can be extremely complex and is often difficult to disentangle experimentally. Using mouse models of bacterial infections, several laboratories have successfully applied genetic approaches to identify novel host genes required for innate immune defense. In addition, a variety of creative bacterial genetic schemes have been developed to identify key bacterial genes involved in triggering or evading host immunity. In cases where both the host and pathogen are amenable to genetic manipulation, a combination of host and pathogen genetic approaches can be used. Focusing on bacterial infections of mice, this review summarizes the benefits and limitations of applying genetic analysis to the study of host-pathogen interactions. In particular, we consider how prokaryotic and eukaryotic genetic strategies can be combined, or "squared," to yield new insights in host-pathogen biology.

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.