Next Article in Journal
What Do We Know about Survival in Skeletally Premature Children Aged 0 to 10 Years with Ewing Sarcoma? A Multicenter 10-Year Follow-Up Study in 60 Patients
Previous Article in Journal
Impact of Cancer Stem Cells on Therapy Resistance in Gastric Cancer
Previous Article in Special Issue
Updated Prognostic Factors in Localized NSCLC
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Multiplex Droplet Digital PCR Assay for Detection of MET and HER2 Genes Amplification in Non-Small Cell Lung Cancer

by
Igor P. Oscorbin
1,*,
Maria A. Smertina
1,
Ksenia A. Pronyaeva
1,2,
Mikhail E. Voskoboev
1,3,
Ulyana A. Boyarskikh
1,
Andrey A. Kechin
1,
Irina A. Demidova
4 and
Maxim L. Filipenko
1
1
Institute of Chemical Biology and Fundamental Medicine, Siberian Branch of the Russian Academy of Sciences (ICBFM SB RAS), 630090 Novosibirsk, Russia
2
V. Zelman Institute for Medicine and Psychology, Novosibirsk State University, 630090 Novosibirsk, Russia
3
Department of Natural Sciences, Novosibirsk State University, 630090 Novosibirsk, Russia
4
Moscow City Oncology Hospital, 143423 Stepanovskoe, Russia
*
Author to whom correspondence should be addressed.
Cancers 2022, 14(6), 1458; https://doi.org/10.3390/cancers14061458
Submission received: 31 January 2022 / Revised: 8 March 2022 / Accepted: 9 March 2022 / Published: 11 March 2022
(This article belongs to the Special Issue Advances in Prognosis and Theranostics of Lung Cancer)

Abstract

:

Simple Summary

Non-small-cell lung cancer (NSCLC) remains one of the most common tumors with a high mortality and morbidity rate. Alterations in HER2 and MET could be a target for anti-tumor drugs or lead to resistance to anti-EGFR therapeutics. Existing methods for detecting HER2 and MET amplifications are time and labor-consuming, and alternative methods are needed. We report the first multiplex droplet digital PCR (ddPCR) assay for the simultaneous detection of MET and HER2 amplification in NSCLC samples. The suitability of qPCR was assessed for the optimization of multiplex ddPCR, and optimal ddPCR conditions were selected. The developed ddPCR was validated on artificial samples with various DNA concentrations and MET and HER2 ratios. Using ddPCR, 436 EGFR-negative NSCLC samples were analyzed, and, among them, five specimens (1.15%) were MET-positive and six samples (1.38%) were HER2-positive. The multiplex ddPCR assay could be used for screening MET and HER2 amplification in NSCLC samples.

Abstract

Non-small-cell lung cancer (NSCLC), a subtype of lung cancer, remains one of the most common tumors with a high mortality and morbidity rate. Numerous targeted drugs were implemented or are now developed for the treatment of NSCLC. Two genes, HER2 and MET, are among targets for these specific therapeutic agents. Alterations in HER2 and MET could lead to primary or acquired resistance to commonly used anti-EGFR drugs. Using current methods for detecting HER2 and MET amplifications is time and labor-consuming; alternative methods are required for HER2 and MET testing. We developed the first multiplex droplet digital PCR assay for the simultaneous detection of MET and HER2 amplification in NSCLC samples. The suitability of qPCR was assessed for the optimization of multiplex ddPCR. The optimal elongation temperature, reference genes for DNA quantification, and amplicon length were selected. The developed ddPCR was validated on control samples with various DNA concentrations and ratios of MET and HER2 genes. Using ddPCR, 436 EGFR-negative NSCLC samples were analyzed. Among the tested samples, five specimens (1.15%) showed a higher ratio of MET, and six samples (1.38%) showed a higher ratio of HER2. The reported multiplex ddPCR assay could be used for the routine screening of MET and HER2 amplification in NSCLC samples.

1. Introduction

Currently, lung cancer remains one of the most common malignancies. In 2020, lung cancer became the second most frequently diagnosed cancer, with 2.2 million estimated new cases (11.4% of all new cancer cases), and remained a leading reason for cancer-associated mortality, with an estimated 1.8 million deaths [1]. Thus, developing new lung cancer diagnostics and treatment approaches is a pressing challenge that can hardly be overestimated. Based on cell morphology and immunohistochemistry, lung cancer is divided into several subtypes: small cell lung cancer (15% of all cases) and non-small-cell lung cancer (NSCLC). NSCLC is further classified into the most common adenocarcinomas (63% of NSCLC), large-cell carcinomas (7% of NSCLC), and squamous cell carcinomas (18% of NSCLC) [2,3]. Over the past few decades, extensive studies revealed the genetic landscape of NSCLC, which allowed the identification of promising targets for therapeutic agents. Tyrosine kinase inhibitors against activated EGFR (epidermal growth factor receptor) [4] and ALK [5] are the most well-known such targeted drugs, widely used in routine clinical practice. Recent advances in the molecular oncology of NSCLC have broadened the spectrum of novel possible targets. Several of these genes—ROS1, BRAF, NTRK1-NTRK3, MET, HER2, RET, and KRAS—are high-potential candidates with several developing or approved inhibitors [6]. Thus, in 2021, FDA (U.S. Food and Drug Administration) approved sotorasib and amivantanab for the treatment of tumors harboring KRAS G12C and 20 exon insertions in EGFR genes, respectively [7,8]. It is believed that the list of available target drugs will be expanded more in the coming years against other mentioned gene products.
Among novel targets for inhibitors, MET (mesenchymal epithelial transition) and HER2 (human epidermal growth factor receptor 2) have been extensively studied. MET and HER2 alterations are also a reason for the resistance to commonly used anti-EGFR therapy. MET is an oncogene that encodes a receptor tyrosine kinase for hepatocyte growth factor (HGF, also known as scatter factor). The binding of HGF to the MET receptor activates downstream signaling pathways, such as RAS/ERK/MAPK and PI3K/AKT, and leads to increased cell motility and survival and enhanced angiogenesis [9,10]. MET overexpression is frequent in NSCLC (35–72%), whereas MET exon 14 skipping and MET amplification are relatively rare events (3–5% and 1–5%, respectively) [11]. The FDA has approved tyrosine kinase inhibitor crizotinib for NSCLC patients with MET exon 14 mutations [12]. For patients with MET-amplified NSCLC, crizotinib, capmatinib, and Sym015-01 are undergoing clinical trials [6]. Capmatinib and Sym015-01 showed potential as a cure for patients with NSCLC harboring MET amplification.
HER2 is an oncogene encoding receptor tyrosine kinase ERBB2 that has no directly activating ligand. The constitutive activation of HER2 promotes cell growth, proliferation, and survival [13]. HER2 is a well-known target for the treatment of breast cancer, with a long history of application of the respective anti-HER2 drugs [14,15,16,17]. HER2 overexpression is found in 2.8–38% of NSCLC cases, with a 1.5–4% frequency of mutations and 1.1–3% frequency of amplification [18,19]. Recently several anti-HER2 drugs have been assessed for patients with NSCLC [6]; the anti-HER2 antibody–drug conjugate trastuzumab–deruxtecan received FDA approval for NSCLC harboring HER2 mutations [20]. Other anti-HER2 therapeutics are undergoing clinical trials.
Normally, for the detection of gene copy amplification events, academicians and clinical pathologists apply immunohistochemistry (IHC), fluorescence in situ hybridization (FISH), next generations sequencing (NGS), and multiplex ligation-dependent probe amplification (MLPA). However, all of these methods have limitations. IHC scoring could be subjective, FISH requires high-quality FFPE samples for testing, and NGS is amenable to contamination. Both techniques, as well as MLPA, are labor and time-consuming, demanding to the level of personnel training, relatively costly, and require a high input of the starting material for the analysis. The latter should be specifically noted; molecular pathologists often receive only biopsy samples from NSCLC patients, as no surgery is available. In addition, liquid biopsy, which is testing circulating tumor DNA in blood plasma, has recently been widely used in clinical practice. Liquid biopsy allows for the real-time monitoring of the appearance of molecular alterations leading to acquired resistance to target therapy. However, such testing requires sensitive and fast approaches.
Recently, various digital PCR techniques have been developed and used for numerous tasks, from somatic mutations and genetic modifications detection to the quantifying of viral genomic material. One of the possible applications of digital PCR is the testing of gene copy number events. Thus, numerous digital PCR-based approaches were reported for gene amplification detection for MET [21,22,23,24] and HER2 [25,26,27,28,29,30] genes. However, most of these techniques are monoplex or duplex and require several separate reactions to analyze several genes. Different approaches for multiplexing digital PCR have been reported to address that problem based on the amplitude of end-point fluorescence or mixing of probes with different fluorophores. To date, no multiplex digital PCR has been published for MET and HER2 gene amplification.
We report the first tetraplex ddPCR-based assay to detect MET and HER2 gene amplification. The developed assay utilizes the amplitude of end-point fluorescence for multiplexing MET and HER2 target genes and RPP30 and RAD51AP2 references genes. After optimization, we validated the ddPCR test on 436 NSCLC samples and proved its suitability for practical usage.

