Comparative genomics between human and animal associated subspecies of the Mycobacterium avium complex: a basis for pathogenicity

Background

A human isolate of Mycobacterium avium subsp. paratuberculosis (M. paratuberculosis 43525) was sequenced and compared genomically to other mycobacterial pathogens. M. paratuberculosis 43525 was recently isolated from a patient with ulcerative colitis and belongs to the M. avium complex, a group known to infect both humans and animals. While M. paratuberculosis is a known pathogen of livestock, there are only 20 human isolates from the last 20 years, therefore we took the opportunity to perform a whole genome comparison between human and animal mycobacterial pathogens. We also compared virulence determinants such as the mycobactin cluster, PE/PPE genes and mammalian cell entry (mce) operons between MAC subspecies that infect animals and those that infect humans. M. tuberculosis was also included in these analyses given its predominant role as a human pathogen.

Results

This genome comparison showed the PE/PPE profile of M. paratuberculosis 43525 to be largely the same as other M. paratuberculosis isolates, except that it had one PPE and one PE_PGRS protein that are only present in human MAC strains and M. tuberculosis. PE/PPE proteins that were unique to M. paratuberculosis 43525, M. avium subsp. hominissuis and a caprine M. paratuberculosis isolate, were also identified. In addition, the mycobactin cluster differed between human and animal isolates and a unique mce operon flanked by two mycobactin genes, mbtA and mbtJ, was identified in all available M. paratuberculosis genomes.

Conclusions

Despite the whole genome comparison placing M. paratuberculosis 43525 as closely related to bovine M. paratuberculosis, key virulence factors were similar to human mycobacterial pathogens. This study highlights key factors of mycobacterial pathogenesis in humans and forms the basis for future functional studies.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1889-2) contains supplementary material, which is available to authorized users.     

Use of Whole Genome Sequencing to Determine the Microevolution of Mycobacterium tuberculosis during an Outbreak

Rationale

Current tools available to study the molecular epidemiology of tuberculosis do not provide information about the directionality and sequence of transmission for tuberculosis cases occurring over a short period of time, such as during an outbreak. Recently, whole genome sequencing has been used to study molecular epidemiology of Mycobacterium tuberculosis over short time periods.

Objective

To describe the microevolution of M. tuberculosis during an outbreak caused by one drug-susceptible strain.

Method and Measurements

We included 9 patients with tuberculosis diagnosed during a period of 22 months, from a population-based study of the molecular epidemiology in San Francisco. Whole genome sequencing was performed using Illumina’s sequencing by synthesis technology. A custom program written in Python was used to determine single nucleotide polymorphisms which were confirmed by PCR product Sanger sequencing.

Main results

We obtained an average of 95.7% (94.1–96.9%) coverage for each isolate and an average fold read depth of 73 (1 to 250). We found 7 single nucleotide polymorphisms among the 9 isolates. The single nucleotide polymorphisms data confirmed all except one known epidemiological link. The outbreak strain resulted in 5 bacterial variants originating from the index case A1 with 0–2 mutations per transmission event that resulted in a secondary case.

Conclusions

Whole genome sequencing analysis from a recent outbreak of tuberculosis enabled us to identify microevolutionary events observable during transmission, to determine 0–2 single nucleotide polymorphisms per transmission event that resulted in a secondary case, and to identify new epidemiologic links in the chain of transmission.     

Distribution of spoligotyping defined genotypic lineages among drug-resistant Mycobacterium tuberculosis complex clinical isolates in Ankara, Turkey

Background

Investigation of genetic heterogeneity and spoligotype-defined lineages of drug-resistant Mycobacterium tuberculosis clinical isolates collected during a three-year period in two university hospitals and National Tuberculosis Reference and Research Laboratory in Ankara, Turkey.

Methods and Findings

