Evaluating the utility of ZNF331 promoter methylation as a prognostic and predictive marker in stage III colon cancer: results from CALGB 89803 (Alliance).

While epigenomic alterations are common in colorectal cancers (CRC), few epigenomic biomarkers that risk-stratify patients have been identified. We thus sought to determine the potential of ZNF331 promoter hypermethylation (mZNF331) as a prognostic and predictive marker in colon cancer. We examined the association of mZNF331 with clinicopathologic features, relapse, survival, and treatment efficacy in patients with stage III colon cancer treated within a randomized adjuvant chemotherapy trial (CALGB/Alliance89803). Residual tumour tissue was available for genomic DNA extraction and methylation analysis for 385 patients. ZNF331 promoter methylation status was determined by bisulphite conversion and fluorescence-based real-time polymerase chain reaction. Kaplan-Meier estimator and Cox proportional hazard models were used to assess the prognostic and predictive role of mZNF331 in this well-annotated dataset, adjusting for clinicopathologic features and standard molecular markers. mZNF331 was observed in 267/385 (69.4%) evaluable cases. Histopathologic features were largely similar between patients with mZNF331 compared to unmethylated ZNF331 (unmZNFF31). There was no significant difference in disease-free or overall survival between patients with mZNF331 versus unmZNF331 colon cancers, even when adjusting for clinicopathologic features and molecular marker status. Similarly, there was no difference in disease-free or overall survival across treatment arms when stratified by ZNF331 methylation status. While ZNF331 promoter hypermethylation is frequently observed in CRC, our current study of a small subset of patients with stage III colon cancer suggests limited applicability as a prognostic marker. Larger studies may provide more insight and clarity into the applicability of mZNF331 as a prognostic and predictive marker.

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.

Germline cis variant determines epigenetic regulation of the anti-cancer drug metabolism gene dihydropyrimidine dehydrogenase (DPYD).

Enhancers are critical for regulating tissue-specific gene expression, and genetic variants within enhancer regions have been suggested to contribute to various cancer-related processes, including therapeutic resistance. However, the precise mechanisms remain elusive. Using a well-defined drug-gene pair, we identified an enhancer region for dihydropyrimidine dehydrogenase (DPD, DPYD gene) expression that is relevant to the metabolism of the anti-cancer drug 5-fluorouracil (5-FU). Using reporter systems, CRISPR genome-edited cell models, and human liver specimens, we demonstrated in vitro and vivo that genotype status for the common germline variant (rs4294451; 27% global minor allele frequency) located within this novel enhancer controls DPYD transcription and alters resistance to 5-FU. The variant genotype increases recruitment of the transcription factor CEBPB to the enhancer and alters the level of direct interactions between the enhancer and DPYD promoter. Our data provide insight into the regulatory mechanisms controlling sensitivity and resistance to 5-FU.

Gene expression biomarkers differentiate overall survival of colorectal cancer upon targeted therapies.

While monoclonal antibody-based targeted therapies have substantially improved progression-free survival in cancer patients, the variability in individual responses poses a significant challenge in patient care. Therefore, identifying cancer subtypes and their associated biomarkers is required for assigning effective treatment. In this study, we integrated genotype and pre-treatment tissue RNA-seq data and identified biomarkers causally associated with the overall survival (OS) of colorectal cancer (CRC) patients treated with either cetuximab or bevacizumab. We performed enrichment analysis for specific consensus molecular subtypes (CMS) of colorectal cancer and evaluated differential expression of identified genes using paired tumor and normal tissue from an external cohort. In addition, we replicated the causal effect of these genes on OS using validation cohort and assessed their association with the Cancer Genome Atlas Program data as an external cohort. One of the replicated findings was WDR62, whose overexpression shortened OS of patients treated with cetuximab. Enrichment of its over expression in CMS1 and low expression in CMS4 suggests that patients with CMS4 subtype may drive greater benefit from cetuximab. In summary, this study highlights the importance of integrating different omics data for identifying promising biomarkers specific to a treatment or a cancer subtype.

