Discovery of stable and variable differences in the Mycobacterium avium subsp. paratuberculosis type I, II, and III genomes by pan-genome microarray analysis

Comparative Genomics of Mycobacterium avium Subspecies Paratuberculosis Sheep Strains

Mycobacterium avium subspecies paratuberculosis (MAP) is the aetiological agent of Johne's disease (JD), a chronic enteritis that causes major losses to the global livestock industry. Further, it has been associated with human Crohn's disease. Several strains of MAP have been identified, the two major groups being sheep strain MAP, which includes the Type I and Type III sub-lineages, and the cattle strain or Type II MAP lineage, of which bison strains are a sub-grouping. Major genotypic, phenotypic and pathogenic variations have been identified in prior comparisons, but the research has predominately focused on cattle strains of MAP. In countries where the sheep industries are more prevalent, however, such as Australia and New Zealand, ovine JD is a substantial burden. An information gap exists regarding the genomic differences between sheep strain sub-lineages and the relevance of Type I and Type III MAP in terms of epidemiology and/or pathogenicity. We therefore investigated sheep MAP isolates from Australia and New Zealand using whole genome sequencing. For additional context, sheep MAP genome datasets were downloaded from the Sequence Read Archive and GenBank. The final dataset contained 18 Type III and 16 Type I isolates and the K10 cattle strain MAP reference genome. Using a pan-genome approach, an updated global phylogeny for sheep MAP from de novo assemblies was produced. When rooted with the K10 cattle reference strain, two distinct clades representing the lineages were apparent. The Australian and New Zealand isolates formed a distinct sub-clade within the type I lineage, while the European type I isolates formed another less closely related group. Within the type III lineage, isolates appeared more genetically diverse and were from a greater number of continents. Querying of the pan-genome and verification using BLAST analysis revealed lineage-specific variations (n = 13) including genes responsible for metabolism and stress responses. The genetic differences identified may represent important epidemiological and virulence traits specific to sheep MAP. This knowledge will potentially contribute to improved vaccine development and control measures for these strains.

Complete Genome Sequence of Ovine Mycobacterium avium subsp. paratuberculosis Strain JIII-386 (MAP-S/type III) and Its Comparison to MAP-S/type I, MAP-C, and M. avium Complex Genomes

Mycobacterium avium (M. a.) subsp. paratuberculosis (MAP) is a worldwide-distributed obligate pathogen in ruminants causing Johne’s disease. Due to a lack of complete subtype III genome sequences, there is not yet conclusive information about genetic differences between strains of cattle (MAP-C, type II) and sheep (MAP-S) type, and especially between MAP-S subtypes I, and III. Here we present the complete, circular genome of MAP-S/type III strain JIII-386 (DE) closed by Nanopore-technology and its comparison with MAP-S/type I closed genome of strain Telford (AUS), MAP-S/type III draft genome of strain S397 (U.S.), twelve closed MAP-C strains, and eight closed M.-a.-complex-strains. Structural comparative alignments revealed clearly the mosaic nature of MAP, emphasized differences between the subtypes and the higher diversity of MAP-S genomes. The comparison of various genomic elements including transposases and genomic islands provide new insights in MAP genomics. MAP type specific phenotypic features may be attributed to genes of known large sequence polymorphisms (LSP^Ss) regions I–IV and deletions #1 and #2, confirmed here, but could also result from identified frameshifts or interruptions of various virulence-associated genes (e.g., mbtC in MAP-S). Comprehensive core and pan genome analysis uncovered unique genes (e.g., cytochromes) and genes probably acquired by horizontal gene transfer in different MAP-types and subtypes, but also emphasized the highly conserved and close relationship, and the complex evolution of M.-a.-strains.

Novel recombinant Mce-truncated protein based ELISA for the diagnosis of Mycobacterium avium subsp. paratuberculosis infection in domestic livestock

Johne’s disease (JD) is an infectious wasting condition of ruminants caused by Mycobacterium avium subsp. paratuberculosis (MAP) in domestic livestock of every country that has been investigated. Controlling JD is problematic due to the lack of sensitive, specific, efficient, and cost-effective diagnostic tests. A major challenge in the development of diagnostics like ELISA is the selection of an ideal antigen/(s) that is pathogen-specific and allows sensitive recognition. Therefore, the purpose of this study was to identify and use Mce-truncated protein-based ELISA assay for the diagnosis of MAP infection with high sensitivity and specificity. In silico epitope prediction by epitope mapping throughout the whole length of MAP2191 protein revealed that C-terminal portion of this protein presented potential T- and B-cell epitopes. Therefore, a novel Mce-truncated protein encoded by the selected region of MAP2191 gene was expressed, purified with Ni-NTA gel matrix and confirmed by SDS PAGE and western blot. A profiling ELISA assay was developed to evaluate sera from MAP infected and non-infected ruminant species for antibodies against Mce-truncated protein to infer the immunogenicity of this protein in the host. Using this Mce protein-based ELISA, 251 goats, 53 sheep, 117 buffaloes, and 33 cattle serum samples were screened and 49.4, 51.0, 69.2, and 54.6% animals, respectively, were found positive. Comparing with i-ELISA, the new Mce-based ELISA kit showed a relatively higher specificity but suffered from slightly reduced sensitivity. Mce-based ELISA excluded apparently false positive results of i-ELISA. Mce protein was found to be antigenic and Mce-ELISA test could be employed as a diagnostic test for JD in domestic livestock in view of the a relatively higher specificity and accuracy. The antigenic potential of Mce antigen can also be exploited for the development of a new vaccine for the control of MAP infection.

