- Research note
- Open Access
Monitoring dynamic cytotoxic chemotherapy response in castration-resistant prostate cancer using plasma cell-free DNA (cfDNA)
BMC Research Notesvolume 12, Article number: 275 (2019)
Cell-free DNA (cfDNA) is an attractive cancer biomarker, as it is thought to reflect a component of the underlying genetic makeup of the tumor and is readily accessible in serial fashion. Because chemotherapy regimens are expected to act rapidly on cancer and cfDNA is cleared from the blood within minutes, we hypothesized that cfDNA would reflect immediate effects of treatment. Here, we developed a method for monitoring long cfDNA fragments, and report dynamic changes in response to cytotoxic chemotherapy.
Peripheral blood was obtained from 15 patients with metastatic castration-resistant prostate cancer (CRPC) immediately before and after cytotoxic chemotherapy infusion. cfDNA was extracted and quantified for long interspersed nuclear elements (LINE1; 297 bp) using qPCR. Targeted deep sequencing was performed to quantify the frequency of mutations in exon 8 of the androgen receptor (AR), a mutational hotspot region in CRPC. Single nucleotide mutations in AR exon 8 were found in 6 subjects (6/15 = 40%). Analytical variability was minimized by pooling independent PCR reactions for each library. In 5 patients, tumor-derived long cfDNA levels were found to change immediately after infusion. Detailed analysis of one subject suggests that cytotoxic chemotherapy can produce rapidly observable effects on cfDNA.
Recent studies have begun to elucidate the relationship between molecular marker status in the blood and therapy response in castration-resistant prostate cancer (CRPC) [1,2,3]. At the genetic level, AR point mutations are associated with resistance to therapies [1, 2, 4,5,6,7]. Clinically, response to therapy is typically evaluated via measurement of prostate-specific antigen (PSA) at regular intervals, with PSA decline rate being shown to inform on overall survival [8, 9]. However, PSA response can take weeks to months to present, and changes may be hard to interpret, as they can reflect both PSA release from cell lysis and production from intact cells. Hence, to develop personalized treatment strategies there is a need to more rapidly monitor patients while receiving cancer treatment and to obtain information reflective of the molecular make-up of their tumors.
Plasma cell-free DNA (cfDNA) can be readily accessed in serial fashion to analyze the dynamic clonal make-up of cancers across tumor evolution and therapy [1, 10,11,12,13,14,15,16,17,18,19,20,21,22]. This is especially valuable in metastatic disease, where sequential access to tissues is not feasible. In prostate cancer, researchers have used cfDNA to track changes in AR mutation status in response to the targeted therapies abiraterone and enzalutamide [1, 2, 19, 20, 23, 24]. Both overall cfDNA levels and specific mutations have been suggested to inform on response to therapy and outcome [21, 22, 25,26,27].
Previous studies analyzed cfDNA after several weeks to months of treatment. The goal of the current study was to look for therapy effects at shorter time frames, just hours after administration of therapy. Because the half-life of cfDNA is < 1 h [28,29,30,31], this approach may offer the potential for monitoring dynamic changes in response to treatment. We therefore developed a targeted approach to track AR mutations in plasma to attempt to correlate with systemic response in CRPC patients at the time of chemotherapy administration.
Materials and methods
Plasma was obtained from patients with metastatic prostate cancer treated with docetaxel-based chemotherapy as part of a correlative clinical trial. Dexamethasone and antiemetics were administered as per standard of care. Serial samples were obtained immediately before, 30–60 min after, and in our case study, 24 h after docetaxel infusion. Baseline patient demographics are presented in Additional file 1.
EDTA tubes of whole blood were collected and processed within 2 h. Plasma was removed after centrifugation at 1600 rpm at 4 °C and stored in 1 ml aliquots at − 80 °C.
cfDNA extraction and quantitation
Plasma samples were thawed on ice and centrifuged at 4 °C to remove residual cells. cfDNA was extracted from 200 µl aliquots using the QIAamp DNA Mini Blood Mini Kit (Qiagen). Samples were eluted in 100 µl and 2 µl were quantified for long interspersed nuclear elements (LINE1, 297 bp amplicon) using qPCR as described previously .
Amplification of the androgen receptor exon 8
Two sequences that spanned exon 8 of the AR were amplified using Phusion polymerase (New England Biolabs) and the following primers: 5′-CCACCTCCTTGTCAACCCTG-3′ and 5′-GATCTCTGCCATCATTTCCGG-3′ for amplicon A1 (138 bp), and 5′-CAAGTCACACATGGTGAGCG-3′ and 5′-GGTTTCCAATGCTTCACTGGGT-3′ for amplicon A2 (128 bp). DNA (1 µl) was amplified as follows: (1) 98 °C 30 s, (2) 30 cycles: 98 °C 10 s, 58 °C 10 s, 72 °C 10 s, (3) 72 °C 5 min. PCR products were purified using Agencourt AMPure XP beads at 1.2× concentration (Beckman Coulter). Reactions were eluted in 10 µl.
Targeted sequencing of androgen receptor exon 8
Amplicon DNA concentrations were determined using the BioAnalyzer with the High Sensitivity DNA Analysis Kit (Agilent Technologies). Libraries were generated with Kapa Hyper Prep Kit (Kapa Biosystems). Barcoded libraries were quantified using the BioAnalyzer, pooled and sequenced using an Illumina MiSeq. For each sample, a pool of 10 independent PCR reactions was analyzed to reduce technical variability. Results from a single patient sample with two single nucleotide variants (SNVs) are demonstrated in Additional file 1, where results from the pooled sample were equivalent to the average of the observed frequencies of each product weighted by its coverage.
To estimate the sensitivity of our approach, samples with known SNV frequencies were sequenced. We amplified variants of the AR using two different plasmids as template, (1) pEGFP-C1-AR  (a gift from Dr. Michael Mancini, Addgene plasmid # 28,235) (AR-LBD WT) and (2) p270 pCMV-CRE-M-AR(LNCaP) (a gift from Dr. Jeffrey Green Addgene plasmid # 12,496) (AR-LBD T878A). Amplicons were generated using the following primers: 5′-CCTGCTCAAGACGCTTCTACC-3′ and 5′-TCACTGGGTGTGGAAATAGATGG-3′. Purified amplicons were diluted to 350 ng/ml final concentration, mixed at various ratios and sequenced. Additional file 2 compares the observed fraction (99.6%, 78.0%, 50.4%,13.3%, 5.9%, 0.4%) with the actual fraction of mutant DNA (100%, 75%, 50%, 10%, 5%, 0%), highlighting the high accuracy of our pipeline in quantifying SNVs (R2 = 0.999) with low absolute error (root mean square RMS = 1.9%).
After demultiplexing and adapter/barcode trimming, reads were mapped to human chromosome X (hg19) using BWA v. 0.7.9 . FixMateInformation and MarkDuplicateWithMateCigar from the package picard v. 1.130 (http://broadinstitute.github.io/picard) were used to ensure and verify correct mate-pair information. Resulting alignments were sorted and indexed using samtools v. 1.3.1  and visually inspected using the Integrated Genome Viewer . Variant positions in the target region were identified via visual inspection and read counts for all observed nucleotides at the variant position were obtained using custom Java code based on the picard library (http://broadinstitute.github.io/picard). The variant frequency was calculated as the number of reads with the observed alternate allele divided by the total number of reads at that position and expressed as a percentage. Statistical analyses were performed using the statistical language R v. 3.2.3 (http://www.R-project.org/) and RStudio v. 0.99 (http://www.rstudio.com/).
cfDNA was isolated from 43 plasma samples from 15 subjects and analyzed with qPCR of LINE1 (297 bp), a method in previous studies used to quantify yields (Additional file 1) [32, 37]. The median (range) LINE1 value observed across all samples was 13.1 (0.4–819.8) ng/ml. In the 17 sample pairs obtained immediately before and after docetaxel administration, LINE1 appeared to decrease after chemotherapy (pre: 13.1 (1.3–143.6) vs. post: 8.5 (0.4–311.1 ng/ml). However, a Wilcoxon signed rank test showed that the change was not significant (V= 69, p = 0.75).
Targeted sequencing of AR exon 8
We amplified two regions spanning exon 8 of the AR, which contains a cluster of AR mutation hotspots reported in CRPC . Amplified product was used for library construction and subsequent massive parallel sequencing (median depth > 1.6 million reads). We processed a total of 34 samples from 15 patients with a mean Q30 of 95.3%, SD 1.9% (Additional file 1), and clinically informative genomic profiling of cfDNA was feasible in all samples.
Detection of genomic alterations in AR exon 8 in patients with CRPC
To identify patients with mutations in the exon 8 of AR, we sequenced cfDNA from 15 patients’ plasma at a single time point. We set our threshold of SNV identification to ≥ 1% frequency. The cut-off was chosen based on our estimate of approximately 100 AR template molecules in each PCR reaction (rough calculation based on LINE1 DNA concentration). Assuming cfDNA is for the most part derived from tumor, we considered 1% a conservative threshold and appropriate for these analyses. Exact frequencies are reported in Additional file 1. Using this threshold, 6 (40%) of the 15 patients analyzed were found to have at least one SNV (Additional file 1). These results are in line with previously reported data [4, 5]. A total of 5 SNVs were identified across all time points, 4 of which were previously reported in CRPC: H875Y and T878A, D891H and Q903H (Table 1, Additional file 1). The exception D880Y has been associated with androgen insensitivity syndrome , but not with cancer. Genomic coordinates and COSMIC database IDs are listed in Additional file 1.
List of patients with at least one AR mutation. Specific SNVs listed in subsequent columns. Numbers indicate the frequency range of SNVs measured across sequential samples for each patient. The detection threshold for valid SNVs was set to frequencies ≥ 1%.
Rapid changes in cfDNA post cytotoxic chemotherapy
Next, we asked whether the taxane-based chemotherapy could affect the dynamics of tumor-derived cfDNA immediately after administration. Therefore, we expanded our sequencing assay to sequential sample draws from patients with detected AR mutations. In 5 patients, we detected SNVs at ≥ 1% frequency in both pre- and post-chemotherapy samples (Fig. 1, Additional file 1). In patient P7, we measured the greatest change in SNV frequency after chemotherapy (61% to 79%). Due to considerable changes in two SNV frequencies in Patient 17, we decided to analyze the case in more detail.
Patient P17 had metastatic prostate cancer initially treated with androgen deprivation therapy (leuprolide acetate and bicalutamide) for 4 months before progression to castration-resistant disease. He was then enrolled in a clinical trial and treated with 4 cycles docetaxel with bevacizumab and everolimus  and responded with a transient 49% decline in PSA followed by rapid clinical progression with PSA elevation and symptomatic bone pain. The patient died from prostate cancer 5 months after docetaxel was initiated.
Figure 2 depicts the patient’s dynamic response across treatment. Prior to the first cycle of treatment, total cfDNA content measured by LINE1 was 34.4 ng/ml at baseline. Targeted deep sequencing of AR exon 8 revealed two independent mutations at 11% (H875Y) and 13% (T878A), respectively. Since the sequencing reads contained one or the other alteration, and AR amplifications and mutations are described as mutually exclusive [19, 41], we conclude that the 2 mutations represent two different clones within the patient’s tumor. Immediately after docetaxel infusion, LINE1 decreased to 0.4 ng/ml, while the frequencies of the two AR mutations increased to 15% (H875Y) and 18% (T878A). After 24 h, LINE1increased to 18.1 ng/ml, representing a recovery of 52%, while mutant fraction values remained stably elevated at 15% (H875Y) and 20% (T878A) (Fig. 2b, c). The response pattern in cycle 2 was similar: baseline LINE1 was 67.4 ng/ml. The two independent mutations were measured 3% (H875Y) and 7% (T878A). Immediately after docetaxel infusion, LINE1 decreased to 8.5 ng/ml, while the SNV frequencies increased to 6% (H875Y) and 17% (T878A). After 24 h, LINE1 increased to 22.8 ng/ml, representing a recovery of 34%, while mutant fraction values increased further to 9% (H875Y) and 31% (T878A) (Fig. 2b, c).
In summary, our assay measured dynamic changes in long cfDNA fragments (297 bp) within hours of repeated cycles of cytotoxic chemotherapy, suggestive of response. Future studies involving more patients are required to determine the prognostic and theragnostic value of these short-term dynamics.
Recent studies highlight the dynamics of cfDNA to elucidate mechanisms of resistance to targeted therapies in CRPC that are typically acquired over a course of several months [1, 2, 19, 20, 23]. Here, we demonstrate the measurement of changes in cfDNA profiles within 1 h of chemotherapy infusion.
In the subject we present here in detail, we observed a rapid decline in circulating cfDNA content measured by qPCR of long LINE1. The reported mean half-life of DNA entering the circulation ranges from minutes up to a few hours before it is degraded and cleared by the kidneys, liver and spleen [28, 29]. These data support our hypothesis that a rapid response to therapy might be reflected in the blood circulation. To answer this question will require additional studies including subjects without chemotherapy treatment to determine normal fluctuation of cfDNA.
Surprisingly, the decline in long LINE1 in our case study coincided with elevated allele frequency of mutant AR. These results suggest that chemotherapy affected cfDNA derived from the tumor differently than from healthy tissue. Furthermore, the two separate mutations increased at varying rates, indicative of different clonal dynamics. Taken together, these results suggest that the heterogeneous response dynamics of different tumor clones can be closely monitored using cfDNA. If confirmed in a larger patient study, this could impact the development of patient-specific drug combination therapies and dosing strategies through real-time monitoring.
We cannot rule out that differences in concomitant medications or intra-patient variations influenced the results. For example, we did not specifically control for glucocorticoid use although high doses dexamethasone are usually given as pre-medication before docetaxel . Similarly, we did not longitudinally sample patients off of chemotherapy. It is possible that some of the variations we observed reflect underlying hour-by-hour changes independent of therapeutic intervention. Finally, the small number of patients analyzed across various courses of treatment is a limitation.
Availability of data and materials
The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.
castration-resistant prostate cancer
long interspersed nuclear elements 1
prostate-specific antigen PSA
root mean square
Carreira S, Romanel A, Goodall J, Grist E, Ferraldeschi R, Miranda S, et al. Tumor clone dynamics in lethal prostate cancer. Sci Transl Med. 2014;6(254):254ra125.
Lallous N, Volik SV, Awrey S, Leblanc E, Tse R, Murillo J, et al. Functional analysis of androgen receptor mutations that confer anti-androgen resistance identified in circulating cell-free DNA from prostate cancer patients. Genome Biol. 2016;17(1):1–15.
Handy CE, Antonarakis ES. Sequencing treatment for castration-resistant prostate cancer. Curr Treat Options Oncol. 2016;17(12):64.
Taplin ME, Bubley GJ, Ko YJ, Small EJ, Upton M, Rajeshkumar B. Selection for androgen receptor mutations in prostate cancers treated with androgen antagonist. Cancer Res. 1999;59(11):2511–5.
Taplin ME, Bubley GJ, Shuster TD, Frantz ME, Spooner AE, Ogata GK. Mutation of the androgen-receptor gene in metastatic androgen-independent prostate cancer. N Engl J Med. 1995;332(21):1393–8.
Culig Z, Hoffmann J, Erdel M, Eder IE, Hobisch A, Hittmair A, et al. Switch from antagonist to agonist of the androgen receptor bicalutamide is associated with prostate tumour progression in a new model system. Br J Cancer. 1999;81:242–51.
Balbas MD, Evans MJ, Hosfield DJ, Wongvipat J, Arora VK, Watson PA. Overcoming mutation-based resistance to antiandrogens with rational drug design. Elife. 2013;2:e00499.
Hanninen M, Venner P, North S. A rapid PSA half-life following docetaxel chemotherapy is associated with improved survival in hormone refractory prostate cancer. Can Urol Assoc J. 2009;3(5):369–74.
Choueiri TK, Xie W, D’Amico AV, Ross RW, Hu JC, Pomerantz M, et al. Time to PSA nadir independently predicts overall survival in metastatic hormone sensitive prostate cancer patients treated with androgen deprivation therapy. Cancer. 2009;115(5):981–7.
Manuscript A, Malignancies H, Bettegowda C, Sausen M, Leary RJ, Kinde I, et al. Detection of circulating tumor DNA in early- and late-stage human malignancies. Sci Transl Med. 2014;6(224):224ra24.
Ulz P, Belic J, Graf R, Auer M, Lafer I, Fischereder K, et al. Whole-genome plasma sequencing reveals focal amplifications as a driving force in metastatic prostate cancer. Nat Commun. 2016;7:12008.
O’Leary B, Hrebien S, Morden JP, Beaney M, Fribbens C, Huang X, et al. Early circulating tumor DNA dynamics and clonal selection with palbociclib and fulvestrant for breast cancer. Nat Commun. 2018;9(1):896.
Murtaza M, Dawson S-J, Tsui DWY, Gale D, Forshew T, Piskorz AM, et al. Non-invasive analysis of acquired resistance to cancer therapy by sequencing of plasma DNA. Nature. 2013;497(7447):108–12.
Murtaza M, Dawson SJ, Pogrebniak K, Rueda OM, Provenzano E, Grant J, et al. Multifocal clonal evolution characterized using circulating tumour DNA in a case of metastatic breast cancer. Nat Commun. 2015;6:8760.
Mok T, Wu Y-L, Lee JS, Yu C-J, Sriuranpong V, Sandoval-Tan J, et al. Detection and dynamic changes of EGFR mutations from circulating tumor DNA as a predictor of survival outcomes in NSCLC patients treated with first-line intercalated erlotinib and chemotherapy. Clin Cancer Res. 2015;21(14):3196–203.
Spindler KLG, Pallisgaard N, Andersen RF, Jakobsen A. Changes in mutational status during third-line treatment for metastatic colorectal cancer—results of consecutive measurement of cell free DNA, KRAS and BRAF in the plasma. Int J Cancer. 2014;135(9):2215–22.
Dawson S-J, Tsui DWY, Murtaza M, Biggs H, Rueda OM, Chin S-F, et al. Analysis of circulating tumor DNA to monitor metastatic breast cancer. N Engl J Med. 2013;368(13):1199–209.
Xi L, Pham T, Payabyab EC, Sherry RM, Rosenberg SA, Raffeld M. Circulating tumor DNA as an early indicator of response to T-cell transfer immunotherapy in metastatic melanoma. Clin Cancer Res. 2016;22:5480–6.
Romanel A, Tandefelt DG, Conteduca V, Jayaram A, Casiraghi N, Wetterskog D. Plasma AR and abiraterone-resistant prostate cancer. Sci Transl Med. 2015;7(312):312re10.
Wyatt AW, Azad AA, Volik SV, Annala M, Beja K, McConeghy B, et al. Genomic alterations in cell-free DNA and enzalutamide resistance in castration-resistant prostate cancer. JAMA Oncol. 2016;2(12):1598–606.
Garcia-Murillas I, Schiavon G, Weigelt B, Ng C, Hrebien S, Cutts RJ, et al. Mutation tracking in circulating tumor DNA predicts relapse in early breast cancer. Sci Transl Med. 2015;7(302):302ra133.
Hendriks RJ, Dijkstra S, Smit FP, Vandersmissen J, Van de Voorde H, Mulders PFA, et al. Epigenetic markers in circulating cell-free DNA as prognostic markers for survival of castration-resistant prostate cancer patients. Prostate. 2018;78(5):336–42.
Azad AA, Volik SV, Wyatt AW, Haegert A, Bihan S, Bell RH. Androgen receptor gene aberrations in circulating cell-free DNA: biomarkers of therapeutic resistance in castration-resistant prostate cancer. Clin Cancer Res. 2015;21(10):2315–24.
Belic J, Graf R, Bauernhofer T, Cherkas Y, Ulz P, Waldispuehl-Geigl J, et al. Genomic alterations in plasma DNA from patients with metastasized prostate cancer receiving abiraterone or enzalutamide. Int J Cancer. 2018;143(5):1236–48.
Kwee S, Song M-A, Cheng I, Loo L, Tiirikainen M. Measurement of circulating cell-free DNA in relation to 18F-fluorocholine PET/CT imaging in chemotherapy-treated advanced prostate cancer. Clin Transl Sci. 2012;5(1):65–70.
Kienel A, Porres D, Heidenreich A, Pfister D. cfDNA as a prognostic marker of response to taxane based chemotherapy in patients with prostate cancer. J Urol. 2015;194(4):966–71.
Conteduca V, Wetterskog D, Sharabiani MTA, Grande E, Fernandez-Perez MP, Jayaram A, et al. Androgen receptor gene status in plasma DNA associates with worse outcome on enzalutamide or abiraterone for castration-resistant prostate cancer: a multi-institution correlative biomarker study. Ann Oncol. 2017;28(7):1508–16.
Emlen W, Mannik M. Effect of DNA size and strandedness on the in vivo clearance and organ localization of DNA. Clin Exp Immunol. 1984;56(1):185–92.
Lo YMD, Zhang J, Leung TN, Lau TK, Chang AMZ, Hjelm NM. Rapid clearance of fetal DNA from maternal plasma. Am J Hum Genet. 1999;64(1):218–24.
Yu SCY, Lee SWY, Jiang P, Leung TY, Chan KCA, Chiu RWK, et al. High-resolution profiling of fetal DNA clearance from maternal plasma by massively parallel sequencing. Clin Chem. 2013;59(8):1228 LP-1237.
Diehl F, Schmidt K, Choti MA, Romans K, Goodman S, Li M, et al. Circulating mutant DNA to assess tumor dynamics. Nat Med. 2008;31(14):985–90.
Sunami E, Vu A-T, Nguyen SL, Giuliano AE, Hoon DSB. Quantification of LINE1 in circulating DNA as a molecular biomarker of breast cancer. Ann N Y Acad Sci. 2008;1137(1):171–4.
Stenoien DL, Cummings CJ, Adams HP, Mancini MG, Patel K, DeMartino GN, et al. Polyglutamine-expanded androgen receptors form aggregates that sequester heat shock proteins, proteasome components and SRC-1, and are suppressed by the HDJ-2 chaperone. Hum Mol Genet. 1999;8(5):731–41.
Li H, Durbin R. Fast and accurate long-read alignment with Burrows-Wheeler transform. Bioinformatics. 2010;26(5):589–95.
Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, et al. The sequence alignment/map format and SAMtools. Bioinformatics. 2009;25(16):2078–9.
Thorvaldsdóttir H, Robinson JT, Mesirov JP. Integrative Genomics Viewer (IGV): high-performance genomics data visualization and exploration. Br Bioinform. 2013;14(2):178–92.
Wei L, Wu W, Han L, Yu W, Du Y. A quantitative analysis of the potential biomarkers of non-small cell lung cancer by circulating cell-free DNA. Oncol Lett. 2018;16:4353–60.
Gottlieb B, Beitel LK, Nadarajah A, Paliouras M, Trifiro M. The androgen receptor gene mutations database: 2012 update. Hum Mutat. 2012;33(5):887–94.
Chávez B, Méndez JP, Ulloa-Aguirre A, Larrea F, Vilchis F. Eight novel mutations of the androgen receptor gene in patients with androgen insensitivity syndrome. J Hum Genet. 2001;46(10):560–5.
Gross ME, Dorff TB, Quinn DI, Diaz PM, Castellanos OO, Agus DB. Safety and efficacy of docetaxel, bevacizumab, and everolimus for castration-resistant prostate cancer (CRPC). Clin Genitourin Cancer. 2018;16(1):e11–21.
Robinson D, Van Allen EM, Wu YMM, Schultz N, Lonigro RJJ, Mosquera JMM, et al. Integrative clinical genomics of advanced prostate cancer. Cell. 2015;161(5):1215–28.
Hesketh PJ, Kris MG, Basch E, Bohlke K, Barbour SY, Clark-Snow RA, et al. Antiemetics: American Society of Clinical Oncology clinical practice guideline update. J Clin Oncol. 2017;35(28):3240–61.
We thank Dr. Michael Mancini for the pEGFP-C1-AR construct and Dr. Jeffrey Green for the p270 pCMV-CRE-M-AR(LNCaP). We thank Olga Castellanos, Brandi Scott, and the clinical and research staff at the USC Westside Cancer Center, the USC Norris Comprehensive Cancer Center and the Next Generation Sequencing Core Facility at the USC Epigenome Center. This work was partially supported by National Cancer Institute Cancer Center Shared Grant award P30CA014089. The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Cancer Institute or the National Institutes of Health who had no influence on the design, analysis, or conclusions of the study. Also thanks to the patients who graciously volunteered to participate in this trial. We would like to express our deepest gratitude to our philanthropic supporters: the Stephenson family, Emmet, Toni and Tessa, for their donation to support this research.
This work was partially supported by National Cancer Institute Cancer Center Shared Grant award P30CA014089. The funding body did not play a role in the design of the study and collection, analysis, or interpretation of data or in writing the manuscript.
Ethics approval and consent to participate
Ethical approval and oversight for prospective sample collection was performed with the consent of all research subjects under auspices of the University of Southern California Institutional Review Board (IRB Protocol number: HS-11-00450). Prospective, voluntary written informed consent was obtained from all subjects prior to data and specimen collection.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.