Comparative genome analysis of fish and human isolates of Mycobacterium marinum

Temporal genome-wide fitness analysis of Mycobacterium marinum during infection reveals the genetic requirement for virulence and survival in amoebae and microglial cells

Tuberculosis remains the most pervasive infectious disease and the recent emergence of drug-resistant strains emphasizes the need for more efficient drug treatments. A key feature of pathogenesis, conserved between the human pathogen Mycobacterium tuberculosis and the model pathogen Mycobacterium marinum, is the metabolic switch to lipid catabolism and altered expression of virulence genes at different stages of infection. This study aims to identify genes involved in sustaining viable intracellular infection. We applied transposon sequencing (Tn-Seq) to M. marinum, an unbiased genome-wide strategy combining saturation insertional mutagenesis and high-throughput sequencing. This approach allowed us to identify the localization and relative abundance of insertions in pools of transposon mutants. Gene essentiality and fitness cost of mutations were quantitatively compared between in vitro growth and different stages of infection in two evolutionary distinct phagocytes, the amoeba Dictyostelium discoideum and the murine BV2 microglial cells. In the M. marinum genome, 57% of TA sites were disrupted and 568 genes (10.2%) were essential, which is comparable to previous Tn-Seq studies on M. tuberculosis and M. bovis. Major pathways involved in the survival of M. marinum during infection of D. discoideum are related to DNA damage repair, lipid and vitamin metabolism, the type VII secretion system (T7SS) ESX-1, and the Mce1 lipid transport system. These pathways, except Mce1 and some glycolytic enzymes, were similarly affected in BV2 cells. These differences suggest subtly distinct nutrient availability or requirement in different host cells despite the known predominant use of lipids in both amoeba and microglial cells.IMPORTANCEThe emergence of biochemically and genetically tractable host model organisms for infection studies holds the promise to accelerate the pace of discoveries related to the evolution of innate immunity and the dissection of conserved mechanisms of cell-autonomous defenses. Here, we have used the genetically and biochemically tractable infection model system Dictyostelium discoideum/Mycobacterium marinum to apply a genome-wide transposon-sequencing experimental strategy to reveal comprehensively which mutations confer a fitness advantage or disadvantage during infection and compare these to a similar experiment performed using the murine microglial BV2 cells as host for M. marinum to identify conservation of virulence pathways between hosts.

Mycobacterial formation of intracellular lipid inclusions is a dynamic process associated with rapid replication

Intracellular lipid inclusions (ILI) are triacylglyceride rich organelles produced by mycobacteria thought to serve as energy reservoirs. It is believed that ILI are formed as a result of a dosR mediated transition from replicative growth to non-replicating persistence (NRP). ILI rich Mycobacterium tuberculosis (Mtb) bacilli have been reported during infection and in sputum, establishing their importance in Mtb pathogenesis. Studies conducted in mycobacteria such as Mycobacterium smegmatis, Mycobacterium abscessus, or lab Mtb strains have demonstrated ILI formation in the presence of hypoxic, nitric oxide, nutrient limitation, or low nitrogen stress, conditions believed to emulate the host environment within which Mtb resides. Here, we show that M. marinum and clinical Mtb isolates make ILI during active replication in axenic culture independent of environmental stressors. By tracking ILI formation dynamics we demonstrate that ILI are quickly formed in the presence of fresh media or exogenous fatty acids but are rapidly depleted while bacteria are still actively replicating. We also show that the cell envelope is an alternate site for neutral lipid accumulation observed during stationary phase. In addition, we screen a panel of 60 clinical isolates and observe variation in ILI production during early log phase growth between and among Mtb lineages. Finally, we show that dosR expression level does not strictly correlate with ILI accumulation in fresh clinical isolates. Taken together, our data provide evidence of an active ILI formation pathway in replicating mycobacteria cultured in the absence of stressors, suggesting a decoupling of ILI formation from NRP.

Genomic Degeneration and Reduction in the Fish Pathogen Mycobacterium shottsii

Mycobacterium shottsii is a dysgonic, nonpigmented mycobacterium originally isolated from diseased striped bass (Morone saxatilis) in the Chesapeake Bay, USA. Genomic analysis reveals that M. shottsii is a Mycobacterium ulcerans/Mycobacterium marinum clade (MuMC) member, but unlike the superficially similar M. pseudoshottsii, also isolated from striped bass, it is not an M. ulcerans ecovar, instead belonging to a transitional group of strains basal to proposed "Aronson" and "M" lineages. Although phylogenetically distinct from the human pathogen M. ulcerans, the M. shottsii genome shows parallel but nonhomologous genomic degeneration, including massive accumulation of pseudogenes accompanied by proliferation of unique insertion sequences (ISMysh01, ISMysh03), large-scale deletions, and genomic reorganization relative to typical M. marinum strains. Coupled with its observed ecological characteristics and loss of chromogenicity, the genomic structure of M. shottsii is suggestive of evolution toward a state of obligate pathogenicity, as observed for other Mycobacterium spp., including M. ulcerans, M. tuberculosis, and M. leprae. IMPORTANCE Morone saxatilis (striped bass) is an ecologically and economically important finfish species on the United States east coast. Mycobacterium shottsii and Mycobacterium pseudoshottsii were originally described in the early 2000s as novel species from outbreaks of visceral and dermal mycobacteriosis in this species. Biochemical and genetic characterization place these species within the Mycobacterium ulcerans/M. marinum clade (MuMC), and M. pseudoshottsii has been proposed as an ecovar of M. ulcerans. Here, we describe the complete genome of M. shottsii, demonstrating that it is clearly not an M. ulcerans ecovar; however, it has undergone parallel genomic modification suggestive of a transition to obligate pathogenicity. As in M. ulcerans, the M. shottsii genome demonstrates widespread pseudogene formation driven by proliferation of insertion sequences, as well as genomic reorganization. This work clarifies the phylogenetic position of M. shottsii relative to other MuMC members and provides insight into processes shaping its genomic structure.