2. Materials and Methods

2.1. Clinical Samples and DNA Extraction

DNA was extracted from FFPE tumor tissue sections using QIAamp DNA FFPE Tissue Kit (Qiagen, Hilden, Germany) following the manufacturer’s protocol. FFPE sections were obtained from 436 NSCLC (lung adenocarcinoma) patients who have been operated on in 14 regional cancer centers across the Siberian and the Far Eastern Federal Districts of Russia.

2.2. Detection of Common EGFR Gene Mutations

Common EGFR mutations were analyzed using the cobas® EGFR Mutation Test v2, Version 1 (Roche, Basel, Switzerland) according to the manufacturer’s instructions. The amplification procedure was conducted using cobas Z480 thermal cycler (Roche, Basel, Switzerland).

2.3. Positive Controls

Plasmids carrying MET and HER2 genes fragments served as positive controls and were used to assess method sensitivity. All control plasmids used in this study were made in Shanghai RealGene Bio-tech, Inc. (Shanghai, China), purified, and sequenced (SB RAS Genomics Core Facility, ICBFM SB RAS, Novosibirsk, Russia). The purified plasmids were linearized using BamHI restriction endonuclease (SibEnzyme, Novosibirsk, Russia) and quantified using NanoDrop Lite A4 spectrophotometer (Thermo Fisher Scientific, Waltham, MA, USA) and droplet digital PCR.

2.4. Fluorescence In Situ Hybridization

Fluorescence in situ hybridization (FISH) was performed using a commercially available assay according to manufacturer recommendations. Briefly, 3 μm tissue section was hybridized overnight with the IVD MET (7q31)/SE7 (D7Z1) probe (Kreatech, Leica Biosystems, Lincolnshire, IL, USA). Normal tissue, including vessels, fibroblasts, lymphocytes, or non-neoplastic lung tissue, served as internal controls. Tumor tissue was entirely scanned for amplification hot spots by using a ×63 objective and appropriate filter sets (ZeissAxioimagerZ2 fluorescent microscope; Carl Zeiss, Oberkochen, Germany). Fifty cancer cells in four random areas were scored and interpreted using criteria reported by Camidge et al. [31].

2.5. Quantitative PCR

Quantitative PCR reactions were performed in 20 µL volume containing 65 mM Tris-HCl, pH 8.9, 24 mM (NH4) 2SO4, 0.05% Tween-20, 3 mM MgSO4, 0.2 mM dNTPs, 900 nM primers, 250 nM TaqMan probe (Table 1), and 1 U of Taq-polymerase. Amplification was carried out in CFX96 Real-Time PCR Detection System (Bio-Rad, Hercules, CA, USA) according to the following program: 95 °C for 15 min followed by 45 cycles of 95 °C for 10 s, and 60 °C for 40 s with a collection of fluorescent signals at FAM and HEX channels. In the case of the elongation temperature gradient, this temperature was varied at the range of 54–64 °C. The acquired data were analyzed with CFX Manager software (Bio-Rad, Hercules, CA, USA) (Figure 1).

2.6. Droplet Digital PCR

The ddPCR assay was performed using the QX200 system (Bio-Rad, Hercules, CA, USA) according to the manufacturer’s recommendations. The reaction mixture in a volume of 20 µL consisted of 1× ddPCR master mix (Bio-Rad, Hercules, CA, USA), primers, probes, and DNA template. For quantification of standard samples, reactions contained 0.9 µM Lac-U/R primers, 0.25 µM Lac-P probe (Table 1), and approximately 104–103 copies of the tested plasmid standard. For the MET and HER2 copy number analysis, concentrations of primers and probes were varied at the range of 450–900 and 125–250 nM, respectively; the exact concentrations are indicated in the respective sections. DNA amount of FFPE samples was in the range of 100–15000 copies per reaction (0.33–49.50 ng/reaction), and all specimens were analyzed in duplicates.
The entire reaction mixtures with 70 µL of droplet generation oil (Bio-Rad, Hercules, CA, USA) were loaded into a disposable plastic cartridge (Bio-Rad, Hercules, CA, USA) and placed in the droplet generator. After processing, the droplets obtained from each sample were transferred to a 96-well PCR plate (Eppendorf, Hamburg, Germany). The amplification was carried out using T100 Thermal Cycler (Bio-Rad, Hercules, CA, USA) according to the program: DNA polymerase activation at 95 °C for 10 min followed by 45 cycles of PCR amplification (94 °C for 30 s and 57 °C for 60 s), and 98 °C for 10 min, 2 °C/s ramp rate at all steps. After PCR, the droplets were counted with the QX200 Droplet Reader. The data obtained were analyzed with QuantaSoft Analysis Pro software (Bio-Rad, Hercules, CA, USA) and curated as explained in Section 3.2.3.

3. Results

3.1. Selection of References, Primers, and Probes

Here, we developed multiplex digital droplet PCR to test MET (7q31.2) and HER2 (17q12) amplification events in FFPE samples from patients with NSCLC. For that purpose, we decided to quantify concentrations of target genes MET and HER2 and two reference genes in a single tetraplex reaction. The multiplexing of ddPCR was based on the amplitude of end-point fluorescence. Thus, four independent amplification reactions were simultaneously performed in a single tube, two probes for two targets were FAM-labeled, and two probes were HEX-labeled. Probes with similar fluorophores were added to the reaction mixes at different concentrations, which resulted in a different end-point fluorescence for each target. Theoretically, up to 16 types of droplet clusters can be observed in such an assay, with various combinations of all four templates. The observed concentrations of targets allow for the calculation of a ratio between each target and the other three loci, and enable us to make a conclusion about the presence of MET and HER2 amplification in the sample.
As candidates for references, several genes were selected based on published works. The criteria for candidate genes were the reported absence of registered copy number amplifications events in cancer cells or previous usage for human genomic DNA quantification. Previously, EIF5B, DCK, RPS27A, and PMM1 genes were used as references to assess HER2, MET, and EGFR CNAs [24,30]. Here, we also tested reference genes on CNA in cancer cells using sequencing data from The Cancer Genome Atlas (TCGA) project. Several candidates met the mentioned requirements and were assessed for suitability in multiplex ddPCR: ALB (albumin, 4q13.3) [32,33], RPP30 (ribonuclease P protein subunit P30, 10q23.31) [33,34], RAD51AP2 (RAD51 associated protein 2, 2p24.2), TEKT3 (tektin 3, 17p12), and PARK7 (Parkinsonism associated deglycase, 1p36.23) [35].
Primers and probes for multiplex ddPCR were designed as described in our previous work [34]. Briefly, OligoAnalyzer software was used to select primers according to the following criteria: annealing temperature in the range of 60–62 °C for primers and 65–67 °C for probes, absence of G at the 5’-terminus of probes, the lengths of primers and probes in the range of 18–27 bases, Gibbs free energy ΔG of self-dimers ≥ −8 kcal/mol, and the amplicon length in the range of 61–100 bp. Sequences of selected oligonucleotides are listed in Table 1.

3.2. Optimization of Multiplex ddPCR Assay

3.2.1. Optimal Elongation Temperature of ddPCR