HER2 Gene Expression Levels Are Predictive and Prognostic in Patients With Metastatic Colorectal Cancer Enrolled in CALGB/SWOG 80405.

PURPOSE: The phase III Cancer and Leukemia Group B (CALGB)/SWOG 80405 trial found no difference in overall survival (OS) in patients with metastatic colorectal cancer receiving first-line chemotherapy in combination with either bevacizumab or cetuximab. We investigated the potential prognostic and predictive value of HER2 amplification and gene expression using next-generation sequencing (NGS) and NanoString data. PATIENTS AND METHODS: Primary tumor DNA from 559 patients was profiled for HER2 amplification by NGS (FoundationOne CDx). Tumor tissue from 925 patients was tested for NanoString gene expression using an 800-gene panel. OS and progression-free survival (PFS) were the time-to-event end points. RESULTS: High HER2 expression (dichotomized at median) was associated with longer PFS (11.6 v 10 months, P = .012) and OS (32 v 25.3 months, P = .033), independent of treatment. An OS benefit for cetuximab versus bevacizumab was observed in the high HER2 expression group (P = .02), whereas a worse PFS for cetuximab was seen in the low-expression group (P = .019). When modeled as a continuous variable, increased HER2 expression was associated with longer OS (hazard ratio [HR], 0.83 [95% CI, 0.75 to 0.93]; adjusted P = .0007) and PFS (HR, 0.82 [95% CI, 0.74 to 0.91]; adjusted P = .0002), reaching a plateau effect after the median. In patients with HER2 expression lower than median, treatment with cetuximab was associated with worse PFS (HR, 1.38 [95% CI, 1.12 to 1.71]; adjusted P = .0027) and OS (HR, 1.28 [95% CI, 1.02 to 1.59]; adjusted P = .03) compared with that with bevacizumab. A significant interaction between HER2 expression and the treatment arm was observed for OS (Pintx = .017), PFS (Pintx = .048), and objective response rate (Pintx = .001). CONCLUSION: HER2 gene expression was prognostic and predictive in CALGB/SWOG 80405. HER2 tumor expression may inform treatment selection for patients with low HER2 favoring bevacizumab- versus cetuximab-based therapies.

Model-Based Prediction of Irinotecan-Induced Grade 4 Neutropenia in Cancer Patients: Influence of Incorporating Germline Genetic Factors in the Model.

Neutropenia is the major dose-limiting toxicity of irinotecan-based therapy. The objective of this study was to assess whether inclusion of germline genetic variants into a population pharmacokinetic/pharmacodynamic model can improve prediction of irinotecan-induced grade 4 neutropenia and identify novel variants of clinical value. A semimechanistic population pharmacokinetic/pharmacodynamic model was used to predict neutrophil response over time in 197 patients receiving irinotecan. Covariate analysis was performed for demographic/clinical factors and 4,781 genetic variants in 84 drug response- and toxicity-related genes to identify covariates associated with neutrophil response. We evaluated the predictive value of the model for grade 4 neutropenia reflecting different clinical scenarios of available data on identified demographic/clinical covariates, baseline and post-treatment absolute neutrophil counts (ANCs), individual pharmacokinetics, and germline genetic variation. Adding 8 genetic identified covariates (rs10929302 (UGT1A1), rs1042482 (DPYD), rs2859101 (HLA-DQB3), rs61754806 (NR3C1), rs9266271 (HLA-B), rs7294 (VKORC1), rs1051713 (ALOX5), and ABCB1 rare variant burden) to a model using only baseline ANCs improved prediction of irinotecan-induced grade 4 neutropenia from area under the receiver operating characteristic curve (AUC-ROC) of 50-64% (95% confidence interval (CI), 54-74%). Individual pharmacokinetics further improved the prediction to 74% (95% CI, 64-84%). When weekly ANC was available, the identified covariates and individual pharmacokinetics yielded no additional contribution to the prediction. The model including only ANCs at baseline and at week 1 achieved an AUC-ROC of 78% (95% CI, 69-88%). Germline DNA genetic variants may contribute to the prediction of irinotecan-induced grade 4 neutropenia when incorporated into a population pharmacokinetic/pharmacodynamic model. This approach is generalizable to drugs that induce neutropenia and ultimately allows for personalized intervention to enhance patient safety.

