Comparative genomic analysis of Mycobacterium avium subspecies obtained from multiple host species

The pan-genome of Mycobacterium avium subsp. paratuberculosis (Map) confirms ancestral lineage and reveals gene rearrangements within Map Type S

BACKGROUND

To date genomic studies on Map have concentrated on Type C strains with only a few Type S strains included for comparison. In this study the entire pan-genome of 261 Map genomes (205 Type C, 52 Type S and 4 Type B) and 7 Mycobacterium avium complex (Mac) genomes were analysed to identify genomic similarities and differences between the strains and provide more insight into the evolutionary relationship within this Mycobacterial species.

RESULTS

Our analysis of the core genome of all the Map isolates identified two distinct lineages, Type S and Type C Map that is consistent with previous phylogenetic studies of Map. Pan-genome analysis revealed that Map has a larger accessory genome than Mycobacterium avium subsp. avium (Maa) and Type C Map has a larger accessory genome than Type S Map. In addition, we found large rearrangements within Type S strains of Map and little to none in Type C and Type B strains. There were 50 core genes identified that were unique to Type S Map and there were no unique core genes identified between Type B and Type C Map strains. In Type C Map we identified an additional CE10 CAZyme class which was identified as an alpha/beta hydrolase and an additional polyketide and non-ribosomal peptide synthetase cluster. Consistent with previous analysis no plasmids and only incomplete prophages were identified in the genomes of Map. There were 45 hypothetical CRISPR elements identified with no associated cas genes.

CONCLUSION

This is the most comprehensive comparison of the genomic content of Map isolates to date and included the closing of eight Map genomes. The analysis revealed that there is greater variation in gene synteny within Type S strains when compared to Type C indicating that the Type C Map strain emerged after Type S. Further analysis of Type C and Type B genomes revealed that they are structurally similar with little to no genetic variation and that Type B Map may be a distinct clade within Type C Map and not a different strain type of Map. The evolutionary lineage of Maa and Map was confirmed as emerging after M. hominissuis.

Vitamin D3 alters macrophage phenotype and endosomal trafficking markers in dairy cattle naturally infected with Mycobacterium avium subsp. paratuberculosis

Macrophages are important host defense cells in ruminant paratuberculosis (Johne’s Disease; JD), a chronic enteritis caused by Mycobacterium avium subsp. paratuberculosis (MAP). Classical macrophage functions of pathogen trafficking, degradation, and antigen presentation are interrupted in mycobacterial infection. Immunologic stimulation by 25-hydroxyvitamin D₃ (25(OH)D₃) and 1,25-dihydroxyvitamin D₃ (1,25(OH)₂D₃) enhances bovine macrophage function. The present study aimed to investigate the role of vitamin D₃ on macrophage phenotype and endosomal trafficking of MAP in monocyte-derived macrophages (MDMs) cultured from JD-, JD+ subclinical, and JD+ clinically infected cattle. MDMs were pre-treated 100 ng/ml 25(OH)D₃ or 4 ng/ml 1,25(OH)₂D₃ and incubated 24 hrs with MAP at 10:1 multiplicity of infection (MOI). In vitro MAP infection upregulated pro-inflammatory (M1) CD80 and downregulated resolution/repair (M2) CD163. Vitamin D₃ generally decreased CD80 and increased CD163 expression. Furthermore, early endosomal marker Rab5 was upregulated 140× across all stages of paratuberculosis infection following in vitro MAP infection; however, Rab5 was reduced in MAP-activated MDMs from JD+ subclinical and JD+ clinical cows compared to healthy controls. Rab7 expression decreased in control and clinical cows following MDM infection with MAP. Both forms of vitamin D₃ reduced Rab5 expression in infected MDMs from JD- control cows, while 1,25(OH)₂D₃ decreased Rab7 expression in JD- and JD+ subclinical animals regardless of MAP infection in vitro. Vitamin D₃ promoted phagocytosis in MDMs from JD- and JD+ clinical cows treated with either vitamin D₃ analog. Results from this study show exogenous vitamin D₃ influences macrophage M1/M2 polarization and Rab GTPase expression within MDM culture.

Global Phylogeny of Mycobacterium avium and Identification of Mutation Hotspots During Niche Adaptation

Mycobacterium avium is separated into four subspecies: M. avium subspecies avium (MAA), M. avium subspecies silvaticum (MAS), M. avium subspecies hominissuis (MAH), and M. avium subspecies paratuberculosis (MAP). Understanding the mechanisms of host and tissue adaptation leading to their clinical significance is vital to reduce the economic, welfare, and public health concerns associated with diseases they may cause in humans and animals. Despite substantial phenotypic diversity, the subspecies nomenclature is controversial due to high genetic similarity. Consequently, a set of 1,230 M. avium genomes was used to generate a phylogeny, investigate SNP hotspots, and identify subspecies-specific genes. Phylogeny reiterated the findings from previous work and established that Mycobacterium avium is a species made up of one highly diverse subspecies, known as MAH, and at least two clonal pathogens, named MAA and MAP. Pan-genomes identified coding sequences unique to each subspecies, and in conjunction with a mapping approach, mutation hotspot regions were revealed compared to the reference genomes for MAA, MAH, and MAP. These subspecies-specific genes may serve as valuable biomarkers, providing a deeper understanding of genetic differences between M. avium subspecies and the virulence mechanisms of mycobacteria. Furthermore, SNP analysis demonstrated common regions between subspecies that have undergone extensive mutations during niche adaptation. The findings provide insights into host and tissue specificity of this genetically conserved but phenotypically diverse species, with the potential to provide new diagnostic targets and epidemiological and therapeutic advances.