Understanding the transmission of Mycobacterium ulcerans: A step towards controlling Buruli ulcer

Mycobacterium ulcerans is the causative agent of Buruli ulcer, a rare but chronic debilitating skin and soft tissue disease found predominantly in West Africa and Southeast Australia. While a moderate body of research has examined the distribution of M. ulcerans, the specific route(s) of transmission of this bacterium remain unknown, hindering control efforts. M. ulcerans is considered an environmental pathogen given it is associated with lentic ecosystems and human-to-human spread is negligible. However, the pathogen is also carried by various mammals and invertebrates, which may serve as key reservoirs and mechanical vectors, respectively. Here, we examine and review recent evidence from these endemic regions on potential transmission pathways, noting differences in findings between Africa and Australia, and summarising the risk and protective factors associated with Buruli ulcer transmission. We also discuss evidence suggesting that environmental disturbance and human population changes precede outbreaks. We note five key research priorities, including adoption of One Health frameworks, to resolve transmission pathways and inform control strategies to reduce the spread of Buruli ulcer.

Buruli ulcer is a debilitating skin and soft tissue disease characterised by large ulcerative wounds that are treated with antibiotics or with adjunctive surgery for advanced cases. Found predominantly in West Africa and Southeast Australia, the causative agent is the environmental bacterial pathogen Mycobacterium ulcerans. Lack of understanding of transmission pathways, combined with the absence of a vaccine, has hindered efforts to control the spread of M. ulcerans. Here, in order to identify probable transmission pathways and inform future studies, we review literature linking M. ulcerans to environmental reservoirs, mammalian hosts, and potential invertebrate vectors. We also summarise factors and behaviours that reduce the risk of developing Buruli ulcer, to inform effective prevention strategies and further shed light on transmission pathways.

Mycobacterium marinum

Mycobacterium marinum is a well-known pathogenic mycobacterium for skin and soft tissue infections and is associated with fishes and water. Among nontuberculous mycobacteria (NTM), it is the leading cause of extrarespiratory human infections worldwide. In addition, there is a specific scientific interest in M. marinum because of its genetic relatedness to Mycobacterium tuberculosis and because experimental infection of M. marinum in fishes mimics tuberculosis pathogenesis. Microbiological characteristics include the fact that it grows in 7 to 14 days with photochromogenic colonies and is difficult to differentiate from Mycobacterium ulcerans and other mycolactone-producing NTM on a molecular basis. The diagnosis is highly suspected by the mode of infection, which is related to the hobby of fishkeeping, professional handling of marine shells, or swimming in nonchlorinated pools. Clinics distinguished skin and soft tissue lesions (typically sporotrichoid or subacute hand nodules) and lesions disseminated to joint and bone, often related with the local use of corticosteroids. In clinical microbiology, microscopy and culture are often negative because growth requires low temperature (30°C) and several weeks to succeed in primary cultivation. The treatment is not standardized, and no randomized control trials have been done. Therapy is a combination of surgery and antimicrobial agents such as cyclines and rifampin, with successful outcome in most of the skin diseases but less frequently in deep tissue infections. Prevention can be useful with hand protection recommendations for professionals and all persons manipulating fishes or fish tank water and use of alcohol disinfection after contact.

Draft Genome Sequence of the Fish Pathogen Mycobacterium pseudoshottsii Strain JCM15466, a Species Closely Related to M. marinum

Mycobacterium pseudoshottsii is a slowly growing photochromogenic mycobacterium and fish pathogen isolated from wild marine fishes. M. pseudoshottsii closely resembles M. marinum, which is a human and animal pathogen. Here, we report the draft genome sequence of M. pseudoshottsii strain JCM15466, originally isolated from striped bass, Morone saxatilis.

Whole genome analyses of marine fish pathogenic isolate, Mycobacterium sp. 012931

