Rapid deletion-based subtyping system for the Manila family of Mycobacterium tuberculosis

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.

Whole genome SNP analysis suggests unique virulence factor differences of the Beijing and Manila families of Mycobacterium tuberculosis found in Hawaii

While tuberculosis (TB) remains a global disease, the WHO estimates that 62% of the incident TB cases in 2016 occurred in the WHO South-East Asia and Western Pacific regions. TB in the Pacific is composed predominantly of two genetic families of Mycobacterium tuberculosis (Mtb): Beijing and Manila. The Manila family is historically under-studied relative to the families that comprise the majority of TB in Europe and North America (e.g. lineage 4), and it remains unclear why this lineage has persisted in Filipino populations despite the predominance of more globally successful Mtb lineages in most of the world. The Beijing family is of particular interest as it is increasingly associated with drug resistance throughout the world. Both of these lineages are important to the State of Hawaii, where they comprise over two-thirds of TB cases. Here, we performed whole genome sequencing on 82 Beijing family, Manila family, and outgroup clinical Mtb isolates from Hawaii to identify lineage-specific SNPs (SNPs found in all isolates from their respective families, and exclusively in those families) in established virulence factor genes. Six non-silent lineage-specific virulence factor SNPs were found in the Beijing family, including mutations in alternative sigma factor sigG and polyketide synthases pks5 and pks7. The Manila family displayed more than eleven non-silent lineage-specific and characteristic virulence factor mutations, including in genes coding for MCE-family protein Mce1B, two mutations in fatty-acid-AMP ligase FadD26, and virulence-regulating transcriptional regulator VirS. This study further identified an ancient clade that shared some virulence factor mutations with the Manila family, and investigated the relationship of those and other “Manila-like” spoligotypes to the Manila family with this SNP dataset. This work identified a set of virulence genes that are worth pursuing to determine potential differences in transmission or virulence displayed by these two Mtb families.

Whole-genome sequencing of the Mycobacterium tuberculosis Manila sublineage results in less clustering and better resolution than mycobacterial interspersed repetitive-unit-variable-number tandem-repeat (MIRU-VNTR) typing and spoligotyping

Mycobacterium tuberculosis isolates of the Manila sublineage are genetically homogeneous. In this study, we used whole-genome sequencing (WGS) to type a collection of 36 M. tuberculosis isolates of the Manila family. WGS enabled the subtyping of these 36 isolates into at least 10 distinct clusters. Our results indicate that WGS is a powerful approach to determining the relatedness of Manila family M. tuberculosis isolates.

Clonal expansion of both modern and ancient genotypes of Mycobacterium tuberculosis in southern Taiwan

We present the first comprehensive analysis of Mycobacterium tuberculosis isolates circulating in the Kaohsiung region of southern Taiwan. The major spoligotypes found in the 224 isolates studied were Beijing lineages (n = 97; 43.3%), EAI lineages (n = 72; 32.1%) and Haarlem lineages (n = 18; 8.0%). By 24 MIRU-VNTR typing, 174 patterns were identified, including 24 clusters of 74 isolates and 150 unique patterns. The combination of spoligotyping and 12-MIRU-VNTR revealed that 129 (57.6%) of the 224 isolates were clustered in 18 genotypes. Moreover, 63.6% (7/11) of infected persons younger than 30 years had a Beijing strain, which could suggest recent spread among younger persons by this family of TB strains in Kaohsiung. Among the 94 Beijing family (SIT1, SIT250 and SIT1674) isolates further analyzed for SNPs by mass spectrometry, the most frequent strain found was ST10 (n = 49; 52%), followed by ST22 (n = 17; 18%) and ST19 (n = 11; 12%). Among the EAI-Manila family isolates analyzed by region deletion-based subtyping, the most frequent strain found was RD type 1 (n = 63; 87.5%), followed by RD type 2 (n = 9; 12.5%). In our previous study, the proportion of modern Beijing strains (52.5%) in northern Taiwan was significantly higher than the proportion of EAI strains (11%). In contrast, in the present study, EAI strains comprised up to 32% of Beijing strains in southern Taiwan. In conclusion, both ‘modern’ (Beijing) and ‘ancient’ (EAI) M. tuberculosis strains are prevalent in the Kaohsiung region, perhaps suggesting that both strains are somehow more adapted to southern Taiwan. It will be interesting to investigate the dynamics of the lineage composition by different selection pressures.