NR2E3 loss disrupts photoreceptor cell maturation and fate in human organoid models of retinal development

While dysfunction and death of light-detecting photoreceptor cells underlie most inherited retinal dystrophies, knowledge of the species-specific details of human rod and cone photoreceptor cell development remains limited. Here, we generate retinal organoids carrying retinal disease-causing variants in NR2E3, as well as isogenic and unrelated controls. Organoids were sampled using single-cell RNA sequencing across the developmental window encompassing photoreceptor specification, emergence, and maturation. Using scRNAseq data, we reconstruct the rod photoreceptor developmental lineage and identify a branchpoint unique to the disease state. We show that the rod-specific transcription factor NR2E3 is required for the proper expression of genes involved in phototransduction, including rhodopsin, which is absent in divergent rods. NR2E3-null rods additionally misexpress several cone-specific phototransduction genes. Using joint multimodal single-cell sequencing, we further identify putative regulatory sites where rod-specific factors act to steer photoreceptor cell development. Finally, we show that rod-committed photoreceptor cells form and persist throughout life in a patient with NR2E3-associated disease. Importantly, these findings are strikingly different than those observed in Nr2e3 rodent models. Together, these data provide a roadmap of human photoreceptor development and leverage patient iPSCs to define the specific roles of rod transcription factors in photoreceptor cell emergence and maturation in health and disease.

Transcriptomic and Chromatin Accessibility Analysis of the Human Macular and Peripheral Retinal Pigment Epithelium at the Single-Cell Level

Some human retinal diseases are characterized by pathology that is restricted to specific cell types and to specific regions of the eye. Several disease entities either selectively affect or spare the macula, the retina region at the center of the posterior pole. Photoreceptor cells in the macula are involved in high-acuity vision and require metabolic support from non-neuronal cell types. Some macular diseases involve the retinal pigment epithelium (RPE), an epithelial cell layer with several metabolic-support functions essential for the overlying photoreceptors. In the current study, the ways in which RPE confers region-specific disease susceptibility were determined by examining heterogeneity within RPE tissue from human donors. RPE nuclei from the macular and peripheral retina were profiled using joint single-nucleus RNA and ATAC sequencing. The expression of several genes differed between macular and peripheral RPE. Region-specific ATAC peaks were found, suggesting regulatory elements used exclusively by macular or peripheral RPE. Across anatomic regions, subpopulations of RPE were identified that appeared to have differential levels of expression of visual cycle genes. Finally, loci associated with age-related macular degeneration were examined for a better understanding of RPE-specific disease phenotypes. These findings showed variations in the regulation of gene expression in the human RPE by region and subpopulation, and provide a source for a better understanding of the molecular basis of macular disease.

NR4A1 deletion promotes pro-angiogenic polarization of macrophages derived from classical monocytes in a mouse model of neovascular age-related macular degeneration

BACKGROUND

Neovascular age-related macular degeneration causes vision loss from destructive angiogenesis, termed choroidal neovascularization (CNV). Cx3cr1-/- mice display alterations in non-classical monocytes and microglia with increased CNV size, suggesting that non-classical monocytes may inhibit CNV formation. NR4A1 is a transcription factor that is necessary for maturation of non-classical monocytes from classical monocytes. While Nr4a1-/- mice are deficient in non-classical monocytes, results are confounded by macrophage hyper-activation. Nr4a1se2/se2 mice lack a transcriptional activator, resulting in non-classical monocyte loss without macrophage hyper-activation.

MAIN BODY

We subjected Nr4a1-/- and Nr4a1se2/se2 mice to the laser-induced CNV model and performed multi-parameter flow cytometry. We found that both models lack non-classical monocytes, but only Nr4a1-/- mice displayed increased CNV area. Additionally, CD11c+ macrophages were increased in Nr4a1-/- mice. Single-cell transcriptomic analysis uncovered that CD11c+ macrophages were enriched from Nr4a1-/- mice and expressed a pro-angiogenic transcriptomic profile that was disparate from prior reports of macrophage hyper-activation.

CONCLUSIONS

These results suggest that non-classical monocytes are dispensable during CNV, and NR4A1 deficiency results in increased recruitment of pro-angiogenic macrophages.

Gene Expression Within a Human Choroidal Neovascular Membrane Using Spatial Transcriptomics

Purpose

Macular neovascularization is a relatively common and potentially visually devastating complication of age-related macular degeneration. In macular neovascularization, pathologic angiogenesis can originate from either the choroid or the retina, but we have limited understanding of how different cell types become dysregulated in this dynamic process.

Methods

To study how gene expression is altered in focal areas of pathology, we performed spatial RNA sequencing on a human donor eye with macular neovascularization as well as a healthy control donor. We performed differential expression to identify genes enriched within the area of macular neovascularization and used deconvolution algorithms to predict the originating cell type of these dysregulated genes.

Results

Within the area of neovascularization, endothelial cells demonstrated increased expression of genes related to Rho family GTPase signaling and integrin signaling. Likewise, VEGF and TGFB1 were identified as potential upstream regulators that could drive the observed gene expression changes produced by endothelial and retinal pigment epithelium cells in the macular neovascularization donor. These spatial gene expression profiles were compared to previous single-cell gene expression experiments in human age-related macular degeneration as well as a model of laser-induced neovascularization in mice. As a secondary aim, we investigated regional gene expression patterns within the macular neural retina and between the macular and peripheral choroid.

Conclusions

Overall, this study spatially analyzes gene expression across the retina, retinal pigment epithelium, and choroid in health and describes a set of candidate molecules that become dysregulated in macular neovascularization.

Multimodal single-cell analysis of nonrandom heteroplasmy distribution in human retinal mitochondrial disease

Variants within the high copy number mitochondrial genome (mtDNA) can disrupt organelle function and lead to severe multisystem disease. The wide range of manifestations observed in patients with mitochondrial disease results from varying fractions of abnormal mtDNA molecules in different cells and tissues, a phenomenon termed heteroplasmy. However, the landscape of heteroplasmy across cell types within tissues and its influence on phenotype expression in affected patients remains largely unexplored. Here, we identify nonrandom distribution of a pathogenic mtDNA variant across a complex tissue using single-cell RNA-Seq, mitochondrial single-cell ATAC sequencing, and multimodal single-cell sequencing. We profiled the transcriptome, chromatin accessibility state, and heteroplasmy in cells from the eyes of a patient with mitochondrial encephalopathy, lactic acidosis, and stroke-like episodes (MELAS) and from healthy control donors. Utilizing the retina as a model for complex multilineage tissues, we found that the proportion of the pathogenic m.3243A>G allele was neither evenly nor randomly distributed across diverse cell types. All neuroectoderm-derived neural cells exhibited a high percentage of the mutant variant. However, a subset of mesoderm-derived lineage, namely the vasculature of the choroid, was near homoplasmic for the WT allele. Gene expression and chromatin accessibility profiles of cell types with high and low proportions of m.3243A>G implicate mTOR signaling in the cellular response to heteroplasmy. We further found by multimodal single-cell sequencing of retinal pigment epithelial cells that a high proportion of the pathogenic mtDNA variant was associated with transcriptionally and morphologically abnormal cells. Together, these findings show the nonrandom nature of mitochondrial variant partitioning in human mitochondrial disease and underscore its implications for mitochondrial disease pathogenesis and treatment.