Efficient method for targeted gene disruption by homologous recombination in Mycobacterium avium subspecie paratuberculosis

Targeted gene disruption by homologous recombination, has been widely used in mycobacterium species to understand the genetic basis of virulence and persistence in the host and to develop efficacious potential live vaccines. However, in slow growing pathogenic mycobacteria as Mycobacterium avium subsp paratuberculosis (MAP), these methods have been inefficient, in part due to the low frequency of legitimate homologous recombination. Another feature of mycobacteria is the low efficiency of transformation; therefore, some years ago, a phage-mediated transduction process was developed to introduce DNA into mycobacteria. This strategy is very efficient, due to the high rate of infection of the phage. This report describes a genetic method for the generation of targeted deletion mutations in MAP by allelic exchange using in vitro-generated specialized transducing mycobacteriophages, which does not require the critical packaging step and that could also be applied to other mycobacteria. We provide a detailed gene deletion methodology and demonstrate the use of this genetic system by deleting the mce4 operon of MAP. Finally, our results showed that the deletion of mce4 in MAP induces triacylglycerol accumulation; alter morphology and aggregation in liquid culture.

Mammalian cell entry operons; novel and major subset candidates for diagnostics with special reference to Mycobacterium avium subspecies paratuberculosis infection

Mammalian cell entry (mce) genes are the components of the mce operon and play a vital role in the entry of Mycobacteria into the mammalian cell and their survival within phagocytes and epithelial cells. Mce operons are present in the DNA of Mycobacteria and translate proteins associated with the invasion and long-term existence of these pathogens in macrophages. The exact mechanism of action of mce genes and their functions are not clear yet. However, with the loss of these genes Mycobacteria lose their pathogenicity. Mycobacterium avium subspecies paratuberculosis (MAP), the etiological agent of Johne’s disease, is the cause of chronic enteritis of animals and significantly affects economic impact on the livestock industry. Since MAP is not inactivated during pasteurization, human population is continuously at the risk of getting exposed to MAP infection through consumption of dairy products. There is need for new candidate genes and/or proteins for developing improved diagnostic assays for the diagnosis of MAP infection and for the control of disease. Increasing evidences showed that expression of mce genes is important for the virulence of MAP. Whole-genome DNA microarray representing MAP revealed that there are 14 large sequence polymorphisms with LSPP12 being the most widely conserved MAP-specific region that included a cluster of six homologs of mce-family involved in lipid metabolism. On the other hand, LSP11 comprising part of mce2 operon was absent in MAP isolates. This review summarizes the advancement of research on mce genes of Mycobacteria with special reference to the MAP infection.

Evaluation of associations between genotypes of Mycobacterium avium subsp. paratuberculsis and presence of intestinal lesions characteristic of paratuberculosis

Mycobacterium avium subsp. paratuberculosis (MAP) is the causative agent of paratuberculosis affecting ruminants worldwide. Depending on the MAP-Type (MAP-C or MAP-S, cattle or sheep type), strains differ in virulence and host preference. There is not yet any strong evidence indicating that individual field strains of the same MAP-subgroup exhibit differences in virulence. The aim of this study was to evaluate a potential association between the genotype of individual field strains belonging to the MAP-C group and the presence of macroscopic intestinal lesions characteristic of paratuberculosis in the infected animals. 88 MAP-C isolates were sampled from clinically healthy cows at slaughter. Cows were grouped as A (n=46) with, and B (n=42) without macroscopic intestinal lesions. Sampled cows from both the A and B groups came from different farms and had a similar age distribution. MAP isolates were characterized by MIRU-VNTR and IS900-RFLP analysis. Resulting genotypes were examined for an association with the presence of macroscopic intestinal lesions characteristic of paratuberculosis. MAP isolates from groups A and B exhibited similar strain diversity: 20 and 18 combined genotypes, altogether 32 genotypes. Six of these genotypes were detected in both groups. Although no association was found between individual combined genotypes and presence of macroscopic intestinal lesions, IS900-RFLP-(BstEII)-Type-C1 (the most common type worldwide) was found more often in group A (p<0.01). The data give only weak indication for the existence of differences in virulence among MAP-cattle type isolates. Differences in the development and severity of lesions may rather depend on unknown host factors or inoculation dose. Virulence properties of IS900-RFLP-(BstEII)-Type-C1 isolates should be examined in more detail.

Comprehensive insights in the Mycobacterium avium subsp. paratuberculosis genome using new WGS data of sheep strain JIII-386 from Germany

Mycobacterium avium (M. a.) subsp. paratuberculosis (MAP)—the etiologic agent of Johne’s disease—affects cattle, sheep, and other ruminants worldwide. To decipher phenotypic differences among sheep and cattle strains (belonging to MAP-S [Type-I/III], respectively, MAP-C [Type-II]), comparative genome analysis needs data from diverse isolates originating from different geographic regions of the world. This study presents the so far best assembled genome of a MAP-S-strain: Sheep isolate JIII-386 from Germany. One newly sequenced cattle isolate (JII-1961, Germany), four published MAP strains of MAP-C and MAP-S from the United States and Australia, and M. a. subsp. hominissuis (MAH) strain 104 were used for assembly improvement and comparisons. All genomes were annotated by BacProt and results compared with NCBI (National Center for Biotechnology Information) annotation. Corresponding protein-coding sequences (CDSs) were detected, but also CDSs that were exclusively determined by either NCBI or BacProt. A new Shine–Dalgarno sequence motif (5′-AGCTGG-3′) was extracted. Novel CDSs including PE-PGRS family protein genes and about 80 noncoding RNAs exhibiting high sequence conservation are presented. Previously found genetic differences between MAP-types are partially revised. Four of ten assumed MAP-S-specific large sequence polymorphism regions (LSP^Ss) are still present in MAP-C strains; new LSP^Ss were identified. Independently of the regional origin of the strains, the number of individual CDSs and single nucleotide variants confirms the strong similarity of MAP-C strains and shows higher diversity among MAP-S strains. This study gives ambiguous results regarding the hypothesis that MAP-S is the evolutionary intermediate between MAH and MAP-C, but it clearly shows a higher similarity of MAP to MAH than to Mycobacterium intracellulare.