Application of Bayesian modeling for diagnostic assays of Mycobacterium avium subsp. paratuberculosis in sheep and goats flocks

Background

This study aimed to screen the sera of goats and sheep from flocks suspected of Mycobacterium avium subsp. paratuberculosis (MAP) infection by a newly standardized Mce-truncated ELISA (Mt-ELISA) kit for the detection of antibodies against MAP. Four diagnostic applied tests were evaluated including Indigenous plate-ELISA (IP-ELISA), Mt-ELISA, fecal Polymerase Chain Reaction (f-PCR) and fecal culture (FC).

Materials and methods

Assuming the absence of a gold standard, latent-class models in a Bayesian framework were used to estimate the diagnostic accuracy of the four tests for MAP.

Results

Mt-ELISA had higher Sensitivity (Se) in sheep (posterior median: 0.68 (95% Probability Interval (PI): 0.43–0.95), while IP-ELISA recorded the highest Se in goats as 0.83 (95% PI, 0.61–0.97). The f-PCR Se estimate slightly differed between species [sheep 0.36 (0.19–0.58), goats 0.19 (0.08–0.35)], while the Se of FC was similar between species [sheep 0.29 (0.15–0.51), goats 0.27 (0.13–0.45)]. The specificity estimates for all tests were high, close to unity, and similar between species.

Conclusion

Overall, the results showed that the Mt-ELISA method can be used for MAP detection in small ruminants’ flocks.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12917-022-03141-7.

Diagnostic Sequences That Distinguish M. avium Subspecies Strains

Over a decade ago Mycobacterium avium subspecies paratuberculosis (Map) specific genes were initially identified in a whole genome context by comparing draft genome sequences of Map strain K-10 with Mycobacterium avium subspecies hominissuis (Mah) strain 104. This resulted in identification of 32 Map specific genes, not including repetitive elements, based on the two-genome comparison. The goal of this study was to define a more complete catalog of M. avium subspecies-specific genes. This is important for obtaining additional diagnostic targets for Johne's disease detection and for understanding the unique biology, evolution and niche adaptation of these organisms. There are now over 28 complete genome sequences representing three M. avium subspecies, including avium (Maa), Mah, and Map. We have conducted a comprehensive comparison of these genomes using two independent pan genomic comparison tools, PanOCT and Roary. This has led to the identification of more than 250 subspecies defining genes common to both analyses. The majority of these genes are arranged in clusters called genomic islands. We further reduced the number of diagnostic targets by excluding sequences having high BLAST similarity to other mycobacterial species recently added to the National Center for Biotechnology Information database. Genes identified as diagnostic following these bioinformatic approaches were further tested by DNA amplification PCR on an additional 20 M. avium subspecies strains. This combined approach confirmed 86 genes as Map-specific, seven as Maa-specific and three as Mah-specific. A single-tube PCR reaction was conducted as a proof of concept method to quickly distinguish M. avium subspecies strains. With these novel data, researchers can classify isolates in their freezers, quickly characterize clinical samples, and functionally analyze these unique genes.

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.

Genetic Diversity Among Mycobacterium avium Subspecies Revealed by Analysis of Complete Genome Sequences

Mycobacterium avium comprises four subspecies that contain both human and veterinary pathogens. At the inception of this study, twenty-eight M. avium genomes had been annotated as RefSeq genomes, facilitating direct comparisons. These genomes represent strains from around the world and provided a unique opportunity to examine genome dynamics in this species. Each genome was confirmed to be classified correctly based on SNP genotyping, nucleotide identity and presence/absence of repetitive elements or other typing methods. The Mycobacterium avium subspecies paratuberculosis (Map) genome size and organization was remarkably consistent, averaging 4.8 Mb with a variance of only 29.6 kb among the 13 strains. Comparing recombination events along with the larger genome size and variance observed among Mycobacterium avium subspecies avium (Maa) and Mycobacterium avium subspecies hominissuis (Mah) strains (collectively termed non-Map) suggests horizontal gene transfer occurs in non-Map, but not in Map strains. Overall, M. avium subspecies could be divided into two major sub-divisions, with the Map type II (bovine strains) clustering tightly on one end of a phylogenetic spectrum and Mah strains clustering more loosely together on the other end. The most evolutionarily distinct Map strain was an ovine strain, designated Telford, which had >1,000 SNPs and showed large rearrangements compared to the bovine type II strains. The Telford strain clustered with Maa strains as an intermediate between Map type II and Mah. SNP analysis and genome organization analyses repeatedly demonstrated the conserved nature of Map versus the mosaic nature of non-Map M. avium strains. Finally, core and pangenomes were developed for Map and non-Map strains. A total of 80% Map genes belonged to the Map core genome, while only 40% of non-Map genes belonged to the non-Map core genome. These genomes provide a more complete and detailed comparison of these subspecies strains as well as a blueprint for how genetic diversity originated.

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.

First molecular epidemiological study of Mycobacterium avium subsp. paratuberculosis in cattle and buffalo from different regions of Brazil

Paratuberculosis is an incurable disease in ruminants with great worldwide economic impact, caused by Mycobacterium avium subsp. paratuberculosis (MAP). The objective of this study was to carry out a study of the molecular epidemiology of the MAP using the restriction enzyme analysis (REA) technique of IS1311 MAP region in biological samples of feces, intestinal tissue, and mesenteric lymph nodes of cattle and buffaloes from six Brazilian states. In total, 109 samples of feces and tissues of cattle and buffaloes were collected from animal paratuberculosis suspected. Twenty-five samples were positive in the detection of the DNA of the IS900 region of MAP and it was possible to type 18 strains in the analysis of the region IS1311, being 100% of them identified as belonging to subtype Bison MAP strain. This is the first epidemiological molecular study of MAP in Brazil. The results indicate that paratuberculosis is widespread in cattle and in buffaloes in several regions of Brazil, and the subtype Bison MAP strain was the only one identified in the samples analyzed in this study, demonstrating the similarity between the strains from different states tested. These results provide the necessary support for the implementation of paratuberculosis control strategies in cattle and buffaloes in Brazil.

Cell wall peptidolipids of Mycobacterium avium: from genetic prediction to exact structure of a nonribosomal peptide

Mycobacteria have a complex cell wall structure that includes many lipids; however, even within a single subspecies of Mycobacterium avium these lipids can differ. Total lipids from an M. avium subsp. paratuberculosis (Map) ovine strain (S-type) contained no identifiable glycopeptidolipids or lipopentapeptide (L5P), yet both lipids are present in other M. avium subspecies. We determined the genetic and phenotypic basis for this difference using sequence analysis as well as biochemical and physico-chemical approaches. This strategy showed that a nonribosomal peptide synthase, encoded by mps1, contains three amino acid specifying modules in ovine strains, compared to five modules in bovine strains (C-type). Sequence analysis predicted these modules would produce the tripeptide Phe-N-Methyl-Val-Ala with a lipid moiety, termed lipotripeptide (L3P). Comprehensive physico-chemical analysis of Map S397 extracts confirmed the structural formula of the native L3P as D-Phe-N-Methyl-L-Val-L-Ala-OMe attached in N-ter to a 20-carbon fatty acid chain. These data demonstrate that S-type strains, which are more adapted in sheep, produce a unique lipid. There is a dose-dependent effect observed for L3P on upregulation of CD25+ CD8 T cells from infected cows, while L5P effects were static. In contrast, L5P demonstrated a significantly stronger induction of CD25+ B cells from infected animals compared to L3P.

Methodological and Clinical Aspects of the Molecular Epidemiology of Mycobacterium tuberculosis and Other Mycobacteria

Molecular typing has revolutionized epidemiological studies of infectious diseases, including those of a mycobacterial etiology. With the advent of fingerprinting techniques, many traditional concepts regarding transmission, infectivity, or pathogenicity of mycobacterial bacilli have been revisited, and their conventional interpretations have been challenged. Since the mid-1990s, when the first typing methods were introduced, a plethora of other modalities have been proposed. So-called molecular epidemiology has become an essential subdiscipline of modern mycobacteriology. It serves as a resource for understanding the key issues in the epidemiology of tuberculosis and other mycobacterial diseases. Among these issues are disclosing sources of infection, quantifying recent transmission, identifying transmission links, discerning reinfection from relapse, tracking the geographic distribution and clonal expansion of specific strains, and exploring the genetic mechanisms underlying specific phenotypic traits, including virulence, organ tropism, transmissibility, or drug resistance. Since genotyping continues to unravel the biology of mycobacteria, it offers enormous promise in the fight against and prevention of the diseases caused by these pathogens. In this review, molecular typing methods for Mycobacterium tuberculosis and nontuberculous mycobacteria elaborated over the last 2 decades are summarized. The relevance of these methods to the epidemiological investigation, diagnosis, evolution, and control of mycobacterial diseases is discussed.

Molecular analysis and MIRU-VNTR typing of Mycobacterium avium subsp. avium, 'hominissuis' and silvaticum strains of veterinary origin

Besides Mycobacterium avium subsp. paratuberculosis (MAP), M. avium subsp. avium (MAA), M. avium subsp. silvaticum (MAS), and 'M. avium subsp. hominissuis' (MAH) are equally important members of M. avium complex, with worldwide distribution and zoonotic potential. Genotypic discrimination is a prerequisite to epidemiological studies which can facilitate disease prevention through revealing infection sources and transmission routes. The primary aim of this study was to identify the genetic diversity within 135 MAA, 62 MAS, and 84 MAH strains isolated from wild and domestic mammals, reptiles and birds. Strains were tested for the presence of large sequence polymorphism LSP(A)17 and were submitted to Mycobacterial interspersed repetitive units-variable-number tandem repeat (MIRU-VNTR) analysis at 8 loci, including MIRU1, 2, 3, and 4, VNTR25, 32, and 259, and MATR9. In 12 strains hsp65 sequence code type was also determined. LSP(A)17 was present only in 19.9% of the strains. All LSP(A)17 positive strains belonged to subspecies MAH. The discriminatory power of the MIRU-VNTR loci set used reached 0.9228. Altogether 54 different genotypes were detected. Within MAH, MAA, and MAS strains 33, 16, and 5 different genotypes were observed. The described genotypes were not restricted to geographic regions or host species, but proved to be subspecies specific. Our knowledge about MAS is limited due to isolation and identification difficulties. This is the first study including a large number of MAS field strains. Our results demonstrate the high diversity of MAH and MAA strains and the relative uniformity of MAS strains.

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.

Molecular identification of Mycobacterium avium subsp. silvaticum by duplex high-resolution melt analysis and subspecies-specific real-time PCR

Accurate identification of mycobacterial species and subspecies is essential to evaluate their significance and to perform epidemiological studies. The subspecies of Mycobacterium avium have different attributes but coincide in their zoonotic potential. Our knowledge about M. avium subsp. silvaticum is limited, since its identification is uncertain. Mycobacterium avium subsp. avium and M. avium subsp. silvaticum can be discriminated from each other based only on phenotypic characteristics, as they have almost identical genome sequences. Here we describe the development of a diagnostic method which enables the molecular identification of M. avium subsp. silvaticum and discrimination from M. avium subsp. avium based on genomic differences in a duplex high-resolution melt and M. avium subsp. silvaticum-specific mismatch real-time PCR. The developed assay was tested on reference strains and 199 field isolates, which were analyzed by phenotypic methods previously. This assay not only identified all 63 M. avium subsp. silvaticum and 138 M. avium subsp. avium strains correctly but also enabled the detection of mixed M. avium subsp. avium-M. avium subsp. silvaticum cultures. This is the first time that such a large panel of strains has been analyzed, and we also report the first isolation of M. avium subsp. silvaticum from red fox, red deer, wild boar, cattle, and badger. This assay is reliable, rapid, simple, inexpensive, and robust. It eliminates the long-existing problem of ambiguous phenotypic identification and opens up the possibility for detailed and comprehensive strain studies.

Comparative functional pan-genome analyses to build connections between genomic dynamics and phenotypic evolution in polycyclic aromatic hydrocarbon metabolism in the genus Mycobacterium

Background

The bacterial genus Mycobacterium is of great interest in the medical and biotechnological fields. Despite a flood of genome sequencing and functional genomics data, significant gaps in knowledge between genome and phenome seriously hinder efforts toward the treatment of mycobacterial diseases and practical biotechnological applications. In this study, we propose the use of systematic, comparative functional pan-genomic analysis to build connections between genomic dynamics and phenotypic evolution in polycyclic aromatic hydrocarbon (PAH) metabolism in the genus Mycobacterium.

Results

Phylogenetic, phenotypic, and genomic information for 27 completely genome-sequenced mycobacteria was systematically integrated to reconstruct a mycobacterial phenotype network (MPN) with a pan-genomic concept at a network level. In the MPN, mycobacterial phenotypes show typical scale-free relationships. PAH degradation is an isolated phenotype with the lowest connection degree, consistent with phylogenetic and environmental isolation of PAH degraders. A series of functional pan-genomic analyses provide conserved and unique types of genomic evidence for strong epistatic and pleiotropic impacts on evolutionary trajectories of the PAH-degrading phenotype. Under strong natural selection, the detailed gene gain/loss patterns from horizontal gene transfer (HGT)/deletion events hypothesize a plausible evolutionary path, an epistasis-based birth and pleiotropy-dependent death, for PAH metabolism in the genus Mycobacterium. This study generated a practical mycobacterial compendium of phenotypic and genomic changes, focusing on the PAH-degrading phenotype, with a pan-genomic perspective of the evolutionary events and the environmental challenges.

Conclusions

Our findings suggest that when selection acts on PAH metabolism, only a small fraction of possible trajectories is likely to be observed, owing mainly to a combination of the ambiguous phenotypic effects of PAHs and the corresponding pleiotropy- and epistasis-dependent evolutionary adaptation. Evolutionary constraints on the selection of trajectories, like those seen in PAH-degrading phenotypes, are likely to apply to the evolution of other phenotypes in the genus Mycobacterium.

Electronic supplementary material

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

Role of host- and pathogen-associated lipids in directing the immune response in mycobacterial infections, with emphasis on Mycobacterium avium subsp. paratuberculosis

Mycobacteria have a complex cell wall with a high lipid content that confers unique advantages for bacterial survival in the hostile host environment, leading to long-term infection. There is a wealth of evidence suggesting the role cell wall-associated lipid antigens play at the host-pathogen interface by contributing to bacterial virulence. One pathway that pathogenic mycobacteria use to subvert host immune pathways to their advantage is host cholesterol/lipid homeostasis. This review focuses on the possible role of pathogen- and host-associated lipids in the survival and persistence of pathogenic mycobacteria with emphasis on Mycobacterium avium subsp. paratuberculosis. We draw upon literature in diverse areas of infectious and metabolic diseases and explain a mechanism by which mycobacterial-induced changes in the host cellular energy state could account for phenomena that are a hallmark of chronic mycobacterial diseases.

Screening of Mycobacterium avium subsp. paratuberculosis mutants for attenuation in a bovine monocyte-derived macrophage model

Vaccination remains a major tool for prevention and progression of Johne's disease, a chronic enteritis of ruminants worldwide. Currently there is only one licensed vaccine within the United States and two vaccines licensed internationally against Johne's disease. All licensed vaccines reduce fecal shedding of Mycobacterium avium subsp. paratuberculosis (MAP) and delay disease progression. However, there are no available vaccines that prevent disease onset. A joint effort by the Johne's Disease Integrated Program (JDIP), a USDA-funded consortium, and USDA—APHIS/VS sought to identify transposon insertion mutant strains as vaccine candidates in part of a three phase study. The focus of the Phase I study was to evaluate MAP mutant attenuation in a well-defined in vitro bovine monocyte-derived macrophage (MDM) model. Attenuation was determined by colony forming unit (CFUs) counts and slope estimates. Based on CFU counts alone, the MDM model did not identify any mutant that significantly differed from the wild-type control, MAP K-10. Slope estimates using mixed models approach identified six mutants as being attenuated. These were enrolled in protection studies involving murine and baby goat vaccination-challenge models. MDM based approach identified trends in attenuation but this did not correlate with protection in a natural host model. These results suggest the need for alternative strategies for Johne's disease vaccine candidate screening and evaluation.

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.

Quantification of Mycobacterium avium subsp. paratuberculosis strains representing distinct genotypes and isolated from domestic and wildlife animal species by use of an automatic liquid culture system

Quantification of 11 clinical strains of Mycobacterium avium subsp. paratuberculosis isolated from domestic (cattle, sheep, and goat) and wildlife (fallow deer, deer, wild boar, and bison) animal species in an automatic liquid culture system (Bactec MGIT 960) was accomplished. The strains were previously isolated and typed using IS1311 PCR followed by restriction endonuclease analysis (PCR-REA) into type C, S, or B. A strain-specific quantification curve was generated for each M. avium subsp. paratuberculosis strain by relating the time to detection in the liquid culture system to the estimated log(10) CFU in each inoculum. According to their growth curves, the tested M. avium subsp. paratuberculosis strains were classified into two distinct groups. The first group included the S-type strain isolated from goat and all the sheep strains with C, S, and B genotypes. A second group contained the C- and B-type strains isolated from cattle, goat, and wildlife animals with the exception of the fallow deer strain. The strains isolated from cattle or sheep showed similar strain-specific standard curves irrespective of their genotype. In contrast, the strains isolated from goat or from wildlife animal species varied in their rates of growth in liquid culture. Universal-standard curves and algorithms for the quantification of each group of strains were generated. In addition, the liquid culture system was compared with a real-time quantitative PCR system for the quantification of the 11 M. avium subsp. paratuberculosis strains. Correlations between the estimated log(10) CFU and M. avium subsp. paratuberculosis DNA copy numbers were very high for all the tested strains (R ≥ 0.9).

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.

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.

Monoclonal Antibodies Bind A SNP-Sensitive Epitope that is Present Uniquely in Mycobacterium avium Subspecies Paratuberculosis

Due to a close genetic relatedness, there is no known antibody that detects Mycobacterium avium subspecies paratuberculosis (MAP), which causes Johne’s disease in cattle and sheep, and does not cross-react with other M. avium subspecies. In the present study, a monoclonal antibody (MAb; 17A12) was identified from mice immunized with a cell membrane fraction of MAP strain K-10. This antibody is 100% specific as it detected a 25-kDa protein in all 29 MAP whole cell lysates, but did not bind to any of the 29 non-paratuberculosis strains tested in immunoblot assays. However, the antibody revealed variable reactivity levels in MAP strains as it detected higher levels in bovine isolates but comparably lower levels in ovine isolates of MAP. In order to identify the target binding protein for 17A12, a lambda phage expression library of MAP genomic fragments was screened with the MAb. Four reactive clones were identified, sequenced and all shown to be overlapping. Further analysis revealed all four clones expressed an unknown protein encoded by a sequence that is not annotated in the K-10 genome and overlapped with MAP3422c on the opposing DNA strand. The epitope of 17A12 was precisely defined to seven amino acids and was used to query the K-10 genome. Similarity searches revealed another protein, encoded by MAP1025, possessed a similar epitope (one-amino acid mismatch) that also reacted strongly to the antibody. A single nucleotide polymorphism (SNP) in MAP1025 was then identified by comparative sequence analysis, which results in a Pro28His change at residue 28, the first amino acid within the 17A12 epitope. This SNP is present in all MAP strains but absent in all non-MAP strains and accounts for the specificity of the 17A12 antibody. This new antibody is the first ever isolated that binds only to the paratuberculosis subspecies of M. avium and opens new possibilities for the specific detection of this significant ruminant pathogen.

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.

Experimental infection of a bovine model with human isolates of Mycobacterium avium subsp. paratuberculosis

Mycobacterium avium subsp. paratuberculosis (Map), the etiologic agent of Johne's disease (JD) in ruminants, has been implicated in the pathogenesis of Crohn's disease (CD) in humans. We developed a bovine ileal cannulation model to facilitate comparison of the immune response to Map and the mechanisms of pathogenesis in cattle and humans. Initial studies showed a T cannula could be maintained for up to a year in calves without inducing inflammation or adversely affecting intestinal function. Map introduced through the cannula established a persistent low level of infection without inflammation. Infection elicited an immune response to Map antigens detectable by flow cytometry. Further studies now show the cannulation model can be used with cows during the later stage of infection, affording access to the target tissue at all stages of infection. The studies also revealed no difference in infectivity or immunogenicity of isolates of Map obtained from cattle or humans with CD. Comparison of the immune response to Map during the early and late stages of infection using PCR, flow cytometry and QRT-PCR, showed the immune response early in the disease process is dominated by CD4 T cells. A CD8 response is delayed but comparable at later stages of infection. Genes for pro-inflammatory cytokines IFN-γ and the recently identified genes encoding IL-17 and IL-22 are up regulated in infected animals. These findings reveal that both human and bovine isolates of Map can establish infection and induce similar immune responses in a bovine model. They also reveal the cytokine responses elicited in cattle are similar to those implicated in CD pathogenesis.

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.

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.

Iron-sparing response of Mycobacterium avium subsp. paratuberculosis is strain dependent

Background

Two genotypically and microbiologically distinct strains of Mycobacterium avium subsp. paratuberculosis (MAP) exist - S and C MAP strains that primarily infect sheep and cattle, respectively. Concentration of iron in the cultivation medium has been suggested as one contributing factor for the observed microbiologic differences. We recently demonstrated that S strains have defective iron storage systems, leading us to propose that these strains might experience iron toxicity when excess iron is provided in the medium. To test this hypothesis, we carried out transcriptional and proteomic profiling of these MAP strains under iron-replete or -deplete conditions.

Results

We first complemented M. smegmatisΔideR with IdeR of C MAP or that derived from S MAP and compared their transcription profiles using M. smegmatis mc²155 microarrays. In the presence of iron, sIdeR repressed expression of bfrA and MAP2073c, a ferritin domain containing protein suggesting that transcriptional control of iron storage may be defective in S strain. We next performed transcriptional and proteomic profiling of the two strain types of MAP under iron-deplete and -replete conditions. Under iron-replete conditions, C strain upregulated iron storage (BfrA), virulence associated (Esx-5 and antigen85 complex), and ribosomal proteins. In striking contrast, S strain downregulated these proteins under iron-replete conditions. iTRAQ (isobaric tag for relative and absolute quantitation) based protein quantitation resulted in the identification of four unannotated proteins. Two of these were upregulated by a C MAP strain in response to iron supplementation. The iron-sparing response to iron limitation was unique to the C strain as evidenced by repression of non-essential iron utilization enzymes (aconitase and succinate dehydrogenase) and upregulation of proteins of essential function (iron transport, [Fe-S] cluster biogenesis and cell division).

Conclusions

Taken together, our study revealed that C and S strains of MAP utilize divergent metabolic pathways to accommodate in vitro iron stress. The knowledge of the metabolic pathways these divergent responses play a role in are important to 1) advance our ability to culture the two different strains of MAP efficiently, 2) aid in diagnosis and control of Johne's disease, and 3) advance our understanding of MAP virulence.