GENE EXPRESSION WITHIN A HUMAN CHOROIDAL NEOVASCULAR MEMBRANE USING SPATIAL TRANSCRIPTOMICS

Macular neovascularization is a relatively common and potentially visually devastating complication of age-related macular degeneration. In macular neovascularization, pathologic angiogenesis can originate from either the choroid or the retina, but we have limited understanding of how different cell types become dysregulated in this dynamic process. In this study, we performed spatial RNA sequencing on a human donor eye with macular neovascularization as well as a healthy control donor. We identified genes enriched within the area of macular neovascularization and used deconvolution algorithms to predict the originating cell type of these dysregulated genes. Within the area of neovascularization, endothelial cells were predicted to increase expression of genes related to Rho family GTPase signaling and integrin signaling. Likewise, VEGF and TGFB1 were identified as potential upstream regulators that could drive the observed gene expression changes produced by endothelial and retinal pigment epithelium cells in the macular neovascularization donor. These spatial gene expression profiles were compared to previous single-cell gene expression experiments in human age-related macular degeneration as well as a model of laser-induced neovascularization in mice. As a secondary aim, we also investigated spatial gene expression patterns within the macular neural retina and between the macular and peripheral choroid. We recapitulated previously described regional-specific gene expression patterns across both tissues. Overall, this study spatially analyzes gene expression across the retina, retinal pigment epithelium, and choroid in health and describes a set of candidate molecules that become dysregulated in macular neovascularization.

Propensity of Patient-Derived iPSCs for Retinal Differentiation: Implications for Autologous Cell Replacement

Prior to use, newly generated induced pluripotent stem cells (iPSC) should be thoroughly validated. While excellent validation and release testing assays designed to evaluate potency, genetic integrity, and sterility exist, they do not have the ability to predict cell type-specific differentiation capacity. Selection of iPSC lines that have limited capacity to produce high-quality transplantable cells, places significant strain on valuable clinical manufacturing resources. The purpose of this study was to determine the degree and root cause of variability in retinal differentiation capacity between cGMP-derived patient iPSC lines. In turn, our goal was to develop a release testing assay that could be used to augment the widely used ScoreCard panel. IPSCs were generated from 15 patients (14-76 years old), differentiated into retinal organoids, and scored based on their retinal differentiation capacity. Despite significant differences in retinal differentiation propensity, RNA-sequencing revealed remarkable similarity between patient-derived iPSC lines prior to differentiation. At 7 days of differentiation, significant differences in gene expression could be detected. Ingenuity pathway analysis revealed perturbations in pathways associated with pluripotency and early cell fate commitment. For example, good and poor producers had noticeably different expressions of OCT4 and SOX2 effector genes. QPCR assays targeting genes identified via RNA sequencing were developed and validated in a masked fashion using iPSCs from 8 independent patients. A subset of 14 genes, which include the retinal cell fate markers RAX, LHX2, VSX2, and SIX6 (all elevated in the good producers), were found to be predictive of retinal differentiation propensity.

Automating iPSC generation to enable autologous photoreceptor cell replacement therapy

BACKGROUND

Inherited retinal degeneration is a leading cause of incurable vision loss in the developed world. While autologous iPSC mediated photoreceptor cell replacement is theoretically possible, the lack of commercially available technologies designed to enable high throughput parallel production of patient specific therapeutics has hindered clinical translation.

METHODS

In this study, we describe the use of the Cell X precision robotic cell culture platform to enable parallel production of clinical grade patient specific iPSCs. The Cell X is housed within an ISO Class 5 cGMP compliant closed aseptic isolator (Biospherix XVivo X2), where all procedures from fibroblast culture to iPSC generation, clonal expansion and retinal differentiation were performed.

RESULTS

Patient iPSCs generated using the Cell X platform were determined to be pluripotent via score card analysis and genetically stable via karyotyping. As determined via immunostaining and confocal microscopy, iPSCs generated using the Cell X platform gave rise to retinal organoids that were indistinguishable from organoids derived from manually generated iPSCs. In addition, at 120 days post-differentiation, single-cell RNA sequencing analysis revealed that cells generated using the Cell X platform were comparable to those generated under manual conditions in a separate laboratory.

CONCLUSION

We have successfully developed a robotic iPSC generation platform and standard operating procedures for production of high-quality photoreceptor precursor cells that are compatible with current good manufacturing practices. This system will enable clinical grade production of iPSCs for autologous retinal cell replacement.

Systems genomics in age-related macular degeneration

Genomic studies in age-related macular degeneration (AMD) have identified genetic variants that account for the majority of AMD risk. An important next step is to understand the functional consequences and downstream effects of the identified AMD-associated genetic variants. Instrumental for this next step are 'omics' technologies, which enable high-throughput characterization and quantification of biological molecules, and subsequent integration of genomics with these omics datasets, a field referred to as systems genomics. Single cell sequencing studies of the retina and choroid demonstrated that the majority of candidate AMD genes identified through genomic studies are expressed in non-neuronal cells, such as the retinal pigment epithelium (RPE), glia, myeloid and choroidal cells, highlighting that many different retinal and choroidal cell types contribute to the pathogenesis of AMD. Expression quantitative trait locus (eQTL) studies in retinal tissue have identified putative causal genes by demonstrating a genetic overlap between gene regulation and AMD risk. Linking genetic data to complement measurements in the systemic circulation has aided in understanding the effect of AMD-associated genetic variants in the complement system, and supports that protein QTL (pQTL) studies in plasma or serum samples may aid in understanding the effect of genetic variants and pinpointing causal genes in AMD. A recent epigenomic study fine-mapped AMD causal variants by determing regulatory regions in RPE cells differentiated from induced pluripotent stem cells (iPSC-RPE). Another approach that is being employed to pinpoint causal AMD genes is to produce synthetic DNA assemblons representing risk and protective haplotypes, which are then delivered to cellular or animal model systems. Pinpointing causal genes and understanding disease mechanisms is crucial for the next step towards clinical translation. Clinical trials targeting proteins encoded by the AMD-associated genomic loci C3, CFB, CFI, CFH, and ARMS2/HTRA1 are currently ongoing, and a phase III clinical trial for C3 inhibition recently showed a modest reduction of lesion growth in geographic atrophy. The EYERISK consortium recently developed a genetic test for AMD that allows genotyping of common and rare variants in AMD-associated genes. Polygenic risk scores (PRS) were applied to quantify AMD genetic risk, and may aid in predicting AMD progression. In conclusion, genomic studies represent a turning point in our exploration of AMD. The results of those studies now serve as a driving force for several clinical trials. Expanding to omics and systems genomics will further decipher function and causality from the associations that have been reported, and will enable the development of therapies that will lessen the burden of AMD.