A total of 95 drug-resistant M. tuberculosis isolates collected from three different centers were included in this study. Susceptibility testing of the isolates to four major antituberculous drugs was performed using proportion method on Löwenstein–Jensen medium and BACTEC 460-TB system. All clinical isolates were typed by using spoligotyping and IS6110-restriction fragment length polymorphism (RFLP) methods. Seventy-three of the 95 (76.8%) drug resistant M. tuberculosis isolates were isoniazid-resistant, 45 (47.4%) were rifampicin-resistant, 32 (33.7%) were streptomycin-resistant and 31 (32.6%) were ethambutol-resistant. The proportion of multidrug-resistant isolates (MDR) was 42.1%. By using spoligotyping, 35 distinct patterns were observed; 75 clinical isolates were grouped in 15 clusters (clustering rate of 79%) and 20 isolates displayed unique patterns. Five of these 20 unique patterns corresponded to orphan patterns in the SITVIT2 database, while 4 shared types containing 8 isolates were newly created. The most prevalent M. tuberculosis lineages were: Haarlem (23/95, 24.2%), ill-defined T superfamily (22/95, 23.2%), the Turkey family (19/95, 20%; previously designated as LAM7-TUR), Beijing (6/95, 6.3%), and Latin-America & Mediterranean (LAM, 5/95 or 5.3%), followed by Manu (3/95, 3.2%) and S (1/95, 1%) lineages. Four of the six Beijing family isolates (66.7%) were MDR. A combination of IS6110-RFLP and spoligotyping reduced the clustering rate from 79% to 11.5% among the drug resistant isolates.

Conclusions

The results obtained showed that ill-defined T, Haarlem, the Turkey family (previously designated as LAM7-TUR family with high phylogeographical specifity for Turkey), Beijing and LAM were predominant lineages observed in almost 80% of the drug-Resistant M. tuberculosis complex clinical isolates in Ankara, Turkey.     

Detecting Novel Genetic Variants Associated with Isoniazid-Resistant Mycobacterium tuberculosis

Background

Isoniazid (INH) is a highly effective antibiotic central for the treatment of Mycobacterium tuberculosis (MTB). INH-resistant MTB clinical isolates are frequently mutated in the katG gene and the inhA promoter region, but 10 to 37% of INH-resistant clinical isolates have no detectable alterations in currently known gene targets associated with INH-resistance. We aimed to identify novel genes associated with INH-resistance in these latter isolates.

Methodology/Principal Findings

INH-resistant clinical isolates of MTB were pre-screened for mutations in the katG, inhA, kasA and ndh genes and the regulatory regions of inhA and ahpC. Twelve INH-resistant isolates with no mutations, and 17 INH-susceptible MTB isolates were subjected to whole genome sequencing. Phylogenetically related variants and synonymous mutations were excluded and further analysis revealed mutations in 60 genes and 4 intergenic regions associated with INH-resistance. Sanger sequencing verification of 45 genes confirmed that mutations in 40 genes were observed only in INH-resistant isolates and not in INH-susceptible isolates. The ratios of non-synonymous to synonymous mutations (dN/dS ratio) for the INH-resistance associated mutations identified in this study were 1.234 for INH-resistant and 0.654 for INH-susceptible isolates, strongly suggesting that these mutations are indeed associated with INH-resistance.

Conclusion

The discovery of novel targets associated with INH-resistance described in this study may potentially be important for the development of improved molecular detection strategies.     

Lipid metabolism and Type VII secretion systems dominate the genome scale virulence profile of Mycobacterium tuberculosis in human dendritic cells

Background

Mycobacterium tuberculosis continues to kill more people than any other bacterium. Although its archetypal host cell is the macrophage, it also enters, and survives within, dendritic cells (DCs). By modulating the behaviour of the DC, M. tuberculosis is able to manipulate the host’s immune response and establish an infection. To identify the M. tuberculosis genes required for survival within DCs we infected primary human DCs with an M. tuberculosis transposon library and identified mutations with a reduced ability to survive.

Results

