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.

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.

gwasurvivr: an R package for genome-wide survival analysis.

SUMMARY: To address the limited software options for performing survival analyses with millions of SNPs, we developed gwasurvivr, an R/Bioconductor package with a simple interface for conducting genome-wide survival analyses using VCF (outputted from Michigan or Sanger imputation servers), IMPUTE2 or PLINK files. To decrease the number of iterations needed for convergence when optimizing the parameter estimates in the Cox model, we modified the R package survival; covariates in the model are first fit without the SNP, and those parameter estimates are used as initial points. We benchmarked gwasurvivr with other software capable of conducting genome-wide survival analysis (genipe, SurvivalGWAS_SV and GWASTools). gwasurvivr is significantly faster and shows better scalability as sample size, number of SNPs and number of covariates increases. AVAILABILITY AND IMPLEMENTATION: gwasurvivr, including source code, documentation and vignette are available at: http://bioconductor.org/packages/gwasurvivr. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

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.