Torch-eCpG: a fast and scalable eQTM mapper for thousands of molecular phenotypes with graphical processing units.

BACKGROUND: Gene expression may be regulated by the DNA methylation of regulatory elements in cis, distal, and trans regions. One method to evaluate the relationship between DNA methylation and gene expression is the mapping of expression quantitative trait methylation (eQTM) loci (also called expression associated CpG loci, eCpG). However, no open-source tools are available to provide eQTM mapping. In addition, eQTM mapping can involve a large number of comparisons which may prevent the analyses due to limitations of computational resources. Here, we describe Torch-eCpG, an open-source tool to perform eQTM mapping that includes an optimized implementation that can use the graphical processing unit (GPU) to reduce runtime. RESULTS: We demonstrate the analyses using the tool are reproducible, up to 18 × faster using the GPU, and scale linearly with increasing methylation loci. CONCLUSIONS: Torch-eCpG is a fast, reliable, and scalable tool to perform eQTM mapping. Source code for Torch-eCpG is available at https://github.com/kordk/torch-ecpg .

Neurodegenerative disease pathways are perturbed in patients with cancer who self-report cognitive changes and anxiety: A pathway impact analysis.

BACKGROUND: Cancer-related cognitive impairment (CRCI) and anxiety co-occur in patients with cancer. Little is known about mechanisms for the co-occurrence of these two symptoms. The purposes of this secondary analysis were to evaluate for perturbed pathways associated with the co-occurrence of self-reported CRCI and anxiety in patients with low versus high levels of these two symptoms and to identify potential mechanisms for the co-occurrence of CRCI and anxiety using biological processes common across any perturbed neurodegenerative disease pathways. METHODS: Patients completed the Attentional Function Index and the Spielberger State-Trait Anxiety Inventory six times over two cycles of chemotherapy. Based on findings from a previous latent profile analysis, patients were grouped into none versus both high levels of these symptoms. Gene expression was quantified, and pathway impact analyses were performed. Signaling pathways for evaluation were defined with the Kyoto Encyclopedia of Genes and Genomes database. RESULTS: A total of 451 patients had data available for analysis. Approximately 85.0% of patients were in the none class and 15.0% were in the both high class. Pathway impact analyses identified five perturbed pathways related to neurodegenerative diseases (i.e., amyotrophic lateral sclerosis, Huntington disease, Parkinson disease, prion disease, and pathways of neurodegeneration-multiple diseases). Apoptosis, mitochondrial dysfunction, oxidative stress, and endoplasmic reticulum stress were common biological processes across these pathways. CONCLUSIONS: This study is the first to describe perturbations in neurodegenerative disease pathways associated with CRCI and anxiety in patients receiving chemotherapy. These findings provide new insights into potential targets for the development of mechanistically based interventions.

Perturbations in inflammatory pathways are associated with shortness of breath profiles in oncology patients receiving chemotherapy.

PURPOSE: One plausible mechanistic hypothesis is the potential contribution of inflammatory mechanisms to shortness of breath. This study was aimed to evaluate for associations between the occurrence of shortness of breath and perturbations in inflammatory pathways. METHODS: Patients with cancer reported the occurrence of shortness of breath six times over two cycles of chemotherapy. Latent class analysis was used to identify subgroups of patients with distinct shortness of breath occurrence profiles (i.e., none (70.5%), decreasing (8.2%), increasing (7.8%), high (13.5%)). Using an extreme phenotype approach, whole transcriptome differential gene expression and pathway impact analyses were performed to evaluate for perturbed signaling pathways associated with shortness of breath between the none and high classes. Two independent samples (RNA-sequencing (n = 293) and microarray (n = 295) methodologies) were evaluated. Fisher's combined probability method was used to combine these results to obtain a global test of the null hypothesis. In addition, an unweighted knowledge network was created using the specific pathway maps to evaluate for interconnections among these pathways. RESULTS: Twenty-nine Kyoto Encyclopedia of Genes and Genomes inflammatory signaling pathways were perturbed. The mitogen-activated protein kinase signaling pathway node had the highest closeness, betweenness, and degree scores. In addition, five common respiratory disease-related pathways, that may share mechanisms with cancer-related shortness of breath, were perturbed. CONCLUSIONS: Findings provide preliminary support for the hypothesis that inflammation contribute to the occurrence of shortness of breath in patients with cancer. In addition, the mechanisms that underlie shortness of breath in oncology patients may be similar to other respiratory diseases.

