Bedside Back to Bench: Building Bridges between Basic and Clinical Genomic Research.

Genome sequencing has revolutionized the diagnosis of genetic diseases. Close collaborations between basic scientists and clinical genomicists are now needed to link genetic variants with disease causation. To facilitate such collaborations, we recommend prioritizing clinically relevant genes for functional studies, developing reference variant-phenotype databases, adopting phenotype description standards, and promoting data sharing.

Common variation in a long non-coding RNA gene modulates variation of circulating TGF-ß2 levels in metastatic colorectal cancer patients (Alliance).

BACKGROUND: Herein, we report results from a genome-wide study conducted to identify protein quantitative trait loci (pQTL) for circulating angiogenic and inflammatory protein markers in patients with metastatic colorectal cancer (mCRC). The study was conducted using genotype, protein marker, and baseline clinical and demographic data from CALGB/SWOG 80405 (Alliance), a randomized phase III study designed to assess outcomes of adding VEGF or EGFR inhibitors to systemic chemotherapy in mCRC patients. Germline DNA derived from blood was genotyped on whole-genome array platforms. The abundance of protein markers was quantified using a multiplex enzyme-linked immunosorbent assay from plasma derived from peripheral venous blood collected at baseline. A robust rank-based method was used to assess the statistical significance of each variant and protein pair against a strict genome-wide level. A given pQTL was tested for validation in two external datasets of prostate (CALGB 90401) and pancreatic cancer (CALGB 80303) patients. Bioinformatics analyses were conducted to further establish biological bases for these findings. RESULTS: The final analysis was carried out based on data from 540,021 common typed genetic variants and 23 protein markers from 869 genetically estimated European patients with mCRC. Correcting for multiple testing, the analysis discovered a novel cis-pQTL in LINC02869, a long non-coding RNA gene, for circulating TGF-ß2 levels (rs11118119; AAF = 0.11; P-value < 1.4e-14). This finding was validated in a cohort of 538 prostate cancer patients from CALGB 90401 (AAF = 0.10, P-value < 3.3e-25). The analysis also validated a cis-pQTL we had previously reported for VEGF-A in advanced pancreatic cancer, and additionally identified trans-pQTLs for VEGF-R3, and cis-pQTLs for CD73. CONCLUSIONS: This study has provided evidence of a novel cis germline genetic variant that regulates circulating TGF-ß2 levels in plasma of patients with advanced mCRC and prostate cancer. Moreover, the validation of previously identified pQTLs for VEGF-A, CD73, and VEGF-R3, potentiates the validity of these associations.

Dihydropyrimidine dehydrogenase gene variants for predicting grade 4-5 fluoropyrimidine-induced toxicity: FUSAFE individual patient data meta-analysis.