Parallel sequencing of the transposon inserts of the surviving mutants identified a large number of genes as being required for optimal intracellular fitness in DCs. Loci whose mutation attenuated intracellular survival included those involved in synthesising cell wall lipids, not only the well-established virulence factors, pDIM and cord factor, but also sulfolipids and PGL, which have not previously been identified as having a direct virulence role in cells. Other attenuated loci included the secretion systems ESX-1, ESX-2 and ESX-4, alongside many PPE genes, implicating a role for ESX-5. In contrast the canonical ESAT-6 family of ESX substrates did not have intra-DC fitness costs suggesting an alternative ESX-1 associated virulence mechanism. With the aid of a gene-nutrient interaction model, metabolic processes such as cholesterol side chain catabolism, nitrate reductase and cysteine-methionine metabolism were also identified as important for survival in DCs.

Conclusion

We conclude that many of the virulence factors required for survival in DC are shared with macrophages, but that survival in DCs also requires several additional functions, such as cysteine-methionine metabolism, PGLs, sulfolipids, ESX systems and PPE genes.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1569-2) contains supplementary material, which is available to authorized users.     

Construction of a virtual Mycobacterium tuberculosis consensus genome and its application to data from a next generation sequencer

Background

Although Mycobacterium tuberculosis isolates are consisted of several different lineages and the epidemiology analyses are usually assessed relative to a particular reference genome, M. tuberculosis H37Rv, which might introduce some biased results. Those analyses are essentially based genome sequence information of M. tuberculosis and could be performed in sillico in theory, with whole genome sequence (WGS) data available in the databases and obtained by next generation sequencers (NGSs). As an approach to establish higher resolution methods for such analyses, whole genome sequences of the M. tuberculosis complexes (MTBCs) strains available on databases were aligned to construct virtual reference genome sequences called the consensus sequence (CS), and evaluated its feasibility in in sillico epidemiological analyses.

Results

The consensus sequence (CS) was successfully constructed and utilized to perform phylogenetic analysis, evaluation of read mapping efficacy, which is crucial for detecting single nucleotide polymorphisms (SNPs), and various MTBC typing methods virtually including spoligotyping, VNTR, Long sequence polymorphism and Beijing typing. SNPs detected based on CS, in comparison with H37Rv, were utilized in concatemer-based phylogenetic analysis to determine their reliability relative to a phylogenetic tree based on whole genome alignment as the gold standard. Statistical comparison of phylogenic trees based on CS with that of H37Rv indicated the former showed always better results that that of later. SNP detection and concatenation with CS was advantageous because the frequency of crucial SNPs distinguishing among strain lineages was higher than those of H37Rv. The number of SNPs detected was lower with the consensus than with the H37Rv sequence, resulting in a significant reduction in computational time. Performance of each virtual typing was satisfactory and accorded with those published when those are available.

Conclusions

These results indicated that virtual CS constructed from genome sequence data is an ideal approach as a reference for MTBC studies.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1368-9) contains supplementary material, which is available to authorized users.     

Fitness of Mycobacterium tuberculosis Strains of the W-Beijing and Non-W-Beijing Genotype

Background

Multidrug resistant tuberculosis (MDR-TB) is a major threat for global tuberculosis control. The W-Beijing Mycobacterium tuberculosis genotype has been associated with drug resistance. Elucidation of the mechanisms underlying this epidemiological finding may have an important role in the control of MDR-TB. The aim of this study was to evaluate the fitness of drug-susceptible and MDR M. tuberculosis strains of the W-Beijing genotype compared with that of Non-W-Beijing strains.

Methodology/Principal Findings

Fitness of M. tuberculosis strains was determined by evaluating the difference in the growth curves obtained in the MGIT960 automated system and assessing the competitive growth capacity between W-Beijing and non-W-Beijing strains. The W-Beijing MDR strains had a significant longer lag phase duration compared to the other strains but did not present a significant fitness cost. When grown in competition they had an advantage only in medium containing 0.1% Tween 80.

Conclusions/Significance