Primary transcriptomes of Mycobacterium avium subsp. paratuberculosis reveal proprietary pathways in tissue and macrophages

Background

Mycobacterium avium subsp. paratuberculosis (MAP) persistently infects intestines and mesenteric lymph nodes leading to a prolonged subclinical disease. The MAP genome sequence was published in 2005, yet its transcriptional organization in natural infection is unknown. While prior research analyzed regulated gene sets utilizing defined, in vitro stress related or advanced surgical methods with various animal species, we investigated the intracellular lifestyle of MAP in the intestines and lymph nodes to understand the MAP pathways that function to govern this persistence.

Results

Our transcriptional analysis shows that 21%, 8% and 3% of the entire MAP genome was represented either inside tissues, macrophages or both, respectively. Transcripts belonging to latency and cell envelope biogenesis were upregulated in the intestinal tissues whereas those belonging to intracellular trafficking and secretion were upregulated inside the macrophages. Transcriptomes of natural infection and in vitro macrophage infection shared genes involved in transcription and inorganic ion transport and metabolism. MAP specific genes within large sequence polymorphisms of ancestral M. avium complex were downregulated exclusively in natural infection.

Conclusions

We have unveiled common and unique MAP pathways associated with persistence, cell wall biogenesis and virulence in naturally infected cow intestines, lymph nodes and in vitro infected macrophages. This dichotomy also suggests that in vitro macrophage models may be insufficient in providing accurate information on the events that transpire during natural infection. This is the first report to examine the primary transcriptome of MAP at the local infection site (i.e. intestinal tissue). Regulatory pathways that govern the lifecycle of MAP appear to be specified by tissue and cell type. While tissues show a "shut-down" of major MAP metabolic genes, infected macrophages upregulate several MAP specific genes along with a putative pathogenicity island responsible for iron acquisition. Many of these regulatory pathways rely on the advanced interplay of host and pathogen and in order to decipher their message, an interactome must be established using a systems biology approach. Identified MAP pathways place current research into direct alignment in meeting the future challenge of creating a MAP-host interactome.

