Bacterial Classification of Fish-Pathogenic Mycobacterium Species by Multigene Phylogenetic Analyses and MALDI Biotyper Identification System

Mycobacterium marinum is difficult to distinguish from other species of Mycobacterium isolated from fish using biochemical methods. Here, we used genetic and proteomic analyses to distinguish three Mycobacterium strains: M. marinum strains MB2 and Europe were isolated from tropical and marine fish in Thailand and Europe, and Mycobacterium sp. 012931 strain was isolated from yellowtail in Japan. In phylogenetic trees based on gyrB, rpoB, and Ag85B genes, Mycobacterium sp. 012931 clustered with M. marinum strains MB2 and Europe, but in trees based on 16S rRNA, hsp65, and Ag85A genes Mycobacterium sp. 012931 did not cluster with the other strains. In proteomic analyses using a Bruker matrix-assisted laser desorption ionization Biotyper, the mass profile of Mycobacterium sp. 012931 differed from the mass profiles of the other two fish M. marinum strains. Therefore, Mycobacterium sp. 012931 is similar to M. marinum but is not the same, suggesting that it could be a subspecies of M. marinum.     

Strain Variation in Mycobacteriummarinum Fish Isolates

A molecular characterization of two Mycobacterium marinum genes, 16S rRNA and hsp65, was carried out with a total of 21 isolates from various species of fish from both marine and freshwater environments of Israel, Europe, and the Far East. The nucleotide sequences of both genes revealed that all M. marinum isolates from fish in Israel belonged to two different strains, one infecting marine (cultured and wild) fish and the other infecting freshwater (cultured) fish. A restriction enzyme map based on the nucleotide sequences of both genes confirmed the divergence of the Israeli marine isolates from the freshwater isolates and differentiated the Israeli isolates from the foreign isolates, with the exception of one of three Greek isolates from marine fish which was identical to the Israeli marine isolates. The second isolate from Greece exhibited a single base alteration in the 16S rRNA sequence, whereas the third isolate was most likely a new Mycobacterium species. Isolates from Denmark and Thailand shared high sequence homology to complete identity with reference strain ATCC 927. Combined analysis of the two gene sequences increased the detection of intraspecific variations and was thus of importance in studying the taxonomy and epidemiology of this aquatic pathogen. Whether the Israeli M. marinum strain infecting marine fish is endemic to the Red Sea and found extremely susceptible hosts in the exotic species imported for aquaculture or rather was accidentally introduced with occasional imports of fingerlings from the Mediterranean Sea could not be determined.     

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.     

Uptake of Carbon Monoxide and Hydrogen at Environmentally Relevant Concentrations by Mycobacteria†

Liquid culture assays revealed a previously unreported capacity for Mycobacterium bovis BCG, M. gordonae, and M. marinum to oxidize CO and for M. smegmatis to consume molecular hydrogen. M. bovis BCG, M. gordonae, M. smegmatis, and M. tuberculosis H37Ra oxidized CO at environmentally relevant concentrations (<50 ppm); H₂ oxidation by M. gordonae and M. smegmatis also occurred at environmentally relevant concentrations (<10 ppm). CO was not consumed by M. avium or M. microti, although the latter appeared to possess CO dehydrogenase (CODH) genes based on PCR results with primers designed for the CODH large subunit, coxL. M. smegmatis and M. gordonae oxidized CO under suboxic (10 and 1% atmospheric oxygen) and anoxic conditions in the presence of nitrate; no oxidation occurred under anoxic conditions without nitrate. Similar results were obtained for H₂ oxidation by M. smegmatis. Phylogenetic analyses of coxL PCR products indicated that mycobacterial sequences form a subclade distinct from that of other bacterial coxL, with limited differentiation among fast- and slow-growing strains.     

Genotypic characteristics of a Mycobacterium sp. isolated from yellowtail Seriola quinqueradiata and striped jack Pseudocaranx dentex in Japan

