Development of a multiplex droplet digital PCR assay for simultaneous detection and quantification of Escherichia coli, E. marmotae, and E. ruysiae in water samples

Escherichia coli are widely used by water quality managers as Fecal Indicator Bacteria, but current quantification methods do not differentiate them from benign, environmental Escherichia species such as E. marmotae (formerly named cryptic clade V) or E. ruysiae (cryptic clades III and IV). Reliable and specific techniques for their identification are required to avoid confounding microbial water quality assessments. To address this, a multiplex droplet digital PCR (ddPCR) assay targeting lipB (E. coli and E. ruysiae) and bglC (E. marmotae) was designed. The ddPCR performance was assessed using in silico analysis; genomic DNA from 40 local, international, and reference strains of target and non-target coliforms; and spiked water samples in a range relevant to water quality managers (1 to 1000 cells/100 mL). Results were compared to an analogous quantitative PCR (qPCR) and the Colilert method. Both PCR assays showed excellent sensitivity with a limit of detection of 0.05 pg/μL and 0.005 pg/μl for ddPCR and qPCR respectively, and of quantification of 0.5 pg/μL of genomic DNA. The ddPCR allowed differentiation and quantification of three Escherichia species per run by amplitude multiplexing and showed a high concordance with concentrations measured by Colilert once proportional bias was accounted for. In silico specificity testing underlined the possibility to further detect and distinguish Escherichia cryptic clade VI. Finally, the applicability of the ddPCR was successfully tested on environmental water samples where E. marmotae and E. ruysiae potentially confound E. coli counts based on the Most Probable Number method, highlighting the utility of this novel ddPCR as an efficient and rapid discriminatory test to improve water quality assessments.

Genomic analysis of canine pneumoviruses and canine respiratory coronavirus from New Zealand

AIMS

To isolate canine respiratory coronavirus (CRCoV) and canine pneumovirus (CnPnV) in cell culture and to compare partial genomic sequences of CRCoV and CnPnV from New Zealand with those from other countries.

METHODS

Oropharyngeal swab samples from dogs affected by canine infectious respiratory disease syndrome that were positive for CnPnV (n = 15) or CRCoV (n = 1) by virus-specific reverse transcriptase quantitative PCR (RT-qPCR) in a previous study comprised the starting material. Virus isolation was performed in HRT-18 cells for CRCoV and RAW 264.7 and Vero cells for CnPnV. The entire sequence of CnPnV G protein (1,266 nucleotides) and most (8,063/9,707 nucleotides) of the 3' region of CRCoV that codes for 10 structural and accessory proteins were amplified and sequenced. The sequences were analysed and compared with other sequences available in GenBank using standard molecular tools including phylogenetic analysis.

RESULTS

Virus isolation was unsuccessful for both CRCoV and CnPnV. Pneumovirus G protein was amplified from 3/15 (20%) samples that were positive for CnPnV RNA by RT-qPCR. Two of these (NZ-048 and NZ-049) were 100% identical to each other, and 90.9% identical to the third one (NZ-007). Based on phylogenetic analysis of the G protein gene, CnPnV NZ-048 and NZ-049 clustered with sequences from the USA, Thailand and Italy in group A, and CnPnV NZ-007 clustered with sequences from the USA in group B. The characteristics of the predicted genes (length, position) and their putative protein products (size, predicted structure, presence of N- and O-glycosylation sites) of the New Zealand CRCoV sequence were consistent with those reported previously, except for the region located between open reading frame (ORF)3 (coding for S protein) and ORF6 (coding for E protein). The New Zealand virus was predicted to encode 5.9 kDa, 27 kDa and 12.7 kDa proteins, which differed from the putative coding capacity of this region reported for CRCoV from other countries.

CONCLUSIONS

This report represents the first characterisation of partial genomic sequences of CRCoV and CnPnV from New Zealand. Our results suggest that the population of CnPnV circulating in New Zealand is not homogeneous, and that the viruses from two clades described overseas are also present here. Limited conclusions can be made based on only one CRCoV sequence, but the putative differences in the coding capacity of New Zealand CRCoV support the previously reported variability of this region. The reasons for such variability and its biological implications need to be further elucidated.

Draft genome sequences of Escherichia spp. isolates from New Zealand environmental sources

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.