Characterization of a novel Lactobacillus species closely related to Lactobacillus johnsonii using a combination of molecular and comparative genomics methods

BACKGROUND

Comparative genomic hybridization (CGH) constitutes a powerful tool for identification and characterization of bacterial strains. In this study we have applied this technique for the characterization of a number of Lactobacillus strains isolated from the intestinal content of rats fed with a diet supplemented with sorbitol.

RESULTS

Phylogenetic analysis based on 16S rRNA gene, recA, pheS, pyrG and tuf sequences identified five bacterial strains isolated from the intestinal content of rats as belonging to the recently described Lactobacillus taiwanensis species. DNA-DNA hybridization experiments confirmed that these five strains are distinct but closely related to Lactobacillus johnsonii and Lactobacillus gasseri. A whole genome DNA microarray designed for the probiotic L. johnsonii strain NCC533 was used for CGH analysis of L. johnsonii ATCC 33200T, L. johnsonii BL261, L. gasseri ATCC 33323T and L. taiwanensis BL263. In these experiments, the fluorescence ratio distributions obtained with L. taiwanensis and L. gasseri showed characteristic inter-species profiles. The percentage of conserved L. johnsonii NCC533 genes was about 83% in the L. johnsonii strains comparisons and decreased to 51% and 47% for L. taiwanensis and L. gasseri, respectively. These results confirmed the separate status of L. taiwanensis from L. johnsonii at the level of species, and also that L. taiwanensis is closer to L. johnsonii than L. gasseri is to L. johnsonii.