After the selection of reference genes and oligonucleotides, designed sets of primers and probes were assessed in both quantitative PCR and droplet digital PCR simultaneously. The main purpose of optimization was proper droplet clustering in ddPCR, where all clusters can be clearly distinguished and no “rain” (droplets with intermediate fluorescence between two clusters) is observed. Currently, qPCR allows us to obtain results faster and with a lower cost than ddPCR. Therefore, we decided to test all primers and probes in qPCR and ddPCR in parallel to determine the suitability of qPCR for ddPCR optimization.
At first, we determined the optimal elongation temperature for all primers and probes. For that, qPCR and ddPCR were performed with similar concentrations of primers and probes—900 and 250 nM, respectively—and DNA from white blood cells of a healthy donor was used as a template. The elongation temperature (Telongation) varied in the range of 54–64 °C. The results for qPCR are presented in Figure 1, and, for ddPCR, in Figure 2.
For all sets of primers and probes, the optimal Telongation for qPCR was in the range of 64-60 °C, as minimal Cq-values and maximal end-point fluorescence were observed at these temperatures, indicating the highest qPCR efficacy. At a lower Telongation, the Cq-value increased for most primers and probes. However, in ddPCR, for most primers and probes, the highest end-point fluorescence was observed at temperatures lower than 58 °C, and a further decrease in Telongation did not decrease the efficacy of ddPCR. No “rain” in ddPCR was observed at optimal temperatures for all primers and probes. Thus, the obtained results showed a discrepancy between the efficacy of qPCR and ddPCR at different Telongation.

3.2.2. Basic and End-Point Fluorescence in qPCR and ddPCR

Simultaneously with the optimal elongation temperature, we compared the basic and end-point fluorescence of primers and probes in qPCR and ddPCR. These two parameters define the suitability of monoplexes for multiplexing, as high, basic, or low end-point fluorescence would hinder the combining of primers and probes in amplitude-based multiplex ddPCR. The results for qPCR are presented in Figure 3a,c, and, for ddPCR, in Figure 3b,d.
The basic fluorescence was different in qPCR and ddPCR (Figure 3a,b). Thus, the fluorescence of PARK7 and RAD51AP2 was the highest in qPCR, whereas, in ddPCR, the fluorescence of TEKT3, HER2, and RPP30 was higher than PARK7 and RAD51AP2. The fluorescence of ALB in qPCR was equal to MET and RPP30, whereas, in ddPCR, it was two times lower compared with the same loci. Therefore, qPCR data regarding the basic level of fluorescence could not be extrapolated directly in ddPCR.
Unlike basic fluorescence, the amplitude of end-point fluorescence is more stable between qPCR and ddPCR (Figure 3c,d). In most cases, this parameter was independent of Telongation in qPCR, with the exception of PARK7 and HER2, where end-point fluorescence was the highest at 64 and 54 °C, respectively. A direct correlation was observed between end-point fluorescence and Telongation, as mentioned above. More importantly, the highest end-point fluorescence was almost the same for different loci, except for MET, whose end-point fluorescence was two to three times lower than for other loci. Notably, the ratio of end-point fluorescence of all loci was the same for qPCR and ddPCR, when MET was the lowest and HER2 was the highest. Taken together, these results indicate that qPCR data cannot be directly used for the prediction of basic fluorescence levels in ddPCR, but are suitable for the forecasting of end-point fluorescence.
Preliminary tests showed that all primers and probes were suitable for ddPCR and can be used in further experiments.

3.2.3. Reference Gene Selection

After the assessment of elongation temperature and end-point fluorescence, primers and probes for different loci were duplexed to select the most favorable combinations of target and references genes. All primers and probes with similar fluorophores were mixed in two possible combinations, where the primers/probe for the target and reference were added to final concentrations of 900/250 nM and 450/125 nM. The results of combined duplexes are presented in Figure 4.
The purpose of duplexing was to achieve clearly distinguishable clusters for different sets of primers and probes. Therefore, for HER2, the best results demonstrated a combination of standard concentration HER2 900 nM primers and a 250 nM probe with RAD51AP2 450 nM primers and a 125 nM probe. With other combinations, clusters for different loci were closer to each other. Based on the same criteria, the best combination with MET was 450 nM primers and a 125 nM probe MET, and 900 nM primers and a 250 nM probe RPP30. However, the end-point fluorescence of MET was low, and the 900 nM primers and 250 nM probe were used in further experiments to facilitate the detection of MET clusters.
Notably, the results of duplexing were in agreement with basic and end-point fluorescence levels of monoplexes. Thus, the low end-point fluorescence of MET led to a poor clustering of duplexes with MET, where the fluorescence amplitude was insufficient for proper cluster identification. However, the observed rain in combinations of HER2 and PARK7 could not be predicted in respective monoplexes, which underlines the need to optimize multiplex ddPCR without skipping steps.

3.2.4. Optimization of MET and HER2 Primers

At the last stage of optimization, we assessed the designed tetraplex ddPCR on a set of 46 DNA specimens purified from NSCLC FFPE tissue blocks. FFPE DNA samples are known to be highly fragmented and chemically modified, which leads to stalling DNA polymerase and a lower efficacy of PCR. Thus, the presence of rain, together with the actual DNA concentration for each target and reference loci, was appraised, and the ratio of concentrations between target and reference genes was calculated. The results of a pilot experiment with FFPE DNA samples are presented at Figure 5.
The designed tetraplex ddPCR showed clearly distinguishable clusters for all loci and was suitable for the analysis of DNA from FFPE samples; rain was moderate and did not hinder DNA quantification (Figure 5a). However, the observed ratio between the target genes (HER2, MET) and reference loci was close to 0.8, far from the expected ratio close to 1 (Figure 5b). The reason was that longer amplicons of MET and HER2 were 100 bp, whereas amplicons of RAD51AP2 and RPP30 were 71 and 64 bp, respectively. To address this emerging issue, primers for target genes were redesigned with lower lengths: 72 bp for HER2 and 77 bp for MET; the probes for target genes remained the same. New sets were validated in monoplexes and demonstrated the same efficacy of amplification (Figure 2). Thus, in further experiments, new primers for HER2 and MET were used, and a close-to-1 ratio between the reference and target genes was observed (Figure 5b). Thus, the issue of underrepresented HER2 and MET was successfully solved, and tetraplex ddPCR with shorter amplicons was used for further work.

3.3. Validation of Multiplex ddPCR Assay

3.3.1. Fresh-Frozen Tissue and FFPE Sample in ddPCR

Currently, the fixation of human tissues in a formalin solution followed by embedding them in paraffin is a common way to preserve a tumor sample for further studies. However, the process of sample fixation results in DNA fragmentation and the chemical modification of genomic DNA, such as the formation of crosslinks, cytosine deamination, and generation of abasic sites. These lesions can lead to DNA polymerase stalling anda reduction in DNA amplification efficacy dyeing PCR. To validate the suitability of the multiplex ddPCR assay for the FFPE samples analysis, we used matched fresh-frozen (FF) and FFPE samples from the single NSCLC tumor. DNA from matched FF and FFPE samples was purified and added in ddPCR with a concentration of 2000 copies per reaction, and all samples were tested in duplicate with the monoplex, duplex, and multiplex ddPCR assay (Figure 6).
The presented results indicate that, despite the degradation of DNA caused by the fixation process, DNA from FFPE samples remains suitable for the analysis. No difference in observed DNA concentrations was noticed for either monoplex, duplex, or multiplex analysis of fresh-frozen and FFPE sample. Thus, the developed multiplex ddPCR assay is suitable for the analysis of damaged DNA from FFPE samples.

3.3.2. Influence of DNA Concentration on the Performance of ddPCR

