Single nucleotide polymorphisms in cell wall biosynthesis-associated genes and phylogeny of Mycobacterium tuberculosis lineages

Using genomics to understand the origin and dispersion of multidrug and extensively drug resistant tuberculosis in Portugal

Portugal is a low incidence country for tuberculosis (TB) disease. Now figuring among TB low incidence countries, it has since the 1990s reported multidrug resistant and extensively drug resistant (XDR) TB cases, driven predominantly by two strain-types: Lisboa3 and Q1. This study describes the largest characterization of the evolutionary trajectory of M/XDR-TB strains in Portugal, spanning a time-period of two decades. By combining whole-genome sequencing and phenotypic susceptibility data for 207 isolates, we report the geospatial patterns of drug resistant TB, particularly the dispersion of Lisboa3 and Q1 clades, which underly 64.2% and 94.0% of all MDR-TB and XDR-TB isolates, respectively. Genomic-based similarity and a phylogenetic analysis revealed multiple clusters (n = 16) reflecting ongoing and uncontrolled recent transmission of M/XDR-TB, predominantly associated with the Lisboa3 and Q1 clades. These clades are now thought to be evolving in a polycentric mode across multiple geographical districts. The inferred evolutionary history is compatible with MDR- and XDR-TB originating in Portugal in the 70's and 80's, respectively, but with subsequent multiple emergence events of MDR and XDR-TB particularly involving the Lisboa3 clade. A SNP barcode was defined for Lisboa3 and Q1 and comparison with a phylogeny of global strain-types (n = 28 385) revealed the presence of Lisboa3 and Q1 strains in Europe, South America and Africa. In summary, Portugal displays an unusual and unique epidemiological setting shaped by >40 years of uncontrolled circulation of two main phylogenetic clades, leading to a sympatric evolutionary trajectory towards XDR-TB with the potential for global reach.

Conservation in gene encoding Mycobacterium tuberculosis antigen Rv2660 and a high predicted population coverage of H56 multistage vaccine in South Africa

H56/AERAS-456+IC31 (H56), composed of two early secretion proteins, Ag85B and ESAT-6, and a latency associated protein, Rv2660, and the IC31 Intercell adjuvant, is a new fusion subunit vaccine candidate designed to induce immunity against both new infection and reactivation of latent tuberculosis infection. Efficacy of subunit vaccines may be affected by the diversity of vaccine antigens among clinical strains and the extent of recognition by the diverse HLA molecules in the recipient population. Although a previous study showed the conservative nature of Ag85B- and ESAT-6-encoding genes, genetic diversity of Rv2660c that encodes RV2660 is largely unknown. The population coverage of H56 as a whole yet remains to be assessed. The present study was conducted to address these important knowledge gaps. DNA sequence analysis of Rv2660c found no variation among 83 of the 84 investigated clinical strains belonging to four genetic lineages. H56 was predicted to have as high as 99.6% population coverage in the South Africa population using the Immune Epitope Database (IEDB) Population Coverage Tool. Further comparison of H56 population coverage between South African Blacks and Caucasians based on the phenotypic frequencies of binding MHC Class I and Class II supertype alleles found that all of the nine MHC-I and six of eight MHC-II human leukocyte antigen (HLA) supertype alleles analyzed were significantly differentially expressed between the two subpopulations. This finding suggests the presence of race-specific functional binding motifs of MHC-I and MHC-II HLA alleles, which, in turn, highlights the importance of including diverse populations in vaccine clinical evaluation. In conclusion, H56 vaccine is predicted to have a promising population coverage in South Africa; this study demonstrates the utility of integrating comparative genomics and bioinformatics in bridging animal and clinical studies of novel TB vaccines.

Aqueous Molecular Dynamics Simulations of the M. tuberculosis Enoyl-ACP Reductase-NADH System and Its Complex with a Substrate Mimic or Diphenyl Ethers Inhibitors