CONCLUSION

Conventional taxonomic analyses and microarray-based CGH analysis have been used for the identification and characterization of the newly species L. taiwanensis. The microarray-based CGH technology has been shown as a remarkable tool for the identification and fine discrimination between phylogenetically close species, and additionally provided insight into the adaptation of the strain L. taiwanensis BL263 to its ecological niche.

Broad conservation of milk utilization genes in Bifidobacterium longum subsp. infantis as revealed by comparative genomic hybridization

Human milk oligosaccharides (HMOs) are the third-largest solid component of milk. Their structural complexity renders them nondigestible to the host but liable to hydrolytic enzymes of the infant colonic microbiota. Bifidobacteria and, frequently, Bifidobacterium longum strains predominate the colonic microbiota of exclusively breast-fed infants. Among the three recognized subspecies of B. longum, B. longum subsp. infantis achieves high levels of cell growth on HMOs and is associated with early colonization of the infant gut. The B. longum subsp. infantis ATCC 15697 genome features five distinct gene clusters with the predicted capacity to bind, cleave, and import milk oligosaccharides. Comparative genomic hybridizations (CGHs) were used to associate genotypic biomarkers among 15 B. longum strains exhibiting various HMO utilization phenotypes and host associations. Multilocus sequence typing provided taxonomic subspecies designations and grouped the strains between B. longum subsp. infantis and B. longum subsp. longum. CGH analysis determined that HMO utilization gene regions are exclusively conserved across all B. longum subsp. infantis strains capable of growth on HMOs and have diverged in B. longum subsp. longum strains that cannot grow on HMOs. These regions contain fucosidases, sialidases, glycosyl hydrolases, ABC transporters, and family 1 solute binding proteins and are likely needed for efficient metabolism of HMOs. Urea metabolism genes and their activity were exclusively conserved in B. longum subsp. infantis. These results imply that the B. longum has at least two distinct subspecies: B. longum subsp. infantis, specialized to utilize milk carbon, and B. longum subsp. longum, specialized for plant-derived carbon metabolism.

