Genomic and phenotypic comparison of two Salmonella Typhimurium strains responsible for consecutive salmonellosis outbreaks in New Zealand.

Salmonella enterica serovar Typhimurium DT160 was the predominant cause of notified human salmonellosis cases in New Zealand from 2000 to 2010, before it was superseded by another S. Typhimurium strain, DT56 variant (DT56v). Whole genome sequencing and phenotypic testing were used to compare 109 DT160 isolates with eight DT56v isolates from New Zealand animal and human sources. Phylogenetic analysis provided evidence that DT160 and DT56v strains were distantly related with an estimated date of common ancestor between 1769 and 1821. The strains replicated at different rates but had similar antimicrobial susceptibility profiles. Both strains were resistant to the phage expressed from the chromosome of the other strain, which may have contributed to the emergence of DT56v. DT160 contained the pSLT virulence plasmid, and the sseJ and sseK2 genes that may have contributed to the higher reported prevalence compared to DT56v. A linear pBSSB1-family plasmid was also found in one of the DT56v isolates, but there was no evidence that this plasmid affected bacterial replication or antimicrobial susceptibility. One of the DT56v isolates was also sequenced using long-read technology and found to contain an uncommon chromosome arrangement for a Typhimurium isolate. This study demonstrates how comparative genomics and phenotypic testing can help identify strain-specific elements and factors that may have influenced the emergence and supersession of bacterial strains of public health importance.

Using whole genome sequencing to investigate transmission in a multi-host system: bovine tuberculosis in New Zealand.

BACKGROUND: Bovine tuberculosis (bTB), caused by Mycobacterium bovis, is an important livestock disease raising public health and economic concerns around the world. In New Zealand, a number of wildlife species are implicated in the spread and persistence of bTB in cattle populations, most notably the brushtail possum (Trichosurus vulpecula). Whole Genome Sequenced (WGS) M. bovis isolates sourced from infected cattle and wildlife across New Zealand were analysed. Bayesian phylogenetic analyses were conducted to estimate the substitution rate of the sampled population and investigate the role of wildlife. In addition, the utility of WGS was examined with a view to these methods being incorporated into routine bTB surveillance. RESULTS: A high rate of exchange was evident between the sampled wildlife and cattle populations but directional estimates of inter-species transmission were sensitive to the sampling strategy employed. A relatively high substitution rate was estimated, this, in combination with a strong spatial signature and a good agreement to previous typing methods, acts to endorse WGS as a typing tool. CONCLUSIONS: In agreement with the current knowledge of bTB in New Zealand, transmission of M. bovis between cattle and wildlife was evident. Without direction, these estimates are less informative but taken in conjunction with the low prevalence of bTB in New Zealand's cattle population it is likely that, currently, wildlife populations are acting as the main bTB reservoir. Wildlife should therefore continue to be targeted if bTB is to be eradicated from New Zealand. WGS will be a considerable aid to bTB eradication by greatly improving the discriminatory power of molecular typing data. The substitution rates estimated here will be an important part of epidemiological investigations using WGS data.

Whole-Genome Analysis of Mycobacterium avium subsp. paratuberculosis IS900 Insertions Reveals Strain Type-Specific Modalities.

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.

Comparative Genomics of Mycobacterium avium Subspecies Paratuberculosis Sheep Strains.

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

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

Mycobacterium avium (M. a.) subsp. paratuberculosis (MAP) is a worldwide-distributed obligate pathogen in ruminants causing Johne's disease. Due to a lack of complete subtype III genome sequences, there is not yet conclusive information about genetic differences between strains of cattle (MAP-C, type II) and sheep (MAP-S) type, and especially between MAP-S subtypes I, and III. Here we present the complete, circular genome of MAP-S/type III strain JIII-386 (DE) closed by Nanopore-technology and its comparison with MAP-S/type I closed genome of strain Telford (AUS), MAP-S/type III draft genome of strain S397 (U.S.), twelve closed MAP-C strains, and eight closed M.-a.-complex-strains. Structural comparative alignments revealed clearly the mosaic nature of MAP, emphasized differences between the subtypes and the higher diversity of MAP-S genomes. The comparison of various genomic elements including transposases and genomic islands provide new insights in MAP genomics. MAP type specific phenotypic features may be attributed to genes of known large sequence polymorphisms (LSPSs) 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.