Retinal Pigment Epithelium-Secreted VEGF-A Induces Alpha-2-Macroglobulin Expression in Endothelial Cells

Alpha-2-macroglobulin (A2M) is a protease inhibitor that regulates extracellular matrix (ECM) stability and turnover. Here, we show that A2M is expressed by endothelial cells (ECs) from human eye choroid. We demonstrate that retinal pigment epithelium (RPE)-conditioned medium induces A2M expression specifically in ECs. Experiments using chemical inhibitors, blocking antibodies, and recombinant proteins revealed a key role of VEGF-A in RPE-mediated A2M induction in ECs. Furthermore, incubation of ECs with RPE-conditioned medium reduces matrix metalloproteinase-2 gelatinase activity of culture supernatants, which is partially restored after A2M knockdown in ECs. We propose that dysfunctional RPE or choroidal blood vessels, as observed in retinal diseases such as age-related macular degeneration, may disrupt the crosstalk mechanism we describe here leading to alterations in the homeostasis of choroidal ECM, Bruch's membrane and visual function.

The Essential Role of the Choriocapillaris in Vision: Novel Insights from Imaging and Molecular Biology

The choriocapillaris, a dense capillary network located at the posterior pole of the eye, is essential for supporting normal vision, supplying nutrients, and removing waste products from photoreceptor cells and the retinal pigment epithelium. The anatomical location, heterogeneity, and homeostatic interactions with surrounding cell types make the choroid complex to study both in vivo and in vitro. Recent advances in single-cell RNA sequencing, in vivo imaging, and in vitro cell modeling are vastly improving our knowledge of the choroid and its role in normal health and in age-related macular degeneration (AMD). Histologically, loss of endothelial cells (ECs) of the choriocapillaris occurs early in AMD concomitant with elevated formation of the membrane attack complex of complement. Advanced imaging has allowed us to visualize early choroidal blood flow changes in AMD in living patients, supporting histological findings of loss of choroidal ECs. Single-cell RNA sequencing is being used to characterize choroidal cell types transcriptionally and discover their altered patterns of gene expression in aging and disease. Advances in induced pluripotent stem cell protocols and 3D cultures will allow us to closely mimic the in vivo microenvironment of the choroid in vitro to better understand the mechanism leading to choriocapillaris loss in AMD.

PRC2 loss drives MPNST metastasis and matrix remodeling

The histone methyltransferase PRC2 plays a complex role in cancer. Malignant peripheral nerve sheath tumors (MPNSTs) are aggressive sarcomas with frequent loss-of-function mutations in PRC2 that are associated with poor outcome. Here, we identify a critical role for PRC2 loss in driving MPNST metastasis. PRC2-dependent metastatic phenotypes include increased collagen-dependent invasion, upregulation of matrix remodeling enzymes, and elevated lung metastasis in orthotopic mouse models. Furthermore, clinical sample analysis determines that PRC2 loss correlates with metastatic disease, increased fibrosis, and decreased survival in MPNST patients. These results may have broad implications for PRC2 function across multiple cancers and provide a strong rationale for investigating potential therapies targeting ECM remodeling enzymes and tumor fibrosis to improve outcomes in MPNST patients.

Choroidal endothelial and macrophage gene expression in atrophic and neovascular macular degeneration

The human choroid is a heterogeneous, highly vascular connective tissue that dysfunctions in age-related macular degeneration (AMD). In this study, we performed single-cell RNA sequencing on 21 human choroids, 11 of which were derived from donors with early atrophic or neovascular AMD. Using this large donor cohort, we identified new gene expression signatures and immunohistochemically characterized discrete populations of resident macrophages, monocytes/inflammatory macrophages and dendritic cells. These three immune populations demonstrated unique expression patterns for AMD genetic risk factors, with dendritic cells possessing the highest expression of the neovascular AMD-associated MMP9 gene. Additionally, we performed trajectory analysis to model transcriptomic changes across the choroidal vasculature, and we identified expression signatures for endothelial cells from choroidal arterioles and venules. Finally, we performed differential expression analysis between control, early atrophic AMD, and neovascular AMD samples, and we observed that early atrophic AMD samples had high expression of SPARCL1, a gene that has been shown to increase in response to endothelial damage. Likewise, neovascular endothelial cells harbored gene expression changes consistent with endothelial cell damage and demonstrated increased expression of the sialomucins CD34 and ENCM, which were also observed at the protein level within neovascular membranes. Overall, this study characterizes the molecular features of new populations of choroidal endothelial cells and mononuclear phagocytes in a large cohort of AMD and control human donors.

Local Factor H production by human choroidal endothelial cells mitigates complement deposition: implications for macular degeneration

Activation of the alternative complement pathway is an initiating event in the pathology of Age-related Macular Degeneration (AMD). Unchecked complement activation leads to the formation of a pro-lytic pore, the Membrane Attack Complex (MAC). MAC deposition is observed on the choriocapillaris of AMD patients and likely causes lysis of choroidal endothelial cells (CECs). Complement factor H (FH, encoded by the gene CFH), is an inhibitor of complement. Both loss of function of FH and reduced choroidal levels of FH have been reported in AMD. It is plausible that reduced local FH availability promotes MAC deposition on CECs. FH is produced primarily in the liver; however, cells including the retinal pigment epithelium can produce FH locally. We hypothesized that CECs produce FH locally to protect against MAC deposition. We aimed to investigate the effect of reduced FH levels in the choroid to determine whether increasing local FH could protect CECs from MAC deposition. We demonstrated that siRNA knockdown of FH (CFH) in human immortalized CECs results in increased MAC deposition. We generated AMD iPSC-derived CECs and found that overexpression of FH protects against MAC deposition. These results suggest that local CEC-produced FH protects against MAC deposition, and that increasing local FH protein may be beneficial in limiting MAC deposition in AMD. This article is protected by copyright. All rights reserved.

CD177 modulates the function and homeostasis of tumor-infiltrating regulatory T cells

Regulatory T (Treg) cells are one of the major immunosuppressive cell types in cancer and a potential target for immunotherapy, but targeting tumor-infiltrating (TI) Treg cells has been challenging. Here, using single-cell RNA sequencing of immune cells from renal clear cell carcinoma (ccRCC) patients, we identify two distinct transcriptional fates for TI Treg cells, Fate-1 and Fate-2. The Fate-1 signature is associated with a poorer prognosis in ccRCC and several other solid cancers. CD177, a cell surface protein normally expressed on neutrophil, is specifically expressed on Fate-1 TI Treg cells in several solid cancer types, but not on other TI or peripheral Treg cells. Mechanistically, blocking CD177 reduces the suppressive activity of Treg cells in vitro, while Treg-specific deletion of Cd177 leads to decreased tumor growth and reduced TI Treg frequency in mice. Our results thus uncover a functional CD177+ TI Treg population that may serve as a target for TI Treg-specific immunotherapy.

