Search

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

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.

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.

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.

Abstract 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

Abstract 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, m4C and m6A, occurring in particular sequence motifs within the genome.

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

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.

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

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 blaSHV−67, blaSHV−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.

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.

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.

The current study evaluated the bacterial diversity of six hot water spring clusters in Sri Lanka by Illumina MiSeq sequencing of the V3–V4 region of the 16S rRNA gene. Bacterial abundance measures and diversity statistics were assessed using QIIME2 metagenomics workflow, and the results were compared according to the region, the water temperature at the surface (36–59 °C), and pH (6.25–8.35). The predominant phyla observed were Proteobacteria, Actinobacteria, Firmicutes, [Thermi], and Cyanobacteria. A low abundance of Bacteroidetes, Chloroflexi, Acidobacteria, TM7, and Spirochaetes was detected in most of the springs. Several important bacterial species such as Deinococcus geothermalis that can tolerate Martian-like conditions, genera such as Legionella and Campylobacter that contain pathogenic species, sulfur metabolizing Desulfovibrio, Desulfatirhabdium, Desulforhabdus, Desulfacinum, Thermodesulfovibrio, Desulfovirga, and Thiobacter species, and several other species with the potential practical industrial application were detected. Several opportunistic human pathogens were detected in the water samples and raised a public health concern about the management of post-bathing water. Based on the Bray Curtis beta diversity metric, the microbial distribution correlated with temperature rather than the geographic distance. This study provides valuable new insights into the bacterial diversity of the hot springs in Sri Lanka. Future research needs to be conducted on industrially important thermophiles identified in this study.

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.

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.

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.

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.

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.

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.

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.

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. Keywords: Tick genome, Pac Bio de novo assembly, Genome annotation, Cattle tick

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.

During 1998–2012, an extended outbreak of Salmonella enterica serovar Typhimurium definitive type 160 (DT160) affected >3,000 humans and killed wild birds in New Zealand. However, the relationship between DT160 within these 2 host groups and the origin of the outbreak are unknown. Whole-genome sequencing was used to compare 109 Salmonella Typhimurium DT160 isolates from sources throughout New Zealand. We provide evidence that DT160 was introduced into New Zealand around 1997 and rapidly propagated throughout the country, becoming more genetically diverse over time. The genetic heterogeneity was evenly distributed across multiple predicted functional protein groups, and we found no evidence of host group differentiation between isolates collected from human, poultry, bovid, and wild bird sources, indicating ongoing transmission between these host groups. Our findings demonstrate how a comparative genomic approach can be used to gain insight into outbreaks, disease transmission, and the evolution of a multihost pathogen after a probable point-source introduction.

Abstract Current culture methods to investigate changes in Escherichia coli community structure are often slow and laborious. Genes such as gnd (6-phosphogluconate dehydrogenase) have a highly variable nucleotide sequence and may provide a target for E. coli microbiome analysis using culture-independent methods. Metabarcoded PCR primers were used to generate separate libraries from calf faecal samples for high throughput sequencing. Although a total of 348 separate gnd sequence types (gSTs) were identified, 188 were likely to be due to sequencing errors. Of the remaining 160 gSTs, 92 did not match those in a database of 319 separate gnd sequences. ‘Animal’ was the main determinant of E. coli diversity with limited impact of sample type or DNA extraction method on intra-host E. coli community variation from faeces and recto-anal mucosal swab samples. This culture-independent study has addressed the difficulties of quantifying bacterial intra-species diversity and revealed that, whilst individual animals may harbour >50 separate E. coli strains, communities are dominated by

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.