Quantifying replication slippage error in Cryptosporidium metabarcoding studies

Genetic variation in Cryptosporidium, a common protozoan gut parasite in humans, is often based on marker genes containing trinucleotide repeats, which differentiate subtypes and track outbreaks. However, repeat regions have high replication slippage rates, making it difficult to discern biological diversity from error. Here, we synthesised Cryptosporidium DNA in clonal plasmid vectors, amplified them in different mock community ratios and sequenced them using next generation sequencing to determine the rate of replication slippage with dada2. Our results indicate that slippage rates increase with the length of the repeat region and can contribute to error rates of up to 20%.

Impact of systemic antimicrobial therapy on the faecal microbiome in symptomatic dairy cows

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.

Species-wide genomics of kākāpō provides tools to accelerate recovery

The kākāpō is a critically endangered, intensively managed, long-lived nocturnal parrot endemic to Aotearoa New Zealand. We generated and analysed whole-genome sequence data for nearly all individuals living in early 2018 (169 individuals) to generate a high-quality species-wide genetic variant callset. We leverage extensive long-term metadata to quantify genome-wide diversity of the species over time and present new approaches using probabilistic programming, combined with a phenotype dataset spanning five decades, to disentangle phenotypic variance into environmental and genetic effects while quantifying uncertainty in small populations. We find associations for growth, disease susceptibility, clutch size and egg fertility within genic regions previously shown to influence these traits in other species. Finally, we generate breeding values to predict phenotype and illustrate that active management over the past 45 years has maintained both genome-wide diversity and diversity in breeding values and, hence, evolutionary potential. We provide new pathways for informing future conservation management decisions for kākāpō, including prioritizing individuals for translocation and monitoring individuals with poor growth or high disease risk. Overall, by explicitly addressing the challenge of the small sample size, we provide a template for the inclusion of genomic data that will be transformational for species recovery efforts around the globe.

Pangenome graphs in infectious disease: a comprehensive genetic variation analysis of Neisseria meningitidis leveraging Oxford Nanopore long reads

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.

Abundant dsRNA picobirnaviruses show little geographic or host association in terrestrial systems

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.

Investigating the genetic components of tuber bruising in a breeding population of tetraploid potatoes

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.

Aristaeella hokkaidonensis gen. nov. sp. nov. and Aristaeella lactis sp. nov., two rumen bacterial species of a novel proposed family, Aristaeellaceae fam. nov

Two strains of Gram-negative, anaerobic, rod-shaped bacteria, from an abundant but uncharacterized rumen bacterial group of the order 'Christensenellales', were phylogenetically and phenotypically characterized. These strains, designated R-7^T and WTE2008^T, shared 98.6-99.0 % sequence identity between their 16S rRNA gene sequences. R-7^T and WTE2008^T clustered together on a distinct branch from other Christensenellaceae strains and had <88.1 % sequence identity to the closest type-strain sequence from Luoshenia tenuis NSJ-44^T. The genome sequences of R-7^T and WTE2008^T had 83.6 % average nucleotide identity to each other, and taxonomic assignment using the Genome Taxonomy Database indicates these are separate species within a novel family of the order 'Christensenellales'. Cells of R-7^T and WTE2008^T lacked any obvious appendages and their cell wall ultra-structures were characteristic of Gram-negative bacteria. The five most abundant cellular fatty acids of both strains were C_16 : 0, C_16 : 0 iso, C_17 : 0 anteiso, C_18 : 0 and C_15 : 0 anteiso. The strains used a wide range of the 23 soluble carbon sources tested, and grew best on cellobiose, but not on sugar-alcohols. Xylan and pectin were fermented by both strains, but not cellulose. Acetate, hydrogen, ethanol and lactate were the major fermentation end products. R-7^T produced considerably more hydrogen than WTE2008^T, which produced more lactate. Based on these analyses, Aristaeellaceae fam. nov. and Aristaeella gen. nov., with type species Aristaeella hokkaidonensis sp. nov., are proposed. Strains R-7^T (=DSM 112795^T=JCM 34733^T) and WTE2008^T (=DSM 112788^T=JCM 34734^T) are the proposed type strains for Aristaeella hokkaidonensis sp. nov. and Aristaeella lactis sp. nov., respectively.

