Development and validation of a deep learning-based microsatellite instability predictor from prostate cancer whole-slide images.

Microsatellite instability-high (MSI-H) is a tumor-agnostic biomarker for immune checkpoint inhibitor therapy. However, MSI status is not routinely tested in prostate cancer, in part due to low prevalence and assay cost. As such, prediction of MSI status from hematoxylin and eosin (H&E) stained whole-slide images (WSIs) could identify prostate cancer patients most likely to benefit from confirmatory testing to evaluate their eligibility for immunotherapy and need for Lynch syndrome testing. Prostate biopsies and surgical resections from prostate cancer patients referred to our institution were analyzed. MSI status was determined by next-generation sequencing. Patients sequenced before a cutoff date formed an algorithm development set (n = 4015, MSI-H 1.8%) and a paired validation set (n = 173, MSI-H 19.7%) that consisted of two serial sections from each sample, one stained and scanned internally and the other at an external site. Patients sequenced after the cutoff date formed a temporally independent validation set (n = 1350, MSI-H 2.3%). Attention-based multiple instance learning models were trained to predict MSI-H from H&E WSIs. The predictor achieved area under the receiver operating characteristic curve values of 0.78 (95% CI [0.69-0.86]), 0.72 (95% CI [0.63-0.81]), and 0.72 (95% CI [0.62-0.82]) on the internally prepared, externally prepared, and temporal validation sets, respectively, showing effective predictability and generalization to both external staining/scanning processes and temporally independent samples. While MSI-H status is significantly correlated with Gleason score, the model remained predictive within each Gleason score subgroup.

Prediction of MET Overexpression in Lung Adenocarcinoma from Hematoxylin and Eosin Images.

Mesenchymal epithelial transition (MET) protein overexpression is a targetable event in non-small cell lung cancer and is the subject of active drug development. Challenges in identifying patients for these therapies include lack of access to validated testing, such as standardized immunohistochemistry assessment, and consumption of valuable tissue for a single gene/protein assay. Development of prescreening algorithms using routinely available digitized hematoxylin and eosin (H&E)-stained slides to predict MET overexpression could promote testing for those who will benefit most. Recent literature reports a positive correlation between MET protein overexpression and RNA expression. In this work, a large database of matched H&E slides and RNA expression data were leveraged to train a weakly supervised model to predict MET RNA overexpression directly from H&E images. This model was evaluated on an independent holdout test set of 300 overexpressed and 289 normal patients, demonstrating a receiver operating characteristic area under curve of 0.70 (95th percentile interval: 0.66 to 0.74) with stable performance characteristics across different patient clinical variables and robust to synthetic noise on the test set. These results suggest that H&E-based predictive models could be useful to prioritize patients for confirmatory testing of MET protein or MET gene expression status.

A functional cellular framework for sex and estrous cycle-dependent gene expression and behavior.

Sex hormones exert a profound influence on gendered behaviors. How individual sex hormone-responsive neuronal populations regulate diverse sex-typical behaviors is unclear. We performed orthogonal, genetically targeted sequencing of four estrogen receptor 1-expressing (Esr1+) populations and identified 1,415 genes expressed differentially between sexes or estrous states. Unique subsets of these genes were distributed across all 137 transcriptomically defined Esr1+ cell types, including estrous stage-specific ones, that comprise the four populations. We used differentially expressed genes labeling single Esr1+ cell types as entry points to functionally characterize two such cell types, BNSTprTac1/Esr1 and VMHvlCckar/Esr1. We observed that these two cell types, but not the other Esr1+ cell types in these populations, are essential for sex recognition in males and mating in females, respectively. Furthermore, VMHvlCckar/Esr1 cell type projections are distinct from those of other VMHvlEsr1 cell types. Together, projection and functional specialization of dimorphic cell types enables sex hormone-responsive populations to regulate diverse social behaviors.

Novel genetic variants associated with mortality after unrelated donor allogeneic hematopoietic cell transplantation.