The amount of DNA in a reaction is one of the crucial variables that define the sensitivity of testing. A low DNA concentration in a sample, i.e., close to dozens of molecules per microliter, leads to the stochastic presence of DNA molecules in the ddPCR reaction and possible false results of the analysis. A high concentration of DNA, i.e., more than 104 molecules per microliter, results in an overload of ddPCR and a greater technical error with incorrect DNA quantification.
Another problem that could arise during the testing of FFPE samples is the ratio of amplified genes differing from whole numbers. Both the percentage of tumor cells and the exact degree of gene amplification can vary in a broad range in the analyzed specimen, leading to various observed gene amplification degrees. In that sense, the assessment of the minimal detectable gene amplification rate becomes necessary to determine the assay’s overall performance.
To assess the effect of the DNA concentration on the ddPCR performance, we prepared a set of samples with a varied ratio of MET and HER2 genes (2×, 4×, 8×) and different DNA concentrations. The samples were prepared by mixing genomic DNA from white blood cells of a healthy donor with linearized plasmids carrying fragments of MET and HER2 genes. The concentration of plasmids was measured with ddPCR and β-lactamase as a target for DNA quantification. Control samples were analyzed using the developed multiplex ddPCR-based assay. DNA was added in the reactions in the range of 100–15,000 copies per reaction. Simultaneously, the wild-type FFPE sample from the NSCLC tumor was titrated in ddPCR in the same range of DNA concentrations. All samples were tested in duplicate (Figure 7).
The copy number ratios were calculated by normalizing an amplicon copy number by the median value for all sample values and both reference values. To calculate the p-value, we evaluated the probability of meeting the normalized value for the same target among other samples with not less (for an amplification) or not more (for a deletion) than the value for a target considered using Equation (1):
p v a l u e = M x x i / x x i N ,
where M x x i / x x i is the number of samples for which the same target has the same or a higher (for amplification) or a lower (for deletion) copy number compared to both reference genes; x is the number of a sample; N is the total number of samples for which the current target has a high or low enough copy number.
For all samples, when 15,000, 3000, or 600 copies of DNA were loaded in reactions, the elevated ratio of MET or HER2 genes was clearly distinguishable. No higher presence of target genes was detected when no respective plasmid was added to the sample. The observed ratios of MET and HER2 genes were also close to the intended values, and calculating the ratios for the target genes was possible. The ratios of reference genes RAD51AP2 and RPP30 were also close to 1. When 100 copies of the template were analyzed, the SD of the ratio estimation increased and 2× and 4× HER2 ratios became indiscernible. However, even in that case, MET or HER2 amplifications were detectable. The titration of a wild-type FFPE NSCLC showed that the ratios of genes become unstable and varied at the range of 0.5–1.5, possibly hindering the CNA detection. Thus, the developed multiplex ddPCR assay can be used for the analysis in the range of DNA amount of 600–15,000 copies per reaction, and for the testing of FFPE samples with a low DNA concentration.

3.3.3. Testing of MET-Positive Sample

To validate the multiplex ddPCR on a real FFPE DNA sample, we used a DNA sample from a patient with NSCLC and found MET amplification. MET amplification was detected using a commercially available IVD MET (7q31)/SE7 (D7Z1) probe in accordance with the manufacturer’s instructions (Figure 8a). Consequently, the sample was analyzed using the developed ddPCR assay (Figure 8b).
MET amplification levels were defined on the basis of the MET-to-CEP7 ratio—low, greater than or equal to 1.8 and less than or equal to 2.2; medium, greater than 2.2 and less than 4.0; and high, greater than or equal to 4.0. The case showed heterogeneous distributions of MET signals, with 50% of cells without signs of amplification and 50% with high amplification levels greater than 5.0.
The multiplex ddPCR assay revealed MET amplification in the tested sample with a ratio of the MET gene to the reference RAD51AP2 and RPP30 genes equal to 3 (p ≤ 0.025). Therefore, the results of the ddPCR assay were confirmed by FISH, giving a MET amplification degree in the range of 3–5. The developed ddPCR was proven to be suitable for MET amplification detection, and we proceeded to test NSCLC DNA samples.

3.4. Testing of Clinical Samples

The validated multiplex ddPCR assay was used to analyze 436 EGFR-negative DNA samples from NSCLC FFPE tissue blocks. The samples were collected from patients from Siberia and the Far East of Russia and blinded before the analysis to hide the personal data of the patients. The EGFR status was determined using commercial qPCR-based kits. For each sample, the ratio of concentrations between target and reference genes was calculated. The obtained results are presented in Figure 9.
Among 436 tested NSCLC samples, five specimens (1.15%) showed a significantly higher ratio of MET to both reference loci (p ≤ 0.025), and six samples (1.38%) demonstrated a higher ratio of HER2 (p ≤ 0.05). The increased ratio of MET was in the range of 1.92–3.22, and the increased ratio of HER2 was in the range of 1.62–2.76. Keeping in mind the 15% common EGFR mutations rate in the studied population, the overall rate of MET and HER2 amplifications was 0.98% and 1.12%, respectively.
It should be noted that, in five samples, the MET ratio was in the range of 0.25–0.41, which can be a result of bad clustering or a low DNA concentration in the respective samples. Another possible explanation is the presence of somatic mutations in primers and probe sites, leading to the observed efficacy of amplification.
Interestingly, a higher ratio of RAD51AP2 (1.38–1.54) and a higher ratio of RPP30 (1.32–1.43) were observed in six and five samples, respectively. Such an increased ratio could be an analysis artifact caused by the selective inhibition of amplification due to the presence of contaminants in the tested samples or insufficient clustering. Another possible explanation is the presence of true amplification events involving reference genes in a small fraction of tumor cells. All samples with abnormal ratios of either target or reference genes were repeatedly analyzed using duplex ddPCR with one target and one reference gene. In all cases, the results of multiplex ddPCR were confirmed and similar ratios were obtained with duplex ddPCR (date not shown).

4. Discussion

The rapid growth of available and promising target drugs challenges modern molecular diagnostics in oncology. Emerging therapeutics demand the development of tests for the detection of respective target molecular alterations. Currently, some anti-MET and anti-HER2 drugs are undergoing clinical trials and require the design of appropriate molecular tests for MET and HER2 amplification. MET [36,37,38,39] and HER2 [40] alterations are also a reason for the resistance to commonly used anti-EGFR therapy. Existing approaches for MET and HER2 testing, such as IHC, FISH, or NGS, could be subjective [41,42], time- and labor-consuming, costly, or, in the case of NGS, are not fully implemented for practice. An additional issue is the low amount of starting material: as in the case of NSLSC, FFPE samples are often prepared from tumor biopsies. Thus, alternative methods are needed to address the mentioned problems. Digital PCR is a sensitive method requiring a low amount of DNA that gives precise information about the DNA concentration. The multiplexing of digital PCR further reduces labor and reagents costs. Therefore, multiplex digital PCR could be used to facilitate the molecular diagnostics of NSCLC. Numerous reports described monoplex or duplex ddPCR for the detection of MET and HER2 amplification in plasma and FFPE samples. However, to our knowledge, no multiplex ddPCR was designed for the simultaneous testing of MET and HER2 copy number alterations (CNAs).
In the present work, we applied the amplitude-based approach for the multiplexing of digital PCR. This approach relies on using probes with the same fluorophore for the detection of different targets. The probes are added to the reaction in different concentrations, resulting in a difference in final fluorescence and clustering. This strategy is widely used to design multiplex ddPCR, and was successfully applied by us to detect BRCA1 gene copy number variations.
To optimize reaction conditions for multiplex ddPCR, we tested the optimal elongation temperature and fluorescence levels of primers and probes for the target MET and HER2 gene and a number of possible reference genes. Simultaneously with digital PCR, primers and probes were tested in quantitative PCR; the latter method is less time-consuming and costly. Thus, we assessed the possibility of extrapolating the results of qPCR on ddPCR for speeding up and lowering the cost of ddPCR optimization.
A direct comparison of qPCR and ddPCR results gave contradicting results. Thus, the decrease in qPCR efficacy at a lower elongation temperature did not correlate with the efficacy of ddPCR. Another important variable, the basic fluorescence level in ddPCR, cannot be predicted using qPCR. The ratio of basic fluorescence between different primers was different in qPCR and ddPCR, resulting in an impossible direct prediction from qPCR data. It should be recalled that primers and probes with a high, basic, and low end-point fluorescence can be obtained using qPCR and do not have to be used for ddPCR, saving time and effort. End-point fluorescence was more stable between qPCR and ddPCR; the ratio of fluorescence between different targets was similar in qPCR and ddPCR. Thus, end-point fluorescence in ddPCR can be forecasted based on qPCR. To sum up, qPCR has a limited predicting ability for the determination of ddPCR optimal conditions. It allows us to predict end-point fluorescence for various targets, but has reservations in the farseeing of the optimal elongation temperature and basic fluorescence. Therefore, the testing of primers and probes in ddPCR cannot be skipped and replaced by qPCR.
Several targets were chosen as two possible optimal references. Such precautions were made to prevent potential issues with a lower PCR efficacy after multiplexing, as primers for different targets can interfere and affect clustering in ddPCR. After testing in monoplexes, primers for different loci were tested in duplex ddPCR. Theoretically, data on the optimal elongation temperature and fluorescence from monoplexes are enough to calculate optimal primer and probe concentrations for multiplexes, and a duplex test is not necessary. However, primers from different sets can be incompatible in a single reaction, as was mentioned above, and this possibility can only be assessed in an actual experiment. Here, results for almost all combinations of reference and target genes were in accordance with the data from monoplexes. However, HER2/PARK7 combinations demonstrated rain in a duplex that was not observed in the respective monoplexes. This notion underlines the need for the prior duplex testing of different primer combinations before combining them in a more complex multiplex.
At the final step of optimization, the resulting tetraplex assay was validated on a small set of real FFPE samples to determine the actual performance ddPCR. It was observed that the target loci ratio was lower than that for the reference gene. The possible reason behind the underestimation was the length of target gene amplicons, which was 30 bp longer than the reference loci. It is well-known that DNA from FFPE samples is highly fragmented and chemically modified. Here, it led to the underestimation of target genes due to a lower concentration of longer DNA template molecules. After introducing shorter amplicons, the ratio for MET and HER2 genes became close to theoretical 1, supporting the hypothesis about the influence of DNA amplicons. This observation highlights the need for a careful primer design for the amplification of FFPE samples when the length of amplicons can significantly affect the actual performance. Commercial assays for different targets should be assessed with special precautions, as the length of amplicons could be unknown. Not matching amplicon lengths can result in biased data, especially for CNA detection.
The developed multiplex ddPCR assay was successfully validated using a MET-positive NSCLC sample and a set of 436 EGFR-negative samples. Thus, ddPCR confirmed the presence of MET amplification in the sample tested by FISH. After correction on the frequency of common EGFR mutations, the observed rates of MET and HER2 amplifications were 0.98% and 1.12%, respectively. These findings are in good agreement with previous studies, demonstrating a similar 1–5% frequency of MET and HER2 CNAs. Several works also reported the coexistence of EGFR mutations with various MET and HER2 alterations; however, those events were rare [19]. For that reason, no EGFR-positive samples were tested here.
Several limitations of the present study should be mentioned. The set of FFPE samples used for the validation of the developed ddPCR assay was relatively small. The insufficient number of tested NSCLC specimens could lead to an erroneous estimation of actual MET and HER2 amplification rates. Only a single sample was validated as MET-positive using gold standard FISH, while a poor correlation between FISH and ddPCR was once reported [43].
Droplet digital PCR also has some reservations. If handled without proper care, ddPCR could be contaminated by amplicons produced in previous runs. Contamination leads to biased results and requires substantial efforts for elimination. Therefore, the ddPCR workflow should be organized in the one-way movement of samples and should prevent any possible contacts of ddPCR reagents with samples and consumables after amplification. Digital PCR is a targeted method that could detect only known mutations, whereas NGS allows for the detection of all alterations in analyzed genomic sites. Expanding the spectrum of novel targets makes NGS promising for the simultaneous detection of all possible mutations with available target drugs. Targeted methods require numerous separate reactions to test many mutations, whereas NGS allows us to analyze all of them in a single run. However, ddPCR remains more sensitive than NGS and can clarify the results of a preliminary screening made by NGS or other methods. A high sensitivity of the analysis can be crucial in the case of liquid biopsy testing, where extremely low amounts of tumor DNA could be present in the studied samples. Therefore, ddPCR could be applied for circulating tumor DNA analysis, taking a sensitive approach for analyzing samples with a low amount of tumor DNA. Monoplex ddPCR assays were successfully validated previously, making it possible to develop more complex multiplex ddPCR techniques.
To sum up, the developed ddPCR-based assay was optimized using both qPCR and ddPCR data and validated on samples with various DNA concentrations and with the MET and HER2 amplification degree. The assay was used to test 436 EGFR-negative NSCLC samples. The observed rates of MET and HER2 amplification were in agreement with previous studies.

