Debashis Sahoo

Single-cell dissection of transcriptional heterogeneity in human colon tumors.

Cancer is often viewed as a caricature of normal developmental processes, but the extent to which its cellular heterogeneity truly recapitulates multilineage differentiation processes of normal tissues remains unknown. Here we implement single-cell PCR gene-expression analysis to dissect the cellular composition of primary human normal colon and colon cancer epithelia. We show that human colon cancer tissues contain distinct cell populations whose transcriptional identities mirror those of the different cellular lineages of normal colon. By creating monoclonal tumor xenografts from injection of a single (n = 1) cell, we demonstrate that the transcriptional diversity of cancer tissues is largely explained by in vivo multilineage differentiation and not only by clonal genetic heterogeneity. Finally, we show that the different gene-expression programs linked to multilineage differentiation are strongly associated with patient survival. We develop two-gene classifier systems (KRT20 versus CA1, MS4A12, CD177, SLC26A3) that predict clinical outcomes with hazard ratios superior to those of pathological grade and comparable to those of microarray-derived multigene expression signatures.

Human bone marrow hematopoietic stem cells are increased in frequency and myeloid-biased with age

In the human hematopoietic system, aging is associated with decreased bone marrow cellularity, decreased adaptive immune system function, and increased incidence of anemia and other hematological disorders and malignancies. Recent studies in mice suggest that changes within the hematopoietic stem cell (HSC) population during aging contribute significantly to the manifestation of these age-associated hematopoietic pathologies. Though the mouse HSC population has been shown to change both quantitatively and functionally with age, changes in the human HSC and progenitor cell populations during aging have been incompletely characterized. To elucidate the properties of an aged human hematopoietic system that may predispose to age-associated hematopoietic dysfunction, we evaluated immunophenotypic HSC and other hematopoietic progenitor populations from healthy, hematologically normal young and elderly human bone marrow samples. We found that aged immunophenotypic human HSC increase in frequency, are less quiescent, and exhibit myeloid-biased differentiation potential compared with young HSC. Gene expression profiling revealed that aged immunophenotypic human HSC transcriptionally up-regulate genes associated with cell cycle, myeloid lineage specification, and myeloid malignancies. These age-associated alterations in the frequency, developmental potential, and gene expression profile of human HSC are similar to those changes observed in mouse HSC, suggesting that hematopoietic aging is an evolutionarily conserved process.

An antibody against SSEA-5 glycan on human pluripotent stem cells enables removal of teratoma-forming cells

An important risk in the clinical application of human pluripotent stem cells (hPSCs), including human embryonic and induced pluripotent stem cells (hESCs and hiPSCs), is teratoma formation by residual undifferentiated cells. We raised a monoclonal antibody against hESCs, designated anti–stage-specific embryonic antigen (SSEA)-5, which binds a previously unidentified antigen highly and specifically expressed on hPSCs—the H type-1 glycan. Separation based on SSEA-5 expression through fluorescence-activated cell sorting (FACS) greatly reduced teratoma-formation potential of heterogeneously differentiated cultures. To ensure complete removal of teratoma-forming cells, we identified additional pluripotency surface markers (PSMs) exhibiting a large dynamic expression range during differentiation: CD9, CD30, CD50, CD90 and CD200. Immunohistochemistry studies of human fetal tissues and bioinformatics analysis of a microarray database revealed that concurrent expression of these markers is both common and specific to hPSCs. Immunodepletion with antibodies against SSEA-5 and two additional PSMs completely removed teratoma-formation potential from incompletely differentiated hESC cultures.

Ly6d marks the earliest stage of B-cell specification and identifies the branchpoint between B-cell and T-cell development