BACKGROUND: Dihydropyrimidine dehydrogenase (DPD) deficiency is the main known cause of life-threatening fluoropyrimidine (FP)-induced toxicities. We conducted a meta-analysis on individual patient data to assess the contribution of deleterious DPYD variants *2A/D949V/*13/HapB3 (recommended by EMA) and clinical factors, for predicting G4-5 toxicity. METHODS: Study eligibility criteria included recruitment of Caucasian patients without DPD-based FP-dose adjustment. Main endpoint was 12-week haematological or digestive G4-5 toxicity. The value of DPYD variants *2A/p.D949V/*13 merged, HapB3, and MIR27A rs895819 was evaluated using multivariable logistic models (AUC). RESULTS: Among 25 eligible studies, complete clinical variables and primary endpoint were available in 15 studies (8733 patients). Twelve-week G4-5 toxicity prevalence was 7.3% (641 events). The clinical model included age, sex, body mass index, schedule of FP-administration, concomitant anticancer drugs. Adding *2A/p.D949V/*13 variants (at least one allele, prevalence 2.2%, OR 9.5 [95%CI 6.7-13.5]) significantly improved the model (p < 0.0001). The addition of HapB3 (prevalence 4.0%, 98.6% heterozygous), in spite of significant association with toxicity (OR 1.8 [95%CI 1.2-2.7]), did not improve the model. MIR27A rs895819 was not associated with toxicity, irrespective of DPYD variants. CONCLUSIONS: FUSAFE meta-analysis highlights the major relevance of DPYD *2A/p.D949V/*13 combined with clinical variables to identify patients at risk of very severe FP-related toxicity.

Pharmacogenetic and clinical risk factors for bevacizumab-related gastrointestinal hemorrhage in prostate cancer patients treated on CALGB 90401 (Alliance).

The objective of this study was to discover clinical and pharmacogenetic factors associated with bevacizumab-related gastrointestinal hemorrhage in Cancer and Leukemia Group B (Alliance) 90401. Patients with metastatic castration-resistant prostate cancer received docetaxel and prednisone ± bevacizumab. Patients were genotyped using Illumina HumanHap610-Quad and assessed using cause-specific risk for association between single nucleotide polymorphisms (SNPs) and gastrointestinal hemorrhage. In 1008 patients, grade 2 or higher gastrointestinal hemorrhage occurred in 9.5% and 3.8% of bevacizumab (n = 503) and placebo (n = 505) treated patients, respectively. Bevacizumab (P < 0.001) and age (P = 0.002) were associated with gastrointestinal hemorrhage. In 616 genetically estimated Europeans (n = 314 bevacizumab and n = 302 placebo treated patients), grade 2 or higher gastrointestinal hemorrhage occurred in 9.6% and 2.0% of patients, respectively. One SNP (rs1478947; HR 6.26; 95% CI 3.19-12.28; P = 9.40 × 10-8) surpassed Bonferroni-corrected significance. Grade 2 or higher gastrointestinal hemorrhage rate was 33.3% and 6.2% in bevacizumab-treated patients with the AA/AG and GG genotypes, versus 2.9% and 1.9% in the placebo arm, respectively. Prospective validation of these findings and functional analyses are needed to better understand the genetic contribution to treatment-related gastrointestinal hemorrhage.

Deceptive learning in histopathology.

AIMS: Deep learning holds immense potential for histopathology, automating tasks that are simple for expert pathologists and revealing novel biology for tasks that were previously considered difficult or impossible to solve by eye alone. However, the extent to which the visual strategies learned by deep learning models in histopathological analysis are trustworthy or not has yet to be systematically analysed. Here, we systematically evaluate deep neural networks (DNNs) trained for histopathological analysis in order to understand if their learned strategies are trustworthy or deceptive. METHODS AND RESULTS: We trained a variety of DNNs on a novel data set of 221 whole-slide images (WSIs) from lung adenocarcinoma patients, and evaluated their effectiveness at (1) molecular profiling of KRAS versus EGFR mutations, (2) determining the primary tissue of a tumour and (3) tumour detection. While DNNs achieved above-chance performance on molecular profiling, they did so by exploiting correlations between histological subtypes and mutations, and failed to generalise to a challenging test set obtained through laser capture microdissection (LCM). In contrast, DNNs learned robust and trustworthy strategies for determining the primary tissue of a tumour as well as detecting and localising tumours in tissue. CONCLUSIONS: Our work demonstrates that DNNs hold immense promise for aiding pathologists in analysing tissue. However, they are also capable of achieving seemingly strong performance by learning deceptive strategies that leverage spurious correlations, and are ultimately unsuitable for research or clinical work. The framework we propose for model evaluation and interpretation is an important step towards developing reliable automated systems for histopathological analysis.

High-throughput screening and genome-wide analyses of 44 anticancer drugs in the 1000 Genomes cell lines reveals an association of the NQO1 gene with the response of multiple anticancer drugs.

Cancer patients exhibit a broad range of inter-individual variability in response and toxicity to widely used anticancer drugs, and genetic variation is a major contributor to this variability. To identify new genes that influence the response of 44 FDA-approved anticancer drug treatments widely used to treat various types of cancer, we conducted high-throughput screening and genome-wide association mapping using 680 lymphoblastoid cell lines from the 1000 Genomes Project. The drug treatments considered in this study represent nine drug classes widely used in the treatment of cancer in addition to the paclitaxel + epirubicin combination therapy commonly used for breast cancer patients. Our genome-wide association study (GWAS) found several significant and suggestive associations. We prioritized consistent associations for functional follow-up using gene-expression analyses. The NAD(P)H quinone dehydrogenase 1 (NQO1) gene was found to be associated with the dose-response of arsenic trioxide, erlotinib, trametinib, and a combination treatment of paclitaxel + epirubicin. NQO1 has previously been shown as a biomarker of epirubicin response, but our results reveal novel associations with these additional treatments. Baseline gene expression of NQO1 was positively correlated with response for 43 of the 44 treatments surveyed. By interrogating the functional mechanisms of this association, the results demonstrate differences in both baseline and drug-exposed induction.

Biomarkers for Response to Anti-PD-1/Anti-PD-L1 Immune Checkpoint Inhibitors: A Large Meta-Analysis.

BACKGROUND: Immune checkpoint inhibitors (ICIs) that block PD-1/PD-L1 have consistently demonstrated durable clinical activity across multiple histologies but have low overall response rates for many cancers-indicating that too few patients benefit from ICIs. Many studies have explored potential predictive biomarkers (eg, PD-1/PD-L1 expression, tumor mutational burden [TMB]), no consensus biomarker has been identified. METHODS: This meta-analysis combined predictive accuracy metrics for various biomarkers, across multiple cancer types, to determine which biomarkers are most accurate for predicting ICI response. Data from 18,792 patients from 100 peer-reviewed studies that evaluated putative biomarkers for response to anti-PD-1/anti- PD-L1 treatment were meta-analyzed using bivariate linear mixed models. Biomarker performance was assessed based on the global area under the receiver operating characteristic curve (AUC) and 95% bootstrap confidence intervals. RESULTS: PD-L1 immunohistochemistry, TMB, and multimodal biomarkers discriminated responders and nonresponders better than random assignment (AUCs >.50). Excluding multimodal biomarkers, these biomarkers correctly classified at least 50% of the responders (sensitivity 95% CIs, >.50). Notably, variation in biomarker performance was observed across cancer types. CONCLUSIONS: Although some biomarkers consistently performed better, heterogeneity in performance was observed across cancer types, and additional research is needed to identify highly accurate and precise biomarkers for widespread clinical use.

A two-stage genome-wide association study to identify novel genetic loci associated with acute radiotherapy toxicity in nasopharyngeal carcinoma.

BACKGROUND: Genetic variants associated with acute side effects of radiotherapy in nasopharyngeal carcinoma (NPC) remain largely unknown. METHODS: We performed a two-stage genome-wide association analysis including a total of 1084 patients, where 319 individuals in the discovery stage were genotyped for 688,783 SNPs using whole genome-wide screening microarray. Significant variants were then validated in an independent cohort of 765 patients using the MassARRAY system. Gene mapping, linkage disequilibrium, genome-wide association analysis, and polygenic risk score were conducted or calculated using FUMA, LDBlockShow, PLINK, and PRSice software programs, respectively. RESULTS: Five SNPs (rs6711678, rs4848597, rs4848598, rs2091255, and rs584547) showed statistical significance after validation. Radiotherapy toxicity was more serious in mutant minor allele carriers of all five SNPs. Stratified analysis further indicated that rs6711678, rs4848597, rs4848598, and rs2091255 correlated with skin toxicity in patients of EBV positive, late stage (III and IV), receiving both concurrent chemoradiotherapy and induction/adjuvant chemotherapy, and with OR values ranging from 1.92 to 2.66. For rs584547, high occurrence of dysphagia was found in A allele carriers in both the discovery (P = 1.27 × 10- 6, OR = 1.55) and validation (P = 0.002, OR = 4.20) cohorts. Furthermore, prediction models integrating both genetic and clinical factors for skin reaction and dysphagia were established. The area under curve (AUC) value of receiver operating characteristic (ROC) curves were 0.657 (skin reaction) and 0.788 (dysphagia). CONCLUSIONS: Rs6711678, rs4848597, rs4848598, and rs2091255 on chromosome 2q14.2 and rs584547 were found to be novel risk loci for skin toxicity and dysphagia in NPC patients receiving radiotherapy. TRIAL REGISTRATION: Chinese Clinical Trial Register (registration number: ChiCTR-OPC-14005257 and CTXY-140007-2).

Bevacizumab-induced hypertension and proteinuria: a genome-wide study of more than 1000 patients.

BACKGROUND: Hypertension and proteinuria are common bevacizumab-induced toxicities. No validated biomarkers are available for identifying patients at risk of these toxicities. METHODS: A genome-wide association study (GWAS) meta-analysis was performed in 1039 bevacizumab-treated patients of European ancestry in four clinical trials (CALGB 40502, 40503, 80303, 90401). Grade ≥2 hypertension and proteinuria were recorded (CTCAE v.3.0). Single-nucleotide polymorphism (SNP)-toxicity associations were determined using a cause-specific Cox model adjusting for age and sex. RESULTS: The most significant SNP associated with hypertension with concordant effect in three out of the four studies (p-value <0.05 for each study) was rs6770663 (A > G) in KCNAB1, with the G allele increasing the risk of hypertension (p-value = 4.16 × 10-6). The effect of the G allele was replicated in ECOG-ACRIN E5103 in 582 patients (p-value = 0.005). The meta-analysis of all five studies for rs6770663 led to p-value = 7.73 × 10-8, close to genome-wide significance. The most significant SNP associated with proteinuria was rs339947 (C > A, between DNAH5 and TRIO), with the A allele increasing the risk of proteinuria (p-value = 1.58 × 10-7). CONCLUSIONS: The results from the largest study of bevacizumab toxicity provide new markers of drug safety for further evaluations. SNP in KCNAB1 validated in an independent dataset provides evidence toward its clinical applicability to predict bevacizumab-induced hypertension. ClinicalTrials.gov Identifier: NCT00785291 (CALGB 40502); NCT00601900 (CALGB 40503); NCT00088894 (CALGB 80303) and NCT00110214 (CALGB 90401).

Race and smoking status associated with paclitaxel drug response in patient-derived lymphoblastoid cell lines.

The use of ex-vivo model systems to provide a level of forecasting for in-vivo characteristics remains an important need for cancer therapeutics. The use of lymphoblastoid cell lines (LCLs) is an attractive approach for pharmacogenomics and toxicogenomics, due to their scalability, efficiency, and cost-effectiveness. There is little data on the impact of demographic or clinical covariates on LCL response to chemotherapy. Paclitaxel sensitivity was determined in LCLs from 93 breast cancer patients from the University of North Carolina Lineberger Comprehensive Cancer Center Breast Cancer Database to test for potential associations and/or confounders in paclitaxel dose-response assays. Measures of paclitaxel cell viability were associated with patient data included treatment regimens, cancer status, demographic and environmental variables, and clinical outcomes. We used multivariate analysis of variance to identify the in-vivo variables associated with ex-vivo dose-response. In this unique dataset that includes both in-vivo and ex-vivo data from breast cancer patients, race (P = 0.0049) and smoking status (P = 0.0050) were found to be significantly associated with ex-vivo dose-response in LCLs. Racial differences in clinical dose-response have been previously described, but the smoking association has not been reported. Our results indicate that in-vivo smoking status can influence ex-vivo dose-response in LCLs, and more precise measures of covariates may allow for more precise forecasting of clinical effect. In addition, understanding the mechanism by which exposure to smoking in-vivo effects ex-vivo dose-response in LCLs may open up new avenues in the quest for better therapeutic prediction.