School of Mathematics and Statistics - Research Publications

Permanent URI for this collection

Search Results

Now showing 1 - 10 of 64
  • Item
    No Preview Available
    Metaphor-A workflow for streamlined assembly and binning of metagenomes.
    Salazar, VW ; Shaban, B ; Quiroga, MDM ; Turnbull, R ; Tescari, E ; Rossetto Marcelino, V ; Verbruggen, H ; Lê Cao, K-A (Oxford University Press (OUP), 2022-12-28)
    Recent advances in bioinformatics and high-throughput sequencing have enabled the large-scale recovery of genomes from metagenomes. This has the potential to bring important insights as researchers can bypass cultivation and analyze genomes sourced directly from environmental samples. There are, however, technical challenges associated with this process, most notably the complexity of computational workflows required to process metagenomic data, which include dozens of bioinformatics software tools, each with their own set of customizable parameters that affect the final output of the workflow. At the core of these workflows are the processes of assembly-combining the short-input reads into longer, contiguous fragments (contigs)-and binning, clustering these contigs into individual genome bins. The limitations of assembly and binning algorithms also pose different challenges depending on the selected strategy to execute them. Both of these processes can be done for each sample separately or by pooling together multiple samples to leverage information from a combination of samples. Here we present Metaphor, a fully automated workflow for genome-resolved metagenomics (GRM). Metaphor differs from existing GRM workflows by offering flexible approaches for the assembly and binning of the input data and by combining multiple binning algorithms with a bin refinement step to achieve high-quality genome bins. Moreover, Metaphor generates reports to evaluate the performance of the workflow. We showcase the functionality of Metaphor on different synthetic datasets and the impact of available assembly and binning strategies on the final results.
  • Item
    Thumbnail Image
    IL11 activates the placental inflammasome to drive preeclampsia
    Menkhorst, E ; Santos, LL ; Zhou, W ; Yang, G ; Winship, AL ; Rainczuk, KE ; Nguyen, P ; Zhang, J-G ; Moore, P ; Williams, M ; Le Cao, K-A ; Mansell, A ; Dimitriadis, E (FRONTIERS MEDIA SA, 2023-05-24)
    INTRODUCTION: Preeclampsia is a life-threatening disorder of pregnancy unique to humans. Interleukin (IL)11 is elevated in serum from pregnancies that subsequently develop early-onset preeclampsia and pharmacological elevation of IL11 in pregnant mice causes the development of early-onset preeclampsia-like features (hypertension, proteinuria, and fetal growth restriction). However, the mechanism by which IL11 drives preeclampsia is unknown. METHOD: Pregnant mice were administered PEGylated (PEG)IL11 or control (PEG) from embryonic day (E)10-16 and the effect on inflammasome activation, systolic blood pressure (during gestation and at 50/90 days post-natal), placental development, and fetal/post-natal pup growth measured. RNAseq analysis was performed on E13 placenta. Human 1st trimester placental villi were treated with IL11 and the effect on inflammasome activation and pyroptosis identified by immunohistochemistry and ELISA. RESULT: PEGIL11 activated the placental inflammasome causing inflammation, fibrosis, and acute and chronic hypertension in wild-type mice. Global and placental-specific loss of the inflammasome adaptor protein Asc and global loss of the Nlrp3 sensor protein prevented PEGIL11-induced fibrosis and hypertension in mice but did not prevent PEGIL11-induced fetal growth restriction or stillbirths. RNA-sequencing and histology identified that PEGIL11 inhibited trophoblast differentiation towards spongiotrophoblast and syncytiotrophoblast lineages in mice and extravillous trophoblast lineages in human placental villi. DISCUSSION: Inhibition of ASC/NLRP3 inflammasome activity could prevent IL11-induced inflammation and fibrosis in various disease states including preeclampsia.
  • Item
    Thumbnail Image
    Randomized phase I trial of antigen- specific tolerizing immunotherapy with peptide/calcitriol liposomes in ACPA plus rheumatoid arthritis
    Sonigra, A ; Nel, HJ ; Wehr, P ; Ramnoruth, N ; Patel, S ; van Schie, KA ; Bladen, MW ; Mehdi, AM ; Tesiram, J ; Talekar, M ; Rossjohn, J ; Reid, HH ; Stuurman, FE ; Roberts, H ; Vecchio, P ; Gourley, I ; Rigby, M ; Becart, S ; Toes, RE ; Scherer, HU ; Cao, K-AL ; Campbell, K ; Thomas, R (AMER SOC CLINICAL INVESTIGATION INC, 2022-10-24)
    BACKGROUNDAntigen-specific regulation of autoimmune disease is a major goal. In seropositive rheumatoid arthritis (RA), T cell help to autoreactive B cells matures the citrullinated (Cit) antigen-specific immune response, generating RA-specific V domain glycosylated anti-Cit protein antibodies (ACPA VDG) before arthritis onset. Low or escalating antigen administration under "sub-immunogenic" conditions favors tolerance. We explored safety, pharmacokinetics, and immunological and clinical effects of s.c. DEN-181, comprising liposomes encapsulating self-peptide collagen II259-273 (CII) and NF-κB inhibitor 1,25-dihydroxycholecalciferol.METHODSA double-blind, placebo-controlled, exploratory, single-ascending-dose, phase I trial assessed the impact of low, medium, and high DEN-181 doses on peripheral blood CII-specific and bystander Cit64vimentin59-71-specific (Cit-Vim-specific) autoreactive T cell responses, cytokines, and ACPA in 17 HLA-DRB1*04:01+ or *01:01+ ACPA+ RA patients on methotrexate.RESULTSDEN-181 was well tolerated. Relative to placebo and normalized to baseline values, Cit-Vim-specific T cells decreased in patients administered medium and high doses of DEN-181. Relative to placebo, percentage of CII-specific programmed cell death 1+ T cells increased within 28 days of DEN-181. Exploratory analysis in DEN-181-treated patients suggested improved RA disease activity was associated with expansion of CII-specific and Cit-Vim-specific T cells; reduction in ACPA VDG, memory B cells, and inflammatory myeloid populations; and enrichment in CCR7+ and naive T cells. Single-cell sequencing identified T cell transcripts associated with tolerogenic TCR signaling and exhaustion after low or medium doses of DEN-181.CONCLUSIONThe safety and immunomodulatory activity of low/medium DEN-181 doses provide rationale to further assess antigen-specific immunomodulatory therapy in ACPA+ RA.TRIAL REGISTRATIONAnzctr.org.au identifier ACTRN12617001482358, updated September 8, 2022.FUNDINGInnovative Medicines Initiative 2 Joint Undertaking (grant agreement 777357), supported by European Union's Horizon 2020 research and innovation programme and European Federation of Pharmaceutical Industries and Associations; Arthritis Queensland; National Health and Medical Research Council (NHMRC) Senior Research Fellowship; and NHMRC grant 2008287.
  • Item
    Thumbnail Image
    Protocol Microbiota DNA isolation, 16S rRNA amplicon sequencing, and bioinformatic analysis for bacterial microbiome profiling of rodent fecal samples
    Love, CJ ; Gubert, C ; Kodikara, S ; Kong, G ; Cao, K-AL ; Hannan, AJ (ELSEVIER, 2022-12-16)
    Fecal samples are frequently used to characterize bacterial populations of the gastrointestinal tract. A protocol is provided to profile gut bacterial populations using rodent fecal samples. We describe the optimal procedures for collecting rodent fecal samples, isolating genomic DNA, 16S rRNA gene V4 region sequencing, and bioinformatic analyses. This protocol includes detailed instructions and example outputs to ensure accurate, reproducible results and data visualization. Comprehensive troubleshooting and limitation sections address technical and statistical issues that may arise when profiling microbiota. For complete details on the use and execution of this protocol, please refer to Gubert et al. (2022).
  • Item
    No Preview Available
    PLSDA-batch: a multivariate framework to correct for batch effects in microbiome data
    Wang, Y ; Cao, K-AL (OXFORD UNIV PRESS, 2023-03-19)
    Microbial communities are highly dynamic and sensitive to changes in the environment. Thus, microbiome data are highly susceptible to batch effects, defined as sources of unwanted variation that are not related to and obscure any factors of interest. Existing batch effect correction methods have been primarily developed for gene expression data. As such, they do not consider the inherent characteristics of microbiome data, including zero inflation, overdispersion and correlation between variables. We introduce new multivariate and non-parametric batch effect correction methods based on Partial Least Squares Discriminant Analysis (PLSDA). PLSDA-batch first estimates treatment and batch variation with latent components, then subtracts batch-associated components from the data. The resulting batch-effect-corrected data can then be input in any downstream statistical analysis. Two variants are proposed to handle unbalanced batch x treatment designs and to avoid overfitting when estimating the components via variable selection. We compare our approaches with popular methods managing batch effects, namely, removeBatchEffect, ComBat and Surrogate Variable Analysis, in simulated and three case studies using various visual and numerical assessments. We show that our three methods lead to competitive performance in removing batch variation while preserving treatment variation, especially for unbalanced batch $\times $ treatment designs. Our downstream analyses show selections of biologically relevant taxa. This work demonstrates that batch effect correction methods can improve microbiome research outputs. Reproducible code and vignettes are available on GitHub.
  • Item
    Thumbnail Image
    MetaGenePipe: An Automated, Portable Pipeline for Contig-based Functional and Taxonomic Analysis
    Shaban, B ; Quiroga, M ; Turnbull, R ; Tescari, E ; Lê Cao, K-A ; Verbruggen, H (The Open Journal, 2023)
    MetaGenePipe (MGP) is an efficient, flexible, portable, and scalable metagenomics pipeline that uses performant bioinformatics software suites and genomic databases to create an accurate taxonomic and functional characterization of the prokaryotic fraction of sequenced microbiomes. Written in the Workflow Definition Language (WDL), MGP produces output that can be explored and interpreted directly, or can be used for downstream analysis. MGP is a pipeline-development best practice tool that uses Singularity for containerization and includes a setup script that downloads the necessary databases for setup. The source code for MGP is freely available and distributed under the Apache 2.0 license.
  • Item
    Thumbnail Image
    Faecal microbiota transplant ameliorates gut dysbiosis and cognitive deficits in Huntington's disease mice
    Gubert, C ; Choo, JM ; Love, CJ ; Kodikara, S ; Masson, BA ; Liew, JJM ; Wang, Y ; Kong, G ; Narayana, VK ; Renoir, T ; Cao, K-AL ; Rogers, GB ; Hannan, AJ (OXFORD UNIV PRESS, 2022-07-04)
    Huntington's disease is a neurodegenerative disorder involving psychiatric, cognitive and motor symptoms. Huntington's disease is caused by a tandem-repeat expansion in the huntingtin gene, which is widely expressed throughout the brain and body, including the gastrointestinal system. There are currently no effective disease-modifying treatments available for this fatal disorder. Despite recent evidence of gut microbiome disruption in preclinical and clinical Huntington's disease, its potential as a target for therapeutic interventions has not been explored. The microbiota-gut-brain axis provides a potential pathway through which changes in the gut could modulate brain function, including cognition. We now show that faecal microbiota transplant (FMT) from wild-type into Huntington's disease mice positively modulates cognitive outcomes, particularly in females. In Huntington's disease male mice, we revealed an inefficiency of FMT engraftment, which is potentially due to the more pronounced changes in the structure, composition and instability of the gut microbial community, and the imbalance in acetate and gut immune profiles found in these mice. This study demonstrates a role for gut microbiome modulation in ameliorating cognitive deficits modelling dementia in Huntington's disease. Our findings pave the way for the development of future therapeutic approaches, including FMT and other forms of gut microbiome modulation, as potential clinical interventions for Huntington's disease.
  • Item
    Thumbnail Image
    Statistical challenges in longitudinal microbiome data analysis
    Kodikara, S ; Ellul, S ; Le Cao, K-A (OXFORD UNIV PRESS, 2022-07-18)
    The microbiome is a complex and dynamic community of microorganisms that co-exist interdependently within an ecosystem, and interact with its host or environment. Longitudinal studies can capture temporal variation within the microbiome to gain mechanistic insights into microbial systems; however, current statistical methods are limited due to the complex and inherent features of the data. We have identified three analytical objectives in longitudinal microbial studies: (1) differential abundance over time and between sample groups, demographic factors or clinical variables of interest; (2) clustering of microorganisms evolving concomitantly across time and (3) network modelling to identify temporal relationships between microorganisms. This review explores the strengths and limitations of current methods to fulfill these objectives, compares different methods in simulation and case studies for objectives (1) and (2), and highlights opportunities for further methodological developments. R tutorials are provided to reproduce the analyses conducted in this review.
  • Item
    Thumbnail Image
    Sincast: a computational framework to predict cell identities in single-cell transcriptomes using bulk atlases as references
    Deng, Y ; Choi, J ; Cao, K-AL (OXFORD UNIV PRESS, 2022-05-13)
    Characterizing the molecular identity of a cell is an essential step in single-cell RNA sequencing (scRNA-seq) data analysis. Numerous tools exist for predicting cell identity using single-cell reference atlases. However, many challenges remain, including correcting for inherent batch effects between reference and query data andinsufficient phenotype data from the reference. One solution is to project single-cell data onto established bulk reference atlases to leverage their rich phenotype information. Sincast is a computational framework to query scRNA-seq data by projection onto bulk reference atlases. Prior to projection, single-cell data are transformed to be directly comparable to bulk data, either with pseudo-bulk aggregation or graph-based imputation to address sparse single-cell expression profiles. Sincast avoids batch effect correction, and cell identity is predicted along a continuum to highlight new cell states not found in the reference atlas. In several case study scenarios, we show that Sincast projects single cells into the correct biological niches in the expression space of the bulk reference atlas. We demonstrate the effectiveness of our imputation approach that was specifically developed for querying scRNA-seq data based on bulk reference atlases. We show that Sincast is an efficient and powerful tool for single-cell profiling that will facilitate downstream analysis of scRNA-seq data.
  • Item
    Thumbnail Image
    Host Traits and Phylogeny Contribute to Shaping Coral-Bacterial Symbioses
    Ricci, F ; Tandon, K ; Black, JR ; Cao, K-AL ; Blackall, LL ; Verbruggen, H ; Raina, J-B (AMER SOC MICROBIOLOGY, 2022-04-26)
    The success of tropical scleractinian corals depends on their ability to establish symbioses with microbial partners. Host phylogeny and traits are known to shape the coral microbiome, but to what extent they affect its composition remains unclear. Here, by using 12 coral species representing the complex and robust clades, we explored the influence of host phylogeny, skeletal architecture, and reproductive mode on the microbiome composition, and further investigated the structure of the tissue and skeleton bacterial communities. Our results show that host phylogeny and traits explained 14% of the tissue and 13% of the skeletal microbiome composition, providing evidence that these predictors contributed to shaping the holobiont in terms of presence and relative abundance of bacterial symbionts. Based on our data, we conclude that host phylogeny affects the presence of specific microbial lineages, reproductive mode predictably influences the microbiome composition, and skeletal architecture works like a filter that affects bacterial relative abundance. We show that the β-diversity of coral tissue and skeleton microbiomes differed, but we found that a large overlapping fraction of bacterial sequences were recovered from both anatomical compartments, supporting the hypothesis that the skeleton can function as a microbial reservoir. Additionally, our analysis of the microbiome structure shows that 99.6% of tissue and 99.7% of skeletal amplicon sequence variants (ASVs) were not consistently present in at least 30% of the samples, suggesting that the coral tissue and skeleton are dominated by rare bacteria. Together, these results provide novel insights into the processes driving coral-bacterial symbioses, along with an improved understanding of the scleractinian microbiome.