Search

Freshwater samples (n = 199) were obtained from 41 sites with contrasting land-uses (avian, low impact, dairy, urban, sheep and beef, and mixed sheep, beef and dairy) and the E. coli phylotype of 3980 isolates (20 per water sample enrichment) was determined. Eight phylotypes were identified with B1 (48.04%), B2 (14.87%) and A (14.79%) the most abundant. Escherichia marmotae (n = 22), and Escherichia ruysiae (n = 1), were rare (0.68%) suggesting that these environmental strains are unlikely to confound water quality assessments. Phylotypes A and B1 were overrepresented in dairy and urban sites (p < 0.0001), whilst B2 were overrepresented in low impact sites (p < 0.0001). Pathogens ((Salmonella, Campylobacter, Cryptosporidium or Giardia) and the presence of diarrhoeagenic E. coli-associated genes (stx and eae) were detected in 89.9% (179/199) samples, including 80.5% (33/41) of samples with putative non-recent faecal inputs. Quantitative PCR to detect microbial source tracking targets from human, ruminant and avian contamination were concordant with land-use type and E. coli phylotype abundance. This study demonstrated that a potential recreational health risk remains where pathogens occurred in water samples with low E. coli concentration, potential non-recent faecal sources, low impact sites and where human, ruminant and avian faecal sources were absent., (© 2024 The Author(s). Environmental Microbiology Reports published by John Wiley & Sons Ltd.)

Escherichia coli is often used as a fecal indicator bacterium for water quality monitoring. We report the draft genome sequences of 500 Escherichia isolates including newly described Escherichia species, namely Escherichia marmotae , Escherichia ruysiae, and Escherichia whittamii , obtained from diverse environmental sources to assist with improved public health risk assessments., Competing Interests: The authors declare no conflict of interest.

Antimicrobial resistance is a global threat to human and animal health, with the misuse and overuse of antimicrobials suggested as the main drivers of resistance. Antimicrobial therapy can alter the bacterial community composition and the faecal resistome in cattle. Little is known about the impact of systemic antimicrobial therapy on the faecal microbiome in dairy cows in the presence of disease. Therefore, this study aimed to assess the impact of systemic antimicrobial therapy on the faecal microbiome in dairy cows in the pastoral farm environment, by analysing faecal samples from cattle impacted by several different clinically-defined conditions and corresponding antimicrobial treatments. Analysis at the individual animal level showed a decrease in bacterial diversity and richness during antimicrobial treatment but, in many cases, the microbiome diversity recovered post-treatment when the cow re-entered the milking herd. Perturbations in the microbiome composition and the ability of the microbiome to recover were specific at the individual animal level, highlighting that the animal is the main driver of variation. Other factors such as disease severity, the type and duration of antimicrobial treatment and changes in environmental factors may also impact the bovine faecal microbiome. AmpC-producing Escherichia coli were isolated from faeces collected during and post-treatment with ceftiofur from one cow while no third-generation cephalosporin resistant E. coli were isolated from the untreated cow samples. This isolation of genetically similar plasmid-mediated AmpC-producing E. coli has implications for the development and dissemination of antibiotic resistant bacteria and supports the reduction in the use of critically important antimicrobials., Competing Interests: The authors have declared that no competing interests exist., (Copyright: © 2024 Collis et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)

Whole genome sequencing has revolutionized infectious disease surveillance for tracking and monitoring the spread and evolution of pathogens. However, using a linear reference genome for genomic analyses may introduce biases, especially when studies are conducted on highly variable bacterial genomes of the same species. Pangenome graphs provide an efficient model for representing and analyzing multiple genomes and their variants as a graph structure that includes all types of variations. In this study, we present a practical bioinformatics pipeline that employs the PanGenome Graph Builder and the Variation Graph toolkit to build pangenomes from assembled genomes, align whole genome sequencing data and call variants against a graph reference. The pangenome graph enables the identification of structural variants, rearrangements, and small variants (e.g., single nucleotide polymorphisms and insertions/deletions) simultaneously. We demonstrate that using a pangenome graph, instead of a single linear reference genome, improves mapping rates and variant calling for both simulated and real datasets of the pathogen Neisseria meningitidis . Overall, pangenome graphs offer a promising approach for comparative genomics and comprehensive genetic variation analysis in infectious disease. Moreover, this innovative pipeline, leveraging pangenome graphs, can bridge variant analysis, genome assembly, population genetics, and evolutionary biology, expanding the reach of genomic understanding and applications., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2023 Yang, Guarracino, Biggs, Black, Ismail, Wold, Merriman, Prins, Garrison and de Ligt.)

Picobirnaviruses are double-stranded RNA viruses known from a wide range of host species and locations but with unknown pathogenicity and host relationships. Here, we examined the diversity of picobirnaviruses from cattle and gorillas within and around Bwindi Impenetrable Forest National Park (BIFNP), Uganda, where wild and domesticated animals and humans live in relatively close contact. We use metagenomic sequencing with bioinformatic analyses to examine genetic diversity. We compared our findings to global Picobirnavirus diversity using clustering-based analyses. Picobirnavirus diversity at Bwindi was high, with 14 near-complete RdRp and 15 capsid protein sequences, and 497 new partial viral sequences recovered from 44 gorilla samples and 664 from 16 cattle samples. Sequences were distributed throughout a phylogenetic tree of globally derived picobirnaviruses. The relationship with Picobirnavirus diversity and host taxonomy follows a similar pattern to the global dataset, generally lacking pattern with either host or geography., Competing Interests: Declaration of Competing Interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2023 The Authors. Published by Elsevier B.V. All rights reserved.)

Cryptosporidiosis is a worldwide diarrheal disease caused by the protozoan Cryptosporidium . The primary symptom is diarrhea, but patients may exhibit different symptoms based on the species of the Cryptosporidium parasite they are infected with. Furthermore, some genotypes within species are more transmissible and apparently virulent than others. The mechanisms underpinning these differences are not understood, and an effective in vitro system for Cryptosporidium culture would help advance our understanding of these differences. Using COLO-680N cells, we employed flow cytometry and microscopy along with the C. parvum -specific antibody Sporo-Glo™ to characterize infected cells 48 h following an infection with C. parvum or C. hominis . The Cryptosporidium parvum -infected cells showed higher levels of signal using Sporo-Glo™ than C. hominis -infected cells, which was likely because Sporo-Glo™ was generated against C. parvum . We found a subset of cells from infected cultures that expressed a novel, dose-dependent auto-fluorescent signal that was detectable across a range of wavelengths. The population of cells that expressed this signal increased proportionately to the multiplicity of infection. The spectral cytometry results confirmed that the signature of this subset of host cells closely matched that of oocysts present in the infectious ecosystem, pointing to a parasitic origin. Present in both C. parvum and C. hominis cultures, we named this Sig M, and due to its distinct profile in cells from both infections, it could be a better marker for assessing Cryptosporidium infection in COLO-680N cells than Sporo-Glo™. We also noted Sig M's impact on Sporo-Glo™ detection as Sporo-Glo™ uses fluoroscein-isothiocynate, which is detected where Sig M also fluoresces. Lastly, we used NanoString nCounter ® analysis to investigate the transcriptomic landscape for the two Cryptosporidium species, assessing the gene expression of 144 host and parasite genes. Despite the host gene expression being at high levels, the levels of putative intracellular Cryptosporidium gene expression were low, with no significant difference from controls, which could be, in part, explained by the abundance of uninfected cells present as determined by both Sporo-Glo™ and Sig M analyses. This study shows for the first time that a natural auto-fluorescent signal, Sig M, linked to Cryptosporidium infection can be detected in infected host cells without any fluorescent labeling strategies and that the COLO-680N cell line and spectral cytometry could be useful tools to advance the understanding of Cryptosporidium infectivity., Competing Interests: JR is employed by the company Flowjoanna Tāpui Ltd., New Zealand, and AH and NH are employed by the company AgResearch Ltd., New Zealand. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2023 Ogbuigwe, Roberts, Knox, Heiser, Pita, Haack, Garcia-Ramirez, Velathanthiri, Biggs, French and Hayman.)

