rpoB Genotypes of Mycobacterium tuberculosis Beijing Family Isolates from East Asian Countries

The 81-bp region of the rpoB gene in 66 Rif(r) Mycobacterium tuberculosis isolates from China, Japan, Korea, and Taiwan was analyzed. Twelve single-nucleotide substitutions in the rpoB gene were detected. The most prevalent mutations were at Ser-531 (52%), Asp-516 (17%), and His-526 (11%). Mutations were not found in seven (11%) of the isolates. Higher mutation rates in 50 Beijing family isolates were found than in other isolates for mutations at Asp-516 (18 and 12.5%, respectively) and His-526 (12 and 6.3%, respectively). The different rates of mutation may reflect the choice of rifamycin analogs

Characterization of the Manila Family of Mycobacterium tuberculosis

Forty-eight Mycobacterium tuberculosis strains were obtained from patients living in metropolitan Manila, Republic of the Philippines. Three molecular typing methods, IS6110 restriction fragment length polymorphism, spoligotyping, and DNA sequencing of the oxyR, gyrA, and katG loci, established that these strains have restricted diversity and are members of a related genetic group of organisms. Comparison of the DNA fingerprint patterns with those in international databases confirmed the uniqueness of this group of isolates, which we designate the Manila family of M. tuberculosis

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

Genomic analyses of the ancestral Manila family of <i>Mycobacterium tuberculosis</i>

<div><p>With its airborne transmission and prolonged latency period, <i>Mycobacterium tuberculosis</i> spreads worldwide as one of the most successful bacterial pathogens and continues to kill millions of people every year. <i>M</i>. <i>tuberculosis</i> lineage 1 is inferred to originate ancestrally based on the presence of the 52-bp TbD1 sequence and analysis of single nucleotide polymorphisms. Previously, we briefly reported the complete genome sequencing of <i>M</i>. <i>tuberculosis</i> strains 96121 and 96075, which belong to the ancient Manila family and modern Beijing family respectively. Here we present the comprehensive genomic analyses of the Manila family in lineage 1 compared to complete genomes in lineages 2–4. Principal component analysis of the presence and absence of CRISPR spacers suggests that Manila isolate 96121 is distinctly distant from lineages 2–4. We further identify a truncated <i>whiB5</i> gene and a putative operon consisting of genes encoding a putative serine/threonine kinase PknH and a putative ABC transporter, which are only found in the genomes of Manila family isolates. Six single nucleotide polymorphisms are uniquely conserved in 38 Manila strains. Moreover, when compared to <i>M</i>. <i>tuberculosis</i> H37Rv, 59 proteins are under positive selection in Manila family isolate 96121 but not in Beijing family isolate 96075. The unique features further serve as biomarkers for Manila strains and may shed light on the limited transmission of this ancestral lineage outside of its Filipino host population.</p></div

Genomic sequencing is required for identification of tuberculosis transmission in Hawaii

Abstract Background Tuberculosis (TB) caused an estimated 1.4 million deaths and 10.4 million new cases globally in 2015. TB rates in the United States continue to steadily decline, yet rates in the State of Hawaii are perennially among the highest in the nation due to a continuous influx of immigrants from the Western Pacific and Asia. TB in Hawaii is composed of a unique distribution of genetic lineages, with the Beijing and Manila families of Mycobacterium tuberculosis (Mtb) comprising over two-thirds of TB cases. Standard fingerprinting methods (spoligotyping plus 24-loci Mycobacterial Interspersed Repetitive Units-Variable Number Tandem Repeats [MIRU-VNTR] fingerprinting) perform poorly when used to identify actual transmission clusters composed of isolates from these two families. Those typing methods typically group isolates from these families into large clusters of non-linked isolates with identical fingerprints. Next-generation whole-genome sequencing (WGS) provides a new tool for molecular epidemiology that can resolve clusters of isolates with identical spoligotyping and MIRU-VNTR fingerprints. Methods We performed WGS and SNP analysis and evaluated epidemiological data to investigate 19 apparent TB transmission clusters in Hawaii from 2003 to 2017 in order to assess WGS’ ability to resolve putative Mtb clusters from the Beijing and Manila families. This project additionally investigated MIRU-VNTR allele prevalence to determine why standard Mtb fingerprinting fails to usefully distinguish actual transmission clusters from these two Mtb families. Results WGS excluded transmission events in seven of these putative clusters, confirmed transmission in eight, and identified both transmission-linked and non-linked isolates in four. For epidemiologically identified clusters, while the sensitivity of MIRU-VNTR fingerprinting for identifying actual transmission clusters was found to be 100%, its specificity was only 28.6% relative to WGS. We identified that the Beijing and Manila families’ significantly lower Shannon evenness of MIRU-VNTR allele distributions than lineage 4 was the cause of standard fingerprinting’s poor performance when identifying transmission in Beijing and Manila family clusters. Conclusions This study demonstrated that WGS is necessary for epidemiological investigation of TB in Hawaii and the Pacific

Strain information list.

<p>Strain information list.</p

CRISPR spacer rearrangements in Manila family 96121, Beijing family 96075, and H37Rv.

<p>Green highlights spacers conserved with other MTBC species but lost in modern Beijing family isolate 96075 and H37Rv. Red highlights spacers specific in Manila family isolate 96121. Blue highlights spacers present in H37Rv but absent in Manila family isolate 96121 and Beijing family isolate 96075.</p

List of proteins under positive selection in Manila family 96121 but not Beijing family 96075.

<p>List of proteins under positive selection in Manila family 96121 but not Beijing family 96075.</p