Genetic diversity of Mycobacterium avium subspecies paratuberculosis and the influence of strain type on infection and pathogenesis: a review

Mycobacterium avium subspecies paratuberculosis (Map) is an important pathogen that causes a chronic, progressive granulomatous enteritis known as Johne's disease or paratuberculosis. The disease is endemic in many parts of the world and responsible for considerable losses to the livestock and associated industries. Diagnosis and control are problematic, due mostly to the long incubation period of the disease when infected animals show no clinical signs and are difficult to detect, and the ability of the organism to survive and persist in the environment. The existence of phenotypically distinct strains of Map has been known since the 1930s but the genetic differentiation of Map strain types has been challenging and only recent technologies have proven sufficiently discriminative for strain comparisons, tracing the sources of infection and epidemiological studies. It is important to understand the differences that exist between Map strains and how they influence both development and transmission of disease. This information is required to develop improved diagnostics and effective vaccines for controlling Johne's disease. Here I review the current classification of Map strain types, the sources of the genetic variability within strains, growth characteristics and epidemiological traits associated with strain type and the influence of strain type on infection and pathogenicity.

Typing of Mycobacterium avium subspecies paratuberculosis isolates from Newfoundland using fragment analysis

Short Sequence Repeat (SSR) typing of Mycobacterium avium subspecies paratuberculosis (Map) isolates is one of the most commonly used method for genotyping this pathogen. Currently used techniques have challenges in analyzing mononucleotide repeats >15 bp, which include some of the Map SSRs. Fragment analysis is a relatively simple technique, which can accurately measure the size of DNA fragments and can be used to calculate the repeat length of the target SSR loci. In the present study, fragment analysis was used to analyze 4 Map SSR loci known to provide sufficient discriminatory power to determine the relationship between Map isolates. Eighty-five Map isolates from 18 animals from the island of Newfoundland were successfully genotyped using fragment analysis. To the best of our knowledge, this is the first report on Map SSR diversity from Newfoundland dairy farms. Previously unreported Map SSR-types or combinations were also identified during the course of the described work. In addition, multiple Map SSR-types were isolated from a single animal in many cases, which is not a common finding.

Disruption of Mycobacterium avium subsp. paratuberculosis-specific genes impairs in vivo fitness

Background

Mycobacterium avium subsp. paratuberculosis (MAP) is an obligate intracellular pathogen that infects many ruminant species. The acquisition of foreign genes via horizontal gene transfer has been postulated to contribute to its pathogenesis, as these genetic elements are absent from its putative ancestor, M. avium subsp. hominissuis (MAH), an environmental organism with lesser pathogenicity. In this study, high-throughput sequencing of MAP transposon libraries were analyzed to qualitatively and quantitatively determine the contribution of individual genes to bacterial survival during infection.

Results

Out of 52384 TA dinucleotides present in the MAP K-10 genome, 12607 had a MycoMarT7 transposon in the input pool, interrupting 2443 of the 4350 genes in the MAP genome (56%). Of 96 genes situated in MAP-specific genomic islands, 82 were disrupted in the input pool, indicating that MAP-specific genomic regions are dispensable for in vitro growth (odds ratio = 0.21). Following 5 independent in vivo infections with this pool of mutants, the correlation between output pools was high for 4 of 5 (R = 0.49 to 0.61) enabling us to define genes whose disruption reproducibly reduced bacterial fitness in vivo. At three different thresholds for reduced fitness in vivo, MAP-specific genes were over-represented in the list of predicted essential genes. We also identified additional genes that were severely depleted after infection, and several of them have orthologues that are essential genes in M. tuberculosis.

Conclusions

This work indicates that the genetic elements required for the in vivo survival of MAP represent a combination of conserved mycobacterial virulence genes and MAP-specific genes acquired via horizontal gene transfer. In addition, the in vitro and in vivo essential genes identified in this study may be further characterized to offer a better understanding of MAP pathogenesis, and potentially contribute to the discovery of novel therapeutic and vaccine targets.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-415) contains supplementary material, which is available to authorized users.

SNP genotyping of animal and human derived isolates of Mycobacterium avium subsp. paratuberculosis