5. Conclusions

Here, we report the first multiplex droplet digital PCR assay for the simultaneous detection of MET and HER2 amplification in NSCLC samples. The suitability of qPCR was assessed for the optimization of multiplex ddPCR. The optimal elongation temperature, reference genes for DNA quantification, and amplicon length were selected. The developed multiplex ddPCR was validated on control samples with various DNA concentrations and ratios of MET and HER2 genes, and on a MET-positive NSCLC sample. Using a ddPCR-based assay, 436 EGFR-negative NSCLC samples were analyzed. Among the tested samples, five specimens (1.15%) showed a significantly higher ratio of MET (p ≤ 0.025), and six samples (1.38%) demonstrated a higher ratio of HER2 (p ≤ 0.05). The increased ratio of MET was in the range of 1.92–3.22 and the increased ratio of HER2 was in the range of 1.62–2.76. The multiplex ddPCR assay reported here could be used for a routine screening of MET and HER2 amplification in NSCLC samples.

Author Contributions

Conceptualization, M.L.F. and I.A.D.; methodology, I.P.O.; software, A.A.K.; validation, M.A.S., K.A.P., M.E.V. and U.A.B.; formal analysis, A.A.K.; investigation, M.A.S., K.A.P., M.E.V. and U.A.B.; resources, M.L.F.; data curation, M.A.S.; writing—original draft preparation, I.P.O.; writing—review and editing, M.L.F.; visualization, I.P.O. and A.A.K.; supervision, M.L.F.; project administration, M.L.F. and I.A.D.; funding acquisition, M.L.F. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Russian State funded budget project 121031300045-2 for ICBFM SB RAS. The APC was funded by 121031300045-2.

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki, and approved by the Local Ethics Committee of ICBFM SB RAS (protocol code 7, date of approval 10 July 2021).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study. Written informed consent has been obtained from the patients to publish this paper.

Data Availability Statement

The data presented in this study are available on request from the corresponding author. The data are not publicly available due to ethical reasons.

Acknowledgments