In Japan, a Mycobacterium marinum‐like mycobacterium was isolated from the yellowtail, Seriola quinqueradiata. The species was identified as M. marinum by a commercial mycobacterial DNA‐DNA hybridization kit. Nevertheless, PCR restriction analysis of the DNA of its RNA polymerase β‐subunit gene definitively showed that this Mycobacterium sp. was M. ulcerans. PCR analysis revealed the genotypic characteristics of M. ulcerans in the Mycobacterium sp., only the mup053 gene sequence being absent, as has been found previously in other piscine mycobacteria such as M. marinum strains DL240490 and DL045 and M. pseudoshottsii. With one exception, this Mycobacterium sp. and M. pseudoshottsii had identical 16S rRNA gene sequences, which is also probably true of M. marinum strains DL240490 and DL045. Similarly, according to comparisons of the 16S rRNA gene, ITS region, and hsp65 gene sequences, this Mycobacterium sp. is more closely related to M. pseudoshottsii than to M. ulcerans or M. marinum. A PCR product of approximately 2000 bp was amplified from region of difference 9 in the Mycobacterium sp. The nucleotide sequence revealed insertion of IS2404, the sequence of which is 1366 bp long. The novel single nucleotide polymorphisms identified in this region distinguished this Mycobacterium sp. from M. marinum strain DL240490 and M. pseudoshottsii. The present findings raise the possibility that these species have a common ancestor. Further studies are required to improve our understanding of the relationship between their geographical origin and genetic diversity.     

Comparative Genomics of Field Isolates of Mycobacterium bovis and M. caprae Provides Evidence for Possible Correlates with Bacterial Viability and Virulence

Mycobacteria of the Mycobacterium tuberculosis complex (MTBC) greatly affect humans and animals worldwide. The life cycle of mycobacteria is complex and the mechanisms resulting in pathogen infection and survival in host cells are not fully understood. Recently, comparative genomics analyses have provided new insights into the evolution and adaptation of the MTBC to survive inside the host. However, most of this information has been obtained using M. tuberculosis but not other members of the MTBC such as M. bovis and M. caprae. In this study, the genome of three M. bovis (MB1, MB3, MB4) and one M. caprae (MB2) field isolates with different lesion score, prevalence and host distribution phenotypes were sequenced. Genome sequence information was used for whole-genome and protein-targeted comparative genomics analysis with the aim of finding correlates with phenotypic variation with potential implications for tuberculosis (TB) disease risk assessment and control. At the whole-genome level the results of the first comparative genomics study of field isolates of M. bovis including M. caprae showed that as previously reported for M. tuberculosis, sequential chromosomal nucleotide substitutions were the main driver of the M. bovis genome evolution. The phylogenetic analysis provided a strong support for the M. bovis/M. caprae clade, but supported M. caprae as a separate species. The comparison of the MB1 and MB4 isolates revealed differences in genome sequence, including gene families that are important for bacterial infection and transmission, thus highlighting differences with functional implications between isolates otherwise classified with the same spoligotype. Strategic protein-targeted analysis using the ESX or type VII secretion system, proteins linking stress response with lipid metabolism, host T cell epitopes of mycobacteria, antigens and peptidoglycan assembly protein identified new genetic markers and candidate vaccine antigens that warrant further study to develop tools to evaluate risks for TB disease caused by M. bovis/M.caprae and for TB control in humans and animals.     

MKAN27435 Is Required for the Biosynthesis of Higher Subclasses of Lipooligosaccharides in Mycobacterium kansasii

Lipooligosaccharides are glycolipids found in the cell wall of many mycobacterial species including the opportunistic pathogen Mycobacterium kansasii. The genome of M. kansasii ATCC12478 contains a cluster with genes orthologous to Mycobacterium marinum LOS biosynthesis genes. To initiate a genetic dissection of this cluster and demonstrate its role in LOS biosynthesis in M. kansasii, we chose MKAN27435, a gene encoding a putative glycosyltransferase. Using Specialized Transduction, a phage-based gene knockout tool previously used to generate null mutants in other mycobacteria, we generated a MKAN27435 null mutant. The mutant strain was found to be defective in the biosynthesis of higher LOS subspecies, viz LOS-IV, LOS-V, LOS-VI and LOS-VII. Additionally, a range of low abundance species were detected in the mutant strain and mass spectroscopic analysis indicated that these were shunt products generated from LOS-III by the addition of up to six molecules of a pentose.     

Discovery and evaluation of candidate sex-determining genes and xenobiotics in the gonads of lake sturgeon (Acipenser fulvescens)

Modern pyrosequencing has the potential to uncover many interesting aspects of genome evolution, even in lineages where genomic resources are scarce. In particular, 454 pyrosequencing of nonmodel species has been used to characterize expressed sequence tags, xenobiotics, gene ontologies, and relative levels of gene expression. Herein, we use pyrosequencing to study the evolution of genes expressed in the gonads of a polyploid fish, the lake sturgeon (Acipenser fulvescens). Using 454 pyrosequencing of transcribed genes, we produced more than 125 MB of sequence data from 473,577 high-quality sequencing reads. Sequences that passed stringent quality control thresholds were assembled into 12,791 male contigs and 32,629 female contigs. Average depth of coverage was 4.2 × for the male assembly and 5.5× for the female assembly. Analytical rarefaction indicates that our assemblies include most of the genes expressed in lake sturgeon gonads. Over 86,700 sequencing reads were assigned gene ontologies, many to general housekeeping genes like protein, RNA, and ion binding genes. We searched specifically for sex determining genes and documented significant sex differences in the expression of two genes involved in animal sex determination, DMRT1 and TRA-1. DMRT1 is the master sex determining gene in birds and in medaka (Oryzias latipes) whereas TRA-1 helps direct sexual differentiation in nematodes. We also searched the lake sturgeon assembly for evidence of xenobiotic organisms that may exist as endosymbionts. Our results suggest that exogenous parasites (trematodes) and pathogens (protozoans) apparently have infected lake sturgeon gonads, and the trematodes have horizontally transferred some genes to the lake sturgeon genome.     

Role of the horizontal gene exchange in evolution of pathogenic Mycobacteria

Background

Mycobacterium tuberculosis is one of the most dangerous human pathogens, the causative agent of tuberculosis. While this pathogen is considered as extremely clonal and resistant to horizontal gene exchange, there are many facts supporting the hypothesis that on the early stages of evolution the development of pathogenicity of ancestral Mtb has started with a horizontal acquisition of virulence factors. Episodes of infections caused by non-tuberculosis Mycobacteria reported worldwide may suggest a potential for new pathogens to appear. If so, what is the role of horizontal gene transfer in this process?

Results

Availing of accessibility of complete genomes sequences of multiple pathogenic, conditionally pathogenic and saprophytic Mycobacteria, a genome comparative study was performed to investigate the distribution of genomic islands among bacteria and identify ontological links between these mobile elements. It was shown that the ancient genomic islands from M. tuberculosis still may be rooted to the pool of mobile genetic vectors distributed among Mycobacteria. A frequent exchange of genes was observed between M. marinum and several saprophytic and conditionally pathogenic species. Among them M. avium was the most promiscuous species acquiring genetic materials from diverse origins.

Conclusions

Recent activation of genetic vectors circulating among Mycobacteria potentially may lead to emergence of new pathogens from environmental and conditionally pathogenic Mycobacteria. The species which require monitoring are M. marinum and M. avium as they eagerly acquire genes from different sources and may become donors of virulence gene cassettes to other micro-organisms.

     

Comparative Genome Analysis Provides Insights into the Pathogenicity of Flavobacterium psychrophilum

Flavobacterium psychrophilum is a fish pathogen in salmonid aquaculture worldwide that causes cold water disease (CWD) and rainbow trout fry syndrome (RTFS). Comparative genome analyses of 11 F. psychrophilum isolates representing temporally and geographically distant populations were used to describe the F. psychrophilum pan-genome and to examine virulence factors, prophages, CRISPR arrays, and genomic islands present in the genomes. Analysis of the genomic DNA sequences were complemented with selected phenotypic characteristics of the strains. The pan genome analysis showed that F. psychrophilum could hold at least 3373 genes, while the core genome contained 1743 genes. On average, 67 new genes were detected for every new genome added to the analysis, indicating that F. psychrophilum possesses an open pan genome. The putative virulence factors were equally distributed among isolates, independent of geographic location, year of isolation and source of isolates. Only one prophage-related sequence was found which corresponded to the previously described prophage 6H, and appeared in 5 out of 11 isolates. CRISPR array analysis revealed two different loci with dissimilar spacer content, which only matched one sequence in the database, the temperate bacteriophage 6H. Genomic Islands (GIs) were identified in F. psychrophilum isolates 950106-1/1 and CSF 259–93, associated with toxins and antibiotic resistance. Finally, phenotypic characterization revealed a high degree of similarity among the strains with respect to biofilm formation and secretion of extracellular enzymes. Global scale dispersion of virulence factors in the genomes and the abilities for biofilm formation, hemolytic activity and secretion of extracellular enzymes among the strains suggested that F. psychrophilum isolates have a similar mode of action on adhesion, colonization and destruction of fish tissues across large spatial and temporal scales of occurrence. Overall, the genomic characterization and phenotypic properties may provide new insights to the mechanisms of pathogenicity in F. psychrophilum.     

Hox Genes in the Parasitic Platyhelminthes Mesocestoides corti, Echinococcus multilocularis, and Schistosoma mansoni: Evidence for a Reduced Hox Complement

Little is known about the Hox gene complement in parasitic platyhelminthes (Neodermata). With the aim of identifying Hox genes in this group we performed two independent strategies: we performed a PCR survey with degenerate primers directed to the Hox homeobox in the cestode Mesocestoides corti, and we searched genomic assemblies of Echinococcus multilocularis and Schistosoma mansoni. We identified two Hox genes in M. corti, seven in E. multilocularis, and nine in S. mansoni (including five previously reported). The affinities of these sequences, and other previously reported Hox sequences from flatworms, were determined according to phylogenetic analysis, presence of characteristic parapeptide sequences, and unusual intron positions. Our results suggest that the last common ancestor of triclads and neodermatans had a Hox gene complement of at least seven genes, and that this was probably derived by gene loss from a larger ancestral Hox complement in lophotrochozoans.     

Rapid Detection and Identification of Nontuberculous Mycobacterial Pathogens in Fish by Using High-Resolution Melting Analysis

Mycobacterial infections in fish are commonly referred to as piscine mycobacteriosis, irrespectively of the specific identity of the causal organism. They usually cause a chronic disease and sometimes may result in high mortalities and severe economic losses. Nearly 20 species of Mycobacterium have been reported to infect fish. Among them, Mycobacterium marinum, M. fortuitum, and M. chelonae are generally considered the major agents responsible for fish mycobacteriosis. As no quick and inexpensive diagnostic test exists, we tested the potential of high-resolution melting analysis (HRMA) to rapidly identify and differentiate several Mycobacterium species involved in fish infections. By analyzing both the melting temperature and melting profile of the 16S-23S rRNA internal transcribed spacer (ITS), we were able to discriminate 12 different species simultaneously. Sensitivity tests conducted on purified M. marinum and M. fortuitum DNA revealed a limit of detection of 10 genome equivalents per reaction. The primers used in this procedure did not lead to any amplification signal with 16 control non-Mycobacterium species, thereby demonstrating their specificity for the genus Mycobacterium.     

Mycobacterium ulcerans Persistence at a Village Water Source of Buruli Ulcer Patients

Buruli ulcer (BU), a neglected tropical disease of the skin and subcutaneous tissue, is caused by Mycobacterium ulcerans and is the third most common mycobacterial disease after tuberculosis and leprosy. While there is a strong association of the occurrence of the disease with stagnant or slow flowing water bodies, the exact mode of transmission of BU is not clear. M. ulcerans has emerged from the environmental fish pathogen M. marinum by acquisition of a virulence plasmid encoding the enzymes required for the production of the cytotoxic macrolide toxin mycolactone, which is a key factor in the pathogenesis of BU. Comparative genomic studies have further shown extensive pseudogene formation and downsizing of the M. ulcerans genome, indicative for an adaptation to a more stable ecological niche. This has raised the question whether this pathogen is still present in water-associated environmental reservoirs. Here we show persistence of M. ulcerans specific DNA sequences over a period of more than two years at a water contact location of BU patients in an endemic village of Cameroon. At defined positions in a shallow water hole used by the villagers for washing and bathing, detritus remained consistently positive for M. ulcerans DNA. The observed mean real-time PCR Ct difference of 1.45 between the insertion sequences IS2606 and IS2404 indicated that lineage 3 M. ulcerans, which cause human disease, persisted in this environment after successful treatment of all local patients. Underwater decaying organic matter may therefore represent a reservoir of M. ulcerans for direct infection of skin lesions or vector-associated transmission.     

Oxidative Stress Response and Characterization of the oxyR-ahpC and furA-katG Loci in Mycobacterium marinum