Mycobacterium avium subsp. paratuberculosis (MAP) is the aetiological agent of Johne's disease (JD), a chronic granulomatous enteritis that affects ruminants worldwide. While the ability of MAP to cause disease in animals is clear, the role of this bacterium in human inflammatory bowel diseases remains unresolved. Previous whole genome sequencing of MAP isolates derived from human and three animal hosts showed that human isolates were genetically similar and showed a close phylogenetic relationship to one bovine isolate. In contrast, other animal derived isolates were more genetically diverse. The present study aimed to investigate the frequency of this human strain across 52 wild-type MAP isolates, collected predominantly from Australia. A Luminex based SNP genotyping approach was utilised to genotype SNPs that had previously been shown to be specific to the human, bovine or ovine isolate types. Fourteen SNPs were initially evaluated across a reference panel of isolates with known genotypes. A subset of seven SNPs was chosen for analysis within the wild-type collection. Of the seven SNPs, three were found to be unique to paediatric human isolates. No wild-type isolates contain these SNP alleles. Interestingly, and in contrast to the paediatric isolates, three additional adult human isolates (derived from adult Crohn's disease patients) also did not contain these SNP alleles. Furthermore we identified two SNPs, which demonstrate extensive polymorphism within the animal-derived MAP isolates. One of which appears unique to ovine and a single camel isolate. From this study we suggest the existence of genetic heterogeneity between human derived MAP isolates, some of which are highly similar to those derived from bovine hosts, but others of which are more divergent.

Genetic diversity and phylogeny of Mycobacterium avium

Mycobacterium avium, one of the species of the M. avium complex (MAC), includes 4 subspecies, i.e., M. avium subsp. hominissuis (MAH), M. avium subsp. avium (MAA), M. avium subsp. silvaticum (MAS) and M. avium subsp. paratuberculosis (MAP), in turn classified into the S (sheep) and C (cattle) types. These subspecies, although closely related, represent distinct organisms, each endowed with specific pathogenetic and host range characteristics, ranging from environmental opportunistic bacteria that cause infections in swine and immunocompromised patients to pathogens of birds and ruminants. The present review summarizes the basic epidemiological and pathological features of the M. avium subspecies, describes the major genomic events responsible of M. avium subspecies diversity (insertion sequences, sequence variations in specific chromosome loci or genes, deletions, duplications and insertions of large genomic regions) and then reconstructs the phylogenetic relationships among the M. avium subspecies.

Novel feature of Mycobacterium avium subsp. paratuberculosis, highlighted by characterization of the heparin-binding hemagglutinin adhesin

Mycobacterium avium subsp. paratuberculosis comprises two genotypically defined groups, known as the cattle (C) and sheep (S) groups. Recent studies have reported phenotypic differences between M. avium subsp. paratuberculosis groups C and S, including growth rates, infectivity for macrophages, and iron metabolism. In this study, we investigated the genotypes and biological properties of the virulence factor heparin-binding hemagglutinin adhesin (HBHA) for both groups. In Mycobacterium tuberculosis, HBHA is a major adhesin involved in mycobacterium-host interactions and extrapulmonary dissemination of infection. To investigate HBHA in M. avium subsp. paratuberculosis, we studied hbhA polymorphisms by fragment analysis using the GeneMapper technology across a large collection of isolates genotyped by mycobacterial interspersed repetitive-unit-variable-number tandem-repeat (MIRU-VNTR) and IS900 restriction fragment length polymorphism (RFLP-IS900) analyses. Furthermore, we analyzed the structure-function relationships of recombinant HBHA proteins of types C and S by heparin-Sepharose chromatography and surface plasmon resonance (SPR) analyses. In silico analysis revealed two forms of HBHA, corresponding to the prototype genomes for the C and S types of M. avium subsp. paratuberculosis. This observation was confirmed using GeneMapper on 85 M. avium subsp. paratuberculosis strains, including 67 strains of type C and 18 strains of type S. We found that HBHAs from all type C strains contain a short C-terminal domain, while those of type S present a long C-terminal domain, similar to that produced by Mycobacterium avium subsp. avium. The purification of recombinant HBHA from M. avium subsp. paratuberculosis of both types by heparin-Sepharose chromatography highlighted a correlation between their affinities for heparin and the lengths of their C-terminal domains, which was confirmed by SPR analysis. Thus, types C and S of M. avium subsp. paratuberculosis may be distinguished by the types of HBHA they produce, which differ in size and adherence properties, thereby providing new evidence that strengthens the genotypic differences between the C and S types of M. avium subsp. paratuberculosis.

Genomic variations associated with attenuation in Mycobacterium avium subsp. paratuberculosis vaccine strains

BACKGROUND

Mycobacterium avium subspecies paratuberculosis (MAP) whole cell vaccines have been widely used tools in the control of Johne's disease in animals despite being unable to provide complete protection. Current vaccine strains derive from stocks created many decades ago; however their genotypes, underlying mechanisms and relative degree of their attenuation are largely unknown.

RESULTS

Using mouse virulence studies we confirm that MAP vaccine strains 316 F, II and 2e have diverse but clearly attenuated survival and persistence characteristics compared with wild type strains. Using a pan genomic microarray we characterise the genomic variations in a panel of vaccine strains sourced from stocks spanning over 40 years of maintenance. We describe multiple genomic variations specific for individual vaccine stocks in both deletion (26-32 Kbp) and tandem duplicated (11-40 Kbp) large variable genomic islands and insertion sequence copy numbers. We show individual differences suitable for diagnostic differentiation between vaccine and wild type genotypes and provide evidence for functionality of some of the deleted MAP-specific genes and their possible relation to attenuation.

CONCLUSIONS

This study shows how culture environments have influenced MAP genome diversity resulting in large tandem genomic duplications, deletions and transposable element activity. In combination with classical selective systematic subculture this has led to fixation of specific MAP genomic alterations in some vaccine strain lineages which link the resulting attenuated phenotypes with deficiencies in high reactive oxygen species handling.

Inter- and intra-subtype genotypic differences that differentiate Mycobacterium avium subspecies paratuberculosis strains

Background