Use of signal quality measurements to gain efficiency in the analysis of cDNA microarray data

This research provides a new way to measure error in microarray data in order to improve gene expression analysis. Microarray data contains many sources of error. In order to glean information about mRNA expression levels, the true signal must first be segregated from noise. This research focuses on the variation that can be captured at the spot level in cDNA microarray images. Variation at other levels, due to differences at the array, dye, and block levels, can be corrected for by a variety of existing normalization procedures. Two signal quality estimates that capture the reliability of each spot printed on a microarray are described. A parametric estimate of within-spot variance, referred to here as sigma(2)(spot), assumes that pixels follow a normal distribution and are spatially correlated. A non-parametric estimate of error, called the mean square prediction error (MSPE), assumes that spots of high quality possess pixels that are similar to their neighbors. This paper will provide a framework to use either spot quality measure in downstream analysis, specifically as weights in regression models. Using these spot quality estimates as weights can result in greater efficiency, in a statistical sense, when modeling microarray data.

Lack of insertional-deletional polymorphism in a collection of Mycobacterium ulcerans isolates from Ghanaian Buruli ulcer patients

Mycobacterium ulcerans causes the devastating infectious skin disease Buruli ulcer and has a monomorphic population structure. The resolution of conventional genetic fingerprinting methods is therefore not sufficient for microepidemiological studies aiming to characterize transmission pathways. In a previous comparative genomic hybridization analysis with a microarray covering part of the M. ulcerans genome, we have found extensive insertional-deletional sequence polymorphisms among M. ulcerans isolates of diverse geographic origins that allowed us to distinguish between strains coming from different continents. Since large numbers of insertion sequences are spread over the genome of African M. ulcerans strains, we reasoned that these may drive large sequence polymorphisms in otherwise clonal local mycobacterial populations. In this study, we used a printed DNA microarray covering the whole genome of the Ghanaian M. ulcerans reference strain Agy99 for comparative genomic hybridization. The assay identified multiple regions of difference when DNA of a Japanese M. ulcerans strain was analyzed. In contrast, not a single insertional-deletional genomic variation was found within a panel of disease isolates coming from an area of Ghana where Buruli ulcer is endemic. These results indicate that, despite the expectations deduced from other mycobacterial pathogens, only analyses of single nucleotide polymorphisms will have the potential to differentiate local populations of M. ulcerans.