The authors greatly appreciate the efforts of Vitali E. Kuznetsov in the technical support of ddPCR experiments. The authors greatly appreciate assistance of the inter-regional non-governmental organization “Society of molecular geneticists in oncology and oncohematology.”

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Sung, H.; Ferlay, J.; Siegel, R.L.; Laversanne, M.; Soerjomataram, I.; Jemal, A.; Bray, F. Global Cancer Statistics 2020: GLOBOCAN Estimates of Incidence and Mortality Worldwide for 36 Cancers in 185 Countries. CA Cancer J. Clin. 2021, 71, 209–249. [Google Scholar] [CrossRef]
  2. Thai, A.A.; Solomon, B.J.; Sequist, L.V.; Gainor, J.F.; Heist, R.S. Lung cancer. Lancet 2021, 398, 535–554. [Google Scholar] [CrossRef]
  3. Gridelli, C.; Rossi, A.; Carbone, D.P.; Guarize, J.; Karachaliou, N.; Mok, T.; Petrella, F.; Spaggiari, L.; Rosell, R. Non-small-cell lung cancer. Nat. Rev. Dis. Primers 2015, 1, 15009. [Google Scholar] [CrossRef]
  4. Lynch, T.J.; Bell, D.W.; Sordella, R.; Gurubhagavatula, S.; Okimoto, R.A.; Brannigan, B.W.; Harris, P.L.; Haserlat, S.M.; Supko, J.G.; Haluska, F.G.; et al. Activating mutations in the epidermal growth factor receptor underlying responsiveness of non-small-cell lung cancer to gefitinib. N. Engl. J. Med. 2004, 350, 2129–2139. [Google Scholar] [CrossRef]
  5. Kwak, E.L.; Bang, Y.-J.; Camidge, D.R.; Shaw, A.T.; Solomon, B.; Maki, R.G.; Ou, S.-H.I.; Dezube, B.J.; Jänne, P.A.; Costa, D.B.; et al. Anaplastic lymphoma kinase inhibition in non-small-cell lung cancer. N. Engl. J. Med. 2010, 363, 1693–1703. [Google Scholar] [CrossRef] [Green Version]
  6. Melosky, B.; Wheatley-Price, P.; Juergens, R.A.; Sacher, A.; Leighl, N.B.; Tsao, M.-S.; Cheema, P.; Snow, S.; Liu, G.; Card, P.B.; et al. The rapidly evolving landscape of novel targeted therapies in advanced non-small cell lung cancer. Lung Cancer 2021, 160, 136–151. [Google Scholar] [CrossRef]
  7. Hong, D.S.; Fakih, M.G.; Strickler, J.H.; Desai, J.; Durm, G.A.; Shapiro, G.I.; Falchook, G.S.; Price, T.J.; Sacher, A.; Denlinger, C.S.; et al. KRAS G12C Inhibition with Sotorasib in Advanced Solid Tumors. N. Engl. J. Med. 2020, 383, 1207–1217. [Google Scholar] [CrossRef]
  8. Park, K.; Haura, E.B.; Leighl, N.B.; Mitchell, P.; Shu, C.A.; Girard, N.; Viteri, S.; Han, J.-Y.; Kim, S.-W.; Lee, C.K.; et al. Amivantamab in EGFR Exon 20 Insertion-Mutated Non-Small-Cell Lung Cancer Progressing on Platinum Chemotherapy: Initial Results From the CHRYSALIS Phase I Study. J. Clin. Oncol. 2021, 39, 3391–3402. [Google Scholar] [CrossRef]
  9. Trusolino, L.; Bertotti, A.; Comoglio, P.M. MET signalling: Principles and functions in development, organ regeneration and cancer. Nat. Rev. Mol. Cell Biol. 2010, 11, 834–848. [Google Scholar] [CrossRef]
  10. Fujino, T.; Suda, K.; Mitsudomi, T. Emerging MET tyrosine kinase inhibitors for the treatment of non-small cell lung cancer. Expert Opin. Emerg. Drugs 2020, 25, 229–249. [Google Scholar] [CrossRef]
  11. Friedlaender, A.; Drilon, A.; Banna, G.L.; Peters, S.; Addeo, A. The METeoric rise of MET in lung cancer. Cancer 2020, 126, 4826–4837. [Google Scholar] [CrossRef]
  12. Malik, S.M.; Maher, V.E.; Bijwaard, K.E.; Becker, R.L.; Zhang, L.; Tang, S.W.; Song, P.; Liu, Q.; Marathe, A.; Gehrke, B.; et al. U.S. Food and Drug Administration approval: Crizotinib for treatment of advanced or metastatic non-small cell lung cancer that is anaplastic lymphoma kinase positive. Clin. Cancer Res. 2014, 20, 2029–2034. [Google Scholar] [CrossRef] [Green Version]
  13. Iqbal, N.; Iqbal, N. Human Epidermal Growth Factor Receptor 2 (HER2) in Cancers: Overexpression and Therapeutic Implications. Mol. Biol. Int. 2014, 2014, 852748. [Google Scholar] [CrossRef]
  14. Albanell, J.; Baselga, J. Trastuzumab, a humanized anti-HER2 monoclonal antibody, for the treatment of breast cancer. Drugs Today (Barc) 1999, 35, 931–946. [Google Scholar]
  15. Blumenthal, G.M.; Scher, N.S.; Cortazar, P.; Chattopadhyay, S.; Tang, S.; Song, P.; Liu, Q.; Ringgold, K.; Pilaro, A.M.; Tilley, A.; et al. First FDA approval of dual anti-HER2 regimen: Pertuzumab in combination with trastuzumab and docetaxel for HER2-positive metastatic breast cancer. Clin. Cancer Res. 2013, 19, 4911–4916. [Google Scholar] [CrossRef] [Green Version]
  16. Royce, M.; Osgood, C.L.; Amatya, A.K.; Fiero, M.H.; Chang, C.J.G.; Ricks, T.K.; Shetty, K.A.; Kraft, J.; Qiu, J.; Song, P.; et al. FDA Approval Summary: Margetuximab plus Chemotherapy for Advanced or Metastatic HER2-Positive Breast Cancer. Clin. Cancer Res. 2021. [Google Scholar] [CrossRef]
  17. Jebbink, M.; de Langen, A.J.; Boelens, M.C.; Monkhorst, K.; Smit, E.F. The force of HER2—A druggable target in NSCLC? Cancer Treat. Rev. 2020, 86, 101996. [Google Scholar] [CrossRef]
  18. Zhao, J.; Xia, Y. Targeting HER2 Alterations in Non–Small-Cell Lung Cancer: A Comprehensive Review. JCO Precis. Oncol. 2020, 4, 411–425. [Google Scholar] [CrossRef]
  19. Nagasaka, M.; Singh, V.; Baca, Y.; Sukari, A.; Kim, C.; Mamdani, H.; Spira, A.I.; Uprety, D.; Bepler, G.; Kim, E.S.; et al. The Effects of HER2 Alterations in EGFR Mutant Non-small Cell Lung Cancer. Clin. Lung Cancer 2022, 23, 52–59. [Google Scholar] [CrossRef]
  20. Li, B.T.; Smit, E.F.; Goto, Y.; Nakagawa, K.; Udagawa, H.; Mazières, J.; Nagasaka, M.; Bazhenova, L.; Saltos, A.N.; Felip, E.; et al. Trastuzumab Deruxtecan in HER2-Mutant Non-Small-Cell Lung Cancer. N. Engl. J. Med. 2022, 386, 241–251. [Google Scholar] [CrossRef]
  21. Zhang, Y.; Tang, E.-T.; Du, Z. Detection of MET Gene Copy Number in Cancer Samples Using the Droplet Digital PCR Method. PLoS ONE 2016, 11, e0146784. [Google Scholar] [CrossRef]
  22. Ishii, H.; Azuma, K.; Sakai, K.; Kawahara, A.; Yamada, K.; Tokito, T.; Okamoto, I.; Nishio, K.; Hoshino, T. Digital PCR analysis of plasma cell-free DNA for non-invasive detection of drug resistance mechanisms in EGFR mutant NSCLC: Correlation with paired tumor samples. Oncotarget 2015, 6, 30850–30858. [Google Scholar] [CrossRef] [Green Version]
  23. Oddo, D.; Siravegna, G.; Gloghini, A.; Vernieri, C.; Mussolin, B.; Morano, F.; Crisafulli, G.; Berenato, R.; Corti, G.; Volpi, C.C.; et al. Emergence of MET hyper-amplification at progression to MET and BRAF inhibition in colorectal cancer. Br. J. Cancer 2017, 117, 347–352. [Google Scholar] [CrossRef]
  24. He, H.-J.; Das, B.; Cleveland, M.H.; Chen, L.; Camalier, C.E.; Liu, L.-C.; Norman, K.L.; Fellowes, A.P.; McEvoy, C.R.; Lund, S.P.; et al. Development and interlaboratory evaluation of a NIST Reference Material RM 8366 for EGFR and MET gene copy number measurements. Clin. Chem. Lab. Med. 2019, 57, 1142–1152. [Google Scholar] [CrossRef]
  25. Chen, W.; Zheng, J.; Wu, C.; Liu, S.; Chen, Y.; Liu, X.; Du, J.; Wang, J. Breast Cancer Subtype Classification Using 4-Plex Droplet Digital PCR. Clin. Chem. 2019, 65, 1051–1059. [Google Scholar] [CrossRef]
  26. Gevensleben, H.; Garcia-Murillas, I.; Graeser, M.K.; Schiavon, G.; Osin, P.; Parton, M.; Smith, I.E.; Ashworth, A.; Turner, N.C. Noninvasive detection of HER2 amplification with plasma DNA digital PCR. Clin. Cancer Res. 2013, 19, 3276–3284. [Google Scholar] [CrossRef] [Green Version]
  27. Belgrader, P.; Tanner, S.C.; Regan, J.F.; Koehler, R.; Hindson, B.J.; Brown, A.S. Droplet Digital PCR Measurement of HER2 Copy Number Alteration in Formalin-Fixed Paraffin-Embedded Breast Carcinoma Tissue. Clin. Chem. 2013, 59, 991–994. [Google Scholar] [CrossRef] [Green Version]
  28. Heredia, N.J.; Belgrader, P.; Wang, S.; Koehler, R.; Regan, J.; Cosman, A.M.; Saxonov, S.; Hindson, B.; Tanner, S.C.; Brown, A.S.; et al. Droplet DigitalTM PCR quantitation of HER2 expression in FFPE breast cancer samples. Methods 2013, 59, S20–S23. [Google Scholar] [CrossRef]
  29. Turner, N.C.; Kingston, B.; Kilburn, L.S.; Kernaghan, S.; Wardley, A.M.; Macpherson, I.R.; Baird, R.D.; Roylance, R.; Stephens, P.; Oikonomidou, O.; et al. Circulating tumour DNA analysis to direct therapy in advanced breast cancer (plasmaMATCH): A multicentre, multicohort, phase 2a, platform trial. Lancet Oncol. 2020, 21, 1296–1308. [Google Scholar] [CrossRef]
  30. He, H.-J.; Almeida, J.L.; Lund, S.P.; Steffen, C.R.; Choquette, S.; Cole, K.D. Development of NIST standard reference material 2373: Genomic DNA standards for HER2 measurements. Biomol. Detect. Quantif. 2016, 8, 1–8. [Google Scholar] [CrossRef] [Green Version]
  31. Camidge, D.R.; Otterson, G.A.; Clark, J.W.; Ignatius Ou, S.-H.; Weiss, J.; Ades, S.; Shapiro, G.I.; Socinski, M.A.; Murphy, D.A.; Conte, U.; et al. Crizotinib in Patients With MET-Amplified NSCLC. J. Thorac. Oncol. 2021, 16, 1017–1029. [Google Scholar] [CrossRef]
  32. Jones, M.; Williams, J.; Gärtner, K.; Phillips, R.; Hurst, J.; Frater, J. Low copy target detection by Droplet Digital PCR through application of a novel open access bioinformatic pipeline, “definetherain”. J. Virol. Methods 2014, 202, 46–53. [Google Scholar] [CrossRef] [Green Version]
  33. Della Starza, I.; Nunes, V.; Cavalli, M.; De Novi, L.A.; Ilari, C.; Apicella, V.; Vitale, A.; Testi, A.M.; Del Giudice, I.; Chiaretti, S.; et al. Comparative analysis between RQ-PCR and digital-droplet-PCR of immunoglobulin/T-cell receptor gene rearrangements to monitor minimal residual disease in acute lymphoblastic leukaemia. Br. J. Haematol. 2016, 174, 541–549. [Google Scholar] [CrossRef] [Green Version]
  34. Oscorbin, I.; Kechin, A.; Boyarskikh, U.; Filipenko, M. Multiplex ddPCR assay for screening copy number variations in BRCA1 gene. Breast Cancer Res. Treat. 2019, 178, 545–555. [Google Scholar] [CrossRef]
  35. Beroukhim, R.; Mermel, C.H.; Porter, D.; Wei, G.; Raychaudhuri, S.; Donovan, J.; Barretina, J.; Boehm, J.S.; Dobson, J.; Urashima, M.; et al. The landscape of somatic copy-number alteration across human cancers. Nature 2010, 463, 899–905. [Google Scholar] [CrossRef]
  36. Ramalingam, S.S.; Cheng, Y.; Zhou, C.; Ohe, Y.; Imamura, F.; Cho, B.C.; Lin, M.-C.; Majem, M.; Shah, R.; Rukazenkov, Y.; et al. Mechanisms of acquired resistance to first-line osimertinib: Preliminary data from the phase III FLAURA study. Ann. Oncol. 2018, 29, viii740. [Google Scholar] [CrossRef]
  37. Papadimitrakopoulou, V.A.; Wu, Y.-L.; Han, J.-Y.; Ahn, M.-J.; Ramalingam, S.S.; John, T.; Okamoto, I.; Yang, J.C.-H.; Bulusu, K.C.; Laus, G.; et al. Analysis of resistance mechanisms to osimertinib in patients with EGFR T790M advanced NSCLC from the AURA3 study. Ann. Oncol. 2018, 29, viii741. [Google Scholar] [CrossRef]
  38. Turke, A.B.; Zejnullahu, K.; Wu, Y.-L.; Song, Y.; Dias-Santagata, D.; Lifshits, E.; Toschi, L.; Rogers, A.; Mok, T.; Sequist, L.; et al. Preexistence and clonal selection of MET amplification in EGFR mutant NSCLC. Cancer Cell 2010, 17, 77–88. [Google Scholar] [CrossRef] [Green Version]
  39. Bean, J.; Brennan, C.; Shih, J.-Y.; Riely, G.; Viale, A.; Wang, L.; Chitale, D.; Motoi, N.; Szoke, J.; Broderick, S.; et al. MET amplification occurs with or without T790M mutations in EGFR mutant lung tumors with acquired resistance to gefitinib or erlotinib. Proc. Natl. Acad. Sci. USA. 2007, 104, 20932–20937. [Google Scholar] [CrossRef] [Green Version]
  40. Yu, H.A.; Arcila, M.E.; Rekhtman, N.; Sima, C.S.; Zakowski, M.F.; Pao, W.; Kris, M.G.; Miller, V.A.; Ladanyi, M.; Riely, G.J. Analysis of Tumor Specimens at the Time of Acquired Resistance to EGFR-TKI Therapy in 155 Patients with EGFR -Mutant Lung Cancers. Clin. Cancer Res. 2013, 19, 2240–2247. [Google Scholar] [CrossRef] [Green Version]
  41. Guo, R.; Berry, L.D.; Aisner, D.L.; Sheren, J.; Boyle, T.; Bunn, P.A.; Johnson, B.E.; Kwiatkowski, D.J.; Drilon, A.; Sholl, L.M.; et al. MET IHC Is a Poor Screen for MET Amplification or MET Exon 14 Mutations in Lung Adenocarcinomas: Data from a Tri-Institutional Cohort of the Lung Cancer Mutation Consortium. J. Thorac. Oncol. 2019, 14, 1666–1671. [Google Scholar] [CrossRef]
  42. Nitta, H.; Kelly, B.D.; Allred, C.; Jewell, S.; Banks, P.; Dennis, E.; Grogan, T.M. The assessment of HER2 status in breast cancer: The past, the present, and the future. Pathol. Int. 2016, 66, 313–324. [Google Scholar] [CrossRef] [Green Version]
  43. Heydt, C.; Becher, A.-K.; Wagener-Ryczek, S.; Ball, M.; Schultheis, A.M.; Schallenberg, S.; Rüsseler, V.; Büttner, R.; Merkelbach-Bruse, S. Comparison of in situ and extraction-based methods for the detection of MET amplifications in solid tumors. Comput. Struct. Biotechnol. J. 2019, 17, 1339–1347. [Google Scholar] [CrossRef]