BACKGROUND: Identification of non-human leukocyte antigen (HLA) genetic risk factors could improve survival after allogeneic blood or marrow transplant (BMT) through matching at additional loci or individualizing risk prediction. We hypothesized that non-HLA loci contributed significantly to 1-year overall survival (OS), disease related mortality (DRM) or transplant related mortality (TRM) after unrelated donor (URD)BMT. METHODS: We performed a genome-wide association study (GWAS) in 2,887 acute myeloid leukemia (AML), myelodysplastic syndrome (MDS) and acute lymphoblastic leukemia (ALL) patients and their ≥8/8 HLA-matched URDs comprising two independent cohorts treated from 2000-2011. FINDINGS: Using meta-analyses of both cohorts, genome-wide significant associations (p < 5 × 10-8) were identified in: recipient genomes with OS at MBNL1 (rs9990017, HR = 1.4, 95% CI 1.24-1.56, p = 3.3 × 10-8) and donor-recipient genotype mismatch with OS at LINC02774 (rs10927108, HR = 1.34, 95% CI 1.21-1.48, p = 2.0 × 10-8); donor genomes with DRM at PCNX4 (rs79076914, HR = 1.7, 95% CI 1.41-2.05, p = 3.15 × 10-8), LINC01194 (rs79498125, HR = 1.86, 95% CI 1.49-2.31, p = 2.84 × 10-8), ARID5B (rs2167710, HR = 1.5, 95% CI 1.31-1.73, p = 6.9 × 10-9) and CT49 (rs32250, HR = 1.44, 95% CI1.26-1.64, p = 2.6 × 10-8); recipient genomes at PILRB with TRM (rs141591562, HR = 2.33, 95% CI 1.74-3.12, p = 1.26 × 10-8) and donor-recipient genotype mismatch between EPGN and MTHF2DL with TRM (rs75868097, HR = 2.66, 95% CI 1.92-3.58, p = 4.6 × 10-9). Results publicly available at https://fuma.ctglab.nl/browse. INTERPRETATION: These data provide the first evidence that non-HLA common genetic variation at novel loci with biochemical function significantly impacts 1-year URD-BMT survival. Our findings have implications for donor selection, could guide treatment strategies and provide individualized risk prediction after future validation and functional studies. FUNDING: This project was funded by grants from the National Institutes of Health, USA.

Genome-Wide Association Analyses Identify Variants in IRF4 Associated With Acute Myeloid Leukemia and Myelodysplastic Syndrome Susceptibility.

The role of common genetic variation in susceptibility to acute myeloid leukemia (AML), and myelodysplastic syndrome (MDS), a group of rare clonal hematologic disorders characterized by dysplastic hematopoiesis and high mortality, remains unclear. We performed AML and MDS genome-wide association studies (GWAS) in the DISCOVeRY-BMT cohorts (2,309 cases and 2,814 controls). Association analysis based on subsets (ASSET) was used to conduct a summary statistics SNP-based analysis of MDS and AML subtypes. For each AML and MDS case and control we used PrediXcan to estimate the component of gene expression determined by their genetic profile and correlate this imputed gene expression level with risk of developing disease in a transcriptome-wide association study (TWAS). ASSET identified an increased risk for de novo AML and MDS (OR = 1.38, 95% CI, 1.26-1.51, Pmeta = 2.8 × 10-12) in patients carrying the T allele at s12203592 in Interferon Regulatory Factor 4 (IRF4), a transcription factor which regulates myeloid and lymphoid hematopoietic differentiation. Our TWAS analyses showed increased IRF4 gene expression is associated with increased risk of de novo AML and MDS (OR = 3.90, 95% CI, 2.36-6.44, Pmeta = 1.0 × 10-7). The identification of IRF4 by both GWAS and TWAS contributes valuable insight on the role of genetic variation in AML and MDS susceptibility.

Prognostic impact of pre-transplant chromosomal aberrations in peripheral blood of patients undergoing unrelated donor hematopoietic cell transplant for acute myeloid leukemia.