CALGB 80802 (Alliance): Impact of Sorafenib with and without Doxorubicin on Hepatitis C Infection in Patients with Advanced Hepatocellular Carcinoma.

Sorafenib blocks nonstructural protein 5A (NS5A)-recruited c-Raf-mediated hepatitis C virus (HCV) replication and gene expression. Release of Raf-1-Ask-1 dimer and inhibition of Raf-1 via sorafenib putatively differ in the presence or absence of doxorubicin. Cancer and Leukemia Group B (CALGB) 80802 (Alliance) randomized phase III trial of doxorubicin plus sorafenib versus sorafenib in patients with advanced hepatocellular carcinoma (HCC), showed no improvement in median overall survival (OS). Whether HCV viral load impacts therapy and whether any correlation between HCV titers and outcome based on HCV was studied. In patients with HCV, HCV titer levels were evaluated at baseline and at multiple postbaseline timepoints until disease progression or treatment discontinuation. HCV titer levels were evaluated in relation to OS and progression-free survival (PFS). Among 53 patients with baseline HCV data, 12 patients had undetectable HCV (HCV-UN). Postbaseline HCV titer levels did not significantly differ between treatment arms. One patient in each arm went from detectable to HCV-UN with greater than 2 log-fold titer levels reduction. Aside from these 2 HCV-UN patients, HCV titers remained stable on treatment. Patients who had HCV-UN at baseline were 3.5 times more likely to progress and/or die from HCC compared with HCV detectable (HR = 3.51; 95% confidence interval: 1.58-7.78; P = 0.002). HCV titer levels remained unchanged, negating any sorafenib impact onto HCV titer levels. Although an overall negative phase III study, patients treated with doxorubicin plus sorafenib and sorafenib only, on CALGB 80802 had worse PFS if HCV-UN. Higher levels of HCV titers at baseline were associated with significantly improved PFS. SIGNIFICANCE: Sorafenib therapy for HCC may impact HCV replication and viral gene expression. In HCV-positive patients accrued to CLAGB 80802 phase III study evaluating the addition of doxorubicin to sorafenib, HCV titer levels were evaluated at baseline and different timepoints. Sorafenib did not impact HCV titer levels. Despite an improved PFS in patients with detectable higher level HCV titers at baseline, no difference in OS was noted.

Predictive value of CDC37 gene expression for targeted therapy in metastatic colorectal cancer.