It was not possible to confirm a selective advantage of W-Beijing strains to grow, except for differences in their resistance to Tween 80. Further studies are needed to elucidate the putative advantage of W-Beijing strains compared to other genotypes.     

Optimizing multiplex SNP-based data analysis for genotyping of Mycobacterium tuberculosis isolates

Background

Multiplex ligation-dependent probe amplification (MLPA) is a powerful tool to identify genomic polymorphisms. We have previously developed a single nucleotide polymorphism (SNP) and large sequence polymorphisms (LSP)-based MLPA assay using a read out on a liquid bead array to screen for 47 genetic markers in the Mycobacterium tuberculosis genome. In our assay we obtain information regarding the Mycobacterium tuberculosis lineage and drug resistance simultaneously. Previously we called the presence or absence of a genotypic marker based on a threshold signal level. Here we present a more elaborate data analysis method to standardize and streamline the interpretation of data generated by MLPA. The new data analysis method also identifies intermediate signals in addition to classification of signals as positive and negative. Intermediate calls can be informative with respect to identifying the simultaneous presence of sensitive and resistant alleles or infection with multiple different Mycobacterium tuberculosis strains.

Results

To validate our analysis method 100 DNA isolates of Mycobacterium tuberculosis extracted from cultured patient material collected at the National TB Reference Laboratory of the National Center for Tuberculosis and Lung Diseases in Tbilisi, Republic of Georgia were tested by MLPA. The data generated were interpreted blindly and then compared to results obtained by reference methods. MLPA profiles containing intermediate calls are flagged for expert review whereas the majority of profiles, not containing intermediate calls, were called automatically. No intermediate signals were identified in 74/100 isolates and in the remaining 26 isolates at least one genetic marker produced an intermediate signal.

Conclusion

Based on excellent agreement with the reference methods we conclude that the new data analysis method performed well. The streamlined data processing and standardized data interpretation allows the comparison of the Mycobacterium tuberculosis MLPA results between different experiments. All together this will facilitate the implementation of the MLPA assay in different settings.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-572) contains supplementary material, which is available to authorized users.     

Genetic diversity of Mycobacterium tuberculosis isolates from Tibetans in Tibet, China

Background

Tuberculosis (TB) is a serious health problem in Tibet where Tibetans are the major ethnic group. Although genotyping of Mycobacterium tuberculosis (M. tuberculosis) isolates is a valuable tool for TB control, our knowledge of population structure of M. tuberculosis circulating in Tibet is limited.

Methodology/Principal Findings

In our study, a total of 576 M. tuberculosis isolates from Tibetans in Tibet, China, were analyzed via spoligotyping and 24-locus MIRU-VNTR. The Beijing genotype was the most prevalent family (90.63%, n = 522). Shared-type (ST) 1 was the most dominant genotype (88.89%, n = 512). We found that there was no association between the Beijing genotype and sex, age and treatment status. In this sample collection, 7 of the 24 MIRU-VNTR loci were highly or moderately discriminative according to their Hunter-Gaston discriminatory index. An informative set of 12 loci had similar discriminatory power with 24 loci set.

Conclusions/Significance

The population structure of M. tuberculosis isolates in Tibetans is homogeneous and dominated by Beijing genotype. The analysis of 24-locus MIRU-VNTR data might be useful to select appropriate VNTR loci for the genotyping of M. tuberculosis.     

First Insights into the Phylogenetic Diversity of Mycobacterium tuberculosis in Nepal

Background

Tuberculosis (TB) is a major public health problem in Nepal. Strain variation in Mycobacterium tuberculosis may influence the outcome of TB infection and disease. To date, the phylogenetic diversity of M. tuberculosis in Nepal is unknown.

Methods and Findings