Background: Tuber bruising in tetraploid potatoes (Solanum tuberosum) is a trait of economic importance, as it affects tubers' fitness for sale. Understanding the genetic components affecting tuber bruising is a key step in developing potato lines with increased resistance to bruising. As the tetraploid setting renders genetic analyses more complex, there is still much to learn about this complex phenotype. Here, we used capture sequencing data on a panel of half-sibling populations from a breeding programme to perform a genome-wide association analysis (GWAS) for tuber bruising. In addition, we collected transcriptomic data to enrich the GWAS results. However, there is currently no satisfactory method to represent both GWAS and transcriptomics analysis results in a single visualisation and to compare them with existing knowledge about the biological system under study., Results: When investigating population structure, we found that the STRUCTURE algorithm yielded greater insights than discriminant analysis of principal components (DAPC). Importantly, we found that markers with the highest (though non-significant) association scores were consistent with previous findings on tuber bruising. In addition, new genomic regions were found to be associated with tuber bruising. The GWAS results were backed by the transcriptomics differential expression analysis. The differential expression notably highlighted for the first time the role of two genes involved in cellular strength and mechanical force sensing in tuber resistance to bruising. We proposed a new visualisation, the HIDECAN plot, to integrate the results from the genomics and transcriptomics analyses, along with previous knowledge about genomic regions and candidate genes associated with the trait., Conclusion: This study offers a unique genome-wide exploration of the genetic components of tuber bruising. The role of genetic components affecting cellular strength and resistance to physical force, as well as mechanosensing mechanisms, was highlighted for the first time in the context of tuber bruising. We showcase the usefulness of genomic data from breeding programmes in identifying genomic regions whose association with the trait of interest merit further investigation. We demonstrate how confidence in these discoveries and their biological relevance can be increased by integrating results from transcriptomics analyses. The newly proposed visualisation provides a clear framework to summarise of both genomics and transcriptomics analyses, and places them in the context of previous knowledge on the trait of interest., (© 2023. The Author(s).)

Accumulation of heavy metals (HMs) in agricultural soil following the application of superphosphate fertilisers seems to induce resistance of soil bacteria to HMs and appears to co-select for resistance to antibiotics (Ab). This study aimed to investigate the selection of co-resistance of soil bacteria to HMs and Ab in uncontaminated soil incubated for 6 weeks at 25 °C in laboratory microcosms spiked with ranges of concentrations of cadmium (Cd), zinc (Zn) and mercury (Hg). Co-selection of HM and Ab resistance was assessed using plate culture on media with a range of HM and Ab concentrations, and pollution-induced community tolerance (PICT) assays. Bacterial diversity was profiled via terminal restriction fragment length polymorphism (TRFLP) assay and 16S rDNA sequencing of genomic DNA isolated from selected microcosms. Based on sequence data, the microbial communities exposed to HMs were found to differ significantly compared to control microcosms with no added HM across a range of taxonomic levels.

Objectives: To assess whether having a pet in the home is a risk factor for community-acquired urinary tract infections associated with extended-spectrum β-lactamase (ESBL)- or AmpC β-lactamase (ACBL)- producing Enterobacterales., Methods: An unmatched case-control study was conducted between August 2015 and September 2017. Cases (n = 141) were people with community-acquired urinary tract infection (UTI) caused by ESBL- or ACBL-producing Enterobacterales. Controls (n = 525) were recruited from the community. A telephone questionnaire on pet ownership and other factors was administered, and associations were assessed using logistic regression., Results: Pet ownership was not associated with ESBL- or ACBL-producing Enterobacterales-related human UTIs. A positive association was observed for recent antimicrobial treatment, travel to Asia in the previous year, and a doctor's visit in the last 6 months. Among isolates with an ESBL-/ACBL-producing phenotype, 126/134 (94%) were Escherichia coli, with sequence type 131 being the most common (47/126)., Conclusions: Companion animals in the home were not found to be associated with ESBL- or ACBL-producing Enterobacterales-related community-acquired UTIs in New Zealand. Risk factors included overseas travel, recent antibiotic use, and doctor visits., Competing Interests: Conflict of interest The authors have no conflicts of interest., (Copyright © 2022 The Author(s). Published by Elsevier Ltd.. All rights reserved.)

Insects that are freeze-tolerant start freezing at high sub-zero temperatures and produce small ice crystals. They do this using ice-nucleating agents that facilitate intercellular ice growth and prevent formation of large crystals where they can damage tissues. In Aotearoa/New Zealand the majority of cold adapted invertebrates studied survive freezing at any time of year, with ice formation beginning in the rich microbiome of the gut. Some freeze-tolerant insects are known to host symbiotic bacteria and/or fungi that produce ice-nucleating agents and we speculate that gut microbes of many New Zealand insects may provide ice-nucleating active compounds that moderate freezing. We consider too the possibility that evolutionary disparate freeze-tolerant insect species share gut microbes that are a source of ice-nucleating agents and so we describe potential transmission pathways of shared gut fauna. Despite more than 30 years of research into the freeze-tolerant mechanisms of Southern Hemisphere insects, the role of exogenous ice-nucleating agents has been neglected. Key traits of three New Zealand freeze-tolerant lineages are considered in light of the supercooling point (temperature of ice crystal formation) of microbial ice-nucleating particles, the initiation site of freezing, and the implications for invertebrate parasites. We outline approaches that could be used to investigate potential sources of ice-nucleating agents in freeze-tolerant insects and the tools employed to study insect microbiomes.