Mycobacterium is a genus within the order Actinomycetales that comprises of a large number of well-characterized species, several of which includes pathogens known to cause serious disease in human and animal. Here, we report the whole genome sequence of Mycobacterium sp. strain 012931 isolated from the marine fish, yellowtail (Seriola quinqueradiata). Mycobacterium sp. 012931 is a fish pathogen causing serious damage to aquaculture farms in Japan. DNA dot plot analysis showed that Mycobacterium sp. 012931 was more closely related to Mycobacterium marinum when compared across several Mycobacterium species. However, little conservation of the gene order was observed between Mycobacterium sp. 012931 and M. marinum genome. The annotated 5,464 genes of Mycobacterium sp. 012931 was classified into 26 subsystems. The insertion/deletion gene analysis shows Mycobacterium sp. 012931 had 643 unique genes that were not found in the M. marinum strains. In the virulence, disease, and defense subsystem, both insertion and deletion genes of Mycobacterium sp. 012931 were associated with the PPE gene cluster of Mycobacteria. Of seven plcB genes in Mycobacterium sp. 012931, plcB_2 and plcB_3 showed low identities with those of M. marinum strains. Therefore, Mycobacterium sp. 012931 has differences on genetic and virulence from M. marinum and may induce different interaction mechanisms between host and pathogen.

Unusual large-scale chromosomal rearrangements in Mycobacterium tuberculosis Beijing B0/W148 cluster isolates

The Mycobacterium tuberculosis (MTB) Beijing family isolates are geographically widespread, and there are examples of Beijing isolates that are hypervirulent and associated with drug resistance. One-fourth of Beijing genotype isolates found in Russia belong to the B0/W148 group. The aim of the present study was to investigate features of these endemic strains on a genomic level. Four Russian clinical isolates of this group were sequenced, and the data obtained was compared with published sequences of various MTB strain genomes, including genome of strain W-148 of the same B0/W148 group. The comparison of the W-148 and H37Rv genomes revealed two independent inversions of large segments of the chromosome. The same inversions were found in one of the studied strains after deep sequencing using both the fragment and mate-paired libraries. Additionally, inversions were confirmed by RFLP hybridization analysis. The discovered rearrangements were verified by PCR in all four newly sequenced strains in the study and in four additional strains of the same Beijing B0/W148 group. The other 32 MTB strains from different phylogenetic lineages were tested and revealed no inversions. We suggest that the initial largest inversion changed the orientation of the three megabase (Mb) segment of the chromosome, and the second one occurred in the previously inverted region and partly restored the orientation of the 2.1 Mb inner segment of the region. This is another remarkable example of genomic rearrangements in the MTB in addition to the recently published of large-scale duplications. The described cases suggest that large-scale genomic rearrangements in the currently circulating MTB isolates may occur more frequently than previously considered, and we hope that further studies will help to determine the exact mechanism of such events.

Transcriptome profiling analysis on whole bodies of microbial challenged Eriocheir sinensis larvae for immune gene identification and SNP development

To study crab immunogenetics of individuals, newly hatched Eriocheir sinensis larvae were stimulated with a mixture of three pathogen strains (Gram-positive bacteria Micrococcus luteus, Gram-negative bacteria Vibrio alginolyticus and fungi Pichia pastoris; 10⁸ cfu·mL^-1). A total of 44,767,566 Illumina clean reads corresponding to 4.52 Gb nucleotides were generated and assembled into 100,252 unigenes (average length: 1,042 bp; range: 201-19,357 bp). 17,097 (26.09%) of 65,535 non-redundant unigenes were annotated in NCBI non-redundant protein (Nr) database. Moreover, 23,188 (35.38%) unigenes were assigned to three Gene Ontology (GO) categories, 15,071 (23.00%) to twenty-six Clusters of orthologous Groups (COG) and 8,574 (13.08%) to six Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways, respectively. Numerous genes were further identified to be associated with multiple immune pathways, including Toll, immune deficiency (IMD), janus kinase (JAK)-signal transducers and activators of transcription (STAT) and mitogen-activated protein kinase (MAPK) pathways. Some of them, such as tumor necrosis factor receptor associated factor 6 (TRAF6), fibroblast growth factor (FGF), protein-tyrosine phosphatase (PTP), JNK-interacting protein 1 (JIP1), were first identified in E. sinensis. TRAF6 was even first discovered in crabs. Additionally, 49,555 single nucleotide polymorphisms (SNPs) were developed from over 13,309 unigenes. This is the first transcriptome report of whole bodies of E. sinensis larvae after immune challenge. Data generated here not only provide detail information to identify novel genes in genome reference-free E. sinensis, but also facilitate our understanding on host immunity and defense mechanism of the crab at whole transcriptome level.

Whole-Genome Sequence of Fish-Pathogenic Mycobacterium sp. Strain 012931, Isolated from Yellowtail (Seriola quinqueradiata)

The genus Mycobacterium comprises a large number of well-characterized species, several of which are human and animal pathogens. Here, we report the whole-genome sequence of Mycobacterium sp. strain 012931, a fish pathogen responsible for huge losses in aquaculture farms in Japan. The strain was isolated from a marine fish, yellowtail (Seriola quinqueradiata).