Human Bone Marrow Organoids for Disease Modeling, Discovery, and Validation of Therapeutic Targets in Hematologic Malignancies.

A lack of models that recapitulate the complexity of human bone marrow has hampered mechanistic studies of normal and malignant hematopoiesis and the validation of novel therapies. Here, we describe a step-wise, directed-differentiation protocol in which organoids are generated from induced pluripotent stem cells committed to mesenchymal, endothelial, and hematopoietic lineages. These 3D structures capture key features of human bone marrow-stroma, lumen-forming sinusoids, and myeloid cells including proplatelet-forming megakaryocytes. The organoids supported the engraftment and survival of cells from patients with blood malignancies, including cancer types notoriously difficult to maintain ex vivo. Fibrosis of the organoid occurred following TGFß stimulation and engraftment with myelofibrosis but not healthy donor-derived cells, validating this platform as a powerful tool for studies of malignant cells and their interactions within a human bone marrow-like milieu. This enabling technology is likely to accelerate the discovery and prioritization of novel targets for bone marrow disorders and blood cancers. SIGNIFICANCE: We present a human bone marrow organoid that supports the growth of primary cells from patients with myeloid and lymphoid blood cancers. This model allows for mechanistic studies of blood cancers in the context of their microenvironment and provides a much-needed ex vivo tool for the prioritization of new therapeutics. See related commentary by Derecka and Crispino, p. 263. This article is highlighted in the In This Issue feature, p. 247.

Processing single-cell RNA-seq datasets using SingCellaR.

Single-cell RNA sequencing has led to unprecedented levels of data complexity. Although several computational platforms are available, performing data analyses for multiple datasets remains a significant challenge. Here, we provide a comprehensive analytical protocol to interrogate multiple datasets on SingCellaR, an analysis package in R. This tool can be applied to general single-cell transcriptome analyses. We demonstrate steps for data analyses and visualization using bespoke pipelines, in conjunction with existing analysis tools to study human hematopoietic stem and progenitor cells. For complete details on the use and execution of this protocol, please refer to Roy et al. (2021).

Transitions in lineage specification and gene regulatory networks in hematopoietic stem/progenitor cells over human development.

Human hematopoiesis is a dynamic process that starts in utero 18-21 days post-conception. Understanding the site- and stage-specific variation in hematopoiesis is important if we are to understand the origin of hematological disorders, many of which occur at specific points in the human lifespan. To unravel how the hematopoietic stem/progenitor cell (HSPC) compartment changes during human ontogeny and the underlying gene regulatory mechanisms, we compare 57,489 HSPCs from 5 different tissues spanning 4 developmental stages through the human lifetime. Single-cell transcriptomic analysis identifies significant site- and developmental stage-specific transitions in cellular architecture and gene regulatory networks. Hematopoietic stem cells show progression from cycling to quiescence and increased inflammatory signaling during ontogeny. We demonstrate the utility of this dataset for understanding aberrant hematopoiesis through comparison to two cancers that present at distinct time points in postnatal life-juvenile myelomonocytic leukemia, a childhood cancer, and myelofibrosis, which classically presents in older adults.

Heterogeneous disease-propagating stem cells in juvenile myelomonocytic leukemia.

Juvenile myelomonocytic leukemia (JMML) is a poor-prognosis childhood leukemia usually caused by RAS-pathway mutations. The cellular hierarchy in JMML is poorly characterized, including the identity of leukemia stem cells (LSCs). FACS and single-cell RNA sequencing reveal marked heterogeneity of JMML hematopoietic stem/progenitor cells (HSPCs), including an aberrant Lin-CD34+CD38-CD90+CD45RA+ population. Single-cell HSPC index-sorting and clonogenic assays show that (1) all somatic mutations can be backtracked to the phenotypic HSC compartment, with RAS-pathway mutations as a "first hit," (2) mutations are acquired with both linear and branching patterns of clonal evolution, and (3) mutant HSPCs are present after allogeneic HSC transplant before molecular/clinical evidence of relapse. Stem cell assays reveal interpatient heterogeneity of JMML LSCs, which are present in, but not confined to, the phenotypic HSC compartment. RNA sequencing of JMML LSC reveals up-regulation of stem cell and fetal genes (HLF, MEIS1, CNN3, VNN2, and HMGA2) and candidate therapeutic targets/biomarkers (MTOR, SLC2A1, and CD96), paving the way for LSC-directed disease monitoring and therapy in this disease.

Cell-intrinsic depletion of Aml1-ETO-expressing pre-leukemic hematopoietic stem cells by K-Ras activating mutation.