Patient derived stem cells for discovery and validation of novel pathogenic variants in inherited retinal disease

Our understanding of inherited retinal disease has benefited immensely from molecular genetic analysis over the past several decades. New technologies that allow for increasingly detailed examination of a patient's DNA have expanded the catalog of genes and specific variants that cause retinal disease. In turn, the identification of pathogenic variants has allowed the development of gene therapies and low-cost, clinically focused genetic testing. Despite this progress, a relatively large fraction (at least 20%) of patients with clinical features suggestive of an inherited retinal disease still do not have a molecular diagnosis today. Variants that are not obviously disruptive to the codon sequence of exons can be difficult to distinguish from the background of benign human genetic variations. Some of these variants exert their pathogenic effect not by altering the primary amino acid sequence, but by modulating gene expression, isoform splicing, or other transcript-level mechanisms. While not discoverable by DNA sequencing methods alone, these variants are excellent targets for studies of the retinal transcriptome. In this review, we present an overview of the current state of pathogenic variant discovery in retinal disease and identify some of the remaining barriers. We also explore the utility of new technologies, specifically patient-derived induced pluripotent stem cell (iPSC)-based modeling, in further expanding the catalog of disease-causing variants using transcriptome-focused methods. Finally, we outline bioinformatic analysis techniques that will allow this new method of variant discovery in retinal disease. As the knowledge gleaned from previous technologies is informing targets for therapies today, we believe that integrating new technologies, such as iPSC-based modeling, into the molecular diagnosis pipeline will enable a new wave of variant discovery and expanded treatment of inherited retinal disease.

Microfluidic processing of stem cells for autologous cell replacement

Autologous photoreceptor cell replacement is one of the most promising approaches currently under development for the treatment of inherited retinal degenerative blindness. Unlike endogenous stem cell populations, induced pluripotent stem cells (iPSCs) can be differentiated into both rod and cone photoreceptors in high numbers, making them ideal for this application. That said, in addition to photoreceptor cells, state of the art retinal differentiation protocols give rise to all of the different cell types of the normal retina, the majority of which are not required and may in fact hinder successful photoreceptor cell replacement. As such, following differentiation photoreceptor cell enrichment will likely be required. In addition, to prevent the newly generated photoreceptor cells from suffering the same fate as the patient's original cells, correction of the patient's disease‐causing genetic mutations will be necessary. In this review we discuss literature pertaining to the use of different cell sorting and transfection approaches with a focus on the development and use of novel next generation microfluidic devices. We will discuss how gold standard strategies have been used, the advantages and disadvantages of each, and how novel microfluidic platforms can be incorporated into the clinical manufacturing pipeline to reduce the complexity, cost, and regulatory burden associated with clinical grade production of photoreceptor cells for autologous cell replacement.

Direct and indirect immune effects of CMP-001, a virus-like particle containing a TLR9 agonist

Background

CMP-001, also known as vidutolimod, is a virus-like particle containing a TLR9 agonist that is showing promise in early clinical trials. Our group previously demonstrated that the immunostimulatory effects of CMP-001 are dependent on an anti-Qβ antibody response which results in opsonization of CMP-001 and uptake by plasmacytoid dendritic cells (pDCs) that then produce interferon (IFN)-α. IFN-α then leads to an antitumor T-cell response that is responsible for the in vivo efficacy of CMP-001. Here, we explore mechanisms by which the initial effects of CMP-001 on pDCs activate other cells that can contribute to development of an antitumor T-cell response.

Methods

Uptake of CMP-001 by various peripheral blood mononuclear cell (PBMC) populations and response to anti-Qβ-coated CMP-001 were evaluated by flow cytometry and single-cell RNA sequencing. Purified monocytes were treated with anti-Qβ-coated CMP-001 or recombinant IFN-α to evaluate direct and secondary effects of anti-Qβ-coated CMP-001 on monocytes.

Results

Monocytes had the highest per cell uptake of anti-Qβ-coated CMP-001 with lower levels of uptake by pDCs and other cell types. Treatment of PBMCs with anti-Qβ-coated CMP-001 induced upregulation of IFN-responsive genes including CXCL10, PDL1, and indoleamine-2,3-dioxygenase (IDO) expression by monocytes. Most of the impact of anti-Qβ-coated CMP-001 on monocytes was indirect and mediated by IFN-α, but uptake of anti-Qβ-coated CMP-001 altered the monocytic response to IFN-α and resulted in enhanced expression of PDL1, IDO, and CD80 and suppressed expression of CXCL10. These changes included an enhanced ability to induce autologous CD4 T-cell proliferation.

Conclusions

Anti-Qβ-coated CMP-001 induces IFN-α production by pDCs which has secondary effects on a variety of cells including monocytes. Uptake of anti-Qβ-coated CMP-001 by monocytes alters their response to IFN-α, resulting in enhanced expression of PDL1, IDO and CD80 and suppressed expression of CXCL10. Despite aspects of an immunosuppressive phenotype, these monocytes demonstrated increased ability to augment autologous CD4 T-cell proliferation. These findings shed light on the complexity of the mechanism of action of anti-Qβ-coated CMP-001 and provide insight into pathways that may be targeted to further enhance the efficacy of this novel approach to immunotherapy.

Human photoreceptor cells from different macular subregions have distinct transcriptional profiles

The human neural retina is a light sensitive tissue with remarkable spatial and cellular organization. Compared with the periphery, the central retina contains more densely packed cone photoreceptor cells with unique morphologies and synaptic wiring. Some regions of the central retina exhibit selective degeneration or preservation in response to retinal disease and the basis for this variation is unknown. In this study, we used both bulk and single-cell RNA sequencing to compare gene expression within concentric regions of the central retina. We identified unique gene expression patterns of foveal cone photoreceptor cells, including many foveal-enriched transcription factors. In addition, we found that the genes RORB1, PPFIA1 and KCNAB2 are differentially spliced in the foveal, parafoveal and macular regions. These results provide a highly detailed spatial characterization of the retinal transcriptome and highlight unique molecular features of different retinal regions.

Mapping the immune environment in clear cell renal carcinoma by single-cell genomics

