Pathogen lineage-based genome-wide association study identified CD53 as susceptible locus in tuberculosis.

Tuberculosis (TB) is known to be affected by host genetic factors. We reported a specific genetic risk factor through a genome-wide association study (GWAS) that focused on young age onset TB. In this study, we further focused on the heterogeneity of Mycobacterium tuberculosis (M. tb) lineages and assessed its possible interaction with age at onset on host genetic factors. We identified the pathogen lineage in 686 Thai TB cases and GWAS stratified by both infected pathogen lineage information and age at onset revealed a genome-wide significant association of one single-nucleotide polymorphism (SNP) on chromosome 1p13, which was specifically associated with non-Beijing lineage-infected old age onset cases (P=2.54E-08, OR=1.74 (95% CI=1.43-2.12)), when we compared them to the population-matched healthy controls. This SNP locates near the CD53 gene, which encodes a leukocyte surface glycoprotein. Interestingly, the expression of CD53 was also correlated with the patients' active TB status. This is the first report of a pathogen lineage-based genome-wide association study. The results suggested that host genetic risk in TB is depended upon the pathogen genetic background and demonstrate the importance of analyzing the interaction between host and pathogen genomes in TB.

Genetic diversity and dynamic distribution of Mycobacterium tuberculosis isolates causing pulmonary and extrapulmonary tuberculosis in Thailand.

This study examined the genetic diversity and dynamicity of circulating Mycobacterium tuberculosis strains in Thailand using nearly neutral molecular markers. The single nucleotide polymorphism (SNP)-based genotypes of 1,414 culture-positive M. tuberculosis isolates from 1,282 pulmonary tuberculosis (PTB) and 132 extrapulmonary TB (EPTB) patients collected from 1995 to 2011 were characterized. Among the eight SNP cluster groups (SCG), SCG2 (44.1%), which included the Beijing (BJ) genotype, and SCG1 (39.4%), an East African Indian genotype, were dominant. Comparisons between the genotypes of M. tuberculosis isolates causing PTB and EPTB in HIV-negative cases revealed similar prevalence trends although genetic diversity was higher in the PTB patients. The identification of 10 reported sequence types (STs) and three novel STs was hypothesized to indicate preferential expansion of the SCG2 genotype, especially the modern BJ ST10 (15.6%) and ancestral BJ ST19 (13.1%). An association between SCG2 and SCG1 genotypes and particular patient age groups implies the existence of different genetic advantages among the bacterial populations. The results revealed that increasing numbers of young patients were infected with M. tuberculosis SCGs 2 and 5, which contrasts with the reduction of the SCG1 genotype. Our results indicate the selection and dissemination of potent M. tuberculosis genotypes in this population. The determination of heterogeneity and dynamic population changes of circulating M. tuberculosis strains in countries using the Mycobacterium bovis BCG (bacillus Calmette-Guérin) vaccine are beneficial for vaccine development and control strategies.

Genome-wide host-pathogen analyses reveal genetic interaction points in tuberculosis disease.

The genetics underlying tuberculosis (TB) pathophysiology are poorly understood. Human genome-wide association studies have failed so far to reveal reproducible susceptibility loci, attributed in part to the influence of the underlying Mycobacterium tuberculosis (Mtb) bacterial genotype on the outcome of the infection. Several studies have found associations of human genetic polymorphisms with Mtb phylo-lineages, but studies analysing genome-genome interactions are needed. By implementing a phylogenetic tree-based Mtb-to-human analysis for 714 TB patients from Thailand, we identify eight putative genetic interaction points (P < 5 × 10-8) including human loci DAP and RIMS3, both linked to the IFNγ cytokine and host immune system, as well as FSTL5, previously associated with susceptibility to TB. Many of the corresponding Mtb markers are lineage specific. The genome-to-genome analysis reveals a complex interactome picture, supports host-pathogen adaptation and co-evolution in TB, and has potential applications to large-scale studies across many TB endemic populations matched for host-pathogen genomic diversity.

Pulsed field gel electrophoresis analysis of Vibrio cholerae isolates in southern Thailand.