Common lymphoid progenitors (CLPs) clonally produce both B- and T-cell lineages, but have little myeloid potential in vivo. However, some studies claim that the upstream lymphoid-primed multipotent progenitor (LMPP) is the thymic seeding population, and suggest that CLPs are primarily B-cell-restricted. To identify surface proteins that distinguish functional CLPs from B-cell progenitors, we used a new computational method of Mining Developmentally Regulated Genes (MiDReG). We identified Ly6d, which divides CLPs into two distinct populations: one that retains full in vivo lymphoid potential and produces more thymocytes at early timepoints than LMPP, and another that behaves essentially as a B-cell progenitor.

MicroRNA-125b expands hematopoietic stem cells and enriches for the lymphoid-balanced and lymphoid-biased subsets

MicroRNAs profoundly impact hematopoietic cells by regulating progenitor cell-fate decisions, as well as mature immune effector function. However to date, microRNAs that regulate hematopoietic stem cell (HSC) function have been less well characterized. Here we show that microRNA-125b (miR-125b) is highly expressed in HSCs and its expression decreases in committed progenitors. Overexpression of miR-125b in mouse HSC enhances their function, demonstrated through serial transplantation of highly purified HSC, and enriches for the previously described Slamf1loCD34− lymphoid-balanced and the Slamf1negCD34− lymphoid-biased cell subsets within the multipotent HSC (CD34-KLS) fraction. Mature peripheral blood cells derived from the miR-125b–overexpressing HSC are skewed toward the lymphoid lineage. Consistent with this observation, miR-125b overexpression significantly increases the number of early B-progenitor cells within the spleen and induces the expansion and enrichment of the lymphoid-balanced and lymphoid-biased HSC subset via an antiapoptotic mechanism, reducing the mRNA expression levels of two proapoptotic targets, Bmf and KLF13. The antiapoptotic effect of miR-125b is more pronounced in the lymphoid-biased HSC subset because of their intrinsic higher baseline levels of apoptosis. These effects of miR-125b are associated with the development of lymphoproliferative disease, marked by expansion of CD8+ T lymphocytes. Taken together, these data reveal that miR-125b regulates HSC survival and can promote lymphoid-fate decisions at the level of the HSC by preferentially expanding lymphoid-balanced and lymphoid-biased HSC.

The genome sequence of the colonial chordate, Botryllus schlosseri

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

Gene Expression Commons: An Open Platform for Absolute Gene Expression Profiling