Oxidative stress response in pathogenic mycobacteria is believed to be of significance for host-pathogen interactions at various stages of infection. It also plays a role in determining the intrinsic susceptibility to isoniazid in mycobacterial species. In this work, we characterized the oxyR-ahpC and furA-katG loci in the nontuberculous pathogen Mycobacterium marinum. In contrast to Mycobacterium smegmatis and like Mycobacterium tuberculosis and Mycobacterium leprae, M. marinum was shown to possess a closely linked and divergently oriented equivalents of the regulator of peroxide stress response oxyR and its subordinate gene ahpC, encoding a homolog of alkyl hydroperoxide reductase. Purified mycobacterial OxyR was found to bind to the oxyR-ahpC promoter region from M. marinum and additional mycobacterial species. Mobility shift DNA binding analyses using OxyR binding sites from several mycobacteria and a panel of in vitro-generated mutants validated the proposed consensus mycobacterial recognition sequence. M. marinum AhpC levels detected by immunoblotting, were increased upon treatment with H₂O₂, in keeping with the presence of a functional OxyR and its binding site within the promoter region of ahpC. In contrast, OxyR did not bind to the sequences upstream of the katG structural gene, and katG expression did not follow the pattern seen with ahpC. Instead, a new open reading frame encoding a homolog of the ferric uptake regulator Fur was identified immediately upstream of katG in M. marinum. The furA-katG linkage and arrangement are ubiquitous in mycobacteria, suggesting the presence of additional regulators of oxidative stress response and potentially explaining the observed differences in ahpC and katG expression. Collectively, these findings broaden our understanding of oxidative stress response in mycobacteria. They also suggest that M. marinum will be useful as a model system for studying the role of oxidative stress response in mycobacterial physiology, intracellular survival, and other host-pathogen interactions associated with mycobacterial diseases.     

Whole-Genome Sequencing and Comparative Analysis of Mycobacterium brisbanense Reveals a Possible Soil Origin and Capability in Fertiliser Synthesis

Mycobacterium brisbanense is a member of Mycobacterium fortuitum third biovariant complex, which includes rapidly growing Mycobacterium spp. that normally inhabit soil, dust and water, and can sometimes cause respiratory tract infections in humans. We present the first whole-genome analysis of M. brisbanense UM_WWY which was isolated from a 70-year-old Malaysian patient. Molecular phylogenetic analyses confirmed the identification of this strain as M. brisbanense and showed that it has an unusually large genome compared with related mycobacteria. The large genome size of M. brisbanense UM_WWY (~7.7Mbp) is consistent with further findings that this strain has a highly variable genome structure that contains many putative horizontally transferred genomic islands and prophage. Comparative analysis showed that M. brisbanense UM_WWY is the only Mycobacterium species that possesses a complete set of genes encoding enzymes involved in the urea cycle, suggesting that this soil bacterium is able to synthesize urea for use as plant fertilizers. It is likely that M. brisbanense UM_WWY is adapted to live in soil as its primary habitat since the genome contains many genes associated with nitrogen metabolism. Nevertheless, a large number of predicted virulence genes were identified in M. brisbanense UM_WWY that are mostly shared with well-studied mycobacterial pathogens such as Mycobacterium tuberculosis and Mycobacterium abscessus. These findings are consistent with the role of M. brisbanense as an opportunistic pathogen of humans. The whole-genome study of UM_WWY has provided the basis for future work of M. brisbanense.     

Nuclear DNA content variation of three Miscanthus species in China

In order to estimate the variation in nuclear genome size in Miscanthus, flow cytometry of nuclei stained by propidium iodide was carried out using 36 populations of three Miscanthus species: M. lutarioriparius, M. sacchariflorus and M. sinensis, which were sampled from cold northern to warm and humid southern and central China, as well as near the sea level in eastern China to mountains in western China. The DNA content of diploid was 4.37 ± 0.02 pg/2C in M. lutarioriparius, 4.37 ± 0.01 pg/2C in M. sacchariflorus, and 5.37 ± 0.03 pg/2C in M. sinensis, respectively. There was no intraspecific variation in the three Miscanthus species at the diploid level, suggesting that the genome size was stable within species and the diverse environments did not induce variation in genome size at the diploid level. However, tetraploid populations were found in M. lutarioriparius and M. sacchariflorus, and their genome sizes were 8.56 and 8.54 pg, respectively, which are lower than expected values (8.74 pg), indicating the genome downsizing after polyploidization in the genus. Our results showed that the plant height of M. lutarioriparius was the highest one among the three species and the species was more closely related to M. sacchariflorus than M. sinensis. The intra-species genomic variation and inter-species differentiation in Miscanthus species provide important genetic and genomic information for the development of Miscanthus, especially for the endemic species, M. lutarioriparius, (together with Miscanthus × giganteus) which are now emerging as a key bio-energy crop because of their high yields and strong adaptability.