Figure 1. Results of qPCR with monoplexes. (a) Cq values; (b) end-point fluorescence after baseline subtraction. X-axis marks elongation temperature, Y-axis marks Cq values (a), or end-point fluorescence (b). Whiskers designate standard deviation in two experiments.
Figure 1. Results of qPCR with monoplexes. (a) Cq values; (b) end-point fluorescence after baseline subtraction. X-axis marks elongation temperature, Y-axis marks Cq values (a), or end-point fluorescence (b). Whiskers designate standard deviation in two experiments.
Cancers 14 01458 g001
Figure 2. Results of ddPCR with monoplexes. Each panel represents results for the set of primers and probes that are designated at the left upper angle. X-axis marks elongation temperature, Y-axis marks end-point fluorescence.
Figure 2. Results of ddPCR with monoplexes. Each panel represents results for the set of primers and probes that are designated at the left upper angle. X-axis marks elongation temperature, Y-axis marks end-point fluorescence.
Cancers 14 01458 g002
Figure 3. Fluorescence of monoplexes in qPCR and ddPCR. (a) Basic fluorescence of qPCR; (b) basic fluorescence of ddPCR; (c) end-point fluorescence after baseline subtraction of qPCR; (d) amplitude of end-point fluorescence ddPCR. Y-axis marks elongation temperature and monoplex, X-axis marks fluorescence in relative fluorescence units. Whiskers designate standard deviation in two experiments.
Figure 3. Fluorescence of monoplexes in qPCR and ddPCR. (a) Basic fluorescence of qPCR; (b) basic fluorescence of ddPCR; (c) end-point fluorescence after baseline subtraction of qPCR; (d) amplitude of end-point fluorescence ddPCR. Y-axis marks elongation temperature and monoplex, X-axis marks fluorescence in relative fluorescence units. Whiskers designate standard deviation in two experiments.
Cancers 14 01458 g003
Figure 4. Results of ddPCR with duplexes. Each panel represents results for the set of primers and probes: (a) HER2-1, ALB, (b) HER2-1, PARK7, (c) HER2-1, RAD51AP2, (d) MET-1, RPP30, (e) MET-1, TEKT3. X-axis represents combinations of primers and probes with different concentrations, Y-axis represents end-point fluorescence. High concentrations of primers/probe correspond to 900/250 nM, low concentrations of primers/probe correspond to 450/125 nM. Clusters are designated at the sides of each panel.
Figure 4. Results of ddPCR with duplexes. Each panel represents results for the set of primers and probes: (a) HER2-1, ALB, (b) HER2-1, PARK7, (c) HER2-1, RAD51AP2, (d) MET-1, RPP30, (e) MET-1, TEKT3. X-axis represents combinations of primers and probes with different concentrations, Y-axis represents end-point fluorescence. High concentrations of primers/probe correspond to 900/250 nM, low concentrations of primers/probe correspond to 450/125 nM. Clusters are designated at the sides of each panel.
Cancers 14 01458 g004
Figure 5. Results of pilot ddPCR with DNA from NSCLC FFPE samples. (a) 2D ddPCR plot. X-axis represents the amplitude of HEX-channel, Y-axis represents amplitude of FAM-channel. Loci are marked near respective clusters. (b) Tukey box-and-whiskers plot, X-axis represents genes, Y-axis represents the ratio of gene concentration to median concentration for both reference loci. Data for long amplicons are presented at the right side of the plot, data for short amplicons are presented at the left side.
Figure 5. Results of pilot ddPCR with DNA from NSCLC FFPE samples. (a) 2D ddPCR plot. X-axis represents the amplitude of HEX-channel, Y-axis represents amplitude of FAM-channel. Loci are marked near respective clusters. (b) Tukey box-and-whiskers plot, X-axis represents genes, Y-axis represents the ratio of gene concentration to median concentration for both reference loci. Data for long amplicons are presented at the right side of the plot, data for short amplicons are presented at the left side.
Cancers 14 01458 g005
Figure 6. Influence of sample type on ddPCR performance. X-axis represents the concentration of genes; Y-axis shows the DNA sample and type of ddPCR—monoplex, duplex, or multiplex. FF—fresh-frozen, FFPE—formalin-fixed, paraffin-embedded samples. Whiskers designate standard deviation in two experiments.
Figure 6. Influence of sample type on ddPCR performance. X-axis represents the concentration of genes; Y-axis shows the DNA sample and type of ddPCR—monoplex, duplex, or multiplex. FF—fresh-frozen, FFPE—formalin-fixed, paraffin-embedded samples. Whiskers designate standard deviation in two experiments.
Cancers 14 01458 g006
Figure 7. Influence of DNA amount on ddPCR performance. (a) 15,000 copies per reaction; (b) 3000 copies per reaction; (c) 600 copies per reaction; (d) 100 copies per reaction; (e) FFPE NSCLC sample titration. Each specimen was analyzed in duplicate. X-axis represents the ratio of target or reference gene concentration to normalized concentration of both reference genes, and Y-axis shows the DNA sample. Whiskers designate standard deviation in two experiments.
Figure 7. Influence of DNA amount on ddPCR performance. (a) 15,000 copies per reaction; (b) 3000 copies per reaction; (c) 600 copies per reaction; (d) 100 copies per reaction; (e) FFPE NSCLC sample titration. Each specimen was analyzed in duplicate. X-axis represents the ratio of target or reference gene concentration to normalized concentration of both reference genes, and Y-axis shows the DNA sample. Whiskers designate standard deviation in two experiments.
Cancers 14 01458 g007aCancers 14 01458 g007b
Figure 8. Results of MET-positive NSCLC sample analysis using FISH and ddPCR. (a) FISH assay. Green dots correspond to D7Z1 probes marking the chromosome 7 centromere, red dots correspond to probes marking the MET gene region at 7q31. (b) Multiplex ddPCR assay. X-axis shows genes, and Y-axis represents the ratio of gene concentration to normalized concentration of both reference genes.
Figure 8. Results of MET-positive NSCLC sample analysis using FISH and ddPCR. (a) FISH assay. Green dots correspond to D7Z1 probes marking the chromosome 7 centromere, red dots correspond to probes marking the MET gene region at 7q31. (b) Multiplex ddPCR assay. X-axis shows genes, and Y-axis represents the ratio of gene concentration to normalized concentration of both reference genes.
Cancers 14 01458 g008
Figure 9. Results of NSCLC FFPE sample analyses by ddPCR. Tukey box-and-whiskers plot, X-axis represents genes, Y-axis represents the ratio of gene concentration to median concentration for both reference loci.
Figure 9. Results of NSCLC FFPE sample analyses by ddPCR. Tukey box-and-whiskers plot, X-axis represents genes, Y-axis represents the ratio of gene concentration to median concentration for both reference loci.
Cancers 14 01458 g009
Table 1. Sequences of primers and probes.
Table 1. Sequences of primers and probes.
Gene, Amplicon LengthOligonucleotide5′-Sequence-3′
HER2 set 1 100 bpHER2-F1CATCTTCCACAAGAACAACCAG
HER2-R1GTCTGAGAGAAGAGGCAGCA
HER2-PFAM-AGACACCAACCGCTCTCGGGC-BHQ1
HER2 set 2 72 bpHER2-F2CAGCTGGCTCTCACACTGA
HER2-R2AAGAGGCAGCAGGGAGG
HER2-PFAM-AGACACCAACCGCTCTCGGGC-BHQ1
MET set 1 100 bpMET-F1AAGTGGGGAACTGATGTGACTTA
MET-R1CAAGATCGTCAACAAAAACAATG
MET-PHEX-CGGACCCAATCATGAGCACTG-BHQ1
MET set 2 77 bpMET-F2CAATGTGAGATGTCTCCAGCAT
MET-R2GGGAACTGATGTGACTTACCCTA
MET-PHEX-CGGACCCAATCATGAGCACTG-BHQ1
ALB 94 bpALB-UGACTTGCCAAGACATATGAAACC
ALB-FAMFAM–TGCTGTGCCGCTGCAGATCC-BHQ1
ALB-RTCCAACAATAAACCTACCACTTTG
PARK7 74 bpPARK7-FCGTTGCAGGCCTGGCTGGAAA
PARK7-RCAAGGCTGGCATCAGGACAAATG
PARK7-PFAM-CCCAGTACAGTGTAGCCGTGATGTGG-BHQ1
RAD51AP2 74 bpRAD51AP2-UTGGTGACTTTTGCCCATATTAG
RAD51AP2-RGTGGTCAACAAAATACGTGCA
RAD51AP2-PFAM-CCTGCTATAGTATCATGGAACGAGG-BHQ1
RPP30 61 bpRPP30-FCAGATTTGGACCTGCGAGCG
RPP30-RGAGCGGCTGTCTCCACAAGT
RP30-HEXHEX-TCTGACCTGAAGGCTCTGCGCG-BHQ2
TEKT3 67 bpTEKT3-FGCGTGTGCCTGACTTACGTTG
TEKT3-RACATCAGACGGTGTCGGCTAC
TEKT3-PHEX-CGACCCTCTCCACTCCGCG-BHQ2
Beta-lactamaseLac-UCGTCGTTTGGTATGGCTTCATTC
Lac-RAGGACCGAAGGAGCTAACCG
Lac-PHEX-CGGTTCCCAACGATCAAGGCGAG-BHQ2
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Oscorbin, I.P.; Smertina, M.A.; Pronyaeva, K.A.; Voskoboev, M.E.; Boyarskikh, U.A.; Kechin, A.A.; Demidova, I.A.; Filipenko, M.L. Multiplex Droplet Digital PCR Assay for Detection of MET and HER2 Genes Amplification in Non-Small Cell Lung Cancer. Cancers 2022, 14, 1458. https://doi.org/10.3390/cancers14061458