Gene expression profiling using microarrays has been limited to comparisons of gene expression between small numbers of samples within individual experiments. However, the unknown and variable sensitivities of each probeset have rendered the absolute expression of any given gene nearly impossible to estimate. We have overcome this limitation by using a very large number (>10,000) of varied microarray data as a common reference, so that statistical attributes of each probeset, such as the dynamic range and threshold between low and high expression, can be reliably discovered through meta-analysis. This strategy is implemented in a web-based platform named “Gene Expression Commons” (https://gexc.stanford.edu/) which contains data of 39 distinct highly purified mouse hematopoietic stem/progenitor/differentiated cell populations covering almost the entire hematopoietic system. Since the Gene Expression Commons is designed as an open platform, investigators can explore the expression level of any gene, search by expression patterns of interest, submit their own microarray data, and design their own working models representing biological relationship among samples.

CDX2 as a Prognostic Biomarker in Stage II and Stage III Colon Cancer

Background The identification of high-risk stage II colon cancers is key to the selection of patients who require adjuvant treatment after surgery. Microarray-based multigene-expression signatures derived from stem cells and progenitor cells hold promise, but they are difficult to use in clinical practice. Methods We used a new bioinformatics approach to search for biomarkers of colon epithelial differentiation across gene-expression arrays and then ranked candidate genes according to the availability of clinical-grade diagnostic assays. With the use of subgroup analysis involving independent and retrospective cohorts of patients with stage II or stage III colon cancer, the top candidate gene was tested for its association with disease-free survival and a benefit from adjuvant chemotherapy. Results The transcription factor CDX2 ranked first in our screening test. A group of 87 of 2115 tumor samples (4.1%) lacked CDX2 expression. In the discovery data set, which included 466 patients, the rate of 5-year disease-free survival was lower among the 32 patients (6.9%) with CDX2-negative colon cancers than among the 434 (93.1%) with CDX2-positive colon cancers (hazard ratio for disease recurrence, 3.44; 95% confidence interval [CI], 1.60 to 7.38; P=0.002). In the validation data set, which included 314 patients, the rate of 5-year disease-free survival was lower among the 38 patients (12.1%) with CDX2 protein-negative colon cancers than among the 276 (87.9%) with CDX2 protein-positive colon cancers (hazard ratio, 2.42; 95% CI, 1.36 to 4.29; P=0.003). In both these groups, these findings were independent of the patients age, sex, and tumor stage and grade. Among patients with stage II cancer, the difference in 5-year disease-free survival was significant both in the discovery data set (49% among 15 patients with CDX2-negative tumors vs. 87% among 191 patients with CDX2-positive tumors, P=0.003) and in the validation data set (51% among 15 patients with CDX2-negative tumors vs. 80% among 106 patients with CDX2-positive tumors, P=0.004). In a pooled database of all patient cohorts, the rate of 5-year disease-free survival was higher among 23 patients with stage II CDX2-negative tumors who were treated with adjuvant chemotherapy than among 25 who were not treated with adjuvant chemotherapy (91% vs. 56%, P=0.006). Conclusions Lack of CDX2 expression identified a subgroup of patients with high-risk stage II colon cancer who appeared to benefit from adjuvant chemotherapy. (Funded by the National Comprehensive Cancer Network, the National Institutes of Health, and others.).

Three differentiation states risk-stratify bladder cancer into distinct subtypes

Current clinical judgment in bladder cancer (BC) relies primarily on pathological stage and grade. We investigated whether a molecular classification of tumor cell differentiation, based on a developmental biology approach, can provide additional prognostic information. Exploiting large preexisting gene-expression databases, we developed a biologically supervised computational model to predict markers that correspond with BC differentiation. To provide mechanistic insight, we assessed relative tumorigenicity and differentiation potential via xenotransplantation. We then correlated the prognostic utility of the identified markers to outcomes within gene expression and formalin-fixed paraffin-embedded (FFPE) tissue datasets. Our data indicate that BC can be subclassified into three subtypes, on the basis of their differentiation states: basal, intermediate, and differentiated, where only the most primitive tumor cell subpopulation within each subtype is capable of generating xenograft tumors and recapitulating downstream populations. We found that keratin 14 (KRT14) marks the most primitive differentiation state that precedes KRT5 and KRT20 expression. Furthermore, KRT14 expression is consistently associated with worse prognosis in both univariate and multivariate analyses. We identify here three distinct BC subtypes on the basis of their differentiation states, each harboring a unique tumor-initiating population.

Hoxb5 marks long-term haematopoietic stem cells and reveals a homogenous perivascular niche

Haematopoietic stem cells (HSCs) are arguably the most extensively characterized tissue stem cells. Since the identification of HSCs by prospective isolation, complex multi-parameter flow cytometric isolation of phenotypic subsets has facilitated studies on many aspects of HSC biology, including self-renewal, differentiation, ageing, niche, and diversity. Here we demonstrate by unbiased multi-step screening, identification of a single gene, homeobox B5 (Hoxb5, also known as Hox-2.1), with expression in the bone marrow that is limited to long-term (LT)-HSCs in mice. Using a mouse single-colour tri-mCherry reporter driven by endogenous Hoxb5 regulation, we show that only the Hoxb5+ HSCs exhibit long-term reconstitution capacity after transplantation in primary transplant recipients and, notably, in secondary recipients. Only 7–35% of various previously defined immunophenotypic HSCs are LT-HSCs. Finally, by in situ imaging of mouse bone marrow, we show that >94% of LT-HSCs (Hoxb5+) are directly attached to VE-cadherin+ cells, implicating the perivascular space as a near-homogenous location of LT-HSCs.