1Department of Biomedical Sciences - BIOMORF,
2Unit of Medical Genetics,
Correspondence Address: Prof. Saverio Alberti, Unit of Medical Genetics, BIOMORF Department of Biomedical Sciences, University of Messina, via Consolare Valeria, 98125 Messina, Italy. E-mail: firstname.lastname@example.org
© The Author(s) 2019. Open Access This article is licensed under a Creative Commons Attribution 4.0 International License (https://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, sharing, adaptation, distribution and reproduction in any medium or format, for any purpose, even commercially, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.
Genetic/genomic profiling at a single-patient level is expected to provide critical information for determining inter-individual drug toxicity and potential efficacy in cancer therapy. A better definition of cancer subtypes at a molecular level, may correspondingly complement such pharmacogenetic and pharmacogenomic approaches, for more effective personalized treatments. Current pharmacogenetic/pharmacogenomic strategies are largely based on the identification of known polymorphisms, thus limiting the discovery of novel or rarer genetic variants. Recent improvements in cost and throughput of next generation sequencing (NGS) are now making whole-genome profiling a plausible alternative for clinical procedures. Beyond classical pharmacogenetic/pharmacogenomic traits for drug metabolism, NGS screening programs of cancer genomes may lead to the identification of novel cancer-driving mutations. These may not only constitute novel therapeutic targets, but also effector determinants for metabolic pathways linked to drug metabolism. An additional advantage is that cancer NGS profiling is now leading to discovering targetable mutations, e.g., in glioblastomas and pancreatic cancers, which were originally discovered in other tumor types, thus allowing for effective repurposing of active drugs already on the market.
Pharmacogenetics, pharmacogenomics, cancer, next-generation sequencing, genomic variants
Advances in cancer treatment over the last decades have more and more relied on identifying therapy-responsive subgroups on the basis of evolving tumor subgroup classifications, as based on pathological and clinical parameters, tumor stage, cancer biomarkers[4,5] and cancer drivers[6-14], including both somatic and hereditary epigenetic changes[15,16]. However, critical limitations remain in our ability to predict therapeutic efficacy and toxicity for tumor subgroups which appear homogeneous by all the above criteria, such target remaining an unmet medical need.
Over the past several years, there has been a major shift in cancer diagnostics from physical examination/in vivo imaging/histopathological analysis to assessment of tumor biomarkers, cancer drivers and targetable genomic mutations. Genetic variation in drug response is well documented and influences treatment efficacy and/or toxicity. Numerous genomic variants are increasingly recognized as both potential therapeutic targets and drug metabolism modifiers[18,19], thus gaining recognition by regulatory agencies [Figure 1, Supplementary Table 1]. The concept of precision/personalized treatment is more and more supported by these complementary strategies, so that clinical decisions are adjusted on the basis of each patient’s genetic background, as complemented by detailed information on the patient’s tumor[20,21].
Figure 1. Pharmacogenetic biomarkers in FDA drug labels over the last decade (www.fda.gov/BiologicsBloodVaccines/DevelopmentApprovalProcess/BiologicalApprovalsbyYear/ucm596371.htm) has undergone an exponential increase. It should be noted that this data-set covers both somatic mutations and hereditary variants, together with genes that indirectly affect pharmacokinetics through drug-drug interactions
Pharmacogenetics and pharmacogenomics constitute a branch of molecular medicine focused on the study of human genome in the context of therapeutic decision-making. Pharmacogenetics is the study of inherited variations in the DNA sequence of known genes, e.g., involved in drug metabolic pathways, which can affect individual responses to drugs, both in terms of therapeutic responses and adverse effects. Pharmacogenomics is the study of acquired and inherited genetic differences in relation to drug response and drug behavior, through a systematic examination of DNA sequence at the level of the entire genome.
Germline mutations may help to predict alterations in drug-metabolism-related genes (e.g., CYP450) and pathways. This can aid in correspondingly adjusting for systemic drug exposure, with the aim of minimizing drug toxicity and improving disease outcome[17,25,26]. Drug’ efficacy and toxicity rely on complex interactions between different gene-encoded metabolic pathways and environmental factors. Highly, complementary approaches are required to dissect such complexity [Table 1].
Discovery strategies for novel pharmacogenetic and pharmacogenomic traits
|Candidate Polymorphism Analysis|
|1. Rapid execution of the assay||1. Polymorphisms need to have strong effects toward the phenotype|
|2. Focus on genes likely involved in treatment response and toxicity||2. It is based on validated knowledge, ad may miss potentially-involved unknown genes|
|3. It ignores de-novo mutations in target genes|
|1. Focus on pathways, downstream the gene(s) of interest, that are highly likely to be involved in the drug action||1. It may miss potentially involved, but still unidentified, signaling cascades|
|2. It highlights whole signaling cascades, for higher sensitivity for genes with smaller phenotypic effects||2. Data analysis is complex given the interplay of multiple interacting genes|
|3. It can identify new polymorphisms or new genes within a given pathway||3. It requires investigating large sample case-series|
|4. More likely to explain inter-individual variation in drug response|
|“Whole genome strategies”|
|1. They provide a complete gene- or protein- expression profile (tumor or individual)||1. The lack of hypothesis-driven analyses may increase the risk of false positives|
|2. They provide information on novel associations||2. Complex data management and analysis procedures are required|
|3. They generate large amounts of data||3. Costs and complexity still high for the clinics|
|4. Useful in predicting tumor response|
Cancer-causing DNA alterations, such as somatic DNA mutations and inherited DNA variants, are not a direct focus of pharmacogenetic and pharmacogenomic studies. However, mutated cancer drivers are becoming more and more actionable targets for therapy. They can also affect key metabolic pathways that may modify drug pharmacokinetics and pharmacodynamics. As such, they can play a key role in pharmacogenetic and pharmacogenomic discovery.
Key examples are DNA-damage response pathways (ATM, CHEK2, BRIP1, BRCA1, BRCA2, PALB2), which are associated with an increase in breast cancer risk. In particular, germline pathogenic variants in BARD1, BRCA1, BRCA2, PALB2 and RAD51D are associated with high risk (odds ratio > 5.0) for triple-negative breast cancer (TNBC)[28,29]. BRCA1, BRCA2, PALB2 variants are also associated with increased risk of ovarian and other cancers[31,32]. Pathogenic variants in BRIP1, RAD51C and TP53 are associated with moderate risk (odds ratio > 2) for TNBC, whereby hereditary pathogenic variants are detected in 12.0% of TNBC.
Corresponding mutations in the BRCA1, BRCA2, PALB2 genes are detected in a fraction of sporadic breast and ovarian cancer. Notably, mutations in BRCA1, BRCA2 are associated to better therapeutic response to PARP inhibitors, making such mutations effective tools for therapy choice. Other, cancer-associated somatic mutations, e.g., the mutations at codons 12 or 13 in the Ki-RAS gene, predict lack of response to EGFR-targeted therapy in colon cancer. Mutations in the EGFR kinase ATP-binding pocket mandate choosing specific tyrosine kinase inhibitors. Of interest, such EGFR mutants also are predictors of response to unrelated therapy.
It should be noted that mutated cancer-driving genes may play a key role as metabolic-pathway modifiers. Examples are mutations in TP53, which differentially impact on metabolic pathways and apoptotic responses, thus modifying the impact of anticancer chemotherapy. Not surprisingly, TP53 mutations are predictive of poor response to therapy. Activation of transcription of c-Myc impacts on main metabolic pathways, on ribosomal biogenesis, and on lipid metabolism. In vivo efficacy of the Bcl-2 antagonist ABT-737 depends on c-Myc activation. Furthermore, c-myc/p53 interactions determine sensitivity of colon cancer to 5-Fluorouracil. Mutations of the PIK3CA gene are predictors of response to AKT and mTOR inhibitors. Hence, taking into account the above caveats, cancer drivers can constitute rather attractive targets for pharmacogenetic investigation.
As whole-genome information is rapidly accumulating [Figure 1], it is likely that more and more genetic factors, or clusters of them, will be discovered that may affect both tumor progression, drug pharmacodynamics and overall response to therapy. Such cross-feeding amid convergent research fields is expected to foster better knowledge about gene-gene interactions vs. therapeutic drug metabolism.
This type of analysis seeks polymorphic DNA sequences within specific genes, known to impact on the pharmacokinetics or pharmacodynamics of a compound. Such information may aid in the selection among different therapeutic strategies. When sufficiently large clinical data are available, they may also aid in drug dosage selection.
A potential functional impact of genes mapping near polymorphic sites can be explored through gene silencing. Gene silencing can be obtained through RNA interference via short hairpin RNA (shRNA; preceding methods relied on siRNA)[44-46]. shRNA have been proposed for use as pharmacological compounds. More recent applications include Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR) technology, which can be utilized to ablate target transcription factors, chromatin-modifying factors, and noncoding RNA, thus providing powerful instruments for gene silencing purposes[48-51].
This strategy covers a wider range than the previous one. It still retains some of the limitation of the candidate-gene approach, namely the need of previous knowledge regarding the biological cascades correlated to a compound action. However, the pathway approach extends such an analysis to larger numbers of related genes, whose altered function may impact on treatment efficacy. This approach may add significant chances to explain inter-individual variation in drug efficacy/toxicity, as it includes evaluation of potential epistatic effects and of the influence of other cis-regulatory elements. These studies may also discover gene-gene interactions as pharmacological targets.
Newly added research strategies in this field are based on the concept of synthetic lethality. A synthetic lethality event occurs when two or more genes are simultaneously perturbed and exposure to a drug results in cellular or organism death/impairment. Large-scale analysis in yeast for interactions among orthologs of human tumor suppressor genes allowed to evaluate thousands of genotype-drug combinations, that were then transferred to cancer cells in vitro. This resulted in the identification of networks of conserved, synthetic lethal interactions. Among them, the interaction of topoisomerases with RAD17 and of checkpoint kinases with BLM were validated by patient survival data. These techniques inform on gene-gene and gene-drug dependencies, taking advantage of better knowledge on control gene pathways, candidate for targeted drug design. PARP inhibitor efficacy in BRCA1 and BRCA2 mutation carriers is limited by inherent and acquired resistance. By using synthetic lethality in combination with PARP inhibitors, BRCA2-mutant cells were found to be dependent on base excision repair, Atr activation and mRNA splicing. Subsequently, FEN1 and APEX2 were identified as BRCA2 synthetic lethal targets.
The previously described approaches suffer from a critical limitation, i.e., they only focus on pre-specified sets of genes (knowledge-based), thereby ignoring potential candidates whose effect has not been linked, yet, to treatment efficacy/toxicity. This limitation can be overcome by global strategies, whereby the focus of the search is extended to the entire genome, transcriptome or proteome of an individual or group of individuals. Since the completion of the human genome project and the early efforts to map human genetic variations, tools were developed to evaluate the relationship between these alterations and drug efficacy in an essentially unbiased fashion. Using high-resolution genotyping platforms, genome-wide association studies (GWAS) were rapidly conducted on nearly all common cancers. Although hundreds of statistically robust risk variants, largely in the form of single-nucleotide polymorphisms (SNP), were identified, each genetic variant was mostly associated with a modest increase in disease risk (relative risk ≤ 1.5). Moreover, as ≈90% of risk variants reside in noncoding introns, causal factors associated to most risk loci have remained elusive. Albeit a useful tool for investigation and discovery of new genes, the use of whole-genome strategies in clinical practice is restricted by high costs, labor intensity and complex interpretation of the data. Thermodynamic limitations add to the high amount of noise associated with this system.
Better knowledge on cancer drivers[7-14], on their use as therapy-target molecules[13,58] and on their interaction with drug efficacy/ toxicity determinants[9,13,14,39,41,59] has reached milestone successes[39,59]. However, current strategies remain largely based on known polymorphisms (from the human genome project), thus limiting the discovery of novel/rarer mutations [Figure 2].
Figure 2. The efficiency of discovery of novel pharmacogenetic/pharmacogenomic biomarkers, modified from. The combined expected impact of genomic variants on disease occurrence and genomic diagnostics is indicated (light green: low; deep green: high)
FDA defines genomic biomarkers as a measurable DNA and/or RNA characteristic that can be used as an indicator of either normal biologic processes, pathogenic phenomena, and/or response to therapeutic or other interventions. Frequently used genomic biomarkers include SNP and DNA segments copy number variations (CNV) (www.fda.gov/ucm/groups/fdagov-public/@fdagov-drugs-gen/documents/document/ucm073162.pdf). Applications of genomic biomarkers in pharmacogenomics and pharmacogenetics are more and more frequent (www.fda.gov/ucm/groups/fdagov-public/@fdagov-drugs-gen/documents/document/ucm073162.pdf). Currently, though, no specific recommendation for the use of genomic biomarkers’ have been established. Somewhat general guidelines for the use of biomarkers can be found at (www.fda.gov/ucm/groups/fdagov-public/@fdagov-drugs-gen/documents/document/ucm628118.pdf).
A global profiling of each individual may critically help identify determinants of predisposition to drug toxicity, as well as to drug efficacy. Improvements in cost and throughput [Figure 3] have recently surpassed competing development of leading microelectronic technologies (Moore’s law), and have made next-generation sequencing (NGS) approaches realistic and feasible [Supplementary Tables 2 and 3].
Figure 3. NGS costs per whole genome sequencing. Inset: number of whole genomes sequenced worldwide, as a compilation from publicly available genome-sequencing projects (as modified from www.genome.gov/images/content/costpergenome_2017.jpg)
Data from publicly available genome-sequencing project data indicate fitting with an exponential growth curve over time (Figure 3, inset; as modified from cdn.static economist.com/sites/default/files/images/2014/11/articles/main/20151101_stc002_595.png). Correspondingly, information on the complexity of cancer states and on tumor progression trajectories are being rapidly obtained[21,61].
Such efforts have already led to a marked increase in the number of targeted therapies approved for the treatment of patients with specific types of malignancies harboring specific types of sequence alterations[Supplementary Table 1].
The field of pharmacogenetics/pharmacogenomics has uncovered an abundance of actionable, clinically relevant markers, through distinct, complementary approaches.
Hotspot panels, examples of which are listed in Supplementary Table 4, target frequently mutated gene stretches (hotspots), which are expected to bear diagnostic/prognostic significance. Such panels have been successfully used to guide cancer treatment, e.g., for lung cancer. Most hotspot panels detect somatic mutations in tumor DNA. Previous comments on the distiction between cancer drivers/therapy targets vs. pharmacogenetic/pharmacogenomic determinants also apply here. As these panels only analyze known, pre-determined areas within specific genes, they are not amenable to the discovery of novel mutations. Even if agreement on guidelines that should be followed in clinical practice while using testing panels has yet to be reached, several consortia have been created to aid in the use of pharmacogenetic data (e.g., Clinical Pharmacogenetics Implementation Consortium)[65-67].
Similar to the hotspot ones, these panels focus on specific genes. However, they include all exons of the targeted genes, so that alterations not explicitly included in hotspot panels can also be investigated. Panels of this type have been developed for melanoma, colon and gastric cancer. While these can be considered an evolution from hotspot panels, they still do not guarantee complete coverage of the target genes (no intronic or regulatory sequences are included) and remain anchored to preexisting knowledge (only genes already known to be involved in drugs’ action are included).
Similarly to the actionable panels, disease-focused panels include sets of genes (exomes) known to be involved in a specific disease. These are mainly used as a screening panel to evaluate the hereditary risk of developing specific tumors, e.g., BRCA-1, -2-driven breast and ovarian cancers. This type of panel is more useful for screening programs and in early/preventive management of cancer.
The difficulty in recruiting adequately large sample groups of homogeneous tumors within realistic time-frames rarely justifies the costs of production of different disease focused panels. This has led life science companies to create larger, more comprehensive panels in which multiple genes that correlate with multiple diseases are included. These panels can include thousands of genes (exons only) which are selected among the ones that are most commonly tested. Like other targeted panels, comprehensive panels depend on pre-existing knowledge on the correlation between genes and the investigated phenotype.
As of November 27th, 2018, ~1677 NGS-based clinical tests for cancer were available across the world (search terms “NGS” AND “Cancer” in the NIH Genetic Testing Registry, website www.ncbi.nlm.nih.gov).
Despite the ever-increasing number of genes that are included in targeted panels, the risk of not including elements that are possibly correlated with the investigated phenotype remains high. To overcome this limitation increasing numbers of researchers resort to whole exome sequencing (WES). This approach allows to identify both known and unknown alterations that occur within coding regions. Some groups estimate that WES may cover ≈85% of disease-related variants. This estimate should be taken with great caution, as approximately 98% of the genome remains uncovered. Such regions include regulatory elements, non-coding mRNA, splice sites, anchor sites to the nuclear membrane, chromatin modifiers, and other regions we simply don’t know enough about.
Whole genome sequencing (WGS) provides a global coverage of the human genome. This approach represents the last and most comprehensive step in genomic analysis. The feasibility of this approach in recent years was aided by the introduction of higher-processivity NGS technologies (higher speed, lower costs). The main limitation of this approach is essentially related to the complexity of data post-analysis, and related statistics, with considerable risk of both false positives and false negatives.
With a coverage comparable to that of WES, RNA sequencing profiles the transcriptome of a tumor. Differently from WES, whose target is the subject’s DNA, RNA provides a dynamic/functional image of a cancer. It has the advantage of detecting the presence of fusion (onco)genes, which are often found in cancer cases[8,68,69]. It may aid clinicians in the choice of drug therapy and basic science for the identification of novel potential drug targets.
This approach is normally used to understand DNA-protein binding interactions, e.g., regulatory protein-binding sites within a DNA region. This also has the capability to detect epigenetic alterations[70,71]. The use of chromatin immunoprecipitation sequencing in cancer therapy is limited, but it has proven useful in stratifying breast cancer subtypes for treatment.
The birth of NGS sequencing can be placed in 2005, when 454 Life Sciences launched their sequencing-by-synthesis tool, as the first instrument of a 2nd generation sequencing technology. Since then, many companies have developed different NGS products, which are usually based on different, proprietary procedures and base-detection chemistry. Current best-ranking technologies, comparative performance and preferential platform applications are presented in Supplementary Tables 2 and 3.
Albeit reference technologies vary widely across NGS platforms, workflows share several steps: (1) DNA extraction; (2) library preparation, with addition of adaptors and barcodes/indexes; (3) template preparation, either by bridge amplification or emulsion PCR; (4) automated sequencing [Supplementary Table 3].
The pioneer Roche/454 (www.sequencing.roche.com) sequencing used pyrosequencing, as based on the detection of pyrophosphate released after nucleotide incorporation in newly synthesized DNA. Currently, this platform is rarely utilized, because of long run times and high costs.
The Supported Oligonucleotide Ligation and Detection (SOLiD) platform (www.thermofisher.com), is a short-read sequencing technology based on ligation. DNA fragments are ligated to oligonucleotide adapters, attached to beads, and clonally amplified by emulsion PCR. Beads with clonally amplified template are immobilized onto a derivatized-glass flow-cell surface, and sequencing is started by annealing a primer oligonucleotide complementary to the adapter. At each cycle the complementary strand is removed and a new sequencing cycle starts at the position n-1 of the template. Sequencing cycles are repeated until each base is sequenced twice. The recovered data from the color space can be translated to letters of DNA bases and the sequence of the DNA fragment is assembled.
The IonTorrent platform (www.thermofisher.com) is based on sequencing by synthesis and the detection is based on solid state pH meters, which measure the hydrogen ions that are released during DNA polymerization.
Illumina (www.emea.illumina.com) currently occupies the most prominent segment of the NGS market. This technology is based on sequencing by synthesis of template-complementary strands and on fluorescence-based detection of reversibly blocked terminator nucleotides.
All technologies above (2nd generation sequencing) require demanding protocols with serial PCR steps that results in an increased time of processing and cost. Moreover, genome regions with high density of repeats are difficult to decipher when using short reads. With the entry in the market of 3rd generation sequencers, several such hurdles have been overcome. To date, popular platforms are the Pacific Biosciences and the Oxford Nanopore technologies.
Pacific Biosciences use single molecule real-time (SMRT) technology. Library preparation leads to a closed circular DNA molecule by ligating an adaptor molecule to both ends of the target DNA molecule to be sequenced. The circular DNA molecule is then loaded into a cell containing 150,000 zeptolitre wells (ZMW)[77,78]. Each ZMW contains a DNA polymerase attached to the bottom and the target DNA fragment for sequencing. During the sequencing reaction, the DNA fragment is labeled by the DNA polymerase with fluorescent nucleotides and the corresponding emitted signal is recorded (www.pacb.com).
The Nanopore technology identifies DNA bases by measuring changes in electric conductivity generated as DNA strands pass through a biological pore. The chemical differences of each base would result, in theory, in detectably altered current flow through the pore. During the library preparation step, fragmented DNA is repaired using a PreCR step. Two adaptors are then added to the DNA, a Y adapter and a hairpin adaptor. A motor protein unzips the double stranded DNA at the Y adapter and feeds the DNA as a single strand through the nanopore (www.nanoporetech.com)[77,78].
The BGI BGISEQ-500 sequencer uses a Probe-Anchor Synthesis (cPAS) and a DNA Nanoballs (DNB) technology. The cPAS chemistry works by incorporating a fluorescent probe into a DNA anchor on the DNB, followed by high-resolution digital imaging. This combination of linear amplification and DNB technology reduces the error rate while enhancing the signal. In addition, the size of the DNB is controlled in such a way that only one DNB is bound per active site (www.bgi.com).
The rapid evolution of NGS technology keeps enhancing throughput and reducing run time and costs. Thus, NGS appears ready to offer opportunities for implementation in clinical procedures. On the other hand, though, the inevitably labor-intensive and time-demanding clinical experimentation is vastly lagging behind, thus fostering ‘home-made’ adoption of genomic tests (www.23andme.com/en-int/), in the virtual absence of experimentally-validated guidelines. Genetic information demands may create further health care disparities because of the high cost of these technologies and even pose health care risks if results are prematurely translated to the consumer market outside of regulatory protection (www.genomeweb.com/dxpgx/fda-warns-consumer-genomics-firms-illumina-selling-unapproved-dx-products#.XBLPfBNKi-U). As such, these practices are devoid of validated medical status (www.nytimes.com/2013/11/26/business/fda-demands-a-halt-to-a-dna-test-kits-marketing.html).
FDA and EMA have established recommendations to adopt pharmacogenetic and pharmacogenomic methods in research and diagnostics. Adequately numbered studies are highly recommended for achieving meaningful screening power, particularly when target mutations are rare (www.fda.gov/ucm/groups/fdagov-public/@fdagov-meddev-gen/documents/document/ucm071075.pdf). Regulatory agencies correspondingly urge to create public databases, that would include global pharmacogenetic data to provide key scientific input to both basic and clinical research (www.fda.gov/ucm/groups/fdagov-public/@fdagov-meddev-gen/documents/document/ucm509837.pdf).
Pharmacogenetic and pharmacogenomic tests will play ever more important roles in efforts to prevent adverse drug reactions. In this case, drug-metabolizing enzymes, drug transporters and drug targets are potential target genes, as they may alter the drug action or metabolism. HLA typing itself may be considered a biomarker of drug response, as it is associated to distinct haplotypes/populations with different prevalence of specific mutations in actionable genes. According to EMA, the evaluation of HLA can be carried out in drug developmental programs, in order to discover new predictive HLA biomarkers. The use of whole exome sequencing for the HLA region is strongly recommended by EMA, and is proposed to become the gold standard for HLA typing (www.ema.europa.eu/documents/scientific-guideline/guideline-good-pharmacogenomic-practice-first-version_en.pdf). Technical recommendations have also been issued, as EMA recommends that the technical predictive value of NGS should be at least 99.9%. In germline genetics, a minimum coverage of > 30× is desirable. A higher one should be pursued if a genetic variant is uncommon (www.ema.europa.eu/documents/scientific-guideline/guideline-good-pharmacogenomic-practice-first-version_en.pdf). It is expected, though, that improvements in current technologies will rapidly superseed these thresholds.
Rapid progress in deciphering cancer genomes is being achieved through ongoing international efforts, including The Cancer Genome Atlas and the International Cancer Genome Consortium. These collaborative efforts aim at sequencing up to 500 clinically well-annotated tumor groups, with the projected generation of an enormous amount of genomic data. To date, the Cancer Genome Atlas has already published initial analyses of glioblastomas and ovarian, colorectal, and breast cancers, identifying a number of deregulated genes that may set the stage for the development of targeted therapies and for associated pharmacogenetic/ pharmacogenomic biomarkers.
The shift to knowledge-driven cancer treatment requires novel classification strategies of cancer endotypes. Improvements in 2nd and 3rd generation NGS [Supplementary Tables 2 and 3] have paved the way for the use of pharmacogenomics testing at a whole-genome level.
Clinical oncology utilization of current knowledge on cancer hereditary genetics (ascopubs.org/toc/edbk/current) largely focuses on highly penetrant genes, which can account for a considerable fraction of a given cancer risk. Main examples comprise mismatch DNA repair genes, for Lynch syndrome, BRCA-1, -2, PALB2 for ovarian, breast and prostate cancer, TP53 for Li-Fraumeni syndrome (www.cancer.gov/about-cancer/causes-prevention/genetics/genetic-testing-fact-sheet#q4) and are extensively treated in other chapters of this journal issue.
Cancer genomics can extensively identify and quantify cancer inherited risk. Tumor DNA sequencing and comparison to the germline genome may identify variants associated with hereditary predisposition to cancer. Such an ability to achieve high-throughput genotyping often surpasses, though, our current ability to interpret and appropriately apply the vast amounts of data that are generated. Although known hereditary cancer susceptibility syndromes are more than one hundred, mutations in high-penetrance genes explain only a fraction of the heritability of human cancers [Figure 2].
GWAS have been conducted on nearly all common cancers. However, given the modest effect size for most risk variants, the clinical utility of genomic profiling for risk stratification based on GWAS data has been limited. A large international consortium study led to the identification of 49 new loci for breast cancer, 26 for prostate cancer, and 8 for ovarian cancer. Extension of such approaches may allow better genetic susceptibility models of cancer risk.
Despite such extensive research, only approximately 30% of familial breast cancer risk is explained by known genetic factors. WES in affected family members from 13 breast cancer families identified two families with mutations in XRCC2, including a protein-truncating change and a probable deleterious missense mutation. Another Fanconi pathway gene, SLX4 was found in only one of ≈700 BRCA-negative breast cancer kindreds. Most recently, a large pooled NGS study focusing on DNA repair pathways identified mutations in the p53-inducible protein phosphatase PPM1D as occurring mosaically in individuals with predisposition to breast and ovarian cancers. SNP-associated risk analysis identified a frameshift mutation in the BRIP1 Fanconi pathway gene, with an odds ratio of 8.1 for ovarian cancer. As indicated above, the identification of genes that are associated with homologous recombination DNA-repair has a role for targeted therapies, such as poly (ADP-ribose) polymerase inhibitors. WGS in families with multiple adenomas and/or colorectal cancer recently identified heterozygous POLE and POLD1 germline variations. GREM1 is the first gene to have been implicated in the genetic etiology of hereditary mixed polyposis syndrome. In prostate cancer, NGS of linkage regions on chromosome 17 helped to identify a mutation in the homeobox gene HOXB13 as a prostate cancer driver and potential therapy target.
NGS screening programs of cancer genomes have surprisingly led to the identification of cancer-driving mutations in tumor types that were essentially unrelated to the ones under investigations. As several driver mutations are therapy targets of drugs already in the clinics, this is rapidly leading to the repurposing of drugs that are already known to be effective against other tumor types. These serendipitous findings are now leading to a research effort that promises to be effective in adding new molecular-target therapies to tumor types often in large need of impactful treatments, such as glioblastomas and pancreatic cancer[21,27,80].
One of the first successes of NGS in this respect was the use of WES to identify risk factors that are linked to pancreatic cancer, such as PALB2 and ATM. Other examples are the E318K MITF variant, which as indicated above associates to melanoma, but was subsequently found to associate to renal cancer; the POLE and POLD1 germline variations, which are associated to colon cancer, however mutations in POLD1 also increase endometrial cancer risk; somatic mutations in BAP1, which were first identified in mesothelioma and uveal and cutaneous melanomas, were subsequently associated to renal cell cancer.
NGS protocols require highly standardized procedures for pre-analytic, analytic, and post-analytic processes. Such requirements follow the guidelines for clinical laboratory tests according to the Clinical Laboratory Improvement Amendments (CLIA) (www.cdc.gov/CLIA).
Although efficient and cost-effective, NGS has the disadvantages of high error rates and short read lengths, enrichment of rare variants, and a large proportion of missing values. A comparison of the accuracy and completeness of variant calling for two commonly used sequencing platforms found that although both technologies achieved a relatively high concordance (88%) for unique single-nucleotide variants, concordance for indel detection was only 27%.
The American College of Medical Genetics and Genomics (ACMG) has correspondingly developed a position statement for the detection of germline mutations by whole exome and genome sequencing and for the validation of NGS methods and platforms, through monitoring NGS testing, data interpretation and reporting. Similar quality assurance guidelines for NGS in diagnostic pathology are being established in Europe. The ACMG and the Association for Molecular Pathology (AMP) have jointly published recent laboratory standards and guidelines.
Accordingly, any new alteration found through NGS (WES, WGS, RNA-Seq) needs to be confirmed through another method. The main method is Sanger sequencing. Other validation methods have been proposed, such as orthogonal NGS, i.e., the parallel sequencing of target genomes using two different NGS platforms[84,85]. As for clinical applications, pre-analytic sample quality, such as tumor content of specimens, DNA integrity and yield, still are key to robust findings and must be emphasized to avoid pre-laboratory errors. Uncertainty remains on best practices regarding identification of indels and of epigenetic changes[86,87], which remain technically demanding for current NGS technologies.
As previously mentioned, storing, managing, analyzing, and interpreting genome-wide data are now rapidly becoming the bottleneck of NGS analysis procedures. Caution is also required in data interpretation, as each distinct NGS procedure has sequencing error biases [Supplementary Table 3]. The occurrence of CNVs, as that of large insertions, causes inherent difficulties (suboptimal sequence coverage and mapping quality) in NGS readings, which relate to the short-read-based nature of most NGS platforms. Recent technological advances aim at tackling this issue. High-homology regions also lead to inappropriate alignment of the reads[84,89], and this is particularly important in pharmacogenetic tests for pharmacokinetics genes, as CYP isoforms share a large sequence similarity. Also, the presence of a too high or too low GC content may affect the accessibility of gene regions. It has been estimated that more than 12% of the human exome contains “difficult or unable to be analyzed” regions for NGS because of the presence of multiple, highly homologous genes.
Another important issue is related to the creation of reference databases. Several efforts have been initiated to reach this aim, such as the National Center for Biotechnology Information-sponsored data repository ClinVar. Other large pharmacogenomic data repositories such as The Pharmacogenomics Knowledge Base (PharmGKB) developed by Stanford University (www.pharmgkb.org) and other data sets, including those available from the National Institutes of Health (NIH) GWAS collection (www.ncbi.nlm.nih.gov/snp), through the use of electronic medical records have made it possible to computationally analyze personal genomes for potential translation of pharmacogenomics into clinical practice.
Neoplasias are heterogeneous diseases that interact with complex microenvironments and batch analyses inevitably “average” such target heterogeneity. Genomic sequencing at the single-cell level bears the potential of identifying both distinct genetic cancer drivers and control networks/ druggable targets in heterogeneous neoplastic populations. Single-cell genomic assays faithfully detect somatic mitochondrial DNA mutations, track cellular relationships and hierarchies, and enable definition of clonal architecture in human cancer.
Combined single-cell RNA sequencing and genotyping can profile distinct subclones of the same tumor. This allows to identify dysregulated transcriptional programs driven by potential drug-targetable genes, with implications for targeted medicine. Single-cell DNA and RNA sequencing during neoadjuvant chemotherapy can identify patients in which treatment leads to clonal extinction vs. those in which clones persist after treatment. Notably, this analysis can reveal pre-existing resistant genotypes that could become pharmacogenomic discovery targets. Improvements of single-cell RNA-sequencing procedures are ongoing to correspondingly enhance the efficiency of such screening procedures.
Corresponding strategies can be applied to single-cell profiling and functional screening of long non-coding RNA. These approaches will critically allow extending the identification of candidate diagnostic control networks and therapeutic targets beyond the protein-coding regions of the genome.
Single cell proteomics by flow cytometry allow to track and analyse signalling events in individual cancer cells and to create signalling network maps in each cell, to identify both common fundamental regulatory themes and population heterogeneity. This can identify pathways that are activated in therapy-resistant cells and can provide predictive, actionable cancer targets. Additional data can be provided by mass cytometry, for high-dimensional, quantitative analysis of the effects of bioactive molecules on cell populations at single-cell resolution. Correlation of proteomic, transcriptomic and mutagenomic profiles with intracellular signaling molecules allows to correlate biological functions, such as metabolism, survival, DNA damage, cell cycle and apoptosis, to provide determination of network states of heterogeneous populations of individual cancer cells. Improvements in bioinformatic single-cell data comparison will be instrumental in taking the greatest advantage of single-cell measurements.
Our current knowledge of rare and common somatic/genetic variants associated with cancer risk, pharmacological treatment and disease outcome has led to significant progress, as well as to a number of challenges associated with the clinical translation of these discoveries. Improved pharmacogenetic and pharmacogenomic knowledge will extend known associations of genetic variabilily with drug responses. Key limitations of traditional approaches, i.e., exclusive focus on already known sets of genes, will most likely be overcome by global strategies, whereby the focus of the search will be extended to the entire genome, transcriptome or proteome of an individual or groups of individuals.
Better characterization of a cancer entity will doubtless contribute to improved therapy. Additional knowledge will be gathered by novel research strategies, such as those based on the concept of synthetic lethality. A synthetic lethality event occurs when two or more genes are simultaneously perturbed and exposure to a drug results in cellular or organism death/impairment. Large-scale analysis in yeast and human cells have already resulted in the identification of actionable networks of conserved, synthetic lethal interactions.
An additional level of cancer pharmacogenomics will be reached by single cell omics analyses. These will obtain information from individual cancer cells, to identify both common regulatory themes as well as population heterogeneity. This is expected to lead to the identification of subpopulations of therapy-resistant cells and of the corresponding control pathways.
A key role is all the analyses above will be provided by next-generation bioinformatic approaches, which will be instrumental in taking the greatest advantage of whole-genome sequencing at population levels and of single-cell measurements.
According to National Comprehensive Cancer Network, NGS-based pharmacogenetic/pharmacogenomic assays still are qualified as second-tier tests. The growing body of pharmacogenomic data may soon allow to go well beyond such a stage and to provide strong means to guide clinicians in the selection of safer and more effective therapies.
We thank Dr. A. Capra for web data analysis during the preparation of the manuscript and Prof. S. Cuzzocrea and L. Chiara for their unfailing support throughout the study.Authors’ contributions
Made substantial contributions to the conception and design of the study: Crisafulli C, Alberti S
Performed analyses of the meta-data from DNA sequence databases: Crisafulli C, Calabrò M, Epasto LM
Performed the comparative profiling of the NGS analyses: Romeo PD
Contributed to data acquisition and took part to the drafting, discussion and revision of the manuscript text: Crisafulli C, Romeo PD, Calabrò M, Epasto LM, Alberti SAvailability of data and materials
Not applicable.Financial support and sponsorship
The support of the Italian Ministry of Development - FESR, of the Italian Ministry of University and Research (Smart Cities and Communities SCN_00558) and of the Policlinico “G. Martino” of Messina is gratefully acknowledged. None of the funding bodies had any role in the manuscript design, in the collection, analysis and interpretation of the data, and in the writing of the manuscript.Conflicts of interest
All authors declared that there are no conflicts of interest.Ethical approval and consent to participate
Not applicable.Consent for publication
© The Author(s) 2019.
1. Relli V, Trerotola M, Guerra E, Alberti S. Distinct lung cancer subtypes associate to distinct drivers of tumor progression. Oncotarget 2018;9:35528-40.DOIPubMedPMC
2. Ambrogi F, Biganzoli E, Querzoli P, Ferretti S, Boracchi P, et al. Molecular subtyping of breast cancer from traditional tumor marker profiles using parallel clustering methods. Clin Cancer Res 2006;12:781-90.DOIPubMed
3. Querzoli P, Pedriali M, Rinaldi R, Lombardi AR, Biganzoli E, et al. Axillary lymph node nanometastases are prognostic factors for disease-free survival and metastatic relapse in breast cancer patients. Clin Cancer Res 2006;12:6696-701.DOIPubMed
4. Querzoli P, Coradini D, Pedriali M, Boracchi P, Ambrogi F, et al. An immunohistochemically positive E-cadherin status is not always predictive for a good prognosis in human breast cancer. Br J Cancer 2010;103:1835-9.DOIPubMedPMC
5. Biganzoli E, Coradini D, Ambrogi F, Garibaldi JM, Lisboa P, et al. p53 status identifies two subgroups of triple-negative breast cancers with distinct biological features. Jpn J Clin Oncol 2011;41:172-9.DOIPubMed
6. Tripaldi R, Stuppia L, Alberti S. Human height genes and cancer. BBA reviews cancer 2013;1836:27-41.DOI
7. Zanna P, Trerotola M, Vacca G, Bonasera V, Palombo B, et al. Trop-1 is a novel cell growth stimulatory molecule that marks early stages of tumor progression. Cancer 2007;110:452-64.DOI
8. Guerra E, Trerotola M, Dell’ Arciprete R, Bonasera V, Palombo B, et al. A bi-cistronic CYCLIN D1-TROP2 mRNA chimera demonstrates a novel oncogenic mechanism in human cancer. Cancer Res 2008;68:8113-21.DOIPubMed
9. Guerra E, Trerotola M, Tripaldi R, Aloisi AL, Simeone P, et al. Trop-2 induces tumor growth through Akt and determines sensitivity to Akt inhibitors. Clin Cancer Res 2016;22:4197-205.DOIPubMed
10. Ambrogi F, Fornili M, Boracchi P, Trerotola M, Relli V, et al. Trop-2 is a determinant of breast cancer survival. PLoS One 2014;9:e96993.DOIPubMedPMC
11. Guerra E, Trerotola M, Aloisi AL, Tripaldi R, Vacca G, et al. The Trop-2 signalling network in cancer growth. Oncogene 2013;32:1594-600.DOIPubMed
12. Trerotola M, Cantanelli P, Guerra E, Tripaldi R, Aloisi AL, et al. Up-regulation of Trop-2 quantitatively stimulates human cancer growth. Oncogene 2013;32:222-33.DOIPubMed
13. De Vita VT, Lawrence TS, Rosenberg SA. De Vita, Hellman & Rosenberg’s Cancer: Principles & Practice of Oncology. Available from: http://www.kubalibri.cz/files/176---DeVita,-Hellman,-and-Rosenberg-s-Cancer---Principles---Practice-of-Oncology.pdf. [Last accessed on 9 Apr 2019].
14. Tomasetti C, Li L, Vogelstein B. Stem cell divisions, somatic mutations, cancer etiology, and cancer prevention. Science 2017;355:1330-4.DOIPubMedPMC
15. Trerotola M, Relli V, Simeone P, Alberti S. Epigenetic inheritance and the missing heritability. Hum Genomics 2015;9:17.DOIPubMedPMC
16. Simeone P, Alberti S. Epigenetic heredity of human height. Physiol Rep 2014;2:e12047.DOIPubMedPMC
17. Patel JN, Wiebe LA, Dunnenberger HM, McLeod HL. Value of supportive care pharmacogenomics in oncology practice. Oncologist 2018;23:956-64.DOIPubMedPMC
18. Cascorbi I, Bruhn O, Werk AN. Challenges in pharmacogenetics. Eur J Clin Pharmacol 2013;69 Suppl 1:17-23.DOIPubMed
19. Garralda E, Paz K, Lopez-Casas PP, Jones S, Katz A, et al. Integrated next-generation sequencing and avatar mouse models for personalized cancer treatment. Clin Cancer Res 2014;20:2476-84.DOIPubMedPMC
20. Gillis NK, Patel JN, Innocenti F. Clinical implementation of germ line cancer pharmacogenetic variants during the next-generation sequencing era. Clin Pharmacol Ther 2014;95:269-80.DOIPubMedPMC
21. Adams DR, Eng CM. Next-generation sequencing to diagnose suspected genetic disorders. N Engl J Med 2018;379:1353-62.DOI
22. Hynes SO, Pang B, James JA, Maxwell P, Salto-Tellez M. Tissue-based next generation sequencing: application in a universal healthcare system. Br J Cancer 2017;116:553-60.DOIPubMedPMC
23. Klotz U. The role of pharmacogenetics in the metabolism of antiepileptic drugs: pharmacokinetic and therapeutic implications. Clin Pharmacokinet 2007;46:271-9.DOIPubMed
24. Evans WE, Relling MV. Pharmacogenomics: translating functional genomics into rational therapeutics. Science 1999;286:487-91.DOIPubMed
25. Pinto N, Cohn SL, Dolan ME. Using germline genomics to individualize pediatric cancer treatments. Clin Cancer Res 2012;18:2791-800.DOIPubMedPMC
26. Crona D, Innocenti F. Can knowledge of germline markers of toxicity optimize dosing and efficacy of cancer therapy? Biomark Med 2012;6:349-62.DOIPubMedPMC
27. Stadler ZK, Schrader KA, Vijai J, Robson ME, Offit K. Cancer Genomics and Inherited Risk. Journal of Clinical Oncology 2014;32:687-98.DOIPubMedPMC
28. Shimelis H, LaDuca H, Hu C, Hart SN, Na J, et al. Triple-negative breast cancer risk genes identified by multigene hereditary cancer panel testing. J Natl Cancer Inst 2018. Epub ahead of print, DOI: 10.1093/jnci/djy106PubMedPMC
29. Tai YC, Domchek S, Parmigiani G, Chen S. Breast cancer risk among male BRCA1 and BRCA2 mutation carriers. J Natl Cancer Inst 2007;99:1811-4.DOIPubMedPMC
30. Mavaddat N, Peock S, Frost D, Ellis S, Platte R, et al. Cancer risks for BRCA1 and BRCA2 mutation carriers: results from prospective analysis of EMBRACE. J Natl Cancer Inst 2013;105:812-22.DOIPubMed
31. Hahn SA, Greenhalf B, Ellis I, Sina-Frey M, Rieder H, et al. BRCA2 germline mutations in familial pancreatic carcinoma. JNCI Cancer Spectrum 2003;95:214-21.DOIPubMed
32. Rafnar T, Sigurjonsdottir GR, Stacey SN, Halldorsson G, Sulem P, et al. Association of BRCA2 K3326* with small cell lung cancer and squamous cell cancer of the skin. J Natl Cancer Inst 2018;110:967-74.DOIPubMedPMC
33. Guerra E, Cimadamore A, Simeone P, Vacca G, Lattanzio R, et al. p53, cathepsin D, Bcl-2 are joint prognostic indicators of breast cancer metastatic spreading. BMC Cancer 2016;16:649.DOIPubMedPMC
34. Smith KL, Adank M, Kauff N, Lafaro K, Boyd J, et al. BRCA mutations in women with ductal carcinoma in situ. Clin Cancer Res 2007;13:4306-10.DOIPubMed
35. Robson M, Im SA, Senkus E, Xu B, Domchek SM, et al. Olaparib for metastatic breast cancer in patients with a germline BRCA mutation. N Engl J Med 2017;377:1700.DOIPubMed
36. Peeters M, Douillard JY, Van Cutsem E, Siena S, Zhang K, et al. Mutant KRAS codon 12 and 13 alleles in patients with metastatic colorectal cancer: assessment as prognostic and predictive biomarkers of response to panitumumab. J Clin Oncol 2013;31:759-65.DOIPubMed
37. Jia Y, Yun CH, Park E, Ercan D, Manuia M, et al. Overcoming EGFR(T790M) and EGFR(C797S) resistance with mutant-selective allosteric inhibitors. Nature 2016;534:129-32.DOIPubMedPMC
38. Lee CK, Wu YL, Ding PN, Lord SJ, Inoue A, et al. Impact of specific epidermal growthfactor receptor (EGFR) mutations and clinical characteristics on outcomes after treatment with EGFR tyrosine kinase inhibitors versus chemotherapy in EGFR-mutant lung cancer: a meta-analysis. J Clin Oncol 2015;33:1958-65.DOIPubMed
39. Sabapathy K, Lane DP. Therapeutic targeting of p53: all mutants are equal, but some mutants are more equal than others. Nat Rev Clin Oncol 2018;15:13-30.DOIPubMed
40. Zirath H, Frenzel A, Oliynyk G, Segerstrom L, Westermark UK, et al. MYC inhibition induces metabolic changes leading to accumulation of lipid droplets in tumor cells. Proc Natl Acad Sci U S A 2013;110:10258-63.DOIPubMedPMC
41. Mason KD, Vandenberg CJ, Scott CL, Wei AH, Cory S, et al. In vivo efficacy of the Bcl-2 antagonist ABT-737 against aggressive Myc-driven lymphomas. Proc Natl Acad Sci U S A 2008;105:17961-6.DOIPubMedPMC
42. Arango D, Corner GA, Wadler S, Catalano PJ, Augenlicht LH. c-myc/p53 interaction determines sensitivity of human colon carcinoma cells to 5-fluorouracil in vitro and in vivo. Cancer Res 2001;61:4910-5.PubMed
43. Janku F, Wheler JJ, Westin SN, Moulder SL, Naing A, et al. PI3K/AKT/mTOR inhibitors in patients with breast and gynecologic malignancies harboring PIK3CA mutations. J Clin Oncol 2012;30:777-82.DOIPubMedPMC
44. Munoz DM, Cassiani PJ, Li L, Billy E, Korn JM, et al. CRISPR screens provide a comprehensive assessment of cancer vulnerabilities but generate false-positive hits for highly amplified genomic regions. Cancer Discov 2016;6:900-13.DOIPubMed
45. Tzelepis K, Koike-Yusa H, De Braekeleer E, Li Y, Metzakopian E, et al. A CRISPR dropout screen identifies genetic vulnerabilities and therapeutic targets in acute myeloid leukemia. Cell Rep 2016;17:1193-205.DOIPubMedPMC
46. Wang T, Yu H, Hughes NW, Liu B, Kendirli A, et al. Gene essentiality profiling reveals gene networks and synthetic lethal interactions with oncogenic ras. Cell 2017;168:890-903.DOIPubMedPMC
47. Shen H, Sun T, Ferrari M. Nanovector delivery of siRNA for cancer therapy. Cancer Gene Ther 2012;19:367-73.DOIPubMedPMC
48. Day FR, Thompson DJ, Helgason H, Chasman DI, Finucane H, et al. Genomic analyses identify hundreds of variants associated with age at menarche and support a role for puberty timing in cancer risk. Nat Genet 2017;49:834-41.DOIPubMedPMC
49. Dzobo K. Epigenomics-guided drug development: recent advances in solving the cancer treatment “jigsaw puzzle”. OMICS 2019;23:70-85.DOIPubMed
50. Ling A, Gruener RF, Fessler J, Huang RS. More than fishing for a cure: the promises and pitfalls of high throughput cancer cell line screens. Pharmacol Ther 2018;191:178-89.DOIPubMed
51. Rubin AJ, Parker KR, Satpathy AT, Qi Y, Wu B, et al. Coupled single-cell CRISPR screening and epigenomic profiling reveals causal gene regulatory networks. Cell 2019;176:361-76.DOIPubMedPMC
52. Nijman SM. Synthetic lethality: general principles, utility and detection using genetic screens in human cells. FEBS Lett 2011;585:1-6.DOIPubMedPMC
53. Srivas R, Shen JP, Yang CC, Sun SM, Li J, et al. A network of conserved synthetic lethal interactions for exploration of precision cancer therapy. Mol Cell 2016;63:514-25.DOIPubMedPMC
54. Dixon SJ, Costanzo M, Baryshnikova A, Andrews B, Boone C. Systematic mapping of genetic interaction networks. Annu Rev Genet 2009;43:601-25.DOIPubMed
55. Mengwasser KE, Adeyemi RO, Leng Y, Choi MY, Clairmont C, et al. Genetic screens reveal FEN1 and APEX2 as BRCA2 synthetic lethal targets. Mol Cell 2019;73:885-99.DOIPubMed
56. Tharp AP, Maffini MV, Hunt PA, VandeVoort CA, Sonnenschein C, et al. Bisphenol A alters the development of the rhesus monkey mammary gland. Proc Natl Acad Sci U S A 2012;109:8190-5.DOIPubMedPMC
57. Carletti E, Guerra E, Alberti S. The forgotten variables of DNA array hybridization. Trends Biotechnol 2006;24:443-8.DOIPubMed
58. Kalemkerian GP, Narula N, Kennedy EB, Biermann WA, Donington J, et al. Molecular testing guideline for the selection of patients with lung cancer for treatment with targeted tyrosine kinase inhibitors: American society of clinical oncology endorsement of the college of American pathologists/international association for the study of lung cancer/association for molecular pathology clinical practice guideline update. J Clin Oncol 2018;36:911-9.DOIPubMed
59. Antonia SJ, Villegas A, Daniel D, Vicente D, Murakami S, et al. Overall survival with durvalumab after chemoradiotherapy in stage III NSCLC. New England Journal of Medicine 2018;379:2342-50.DOIPubMed
60. Ortega VE, Meyers DA. Pharmacogenetics: implications of race and ethnicity on defining genetic profiles for personalized medicine. J Allergy Clin Immunol 2014;133:16-26.DOIPubMedPMC
61. Rabbani B, Nakaoka H, Akhondzadeh S, Tekin M, Mahdieh N. Next generation sequencing: implications in personalized medicine and pharmacogenomics. Mol Biosyst 2016;12:1818-30.DOIPubMed
62. Garraway LA. Genomics-driven oncology: framework for an emerging paradigm. J Clin Oncol 2013;31:1806-14.DOIPubMed
63. Dong L, Wang W, Li A, Kansal R, Chen Y, et al. Clinical next generation sequencing for precision medicine in cancer. Curr Genomics 2015;16:253-63.DOIPubMedPMC
64. Filipski KK, Mechanic LE, Long R, Freedman AN. Pharmacogenomics in oncology care. Front Genet 2014;5:73.DOIPubMedPMC
65. Caudle KE, Klein TE, Hoffman JM, Muller DJ, Whirl-Carrillo M, et al. Incorporation of pharmacogenomics into routine clinical practice: the Clinical Pharmacogenetics Implementation Consortium (CPIC) guideline development process. Curr Drug Metab 2014;15:209-17.DOIPubMedPMC
66. Relling MV, Schwab M, Whirl-Carrillo M, Suarez-Kurtz G, Pui CH, et al. Clinical Pharmacogenetics Implementation Consortium (CPIC) guideline for thiopurine dosing based on TPMT and NUDT15 genotypes: 2018 update. Clin Pharmacol Ther 2018. Available from: https://ascpt.onlinelibrary.wiley.com/doi/full/10.1002/cpt.1304. [Last accessed on 9 Apr 2019].
67. Luzum JA, Pakyz RE, Elsey AR, Haidar CE, Peterson JF, et al. The pharmacogenomics research network translational pharmacogenetics program: outcomes and metrics of pharmacogenetic implementations across diverse healthcare systems. Clin Pharmacol Ther 2017;102:502-10.DOIPubMedPMC
68. Romani A, Guerra M, Trerotola M, Alberti S. Detection and analysis of spliced chimeric mRNAs in sequence databanks. Nucleic Acids Res 2003;31:1-8.
69. Plebani R, Oliver GR, Trerotola M, Guerra E, Cantanelli P, et al. Long-range transcriptome sequencing reveals cancer cell growth regulatory chimeric mRNA. Neoplasia 2012;14:1087-96.DOIPubMedPMC
70. Alberti S, Herzenberg LA. DNA methylation prevents transfection of genes for specific surface antigens. Proc Natl Acad Sci USA 1988;85:8391-4.PubMedPMC
71. Alberti S, Nutini M, Herzenberg LA. DNA methylation prevents the amplification of TROP1, a tumor associated cell surface antigen gene. Proc Natl Acad Sci USA 1994;91:5833-7.PubMedPMC
72. Nevedomskaya E, Wessels L, Zwart W. Genome-wide epigenetic profiling of breast cancer tumors treated with aromatase inhibitors. Genom Data 2014;2:195-8.DOIPubMedPMC
73. Margulies M, Egholm M, Altman WE, Attiya S, Bader JS, et al. Genome sequencing in microfabricated high-density picolitre reactors. Nature 2005;437:376-80.DOIPubMedPMC
74. Voelkerding KV, Dames SA, Durtschi JD. Next-generation sequencing: from basic research to diagnostics. Clin Chem 2009;55:641-58.DOIPubMed
75. Kchouk M, Gibrat JF, Elloumi M. Generations of sequencing technologies: from first to next generation. Biol Med 2017;9:8.DOI
76. Besser J, Carleton HA, Gerner-Smidt P, Lindsey RL, Trees E. Next-generation sequencing technologies and their application to the study and control of bacterial infections. Clin Microbiol Infect 2018;24:335-41.DOIPubMedPMC
77. Norris AL, Workman RE, Fan Y, Eshleman JR, Timp W. Nanopore sequencing detects structural variants in cancer. Cancer Biol Ther 2016;17:246-53.DOIPubMedPMC
78. Cao Y, Fanning S, Proos S, Jordan K, Srikumar S. A review on the applications of next generation sequencing technologies as applied to food-related microbiome studies. Front Microbiol 2017;8:1829.DOIPubMedPMC
79. Simeone P, Alberti S. RE: HABP2 G534E mutation in familial nonmedullary thyroid cancer. J Natl Cancer Inst 2016;108.DOIPubMedPMC
80. Herbst RS, Morgensztern D, Boshoff C. The biology and management of non-small cell lung cancer. Nature 2018;553:446-54.DOIPubMed
81. Rehm HL, Bale SJ, Bayrak-Toydemir P, Berg JS, Brown KK, et al. ACMG clinical laboratory standards for next-generation sequencing. Genet Med 2013;15:733-47.DOIPubMedPMC
82. Dubbink HJ, Deans ZC, Tops BB, van Kemenade FJ, Koljenovic S, et al. Next generation diagnostic molecular pathology: critical appraisal of quality assurance in Europe. Mol Oncol 2014;8:830-9.DOIPubMedPMC
83. Li MM, Datto M, Duncavage EJ, Kulkarni S, Lindeman NI, et al. Standards and guidelines for the interpretation and reporting of sequence variants in cancer: a joint consensus recommendation of the association for molecular pathology, American society of clinical oncology, and College of American pathologists. J Mol Diagn 2017;19:4-23.DOIPubMedPMC
84. Schwarz UI, Gulilat M, Kim RB. The role of next-generation sequencing in pharmacogenetics and pharmacogenomics. Cold Spring Harb Perspect Med 2019. Epub ahead of print [PMID: 29844222 DOI: 10.1101/cshperspect.a033027]PubMed
85. Chennagiri N, White EJ, Frieden A, Lopez E, Lieber DS, et al. Orthogonal NGS for high throughput clinical diagnostics. Sci Rep 2016;6:24650.DOIPubMedPMC
86. Nasr AF, Nutini M, Palombo B, Guerra E, Alberti S. Mutations ofTP53 induce loss of DNA methylation and amplification of the TROP1 gene. Oncogene 2003;22:1668-77.DOIPubMed
87. Li Y, Melnikov AA, Levenson V, Guerra E, Simeone P, et al. A seven-gene CpG-island methylation panel predicts breast cancer progression. BMC Cancer 2015;15:417.DOIPubMedPMC
88. Lauschke VM, Ingelman-Sundberg M. Requirements for comprehensive pharmacogenetic genotyping platforms. Pharmacogenomics 2016;17:917-24.DOIPubMed
89. Mandelker D, Schmidt RJ, Ankala A, McDonald Gibson K, Bowser M, et al. Navigating highly homologous genes in a molecular diagnostic setting: a resource for clinical next-generation sequencing. Genet Med 2016;18:1282-9.DOIPubMed
90. Rasmussen-Torvik LJ, Stallings SC, Gordon AS, Almoguera B, Basford MA, et al. Design and anticipated outcomes of the eMERGE-PGx project: a multicenter pilot for preemptive pharmacogenomics in electronic health record systems. Clin Pharmacol Ther 2014;96:482-9.DOIPubMedPMC
91. Chu WK, Edge P, Lee HS, Bansal V, Bafna V, et al. Ultraaccurate genome sequencing and haplotyping of single human cells. Proceedings of the National Academy of Sciences 2017;114:12512-7.DOIPubMedPMC
92. Kim C, Gao R, Sei E, Brandt R, Hartman J, et al. Chemoresistance evolution in triple-negative breast cancer delineated by single-cell sequencing. Cell 2018;173:879-93.DOIPubMedPMC
93. Zhou J, Xu J, Zhang L, Liu S, Ma Y, et al. Combined single-cell profiling of lncRNAs and functional screening reveals that H19 is pivotal for embryonic hematopoietic stem cell development. Cell Stem Cell 2019;24:285-98.DOIPubMed
94. Irish JM, Kotecha N, Nolan GP. Mapping normal and cancer cell signalling networks: towards single-cell proteomics. Nat Rev Cancer 2006;6:146-55.DOIPubMed
95. Treutlein B, Camp JG. Evaluation of Wu et al.: strengthening single-cell bioinformatic comparisons for improving organoid differentiation. Cell Stem Cell 2019;24:343-4.DOIPubMed
Jun Koseki et al., Cancer Drug Resistance, 2019
Miki Mori et al., Cancer Drug Resistance, 2019
Jan Willem van den Berg et al., Mini-invasive Surgery, 2019
Ilaria Fiordoliva et al., Journal of Cancer Metastasis and Treatment, 2019
Pooja Murthy et al., Cancer Drug Resistance, 2019
Jose E. Belizario et al., Cancer Drug Resistance, 2019
Elizabeth Harrod et al., Hepatoma Research, 2019
Yukio Nishimura et al., Cancer Drug Resistance, 2019
Jin-Liang Zhang et al., Hepatoma Research, 2019
Gerhard Hamilton et al., Cancer Drug Resistance, 2019