Co-Selection of Bacterial Metal and Antibiotic Resistance in Soil Laboratory Microcosms

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.

Extended-spectrum β-lactamase- and AmpC β-lactamase-producing Enterobacterales associated with urinary tract infections in the New Zealand community: a case-control study

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.

Insect Freeze-Tolerance Downunder: The Microbial Connection

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.

Urinary MicroRNA Analysis Indicates an Epigenetic Regulation of Chronic Kidney Disease of Unknown Etiology in Sri Lanka

BACKGROUND

Chronic kidney disease of unknown etiology (CKDu) is reported among male paddy farmers in the dry zone of Sri Lanka. The exact cause of this disease remains undetermined. Genetic susceptibility is identified as a major risk factor for CKDu Objectives: In this study, small urinary RNAs were characterized in CKDu patients, healthy endemic and non-endemic controls. Differently expressed urinary miRNAs and their associated pathways were identified in the study population.

METHODS

Healthy and diseased male volunteers (n = 9) were recruited from Girandurukotte (endemic) and Mawanella (non-endemic) districts. Urinary small RNAs were purified and sequenced using Illumina MiSeqTM. The sequence trace files were assembled and analyzed. Differentially ex-pressed miRNAs among these three groups were identified and pathway analysis was conducted.

RESULTS

The urine samples contained 130,623 sequence reads identified as non-coding RNAs, PIWI-interacting RNAs (piRNA), and miRNAs. Approximately four percent of the total small RNA reads represented miRNA, and 29% represented piRNA. A total of 409 miRNA species were ex-pressed in urine. Interestingly, both diseased and endemic controls population showed significantly low expression of miRNA and piRNA. Regardless of the health status, the endemic population ex-pressed significantly low levels of miR-10a, miR-21, miR-148a, and miR-30a which have been linked with several environmental toxins Conclusion: Significant downregulation of miRNA and piRNA expression in both diseased and healthy endemic samples indicates an epigenetic regulation of CKDu involving genetic and environmental interaction. Further studies of specific miRNA species are required to develop a miRNA panel to identify individuals susceptible to CKDu.

Mass Spectrometry-Based Lipidomics of Brown Adipose Tissue and Plasma of New-Born Lambs Subjected to Short-Term Cold Exposure

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.

Prevalence and distribution of extended-spectrum β-lactamase and AmpC-producing Escherichia coli in two New Zealand dairy farm environments

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.

Genetic characterization of hypervirulent Klebsiella pneumoniae responsible for acute death in captive marmosets

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.

A large chromosomal inversion affects antimicrobial sensitivity of Escherichia coli to sodium deoxycholate

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 IS5 element. Besides lowering deoxycholate susceptibility, the IS5-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.

Uncovering the genetic diversity of Giardia intestinalis in isolates from outbreaks in New Zealand

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.

Graphical Abstract

Supplementary Information

The online version contains supplementary material available at 10.1186/s40249-022-00969-x.

Genomic adaptations of Campylobacter jejuni to long-term human colonization

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.

Aliarcobacter, Halarcobacter, Malaciobacter, Pseudarcobacter and Poseidonibacter are later synonyms of Arcobacter: transfer of Poseidonibacter parvus, Poseidonibacter antarcticus, 'Halarcobacter arenosus', and 'Aliarcobacter vitoriensis' to Arcobacter as Arcobacter parvus comb. nov., Arcobacter antarcticus comb. nov., Arcobacter arenosus comb. nov. and Arcobacter vitoriensis comb. nov

