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.     

Comparative Genome Analysis of Fish and Human Isolates of Mycobacterium marinum

Mycobacterium marinum is a major causative agent of mycobacteriosis in fish that has a broad range of hosts, including in human isolates. So far, genomic analyses have focused on the human isolate. Here, we compared the draft genome sequences of two strains of M. marinum isolated from fish (MB2 and Europe) with the M. marinum M isolated from humans. M. marinum MB2 and Europe have single, circular chromosomes of 6,134,389 and 6,029,340 bp, and average G + C contents of 65.7 and 65.5 %, respectively. A total of 5,464 coding DNA sequences were annotated in both M. marinum MB2 and Europe genome. Dot plot analyses showed that M. marinum MB2 and Europe were closer to M. marinum M when compared to three other Mycobacterium species. The insertion/deletion gene analysis showed that M. marinum MB2 and Europe contained 342 and 487 genes that were not found in M. marinum M, and lacked 625 and 776 genes found in M. marinum M, respectively. Most of the inserted and deleted genes were classified in the fatty acid, lipid, and isoprenoid subsystem and the virulence, disease, and defense subsystem. Therefore, these results provide insights into the genomic diversity associated with variable hosts and pathogens.     

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.     

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.     

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.     

Mycobacterium marinum produces distinct mycobactin and carboxymycobactin siderophores to promote growth in broth and phagocytes

Mycobacterium marinum is a model organism for pathogenic Mycobacterium species, including Mycobacterium tuberculosis, the causative agent of tuberculosis. These pathogens enter phagocytes and replicate within the Mycobacterium‐containing vacuole, possibly followed by vacuole exit and growth in the host cell cytosol. Mycobacteria release siderophores called mycobactins to scavenge iron, an essential yet poorly soluble and available micronutrient. To investigate the role of M. marinum mycobactins, we purified by organic solvent extraction and identified by mass spectrometry the lipid‐bound mycobactin (MBT) and the water‐soluble variant carboxymycobactin (cMBT). Moreover, we generated by specialised phage transduction a defined M. marinum ΔmbtB deletion mutant predicted to be defective for mycobactin production. The M. marinum ΔmbtB mutant strain showed a severe growth defect in broth and phagocytes, which was partially complemented by supplying the mbtB gene on a plasmid. Furthermore, purified Fe‐MBT or Fe‐cMBT improved the growth of wild type as well as ΔmbtB mutant bacteria on minimal plates, but only Fe‐cMBT promoted the growth of wild‐type M. marinum during phagocyte infection. Finally, the intracellular growth of M. marinum ΔmbtB in Acanthamoeba castellanii amoebae was restored by coinfection with wild‐type bacteria. Our study identifies and characterises the M. marinum MBT and cMBT siderophores and reveals the requirement of mycobactins for extra‐ and intracellular growth of the pathogen.     

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.

     

Systematic Phenotyping of a Large-Scale Candida glabrata Deletion Collection Reveals Novel Antifungal Tolerance Genes

The opportunistic fungal pathogen Candida glabrata is a frequent cause of candidiasis, causing infections ranging from superficial to life-threatening disseminated disease. The inherent tolerance of C. glabrata to azole drugs makes this pathogen a serious clinical threat. To identify novel genes implicated in antifungal drug tolerance, we have constructed a large-scale C. glabrata deletion library consisting of 619 unique, individually bar-coded mutant strains, each lacking one specific gene, all together representing almost 12% of the genome. Functional analysis of this library in a series of phenotypic and fitness assays identified numerous genes required for growth of C. glabrata under normal or specific stress conditions, as well as a number of novel genes involved in tolerance to clinically important antifungal drugs such as azoles and echinocandins. We identified 38 deletion strains displaying strongly increased susceptibility to caspofungin, 28 of which encoding proteins that have not previously been linked to echinocandin tolerance. Our results demonstrate the potential of the C. glabrata mutant collection as a valuable resource in functional genomics studies of this important fungal pathogen of humans, and to facilitate the identification of putative novel antifungal drug target and virulence genes.     