Complete Genome Sequence of the Telford Type S Strain of Mycobacterium avium subsp. paratuberculosis.

Mycobacterium avium subsp. paratuberculosis is the causative agent of Johne's disease (JD). Here, we report the complete genome sequence of Telford 9.2, a well-characterized representative strain of the M. avium subsp. paratuberculosis S subtype that is endemic in New Zealand and Australian sheep.

Pathology and molecular epidemiology of Mycobacterium pinnipedii tuberculosis in native New Zealand marine mammals.

Mycobacterium pinnipedii causes tuberculosis in a number of pinniped species, and transmission to cattle and humans has been reported. The aims of this study were to: characterize the pathology and prevalence of tuberculosis in New Zealand marine mammals; use molecular diagnostic methods to confirm and type the causal agent; and to explore relationships between type and host characteristics. Tuberculosis was diagnosed in 30 pinnipeds and one cetacean. Most affected pinnipeds had involvement of the pulmonary system, supporting inhalation as the most common route of infection, although ingestion was a possible route in the cetacean. PCR for the RD2 gene confirmed M. pinnipedii as the causal agent in 23/31 (74%) cases (22 using DNA from cultured organisms, and one using DNA from formalin-fixed paraffin-embedded (FFPE) tissue), including the first published report in a cetacean. RD2 PCR results were compared for 22 cases where both cultured organisms and FFPE tissues were available, with successful identification of M. pinnipedii in 7/22 (31.8%). In cases with moderate to large numbers of acid-fast bacilli, RD2 PCR on FFPE tissue provided a rapid, inexpensive method for confirming M. pinnipedii infection without the need for culture. VNTR typing distinguished New Zealand M. pinnipedii isolates from M. pinnipedii isolated from Australian pinnipeds and from common types of M. bovis in New Zealand. Most (16/18) M. pinnipedii isolates from New Zealand sea lions were one of two common VNTR types whereas the cetacean isolate was a type detected previously in New Zealand cattle.

Membrane segment organization in the stator complex of the flagellar motor: implications for proton flow and proton-induced conformational change.

MotA and MotB are membrane proteins that form the stator of the bacterial flagellar motor. Each motor contains several MotA 4MotB 2 complexes, which function independently to conduct protons across the membrane and couple proton flow to rotation. The mechanism of rotation is not understood in detail but is thought to involve conformational changes in the stator complexes driven by proton association/dissociation at a critical Asp residue of MotB (Asp 32 in the protein of Escherichia coli). MotA has four membrane segments and MotB has one. Previous studies using targeted disulfide cross-linking showed that the membrane segments of the two MotB subunits are together at the center of the complex, surrounded by the TM3 and TM4 segments of the four MotA subunits. Here, the cross-linking studies are extended to TM1 and TM2 of MotA, using Cys residues introduced in several positions in the segments. The observed patterns of disulfide cross-linking indicate that the TM2 segment is positioned between segments TM3 and TM4 of the same subunit, where it could contribute to the proton-channel-forming part of the structure. TM1 is at the interface between TM4 of its own subunit and the TM3 segment of another subunit, where it could stabilize the complex. A structural model based on the cross-linking results shows unobstructed pathways reaching from the periplasm to the Asp 32 residues near the inner ends of the MotB segments. The model indicates a close proximity for certain conserved, functionally important residues. The results are used to develop an explicit model for the proton-induced conformational change in the stator.

Whole Genome Sequencing for Determining the Source of Mycobacterium bovis Infections in Livestock Herds and Wildlife in New Zealand.

The ability to DNA fingerprint Mycobacterium bovis isolates helped to define the role of wildlife in the persistence of bovine tuberculosis in New Zealand. DNA fingerprinting results currently help to guide wildlife control measures and also aid in tracing the source of infections that result from movement of livestock. During the last 5 years we have developed the ability to distinguish New Zealand (NZ) M. bovis isolates by comparing the sequences of whole genome sequenced (WGS) M. bovis samples. WGS provides much higher resolution than our other established typing methods and greatly improves the definition of the regional localization of NZ M. bovis types. Three outbreak investigations are described and results demonstrate how WGS analysis has led to the confirmation of epidemiological sourcing of infection, to better definition of new sources of infection by ruling out other possible sources, and has revealed probable wildlife infection in an area considered to be free of infected wildlife. The routine use of WGS analyses for sourcing new M. bovis infections will be an important component of the strategy employed to eradicate bovine TB from NZ livestock and wildlife.

