Bacterial artificial chromosome fingerprint arrays for the differentiation of transcriptomic differences in mycobacteria

Although microarray technology has become more widespread as a discovery tool for bacterial pathogenesis, it remains a method available only to laboratories with access to expensive equipment and costly analysis software. Mycobacterium tuberculosis, the causative agent for tuberculosis (TB), afflicts one-third of the global population, and kills between 2 and 3 million people per year. While the majority of cases of TB occur in developing areas of the world, facilities in these regions may not be able to support microarray analysis. Additionally, a major limitation of microarrays is that only genes on the array are being assayed. With acquired virulence and drug resistance in microbes, a method less dependent on a predetermined list of gene targets is advantageous. We present a method of expression analysis based on bacterial artificial chromosomes (BACs) that can be applied with standard laboratory equipment and free analysis software. This technique, bacterial artificial chromosome fingerprint arrays (BACFA), was developed and utilised to identify expression differences between intracellular strains of M. tuberculosis, one virulent (H37Rv) and one attenuated (H37Ra). Southern blots of restriction-enzyme digested BAC fragments were sequentially hybridised with strain-specific cDNA probes to generate expression profiles that were used to isolate expression differences in broth grown and intracellular bacteria. Repeat comparisons of intracellular profiles via BACFA identified genomic regions differentially expressed by the two strains. Quantitative real-time PCR was used to assess the genes located in these fragments in order to confirm or deny the differential regulation of genes. In total, we identified six genes that were differentially regulated between strains inside the host cell (pks2, aceE, Rv1571, and frdBCD). We report that BACFA is an effective technique in the expression analysis of bacteria and can be considered complementary to the high-throughput analysis offered by microarrays.     

Functional and comparative genomics of pathogenic bacteria

Microarray expression profiling and the development of data-mining tools and new statistical instruments affords an unprecedented opportunity for the genome-scale study of bacterial pathogenicity. Expression profiles obtained from bacteria grown in media simulating host microenvironments yield a portrait of interacting metabolic pathways and multistage developmental programs and disclose regulatory networks. The analysis of closely related strains and species by microarray-based comparative genomics provides a measure of genetic variability within natural populations and identifies crucial differences between pathogen and commensal. In the near future, the combined use of bacterial and host microarrays to study the same infected tissue will reveal the host-pathogen dialogue in a gene-by-gene and site- and time-specific manner. This review discusses the use of microarray-based expression profiling to identify genes of pathogenic bacteria that are differentially regulated in response to host-specific signals. Additionally, the review describes the application of microarray methods to disclose differences in gene content between taxonomically related strains that vary with respect to pathogenic phenotype.     

Signature patterns revealed by microarray analyses of mice infected with influenza virus A and Streptococcus pneumoniae

We used cDNA microarrays to identify differentially expressed genes in mice in response to infections with influenza virus A/PR/8/34 (H1N1) and Streptococcus pneumoniae. Expression microarray analysis showed up-regulation and down-regulation of many genes involved in the defense, inflammatory response and intracellular signaling pathways including chemokine, apoptosis, MAPK, Notch, Jak-STAT, T-cell receptor and complement and coagulation cascades. We have revealed signature patterns of gene expression in mice infected with two different classes of pathogens: influenza virus A and S. pneumoniae. Quantitative real-time RT-PCR results confirmed microarray results for most of the genes tested. These studies document clear differences in gene expression profiles between mice infected with influenza virus A and S. pneumoniae. Identification of genes that are differentially expressed after respiratory infections can provide insights into the mechanisms by which the host interacts with different pathogens, useful information about stage of diseases and selection of suitable targets for early diagnosis and treatments. The advantage of this novel approach is that the detection of pathogens is based on the differences in host gene expression profiles in response to different pathogens instead of detecting pathogens directly.     

Screening and comparison of differentially expressed genes between one MDR-TB strain and the virulent M. tuberculosis H37Rv

To screen and compare the differentially expressed genes between one MDR-TB strain separated from one child patient and the virulent Mycobacterium tuberculosis H37Rv, suppression subtractive hybridization (SSH) technology was used to build a library of cDNAs that were differentially expressed in the MDR and H37Rv. From this cDNA library, genes that were expressed in the MDR-TB but not in the H37Rv were selected for gene sequencing and homology analysis; 113 positive clones were obtained, their cDNA fragments were sequenced, and homology analysis was performed. Four novel sequences were identified. The results provide a partial list of genes differentially expressed in MDR-TB and four novel genes were found. Identification of these genes may contribute to our understanding of MDR-TB development.     

Rv2190c, an NlpC/P60 Family Protein, Is Required for Full Virulence of Mycobacterium tuberculosis

Mycobacterium tuberculosis, the etiologic agent of tuberculosis (TB) possesses at least five genes predicted to encode proteins with NlpC/P60 hydrolase domains, including the relatively uncharacterized Rv2190c. As NlpC/P60 domain-containing proteins are associated with diverse roles in bacterial physiology, our objective was to characterize Rv2190c in M. tuberculosis growth and virulence. Our data indicate that lack of Rv2190c is associated with impaired growth, both in vitro and during an in vivo mouse model of TB. These growth defects are associated with altered colony morphology and phthiocerol dimycocerosate levels, indicating that Rv2190c is involved in cell wall maintenance and composition. In addition, we have demonstrated that Rv2190c is expressed during active growth phase and that its protein product is immunogenic during infection. Our findings have significant implications, both for better understanding the role of Rv2190c in M. tuberculosis biology and also for translational developments.     

