Computational pharmacogenotype extraction from clinical next-generation sequencing.

Background: Next-generation sequencing (NGS), including whole genome sequencing (WGS) and whole exome sequencing (WES), is increasingly being used for clinic care. While NGS data have the potential to be repurposed to support clinical pharmacogenomics (PGx), current computational approaches have not been widely validated using clinical data. In this study, we assessed the accuracy of the Aldy computational method to extract PGx genotypes from WGS and WES data for 14 and 13 major pharmacogenes, respectively. Methods: Germline DNA was isolated from whole blood samples collected for 264 patients seen at our institutional molecular solid tumor board. DNA was used for panel-based genotyping within our institutional Clinical Laboratory Improvement Amendments- (CLIA-) certified PGx laboratory. DNA was also sent to other CLIA-certified commercial laboratories for clinical WGS or WES. Aldy v3.3 and v4.4 were used to extract PGx genotypes from these NGS data, and results were compared to the panel-based genotyping reference standard that contained 45 star allele-defining variants within CYP2B6, CYP2C8, CYP2C9, CYP2C19, CYP2D6, CYP3A4, CYP3A5, CYP4F2, DPYD, G6PD, NUDT15, SLCO1B1, TPMT, and VKORC1. Results: Mean WGS read depth was >30x for all variant regions except for G6PD (average read depth was 29 reads), and mean WES read depth was >30x for all variant regions. For 94 patients with WGS, Aldy v3.3 diplotype calls were concordant with those from the genotyping reference standard in 99.5% of cases when excluding diplotypes with additional major star alleles not tested by targeted genotyping, ambiguous phasing, and CYP2D6 hybrid alleles. Aldy v3.3 identified 15 additional clinically actionable star alleles not covered by genotyping within CYP2B6, CYP2C19, DPYD, SLCO1B1, and NUDT15. Within the WGS cohort, Aldy v4.4 diplotype calls were concordant with those from genotyping in 99.7% of cases. When excluding patients with CYP2D6 copy number variation, all Aldy v4.4 diplotype calls except for one CYP3A4 diplotype call were concordant with genotyping for 161 patients in the WES cohort. Conclusion: Aldy v3.3 and v4.4 called diplotypes for major pharmacogenes from clinical WES and WGS data with >99% accuracy. These findings support the use of Aldy to repurpose clinical NGS data to inform clinical PGx.

Initial Phase I Safety Study of Gedatolisib plus Cofetuzumab Pelidotin for Patients with Metastatic Triple-Negative Breast Cancer.

PURPOSE: The PI3K pathway is dysregulated in the majority of triple-negative breast cancers (TNBC), yet single-agent inhibition of PI3K has been ineffective in TNBC. PI3K inhibition leads to an immediate compensatory upregulation of the Wnt pathway. Dual targeting of both pathways is highly synergistic against TNBC models in vitro and in vivo. We initiated a phase I clinical trial combining gedatolisib, a pan-class I isoform PI3K/mTOR inhibitor, and cofetuzumab pelidotin, an antibody-drug conjugate against the cell-surface PTK7 protein (Wnt pathway coreceptor) with an auristatin payload. PATIENTS AND METHODS: Participants (pt) had metastatic TNBC or estrogen receptor (ER) low (ER and PgR < 5%, HER2-negative) breast cancer, and had received at least one prior chemotherapy for advanced disease. The primary objective was safety. Secondary endpoints included overall response rate (ORR), clinical benefit at 18 weeks (CB18), progression-free survival (PFS), and correlative analyses. RESULTS: A total of 18 pts were enrolled in three dose cohorts: gedatolisib 110 mg weekly + cofetuzumab pelidotin 1.4 mg/kg every 3 weeks (n = 4), 180 mg + 1.4 mg/kg (n = 3), and 180 mg + 2.8 mg/kg (n = 11). Nausea, anorexia, fatigue, and mucositis were common but rarely reached ≥grade 3 severity. Myelosuppression was uncommon. ORR was 16.7% (3/18). An additional 3 pts had stable disease (of these 2 had stable disease for >18 weeks); CB18 was 27.8%. Median PFS was 2.0 months (95% confidence interval for PFS: 1.2-6.2). Pts with clinical benefit were enriched with genomic alterations in the PI3K and PTK7 pathways. CONCLUSIONS: The combination of gedatolisib + cofetuzumab pelidotin was well tolerated and demonstrated promising clinical activity. Further investigation of this drug combination in metastatic TNBC is warranted.