This paper re-examines the taxonomic positions of recently described Poseidonibacter (P. parvum and P. antarcticus), Aliarcobacter ('Al. vitoriensis'), Halarcobacter ('H. arenosus') and Arcobacter (A. caeni, A. lacus) species, and other species proposed to represent novel genera highly related to the genus Arcobacter. Phylogenomic and several overall genome relatedness indices (OGRIs) were applied to a total of 118 representative genomes for this purpose. Phylogenomic analyses demonstrated the Arcobacter clade to be distinct from other Epsilonproteobacteria, clearly defined and containing closely related species. Aliarcobacter butzleri and Malaciobacter pacificus did not cluster with other members of these proposed genera, indicating incoherence of these genera. Every OGRI measure applied indicated a high level of relatedness among all Arcobacter clade species, including the recently described taxa studied here, and substantially lower between type species representatives for other Epsilonproteobacteria. Where published guidelines were available, OGRI values for Arcobacter clade species were either unsupportive of division into other genera or were at the lowest boundary range (for average amino acid identity). We propose that Aliarcobacter, Halarcobacter, Malaciobacter, Pseudarcobacter, Poseidonibacter and Arcobacter sensu stricto be considered members of a single genus, Arcobacter, and subsequently transfer P. parvum, P. antarcticus, 'Al. vitoriensis' and 'H. arenosus' to Arcobacter as Arcobacter parvum comb. nov., Arcobacter antarcticus comb. nov., Arcobacter vitoriensis comb. nov. and Arcobacter arenosus comb. nov.

Genomic Profiling of Mycobacterium tuberculosis Strains, Myanmar

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.

Whole-genome sequencing and ad hoc shared genome analysis of Staphylococcus aureus isolates from a New Zealand primary school

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.

Resilience of Faecal Microbiota in Stabled Thoroughbred Horses Following Abrupt Dietary Transition between Freshly Cut Pasture and Three Forage-Based Diets

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.

A molecular survey of canine respiratory viruses in New Zealand

AIMS

The aim of this study was to identify viruses associated with canine infectious respiratory disease syndrome (CIRDS) among a population of New Zealand dogs.

METHODS

Convenience samples of oropharyngeal swabs were collected from 116 dogs, including 56 CIRDS-affected and 60 healthy dogs from various locations in New Zealand between March 2014 and February 2016. Pooled samples from CIRDS-affected (n = 50) and from healthy (n = 50) dogs were tested for the presence of canine respiratory viruses using next generation sequencing (NGS). Individual samples (n = 116) were then tested by quantitative PCR (qPCR) and reverse transcriptase qPCR (RT-qPCR) for specific viruses. Groups were compared using Fisher's exact or χ² tests. The effect of explanatory variables (age, sex, type of household, presence of viral infection) on the response variable (CIRDS-affected or not) was tested using RR.

RESULTS

Canine pneumovirus (CnPnV), canine respiratory coronavirus (CRCoV), canine herpesvirus-1 (CHV-1), canine picornavirus and influenza C virus sequences were identified by NGS in the pooled sample from CIRDS-affected but not healthy dogs. At least one virus was detected by qPCR/RT-qPCR in 20/56 (36%) samples from CIRDS dogs and in 23/60 (38%) samples from healthy dogs (p = 0.84). CIRDS-affected dogs were most commonly positive for CnPnV (14/56, 25%) followed by canine adenovirus-2 (CAdV-2, 5/56, 9%), canine parainfluenza virus (CpiV) and CHV-1 (2/56, 4% each), and CRCoV (1/56, 2%). Only CnPnV (17/60, 28%) and CAdV-2 (14/60, 23%) were identified in samples from healthy dogs, and CAdV-2 was more likely to be detected healthy than diseased dogs (RR 0.38; 95% CI = 0.15-0.99; p = 0.045).

CONCLUSIONS

The frequency of detection of viruses traditionally linked to CIRDS (CAdV-2 and CPiV) among diseased dogs was low. This suggests that other pathogens are likely to have contributed to development of CIRDS among sampled dogs. Our data represent the first detection of CnPnV in New Zealand, but the role of this virus in CIRDS remains unclear. On-going monitoring of canine respiratory pathogens by NGS would be beneficial, as it allows rapid detection of novel viruses that may be introduced to the New Zealand canine population in the future. Such monitoring could be done using pooled samples to minimise costs.

CLINICAL RELEVANCE

Testing for novel respiratory viruses such as CnPnV and CRCoV should be considered in all routine laboratory investigations of CIRDS cases, particularly in dogs vaccinated with currently available kennel cough vaccines.

One dog's waste is another dog's wealth: A pilot study of fecal microbiota transplantation in dogs with acute hemorrhagic diarrhea syndrome

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.

Novel mechanisms of TolC-independent decreased bile-salt susceptibility in Escherichia coli