Deregulated microRNAs in CD4+ T cells from individuals with latent tuberculosis versus active tuberculosis

The mechanisms of latent tuberculosis (TB) infection remain elusive. Roles of microRNA (miRNA) have been highlighted in pathogen–host interactions recently. To identify miRNAs involved in the immune response to TB, expression profiles of miRNAs in CD4⁺ T cells from patients with latent TB, active TB and healthy controls were investigated by microarray assay and validated by RT‐qPCR. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis were used to analyse the significant functions and involvement in signalling pathways of the differentially expressed miRNAs. To identify potential target genes for miR‐29, interferon‐γ (IFN‐γ) mRNA expression was measured by RT‐qPCR. Our results showed that 27 miRNAs were deregulated among the three groups. RT‐qPCR results were generally consistent with the microarray data. We observed an inverse correlation between miR‐29 level and IFN‐γ mRNA expression in CD4⁺ T cells. GO and KEGG pathway analysis showed that the possible target genes of deregulated miRNAs were significantly enriched in mitogen‐activated protein kinase signalling pathway, focal adhesion and extracellular matrix receptor interaction, which might be involved in the transition from latent to active TB. In all, for the first time, our study revealed that some miRNAs in CD4⁺ T cells were altered in latent and active TB. Function and pathway analysis highlighted the possible involvement of miRNA‐deregulated mRNAs in TB. The study might help to improve understanding of the relationship between miRNAs in CD4⁺ T cells and TB, and laid an important foundation for further identification of the underlying mechanisms of latent TB infection and its reactivation.     

Novel serological biomarker panel using protein microarray can distinguish active TB from latent TB infection

BackgroundRapid laboratory technologies which can effectively distinguish active tuberculosis (ATB) from controls and latent tuberculosis infection (LTBI) are lacked.The objective of this study is to explore MTB biomarkers in serum that can distinguish ATB from LTBI.MethodsWe constructed a tuberculosis protein microarray containing 64 MTB associated antigens. We then used this microarray to screen 180 serum samples, from patients with ATB and LTBI, and healthy volunteer controls. Both SAM (Significance analysis of microarrays) and ROC curve analysis were used to identify the differentially recognized biomarkers between groups. Extra 300 serum samples from patients with ATB and LTBI, and healthy volunteer controls were employed to validate the identified biomarkers using ELISA-based method.ResultsAccording to the results, the best biomarker combinations of 4 proteins (Rv1860, RV3881c, Rv2031c and Rv3803c) were selected. The biomarker panel containing these 4 proteins has reached a sensitivity of 93.3% and specificity of 97.7% for distinguishing ATB from LTBI, and a sensitivity of 86% and specificity of 97.6% for distinguishing ATB from HC.ConclusionThe biomarker combination in this study has high sensitivity and specificity in distinguishing ATB from LTBI, suggesting it is worthy for further validation in more clinical samples.     

Requirement of gene fadD33 for the growth of Mycobacterium tuberculosis in a hepatocyte cell line

Gene fadD33 of Mycobacterium tuberculosis, one of the 36 homologues of gene fadD of Escherichia coli identified in the M. tuberculosis genome, predictively encodes an acyl-CoA synthase, an enzyme involved in fatty acids metabolism. The gene is underexpressed in the attenuated strain M. tuberculosis H37Ra relative to virulent H37Rv and plays a role in M. tuberculosis virulence in BALB/c mice by supporting mycobacterial replication in the liver. In the present paper, we investigated the role of fadD33 expression in bacterial growth within the hepatocyte cell line HepG2, as well as in human monocyte-derived THP-1 cells and peripheral blood mononuclear cells. M. tuberculosis H37Rv proved able to grow within HepG2 cells, while the intracellular replication of M. tuberculosis H37Ra was markedly impaired; complementation of strain H37Ra with gene fadD33 restored its replication to the levels of H37Rv. Moreover, disruption of gene fadD33 by allelic exchange mutagenesis reduced the intracellular growth of M. tuberculosis H37Rv, and complementation of the fadD33-disrupted mutant with gene fadD33 restored bacterial replication. Conversely, fadD33 expression proved unable to influence M. tuberculosis growth in human phagocytes, as fadD33-disrupted M. tuberculosis H37Rv mutant, as well as fadD33-complemented M. tuberculosis H37Ra, grew within THP-1 cells and peripheral monocytes basically at the same rates as parent H37Rv and H37Ra strains. The results of these experiments indicate that gene fadD33 expression confers growth advantage to M. tuberculosis in immortalized hepatocytes, but not in macrophages, thus emphasizing the importance of fadD33 in liver-specific replication of M. tuberculosis.     

Microarrays for public health: genomic epidemiology of tuberculosis