Venetoclax and dinaciclib elicit synergistic preclinical efficacy against hypodiploid acute lymphoblastic leukemia.

Hypodiploid acute lymphoblastic leukemia (ALL) is an aggressive blood cancer with a poor prognosis despite intensive chemotherapy or stem cell transplant. Children and adolescents with positive end-of-induction minimal residual disease have an overall survival lower than 30%. However, data regarding therapeutic alternatives for this disease is nearly nonexistent, emphasizing the critical need for new or adjunctive therapies that can improve outcomes. We previously reported on the therapeutic efficacy of venetoclax (ABT-199) in hypodiploid B-lineage ALL but with limitations as monotherapy. In this study, we set out to identify drugs enhancing the anti-leukemic effect of venetoclax in hypodiploid ALL. Using a highthroughput drug screen, we identified dinaciclib, a cyclin-dependent kinase inhibitor that worked synergistically with venetoclax to induce cell death in hypodiploid cell lines. This combination eradicated leukemic blasts within hypodiploid ALL patient-derived xenografts mice with low off-target toxicity. Our findings suggest that dual inhibition of BCL-2 (venetoclax) and CDK9/MCL-1 (dinaciclib) is a promising therapeutic approach in hypodiploid ALL, warranting further investigation to inform clinical trials in this high-risk patient population.

Prediction of morning fatigue severity in outpatients receiving chemotherapy: less may still be more.

INTRODUCTION: Fatigue is the most common and debilitating symptom experienced by cancer patients undergoing chemotherapy (CTX). Prediction of symptom severity can assist clinicians to identify high-risk patients and provide education to decrease symptom severity. The purpose of this study was to predict the severity of morning fatigue in the week following the administration of CTX. METHODS: Outpatients (n = 1217) completed questionnaires 1 week prior to and 1 week following administration of CTX. Morning fatigue was measured using the Lee Fatigue Scale (LFS). Separate prediction models for morning fatigue severity were created using 157 demographic, clinical, symptom, and psychosocial adjustment characteristics and either morning fatigue scores or individual fatigue item scores. Prediction models were created using two regression and five machine learning approaches. RESULTS: Elastic net models provided the best fit across all models. For the EN model using individual LFS item scores, two of the 13 individual LFS items (i.e., "worn out," "exhausted") were the strongest predictors. CONCLUSIONS: This study is the first to use machine learning techniques to accurately predict the severity of morning fatigue from prior to through the week following the administration of CTX using total and individual item scores from the Lee Fatigue Scale (LFS). Our findings suggest that the language used to assess clinical fatigue in oncology patients is important and that two simple questions may be used to predict morning fatigue severity.

Next-Generation Sequencing of Retinoblastoma Identifies Pathogenic Alterations beyond RB1 Inactivation That Correlate with Aggressive Histopathologic Features.

PURPOSE: To determine the usefulness of a comprehensive, targeted-capture next-generation sequencing (NGS) assay for the clinical management of children undergoing enucleation for retinoblastoma. DESIGN: Cohort study. PARTICIPANTS: Thirty-two children with retinoblastoma. METHODS: We performed targeted NGS using the UCSF500 Cancer Panel (University of California, San Francisco, San Francisco, CA) on formalin-fixed, paraffin-embedded tumor tissue along with constitutional DNA isolated from peripheral blood, buccal swab, or uninvolved optic nerve. Peripheral blood samples were also sent to a commercial laboratory for germline RB1 mutation testing. MAIN OUTCOME MEASURES: Presence or absence of germline RB1 mutation or deletion, tumor genetic profile, and association of genetic alterations with clinicopathologic features. RESULTS: Germline mutation or deletion of the RB1 gene was identified in all children with bilateral retinoblastoma (n = 12), and these NGS results were 100% concordant with commercial germline RB1 mutation analysis. In tumor tissue tested with NGS, biallelic inactivation of RB1 was identified in 28 tumors and focal MYCN amplification was identified in 4 tumors (2 with wild-type RB1 and 2 with biallelic RB1 inactivation). Additional likely pathogenic alterations beyond RB1 were identified in 13 tumors (41%), several of which have not been reported previously in retinoblastoma. These included focal amplifications of MDM4 and RAF1, as well as damaging mutations involving BCOR, ARID1A, MGA, FAT1, and ATRX. The presence of additional likely pathogenetic mutations beyond RB1 inactivation was associated with aggressive histopathologic features, including higher histologic grade and anaplasia, and also with both unilateral and sporadic disease. CONCLUSIONS: Comprehensive NGS analysis reliably detects relevant mutations, amplifications, and chromosomal copy number changes in retinoblastoma. The presence of genetic alterations beyond RB1 inactivation correlates with aggressive histopathologic features.