Molecular dynamics (MD) simulations of 12 aqueous systems of the NADH-dependent enoyl-ACP reductase from Mycobacterium tuberculosis (InhA) were carried out for up to 20–40 ns using the GROMACS 4.5 package. Simulations of the holoenzyme, holoenzyme-substrate, and 10 holoenzyme-inhibitor complexes were conducted in order to gain more insight about the secondary structure motifs of the InhA substrate-binding pocket. We monitored the lifetime of the main intermolecular interactions: hydrogen bonds and hydrophobic contacts. Our MD simulations demonstrate the importance of evaluating the conformational changes that occur close to the active site of the enzyme-cofactor complex before and after binding of the ligand and the influence of the water molecules. Moreover, the protein-inhibitor total steric (E_LJ) and electrostatic (E_C) interaction energies, related to Gly96 and Tyr158, are able to explain 80% of the biological response variance according to the best linear equation, pK_i = 7.772 − 0.1885 × Gly96 + 0.0517 × Tyr158 (R² = 0.80; n = 10), where interactions with Gly96, mainly electrostatic, increase the biological response, while those with Tyr158 decrease. These results will help to understand the structure-activity relationships and to design new and more potent anti-TB drugs.

Strain classification of Mycobacterium tuberculosis isolates in Brazil based on genotypes obtained by spoligotyping, mycobacterial interspersed repetitive unit typing and the presence of large sequence and single nucleotide polymorphism

Rio de Janeiro is endemic for tuberculosis (TB) and presents the second largest prevalence of the disease in Brazil. Here, we present the bacterial population structure of 218 isolates of Mycobacterium tuberculosis, derived from 186 patients that were diagnosed between January 2008 and December 2009. Genotypes were generated by means of spoligotyping, 24 MIRU-VNTR typing and presence of fbpC¹⁰³, RD^Rio and RD174. The results confirmed earlier data that predominant genotypes in Rio de Janeiro are those of the Euro American Lineages (99%). However, we observed differences between the classification by spoligotyping when comparing to that of 24 MIRU-VNTR typing, being respectively 43.6% vs. 62.4% of LAM, 34.9% vs. 9.6% of T and 18.3% vs. 21.5% of Haarlem. Among isolates classified as LAM by MIRU typing, 28.0% did not present the characteristic spoligotype profile with absence of spacers 21 to 24 and 32 to 36 and we designated these conveniently as “LAM-like”, 79.3% of these presenting the LAM-specific SNP fbpC¹⁰³. The frequency of RD^Rio and RD174 in the LAM strains, as defined both by spoligotyping and 24 MIRU-VNTR loci, were respectively 11% and 15.4%, demonstrating that RD174 is not always a marker for LAM/RD^Rio strains. We conclude that, although spoligotyping alone is a tool for classification of strains of the Euro-American lineage, when combined with MIRU-VNTRs, SNPs and RD typing, it leads to a much better understanding of the bacterial population structure and phylogenetic relationships among strains of M. tuberculosis in regions with high incidence of TB.

Phylogenetic and population genetic analysis of Salmonella enterica subsp. enterica serovar Infantis strains isolated in Japan using whole genome sequence data

Salmonella enterica subsp. enterica serovar Infantis has been reported to be carried asymptomatically in chickens and isolated from some human patients with diarrhea. The aim of this study was to investigate the phylogeny of S. Infantis strains isolated in Japan from chicken meat, chicken egg shells, environmental samples from a grading and packing center for chicken eggs (GP center), diarrhea patients, and asymptomatic carriers based on whole genome sequence data. The S. Infantis strains were in five clusters in a phylogenetic tree reconstructed by the maximum likelihood method. The five clusters were confirmed by neighbor-net and Bayesian cluster analyses. Two of the five clusters formed a group containing all of the strains isolated from chicken meat and some of the strains isolated from diarrhea patients and asymptomatic carriers. The median-joining network reconstructed in this study showed that strains in one of these two clusters diverged from one node with similar relatively short branches, suggesting clonal dissemination of these strains. The other three clusters formed a group containing all of the strains isolated from chicken egg shells and the GP center, and the remaining strains from diarrhea patients and asymptomatic carriers. Interestingly, strains isolated from patients did not cluster in only one group, indicating that none of the S. Infantis strains in this study had significantly higher human pathogenicity. The population genetic analyses in this study showed the separation of the five clusters into two groups was concordant with the sources where the strains in the clusters were isolated. These results suggested that evolutionary groups with higher hierarchy than the clusters identified in this study may be present, although such groups could not be determined by phylogenetic, neighbor-net, and Bayesian analyses in this study. Determination of higher level S. Infantis evolutionary groups should be investigated using other types of genetic markers.

Combined species identification, genotyping, and drug resistance detection of Mycobacterium tuberculosis cultures by MLPA on a bead-based array