Analytical Validation of a Computational Method for Pharmacogenetic Genotyping from Clinical Whole Exome Sequencing.

Germline whole exome sequencing from molecular tumor boards has the potential to be repurposed to support clinical pharmacogenomics. However, accurately calling pharmacogenomics-relevant genotypes from exome sequencing data remains challenging. Accordingly, this study assessed the analytical validity of the computational tool, Aldy, in calling pharmacogenomics-relevant genotypes from exome sequencing data for 13 major pharmacogenes. Germline DNA from whole blood was obtained for 164 subjects seen at an institutional molecular solid tumor board. All subjects had whole exome sequencing from Ashion Analytics and panel-based genotyping from an institutional pharmacogenomics laboratory. Aldy version 3.3 was operationalized on the LifeOmic Precision Health Cloud with copy number fixed to two copies per gene. Aldy results were compared with those from genotyping for 56 star allele-defining variants within CYP2B6, CYP2C8, CYP2C9, CYP2C19, CYP2D6, CYP3A4, CYP3A5, CYP4F2, DPYD, G6PD, NUDT15, SLCO1B1, and TPMT. Read depth was >100× for all variants except CYP3A4∗22. For 75 subjects in the validation cohort, all 3393 Aldy variant calls were concordant with genotyping. Aldy calls for 736 diplotypes containing alleles assessed by both platforms were also concordant. Aldy identified additional star alleles not covered by targeted genotyping for 139 diplotypes. Aldy accurately called variants and diplotypes for 13 major pharmacogenes, except for CYP2D6 variants involving copy number variations, thus allowing repurposing of whole exome sequencing to support clinical pharmacogenomics.

Potential Utility of Pre-Emptive Germline Pharmacogenetics in Breast Cancer.

Patients with breast cancer often receive many drugs to manage the cancer, side effects associated with cancer treatment, and co-morbidities (i.e., polypharmacy). Drug-drug and drug-gene interactions contribute to the risk of adverse events (AEs), which could lead to non-adherence and reduced efficacy. Here we investigated several well-characterized inherited (germline) pharmacogenetic (PGx) targets in 225 patients with breast cancer. All relevant clinical, pharmaceutical, and PGx diplotype data were aggregated into a single unifying informatics platform to enable an exploratory analysis of the cohort and to evaluate pharmacy ordering patterns. Of the drugs recorded, there were 38 for which high levels of evidence for clinical actionability with PGx was available from the US FDA and/or the Clinical Pharmacogenetics Implementation Consortium (CPIC). These data were associated with 10 pharmacogenes: DPYD, CYP2C9, CYP2C19, CYP2D6, CYP3A5, CYP4F2, G6PD, MT-RNR1, SLCO1B1, and VKORC1. All patients were taking at least one of the 38 drugs and had inherited at least one actionable PGx variant that would have informed prescribing decisions if this information had been available pre-emptively. The non-cancer drugs with PGx implications that were common (prescribed to at least one-third of patients) included anti-depressants, anti-infectives, non-steroidal anti-inflammatory drugs, opioids, and proton pump inhibitors. Based on these results, we conclude that pre-emptive PGx testing may benefit patients with breast cancer by informing drug and dose selection to maximize efficacy and minimize AEs.

Molecular dissection of CRC primary tumors and their matched liver metastases reveals critical role of immune microenvironment, EMT and angiogenesis in cancer metastasis.

Metastasis is the primary cause of cancer mortality. The primary tumors of colorectal cancer (CRC) often metastasize to the liver. In this study, we have collected 122 samples from 45 CRC patients. Among them, 32 patients have primary tumors, adjacent normal tissues, and matched liver metastases. Thirteen patients have primary tumors without distant metastasis and matched normal tissues. Characterization of these samples was conducted by whole-exome and RNA sequencing and SNP6.0 analysis. Our results revealed no significant difference in genetic alterations including common oncogenic mutations, whole genome mutations and copy number variations between primary and metastatic tumors. We then assembled gene co-expression networks and identified metastasis-correlated gene networks of immune-suppression, epithelial-mesenchymal transition (EMT) and angiogenesis as the key events and potentially synergistic drivers associated with CRC metastasis. Further independent cohort validation using published datasets has verified that these specific gene networks are up regulated throughout the tumor progression. The gene networks of EMT, angiogenesis, immune-suppression and T cell exhaustion are closely correlated with the poor patient outcome and intrinsic anti-PD-1 resistance. These results offer insights of combinational strategy for the treatment of metastatic CRC.