Identification and functional characterization of the iron-dependent regulator (IdeR) of Mycobacterium avium subsp. paratuberculosis

Mycobacterium avium subspecies paratuberculosis (MAP), the causative agent of Johne's disease in cattle and sheep, has unique iron requirements in that it is mycobactin-dependent for cultivation in vitro. The iron-dependent regulator (IdeR) is a well-characterized global regulator responsible for maintaining iron homeostasis in Mycobacterium tuberculosis (MTB). We identified an orthologous segment in the MAP genome, MAP2827, with >93 % amino acid identity to MTB IdeR. Electrophoretic mobility shift assays and DNase protection assays confirmed that MAP2827 binds the 19 bp consensus motif (iron box) on the MAP genome. Sequencing of MAP2827 from multiple isolates revealed a non-synonymous change (R91G) exclusive to sheep strains. Reporter gene assays and quantitative real-time RT-PCR assays in two diverse MAP strains and in an ideR deletion mutant of M. smegmatis (mc(2)155) suggested that both sheep MAP IdeR (sIdeR) and cattle MAP IdeR (cIdeR) repress mbtB transcription at high iron concentrations and relieve repression at low iron concentrations. On the other hand, bfrA (an iron storage gene) was upregulated by cIdeR when presented with MTB or the cattle MAP bfrA promoter, and was downregulated by sIdeR in the presence of MTB, or sheep or cattle MAP bfrA promoters, at high iron concentrations. The differential iron regulatory mechanisms between IdeR-regulated genes across strains may contribute to the differential growth or pathogenic characteristics of sheep and cattle MAP strains. Taken together, our study provides a possible reason for mycobactin dependency and suggests strong implications in the differential iron acquisition and storage mechanisms in MAP.