We analyzed 261 M. tuberculosis isolates recovered from pulmonary TB patients recruited between August 2009 and August 2010 in Nepal. M. tuberculosis lineages were determined by single nucleotide polymorphisms (SNP) typing and spoligotyping. Drug resistance was determined by sequencing the hot spot regions of the relevant target genes. Overall, 164 (62.8%) TB patients were new, and 97 (37.2%) were previously treated. Any drug resistance was detected in 50 (19.2%) isolates, and 16 (6.1%) were multidrug-resistant. The most frequent M. tuberculosis lineage was Lineage 3 (CAS/Delhi) with 106 isolates (40.6%), followed by Lineage 2 (East-Asian lineage, includes Beijing genotype) with 84 isolates (32.2%), Lineage 4 (Euro-American lineage) with 41 (15.7%) isolates, and Lineage 1 (Indo-Oceanic lineage) with 30 isolates (11.5%). Based on spoligotyping, we found 45 different spoligotyping patterns that were previously described. The Beijing (83 isolates, 31.8%) and CAS spoligotype (52, 19.9%) were the dominant spoligotypes. A total of 36 (13.8%) isolates could not be assigned to any known spoligotyping pattern. Lineage 2 was associated with female sex (adjusted odds ratio [aOR] 2.58, 95% confidence interval [95% CI] 1.42–4.67, p = 0.002), and any drug resistance (aOR 2.79; 95% CI 1.43–5.45; p = 0.002). We found no evidence for an association of Lineage 2 with age or BCG vaccination status.

Conclusions

We found a large genetic diversity of M. tuberculosis in Nepal with representation of all four major lineages. Lineages 3 and 2 were dominating. Lineage 2 was associated with clinical characteristics. This study fills an important gap on the map of the M. tuberculosis genetic diversity in the Asian region.     

Identification of T-Cell Antigens Specific for Latent Mycobacterium Tuberculosis Infection

Background

T-cell responses against dormancy-, resuscitation-, and reactivation-associated antigens of Mycobacterium tuberculosis are candidate biomarkers of latent infection in humans.

Methodology/Principal Findings

We established an assay based on two rounds of in vitro restimulation and intracellular cytokine analysis that detects T-cell responses to antigens expressed during latent M. tuberculosis infection. Comparison between active pulmonary tuberculosis (TB) patients and healthy latently M. tuberculosis-infected donors (LTBI) revealed significantly higher T-cell responses against 7 of 35 tested M. tuberculosis latency-associated antigens in LTBI. Notably, T cells specific for Rv3407 were exclusively detected in LTBI but not in TB patients. The T-cell IFNγ response against Rv3407 in individual donors was the most influential factor in discrimination analysis that classified TB patients and LTBI with 83% accuracy using cross-validation. Rv3407 peptide pool stimulations revealed distinct candidate epitopes in four LTBI.

Conclusions

Our findings further support the hypothesis that the latency-associated antigens can be exploited as biomarkers for LTBI.     

Genomic Diversity among Drug Sensitive and Multidrug Resistant Isolates of Mycobacterium tuberculosis with Identical DNA Fingerprints

Background

Mycobacterium tuberculosis complex (MTBC), the causative agent of tuberculosis (TB), is characterized by low sequence diversity making this bacterium one of the classical examples of a genetically monomorphic pathogen. Because of this limited DNA sequence variation, routine genotyping of clinical MTBC isolates for epidemiological purposes relies on highly discriminatory DNA fingerprinting methods based on mobile and repetitive genetic elements. According to the standard view, isolates exhibiting the same fingerprinting pattern are considered direct progeny of the same bacterial clone, and most likely reflect ongoing transmission or disease relapse within individual patients.

Methodology/Principal Findings

Here we further investigated this assumption and used massively parallel whole-genome sequencing to compare one drug-susceptible (K-1) and one multidrug resistant (MDR) isolate (K-2) of a rapidly spreading M. tuberculosis Beijing genotype clone from a high incidence region (Karakalpakstan, Uzbekistan). Both isolates shared the same IS6110 RFLP pattern and the same allele at 23 out of 24 MIRU-VNTR loci.We generated 23.9 million (K-1) and 33.0 million (K-2) paired 50 bp purity filtered reads corresponding to a mean coverage of 483.5 fold and 656.1 fold respectively. Compared with the laboratory strain H37Rv both Beijing isolates shared 1,209 SNPs. The two Beijing isolates differed by 130 SNPs and one large deletion. The susceptible isolate had 55 specific SNPs, while the MDR variant had 75 specific SNPs, including the five known resistance-conferring mutations.