Genomic characterization of intrinsic and acquired resistance to cetuximab in colorectal cancer patients.

Anti-EGFR antibodies are effective in therapies for late-stage colorectal cancer (CRC); however, many tumours are unresponsive or develop resistance. We performed genomic analysis of intrinsic and acquired resistance to anti-EGFR therapy in prospectively collected tumour samples from 25 CRC patients receiving cetuximab (an EGFR inhibitor). Of 25 CRC patients, 13 displayed intrinsic resistance to cetuximab; 12 were intrinsically sensitive. We obtained six re-biopsy samples at acquired resistance from the intrinsically sensitive patients. NCOA4-RET and LMNA-NTRK1 fusions and NRG1 and GNAS amplifications were found in intrinsic-resistant patients. In cetuximab-sensitive patients, we found KRAS K117N and A146T mutations in addition to BRAF V600E, AKT1 E17K, PIK3CA E542K, and FGFR1 or ERBB2 amplifications. The comparison between baseline and acquired-resistant tumours revealed an extreme shift in variant allele frequency of somatic variants, suggesting that cetuximab exposure dramatically selected for rare resistant subclones that were initially undetectable. There was also an increase in epithelial-to-mesenchymal transition at acquired resistance, with a reduction in the immune infiltrate. Furthermore, characterization of an acquired-resistant, patient-derived cell line showed that PI3K/mTOR inhibition could rescue cetuximab resistance. Thus, we uncovered novel genomic alterations that elucidate the mechanisms of sensitivity and resistance to anti-EGFR therapy in metastatic CRC patients.

Genomic Aberrations that Activate D-type Cyclins Are Associated with Enhanced Sensitivity to the CDK4 and CDK6 Inhibitor Abemaciclib.

Most cancers preserve functional retinoblastoma (Rb) and may, therefore, respond to inhibition of D-cyclin-dependent Rb kinases, CDK4 and CDK6. To date, CDK4/6 inhibitors have shown promising clinical activity in breast cancer and lymphomas, but it is not clear which additional Rb-positive cancers might benefit from these agents. No systematic survey to compare relative sensitivities across tumor types and define molecular determinants of response has been described. We report a subset of cancers highly sensitive to CDK4/6 inhibition and characterized by various genomic aberrations known to elevate D-cyclin levels and describe a recurrent CCND1 3'UTR mutation associated with increased expression in endometrial cancer. The results suggest multiple additional classes of cancer that may benefit from CDK4/6-inhibiting drugs such as abemaciclib.

New perspectives on the biology of fragile X syndrome.

Fragile X syndrome (FXS) is a trinucleotide repeat disorder caused by a CGG repeat expansion in FMR1, and loss of its protein product FMRP. Recent studies have provided increased support for the role of FMRP in translational repression via ribosomal stalling and the microRNA pathway. In neurons, particular focus has been placed on identifying the signaling pathways such as PI3K and mTOR downstream of group 1 metabotropic glutamate receptors (mGluR1/5) that regulate FMRP. New evidence also suggests that loss of FMRP causes presynaptic dysfunction and abnormal adult neurogenesis. In addition, studies on FXS stem cells especially induced pluripotent stem (iPS) cells and new sequencing efforts hold out promise for deeper understanding of the silencing process and mutation spectrum of FMR1.

Molecular mechanisms of fragile X syndrome: a twenty-year perspective.

Fragile X syndrome (FXS) is a common form of inherited intellectual disability and is one of the leading known causes of autism. The mutation responsible for FXS is a large expansion of the trinucleotide CGG repeat in the 5' untranslated region of the X-linked gene FMR1. This expansion leads to DNA methylation of FMR1 and to transcriptional silencing, which results in the absence of the gene product, FMRP, a selective messenger RNA (mRNA)-binding protein that regulates the translation of a subset of dendritic mRNAs. FMRP is critical for mGluR (metabotropic glutamate receptor)-dependent long-term depression, as well as for other forms of synaptic plasticity; its absence causes excessive and persistent protein synthesis in postsynaptic dendrites and dysregulated synaptic function. Studies continue to refine our understanding of FMRP's role in synaptic plasticity and to uncover new functions of this protein, which have illuminated therapeutic approaches for FXS.

Signatures of founder effects, admixture, and selection in the Ashkenazi Jewish population.