Optical mapping of the Mycobacterium avium subspecies paratuberculosis genome

BACKGROUND

Infection of cattle with Mycobacterium avium subspecies paratuberculosis (M. ap) causes severe economic losses to the dairy industry in the USA and worldwide. In an effort to better examine diversity among M. ap strains, we used optical mapping to profile genomic variations between strains of M. ap K-10 (sequenced strain) and M. ap ATCC 19698 (type strain).

RESULTS

The assembled physical restriction map of M. ap ATCC 19698 showed a genome size of 4,839 kb compared to the sequenced K-10 genome of 4,830 kb. Interestingly, alignment of the optical map of the M. ap ATCC 19698 genome to the complete M. ap K-10 genome sequence revealed a 648-kb inversion around the origin of replication. However, Southern blotting, PCR amplification and sequencing analyses of the inverted region revealed that the genome of M. ap K-10 differs from the published sequence in the region starting from 4,197,080 bp to 11,150 bp, spanning the origin of replication. Additionally, two new copies of the coding sequences > 99.8% were identified, identical to the MAP0849c and MAP0850c genes located immediately downstream of the MAP3758c gene.

CONCLUSION

The optical map of M. ap ATCC 19698 clearly indicated the miss-assembly of the sequenced genome of M. ap K-10. Moreover, it identified 2 new genes in M. ap K-10 genome. This analysis strongly advocates for the utility of physical mapping protocols to complement genome sequencing projects.

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

Mycobacterium avium subsp. paratuberculosis is an important animal pathogen widely disseminated in the environment that has also been associated with Crohn's disease in humans. Three M. avium subsp. paratuberculosis genomotypes are recognized, but genomic differences have not been fully described. To further investigate these potential differences, a 60-mer oligonucleotide microarray (designated the MAPAC array), based on the combined genomes of M. avium subsp. paratuberculosis (strain K-10) and Mycobacterium avium subsp. hominissuis (strain 104), was designed and validated. By use of a test panel of defined M. avium subsp. paratuberculosis strains, the MAPAC array was able to identify a set of large sequence polymorphisms (LSPs) diagnostic for each of the three major M. avium subsp. paratuberculosis types. M. avium subsp. paratuberculosis type II strains contained a smaller genomic complement than M. avium subsp. paratuberculosis type I and M. avium subsp. paratuberculosis type III genomotypes, which included a set of genomic regions also found in M. avium subsp. hominissuis 104. Specific PCRs for genes within LSPs that differentiated M. avium subsp. paratuberculosis types were devised and shown to accurately screen a panel (n = 78) of M. avium subsp. paratuberculosis strains. Analysis of insertion/deletion region INDEL12 showed deletion events causing a reduction in the complement of mycobacterial cell entry genes in M. avium subsp. paratuberculosis type II strains and significantly altering the coding of a major immunologic protein (MPT64) associated with persistence and granuloma formation. Analysis of MAPAC data also identified signal variations in several genomic regions, termed variable genomic islands (vGIs), suggestive of transient duplication/deletion events. vGIs contained significantly low GC% and were immediately flanked by insertion sequences, integrases, or short inverted repeat sequences. Quantitative PCR demonstrated that variation in vGI signals could be associated with colony growth rate and morphology.

Insertion and deletion events that define the pathogen Mycobacterium avium subsp. paratuberculosis

Mycobacterium avium comprises genetically related yet phenotypically distinct subspecies. Consistent with their common origin, whole-genome sequence comparisons have revealed extensive synteny among M. avium organisms. However, the sequenced strains also display numerous regions of heterogeneity that likely contribute to the diversity of the individual subspecies. Starting from a phylogenetic framework derived by multilocus sequence analysis, we examined the distribution of 25 large sequence polymorphisms across a panel of genetically defined M. avium strains. This distribution was most variable among M. avium subsp. hominissuis isolates. In contrast, M. avium subsp. paratuberculosis strains exhibited a characteristic profile, with all isolates containing a set of genomic insertions absent from other M. avium strains. The emergence of the pathogen from its putative M. avium subsp. hominissuis ancestor entailed the acquisition of approximately 125 kb of novel genetic material, followed by a second phase, characterized by reductive genomics. One genomic deletion is common to all isolates while additional deletions distinguish two major lineages of M. avium subsp. paratuberculosis. For the average strain, these losses total at least 38 kb (sheep lineage) to 90 kb (cattle lineage). This biphasic pattern of evolution, characterized by chromosomal gene acquisition with subsequent gene loss, describes the emergence of M. avium subsp. paratuberculosis and may serve as a general model for the origin of pathogenic mycobacteria.