Bile salts, including sodium deoxycholate (DOC), are secreted into the intestine to aid fat digestion and contribute to antimicrobial protection. Gram-negative pathogens such as Escherichia coli, however, are highly resistant to DOC, using multiple mechanisms of which the multidrug efflux pump AcrAB-TolC is the dominant one. Given that TolC-mediated efflux masks the interaction of DOC with potential targets, we sought to identify those targets by identifying genes whose mutations cause an increase in the MIC to DOC relative to the ∆tolC parental strain, that lacks TolC-associated functional efflux pumps. Using a mutant screen, we isolated twenty independent spontaneous mutants that had a higher MICDOC than the E. coli parental ∆tolC strain. Whole genome sequencing of these mutants mapped most mutations to the ptsI or cyaA gene. Analysis of knock-out mutants and complementation showed that elimination of PtsI, a component of the carbohydrate phosphotransferase system, or one of the two key proteins involved in cAMP synthesis and signaling, adenylate cyclase (CyaA) or cAMP receptor protein (Crp) causes low-level increased resistance of a ∆tolC E. coli strain to DOC.

Genetic characterisation of Campylobacter concisus: Strategies for improved genomospecies discrimination

Although at least two genetically distinct groups, or genomospecies, have been well documented for Campylobacter concisus, no phenotype has yet been identified for their differentiation and thus formal description as separate species. C. concisus has been isolated from a variety of sites in the human body, including saliva and stool samples from both healthy and diarrhoeic individuals. We evaluated the ability of a range of whole genome-based tools to distinguish between the two C. concisus genomospecies (GS) using a collection of 190 C. concisus genomes. Nine genomes from related Campylobacter species were included in some analyses to provide context. Analyses incorporating sequence analysis of multiple ribosomal genes generated similar levels of C. concisus GS discrimination as genome-wide comparisons. The C. concisus genomes formed two groups; GS1 represented by ATCC 33237^T and GS2 by CCUG 19995. The two C. concisus GS were separated from the nine genomes of related species. GS1 and GS2 also differed in G+C content with medians of 37.56% and 39.51%, respectively. The groups are consistent with previously established GS and are supported by DNA reassociation results. Average Nucleotide Identity using MUMmer (ANIm) and Genome BLAST Distance Phylogeny generated in silico DNA-DNA hybridisation (isDDH) (against ATCC 33237^T and CCUG 19995), plus G+C content provides cluster-independent GS discrimination suitable for routine use. Pan-genomic analysis identified genes specific to GS1 and GS2. WGS data and genomic species identification methods support the existence of two GS within C. concisus. These data provide genome-level metrics for strain identification to genomospecies level.

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.

Carriage of Extended-Spectrum-Beta-Lactamase- and AmpC Beta-Lactamase-Producing Escherichia coli Strains from Humans and Pets in the Same Households

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. coliIMPORTANCEEnterobacteriaceae 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.

Does land use affect pathogen presence in New Zealand drinking water supplies?

Four microbes (Campylobacter spp., Escherichia coli, Cryptosporidium spp. and Giardia spp.) were monitored in 16 waterways that supply public drinking water for 13 New Zealand towns and cities. Over 500 samples were collected from the abstraction point at each study site every three months between 2009 and 2019. The waterways represent a range from small to large, free flowing to reservoir impoundments, draining catchments of entirely native vegetation to those dominated by pastoral agriculture. We used machine learning algorithms to explore the relative contribution of land use, catchment geology, vegetation, topography, and water quality characteristics of the catchment to determining the abundance and/or presence of each microbe. Sites on rivers draining predominantly agricultural catchments, the Waikato River, Oroua River and Waiorohi Stream had all four microbes present, often in high numbers, throughout the sampling interval. Other sites, such as the Hutt River and Big Huia Creek in Wellington which drain catchments of native vegetation, never had pathogenic microbes detected, or unsafe levels of E. coli. Boosted Regression Tree models could predict abundances and presence/absence of all four microbes with good precision using a wide range of potential environmental predictors covering land use, geology, vegetation, topography, and nutrient concentrations. Models were more accurate for protozoa than bacteria but did not differ markedly in their ability to predict abundance or presence/absence. Environmental drivers of microbe abundance or presence/absence also differed depending on whether the microbe was protozoan or bacterial. Protozoa were more prevalent in waterways with lower water quality, higher numbers of ruminants in the catchment, and in September and December. Bacteria were more abundant with higher rainfall, saturated soils, and catchments with greater than 35% of the land in agriculture. Although modern water treatment protocols will usually remove many pathogens from drinking water, several recent outbreaks of waterborne disease due to treatment failures, have highlighted the need to manage water supplies on multiple fronts. This research has identified potential catchment level variables, and thresholds, that could be better managed to reduce the potential for pathogens to enter drinking water supplies.