Mycobacterium avium subspecies paratuberculosis (Map) is the aetiological agent of Johne’s disease or paratuberculosis and is included within the Mycobacterium avium complex (MAC). Map strains are of two major types often referred to as ‘Sheep’ or ‘S-type’ and ‘Cattle’ or ‘C-type’. With the advent of more discriminatory typing techniques it has been possible to further classify the S-type strains into two groups referred to as Type I and Type III. This study was undertaken to genotype a large panel of S-type small ruminant isolates from different hosts and geographical origins and to compare them with a large panel of well documented C-type isolates to assess the genetic diversity of these strain types. Methods used included Mycobacterial Interspersed Repetitive Units - Variable-Number Tandem Repeat analysis (MIRU-VNTR), analysis of Large Sequence Polymorphisms by PCR (LSP analysis), Single Nucleotide Polymorphism (SNP) analysis of gyr genes, Pulsed-Field Gel Electrophoresis (PFGE) and Restriction Fragment Length Polymorphism analysis coupled with hybridization to IS900 (IS900-RFLP) analysis.

Results

The presence of LSP^A4 and absence of LSP^A20 was confirmed in all 24 Map S-type strains analysed. SNPs within the gyr genes divided the S-type strains into types I and III. Twenty four PFGE multiplex profiles and eleven different IS900-RFLP profiles were identified among the S-type isolates, some of them not previously published. Both PFGE and IS900-RFLP segregated the S-type strains into types I and III and the results concurred with those of the gyr SNP analysis. Nine MIRU-VNTR genotypes were identified in these isolates. MIRU-VNTR analysis differentiated Map strains from other members of Mycobacterium avium Complex, and Map S-type from C-type but not type I from III. Pigmented Map isolates were found of type I or III.

Conclusion

This is the largest panel of S-type strains investigated to date. The S-type strains could be further divided into two subtypes, I and III by some of the typing techniques (IS900-RFLP, PFGE and SNP analysis of the gyr genes). MIRU-VNTR did not divide the strains into the subtypes I and III but did detect genetic differences between isolates within each of the subtypes. Pigmentation is not exclusively associated with type I strains.

Maize (Zea mays L.) genome diversity as revealed by RNA-sequencing

Maize is rich in genetic and phenotypic diversity. Understanding the sequence, structural, and expression variation that contributes to phenotypic diversity would facilitate more efficient varietal improvement. RNA based sequencing (RNA-seq) is a powerful approach for transcriptional analysis, assessing sequence variation, and identifying novel transcript sequences, particularly in large, complex, repetitive genomes such as maize. In this study, we sequenced RNA from whole seedlings of 21 maize inbred lines representing diverse North American and exotic germplasm. Single nucleotide polymorphism (SNP) detection identified 351,710 polymorphic loci distributed throughout the genome covering 22,830 annotated genes. Tight clustering of two distinct heterotic groups and exotic lines was evident using these SNPs as genetic markers. Transcript abundance analysis revealed minimal variation in the total number of genes expressed across these 21 lines (57.1% to 66.0%). However, the transcribed gene set among the 21 lines varied, with 48.7% expressed in all of the lines, 27.9% expressed in one to 20 lines, and 23.4% expressed in none of the lines. De novo assembly of RNA-seq reads that did not map to the reference B73 genome sequence revealed 1,321 high confidence novel transcripts, of which, 564 loci were present in all 21 lines, including B73, and 757 loci were restricted to a subset of the lines. RT-PCR validation demonstrated 87.5% concordance with the computational prediction of these expressed novel transcripts. Intriguingly, 145 of the novel de novo assembled loci were present in lines from only one of the two heterotic groups consistent with the hypothesis that, in addition to sequence polymorphisms and transcript abundance, transcript presence/absence variation is present and, thereby, may be a mechanism contributing to the genetic basis of heterosis.

A 16 kb naturally occurring genomic deletion including mce and PPE genes in Mycobacterium avium subspecies paratuberculosis isolates from goats with Johne's disease

In this study we characterise the genomic and transcriptomic variability of a natural deletion strain of Mycobacterium avium subspecies paratuberculosis (MAP) prevalent in Spanish Guadarrama goats. Using a pan-genome microarray including MAP and M. avium subspecies hominissuis 104 genomes (MAPAC) we demonstrate the genotype to be MAP Type II with a single deletion of 19 contiguous ORFs (16 kb) including a complete mammalian cell entry (mce7_1) operon and adjacent proline-glutamic acid (PE)/proline-proline-glutamic acid (PPE) genes. A deletion specific PCR test was developed and a subsequent screening identified four goat herds infected with the variant strain. Each was located in central Spain and showed epidemiological links suggestive of transmission between herds. A majority of animals infected with the variant manifested a paucibacillary form of the disease. Comparisons between virulent complete genome compliment strains isolated from multibacillary diseased goats and the MAP variant strain during entry into activated macrophages demonstrated an increased sensitivity in the variant to intracellular killing in human and ovine macrophages. As PPE and mce genes are associated with mycobacterial virulence and pathogenesis we investigated the interplay of these gene sets during cell entry using the MAPAC array. This showed significant differential transcriptome profiles compared to full genome complement MAP controls that included changes in other undeleted mce operons and PE/PPE genes, esx-like signalling operons and stress response/fatty acid metabolism pathways. This strain represents the first report of a MAP Type II genotype with significant natural genomic deletions which remains able to cause disease and is transmissible in goats.

Genome sequencing of ovine isolates of Mycobacterium avium subspecies paratuberculosis offers insights into host association