Fatigue is a common adverse effect of external beam radiation therapy in cancer patients. Mechanisms causing radiation fatigue remain unclear, although linkage to skin irradiation has been suggested. β-Endorphin, an endogenous opioid, is synthesized in skin following genotoxic ultraviolet irradiation and acts systemically, producing addiction. Exogenous opiates with the same receptor activity as β-endorphin can cause fatigue. Using rodent models of radiation therapy, exposing tails and sparing vital organs, we tested whether skin-derived β-endorphin contributes to radiation-induced fatigue. Over a 6-week radiation regimen, plasma β-endorphin increased in rats, paralleled by opiate phenotypes (elevated pain thresholds, Straub tail) and fatigue-like behavior, which was reversed in animals treated by the opiate antagonist naloxone. Mechanistically, all these phenotypes were blocked by opiate antagonist treatment and were undetected in either β-endorphin knockout mice or mice lacking keratinocyte p53 expression. These findings implicate skin-derived β-endorphin in systemic effects of radiation therapy. Opioid antagonism may warrant testing in humans as treatment or prevention of radiation-induced fatigue.

During the early postnatal period, lambs have the ability to thermoregulate body temperature via non-shivering thermogenesis through brown adipose tissue (BAT), which soon after birth begins to transform into white adipose tissue. An RNA seq approach was used to characterize the transcriptome of BAT and thyroid tissue in newborn lambs exposed to cold conditions. Fifteen newborn Romney lambs were selected and divided into three groups: group 1 ( n = 3) was a control, and groups 2 and 3 ( n = 6 each) were kept indoors for two days at an ambient temperature (20-22 °C) or at a cold temperature (4 °C), respectively. Sequencing was performed using a paired-end strategy through the BGISEQ-500 platform, followed by the identification of differentially expressed genes using DESeq2 and an enrichment analysis by g:Profiler. This study provides an in-depth expression network of the main characters involved in the thermogenesis and fat-whitening mechanisms that take place in the newborn lamb. Data revealed no significant differential expression of key thermogenic factors such as uncoupling protein 1, suggesting that the heat production peak under cold exposure might occur so rapidly and in such an immediate way that it may seem undetectable in BAT by day three of life. Moreover, these changes in expression might indicate the start of the whitening process of the adipose tissue, concluding the non-shivering thermogenesis period.

During cold exposure, brown adipose tissue (BAT) holds the key mechanism in the generation of heat, thus inducing thermogenic adaptation in response to cooler environmental changes. This process can lead to a major lipidome remodelling in BAT, where the increase in abundance of many lipid classes plays a significant role in the thermogenic mechanisms for heat production. This study aimed to identify different types of lipids, through liquid chromatography-mass spectrometry (LC-MS), in BAT and plasma during a short-term cold challenge (2-days), or not, in new-born lambs. Fifteen new-born Romney lambs were selected randomly and divided into three groups: Group 1 ( n = 3) with BAT and plasma obtained within 24 h after birth, as a control; Group 2 ( n = 6) kept indoors for two days at an ambient temperature (20-22 °C) and Group 3 ( n = 6) kept indoors for two days at a cold temperature (4 °C). Significant differences in lipid composition of many lipid categories (such as glycerolipids, glycerophospholipids, sphingolipids and sterol lipids) were observed in BAT and plasma under cold conditions, compared with ambient conditions. Data obtained from the present study suggest that short-term cold exposure induces profound changes in BAT and plasma lipidome composition of new-born lambs, which may enhance lipid metabolism via BAT thermogenic activation and adipocyte survival during cold adaptation. Further analysis on the roles of these lipid changes, validation of potential biomarkers for BAT activity, such as LPC 18:1 and PC 35:6, should contribute to the improvement of new-born lamb survival. Collectively, these observations help broaden the knowledge on the variations of lipid composition during cold exposure.

Antimicrobial resistance (AMR) is a global threat to human and animal health, with the misuse and overuse of antimicrobials being suggested as the main driver of resistance. In a global context, New Zealand (NZ) is a relatively low user of antimicrobials in animal production. However, the role antimicrobial usage on pasture-based dairy farms, such as those in NZ, plays in driving the spread of AMR within the dairy farm environment remains equivocal. Culture-based methods were used to determine the prevalence and distribution of extended-spectrum β-lactamase (ESBL)- and AmpC-producing Escherichia coli from farm environmental samples collected over a 15-month period from two NZ dairy farms with contrasting management practices. Whole genome sequencing was utilised to understand the genomic epidemiology and antimicrobial resistance gene repertoire of a subset of third-generation cephalosporin resistant E. coli isolated in this study. There was a low sample level prevalence of ESBL-producing E. coli (faeces 1.7%; farm dairy effluent, 6.7% from Dairy 4 and none from Dairy 1) but AmpC-producing E. coli were more frequently isolated across both farms (faeces 3.3% and 8.3%; farm dairy effluent 38.4%, 6.7% from Dairy 1 and Dairy 4, respectively). ESBL- and AmpC-producing E. coli were isolated from faeces and farm dairy effluent in spring and summer, during months with varying levels of antimicrobial use, but no ESBL- or AmpC-producing E. coli were isolated from bulk tank milk or soil from recently grazed paddocks. Hybrid assemblies using short- and long-read sequence data from a subset of ESBL- and AmpC-producing E. coli enabled the assembly and annotation of nine plasmids from six E. coli , including one plasmid co-harbouring 12 antimicrobial resistance genes. ESBL-producing E. coli were infrequently identified from faeces and farm dairy effluent on the two NZ dairy farms, suggesting they are present at a low prevalence on these farms. Plasmids harbouring several antimicrobial resistance genes were identified, and bacteria carrying such plasmids are a concern for both animal and public health. AMR is a burden for human, animal and environmental health and requires a holistic "One Health" approach to address., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2022 Collis, Biggs, Burgess, Midwinter, Brightwell and Cookson.)

Klebsiella pneumoniae is a Gram-negative bacterium implicated as the causative pathogen in several medical health issues with different strains causing different pathologies including pneumonia, bloodstream infections, meningitis and infections from wounds or surgery. In this study, four captive African marmosets housed in Thailand were found dead. Necropsy and histology revealed congestion of hearts, kidneys and adrenal glands. Twenty-four bacterial isolates were obtained from these four animals with all isolates yielding identical phenotypes indicative of K. pneumoniae based on classical identification schema. All the isolates show the susceptibility to amikacin, cephalexin, doxycycline, gentamicin, and enrofloxacin with intermediate susceptibility to amoxycillin/clavulanic acid. One isolate (20P167W) was chosen for genome analysis and determined to belong to sequence type 65 (ST65). The genome of 20P167W possessed multiple virulence genes including mrk gene cluster and iro and iuc gene cluster (salmochelin and aerobactin, respectively) as well as multiple antibiotic resistance genes including bla SHV -67 , bla SHV -11 , oqxA, oqxB , and fosA genes resembling those found in human isolates; this isolate has a close genetic relationship with isolates from humans in Ireland, but not from Thailand and California sea lions. Phylogenetic studies using SNP show that there was no relation between genetic and geographic distributions of all known strains typing ST65, suggesting that ST65 strains may spread worldwide through multiple international transmission events rather than by local expansions in humans and/or animals. We also predict that K. pneumoniae ST65 has an ability to acquire genetic mobile element from other bacteria, which would allow Klebsiella to become an even greater public health concern., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2022 Pinpimai, Banlunara, Roe, Dittmer, Biggs, Tantilertcharoen, Chankow, Bunpapong, Boonkam and Pirarat.)