Conclusions

Our results suggest that M. tuberculosis isolates exhibiting identical DNA fingerprinting patterns can harbour substantial genomic diversity. Because this heterogeneity is not captured by traditional genotyping of MTBC, some aspects of the transmission dynamics of tuberculosis could be missed or misinterpreted. Furthermore, a valid differentiation between disease relapse and exogenous reinfection might be impossible using standard genotyping tools if the overall diversity of circulating clones is limited. These findings have important implications for clinical trials of new anti-tuberculosis drugs.     

Microscopic-Observation Drug-Susceptibility Assay for the Diagnosis of Drug-Resistant Tuberculosis in Harare,Zimbabwe

Introduction

Limited data exist on use of the microscopic-observation drug-susceptibility (MODS) assay among persons suspected of MDR-TB living in high HIV-prevalence settings.

Methods

We retrospectively reviewed available clinical and drug susceptibility data for drug-resistant TB suspects referred for culture and drug-susceptibility testing between April 1, 2011 and March 1, 2012. The diagnostic accuracy of MODS was estimated against a reference standard including Löwenstein-Jensen (LJ) media and manual liquid (BACTEC MGIT) culture. The accuracy of MODS drug-susceptibility testing (DST) was assessed against a reference standard absolute concentration method.

Results

One hundred thirty-eight sputum samples were collected from 99 drug-resistant TB suspects; in addition, six previously cultured MDR isolates were included for assessment of DST accuracy. Among persons with known HIV infection status, 39/59 (66%) were HIV-infected. Eighty-six percent of patients had a history of prior TB treatment, and 80% of individuals were on antituberculous treatment at the time of sample collection. M. tuberculosis was identified by reference standard culture among 34/98 (35%) MDR-TB suspects. Overall MODS sensitivity for M. tuberculosis detection was 85% (95% CI, 69–95%) and specificity was 93% (95% CI, 84–98%); diagnostic accuracy did not significantly differ by HIV infection status. Median time to positivity was significantly shorter for MODS (7 days; IQR 7–15 days) than MGIT (12 days; IQR 6–16 days) or LJ (28 days; IQR 21–35 days; p<0.001). Of 33 specimens with concurrent DST results, sensitivity of the MODS assay for detection of resistance to isoniazid, rifampin, and MDR-TB was 88% (95% CI, 68–97%), 96% (95% CI, 79–100%), and 91% (95% CI, 72–99%), respectively; specificity was 89% (95% CI, 52–100%), 89% (95% CI, 52–100%), and 90% (95% CI, 56–100%), respectively.

Conclusion

In a high HIV-prevalence setting, MODS diagnosed TB and drug-resistant TB with high sensitivity and shorter turnaround time compared with standard culture and DST methods.     

Disruption of Mycobacterium avium subsp. paratuberculosis-specific genes impairs in vivo fitness

Background

Mycobacterium avium subsp. paratuberculosis (MAP) is an obligate intracellular pathogen that infects many ruminant species. The acquisition of foreign genes via horizontal gene transfer has been postulated to contribute to its pathogenesis, as these genetic elements are absent from its putative ancestor, M. avium subsp. hominissuis (MAH), an environmental organism with lesser pathogenicity. In this study, high-throughput sequencing of MAP transposon libraries were analyzed to qualitatively and quantitatively determine the contribution of individual genes to bacterial survival during infection.

Results