Comparative genomic analysis of pyrene-degrading <Emphasis Type="Italic">Mycobacterium</Emphasis> species: Genomic islands and ring-hydroxylating dioxygenases involved in pyrene degradation

The genome sequences of two pyrene-degrading bacterial strains of Mycobacterium spp. PYR10 and PYR15, isolated from the estuarine wetland of the Han river, South Korea, were determined using the PacBio RS II sequencing platform. The complete genome of strain PYR15 was 6,037,017 bp in length with a GC content of 66.5%, and contained 5,933 protein-coding genes. The genome of strain PYR10 was 5,999,427 bp in length with a GC content of 67.7%, and contained 5,767 protein-coding genes. Based on the average nucleotide identity values, these strains were designated as M. gilvum PYR10 and M. pallens PYR15. A genomic comparison of these pyrene-degrading Mycobacterium strains with pyrene-non-degrading strains revealed that the genomes of pyrene-degrading strains possessed similar repertoires of ringhydroxylating dioxygenases (RHDs), including the pyrenehydroxylating dioxygenases encoded by nidA and nidA3, which could be readily distinguished from those of pyrenenon-degraders. Furthermore, genomic islands, containing catabolic gene clusters, were shared only among the pyrenedegrading Mycobacterium strains and these gene clusters contained RHD genes, including nidAB and nidA3B3. Our genome data should facilitate further studies on the evolution of the polycyclic aromatic hydrocarbon-degradation pathways in the genus Mycobacterium.     

A Novel ESX-1 Locus Reveals that Surface-Associated ESX-1 Substrates Mediate Virulence in Mycobacterium marinum

EsxA (ESAT-6) and EsxB (CFP-10) are virulence factors exported by the ESX-1 system in mycobacterial pathogens. In Mycobacterium marinum, an established model for ESX-1 secretion in Mycobacterium tuberculosis, genes required for ESX-1 export reside at the extended region of difference 1 (RD1) locus. In this study, a novel locus required for ESX-1 export in M. marinum was identified outside the RD1 locus. An M. marinum strain bearing a transposon-insertion between the MMAR_1663 and MMAR_1664 genes exhibited smooth-colony morphology, was deficient for ESX-1 export, was nonhemolytic, and was attenuated for virulence. Genetic complementation revealed a restoration of colony morphology and a partial restoration of virulence in cell culture models. Yet hemolysis and the export of ESX-1 substrates into the bacteriological medium in vitro as measured by both immunoblotting and quantitative proteomics were not restored. We show that genetic complementation of the transposon insertion strain partially restored the translocation of EsxA and EsxB to the mycobacterial cell surface. Our findings indicate that the export of EsxA and EsxB to the cell surface, rather than secretion into the bacteriological medium, correlates with virulence in M. marinum. Together, these findings not only expand the known genetic loci required for ESX-1 secretion in M. marinum but also provide an explanation for the observed disparity between in vitro ESX-1 export and virulence.     

Insights from the comparative genome analysis of natural rubber degrading Nocardia species

Nocardia are known to be a facultative human pathogen and can cause infection in immune compromised patients. Though the details research on the virulence factors of Nocardia are scanty but numerous genes that code such factors were reported from different species of Nocardia. Despite of the presence of several virulence factors, species of this genus have been shown to have role in remediation of many toxic and hazardous materials from the environment. In this study, genome sequences of rubber degrading Nocardia sp. BSTN01 and N.nova SH22a have been analyzed to locate the potential virulence genes. Also, the genomes of facultative pathogenic Nocardia like, N.africana, N. brasiliensis, N. kruczakiae, N. transvalensis and N. veterana have been analyzed to find the gene encoding latex clearing protein (Lcp), a rubber oxygenase enzyme of Gram-positive action bacteria. The study provides an insight about the potentiality of rubberdegrading Nocardia species to emerge as future human pathogens and also the probability of a serious concern if the studied facultative pathogens of Nocardia like N. africana, N. brasiliensis, N. kruczakiae, N. transvalensis and N. veterana are capable of degrading rubber, a regularly used material in clinics. Moreover, use of such possible pathogenic strains for their known role in bioremediation of rubber waste from the environment might be deleterious.     