Resistance to antimicrobials is normally caused by mutations in the drug targets or genes involved in antimicrobial activation or expulsion. Here we show that an Escherichia coli strain, named DOC14, selected for increased resistance to the bile salt sodium deoxycholate, has no mutations in any ORF, but instead has a 2.1 Mb chromosomal inversion. The breakpoints of the inversion are two inverted copies of an IS 5 element. Besides lowering deoxycholate susceptibility, the IS 5 -mediated chromosomal inversion in the DOC14 mutant was found to increase bacterial survival upon exposure to ampicillin and vancomycin, and sensitize the cell to ciprofloxacin and meropenem, but does not affect bacterial growth or cell morphology in a rich medium in the absence of antibacterial molecules. Overall, our findings support the notion that a large chromosomal inversion can benefit bacterial cells under certain conditions, contributing to genetic variability available for selection during evolution. The DOC14 mutant paired with its isogenic parental strain form a useful model as bacterial ancestors in evolution experiments to study how a large chromosomal inversion influences the evolutionary trajectory in response to various environmental stressors.

Generic Escherichia coli is commonly used as an indicator of fecal contamination to assess water quality and human health risk. Where measured E. coli exceedances occur, the presence of other pathogenic microorganisms, such as Shiga toxin-producing E. coli (STEC), is assumed, but confirmatory data are lacking. Putative E. coli isolates ( n  = 709) were isolated from water, sediment, soil, periphyton, and feces samples ( n  = 189) from five sites representing native forest and agricultural environments. Ten E. coli isolates (1.41%) were stx 2 positive, 19 (2.7%) were eae positive, and stx 1 -positive isolates were absent. At the sample level, stx 2 -positive E. coli (5 of 189, 2.6%) and eae -positive isolates (16 of 189, 8.5%) were rare. Using real-time PCR, these STEC-associated virulence factors were determined to be more prevalent in sample enrichments ( stx 1 , 23.9%; stx 2 , 31.4%; eae , 53.7%) and positively correlated with generic E. coli isolate numbers ( P  < 0.05) determined using culture-based methods. Whole-genome sequencing (WGS) was undertaken on a subset of 238 isolates with assemblies representing seven E. coli phylogroups (A, B1, B2, C, D, E, and F), 22 Escherichia marmotae isolates, and 1 Escherichia ruysiae isolate. Virulence factors, including those from extraintestinal pathogenic E. coli, were extremely diverse in isolates from the different locations and were more common in phylogroup B2. Analysis of the virulome from WGS data permitted the identification of gene repertoires that may be involved in environmental fitness and broadly align with phylogroup. Although recovery of STEC isolates was low, our molecular data indicate that they are likely to be widely present in environmental samples containing diverse E. coli phylogroups. IMPORTANCE This study takes a systematic sampling approach to assess the public health risk of Escherichia coli recovered from freshwater sites within forest and farmland. The New Zealand landscape is dominated by livestock farming, and previous work has demonstrated that "recreational exposure to water" is a risk factor for human infection by Shiga toxin-producing Escherichia coli (STEC). Though STEC isolates were rarely isolated from water samples, STEC-associated virulence factors were identified more commonly from water sample culture enrichments and were associated with increased generic E. coli concentrations. Whole-genome sequencing data from both E. coli and newly described Escherichia spp. demonstrated the presence of virulence factors from E. coli pathotypes, including extraintestinal pathogenic E. coli. This has significance for understanding and interpreting the potential health risk from E. coli where water quality is poor and suggests a role of virulence factors in survival and persistence of E. coli and Escherichia spp.

Background: Giardia intestinalis is one of the most common causes of diarrhoea worldwide. Molecular techniques have greatly improved our understanding of the taxonomy and epidemiology of this parasite. Co-infection with mixed (sub-) assemblages has been reported, however, Sanger sequencing is sometimes unable to identify shared subtypes between samples involved in the same epidemiologically linked event, due to samples showing multiple dominant subtypes within the same outbreak. Here, we aimed to use a metabarcoding approach to uncover the genetic diversity within samples from sporadic and outbreak cases of giardiasis to characterise the subtype diversity, and determine if there are common sequences shared by epidemiologically linked cases that are missed by Sanger sequencing., Methods: We built a database with 1109 unique glutamate dehydrogenase (gdh) locus sequences covering most of the assemblages of G. intestinalis and used gdh metabarcoding to analyse 16 samples from sporadic and outbreak cases of giardiasis that occurred in New Zealand between 2010 and 2018., Results: There is considerable diversity of subtypes of G. intestinalis present in each sample. The utilisation of metabarcoding enabled the identification of shared subtypes between samples from the same outbreak. Multiple variants were identified in 13 of 16 samples, with Assemblage B variants most common, and Assemblages E and A present in mixed infections., Conclusions: This study showed that G. intestinalis infections in humans are frequently mixed, with multiple subtypes present in each host. Shared sequences among epidemiologically linked cases not identified through Sanger sequencing were detected. Considering the variation in symptoms observed in cases of giardiasis, and the potential link between symptoms and (sub-) assemblages, the frequency of mixed infections could have implications for our understanding of host-pathogen interactions., (© 2022. The Author(s).)

Legionella longbeachae is an environmental bacterium that is the most clinically significant Legionella species in New Zealand (NZ), causing around two-thirds of all notified cases of Legionnaires' disease. Here we report the sequencing and analysis of the geo-temporal genetic diversity of 54 L. longbeachae serogroup 1 (sg1) clinical isolates, derived from cases from around NZ over a 22-year period, including one complete genome and its associated methylome. The 54 sg1 isolates belonged to two main clades that last shared a common ancestor between 95 BCE and 1694 CE. There was diversity at the genome-structural level, with large-scale arrangements occurring in some regions of the chromosome and evidence of extensive chromosomal and plasmid recombination. This includes the presence of plasmids derived from recombination and horizontal gene transfer between various Legionella species, indicating there has been both intra- and inter-species gene flow. However, because similar plasmids were found among isolates within each clade, plasmid recombination events may pre-empt the emergence of new L. longbeachae strains. Our complete NZ reference genome consisted of a 4.1 Mb chromosome and a 108 kb plasmid. The genome was highly methylated with two known epigenetic modifications, m 4 C and m 6 A, occurring in particular sequence motifs within the genome., (© 2022. The Author(s).)