Somatic mutations in acute myeloid leukemia are acquired sequentially and hierarchically. First, pre-leukemic mutations, such as t(8;21) that encodes AML1-ETO, are acquired within the hematopoietic stem cell (HSC) compartment, while signaling pathway mutations, including KRAS activating mutations, are late events acquired during transformation of leukemic progenitor cells and are rarely detectable in HSC. This raises the possibility that signaling pathway mutations are detrimental to clonal expansion of pre-leukemic HSC. To address this hypothesis, we used conditional genetics to introduce Aml1-ETO and K-RasG12D into murine HSC, either individually or in combination. In the absence of activated Ras, Aml1-ETO-expressing HSC conferred a competitive advantage. However, activated K-Ras had a marked detrimental effect on Aml1-ETO-expressing HSC, leading to loss of both phenotypic and functional HSC. Cell cycle analysis revealed a loss of quiescence in HSC co-expressing Aml1-ETO and K-RasG12D, accompanied by an enrichment in E2F and Myc target gene expression and depletion of HSC self-renewal-associated gene expression. These findings provide a mechanistic basis for the observed absence of KRAS signaling mutations in the pre-malignant HSC compartment.

Association of BRCA1- and BRCA2-deficiency with mutation burden, expression of PD-L1/PD-1, immune infiltrates, and T cell-inflamed signature in breast cancer.

Immune checkpoint inhibitors have demonstrated effective anti-tumour response in cancer types with high mutation burden (e.g. melanoma) and in subset of cancers with features of genomic instability (e.g. mismatch-repair deficiency). One possible explanation for this effect is the increased expression of immune checkpoint molecules and pre-existing adaptive immune response in these cancers. Given that BRCA1 and BRCA2 are integral in maintaining genomic integrity, we hypothesise that the inactivation of these genes may give rise to breast cancers with such immunogenic phenotype. Therefore, using two large series of publicly available breast cancer datasets, namely that from The Cancer Genome Atlas and Wellcome Trust Institute, we sought to investigate the association between BRCA1- and BRCA2-deficiency with features of genomic instability, expression of PD-L1 and PD-1, landscape of inferred tumour-infiltrating immune cells, and T-cell inflamed signature in breast cancers. Here, we report that BRCA1 and BRCA2-deficient breast cancers were associated with features of genomic instability including increased mutation burden. Interestingly, BRCA1-, but not BRCA2-, deficient breast cancers were associated with increased expression of PD-L1 and PD-1, higher abundance of tumour-infiltrating immune cells, and enrichment of T cell-inflamed signature. The differences in immunophenotype between BRCA1- and BRCA2-deficient breast cancers can be attributed, in part, to PTEN gene mutation. Therefore, features of genomic instability such as that mediated by BRCA1- and BRCA2- deficiency in breast cancer were necessary, but not always sufficient, for yielding T cell-inflamed tumour microenvironment, and by extension, predicting clinical benefit from immunotherapy.

Prevalence of BRCA1 and BRCA2 pathogenic variants in a large, unselected breast cancer cohort.

Breast cancer patients with BRCA1/2-driven tumors may benefit from targeted therapy. It is not clear whether current BRCA screening guidelines are effective at identifying these patients. The purpose of our study was to evaluate the prevalence of inherited BRCA1/2 pathogenic variants in a large, clinically representative breast cancer cohort and to estimate the proportion of BRCA1/2 carriers not detected by selectively screening individuals with the highest probability of being carriers according to current clinical guidelines. The study included 5,122 unselected Swedish breast cancer patients diagnosed from 2001 to 2008. Target sequence enrichment (48.48 Fluidigm Access Arrays) and sequencing were performed (Illumina Hi-Seq 2,500 instrument, v4 chemistry). Differences in patient and tumor characteristics of BRCA1/2 carriers who were already identified as part of clinical BRCA1/2 testing routines and additional BRCA1/2 carriers found by sequencing the entire study population were compared using logistic regression models. Ninety-two of 5,099 patients with valid variant calls were identified as BRCA1/2 carriers by screening all study participants (1.8%). Only 416 study participants (8.2%) were screened as part of clinical practice, but this identified 35 out of 92 carriers (38.0%). Clinically identified carriers were younger, less likely postmenopausal and more likely to be associated with familiar ovarian cancer compared to the additional carriers identified by screening all patients. More BRCA2 (34/42, 81.0%) than BRCA1 carriers (23/50, 46%) were missed by clinical screening. In conclusion, BRCA1/2 mutation prevalence in unselected breast cancer patients was 1.8%. Six in ten BRCA carriers were not detected by selective clinical screening of individuals.

Differential Burden of Rare and Common Variants on Tumor Characteristics, Survival, and Mode of Detection in Breast Cancer.