Maternal BMI at the start of pregnancy and offspring epigenome-wide DNA methylation: findings from the pregnancy and childhood epigenetics (PACE) consortium.

Pre-pregnancy maternal obesity is associated with adverse offspring outcomes at birth and later in life. Individual studies have shown that epigenetic modifications such as DNA methylation could contribute. Within the Pregnancy and Childhood Epigenetics (PACE) Consortium, we meta-analysed the association between pre-pregnancy maternal BMI and methylation at over 450,000 sites in newborn blood DNA, across 19 cohorts (9,340 mother-newborn pairs). We attempted to infer causality by comparing the effects of maternal versus paternal BMI and incorporating genetic variation. In four additional cohorts (1,817 mother-child pairs), we meta-analysed the association between maternal BMI at the start of pregnancy and blood methylation in adolescents. In newborns, maternal BMI was associated with small (<0.2% per BMI unit (1 kg/m2), P < 1.06 × 10-7) methylation variation at 9,044 sites throughout the genome. Adjustment for estimated cell proportions greatly attenuated the number of significant CpGs to 104, including 86 sites common to the unadjusted model. At 72/86 sites, the direction of the association was the same in newborns and adolescents, suggesting persistence of signals. However, we found evidence for acausal intrauterine effect of maternal BMI on newborn methylation at just 8/86 sites. In conclusion, this well-powered analysis identified robust associations between maternal adiposity and variations in newborn blood DNA methylation, but these small effects may be better explained by genetic or lifestyle factors than a causal intrauterine mechanism. This highlights the need for large-scale collaborative approaches and the application of causal inference techniques in epigenetic epidemiology.

Epigenetic remodeling in B-cell acute lymphoblastic leukemia occurs in two tracks and employs embryonic stem cell-like signatures.

We investigated DNA methylomes of pediatric B-cell acute lymphoblastic leukemias (B-ALLs) using whole-genome bisulfite sequencing and high-definition microarrays, along with RNA expression profiles. Epigenetic alteration of B-ALLs occurred in two tracks: de novo methylation of small functional compartments and demethylation of large inter-compartmental backbones. The deviations were exaggerated in lamina-associated domains, with differences corresponding to methylation clusters and/or cytogenetic groups. Our data also suggested a pivotal role of polycomb and CTBP2 in de novo methylation, which may be traced back to bivalency status of embryonic stem cells. Driven by these potent epigenetic modulations, suppression of polycomb target genes was observed along with disruption of developmental fate and cell cycle and mismatch repair pathways and altered activities of key upstream regulators.

Genetic variants of Adam17 differentially regulate TGFß signaling to modify vascular pathology in mice and humans.

Outcome of TGFß1 signaling is context dependent and differs between individuals due to germ-line genetic variation. To explore innate genetic variants that determine differential outcome of reduced TGFß1 signaling, we dissected the modifier locus Tgfbm3, on mouse chromosome 12. On a NIH/OlaHsd genetic background, the Tgfbm3b(C57) haplotype suppresses prenatal lethality of Tgfb1(-/-) embryos and enhances nuclear accumulation of mothers against decapentaplegic homolog 2 (Smad2) in embryonic cells. Amino acid polymorphisms within a disintegrin and metalloprotease 17 (Adam17) can account, at least in part, for this Tgfbm3b effect. ADAM17 is known to down-regulate Smad2 signaling by shedding the extracellular domain of TGFßRI, and we show that the C57 variant is hypomorphic for down-regulation of Smad2/3-driven transcription. Genetic variation at Tgfbm3 or pharmacological inhibition of ADAM17, modulates postnatal circulating endothelial progenitor cell (CEPC) numbers via effects on TGFßRI activity. Because CEPC numbers correlate with angiogenic potential, this suggests that variant Adam17 is an innate modifier of adult angiogenesis, acting through TGFßR1. To determine whether human ADAM17 is also polymorphic and interacts with TGFß signaling in human vascular disease, we investigated hereditary hemorrhagic telangiectasia (HHT), which is caused by mutations in TGFß/bone morphogenetic protein receptor genes, ENG, encoding endoglin (HHT1), or ACVRL1 encoding ALK1 (HHT2), and considered a disease of excessive abnormal angiogenesis. HHT manifests highly variable incidence and severity of clinical features, ranging from small mucocutaneous telangiectases to life-threatening visceral and cerebral arteriovenous malformations (AVMs). We show that ADAM17 SNPs associate with the presence of pulmonary AVM in HHT1 but not HHT2, indicating genetic variation in ADAM17 can potentiate a TGFß-regulated vascular disease.