sismonr: simulation of in silico multi-omic networks with adjustable ploidy and post-transcriptional regulation in R

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.

Comparison of the Pathogenic Potential of Campylobacter jejuni, C. upsaliensis and C. helveticus and Limitations of Using Larvae of Galleria mellonella as an Infection Model

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.

Phylogeographic Analysis Reveals Multiple International transmission Events Have Driven the Global Emergence of Escherichia coli O157:H7

BACKGROUND

Shiga toxin-producing Escherchia coli (STEC) O157:H7 is a zoonotic pathogen that causes numerous food and waterborne disease outbreaks. It is globally distributed, but its origin and the temporal sequence of its geographical spread are unknown.

METHODS

We analyzed whole-genome sequencing data of 757 isolates from 4 continents, and performed a pan-genome analysis to identify the core genome and, from this, extracted single-nucleotide polymorphisms. A timed phylogeographic analysis was performed on a subset of the isolates to investigate its worldwide spread.

RESULTS

The common ancestor of this set of isolates occurred around 1890 (1845-1925) and originated from the Netherlands. Phylogeographic analysis identified 34 major transmission events. The earliest were predominantly intercontinental, moving from Europe to Australia around 1937 (1909-1958), to the United States in 1941 (1921-1962), to Canada in 1960 (1943-1979), and from Australia to New Zealand in 1966 (1943-1982). This pre-dates the first reported human case of E. coli O157:H7, which was in 1975 from the United States.

CONCLUSIONS

Inter- and intra-continental transmission events have resulted in the current international distribution of E. coli O157:H7, and it is likely that these events were facilitated by animal movements (eg, Holstein Friesian cattle). These findings will inform policy on action that is crucial to reduce the further spread of E. coli O157:H7 and other (emerging) STEC strains globally.

Microbial contamination in drinking water at public outdoor recreation facilities in New Zealand

AIM

The aim of our study was to assess the presence and risk of waterborne pathogens in the drinking water of outdoor facilities in New Zealand and track potential sources of microbial contamination in water sources.

METHODS AND RESULTS

A serial cross-sectional study with a risk-based sample collection strategy was conducted at 15 public campgrounds over two summer seasons (2011-2012 and 2012-2013). Drinking water supplied to these campgrounds was not compliant with national standards, based on Escherichia coli as an indicator organism, in more than half of the sampling occasions. Campylobacter contamination of drinking water at the campgrounds was likely to be of wild bird origin. Faecal samples from rails (pukeko and weka) were 35 times more likely to return a Campylobacter-positive result compared to passerines. Water treatment using ultraviolet (UV) irradiation or a combination of filtration and UV irradiation or chemicals was more likely to result in water that was compliant with the national standards than water from a tap without any treatment. The use of filters alone was not associated with the likelihood of compliance.

CONCLUSIONS

Providing microbiologically safe drinking water at outdoor recreational facilities is imperative to avoid gastroenteritis outbreaks. This requires an in-depth understanding of potential sources of contamination in drinking water sources and the installation of adequate water treatment facilities.

SIGNIFICANCE AND IMPACT OF THE STUDY

Our study provides evidence that drinking water without treatment or filter-only treatment in public campgrounds is unlikely to comply with national standards for human consumption and extra water treatment measures such as UV irradiation or chemical treatment are needed.

A critical rebuttal of the proposed division of the genus Arcobacter into six genera using comparative genomic, phylogenetic, and phenotypic criteria

