BackgroundAssays of the abundance of immune cell populations in the tumor microenvironment promise to inform immune oncology research and the choice of immunotherapy for individual patients. We propose to measure the intratumoral abundance of various immune cell populations with gene expression. In contrast to IHC and flow cytometry, gene expression assays yield high information content from a clinically practical workflow. Previous studies of gene expression in purified immune cells have reported hundreds of genes showing enrichment in a single cell type, but the utility of these genes in tumor samples is unknown. We use co-expression patterns in large tumor gene expression datasets to evaluate previously reported candidate cell type marker genes lists, eliminate numerous false positives and identify a subset of high confidence marker genes.MethodsUsing a novel statistical tool, we use co-expression patterns in 9986 samples from The Cancer Genome Atlas (TCGA) to evaluate previously reported cell type marker genes. We compare immune cell scores derived from these genes to measurements from flow cytometry and immunohistochemistry. We characterize the reproducibility of our cell scores in replicate runs of RNA extracted from FFPE tumor tissue.ResultsWe identify a list of 60 marker genes whose expression levels measure 14 immune cell populations. Cell type scores calculated from these genes are concordant with flow cytometry and IHC readings, show high reproducibility in replicate RNA samples from FFPE tissue and enable detailed analyses of the anti-tumor immune response in TCGA. In an immunotherapy dataset, they separate responders and non-responders early on therapy and provide an intricate picture of the effects of checkpoint inhibition. Most genes previously reported to be enriched in a single cell type have co-expression patterns inconsistent with cell type specificity.ConclusionsDue to their concise gene set, computational simplicity and utility in tumor samples, these cell type gene signatures may be useful in future discovery research and clinical trials to understand how tumors and therapeutic intervention shape the immune response.Electronic supplementary materialThe online version of this article (doi:10.1186/s40425-017-0215-8) contains supplementary material, which is available to authorized users.

Defects in cardiac valve morphogenesis and septation of the heart chambers constitute some of the most common human congenital abnormalities. Some of these defects originate from errors in atrioventricular (AV) endocardial cushion development. Although this process is being extensively studied in mouse and chick, the zebrafish system presents several advantages over these models, including the ability to carry out forward genetic screens and study vertebrate gene function at the single cell level. In this paper, we analyze the cellular and subcellular architecture of the zebrafish heart during stages of AV cushion and valve development and gain an unprecedented level of resolution into this process. We find that endocardial cells in the AV canal differentiate morphologically before the onset of epithelial to mesenchymal transformation, thereby defining a previously unappreciated step during AV valve formation. We use a combination of novel transgenic lines and fluorescent immunohistochemistry to analyze further the role of various genetic (Notch and Calcineurin signaling) and epigenetic (heart function)pathways in this process. In addition, from a large-scale forward genetic screen we identified 55 mutants, defining 48 different genes, that exhibit defects in discrete stages of AV cushion development. This collection of mutants provides a unique set of tools to further our understanding of the genetic basis of cell behavior and differentiation during AV valve development.

While many factors that modulate the morphogenesis and patterning of the embryonic heart have been identified, relatively little is known about the molecular events that regulate the differentiation of progenitor cells fated to form the myocardium. Here, we show that zebrafish grinch (grn) mutants form a reduced number of myocardial progenitor cells, which results in a profound deficit in cardiomyocyte numbers in the most severe cases. We show that grn encodes the G protein-coupled receptor (GPCR) Agtrl1b, a known regulator of adult cardiovascular physiology. Ectopic expression of Apelin, an Agtrl1b ligand, results in the complete absence of cardiomyocytes. Data from transplantation and transgenic approaches indicate that Agtrl1 signaling plays a cell-autonomous role in myocardial specification, with activity being required coincident with the onset of gastrulation movements. These results support a model in which agtrl1b regulates the migration of cells fated to form myocardial progenitors.

In zebrafish embryos, the nascent embryonic shield first appears as a thickening in the germ ring of the mid-epiboly blastoderm. This site defines the dorsal side of the developing embryo. In this paper, we report that the site of embryonic axis formation is marked earlier at the late-blastula stage by the appearance of a cluster of cells with unique endocytic activities. This cluster of cells is composed of enveloping layer epithelial cells and one to two layers of underlying deep cells. Unlike other marginal blastomeres, cells in this cluster do not participate in involution as the blastoderm undergoes epiboly. These noninvoluting endocytic marginal (NEM) cells can be selectively labeled by applying membrane impermeant fluorescent probes to pre-epiboly and mid-epiboly embryos. During embryonic shield formation, deep cells in the NEM cell cluster rearrange and are displaced forward to the leading edge of the blastoderm. As deep NEM cells move into this location, they become a group of cells known as "forerunner cells." Between 60%- and 80%-epiboly, the forerunner cells coalesce into a coherent cell cluster that forms a wedge-shaped cap at the leading edge of the blastoderm. During embryonic axis formation, deep cells migrate and converge toward the embryonic midline, which is defined by the center of the forerunner cell cluster. At approximately 90% epiboly, the forerunner cell cluster becomes overlapped by the constricting germ ring. At tailbud stage, forerunner cells form the dorsal roof of Kupffer's vesicle, which is located ventral to the nascent chordoneural hinge. On the basis of previous grafting studies and known dorsal gene expression patterns, we discuss possible roles that the NEM/forerunner cell cluster may play in teleost axis formation.

Formation of a functional vasculature during mammalian development is essential for embryonic survival. In addition, imbalance in blood vessel growth contributes to the pathogenesis of numerous disorders. Most of our understanding of vascular development and blood vessel growth comes from investigating the Vegf signaling pathway as well as the recent observation that molecules involved in axon guidance also regulate vascular patterning. In order to take an unbiased, yet focused, approach to identify novel genes regulating vascular development, we performed a three-step ENU mutagenesis screen in zebrafish. We first screened live embryos visually, evaluating blood flow in the main trunk vessels, which form by vasculogenesis, and the intersomitic vessels, which form by angiogenesis. Embryos that displayed reduced or absent circulation were fixed and stained for endogenous alkaline phosphatase activity to reveal blood vessel morphology. All putative mutants were then crossed into the Tg(flk1:EGFP)(s843) transgenic background to facilitate detailed examination of endothelial cells in live and fixed embryos. We screened 4015 genomes and identified 30 mutations affecting various aspects of vascular development. Specifically, we identified 3 genes (or loci) that regulate the specification and/or differentiation of endothelial cells, 8 genes that regulate vascular tube and lumen formation, 8 genes that regulate vascular patterning, and 11 genes that regulate vascular remodeling, integrity and maintenance. Only 4 of these genes had previously been associated with vascular development in zebrafish illustrating the value of this focused screen. The analysis of the newly defined loci should lead to a greater understanding of vascular development and possibly provide new drug targets to treat the numerous pathologies associated with dysregulated blood vessel growth.