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

Characterisation of IS1311 in Mycobacterium avium subspecies paratuberculosis genomes: Typing, continental clustering, microbial evolution and host adaptation

Johne's disease (JD), caused by Mycobacterium avium subspecies paratuberculosis (MAP) is a global burden for livestock producers and has an association with Crohn's disease in humans. Within MAP there are two major lineages, S/Type I/TypeIII and C/Type II, that vary in phenotype including culturability, host preference and virulence. These lineages have been identified using the IS1311 element, which contains a conserved, single nucleotide polymorphism. IS1311 and the closely related IS1245 element belong to the IS256 family of insertion sequences, are dispersed throughout M. avium taxa but remain poorly characterised. To investigate the distribution and diversity of IS1311 in MAP, 805 MAP genomes were collated from public databases. IS1245 was absent, while IS1311 sequence, copy number and insertion loci were conserved between MAP S lineages and varied within the MAP C lineage. One locus was specific to the S strains, which contained nine IS1311 copies. In contrast, C strains contained either seven or eight IS1311 loci. Most insertion loci were associated with the boundaries of homologous regions that had undergone genome rearrangement between the MAP lineages, suggesting that this sequence may be a driver of recombination. Phylogenomic geographic clustering of MAP subtypes was demonstrated for the first time, at continental scale, and indicated that there may have been recent MAP transmission between Europe and North America, in contrast to Australia where importation of live ruminants is generally prohibited. This investigation confirmed the utility of IS1311 typing in epidemiological studies and resolved anomalies in past studies. The results shed light on potential mechanisms of niche/host adaptation, virulence of MAP and global transmission dynamics.

Molecular Diversity of Mycobacterium avium subsp. paratuberculosis in Four Dairy Goat Herds from Thuringia (Germany)

This study investigated the intra- and inter-herd diversity of Mycobacterium avium subsp. paratuberculosis (MAP) isolates from four goat herds in Thuringia (Germany) that were affected by paratuberculosis for several years. The main focus was on the characterization and distribution of genotypes among animals and the environment of goat herd 1. This study included 196 isolates from the feces of 121 infected goats, various tissues from 13 clinically diseased goats, 29 environmental samples from herd 1, and additionally, 22 isolates of different origin from herds 2 to 4. The isolates, sampled between 2018 and 2022, were genotyped using short-sequence-repeat (SSR) analysis, mycobacterial-interspersed repetitive units-variable-number tandem repeat (MIRU-VNTR) analysis, and a single nucleotide polymorphism (SNP)-based assay for phylogenetic grouping. All the isolates belonged to the MAP-C group. In herd 1, one predominant genotype was determined, while two other genotypes were identified very rarely and only in fecal and environmental samples. One of three further genotypes was found in each of herds 2 to 4. The assignment of genotypes to different phylogenetic clades suggested six different infection strains. The results indicated no epidemiological links between the examined herds. Based on the current MAP genotyping data from Germany, possible sources of infection are MAP-contaminated barns previously used by infected cattle and the purchase of sub-clinically infected goats.

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.

More insights about genomic population structure of Mycobacterium avium subspecies paratuberculosis (Map) from multiple hosts in west and central provinces of Iran using a boosted genotyping approach