Comparison of the BBL mycobacteria growth indicator tube, the BACTEC 12B, and solid media for the isolation of Mycobacterium bovis.

We compared different methods for their ability to isolate Mycobacterium bovis from tissue samples from animals with lesions resembling bovine tuberculosis. In the first trial, M. bovis was isolated from 86 of 200 tissue samples that were cultured using 2 liquid media, BACTEC 12B and BBL mycobacteria growth indicator tube (MGIT), and a solid medium, Middlebrook 7H11 supplemented with pyruvate (7H11P). M. bovis was isolated from 2 samples with MGIT but not BACTEC 12B. M. bovis was isolated from 9 samples with BACTEC but not MGIT; these 9 samples came from the North Canterbury/Marlborough region of New Zealand. The proportion of tissues from which M. bovis was isolated with BACTEC 12B or MGIT and the mean time for isolation was different for samples from the North Canterbury/Marlborough region but not the rest of New Zealand. In the second trial, M. bovis was isolated from 401 of 1,033 tissues that were cultured using MGIT, Middlebrook 7H9 broth, or solid 7H11P. The proportion of isolates of M. bovis and the mean time for their isolation with MGIT was different for the North Canterbury/Marlborough and the rest of New Zealand. The reason for this difference was not determined but may be related to the genotypes present in this region. Genotyping using variable number tandem repeats (VNTRs) of 197 isolates of M. bovis revealed that the 44 isolates from North Canterbury/Marlborough were represented by 2 closely related VNTR types that were not found in 153 isolates from the remainder of New Zealand.

Experimental infection of New Zealand Merino sheep with a suspension of Mycobacterium avium subspecies paratuberculosis (Map) strain Telford: Kinetics of the immune response, histopathology and Map culture.

A long-term study was undertaken to monitor immune responses, faecal cultures and clinical disease in sheep experimentally infected with Mycobacterium avium subspecies paratuberculosis (Map) strain Telford. New Zealand Merino lambs (N=56) were challenged with three oral doses of Map suspension. The lambs were weighed and faecal and blood samples obtained at different time-points. At 63 weeks post-challenge, surviving sheep were euthanised and samples of liver, ileo-caecal valve and mesenteric lymph node were collected for histopathology and Map culture. High IFN-γ and antibody responses were evident as early as 8 weeks post-C1 which persisted until the end of the trial. Approximately 92% of the sheep shed Map in faeces at 36 weeks post-challenge, with the prevalence decreasing to around 40% at the end of the trial. Thirteen sheep progressively lost weight and were euthanised between weeks 32 and 58 post-challenge. Nearly 58% of surviving sheep exhibited histo-pathological lesions in at least one of the three tissues sampled, while 42% showed acid-fast bacilli in at least one tissue. A positive Map culture in at least one tissue was obtained from approximately 85% of sheep. These results indicate that the three doses of Map challenge were highly effective in establishing Johne's disease in NZ Merino lambs.

Molecular epidemiology of Mycobacterium avium subsp. paratuberculosis isolated from sheep, cattle and deer on New Zealand pastoral farms.