Background: Campylobacter is a genus of bacteria that has been isolated from the gastrointestinal tract of humans and animals, and the environments they inhabit around the world. Campylobacter adapt to new environments by changes in their gene content and expression, but little is known about how they adapt to long-term human colonization. In this study, the genomes of 31 isolates from a New Zealand patient and 22 isolates from a United Kingdom patient belonging to Campylobacter jejuni sequence type 45 (ST45) were compared with 209 ST45 genomes from other sources to identify the mechanisms by which Campylobacter adapts to long-term human colonization. In addition, the New Zealand patient had their microbiota investigated using 16S rRNA metabarcoding, and their level of inflammation and immunosuppression analyzed using biochemical tests, to determine how Campylobacter adapts to a changing gastrointestinal tract., Results: There was some evidence that long-term colonization led to genome degradation, but more evidence that Campylobacter adapted through the accumulation of non-synonymous single nucleotide polymorphisms (SNPs) and frameshifts in genes involved in cell motility, signal transduction and the major outer membrane protein (MOMP). The New Zealand patient also displayed considerable variation in their microbiome, inflammation and immunosuppression over five months, and the Campylobacter collected from this patient could be divided into two subpopulations, the proportion of which correlated with the amount of gastrointestinal inflammation., Conclusions: This study demonstrates how genomics, phylogenetics, 16S rRNA metabarcoding and biochemical markers can provide insight into how Campylobacter adapts to changing environments within human hosts. This study also demonstrates that long-term human colonization selects for changes in Campylobacter genes involved in cell motility, signal transduction and the MOMP; and that genetically distinct subpopulations of Campylobacter evolve to adapt to the changing gastrointestinal environment., (© 2021. The Author(s).)

Multidrug resistance is a major threat to global elimination of tuberculosis (TB). We performed phenotypic drug-susceptibility testing and whole-genome sequencing for 309 isolates from 342 consecutive patients who were given a diagnosis of TB in Yangon, Myanmar, during July 2016‒June 2018. We identified isolates by using the GeneXpert platform to evaluate drug-resistance profiles. A total of 191 (62%) of 309 isolates had rifampin resistance; 168 (88%) of these rifampin-resistant isolates were not genomically related, indicating the repeated emergence of resistance in the population, rather than extensive local transmission. We did not detect resistance mutations to new oral drugs, including bedaquiline and pretomanid. The current GeneXpert MTB/RIF system needs to be modified by using the newly launched Xpert MTB/XDR cartridge or line-probe assay. Introducing new oral drugs to replace those currently used in treatment regimens for multidrug-resistant TB will also be useful for treating TB in Myanmar.

Epidemiological studies of communicable diseases increasingly use large whole-genome sequencing (WGS) datasets to explore the transmission of pathogens. It is important to obtain an initial overview of datasets and identify closely related isolates, but this can be challenging with large numbers of isolates and imperfect sequencing. We used an ad hoc whole-genome multi locus sequence typing method to summarise data from a longitudinal study of Staphylococcus aureus in a primary school in New Zealand. Each pair of isolates was compared and the number of genes where alleles differed between isolates was tallied to produce a matrix of "allelic differences". We plotted histograms of the number of allelic differences between isolates for: all isolate pairs; pairs of isolates from different individuals; and pairs of isolates from the same individual. 340 sequenced isolates were included, and the ad hoc shared genome contained 445 genes. There were between 0 and 420 allelic differences between isolate pairs and the majority of pairs had more than 260 allelic differences. We found many genetically closely related S. aureus isolates from single individuals and a smaller number of closely-related isolates from separate individuals. Multiple S. aureus isolates from the same individual were usually very closely related or identical over the ad hoc shared genome. Siblings carried genetically similar, but not identical isolates. An ad hoc shared genome approach to WGS analysis can accommodate imperfect sequencing of the included isolates, and can provide insights into relationships between isolates in epidemiological studies with large WGS datasets containing diverse isolates., (© 2021. The Author(s).)

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., (Copyright © 2021 Elsevier GmbH. All rights reserved.)

The management of competition horses in New Zealand often involves rotations of short periods of stall confinement and concentrate feeding, with periods of time at pasture. Under these systems, horses may undergo abrupt dietary changes, with the incorporation of grains or concentrate feeds to the diet to meet performance needs, or sudden changes in the type of forage fed in response to a lack of fresh or conserved forage. Abrupt changes in dietary management are a risk factor for gastrointestinal (GI) disturbances, potentially due to the negative effects observed on the population of GI microbiota. In the present study, the faecal microbiota of horses was investigated to determine how quickly the bacterial communities; (1) responded to dietary change, and (2) stabilised following abrupt dietary transition. Six Thoroughbred mares were stabled for six weeks, consuming freshly cut pasture (weeks 1, 3 and 5), before being abruptly transitioned to conserved forage-based diets, both offered ad libitum. Intestinal markers were administered to measure digesta transit time immediately before each diet change. The conserved forage-based diets were fed according to a 3 × 3 Latin square design (weeks 2, 4 and 6), and comprised a chopped ensiled forage fed exclusively (Diet FE) or with whole oats (Diet FE + O), and perennial ryegrass hay fed with whole oats (Diet H + O). Faecal samples were collected at regular intervals from each horse following the diet changes. High throughput 16S rRNA gene sequencing was used to evaluate the faecal microbiota. There were significant differences in alpha diversity across diets ( p < 0.001), and a significant effect of diet on the beta diversity (ANOSIM, p = 0.001), with clustering of samples observed by diet group. There were differences in the bacterial phyla across diets ( p < 0.003), with the highest relative abundances observed for Firmicutes (62-64%) in the two diets containing chopped ensiled forage, Bacteroidetes (32-38%) in the pasture diets, and Spirochaetes (17%) in the diet containing hay. Major changes in relative abundances of faecal bacteria appeared to correspond with the cumulative percentage of intestinal markers retrieved in the faeces as the increasing amounts of digesta from each new diet transited the animals. A stable faecal microbiota profile was observed in the samples from 96 h after abrupt transition to the treatment diets containing ensiled chopped forage. The present study confirmed that the diversity and community structure of the faecal bacteria in horses is diet-specific and resilient following dietary transition and emphasised the need to have modern horse feeding management that reflects the ecological niche, particularly by incorporating large proportions of forage into equine diets.

Seasonal variation in the faecal microbiota of forage-fed horses was investigated over a 12-month period to determine whether the bacterial diversity fluctuated over time. Horses ( n = 10) were maintained on pasture for one year, with hay supplemented from June to October. At monthly intervals, data were recorded on pasture availability and climate (collected continuously and averaged on monthly basis), pasture and hay samples were collected for nutrient analysis, and faecal samples were collected from all horses to investigate the diversity of faecal microbiota using next-generation sequencing on the Illumina MiSeq platform. The alpha diversity of bacterial genera was high in all samples ( n = 118), with significantly higher Simpson's ( p < 0.001) and Shannon-Wiener ( p < 0.001) diversity indices observed during the months when horses were kept exclusively on pasture compared to the months when pasture was supplemented with hay. There were significant effects of diet, season, and month (ANOSIM, p < 0.01 for each comparison) on the beta diversity of bacterial genera identified in the faeces. While there was some inter-horse variation, hierarchical clustering of beta diversity indices showed separate clades originating for samples obtained during May, June, and July (late-autumn to winter period), and January, February, and March (a period of drought), with a strong association between bacterial taxa and specific nutrients (dry matter, protein, and structural carbohydrates) and climate variables (rainfall and temperature). Our study supports the hypothesis that the diversity and community structure of the faecal microbiota of horses kept on pasture varied over a 12-month period, and this variation reflects changes in the nutrient composition of the pasture, which in turn is influenced by climatic conditions. The findings of this study may have implications for grazing management and the preparation of conserved forages for those horses susceptible to perturbations of the hindgut microbiota.