In response to a large local school-based outbreak of tuberculosis, we have been evaluating the utility of microarray bacterial genomic analysis in outbreak management. After initial comparison of the isolate from the index case with Mycobacterium tuberculosis H37Rv, it was possible to design robust PCRs directed towards strain-specific deletions. Rapid PCR analysis of isolates proved valuable in determining whether or not other isolates were compatible with the outbreak strain and further microarray studies revealed genetic markers that could be used to discriminate between locally circulating strains.We suggest that this approach forms the basis for developing rapid local genotyping schemes applicable to M. tuberculosis and that application to other pathogens warrants consideration.     

Identification of Tuberculosis Susceptibility Genes with Human Macrophage Gene Expression Profiles

Although host genetics influences susceptibility to tuberculosis (TB), few genes determining disease outcome have been identified. We hypothesized that macrophages from individuals with different clinical manifestations of Mycobacterium tuberculosis (Mtb) infection would have distinct gene expression profiles and that polymorphisms in these genes may also be associated with susceptibility to TB. We measured gene expression levels of >38,500 genes from ex vivo Mtb-stimulated macrophages in 12 subjects with 3 clinical phenotypes: latent, pulmonary, and meningeal TB (n = 4 per group). After identifying differentially expressed genes, we confirmed these results in 34 additional subjects by real-time PCR. We also used a case-control study design to examine whether polymorphisms in differentially regulated genes were associated with susceptibility to these different clinical forms of TB. We compared gene expression profiles in Mtb-stimulated and unstimulated macrophages and identified 1,608 and 199 genes that were differentially expressed by >2- and >5-fold, respectively. In an independent sample set of 34 individuals and a subset of highly regulated genes, 90% of the microarray results were confirmed by RT-PCR, including expression levels of CCL1, which distinguished the 3 clinical groups. Furthermore, 6 single nucleotide polymorphisms (SNPs) in CCL1 were found to be associated with TB in a case-control genetic association study with 273 TB cases and 188 controls. To our knowledge, this is the first identification of CCL1 as a gene involved in host susceptibility to TB and the first study to combine microarray and DNA polymorphism studies to identify genes associated with TB susceptibility. These results suggest that genome-wide studies can provide an unbiased method to identify critical macrophage response genes that are associated with different clinical outcomes and that variation in innate immune response genes regulate susceptibility to TB.     

Rv1818c-encoded PE_PGRS protein of Mycobacterium tuberculosis is surface exposed and influences bacterial cell structure

Identification of the novel PE multigene family was an unexpected finding of the genomic sequencing of Mycobacterium tuberculosis. Presently, the biological role of the PE and PE_PGRS proteins encoded by this unique family of mycobacterial genes remains unknown. In this report, a representative PE_PGRS gene (Rv1818c/PE_PGRS33) was selected to investigate the role of these proteins. Cell fractionation studies and fluorescence analysis of recombinant strains of Mycobacterium smegmatis and M. tuberculosis expressing green fluorescent protein (GFP)-tagged proteins indicated that the Rv1818c gene product localized in the mycobacterial cell wall, mostly at the bacterial cell poles, where it is exposed to the extracellular milieu. Further analysis of this PE_PGRS protein showed that the PE domain is necessary for subcellular localization. In addition, the PGRS domain, but not PE, affects bacterial shape and colony morphology when Rv1818c is overexpressed in M. smegmatis and M. tuberculosis. Taken together, the results indicate that PE_PGRS and PE proteins can be associated with the mycobacterial cell wall and influence cellular structure as well as the formation of mycobacterial colonies. Regulated expression of PE genes could have implications for the survival and pathogenesis of mycobacteria within the human host and in other environmental niches.     

3′-End Sequencing for Expression Quantification (3SEQ) from Archival Tumor Samples

Gene expression microarrays are the most widely used technique for genome-wide expression profiling. However, microarrays do not perform well on formalin fixed paraffin embedded tissue (FFPET). Consequently, microarrays cannot be effectively utilized to perform gene expression profiling on the vast majority of archival tumor samples. To address this limitation of gene expression microarrays, we designed a novel procedure (3′-end sequencing for expression quantification (3SEQ)) for gene expression profiling from FFPET using next-generation sequencing. We performed gene expression profiling by 3SEQ and microarray on both frozen tissue and FFPET from two soft tissue tumors (desmoid type fibromatosis (DTF) and solitary fibrous tumor (SFT)) (total n = 23 samples, which were each profiled by at least one of the four platform-tissue preparation combinations). Analysis of 3SEQ data revealed many genes differentially expressed between the tumor types (FDR<0.01) on both the frozen tissue (∼9.6K genes) and FFPET (∼8.1K genes). Analysis of microarray data from frozen tissue revealed fewer differentially expressed genes (∼4.64K), and analysis of microarray data on FFPET revealed very few (69) differentially expressed genes. Functional gene set analysis of 3SEQ data from both frozen tissue and FFPET identified biological pathways known to be important in DTF and SFT pathogenesis and suggested several additional candidate oncogenic pathways in these tumors. These findings demonstrate that 3SEQ is an effective technique for gene expression profiling from archival tumor samples and may facilitate significant advances in translational cancer research.