The present study aimed to describe the molecular diversity of Mycobacterium avium subsp. paratuberculosis (MAP) isolates obtained from sheep, cattle (beef and dairy) and deer farms in New Zealand. A total of 206 independent MAP isolates (15 beef cattle, 89 dairy cattle, 35 deer, 67 sheep) were sourced from 172 species-mobs (15 beef cattle, 66 dairy cattle, 31 deer, 60 sheep). Seventeen subtypes were identified, using a combination of variable number of tandem repeats (VNTR) and short sequence repeat (SSR) methods. Rarefaction analysis, analysis of molecular variance (AMOVA), Fst pairwise comparisons and proportional similarity index (PSI) were used to describe subtype population richness, genetic structure and potential associations between livestock sectors and New Zealand two main islands (North and South). The rarefaction analysis suggests a significantly higher subtype richness in dairy cattle herds when compared to the other livestock sectors. AMOVA results indicate that the main source of subtype variation is attributable to the livestock sector from which samples were sourced suggesting that subtypes are generally sector-specific. The pairwise Fst results were similar, with low Fst values for island differences within a livestock sector when compared to between sector analyses, representing a low subtype differentiation between islands. However, for a given island, potential associations were seen between dominant subtypes and specific livestock sectors. Three subtypes accounted for 76% of the isolates. The most common of these was isolated from sheep and beef cattle in the North Island, the second most frequent subtype was mainly isolated from dairy cattle (either island), while the third most common subtype was associated with deer farmed in the South Island. The PSI analysis suggests similarities in subtypes sourced from sheep and beef cattle. This contrasted with the isolates sourced from other livestock sectors, which tended to present sector-specific subtypes. Sheep and beef cattle were mainly infected with MAP Type I, while dairy cattle and deer were almost exclusively infected with MAP Type II. However, when beef cattle and deer were both present at farm level, they harboured similar subtypes. This study indicates that cross-species transmission of MAP occurs on New Zealand farms although close contact between species appears to be required, as in the case of sheep and beef cattle which are commonly grazed together in New Zealand.

The seal tuberculosis agent, Mycobacterium pinnipedii, infects domestic cattle in New Zealand: epidemiologic factors and DNA strain typing.

The fur seal (Arctocephalus forsteri), which is abundant in coastal areas of New Zealand, harbors several zoonotic pathogens, including Mycobacterium pinnipedii, a member of the Mycobacterium tuberculosis complex. We describe the microbiology and epidemiology of seven cases of M. pinnipedii infection in beef cattle (Bos primigenius) in coastal areas of New Zealand in 1991-2011. Epidemiologic factors were analyzed on six case farms and a telephone survey of 55 neighboring farms. A DNA-strain typing, using analysis of variable number tandem repeats and the direct repeats (VNTR/DR) of those isolates, was used to compare them to M. bovis isolates commonly found in New Zealand cattle and wildlife. In all cases of M. pinnipedii in cattle, only one animal in the herd was found to be infected. In six of seven cases, the lesions were in the thoracic lymph nodes, indicating a likely aerosol pathway. The lack of multiple cases within a herd suggests that cow-to-cow transmission is uncommon, if it occurs at all. There was no significant difference between case and control farms in distance to sea, herd size, herd type, or farming practice. The odds ratio for access to the beach for cattle on the Chatham Islands was significantly higher than it was for farms on the mainland coastal areas (odds ratio [OR] = 3.6, 95% CI = 1.1-11.4) Likewise, the odds ratio for acquiring tuberculosis was increased when farmers had seen seals on the property (OR =  9, 95% CI = 1.4-56.1 ). In all case farms, cattle had access to seals by beach grazing areas or waterways connecting directly with the ocean. The VNTR/DR typing of the isolates showed some variation in the M. pinnipedii isolates, with only two being identical; all isolates were easily distinguishable from M. bovis isolates.

Immune responses associated with progression and control of infection in calves experimentally challenged with Mycobacterium avium subsp. paratuberculosis.

This study examined the immune responses related to the infection, progression and control of Mycobacterium avium subsp. paratuberculosis (MAP) infection in calves. Twenty calves were challenged orally with MAP and 11 non-challenged calves served as controls. Approximately half the calves from each group were sacrificed at either 7 or 15 months post-challenge (PC). The majority of the challenged calves (19/20) shed MAP in feces 2-4 months PC, but thereafter fecal shedding reduced markedly. The severity of infection was reduced at 15 months PC compared to that at 7 months PC as seen from a significantly lower isolation of MAP from tissues and lower lesion scores (P<0.05). In addition, there was a reduction in the upregulation of gene expression of gamma interferon, interleukin-10 (IL-10) and inducible nitric oxide synthase from the antigen-stimulated mesenteric lymph node (MLN) cultures of the challenged calves. No evidence of infection was detected in the control calves. The severity of the infection in individual calves at 15 months PC as indicated from the number of tissue culture positive sites, was negatively related to IL-10 released from antigen-stimulated peripheral blood mononuclear cells (P<0.05). Collectively the data indicated that the severity of the MAP infection was reduced in the calves at 15 months PC and in a specific time period during infection, IL-10 may play a role in reducing the severity of this disease.