Spectrum of heart diseases in a new cardiac service in Nigeria: An echocardiographic study of 1441 subjects in Abeokuta

Background

The recent determination of the complete nucleotide sequence of several Mycobacterium tuberculosis (MTB) genomes allows the use of comparative genomics as a tool for dissecting the nature and consequence of genetic variability within this species. The multiple alignment of the genomes of clinical strains (CDC1551, F11, Haarlem and C), along with the genomes of laboratory strains (H37Rv and H37Ra), provides new insights on the mechanisms of adaptation of this bacterium to the human host.

Findings

The genetic variation found in six M. tuberculosis strains does not involve significant genomic rearrangements. Most of the variation results from deletion and transposition events preferentially associated with insertion sequences and genes of the PE/PPE family but not with genes implicated in virulence. Using a Perl-based software islandsanalyser, which creates a representation of the genetic variation in the genome, we identified differences in the patterns of distribution and frequency of the polymorphisms across the genome. The identification of genes displaying strain-specific polymorphisms and the extrapolation of the number of strain-specific polymorphisms to an unlimited number of genomes indicates that the different strains contain a limited number of unique polymorphisms.

Conclusion

The comparison of multiple genomes demonstrates that the M. tuberculosis genome is currently undergoing an active process of gene decay, analogous to the adaptation process of obligate bacterial symbionts. This observation opens new perspectives into the evolution and the understanding of the pathogenesis of this bacterium.     

Vaccine efficacy of Mycobacterium bovis BCG against Mycobacterium sp. infection in amberjack Seriola dumerili

Mycobacteriosis, caused by the intracellular parasitism Mycobacterium sp., causes economic damages to aquaculture production in Japan, particularly in seriola fish production. Antibiotics are not effective against Mycobacterium sp. and so a potent vaccine is needed. We previously reported that BCG vaccine (Mycobacterium bovis BCG) induces adaptive immunity against Mycobacterium sp. in Japanese flounder, Paralichthys olivaceus. In a phylogenetic tree, the genes for a major antigen, the Ag85 complex, in Mycobacterium sp. TUMSAT-Msp001 are closely related to homologues in Mycobacterium ulcerans. M. bovis BCG was detected until 7 days post-injection at the injection site (muscle) and 28 days post-vaccination in spleen. Cumulative mortality of amberjack, Seriola dumerili vaccinated intramuscularly (i.m.) and intraperitoneally (i.p.) with M. bovis BCG was 32.3% and 59.5% respectively, at 24 days post-infection of Mycobacterium sp., compared to 97.8% in PBS-injected fish. The bacterial counts of Mycobacterium sp. in spleen of both i.m.-and i.p.-vaccinated fish (6.2 × 10³ and 1.3 × 10⁴ CFU/mg tissue, respectively) at 20 days post-infection were significantly lower (P < 0.01) than those of PBS-injected fish (8.0 × 10⁶ CFU/mg). Furthermore, Immersion challenge with Mycobacterium sp. TUMSAT Msp-001 showed 50% RPS value in BCG i.m.-vaccinated fish at the end of the experiment. These results support our previous study using Japanese flounder and suggest that BCG vaccine is also effective against Mycobacterium sp. infection in amberjack.     

Non Mycobacterial Virulence Genes in the Genome of the Emerging Pathogen Mycobacterium abscessus