The proposal to restructure the genus Arcobacter into six distinct genera was critically examined using: comparative analyses of up to 80 Epsilonproteobacterial genome sequences (including 26 arcobacters); phylogenetic analyses of three housekeeping genes and also 342 core genes; and phenotypic criteria. Genome sequences were analysed with tools to calculate Percentage of Conserved Proteins, Average Amino-acid Identity, BLAST-based Average Nucleotide Identity, in silico DNA-DNA hybridisation values, genome-wide Average Nucleotide Identity, Alignment Fractions and G+C percentages. Genome analyses revealed the genus Arcobacter sensu lato to be relatively homogenous, and phylogenetic analyses clearly distinguished the group from other Epsilonproteobacteria. Genomic distinction of the genera proposed by Pérez-Cataluña et al. [2018] was not supported by any of the measures used and a subsequent risk of strain misidentification clearly identified. Similarly, phenotypic analyses supported the delineation of Arcobacter sensu lato but did not justify the position of the proposed novel genera. The present polyphasic taxonomic study strongly supports the continuance of the classification of "aerotolerant campylobacters" as Arcobacter and refutes the proposed genus-level subdivision of Pérez-Cataluña et al. [2018].

Genomic epidemiology and carbon metabolism of Escherichia coli serogroup O145 reflect contrasting phylogenies

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.

Genomic Analysis of Fluoroquinolone- and Tetracycline-Resistant Campylobacter jejuni Sequence Type 6964 in Humans and Poultry, New Zealand, 2014-2016

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.

Use of Genomics to Investigate Historical Importation of Shiga Toxin-Producing Escherichia coli Serogroup O26 and Nontoxigenic Variants into New Zealand

Shiga toxin-producing Escherichia coli serogroup O26 is an important public health pathogen. Phylogenetic bacterial lineages in a country can be associated with the level and timing of international imports of live cattle, the main reservoir. We sequenced the genomes of 152 E. coli O26 isolates from New Zealand and compared them with 252 E. coli O26 genomes from 14 other countries. Gene variation among isolates from humans, animals, and food was strongly associated with country of origin and stx toxin profile but not isolation source. Time of origin estimates indicate serogroup O26 sequence type 21 was introduced at least 3 times into New Zealand from the 1920s to the 1980s, whereas nonvirulent O26 sequence type 29 strains were introduced during the early 2000s. New Zealand's remarkably fewer introductions of Shiga toxin-producing Escherichia coli O26 compared with other countries (such as Japan) might be related to patterns of trade in live cattle.

Long-term Colonization by Campylobacter jejuni Within a Human Host: Evolution, Antimicrobial Resistance, and Adaptation

Background

Campylobacteriosis is inflammation of the gastrointestinal tract as a result of Campylobacter infection. Most campylobacteriosis cases are acute and self-limiting, with Campylobacter excretion ceasing a few weeks after symptoms cease. We identified a patient with fecal specimens positive for Campylobacter jejuni (ST45) intermittently during a 10-year period.

Methods

Sixteen Campylobacter isolates were collected from the patient during 2006-2016. The isolates' genomes were sequenced to determine their relatedness, and their antimicrobial susceptibility patterns and motility were measured to determine the effects of antibiotic therapy and long-term excretion on the Campylobacter population.

Results

Phylogenetic analyses estimated that the isolates shared a date of common ancestor between 1998 and 2006, coinciding with the onset of symptoms for the patient. Genomic analysis identified selection for changes in motility, and antimicrobial susceptibility testing suggested that the Campylobacter population developed resistance to several antibiotics coinciding with periods of antibiotic therapy.

Conclusions

The patient was consistently colonized with organisms from a Campylobacter population that adapted to the internal environment of the patient. Genomic and phylogenetic analyses can give insight into a patient's infection history and the effect of antimicrobial treatment on Campylobacter populations in this unusual situation of long-term colonization of an individual.

Genomic correlates of extraintestinal infection are linked with changes in cell morphology in Campylobacter jejuni

Campylobacter jejuni is the most common cause of bacterial diarrheal disease in the world. Clinical outcomes of infection can range from asymptomatic infection to life-threatening extraintestinal infections. This variability in outcomes for infected patients has raised questions as to whether genetic differences between C. jejuni isolates contribute to their likelihood of causing severe disease. In this study, we compare the genomes of ten C. jejuni isolates that were implicated in extraintestinal infections with reference gastrointestinal isolates, in order to identify unusual patterns of sequence variation associated with infection outcome. We identified a collection of genes that display a higher burden of uncommon mutations in invasive isolates compared with gastrointestinal close relatives, including some that have been previously linked to virulence and invasiveness in C. jejuni. Among the top genes identified were mreB and pgp1, which are both involved in determining cell shape. Electron microscopy confirmed morphological differences in isolates carrying unusual sequence variants of these genes, indicating a possible relationship between extraintestinal infection and changes in cell morphology.