Forty isolates of V. cholorae O1, O139 and non-O1/non-O139 collected from outbreaks in Songkhla and Phuket Provinces of southern Thailand during 1999-2001 and sporadic cases from different regions of Thailand during 1993-2002 were characterized using pulsed field gel electrophoresis (PFGE). Digestion of chromosomal DNA of the V cholerae isolates with restriction endonuclease NotI, followed by PFGE, generated 10 distinct restriction endonuclease analysis patterns consisting of 8 to 13 bands, ranging in size from 78 to 394 kb. PFGE patterns of O1 Inaba strains from the outbreak in Songkhla were identical (P1) except one isolate (P3). The O1 Inaba outbreak strains from Phuket in the same period belonged to P2 pattern, whereas the O1 Ogawa strain from the outbreak in Phuket isolated in 1999 was of P7 pattern. These patterns of O1 Inaba and Ogawa strains were slightly different suggesting that the isolates were epidemiologically related and therefore the outbreaks were likely due to the same V cholerae clone. Isolates of V cholerae O1 Inaba from sporadic cases in the neighboring area (e.g., Pattani Province) in a similar period of time of the outbreak in Songkhla Province had very similar patterns, with only one single band different from those of the outbreak isolates. This indicates that the Inaba strains isolated from Songkhla Province during the 2001 cholera outbreak belonging to P1 pattern had not spread to other regions in 2001 and 2002. On the otherhand, the sporadic isolates collected from other regions of Thailand were quite distinct from the outbreak isolates in Songkhla Province, especially those from Chaiyaphum and Chaing Mai Provinces, which belonged to P5 and P6 pattern, respectively. Isolates of V cholerae O139 and non-O1/non-O139 gave different patterns from that of V. cholerae O1. This study shows that the PFGE technique is markedly advantageous in distinguishing strains of V cholerae isolates leading to insightful detailed charateristics of these isolates in Thailand.

Identification and in silico functional prediction of lineage-specific SNPs distributed in DosR-related proteins and resuscitation-promoting factor proteins of Mycobacterium tuberculosis.

One-third of the world population is infected by Mycobacterium tuberculosis, which may persist in the latent or dormant state. Bacteria can shift to dormancy when encountering harsh conditions such as low oxygen, nutrient starvation, high acidity and host immune defenses. Genes related to the dormancy survival regulator (DosR) regulon are responsible for the inhibition of aerobic respiration and replication, which is required to enter dormancy. Conversely, resuscitation-promoting factor (rpf) proteins participate in reactivation from dormancy and the development of active tuberculosis (TB). Many DosR regulon and rpf proteins are immunodominant T cell antigens that are highly expressed in latent TB infection. They could serve as TB vaccine candidates and be used for diagnostic development. We explored the genetic polymorphisms of 50 DosR-related genes and 5 rpf genes among 1,170 previously sequenced clinical M. tuberculosis genomes. Forty-three lineage- or sublineage-specific nonsynonymous single nucleotide polymorphisms (nsSNPs) were identified. Ten nsSNPs were specific to all Mtb isolates belonging to lineage 1 (L1). Two common sublineages, the Beijing family (L2.2) and EAI2 (L1.2.1), differed at as many as 26 lineage- or sublineage-specific SNPs. DosR regulon genes related to membrane proteins and the rpf family possessed mean dN/dS ratios greater than one, suggesting that they are under positive selection. Although the T cell epitope regions of DosR-related and rpf antigens were quite conserved, we found that the epitopes in L1 had higher rates of genetic polymorphisms than the other lineages. Some mutations in immunogenic epitopes of the antigens were specific to particular M. tuberculosis lineages. Therefore, the genetic diversity of the DosR regulon and rpf proteins might impact the adaptation of M. tuberculosis to the dormant state and the immunogenicity of latency antigens, which warrants further investigation.

Homoplastic single nucleotide polymorphisms contributed to phenotypic diversity in Mycobacterium tuberculosis.

Homoplastic mutations are mutations independently occurring in different clades of an organism. The homoplastic changes may be a result of convergence evolution due to selective pressures. Reports on the analysis of homoplastic mutations in Mycobacterium tuberculosis have been limited. Here we characterized the distribution of homoplastic single nucleotide polymorphisms (SNPs) among genomes of 1,170 clinical M. tuberculosis isolates. They were present in all functional categories of genes, with pe/ppe gene family having the highest ratio of homoplastic SNPs compared to the total SNPs identified in the same functional category. Among the pe/ppe genes, the homoplastic SNPs were common in a relatively small number of homologous genes, including ppe18, the protein of which is a component of a promising candidate vaccine, M72/AS01E. The homoplastic SNPs in ppe18 were particularly common among M. tuberculosis Lineage 1 isolates, suggesting the need for caution in extrapolating the results of the vaccine trial to the population where L1 is endemic in Asia. As expected, homoplastic SNPs strongly associated with drug resistance. Most of these mutations are already well known. However, a number of novel mutations associated with streptomycin resistance were identified, which warrants further investigation. A SNP in the intergenic region upstream of Rv0079 (DATIN) was experimentally shown to increase transcriptional activity of the downstream gene, suggesting that intergenic homoplastic SNPs should have effects on the physiology of the bacterial cells. Our study highlights the potential of homoplastic mutations to produce phenotypic changes. Under selective pressure and during interaction with the host, homoplastic mutations may confer advantages to M. tuberculosis and deserve further characterization.