Clear cell renal cell carcinoma (ccRCC) is one of the most immunologically distinct tumor types due to high response rate to immunotherapies, despite low tumor mutational burden. To characterize the tumor immune microenvironment of ccRCC, we applied single-cell-RNA sequencing (SCRS) along with T-cell-receptor (TCR) sequencing to map the transcriptomic heterogeneity of 25,688 individual CD45+ lymphoid and myeloid cells in matched tumor and blood from three patients with ccRCC. We also included 11,367 immune cells from four other individuals derived from the kidney and peripheral blood to facilitate the identification and assessment of ccRCC-specific differences. There is an overall increase in CD8+ T-cell and macrophage populations in tumor-infiltrated immune cells compared to normal renal tissue. We further demonstrate the divergent cell transcriptional states for tumor-infiltrating CD8+ T cells and identify a MKI67 + proliferative subpopulation being a potential culprit for the progression of ccRCC. Using the SCRS gene expression, we found preferential prediction of clinical outcomes and pathological diseases by subcluster assignment. With further characterization and functional validation, our findings may reveal certain subpopulations of immune cells amenable to therapeutic intervention.

Abstract B08: Identification and characterization of Rho family GTPases as drivers of drug resistance in BRAFV600 mutant melanoma

The serine/threonine protein kinase BRAF is mutated in approximately 50% of cutaneous melanomas, leading to hyperactivation of the MAPK/ERK pathway. The most common mutations, BRAFV600, can be targeted by selective kinase inhibitors, such as vemurafenib. Although initial clinical response to BRAF inhibition (BRAFi) is encouraging, 90% of patients develop drug resistance within a few months. Drug resistance can be delayed, but not prevented, by combining BRAFi with an MEK inhibitor (MEKi), such as cobimetinib. While some resistance mechanisms are known, disease progression on drug cannot be explained in all patients. We performed a gain-of-function mutagenesis screen utilizing the Sleeping Beauty transposon system to identify novel drivers of resistance in BRAFV600E mutant melanoma cells sensitive to current therapies. We chose four of the top candidates from our screen and validated the ability to drive resistance to both vemurafenib and vemurafenib-cobimetinib combination treatment in multiple melanoma cell lines. In an effort to determine the broader role of candidate vemurafenib-resistance drivers, we conducted an additional in vivo mutagenesis screen, of which genetic analysis is ongoing. Our initial cell-based screen identified two members of the Dbl family of guanine nucleotide exchange factors (GEFs), VAV1 and MCF2, as candidate drivers of vemurafenib resistance. A375 melanoma cells overexpressing VAV1 or MCF2 maintain significant growth under vemurafenib treatment, while control cells do not. Functional tests of VAV1 and MCF2 identified that the active form of two Rho family members, RAC1 and CDC42, increases following treatment with vemurafenib, suggesting a PAK-mediated pathway of resistance. In addition, all candidates that were tested elevated ERK signaling in the presence of vemurafenib. Many of the extracellular signaling pathways known to drive increased vemurafenib resistance activate Rho signaling. Our results suggest that Dbl family members may play an important role in this process. Understanding how Rho activation occurs and its consequences for drug resistance in melanoma will provide critical insights into the design and validation of future targeted therapies.

Citation Format: Jacob L. Schillo, Charlotte R. Feddersen, Afshin Varzavand, Hayley R. Vaughn, Lexy S. Wadsworth, Andrew P. Voigt, Eliot Y. Zhu, Jesse D. Riordan, Christopher S. Stipp, Adam J. Dupuy. Identification and characterization of Rho family GTPases as drivers of drug resistance in BRAFV600 mutant melanoma [abstract]. In: Proceedings of the AACR Special Conference on Melanoma: From Biology to Target; 2019 Jan 15-18; Houston, TX. Philadelphia (PA): AACR; Cancer Res 2020;80(19 Suppl):Abstract nr B08.

Stepwise differentiation and functional characterization of human induced pluripotent stem cell-derived choroidal endothelial cells

BACKGROUND

Endothelial cells (ECs) are essential regulators of the vasculature, lining arteries, veins, and capillary beds. While all ECs share a number of structural and molecular features, heterogeneity exists depending on their resident tissue. ECs lining the choriocapillaris in the human eye are lost early in the pathogenesis of age-related macular degeneration (AMD), a common and devastating form of vision loss. In order to study the mechanisms leading to choroidal endothelial cell (CEC) loss and to develop reagents for repairing the choroid, a reproducible in vitro model, which closely mimic CECs, is needed. While a number of protocols have been published to direct induced pluripotent stem cells (iPSCs) into ECs, the goal of this study was to develop methods to differentiate iPSCs into ECs resembling those found in the human choriocapillaris specifically.

METHODS

We transduced human iPSCs with a CDH5p-GFP-ZEO lentiviral vector and selected for transduced iPSCs using blasticidin. We generated embryoid bodies (EBs) from expanded iPSC colonies and transitioned from mTESR™1 to EC media. One day post-EB formation, we induced mesoderm fate commitment via addition of BMP-4, activin A, and FGF-2. On day 5, EBs were adhered to Matrigel-coated plates in EC media containing vascular endothelial cell growth factor (VEGF) and connective tissue growth factor (CTGF) to promote CEC differentiation. On day 14, we selected for CECs using either zeocin resistance or anti-CD31 MACS beads. We expanded CECs post-selection and performed immunocytochemical analysis of CD31, carbonic anhydrase IV (CA4), and RGCC; tube formation assays; and transmission electron microscopy to access vascular function.

RESULTS

We report a detailed protocol whereby we direct iPSC differentiation toward mesoderm and utilize CTGF to specify CECs. The CDH5p-GFP-ZEO lentiviral vector facilitated the selection of iPSC-derived ECs that label with antibodies directed against CD31, CA4, and RGCC; form vascular tubes in vitro; and migrate into empty choroidal vessels. CECs selected using either antibiotic selection or CD31 MACS beads showed similar characteristics, thereby making this protocol easily reproducible with or without lentiviral vectors.

CONCLUSION

ECs generated following this protocol exhibit functional and biochemical characteristics of CECs. This protocol will be useful for developing in vitro models toward understanding the mechanisms of CEC loss early in AMD.

Spectacle: An interactive resource for ocular single-cell RNA sequencing data analysis

Single-cell RNA sequencing has revolutionized ocular gene expression studies. This technology has enabled researchers to identify expression signatures for rare cell types and characterize how gene expression changes across biological conditions, such as topographic region or disease status. However, sharing single-cell RNA sequencing results remains a major obstacle, particular for individuals without a computational background. To address these limitations, we developed Spectacle, an interactive web-based resource for exploring previously published single-cell RNA sequencing data from ocular studies. Spectacle is powered by a locally developed R package, cellcuratoR, which utilizes the Shiny framework in R to generate interactive visualizations for single-cell expression data. Spectacle contains five pre-processed ocular single-cell RNA sequencing data sets and is accessible via the web at OcularGeneExpression.org/singlecell. With Spectacle, users can interactively identify which cell types express a gene of interest, detect transcriptomic subpopulations within a cell type, and perform highly flexible differential expression analyses. The freely-available Spectacle system reduces the bioinformatic barrier for interacting with rich single-cell RNA sequencing studies from ocular tissues, making it easy to quickly identify cell types that express a gene of interest.

Label-free microfluidic enrichment of photoreceptor cells