BACKGROUND: CDC37 is a key determinant of client kinase recruitment to the HSP90 chaperoning system. We hypothesized that kinase-specific dependency on CDC37 alters the efficacy of targeted therapies for metastatic colorectal cancer (mCRC). MATERIAL AND METHODS: Two independent mCRC cohorts were analyzed to compare the survival outcomes between CDC37-high and CDC37-low patients (stratified by the median cutoff values): the CALGB/SWOG 80405 trial (226 and 207 patients receiving first-line bevacizumab- and cetuximab-containing chemotherapies, respectively) and Japanese retrospective (50 refractory patients receiving regorafenib) cohorts. A dataset of specimens submitted to a commercial CLIA-certified laboratory was utilized to characterize molecular profiles of CDC37-high (top quartile, N = 5055) and CDC37-low (bottom quartile, N = 5055) CRCs. RESULTS: In the bevacizumab-treated group, CDC37-high patients showed significantly better progression-free survival (PFS) (median 13.3 vs 9.6 months, hazard ratio [HR] 0.59, 95% confidence interval [CI] 0.44-0.79, p < 0.01) than CDC37-low patients. In the cetuximab-treated group, CDC37-high and CDC37-low patients had similar outcomes. In the regorafenib-treated group, CDC37-high patients showed significantly better overall survival (median 11.3 vs 6.0 months, HR 0.24, 95% CI 0.11-0.54, p < 0.01) and PFS (median 3.5 vs 1.9 months, HR 0.51, 95% CI 0.28-0.94, p = 0.03). Comprehensive molecular profiling revealed that CDC37-high CRCs were associated with higher VEGFA, FLT1, and KDR expressions and activated hypoxia signature. CONCLUSIONS: CDC37-high mCRC patients derived more benefit from anti-VEGF therapies, including bevacizumab and regorafenib, but not from cetuximab. Molecular profiles suggested that such tumors were dependent on angiogenesis-relating pathways.

CCR5 and CCL5 gene expression in colorectal cancer: comprehensive profiling and clinical value.

BACKGROUND: The C-C motif chemokine receptor 5 (CCR5)/C-C motif chemokine ligand 5 (CCL5) axis plays a major role in colorectal cancer (CRC). We aimed to characterize the molecular features associated with CCR5/CCL5 expression in CRC and to determine whether CCR5/CCL5 levels could impact treatment outcomes. METHODS: 7604 CRCs tested with NextGen Sequencing on DNA and RNA were analyzed. Molecular features were evaluated according to CCR5 and CCL5 tumor gene expression quartiles. The impact on treatment outcomes was assessed in two cohorts, including 6341 real-world patients and 429 patients from the Cancer and Leukemia Group B (CALGB)/SWOG 80405 trial. RESULTS: CCR5/CCL5 expression was higher in right-sided versus left-sided tumors, and positively associated with consensus molecular subtypes 1 and 4. Higher CCR5/CCL5 expression was associated with higher tumor mutational burden, deficiency in mismatch repair and programmed cell death ligand 1 (PD-L1) levels. Additionally, high CCR5/CCL5 were associated with higher immune cell infiltration in the tumor microenvironment (TME) of MMR proficient tumors. Ingenuity pathway analysis revealed upregulation of the programmed cell death protein 1 (PD-1)/PD-L1 cancer immunotherapy pathway, phosphatase and tensin homolog (PTEN) and peroxisome proliferator-activated receptors (PPAR) signaling, and cytotoxic T-lymphocyte antigen 4 (CTLA-4) signaling in cytotoxic T lymphocytes, whereas several inflammation-related pathways were downregulated. Low CCR5/CCL5 expression was associated with increased benefit from cetuximab-FOLFOX treatment in the CALGB/SWOG 80405 trial, where significant treatment interaction was observed with biologic agents and chemotherapy backbone. CONCLUSIONS: Our data show a strong association between CCR5/CCL5 gene expression and distinct molecular features, gene expression profiles, TME cell infiltration, and treatment benefit in CRC. Targeting the CCR5/CCL5 axis may have clinical applications in selected CRC subgroups and may play a key role in developing and deploying strategies to modulate the immune TME for CRC treatment.

Germline cis variant determines epigenetic regulation of the anti-cancer drug metabolism gene dihydropyrimidine dehydrogenase (DPYD).

Enhancers are critical for regulating tissue-specific gene expression, and genetic variants within enhancer regions have been suggested to contribute to various cancer-related processes, including therapeutic resistance. However, the precise mechanisms remain elusive. Using a well-defined drug-gene pair, we identified an enhancer region for dihydropyrimidine dehydrogenase (DPD, DPYD gene) expression that is relevant to the metabolism of the anti-cancer drug 5-fluorouracil (5-FU). Using reporter systems, CRISPR genome edited cell models, and human liver specimens, we demonstrated in vitro and vivo that genotype status for the common germline variant (rs4294451; 27% global minor allele frequency) located within this novel enhancer controls DPYD transcription and alters resistance to 5-FU. The variant genotype increases recruitment of the transcription factor CEBPB to the enhancer and alters the level of direct interactions between the enhancer and DPYD promoter. Our data provide insight into the regulatory mechanisms controlling sensitivity and resistance to 5-FU.