Associations between circulating carotenoids, genomic instability and the risk of high-grade prostate cancer.

BACKGROUND: Carotenoids are a class of nutrients with antioxidant properties that have been purported to protect against cancer. However, the reported associations between carotenoids and prostate cancer have been heterogeneous and lacking data on interactions with nucleotide sequence variations and genomic biomarkers. OBJECTIVE: To examine the associations between carotenoid levels and the risk of high-grade prostate cancer, also considering antioxidant-related genes and tumor instability. METHODS: We measured plasma levels of carotenoids and genotyped 20 single nucleotide polymorphisms (SNP) in SOD1, SOD2, SOD3, XRCC1, and OGG1 among 559 men with non-metastatic prostate cancer undergoing radical prostatectomy. We performed copy number analysis in a subset of these men (n = 67) to study tumor instability assessed as Fraction of the Genome Altered (FGA). We examined associations between carotenoids, genotypes, tumor instability and risk of high-grade prostate cancer (Gleason grade ≥ 4 + 3) using logistic and linear regression. RESULTS: Circulating carotenoid levels were inversely associated with the risk of high-grade prostate cancer; odds ratios (OR) and 95% confidence intervals (CI) comparing highest versus lowest quartiles were: 0.34 (95% CI: 0.18-0.66) for α-carotene, 0.31 (95% CI: 0.15-0.63) for ß-carotene, 0.55 (0.28-1.08) for lycopene and 0.37 (0.18-0.75) for total carotenoids. SNPs rs25489 in XRCC1, rs699473 in SOD3 and rs1052133 in OGG1 modified these associations for α-carotene, ß-carotene and lycopene, respectively (P ≤ 0.05). The proportion of men with a high degree of FGA increased with Gleason Score (P < 0.001). Among men with Gleason score ≤ 3 + 4, higher lycopene levels were associated with lower FGA (P = 0.04). CONCLUSION: Circulating carotenoids at diagnosis, particularly among men carrying specific somatic variations, were inversely associated with risk of high-grade prostate cancer. In exploratory analyses, higher lycopene level was associated with less genomic instability among men with low-grade disease which is novel and supports the hypothesis that lycopene may inhibit progression of prostate cancer early in its natural history.

Serial analysis of 3D H-1 MRSI for patients with newly diagnosed GBM treated with combination therapy that includes bevacizumab.

Interpretation of changes in the T1- and T2-weighted MR images from patients with newly diagnosed glioblastoma (GBM) treated with standard of care in conjunction with anti-angiogenic agents is complicated by pseudoprogression and pseudoresponse. The hypothesis being tested in this study was that 3D H-1 magnetic resonance spectroscopic imaging (MRSI) provides estimates of levels of choline, creatine, N-acetylaspartate (NAA), lactate and lipid that change in response to treatment and that metrics describing these characteristics are associated with survival. Thirty-one patients with newly diagnosed GBM and being treated with radiation therapy (RT), temozolomide, erlotinib and bevacizumab were recruited to receive serial MR scans that included 3-D lactate edited MRSI at baseline, mid-RT, post-RT and at specific follow-up time points. The data were processed to provide estimates of metrics representing changes in metabolite levels relative to normal appearing brain. Cox proportional hazards analysis was applied to examine the relationship of these parameters with progression free survival (PFS) and overall survival (OS). There were significant reductions in parameters that describe relative levels of choline to NAA and creatine, indicating that the treatment caused a decrease in tumor cellularity. Changes in the levels of lactate and lipid relative to the NAA from contralateral brain were consistent with vascular normalization. Metabolic parameters from the first serial follow-up scan were associated with PFS and OS, when accounting for age and extent of resection. Integrating metabolic parameters into the assessment of patients with newly diagnosed GBM receiving therapies that include anti-angiogenic agents may be helpful for tracking changes in tumor burden, resolving ambiguities in anatomic images caused by non-specific treatment effects and for predicting outcome.