Background

The genome of Mycobacterium avium subspecies paratuberculosis (MAP) is remarkably homogeneous among the genomes of bovine, human and wildlife isolates. However, previous work in our laboratories with the bovine K-10 strain has revealed substantial differences compared to sheep isolates. To systematically characterize all genomic differences that may be associated with the specific hosts, we sequenced the genomes of three U.S. sheep isolates and also obtained an optical map.

Results

Our analysis of one of the isolates, MAP S397, revealed a genome 4.8 Mb in size with 4,700 open reading frames (ORFs). Comparative analysis of the MAP S397 isolate showed it acquired approximately 10 large sequence regions that are shared with the human M. avium subsp. hominissuis strain 104 and lost 2 large regions that are present in the bovine strain. In addition, optical mapping defined the presence of 7 large inversions between the bovine and ovine genomes (~ 2.36 Mb). Whole-genome sequencing of 2 additional sheep strains of MAP (JTC1074 and JTC7565) further confirmed genomic homogeneity of the sheep isolates despite the presence of polymorphisms on the nucleotide level.

Conclusions

Comparative sequence analysis employed here provided a better understanding of the host association, evolution of members of the M. avium complex and could help in deciphering the phenotypic differences observed among sheep and cattle strains of MAP. A similar approach based on whole-genome sequencing combined with optical mapping could be employed to examine closely related pathogens. We propose an evolutionary scenario for M. avium complex strains based on these genome sequences.

Genome-wide analysis of the emerging infection with Mycobacterium avium subspecies paratuberculosis in the Arabian camels (Camelus dromedarius)

Mycobacterium avium subspecies paratuberculosis (M. ap) is the causative agent of paratuberculosis or Johne's disease (JD) in herbivores with potential involvement in cases of Crohn's disease in humans. JD is spread worldwide and is economically important for both beef and dairy industries. Generally, pathogenic ovine strains (M. ap-S) are mainly found in sheep while bovine strains (M. ap-C) infect other ruminants (e.g. cattle, goat, deer), as well as sheep. In an effort to characterize this emerging infection in dromedary/Arabian camels, we successfully cultured M. ap from several samples collected from infected camels suffering from chronic, intermittent diarrhea suggestive of JD. Gene-based typing of isolates indicated that all isolates belong to sheep lineage of strains of M. ap (M. ap-S), suggesting a putative transmission from infected sheep herds. Screening sheep and goat herds associated with camels identified the circulation of this type in sheep but not goats. The current genome-wide analysis recognizes these camel isolates as a sub-lineage of the sheep strain with a significant number of single nucleotide polymorphisms (SNPs) between sheep and camel isolates (∼1000 SNPs). Such polymorphism could represent geographical differences among isolates or host adaptation of M. ap during camel infection. To our knowledge, this is the first attempt to examine the genomic basis of this emerging infection in camels with implications on the evolution of this important pathogen. The sequenced genomes of M. ap isolates from camels will further assist our efforts to understand JD pathogenesis and the dynamic of disease transmission across animal species.

Proteome-determined type-specific proteins of Mycobacterium avium subspecies paratuberculosis

Mycobacterium avium subspecies paratuberculosis (M. a. paratuberculosis) is a pathogen of ruminants, causing paratuberculosis (characterized by severe emaciation). The disease is endemic in many countries including the UK and places a severe economic burden on the global livestock industry. Two types of M. a. paratuberculosis can be classified by pulsed-field electrophoresis (I/III and II), which are phenotypically distinct and appear to have different host preferences. Proteomes of Type I and Type II M. a. paratuberculosis were analyzed by 2-D gel electrophoresis to determine if any significant differences existed between the subtypes. Seven different strains of Type I and 18 strains of Type II were analyzed and compared to detect type-specific differences. These 'type-specific' differences existed regardless of growth phase and were also exhibited in cells isolated directly from pathogenic lesions. Twenty-three spots predominated on the Type I profile, from which 17 proteins were identified. Twenty-one spots predominated on the Type II profile, from which 16 proteins were identified. None of the proteins identified as differentially represented on the profiles of Type I or Type II corresponded to open reading frames of the defining genomic regions as previously described for the Type I (sheep) and Type II (cattle). Sequence polymorphisms existing in Type I and II strains were identified in some open reading frames or regulatory regions of genes that correspond to proteins expressed in a type-specific fashion. The consequence of these is discussed in relation to protein expression and their impact on the type phenotype is discussed.

Exploring the zoonotic potential of Mycobacterium avium subspecies paratuberculosis through comparative genomics

A comparative genomics approach was utilised to compare the genomes of Mycobacterium avium subspecies paratuberculosis (MAP) isolated from early onset paediatric Crohn's disease (CD) patients as well as Johne's diseased animals. Draft genome sequences were produced for MAP isolates derived from four CD patients, one ulcerative colitis (UC) patient, and two non-inflammatory bowel disease (IBD) control individuals using Illumina sequencing, complemented by comparative genome hybridisation (CGH). MAP isolates derived from two bovine and one ovine host were also subjected to whole genome sequencing and CGH. All seven human derived MAP isolates were highly genetically similar and clustered together with one bovine type isolate following phylogenetic analysis. Three other sequenced isolates (including the reference bovine derived isolate K10) were genetically distinct. The human isolates contained two large tandem duplications, the organisations of which were confirmed by PCR. Designated vGI-17 and vGI-18 these duplications spanned 63 and 109 open reading frames, respectively. PCR screening of over 30 additional MAP isolates (3 human derived, 27 animal derived and one environmental isolate) confirmed that vGI-17 and vGI-18 are common across many isolates. Quantitative real-time PCR of vGI-17 demonstrated that the proportion of cells containing the vGI-17 duplication varied between 0.01 to 15% amongst isolates with human isolates containing a higher proportion of vGI-17 compared to most animal isolates. These findings suggest these duplications are transient genomic rearrangements. We hypothesise that the over-representation of vGI-17 in human derived MAP strains may enhance their ability to infect or persist within a human host by increasing genome redundancy and conferring crude regulation of protein expression across biologically important regions.