The Ashkenazi Jewish (AJ) population has long been viewed as a genetic isolate, yet it is still unclear how population bottlenecks, admixture, or positive selection contribute to its genetic structure. Here we analyzed a large AJ cohort and found higher linkage disequilibrium (LD) and identity-by-descent relative to Europeans, as expected for an isolate. However, paradoxically we also found higher genetic diversity, a sign of an older or more admixed population but not of a long-term isolate. Recent reports have reaffirmed that the AJ population has a common Middle Eastern origin with other Jewish Diaspora populations, but also suggest that the AJ population, compared with other Jews, has had the most European admixture. Our analysis indeed revealed higher European admixture than predicted from previous Y-chromosome analyses. Moreover, we also show that admixture directly correlates with high LD, suggesting that admixture has increased both genetic diversity and LD in the AJ population. Additionally, we applied extended haplotype tests to determine whether positive selection can account for the level of AJ-prevalent diseases. We identified genomic regions under selection that account for lactose and alcohol tolerance, and although we found evidence for positive selection at some AJ-prevalent disease loci, the higher incidence of the majority of these diseases is likely the result of genetic drift following a bottleneck. Thus, the AJ population shows evidence of past founding events; however, admixture and selection have also strongly influenced its current genetic makeup.

Identification of novel FMR1 variants by massively parallel sequencing in developmentally delayed males.

Fragile X syndrome (FXS), the most common inherited form of developmental delay, is typically caused by CGG-repeat expansion in FMR1. However, little attention has been paid to sequence variants in FMR1. Through the use of pooled-template massively parallel sequencing, we identified 130 novel FMR1 sequence variants in a population of 963 developmentally delayed males without CGG-repeat expansion mutations. Among these, we identified a novel missense change, p.R138Q, which alters a conserved residue in the nuclear localization signal of FMRP. We have also identified three promoter mutations in this population, all of which significantly reduce in vitro levels of FMR1 transcription. Additionally, we identified 10 noncoding variants of possible functional significance in the introns and 3'-untranslated region of FMR1, including two predicted splice site mutations. These findings greatly expand the catalog of known FMR1 sequence variants and suggest that FMR1 sequence variants may represent an important cause of developmental delay.

Identification of small molecules rescuing fragile X syndrome phenotypes in Drosophila.

Fragile X syndrome is caused by the functional loss of the fragile X mental retardation 1 (FMR1) gene. Deletion of the FMR1 ortholog in Drosophila melanogaster (Fmr1) recapitulates many phenotypes associated with fragile X syndrome. We have discovered that Fmr1 mutant Drosophila die during development when reared on food containing increased levels of glutamate, which is consistent with the theory that FMR1 loss results in excess glutamate signaling. Using this lethal phenotype, we screened a chemical library of 2,000 compounds and identified nine molecules that rescued the lethality, including three that implicate the GABAergic inhibitory pathway. Indeed, GABA treatment rescued several known Fmr1 mutant phenotypes in flies, including mushroom bodies defects, excess Futsch translation and abnormal male courtship behavior. These data are consistent with GABAergic inhibition of the enhanced excitatory pathway in fragile X syndrome. In addition, our screen reveals that the muscarinic cholinergic receptors may have a role in fragile X syndrome in parallel to the GABAergic pathway. These results point to potential therapeutic approaches for treating fragile X syndrome.

Regulation of angiotensin II-induced G protein signaling by phosducin-like protein.

Phosducin-like protein (PhLP) is a broadly expressed member of the phosducin (Pd) family of G protein betagamma subunit (Gbetagamma)-binding proteins. Though PhLP has been shown to bind Gbetagamma in vitro, little is known about its physiological function. In the present study, the effect of PhLP on angiotensin II (Ang II) signaling was measured in Chinese hamster ovary cells expressing the type 1 Ang II receptor and various amounts of PhLP. Up to 3.6-fold overexpression of PhLP had no effect on Ang II-stimulated inositol trisphosphate (IP(3)) formation, whereas further increases caused an abrupt decrease in IP(3) production with half-maximal inhibition occurring at 6-fold PhLP overexpression. This threshold level for inhibition corresponds to the cellular concentration of cytosolic chaperonin complex, a recently described binding partner that preferentially binds PhLP over Gbetagamma. Results of pertussis toxin sensitivity, GTPgammaS binding, and immunoprecipitation experiments suggest that PhLP inhibits phospholipase Cbeta activation by dual mechanisms: (i) steric blockage of Gbetagamma activation of PLCbeta and (ii) interference with Gbetagamma-dependent cycling of G(q)alpha by the receptor. These results suggest that G protein signaling may be regulated through controlling the cellular concentration of free PhLP by inducing its expression or by regulating its binding to the chaperonin.