Transcription and methylation analyses of preleukemic promyelocytes indicate a dual role for PML/RARA in leukemia initiation.

Acute promyelocytic leukemia is an aggressive malignancy characterized by the accumulation of promyelocytes in the bone marrow. PML/RARA is the primary abnormality implicated in this pathology, but the mechanisms by which this chimeric fusion protein initiates disease are incompletely understood. Identifying PML/RARA targets in vivo is critical for comprehending the road to pathogenesis. Utilizing a novel sorting strategy, we isolated highly purified promyelocyte populations from normal and young preleukemic animals, carried out microarray and methylation profiling analyses, and compared the results from the two groups of animals. Surprisingly, in the absence of secondary lesions, PML/RARA had an overall limited impact on both the transcriptome and methylome. Of interest, we did identify down-regulation of secondary and tertiary granule genes as the first step engaging the myeloid maturation block. Although initially not sufficient to arrest terminal granulopoiesis in vivo, such alterations set the stage for the later, complete differentiation block seen in leukemia. Further, gene set enrichment analysis revealed that PML/RARA promyelocytes exhibit a subtle increase in expression of cell cycle genes, and we show that this leads to both increased proliferation of these cells and expansion of the promyelocyte compartment. Importantly, this proliferation signature was absent from the poorly leukemogenic p50/RARA fusion model, implying a critical role for PML in the altered cell-cycle kinetics and ability to initiate leukemia. Thus, our findings challenge the predominant model in the field and we propose that PML/RARA initiates leukemia by subtly shifting cell fate decisions within the promyelocyte compartment.

PTPRG inhibition by DNA methylation and cooperation with RAS gene activation in childhood acute lymphoblastic leukemia.

While the cytogenetic and genetic characteristics of childhood acute lymphoblastic leukemias (ALL) are well studied, less clearly understood are the contributing epigenetic mechanisms that influence the leukemia phenotype. Our previous studies and others identified gene mutation (RAS) and DNA methylation (FHIT) to be associated with the most common cytogenetic subgroup of childhood ALL, high hyperdiploidy (having five more chromosomes). We screened DNA methylation profiles, using a genome-wide high-dimension platform of 166 childhood ALLs and 6 normal pre-B cell samples and observed a strong association of DNA methylation status at the PTPRG locus in human samples with levels of PTPRG gene expression as well as with RAS gene mutation status. In the 293 cell line, we found that PTPRG expression induces dephosphorylation of ERK, a downstream RAS target that may be critical for mutant RAS-induced cell growth. In addition, PTPRG expression is upregulated by RAS activation under DNA hypomethylating conditions. An element within the PTPRG promoter is bound by the RAS-responsive transcription factor RREB1, also under hypomethylating conditions. In conclusion, we provide evidence that DNA methylation of the PTPRG gene is a complementary event in oncogenesis induced by RAS mutations. Evidence for additional roles for PTPR family member genes is also suggested. This provides a potential therapeutic target for RAS-related leukemias as well as insight into childhood ALL etiology and pathophysiology.

Diversity of antigen-specific responses induced in vivo with CTLA-4 blockade in prostate cancer patients.

CTLA-4 is a surface receptor on activated T cells that delivers an inhibitory signal, serving as an immune checkpoint. Treatment with anti-CTLA-4 Abs can induce clinical responses to different malignancies, but the nature of the induced Ag-specific recognition is largely unknown. Using microarrays spotted with >8000 human proteins, we assessed the diversity of Ab responses modulated by treatment with CTLA-4 blockade and GM-CSF. We find that advanced prostate cancer patients who clinically respond to treatment also develop enhanced Ab responses to a higher number of Ags than nonresponders. These induced Ab responses targeted Ags to which preexisting Abs are more likely to be present in the clinical responders compared with nonresponders. The majority of Ab responses are patient-specific, but immune responses against Ags shared among clinical responders are also detected. One of these shared Ags is PAK6, which is expressed in prostate cancer and to which CD4(+) T cell responses were also induced. Moreover, immunization with PAK6 can be both immunogenic and protective in mouse tumor models. These results demonstrate that immune checkpoint blockade modulates Ag-specific responses to both individualized and shared Ags, some of which can mediate anti-tumor responses.