Members of the Clostridiales R-7 group are abundant bacterial residents of the rumen microbiome; however, they are poorly characterized. We report the complete genome sequences of three members of the R-7 group, FE2010, FE2011, and XBB3002, isolated from the ruminal contents of pasture-grazed dairy cows in New Zealand.

Extended-spectrum beta-lactamase (ESBL)-producing Enterobacteriaceae are considered a critical priority by the World Health Organization. Presented here are two genome sequences of Escherichia coli strains isolated from New Zealand freshwater. The genome sequences' mean size was 5.2 Mb, with a mean of 4,848 coding sequences. Both genomes carried the ESBL bla CTX-M gene.

Cattle are asymptomatic carriers of Shiga toxin-producing Escherichia coli (STEC) strains that can cause serious illness or death in humans. In New Zealand, contact with cattle feces and living near cattle populations are known risk factors for human STEC infection. Contamination of fresh meat with STEC strains also leads to the potential for rejection of consignments by importing countries. We used a combination of PCR/matrix-assisted laser desorption ionization-time of flight mass spectrometry (MALDI-TOF) and whole-genome sequencing (WGS) to evaluate the presence and transmission of STEC on farms and in processing plants to better understand the potential pathways for human exposure and thus mitigate risk. Animal and environmental samples ( n  = 2,580) were collected from six farms and three meat processing plants in New Zealand during multiple sampling sessions in spring of 2015 and 2016. PCR/MALDI-TOF analysis revealed that 6.2% were positive for "Top 7" STEC. Top 7 STEC strains were identified in all sample sources ( n  = 17) tested. A marked increase in Top 7 STEC prevalence was observed between calf hides on farm (6.3% prevalence) and calf hides at processing plants (25.1% prevalence). Whole-genome sequencing was performed on Top 7 STEC bacterial isolates ( n  = 40). Analysis of STEC O26 ( n  = 25 isolates) revealed relatively low genetic diversity on individual farms, consistent with the presence of a resident strain disseminated within the farm environment. Public health efforts should focus on minimizing human contact with fecal material on farms and during handling, transport, and slaughter of calves. Meat processing plants should focus on minimizing cross-contamination between the hides of calves in a cohort during transport, lairage, and slaughter. IMPORTANCE Cattle are asymptomatic carriers of Shiga toxin-producing E. coli (STEC) strains, which can cause serious illness or death in humans. Contact with cattle feces and living near cattle are known risk factors for human STEC infection. This study evaluated STEC carriage in young calves and the farm environment with an in-depth evaluation of six farms and three meat processing plants over 2 years. An advanced molecular detection method and whole-genome sequencing were used to provide a detailed evaluation of the transmission of STEC both within and between farms. The study revealed widespread STEC contamination within the farm environment, but no evidence of recent spread between farms. Contamination of young dairy calf hides increased following transport and holding at meat processing plants. The elimination of STEC in farm environments may be very difficult given the multiple transmission routes; interventions should be targeted at decreasing fecal contamination of calf hides during transport, lairage, and processing., (Copyright © 2021 Browne et al.)

Canine acute hemorrhagic diarrhea syndrome (AHDS) has been associated in some studies with Clostridioides perfringens overgrowth and toxin-mediated necrosis of the intestinal mucosa. We aimed to determine the effect of a single fecal microbiota transplantation (FMT) on clinical scores and fecal microbiomes of 1 and 7 dogs with AHDS from New Zealand and South Africa. We hypothesized that FMT would improve AHDS clinical scores and increase microbiota alpha-diversity and short-chain fatty acid (SCFA)-producing microbial communities' abundances in dogs with AHDS after FMT. We sequenced the V3-V4 region of the 16S-rRNA gene in the feces of AHDS FMT-recipients and sham-treated control dogs, and their healthy donors at admission, discharge, and 30 days post-discharge. There were no significant differences in median AHDS clinical scores between FMT-recipients and sham-treated controls at admission or discharge (P = 0.22, P = 0.41). At admission, the Shannon diversity index (SDI) was lower in AHDS dogs than healthy donors (P = 0.002). The SDI did not change from admission to 30 days in sham-treated dogs yet increased in FMT-recipients from admission to discharge (P = 0.04) to levels not different than donors (P = 0.33) but significantly higher than sham-treated controls (P = 0.002). At 30 days, the SDI did not differ between FMT recipients, sham-treated controls, and donors (P = 0.88). Principal coordinate analysis of the Bray-Curtis index separated post-FMT and donor dogs from pre-FMT and sham-treated dogs (P = 0.009) because of increased SCFA-producing genera's abundances after FMT. A single co-abundance subnetwork contained many of the same OTUs found to be differentially abundant in FMT-recipients, and the abundance of this module was increased in FMT-recipients at discharge and 30 days, compared to sham-treated controls. We conclude in this small pilot study FMT did not have any clinical benefit. A single FMT procedure has the potential to increase bacterial communities of SCFA-producing genera important for intestinal health up to 30 days post-FMT., Competing Interests: Dr. Burchell’s current affiliation did not alter our adherence to PLOS ONE policies on sharing data and materials. None of the authors have or had any conflict of interest with regards to conception, funding or execution of this study, analysis and interpretation of the results, or writing the manuscript of this study.

Draft genomes of seven extended-spectrum β-lactamase (ESBL)-producing Escherichia coli strains recovered from New Zealand waterways are described. The mean genome size was 5.1 Mb, with 4,724 coding sequences. All genomes contained the ESBL gene bla CTX-M , and one carried a plasmid-mediated AmpC gene, bla CMY-2 A multidrug-resistant genotype was detected in three isolates., (Copyright © 2021 Burgess et al.)

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., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2021 Mizzi, Timms, Price-Carter, Gautam, Whittington, Heuer, Biggs and Plain.)

