MARVEL: an integrated alternative splicing analysis platform for single-cell RNA sequencing data.

Alternative splicing is an important source of heterogeneity underlying gene expression between individual cells but remains an understudied area due to the paucity of computational tools to analyze splicing dynamics at single-cell resolution. Here, we present MARVEL, a comprehensive R package for single-cell splicing analysis applicable to RNA sequencing generated from the plate- and droplet-based methods. We performed extensive benchmarking of MARVEL against available tools and demonstrated its utility by analyzing multiple publicly available datasets in diverse cell types, including in disease. MARVEL enables systematic and integrated splicing and gene expression analysis of single cells to characterize the splicing landscape and reveal biological insights.

KIR genotype and haplotype frequencies in the multi-ethnic population of Malaysia.

Killer cell immunoglobulin-like receptors (KIR) genotype and haplotype frequencies have been reported to vary distinctly between populations, which in turn contributes to variation in the alloreactivity of natural killer (NK) cells. Utilizing the diverse KIR genes to identify suitable transplant donors would prove challenging in multi-ethnic countries, even more in resource-limited countries where KIR genotyping has not been established. In this study, we determined the KIR genotypes from 124 unrelated Malaysians consisting of the Malays, Chinese, Indians, and aboriginal people through polymerase chain reaction sequence-specific primer (PCR-SSP) genotyping and employing an expectation-maximization (EM) algorithm to assign haplotypes based on pre-established reference haplotypes. A total of 27 distinct KIR haplotypes were discerned with higher frequencies of haplotype A (55.2%) than haplotype B (44.8%). The most frequent haplotypes were cA01:tA01 (55.2%), cB01:tB01 (18.1%), and cB02:tA01 (13.3%), while the least frequent haplotypes were cB03:tB01 (1.2%), cB04:tB03 (0.4%), and cB03:tA01 (0.4%). Several haplotypes were identified to be unique to a specific ethnic group. The genotype with the highest frequency was genotype AB (71.8%), followed by AA (19.4%), and BB (8.9%). The Indians exhibited the lowest genotype AA but the highest genotype BB, whereas genotype BB was absent in the aboriginal people. Despite the limitations, the genotype and haplotypes in the Malaysian population were successfully highlighted. The identification of ethnic-specific KIR genotypes and haplotypes provides the first step to utilizing KIR in identifying suitable transplant donors to further improve the transplant outcome in the Malaysian population.

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.

VALERIE: Visual-based inspection of alternative splicing events at single-cell resolution.

We present VALERIE (Visualising alternative splicing events from single-cell ribonucleic acid-sequencing experiments), an R package for visualising alternative splicing events at single-cell resolution. To explore any given specified genomic region, corresponding to an alternative splicing event, VALERIE generates an ensemble of informative plots to visualise cell-to-cell heterogeneity of alternative splicing profiles across single cells and performs statistical tests to compare percent spliced-in (PSI) values across the user-defined groups of cells. Among the features available, VALERIE displays PSI values, in lieu of read coverage, which is more suitable for representing alternative splicing profiles for a large number of samples typically generated by single-cell RNA-sequencing experiments. VALERIE is available on the Comprehensive R Archive Network (CRAN): https://cran.r-project.org/web/packages/VALERIE/index.html.

Technological advances and computational approaches for alternative splicing analysis in single cells.

Alternative splicing of RNAs generates isoform diversity, resulting in different proteins that are necessary for maintaining cellular function and identity. The discovery of alternative splicing has been revolutionized by next-generation transcriptomic sequencing mainly using bulk RNA-sequencing, which has unravelled RNA splicing and mis-splicing of normal cells under steady-state and stress conditions. Single-cell RNA-sequencing studies have focused on gene-level expression analysis and revealed gene expression signatures distinguishable between different cellular types. Single-cell alternative splicing is an emerging area of research with the promise to reveal transcriptomic dynamics invisible to bulk- and gene-level analysis. In this review, we will discuss the technological advances for single-cell alternative splicing analysis, computational strategies for isoform detection and quantitation in single cells, and current applications of single-cell alternative splicing analysis and its potential future contributions to personalized medicine.

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.

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.

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.

Repurposing Pentoxifylline for the Treatment of Fibrosis: An Overview.

Fibrosis is a potentially debilitating disease with high morbidity rates. It is estimated that half of all deaths that occur in the USA are attributed to fibrotic disorders. Fibrotic disorders are characterized primarily by disruption in the extracellular matrix deposition and breakdown equilibrium, leading to the accumulation of excessive amounts of extracellular matrix. Given the potentially high prevalence of fibrosis and the paucity of agents currently available for the treatment of this disease, there is an urgent need for the identification of drugs that can be utilized to treat the disease. Pentoxifylline is a methylxanthine derivative that is currently approved for the treatment of vascular diseases, in particular, claudication. Pentoxifylline has three main properties: improving the rheological properties of blood, anti-inflammatory, and antioxidative. Recently, the effectiveness of pentoxifylline in the treatment of fibrosis via attenuating and reversing fibrotic lesions has been demonstrated in several clinical trials and animal studies. As a result of the limited availability of antifibrotic agents in the long-term treatment of fibrosis that can attenuate and even reverse fibrotic lesions effectively, it would be of particular importance to consider the potential clinical utility of pentoxifylline in the treatment of fibrosis. Thus, this paper discusses the evolving roles of pentoxifylline in the treatment of different types of fibrosis.

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.

Evaluation of germline BRCA1 and BRCA2 mutations in a multi-ethnic Asian cohort of ovarian cancer patients.

OBJECTIVE: Despite the discovery of breast and ovarian cancer predisposition genes BRCA1 and BRCA2 more than two decades ago, almost all the available data relate to women of European ancestry, with only a handful of studies in Asian populations. In this study, we determined the frequency of germline alterations in BRCA1 and BRCA2 in ovarian cancer patients from a multi-ethnic cross-sectional cohort of Asian ovarian cancer patients from Malaysia. METHODS: From October 2008 to February 2015, we established a hospital-based cohort of ovarian cancer patients and the germline status of all 218 women with invasive epithelial ovarian cancer was tested using targeted amplification and sequencing of the intron-exon junctions and exonic sequences of BRCA1, BRCA2, PALB2 and TP53. RESULTS: BRCA1 and BRCA2 mutations were found in 8% (17 cases) and 3% (7 cases) of the ovarian cancer patients, respectively. Mutation carriers were diagnosed at a similar age to non-carriers, but were more likely to be Indian, have serous ovarian cancer, and have more relatives with breast or ovarian cancer. Nonetheless, 42% (10/24) of mutation carriers did not have any family history of breast or ovarian cancer and offering genetic counselling and genetic testing only to women with family history would mean that 35% (6/17) of BRCA1 mutation carriers and 57% (4/7) of BRCA2 mutation carriers would not be offered genetic testing. CONCLUSIONS: Our data suggest that, similar to Caucasians, a significant proportion of Asian ovarian cancer was attributed to germline mutations in BRCA1 and to a lesser extent in BRCA2.