Amplified fragment length polymorphism reveals specific epigenetic distinctions between Mycobacterium avium subspecies paratuberculosis isolates of various isolation types

Amplified fragment length polymorphism (AFLP) was employed as a genetic analysis tool for the study of the genetic relatedness of Mycobacterium avium subsp. paratuberculosis isolates harvested from bovine fecal samples and from bovine or human tissues. This analysis revealed genetic differences between these two isolate types that were confirmed through cluster analysis. Dendrogram analysis separated these two isolate types based on the isolation scheme (tissue-associated versus fecal M. avium subsp. paratuberculosis isolates). Further sequence analysis of unique genetic regions from each isolation type revealed no genetic sequence differences. However, Clustal DNA alignments identified AFLP restriction enzyme sites that were undigested in the tissue-associated isolates. AFLP analysis also disclosed that the same AFLP restriction sites were digested in all of the fecal isolates. Sequence analysis further revealed a consensus sequence upstream of the undigested restriction sites for possible methyltransferase recognition in the tissue-associated M. avium subsp. paratuberculosis isolates.

Mycobacterium avium subsp. paratuberculosis in dairy products, meat, and drinking water

Mycobacterium avium subsp. paratuberculosis (Map) is the cause of Johne's disease, a chronic infection of the gut, in ruminant animals that provide milk and/or meat for human consumption. Map also may be involved in Crohn's disease and type 1 diabetes in humans. Although the role of Map in human diseases has not been established, minimizing the exposure of humans to the organism is considered desirable as a precautionary measure. Infected animals can shed Map in feces and milk, and the organism can become disseminated in tissues remote from the gut and its associated lymph nodes. The presence of at least some Map in raw milk and meat and in natural waters is likely, but the numbers of Map in those foods and waters should be reduced through cooking or purification. The available information relating to Map in milk and dairy products, meats, and drinking water is reviewed here for assessment of the risks of exposure to Map from consumption of such foods and water.

Comparative genomics of two newly isolated Dehalococcoides strains and an enrichment using a genus microarray

Comparative genomics of Dehalococcoides strains and an enrichment were performed using a microarray targeting genes from all available sequenced genomes of the Dehalococcoides genus. The microarray was designed with 4305 probe sets to target 98.6% of the open-reading frames from strains 195, CBDB1, BAV1 and VS. The microarrays were validated and applied to query the genomes of two recently isolated Dehalococcoides strains, ANAS1 and ANAS2, and their enrichment source (ANAS) to understand the genome-physiology relationships. Strains ANAS1 and ANAS2 can both couple the reduction of trichloroethene, cis-dichloroethene (DCE) and 1,1-DCE, but not tetrachloroethene and trans-DCE with growth, whereas only strain ANAS2 couples vinyl chloride reduction to growth. Comparative genomic analysis showed that the genomes of both strains are similar to each other and to strain 195, except for genes that are within the previously defined integrated elements or high-plasticity regions. Combined results of the two isolates closely matched the results obtained using genomic DNA of the ANAS enrichment. The genome similarities, together with the distinct chlorinated ethene usage of strains ANAS1, ANAS2 and 195 demonstrate that closely phylogenetically related strains can be physiologically different. This incongruence between physiology and core genome phylogeny seems to be related to the presence of distinct reductive dehalogenase-encoding genes with assigned chlorinated ethene functions (pceA, tceA in strain 195; tceA in strain ANAS1; vcrA in strain ANAS2). Overall, the microarrays are a valuable high-throughput tool for comparative genomics of unsequenced Dehalococcoides-containing samples to provide insights into their gene content and dechlorination functions.

Genomic and transcriptomic studies in Mycobacterium avium subspecies paratuberculosis

Microarray technology is an important tool in functional genomic research. It has enabled a deeper analysis of genomic diversity among bacteria belonging to the Mycobacterium avium complex (MAC). In addition, the expression of thousands of genes can be studied simultaneously in a single experiment. With the complete genome sequence of a bovine isolate of M. avium subspecies paratuberculosis, and the independent construction of DNA microarrays in our laboratories, transcriptomic studies for this veterinary pathogen are now possible. Furthermore, the bovine genome sequence project is completed and bovine arrays have been developed to examine host responses to infection with M. avium subsp. paratuberculosis. Collectively, genomic and transcriptomic data has yielded novel insights surrounding the genetic regulation and biology of Johne's disease.

Probing the pan-genome of Listeria monocytogenes: new insights into intraspecific niche expansion and genomic diversification

Background

Bacterial pathogens often show significant intraspecific variations in ecological fitness, host preference and pathogenic potential to cause infectious disease. The species of Listeria monocytogenes, a facultative intracellular pathogen and the causative agent of human listeriosis, consists of at least three distinct genetic lineages. Two of these lineages predominantly cause human sporadic and epidemic infections, whereas the third lineage has never been implicated in human disease outbreaks despite its overall conservation of many known virulence factors.