DNA Mutational Profiling in Patients With Colorectal Cancer Treated With Standard of Care Reveals Differences in Outcome and Racial Distribution of Mutations.

PURPOSE: CALGB (Alliance)/SWOG 80405 was a randomized phase III trial that in first-line patients with metastatic colorectal cancer (mCRC) treated with bevacizumab or cetuximab with chemotherapy. We aimed to discover novel mutated genes associated with prognosis and differential response to therapy with the biologics. METHODS: Primary tumor DNA from 548 patients was sequenced using FoundationOne. The effect of mutated genes and mutations on overall survival (OS) was tested adjusting for microsatellite instability status, BRAF V600E, all RAS mutations, arm, sex, and age. RESULTS: The median number (lower-upper quartile) of mutated genes was 5 (3-7), 5 (3-6) in microsatellite stable and 12.5 (4.5-32) in microsatellite instability-high tumors. Mutated KRAS and APC were more frequent in Black (53% and 85%) than White (27% and 65%, respectively) patients while BRAF V600E was less frequent in Black (5%) than White (14%) patients. The median OS in patients with BRAF non-V600E (2.2% of patients) was 31.9 months (95% CI, 15.1 to not applicable [NA]) similar to that of BRAF wild-type (WT) patients (31.2 months [95% CI, 29.0 to 33.9]). Mutated LRP1B (10.7% of patients) was associated with improved OS compared with WT LRP1B (hazard ratio, 0.57 [95% CI, 0.40 to 0.80]). RNF43 (5.6% of patients) interacted with treatment arms as, in the cetuximab arm, patients with mutated RNF43 had a median OS of 11.5 (95% CI, 10.8 to NA) months compared with 30.1 (95% CI, 24.9 to 35.3) months in patients with WT RNF43, whereas in the bevacizumab arm, patients with mutated RNF43 had a median OS of 25.0 (95% CI, 14.2 to NA) months compared with 31.3 (95% CI, 29.0 to 34.3) months in patients with WT RNF43. CONCLUSION: These results can provide new tools to predict patient outcome and improve therapeutic decisions and trial participation in patient minorities. The molecular alterations identified in this study may direct biomarker-driven studies.

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.

Sequencing of genes of drug response in tumor DNA and implications for precision medicine in cancer patients.

Tumor DNA sequencing is becoming standard-of-care for patient treatment decisions. We evaluated genotype concordance between tumor DNA and genomic DNA from blood and catalogued functional effects of somatic mutations in 21 drug response genes in 752 solid tumor patients. Using a threshold of 10% difference between tumor and blood DNA variant allele fraction (VAF), concordance for heterogenous genotype calls was 78% and increased to 97.5% using a 30% VAF threshold. Somatic mutations were observed in all 21 drug response genes, and 44% of patients had at least one somatic mutation in these genes. In tumor DNA, eight patients had a frameshift mutation in CYP2C8, which metabolizes taxanes. Overall, somatic copy number losses were more frequent than gains, including for CYP2C19 and CYP2D6 which had the most frequent copy number losses. However, copy number gains in TPMT were more than four times as common as losses. Seven % of patients had copy number gains in ABCB1, a multidrug resistance transporter of anti-cancer agents. These results demonstrate tumor-only DNA sequencing might not be reliable to call germline genotypes of drug response variants.

Identification of Reduced ERAP2 Expression and a Novel HLA Allele as Components of a Risk Score for Susceptibility to Liver Injury Due to Amoxicillin-Clavulanate.