Risk factors associated with large clusters of tuberculosis patients determined by whole-genome sequencing in a high-tuberculosis-burden country.

Whole-genome sequencing (WGS) analysis has great discriminative power for detecting similar molecular fingerprints of suspected tuberculosis (TB) clusters. The proportion of TB cases within clusters and the associated risk factors are important epidemiological parameters guiding appropriate outbreak control strategies in endemic settings. We conducted a hospital-based TB case-cohort study between 2003 and 2011 in the northernmost province of Thailand. We identified TB clusters by Mycobacterium tuberculosis WGS and analysed the risks of TB clustering and the characteristics of large clusters compared with small clusters. Among 1146 TB isolates, we identified 77 clusters with 251 isolates defined by a 5-single-nucleotide variant (SNV) cutoff and 112 clusters with 431 isolates defined by a 12-SNV cutoff. Twelve large clusters with 6 isolates or more in each cluster were identified by a 12-SNV cutoff. Sublineage 2.2.1 (both Ancestral and Modern) strains and imprisonment were independently associated with large clusters. Furthermore, although large clusters of Lineage 2.2.1/Ancestral strains included a high number of prisoners, Lineage 2.2.1/Modern strain clusters were only associated with treatment failures and drug resistance. Heterogeneity among lineage strains was observed with respect to large-cluster characteristics. Patients with an increased TB-transmission tendency should be priority targets for contact investigations and outbreak interventions to prevent ongoing transmission.

A novel Ancestral Beijing sublineage of Mycobacterium tuberculosis suggests the transition site to Modern Beijing sublineages.

Global Mycobacterium tuberculosis population comprises 7 major lineages. The Beijing strains, particularly the ones classified as Modern groups, have been found worldwide, frequently associated with drug resistance, younger ages, outbreaks and appear to be expanding. Here, we report analysis of whole genome sequences of 1170 M. tuberculosis isolates together with their patient profiles. Our samples belonged to Lineage 1-4 (L1-L4) with those of L1 and L2 being equally dominant. Phylogenetic analysis revealed several new or rare sublineages. Differential associations between sublineages of M. tuberculosis and patient profiles, including ages, ethnicity, HIV (human immunodeficiency virus) infection and drug resistance were demonstrated. The Ancestral Beijing strains and some sublineages of L4 were associated with ethnic minorities while L1 was more common in Thais. L2.2.1.Ancestral 4 surprisingly had a mutation that is typical of the Modern Beijing sublineages and was common in Akha and Lahu tribes who have migrated from Southern China in the last century. This may indicate that the evolutionary transition from the Ancestral to Modern Beijing sublineages might be gradual and occur in Southern China, where the presence of multiple ethnic groups might have allowed for the circulations of various co-evolving sublineages which ultimately lead to the emergence of the Modern Beijing strains.

Evidence for Host-Bacterial Co-evolution via Genome Sequence Analysis of 480 Thai Mycobacterium tuberculosis Lineage 1 Isolates.

Tuberculosis presents a global health challenge. Mycobacterium tuberculosis is divided into several lineages, each with a different geographical distribution. M. tuberculosis lineage 1 (L1) is common in the high-burden areas in East Africa and Southeast Asia. Although the founder effect contributes significantly to the phylogeographic profile, co-evolution between the host and M. tuberculosis may also play a role. Here, we reported the genomic analysis of 480 L1 isolates from patients in northern Thailand. The studied bacterial population was genetically diverse, allowing the identification of a total of 18 sublineages distributed into three major clades. The majority of isolates belonged to L1.1 followed by L1.2.1 and L1.2.2. Comparison of the single nucleotide variant (SNV) phylogenetic tree and the clades defined by spoligotyping revealed some monophyletic clades representing EAI2_MNL, EAI2_NTM and EAI6_BGD1 spoligotypes. Our work demonstrates that ambiguity in spoligotype assignment could be partially resolved if the entire DR region is investigated. Using the information to map L1 diversity across Southeast Asia highlighted differences in the dominant strain-types in each individual country, despite extensive interactions between populations over time. This finding supported the hypothesis that there is co-evolution between the bacteria and the host, and have implications for tuberculosis disease control.

Evolution of some variable-number tandem repeat loci among a group of Beijing strains of Mycobacterium tuberculosis.

The patterns of variable-number tandem repeats (VNTR) genotypes of a clonal group of Mycobacterium tuberculosis Beijing isolates with very similar IS6110-restriction fragment length polymorphism (RFLP) patterns were studied. Differences between VNTR were mostly by a single repeat unit. However, a multiple-unit change also occurred. This suggests that a mechanism other than the slipped-strand mispairing might be responsible for the instability of VNTR in M. tuberculosis as well. This finding is useful for inferring phylogenetics of M. tuberculosis based on the VNTR genotypes.