Comparing genomes within the species Mycobacterium tuberculosis

Diversification of gene content in the Mycobacterium tuberculosis complex is determined by phylogenetic and ecological signatures

We analyzed the pan-genome and gene content modulation of the most diverse genome data set of the Mycobacterium tuberculosis complex (MTBC) gathered to date. The closed pan-genome of the MTBC was characterized by reduced accessory and strain-specific genomes, compatible with its clonal nature. However, significantly fewer gene families were shared between MTBC genomes as their phylogenetic distance increased. This effect was only observed in inter-species comparisons, not within-species, which suggests that species-specific ecological characteristics are associated with changes in gene content. Gene loss, resulting from genomic deletions and pseudogenization, was found to drive the variation in gene content. This gene erosion differed among MTBC species and lineages, even within M. tuberculosis, where L2 showed more gene loss than L4. We also show that phylogenetic proximity is not always a good proxy for gene content relatedness in the MTBC, as the gene repertoire of Mycobacterium africanum L6 deviated from its expected phylogenetic niche conservatism. Gene disruptions of virulence factors, represented by pseudogene annotations, are mostly not conserved, being poor predictors of MTBC ecotypes. Each MTBC ecotype carries its own accessory genome, likely influenced by distinct selective pressures such as host and geography. It is important to investigate how gene loss confer new adaptive traits to MTBC strains; the detected heterogeneous gene loss poses a significant challenge in elucidating genetic factors responsible for the diverse phenotypes observed in the MTBC. By detailing specific gene losses, our study serves as a resource for researchers studying the MTBC phenotypes and their immune evasion strategies.IMPORTANCEIn this study, we analyzed the gene content of different ecotypes of the Mycobacterium tuberculosis complex (MTBC), the pathogens of tuberculosis. We found that changes in their gene content are associated with their ecological features, such as host preference. Gene loss was identified as the primary driver of these changes, which can vary even among different strains of the same ecotype. Our study also revealed that the gene content relatedness of these bacteria does not always mirror their evolutionary relationships. In addition, some genes of virulence can be variably lost among strains of the same MTBC ecotype, likely helping them to evade the immune system. Overall, our study highlights the importance of understanding how gene loss can lead to new adaptations in these bacteria and how different selective pressures may influence their genetic makeup.

Mapping the phylogeny and lineage history of geographically distinct BCG vaccine strains

The bacillus Calmette-Guérin (BCG) vaccine has been in use for prevention of tuberculosis for over a century. It remains the only widely available tuberculosis vaccine and its protective efficacy has varied across geographical regions. Since it was developed, the BCG vaccine strain has been shared across different laboratories around the world, where use of differing culture methods has resulted in genetically distinct strains over time. Whilst differing BCG vaccine efficacy around the world is well documented, and the reasons for this may be multifactorial, it has been hypothesized that genetic differences in BCG vaccine strains contribute to this variation. Isolates from an historic archive of lyophilized BCG strains were regrown, DNA was extracted and then whole-genome sequenced using Oxford Nanopore Technologies. The resulting whole-genome data were plotted on a phylogenetic tree and analysed to identify the presence or absence of regions of difference (RDs) and single-nucleotide polymorphisms (SNPs) relating to virulence, growth and cell wall structure. Of 50 strains available, 36 were revived in culture and 39 were sequenced. Morphology differed between the strains distributed before and after 1934. There was phylogenetic association amongst certain geographically classified strains, most notably BCG-Russia, BCG-Japan and BCG-Danish. RD2, RD171 and RD713 deletions were associated with late strains (seeded after 1927). When mapped to BCG-Pasteur 1172, the SNPs in sigK, plaA, mmaA3 and eccC5 were associated with early strains. Whilst BCG-Russia, BCG-Japan and BCG-Danish showed strong geographical isolate clustering, the late strains, including BCG-Pasteur, showed more variation. A wide range of SNPs were seen within geographically classified strains, and as much intra-strain variation as between-strain variation was seen. The date of distribution from the original Pasteur laboratory (early pre-1927 or late post-1927) gave the strongest association with genetic differences in regions of difference and virulence-related SNPs, which agrees with the previous literature.

Using Omics to Study Leprosy, Tuberculosis, and Other Mycobacterial Diseases

Mycobacteria are members of the Actinomycetales order, and they are classified into one family, Mycobacteriaceae. More than 20 mycobacterial species cause disease in humans. The Mycobacterium group, called the Mycobacterium tuberculosis complex (MTBC), has nine closely related species that cause tuberculosis in animals and humans. TB can be detected worldwide and one-fourth of the world’s population is contaminated with tuberculosis. According to the WHO, about two million dies from it, and more than nine million people are newly infected with TB each year. Mycobacterium tuberculosis (M. tuberculosis) is the most potential causative agent of tuberculosis and prompts enormous mortality and morbidity worldwide due to the incompletely understood pathogenesis of human tuberculosis. Moreover, modern diagnostic approaches for human tuberculosis are inefficient and have many lacks, while MTBC species can modulate host immune response and escape host immune attacks to sustain in the human body. “Multi-omics” strategies such as genomics, transcriptomics, proteomics, metabolomics, and deep sequencing technologies could be a comprehensive strategy to investigate the pathogenesis of mycobacterial species in humans and offer significant discovery to find out biomarkers at the early stage of disease in the host. Thus, in this review, we attempt to understand an overview of the mission of “omics” approaches in mycobacterial pathogenesis, including tuberculosis, leprosy, and other mycobacterial diseases.

A Comprehensive Map of Mycobacterium tuberculosis Complex Regions of Difference