BACKGROUND & AIMS: Drug-induced liver injury (DILI) due to amoxicillin-clavulanate (AC) has been associated with HLA-A∗02:01, HLA-DRB1∗15:01, and rs2476601, a missense variant in PTPN22. The aim of this study was to identify novel risk factors for AC-DILI and to construct a genetic risk score (GRS). METHODS: Transcriptome-wide association study and genome-wide association study analyses were performed on 444 AC-DILI cases and 10,397 population-based controls of European descent. Associations were confirmed in a validation cohort (n = 133 cases and 17,836 population-based controls). Discovery and validation AC-DILI cases were also compared with 1358 and 403 non-AC-DILI cases. RESULTS: Transcriptome-wide association study revealed a significant association of AC-DILI risk with reduced liver expression of ERAP2 (P = 3.7 × 10-7), coding for an aminopeptidase involved in antigen presentation. The lead eQTL single nucleotide polymorphism, rs1363907 (G), was associated with AC-DILI risk in the discovery (odds ratio [OR], 1.68; 95% CI, 1.23-1.66; P = 1.7 × 10-7) and validation cohorts (OR, 1.2; 95% CI, 1.04-2.05; P = .03), following a recessive model. We also identified HLA-B∗15:18 as a novel AC-DILI risk factor in both discovery (OR, 4.19; 95% CI, 2.09-8.36; P = 4.9 × 10-5) and validation (OR, 7.78; 95% CI, 2.75-21.99; P = .0001) cohorts. GRS, incorporating rs1363907, rs2476601, HLA-B∗15:18, HLA-A∗02:01, and HLA-DRB1∗15:01, was highly predictive of AC-DILI risk when cases were analyzed against both general population and non-AC-DILI control cohorts. GRS was the most significant predictor in a regression model containing known AC-DILI clinical risk characteristics and significantly improved the predictive model. CONCLUSIONS: We identified novel associations of AC-DILI risk with ERAP2 low expression and with HLA-B∗15:18. GRS based on the 5 risk variants may assist AC-DILI causality assessment and risk management.

Acquired Genomic Alterations on First-Line Chemotherapy With Cetuximab in Advanced Colorectal Cancer: Circulating Tumor DNA Analysis of the CALGB/SWOG-80405 Trial (Alliance).

Clinical trials frequently include multiple end points that mature at different times. The initial report, typically based on the primary end point, may be published when key planned co-primary or secondary analyses are not yet available. Clinical Trial Updates provide an opportunity to disseminate additional results from studies, published in JCO or elsewhere, for which the primary end point has already been reported.Acquired genomic alterations (Acq-GAs), specifically RAS, BRAF, and EGFR-ectodomain mutations and ERBB2 and MET amplifications, are recognized as major mechanisms of resistance to later-line anti-EGFR-antibody therapy in metastatic colorectal cancer (mCRC). However, data regarding emergence of these Acq-GAs under the selective pressure of first-line anti-EGFR-chemotherapy are lacking. We performed next-generation sequencing (Guardant360) on circulating tumor DNA obtained from paired plasma samples (pretreatment and postprogression) from the CALGB/SWOG-80405 trial, which randomly assigned patients with mCRC between first-line chemotherapy with cetuximab (anti-EGFR-chemotherapy) or bevacizumab (anti-VEGF-chemotherapy). The primary objective was to determine the prevalence of Acq-GAs on anti-EGFR-chemotherapy and compare this to the prevalence with anti-VEGF-chemotherapy on trial and pooled estimates (N = 292) seen with later-line anti-EGFR-antibody therapy as reported in the literature. Among the 61 patients on anti-EGFR-chemotherapy, only four (6.6%) developed ≥ 1 Acq-GAs of interest compared with 10.1% (7) on anti-VEGF-chemotherapy (odds ratio, 0.62; 95% CI, 0.20 to 2.11) and 62.0% on anti-EGFR-antibody therapy in later lines (odds ratio, 0.09; 95% CI, 0.03 to 0.23). Acq-GAs, classically associated with anti-EGFR-antibody resistance in later lines (RAS, BRAF, and EGFR-ectodomain mutations; ERBB2 and MET amplifications), were rare with up-front use of anti-EGFR-chemotherapy indicating divergent resistance mechanisms. These findings have critical translational relevance to timing and value of circulating tumor DNA-guided anti-EGFR rechallenge in patients with mCRC, especially those treated with anti-EGFR therapy upfront.

