Patent Application
20240279755
The text of the computer readable sequence listing filed herewith, titled “39554-601_SQL_ST25”, created Jun. 9, 2022, having a file size of 2,873 bytes, is hereby incorporated by reference in its entirety.
The present disclosure relates to materials and methods for measuring human papillomavirus (HPV) circulating tumor DNA (ctDNA). In particular, the disclosure provides a droplet digital PCR-based assay for HPV16 ctDNA, and the use thereof for predicting response to treatment in patients having HPV positive head and neck squamous cell carcinoma (HNSCC), such as oropharyngeal squamous cell carcinoma (OPSSC).
Head and neck squamous cell carcinomas (HNSCC) constitute 3-5% of all malignancies worldwide and there are approximately 600,000 newly diagnosed cases annually.1,2 The majority of patients with HNSCC present with locally and/or regionally advanced disease at diagnosis. Despite the use of combined modality treatment, a significant proportion of patients develop unresectable recurrent or metastatic disease (R/M) HNSCC, for which treatment is palliative and consists of systemic therapy.3-5 Even with the development of novel treatment regimens including immunotherapy, the median survival for patients with R/M HNSCC is typically less than one year.6 There is an increasing incidence of human papillomavirus associated oropharyngeal squamous cell carcinoma (HPV+OPSCC) and this represents one of only four cancers increasing in incidence in the United States.7
A number of prognostic clinical factors have been identified for patients with recurrent/metastatic (R/M) HNSCC including HPV status8. However, the only validated predictive biomarker is PD-L1, allowing for identification of patients who will have survival benefit with pembrolizumab versus chemotherapy as a first line treatment. Despite this, treatment response rates with currently available therapeutics remain low and no biomarker has been validated for dynamic response assessment or prediction of treatment benefit including PD-L19-11. Given this, efficacy of therapy can only able to be assessed after 9-12 weeks of treatment since radiographic responses are often delayed.6,9,11,12 As a result of the limited response rate, many patients experience significant progression leading to airway compromise and impaired functional status, rendering them unsuitable for further cancer directed therapy. A method is needed to rapidly identify those not benefiting from therapy, as well as to avoid progressive disease on a futile therapy and enable an earlier switch to a potentially efficacious treatment. A clear need exists for treatment individualization in R/M HNSCC, which necessitates identification of predictive biomarkers that can inform about response to therapy in real-time. Additionally, such an assay, may have utility for the detection of early recurrence for previously untreated patients, or for detection of response to therapy in early treatment settings.
In some aspects, provided herein are methods of detecting HPV16 circulating tumor DNA (ctDNA). In some embodiments, methods for detecting HPV16 ctDNA comprise contacting a sample comprising cell free DNA with a set of oligonucleotides comprising a forward primer having at least 90% sequence identity to SEQ ID NO: 1, a reverse primer having at least 90% sequence identity to SEQ ID NO: 2; and a probe comprising a detectable label. The forward primer and the reverse primer anneal to a target HPV16 sequence. The methods further comprise amplifying the target HPV16 sequence, if present in the sample, and detecting a signal from the detectable label. Detection of the signal indicates the presence of the target HPV16 sequence in the sample. In some embodiments, the probe comprises a fluorescent label.
In some aspects, methods of detecting HPV16 ctDNA comprise contacting a sample comprising cell free DNA with a set of oligonucleotides comprising a forward primer, a reverse primer, and a probe comprising a detectable label. The methods further comprise amplifying the target HPV16 sequence, if present in the sample, comprising the sequence of SEQ ID NO: 4 or a portion thereof comprising at least 15 continuous bases of SEQ ID NO:4. The methods further comprise detecting a signal from the detectable label, wherein detection of the signal indicates the presence of the target HPV16 sequence in the sample.
For any of the methods described herein, amplifying the target HPV16 sequence may be performed by quantitative PCR (qPCR). For example, the qPCR may be digital PCR. In some embodiments, the PCR is droplet digital PCR.
In some aspects, provided herein are methods of quantifying HPV16 circulating tumor DNA (ctDNA). In some embodiments, methods of quantifying HPV16 ctDNA comprise contacting a sample comprising cell free DNA with a set of oligonucleotides comprising a forward primer having at least 90% sequence identity to SEQ ID NO: 1, a reverse primer having at least 90% sequence identity to SEQ ID NO: 2; and a probe comprising a detectable label. The forward primer and the reverse primer anneal to a target HPV16 sequence. In some embodiments, the methods further comprise amplifying the target HPV16 sequence, if present in the sample, by droplet digital PCR. In some embodiments, the methods further comprise detecting a signal from the detectable label, wherein detection of the signal indicates the presence of the target HPV16 sequence. The methods further comprise quantifying the amount of the target HPV16 sequence present in the sample.
In some embodiments, methods for quantifying HPV16 circulating tumor DNA comprise contacting a sample comprising cell free DNA with a set of oligonucleotides comprising a forward primer, a reverse primer, and a probe comprising a detectable label, wherein the forward primer and the reverse primer anneal to a target HPV16 sequence. The methods further comprise amplifying the target HPV16 sequence, if present in the sample, by droplet digital PCR. In some embodiments, the target HPV16 sequence comprises the sequence of SEQ ID NO: 4 or a portion thereof comprising at least 15 contiguous bases of SEQ ID NO:4. The method further comprises detecting a signal from the detectable label, wherein detection of the signal indicates the presence of the target HPV16 sequence, and quantifying the amount of the target HPV16 sequence present in the sample.
For any of the methods described herein, the probe may be an oligonucleotide having at least 90% sequence identity to SEQ ID NO: 3. The cell free DNA may be isolated form any suitable sample. In some embodiments, the cell free DNA is isolated from a plasma sample. In some embodiments, the cell free DNA is isolated from a saliva sample.
For any of the methods described herein, the forward primer may comprise the nucleotide sequence of SEQ ID NO: 1, the reverse primer may comprise the nucleotide sequence of SEQ ID NO: 2, and the probe may comprise the nucleotide sequence of SEQ ID NO: 3.
In some aspects, provided herein is a set of oligonucleotides for amplifying and detecting a target HPV16 sequence. In some embodiments, the set of oligonucleotides comprises a forward primer having at least 90% sequence identity to SEQ ID NO: 1; a reverse primer having at least 90% sequence identity to SEQ ID NO: 2; and a probe having at least 90% sequence identity to SEQ ID NO: 3.
In some embodiments, the forward primer comprises the nucleotide sequence of SEQ ID NO: 1. In some embodiments, the reverse primer comprises the nucleotide sequence of SEQ ID NO: 2. In some embodiments, the probe comprises the nucleotide sequence of SEQ ID NO: 3. In some embodiments, the probe comprises a detectable label. In some embodiments, the detectable label is a fluorescent label.
In some embodiments, the set of oligonucleotides comprises a forward primer and a reverse primer that generates an amplicon comprising the sequence of SEQ ID NO: 4 or a portion thereof comprising at least 15 contiguous bases of SEQ ID NO:4; and a probe that hybridizes to the amplicon. In some embodiments, the probe comprises a detectable label. In some embodiments, the detectable label is a fluorescent label.
In some aspects, provided herein are methods of predicting response to treatment in a subject with HPV positive head and neck squamous cell carcinoma (HPV+HNSCC). In some embodiments, the method comprises measuring a baseline level HPV16 ctDNA in the subject, measuring a follow on level of HPV16 ctDNA following one or more treatment sessions in the subject, and predicting response to treatment in the subject. In some embodiments, predicting response comprises predicting a positive response to treatment in the subject when a change from the baseline level to the follow on level is below a threshold value. In some embodiments, predicting response comprises predicting a negative response to treatment in the subject when a change from the baseline level to the follow on level is above a threshold value.
In some embodiments, the baseline level is measured prior to any treatment in the subject. In some embodiments, the baseline level is measured following a first treatment cycle and the follow on level is measured following a second treatment cycle.
In some embodiments, measuring the baseline level and/or measuring the follow on level may be performed using a suitable method as described herein. In some embodiments, measuring the baseline level and measuring the follow on level comprises obtaining a sample comprising cell free DNA from the subject, contacting the sample with a set of oligonucleotides comprising a forward primer having at least 90% sequence identity to SEQ ID NO: 1, a reverse primer having at least 90% sequence identity to SEQ ID NO: 2; and a probe comprising a detectable label, wherein the forward primer and the reverse primer anneal to a target HPV16 sequence. In some embodiments, the method further comprises amplifying the target HPV16 sequence, if present in the sample, by droplet digital PCR. In some embodiments, the method further comprises detecting a signal from the detectable label, wherein detection of the signal indicates the presence of the target HPV16 sequence; and quantifying the amount of the target HPV16 present in the sample.
In some embodiments, measuring the baseline level and measuring the follow on level comprises obtaining a sample comprising cell free DNA from the subject, and contacting the sample with a set of oligonucleotides comprising a forward primer, a reverse primer, and a probe comprising a detectable label, wherein the forward primer and the reverse primer anneal to a target HPV16 sequence. In some embodiments, the method further comprises amplifying the target HPV16 sequence, if present in the sample, by droplet digital PCR. The target HPV16 sequences may comprise the sequence of SEQ ID NO: 4 or a portion thereof comprising at least 15 contiguous bases of SEQ ID NO:4. In some embodiments, the method further comprises detecting a signal from the detectable label, wherein detection of the signal indicates the presence of the target HPV16 sequence, and quantifying the amount of the target HPV16 present in the sample.
For any of the methods of predicting response described herein, the probe may comprise a fluorescent label. The probe may be an oligonucleotide having at least 90% sequence identity to SEQ ID NO: 3. The cell free DNA may be isolated from any suitable sample, including a plasma sample or a saliva sample obtained from the subject. In some embodiments, the forward primer comprises the nucleotide sequence of SEQ ID NO: 1, the reverse primer comprises the nucleotide sequence of SEQ ID NO: 2, and the probe comprises the nucleotide sequence of SEQ ID NO: 3.
In some embodiments, the threshold value is 50%. In some embodiments, the threshold value is 60%.
A positive response to treatment may indicate a complete response to treatment, a partial response to treatment, or a stable disease state in the subject. A negative response to treatment may indicate disease progression in the subject. In some embodiments, the treatment is immunotherapy. In some embodiments, the HNSCC is oropharyngeal squamous cell carcinoma (OPSCC).
In some aspects, provided herein are methods of treating a subject with HPV positive head and neck squamous cell carcinoma (HPV+HNSCC). In some embodiments, the method comprises measuring a baseline level HPV16 ctDNA in the subject, measuring a follow on level of HPV16 ctDNA following one or more treatment sessions in the subject, and modulating therapy based upon a change from baseline to the follow on level. In some embodiments, modulating therapy comprises increasing the dose and/or frequency of therapy or providing an alternative or additional therapy to the subject when the change from baseline to the follow on level is above a threshold value. In some embodiments, the baseline level is measured prior to any treatment in the subject. In some embodiments, the baseline level is measured following a first treatment cycle and the follow on level is measured following a second treatment cycle.
In some embodiments, measuring the baseline level and measuring the follow on level comprises obtaining a sample comprising cell free DNA from the subject, and contacting the sample with a set of oligonucleotides comprising a forward primer having at least 90% sequence identity to SEQ ID NO: 1, a reverse primer having at least 90% sequence identity to SEQ ID NO: 2; and a probe comprising a detectable label. In some embodiments, the forward primer and the reverse primer anneal to a target HPV16 sequence. In some embodiments, the method further comprises amplifying the target HPV16 sequence, if present in the sample, by droplet digital PCR, detecting a signal from the detectable label, wherein detection of the signal indicates the presence of the target HPV16 sequence, and quantifying the amount of the target HPV16 present in the sample.
In some embodiments, measuring the baseline level and measuring the follow on level comprises obtaining a sample comprising cell free DNA from the subject, and contacting the sample with a set of oligonucleotides comprising a forward primer, a reverse primer, and a probe comprising a detectable label. The forward primer and the reverse primer anneal to a target HPV16 sequence. In some embodiments, the method further comprises amplifying the target HPV16 sequence, if present in the sample, by droplet digital PCR. In some embodiments, the target HPV16 sequences comprises the sequence of SEQ ID NO: 4 or a portion thereof comprising at least 15 contiguous bases of SEQ ID NO:4. In some embodiments, the method further comprises detecting a signal from the detectable label, wherein detection of the signal indicates the presence of the target HPV16 sequence. In some embodiments, the method further comprises quantifying the amount of the target HPV16 present in the sample.
For any of the methods of measuring the baseline and/or follow on level, the probe may comprise a fluorescent label. The probe may be an oligonucleotide having at least 90% sequence identity to SEQ ID NO: 3. The cell free DNA may be isolated from a suitable sample, including plasma or saliva obtained from the subject.
In some embodiments, the forward primer comprises the nucleotide sequence of SEQ ID NO: 1, the reverse primer comprises the nucleotide sequence of SEQ ID NO: 2, and the probe comprises the nucleotide sequence of SEQ ID NO: 3. In some embodiments, the threshold value is 50%. In some embodiments, the threshold value is 60%. In some embodiments, the treatment is immunotherapy. In some embodiments, the HNSCC is oropharyngeal squamous cell carcinoma (OPSCC).
Although any methods and materials similar or equivalent to those described herein can be used in the practice or testing of embodiments described herein, some preferred methods, compositions, devices, and materials are described herein. However, before the present materials and methods are described, it is to be understood that this invention is not limited to the particular molecules, compositions, methodologies, or protocols herein described, as these may vary in accordance with routine experimentation and optimization. It is also to be understood that the terminology used in the description is for the purpose of describing the particular versions or embodiments only, and is not intended to limit the scope of the embodiments described herein.
Unless otherwise defined, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this invention belongs. However, in case of conflict, the present specification, including definitions, will control. Accordingly, in the context of the embodiments described herein, the following definitions apply.
As used herein and in the appended claims, the singular forms “a”, “an” and “the” include plural reference unless the context clearly dictates otherwise.
As used herein, the term “comprise” and linguistic variations thereof denote the presence of recited feature(s), element(s), method step(s), etc. without the exclusion of the presence of additional feature(s), element(s), method step(s), etc. Conversely, the term “consisting of” and linguistic variations thereof, denotes the presence of recited feature(s), element(s), method step(s), etc. and excludes any unrecited feature(s), element(s), method step(s), etc., except for ordinarily-associated impurities. The phrase “consisting essentially of” denotes the recited feature(s), element(s), method step(s), etc. and any additional feature(s), element(s), method step(s), etc. that do not materially affect the basic nature of the composition, system, or method. Many embodiments herein are described using open “comprising” language. Such embodiments encompass multiple closed “consisting of” and/or “consisting essentially of” embodiments, which may alternatively be claimed or described using such language.
An “amplicon” refers to a nucleic acid fragment formed as a product of natural or artificial amplification events or techniques. For example, an amplicon can be produced by PCR
As used herein, the terms “subject” and “patient” are used interchangeably herein and refer to both human and nonhuman animals. The term “nonhuman animals” includes all vertebrates, e.g., mammals and non-mammals, such as nonhuman primates, sheep, dogs, cats, horses, cows, chickens, amphibians, reptiles, and the like. In some embodiments, the subject is a human. In some embodiments, the subject is a human. In particular embodiments, the subject may be male. In other embodiments, the subject may be female. In some embodiments, the subject is suffering from cancer.
As used herein, “treat”, “treating”, “treatment”, and variations thereof refer to the clinical intervention made in response to a disease, disorder or physiological condition manifested by a patient or to which a patient may be susceptible. The aim of treatment includes the alleviation or prevention of symptoms, slowing or stopping the progression or worsening of a disease, disorder, or condition and/or the remission of the disease, disorder, or condition. In some embodiments, treating a cancer refers to the management and care of the subject for combating and reducing one or more symptoms of the cancer. For example, treating cancer may reduce tumor burden (e.g. reduce the size of one or more tumors in the subject afflicted with cancer and/or reduce the overall number of tumors in the subject afflicted with cancer). Treating a cancer may reduce or completely eliminate the cancer (e.g. completely eliminate the tumor) in the subject.
Solid tumors are known to shed circulating tumor DNA (ctDNA) which might be detectable bodily fluids including saliva, urine, and plasma.13-15 Plasma ctDNA analyses allow for non-invasive longitudinal monitoring of tumor specific genomic alterations. To date a comprehensive prospective analysis of ctDNA assay characteristics and predictive efficacy has not been performed in head and neck cancer.
Described herein is the development of an HPV16 ctDNA test that offers a precise assay for detection of HPV+OPSCC. Moreover, the assay described herein finds use to predict progressive disease prior to radiographic imaging in patients undergoing treatment. In some embodiments, a highly sensitive and specific droplet digital PCR (ddPCR) assay for absolute quantification of HPV ctDNA from plasma specimens is described herein.
For any of the methods described herein, any suitable sample type may be used. The sample may be obtained from the subject and subsequently used for any of the methods described herein. In some embodiments, the sample is obtained from the subject and circulating tumor DNA is isolated from the sample for use in the methods described herein. Suitable samples include fluids (e.g. urine, blood, blood products, sputum, saliva, etc.), solids, tissues, and gases. In some embodiments, the sample is blood (e.g. whole blood) or a blood product such as plasma, serum, and the like. In some embodiments, the sample is a plasma sample. In some embodiments, the sample is a sample obtained from the mouth of the subject. For example, the sample may be a saliva sample. The term “saliva” or “saliva sample” is meant to include any sample containing saliva from the subject, including spit, an oral swab or sponge sample, a mouthwash rinse sample, etc.
Methods of Detecting and/or Quantifying HPV ctDNA
In some aspects, provided herein are methods of detecting HPV circulating tumor DNA (ctDNA). In some aspects, provided herein are methods for detecting HPV subtype 16 (HPV16) ctDNA. In some embodiments, the methods comprise contacting a sample comprising cell free DNA with a set of oligonucleotides. In some embodiments, the set of oligonucleotides comprises a forward primer, a reverse primer, and a probe.
In some embodiments, the set of oligonucleotides generate an amplicon comprising the sequence of: CGAAACCGGTTAGTATAAAAGCAGACATTTTATGCACCAAAAGAG AACTGCAATGTTTCAGGACCCACAGGAGCGAC (SEQ ID NO: 4), or a portion thereof having at least 15 contiguous bases of SEQ ID NO:4 (e.g., at least 15, 18, 20, 25, 30, etc.). In some embodiments, the set of oligonucleotides comprise a forward primer and a reverse primer that generate an amplicon of SEQ ID NO: 4 or a portion thereof, and the probe comprises an oligonucleotide capable of hybridizing to the amplicon.
In some embodiments, the forward primer has at least 90% sequence identity to CGAAACCGGTTAGTATAAAAGCAG (SEQ ID NO: 1). For example, the forward primer may have at least 90%, at least 91%, at least 92%, at least 93%, at least 94%, at least 95%, at least 96%, at least 97%, at least 98%, at least 99%, or 100% sequence identity to SEQ ID NO: 1. In some embodiments, the reverse primer has at least 90% sequence identity to GTCGCTCCTGTGGGTCCT (SEQ ID NO: 2). For example, the reverse primer may have at least 90%, at least 91%, at least 92%, at least 93%, at least 94%, at least 95%, at least 96%, at least 97%, at least 98%, at least 99%, or 100% sequence identity to SEQ ID NO: 2.
In some embodiments, the probe comprises a detectable label. Any suitable detectable label may be used. In some embodiments, the detectable label is a fluorescent label. Suitable fluorescent labels include, for example, FAM (5- or 6-carboxyfluorescein), VIC, NED, Fluorescein, FITC, IRD-700/800, CY3, CY5, CY3.5, CY5.5, HEX, TET, TAMRA, JOE, ROX, BODIPY TMR, Oregon Green, Rhodamine Green, Rhodamine Red, Texas Red, Yakima Yellow, Alexa Fluor PET, Biosearch Blue™, Marina Blue™, Bothell Blue®, Alexa Fluor®, 350 FAM™ SYBR™ Green 1, Fluorescein, EvaGreen™, Alexa Fluor® 488 JOE™, VIC™ HEX™ TET™ CAL Fluor® Gold 540, Yakima Yellow®, ROX™, CAL Fluor® Red 610, Cy3.5™, Texas Red®, Alexa Fluor® 0.568 CyS™, Quasar™ 670, LightCycler Red640@, Alexa Fluor 633 Quasar™ 705, LightCycler Red705®. Alexa Fluor® 680, SYTO™ 9, LC Green®, LC Green® Plus+, and EvaGreen™.
In some embodiments, the probe comprises an oligonucleotide having at least 90% sequence identity to CATTTTATGCACCAAAAGAGAACTGCAATGTTTC (SEQ ID NO: 3). For example, the probe may comprise an oligonucleotide having at least 90%, at least 91%, at least 92%, at least 93%, at least 94%, at least 95%, at least 96%, at least 97%, at least 98%, at least 99%, or 100% sequence identity to SEQ ID NO: 3.
The forward primer and the reverse primer anneal to a target HPV16 sequence (e.g. a target region within the HPV16 genome). HPVs, including HPV16, contain a 7.9-kb circular double-stranded DNA genome that consists of four parts: an early region (E1, 2, 4-7 genes), a late region (L1, 2 genes), a long control region (LCR), and a small, highly variable, non-coding region (NCR) between E5 and L2. The early region and late region of HPV16 are shown in
In some embodiments, the method comprises contacting a sample comprising cell free DNA with a set of oligonucleotides comprising a forward primer having at least 90% sequence identity to SEQ ID NO: 1, a reverse primer having at least 90% sequence identity to SEQ ID NO: 2; and a probe comprising an oligonucleotide having at least 90% sequence identity to SEQ ID NO: 3 and a detectable label. In some embodiments, the forward primer comprises the nucleotide sequence of SEQ ID NO: 1, the reverse primer comprises the nucleotide sequence of SEQ ID NO: 2, and the probe comprises the nucleotide sequence of SEQ ID NO: 3.
The methods further comprise amplifying the target HPV16 sequence, if present in the sample. The methods further comprise detecting a signal from the detectable label, which is indicative of the presence of the target HPV16, if present in the sample. For example, amplification may be performed by polymerase chain reaction (PCR). In some embodiments, the PCR is a quantitative PCR (qPCR). For example, amplification of the target sequence may be performed using a digital PCR technique, such as a dPCR technique selected from droplet digital PCR (ddPCR), BEAMing (beads, emulsion, amplification, and magnetic), and microfluidic chips.
As used herein, “digital PCR” refers to an assay that provides an end-point measurement that provides the ability to quantify nucleic acids without the use of standard curves, as is used in real-time PCR. In a typical digital PCR experiment, the sample is randomly distributed into discrete partitions, such that some contain no nucleic acid template and others contain one or more template copies. The partitions are amplified to the terminal plateau phase of PCR (or end-point) and then read to determine the fraction of positive partitions. If the partitions are of uniform volume, the number of target DNA molecules present may be calculated from the fraction of positive end-point reactions using Poisson statistics, according to the following equation.
wherein λ is the average number of target DNA molecules per replicate reaction and p is the fraction of positive end-point reactions. From k, together with the volume of each replicate PCR and the total number of replicates analyzed, an estimate of the absolute target DNA concentration is calculated. Digital PCR includes a variety of formats, including droplet digital PCR, BEAMing (beads, emulsion, amplification, and magnetic), and microfluidic chips.
In some embodiments, the dPCR is droplet digital PCR. “Droplet digital PCR” (ddPCR) refers to a digital PCR assay that measures absolute quantities by counting nucleic acid molecules encapsulated in discrete, volumetrically defined, water-in-oil droplet partitions that support PCR amplification (Hinson et al., 2011, Anal. Chem. 83:8604-8610, Pinheiro et al., 2012. Anal. Chem. 84:1003-1011). A single ddPCR reaction may be comprised of at least 20,000 partitioned droplets per well. A “droplet” or “water-in-oil droplet” refers to an individual partition of the droplet digital PCR assay. A droplet supports PCR amplification of template molecule(s) using homogenous assay chemistries and workflows similar to those widely used for real-time PCR applications (Hinson et al., 2011, Anal. Chem. 83.8604-8610; Pinheiro et al., 2012, Anal. Chem. 84:1003-1011).
Droplet digital PCR may be performed using any platform that performs a digital PCR assay that measures absolute quantities by counting nucleic acid molecules encapsulated in discrete, volumetrically defined, water-in-oil droplet partitions that support PCR amplification. The strategy for droplet digital PCR may be summarized as follows: a sample is diluted and partitioned into thousands to millions of separate reaction chambers (water-in-oil droplets) so that each contains one or no copies of the nucleic acid molecule of interest. The number of “positive” droplets detected, which contain the target amplicon (i.e., nucleic acid molecule of interest), versus the number of “negative” droplets, which do not contain the target amplicon (i.e., nucleic acid molecule of interest), may be used to determine the number of copies of the nucleic acid molecule of interest that were in the original sample. Examples of droplet digital PCR systems include the QX100™ Droplet Digital PCR System by Bio-Rad, which partitions samples containing nucleic acid template into 20,000 nanoliter-sized droplets; the QX200™ Droplet Digital PCR System by Bio-Rad; and the RainDrop™ digital PCR system by RainDance, which partitions samples containing nucleic acid template into 1,000,000 to 10,000,000 picoliter-sized droplets.
In some embodiments, the method further comprises quantifying HPV16 ctDNA in the sample.
In some aspects, provided herein is a method for quantifying HPV16 ct DNA in a sample. The method may be substantially the same as the above-described methods for detecting HPV16 ctDNA in the sample, and further comprise quantifying the HPV16 ctDNA, if present in the sample. In some embodiments, the method comprises contacting a sample comprising cell free DNA with a set of oligonucleotides comprising a forward primer having at least 90% sequence identity to SEQ ID NO: 1, a reverse primer having at least 90% sequence identity to SEQ ID NO: 2; and a probe comprising a detectable label. In some embodiments, the method further comprises amplifying the target HPV16 sequence, if present in the sample, by droplet digital PCR.
In some embodiments, the target HPV16 sequence comprises SEQ ID NO. 4 or a portion thereof having at least 15 contiguous bases of SEQ ID NO:4. For example, in some embodiments the method comprises contacting a sample comprising cell free DNA with a set of oligonucleotides capable of generating an amplicon comprising SEQ ID NO: 4 or a portion thereof. For example, in some embodiments the set of oligonucleotides comprise a forward primer and a reverse primer that generates an amplicon comprising SEQ ID NO: 4 or a portion thereof, and a probe that hybridizes to the amplicon.
The method further comprises detecting a signal from the detectable label, wherein detection of the signal indicates the presence of the target HPV16 sequence, and quantifying the amount of the target HPV16 present in the sample. Quantification may be performed, for example, by determining the number of target DNA molecules present. The number of target DNA molecules present may be calculated from the fraction of positive end-point reactions using Poisson statistics, as described above.
The methods of detecting and/or quantifying HPV16 ctDNA may be used in a variety of diagnostic or predictive methods. In some embodiments, the methods of detecting or the methods of quantifying HPV16 ctDNA may be used in methods of determining whether a subject has an HPV positive cancer. For example, the methods may be used to determine whether a subject has an HPV positive cervical cancer, head and neck cancer (e.g. oropharyngeal cancer), anal cancer, penile cancer, vaginal cancer, or vulvar cancer. For example, the methods may be used to determine whether a subject has an HPV positive head and neck cancer, such as HPV positive oropharyngeal cancer.
In some embodiments, the methods for detecting or the methods for quantifying HPV16 ctDNA may be used in methods of predicting patient response to therapy in a cancer. Such methods are described in further detail below.
In some aspects, provided herein are methods of predicting response to treatment in a subject. In some embodiments, provided herein are methods for predicting response to treatment in a subject suffering from cancer. The treatment may be any suitable treatment/therapy for cancer. Suitable therapies include, for example, surgery, chemotherapy, radiation therapy, targeted therapy, immunotherapy, stem cell transplant therapy, bone marrow transplant therapy, or hormone therapy. In some embodiments, the treatment is immunotherapy. As used herein, the term “immunotherapy” refers to any type of cancer treatment that helps the immune system fight cancer. In some embodiments, “immunotherapy” refers to treatment with an immune checkpoint inhibitor, T-cell transfer therapy, monoclonal antibodies, treatment vaccines, and/or immune modulators. The cancer may be any cancer type. In some embodiments, the cancer is an HPV positive cancer. For example, the cancer may be an HPV positive cancer selected from cervical cancer, head and neck cancer (e.g. oropharyngeal cancer), anal cancer, penile cancer, vaginal cancer, or vulvar cancer. For example, in some embodiments provided herein is a method of predicting response immunotherapy in a subject suffering from HPV positive head and neck squamous cell carcinoma. For example, the HPV positive head and neck squamous cell carcinoma maybe HPV positive oropharyngeal squamous cell carcinoma.
In some embodiments, the method of predicting response comprises measuring a baseline level HPV16 ctDNA in the subject. The baseline level of HPV16 ctDNA may be measured prior to any treatment in the subject. Alternatively, the baseline level of HPV16 ctDNA may be measured following one or more treatment sessions in the subject. For example, the cancer therapy (e.g., immunotherapy) may be provided to the subject over the course of one or more cycles. Each cycle may comprise a single immunotherapy session (e.g. dose) or multiple sessions (e.g. multiple doses, over the course of multiple days). In some embodiments, the baseline level is measured after a therapy (e.g. an immunotherapy) cycle.
The method further comprises measuring a follow on level of HPV16 ctDNA following one or more treatment sessions in the subject. For example, a baseline level of HPV16 ctDNA may be measured prior to any therapy, and a follow on level may be measured following a therapy cycle. Alternatively, a baseline level of HPV16 ctDNA may be measured after a first therapy cycle, and a follow on level may be measured following a second therapy cycle.
The baseline level and/or follow on level of HPV16 ctDNA may be measured using a suitable method described herein. In some embodiments, measuring the baseline level and measuring the follow on level comprises obtaining a sample comprising cell free DNA from the subject. In some embodiments, the method further comprises contacting the sample with a set of oligonucleotides comprising a forward primer having at least 90% sequence identity to SEQ ID NO: 1, a reverse primer having at least 90% sequence identity to SEQ ID NO: 2; and a probe comprising a detectable label. In some embodiments, the method comprises contacting a sample comprising cell free DNA with a set of oligonucleotides capable of generating an amplicon comprising SEQ ID NO: 4 or a portion thereof.
In some embodiments, the detectable label is a fluorescent label. Suitable fluorescent labels are described above. In some embodiments, the probe is an is an oligonucleotide having at least 90% sequence identity to SEQ ID NO: 3. In some embodiments, forward primer and the reverse primer anneal to a target HPV16 sequence. The method may further comprise amplifying the target HPV16 sequence, if present in the sample, by droplet digital PCR. In some embodiments, the method further comprises detecting a signal from the detectable label, wherein detection of the signal indicates the presence of the target HPV16 sequence. In some embodiments, the method further comprises quantifying the amount of the target HPV16 present in the sample.
The cell free DNA may be isolated from plasma obtained from the subject. For example, plasma may be obtained and cell free DNA may be isolated using any suitable method. In some embodiments, cell free DNA is isolated using a commercially available kit.
The method further involves predicting response to therapy in the subject based upon HPV16 ctDNA levels. In some embodiments, the method comprises predicting a positive response to treatment in the subject when a change from the baseline level to the follow on level is below a threshold value. In some embodiments, the method comprises predicting a negative response to treatment in the subject when a change from the baseline level to the follow on level is above a threshold value. In some embodiments, the threshold value is 50%. In some embodiments, the threshold value is 60%. For example, the method may comprise predicting a positive response to treatment when there is less than a 60% increase in HPV16 ctDNA levels from the baseline level to the follow on level. Alternatively, the method may comprise predicting a negative response to treatment when there a 60% or higher increase in HPV16 ctDNA levels from the baseline level to the follow on level. A positive response to treatment may indicate a complete response to treatment, a partial response to treatment, or a stable disease state in the subject. For example, a positive response may indicate no metastasis, no increase in tumor size, and/or no increase in total number of tumors observed in the subject. A negative response may indicate progression of the disease, such as increased total number of tumors and/or increased tumor size.
In some aspects, provided herein are methods of treating cancer in a subject. The subject may be undergoing or about to begin a suitable treatment/therapy for cancer. Suitable therapies include, for example, surgery, chemotherapy, radiation therapy, targeted therapy, immunotherapy, stem cell transplant therapy, bone marrow transplant therapy, or hormone therapy. In some embodiments, the treatment is immunotherapy. As used herein, the term “immunotherapy” refers to any type of cancer treatment that helps the immune system fight cancer. In some embodiments, “immunotherapy” refers to treatment with an immune checkpoint inhibitor, T-cell transfer therapy, monoclonal antibodies, treatment vaccines, and/or immune modulators. The cancer may be any cancer type. In some embodiments, the cancer is an HPV positive cancer. For example, the cancer may be an HPV positive cancer selected from cervical cancer, head and neck cancer (e.g. oropharyngeal cancer), anal cancer, penile cancer, vaginal cancer, or vulvar cancer. For example, in some embodiments provided herein is a method of treating a subject suffering from HPV positive head and neck squamous cell carcinoma. For example, the HPV positive head and neck squamous cell carcinoma maybe HPV positive oropharyngeal squamous cell carcinoma.
In some embodiments, the method of treating the subject comprises measuring a baseline level HPV16 ctDNA in the subject. The baseline level of HPV16 ctDNA may be measured prior to any treatment in the subject. Alternatively, the baseline level of HPV16 ctDNA may be measured following one or more treatment sessions in the subject. For example, the cancer therapy (e.g., immunotherapy) may be provided to the subject over the course of one or more cycles. Each cycle may comprise a single immunotherapy session (e.g. dose) or multiple sessions (e.g. multiple doses, over the course of multiple days). In some embodiments, the baseline level is measured after a therapy (e.g. an immunotherapy) cycle. The method further comprises measuring a follow on level of HPV16 ctDNA following one or more treatment sessions in the subject. For example, a baseline level of HPV16 ctDNA may be measured prior to any therapy, and a follow on level may be measured following a therapy cycle. Alternatively, a baseline level of HPV16 ctDNA may be measured after a first therapy cycle, and a follow on level may be measured following a second therapy cycle.
The baseline level and/or follow on level of HPV16 ctDNA may be measured using any suitable method. In some embodiments, the baseline level and/or follow on level is measured using a method described herein. In some embodiments, measuring the baseline level and measuring the follow on level comprises obtaining a sample comprising cell free DNA from the subject. In some embodiments, the method further comprises contacting the sample with a set of oligonucleotides comprising a forward primer having at least 90% sequence identity to SEQ ID NO: 1, a reverse primer having at least 90% sequence identity to SEQ ID NO: 2; and a probe comprising a detectable label. In some embodiments, the method comprises contacting a sample comprising cell free DNA with a set of oligonucleotides capable of generating an amplicon comprising SEQ ID NO: 4 or a portion thereof.
In some embodiments, the probe comprises a detectable label. In some embodiments, the detectable label is a fluorescent label. Suitable fluorescent labels are described above. In some embodiments, the probe is an is an oligonucleotide having at least 90% sequence identity to SEQ ID NO: 3. In some embodiments, forward primer and the reverse primer anneal to a target HPV16 sequence. The method may further comprise amplifying the target HPV16 sequence, if present in the sample, by droplet digital PCR. In some embodiments, the method further comprises detecting a signal from the detectable label, wherein detection of the signal indicates the presence of the target HPV16 sequence. In some embodiments, the method further comprises quantifying the amount of the target HPV16 present in the sample.
The cell free DNA may be isolated from plasma obtained from the subject. For example, plasma may be obtained and cell free DNA may be isolated using any suitable method. In some embodiments, cell free DNA is isolated using a commercially available kit.
The method further involves modulating therapy in the subject based upon HPV16 ctDNA levels. In some embodiments, the method comprises increasing the dose and/or frequency of the therapy and/or providing an alternative or additional therapy to the subject when the change from baseline to the follow on level is above a threshold value. For example, the method may comprise increasing the dose and/or frequency of the therapy that has already been provided to the subject (e.g. the therapy provided to the subject prior to measuring the follow on level of HPV16 ctDNA). For example, the method may comprise increasing the frequency of therapy. As another example, the method may comprise increasing the dose of therapy. Increasing the dose of therapy may refer to increasing the amount given per individual dose and/or increasing the total number of doses given to the subject. In some embodiments, the method may comprise providing an alternative therapy to the subject. For example, the method may comprise switching from the first form of therapy to a second, different form a therapy. In some embodiments, the method may comprise providing an additional therapy to the subject.
In some embodiments, the threshold value is 50%. In some embodiments, the threshold value is 60%. For example, in some embodiments the method comprises increasing the dose and/or frequency of the therapy and/or providing an alternative or additional therapy to the subject when the change from the baseline level to the follow on level is above 50% (e.g. at least 50%, at least 55%, at least 60%, at least 65%, at least 70%, at least 75%, etc.).
In some aspects, provided herein is a set of oligonucleotides. The set of oligonucleotides may be used in a method of amplifying and detecting a target HPV16 sequence, such as a target HPV16 ctDNA sequence. The set of oligonucleotides comprises a forward primer having at least 90% sequence identity (e.g. at least 90%, at least 91%, at least 92%, at least 93%, at least 94%, at least 95%, at least 96%, at least 97%, at least 98%, at least 99%, or 100%) to SEQ ID NO: 1, a reverse primer having at least 90% sequence identity (e.g at least 90%, at least 91%, at least 92%, at least 93%, at least 94%, at least 95%, at least 96%, at least 97%, at least 98%, at least 99%, or 100%) to SEQ ID NO: 2, and a probe having at least 90% sequence identity (e.g. at least 90%, at least 91%, at least 92%, at least 93%, at least 94%, at least 95%, at least 96%, at least 97%, at least 98%, at least 99%, or 100%) to SEQ ID NO: 3.
In some embodiments, the forward primer comprises the nucleotide sequence of SEQ ID NO: 1. In some embodiments, the reverse primer comprises the nucleotide sequence of SEQ ID NO: 2. In some embodiments, the probe comprises the nucleotide sequence of SEQ ID NO: 3. In some embodiments, the forward primer comprises the nucleotide sequence of SEQ ID NO: 1, the reverse primer comprises the nucleotide sequence of SEQ ID NO: 2, and the probe comprises the nucleotide sequence of SEQ ID NO: 3. In some embodiments, the probe comprises a detectable label. In some embodiments, the probe comprises a fluorescent label. Suitable fluorescent labels include, for example, FAM (5- or 6-carboxyfluorescein), VIC, NED, Fluorescein, FITC, IRD-700/800, CY3, CY5, CY3.5, CY5.5, HEX, TET, TAMRA, JOE, ROX, BODIPY TMR, Oregon Green, Rhodamine Green, Rhodamine Red, Texas Red, Yakima Yellow, Alexa Fluor PET, Biosearch Blue™, Marina Blue®, Bothell Blue®, Alexa Fluor®, 350 FAM™, SYBR® Green 1, Fluorescein, EvaGreen™, Alexa Fluor® 488 JOE™, VIC™ HEX™ TET™ CAL Fluor® Gold 540, Yakima Yellow®, ROX™, CAL Fluor® Red 610, Cy3.5™, Texas Red®, Alexa Fluor® 0.568 Cy5™, Quasar™ 670, LightCycler Red640®, Alexa Fluor 633 Quasar™ 705, LightCycler Red705®, Alexa Fluor® 680, SYTO® 9, LC Green®, LC Green® Plus+, and EvaGreen™
In some embodiments, the set of oligonucleotides generate an amplicon comprising SEQ ID NO: 4 or a portion thereof comprising at least 15 contiguous bases of SEQ ID NO:4. For example, in some embodiments the set of oligonucleotides comprise a forward primer and a reverse primer that generates an amplicon comprising SEQ ID NO: 4 or a portion thereof, and a probe that hybridizes to the amplicon.
In some aspects, provided herein are kits. The kits may comprise a set of oligonucleotides, as described herein. The kits may further comprise reagents necessary to purify, isolate, detect and/or quantify HPV16 ctDNA. For example, the kits may further comprise PCR reagents, including buffers and enzymes (e.g. dUTPs, dNTPs). The kit may further comprise control samples, if needed. The kit may comprise solid surfaces (e.g., beads) comprising capture reagents (e.g., oligonucleotides) specific for target ctDNA. The kit can further include containers for holding or storing a sample (e.g., a container or cartridge for a plasma sample, a container for a cell free DNA sample, etc.). The kit can also include one or more instrument for assisting with obtaining a test sample, such as a syringe. Where appropriate, the kit optionally also can contain reaction vessels, mixing vessels, and other components that facilitate the preparation of reagents (e.g. a container for mixing reagents for PCR).
The kit may further comprise instructions for use of the kit. Instructions included in kits can be affixed to packaging material or can be included as a package insert, or can be viewed or downloaded from a particular website that is recited as part of the kit packaging or inserted materials. While the instructions are typically written or printed materials they are not limited to such. Any medium capable of storing such instructions and communicating them to an end user is contemplated by this disclosure. Such media include, but are not limited to, electronic storage media (e.g., magnetic discs, tapes, cartridges, chips), optical media (e.g., CD ROM), and the like. As used herein, the term “instructions” can include the address of an internet site that provides the instructions.
HPV16 ctDNA Assay Optimization
To develop a novel HPV16 ctDNA assay, regions of the HPV16 genome that are retained throughout the course of HPV+HNSCC cancer progression were characterized. A detailed characterization of the complete HPV genome in tumor metastases using a focused enrichment analysis has not previously been accomplished using sequencing-based approaches. Therefore, to provide comprehensive empirical evidence for regions of the HPV genome that are retained in metastatic lesions and could be the focus of the assay, targeted capture based sequencing of biopsy specimens collected from a prospective patient cohort was performed (
In parallel, a bioinformatics-based approach was implemented to define local specificity of potential amplicon regions by comparing the HPV16 and HPV18 genomes. These two genomes were selected as they represent the highest and second highest HPV types in HPV+ cancers,31 respectively. An 80 nucleotide sliding window was used to approximate the size of a small primer-probe amplicon-based assay (2×25 nt primers+30 nt probe), given that most abundant cell free DNA expected to be accessible to standard ddPCR assay-design is in the range of 73-165 nts.32 This demonstrated that the HPV16 and HPV18 genomes have minimal pairwise alignment at 80 nucleotide resolution with a few exceptions (
For optimal primer and probe design, the assay was focused on a region of the HPV16 genome with GC-content around 45%, based on data showing the impact of GC-content on polymerase amplification ability of small amplicons.33 Again using an 80 nt sliding window, three potential regions of the focus region for assay design were identified (
Analytic Validation of the HPV16 ctDNA ddPCR Assay
To determine the limit of detection for the assay, serial dilutions of synthetic HPV16 DNA containing the 77 bp amplicon specific to assay V9 were used (
To then confirm specificity of the ddPCR assay, assay results on genomic DNA isolated from a cohort p16+OPSCC tumors were compared to a series of melanoma tumor tissue specimens as HPV-negative controls. This confirmed the high specificity of the ddPCR assay for HPV16+ tumor tissue (
Clinical Performance of HPV ctDNA Assay
Sixteen patients with p16 positive tumors were enrolled between October 2017 and April 2019 in which 102 distinct plasma sample were collected. The primers described herein were able to detect HPV16 ctDNA in baseline samples from 12 of these patients (12/16, 75%), all of whom had HPV type 16. Baseline characteristics of the 12 subjects are summarized in Table 1. Blood was collected with Paxgene tubes for all time points in 10 patients and Streck Cell-Free BCT tubes for all time points in 2 patients. Of these 12 subjects, 6 had ctDNA levels analyzed during two distinct treatment courses representing a total of 18 distinct treatment courses (85 distinct samples). There was an average of 5 samples collected per treatment regimen (range: 2-17). Other than one patient missed plasma collection on Cycle 2, Day 1, there were no missing samples. Sample collection is outlined in
16 of 12 subjects had ctDNA collected during two distinct treatment courses representing a total of 18 distinct treatment courses
2Performed on Cycle 3, Day 1
3Performed on Cycle 4, Day 1
Change in Plasma HPV16 ctDNA Levels Correlates with Radiographically-Determined Treatment Response
The concordance of detected radiographic progression was assessed at the time of standard of care re-imaging with an increase in HPV16 ctDNA observed in blood collected synchronously at the time of re-imaging. Thus, the change in HPV16 ctDNA from baseline to the time of restaging imaging at post-cycle 2 therapy or post-cycle 3 therapy time points was evaluated (
Next, the longitudinal pattern of HPV16 ctDNA within each patient in relationship to their clinical history including treatment modality, symptoms of progression, and imaging findings was evaluated. A sample of patient HPV16 ctDNA levels and histories over time are shown in
Given these observations, it was assessed if early changes in HPV16 ctDNA were predictive of future radiographic response to therapy. To test this, the change between baseline HPV16 ctDNA and that drawn prior to cycle 2 was investigated (
Described herein is the development and detailed analytical validation of a high-precision ddPCR assay to quantify plasma HPV16 ctDNA, and the demonstration through a prospective study that results of radiographic assessment and simultaneously collected HPV16 ctDNA were highly concordant. Furthermore, analysis suggests change in HPV16 ctDNA after just one cycle of treatment is predictive of radiographic progression. Specifically, the data suggests that a <60% increase (or any decrease) in HPV16 ctDNA is associated with a favorable response to therapy (complete response, partial response, or stable disease), whereas patients with ≥60% increase in HPV16 ctDNA have disease progression and presumably do not derive any benefit from therapy. Finally, these data demonstrate that the percent change in HPV16 ctDNA after a single cycle of treatment correlates with that obtained synchronous to re-staging imaging.
Despite introduction of immunotherapeutic agents, survival for patients with R/M HNSCC remains poor.6 Initial response to therapy can have a significant favorable impact on survival for some patients,37 however a clear understanding of which patients will respond to each therapy is lacking. Additionally, futile systemic treatments can have significant toxic side effects. Currently established tissue biomarkers, including PD-L1, could theoretically be used to guide initial treatment selection based on historical outcomes but require invasive biopsy. In OPSCC, biopsies are often challenging due to the nature of the anatomy involved. Circulating biomarkers, on the other hand, allow for dynamic quantification of real-time treatment response obtained via minimally invasive approaches. The use of plasma ctDNA to guide treatment decisions is increasing in the oncology clinical setting.38 For example, PCR-based ctDNA assays for EGFR genotyping in non-small cell lung cancer, and for KRAS genotyping in colorectal cancer have demonstrated clinical validity and thus have received regulatory approval in the United States and Europe.39-41 Little investigation has been performed to date in R/M HPV+OPSCC. This work demonstrated that patients with poor outcomes had higher plasma levels of ctDNA and that patients who responded to treatment were seen to have decrease in plasma HPV16 ctDNA levels whereas patients with progressive disease had increases in plasma HPV16 ctDNA.
The data demonstrate a role for HPV16 ctDNA for early detection of progressive disease in R/M HPV+OPSCC. Accordingly, this assay finds use in guiding treatment decisions regarding efficacy of therapy, thereby potentially switching patients from futile therapy to more efficacious therapy more promptly. Furthermore, HPV16 ctDNA serves as a tool to distinguish clinical benefit and progressive disease on imaging in patients treated with immunotherapy given the phenomenon of pseudo-progression. The assay herein allows for a paradigm-defining approach to head and neck cancer in which systemic therapy is switched after one cycle based on biomarker indications of treatment failure, rather than delaying that intervention until radiographic imaging demonstrates disease progression. This data demonstrates that changes in HPV16 ctDNA is predictive of progressive disease in patients with R/M HNSCC receiving systemic therapy. Changes in HPV16 ctDNA precede radiographic response and thus find use as an early predictive biomarker to guide treatment decisions, improving survival and sparing toxicity.
Development of HPV ctDNA Assay
Droplet Digital PCR (ddPCR) Assay for HPV16 ctDNA
A Taqman probe-based ddPCR assay for the most common HPV subtype, HPV16, was developed as described below. Each ddPCR sample reaction contained 7 uL of DNA with 1 uL HPV target assay mix (forward and reverse primers and TaqMan FAM probe), 1 uL RPP30 assay mix (with HEX probe, Assay ID: dHsaCP2500350), 1 uL nuclease free water, and 10 uL ddPCR Super Mix (no dUTP) (#186-3024, Bio-Rad, Hercules, CA). A ddPCR assay for the RPP30 control gene was used to assess sample quality (Bio-Rad, Hercules, CA). Each sample assay was run in duplicate. UM-SCC-104 (HPV16 positive) and UM-SCC-105 (HPV18 positive) cell line DNA were run to serve as positive and negative controls.
Droplet generation was performed using the QX200 Droplet Digital PCR System (Bio-Rad, Hercules, CA). PCR was run with the following thermocycler conditions: enzyme activation for 10 minutes at 95° C. for 1 cycle, denaturation at 94° C. for 30 seconds and annealing/extension at 60° C. for 60 seconds for 40 cycles at a 2° C./second ramp rate, enzyme deactivation for 10 minutes at 98° C. for 1 cycle, followed by a hold at 4° C. Droplets were read using QuantaSoft Software (Bio-Rad, Hercules, CA). Standardized thresholds for positive and negative droplets were used across samples. Samples were considered positive if they had 2 or more positive droplets at the same amplitude as the positive cell line control. The Bio-Rad system quantifies DNA in copies per uL of the 20 uL PCR reaction, which are the values reported here.
To determine the Limit of detection (LoD), a synthetic 87 bp dsDNA fragment containing the 77 bp HPV16 E6 amplicon from assay V9 (5′ AATGCCGAAACCGGTTAGTATAAAAGC AGACATTTTATGCACCAAAAGAGAACTGCAATGTTTCAGGACCCACAGGAGCGACT CACT 3′) or genomic DNA from UM-SCC-104 cell line was tested in a background human genomic DNA matrix (200,000 genome equivalents (GEs) per 20 ul ddPCR reaction). The 77 bp amplicon (SEQ ID NO: 4) is underlined above. In order to generate background matrix DNA, human genomic DNA (#PR-G3041, Promega) was incubated with 1 unit HindIII enzyme (#R0104S, New England Biolabs) per microgram of DNA at 37° C. for 4 hours followed by enzyme inactivation at 80° C. for 20 min. A 2-fold dilution series (2000 copies or GEs per ul to 4 copies or GEs per ul) was tested in triplicates using ddPCR and LoD was calculated using the formula:
LoB was determined per 20 uL reaction using 200,000 diploid human GEs (HindIII digested) as a non-HPV background matrix (n=15). With UM-SCC-104 cell line, LoD=0.3 copies/20 ul reaction based on standard deviation of 0.2 for the lowest concentration sample (4 GEs/ul); and with synthetic DNA (containing the 77 bp amplicon sequence), LoD=4.2 copies/20 ul reaction based on standard deviation of 2.54 for the lowest concentration sample (4 copies/ul). The Coefficient of Variation (CV) was calculated using the formula % CV=(Standard Deviation/Mean)*100.
Characterizing the Precision of the HPV16 ctDNA ddPCR Assay
One HPV16 ctDNA sample with low counts and one sample with high counts were selected to run reproducibility testing. ctDNA was extracted from 3×2 mL aliquots and each diluted 1:40 as described below. From this, 12 reaction assays were created for each sample. Droplet generation, PCR amplification, and droplet reading were performed as indicated above. HPV16 ctDNA counts were plotted. Means with standard deviations were calculated for the low and high sample. Standard curves were created to determine expected standard deviations for any given HPV16 ctDNA count.
The genomic sequence for HPV16 (K02718.1) and HPV18 (AY262282.1) were downloaded from Genbank. One nucleotide-offset sliding windows of 18, 24, 34 and 77 nt were extracted from the HPV16 reference sequence, corresponding to the lengths of the reverse primer, forward primer, probe and amplicon sequence, respectively. Pairwise alignment scores were calculated between each of the HPV16 k-mer window sequences and the HPV18 genome using the pairwise alignment function (R package, Biostrings v2.52.0) with default settings for the “global-local” alignment type and scoreOnly=TRUE, which returns the maximum alignment score for each k-mer window sequence. The maximum alignment score for each k-mer window was plotted based on the start position of each window. Using the same settings, pairwise alignment scores were also calculated between the HPV18 genome and the HPV16 primer, probe and amplicon sequences, with 0-4 mutations (substituting “N” bases) introduced at random locations. Plots were generated using the R packages, ggplot2 (v3.2.1), cowplot (v1.0.0) and ggforce (v0.3.2).
A prospective cohort was enrolled as part of a University of Michigan Institutional Review Board-approved, observational longitudinal biospecimen collection study (UM IRB 00042189). All patients provided written informed consent. Patients ≥18 years old with histologically documented HPV+OPSCC undergoing systemic therapy were eligible. p16 was determined by tumor immunohistochemistry and was used as a surrogate marker for HPV status of a patient's cancer. All patients were required to have the presence of measurable disease by CT scan, or cutaneous lesions ≥10 mm not assessable on imaging but present on physical exam. Adequate hematopoietic and renal function were required and defined as hemoglobin ≥9.0 g/dL and serum.
Enrolled patients had baseline collection of clinical data on demographics, disease characteristics, treatment plan, laboratory studies (CBC with differential, comprehensive metabolic profile) and radiographic staging studies. Pre-treatment tumor tissue was obtained for next generation sequencing, whenever available. Duration of a treatment cycle was defined based on clinical documentation. Blood was collected at each visit during the treatment course as well as synchronous with radiographic imaging until disease progression, patient withdrawal of consent, or at investigator discretion. 10 mL of blood was collected from patients in sterile PreAnalytix Paxgene tubes (PreAnalytiX GmbH, Switzerland) or Streck Cell-Free DNA BCT tubes (Streck Inc., Omaha, NE, USA) with the same type of tube used for serial collection within a given participant. Previous publications of ctDNA analysis have demonstrated equivalence between these collection tubes.27 Tubes were held at room temperature and processed as per manufacturer's recommendations within 7 days of collections as described in detail below. Re-staging imaging (after 2 or 3 cycles) was obtained at the frequency as determined by the discretion of the treating physician. Timing of the imaging and corresponding results were recorded in the database. Radiologic response was assessed according to RECIST v1.1. In the case of treatment continued despite RECIST-assessed progressive disease, due to suspected pseudoprogression, response was confirmed on subsequent scans. If a patient experienced progression and was started on a subsequent therapy, they were re-enrolled with identical approach of data and specimen collection.
DNA Isolation from Tumor Tissue
Two tumor tissue blocks from patients with R/M HPV+OPSCC from the prospective cohort were available with sufficient material for NGS and eight blocks from primary HPV+OPSCC and 7 with melanoma blocks were available with sufficient material for ddPCR based tumor analysis. DNA extraction from FFPE specimens was performed as previously described.[1, 2] Areas of tumor were identified by a board-certified head and neck pathologist (J.B.M). Tumor and adjacent normal tissue cores were collected and genomic DNA was obtained using a Qiagen Allprep DNA/RNA FFPE kit (Qiagen, Hilden, Germany) and quantified using Qubit as previously described.[2] For ddPCR assays, RPP30 Bio-Rad assay was used to measure RPP30 reference gene and ensure ample DNA being assessed in the case of HPV positive and HPV negative tumors (data not shown).
DNA from tumor tissue was submitted to the University of Michigan Advanced Genomics Core for targeted capture sequencing using the DNA Thruplex kit (Takara Biosciences). Targeted capture was performed using a custom designed probe panel from Nextera with high density coverage of HPV16 which is described in O'Leary et. al 2020.[3] Following library preparation and capture, the samples were sequenced on an Illumina HiSEQ4000, respectively, with paired-end 150 nt reads. Data was de-multiplexed and FastQ files were generated. Samtools depth (samtools/1.9) was used to compute the read depth at each position across the HPV16 genome. The input is a BAM file derived from HPViewer with short reads aligned to the HPV genome. The unmapped, secondary, QC-fail and duplicate reads were skipped from the analysis. Only reads with base quality greater than 5 and mapping quality greater than 20 were counted, plots were created with Excel.
DNA Extraction from Plasma
Plasma was isolated in two sequential centrifugation steps per manufacturer recommendations. Briefly, PaxGene tubes were centrifuged at 1900×g for 15 minutes; plasma layer was transferred into a 15 ml conical tube and centrifuged a second time again at 1900×g for 10 minutes. Streck tubes were centrifuged at 300×g for 20 minutes; plasma layer was transferred into a 15 mL conical tube and centrifuged a second time at 4600×g (maximum possible speed) for 10 minutes. Cell free DNA was isolated from a 2 mL aliquot of plasma using QIAamp® MiniElute® ccfDNA Mini Kit (Qiagen, Germantown, MD) and diluted 1:40 according to manufacturer instructions. Genomic DNA from cell lines UM-SCC-104 (HPV16 positive), UM-SCC-105 (HPV18 positive), and UM-SCC-85 (HPV negative)28 was extracted using the Promega Wizard DNA Purification Kit protocol (Promega A1120, Madison WI). Cell line DNA was used for assay optimization and as positive controls for plasma ctDNA assays.
The primary aim of the prospective pilot cohort was to assess the HPV16 ctDNA assay in a diverse patient population and assess for 1) concordance with imaging results and 2) potential signal of association of early dynamics with radiographic progression. Treatment responses, as determined by radiology, after 2-3 cycles of treatment, based on treatment protocol, were collected from the medical record. For each patient, HPV16 ctDNA counts/mL were plotted over time along with treatment histories. Percent change in HPV16 ctDNA level was calculated relative to a patient's own baseline ((time point level-baseline level)/baseline level) and assessed at two clinical episodes. The HPV16 ctDNA level closest to the radiographic assessment was chosen to test association of HPV16 ctDNA with synchronous imaging results whereas the level on day 1 of cycle 2 of therapy was chosen to explore the ability of HPV16 ctDNA to predict radiographic progression before the standard timeline of assessment. Percent change in HPV16 ctDNA across radiographic response groups were compared with nonparametric Wilcoxon, Kruskal-Wallis, and median one-way analysis tests. ROC curve analysis was performed for the prediction of progressive disease at time of imaging. From the ROC curve, area under the curves were extracted (AUC) as well as the optimal cutpoint using Youden index to maximize the sum of sensitivity and specificity for predicting progression. The optimal cutpoint was estimated from the full sample and in n=500 bootstrapped samples. Statistical analysis was performed using SPSS (v27) and R packages, ggplot2 (v3.2.1), cutpointR (v1.1.0).
For HPV+OPSCC, a large number of clinical trials are attempting to advance new treatment de-escalation strategies that maintain excellent survival outcomes while minimizing long-term toxicities.2-6 However, these trials carry the risk of de-escalating patients that would have otherwise benefited from standard treatment; and, the identification of such patients, with ctDNA or other biomarkers, remains of paramount importance. Retrospective data with small patient cohorts suggest that HPV ctDNA is difficult to detect in patients with early stage disease (26% of patients with T1, 48% of patients with T2 disease)7,8, levels correlate with disease burden,9,10 rapid clearance after chemoradiation (CRT) is associated with improved locoregional control11, and levels increase at the time of recurrence9,12. Importantly, very limited data has been published from low risk OPSCC patients and the published studies have few failures, preventing the rational design of HPV ctDNA trials in this setting.
A phase II clinical trial of low-risk HPV+OPSCC was completed to determine whether treatment stratification by neck dissection, to more accurately pathologically stage patients, could minimize treatment modalities and thereby improve quality of life. The trial cohort therefore provided the unique opportunity to prospectively evaluate HPV ctDNA clearance following various modes of therapy in low-risk disease and to test the utility of HPV ctDNA for early detection of recurrence during surveillance.
Twenty-seven T1 and T2 patients with limited neck disease had plasma drawn before treatment and, 35.7%(5/14) of T1, and 69%(9/13) of T2 patients had detectable HPV ctDNA at baseline. Baseline HPV ctDNA levels were significantly higher in patients with T2 disease compared with T1 disease (averages of 2,467 vs 130.7 copies per milliliter plasma, respectively, P<0.05). There were no differences in baseline plasma HPV ctDNA counts based on N classification or presence of extracapsular extension. All patients had complete clearance of plasma HPV ctDNA at 4 weeks post treatment regardless of modality (
In total, 81 post-treatment plasma specimens from patients in this cohort were evaluated with a median follow up duration of 33.3 months (range: 0.56-39.6 months). One patient had detectable ctDNA at baseline that became undetectable 4 weeks post-surgery. Then, HPV ctDNA was detected six months post-treatment and persisted through the next 6 time points (
Ongoing de-escalation strategies use various modalities to allow for de-intensified adjuvant treatment3-6. While this trial utilized the pathologic findings from neck dissection as a clinical biomarker to guide primary treatment modality, leveraging the cohort to gain insight into the clearance of HPV ctDNA in early stage disease was possible. The data demonstrate that all patients had undetectable plasma HPV ctDNA within one month of completion of treatment regardless of treatment modality and future recurrence. However, one of these patients developed recurrent disease, in spite of undetectable HPV ctDNA in plasma. These data are consistent with current working mathematical models of ctDNA shedding,13 which suggest that tumor size, death and other shedding characteristics may limit the amount of ctDNA in circulation to below detectable levels with current technology. Despite this limitation, recent published data in a variety of settings suggest that plasma HPV ctDNA may be a complementary biomarker for identification of residual or recurrent disease.7-9,11,12 Herein, data is provided showing that HPV ctDNA was detectable 18 months prior to overt clinical or radiographic recurrence.
Thirty-four patients with low-risk HPV+OPSCC (T1-3, N0-N1) were enrolled (NCT02784288). Patients underwent upfront neck dissection with stratification into one of three treatment groups based on post-operative pathology. Patients with a single lymph node <6 cm with no extracapsular extension (ECE), perineural (PNI), or perivascular invasion (PVI) underwent transoral surgery. Patients with two or more positive nodes without adverse features (ECE or positive margins) or a single node in which PVI or PNI was noted were treated with radiotherapy (RT) alone. Patients with ECE underwent CRT. Plasma was isolated from blood collected in Streck tubes from 27 patients prior to the neck dissection, one-month post-treatment (surgery, RT, or CRT), and every 3 months thereafter.
CfDNA was isolated following the QIAamp MinElute ccfDNA Mini Kit protocol. A digital PCR ctDNA assay was performed to amplify a 77 bp region of HPV16_E6 with forward 5′-CGAAACCGGTTAGTATAAAAGCAG-3′, reverse 5′GTCGCTCCTGTGGGTCCT-3′ and probe 5′-FAM-CATTTTATGCACCAAAAGAGAACTGCAATGTTTC-MGBNFQ-3′. Each reaction assay contained 10 μL of 2×dPCR Supermix for Probes (No dUTP), 0.9 μmol/L of respective primers, 0.25 μmol/L of respective probes, and 5 to 50 ng of cfDNA in a final volume of 20 μL. dPCR assays were performed on the QX100 and/or QX200 platforms outfitted with an automated droplet generator (Bio-Rad). An RPP30-HEX reference was used to assess quality (Assay-ID: dHsaCP2500350, Bio-Rad).
9 Hanna G J, Supplee J G, Kuang Y, Mahmood U, Lau C J, Haddad R I, et al. Plasma HPV cell-free DNA monitoring in advanced HPV-associated oropharyngeal cancer. Ann Oncol 2018; 29:1980-6. https://doi.org/10.1093/annonc/mdy251.
It is understood that the foregoing detailed description and accompanying examples are merely illustrative and are not to be taken as limitations upon the scope of the disclosure, which is defined solely by the appended claims and their equivalents.
Various changes and modifications to the disclosed embodiments will be apparent to those skilled in the art.
Any patents and publications referenced herein are herein incorporated by reference in their entireties.
This application claims priority to U.S. Provisional Patent Application No. 63/208,736, filed Jun. 9, 2021, the entire contents of which are incorporated herein by reference for all purposes.
| Filing Document | Filing Date | Country | Kind |
|---|---|---|---|
| PCT/US2022/032766 | 6/9/2022 | WO |
| Number | Date | Country | |
|---|---|---|---|
| 63208736 | Jun 2021 | US |