Inherited retinal degenerative disorders such as retinitis pigmentosa and Usher syndrome are characterized by progressive death of photoreceptor cells. To restore vision to patients blinded by these diseases, a stem cell-based photoreceptor cell replacement strategy will likely be required. Although retinal stem cell differentiation protocols suitable for generating photoreceptor cells exist, they often yield a rather heterogenous mixture of cell types. To enrich the donor cell population for one or a few cell types, scientists have traditionally relied upon the use of antibody-based selection approaches. However, these strategies are quite labor intensive and require animal derived reagents and equipment that are not well suited to current good manufacturing practices (cGMP). The purpose of this study was to develop and evaluate a microfluidic cell sorting device capable of exploiting the physical and mechanical differences between retinal cell types to enrich specific donor cell populations such as Retinal Pigment Epithelial (RPE) cells and photoreceptor cells. Using this device, we were able to separate a mixture of RPE and iPSC-derived photoreceptor precursor cell lines into two substantially enriched fractions. The enrichment factor of the RPE fraction was 2 and that of the photoreceptor precursor cell fraction was 2.7. Similarly, when human retina, obtained from 3 independent donors, was dissociated and passed through the sorting device, the heterogeneous mixture could be reliably sorted into RPE and photoreceptor cell rich fractions. In summary, microfluidic cell sorting is a promising approach for antibody free enrichment of retinal cell populations.

Bulk and single-cell gene expression analyses reveal aging human choriocapillaris has pro-inflammatory phenotype

The human choroidal vasculature is subject to age-related structural and gene expression changes implicated in age-related macular degeneration (AMD). In this study, we performed both bulk and single-cell RNA sequencing on infant (n = 4 for bulk experiments, n = 2 for single-cell experiments) and adult (n = 13 for bulk experiments, n = 6 for single-cell experiments) human donors to characterize how choroidal gene expression changes with age. Differential expression analysis revealed that aged choroidal samples were enriched in genes encoding pro-inflammatory transcription factors and leukocyte transendothelial cell migration adhesion proteins. Such genes were observed to be differentially expressed specifically within choroidal endothelial cells at the single-cell level. Immunohistochemistry experiments support transcriptional findings that CD34 is elevated in infant choriocapillaris endothelial cells while ICAM-1 is enriched in adults. These results suggest several potential drivers of the pro-inflammatory vascular phenotype observed with advancing age.

Single-Cell RNA Sequencing in Human Retinal Degeneration Reveals Distinct Glial Cell Populations

Degenerative diseases affecting retinal photoreceptor cells have numerous etiologies and clinical presentations. We clinically and molecularly studied the retina of a 70-year-old patient with retinal degeneration attributed to autoimmune retinopathy. The patient was followed for 19 years for progressive peripheral visual field loss and pigmentary changes. Single-cell RNA sequencing was performed on foveal and peripheral retina from this patient and four control patients, and cell-specific gene expression differences were identified between healthy and degenerating retina. Distinct populations of glial cells, including astrocytes and Müller cells, were identified in the tissue from the retinal degeneration patient. The glial cell populations demonstrated an expression profile consistent with reactive gliosis. This report provides evidence that glial cells have a distinct transcriptome in the setting of human retinal degeneration and represents a complementary clinical and molecular investigation of a case of progressive retinal disease.

Src-Dependent DBL Family Members Drive Resistance to Vemurafenib in Human Melanoma

The use of selective BRAF inhibitors (BRAFi) has produced remarkable outcomes for patients with advanced cutaneous melanoma harboring a BRAFV600E mutation. Unfortunately, the majority of patients eventually develop drug-resistant disease. We employed a genetic screening approach to identify gain-of-function mechanisms of BRAFi resistance in two independent melanoma cell lines. Our screens identified both known and unappreciated drivers of BRAFi resistance, including multiple members of the DBL family. Mechanistic studies identified a DBL/RAC1/PAK signaling axis capable of driving resistance to both current and next-generation BRAFis. However, we show that the SRC inhibitor, saracatinib, can block the DBL-driven resistance. Our work highlights the utility of our straightforward genetic screening method in identifying new drug combinations to combat acquired BRAFi resistance. SIGNIFICANCE: A simple, rapid, and flexible genetic screening approach identifies genes that drive resistance to MAPK inhibitors when overexpressed in human melanoma cells.

Deciphering the Significance of CD56 Expression in Pediatric Acute Myeloid Leukemia: A Report from the Children's Oncology Group

BACKGROUND

In patients with acute myeloid leukemia (AML), CD56 expression has been associated with adverse clinical outcome. We reported on a phenotype associated with very poor prognosis (RAM) in children enrolled in the Children's Oncology Group trial AAML0531 (Brodersen et al. Leukemia 30 (2016) 2077-2080). RAM is also characterized in part by high-intensity expression of the CD56 antigen. Herein, we investigate underlying biological and clinical differences among CD56-positive AMLs for patients in AAML0531.

METHODS

For 769 newly diagnosed pediatric patients with de novo AML enrolled in AAML0531, bone marrow specimens were submitted for flow cytometric analysis. For each patient, an immunophenotypic expression profile (IEP) was defined by mean fluorescent intensities of assayed surface antigens. Unsupervised hierarchical clustering analysis (HCA) was completed to group patients with similar immunophenotypes. Clusters were then evaluated for CD56 expression. Principal component analysis (PCA) was subsequently applied to determine whether CD56-positive patient groups were nonoverlapping.

RESULTS

HCA of IEPs revealed three unique phenotypic clusters of patients with CD56-positive AML, and PCA showed that these three cohorts are distinct. Cohort 1 (N = 77) showed a prevalence of t(8;21) patients (72%), Cohort 2 (N = 52) a prevalence of 11q23 patients (69%), and Cohort 3 (RAM) (N = 16) a prevalence of patients with co-occurrence of the CBFA2T3-GLIS2 fusion transcript (63%). The 5-year event-free survival (EFS) for Cohorts 1, 2, and 3 were 69, 39, and 19%, respectively.

CONCLUSIONS

When leukemia is considered by its multidimensional immunophenotype and not by the expression of a single antigen, correlations are seen between genotype and there are significant differences in patient outcomes. © 2019 International Clinical Cytometry Society.

Molecular characterization of foveal versus peripheral human retina by single-cell RNA sequencing