Gene signatures derived from transcriptomic-causal networks stratified colorectal cancer patients for effective targeted therapy.

Predictive and prognostic gene signatures derived from interconnectivity among genes can tailor clinical care to patients in cancer treatment. We identified gene interconnectivity as the transcriptomic-causal network by integrating germline genotyping and tumor RNA-seq data from 1,165 patients with metastatic colorectal cancer (CRC). The patients were enrolled in a clinical trial with randomized treatment, either cetuximab or bevacizumab in combination with chemotherapy. We linked the network to overall survival (OS) and detected novel biomarkers by controlling for confounding genes. Our data-driven approach discerned sets of genes, each set collectively stratify patients based on OS. Two signatures under the cetuximab treatment were related to wound healing and macrophages. The signature under the bevacizumab treatment was related to cytotoxicity and we replicated its effect on OS using an external cohort. We also showed that the genes influencing OS within the signatures are downregulated in CRC tumor vs. normal tissue using another external cohort. Furthermore, the corresponding proteins encoded by the genes within the signatures interact each other and are functionally related. In conclusion, this study identified a group of genes that collectively stratified patients based on OS and uncovered promising novel prognostic biomarkers for personalized treatment of CRC using transcriptomic causal networks.

A resource for integrated genomic analysis of the human liver.

In this study, we generated whole-transcriptome RNA-Seq from n = 192 genotyped liver samples and used these data with existing data from the GTEx Project (RNA-Seq) and previous liver eQTL (microarray) studies to create an enhanced transcriptomic sequence resource in the human liver. Analyses of genotype-expression associations show pronounced enrichment of associations with genes of drug response. The associations are primarily consistent across the two RNA-Seq datasets, with some modest variation, indicating the importance of obtaining multiple datasets to produce a robust resource. We further used an empirical Bayesian model to compare eQTL patterns in liver and an additional 20 GTEx tissues, finding that MHC genes, and especially class II genes, are enriched for liver-specific eQTL patterns. To illustrate the utility of the resource to augment GWAS analysis with small sample sizes, we developed a novel meta-analysis technique to combine several liver eQTL data sources. We also illustrate its application using a transcriptome-enhanced re-analysis of a study of neutropenia in pancreatic cancer patients. The associations of genotype with liver expression, including splice variation and its genetic associations, are made available in a searchable genome browser.

ARID1A mutations confer intrinsic and acquired resistance to cetuximab treatment in colorectal cancer.

Most colorectal (CRC) tumors are dependent on EGFR/KRAS/BRAF/MAPK signaling activation. ARID1A is an epigenetic regulator mutated in approximately 5% of non-hypermutated CRC tumors. Here we show that anti-EGFR but not anti-VEGF treatment enriches for emerging ARID1A mutations in CRC patients. In addition, we find that patients with ARID1A mutations, at baseline, are associated with worse outcome when treated with cetuximab- but not bevacizumab-containing therapies; thus, this suggests that ARID1A mutations may provide both an acquired and intrinsic mechanism of resistance to anti-EGFR therapies. We find that, ARID1A and EGFR-pathway genetic alterations are mutually exclusive across lung and colorectal cancers, further supporting a functional connection between these pathways. Our results not only suggest that ARID1A could be potentially used as a predictive biomarker for cetuximab treatment decisions but also provide a rationale for exploring therapeutic MAPK inhibition in an unexpected but genetically defined segment of CRC patients.