Gene Regulatory Networks: A Primer in Biological Processes and Statistical Modelling

Modelling gene regulatory networks requires not only a thorough understanding of the biological system depicted, but also the ability to accurately represent this system from a mathematical perspective. Throughout this chapter, we aim to familiarize the reader with the biological processes and molecular factors at play in the process of gene expression regulation. We first describe the different interactions controlling each step of the expression process, from transcription to mRNA and protein decay. In the second section, we provide statistical tools to accurately represent this biological complexity in the form of mathematical models. Among other considerations, we discuss the topological properties of biological networks, the application of deterministic and stochastic frameworks, and the quantitative modelling of regulation. We particularly focus on the use of such models for the simulation of expression data that can serve as a benchmark for the testing of network inference algorithms.

Isolation of emerging Campylobacter species in working farm dogs and their frozen home-killed raw meat diets

We applied 7 culture methods to 50 working farm dog fecal samples and 6 methods to 50 frozen home-killed raw meat diet samples to optimize recovery of a wide range of Campylobacter spp. Culture methods combined filtration, enrichment broths, and agars at 37°C and 42°C in conventional and hydrogen-enriched microaerobic atmospheres. Overall, a prevalence of 62% (31 of 50) and 6% (3 of 50) was detected in dog and meat samples, respectively, based on Campylobacter genus PCR. A total of 356 Campylobacter spp. isolates were recovered from dogs, with successful isolation by individual methods ranging from 2 to 25 dogs. The species detected most commonly were C. upsaliensis and C. jejuni, and less commonly C. coli and C. lari. Species isolated that are rarely reported from dogs included C. rectus, C. lari subsp. concheus, C. volucris, and Helicobacter winghamensis. Six isolates from dogs positive by Campylobacter genus PCR were confirmed, using 16S rRNA sequencing, as Arcobacter cryaerophilus (1) and Arcobacter butzleri (5). C. jejuni multi-locus sequence typing results revealed a diversity of sequence types in working dogs, with several uncommonly reported from other C. jejuni sources in New Zealand. Overall, 20 isolates from 3 meat samples were positive by Campylobacter genus PCR; 1 meat sample was positive for C. jejuni, 1 for C. rectus, and 1 isolate was subsequently identified as A. butzleri. The method using Campylobacter enrichment broth in a hydrogen-enriched environment on nonselective agar resulted in significantly reduced recovery of Campylobacter spp. from both sample types.

Extended-Spectrum Beta-Lactamase-Producing Enterobacteriaceae in Dairy Farm Environments: A New Zealand Perspective

Antimicrobial resistance (AMR) is a global issue for both human and animal health. Infections caused by antimicrobial-resistant bacteria present treatment option challenges and are often associated with heightened severity of infection. Antimicrobial use (AMU) in human and animal health is a main driver for the development of antimicrobial-resistant bacteria. Increasing levels of AMU and the development and spread of AMR in food-producing animals, especially in poultry and swine production, has been identified as a food safety risk, but dairy production systems have been less studied. A number of farm management practices may impact on animal disease and as a result can influence the use of antimicrobials and subsequently AMR prevalence. However, this relationship is multifactorial and complex. Several AMR transmission pathways between dairy cattle, the environment, and humans have been proposed, including contact with manure-contaminated pastures, direct contact, or through the food chain from contaminated animal-derived products. The World Health Organization has defined a priority list for selected bacterial pathogens of concern to human health according to 10 criteria relating to health and AMR. This list includes human pathogens such as the extended-spectrum beta-lactamase-producing Enterobacteriaceae (ESBL-E), which can be associated with dairy cattle, their environment, as well as animal-derived food products. ESBL-E represent a potential risk to human and animal health and an emerging food safety concern. This review addresses two areas; first, the current understanding of the role of dairy farming in the prevalence and spread of AMR is considered, highlighting research gaps using ESBL-E as an exemplar; and second, a New Zealand perspective is taken to examine how farm management practices may contribute to on-farm AMU and AMR in dairy cattle.