Knowledge of immune cell phenotypes in the tumor microenvironment is essential for understanding mechanisms of cancer progression and immunotherapy response. We profiled 45,000 immune cells from eight breast carcinomas, as well as matched normal breast tissue, blood, and lymph nodes, using single-cell RNA-seq. We developed a preprocessing pipeline, SEQC, and a Bayesian clustering and normalization method, Biscuit, to address computational challenges inherent to single-cell data. Despite significant similarity between normal and tumor tissue-resident immune cells, we observed continuous phenotypic expansions specific to the tumor microenvironment. Analysis of paired single-cell RNA and T cell receptor (TCR) sequencing data from 27,000 additional T cells revealed the combinatorial impact of TCR utilization on phenotypic diversity. Our results support a model of continuous activation in T cells and do not comport with the macrophage polarization model in cancer. Our results have important implications for characterizing tumor-infiltrating immune cells.

We have taken the first steps towards a complete reconstruction of the Mycobacterium tuberculosis regulatory network based on ChIP-Seq and combined this reconstruction with system-wide profiling of messenger RNAs, proteins, metabolites and lipids during hypoxia and re-aeration. Adaptations to hypoxia are thought to have a prominent role in M. tuberculosis pathogenesis. Using ChIP-Seq combined with expression data from the induction of the same factors, we have reconstructed a draft regulatory network based on 50 transcription factors. This network model revealed a direct interconnection between the hypoxic response, lipid catabolism, lipid anabolism and the production of cell wall lipids. As a validation of this model, in response to oxygen availability we observe substantial alterations in lipid content and changes in gene expression and metabolites in corresponding metabolic pathways. The regulatory network reveals transcription factors underlying these changes, allows us to computationally predict expression changes, and indicates that Rv0081 is a regulatory hub.

Tissue damage increases cancer risk through poorly understood mechanisms 1 . In the pancreas, pancreatitis associated with tissue injury collaborates with activating mutations in the Kras oncogene to dramatically accelerate the formation of early neoplastic lesions and ultimately pancreatic cancer 2 , 3 . By integrating genomics, single-cell chromatin assays and spatiotemporally-controlled functional perturbations in autochthonous mouse models, we show that the combination of Kras mutation and tissue damage promotes a unique chromatin state in the pancreatic epithelium that distinguishes neoplastic transformation from normal regeneration and is selected for throughout malignant evolution. This cancer-associated epigenetic state emerges within 48 hours of pancreatic injury, and involves an acinar-to-neoplasia ‘chromatin switch’ that contributes to the early dysregulation of genes defining human pancreatic cancer. Among the genes most rapidly activated upon tissue damage in the pre-malignant pancreatic epithelium is the alarmin cytokine IL-33, which cooperates with mutant Kras in unleashing the epigenetic remodeling program of early neoplasia and neoplastic transformation in the absence of injury. Collectively, our study demonstrates how gene-environment interactions can rapidly produce gene regulatory programs that dictate early neoplastic commitment and provides a molecular framework for understanding the interplay between genetics and environmental cues in cancer initiation.

Transcriptional regulators, including the cohesin complex member STAG2, are recurrently mutated in cancer. The role of STAG2 in gene regulation, hematopoiesis, and tumor suppression remains unresolved. We show that Stag2 deletion in hematopoietic stem and progenitor cells (HSPCs) results in altered hematopoietic function, increased selfrenewal, and impaired differentiation. Chromatin immunoprecipitation (ChIP) sequencing revealed that, although Stag2 and Stag1 bind a shared set of genomic loci, a component of Stag2 binding sites is unoccupied by Stag1, even in Stag2-deficient HSPCs. Although concurrent loss of Stag2 and Stag1 abrogated hematopoiesis, Stag2 loss alone decreased chromatin accessibility and transcription of lineage-specification genes, including Ebf1 and Pax5, leading to increased self-renewal and reduced HSPC commitment to the B cell lineage. Our data illustrate a role for Stag2 in transformation and transcriptional dysregulation distinct from its shared role with Stag1 in chromosomal segregation.

The prevalence of myopia in schoolchildren in Northeastern Iran is considerably lower than that of East Asian populations, but similar to many other populations, including South Africa, Chile and other countries of the Middle East. Whilst comparisons with other studies show that the prevalence of hyperopia and astigmatism in Northeastern Iran is higher than that of some countries, it is lower compared with others.