To improve risk stratification and treatment decisions for patients with acute myeloid leukemia (AML) undergoing hematopoietic cell transplantation (HCT). We used SNP-array data from the DISCOVeRY-BMT study to detect chromosomal aberrations in pre-HCT peripheral blood (collected 2-4 weeks before the administration of conditioning regimen) from 1974 AML patients who received HCT between 2000 and 2011. All aberrations detected in ≥ 10 patients were tested for their association with overall survival (OS), separately by remission status, using the Kaplan-Meier estimator. Cox regression models were used for multivariable analyses. Follow-up was through January 2019. We identified 701 unique chromosomal aberrations in 285 patients (7% of 1438 in complete remission (CR) and 36% of 536 not in CR). Copy-neutral loss-of-heterozygosity (CNLOH) in chr17p in CR patients (3-year OS = 20% vs. 50%, with and without chr17p CNLOH, p = 0.0002), and chr13q in patients not in CR (3-year OS = 4% vs. 26%, with and without chr13q CNLOH, p < 0.0001) are risk factors for poor survival. Models adjusted for clinical factors showed approximately three-fold excess risk of post-HCT mortality with chr17p CNLOH in CR patients (hazard ratio, HR = 3.39, 95% confidence interval CI 1.74-6.60, p = 0.0003), or chr13q CNLOH in patients not in CR (HR = 2.68, 95% CI 1.75-4.09, p < 0.0001). The observed mortality was mostly driven by post-HCT relapse (HR = 2.47, 95% CI 1.01-6.02, p = 0.047 for chr17p CNLOH in CR patients, and HR = 2.58, 95% CI 1.63-4.08, p < 0.0001 for chr13q CNLOH in patients not in CR. Pre-transplant CNLOH in chr13q or chr17p predicts risk of poor outcomes after unrelated donor HCT in AML patients. A large prospective study is warranted to validate the results and evaluate novel strategies to improve survival in those patients.

Multiple functional variants in the IL1RL1 region are pretransplant markers for risk of GVHD and infection deaths.

Graft-versus-host disease (GVHD) and infections are the 2 main causes of death without relapse after allogeneic hematopoietic cell transplantation (HCT). Elevated soluble serum simulation-2 (sST2), the product of IL1RL1 in plasma/serum post-HCT, is a validated GVHD biomarker. Hundreds of SNPs at 2q12.1 have been shown to be strongly associated with sST2 concentrations in healthy populations. We therefore hypothesized that the donor genetic variants in IL1RL1 correlate with sST2 protein levels associated with patient survival outcomes after HCT. We used DISCOVeRY-BMT (Determining the Influence of Susceptibility Conveying Variants Related to 1-Year Mortality after Blood and Marrow Transplantation), a genomic study of >3000 donor-recipient pairs, to inform our hypothesis. We first measured pre-HCT plasma/serum sST2 levels in a subset of DISCOVeRY-BMT donors (n = 757) and tested the association of donor sST2 levels with donor single nucleotide polymorphisms (SNPs) in the 2q12.1 region. Donor SNPs associated with sST2 levels were then tested for association with recipient death caused by acute GVHD (aGVHD)-, infection-, and transplant-related mortality in cohorts 1 and 2. Meta-analyses of cohorts 1 and 2 were performed using fixed-effects inverse variance weighting, and P values were corrected for multiple comparisons. Donor risk alleles in rs22441131 (P meta = .00026) and rs2310241 (P meta = .00033) increased the cumulative incidence of aGVHD death up to fourfold and were associated with high sST2 levels. Donor risk alleles at rs4851601 (P meta = 9.7 × 10-7), rs13019803 (P meta = 8.9 × 10-6), and rs13015714 (P meta = 5.3 × 10-4) increased cumulative incidence of infection death to almost sevenfold and were associated with low sST2 levels. These functional variants are biomarkers of infection or aGVHD death and could facilitate donor selection, prophylaxis, and a conditioning regimen to reduce post-HCT mortality.

The miR-96 and RARγ signaling axis governs androgen signaling and prostate cancer progression.

Expression levels of retinoic acid receptor gamma (NR1B3/RARG, encodes RARγ) are commonly reduced in prostate cancer (PCa). Therefore, we sought to establish the cellular and gene regulatory consequences of reduced RARγ expression, and determine RARγ regulatory mechanisms. RARG shRNA approaches in non-malignant (RWPE-1 and HPr1-AR) and malignant (LNCaP) prostate models revealed that reducing RARγ levels, rather than adding exogenous retinoid ligand, had the greatest impact on prostate cell viability and gene expression. ChIP-Seq defined the RARγ cistrome, which was significantly enriched at active enhancers associated with AR binding sites. Reflecting a significant genomic role for RARγ to regulate androgen signaling, RARγ knockdown in HPr1-AR cells significantly regulated the magnitude of the AR transcriptome. RARγ downregulation was explained by increased miR-96 in PCa cell and mouse models, and TCGA PCa cohorts. Biochemical approaches confirmed that miR-96 directly regulated RARγ expression and function. Capture of the miR-96 targetome by biotin-miR-96 identified that RARγ and a number of RARγ interacting co-factors including TACC1 were all targeted by miR-96, and expression of these genes were prominently altered, positively and negatively, in the TCGA-PRAD cohort. Differential gene expression analyses between tumors in the TCGA-PRAD cohort with lower quartile expression levels of RARG and TACC1 and upper quartile miR-96, compared to the reverse, identified a gene network including several RARγ target genes (e.g., SOX15) that significantly associated with worse disease-free survival (hazard ratio 2.23, 95% CI 1.58 to 2.88, p = 0.015). In summary, miR-96 targets a RARγ network to govern AR signaling, PCa progression and disease outcome.

Genome-wide meta-analyses identifies novel taxane-induced peripheral neuropathy-associated loci.

OBJECTIVE: Taxane containing chemotherapy extends survival for breast cancer patients. However, taxane-induced peripheral neuropathy (TIPN) cannot be predicted, prevented or effectively treated. Using genome-wide analyses, we sought to identify common risk variants for TIPN. PATIENTS AND METHODS: Women with high-risk breast cancer enrolled in SWOG 0221 were genotyped using the Illumina 1M chip. Genome-wide analyses were performed in relation to ≥grade 3 Common Terminology Criteria for Adverse Events (CTCAE) neuropathy in European and African Americans. Data were meta-analyzed with GW associations of CTCAE ≥grade 3 versus

Replication and validation of genetic polymorphisms associated with survival after allogeneic blood or marrow transplant.

Multiple candidate gene-association studies of non-HLA single-nucleotide polymorphisms (SNPs) and outcomes after blood or marrow transplant (BMT) have been conducted. We identified 70 publications reporting 45 SNPs in 36 genes significantly associated with disease-related mortality, progression-free survival, transplant-related mortality, and/or overall survival after BMT. Replication and validation of these SNP associations were performed using DISCOVeRY-BMT (Determining the Influence of Susceptibility COnveying Variants Related to one-Year mortality after BMT), a well-powered genome-wide association study consisting of 2 cohorts, totaling 2888 BMT recipients with acute myeloid leukemia, acute lymphoblastic leukemia, or myelodysplastic syndrome, and their HLA-matched unrelated donors, reported to the Center for International Blood and Marrow Transplant Research. Gene-based tests were used to assess the aggregate effect of SNPs on outcome. None of the previously reported significant SNPs replicated at P < .05 in DISCOVeRY-BMT. Validation analyses showed association with one previously reported donor SNP at P < .05 and survival; more associations would be anticipated by chance alone. No gene-based tests were significant at P < .05. Functional annotation with publicly available data shows these candidate SNPs most likely do not have biochemical function; only 13% of candidate SNPs correlate with gene expression or are predicted to impact transcription factor binding. Of these, half do not impact the candidate gene of interest; the other half correlate with expression of multiple genes. These findings emphasize the peril of pursing candidate approaches and the importance of adequately powered tests of unbiased genome-wide associations with BMT clinical outcomes given the ultimate goal of improving patient outcomes.

Genetic association with B-cell acute lymphoblastic leukemia in allogeneic transplant patients differs by age and sex.

The incidence and mortality rates of B-cell acute lymphoblastic leukemia (B-ALL) differ by age and sex. To determine if inherited genetic susceptibility contributes to these differences we performed 2 genome-wide association studies (GWAS) by age, sex, and subtype and subsequent meta-analyses. The GWAS included 446 B-ALL cases, and 3027 healthy unrelated blood and marrow transplant (BMT) donors as controls from the Determining the Influence of Susceptibility Conveying Variants Related to One-Year Mortality after BMT (DISCOVeRY-BMT) study. We identified 1 novel variant, rs189434316, significantly associated with odds of normal cytogenetic B-ALL (odds ratio from meta-analysis [ORmeta] = 3.7; 95% confidence interval [CI], 2.5, 6.2; P value from meta-analysis [Pmeta] = 6.0 × 10-9). The previously reported pediatric B-ALL GWAS variant, rs11980379 (IKZF1), replicated in B-ALL pediatric patients (ORmeta = 2.3; 95% CI, 1.5, 3.7; Pmeta = 1.0 × 10-9), with evidence of heterogeneity (P = .02) between males and females. Sex differences in single-nucleotide polymorphism effect were seen in those >15 years (OR = 1.7; 95% CI, 1.4, 2.2, PMales = 6.38 × 10-6/OR = 1.1; 95% CI, 0.8, 1.5; PFemales = .6) but not ≤15 years (OR = 2.3; 95% CI, 1.4, 3.8; PMales = .0007/OR = 1.9; 95% CI, 1.2, 3.2; PFemales = .007). The latter association replicated in independent pediatric B-ALL cohorts. A previously identified adolescent and young-adult onset ALL-associated variant in GATA3 is associated with B-ALL risk in those >40 years. Our findings provide more evidence of the influence of genetics on B-ALL age of onset and we have shown the first evidence that IKZF1 associations with B-ALL may be sex and age specific.

Scalable microfluidics for single-cell RNA printing and sequencing.

Many important biological questions demand single-cell transcriptomics on a large scale. Hence, new tools are urgently needed for efficient, inexpensive manipulation of RNA from individual cells. We report a simple platform for trapping single-cell lysates in sealed, picoliter microwells capable of printing RNA on glass or capturing RNA on beads. We then develop a scalable technology for genome-wide, single-cell RNA-Seq. Our device generates pooled libraries from hundreds of individual cells with consumable costs of $0.10-$0.20 per cell and includes five lanes for simultaneous experiments. We anticipate that this system will serve as a general platform for single-cell imaging and sequencing.