Out of 52384 TA dinucleotides present in the MAP K-10 genome, 12607 had a MycoMarT7 transposon in the input pool, interrupting 2443 of the 4350 genes in the MAP genome (56%). Of 96 genes situated in MAP-specific genomic islands, 82 were disrupted in the input pool, indicating that MAP-specific genomic regions are dispensable for in vitro growth (odds ratio = 0.21). Following 5 independent in vivo infections with this pool of mutants, the correlation between output pools was high for 4 of 5 (R = 0.49 to 0.61) enabling us to define genes whose disruption reproducibly reduced bacterial fitness in vivo. At three different thresholds for reduced fitness in vivo, MAP-specific genes were over-represented in the list of predicted essential genes. We also identified additional genes that were severely depleted after infection, and several of them have orthologues that are essential genes in M. tuberculosis.

Conclusions

This work indicates that the genetic elements required for the in vivo survival of MAP represent a combination of conserved mycobacterial virulence genes and MAP-specific genes acquired via horizontal gene transfer. In addition, the in vitro and in vivo essential genes identified in this study may be further characterized to offer a better understanding of MAP pathogenesis, and potentially contribute to the discovery of novel therapeutic and vaccine targets.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-415) contains supplementary material, which is available to authorized users.     

Molecular Epidemiology and Clinical Characteristics of Drug-Resistant Mycobacterium tuberculosis in a Tuberculosis Referral Hospital in China

Background

Despite the large number of drug-resistant tuberculosis (TB) cases in China, few studies have comprehensively analyzed the drug resistance-associated gene mutations and genotypes in relation to the clinical characteristics of M. tuberculosis (Mtb) isolates.

Methodology/Principal Findings

We thus analyzed the phenotypic and genotypic drug resistance profiles of 115 Mtb clinical isolates recovered from a tuberculosis referral hospital in Beijing, China. We also performed genotyping by 28 loci MIRU-VNTR analysis. Socio-demographic and clinical data were retrieved from medical records and analyzed. In total, 78 types of mutations (including 42 previously reported and 36 newly identified ones) were identified in 115 Mtb clinical isolates. There was significant correlation between phenotypic and genotypic drug resistance rates for first-line anti-TB drugs (P<0.001). Genotyping revealed 101 MIRU-VNTR types, with 20 isolates (17.4%) being clustered and 95 isolates (82.6%) having unique genotypes. Higher proportion of re-treatment cases was observed among patients with clustered isolates than those with unique MIRU-VNTR genotypes (75.0% vs. 41.1%). Moreover, clinical epidemiological links were identified among patients infected by Mtb strains belonging to the same clusters, suggesting a potential of transmission among patients.

Conclusions/Significance

Our study provided information on novel potential drug resistance-associated mutations in Mtb. In addition, the genotyping data from our study suggested that enforcement of the implementation of genotyping in diagnostic routines would provide important information for better monitor and control of TB transmission.     

Combination of Cytokine Responses Indicative of Latent TB and Active TB in Malawian Adults

Background

An IFN-γ response to M. tuberculosis-specific antigens is an effective biomarker for M. tuberculosis infection but it cannot discriminate between latent TB infection and active TB disease. Combining a number of cytokine/chemokine responses to M. tuberculosis antigens may enable differentiation of latent TB from active disease.

Methods

Asymptomatic recently-exposed individuals (spouses of TB patients) were recruited and tuberculin skin tested, bled and followed-up for two years. Culture supernatants, from a six-day culture of diluted whole blood samples stimulated with M. tuberculosis-derived PPD or ESAT-6, were measured for IFN-γ, IL-10, IL-13, IL-17, TNF-α and CXCL10 using cytokine ELISAs. In addition, 15 patients with sputum smear-positive pulmonary TB were recruited and tested.

Results

Spouses with positive IFN-γ responses to M. tuberculosis ESAT-6 (>62.5 pg/mL) and TB patients showed high production of IL-17, CXCL10 and TNF-α. Higher production of IL-10 and IL-17 in response to ESAT-6 was observed in the spouses compared with TB patients while the ratios of IFN-γ/IL-10 and IFN-γ/IL-17 in response to M. tuberculosis-derived PPD were significantly higher in TB patients compared with the spouses. Tuberculin skin test results did not correlate with cytokine responses.