Extended-spectrum-beta-lactamase (ESBL)- or AmpC beta-lactamase (ACBL)-producing Escherichia coli bacteria are the most common cause of community-acquired multidrug-resistant urinary tract infections (UTIs) in New Zealand. The carriage of antimicrobial-resistant bacteria has been found in both people and pets from the same household; thus, the home environment may be a place where antimicrobial-resistant bacteria are shared between humans and pets. In this study, we sought to determine whether members (pets and people) of the households of human index cases with a UTI caused by an ESBL- or ACBL-producing E. coli strain also carried an ESBL- or ACBL-producing Enterobacteriaceae strain and, if so, whether it was a clonal match to the index case clinical strain. Index cases with a community-acquired UTI were recruited based on antimicrobial susceptibility testing of urine isolates. Fecal samples were collected from 18 non-index case people and 36 pets across 27 households. Eleven of the 27 households screened had non-index case household members (8/18 people and 5/36 animals) positive for ESBL- and/or ACBL-producing E. coli strains. Whole-genome sequence analysis of 125 E. coli isolates (including the clinical urine isolates) from these 11 households showed that within seven households, the same strain of ESBL-/ACBL-producing E. coli was cultured from both the index case and another person (5/11 households) or pet dog (2/11 households). These results suggest that transmission within the household may contribute to the community spread of ESBL- or ACBL-producing E. coli IMPORTANCE Enterobacteriaceae that produce extended-spectrum beta-lactamases (ESBLs) and AmpC beta-lactamases (ACBLs) are important pathogens and can cause community-acquired illnesses, such as urinary tract infections (UTIs). Fecal carriage of these resistant bacteria by companion animals may pose a risk for transmission to humans. Our work evaluated the sharing of ESBL- and ACBL-producing E. coli isolates between humans and companion animals. We found that in some households, dogs carried the same strain of ESBL-producing E. coli as the household member with a UTI. This suggests that transmission events between humans and animals (or vice versa) are likely occurring within the home environment and, therefore, the community as a whole. This is significant from a health perspective, when considering measures to minimize community transmission, and highlights that in order to manage community spread, we need to consider interventions at the household level., (Copyright © 2020 Toombs-Ruane et al.)

Campylobacter enteritis in humans is primarily associated with C. jejuni/coli infection. Other species cause campylobacteriosis relatively infrequently; while this could be attributed to bias in diagnostic methods, the pathogenicity of non- jejuni/coli Campylobacter spp. such as C. upsaliensis and C. helveticus (isolated from dogs and cats) is uncertain. Galleria mellonella larvae are suitable models of the mammalian innate immune system and have been applied to C. jejuni studies. This study compared the pathogenicity of C. jejuni , C. upsaliensis , and C. helveticus isolates. Larvae inoculated with either C. upsaliensis or C. helveticus showed significantly higher survival than those inoculated with C. jejuni . All three Campylobacter species induced indistinguishable histopathological changes in the larvae. C. jejuni could be isolated from inoculated larvae up to eight days post-inoculation whereas C. upsaliensis and C. helveticus could only be isolated in the first two days. There was a significant variation in the hazard rate between batches of larvae, in Campylobacter strains, and in biological replicates as random effects, and in species and bacterial dose as fixed effects. The Galleria model is applicable to other Campylobacter spp. as well as C. jejuni, but may be subject to significant variation with all Campylobacter species. While C. upsaliensis and C. helveticus cannot be considered non-pathogenic, they are significantly less pathogenic than C. jejuni .

Shiga toxin-producing Escherichia coli (STEC) are a leading cause of foodborne outbreaks of human disease, but they reside harmlessly as an asymptomatic commensal in the ruminant gut. STEC serogroup O145 are difficult to isolate as routine diagnostic methods are unable to distinguish non-O157 serogroups due to their heterogeneous metabolic characteristics, resulting in under-reporting which is likely to conceal their true prevalence. In light of these deficiencies, the purpose of this study was a twofold approach to investigate enhanced STEC O145 diagnostic culture-based methods: firstly, to use a genomic epidemiology approach to understand the genetic diversity and population structure of serogroup O145 at both a local (New Zealand) (n = 47) and global scale (n = 75) and, secondly, to identify metabolic characteristics that will help the development of a differential media for this serogroup. Analysis of a subset of E. coli serogroup O145 strains demonstrated considerable diversity in carbon utilisation, which varied in association with eae subtype and sequence type. Several carbon substrates, such as D-serine and D-malic acid, were utilised by the majority of serogroup O145 strains, which, when coupled with current molecular and culture-based methods, could aid in the identification of presumptive E. coli serogroup O145 isolates. These carbon substrates warrant subsequent testing with additional serogroup O145 strains and non-O145 strains. Serogroup O145 strains displayed extensive genetic heterogeneity that was correlated with sequence type and eae subtype, suggesting these genetic markers are good indicators for distinct E. coli phylogenetic lineages. Pangenome analysis identified a core of 3,036 genes and an open pangenome of >14,000 genes, which is consistent with the identification of distinct phylogenetic lineages. Overall, this study highlighted the phenotypic and genotypic heterogeneity within E. coli serogroup O145, suggesting that the development of a differential media targeting this serogroup will be challenging., Competing Interests: Two of the authors (ALC and GB) are employed by a commercial company, AgResearch Ltd. This does not alter our adherence to PLOS ONE policies on sharing data and materials. The remaining authors have declared that no competing interests exist.

