Aaron Newman
Associate Professor of Biomedical Data Science
Department of Biomedical Data Science
Web page: https://anlab.stanford.edu/
Academic Appointments
-
Associate Professor, Department of Biomedical Data Science
-
Member, Bio-X
-
Member, Stanford Cancer Institute
Honors & Awards
-
3rd Annual Daisy Ho Memorial Award, Department of Experimental Therapeutics, MD Anderson Cancer Center (2023)
-
CZ Biohub Investigator, Chan Zuckerberg Biohub (2022 - 2027)
-
Baxter Faculty Scholar Award, Donald E. and Delia B. Baxter Foundation (2020)
-
Interdisciplinary Initiatives Seed Grants Program Award, Stanford Bio-X (2018 - 2020)
-
Stem Cell Research Award, Stinehart/Reed Foundation (2018 - 2020)
-
K99/R00 Pathway to Independence Award, NIH/NCI (2015 - 2020)
-
Visionary Postdoctoral Fellowship, Dept. of Defense (2012 - 2015)
-
NIH T32 Cancer Biology Training Grant, Stanford University (2012)
-
Siebel Fellow, The Siebel Stem Cell Institute (Stanford/UC Berkeley) (2011-2015)
-
Graduate Dissertation Fellowship Award, UCSB (2009)
Boards, Advisory Committees, Professional Organizations
-
Member, International Society for Computational Biology (ISCB) (2008 - Present)
-
Associate Member, American Association for Cancer Research (AACR) (2012 - Present)
Professional Education
-
PhD, University of California, Santa Barbara, Biomolecular Science and Engineering Program (2010)
Current Research and Scholarly Interests
Our group combines computational and experimental techniques to study the cellular organization of complex tissues, with a focus on determining the phenotypic diversity and clinical significance of tumor cell subsets and their surrounding microenvironments. We have a particular interest in developing innovative data science tools that illuminate the cellular hierarchies and stromal elements that underlie tumor initiation, progression, and response to therapy. As part of this focus, we develop new algorithms to resolve cellular states and multicellular communities, tumor developmental hierarchies, and single-cell spatial relationships from genomic profiles of clinical biospecimens. Key results are further explored experimentally, both in our lab and through collaboration, with the goal of translating promising findings into the clinic.
As a member of the Department of Biomedical Data Science and the Institute for Stem Cell Biology and Regenerative Medicine, and as an affiliate of graduate programs in Biomedical Informatics, Cancer Biology, and Immunology, we are also interested in the development of impactful biomedical data science tools in areas beyond our immediate research focus, including developmental biology, regenerative medicine, and systems immunology.
2024-25 Courses
- Bioinformatics for Stem Cell and Cancer Biology
BIODS 205, STEMREM 205 (Win) - Workshop in Biostatistics
BIODS 260B, STATS 260B (Win) -
Independent Studies (8)
- Directed Investigation
BIOE 392 (Aut, Win, Spr, Sum) - Directed Study
BIOE 391 (Aut, Win, Spr, Sum) - Graduate Research
CBIO 399 (Aut, Win, Spr, Sum) - Graduate Research
IMMUNOL 399 (Aut, Win, Spr, Sum) - Graduate Research
STEMREM 399 (Aut, Win, Spr, Sum) - Medical Scholars Research
MED 370 (Aut, Sum) - Medical Scholars Research
STEMREM 370 (Aut, Win, Spr, Sum) - Undergraduate Research
STEMREM 199 (Aut, Win, Spr, Sum)
- Directed Investigation
-
Prior Year Courses
2023-24 Courses
- Bioinformatics for Stem Cell and Cancer Biology
BIODS 205, STEMREM 205 (Win) - Workshop in Biostatistics
BIODS 260B, STATS 260B (Win)
2022-23 Courses
- Bioinformatics for Stem Cell and Cancer Biology
BIODS 205, STEMREM 205 (Win) - Workshop in Biostatistics
BIODS 260B, STATS 260B (Win)
2021-22 Courses
- Bioinformatics for Stem Cell and Cancer Biology
BIODS 205, STEMREM 205 (Win) - Workshop in Biostatistics
BIODS 260B, STATS 260B (Win)
- Bioinformatics for Stem Cell and Cancer Biology
Stanford Advisees
-
Med Scholar Project Advisor
Samuel Ho -
Doctoral Dissertation Reader (AC)
Jessica Arozqueta Basurto, Quenton Bubb -
Postdoctoral Faculty Sponsor
Wubing Zhang -
Doctoral Dissertation Advisor (AC)
Jeremy D'Silva, Rachel Gleyzer, Sharon Loa, Janella Schwab Lizarraga -
Doctoral (Program)
Minji Kang -
Postdoctoral Research Mentor
Wubing Zhang
Graduate and Fellowship Programs
-
Biomedical Data Science (Phd Program)
All Publications
-
Profiling cell identity and tissue architecture with single-cell and spatial transcriptomics.
Nature reviews. Molecular cell biology
2024
Abstract
Single-cell transcriptomics has broadened our understanding of cellular diversity and gene expression dynamics in healthy and diseased tissues. Recently, spatial transcriptomics has emerged as a tool to contextualize single cells in multicellular neighbourhoods and to identify spatially recurrent phenotypes, or ecotypes. These technologies have generated vast datasets with targeted-transcriptome and whole-transcriptome profiles of hundreds to millions of cells. Such data have provided new insights into developmental hierarchies, cellular plasticity and diverse tissue microenvironments, and spurred a burst of innovation in computational methods for single-cell analysis. In this Review, we discuss recent advancements, ongoing challenges and prospects in identifying and characterizing cell states and multicellular neighbourhoods. We discuss recent progress in sample processing, data integration, identification of subtle cell states, trajectory modelling, deconvolution and spatial analysis. Furthermore, we discuss the increasing application of deep learning, including foundation models, in analysing single-cell and spatial transcriptomics data. Finally, we discuss recent applications of these tools in the fields of stem cell biology, immunology, and tumour biology, and the future of single-cell and spatial transcriptomics in biological research and its translation to the clinic.
View details for DOI 10.1038/s41580-024-00768-2
View details for PubMedID 39169166
View details for PubMedCentralID 9888017
-
High-resolution alignment of single-cell and spatial transcriptomes with CytoSPACE.
Nature biotechnology
2023
Abstract
Recent studies have emphasized the importance of single-cell spatial biology, yet available assays for spatial transcriptomics have limited gene recovery or low spatial resolution. Here we introduce CytoSPACE, an optimization method for mapping individual cells from a single-cell RNA sequencing atlas to spatial expression profiles. Across diverse platforms and tissue types, we show that CytoSPACE outperforms previous methods with respect to noise tolerance and accuracy, enabling tissue cartography at single-cell resolution.
View details for DOI 10.1038/s41587-023-01697-9
View details for PubMedID 36879008
View details for PubMedCentralID 6132072
-
T cell characteristics associated with toxicity to immune checkpoint blockade in patients with melanoma.
Nature medicine
2022
Abstract
Severe immune-related adverse events (irAEs) occur in up to 60% of patients with melanoma treated with immune checkpoint inhibitors (ICIs). However, it is unknown whether a common baseline immunological state precedes irAE development. Here we applied mass cytometry by time of flight, single-cell RNA sequencing, single-cell V(D)J sequencing, bulk RNA sequencing and bulk T cell receptor (TCR) sequencing to study peripheral blood samples from patients with melanoma treated with anti-PD-1 monotherapy or anti-PD-1 and anti-CTLA-4 combination ICIs. By analyzing 93 pre- and early on-ICI blood samples and 3 patient cohorts (n = 27, 26 and 18), we found that 2 pretreatment factors in circulation-activated CD4 memory T cell abundance and TCR diversity-are associated with severe irAE development regardless of organ system involvement. We also explored on-treatment changes in TCR clonality among patients receiving combination therapy and linked our findings to the severity and timing of irAE onset. These results demonstrate circulating T cell characteristics associated with ICI-induced toxicity, with implications for improved diagnostics and clinical management.
View details for DOI 10.1038/s41591-021-01623-z
View details for PubMedID 35027754
-
Atlas of clinically distinct cell states and ecosystems across human solid tumors.
Cell
2021
Abstract
Determining how cells vary with their local signaling environment and organize into distinct cellular communities is critical for understanding processes as diverse as development, aging, and cancer. Here we introduce EcoTyper, a machine learning framework for large-scale identification and validation of cell states and multicellular communities from bulk, single-cell, and spatially resolved gene expression data. When applied to 12 major cell lineages across 16 types of human carcinoma, EcoTyper identified 69 transcriptionally defined cell states. Most states were specific to neoplastic tissue, ubiquitous across tumor types, and significantly prognostic. By analyzing cell-state co-occurrence patterns, we discovered ten clinically distinct multicellular communities with unexpectedly strong conservation, including three with myeloid and stromal elements linked to adverse survival, one enriched in normal tissue, and two associated with early cancer development. This study elucidates fundamental units of cellular organization in human carcinoma and provides a framework for large-scale profiling of cellular ecosystems in any tissue.
View details for DOI 10.1016/j.cell.2021.09.014
View details for PubMedID 34597583
-
Single-cell transcriptional diversity is a hallmark of developmental potential.
Science (New York, N.Y.)
2020; 367 (6476): 405–11
Abstract
Single-cell RNA sequencing (scRNA-seq) is a powerful approach for reconstructing cellular differentiation trajectories. However, inferring both the state and direction of differentiation is challenging. Here, we demonstrate a simple, yet robust, determinant of developmental potential-the number of expressed genes per cell-and leverage this measure of transcriptional diversity to develop a computational framework (CytoTRACE) for predicting differentiation states from scRNA-seq data. When applied to diverse tissue types and organisms, CytoTRACE outperformed previous methods and nearly 19,000 annotated gene sets for resolving 52 experimentally determined developmental trajectories. Additionally, it facilitated the identification of quiescent stem cells and revealed genes that contribute to breast tumorigenesis. This study thus establishes a key RNA-based feature of developmental potential and a platform for delineation of cellular hierarchies.
View details for DOI 10.1126/science.aax0249
View details for PubMedID 31974247
-
Determining cell type abundance and expression from bulk tissues with digital cytometry.
Nature biotechnology
2019
Abstract
Single-cell RNA-sequencing has emerged as a powerful technique for characterizing cellular heterogeneity, but it is currently impractical on large sample cohorts and cannot be applied to fixed specimens collected as part of routine clinical care. We previously developed an approach for digital cytometry, called CIBERSORT, that enables estimation of cell type abundances from bulk tissue transcriptomes. We now introduce CIBERSORTx, a machine learning method that extends this framework to infer cell-type-specific gene expression profiles without physical cell isolation. By minimizing platform-specific variation, CIBERSORTx also allows the use of single-cell RNA-sequencing data for large-scale tissue dissection. We evaluated the utility of CIBERSORTx in multiple tumor types, including melanoma, where single-cell reference profiles were used to dissect bulk clinical specimens, revealing cell-type-specific phenotypic states linked to distinct driver mutations and response to immune checkpoint blockade. We anticipate that digital cytometry will augment single-cell profiling efforts, enabling cost-effective, high-throughput tissue characterization without the need for antibodies, disaggregation or viable cells.
View details for PubMedID 31061481
-
Integrated digital error suppression for improved detection of circulating tumor DNA
NATURE BIOTECHNOLOGY
2016
View details for DOI 10.1038/nbt.3520
-
Robust enumeration of cell subsets from tissue expression profiles
NATURE METHODS
2015
View details for DOI 10.1038/nmeth.3337
-
The prognostic landscape of genes and infiltrating immune cells across human cancers
NATURE MEDICINE
2015
View details for DOI 10.1038/nm.3909
-
An ultrasensitive method for quantitating circulating tumor DNA with broad patient coverage
NATURE MEDICINE
2014; 20 (5): 552-558
Abstract
Circulating tumor DNA (ctDNA) is a promising biomarker for noninvasive assessment of cancer burden, but existing ctDNA detection methods have insufficient sensitivity or patient coverage for broad clinical applicability. Here we introduce cancer personalized profiling by deep sequencing (CAPP-Seq), an economical and ultrasensitive method for quantifying ctDNA. We implemented CAPP-Seq for non-small-cell lung cancer (NSCLC) with a design covering multiple classes of somatic alterations that identified mutations in >95% of tumors. We detected ctDNA in 100% of patients with stage II-IV NSCLC and in 50% of patients with stage I, with 96% specificity for mutant allele fractions down to ∼0.02%. Levels of ctDNA were highly correlated with tumor volume and distinguished between residual disease and treatment-related imaging changes, and measurement of ctDNA levels allowed for earlier response assessment than radiographic approaches. Finally, we evaluated biopsy-free tumor screening and genotyping with CAPP-Seq. We envision that CAPP-Seq could be routinely applied clinically to detect and monitor diverse malignancies, thus facilitating personalized cancer therapy.
View details for DOI 10.1038/nm.3519
View details for Web of Science ID 000335710700028
-
Identification of a colonial chordate histocompatibility gene
SCIENCE
2013; 341 (6144): 384-387
View details for DOI 10.1126/science.1238036
-
Lab-Specific Gene Expression Signatures in Pluripotent Stem Cells
CELL STEM CELL
2010; 7 (2): 258-262
Abstract
Pluripotent stem cells derived from both embryonic and reprogrammed somatic cells have significant potential for human regenerative medicine. Despite similarities in developmental potential, however, several groups have found fundamental differences between embryonic stem cell (ESC) and induced-pluripotent stem cell (iPSC) lines that may have important implications for iPSC-based medical therapies. Using an unsupervised clustering algorithm, we further studied the genetic homogeneity of iPSC and ESC lines by reanalyzing microarray gene expression data from seven different laboratories. Unexpectedly, this analysis revealed a strong correlation between gene expression signatures and specific laboratories in both ESC and iPSC lines. Nearly one-third of the genes with lab-specific expression signatures are also differentially expressed between ESCs and iPSCs. These data are consistent with the hypothesis that in vitro microenvironmental context differentially impacts the gene expression signatures of both iPSCs and ESCs.
View details for DOI 10.1016/j.stem.2010.06.016
View details for Web of Science ID 000281107400017
View details for PubMedID 20682451
-
Community assessment of methods to deconvolve cellular composition from bulk gene expression.
Nature communications
2024; 15 (1): 7362
Abstract
We evaluate deconvolution methods, which infer levels of immune infiltration from bulk expression of tumor samples, through a community-wide DREAM Challenge. We assess six published and 22 community-contributed methods using in vitro and in silico transcriptional profiles of admixed cancer and healthy immune cells. Several published methods predict most cell types well, though they either were not trained to evaluate all functional CD8+ T cell states or do so with low accuracy. Several community-contributed methods address this gap, including a deep learning-based approach, whose strong performance establishes the applicability of this paradigm to deconvolution. Despite being developed largely using immune cells from healthy tissues, deconvolution methods predict levels of tumor-derived immune cells well. Our admixed and purified transcriptional profiles will be a valuable resource for developing deconvolution methods, including in response to common challenges we observe across methods, such as sensitive identification of functional CD4+ T cell states.
View details for DOI 10.1038/s41467-024-50618-0
View details for PubMedID 39191725
View details for PubMedCentralID 6167550
-
Quantification of cerebrospinal fluid tumor DNA in lung cancer patients with suspected leptomeningeal carcinomatosis.
NPJ precision oncology
2024; 8 (1): 121
Abstract
Cerebrospinal fluid tumor-derived DNA (CSF-tDNA) analysis is a promising approach for monitoring the neoplastic processes of the central nervous system. We applied a lung cancer-specific sequencing panel (CAPP-Seq) to 81 CSF, blood, and tissue samples from 24 lung cancer patients who underwent lumbar puncture (LP) for suspected leptomeningeal disease (LMD). A subset of the cohort (N = 12) participated in a prospective trial of osimertinib for refractory LMD in which serial LPs were performed before and during treatment. CSF-tDNA variant allele fractions (VAFs) were significantly higher than plasma circulating tumor DNA (ctDNA) VAFs (median CSF-tDNA, 32.7%; median plasma ctDNA, 1.8%; P < 0.0001). Concentrations of tumor DNA in CSF and plasma were positively correlated (Spearman's ρ, 0.45; P = 0.03). For LMD diagnosis, cytology was 81.8% sensitive and CSF-tDNA was 91.7% sensitive. CSF-tDNA was also strongly prognostic for overall survival (HR = 7.1; P = 0.02). Among patients with progression on targeted therapy, resistance mutations, such as EGFR T790M and MET amplification, were common in peripheral blood but were rare in time-matched CSF, indicating differences in resistance mechanisms based on the anatomic compartment. In the osimertinib cohort, patients with CNS progression had increased CSF-tDNA VAFs at follow-up LP. Post-osimertinib CSF-tDNA VAF was strongly prognostic for CNS progression (HR = 6.2, P = 0.009). Detection of CSF-tDNA in lung cancer patients with suspected LMD is feasible and may have clinical utility. CSF-tDNA improves the sensitivity of LMD diagnosis, enables improved prognostication, and drives therapeutic strategies that account for spatial heterogeneity in resistance mechanisms.
View details for DOI 10.1038/s41698-024-00582-1
View details for PubMedID 38806586
View details for PubMedCentralID 5641214
-
Spatially Segregated Macrophage Populations Predict Distinct Outcomes In Colon Cancer.
Cancer discovery
2024
Abstract
Tumor-associated macrophages are transcriptionally heterogeneous, but the spatial distribution and cell interactions that shape macrophage tissue roles remain poorly characterized. Here, we spatially resolve five distinct human macrophage populations in normal and malignant human breast and colon tissue and reveal their cellular associations. This spatial map reveals that distinct macrophage populations reside in spatially segregated micro-environmental niches with conserved cellular compositions that are repeated across healthy and diseased tissue. We show that IL4I1+ macrophages phagocytose dying cells in areas with high cell turnover and predict good outcome in colon cancer. In contrast, SPP1+ macrophages are enriched in hypoxic and necrotic tumor regions and portend worse outcome in colon cancer. A subset of FOLR2+ macrophages is embedded in plasma cell niches. NLRP3+ macrophages co-localize with neutrophils and activate an inflammasome in tumors. Our findings indicate that a limited number of unique human macrophage niches function as fundamental building blocks in tissue.
View details for DOI 10.1158/2159-8290.CD-23-1300
View details for PubMedID 38552005
-
Sarcoma microenvironment cell states and ecosystems are associated with prognosis and predict response to immunotherapy.
Nature cancer
2024
Abstract
Characterization of the diverse malignant and stromal cell states that make up soft tissue sarcomas and their correlation with patient outcomes has proven difficult using fixed clinical specimens. Here, we employed EcoTyper, a machine-learning framework, to identify the fundamental cell states and cellular ecosystems that make up sarcomas on a large scale using bulk transcriptomes with clinical annotations. We identified and validated 23 sarcoma-specific, transcriptionally defined cell states, many of which were highly prognostic of patient outcomes across independent datasets. We discovered three conserved cellular communities or ecotypes associated with underlying genomic alterations and distinct clinical outcomes. We show that one ecotype defined by tumor-associated macrophages and epithelial-like malignant cells predicts response to immune-checkpoint inhibition but not chemotherapy and validate our findings in an independent cohort. Our results may enable identification of patients with soft tissue sarcomas who could benefit from immunotherapy and help develop new therapeutic strategies.
View details for DOI 10.1038/s43018-024-00743-y
View details for PubMedID 38429415
View details for PubMedCentralID 4486342
-
Single Cell Spatial Biology for Precision Cancer Medicine.
Pacific Symposium on Biocomputing. Pacific Symposium on Biocomputing
2023; 28: 549-553
Abstract
In cancer, complex ecosystems of interacting cell types play fundamental roles in tumor development, progression, and response to therapy. However, the cellular organization, community structure, and spatially defined microenvironments of human tumors remain poorly understood. With the emergence of new technologies for high-throughput spatial profiling of complex tissue specimens, it is now possible to identify clinically significant spatial features with high granularity. In this PSB workshop, we will highlight recent advances in this area and explore how single cell spatial profiling can advance precision cancer medicine.
View details for PubMedID 36541010
-
Single-cell mapping identifies MSI+ cells as a common origin for diverse subtypes of pancreatic cancer.
Cancer cell
2023
Abstract
Identifying the cells from which cancers arise is critical for understanding the molecular underpinnings of tumor evolution. To determine whether stem/progenitor cells can serve as cells of origin, we created a Msi2-CreERT2 knock-in mouse. When crossed to CAG-LSL-MycT58A mice, Msi2-CreERT2 mice developed multiple pancreatic cancer subtypes: ductal, acinar, adenosquamous, and rare anaplastic tumors. Combining single-cell genomics with computational analysis of developmental states and lineage trajectories, we demonstrate that MYC preferentially triggers transformation of the most immature MSI2+ pancreas cells into multi-lineage pre-cancer cells. These pre-cancer cells subsequently diverge to establish pancreatic cancer subtypes by activating distinct transcriptional programs and large-scale genomic changes, and enforced expression of specific signals like Ras can redirect subtype specification. This study shows that multiple pancreatic cancer subtypes can arise from a common pool of MSI2+ cells and provides a powerful model to understand and control the programs that shape divergent fates in pancreatic cancer.
View details for DOI 10.1016/j.ccell.2023.09.008
View details for PubMedID 37802055
-
Diversity of intratumoral regulatory T cells in B-cell non-Hodgkin lymphoma.
Blood advances
2023
Abstract
Tumor-infiltrating regulatory T cells (Tregs) contribute to an immunosuppressive tumor microenvironment. Despite extensive studies, the prognostic impact of tumor-infiltrating Tregs in B-cell non-Hodgkin lymphomas (B-NHL) remains unclear. Emerging studies suggest substantial heterogeneity in the phenotypes and suppressive capacities of Tregs, emphasizing the importance of understanding Treg diversity and the need for additional markers to identify highly suppressive Tregs. Here, we applied single-cell RNA sequencing and TCR sequencing combined with high-dimensional cytometry to decipher the heterogeneity of intratumoral Tregs in diffuse large B-cell lymphoma (DLBCL) and follicular lymphoma (FL), compared with non-malignant tonsillar tissue. We identified three distinct transcriptional states of Tregs; resting, activated and unconventional LAG3+FOXP3- Tregs. Activated Tregs were enriched in B-NHL tumors, co-expressed several checkpoint receptors and had stronger immunosuppressive activity compared with resting Tregs. In FL, activated Tregs were found in closer proximity to CD4+ and CD8+ T cells than other cell types. Furthermore, we used a computational approach to develop unique gene signature matrices which were used to enumerate each Treg subset in cohorts with bulk gene expression data. In two independent FL cohorts, activated Tregs was the major subset, and high abundance was associated with adverse outcome. This study demonstrates that Tregs infiltrating NHL tumors are transcriptionally and functionally diverse. Highly immunosuppressive activated Tregs were enriched in tumor tissue but absent in peripheral blood. Our data suggest that a deeper understanding of Treg heterogeneity in B-NHL could open new paths for rational drug design, facilitating selective targeting to improve anti-tumor immunity.
View details for DOI 10.1182/bloodadvances.2023010158
View details for PubMedID 37695745
-
CD4 T cells and toxicity from immune checkpoint blockade.
Immunological reviews
2023
Abstract
Immune-related toxicities, otherwise known as immune-related adverse events (irAEs), occur in a substantial fraction of cancer patients treated with immune checkpoint inhibitors (ICIs). Ranging from asymptomatic to life-threatening, ICI-induced irAEs can result in hospital admission, high-dose corticosteroid treatment, ICI discontinuation, and in some cases, death. A deeper understanding of the factors underpinning severe irAE development will be essential for improved irAE prediction and prevention, toward maximizing the benefits and safety profiles of ICIs. In recent work, we applied mass cytometry, single-cell RNA sequencing, single-cell V(D)J sequencing, bulk RNA sequencing, and bulk T-cell receptor (TCR) sequencing to identify pretreatment determinants of severe irAE development in patients with advanced melanoma. Across 71 patients separated into three cohorts, we found that two baseline features in circulation-elevated activated CD4 effector memory T-cell abundance and TCR diversity-are associated with severe irAE development, independent of the affected organ system within 3 months of ICI treatment initiation. Here, we provide an extended perspective on this work, synthesize and discuss related literature, and summarize practical considerations for clinical translation. Collectively, these findings lay a foundation for data-driven and mechanistic insights into irAE development, with the potential to reduce ICI morbidity and mortality in the future.
View details for DOI 10.1111/imr.13248
View details for PubMedID 37491734
-
p53 governs an AT1 differentiation programme in lung cancer suppression.
Nature
2023
Abstract
Lung cancer is the leading cause of cancer deaths worldwide1. Mutations in the tumour suppressor gene TP53 occur in 50% of lung adenocarcinomas (LUADs) and are linked to poor prognosis1-4, but how p53 suppresses LUAD development remains enigmatic. We show here that p53 suppresses LUAD by governing cell state, specifically by promoting alveolar type 1 (AT1) differentiation. Using mice that express oncogenic Kras and null, wild-type or hypermorphic Trp53 alleles in alveolar type 2 (AT2) cells, we observed graded effects of p53 on LUAD initiation and progression. RNA sequencing and ATAC sequencing of LUAD cells uncovered a p53-induced AT1 differentiation programme during tumour suppression in vivo through direct DNA binding, chromatin remodelling and induction of genes characteristic of AT1 cells. Single-cell transcriptomics analyses revealed that during LUAD evolution, p53 promotes AT1 differentiation through action in a transitional cell state analogous to a transient intermediary seen during AT2-to-AT1 cell differentiation in alveolar injury repair. Notably, p53 inactivation results in the inappropriate persistence of these transitional cancer cells accompanied by upregulated growth signalling and divergence from lung lineage identity, characteristics associated with LUAD progression. Analysis of Trp53 wild-type and Trp53-null mice showed that p53 also directs alveolar regeneration after injury by regulating AT2 cell self-renewal and promoting transitional cell differentiation into AT1 cells. Collectively, these findings illuminate mechanisms of p53-mediated LUAD suppression, in which p53 governs alveolar differentiation, and suggest that tumour suppression reflects a fundamental role of p53 in orchestrating tissue repair after injury.
View details for DOI 10.1038/s41586-023-06253-8
View details for PubMedID 37468633
View details for PubMedCentralID 4231481
-
Profiling Cellular Ecosystems at Single-Cell Resolution and at Scale with EcoTyper.
Methods in molecular biology (Clifton, N.J.)
2023; 2629: 43-71
Abstract
Tissues are composed of diverse cell types and cellular states that organize into distinct ecosystems with specialized functions. EcoTyper is a collection of machine learning tools for the large-scale delineation of cellular ecosystems and their constituent cell states from bulk, single-cell, and spatially resolved gene expression data. In this chapter, we provide a primer on EcoTyper and demonstrate its use for the discovery and recovery of cell states and ecosystems from healthy and diseased tissue specimens.
View details for DOI 10.1007/978-1-0716-2986-4_4
View details for PubMedID 36929073
View details for PubMedCentralID 9067608
-
Multiomic analysis for optimization of combined focal and immunotherapy protocols in murine pancreatic cancer.
Theranostics
2022; 12 (18): 7884-7902
Abstract
Background: Although combination immunotherapies incorporating local and systemic components have shown promising results in treating solid tumors, varied tumor microenvironments (TMEs) can impact immunotherapeutic efficacy. Method: We designed and evaluated treatment strategies for breast and pancreatic cancer combining magnetic resonance-guided focused ultrasound (MRgFUS) ablation and antibody therapies. With a combination of single-cell sequencing, spectral flow cytometry, and histological analyses, we profiled an immune-suppressed KPC (Kras+/LSL-G12D; Trp53+/LSL-R172H; Pdx1-Cre) pancreatic adenocarcinoma (MT4) model and a dense epithelial neu deletion (NDL) HER2+ mammary adenocarcinoma model with a greater fraction of lymphocytes, natural killer cells and activated dendritic cells. We then performed gene ontology analysis, spectral and digital cytometry to assess the immune response to combination immunotherapies and correlation with survival studies. Result: Based on gene ontology analysis, adding ablation to immunotherapy enriched immune cell migration pathways in the pancreatic cancer model and extensively enriched wound healing pathways in the breast cancer model. With CIBERSORTx digital cytometry, aCD40 + aPD-1 immunotherapy combinations enhanced dendritic cell activation in both models. In the MT4 TME, adding the combination of aCD40 antibody and checkpoint inhibitors (aPD-1 and aCTLA-4) with ablation was synergistic, increasing activated natural killer cells and T cells in distant tumors. Furthermore, ablation with immunotherapy upregulated critical Ly6c myeloid remodeling phenotypes that enhance T-cell effector function and increased granzyme and protease encoding genes by as much as 100-fold. Ablation combined with immunotherapy then extended survival in the MT4 model to a greater extent than immunotherapy alone. Conclusion: In summary, TME profiling informed a successful multicomponent treatment protocol incorporating ablation and facilitated differentiation of TMEs in which ablation is most effective.
View details for DOI 10.7150/thno.73218
View details for PubMedID 36451859
View details for PubMedCentralID PMC9706583
-
Identification of a minority population of LMO2+ breast cancer cells that integrate into the vasculature and initiate metastasis.
Science advances
2022; 8 (45): eabm3548
Abstract
Metastasis is responsible for most breast cancer-related deaths; however, identifying the cellular determinants of metastasis has remained challenging. Here, we identified a minority population of immature THY1+/VEGFA+ tumor epithelial cells in human breast tumor biopsies that display angiogenic features and are marked by the expression of the oncogene, LMO2. Higher abundance of LMO2+ basal cells correlated with tumor endothelial content and predicted poor distant recurrence-free survival in patients. Using MMTV-PyMT/Lmo2CreERT2 mice, we demonstrated that Lmo2 lineage-traced cells integrate into the vasculature and have a higher propensity to metastasize. LMO2 knockdown in human breast tumors reduced lung metastasis by impairing intravasation, leading to a reduced frequency of circulating tumor cells. Mechanistically, we find that LMO2 binds to STAT3 and is required for STAT3 activation by tumor necrosis factor-α and interleukin-6. Collectively, our study identifies a population of metastasis-initiating cells with angiogenic features and establishes the LMO2-STAT3 signaling axis as a therapeutic target in breast cancer metastasis.
View details for DOI 10.1126/sciadv.abm3548
View details for PubMedID 36351009
-
Inferring gene expression from cell-free DNA fragmentation profiles.
Nature biotechnology
2022
Abstract
Profiling of circulating tumor DNA (ctDNA) in the bloodstream shows promise for noninvasive cancer detection. Chromatin fragmentation features have previously been explored to infer gene expression profiles from cell-free DNA (cfDNA), but current fragmentomic methods require high concentrations of tumor-derived DNA and provide limited resolution. Here we describe promoter fragmentation entropy as an epigenomic cfDNA feature that predicts RNA expression levels at individual genes. We developed 'epigenetic expression inference from cell-free DNA-sequencing' (EPIC-seq), a method that uses targeted sequencing of promoters of genes of interest. Profiling 329 blood samples from 201 patients with cancer and 87 healthy adults, we demonstrate classification of subtypes of lung carcinoma and diffuse large B cell lymphoma. Applying EPIC-seq to serial blood samples from patients treated with PD-(L)1 immune-checkpoint inhibitors, we show that gene expression profiles inferred by EPIC-seq are correlated with clinical response. Our results indicate that EPIC-seq could enable noninvasive, high-throughput tissue-of-origin characterization with diagnostic, prognostic and therapeutic potential.
View details for DOI 10.1038/s41587-022-01222-4
View details for PubMedID 35361996
-
Integrative molecular and clinical profiling of acral melanoma links focal amplification of 22q11.21 to metastasis.
Nature communications
2022; 13 (1): 898
Abstract
Acral melanoma, the most common melanoma subtype among non-White individuals, is associated with poor prognosis. However, its key molecular drivers remain obscure. Here, we perform integrative genomic and clinical profiling of acral melanomas from 104 patients treated in North America (n = 37) or China (n = 67). We find that recurrent, late-arising focal amplifications of cytoband 22q11.21 are a leading determinant of inferior survival, strongly associated with metastasis, and linked to downregulation of immunomodulatory genes associated with response to immune checkpoint blockade. Unexpectedly, LZTR1 - a known tumor suppressor in other cancers - is a key candidate oncogene in this cytoband. Silencing of LZTR1 in melanoma cell lines causes apoptotic cell death independent of major hotspot mutations or melanoma subtypes. Conversely, overexpression of LZTR1 in normal human melanocytes initiates processes associated with metastasis, including anchorage-independent growth, formation of spheroids, and an increase in MAPK and SRC activities. Our results provide insights into the etiology of acral melanoma and implicate LZTR1 as a key tumor promoter and therapeutic target.
View details for DOI 10.1038/s41467-022-28566-4
View details for PubMedID 35197475
-
Integrated spatial multiomics reveals fibroblast fate during tissue repair.
Proceedings of the National Academy of Sciences of the United States of America
2021; 118 (41)
Abstract
In the skin, tissue injury results in fibrosis in the form of scars composed of dense extracellular matrix deposited by fibroblasts. The therapeutic goal of regenerative wound healing has remained elusive, in part because principles of fibroblast programming and adaptive response to injury remain incompletely understood. Here, we present a multimodal -omics platform for the comprehensive study of cell populations in complex tissue, which has allowed us to characterize the cells involved in wound healing across both time and space. We employ a stented wound model that recapitulates human tissue repair kinetics and multiple Rainbow transgenic lines to precisely track fibroblast fate during the physiologic response to skin injury. Through integrated analysis of single cell chromatin landscapes and gene expression states, coupled with spatial transcriptomic profiling, we are able to impute fibroblast epigenomes with temporospatial resolution. This has allowed us to reveal potential mechanisms controlling fibroblast fate during migration, proliferation, and differentiation following skin injury, and thereby reexamine the canonical phases of wound healing. These findings have broad implications for the study of tissue repair in complex organ systems.
View details for DOI 10.1073/pnas.2110025118
View details for PubMedID 34620713
-
Landscape of innate lymphoid cells in human head and neck cancer reveals divergent NK cell states in the tumor microenvironment.
Proceedings of the National Academy of Sciences of the United States of America
2021; 118 (28)
Abstract
Natural killer (NK) cells comprise one subset of the innate lymphoid cell (ILC) family. Despite reported antitumor functions of NK cells, their tangible contribution to tumor control in humans remains controversial. This is due to incomplete understanding of the NK cell states within the tumor microenvironment (TME). Here, we demonstrate that peripheral circulating NK cells differentiate down two divergent pathways within the TME, resulting in different end states. One resembles intraepithelial ILC1s (ieILC1) and possesses potent in vivo antitumor activity. The other expresses genes associated with immune hyporesponsiveness and has poor antitumor functional capacity. Interleukin-15 (IL-15) and direct contact between the tumor cells and NK cells are required for the differentiation into CD49a+CD103+ cells, resembling ieILC1s. These data explain the similarity between ieILC1s and tissue-resident NK cells, provide insight into the origin of ieILC1s, and identify the ieILC1-like cell state within the TME to be the NK cell phenotype with the greatest antitumor activity. Because the proportions of the different ILC states vary between tumors, these findings provide a resource for the clinical study of innate immune responses against tumors and the design of novel therapy.
View details for DOI 10.1073/pnas.2101169118
View details for PubMedID 34244432
-
The landscape of tumor cell states and ecosystems in diffuse large B cell lymphoma.
Cancer cell
2021
Abstract
Biological heterogeneity in diffuse large B cell lymphoma (DLBCL) is partly driven by cell-of-origin subtypes and associated genomic lesions, but also by diverse cell types and cell states in the tumor microenvironment (TME). However, dissecting these cell states and their clinical relevance at scale remains challenging. Here, we implemented EcoTyper, a machine-learning framework integrating transcriptome deconvolution and single-cell RNA sequencing, to characterize clinically relevant DLBCL cell states and ecosystems. Using this approach, we identified five cell states of malignant B cells that vary in prognostic associations and differentiation status. We also identified striking variation in cell states for 12 other lineages comprising the TME and forming cell state interactions in stereotyped ecosystems. While cell-of-origin subtypes have distinct TME composition, DLBCL ecosystems capture clinical heterogeneity within existing subtypes and extend beyond cell-of-origin and genotypic classes. These results resolve the DLBCL microenvironment at systems-level resolution and identify opportunities for therapeutic targeting (https://ecotyper.stanford.edu/lymphoma).
View details for DOI 10.1016/j.ccell.2021.08.011
View details for PubMedID 34597589
-
LEFTY1 Is a Dual-SMAD Inhibitor that Promotes Mammary Progenitor Growth and Tumorigenesis.
Cell stem cell
2020
Abstract
SMAD pathways govern epithelial proliferation, and transforming growth factor beta (TGF-beta and BMP signaling through SMAD members has distinct effects on mammary development and homeostasis. Here, we show that LEFTY1, a secreted inhibitor of NODAL/SMAD2 signaling, is produced by mammary progenitor cells and, concomitantly, suppresses SMAD2 and SMAD5 signaling to promote long-term proliferation of normal and malignant mammary epithelial cells. In contrast, BMP7, a NODAL antagonist with context-dependent functions, is produced by basal cells and restrains progenitor cell proliferation. In normal mouse epithelium, LEFTY1 expression in a subset of luminal cells and rare basal cells opposes BMP7 to promote ductal branching. LEFTY1 binds BMPR2 to suppress BMP7-induced activation of SMAD5, and this LEFTY1-BMPR2 interaction is specific to tumor-initiating cells in triple-negative breast cancer xenografts that rely on LEFTY1 for growth. These results suggest that LEFTY1 is an endogenous dual-SMAD inhibitor and that suppressing its function may represent a therapeutic vulnerability in breast cancer.
View details for DOI 10.1016/j.stem.2020.06.017
View details for PubMedID 32693087
-
Analytical validation of digital cytometry (iSort) for leukocyte enumeration using stored blood.
AMER SOC CLINICAL ONCOLOGY. 2020
View details for Web of Science ID 000560368301434
-
Analytical validation of iSort digital cytometry for leukocyte enumeration in clinical tumor specimens.
LIPPINCOTT WILLIAMS & WILKINS. 2020
View details for Web of Science ID 000560368305213
-
Profiling Cell Type Abundance and Expression in Bulk Tissues with CIBERSORTx.
Methods in molecular biology (Clifton, N.J.)
2020; 2117: 135–57
Abstract
CIBERSORTx is a suite of machine learning tools for the assessment of cellular abundance and cell type-specific gene expression patterns from bulk tissue transcriptome profiles. With this framework, single-cell or bulk-sorted RNA sequencing data can be used to learn molecular signatures of distinct cell types from a small collection of biospecimens. These signatures can then be repeatedly applied to characterize cellular heterogeneity from bulk tissue transcriptomes without physical cell isolation. In this chapter, we provide a detailed primer on CIBERSORTx and demonstrate its capabilities for high-throughput profiling of cell types and cellular states in normal and neoplastic tissues.
View details for DOI 10.1007/978-1-0716-0301-7_7
View details for PubMedID 31960376
-
Noninvasive Early Identification of Therapeutic Benefit from Immune Checkpoint Inhibition.
Cell
2020
Abstract
Although treatment of non-small cell lung cancer (NSCLC) with immune checkpoint inhibitors (ICIs) can produce remarkably durable responses, most patients develop early disease progression. Furthermore, initial response assessment by conventional imaging is often unable to identify which patients will achieve durable clinical benefit (DCB). Here, we demonstrate that pre-treatment circulating tumor DNA (ctDNA) and peripheral CD8 T cell levels are independently associated with DCB. We further show that ctDNA dynamics after a single infusion can aid in identification of patients who will achieve DCB. Integrating these determinants, we developed and validated an entirely noninvasive multiparameter assay (DIREct-On, Durable Immunotherapy Response Estimation by immune profiling and ctDNA-On-treatment) that robustly predicts which patients will achieve DCB with higher accuracy than any individual feature. Taken together, these results demonstrate that integrated ctDNA and circulating immune cell profiling can provide accurate, noninvasive, and early forecasting of ultimate outcomes for NSCLC patients receiving ICIs.
View details for DOI 10.1016/j.cell.2020.09.001
View details for PubMedID 33007267
-
Multiomic single cell analysis of normal human bone marrow identifies a unique stem and progenitor population that expands in AML
Proceedings of the Annual Meeting of the American Association for Cancer Research 2020
2020
View details for DOI 10.1158/1538-7445.AM2020-3779
-
Molecular and immunological signatures are related to clinical benefit from treatment with Vocimagene amiretrorepvec (Toca 511) and 5-fluorocytosine (Toca FC) in patients with glioma.
Clinical cancer research : an official journal of the American Association for Cancer Research
2020
Abstract
High grade gliomas are central nervous system tumors with poor prognoses and limited treatment options. Vocimagene amiretrorepvec (Toca 511) is a retroviral replicating vector encoding cytosine deaminase, which converts extended-release 5-fluorocytosine (Toca FC) into the anti-cancer agent 5-fluorouracil. According to preclinical studies, this therapy kills cancer cells and immunosuppressive myeloid cells in the tumor microenvironment, leading to T-cell mediated anti-tumor immune activity. Therefore, we sought to elucidate this immune-related mechanism of action in humans, and to investigate potential molecular and immunological indicators of clinical benefit from therapy.In a Phase I clinical trial (NCT01470794), recurrent high grade glioma patients treated with Toca 511 and Toca FC showed improved survival relative to historical controls, and some had durable complete responses to therapy. As part of this trial, we performed whole exome DNA sequencing, RNA sequencing and multiplex digital ELISA measurements on tumor and blood samples.Genetic analyses suggest mutations, copy number variations and neoantigens are linked to survival. Quantities of tumor immune infiltrates estimated by transcript abundance may potentially predict clinical outcomes. Peak values of cytokines in peripheral blood samples collected during and after therapy could indicate response.These results support an immune-related mechanism of action for Toca 511 and Toca FC, and suggest that molecular and immunological signatures are related to clinical benefit from treatment.
View details for DOI 10.1158/1078-0432.CCR-20-0536
View details for PubMedID 32816892
-
Atlas of clinically-distinct cell states and cellular ecosystems across human solid tumors
BMC. 2019
View details for Web of Science ID 000496473200425
-
The Immune Landscape of Cancer.
Immunity
2019; 51 (2): 411–12
View details for DOI 10.1016/j.immuni.2019.08.004
View details for PubMedID 31433971
-
Dynamic Risk Profiling Using Serial Tumor Biomarkers for Personalized Outcome Prediction.
Cell
2019
Abstract
Accurate prediction of long-term outcomes remains a challenge in the care of cancer patients. Due to the difficulty of serial tumor sampling, previous prediction tools have focused on pretreatment factors. However, emerging non-invasive diagnostics have increased opportunities for serial tumor assessments. We describe the Continuous Individualized Risk Index (CIRI), a method to dynamically determine outcome probabilities for individual patients utilizing risk predictors acquired over time. Similar to "win probability" models in other fields, CIRI provides a real-time probability by integrating risk assessments throughout a patient's course. Applying CIRI to patients with diffuse large B cell lymphoma, we demonstrate improved outcome prediction compared to conventional risk models. We demonstrate CIRI's broader utility in analogous models of chronic lymphocytic leukemia and breast adenocarcinoma and perform a proof-of-concept analysis demonstrating how CIRI could be used to develop predictive biomarkers for therapy selection. We envision thatdynamic risk assessment will facilitate personalized medicine and enable innovative therapeutic paradigms.
View details for DOI 10.1016/j.cell.2019.06.011
View details for PubMedID 31280963
-
A tumor deconvolution DREAM Challenge: Inferring immune infiltration from bulk gene expression data
AMER ASSOC CANCER RESEARCH. 2019
View details for DOI 10.1158/1538-7445.AM2019-1690
View details for Web of Science ID 000488129904217
-
A functional subset of CD8+ T cells during chronic exhaustion is defined by SIRPalpha expression.
Nature communications
2019; 10 (1): 794
Abstract
Prolonged exposure of CD8+ T cells to antigenic stimulation, as in chronic viral infections, leads to a state of diminished function termed exhaustion. We now demonstrate that even during exhaustion there is a subset of functional CD8+ T cells defined by surface expression of SIRPalpha, a protein not previously reported on lymphocytes. On SIRPalpha+ CD8+ T cells, expression of co-inhibitory receptors is counterbalanced by expression of co-stimulatory receptors and it is only SIRPalpha+ cells that actively proliferate, transcribe IFNgamma and show cytolytic activity. Furthermore, target cells that express the ligand for SIRPalpha, CD47, are more susceptible to CD8+ T cell-killing in vivo. SIRPalpha+ CD8+ T cells are evident in mice infected with Friend retrovirus, LCMV Clone 13, and in patients with chronic HCV infections. Furthermore, therapeutic blockade of PD-L1 to reinvigorate CD8+ T cells during chronic infection expands the cytotoxic subset of SIRPalpha+ CD8+ T cells.
View details for PubMedID 30770827
-
Reply to J. Wang et al.
Journal of clinical oncology : official journal of the American Society of Clinical Oncology
2019: JCO1801907
View details for PubMedID 30753108
-
Spatial mapping of the immune microenvironment in primary triple-negative breast cancer (TNBC) and association with neoadjuvant therapy response
AMER ASSOC CANCER RESEARCH. 2019
View details for DOI 10.1158/1538-7445.SABCS18-PD5-12
View details for Web of Science ID 000478677003010
-
Computational approaches for characterizing the tumor immune microenvironment.
Immunology
2019
Abstract
Recent advances in high-throughput molecular profiling technologies and multiplexed imaging platforms have revolutionized our ability to characterize the tumor immune microenvironment. As a result, studies of tumor-associated immune cells increasingly involve complex datasets that require sophisticated methods of computational analysis. In this review, we present an overview of key assays and related bioinformatics tools for analyzing the tumor-associated immune system in bulk tissues and at the single-cell level. In parallel, we describe how data science strategies and novel technologies have advanced tumor immunology and opened the door for new opportunities to exploit host immunity to improve cancer clinical outcomes. This article is protected by copyright. All rights reserved.
View details for DOI 10.1111/imm.13101
View details for PubMedID 31347163
-
Functional significance of U2AF1 S34F mutations in lung adenocarcinomas.
Nature communications
2019; 10 (1): 5712
Abstract
The functional role of U2AF1 mutations in lung adenocarcinomas (LUADs) remains incompletely understood. Here, we report a significant co-occurrence of U2AF1 S34F mutations with ROS1 translocations in LUADs. To characterize this interaction, we profiled effects of S34F on the transcriptome-wide distribution of RNA binding and alternative splicing in cells harboring the ROS1 translocation. Compared to its wild-type counterpart, U2AF1 S34F preferentially binds and modulates splicing of introns containing CAG trinucleotides at their 3' splice junctions. The presence of S34F caused a shift in cross-linking at 3' splice sites, which was significantly associated with alternative splicing of skipped exons. U2AF1 S34F induced expression of genes involved in the epithelial-mesenchymal transition (EMT) and increased tumor cell invasion. Finally, S34F increased splicing of the long over the short SLC34A2-ROS1 isoform, which was also associated with enhanced invasiveness. Taken together, our results suggest a mechanistic interaction between mutant U2AF1 and ROS1 in LUAD.
View details for DOI 10.1038/s41467-019-13392-y
View details for PubMedID 31836708
-
Circulating tumor DNA analysis for detection of minimal residual disease after chemoradiotherapy for localized esophageal cancer.
Gastroenterology
2019
Abstract
Biomarkers are needed to identify patients at risk of tumor progression following chemoradiotherapy for localized esophageal cancer. These could improve identification of patients at risk for cancer progression and selection of therapy.We performed deep sequencing (CAPP-Seq) analyses of plasma cell-free DNA collected from 45 patients before and after chemoradiotherapy for esophageal cancer, as well as DNA from leukocytes, and fixed esophageal tumor biopsies collected during esophagogastroduodenoscopy. Patients were treated from May 2010 through October 2015; 23 patients subsequently underwent esophagectomy and 22 did not undergo surgery. We also sequenced DNA from blood samples from 40 healthy individuals (controls). We analyzed 802 regions of 607 genes for single-nucleotide variants previously associated with esophageal adenocarcinoma or squamous cell carcinoma. Patients underwent imaging analyses 6-8 weeks after chemoradiotherapy and were followed for 5 years. Our primary aim was to determine whether detection of circulating tumor DNA (ctDNA) following chemoradiotherapy is associated with risk of tumor progression (growth of local, regional, or distant tumors, detected by imaging or biopsy).The median proportion of tumor-derived DNA in total cell-free DNA before treatment was 0.07%, indicating that ultrasensitive assays are needed for quantification and analysis of ctDNA from localized esophageal tumors. Detection of ctDNA following chemoradiotherapy was associated with tumor progression (hazard ratio, 18.7; P<.0001), formation of distant metastases (hazard ratio, 32.1; P<.0001), and shorter disease-specific survival times (hazard ratio, 23.1; P<.0001). A higher proportion of patients with tumor progression had new mutations detected in plasma samples collected after chemoradiotherapy than patients without progression (P=.03). Detection of ctDNA after chemoradiotherapy preceded radiographic evidence of tumor progression by an average of 2.8 months. Among patients who received chemoradiotherapy without surgery, combined ctDNA and metabolic imaging analysis predicted progression in 100% of patients with tumor progression, compared with 71% for only ctDNA detection and 57% for only metabolic imaging analysis (P<.001 for comparison of either technique to combined analysis).In an analysis of cell-free DNA in blood samples from patients who underwent chemoradiotherapy for esophageal cancer, detection of ctDNA was associated with tumor progression, metastasis, and disease-specific survival. Analysis of ctDNA might be used to identify patients at highest risk for tumor progression.
View details for DOI 10.1053/j.gastro.2019.10.039
View details for PubMedID 31711920
-
Method of Isolating and Transplanting the Hematopoietic Stem Cell with Its Microenvironment Which Improves Functional Hematopoietic Engraftment
ELSEVIER SCIENCE INC. 2018: E224
View details for DOI 10.1016/j.jamcollsurg.2018.08.605
View details for Web of Science ID 000447772500541
-
Circulating tumor DNA (ctDNA) in B-cell lymphoma
WILEY. 2018: 16–17
View details for Web of Science ID 000444944200019
-
Circulating Tumor DNA Measurements As Early Outcome Predictors in Diffuse Large B-Cell Lymphoma.
Journal of clinical oncology : official journal of the American Society of Clinical Oncology
2018: JCO2018785246
Abstract
Purpose Outcomes for patients with diffuse large B-cell lymphoma remain heterogeneous, with existing methods failing to consistently predict treatment failure. We examined the additional prognostic value of circulating tumor DNA (ctDNA) before and during therapy for predicting patient outcomes. Patients and Methods We studied the dynamics of ctDNA from 217 patients treated at six centers, using a training and validation framework. We densely characterized early ctDNA dynamics during therapy using cancer personalized profiling by deep sequencing to define response-associated thresholds within a discovery set. These thresholds were assessed in two independent validation sets. Finally, we assessed the prognostic value of ctDNA in the context of established risk factors, including the International Prognostic Index and interim positron emission tomography/computed tomography scans. Results Before therapy, ctDNA was detectable in 98% of patients; pretreatment levels were prognostic in both front-line and salvage settings. In the discovery set, ctDNA levels changed rapidly, with a 2-log decrease after one cycle (early molecular response [EMR]) and a 2.5-log decrease after two cycles (major molecular response [MMR]) stratifying outcomes. In the first validation set, patients receiving front-line therapy achieving EMR or MMR had superior outcomes at 24 months (EMR: EFS, 83% v 50%; P = .0015; MMR: EFS, 82% v 46%; P < .001). EMR also predicted superior 24-month outcomes in patients receiving salvage therapy in the first validation set (EFS, 100% v 13%; P = .011). The prognostic value of EMR and MMR was further confirmed in the second validation set. In multivariable analyses including International Prognostic Index and interim positron emission tomography/computed tomography scans across both cohorts, molecular response was independently prognostic of outcomes, including event-free and overall survival. Conclusion Pretreatment ctDNA levels and molecular responses are independently prognostic of outcomes in aggressive lymphomas. These risk factors could potentially guide future personalized risk-directed approaches.
View details for PubMedID 30125215
-
Circulating Tumor DNA Quantitation for Early Response Assessment of Immune Checkpoint Inhibitors for Metastatic Non-Small Cell Lung Cancer
ELSEVIER SCIENCE INC. 2018: E1–E2
View details for Web of Science ID 000432447200003
-
Genomic Feature Selection by Coverage Design Optimization.
Journal of applied statistics
2018; 45 (14): 2658-2676
Abstract
We introduce a novel data reduction technique whereby we select a subset of tiles to "cover" maximally events of interest in large-scale biological datasets (e.g., genetic mutations), while minimizing the number of tiles. A tile is a genomic unit capturing one or more biological events, such as a sequence of base pairs that can be sequenced and observed simultaneously. The goal is to reduce significantly the number of tiles considered to those with areas of dense events in a cohort, thus saving on cost and enhancing interpretability. However, the reduction should not come at the cost of too much information, allowing for sensible statistical analysis after its application. We envisage application of our methods to a variety of high throughput data types, particularly those produced by next generation sequencing (NGS) experiments. The procedure is cast as a convex optimization problem, which is presented, along with methods of its solution. The method is demonstrated on a large dataset of somatic mutations spanning 5000+ patients, each having one of 29 cancer types. Applied to these data, our method dramatically reduces the number of gene locations required for broad coverage of patients and their mutations, giving subject specialists a more easily interpretable snapshot of recurrent mutational profiles in these cancers. The locations identified coincide with previously identified cancer genes. Finally, despite considerable data reduction, we show that our covering designs preserve the cancer discrimination ability of multinomial logistic regression models trained on all of the locations (> 1M).
View details for DOI 10.1080/02664763.2018.1432577
View details for PubMedID 30294060
View details for PubMedCentralID PMC6173524
-
Early B cell changes predict autoimmunity following combination immune checkpoint blockade.
The Journal of clinical investigation
2018
Abstract
Combination checkpoint blockade (CCB) targeting inhibitory CTLA4 and PD1 receptors holds promise for cancer therapy. Immune-related adverse events (IRAEs) remain a major obstacle for the optimal application of CCB in cancer. Here, we analyzed B cell changes in patients with melanoma following treatment with either anti-CTLA4 or anti-PD1, or in combination. CCB therapy led to changes in circulating B cells that were detectable after the first cycle of therapy and characterized by a decline in circulating B cells and an increase in CD21lo B cells and plasmablasts. PD1 expression was higher in the CD21lo B cells, and B cell receptor sequencing of these cells demonstrated greater clonality and a higher frequency of clones compared with CD21hi cells. CCB induced proliferation in the CD21lo compartment, and single-cell RNA sequencing identified B cell activation in cells with genomic profiles of CD21lo B cells in vivo. Increased clonality of circulating B cells following CCB occurred in some patients. Treatment-induced changes in B cells preceded and correlated with both the frequency and timing of IRAEs. Patients with early B cell changes experienced higher rates of grade 3 or higher IRAEs 6 months after CCB. Thus, early changes in B cells following CCB may identify patients who are at increased risk of IRAEs, and preemptive strategies targeting B cells may reduce toxicities in these patients.
View details for PubMedID 29309048
-
Complex mammalian-like haematopoietic system found in a colonial chordate.
Nature
2018
Abstract
Haematopoiesis is an essential process that evolved in multicellular animals. At the heart of this process are haematopoietic stem cells (HSCs), which are multipotent and self-renewing, and generate the entire repertoire of blood and immune cells throughout an animal's life1. Although there have been comprehensive studies on self-renewal, differentiation, physiological regulation and niche occupation in vertebrate HSCs, relatively little is known about the evolutionary origin and niches of these cells. Here we describe the haematopoietic system of Botryllus schlosseri, a colonial tunicate that has a vasculature and circulating blood cells, and interesting stem-cell biology and immunity characteristics2-8. Self-recognition between genetically compatible B. schlosseri colonies leads to the formation of natural parabionts with shared circulation, whereas incompatible colonies reject each other3,4,7. Using flow cytometry, whole-transcriptome sequencing of defined cell populations and diverse functional assays, we identify HSCs, progenitors, immune effector cells and an HSC niche, and demonstrate that self-recognition inhibits allospecific cytotoxic reactions. Our results show that HSC and myeloid lineage immune cells emerged in a common ancestor of tunicates and vertebrates, and also suggest that haematopoietic bone marrow and the B. schlosseri endostyle niche evolved from a common origin.
View details for PubMedID 30518860
-
Circulating tumor DNA levels correlate with response to treatment in LMS patients
AMER ASSOC CANCER RESEARCH. 2018: 38–39
View details for Web of Science ID 000422882000043
-
Combination approach for detecting different types of alterations in circulating tumor DNA in leiomyosarcoma.
Clinical cancer research : an official journal of the American Association for Cancer Research
2018
Abstract
The clinical utility of circulating tumor DNA (ctDNA) monitoring has been shown in tumors that harbor highly recurrent mutations. Leiomyosarcoma (LMS) represents a type of tumor with a wide spectrum of heterogeneous genomic abnormalities; thus, targeting hotspot mutations or a narrow genomic region for ctDNA detection may not be practical. Here we demonstrate a combinatorial approach that integrates different sequencing protocols for the orthogonal detection of single nucleotide variants (SNVs), small indels and copy number alterations (CNAs) in ctDNA.We employed Cancer Personalized Profiling by deep Sequencing (CAPP-Seq) for the analysis of SNVs and indels, together with a genome-wide interrogation of CNAs by Genome Representation Profiling (GRP). We profiled 28 longitudinal plasma samples and 25 tumor specimens from 7 patients with LMS.We detected ctDNA in 6 of 7 of these patients with >98% specificity for mutant allele fractions down to a level of 0.01%. We show that results from CAPP-Seq and GRP are highly concordant, and the combination of these methods allows for more comprehensive monitoring of ctDNA by profiling a wide spectrum of tumor-specific markers. By analyzing multiple tumor specimens in individual patients obtained from different sites and at different times during treatment, we observed clonal evolution of these tumors that was reflected by ctDNA profiles.Our strategy allows for a comprehensive monitoring of a broad spectrum of tumor-specific markers in plasma. Our approach may be clinically useful not only in LMS but also in other tumor types that lack recurrent genomic alterations.
View details for PubMedID 29463554
-
Profiling Tumor Infiltrating Immune Cells with CIBERSORT.
Methods in molecular biology (Clifton, N.J.)
2018; 1711: 243–59
Abstract
Tumor infiltrating leukocytes (TILs) are an integral component of the tumor microenvironment and have been found to correlate with prognosis and response to therapy. Methods to enumerate immune subsets such as immunohistochemistry or flow cytometry suffer from limitations in phenotypic markers and can be challenging to practically implement and standardize. An alternative approach is to acquire aggregative high dimensional data from cellular mixtures and to subsequently infer the cellular components computationally. We recently described CIBERSORT, a versatile computational method for quantifying cell fractions from bulk tissue gene expression profiles (GEPs). Combining support vector regression with prior knowledge of expression profiles from purified leukocyte subsets, CIBERSORT can accurately estimate the immune composition of a tumor biopsy. In this chapter, we provide a primer on the CIBERSORT method and illustrate its use for characterizing TILs in tumor samples profiled by microarray or RNA-Seq.
View details for PubMedID 29344893
-
Macrophage infiltration and genetic landscape of undifferentiated uterine sarcomas.
JCI insight
2017; 2 (11)
Abstract
Endometrial stromal tumors include translocation-associated low- and high-grade endometrial stromal sarcomas (ESS) and highly malignant undifferentiated uterine sarcomas (UUS). UUS is considered a poorly defined group of aggressive tumors and is often seen as a diagnosis of exclusion after ESS and leiomyosarcoma (LMS) have been ruled out. We performed a comprehensive analysis of gene expression, copy number variation, point mutations, and immune cell infiltrates in the largest series to date of all major types of uterine sarcomas to shed light on the biology of UUS and to identify potential novel therapeutic targets. We show that UUS tumors have a distinct molecular profile from LMS and ESS. Gene expression and immunohistochemical analyses revealed the presence of high numbers of tumor-associated macrophages (TAMs) in UUS, which makes UUS patients suitable candidates for therapies targeting TAMs. Our results show a high genomic instability of UUS and downregulation of several TP53-mediated tumor suppressor genes, such as NDN, CDH11, and NDRG4. Moreover, we demonstrate that UUS carry somatic mutations in several oncogenes and tumor suppressor genes implicated in RAS/PI3K/AKT/mTOR, ERBB3, and Hedgehog signaling.
View details for DOI 10.1172/jci.insight.94033
View details for PubMedID 28570276
-
Antigen presentation profiling reveals recognition of lymphoma immunoglobulin neoantigens
NATURE
2017; 543 (7647): 723-?
Abstract
Cancer somatic mutations can generate neoantigens that distinguish malignant from normal cells. However, the personalized identification and validation of neoantigens remains a major challenge. Here we discover neoantigens in human mantle-cell lymphomas by using an integrated genomic and proteomic strategy that interrogates tumour antigen peptides presented by major histocompatibility complex (MHC) class I and class II molecules. We applied this approach to systematically characterize MHC ligands from 17 patients. Remarkably, all discovered neoantigenic peptides were exclusively derived from the lymphoma immunoglobulin heavy- or light-chain variable regions. Although we identified MHC presentation of private polymorphic germline alleles, no mutated peptides were recovered from non-immunoglobulin somatically mutated genes. Somatic mutations within the immunoglobulin variable region were almost exclusively presented by MHC class II. We isolated circulating CD4(+) T cells specific for immunoglobulin-derived neoantigens and found these cells could mediate killing of autologous lymphoma cells. These results demonstrate that an integrative approach combining MHC isolation, peptide identification, and exome sequencing is an effective platform to uncover tumour neoantigens. Application of this strategy to human lymphoma implicates immunoglobulin neoantigens as targets for lymphoma immunotherapy.
View details for DOI 10.1038/nature21433
View details for PubMedID 28329770
-
Data normalization considerations for digital tumor dissection.
Genome biology
2017; 18 (1): 128
Abstract
In a recently published article in Genome Biology, Li and colleagues introduced TIMER, a gene expression deconvolution approach for studying tumor-infiltrating leukocytes (TILs) in 23 cancer types profiled by The Cancer Genome Atlas. Methods to characterize TIL biology are increasingly important, and the authors offer several arguments in favor of their strategy. Several of these claims warrant further discussion and highlight the critical importance of data normalization in gene expression deconvolution applications.Please see related Li et al correspondence: www.dx.doi.org/10.1186/s13059-017-1256-5 and Zheng correspondence: www.dx.doi.org/10.1186/s13059-017-1258-3.
View details for PubMedID 28679399
-
Targeted chromatin ligation, a robust epigenetic profiling technique for small cell numbers.
Nucleic acids research
2017; 45 (17): e153
Abstract
The complexity and inefficiency of chromatin immunoprecipitation strategies restrict their sensitivity and application when examining rare cell populations. We developed a new technique that replaces immunoprecipitation with a simplified chromatin fragmentation and proximity ligation step that eliminates bead purification and washing steps. We present a simple single tube proximity ligation technique, targeted chromatin ligation, that captures histone modification patterns with only 200 cells. Our technique eliminates loss of material and sensitivity due to multiple inefficient steps, while simplifying the workflow to enhance sensitivity and create the potential for novel applications.
View details for PubMedID 28973448
-
Early detection of molecular residual disease in localized lung cancer by circulating tumor DNA profiling.
Cancer discovery
2017
Abstract
Identifying molecular residual disease (MRD) after treatment of localized lung cancer could facilitate early intervention and personalization of adjuvant therapies. Here we apply Cancer Personalized Profiling by Deep Sequencing (CAPP-Seq) circulating tumor DNA (ctDNA) analysis to 255 samples from 40 patients treated with curative intent for stage I-III lung cancer and 54 healthy adults. In 94% of evaluable patients experiencing recurrence, ctDNA was detectable in the first post-treatment blood sample, indicating reliable identification of MRD. Post-treatment ctDNA detection preceded radiographic progression in 72% of patients by a median of 5.2 months and 53% of patients harbored ctDNA mutation profiles associated with favorable responses to tyrosine kinase inhibitors or immune checkpoint blockade. Collectively, these results indicate that ctDNA MRD in lung cancer patients can be accurately detected using CAPP-Seq and may allow personalized adjuvant treatment while disease burden is lowest.
View details for PubMedID 28899864
-
Distinct biological subtypes and patterns of genome evolution in lymphoma revealed by circulating tumor DNA
SCIENCE TRANSLATIONAL MEDICINE
2016; 8 (364)
Abstract
Patients with diffuse large B cell lymphoma (DLBCL) exhibit marked diversity in tumor behavior and outcomes, yet the identification of poor-risk groups remains challenging. In addition, the biology underlying these differences is incompletely understood. We hypothesized that characterization of mutational heterogeneity and genomic evolution using circulating tumor DNA (ctDNA) profiling could reveal molecular determinants of adverse outcomes. To address this hypothesis, we applied cancer personalized profiling by deep sequencing (CAPP-Seq) analysis to tumor biopsies and cell-free DNA samples from 92 lymphoma patients and 24 healthy subjects. At diagnosis, the amount of ctDNA was found to strongly correlate with clinical indices and was independently predictive of patient outcomes. We demonstrate that ctDNA genotyping can classify transcriptionally defined tumor subtypes, including DLBCL cell of origin, directly from plasma. By simultaneously tracking multiple somatic mutations in ctDNA, our approach outperformed immunoglobulin sequencing and radiographic imaging for the detection of minimal residual disease and facilitated noninvasive identification of emergent resistance mutations to targeted therapies. In addition, we identified distinct patterns of clonal evolution distinguishing indolent follicular lymphomas from those that transformed into DLBCL, allowing for potential noninvasive prediction of histological transformation. Collectively, our results demonstrate that ctDNA analysis reveals biological factors that underlie lymphoma clinical outcomes and could facilitate individualized therapy.
View details for DOI 10.1126/scitranslmed.aai8545
View details for PubMedID 27831904
-
Role of KEAP1/NRF2 and TP53 Mutations in Lung Squamous Cell Carcinoma Development and Radiation Resistance.
Cancer discovery
2016
Abstract
Lung squamous cell carcinoma (LSCC) pathogenesis remains incompletely understood, and biomarkers predicting treatment response remain lacking. Here, we describe novel murine LSCC models driven by loss of Trp53 and Keap1, both of which are frequently mutated in human LSCCs. Homozygous inactivation of Keap1 or Trp53 promoted airway basal stem cell (ABSC) self-renewal, suggesting that mutations in these genes lead to expansion of mutant stem cell clones. Deletion of Trp53 and Keap1 in ABSCs, but not more differentiated tracheal cells, produced tumors recapitulating histologic and molecular features of human LSCCs, indicating that they represent the likely cell of origin in this model. Deletion of Keap1 promoted tumor aggressiveness, metastasis, and resistance to oxidative stress and radiotherapy (RT). KEAP1/NRF2 mutation status predicted risk of local recurrence after RT in patients with non-small lung cancer (NSCLC) and could be noninvasively identified in circulating tumor DNA. Thus, KEAP1/NRF2 mutations could serve as predictive biomarkers for personalization of therapeutic strategies for NSCLCs.We developed an LSCC mouse model involving Trp53 and Keap1, which are frequently mutated in human LSCCs. In this model, ABSCs are the cell of origin of these tumors. KEAP1/NRF2 mutations increase radioresistance and predict local tumor recurrence in radiotherapy patients. Our findings are of potential clinical relevance and could lead to personalized treatment strategies for tumors with KEAP1/NRF2 mutations. Cancer Discov; 7(1); 86-101. ©2016 AACR.This article is highlighted in the In This Issue feature, p. 1.
View details for PubMedID 27663899
-
High-throughput genomic profiling of tumor-infiltrating leukocytes.
Current opinion in immunology
2016; 41: 77-84
Abstract
Tumors are complex ecosystems comprised of diverse cell types including malignant cells, mesenchymal cells, and tumor-infiltrating leukocytes (TILs). While TILs are well known to play important roles in many aspects of cancer biology, recent developments in immuno-oncology have spurred considerable interest in TILs, particularly in relation to their optimal engagement by emerging immunotherapies. Traditionally, the enumeration of TIL phenotypic diversity and composition in solid tumors has relied on resolving single cells by flow cytometry and immunohistochemical methods. However, advances in genome-wide technologies and computational methods are now allowing TILs to be profiled with increasingly high resolution and accuracy directly from RNA mixtures of bulk tumor samples. In this review, we highlight recent progress in the development of in silico tumor dissection methods, and illustrate examples of how these strategies can be applied to characterize TILs in human tumors to facilitate personalized cancer therapy.
View details for DOI 10.1016/j.coi.2016.06.006
View details for PubMedID 27372732
-
Circulating tumour DNA profiling reveals heterogeneity of EGFR inhibitor resistance mechanisms in lung cancer patients
NATURE COMMUNICATIONS
2016; 7
Abstract
Circulating tumour DNA (ctDNA) analysis facilitates studies of tumour heterogeneity. Here we employ CAPP-Seq ctDNA analysis to study resistance mechanisms in 43 non-small cell lung cancer (NSCLC) patients treated with the third-generation epidermal growth factor receptor (EGFR) inhibitor rociletinib. We observe multiple resistance mechanisms in 46% of patients after treatment with first-line inhibitors, indicating frequent intra-patient heterogeneity. Rociletinib resistance recurrently involves MET, EGFR, PIK3CA, ERRB2, KRAS and RB1. We describe a novel EGFR L798I mutation and find that EGFR C797S, which arises in ∼33% of patients after osimertinib treatment, occurs in <3% after rociletinib. Increased MET copy number is the most frequent rociletinib resistance mechanism in this cohort and patients with multiple pre-existing mechanisms (T790M and MET) experience inferior responses. Similarly, rociletinib-resistant xenografts develop MET amplification that can be overcome with the MET inhibitor crizotinib. These results underscore the importance of tumour heterogeneity in NSCLC and the utility of ctDNA-based resistance mechanism assessment.
View details for DOI 10.1038/ncomms11815
View details for PubMedID 27283993
-
Identification of tumorigenic cells and therapeutic targets in pancreatic neuroendocrine tumors
PROCEEDINGS OF THE NATIONAL ACADEMY OF SCIENCES OF THE UNITED STATES OF AMERICA
2016; 113 (16): 4464-4469
Abstract
Pancreatic neuroendocrine tumors (PanNETs) are a type of pancreatic cancer with limited therapeutic options. Consequently, most patients with advanced disease die from tumor progression. Current evidence indicates that a subset of cancer cells is responsible for tumor development, metastasis, and recurrence, and targeting these tumor-initiating cells is necessary to eradicate tumors. However, tumor-initiating cells and the biological processes that promote pathogenesis remain largely uncharacterized in PanNETs. Here we profile primary and metastatic tumors from an index patient and demonstrate that MET proto-oncogene activation is important for tumor growth in PanNET xenograft models. We identify a highly tumorigenic cell population within several independent surgically acquired PanNETs characterized by increased cell-surface protein CD90 expression and aldehyde dehydrogenase A1 (ALDHA1) activity, and provide in vitro and in vivo evidence for their stem-like properties. We performed proteomic profiling of 332 antigens in two cell lines and four primary tumors, and showed that CD47, a cell-surface protein that acts as a "don't eat me" signal co-opted by cancers to evade innate immune surveillance, is ubiquitously expressed. Moreover, CD47 coexpresses with MET and is enriched in CD90(hi)cells. Furthermore, blocking CD47 signaling promotes engulfment of tumor cells by macrophages in vitro and inhibits xenograft tumor growth, prevents metastases, and prolongs survival in vivo.
View details for DOI 10.1073/pnas.1600007113
View details for PubMedID 27035983
-
Skin fibrosis. Identification and isolation of a dermal lineage with intrinsic fibrogenic potential.
Science
2015; 348 (6232)
Abstract
Dermal fibroblasts represent a heterogeneous population of cells with diverse features that remain largely undefined. We reveal the presence of at least two fibroblast lineages in murine dorsal skin. Lineage tracing and transplantation assays demonstrate that a single fibroblast lineage is responsible for the bulk of connective tissue deposition during embryonic development, cutaneous wound healing, radiation fibrosis, and cancer stroma formation. Lineage-specific cell ablation leads to diminished connective tissue deposition in wounds and reduces melanoma growth. Using flow cytometry, we identify CD26/DPP4 as a surface marker that allows isolation of this lineage. Small molecule-based inhibition of CD26/DPP4 enzymatic activity during wound healing results in diminished cutaneous scarring. Identification and isolation of these lineages hold promise for translational medicine aimed at in vivo modulation of fibrogenic behavior.
View details for DOI 10.1126/science.aaa2151
View details for PubMedID 25883361
View details for PubMedCentralID PMC5088503
-
Potential clinical utility of ultrasensitive circulating tumor DNA detection with CAPP-Seq.
Expert review of molecular diagnostics
2015: 1–5
Abstract
Tumors continually shed DNA into the circulation, where it can be noninvasively accessed. The ability to accurately detect circulating tumor DNA (ctDNA) could significantly impact the management of patients with nearly every cancer type. Quantitation of ctDNA could allow objective response assessment, detection of minimal residual disease and noninvasive tumor genotyping. The latter application overcomes the barriers currently limiting repeated tumor tissue sampling during therapy. Recent technical advancements have improved upon the sensitivity, specificity and feasibility of ctDNA detection and promise to enable innovative clinical applications. Here, we focus on the potential clinical utility of ctDNA analysis using CAncer Personalized Profiling by deep Sequencing (CAPP-Seq), a novel next-generation sequencing-based approach for ultrasensitive ctDNA detection. Applications of CAPP-Seq for the personalization of cancer detection and therapy are discussed.
View details for DOI 10.1586/14737159.2015.1019476
View details for PubMedID 25773944
-
Large-Scale and Comprehensive Immune Profiling and Functional Analysis of Normal Human Aging.
PloS one
2015; 10 (7)
Abstract
While many age-associated immune changes have been reported, a comprehensive set of metrics of immune aging is lacking. Here we report data from 243 healthy adults aged 40-97, for whom we measured clinical and functional parameters, serum cytokines, cytokines and gene expression in stimulated and unstimulated PBMC, PBMC phenotypes, and cytokine-stimulated pSTAT signaling in whole blood. Although highly heterogeneous across individuals, many of these assays revealed trends by age, sex, and CMV status, to greater or lesser degrees. Age, then sex and CMV status, showed the greatest impact on the immune system, as measured by the percentage of assay readouts with significant differences. An elastic net regression model could optimally predict age with 14 analytes from different assays. This reinforces the importance of multivariate analysis for defining a healthy immune system. These data provide a reference for others measuring immune parameters in older people.
View details for DOI 10.1371/journal.pone.0133627
View details for PubMedID 26197454
-
In Vivo clonal analysis reveals lineage-restricted progenitor characteristics in Mammalian kidney development, maintenance, and regeneration.
Cell reports
2014; 7 (4): 1270-1283
Abstract
The mechanism and magnitude by which the mammalian kidney generates and maintains its proximal tubules, distal tubules, and collecting ducts remain controversial. Here, we use long-term in vivo genetic lineage tracing and clonal analysis of individual cells from kidneys undergoing development, maintenance, and regeneration. We show that the adult mammalian kidney undergoes continuous tubulogenesis via expansions of fate-restricted clones. Kidneys recovering from damage undergo tubulogenesis through expansions of clones with segment-specific borders, and renal spheres developing in vitro from individual cells maintain distinct, segment-specific fates. Analysis of mice derived by transfer of color-marked embryonic stem cells (ESCs) into uncolored blastocysts demonstrates that nephrons are polyclonal, developing from expansions of singly fated clones. Finally, we show that adult renal clones are derived from Wnt-responsive precursors, and their tracing in vivo generates tubules that are segment specific. Collectively, these analyses demonstrate that fate-restricted precursors functioning as unipotent progenitors continuously maintain and self-preserve the mouse kidney throughout life.
View details for DOI 10.1016/j.celrep.2014.04.018
View details for PubMedID 24835991
View details for PubMedCentralID PMC4425291
-
Efficient Selection of Biomineralizing DNA Aptamers Using Deep Sequencing and Population Clustering
ACS Nano
2014; 8 (1): 387-395
View details for DOI 10.1021/nn404448s
-
FACTERA: a practical method for the discovery of genomic rearrangements at breakpoint resolution
BIOINFORMATICS
2014
View details for DOI 10.1093/bioinformatics/btu549
-
Identifying Stem Cell Gene Expression Patterns and Phenotypic Networks with AutoSOME.
Methods in molecular biology (Clifton, N.J.)
2014; 1150: 115-130
Abstract
Stem cells have the unique property of differentiation and self-renewal and play critical roles in normal development, tissue repair, and disease. To promote systems-wide analysis of cells and tissues, we developed AutoSOME, a machine-learning method for identifying coordinated gene expression patterns and correlated cellular phenotypes in whole-transcriptome data, without prior knowledge of cluster number or structure. Here, we present a facile primer demonstrating the use of AutoSOME for identification and characterization of stem cell gene expression signatures and for visualization of transcriptome networks using Cytoscape. This protocol should serve as a general foundation for gene expression cluster analysis of stem cells, with applications for studying pluripotency, multi-lineage potential, and neoplastic disease.
View details for DOI 10.1007/978-1-4939-0512-6_6
View details for PubMedID 24743993
-
The genome sequence of the colonial chordate, Botryllus schlosseri.
eLife
2013; 2
Abstract
Botryllus schlosseri is a colonial urochordate that follows the chordate plan of development following sexual reproduction, but invokes a stem cell-mediated budding program during subsequent rounds of asexual reproduction. As urochordates are considered to be the closest living invertebrate relatives of vertebrates, they are ideal subjects for whole genome sequence analyses. Using a novel method for high-throughput sequencing of eukaryotic genomes, we sequenced and assembled 580 Mbp of the B. schlosseri genome. The genome assembly is comprised of nearly 14,000 intron-containing predicted genes, and 13,500 intron-less predicted genes, 40% of which could be confidently parceled into 13 (of 16 haploid) chromosomes. A comparison of homologous genes between B. schlosseri and other diverse taxonomic groups revealed genomic events underlying the evolution of vertebrates and lymphoid-mediated immunity. The B. schlosseri genome is a community resource for studying alternative modes of reproduction, natural transplantation reactions, and stem cell-mediated regeneration. DOI:http://dx.doi.org/10.7554/eLife.00569.001.
View details for DOI 10.7554/eLife.00569
View details for PubMedID 23840927
-
Systems-level analysis of age-related macular degeneration reveals global biomarkers and phenotype-specific functional networks
GENOME MEDICINE
2012; 4
Abstract
Please see related commentary: http://www.biomedcentral.com/1741-7015/10/21/abstractAge-related macular degeneration (AMD) is a leading cause of blindness that affects the central region of the retinal pigmented epithelium (RPE), choroid, and neural retina. Initially characterized by an accumulation of sub-RPE deposits, AMD leads to progressive retinal degeneration, and in advanced cases, irreversible vision loss. Although genetic analysis, animal models, and cell culture systems have yielded important insights into AMD, the molecular pathways underlying AMD's onset and progression remain poorly delineated. We sought to better understand the molecular underpinnings of this devastating disease by performing the first comparative transcriptome analysis of AMD and normal human donor eyes.RPE-choroid and retina tissue samples were obtained from a common cohort of 31 normal, 26 AMD, and 11 potential pre-AMD human donor eyes. Transcriptome profiles were generated for macular and extramacular regions, and statistical and bioinformatic methods were employed to identify disease-associated gene signatures and functionally enriched protein association networks. Selected genes of high significance were validated using an independent donor cohort.We identified over 50 annotated genes enriched in cell-mediated immune responses that are globally over-expressed in RPE-choroid AMD phenotypes. Using a machine learning model and a second donor cohort, we show that the top 20 global genes are predictive of AMD clinical diagnosis. We also discovered functionally enriched gene sets in the RPE-choroid that delineate the advanced AMD phenotypes, neovascular AMD and geographic atrophy. Moreover, we identified a graded increase of transcript levels in the retina related to wound response, complement cascade, and neurogenesis that strongly correlates with decreased levels of phototransduction transcripts and increased AMD severity. Based on our findings, we assembled protein-protein interactomes that highlight functional networks likely to be involved in AMD pathogenesis.We discovered new global biomarkers and gene expression signatures of AMD. These results are consistent with a model whereby cell-based inflammatory responses represent a central feature of AMD etiology, and depending on genetics, environment, or stochastic factors, may give rise to the advanced AMD phenotypes characterized by angiogenesis and/or cell death. Genes regulating these immunological activities, along with numerous other genes identified here, represent promising new targets for AMD-directed therapeutics and diagnostics.
View details for DOI 10.1186/PREACCEPT-1418491035586234
View details for Web of Science ID 000314566500002
View details for PubMedID 22364233
-
A proteomic approach for the identification of novel lysine methyltransferase substrates
EPIGENETICS & CHROMATIN
2011; 4
Abstract
Signaling via protein lysine methylation has been proposed to play a central role in the regulation of many physiologic and pathologic programs. In contrast to other post-translational modifications such as phosphorylation, proteome-wide approaches to investigate lysine methylation networks do not exist.In the current study, we used the ProtoArray® platform, containing over 9,500 human proteins, and developed and optimized a system for proteome-wide identification of novel methylation events catalyzed by the protein lysine methyltransferase (PKMT) SETD6. This enzyme had previously been shown to methylate the transcription factor RelA, but it was not known whether SETD6 had other substrates. By using two independent detection approaches, we identified novel candidate substrates for SETD6, and verified that all targets tested in vitro and in cells were genuine substrates.We describe a novel proteome-wide methodology for the identification of new PKMT substrates. This technological advance may lead to a better understanding of the enzymatic activity and substrate specificity of the large number (more than 50) PKMTs present in the human proteome, most of which are uncharacterized.
View details for DOI 10.1186/1756-8935-4-19
View details for PubMedID 22024134
-
Global Analysis of Proline-Rich Tandem Repeat Proteins Reveals Broad Phylogenetic Diversity in Plant Secretomes
PLOS ONE
2011; 6 (8)
Abstract
Cell walls, constructed by precisely choreographed changes in the plant secretome, play critical roles in plant cell physiology and development. Along with structural polysaccharides, secreted proline-rich Tandem Repeat Proteins (TRPs) are important for cell wall function, yet the evolutionary diversity of these structural TRPs remains virtually unexplored. Using a systems-level computational approach to analyze taxonomically diverse plant sequence data, we identified 31 distinct Pro-rich TRP classes targeted for secretion. This analysis expands upon the known phylogenetic diversity of extensins, the most widely studied class of wall structural proteins, and demonstrates that extensins evolved before plant vascularization. Our results also show that most Pro-rich TRP classes have unexpectedly restricted evolutionary distributions, revealing considerable differences in plant secretome signatures that define unexplored diversity.
View details for DOI 10.1371/journal.pone.0023167
View details for Web of Science ID 000293511900032
View details for PubMedID 21829715
-
clusterMaker: a multi-algorithm clustering plugin for Cytoscape
BMC BIOINFORMATICS
2011; 12: 436
View details for DOI 10.1186/1471-2105-12-436
-
AutoSOME: a clustering method for identifying gene expression modules without prior knowledge of cluster number
BMC BIOINFORMATICS
2010; 11
Abstract
Clustering the information content of large high-dimensional gene expression datasets has widespread application in "omics" biology. Unfortunately, the underlying structure of these natural datasets is often fuzzy, and the computational identification of data clusters generally requires knowledge about cluster number and geometry.We integrated strategies from machine learning, cartography, and graph theory into a new informatics method for automatically clustering self-organizing map ensembles of high-dimensional data. Our new method, called AutoSOME, readily identifies discrete and fuzzy data clusters without prior knowledge of cluster number or structure in diverse datasets including whole genome microarray data. Visualization of AutoSOME output using network diagrams and differential heat maps reveals unexpected variation among well-characterized cancer cell lines. Co-expression analysis of data from human embryonic and induced pluripotent stem cells using AutoSOME identifies >3400 up-regulated genes associated with pluripotency, and indicates that a recently identified protein-protein interaction network characterizing pluripotency was underestimated by a factor of four.By effectively extracting important information from high-dimensional microarray data without prior knowledge or the need for data filtration, AutoSOME can yield systems-level insights from whole genome microarray expression studies. Due to its generality, this new method should also have practical utility for a variety of data-intensive applications, including the results of deep sequencing experiments. AutoSOME is available for download at http://jimcooperlab.mcdb.ucsb.edu/autosome webcite.
View details for DOI 10.1186/1471-2105-11-117
View details for Web of Science ID 000276296100002
View details for PubMedID 20202218
-
XSTREAM: A practical algorithm for identification and architecture modeling of tandem repeats in protein sequences
BMC BIOINFORMATICS
2007; 8
Abstract
Biological sequence repeats arranged in tandem patterns are widespread in DNA and proteins. While many software tools have been designed to detect DNA tandem repeats (TRs), useful algorithms for identifying protein TRs with varied levels of degeneracy are still needed.To address limitations of current repeat identification methods, and to provide an efficient and flexible algorithm for the detection and analysis of TRs in protein sequences, we designed and implemented a new computational method called XSTREAM. Running time tests confirm the practicality of XSTREAM for analyses of multi-genome datasets. Each of the key capabilities of XSTREAM (e.g., merging, nesting, long-period detection, and TR architecture modeling) are demonstrated using anecdotal examples, and the utility of XSTREAM for identifying TR proteins was validated using data from a recently published paper.We show that XSTREAM is a practical and valuable tool for TR detection in protein and nucleotide sequences at the multi-genome scale, and an effective tool for modeling TR domains with diverse architectures and varied levels of degeneracy. Because of these useful features, XSTREAM has significant potential for the discovery of naturally-evolved modular proteins with applications for engineering novel biostructural and biomimetic materials, and identifying new vaccine and diagnostic targets.
View details for DOI 10.1186/1471-2105-8-382
View details for Web of Science ID 000252936900001
View details for PubMedID 17931424