Comparison of 45 variable number tandem repeat (VNTR) and two direct repeat (DR) assays to restriction endonuclease analysis for typing isolates of Mycobacterium bovis.

Restriction endonuclease analysis (REA), developed 25 years ago for genotyping Mycobacterium bovis strains, is an important tool for bovine tuberculosis control in New Zealand. While REA gives excellent discrimination, it is technically difficult to perform compared to PCR-based typing systems which are faster and simpler to operate. Genotyping of M. bovis by the use of variable number tandem repeat loci (VNTR) and spoligotyping, either alone or together, has now become the preferred approach for typing M. bovis. Here, we evaluated the widest range of VNTR loci yet investigated for M. bovis, including two VNTR loci not previously studied, one of which (4155) had particular utility for characterizing New Zealand isolates. VNTR typing provided substantial geographical resolution of 26 of the most commonly found REA types and this was improved by the addition of two PCR assays based on parts of the direct repeat (DR) locus. Overall, 68 REA types of M. bovis common in New Zealand were discriminated into 33 VNTR/DR groups by using a minimum of nine VNTR and two DR assays. These 11 VNTR/DR assays concorded for three isolates each of 45 of the REA types but showed some variation with at least one of the VNTR/DR assays for the remaining 23 REA types. Major differences were found in allelic variation of some VNTRs between isolates from New Zealand and other countries, emphasizing the importance of adapting M. bovis typing systems to suit individual countries.

Polyphosphate kinase protects Salmonella enterica from weak organic acid stress.

Mutants of Salmonella enterica lacking polyphosphate kinase (ppk) grow poorly in the presence of the weak organic acids acetate, propionate, and benzoate. This sensitivity is corrected by methionine and seems to result from destabilization of MetA (homoserine transsuccinylase), the first enzyme in methionine biosynthesis. The MetA protein is known to be sensitive to thermal inactivation, and ppk mutants are more sensitive to heat-induced methionine auxotrophy. Peroxide increases the sensitivity of ppk mutants to both heat and acid and may oxidatively damage (carbonylate) destabilized MetA. While acid appears to impair methionine biosynthesis, it leads to derepression of MetA and may inhibit growth by causing toxic accumulation of denatured protein. This is supported by the observation that the overexpression of MetA in ppk mutants causes acid sensitivity that is not corrected by methionine. We propose that polyphosphate acts as a chemical chaperone that helps refold MetA and/or may stimulate proteolysis of toxic denatured protein. The instability of MetA protein may provide a metabolic fuse that blocks growth under conditions that denature proteins; the sensitivity of this fuse is modulated by polyphosphate.

Initial disulfide formation steps in the folding of an omega-conotoxin.

To determine whether the native disulfides of omega-conotoxins are preferentially stabilized early in the folding of these small proteins, the rates and equilibria for disulfide formation were measured for three analogues of omega-conotoxin MVIIA. In each analogue, one of the three pairs of disulfide-bonded Cys residues was replaced with Ala residues, leaving four Cys residues that can form six intermediates with one disulfide and three species with two disulfides. For each analogue, all of the disulfide-bonded species were identified, and the equilibrium constants for forming the individual species via exchange with oxidized and reduced glutathione were measured. These equilibrium constants represent effective concentrations of the Cys thiols and ranged from 0.01 to 0.4 M in the fully reduced protein. There was little or no preference for forming the native disulfides, and the equilibria for forming the first and second disulfides decreased only slightly upon the addition of 8 M urea. The data for the four-Cys analogues, together with equilibrium data for the six-Cys form, were also used to estimate effective concentrations for forming a third disulfide once two native disulfides are present. These effective concentrations were approximately 100 and 10 M in the presence of 0 and 8 M urea, respectively. The results indicate that there is little or no preferential formation of native interactions in the folding of these molecules until two disulfides have formed, after which there is a high degree of cooperativity among the native interactions.