Conclusions

CXCL10 and TNF-α may be used as adjunct markers alongside an IFN-γ release assay to diagnose M. tuberculosis infection, and IL-17 and IL-10 production may differentiate individuals with LTBI from active TB.     

A single-step sequencing method for the identification of Mycobacterium tuberculosis complex species

Background

The Mycobacterium tuberculosis complex (MTC) comprises closely related species responsible for strictly human and zoonotic tuberculosis. Accurate species determination is useful for the identification of outbreaks and epidemiological links. Mycobacterium africanum and Mycobacterium canettii are typically restricted to Africa and M. bovis is a re-emerging pathogen. Identification of these species is difficult and expensive.

Methodology/Principal Findings

The Exact Tandem Repeat D (ETR-D; alias Mycobacterial Interspersed Repetitive Unit 4) was sequenced in MTC species type strains and 110 clinical isolates, in parallel to reference polyphasic identification based on phenotype profiling and sequencing of pncA, oxyR, hsp65, gyrB genes and the major polymorphism tandem repeat. Inclusion of M. tuberculosis isolates in the expanding, antibiotic-resistant Beijing clone was determined by Rv0927c gene sequencing. The ETR-D (780-bp) sequence unambiguously identified MTC species type strain except M. pinnipedii and M. microti thanks to six single nucleotide polymorphisms, variable numbers (1–7 copies) of the tandem repeat and two deletions/insertions. The ETR-D sequencing agreed with phenotypic identification in 107/110 clinical isolates and with reference polyphasic molecular identification in all isolates, comprising 98 M. tuberculosis, 5 M. bovis BCG type, 5 M. canettii, and 2 M. africanum. For M. tuberculosis isolates, the ETR-D sequence was not significantly associated with the Beijing clone.

Conclusions/Significance

ETR-D sequencing allowed accurate, single-step identification of the MTC at the species level. It circumvented the current expensive, time-consuming polyphasic approach. It could be used to depict epidemiology of zoonotic and human tuberculosis, especially in African countries where several MTC species are emerging.     

Transmission of Mycobacterium tuberculosis between Farmers and Cattle in Central Ethiopia

Background

Transmission of Mycobacterium tuberculosis (M. tuberculosis) complex could be possible between farmers and their cattle in Ethiopia.

Methodology/Principal Findings

A study was conducted in mixed type multi-purposes cattle raising region of Ethiopia on 287 households (146 households with case of pulmonary tuberculosis (TB) and 141 free of TB) and 287 herds consisting of 2,033 cattle belonging to these households to evaluate transmission of TB between cattle and farmers. Interview, bacteriological examinations and molecular typing were used for human subjects while comparative intradermal tuberculin (CIDT) test, post mortem and bacteriological examinations, and molecular typing were used for animal studies. Herd prevalence of CIDT reactors was 9.4% and was higher (p<0.01) in herds owned by households with TB than in herds owned by TB free households. Animal prevalence was 1.8% and also higher (p<0.01) in cattle owned by households with TB case than in those owned by TB free households. All mycobacteria (141) isolated from farmers were M. tuberculosis, while only five of the 16 isolates from cattle were members of the M. tuberculosis complex (MTC) while the remaining 11 were members of non-tuberculosis mycobacteria (NTM). Further speciation of the five MTC isolates showed that three of the isolates were M. bovis (strain SB1176), while the remaining two were M. tuberculosis strains (SIT149 and SIT53). Pathology scoring method described by “Vordermeier et al. (2002)” was applied and the average severity of pathology in two cattle infected with M. bovis, in 11 infected with NTM and two infected with M. tuberculosis were 5.5, 2.1 and 0.5, respectively.

Conclusions/Significance

The results showed that transmission of TB from farmers to cattle by the airborne route sensitizes the cows but rarely leads to TB. Similarly, low transmission of M. bovis between farmers and their cattle was found, suggesting requirement of ingestion of contaminated milk from cows with tuberculous mastitis.