Genetic variants that increase breast cancer risk can be rare or common. This study tests whether the genetic risk stratification of breast cancer by rare and common variants in established loci can discriminate tumors with different biology, patient survival, and mode of detection. Multinomial logistic regression tested associations between genetic risk load [protein-truncating variant (PTV) carriership in 31 breast cancer predisposition genes-or polygenic risk score (PRS) using 162 single-nucleotide polymorphisms], tumor characteristics, and mode of detection (OR). Ten-year breast cancer-specific survival (HR) was estimated using Cox regression models. In this unselected cohort of 5,099 patients with breast cancer diagnosed in Sweden between 2001 and 2008, PTV carriers (n = 597) were younger and associated with more aggressive tumor phenotypes (ER-negative, large size, high grade, high proliferation, luminal B, and basal-like subtype) and worse outcome (HR, 1.65; 1.16-2.36) than noncarriers. After excluding 92 BRCA1/2 carriers, PTV carriership remained associated with high grade and worse survival (HR, 1.76; 1.21-2.56). In 5,007 BRCA1/2 noncarriers, higher PRS was associated with less aggressive tumor characteristics (ER-positive, PR-positive, small size, low grade, low proliferation, and luminal A subtype). Among patients with low mammographic density (<25%), non-BRCA1/2 PTV carriers were more often interval than screen-detected breast cancer (OR, 1.89; 1.12-3.21) than noncarriers. In contrast, higher PRS was associated with lower risk of interval compared with screen-detected cancer (OR, 0.77; 0.64-0.93) in women with low mammographic density. These findings suggest that rare and common breast cancer susceptibility loci are differentially associated with tumor characteristics, survival, and mode of detection.Significance: These findings offer the potential to improve screening practices for breast cancer by providing a deeper understanding of how risk variants affect disease progression and mode of detection. Cancer Res; 78(21); 6329-38. ©2018 AACR.

Inherited mutations in BRCA1 and BRCA2 in an unselected multiethnic cohort of Asian patients with breast cancer and healthy controls from Malaysia.

BACKGROUND: Genetic testing for BRCA1 and BRCA2 is offered typically to selected women based on age of onset and family history of cancer. However, current internationally accepted genetic testing referral guidelines are built mostly on data from cancer genetics clinics in women of European descent. To evaluate the appropriateness of such guidelines in Asians, we have determined the prevalence of germ line variants in an unselected cohort of Asian patients with breast cancer and healthy controls. METHODS: Germ line DNA from a hospital-based study of 2575 unselected patients with breast cancer and 2809 healthy controls were subjected to amplicon-based targeted sequencing of exonic and proximal splice site junction regions of BRCA1 and BRCA2 using the Fluidigm Access Array system, with sequencing conducted on a Illumina HiSeq2500 platform. Variant calling was performed with GATK UnifiedGenotyper and were validated by Sanger sequencing. RESULTS: Fifty-five (2.1%) BRCA1 and 66 (2.6%) BRCA2 deleterious mutations were identified among patients with breast cancer and five (0.18%) BRCA1 and six (0.21%) BRCA2 mutations among controls. One thousand one hundred and eighty-six (46%) patients and 97 (80%) carriers fulfilled the National Comprehensive Cancer Network guidelines for genetic testing. CONCLUSION: Five per cent of unselected Asian patients with breast cancer carry deleterious variants in BRCA1 or BRCA2. While current referral guidelines identified the majority of carriers, one in two patients would be referred for genetic services. Given that such services are largely unavailable in majority of low-resource settings in Asia, our study highlights the need for more efficient guidelines to identify at-risk individuals in Asia.

Germline APOBEC3B deletion is associated with breast cancer risk in an Asian multi-ethnic cohort and with immune cell presentation.

BACKGROUND: APOBEC3B is a cytosine deaminase implicated in immune response to viral infection, cancer predisposition and carcinogenesis. Germline APOBEC3B deletion is more common in East Asian women and confers a modest risk to breast cancer in both East Asian and Caucasian women. Analysis of tumour samples from women of European descent has shown that germline APOBEC3B deletion is associated with an increased propensity to develop somatic mutations and with an enrichment for immune response-related gene sets. However, this has not been examined in Asian tumour samples, where population differences in genetic and dietary factors may have an impact on the immune system. METHODS: In this study, we determined the prevalence of germline APOBEC3B deletion and its association with breast cancer risk in a cross-sectional hospital-based Asian multi-ethnic cohort of 1451 cases and 1442 controls from Malaysia. We compared gene expression profiles of breast cancers arising from APOBEC3B deletion carriers and non-carriers using microarray analyses. Finally, we characterised the overall abundance of tumour-infiltrating immune cells in breast cancers from TCGA and METABRIC using ESTIMATE and relative frequency of 22 immune cell subsets in breast cancers from METABRIC using CIBERSORT. RESULTS: The minor allelic frequency of APOBEC3B deletion was estimated to be 0.35, 0.42 and 0.16 in female populations of Chinese, Malay and Indian descent, respectively, and that germline APOBEC3B deletion was associated with breast cancer risk with odds ratios of 1.23 (95 % CI: [1.05, 1.44]) for one-copy deletion and 1.38 (95 % CI: [1.10, 1.74]) for two-copy deletion compared to women with no deletion. Germline APOBEC3B deletion was not associated with any clinicopathologic features or the expression of any APOBEC family members but was associated with immune response-related gene sets (FDR q values < 0.05). Analysis of breast cancers from METABRIC revealed breast cancers from APOBEC3B deletion carriers to have significantly higher abundance of tumour-infiltrating immune cells (P < 0.001). CONCLUSIONS: Taken together, our data suggests that tumour-infiltrating immune cells may be an important feature of breast cancers arising in women with APOBEC3B germline deletion, and that this may be of particular interest in Asian women where the germline deletion is more common.