Multiomic single cell sequencing identifies stemlike nature of mixed phenotype acute leukemia.

Despite recent work linking mixed phenotype acute leukemia (MPAL) to certain genetic lesions, specific driver mutations remain undefined for a significant proportion of patients and no genetic subtype is predictive of clinical outcomes. Moreover, therapeutic strategy for MPAL remains unclear, and prognosis is overall poor. We performed multiomic single cell profiling of 14 newly diagnosed adult MPAL patients to characterize the inter- and intra-tumoral transcriptional, immunophenotypic, and genetic landscapes of MPAL. We show that neither genetic profile nor transcriptome reliably correlate with specific MPAL immunophenotypes. Despite this, we find that MPAL blasts express a shared stem cell-like transcriptional profile indicative of high differentiation potential. Patients with the highest differentiation potential demonstrate inferior survival in our dataset. A gene set score, MPAL95, derived from genes highly enriched in the most stem-like MPAL cells, is applicable to bulk RNA sequencing data and is predictive of survival in an independent patient cohort, suggesting a potential strategy for clinical risk stratification.

Mutations in GNA11 in uveal melanoma.

BACKGROUND: Uveal melanoma is the most common intraocular cancer. There are no effective therapies for metastatic disease. Mutations in GNAQ, the gene encoding an alpha subunit of heterotrimeric G proteins, are found in 40% of uveal melanomas. METHODS: We sequenced exon 5 of GNAQ and GNA11, a paralogue of GNAQ, in 713 melanocytic neoplasms of different types (186 uveal melanomas, 139 blue nevi, 106 other nevi, and 282 other melanomas). We sequenced exon 4 of GNAQ and GNA11 in 453 of these samples and in all coding exons of GNAQ and GNA11 in 97 uveal melanomas and 45 blue nevi. RESULTS: We found somatic mutations in exon 5 (affecting Q209) and in exon 4 (affecting R183) in both GNA11 and GNAQ, in a mutually exclusive pattern. Mutations affecting Q209 in GNA11 were present in 7% of blue nevi, 32% of primary uveal melanomas, and 57% of uveal melanoma metastases. In contrast, we observed Q209 mutations in GNAQ in 55% of blue nevi, 45% of uveal melanomas, and 22% of uveal melanoma metastases. Mutations affecting R183 in either GNAQ or GNA11 were less prevalent (2% of blue nevi and 6% of uveal melanomas) than the Q209 mutations. Mutations in GNA11 induced spontaneously metastasizing tumors in a mouse model and activated the mitogen-activated protein kinase pathway. CONCLUSIONS: Of the uveal melanomas we analyzed, 83% had somatic mutations in GNAQ or GNA11. Constitutive activation of the pathway involving these two genes appears to be a major contributor to the development of uveal melanoma. (Funded by the National Institutes of Health and others.).

Germline variation controls the architecture of somatic alterations in tumors.

Studies have suggested that somatic events in tumors can depend on an individual's constitutional genotype. We used squamous cell carcinomas (SCC) of the skin, which arise in high multiplicity in organ transplant recipients, as a model to compare the pattern of somatic alterations within and across individuals. Specifically, we performed array comparative genomic hybridization on 104 tumors from 25 unrelated individuals who each had three or more independently arisen SCCs and compared the profiles occurring within patients to profiles of tumors across a larger set of 135 patients. In general, chromosomal aberrations in SCCs were more similar within than across individuals (two-sided exact-test p-value<1x10(-7)), consistent with the notion that the genetic background was affecting the pattern of somatic changes. To further test this possibility, we performed allele-specific imbalance studies using microsatellite markers mapping to 14 frequently aberrant regions of multiple independent tumors from 65 patients. We identified nine loci which show evidence of preferential allelic imbalance. One of these loci, 8q24, corresponded to a region in which multiple single nucleotide polymorphisms have been associated with increased cancer risk in genome-wide association studies (GWAS). We tested three implicated variants and identified one, rs13281615, with evidence of allele-specific imbalance (p-value=0.012). The finding of an independently identified cancer susceptibility allele with allele-specific imbalance in a genomic region affected by recurrent DNA copy number changes suggest that it may also harbor risk alleles for SCC. Together these data provide strong evidence that the genetic background is a key driver of somatic events in cancer, opening an opportunity to expand this approach to identify cancer risk alleles.