Mycobacterium tuberculosis complex (MTBC) species are classic examples of genetically monomorphic microorganisms due to their low genetic variability. Whole-genome sequencing made it possible to describe both the main species within the complex and M. tuberculosis lineages and sublineages. This differentiation is based on single nucleotide polymorphisms (SNPs) and large sequence polymorphisms in the so-called regions of difference (RDs). Although a number of studies have been performed to elucidate RD localizations, their distribution among MTBC species, and their role in the bacterial life cycle, there are some inconsistencies and ambiguities in the localization of RDs in different members of the complex. To address this issue, we conducted a thorough search for all possible deletions in the WGS data collection comprising 721 samples representing the full MTBC diversity. Discovered deletions were compared with a list of all previously described RDs. As with the SNP-based analysis, we confirmed the specificities of 79 regions at the species, lineage, or sublineage level, 17 of which are described for the first time. We also present RDscan (https://github.com/dbespiatykh/RDscan), an open-source workflow, which detects deletions from short-read sequencing data and correlates the results with high-specificity RDs, curated in this study. Testing of the workflow on a collection comprising ∼7,000 samples showed a high specificity of the found RDs. This study provides novel details that can contribute to a better understanding of the species differentiation within the MTBC and can help to determine how individual clusters evolve within various MTBC species.

IMPORTANCE Reductive genome evolution is one of the most important and intriguing adaptation strategies of different living organisms to their environment. Mycobacterium offers several notorious examples of either naturally reduced (Mycobacterium leprae) or laboratory-reduced (Mycobacterium bovis BCG) genomes. Mycobacterium tuberculosis complex has its phylogeny unambiguously framed by large sequence polymorphisms that present unidirectional unique event changes. In the present study, we curated all known regions of difference and analyzed both Mycobacterium tuberculosis and animal-adapted MTBC species. For 79 loci, we have shown a relationship with phylogenetic units, which can serve as a marker for diagnosing or studying biological effects. Moreover, intersections were found for some loci, which may indicate the nonrandomness of these processes and the involvement of these regions in the adaptation of bacteria to external conditions.

Mycobacterium tuberculosis complex lineage 5 exhibits high levels of within-lineage genomic diversity and differing gene content compared to the type strain H37Rv

Pathogens of the Mycobacterium tuberculosis complex (MTBC) are considered to be monomorphic, with little gene content variation between strains. Nevertheless, several genotypic and phenotypic factors separate strains of the different MTBC lineages (L), especially L5 and L6 (traditionally termed Mycobacterium africanum) strains, from each other. However, this genome variability and gene content, especially of L5 strains, has not been fully explored and may be important for pathobiology and current approaches for genomic analysis of MTBC strains, including transmission studies. By comparing the genomes of 355 L5 clinical strains (including 3 complete genomes and 352 Illumina whole-genome sequenced isolates) to each other and to H37Rv, we identified multiple genes that were differentially present or absent between H37Rv and L5 strains. Additionally, considerable gene content variability was found across L5 strains, including a split in the L5.3 sub-lineage into L5.3.1 and L5.3.2. These gene content differences had a small knock-on effect on transmission cluster estimation, with clustering rates influenced by the selected reference genome, and with potential overestimation of recent transmission when using H37Rv as the reference genome. We conclude that full capture of the gene diversity, especially high-resolution outbreak analysis, requires a variation of the single H37Rv-centric reference genome mapping approach currently used in most whole-genome sequencing data analysis pipelines. Moreover, the high within-lineage gene content variability suggests that the pan-genome of M. tuberculosis is at least several kilobases larger than previously thought, implying that a concatenated or reference-free genome assembly (de novo) approach may be needed for particular questions.

Human major infections: Tuberculosis, treponematoses, leprosy-A paleopathological perspective of their evolution

The key to evolution is reproduction. Pathogens can either kill the human host or can invade the host without causing death, thus ensuring their own survival, reproduction and spread. Tuberculosis, treponematoses and leprosy are widespread chronic infectious diseases whereby the host is not immediately killed. These diseases are examples of the co-evolution of host and pathogen. They can be well studied as the paleopathological record is extensive, spanning over 200 human generations. The paleopathology of each disease has been well documented in the form of published synthetic analyses recording each known case and case frequencies in the samples they were derived from. Here the data from these synthetic analyses were re-analysed to show changes in the prevalence of each disease over time. A total of 69,379 skeletons are included in this study. There was ultimately a decline in the prevalence of each disease over time, this decline was statistically significant (Chi-squared, p<0.001). A trend may start with the increase in the disease’s prevalence before the prevalence declines, in tuberculosis the decline is monotonic. Increase in skeletal changes resulting from the respective diseases appears in the initial period of host-disease contact, followed by a decline resulting from co-adaptation that is mutually beneficial for the disease (spread and maintenance of pathogen) and host (less pathological reactions to the infection). Eventually either the host may become immune or tolerant, or the pathogen tends to be commensalic rather than parasitic.

Transmission phenotype of Mycobacterium tuberculosis strains is mechanistically linked to induction of distinct pulmonary pathology

In a study of household contacts (HHC), households were categorized into High (HT) and Low (LT) transmission groups based on the proportion of HHC with a positive tuberculin skin test. The Mycobacterium tuberculosis (Mtb) strains from HT and LT index cases of the households were designated Mtb-HT and Mtb-LT, respectively. We found that C3HeB/FeJ mice infected with Mtb-LT strains exhibited significantly higher bacterial burden compared to Mtb-HT strains and also developed diffused inflammatory lung pathology. In stark contrast, a significant number of mice infected with Mtb-HT strains developed caseating granulomas, a lesion type with high potential to cavitate. None of the Mtb-HT infected animals developed diffused inflammatory lung pathology. A link was observed between increased in vitro replication of Mtb-LT strains and their ability to induce significantly high lipid droplet formation in macrophages. These results support that distinct early interactions of Mtb-HT and Mtb-LT strains with macrophages and subsequent differential trajectories in pathological disease may be the mechanism underlying their transmission potential.

Identification of Mycobacterial Ribosomal Proteins as Targets for CD4+ T Cells That Enhance Protective Immunity in Tuberculosis

Mycobacterium tuberculosis remains a threat to global health, and a more efficacious vaccine is needed to prevent disease caused by M. tuberculosis We previously reported that the mycobacterial ribosome is a major target of CD4⁺ T cells in mice immunized with a genetically modified Mycobacterium smegmatis strain (IKEPLUS) but not in mice immunized with Mycobacterium bovis BCG. Two specific ribosomal proteins, RplJ and RpsA, were identified as cross-reactive targets of M. tuberculosis, but the breadth of the CD4⁺ T cell response to M. tuberculosis ribosomes was not determined. In the present study, a library of M. tuberculosis ribosomal proteins and in silico-predicted peptide libraries were used to screen CD4⁺ T cell responses in IKEPLUS-immunized mice. This identified 24 out of 57 M. tuberculosis ribosomal proteins distributed over both large and small ribosome subunits as specific CD4⁺ T cell targets. Although BCG did not induce detectable responses against ribosomal proteins or peptide epitopes, the M. tuberculosis ribosomal protein RplJ produced a robust and multifunctional Th1-like CD4⁺ T cell population when administered as a booster vaccine to previously BCG-primed mice. Boosting of BCG-primed immunity with the M. tuberculosis RplJ protein led to significantly reduced lung pathology compared to that in BCG-immunized animals and reductions in the bacterial burdens in the mediastinal lymph node compared to those in naive and standard BCG-vaccinated mice. These results identify the mycobacterial ribosome as a potential source of cryptic or subdominant antigenic targets of protective CD4⁺ T cell responses and suggest that supplementing BCG with ribosomal antigens may enhance protective vaccination against M. tuberculosis.

RD5-mediated lack of PE_PGRS and PPE-MPTR export in BCG vaccine strains results in strong reduction of antigenic repertoire but little impact on protection

Tuberculosis is the deadliest infectious disease worldwide. Although the BCG vaccine is widely used, it does not efficiently protect against pulmonary tuberculosis and an improved tuberculosis vaccine is therefore urgently needed. Mycobacterium tuberculosis uses different ESX/Type VII secretion (T7S) systems to transport proteins important for virulence and host immune responses. We recently reported that secretion of T7S substrates belonging to the mycobacteria-specific Pro-Glu (PE) and Pro-Pro-Glu (PPE) proteins of the PGRS (polymorphic GC-rich sequences) and MPTR (major polymorphic tandem repeat) subfamilies required both a functional ESX-5 system and a functional PPE38/71 protein for secretion. Inactivation of ppe38/71 and the resulting loss of PE_PGRS/PPE-MPTR secretion were linked to increased virulence of M. tuberculosis strains. Here, we show that a predicted total of 89 PE_PGRS/PPE-MPTR surface proteins are not exported by certain animal-adapted strains of the M. tuberculosis complex including M. bovis. This Δppe38/71-associated secretion defect therefore also occurs in the M. bovis-derived tuberculosis vaccine BCG and could be partially restored by introduction of the M. tuberculosis ppe38-locus. Epitope mapping of the PPE-MPTR protein PPE10, further allowed us to monitor T-cell responses in splenocytes from BCG/M. tuberculosis immunized mice, confirming the dependence of PPE10-specific immune-induction on ESX-5/PPE38-mediated secretion. Restoration of PE_PGRS/PPE-MPTR secretion in recombinant BCG neither altered global antigenic presentation or activation of innate immune cells, nor protective efficacy in two different mouse vaccination-infection models. This unexpected finding stimulates a reassessment of the immunomodulatory properties of PE_PGRS/PPE-MPTR proteins, some of which are contained in vaccine formulations currently in clinical evaluation.

One of the major findings of the pioneering Mycobacterium tuberculosis H37Rv genome sequencing project was the identification of the highly abundant PE and PPE proteins, named after their N-terminal motifs Pro–Glu (PE) or Pro–Pro–Glu (PPE). Within the 20 years of research since then, many claims were made that PE/PPE proteins, including the two large subgroups encoded by repetitive sequences with very high GC content (PE_PGRS and PPE-MPTR families), are exported to the bacterial surface or beyond, and show broad immunomodulatory impact on host-pathogen interaction. We thus screened strains from different branches of the M. tuberculosis complex, including the attenuated Mycobacterium bovis BCG vaccine strains, for their capacity to export PE_PGRS/PPE-MPTR proteins. Strikingly, we found that BCG strains were unable to export the plethora of PE_PGRS/PPE-MPTR proteins due to the absence of the region of difference RD5, which in M. tuberculosis encodes PPE38, required for PE_PGRS/PPE-MPTR export. Surprisingly, the restoration of PE_PGRS/PPE-MPTR export by genetic complementation in recombinant BCG did not result in immunomodulatory changes or altered protection in mouse models. Our results thus put into perspective the numerous reports on virulence-associated immunomodulatory impact of individual PE_PGRS and PPE-MPTR proteins and open novel questions on their biological function(s).

Complete Genome Sequencing of Mycobacterium bovis SP38 and Comparative Genomics of Mycobacterium bovis and M. tuberculosis Strains

Mycobacterium bovis causes bovine tuberculosis and is the main organism responsible for zoonotic tuberculosis in humans. We performed the sequencing, assembly and annotation of a Brazilian strain of M. bovis named SP38, and performed comparative genomics of M. bovis genomes deposited in GenBank. M. bovis SP38 has a traditional tuberculous mycobacterium genome of 4,347,648 bp, with 65.5% GC, and 4,216 genes. The majority of CDSs (2,805, 69.3%) have predictive function, while 1,206 (30.07%) are hypothetical. For comparative analysis, 31 M. bovis, 32 M. bovis BCG, and 23 Mycobacterium tuberculosis genomes available in GenBank were selected. M. bovis RDs (regions of difference) and Clonal Complexes (CC) were identified in silico. Genome dynamics of bacterial groups were analyzed by gene orthology and polymorphic sites identification. M. bovis polymorphic sites were used to construct a phylogenetic tree. Our RD analyses resulted in the exclusion of three genomes, mistakenly annotated as virulent M. bovis. M. bovis SP38 along with strain 35 represent the first report of CC European 2 in Brazil, whereas two other M. bovis strains failed to be classified within current CC. Results of M. bovis orthologous genes analysis suggest a process of genome remodeling through genomic decay and gene duplication. Quantification, pairwise comparisons and distribution analyses of polymorphic sites demonstrate greater genetic variability of M. tuberculosis when compared to M. bovis and M. bovis BCG (p ≤ 0.05), indicating that currently defined M. tuberculosis lineages are more genetically diverse than M. bovis CC and animal-adapted MTC (M. tuberculosis Complex) species. As expected, polymorphic sites annotation shows that M. bovis BCG are subjected to different evolutionary pressures when compared to virulent mycobacteria. Lastly, M. bovis phylogeny indicates that polymorphic sites may be used as markers of M. bovis lineages in association with CC. Our findings highlight the need to better understand host-pathogen co-evolution in genetically homogeneous and/or diverse host populations, considering the fact that M. bovis has a broader host range when compared to M. tuberculosis. Also, the identification of M. bovis genomes not classified within CC indicates that the diversity of M. bovis lineages may be larger than previously thought or that current classification should be reviewed.

Modelling Prokaryote Gene Content

The patchy distribution of genes across the prokaryotes may be caused by multiple gene losses or lateral transfer. Probabilistic models of gene gain and loss are needed to distinguish between these possibilities. Existing models allow only single genes to be gained and lost, despite the empirical evidence for multi-gene events. We compare birth-death models (currently the only widely-used models, in which only one gene can be gained or lost at a time) to blocks models (allowing gain and loss of multiple genes within a family). We analyze two pairs of genomes: two E. coli strains, and the distantly-related Archaeoglobus fulgidus (archaea) and Bacillus subtilis (gram positive bacteria). Blocks models describe the data much better than birth-death models. Our models suggest that lateral transfers of multiple genes from the same family are rare (although transfers of single genes are probably common). For both pairs, the estimated median time that a gene will remain in the genome is not much greater than the time separating the common ancestors of the archaea and bacteria. Deep phylogenetic reconstruction from sequence data will therefore depend on choosing genes likely to remain in the genome for a long time. Phylogenies based on the blocks model are more biologically plausible than phylogenies based on the birth-death model.

Human leukocyte antigen A, B and Hepatitis B infection outcome: A meta-analysis

AIM

To evaluate the association between HLA A, B with chronic Hepatitis B by comprehensive meta-analysis.

METHODS

We searched PubMed and Cochrane databases, identified relevant studies, evaluated these for quality by New Castle Ottawa scale (NOS) and further analyzed the qualified data sets. Heterogeneity analyses were performed by Cochrane's Q test and I² tests. Pooled Odds ratio (OR) & 95% Confidence Interval (CI) were obtained by fixed effects, using Mantel-Haenszel's method for homogenous studies, and by using DerSimonian and Laird's method for heterogenous studies. Publication bias was determined by the Beggs test and Eggers test and all tests were two tailed to evaluate their significance.

RESULTS

The meta-analyses on 1652 healthy controls and 659 Chronic Hepatitis B (CHB) patients from 8 studies from various continents revealed a HLA B*07 (p value of Odds ratio (pOR)=0.004; OR Fixed effects=0.480 with 95%CI 0.290-0.794) and B*58 (pOR=0.029; OR Fixed Effects=0.020 with 95%CI 0.381-0.949) associated protection for CHB. The identified HLA B*35 associated risk (pOR 0.009; OR Fixed effect 1.445; 95% confidence interval 1.094-1.907) however did not stand the test of random effect model.

CONCLUSION

While HLA B*07 and B*58 are protective against CHB. The HLA B*35 associated marginal risk need to be further validated in well-designed global study on larger cohorts, considering the population, ethnic, epidemiological and HLA diversity at the sequence level: these may throw further light to utilize these markers in predictive medicine.

Methodological and Clinical Aspects of the Molecular Epidemiology of Mycobacterium tuberculosis and Other Mycobacteria

Molecular typing has revolutionized epidemiological studies of infectious diseases, including those of a mycobacterial etiology. With the advent of fingerprinting techniques, many traditional concepts regarding transmission, infectivity, or pathogenicity of mycobacterial bacilli have been revisited, and their conventional interpretations have been challenged. Since the mid-1990s, when the first typing methods were introduced, a plethora of other modalities have been proposed. So-called molecular epidemiology has become an essential subdiscipline of modern mycobacteriology. It serves as a resource for understanding the key issues in the epidemiology of tuberculosis and other mycobacterial diseases. Among these issues are disclosing sources of infection, quantifying recent transmission, identifying transmission links, discerning reinfection from relapse, tracking the geographic distribution and clonal expansion of specific strains, and exploring the genetic mechanisms underlying specific phenotypic traits, including virulence, organ tropism, transmissibility, or drug resistance. Since genotyping continues to unravel the biology of mycobacteria, it offers enormous promise in the fight against and prevention of the diseases caused by these pathogens. In this review, molecular typing methods for Mycobacterium tuberculosis and nontuberculous mycobacteria elaborated over the last 2 decades are summarized. The relevance of these methods to the epidemiological investigation, diagnosis, evolution, and control of mycobacterial diseases is discussed.

Eukaryotic pathways targeted by the type III secretion system effector protein, BipC, involved in the intracellular lifecycle of Burkholderia pseudomallei

Burkholderia pseudomallei, the etiological agent for melioidosis, is known to secrete a type III secretion system (TTSS) protein into the host’s internal milieu. One of the TTSS effector protein, BipC, has been shown to play an important role in the B. pseudomallei pathogenesis. To identify the host response profile that was directly or indirectly regulated by this protein, genome-wide transcriptome approach was used to examine the gene expression profiles of infected mice. The transcriptome analysis of the liver and spleen revealed that a total of approximately 1,000 genes were transcriptionally affected by BipC. Genes involved in bacterial invasion, regulation of actin cytoskeleton, and MAPK signalling pathway were over-expressed and may be specifically regulated by BipC in vivo. These results suggest that BipC mainly targets pathways related to the cellular processes which could modulate the cellular trafficking processes. The host transcriptional response exhibited remarkable differences with and without the presence of the BipC protein. Overall, the detailed picture of this study provides new insights that BipC may have evolved to efficiently manipulate host-cell pathways which is crucial in the intracellular lifecycle of B. pseudomallei.

Identification of Mycobacterium tuberculosis in the cerebrospinal fluid of patients with meningitis using nested PCR

Tuberculous meningitis (TBM) is the most severe form of tuberculosis. It is caused by Mycobacterium tuberculosis (M. tuberculosis; MT) and it is very difficult to diagnose. The symptoms are similar to other infectious neurological diseases, such as neurocysticercosis, neuroborreliosis, or herpes viral infection. The aim of this study was to identify tuberculosis (TB) in cases of meningitis with clinical and laboratory evidence suggestive of TBM, and to confirm our findings with molecular tests for TB infection. We recruited patients with neurological symptoms who were examined at the neurology services of Hospitals of Instituto Mexicano del Seguro Social (IMSS) in Mexico City. A total of 144 consecutive patients with suggestive infectious meningitis were initially included; 94 cases of meningitis with clinical and laboratory evidence suggestive of TBM were included, but only 50 of these cases fulfilled the criteria for probable TBM. As the controls, we included 50 cases of meningitis with clinical and laboratory evidence suggestive of non-TBM. Cerebrospinal fluid (CSF) was collected from all 100 patients (cases and controls) and tested for TB by multiplex and nested PCR analyses. Nested PCR detected 0.1 fg of M. tuberculosis DNA. TB infection was confirmed with molecular tests in 49 patients from the 50 cases suggestive of TBM and in 1 of the 50 non-TBM cases. The analysis exhibited a sensitivity of 98.0%, a specificity of 92.0%, a positive predictive value of 88.0% and a negative predictive value of 98.0%. The use CSF for the analyses proved to be effective for the rapid diagnosis of TBM using a developed system of multiplex and nested PCR analyses in patients presenting neurological symptoms.

Not all the infected develop the disease - A "Lotus and Cactus" model

The immunogenetic dictum "not all the infected develop the disease" can best be explained by a "Lotus and Cactus" model. Lotuses grow in ponds and cacti in deserts: analogously, we can say that tubercle patient's lung (genetic makeup) functions as an ideal 'broth' for Mycobacterium tuberculosis (M.tb) germs to grow, but not the lungs of an endemic control. HLA association studies from Europe to Asia since 1983 till date, have shown a persistent HLA DR2 (15) association. Further, HLA DR2 and non-DR2 endemic controls showed disparate patterns of immune responses and gene expressions. The host and pathogen MHC diversities, Th1-Th2 paradigm and cytokine circuits all may play a crucial role in TB susceptibility. It is possible to decipher the protective immunity by controlling the known confounders - epidemiological, demographic, socio-biological and also host and pathogen diversities. This has become significant with our understanding on the 'out of Africa' migration and neolithic co-dispersal of M.tb with modern human. Divergence and expansion of various MHCs (eg HLA-DRB1*15, HLA-B*57) and non-MHC alleles in various continents might be responsible for the skewed transmission and distribution of the infectious diseases around the globe. The 'Lotus and Cactus' model proposed here exemplifies this. A holistic genetic epidemiology approach employing modern tools is the need of the hour to better understand infectious disease susceptibility.

Solid and Suspension Microarrays for Microbial Diagnostics

Advancements in molecular technologies have provided new platforms that are being increasingly adopted for use in the clinical microbiology laboratory. Among these, microarray methods are particularly well suited for diagnostics as they allow multiplexing, or the ability to test for multiple targets simultaneously from the same specimen. Microarray technologies commonly used for the detection and identification of microbial targets include solid-state microarrays, electronic microarrays and bead suspension microarrays. Microarray methods have been applied to microbial detection, genotyping and antimicrobial resistance gene detection. Microarrays can offer a panel approach to diagnose specific patient presentations, such as respiratory or gastrointestinal infections, and can discriminate isolates by genotype for tracking epidemiology and outbreak investigations. And, as more information has become available on specific genes and pathways involved in antimicrobial resistance, we are beginning to be able to predict susceptibility patterns based on sequence detection for particular organisms. With further advances in automated microarray processing methods and genotype-phenotype prediction algorithms, these tests will become even more useful as an adjunct or replacement for conventional antimicrobial susceptibility testing, allowing for more rapid selection of targeted therapy for infectious diseases.

Inter-individual differences in the gene content of human gut bacterial species

BACKGROUND

Gene content differences in human gut microbes can lead to inter-individual phenotypic variations such as digestive capacity. It is unclear whether gene content variation is caused by differences in microbial species composition or by the presence of different strains of the same species; the extent of gene content variation in the latter is unknown. Unlike pan-genome studies of cultivable strains, the use of metagenomic data can provide an unbiased view of structural variation of gut bacterial strains by measuring them in their natural habitats, the gut of each individual in this case, representing native boundaries between gut bacterial populations. We analyzed publicly available metagenomic data from fecal samples to characterize inter-individual variation in gut bacterial species.

RESULTS

A comparison of 11 abundant gut bacterial species showed that the gene content of strains from the same species differed, on average, by 13% between individuals. This number is based on gene deletions only and represents a lower limit, yet the variation is already in a similar range as observed between completely sequenced strains of cultivable species. We show that accessory genes that differ considerably between individuals can encode important functions, such as polysaccharide utilization and capsular polysaccharide synthesis loci.

CONCLUSION

Metagenomics can yield insights into gene content variation of strains in complex communities, which cannot be predicted by phylogenetic marker genes alone. The large degree of inter-individual variability in gene content implies that strain resolution must be considered in order to fully assess the functional potential of an individual's human gut microbiome.

Overview and phylogeny of Mycobacterium tuberculosis complex organisms: implications for diagnostics and legislation of bovine tuberculosis

Members of the Mycobacterium tuberculosis complex (MTBC) cause a serious disease with similar pathology, tuberculosis; in this review, bovine tuberculosis will be considered as disease caused by any member of the MTBC in bovids. Bovine tuberculosis is responsible for significant economic loss due to costly eradication programs and trade limitations and poses a threat to both endangered and protected species as well as to public health. We here give an overview on all members of the MTBC, focusing on their isolation from different animal hosts. We also review the recent advances made in elucidating the evolutionary and phylogenetic relationships of members of the MTBC. Because the nomenclature of the MTBC is controversial, its members have been considered species, subspecies or ecotypes, this review discusses the possible implications for diagnostics and the legal consequences of naming of new species.

Current methods in the molecular typing of Mycobacterium tuberculosis and other mycobacteria

In the epidemiology of tuberculosis (TB) and nontuberculous mycobacterial (NTM) diseases, as in all infectious diseases, the key issue is to define the source of infection and to disclose its routes of transmission and dissemination in the environment. For this to be accomplished, the ability of discerning and tracking individual Mycobacterium strains is of critical importance. Molecular typing methods have greatly improved our understanding of the biology of mycobacteria and provide powerful tools to combat the diseases caused by these pathogens. The utility of various typing methods depends on the Mycobacterium species under investigation as well as on the research question. For tuberculosis, different methods have different roles in phylogenetic analyses and person-to-person transmission studies. In NTM diseases, most investigations involve the search for environmental sources or phylogenetic relationships. Here, too, the type of setting determines which methodology is most suitable. Within this review, we summarize currently available molecular methods for strain typing of M. tuberculosis and some NTM species, most commonly associated with human disease. For the various methods, technical practicalities as well as discriminatory power and accomplishments are reviewed.

Molecular characterization of Mycobacterium tuberculosis isolates from North Indian patients with extrapulmonary tuberculosis

Genotypic studies are important to understand the molecular epidemiology and transmission routes of Mycobacterium tuberculosis. In the first and largest study from India, spoligotyping and 24 loci mycobacterial interspersed repetitive units (MIRU) were performed to find genetic profiles of 125 M. tuberculosis strains isolated from patients with extrapulmonary tuberculosis (EPTB) and their drug susceptibility test was performed using BACTEC-MGIT 960. Spoligotyping results were compared with the world Spoligotyping Database of Institute Pasteur de Guadeloupe (SpolDB4). The spoligotyping results showed that 110 (88%) displayed known patterns while 15 (12%) isolates had no matching database. Predominant spoligotypes belonged to CAS family (57.27%). The largest clade comprised of 38 isolates belonging to the CAS1_DEL lineage. Though there was no significant association between specific mycobacterial lineage and extrapulmonary site, a significantly high (p < 0.001) number of Beijing type isolates (28.6%) were isolated from bone and joint samples as compared to cerebrospinal fluid (5%). There was a significant association between Beijing family isolates and multi-drug-resistance, while all MANU genotypes were pan-drug sensitive. The CAS family lineage was most prevalent genotype in the EPTB cases in our population.

Genetic diversity in Mycobacterium tuberculosis

Recent years have witnessed an increased appreciation of the extent and relevance of strain-to-strain variation in Mycobacterium tuberculosis. This paradigm shift can largely be attributed to an improved understanding of the global population structure of this organism, and to the realisation that the various members of the M. tuberculosis complex (MTBC) harbour more genetic diversity than previously realised. Moreover, many studies using experimental models of infection have demonstrated that MTBC diversity translates into significant differences in immunogenecity and virulence . However, linking these experimental phenotypes to relevant clinical phenotypes has been difficult, and to date, largely unsuccessful. Nevertheless, emerging high-throughput technologies, in particular next-generation sequencing , offer new opportunities, and have already lead to important new insights. Given the complexity of the host-pathogen interaction in tuberculosis, systems approaches will be key to define the role of MTBC diversity in the fight against one of humankind's most important pathogens.

Molecular detection of nontuberculous mycobacteria: advantages and limits of a broad-range sequencing approach

The isolation of nontuberculous mycobacteria (NTM) from clinical specimens has become very common in recent years. Such organisms are typically environmental and occasionally pathogenic for humans and animals. Standard diagnosis of mycobacterial infections relies on direct examination and culture. However, molecular tools are now available which allow quicker and more accurate diagnosis. Detection of NTM can be performed directly from clinical samples, although identification is mostly carried out after isolation. Sequencing of genomic targets (such as 16S rRNA, ITS, rpoB or hsp65) allows accurate and rapid identification, but has some technical limitations. A brief summary of the molecular methods available for NTM identification and a discussion of the problems associated with the use of sequencing analysis together with a description of available algorithms for NTM identification are the major objectives of this review.

Use of mouse models to study the variability in virulence associated with specific genotypic lineages of Mycobacterium tuberculosis

The host response against Mycobacterium tuberculosis show a wide spectrum of clinical manifestations in those patients who fail to control the infection. The course of the infection and its epidemiological consequences depend upon a complex interplay of host, environmental and bacterial factors. Experimental animal models have helped to define the influence of bacterial genetic diversity on virulence and on the immune response that is induced. For this purpose, experimental animals such as mice, guinea pigs and rabbits have been infected with selected clinical isolates obtained from outbreaks or from clinical epidemiology settings. Here we review the contribution of mouse models to defining the variability in virulence and immune response in relation to mycobacterial genetic diversity. Low dose aerosol infection in C57Bl mice or high dose intratracheal infection in BALB/c mice have demonstrated wide variability in virulence and immune responses induced by different bacterial genotypes, and each genotype has different phenotypes, with high and low virulence variants. In general, these studies have shown that high prevalent strains from big clusters are more virulent than low prevalent sporadic clinical isolates, and highly virulent strains induce non-protective immune responses with some correlation with clinical-epidemiological data. In the future selected strains from these types of studies should be analyzed with molecular technologies. These kinds of study will contribute to the identification of mycobacterial genes associated with virulence and immunogenicity.

Importance of confirming data on the in vivo efficacy of novel antibacterial drug regimens against various strains of Mycobacterium tuberculosis

In preclinical testing of antituberculosis drugs, laboratory-adapted strains of Mycobacterium tuberculosis are usually used both for in vitro and in vivo studies. However, it is unknown whether the heterogeneity of M. tuberculosis stocks used by various laboratories can result in different outcomes in tests of antituberculosis drug regimens in animal infection models. In head-to-head studies, we investigated whether bactericidal efficacy results in BALB/c mice infected by inhalation with the laboratory-adapted strains H37Rv and Erdman differ from each other and from those obtained with clinical tuberculosis strains. Treatment of mice consisted of dual and triple drug combinations of isoniazid (H), rifampin (R), and pyrazinamide (Z). The results showed that not all strains gave the same in vivo efficacy results for the drug combinations tested. Moreover, the ranking of HRZ and RZ efficacy results was not the same for the two H37Rv strains evaluated. The magnitude of this strain difference also varied between experiments, emphasizing the risk of drawing firm conclusions for human trials based on single animal studies. The results also confirmed that the antagonism seen within the standard HRZ regimen by some investigators appears to be an M. tuberculosis strain-specific phenomenon. In conclusion, the specific identity of M. tuberculosis strain used was found to be an important variable that can change the apparent outcome of in vivo efficacy studies in mice. We highly recommend confirmation of efficacy results in late preclinical testing against a different M. tuberculosis strain than the one used in the initial mouse efficacy study, thereby increasing confidence to advance potent drug regimens to clinical trials.

Strain diversity, epistasis and the evolution of drug resistance in Mycobacterium tuberculosis

Mycobacterium tuberculosis harbours little DNA sequence diversity compared with other bacteria. However, there is mounting evidence that strain-to-strain variation in this organism has been underestimated. We review our current understanding of the genetic diversity among M. tuberculosis clinical strains and discuss the relevance of this diversity for the ongoing global epidemics of drug-resistant tuberculosis. Based on findings in other bacteria, we propose that epistatic interactions between pre-existing differences in strain genetic background, acquired drug-resistance-conferring mutations and compensatory changes could play a role in the emergence and spread of drug-resistant M. tuberculosis.

Genomic stability over 9 years of an isoniazid resistant Mycobacterium tuberculosis outbreak strain in Sweden

In molecular epidemiological studies of drug resistant Mycobacterium tuberculosis (TB) in Sweden a large outbreak of an isoniazid resistant strain was identified, involving 115 patients, mainly from the Horn of Africa. During the outbreak period, the genomic pattern of the outbreak strain has stayed virtually unchanged with regard to drug resistance, IS6110 restriction fragment length polymorphism and spoligotyping patterns. Here we present the complete genome sequence analyses of the index isolate and two isolates sampled nine years after the index case as well as experimental data on the virulence of this outbreak strain. Even though the strain has been present in the community for nine years and passaged between patients at least five times in-between the isolates, we only found four single nucleotide polymorphisms in one of the later isolates and a small (4 amino acids) deletion in the other compared to the index isolate. In contrast to many other evolutionarily successful outbreak lineages (e.g. the Beijing lineage) this outbreak strain appears to be genetically very stable yet evolutionarily successful in a low endemic country such as Sweden. These findings further illustrate that the rate of genomic variation in TB can be highly strain dependent, something that can have important implications for epidemiological studies as well as development of resistance.

Anchor-based whole genome phylogeny (ABWGP): a tool for inferring evolutionary relationship among closely related microorganisms [corrected]

Phenotypic behavior of a group of organisms can be studied using a range of molecular evolutionary tools that help to determine evolutionary relationships. Traditionally a gene or a set of gene sequences was used for generating phylogenetic trees. Incomplete evolutionary information in few selected genes causes problems in phylogenetic tree construction. Whole genomes are used as remedy. Now, the task is to identify the suitable parameters to extract the hidden information from whole genome sequences that truly represent evolutionary information. In this study we explored a random anchor (a stretch of 100 nucleotides) based approach (ABWGP) for finding distance between any two genomes, and used the distance estimates to compute evolutionary trees. A number of strains and species of Mycobacteria were used for this study. Anchor-derived parameters, such as cumulative normalized score, anchor order and indels were computed in a pair-wise manner, and the scores were used to compute distance/phylogenetic trees. The strength of branching was determined by bootstrap analysis. The terminal branches are clearly discernable using the distance estimates described here. In general, different measures gave similar trees except the trees based on indels. Overall the tree topology reflected the known biology of the organisms. This was also true for different strains of Escherichia coli. A new whole genome-based approach has been described here for studying evolutionary relationships among bacterial strains and species.

A survey of nucleotide cyclases in actinobacteria: unique domain organization and expansion of the class III cyclase family in Mycobacterium tuberculosis

Cyclic nucleotides are well-known second messengers involved in the regulation of important metabolic pathways or virulence factors. There are six different classes of nucleotide cyclases that can accomplish the task of generating cAMP, and four of these are restricted to the prokaryotes. The role of cAMP has been implicated in the virulence and regulation of secondary metabolites in the phylum Actinobacteria, which contains important pathogens, such as Mycobacterium tuberculosis, M. leprae, M. bovis and Corynebacterium, and industrial organisms from the genus Streptomyces. We have analysed the actinobacterial genome sequences found in current databases for the presence of different classes of nucleotide cyclases, and find that only class III cyclases are present in these organisms. Importantly, prominent members such as M. tuberculosis and M. leprae have 17 and 4 class III cyclases, respectively, encoded in their genomes, some of which display interesting domain fusions seen for the first time. In addition, a pseudogene corresponding to a cyclase from M. avium has been identified as the only cyclase pseudogene in M. tuberculosis and M. bovis. The Corynebacterium and Streptomyces genomes encode only a single adenylyl cyclase each, both of which have corresponding orthologues in M. tuberculosis. A clustering of the cyclase domains in Actinobacteria reveals the presence of typical eukaryote-like, fungi-like and other bacteria-like class III cyclase sequences within this phylum, suggesting that these proteins may have significant roles to play in this important group of organisms.

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.

A critical assessment of cross-species detection of gene duplicates using comparative genomic hybridization

Background

Comparison of genomic DNA among closely related strains or species is a powerful approach for identifying variation in evolutionary processes. One potent source of genomic variation is gene duplication, which is prevalent among individuals and species. Array comparative genomic hybridization (aCGH) has been successfully utilized to detect this variation among lineages. Here, beyond the demonstration that gene duplicates among species can be quantified with aCGH, we consider the effect of sequence divergence on the ability to detect gene duplicates.

Results

Using the X chromosome genomic content difference between male D. melanogaster and female D. yakuba and D. simulans, we describe a decrease in the ability to accurately measure genomic content (copy number) for orthologs that are only 90% identical. We demonstrate that genome characteristics (e.g. chromatin environment and non-orthologous sequence similarity) can also affect the ability to accurately measure genomic content. We describe a normalization strategy and statistical criteria to be used for the identification of gene duplicates among any species group for which an array platform is available from a closely related species.

Conclusions

Array CGH can be used to effectively identify gene duplication and genome content; however, certain biases are present due to sequence divergence and other genome characteristics resulting from the divergence between lineages. Highly conserved gene duplicates will be more readily recovered by aCGH. Duplicates that have been retained for a selective advantage due to directional selection acting on many loci in one or both gene copies are likely to be under-represented. The results of this study should inform the interpretation of both previously published and future work that employs this powerful technique.

Using comparative genomic hybridization to survey genomic sequence divergence across species: a proof-of-concept from Drosophila

Background

Genome-wide analysis of sequence divergence among species offers profound insights into the evolutionary processes that shape lineages. When full-genome sequencing is not feasible for a broad comparative study, we propose the use of array-based comparative genomic hybridization (aCGH) in order to identify orthologous genes with high sequence divergence. Here we discuss experimental design, statistical power, success rate, sources of variation and potential confounding factors. We used a spotted PCR product microarray platform from Drosophila melanogaster to assess sequence divergence on a gene-by-gene basis in three fully sequenced heterologous species (D. sechellia, D. simulans, and D. yakuba). Because complete genome assemblies are available for these species this study presents a powerful test for the use of aCGH as a tool to measure sequence divergence.

Results

We found a consistent and linear relationship between hybridization ratio and sequence divergence of the sample to the platform species. At higher levels of sequence divergence (< 92% sequence identity to D. melanogaster) ~84% of features had significantly less hybridization to the array in the heterologous species than the platform species, and thus could be identified as "diverged". At lower levels of divergence (≥ 97% identity), only 13% of genes were identified as diverged. While ~40% of the variation in hybridization ratio can be accounted for by variation in sequence identity of the heterologous sample relative to D. melanogaster, other individual characteristics of the DNA sequences, such as GC content, also contribute to variation in hybridization ratio, as does technical variation.

Conclusions

Here we demonstrate that aCGH can accurately be used as a proxy to estimate genome-wide divergence, thus providing an efficient way to evaluate how evolutionary processes and genomic architecture can shape species diversity in non-model systems. Given the increased number of species for which microarray platforms are available, comparative studies can be conducted for many interesting lineages in order to identify highly diverged genes that may be the target of natural selection.

Does M. tuberculosis genomic diversity explain disease diversity?

The outcome of tuberculosis infection and disease is highly variable. This variation has been attributed primarily to host and environmental factors, but better understanding of the global genomic diversity in the M. tuberculosis complex (MTBC) suggests that bacterial factors could also be involved. Review of nearly 100 published reports shows that MTBC strains differ in their virulence and immunogenicity in experimental models, but whether this phenotypic variation plays a role in human disease remains unclear. Given the complex interactions between the host, the pathogen and the environment, linking MTBC genotypic diversity to experimental and clinical phenotypes requires an integrated systems epidemiology approach embedded in a robust evolutionary framework.

Effect of serial subculturing on the genetic composition and cytotoxic activity of Mycobacterium tuberculosis

Continuous subculture has been observed to produce changes in the virulence of micro-organisms, e.g. rabies virus, poliovirus and Mycobacterium bovis BCG. The latter has been used as a vaccine for tuberculosis for the last 100 years; however, in some instances its efficacy has been observed to be very low. In order to determine whether similar changes can be produced in Mycobacterium tuberculosis, we selected four isolates, M. tuberculosis H37Rv, a Beijing strain (DR-689), and two more isolates with deletion of the phospholipase C locus (plcA-plcB-plcC ), and subjected them to serial culturing on Middlebrook 7H9 medium, with or without ox bile. After 100 passages, we performed RFLP-IS6110 analysis to determine whether genomic changes were produced. We also checked their genomic composition by microarray analysis. Changes in virulence were studied by measuring the cytotoxic effect of parental and subcultured isolates on a THP-1 macrophage monolayer. The most visible change was the change of position of an IS6110 band of approximately 1400 bp to approximately 1600 bp in the Beijing isolate subcultured in the ox bile medium. Analysis by microarray and PCR confirmation did not reveal any genomic changes. Cytotoxic activity was decreased in the isolates at levels close to that of BCG, and more consistently in those subcultured in the presence of ox bile.

Polymorphic exact tandem repeat A (PETRA): a newly defined lineage of mycobacterium tuberculosis in israel originating predominantly in Sub-Saharan Africa

As part of the Israel National Program for Prevention and Control of Tuberculosis, the molecular epidemiology of new tuberculosis cases is monitored. Prospective screening showed that about 20% of all new cases of culture-positive tuberculosis (43 of 222) in Israel in the year 2008 were caused by certain Mycobacterium tuberculosis strains of the central Asian (CAS) spoligotype lineage. The identity and similarity of these strains by mycobacterial interspersed repetitive-unit-variable-number tandem-repeat (MIRU-VNTR) typing form a lineage we call PETRA for polymorphic at locus ETR A. The name PETRA was given to 79 strains we have found since the year 2000, because the largest number of strains with MIRU-VNTR profiles identical other than at locus A formed three groups, including 5 of 10 strains that had deleted the ETR A region from their genomes. No PETRA strain was found to be multiple drug resistant (resistant to both isoniazid and rifampin [rifampicin]). Most patients (75% [58 of 77 patients of known origin]) infected with PETRA were of sub-Saharan African origins. The genotypes associated with the 79 PETRA lineage strains presented in this paper suggest that the PETRA lineage is a large, major contributor to new tuberculosis cases in Israel.

The past and future of tuberculosis research

Renewed efforts in tuberculosis (TB) research have led to important new insights into the biology and epidemiology of this devastating disease. Yet, in the face of the modern epidemics of HIV/AIDS, diabetes, and multidrug resistance--all of which contribute to susceptibility to TB--global control of the disease will remain a formidable challenge for years to come. New high-throughput genomics technologies are already contributing to studies of TB's epidemiology, comparative genomics, evolution, and host-pathogen interaction. We argue here, however, that new multidisciplinary approaches--especially the integration of epidemiology with systems biology in what we call "systems epidemiology"--will be required to eliminate TB.

Lack of insertional-deletional polymorphism in a collection of Mycobacterium ulcerans isolates from Ghanaian Buruli ulcer patients

Mycobacterium ulcerans causes the devastating infectious skin disease Buruli ulcer and has a monomorphic population structure. The resolution of conventional genetic fingerprinting methods is therefore not sufficient for microepidemiological studies aiming to characterize transmission pathways. In a previous comparative genomic hybridization analysis with a microarray covering part of the M. ulcerans genome, we have found extensive insertional-deletional sequence polymorphisms among M. ulcerans isolates of diverse geographic origins that allowed us to distinguish between strains coming from different continents. Since large numbers of insertion sequences are spread over the genome of African M. ulcerans strains, we reasoned that these may drive large sequence polymorphisms in otherwise clonal local mycobacterial populations. In this study, we used a printed DNA microarray covering the whole genome of the Ghanaian M. ulcerans reference strain Agy99 for comparative genomic hybridization. The assay identified multiple regions of difference when DNA of a Japanese M. ulcerans strain was analyzed. In contrast, not a single insertional-deletional genomic variation was found within a panel of disease isolates coming from an area of Ghana where Buruli ulcer is endemic. These results indicate that, despite the expectations deduced from other mycobacterial pathogens, only analyses of single nucleotide polymorphisms will have the potential to differentiate local populations of M. ulcerans.

MyBASE: a database for genome polymorphism and gene function studies of Mycobacterium

Background

Mycobacterial pathogens are a major threat to humans. With the increasing availability of functional genomic data, research on mycobacterial pathogenesis and subsequent control strategies will be greatly accelerated. It has been suggested that genome polymorphisms, namely large sequence polymorphisms, can influence the pathogenicity of different mycobacterial strains. However, there is currently no database dedicated to mycobacterial genome polymorphisms with functional interpretations.

Description

We have developed a mycobacterial database (MyBASE) housing genome polymorphism data and gene functions to provide the mycobacterial research community with a useful information resource and analysis platform. Whole genome comparison data produced by our lab and the novel genome polymorphisms identified were deposited into MyBASE. Extensive literature review of genome polymorphism data, mainly large sequence polymorphisms (LSPs), operon predictions and curated annotations of virulence and essentiality of mycobacterial genes are unique features of MyBASE. Large-scale genomic data integration from public resources makes MyBASE a comprehensive data warehouse useful for current research. All data is cross-linked and can be graphically viewed via a toolbox in MyBASE.

Conclusion

As an integrated platform focused on the collection of experimental data from our own lab and published literature, MyBASE will facilitate analysis of genome structure and polymorphisms, which will provide insight into genome evolution. Importantly, the database will also facilitate the comparison of virulence factors among various mycobacterial strains. MyBASE is freely accessible via http://mybase.psych.ac.cn.

Defining mycobacteria: Shared and specific genome features for different lifestyles

During the last decade, the combination of rapid whole genome sequencing capabilities, application of genetic and computational tools, and establishment of model systems for the study of a range of species for a spectrum of biological questions has enhanced our cumulative knowledge of mycobacteria in terms of their growth properties and requirements. The adaption of the corynebacterial surrogate system has simplified the study of cell wall biosynthetic machinery common to actinobacteria. Comparative genomics supported by experimentation reveals that superimposed on a common core of 'mycobacterial' gene set, pathogenic mycobacteria are endowed with multiple copies of several protein families that encode novel secretion and transport systems such as mce and esx; immunomodulators named PE/PPE proteins, and polyketide synthases for synthesis of complex lipids. The precise timing of expression, engagement and interactions involving one or more of these redundant proteins in their host environments likely play a role in the definition and differentiation of species and their disease phenotypes. Besides these, only a few species specific 'virulence' factors i.e., macromolecules have been discovered. Other subtleties may also arise from modifications of shared macromolecules. In contrast, to cope with the broad and changing growth conditions, their saprophytic relatives have larger genomes, in which the excess coding capacity is dedicated to transcriptional regulators, transporters for nutrients and toxic metabolites, biosynthesis of secondary metabolites and catabolic pathways. In this review, we present a sampling of the tools and techniques that are being implemented to tease apart aspects of physiology, phylogeny, ecology and pathology and illustrate the dominant genomic characteristics of representative species. The investigation of clinical isolates, natural disease states and discovery of new diagnostics, vaccines and drugs for existing and emerging mycobacterial diseases, particularly for multidrug resistant strains are the challenges in the coming decades.

Molecular Epidemiology

Molecular epidemiology is now an established discipline in epidemiology.(1) It is the contemporary stage in the evolution of laboratory-based epidemiology that may have begun with the discovery in the late 1800s of ways to differentiate bacterial organisms by pure culture in artificial media.(2) Molecular epidemiology uses new molecular biology tools to address questions difficult or not possible to address by old laboratory tools. Just as statistical tools have become indispensable in epidemiological investigations and interpretations of epidemiologic data, molecular biology tools today have come to elucidate epidemiologic features of diseases that cannot be easily characterized by conventional techniques. Applied to infectious diseases, molecular biology methods have also come to challenge our traditional notions about the epidemiology of these diseases and have engendered novel opportunities for their prevention and control. This chapter will (1) review definitions commonly used in molecular epidemiology, (2) present an overview of molecular biology methods used to study infectious disease epidemiology, and (3) describe examples of the types of epidemiologic problems that can be addressed by molecular biology techniques, highlighting new concepts that emerged in the process of applying this approach to study bacterial infectious diseases.

Comparative genomics of vesicomyid clam (Bivalvia: Mollusca) chemosynthetic symbionts

BACKGROUND

The Vesicomyidae (Bivalvia: Mollusca) are a family of clams that form symbioses with chemosynthetic gamma-proteobacteria. They exist in environments such as hydrothermal vents and cold seeps and have a reduced gut and feeding groove, indicating a large dependence on their endosymbionts for nutrition. Recently, two vesicomyid symbiont genomes were sequenced, illuminating the possible nutritional contributions of the symbiont to the host and making genome-wide evolutionary analyses possible.

RESULTS

To examine the genomic evolution of the vesicomyid symbionts, a comparative genomics framework, including the existing genomic data combined with heterologous microarray hybridization results, was used to analyze conserved gene content in four vesicomyid symbiont genomes. These four symbionts were chosen to include a broad phylogenetic sampling of the vesicomyid symbionts and represent distinct chemosynthetic environments: cold seeps and hydrothermal vents.

CONCLUSION

The results of this comparative genomics analysis emphasize the importance of the symbionts' chemoautotrophic metabolism within their hosts. The fact that these symbionts appear to be metabolically capable autotrophs underscores the extent to which the host depends on them for nutrition and reveals the key to invertebrate colonization of these challenging environments.

Does the lipR gene of tubercle bacilli have a role in tuberculosis transmission and pathogenesis?

Mycobacterium tuberculosis lipases, a diverse class of enzymes involved in lipid metabolism, may have an important role in tuberculosis (TB) pathogenesis. We explored the association of large sequence polymorphism (LSP) in one of the M. tuberculosis lipase-encoding genes, lipR (Rv3084), with patient characteristics using a population-based sample of clinical isolates to elucidate the potential role of lipR in TB pathogenesis. LSP in lipR was found in 104 (15.6%) of 665 isolates, of which 96% belonged to principal genetic group 3. When linkage by molecular type and epidemiologic evidence were compared, molecularly clustered cases infected with a lipR LSP isolate were more often epidemiologically linked than clustered cases infected with a lipR wild-type isolate. Further epidemiologic and functional studies are necessary to determine if the association between this lipR LSP and recent transmission we identified in this population reflects a functional role of lipR in TB transmission and pathogenesis or other unidentified mechanisms.