AMA Style

Oscorbin IP, Smertina MA, Pronyaeva KA, Voskoboev ME, Boyarskikh UA, Kechin AA, Demidova IA, Filipenko ML. Multiplex Droplet Digital PCR Assay for Detection of MET and HER2 Genes Amplification in Non-Small Cell Lung Cancer. Cancers. 2022; 14(6):1458. https://doi.org/10.3390/cancers14061458

Chicago/Turabian Style

Oscorbin, Igor P., Maria A. Smertina, Ksenia A. Pronyaeva, Mikhail E. Voskoboev, Ulyana A. Boyarskikh, Andrey A. Kechin, Irina A. Demidova, and Maxim L. Filipenko. 2022. "Multiplex Droplet Digital PCR Assay for Detection of MET and HER2 Genes Amplification in Non-Small Cell Lung Cancer" Cancers 14, no. 6: 1458. https://doi.org/10.3390/cancers14061458

APA Style

Oscorbin, I. P., Smertina, M. A., Pronyaeva, K. A., Voskoboev, M. E., Boyarskikh, U. A., Kechin, A. A., Demidova, I. A., & Filipenko, M. L. (2022). Multiplex Droplet Digital PCR Assay for Detection of MET and HER2 Genes Amplification in Non-Small Cell Lung Cancer. Cancers, 14(6), 1458. https://doi.org/10.3390/cancers14061458

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop