Browsing by Author "Biggs PJ"

Now showing 1 - 20 of 32
  • Thumbnail Image
    Item
    A large chromosomal inversion affects antimicrobial sensitivity of Escherichia coli to sodium deoxycholate
    (Microbiology Society, 2022-08-12) Le VVH; León-Quezada RI; Biggs PJ; Rakonjac J
    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.
  • Thumbnail Image
    Item
    A novel, stain-free, natural auto-fluorescent signal, Sig M, identified from cytometric and transcriptomic analysis of infectivity of Cryptosporidium hominis and Cryptosporidium parvum.
    (Frontiers Media S.A., 2023-05-22) Ogbuigwe P; Roberts JM; Knox MA; Heiser A; Pita A; Haack NA; Garcia-Ramirez JC; Velathanthiri N; Biggs PJ; French NP; Hayman DTS; Xu R
    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.
  • Thumbnail Image
    Item
    Abundant dsRNA picobirnaviruses show little geographic or host association in terrestrial systems.
    (Elsevier, 2023-08) Knox MA; Wierenga J; Biggs PJ; Gedye K; Almeida V; Hall R; Kalema-Zikusoka G; Rubanga S; Ngabirano A; Valdivia-Granda W; Hayman DTS
    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.
  • Thumbnail Image
    Item
    Aristaeella hokkaidonensis gen. nov. sp. nov. and Aristaeella lactis sp. nov., two rumen bacterial species of a novel proposed family, Aristaeellaceae fam. nov.
    (Microbiology Society, 2023-05-12) Mahoney-Kurpe SC; Palevich N; Noel SJ; Gagic D; Biggs PJ; Soni P; Reid PM; Koike S; Kobayashi Y; Janssen PH; Attwood GT; Moon CD
    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-7T and WTE2008T, shared 98.6-99.0 % sequence identity between their 16S rRNA gene sequences. R-7T and WTE2008T 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-44T. The genome sequences of R-7T and WTE2008T 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-7T and WTE2008T 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 C16 : 0, C16 : 0 iso, C17 : 0 anteiso, C18 : 0 and C15 : 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-7T produced considerably more hydrogen than WTE2008T, 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-7T (=DSM 112795T=JCM 34733T) and WTE2008T (=DSM 112788T=JCM 34734T) are the proposed type strains for Aristaeella hokkaidonensis sp. nov. and Aristaeella lactis sp. nov., respectively.
  • Thumbnail Image
    Item
    Assessing antimicrobial resistance in pasture-based dairy farms: a 15-month surveillance study in New Zealand.
    (American Society for Microbiology, 2024-10-23) Collis RM; Biggs PJ; Burgess SA; Midwinter AC; Liu J; Brightwell G; Cookson AL
    Antimicrobial resistance is a global public and animal health concern. Antimicrobial resistance genes (ARGs) have been detected in dairy farm environments globally; however, few longitudinal studies have utilized shotgun metagenomics for ARG surveillance in pasture-based systems. This 15-month study aimed to undertake a baseline survey using shotgun metagenomics to assess the relative abundance and diversity of ARGs in two pasture-based dairy farm environments in New Zealand with different management practices. There was no statistically significant difference in overall ARG relative abundance between the two dairy farms (P = 0.321) during the study period. Compared with overseas data, the relative abundance of ARG copies per 16S rRNA gene in feces (0.08-0.17), effluent (0.03-0.37), soil (0.20-0.63), and bulk tank milk (0.0-0.12) samples was low. Models comparing the presence or absence of resistance classes found in >10% of all feces, effluent, and soil samples demonstrated no statistically significant associations (P > 0.05) with "season," and only multi-metal (P = 0.020) and tetracycline (P = 0.0003) resistance were significant at the "farm" level. Effluent samples harbored the most diverse ARGs, some with a recognized public health risk, whereas soil samples had the highest ARG relative abundance but without recognized health risks. This highlights the importance of considering the genomic context and risk of ARGs in metagenomic data sets. This study suggests that antimicrobial resistance on pasture-based dairy farms is low and provides essential baseline ARG surveillance data for such farming systems. IMPORTANCE: Antimicrobial resistance is a global threat to human and animal health. Despite the detection of antimicrobial resistance genes (ARGs) in dairy farm environments globally, longitudinal surveillance in pasture-based systems remains limited. This study assessed the relative abundance and diversity of ARGs in two New Zealand dairy farms with different management practices and provided important baseline ARG surveillance data on pasture-based dairy farms. The overall ARG relative abundance on these two farms was low, which provides further evidence for consumers of the safety of New Zealand's export products. Effluent samples harbored the most diverse range of ARGs, some of which were classified with a recognized risk to public health, whereas soil samples had the highest ARG relative abundance; however, the soil ARGs were not classified with a recognized public health risk. This emphasizes the need to consider genomic context and risk as well as ARG relative abundance in resistome studies.
  • Thumbnail Image
    Item
    Comparative Genomics of Mycobacterium avium Subspecies Paratuberculosis Sheep Strains
    (Frontiers Media S.A., 2021-02-15) Mizzi R; Timms VJ; Price-Carter ML; Gautam M; Whittington R; Heuer C; Biggs PJ; Plain KM; Salgado M
    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.
  • Thumbnail Image
    Item
    Complete Genome Sequences of Three Clostridiales R-7 Group Strains Isolated from the Bovine Rumen in New Zealand
    (American Society for Microbiology, 2021-07-01) Mahoney-Kurpe SC; Palevich N; Noel SJ; Kumar S; Gagic D; Biggs PJ; Janssen PH; Attwood GT; Moon CD; Putonti C
    Members of the Clostridiales R-7 group are abundant bacterial residents of the rumen microbiome; however, they are poorly characterized. We report the complete genome sequences of three members of the R-7 group, FE2010, FE2011, and XBB3002, isolated from the ruminal contents of pasture-grazed dairy cows in New Zealand.
  • Thumbnail Image
    Item
    Development of a multiplex droplet digital PCR assay for simultaneous detection and quantification of Escherichia coli, E. marmotae, and E. Ruysiae in water samples
    (Elsevier BV, 2024-03-01) Moinet M; Collis RM; Rogers L; Devane ML; Biggs PJ; Stott R; Marshall J; Muirhead R; Cookson AL
    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.
  • Thumbnail Image
    Item
    Draft genome sequences of Escherichia spp. isolates from New Zealand environmental sources.
    (American Society for Microbiology, 2024-03-12) Biggs PJ; Moinet M; Rogers LE; Devane M; Muirhead R; Stott R; Marshall JC; Cookson AL; Dennehy JJ
    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.
  • Thumbnail Image
    Item
    Extensive epigenetic modification with large-scale chromosomal and plasmid recombination characterise the Legionella longbeachae serogroup 1 genome
    (Springer Nature Limited, 2022-04-06) Slow S; Anderson T; Murdoch DR; Bloomfield S; Winter D; Biggs PJ
    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.
  • Thumbnail Image
    Item
    Genomic adaptations of Campylobacter jejuni to long-term human colonization
    (BioMed Central Ltd, 2021-12-10) Bloomfield SJ; Midwinter AC; Biggs PJ; French NP; Marshall JC; Hayman DTS; Carter PE; Mather AE; Fayaz A; Thornley C; Kelly DJ; Benschop J
    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.
  • Thumbnail Image
    Item
    Genomic analysis of canine pneumoviruses and canine respiratory coronavirus from New Zealand.
    (Taylor and Francis Group, 2024-07-01) Dunowska M; More GD; Biggs PJ; Cave NJ
    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.
  • Thumbnail Image
    Item
    Genomic Profiling of Mycobacterium tuberculosis Strains, Myanmar
    (Centers for Disease Control and Prevention, 2021-11) Aung HL; Nyunt WW; Fong Y; Biggs PJ; Winkworth RC; Lockhart PJ; Yeo TW; Hill PC; Cook GM; Aung ST
    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.
  • Thumbnail Image
    Item
    Insect Freeze-Tolerance Downunder: The Microbial Connection
    (MDPI (Basel, Switzerland), 2023-01-13) Morgan-Richards M; Marshall CJ; Biggs PJ; Trewick SA; Hoffmann KH
    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.
  • Thumbnail Image
    Item
    Investigating the genetic components of tuber bruising in a breeding population of tetraploid potatoes
    (BioMed Central Ltd, 2023-05-05) Angelin-Bonnet O; Thomson S; Vignes M; Biggs PJ; Monaghan K; Bloomer R; Wright K; Baldwin S
    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.
  • Thumbnail Image
    Item
    Lost in the Forest
    (Cold Spring Harbor Laboratory, 2022) Smith HL; Biggs PJ; French NP; Smith ANH; Marshall JC
    To date, there remains no satisfactory solution for absent levels in random forest models. Absent levels are levels of a predictor variable encountered during prediction for which no explicit rule exists. Imposing an order on nominal predictors allows absent levels to be integrated and used for prediction. The ordering of predictors has traditionally been via class probabilities with absent levels designated the lowest order. Using a combination of simulated data and pathogen source-attribution models using whole-genome sequencing data, we examine how the method of ordering predictors with absent levels can (i) systematically bias a model, and (ii) affect the out-of-bag error rate. We show that the traditional approach is systematically biased and underestimates out-of-bag error rates, and that this bias is resolved by ordering absent levels according to the a priori hypothesis of equal class probability. We present a novel method of ordering predictors via principal coordinates analysis (PCO) which capitalizes on the similarity between pairs of predictor levels. Absent levels are designated an order according to their similarity to each of the other levels in the training data. We show that the PCO method performs at least as well as the traditional approach of ordering and is not biased.
  • Thumbnail Image
    Item
    Lost in the Forest: Encoding categorical variables and the absent levels problem
    (Springer Nature, 2024-04-10) Smith HL; Biggs PJ; French NP; Smith ANH; Marshall JC; Gama J
    Levels of a predictor variable that are absent when a classification tree is grown can not be subject to an explicit splitting rule. This is an issue if these absent levels are present in a new observation for prediction. To date, there remains no satisfactory solution for absent levels in random forest models. Unlike missing data, absent levels are fully observed and known. Ordinal encoding of predictors allows absent levels to be integrated and used for prediction. Using a case study on source attribution of Campylobacter species using whole genome sequencing (WGS) data as predictors, we examine how target-agnostic versus target-based encoding of predictor variables with absent levels affects the accuracy of random forest models. We show that a target-based encoding approach using class probabilities, with absent levels designated the highest rank, is systematically biased, and that this bias is resolved by encoding absent levels according to the a priori hypothesis of equal class probability. We present a novel method of ordinal encoding predictors via principal coordinates analysis (PCO) which capitalizes on the similarity between pairs of predictor levels. Absent levels are encoded according to their similarity to each of the other levels in the training data. We show that the PCO-encoding method performs at least as well as the target-based approach and is not biased.
  • Thumbnail Image
    Item
    Mass Spectrometry-Based Lipidomics of Brown Adipose Tissue and Plasma of New-Born Lambs Subjected to Short-Term Cold Exposure
    (MDPI (Basel, Switzerland), 2022-10-14) Graña-Baumgartner A; Dukkipati VSR; Biggs PJ; Kenyon PR; Blair HT; Lopez-Villalobos N; Ross AB; Czauderna M
    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.
  • Thumbnail Image
    Item
    Microbial diversity and ecology of geothermal springs in the high-grade metamorphic terrain of Sri Lanka
    (Elsevier Ltd, 2022-04) Rupasinghe R; Amarasena S; Wickramarathna S; Biggs PJ; Chandrajith R; Wickramasinghe S
    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.
  • Thumbnail Image
    Item
    One dog's waste is another dog's wealth: A pilot study of fecal microbiota transplantation in dogs with acute hemorrhagic diarrhea syndrome
    (PLOS, 2021-04-19) Gal A; Barko PC; Biggs PJ; Gedye KR; Midwinter AC; Williams DA; Burchell RK; Pazzi P; Carbonero F
    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.