To investigate the population genetic of Mycobacterium avium subsp. paratuberculosis (Map) in Iran, Mycobacterial Interspersed Repetitive Units (MIRUs) and Multi Locus Short Sequence Repeat (MLSSR) system were employed. Numerous genotypes by MIRU (N = 11) and MLSSR (N = 9) methods bearing discriminatory indices of 0.90 and 0.79 respectively, were obtained. Browsing the INRA-Nouzilly list (http://mac-inmv.tours.inra.fr/) detected 3 of the found patterns as new types. Some loci either MIRU-VNTR or SSR proved more polymorphic and therefore are recommended to be applied in priority for strain typing in the Iranian environment. While identical MIRU-VNTR or MLSSR patterns were detected among different conspecifics and geographical locations, dissimilar types were also observed at the same farms an indication of coexistence of Map strains within one herd. We suggest extension of the genotyping work described here to include more endogenous isolates in order to better analysis of transmission and virulence in epidemiology and control of paratuberculosis.

Can Genome Sequencing Coupled to Flux Balance Analyses Offer Precision Guidance for Industrial Strain Development? The Lessons from Carbon Trafficking in Corynebacterium glutamicum ATCC 21573

Systems biology tools offer new prospects for industrial strain selection. For bacteria that are significant for industrial applications, whole-genome sequencing coupled to flux balance analysis (FBA) can help unpack the complex relationships between genome mutations and carbon trafficking. This work investigates the l-tyrosine (l-Tyr) overproducing model system Corynebacterium glutamicum ATCC 21573 with an eye to more rational and precision strain development. Using genome-wide mutational analysis of C. glutamicum, we identified 27,611 single nucleotide polymorphisms and 479 insertion/deletion mutations. Mutations in the carbon uptake machinery have led to phosphotransferase system-independent routes as corroborated with FBA. Mutations within the central carbon metabolism of C. glutamicum impaired the carbon flux, as evidenced by the lower growth rate. The entry to and flow through the tricarboxylic acid cycle was affected by mutations in pyruvate and α-ketoglutarate dehydrogenase complexes, citrate synthase, and isocitrate dehydrogenase. FBA indicated that the estimated flux through the shikimate pathway became larger as the l-Tyr production rate increased. In addition, protocatechuate export was probabilistically impossible, which could have contributed to the l-Tyr accumulation. Interestingly, aroG and cg0975, which have received previous attention for aromatic amino acid overproduction, were not mutated. From the branch point molecule, prephenate, the change in the promoter region of pheA could be an influential contributor. In summary, we suggest that genome sequencing coupled with FBA is well poised to offer rational guidance for industrial strain development, as evidenced by these findings on carbon trafficking in C. glutamicum ATCC 21573.

Improved DNA Amplification of the Hallmark IS900 Element in Mycobacterium avium subsp. paratuberculosis: a Reexamination Based on Whole-Genome Sequence Analysis

Amplification of the IS900 multicopy element is a hallmark nucleic acid-based diagnostic test for Mycobacterium avium subsp. paratuberculosis, which causes Johne's disease in ruminants. This assay is frequently used to determine the presence of the bacterium in feces of infected cattle and sheep. Two IS900 primer sets developed in the 1990s were widely used for decades, and their use has continued in current studies. However, these primers were developed prior to the availability of complete genome sequences. Recent sequence analysis of the binding locations for one primer pair (P90/P91) identified errors and binding inefficiencies that can be easily corrected to further increase detection sensitivity. The P90 primer is missing two nucleotides that should be present near the 3' end, and it does not bind all copies of IS900 due to 5' deletions at some IS900 loci. These IS900 primer pairs, along with newly developed primers, were tested by real-time PCR on purified genomic DNA to determine which primer set performed the best and how primer design errors affect amplification efficiencies. The newly designed PCR primer set (JB5) showed increased sensitivity by two to three quantification cycles using purified genomic DNA and was similar in efficiency to 150C/921. These tests were extended using DNA from feces and tissues of infected cows, which showed similar results. Finally, a 167-bp partial duplication of IS900 was found in type I strains. Although P90 and P91 primers successfully amplify M. avium subsp. paratuberculosis DNA, their use should be discontinued in favor of more efficient primer pairs in future studies. IMPORTANCE This study is an example of how applied genomic analysis can aid diagnostic test improvements. Detection of Mycobacterium avium subsp. paratuberculosis infection of livestock prior to the appearance of clinical disease signs is very difficult but essential for identifying animals shedding the bacterium to prevent transmission of Johne's disease. Total M. avium subsp. paratuberculosis quantity in the feces as determined by real-time PCR (qPCR) using the IS900 target indicates bacterial shedding status and potential for transmission of the pathogen. However, legacy primers designed prior to the availability of complete genome sequences that are used in these tests to detect M. avium subsp. paratuberculosis were based on data from only a single copy of IS900 and not considering all copies collectively as a group. This approach resulted in primer design errors which can be easily corrected to improve test sensitivities. We tested original primers that contain these errors and their corrected versions by qPCR and showed improved sensitivity on purified genomic DNA as well as fecal and tissue samples. These findings may help detect the organism from environmental samples on farms where sensitivity is currently lacking.

Strain diversity in Mycobacterium avium subsp. paratuberculosis-positive bovine fecal samples collected in Switzerland

Paratuberculosis or Johne's disease is a chronic intestinal disease in domestic and wild ruminants. It affects global dairy economy and is caused by Mycobacterium avium subsp. paratuberculosis (MAP). The objective of this study was to analyze strain diversity in MAP-positive fecal samples by using a particular single nucleotide polymorphism (SNP) distinguishing between cattle (C-) and sheep (S-) type MAP and analysis of SNPs within gyrA and gyrB genes differentiating between Types I, II, and III. Moreover, mycobacterial interspersed repetitive unit and variable-number tandem repeat (MIRU-VNTR) analysis using eight established loci was performed. A total of 90 fecal samples from diseased animals presenting diarrhea and/or weight loss, originating from 59 bovine herds across 16 cantons of Switzerland were screened by PCR for the MAP-specific F57 and IS900 genes and were further subtyped. 96.7% and 3.3% of the samples contained C- and S-type MAP, respectively. Ten INRA Nouzilly MIRU-VNTR (INMV) profiles, with a discriminatory index of 0.802, calculated based on 65 epidemiological independent genotypes, were detected: INMV 1 (33.8%), INMV 2 (23.1%), INMV 6 (16.9%), INMV 9 (9.2%), INMV 116 (4.6%), INMV 3 (3.1%), INMV 5 (3.1%) and INMV 72 (1.5%), including two novel INMV profiles, namely INMV 253 (3.1%; S-type III) and INMV 252 (1.5%; C-type). INMV 1, INMV 2, and INMV 6 comprised almost 75% of the F57- and IS900-positive samples. Typing data from 11 herds suggest that there are some herds with intra-herd diversity of genotypes. The results of this study indicate a heterogeneity of MAP in Switzerland.

Agronomic efficiency and genome mining analysis of the wheat-biostimulant rhizospheric bacterium Pseudomonas pergaminensis sp. nov. strain 1008T

Pseudomonas sp. strain 1008 was isolated from the rhizosphere of field grown wheat plants at the tillering stage in an agricultural plot near Pergamino city, Argentina. Based on its in vitro phosphate solubilizing capacity and the production of IAA, strain 1008 was formulated as an inoculant for bacterization of wheat seeds and subjected to multiple field assays within the period 2010-2017. Pseudomonas sp. strain 1008 showed a robust positive impact on the grain yield (+8% on average) across a number of campaigns, soil properties, seed genotypes, and with no significant influence of the simultaneous seed treatment with a fungicide, strongly supporting the use of this biostimulant bacterium as an agricultural input for promoting the yield of wheat. Full genome sequencing revealed that strain 1008 has the capacity to access a number of sources of inorganic and organic phosphorus, to compete for iron scavenging, to produce auxin, 2,3-butanediol and acetoin, and to metabolize GABA. Additionally, the genome of strain 1008 harbors several loci related to rhizosphere competitiveness, but it is devoid of biosynthetic gene clusters for production of typical secondary metabolites of biocontrol representatives of the Pseudomonas genus. Finally, the phylogenomic, phenotypic, and chemotaxonomic comparative analysis of strain 1008 with related taxa strongly suggests that this wheat rhizospheric biostimulant isolate is a representative of a novel species within the genus Pseudomonas, for which the name Pseudomonas pergaminensis sp. nov. (type strain 1008T = DSM 113453T = ATCC TSD-287T) is proposed.

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.

Activity-Based Protein Profiling for the Identification of Novel Carbohydrate-Active Enzymes Involved in Xylan Degradation in the Hyperthermophilic Euryarchaeon Thermococcus sp. Strain 2319x1E

Activity-based protein profiling (ABPP) has so far scarcely been applied in Archaea in general and, especially, in extremophilic organisms. We herein isolated a novel Thermococcus strain designated sp. strain 2319x1E derived from the same enrichment culture as the recently reported Thermococcus sp. strain 2319x1. Both strains are able to grow with xylan as the sole carbon and energy source, and for Thermococcus sp. strain 2319x1E (optimal growth at 85°C, pH 6–7), the induction of xylanolytic activity in the presence of xylan was demonstrated. Since the solely sequence-based identification of xylanolytic enzymes is hardly possible, we established a complementary approach by conducting comparative full proteome analysis in combination with ABPP using α- or β-glycosidase selective probes and subsequent mass spectrometry (MS)-based analysis. This complementary proteomics approach in combination with recombinant protein expression and classical enzyme characterization enabled the identification of a novel bifunctional maltose-forming α-amylase and deacetylase (EGDIFPOO_00674) belonging to the GH57 family and a promiscuous β-glycosidase (EGIDFPOO_00532) with β-xylosidase activity. We thereby further substantiated the general applicability of ABPP in archaea and expanded the ABPP repertoire for the identification of glycoside hydrolases in hyperthermophiles.

Serological investigation and genotyping of Mycobacterium avium subsp. paratuberculosis in sheep and goats in Inner Mongolia, China

Paratuberculosis a contagious and chronic disease in domestic and wild ruminants, is caused by Mycobacterium avium subspecies paratuberculosis (MAP). Typical clinical signs include intractable diarrhea, progressive emaciation, proliferative enteropathy, and mesenteric lymphadenitis. Paratuberculosis is endemic to many parts of the world and responsible for considerable economic losses. In this study, different types of paratuberculosis and MAP in sheep and goats were investigated in Inner Mongolia, a northern province in China contiguous with two countries and eight other provinces. A total of 4434 serum samples were collected from six cities in the western, central, and eastern regions of Inner Mongolia and analyzed using the ELISA test. In addition, tissue samples were collected from seven animals that were suspected to be infected with MAP. Finally, these tissues samples were analyzed by histopathological examination followed by polymerase chain reaction (PCR), IS1311 PCR-restriction enzyme analysis (PCR-REA), and a sequence analysis of five genes. Among all 4434 ruminant serum samples collected from the six cities in the western, central, and eastern regions of Inner Mongolia, 7.60% (337/4434) measured positive for the MAP antibody. The proportions of positive MAP antibody results for serum samples collected in the western, central, and eastern regions were 5.10% (105/2058), 6.63% (85/1282), and 13.44% (147/1094), respectively. For the seven suspected infected animals selected from the herd with the highest rate of positivity, the gross pathology and histopathology of the necropsied animals were found to be consistent with the pathological features of paratuberculosis. The PCR analysis further confirmed the diagnosis of paratuberculosis. The rest of the results demonstrated that herds of sheep and goats in Inner Mongolia were infected with both MAP type II and type III. To the best of our knowledge, this is the first study of the two subtypes of MAP strains in sheep and goats in Inner Mongolia.

Whole-Genome Analysis of Mycobacterium avium subsp. paratuberculosis IS900 -

Mycobacterium avium subsp. paratuberculosis (Map) is the etiological agent of Johne’s disease in ruminants. The IS900 insertion sequence (IS) has been used widely as an epidemiological marker and target for PCR diagnosis. Updated DNA sequencing technologies have led to a rapid increase in available Map genomes, which makes it possible to analyze the distribution of IS900 in this slow-growing bacterium. The objective of this study is to characterize the distribution of the IS900 element and how it affects genomic evolution and gene function of Map. A secondary goal is to develop automated in silico restriction fragment length polymorphism (RFLP) analysis using IS900. Complete genomes from the major phylogenetic lineages known as C-type and S-type (including subtypes I and III), were chosen to represent the genetic diversity of Map. IS900 elements were located in these genomes using BLAST software and the relevant fragments extracted. An in silico RFLP analysis using the BstEII restriction site was performed to obtain exact sizes of the DNA fragments carrying a copy of IS900 and the resulting RFLP profiles were analyzed and compared by digital visualization of the separated restriction fragments. The program developed for this study allowed automated localization of IS900 sequences to identify their position within each genome along with the exact number of copies per genome. The number of IS900 copies ranged from 16 in the C-type isolate to 22 in the S-type subtype I isolate. A loci-by-loci sequence alignment of all IS900 copies within the three genomes revealed new sequence polymorphisms that define three sequevars distinguishing the subtypes. Nine IS900 insertion site locations were conserved across all genomes studied while smaller subsets were unique to a particular lineage. Preferential insertion motif sequences were identified for IS900 along with genes bordering all IS900 insertions. Rarely did IS900 insert within coding sequences as only three genes were disrupted in this way. This study makes it possible to automate IS900 distribution in Map genomes to enrich knowledge on the distribution dynamics of this IS for epidemiological purposes, for understanding Map evolution and for studying the biological implications of IS900 insertions.

Complete Genome Sequence of a Type III Ovine Strain of Mycobacterium avium subsp. paratuberculosis

The complete genome sequence of a type III strain of Mycobacterium avium subsp. paratuberculosis was determined. The genome size for this pathogen of sheep is 4,895,755 bp with no plasmid DNA. The chromosome contains 19 copies of the hallmark IS900 element, which is routinely used to identify this subspecies.

The complete genome sequence of a type III strain of Mycobacterium avium subsp. paratuberculosis was determined. The genome size for this pathogen of sheep is 4,895,755 bp with no plasmid DNA. The chromosome contains 19 copies of the hallmark IS900 element, which is routinely used to identify this subspecies.