The human retina is a complex tissue responsible for detecting photons of light and converting information from these photons into the neurochemical signals interpreted as vision. Such visual signaling not only requires sophisticated interactions between multiple classes of neurons, but also spatially-dependent molecular specialization of individual cell types. In this study, we performed single-cell RNA sequencing on neural retina isolated from both the fovea and peripheral retina in three human donors. We recovered a total of 8,217 cells, with 3,578 cells originating from the fovea and 4,639 cells originating from the periphery. Expression profiles for all major retinal cell types were compiled, and differential expression analysis was performed between cells of foveal versus peripheral origin. Globally, mRNA for the serum iron binding protein transferrin (TF), which has been associated with age-related macular degeneration pathogenesis, was enriched in peripheral samples. Cone photoreceptor cells were of particular interest and formed two predominant clusters based on gene expression. One cone cluster had 96% of cells originating from foveal samples, while the second cone cluster consisted exclusively of peripherally isolated cells. A total of 148 genes were differentially expressed between cones from the fovea versus periphery. Interestingly, peripheral cones were enriched for the gene encoding Beta-Carotene Oxygenase 2 (BCO2). A relative deficiency of this enzyme may account for the accumulation of carotenoids responsible for yellow pigment deposition within the macula. Overall, this data set provides rich expression profiles of the major human retinal cell types and highlights transcriptomic features that distinguish foveal and peripheral cells.

A simplified transposon mutagenesis method to perform phenotypic forward genetic screens in cultured cells

BACKGROUND

The introduction of genome-wide shRNA and CRISPR libraries has facilitated cell-based screens to identify loss-of-function mutations associated with a phenotype of interest. Approaches to perform analogous gain-of-function screens are less common, although some reports have utilized arrayed viral expression libraries or the CRISPR activation system. However, a variety of technical and logistical challenges make these approaches difficult for many labs to execute. In addition, genome-wide shRNA or CRISPR libraries typically contain of hundreds of thousands of individual engineered elements, and the associated complexity creates issues with replication and reproducibility for these methods.

RESULTS

Here we describe a simple, reproducible approach using the SB transposon system to perform phenotypic cell-based genetic screens. This approach employs only three plasmids to perform unbiased, whole-genome transposon mutagenesis. We also describe a ligation-mediated PCR method that can be used in conjunction with the included software tools to map raw sequence data, identify candidate genes associated with phenotypes of interest, and predict the impact of recurrent transposon insertions on candidate gene function. Finally, we demonstrate the high reproducibility of our approach by having three individuals perform independent replicates of a mutagenesis screen to identify drivers of vemurafenib resistance in cultured melanoma cells.

CONCLUSIONS

Collectively, our work establishes a facile, adaptable method that can be performed by labs of any size to perform robust, genome-wide screens to identify genes that influence phenotypes of interest.

Single-Cell Profiling of Cutaneous T-Cell Lymphoma Reveals Underlying Heterogeneity Associated with Disease Progression

PURPOSE

Cutaneous T-cell lymphomas (CTCL), encompassing a spectrum of T-cell lymphoproliferative disorders involving the skin, have collectively increased in incidence over the last 40 years. Sézary syndrome is an aggressive form of CTCL characterized by significant presence of malignant cells in both the blood and skin. The guarded prognosis for Sézary syndrome reflects a lack of reliably effective therapy, due, in part, to an incomplete understanding of disease pathogenesis.

EXPERIMENTAL DESIGN

Using single-cell sequencing of RNA and the machine-learning reverse graph embedding approach in the Monocle package, we defined a model featuring distinct transcriptomic states within Sézary syndrome. Gene expression used to differentiate the unique transcriptional states were further used to develop a boosted tree classification for early versus late CTCL disease.

RESULTS

Our analysis showed the involvement of FOXP3 + malignant T cells during clonal evolution, transitioning from FOXP3 + T cells to GATA3 + or IKZF2 + (HELIOS) tumor cells. Transcriptomic diversities in a clonal tumor can be used to predict disease stage, and we were able to characterize a gene signature that predicts disease stage with close to 80% accuracy. FOXP3 was found to be the most important factor to predict early disease in CTCL, along with another 19 genes used to predict CTCL stage.

CONCLUSIONS

This work offers insight into the heterogeneity of Sézary syndrome, providing better understanding of the transcriptomic diversities within a clonal tumor. This transcriptional heterogeneity can predict tumor stage and thereby offer guidance for therapy.

TRGAted: A web tool for survival analysis using protein data in the Cancer Genome Atlas

Reverse-phase protein arrays (RPPAs) are a highthroughput approach to protein quantification utilizing antibody-based micro-to-nano scale dot blot. Within the Cancer Genome Atlas (TCGA), RPPAs were used to quantify over 200 proteins in 8,167 tumor and metastatic samples. Protein-level data has particular advantages in assessing putative prognostic or therapeutic targets in tumors. However, many of the available pipelines do not allow for the partitioning of clinical and RPPA information to make meaningful conclusions. We developed a cloud-based application, TRGAted to enable researchers to better examine patient survival based on single or multiple proteins across 31 cancer types in the TCGA. TRGAted contains up-to-date overall survival, disease-specific survival, disease-free interval and progression-free interval information. Furthermore, survival information for primary tumor samples can be stratified based on gender, age, tumor stage, histological type, and subtype, allowing for highly adaptive and intuitive user experience. The code and processed data are open sourced and available on github and contains a tutorial built into the application for assisting users.

TRGAted: A web tool for survival analysis using protein data in the Cancer Genome Atlas

Reverse-phase protein arrays (RPPAs) are a highthroughput approach to protein quantification utilizing antibody-based micro-to-nano scale dot blot. Within the Cancer Genome Atlas (TCGA), RPPAs were used to quantify over 200 proteins in 8,167 tumor and metastatic samples. Protein-level data has particular advantages in assessing putative prognostic or therapeutic targets in tumors. However, many of the available pipelines do not allow for the partitioning of clinical and RPPA information to make meaningful conclusions. We developed a cloud-based application, TRGAted to enable researchers to better examine patient survival based on single or multiple proteins across 31 cancer types in the TCGA. TRGAted contains up-to-date overall survival, disease-specific survival, disease-free interval and progression-free interval information. Furthermore, survival information for primary tumor samples can be stratified based on gender, age, tumor stage, histological type, and subtype, allowing for highly adaptive and intuitive user experience. The code and processed data are open sourced and available on github and contains a tutorial built into the application for assisting users.

Abstract LB-096: A gain-of-function screen identifies drivers of drug resistance in BRAFV600E melanoma cell lines

Melanoma is the deadliest form of skin cancer, and until recently patients had few treatment options upon diagnosis. Fortunately, this lack of treatment options changed with the discovery of vemurafenib, a highly selective kinase inhibitor that specifically targets the BRAFV600E mutant protein, found in ~50% of all melanoma cases. Initially, the kinase inhibitor provided complete or partial response to over 50% of patients and prolonged progression-free survival. However, the majority of patients relapse once tumors acquire resistance to vemurafenib. Some of the genetic resistance mechanisms are known, including BRAF amplification, expression of truncated BRAF, and RAS mutation. However, a significant portion of vemurafenib resistance cases is unexplained. While drug resistance can be delayed by combining vemurafenib with cobimetinib, a MEK inhibitor, mechanisms of resistance to this drug combination are not well characterized.