Summary: We present sismonr, an R package for an integral generation and simulation of in silico biological systems. The package generates gene regulatory networks, which include protein-coding and non-coding genes along with different transcriptional and post-transcriptional regulations. The effect of genetic mutations on the system behaviour is accounted for via the simulation of genetically different in silico individuals. The ploidy of the system is not restricted to the usual haploid or diploid situations but can be defined by the user to higher ploidies. A choice of stochastic simulation algorithms allows us to simulate the expression profiles of the genes in the in silico system. We illustrate the use of sismonr by simulating the anthocyanin biosynthesis regulation pathway for three genetically distinct in silico plants., Availability and Implementation: The sismonr package is implemented in R and Julia and is publicly available on the CRAN repository (https://CRAN.R-project.org/package=sismonr). A detailed tutorial is available from GitHub at https://oliviaab.github.io/sismonr/., Supplementary Information: Supplementary data are available at Bioinformatics online., (© The Author(s) 2020. Published by Oxford University Press. All rights reserved. For permissions, please e-mail: journals.permissions@oup.com.)

The draft genome sequences for eight isolates of Campylobacter helveticus isolated from companion animals are described and compared with that of the type strain. On average, the genomes are 1,825,025 bp long and have a GC content of 34.4% and 1,885 coding DNA sequences (CDSs). CRISPRs were detected in only one isolate and phages in none., (Copyright © 2020 Bojanić et al.)

Bacteria within the digestive tract of adult honey bees are likely to play a key role in the digestion of sugar-rich foods. However, the influence of diet on honey bee gut bacteria is not well understood. During periods of low floral abundance, beekeepers often supplement the natural sources of carbohydrate that honey bees collect, such as nectar, with various forms of carbohydrates such as sucrose (a disaccharide) and invert sugar (a mixture of the monosaccharides glucose and fructose). We compared the effect of these sugar supplements on the relative abundance of bacteria in the gut of bees by feeding bees from a single colony, two natural diets: mānuka honey, a monofloral honey with known antibacterial properties, and a hive diet; and artificial diets of invert sugar, sucrose solution, and sucrose solutions containing synthesised compounds associated with the antibacterial properties of mānuka honey. 16S ribosomal RNA (rRNA)-based sequencing showed that dietary regimes containing mānuka honey, sucrose and invert sugar did not alter the relative abundance of dominant core bacteria after 6 days of being fed these diets. However, sucrose-rich diets increased the relative abundances of three sub-dominant core bacteria, Rhizobiaceae, Acetobacteraceae, and Lactobacillus kunkeei, and decreased the relative abundance of Frischella perrara, all which significantly altered the bacterial composition. Acetogenic bacteria from the Rhizobiaceae and Acetobacteraceae families increased two- to five-fold when bees were fed sucrose. These results suggest that sucrose fuels the proliferation of specific low abundance primary sucrose-feeders, which metabolise sugars into monosaccharides, and then to acetate., Competing Interests: We have the following interests as this study received funding from New Zealand Apples & Pears Inc, Honey New Zealand (International) Limited, seven regional New Zealand Beekeeping clubs, Beeline Apiaries, Arataki Honey, Hawke’s Bay, Galatea Apiaries, and NZ Beeswax Ltd. The invert sugar used in the research was supplied free of charge by Natural Sugars NZ. There are no patents, products in development or marketed products to declare. This does not alter our adherence to the PLOS ONE policies on sharing data and materials.

In 2014, antimicrobial drug-resistant Campylobacter jejuni sequence type 6964 emerged contemporaneously in poultry from 3 supply companies in the North Island of New Zealand and as a major cause of campylobacteriosis in humans in New Zealand. This lineage, not previously identified in New Zealand, was resistant to tetracycline and fluoroquinolones. Genomic analysis revealed divergence into 2 major clades; both clades were associated with human infection, 1 with poultry companies A and B and the other with company C. Accessory genome evolution was associated with a plasmid, phage insertions, and natural transformation. We hypothesize that the tetO gene and a phage were inserted into the chromosome after conjugation, leaving a remnant plasmid that was lost from isolates from company C. The emergence and rapid spread of a resistant clone of C. jejuni in New Zealand, coupled with evolutionary change in the accessory genome, demonstrate the need for ongoing Campylobacter surveillance among poultry and humans.

Draft genomes of five Campylobacter isolates recovered from New Zealand brushtail possums are described. Genome sizes ranged from 1.591 Mbp to 1.594 Mbp, with G+C contents of 29.9% to 29.95%. Comparison to Australian Campylobacter 16S rRNA gene sequences suggests that the species may be common to possums., (Copyright © 2019 Wilkinson et al.)

The global spread of multidrug-resistant enterobacteria warrants new strategies to combat these pathogens. One possible approach is the reconsideration of "old" antimicrobials, which remain effective after decades of use. Synthetic 5-nitrofurans such as furazolidone, nitrofurantoin, and nitrofurazone are such a class of antimicrobial drugs. Recent epidemiological data showed a very low prevalence of resistance to this antimicrobial class among clinical Escherichia coli isolates in various parts of the world, forecasting the increasing importance of its uses to battle antibiotic-resistant enterobacteria. However, although they have had a long history of clinical use, a detailed understanding of the 5-nitrofurans' mechanisms of action remains limited. Nitrofurans are known as prodrugs that are activated in E. coli by reduction catalyzed by two redundant nitroreductases, NfsA and NfsB. Furazolidone, nevertheless, retains relatively significant antibacterial activity in the nitroreductase-deficient Δ nfsA Δ nfsB E. coli strain, indicating the presence of additional activating enzymes and/or antibacterial activity of the unreduced form. Using genome sequencing, genetic, biochemical, and bioinformatic approaches, we discovered a novel 5-nitrofuran-activating enzyme, AhpF, in E. coli The discovery of a new nitrofuran-reducing enzyme opens new avenues for overcoming 5-nitrofuran resistance, such as designing nitrofuran analogues with higher affinity for AhpF or screening for adjuvants that enhance AhpF expression., (Copyright © 2019 American Society for Microbiology.)

The longhorned tick, Haemaphysalis longicornis , feeds upon a wide range of bird and mammalian hosts. Mammalian hosts include cattle, deer, sheep, goats, humans, and horses. This tick is known to transmit a number of pathogens causing tick-borne diseases, and was the vector of a recent serious outbreak of oriental theileriosis in New Zealand. A New Zealand-USA consortium was established to sequence, assemble, and annotate the genome of this tick, using ticks obtained from New Zealand's North Island. In New Zealand, the tick is considered exclusively parthenogenetic and this trait was deemed useful for genome assembly. Very high molecular weight genomic DNA was sequenced on the Illumina HiSeq4000 and the long-read Pac Bio Sequel platforms. Twenty-eight SMRT cells produced a total of 21.3 million reads which were assembled with Canu on a reserved supercomputer node with access to 12 TB of RAM, running continuously for over 24 days. The final assembly dataset consisted of 34,211 contigs with an average contig length of 215,205 bp. The quality of the annotated genome was assessed by BUSCO analysis, an approach that provides quantitative measures for the quality of an assembled genome. Over 95% of the BUSCO gene set was found in the assembled genome. Only 48 of the 1066 BUSCO genes were missing and only 9 were present in a fragmented condition. The raw sequencing reads and the assembled contigs/scaffolds are archived at the National Center for Biotechnology Information., (© 2019 The Author(s).)

The use of culture methods to detect Escherichia coli diversity does not provide sufficient resolution to identify strains present at low levels. Here, we target the hypervariable gnd gene and describe a database containing 534 distinct partial gnd sequences and associated O groups for use with culture-independent E. coli community analysis., (Copyright © 2019 Cookson et al.)

Colorectal cancer is a heterogenous and mostly sporadic disease, the development of which is associated with microbial dysbiosis. Recent advances in subtype classification have successfully stratified the disease using molecular profiling. To understand potential relationships between molecular mechanisms differentiating the subtypes of colorectal cancer and composition of gut microbial community, we classified a set of 34 tumour samples into molecular subtypes using RNA-sequencing gene expression profiles and determined relative abundances of bacterial taxonomic groups. To identify bacterial community composition, 16S rRNA amplicon metabarcoding was used as well as whole genome metagenomics of the non-human part of RNA-sequencing data. The generated data expands the collection of the data sources related to the disease and connects molecular aspects of the cancer with environmental impact of microbial community.

Campylobacter spp. are frequently found associated with the avian intestinal tract. Most are commensals, but some can cause human campylobacteriosis. Here, we report the draft genome sequences of three strains of a novel Campylobacter sp. isolated from urban birds and a rural river in New Zealand., (Copyright © 2019 Wilkinson et al.)

Escherichia coli P50 is a canine uropathogenic isolate sampled in the Wellington region of New Zealand. We report the draft genome sequence of this isolate, which contains characteristic virulence genes for urinary tract infections and is predicted to be capable of causing human infections.