Results

Here we compare the genomes of 26 L. monocytogenes strains representing the three lineages based on both in silico comparative genomic analysis and high-density, pan-genomic DNA array hybridizations. We uncover 86 genes and 8 small regulatory RNAs that likely make L. monocytogenes lineages differ in carbohydrate utilization and stress resistance during their residence in natural habitats and passage through the host gastrointestinal tract. We also identify 2,330 to 2,456 core genes that define this species along with an open pan-genome pool that contains more than 4,052 genes. Phylogenomic reconstructions based on 3,560 homologous groups allowed robust estimation of phylogenetic relatedness among L. monocytogenes strains.

Conclusions

Our pan-genome approach enables accurate co-analysis of DNA sequence and hybridization array data for both core gene estimation and phylogenomics. Application of our method to the pan-genome of L. monocytogenes sheds new insights into the intraspecific niche expansion and evolution of this important foodborne pathogen.

PanCGHweb: a web tool for genotype calling in pangenome CGH data

Summary: A pangenome is the total of genes present in strains of the same species. Pangenome microarrays allow determining the genomic content of bacterial strains more accurately than conventional comparative genome hybridization microarrays. PanCGHweb is the first tool that effectively calls genotype based on pangenome microarray data.

Availability: PanCGHweb, the web tool is accessible from: http://bamics2.cmbi.ru.nl/websoftware/pancgh/

Contact:sacha.vanhijum@nizo.nl

Rapid identification and differentiation of Mycobacterium avium subspecies paratuberculosis types by use of real-time PCR and high-resolution melt analysis of the MAP1506 locus

High-resolution melt (HRM) analysis can identify sequence polymorphisms by comparing the melting curves of amplicons generated by real-time PCR amplification. We describe the application of this technique to identify Mycobacterium avium subspecies paratuberculosis types I, II, and III. The HRM approach was based on type-specific nucleotide sequences in MAP1506, a member of the PPE (proline-proline-glutamic acid) gene family.

Molecular characterization of Mycobacterium avium subspecies paratuberculosis Types II and III isolates by a combination of MIRU-VNTR loci

Mycobacterial interspersed repetitive units and variable number tandem repeats typing (MIRU-VNTR) is a useful technique that has been recently applied to characterize members of the Mycobacterium avium complex (MAC). The aim of this study was to examine the genetic variability among a collection of Spanish M. avium subspecies paratuberculosis (M. a. paratuberculosis) isolates with a combination of MIRU-VNTR loci. For this purpose we tested six MIRU-VNTR loci (MIRU-2, MIRU-3, VNTR-25, VNTR-32, VNTR-292 and VNTR-259) in 70 M. a. paratuberculosis isolates of Types II and III that were recovered from 22 Spanish localities during a nine-year period (1998-2007). The combination of five loci (MIRU-2, MIRU-3, VNTR-25, VNTR-32 and VNTR-259) enabled the differentiation of 12 allelic profiles, with a resulting Hunter and Gaston discriminatory index (HGDI) of 0.84. Moreover, we obtained MIRU-VNTR patterns that were unique for each of the M. a. paratuberculosis types analyzed (II and III); other patterns were host-related or restricted to geographic areas. Therefore, this MIRU-VNTR approach could be a useful sub-typing molecular tool in order to get a better sense of the epidemiology of Johne's disease.

A rhodanine agent active against non-replicating intracellular Mycobacterium avium subspecies paratuberculosis

Background

Antibiotic therapy targeting chronic mycobacterial disease is often ineffective due to problems with the emergence of drug resistance and non-replicating persistent intracellular antibiotic resistant phenotypes. Strategies which include agents able to enhance host cell killing mechanisms could represent an alternative to conventional methods with the potential for host clearance if active against dormant phenotypes. Investigations of agents with potential activity against non-replicating mycobacteria however are restricted due to a need for assays that can assess bacterial viability without having to culture.

Results

This study describes the development and use of a pre16S ribosomal gene RNA/DNA ratio viability assay which is independent of the need for culture, supported by a novel thin layer accelerated mycobacterial colony forming method for determining viability and culturability of MAP in intracellular environments. We describe the use of these tools to demonstrate intracellular killing activity of a novel rhodanine agent (D157070) against the intracellular pathogen Mycobacterium avium subspecies paratuberculosis (MAP) and show that the culturability of MAP decreases relative to its viability on intracellular entry suggesting the induction of a non-culturable phenotype. We further demonstrate that D157070, although having no direct activity against the culturability of extracellular MAP, can bind to cultured MAP cells and has significant influence on the MAP transcriptome, particularly with respect of δ^L associated genes. D157070 is shown to be taken up by bovine and human cells and able to enhance host cell killing, as measured by significant decreases in both culturability and viability of intracellular MAP.

Conclusions

This work suggests that pre16srRNA gene ratios represent a viable method for studying MAP viability. In addition, the rhodanine agent D157070 tested is non-toxic and enhances cell killing activity against both growing and latent MAP phenotypes.

Single nucleotide polymorphisms in the IS900 sequence of Mycobacterium avium subsp. paratuberculosis are strain type specific

Insertion sequence IS900 is used as a target for the identification of Mycobacterium avium subsp. paratuberculosis. Previous reports have revealed single nucleotide polymorphisms within IS900. This study, which analyzed the IS900 sequences of a panel of isolates representing M. avium subsp. paratuberculosis strain types I, II, and III, revealed conserved type-specific polymorphisms that could be utilized as a tool for diagnostic and epidemiological purposes.