We performed a gain-of-function screen to identify drivers of drug resistance to vemurafenib or vemurafenib combined with cobimetinib. Briefly, A375 melanoma cells were engineered to express the hyperactive transposase SB100. These cells were then subsequently transfected with the mutagenic T2Onc3 transposon to produce a population of mutagenized A375 cells. Three weeks of vemurafenib treatment selected for cells with the acquired ability to grow in the presence of drug. Candidate driver mutations were then identified by profiling the sites of transposon insertion in drug-resistant cells to find genes that were recurrently over-expressed by transposon insertions. This led to the discovery of five genes that act as drivers of vemurafenib resistance when overexpressed in the human melanoma cell line A375. Candidate genes were then validated in multiple melanoma cell lines by assessing cell growth and viability in varying concentrations of vemurafenib.

An analogous screen was performed in A375 to identify drivers of resistance to the combination of vemurafenib and cobimetinib. In addition, we are assessing the diversity of drug resistance mechanisms by performing mutagenesis in three additional human melanoma cell lines.

The characterization of the novel drug resistance pathways identified here will improve treatment outcomes for patients with BRAFV600E mutant melanoma.

Citation Format: Charlotte R. Feddersen, Jacob L. Schillo, Hayley R. Vaughn, Andrew P. Voigt, Eliot Y. Zhu, Lexy S. Wadsworth, Christopher S. Stipp, Adam J. Dupuy. A gain-of-function screen identifies drivers of drug resistance in BRAFV600E melanoma cell lines [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2018; 2018 Apr 14-18; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2018;78(13 Suppl):Abstract nr LB-096.

Phenotype in combination with genotype improves outcome prediction in acute myeloid leukemia: a report from Children's Oncology Group protocol AAML0531

Diagnostic biomarkers can be used to determine relapse risk in acute myeloid leukemia, and certain genetic aberrancies have prognostic relevance. A diagnostic immunophenotypic expression profile, which quantifies the amounts of distinct gene products, not just their presence or absence, was established in order to improve outcome prediction for patients with acute myeloid leukemia. The immunophenotypic expression profile, which defines each patient’s leukemia as a location in 15-dimensional space, was generated for 769 patients enrolled in the Children’s Oncology Group AAML0531 protocol. Unsupervised hierarchical clustering grouped patients with similar immunophenotypic expression profiles into eleven patient cohorts, demonstrating high associations among phenotype, genotype, morphology, and outcome. Of 95 patients with inv(16), 79% segregated in Cluster A. Of 109 patients with t(8;21), 92% segregated in Clusters A and B. Of 152 patients with 11q23 alterations, 78% segregated in Clusters D, E, F, G, or H. For both inv(16) and 11q23 abnormalities, differential phenotypic expression identified patient groups with different survival characteristics (P<0.05). Clinical outcome analysis revealed that Cluster B (predominantly t(8;21)) was associated with favorable outcome (P<0.001) and Clusters E, G, H, and K were associated with adverse outcomes (P<0.05). Multivariable regression analysis revealed that Clusters E, G, H, and K were independently associated with worse survival (P range <0.001 to 0.008). The Children’s Oncology Group AAML0531 trial: clinicaltrials.gov Identifier: 00372593.

Consistent quantitative gene product expression: #2. Antigen intensities on bone marrow cells are invariant between individuals

Five reference populations in bone marrow specimens were identified by flow cytometry using specific combinations of reagents in order define the variation of gene product expression intensities both within and between individuals. Mature lymphocytes, uncommitted progenitor cells, promyelocytes, mature monocytes and mature neutrophils can be reproducibly identified as distinct clusters of events in heterogeneous, maturing bone marrow specimens. Support Vector Machines were used to identify the reference populations in order to reduce subjective bias in manually defining boundaries of these populations since they were not discretely separated from the remainder of the cells. Reference populations were identified in 50 randomly selected bone marrow aspirates obtained over a period spanning 3 years and 6 months from pediatric patients following chemotherapy for acute myeloid leukemia (AML). The quantitative expression of gene products (cell surface antigens) and light scattering characteristics on these stressed specimens were demonstrated to be tightly regulated both within individuals and between individuals. Within an individual most gene products (CD45, CD34, CD14, CD16, CD64, CD33) demonstrated limited variability with a standard deviation of <0.20 log units while CD13 and CD36 exhibited broader variation >0.25 log units. Surprisingly, with the exception of CD33, the variation of the mean intensities of each antigen between individuals was even less than the variation within an individual. These data confirm that the amounts of gene products expressed on normal developing cells are highly regulated but differ in intensities between different lineages and during the maturational pathway of those lineages. The amounts of gene products expressed at specific stages of development of each lineage are a biologic constant with minimal variation within or between individuals. © 2016 The Authors. Cytometry Part A Published by Wiley Periodicals, Inc. on behalf of ISAC.

Consistent quantitative gene product expression: #3. Invariance with age

The quantitative expression of cell surface antigens and light scattering properties of five cellular reference populations in stressed bone marrow specimens were compared between pediatric and adult patients treated for acute myeloid leukemia (AML). The mean intensity of each antigen as well as the within patient and between patient variability showed striking consistency between the two different age groups. The only difference between the groups of specimens was the proportion of progenitor cells in the adult cohort averaged less than three times the proportion in the pediatric cohort. These data show that the amounts of gene products expressed on bone marrow cells are invariant with age. © 2016 The Authors. Cytometry Part A Published by Wiley Periodicals, Inc. on behalf of ISAC.

Consistent quantitative gene product expression: #1. Automated identification of regenerating bone marrow cell populations using support vector machines

Identification and quantification of maturing hematopoietic cell populations in flow cytometry data sets is a complex and sometimes irreproducible step in data analysis. Supervised machine learning algorithms present promise to automatically classify cells into populations, reducing subjective bias in data analysis. We describe the use of support vector machines (SVMs), a supervised algorithm, to reproducibly identify two distinctly different populations of normal hematopoietic cells, mature lymphocytes and uncommitted progenitor cells, in the challenging setting of pediatric bone marrow specimens obtained 1 month after chemotherapy. Four‐color flow cytometry data were collected on a FACS Calibur for 77 randomly selected postchemotherapy pediatric patients enrolled on the Children's Oncology Group clinical trial AAML1031. These patients demonstrated no evidence of detectable residual disease and were divided into training (n = 27) and testing (n = 50) cohorts. SVMs were trained to identify mature lymphocytes and uncommitted progenitor cells in the training cohort before independent evaluation of prediction efficiency in the testing cohort. Both SVMs demonstrated high predictive performance (lymphocyte SVM: sensitivity >0.99, specificity >0.99; uncommitted progenitor cell SVM: sensitivity = 0.94, specificity >0.99) and closely mirrored manual cell classifications by two expert‐analysts. SVMs present an efficient, automated methodology for identifying normal cell populations even in stressed bone marrows, replicating the performance of an expert while reducing the intrinsic bias of gating procedures between multiple analysts. © 2016 The Authors. Cytometry Part A published by Wiley Periodicals, Inc. on behalf of ISAC.