Mycobacterium abscessus is an emerging rapidly growing mycobacterium (RGM) causing a pseudotuberculous lung disease to which patients with cystic fibrosis (CF) are particularly susceptible. We report here its complete genome sequence. The genome of M. abscessus (CIP 104536T) consists of a 5,067,172-bp circular chromosome including 4920 predicted coding sequences (CDS), an 81-kb full-length prophage and 5 IS elements, and a 23-kb mercury resistance plasmid almost identical to pMM23 from Mycobacterium marinum. The chromosome encodes many virulence proteins and virulence protein families absent or present in only small numbers in the model RGM species Mycobacterium smegmatis. Many of these proteins are encoded by genes belonging to a “mycobacterial” gene pool (e.g. PE and PPE proteins, MCE and YrbE proteins, lipoprotein LpqH precursors). However, many others (e.g. phospholipase C, MgtC, MsrA, ABC Fe(3+) transporter) appear to have been horizontally acquired from distantly related environmental bacteria with a high G+C content, mostly actinobacteria (e.g. Rhodococcus sp., Streptomyces sp.) and pseudomonads. We also identified several metabolic regions acquired from actinobacteria and pseudomonads (relating to phenazine biosynthesis, homogentisate catabolism, phenylacetic acid degradation, DNA degradation) not present in the M. smegmatis genome. Many of the “non mycobacterial” factors detected in M. abscessus are also present in two of the pathogens most frequently isolated from CF patients, Pseudomonas aeruginosa and Burkholderia cepacia. This study elucidates the genetic basis of the unique pathogenicity of M. abscessus among RGM, and raises the question of similar mechanisms of pathogenicity shared by unrelated organisms in CF patients.     

Severe inhibition of lipooligosaccharide synthesis induces TLR2-dependent elimination of <Emphasis Type="Italic">Mycobacterium marinum</Emphasis> from THP1-derived macrophages

Background

Although mycobacterial glycolipids are among the first-line molecules involved in host–pathogen interactions, their contribution in virulence remains incomplete. Mycobacterium marinum is a waterborne pathogen of fish and other ectotherms, closely related to Mycobacterium tuberculosis. Since it causes tuberculosis-like systemic infection it is widely used as a model organism for studying the pathogenesis of tuberculosis. It is also an occasional opportunistic human pathogen. The M. marinum surface-exposed lipooligosaccharides (LOS) are immunogenic molecules that participate in the early interactions with macrophages and modulate the host immune system. Four major LOS species, designated LOS-I to LOS-IV, have been identified and characterized in M. marinum. Herein, we investigated the interactions between a panel of defined M. marinum LOS mutants that exhibited various degrees of truncation in the LOS structure, and human-derived THP-1 macrophages to address the potential of LOSs to act as pro- or avirulence factors.

Results

A moderately truncated LOS structure did not interfere with M. marinum invasion. However, a deeper shortening of the LOS structure was associated with increased entry of M. marinum into host cells and increased elimination of the bacilli by the macrophages. These effects were dependent on Toll-like receptor 2.

Conclusion

We provide the first evidence that LOSs inhibit the interaction between mycobacterial cell wall ligands and appropriate macrophage pattern recognition receptors, affecting uptake and elimination of the bacteria by host phagocytes.

     

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.     

Identification and Characterization of Putative Virulence Genes and Gene Clusters in Aeromonas hydrophila PPD134/91

Aeromonas hydrophila is a gram-negative opportunistic pathogen of animals and humans. The pathogenesis of A. hydrophila is multifactorial. Genomic subtraction and markers of genomic islands (GIs) were used to identify putative virulence genes in A. hydrophila PPD134/91. Two rounds of genomic subtraction led to the identification of 22 unique DNA fragments encoding 19 putative virulence factors and seven new open reading frames, which are commonly present in the eight virulence strains examined. In addition, four GIs were found, including O-antigen, capsule, phage-associated, and type III secretion system (TTSS) gene clusters. These putative virulence genes and gene clusters were positioned on a physical map of A. hydrophila PPD134/91 to determine their genetic organization in this bacterium. Further in vivo study of insertion and deletion mutants showed that the TTSS may be one of the important virulence factors in A. hydrophila pathogenesis. Furthermore, deletions of multiple virulence factors such as S-layer, serine protease, and metalloprotease also increased the 50% lethal dose to the same level as the TTSS mutation (about 1 log) in a blue gourami infection model. This observation sheds light on the multifactorial and concerted nature of pathogenicity in A. hydrophila. The large number of putative virulence genes identified in this study will form the basis for further investigation of this emerging pathogen and help to develop effective vaccines, diagnostics, and novel therapeutics.     