The population structure of Mycobacterium tuberculosis is typically clonal therefore genotypic lineages can be unequivocally identified by characteristic markers such as mutations or genomic deletions. In addition, drug resistance is mainly mediated by mutations. These issues make multiplexed detection of selected mutations potentially a very powerful tool to characterise Mycobacterium tuberculosis. We used Multiplex Ligation-dependent Probe Amplification (MLPA) to screen for dispersed mutations, which can be successfully applied to Mycobacterium tuberculosis as was previously shown. Here we selected 47 discriminative and informative markers and designed MLPA probes accordingly to allow analysis with a liquid bead array and robust reader (Luminex MAGPIX technology). To validate the bead-based MLPA, we screened a panel of 88 selected strains, previously characterised by other methods with the developed multiplex assay using automated positive and negative calling. In total 3059 characteristics were screened and 3034 (99.2%) were consistent with previous molecular characterizations, of which 2056 (67.2%) were directly supported by other molecular methods, and 978 (32.0%) were consistent with but not directly supported by previous molecular characterizations. Results directly conflicting or inconsistent with previous methods, were obtained for 25 (0.8%) of the characteristics tested. Here we report the validation of the bead-based MLPA and demonstrate its potential to simultaneously identify a range of drug resistance markers, discriminate the species within the Mycobacterium tuberculosis complex, determine the genetic lineage and detect and identify the clinically most relevant non-tuberculous mycobacterial species. The detection of multiple genetic markers in clinically derived Mycobacterium tuberculosis strains with a multiplex assay could reduce the number of TB-dedicated screening methods needed for full characterization. Additionally, as a proportion of the markers screened are specific to certain Mycobacterium tuberculosis lineages each profile can be checked for internal consistency. Strain characterization can allow selection of appropriate treatment and thereby improve treatment outcome and patient management.

Prokaryotic systematics in the genomics era

As an essential and basic biological discipline, prokaryotic systematics is entering the era of genomics. This paradigmatic shift is significant not only for understanding molecular phylogeny at the whole genome level but also in revealing the genetic or epigenetic basis that accounts for the phenotypic criteria used to classify and identify species. These developments provide an opportunity and a challenge for systematists to reanalyze the molecular mechanisms underlying the taxonomic characteristics of prokaryotes by drawing the knowledge from studies of genomics and/or functional genomics employing platform technologies and related bioinformatics tools. It is expected that taxonomic books, such as Bergey's Manual of Systematic Bacteriology may evolve into a systematics library indexed by phylogenomic information with an comprehensive understanding of prokaryotic speciation and associated increasing knowledge of biological phenomena.

Genotypic analysis of genes associated with transmission and drug resistance in the Beijing lineage of Mycobacterium tuberculosis

The Beijing genotype of Mycobacterium tuberculosis is an endemic lineage in East Asia that has disseminated worldwide. It is a major health concern, as it is geographically widespread and is considered to be hypervirulent. To elucidate its genetic diversity in Taiwan, phylogenetic reconstruction was performed using 338 M. tuberculosis Beijing family clinical isolates. Region-of-difference analysis revealed the strains from Taiwan to be distributed among six subgroups of a phylogenetic tree. Synonymous single nucleotide polymorphisms at 10 chromosomal positions were also analysed. Among the 338 isolates analysed for single-nucleotide polymorphisms by using mass spectrometry, the most frequent strain found was ST10 (53.3%), followed by ST19 (14.8%) and ST22 (14.5%). Tests of drug resistance showed that the sublineages ST10, ST19 and ST26 were over-represented in the multidrug-resistant population. The presence of mutations in putative genes coding for DNA repair enzymes, which could confer a mutator phenotype to facilitate spreading of the pathogen, did not demonstrate an association with multidrug resistance. Therefore, the DNA repair genes may be involved in transmission but not in drug resistance.

Single-nucleotide polymorphism in the fadD28 gene as a genetic marker for East Asia Lineage Mycobacterium tuberculosis

Mycobacterium tuberculosis isolates with a region of difference 105 (RD105) deletion, mainly Beijing family spoligotypes, were phylogenetically grouped into the East Asia lineage. We identified a single nucleotide polymorphism in codon 507, ATC to ATT, of the fadD28 gene as a robust marker and developed a rapid assay for East Asia lineage M. tuberculosis.