School of Mathematics and Statistics - Research Publications
Permanent URI for this collection
Search Results
1 - 10 of 11
-
ItemNo Preview AvailableMolecular profiling reveals features of clinical immunity and immunosuppression in asymptomatic P. falciparum malaria(WILEY, 2022-04-01)Clinical immunity to P. falciparum malaria is non-sterilizing, with adults often experiencing asymptomatic infection. Historically, asymptomatic malaria has been viewed as beneficial and required to help maintain clinical immunity. Emerging views suggest that these infections are detrimental and constitute a parasite reservoir that perpetuates transmission. To define the impact of asymptomatic malaria, we pursued a systems approach integrating antibody responses, mass cytometry, and transcriptional profiling of individuals experiencing symptomatic and asymptomatic P. falciparum infection. Defined populations of classical and atypical memory B cells and a TH2 cell bias were associated with reduced risk of clinical malaria. Despite these protective responses, asymptomatic malaria featured an immunosuppressive transcriptional signature with upregulation of pathways involved in the inhibition of T-cell function, and CTLA-4 as a predicted regulator in these processes. As proof of concept, we demonstrated a role for CTLA-4 in the development of asymptomatic parasitemia in infection models. The results suggest that asymptomatic malaria is not innocuous and might not support the induction of immune processes to fully control parasitemia or efficiently respond to malaria vaccines.
-
ItemNo Preview AvailableTargeting histone acetylation dynamics and oncogenic transcription by catalytic P300/CBP inhibition(CELL PRESS, 2021-05-20)To separate causal effects of histone acetylation on chromatin accessibility and transcriptional output, we used integrated epigenomic and transcriptomic analyses following acute inhibition of major cellular lysine acetyltransferases P300 and CBP in hematological malignancies. We found that catalytic P300/CBP inhibition dynamically perturbs steady-state acetylation kinetics and suppresses oncogenic transcriptional networks in the absence of changes to chromatin accessibility. CRISPR-Cas9 screening identified NCOR1 and HDAC3 transcriptional co-repressors as the principal antagonists of P300/CBP by counteracting acetylation turnover kinetics. Finally, deacetylation of H3K27 provides nucleation sites for reciprocal methylation switching, a feature that can be exploited therapeutically by concomitant KDM6A and P300/CBP inhibition. Overall, this study indicates that the steady-state histone acetylation-methylation equilibrium functions as a molecular rheostat governing cellular transcription that is amenable to therapeutic exploitation as an anti-cancer regimen.
-
ItemThe concerted change in the distribution of cell cycle phases and zone composition in germinal centers is regulated by IL-21(NATURE PORTFOLIO, 2021-12-09)Humoral immune responses require germinal centres (GC) for antibody affinity maturation. Within GC, B cell proliferation and mutation are segregated from affinity-based positive selection in the dark zone (DZ) and light zone (LZ) substructures, respectively. While IL-21 is known to be important in affinity maturation and GC maintenance, here we show it is required for both establishing normal zone representation and preventing the accumulation of cells in the G1 cell cycle stage in the GC LZ. Cell cycle progression of DZ B cells is unaffected by IL-21 availability, as is the zone phenotype of the most highly proliferative GC B cells. Collectively, this study characterises the development of GC zones as a function of time and B cell proliferation and identifies IL-21 as an important regulator of these processes. These data help explain the requirement for IL-21 in normal antibody affinity maturation.
-
ItemR code and downstream analysis objects for the scRNA-seq atlas of normal and tumorigenic human breast tissue(NATURE PORTFOLIO, 2022-03-23)Breast cancer is a common and highly heterogeneous disease. Understanding cellular diversity in the mammary gland and its surrounding micro-environment across different states can provide insight into cancer development in the human breast. Recently, we published a large-scale single-cell RNA expression atlas of the human breast spanning normal, preneoplastic and tumorigenic states. Single-cell expression profiles of nearly 430,000 cells were obtained from 69 distinct surgical tissue specimens from 55 patients. This article extends the study by providing quality filtering thresholds, downstream processed R data objects, complete cell annotation and R code to reproduce all the analyses. Data quality assessment measures are presented and details are provided for all the bioinformatic analyses that produced results described in the study.
-
ItemCalling differentially methylated regions from whole genome bisulphite sequencing with DMRcate(OXFORD UNIV PRESS, 2021-07-28)Whole genome bisulphite sequencing (WGBS) permits the genome-wide study of single molecule methylation patterns. One of the key goals of mammalian cell-type identity studies, in both normal differentiation and disease, is to locate differential methylation patterns across the genome. We discuss the most desirable characteristics for DML (differentially methylated locus) and DMR (differentially methylated region) detection tools in a genome-wide context and choose a set of statistical methods that fully or partially satisfy these considerations to compare for benchmarking. Our data simulation strategy is both biologically informed-employing distribution parameters derived from large-scale consortium datasets-and thorough. We report DML detection ability with respect to coverage, group methylation difference, sample size, variability and covariate size, both marginally and jointly, and exhaustively with respect to parameter combination. We also benchmark these methods on FDR control and computational time. We use this result to backend and introduce an expanded version of DMRcate: an existing DMR detection tool for microarray data that we have extended to now call DMRs from WGBS data. We compare DMRcate to a set of alternative DMR callers using a similarly realistic simulation strategy. We find DMRcate and RADmeth are the best predictors of DMRs, and conclusively find DMRcate the fastest.
-
ItemIn vivo genome-editing screen identifies tumor suppressor genes that cooperate with Trp53 loss during mammary tumorigenesis(WILEY, 2022-01-26)Breast cancer is a heterogeneous disease that comprises multiple histological and molecular subtypes. To gain insight into mutations that drive breast tumorigenesis, we describe a pipeline for the identification and validation of tumor suppressor genes. Based on an in vivo genome-wide CRISPR/Cas9 screen in Trp53+/- heterozygous mice, we identified tumor suppressor genes that included the scaffold protein Axin1, the protein kinase A regulatory subunit gene Prkar1a, as well as the proof-of-concept genes Pten, Nf1, and Trp53 itself. Ex vivo editing of primary mammary epithelial organoids was performed to further interrogate the roles of Axin1 and Prkar1a. Increased proliferation and profound changes in mammary organoid morphology were observed for Axin1/Trp53 and Prkar1a/Trp53 double mutants compared to Pten/Trp53 double mutants. Furthermore, direct in vivo genome editing via intraductal injection of lentiviruses engineered to express dual short-guide RNAs revealed that mutagenesis of Trp53 and either Prkar1a, Axin1, or Pten markedly accelerated tumor development compared to Trp53-only mutants. This proof-of-principle study highlights the application of in vivo CRISPR/Cas9 editing for uncovering cooperativity between defects in tumor suppressor genes that elicit mammary tumorigenesis.
-
ItemSingle cell transcriptome atlas of mouse mammary epithelial cells across development(BMC, 2021-06-29)BACKGROUND: Heterogeneity within the mouse mammary epithelium and potential lineage relationships have been recently explored by single-cell RNA profiling. To further understand how cellular diversity changes during mammary ontogeny, we profiled single cells from nine different developmental stages spanning late embryogenesis, early postnatal, prepuberty, adult, mid-pregnancy, late-pregnancy, and post-involution, as well as the transcriptomes of micro-dissected terminal end buds (TEBs) and subtending ducts during puberty. METHODS: The single cell transcriptomes of 132,599 mammary epithelial cells from 9 different developmental stages were determined on the 10x Genomics Chromium platform, and integrative analyses were performed to compare specific time points. RESULTS: The mammary rudiment at E18.5 closely aligned with the basal lineage, while prepubertal epithelial cells exhibited lineage segregation but to a less differentiated state than their adult counterparts. Comparison of micro-dissected TEBs versus ducts showed that luminal cells within TEBs harbored intermediate expression profiles. Ductal basal cells exhibited increased chromatin accessibility of luminal genes compared to their TEB counterparts suggesting that lineage-specific chromatin is established within the subtending ducts during puberty. An integrative analysis of five stages spanning the pregnancy cycle revealed distinct stage-specific profiles and the presence of cycling basal, mixed-lineage, and 'late' alveolar intermediates in pregnancy. Moreover, a number of intermediates were uncovered along the basal-luminal progenitor cell axis, suggesting a continuum of alveolar-restricted progenitor states. CONCLUSIONS: This extended single cell transcriptome atlas of mouse mammary epithelial cells provides the most complete coverage for mammary epithelial cells during morphogenesis to date. Together with chromatin accessibility analysis of TEB structures, it represents a valuable framework for understanding developmental decisions within the mouse mammary gland.
-
ItemA guide to creating design matrices for gene expression experiments.(F1000 Research Ltd, 2020)Differential expression analysis of genomic data types, such as RNA-sequencing experiments, use linear models to determine the size and direction of the changes in gene expression. For RNA-sequencing, there are several established software packages for this purpose accompanied with analysis pipelines that are well described. However, there are two crucial steps in the analysis process that can be a stumbling block for many -- the set up an appropriate model via design matrices and the set up of comparisons of interest via contrast matrices. These steps are particularly troublesome because an extensive catalogue for design and contrast matrices does not currently exist. One would usually search for example case studies across different platforms and mix and match the advice from those sources to suit the dataset they have at hand. This article guides the reader through the basics of how to set up design and contrast matrices. We take a practical approach by providing code and graphical representation of each case study, starting with simpler examples (e.g. models with a single explanatory variable) and move onto more complex ones (e.g. interaction models, mixed effects models, higher order time series and cyclical models). Although our work has been written specifically with a limma-style pipeline in mind, most of it is also applicable to other software packages for differential expression analysis, and the ideas covered can be adapted to data analysis of other high-throughput technologies. Where appropriate, we explain the interpretation and differences between models to aid readers in their own model choices. Unnecessary jargon and theory is omitted where possible so that our work is accessible to a wide audience of readers, from beginners to those with experience in genomics data analysis.
-
ItemNo Preview AvailableGermline heterozygous mutations in Nxf1 perturb RNA metabolism and trigger thrombocytopenia and lymphopenia in mice(ELSEVIER, 2020-04-14)In eukaryotic cells, messenger RNA (mRNA) molecules are exported from the nucleus to the cytoplasm, where they are translated. The highly conserved protein nuclear RNA export factor1 (Nxf1) is an important mediator of this process. Although studies in yeast and in human cell lines have shed light on the biochemical mechanisms of Nxf1 function, its contribution to mammalian physiology is less clear. Several groups have identified recurrent NXF1 mutations in chronic lymphocytic leukemia (CLL), placing it alongside several RNA-metabolism factors (including SF3B1, XPO, RPS15) whose dysregulation is thought to contribute to CLL pathogenesis. We report here an allelic series of germline point mutations in murine Nxf1. Mice heterozygous for these loss-of-function Nxf1 mutations exhibit thrombocytopenia and lymphopenia, together with milder hematological defects. This is primarily caused by cell-intrinsic defects in the survival of platelets and peripheral lymphocytes, which are sensitized to intrinsic apoptosis. In contrast, Nxf1 mutations have almost no effect on red blood cell homeostasis. Comparative transcriptome analysis of platelets, lymphocytes, and erythrocytes from Nxf1-mutant mice shows that, in response to impaired Nxf1 function, the cytoplasmic representation of transcripts encoding regulators of RNA metabolism is altered in a unique, lineage-specific way. Thus, blood cell lineages exhibit differential requirements for Nxf1-mediated global mRNA export.
-
ItemTargeting triple-negative breast cancers with the Smac-mimetic birinapant(Springer Nature, 2020-04-27)Smac mimetics target inhibitor of apoptosis (IAP) proteins, thereby suppressing their function to facilitate tumor cell death. Here we have evaluated the efficacy of the preclinical Smac-mimetic compound A and the clinical lead birinapant on breast cancer cells. Both exhibited potent in vitro activity in triple-negative breast cancer (TNBC) cells, including those from patient-derived xenograft (PDX) models. Birinapant was further studied using in vivo PDX models of TNBC and estrogen receptor-positive (ER+) breast cancer. Birinapant exhibited single agent activity in all TNBC PDX models and augmented response to docetaxel, the latter through induction of TNF. Transcriptomic analysis of TCGA datasets revealed that genes encoding mediators of Smac-mimetic-induced cell death were expressed at higher levels in TNBC compared with ER+ breast cancer, resulting in a molecular signature associated with responsiveness to Smac mimetics. In addition, the cell death complex was preferentially formed in TNBCs versus ER+ cells in response to Smac mimetics. Taken together, our findings provide a rationale for prospectively selecting patients whose breast tumors contain a competent death receptor signaling pathway for the further evaluation of birinapant in the clinic.