Genomic and proteomic analyses of the coral pathogen Vibrio coralliilyticus reveal a diverse virulence repertoire

Vibrio coralliilyticus has been implicated as an important pathogen of coral species worldwide. In this study, the nearly complete genome of Vibrio coralliilyticus strain P1 (LMG23696) was sequenced and proteases implicated in virulence of the strain were specifically investigated. The genome sequence of P1 (5 513 256 bp in size) consisted of 5222 coding sequences and 58 RNA genes (53 tRNAs and at least 5 rRNAs). Seventeen metalloprotease and effector (vgrG, hlyA and hcp) genes were identified in the genome and expressed proteases were also detected in the secretome of P1. As the VcpA zinc-metalloprotease has been considered an important virulence factor of V. coralliilyticus, a vcpA deletion mutant was constructed to evaluate the effect of this gene in animal pathogenesis. Both wild-type and mutant (ΔvcpA) strains exhibited similar virulence characteristics that resulted in high mortality in Artemia and Drosophila pathogenicity bioassays and strong photosystem II inactivation of the coral dinoflagellate endosymbiont (Symbiodinium). In contrast, the ΔvcpA mutant demonstrated higher hemolytic activity and secreted 18 proteins not secreted by the wild type. These proteins included four types of metalloproteases, a chitinase, a hemolysin-related protein RbmC, the Hcp protein and 12 hypothetical proteins. Overall, the results of this study indicate that V. coralliilyticus strain P1 has a diverse virulence repertoire that possibly enables this bacterium to be an efficient animal pathogen.     

Whole-genome analysis of Fusarium graminearum insertional mutants identifies virulence associated genes and unmasks untagged chromosomal deletions

Background

Identifying pathogen virulence genes required to cause disease is crucial to understand the mechanisms underlying the pathogenic process. Plasmid insertion mutagenesis of fungal protoplasts is frequently used for this purpose in filamentous ascomycetes. Post transformation, the mutant population is screened for loss of virulence to a specific plant or animal host. Identifying the insertion event has previously met with varying degrees of success, from a cleanly disrupted gene with minimal deletion of nucleotides at the insertion point to multiple-copy insertion events and large deletions of chromosomal regions. Currently, extensive mutant collections exist in laboratories globally where it was hitherto impossible to identify all the affected genes.

Results

We used a whole-genome sequencing (WGS) approach using Illumina HiSeq 2000 technology to investigate DNA tag insertion points and chromosomal deletion events in mutagenised, reduced virulence F. graminearum isolates identified in disease tests on wheat (Triticum aestivum). We developed the FindInsertSeq workflow to localise the DNA tag insertions to the nucleotide level. The workflow was tested using four mutants showing evidence of single and multi-copy insertions in DNA blot analysis. FindInsertSeq was able to identify both single and multi-copy concatenation insertion sites. By comparing sequencing coverage, unexpected molecular recombination events such as large tagged and untagged chromosomal deletions, and DNA amplification were observed in three of the analysed mutants. A random data sampling approach revealed the minimum genome coverage required to survey the F. graminearum genome for alterations.

Conclusions

This study demonstrates that whole-genome re-sequencing to 22x fold genome coverage is an efficient tool to characterise single and multi-copy insertion mutants in the filamentous ascomycete Fusarium graminearum. In some cases insertion events are accompanied with large untagged chromosomal deletions while in other cases a straight-forward insertion event could be confirmed. The FindInsertSeq analysis workflow presented in this study enables researchers to efficiently characterise insertion and deletion mutants.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1412-9) contains supplementary material, which is available to authorized users.