Torch-eCpG: A fast and scalable eQTM mapper for thousands of molecular phenotypes with graphical processing units.

Background: Gene expression may be regulated by the DNA methylation of regulatory elements in cis, distal, and trans regions. One method to evaluate the relationship between DNA methylation and gene expression is the mapping of expression quantitative trait methylation (eQTM) loci (also called expression associated CpG loci, eCpG). However, no open-source tools are available to provide eQTM mapping. In addition, eQTM mapping can involve a large number of comparisons which may prevent the analyses due to limitations of computational resources. Here, we describe Torch-eCpG, an open-source tool to perform eQTM mapping that includes an optimized implementation that can use the graphical processing unit (GPU) to reduce runtime. Results: We demonstrate the analyses using the tool are reproducible, up to 18x faster using the GPU, and scale linearly with increasing methylation loci. Conclusions: Torch-eCpG is a fast, reliable, and scalable tool to perform eQTM mapping. Source code for Torch-eCpG is available at https://github.com/kordk/torch-ecpg.

Therapeutic implications of transcriptomics in head and neck cancer patient-derived xenografts.

There are currently no clinical strategies utilizing tumor gene expression to inform therapeutic selection for patients with head and neck squamous cell carcinoma (HNSCC). One of the challenges in developing predictive biomarkers is the limited characterization of preclinical HNSCC models. Patient-derived xenografts (PDXs) are increasingly recognized as translationally relevant preclinical avatars for human tumors; however, the overall transcriptomic concordance of HNSCC PDXs with primary human HNSCC is understudied, especially in human papillomavirus-associated (HPV+) disease. Here, we characterized 64 HNSCC PDXs (16 HPV+ and 48 HPV-) at the transcriptomic level using RNA-sequencing. The range of human-specific reads per PDX varied from 64.6%-96.5%, with a comparison of the most differentially expressed genes before and after removal of mouse transcripts revealing no significant benefit to filtering out mouse mRNA reads in this cohort. We demonstrate that four previously established HNSCC molecular subtypes found in The Cancer Genome Atlas (TCGA) are also clearly recapitulated in HNSCC PDXs. Unsupervised hierarchical clustering yielded a striking natural division of HNSCC PDXs by HPV status, with C19orf57 (BRME1), a gene previously correlated with positive response to cisplatin in cervical cancer, among the most significantly differentially expressed genes between HPV+ and HPV- PDXs. In vivo experiments demonstrated a possible relationship between increased C19orf57 expression and superior anti-tumor responses of PDXs to cisplatin, which should be investigated further. These findings highlight the value of PDXs as models for HPV+ and HPV- HNSCC, providing a resource for future discovery of predictive biomarkers to guide treatment selection in HNSCC.

Multiomic Single Cell Sequencing Identifies Stemlike Nature of Mixed Phenotype Acute Leukemia and Provides Novel Risk Stratification.

Mixed phenotype acute leukemia (MPAL) is a leukemia whose biologic drivers are poorly understood, therapeutic strategy remains unclear, and prognosis is poor. We performed multiomic single cell (SC) profiling of 14 newly diagnosed adult MPAL patients to characterize the immunophenotypic, genetic, and transcriptional landscapes of MPAL. We show that neither genetic profile nor transcriptome reliably correlate with specific MPAL immunophenotypes. However, progressive acquisition of mutations is associated with increased expression of immunophenotypic markers of immaturity. Using SC transcriptional profiling, we find that MPAL blasts express a stem cell-like transcriptional profile distinct from other acute leukemias and indicative of high differentiation potential. Further, patients with the highest differentiation potential demonstrated inferior survival in our dataset. A gene set score, MPAL95, derived from genes highly enriched in this cohort, is applicable to bulk RNA sequencing data and was predictive of survival in an independent patient cohort, suggesting utility for clinical risk stratification.