Identification and single-base gene-editing functional validation of a cis-EPO variant as a genetic predictor for EPO-increasing therapies.

Hypoxia-inducible factor prolyl hydroxylase inhibitors (HIF-PHIs) are currently under clinical development for treating anemia in chronic kidney disease (CKD), but it is important to monitor their cardiovascular safety. Genetic variants can be used as predictors to help inform the potential risk of adverse effects associated with drug treatments. We therefore aimed to use human genetics to help assess the risk of adverse cardiovascular events associated with therapeutically altered EPO levels to help inform clinical trials studying the safety of HIF-PHIs. By performing a genome-wide association meta-analysis of EPO (n = 6,127), we identified a cis-EPO variant (rs1617640) lying in the EPO promoter region. We validated this variant as most likely causal in controlling EPO levels by using genetic and functional approaches, including single-base gene editing. Using this variant as a partial predictor for therapeutic modulation of EPO and large genome-wide association data in Mendelian randomization tests, we found no evidence (at p < 0.05) that genetically predicted long-term rises in endogenous EPO, equivalent to a 2.2-unit increase, increased risk of coronary artery disease (CAD, OR [95% CI] = 1.01 [0.93, 1.07]), myocardial infarction (MI, OR [95% CI] = 0.99 [0.87, 1.15]), or stroke (OR [95% CI] = 0.97 [0.87, 1.07]). We could exclude increased odds of 1.15 for cardiovascular disease for a 2.2-unit EPO increase. A combination of genetic and functional studies provides a powerful approach to investigate the potential therapeutic profile of EPO-increasing therapies for treating anemia in CKD.

Polygenic Risk Scores for Blood Pressure to Assess the Risk of Severe Bevacizumab-Induced Hypertension in Cancer Patients (Alliance).

Hypertension is a common bevacizumab-induced toxicity. No markers are available to predict patients at risk of developing hypertension. We hypothesized that genetic risk of essential hypertension, as measured by a blood pressure polygenic risk score (PRS), would be associated with risk of severe bevacizumab-induced hypertension. PRSs were calculated for 1,027 bevacizumab-treated patients of European descent with cancer from four clinical trials (Alliance for Clinical Trials in Oncology (Alliance) / Cancer and Leukemia Group B (CALGB) 80303, 40503, 90401, 40502) using summary systolic blood pressure (SBP) and diastolic blood pressure (DBP) genome-wide association results obtained from 757,601 individuals of European descent. The association between PRS and grade 3 bevacizumab-induced hypertension (Common Toxicity Criteria for Adverse Events version 3) in each trial was performed by multivariable logistic regression. Fixed-effect meta-analyses odds ratios (ORs) per standard deviation (SD) of the association of PRS (quantitative) and hypertension across trials were estimated by inverse-variance weighting. PRSs were additionally stratified into quintiles, with the bottom quintile as the referent group. The OR of the association between hypertension and each quintile vs. the referent group was determined by logistic regression. The most significant PRS (quantitative)-hypertension association included up to 67 single-nucleotide variants (SNPs) associated with SBP (P = 0.0077, OR per SD = 1.31, 95% confidence interval (CI), 1.07-1.60), and up to 53 SNPs associated with DBP (P = 0.0209, OR per SD = 1.27, 95% CI, 1.04-1.56). Patients in the top quintile had a higher risk of developing bevacizumab-induced hypertension compared with patients in the bottom quintile using SNPs associated with SBP (P = 4.75 × 10-4 , OR = 3.72, 95% CI, 1.84-8.16) and DBP (P = 0.076, OR = 1.83, 95% CI, 0.95-3.64). Genetic variants associated with essential hypertension, mainly SBP, increase the risk of severe bevacizumab-induced hypertension.