Clonal Origin and Lineage Ambiguity in Mixed Neuroendocrine Carcinoma of the Uterine Cervix.

Small-cell neuroendocrine carcinoma (SCNEC) of the cervix is a rare disease characterized by a high incidence of mixed tumors with other types of cancer. The mechanism underlying this mixed phenotype is not well understood. This study established a panel of organoid lines from patients with SCNEC of the cervix and ultimately focused on one line, which retained a mixed tumor phenotype, both in vitro and in vivo. Histologically, both organoids and xenograft tumors showed distinct differentiation into either SCNEC or adenocarcinoma in some regions and ambiguous differentiation in others. Tracking single cells indicated the existence of cells with bipotential differentiation toward SCNEC and adenocarcinomas. Single-cell transcriptional analysis identified three distinct clusters: SCNEC-like, adenocarcinoma-like, and a cluster lacking specific differentiation markers. The expression of neuroendocrine markers was enriched in the SCNEC-like cluster but not exclusively. Human papillomavirus 18 E6 was enriched in the SCNEC-like cluster, which showed higher proliferation and lower levels of the p53 pathway. After treatment with anticancer drugs, the expression of adenocarcinoma markers increased, whereas that of SCNEC decreased. Using a reporter system for keratin 19 expression, changes in the differentiation of each cell were shown to be associated with the shift in differentiation induced by drug treatment. These data suggest that mixed SCNEC/cervical tumors have a clonal origin and are characterized by an ambiguous and flexible differentiation state.

Analytical Performance of a Highly Sensitive System to Detect Gene Variants Using Next-Generation Sequencing for Lung Cancer Companion Diagnostics.

The recent increase in the number of molecular targeted agents for lung cancer has led to the demand for the simultaneous testing of multiple genes. Although gene panels using next-generation sequencing (NGS) are ideal, conventional panels require a high tumor content, and biopsy samples often do not meet this requirement. We developed a new NGS panel, called compact panel, characterized by high sensitivity, with detection limits for mutations of 0.14%, 0.20%, 0.48%, 0.24%, and 0.20% for EGFR exon 19 deletion, L858R, T790M, BRAF V600E, and KRAS G12C, respectively. Mutation detection also had a high quantitative ability, with correlation coefficients ranging from 0.966 to 0.992. The threshold for fusion detection was 1%. The panel exhibited good concordance with the approved tests. The identity rates were as follows: EGFR positive, 100% (95% confidence interval, 95.5-100); EGFR negative, 90.9 (82.2-96.3); BRAF positive, 100 (59.0-100); BRAF negative, 100 (94.9-100); KRAS G12C positive, 100 (92.7-100); KRAS G12C negative, 100 (93.0-100); ALK positive, 96.7 (83.8-99.9); ALK negative, 98.4 (97.2-99.2); ROS1 positive, 100 (66.4-100); ROS1 negative, 99.0 (94.6-100); MET positive, 98.0 (89.0-99.9); MET negative 100 (92.8-100); RET positive, 93.8 (69.8-100); RET negative, 100 (94.9-100). The analytical performance showed that the panel could handle various types of biopsy samples obtained by routine clinical practice without requiring strict pathological monitoring, as in the case of conventional NGS panels.

Quantitative detection of ALK fusion breakpoints in plasma cell-free DNA from patients with non-small cell lung cancer using PCR-based target sequencing with a tiling primer set and two-step mapping/alignment.

BACKGROUND: Tyrosine kinase inhibitors targeted to anaplastic lymphoma kinase (ALK) have been demonstrated to be effective for lung cancer patients with an ALK fusion gene. Application of liquid biopsy, i.e., detection and quantitation of the fusion product in plasma cell-free DNA (cfDNA), could improve clinical practice. To detect ALK fusions, because fusion breakpoints occur somewhere in intron 19 of the ALK gene, sequencing of the entire intron is required to locate breakpoints. RESULTS: We constructed a target sequencing system using an adapter and a set of primers that cover the entire ALK intron 19. This system can amplify fragments, including breakpoints, regardless of fusion partners. The data analysis pipeline firstly detected fusions by alignment to selected target sequences, and then quantitated the fusion alleles aligning to the identified breakpoint sequences. Performance was validated using 20 cfDNA samples from ALK-positive non-small cell lung cancer patients and samples from 10 healthy volunteers. Sensitivity and specificity were 50 and 100%, respectively. CONCLUSIONS: We demonstrated that PCR-based target sequencing using a tiling primer set and two-step mapping/alignment quantitatively detected ALK fusions in cfDNA from lung cancer patients. The system offers an alternative to existing approaches based on hybridization capture.

Recent Advances in Liquid Biopsy in Precision Oncology Research.

Liquid biopsy is a minimally invasive test for cancer genetic status based on circulating tumor DNA (ctDNA), circulating tumor cells, or other tumor-derived materials in blood plasma. Although the minimal invasiveness and time resolution are attractive features of liquid biopsy, the limited amount of ctDNA in plasma poses problems. Recent developments in digital PCR and next-generation sequencing (NGS)-based technology have improved the accuracy of liquid biopsy. In particular, molecular barcoding technology in NGS-based methods, i.e., tagging of molecular barcodes to cell-free DNA before amplification, reduces technical errors by validating the consensus of sequences originating from a single molecule, leading to marked improvement of the accuracy and detection limit. However, substitutions caused by DNA damage and somatic mutations originating from normal cells are still obstacles to the sensitive detection of mutations on ctDNA. Since there have been only a few clinical applications, a deeper understanding of ctDNA biology and more advanced analytical technology are needed for the practical application of liquid biopsy.

High-throughput screening in colorectal cancer tissue-originated spheroids.

Patient-derived cancer organoid culture is an important live material that reflects clinical heterogeneity. However, the limited amount of organoids available for each case as well as the considerable amount of time and cost to expand in vitro makes it impractical to perform high-throughput drug screening using organoid cultures from multiple patients. Here, we report an advanced system for the high-throughput screening of 2427 drugs using the cancer tissue-originated spheroid (CTOS) method. In this system, we apply the CTOS method in an ex vivo platform from xenograft tumors, using machines to handle CTOS and reagents, and testing a CTOS reference panel of multiple CTOS lines for the hit drugs. CTOS passages in xenograft tumors resulted in minimal changes of morphological and genomic status, and xenograft tumor generation efficiently expanded the number of CTOS to evaluate multiple drugs. Our panel of colorectal cancer CTOS lines exhibited diverse sensitivities to the hit compounds, demonstrating the usefulness of this system for investigating highly heterogeneous disease.

Selective identification of somatic mutations in pancreatic cancer cells through a combination of next-generation sequencing of plasma DNA using molecular barcodes and a bioinformatic variant filter.

The accuracy of next-generation sequencing (NGS) for detecting tumor-specific mutations in plasma DNA is hindered by errors introduced during PCR/sequencing, base substitutions caused by DNA damage, and pre-existing mutations in normal cells that are present at a low frequency. Here, we performed NGS of genes related to pancreatic cancer (comprising 2.8 kb of genomic DNA) in plasma DNA (average 4.5 ng) using molecular barcodes. The average number of sequenced molecules was 900, and the sequencing depth per molecule was 100 or more. We also developed a bioinformatic variant filter, called CV78, to remove variants that were not considered to be tumor-specific, i.e., those that are either absent or occur at low frequencies in the Catalogue of Somatic Mutations in Cancer database. In a cohort comprising 57 pancreatic cancer patients and 12 healthy individuals, sequencing initially identified variants in 31 (54%) and 5 (42%), respectively, whereas after applying the CV78 filter, 19 (33%) and zero were variant-positive. In a validation cohort consisting of 86 patients with pancreatic cancer and 20 patients with intraductal papillary mucinous neoplasm (IPMN), 62 (72%) with pancreatic cancer patients and 10 (50%) IPMN patients were initially variant positive. After CV78 filtering, these values were reduced to 32 (37%) and 1 (5%), respectively. The variant allele frequency of filtered variants in plasma ranged from 0.25% to 76.1%. Therefore, combining NGS and molecular barcodes with subsequent filtering is likely to eliminate most non-tumor-specific mutations.

HLA genotyping by next-generation sequencing of complementary DNA.

BACKGROUND: Genotyping of the human leucocyte antigen (HLA) is indispensable for various medical treatments. However, unambiguous genotyping is technically challenging due to high polymorphism of the corresponding genomic region. Next-generation sequencing is changing the landscape of genotyping. In addition to high throughput of data, its additional advantage is that DNA templates are derived from single molecules, which is a strong merit for the phasing problem. Although most currently developed technologies use genomic DNA, use of cDNA could enable genotyping with reduced costs in data production and analysis. We thus developed an HLA genotyping system based on next-generation sequencing of cDNA. METHODS: Each HLA gene was divided into 3 or 4 target regions subjected to PCR amplification and subsequent sequencing with Ion Torrent PGM. The sequence data were then subjected to an automated analysis. The principle of the analysis was to construct candidate sequences generated from all possible combinations of variable bases and arrange them in decreasing order of the number of reads. Upon collecting candidate sequences from all target regions, 2 haplotypes were usually assigned. Cases not assigned 2 haplotypes were forwarded to 4 additional processes: selection of candidate sequences applying more stringent criteria, removal of artificial haplotypes, selection of candidate sequences with a relaxed threshold for sequence matching, and countermeasure for incomplete sequences in the HLA database. RESULTS: The genotyping system was evaluated using 30 samples; the overall accuracy was 97.0% at the field 3 level and 98.3% at the G group level. With one sample, genotyping of DPB1 was not completed due to short read size. We then developed a method for complete sequencing of individual molecules of the DPB1 gene, using the molecular barcode technology. CONCLUSION: The performance of the automatic genotyping system was comparable to that of systems developed in previous studies. Thus, next-generation sequencing of cDNA is a viable option for HLA genotyping.

Homozygous inactivation of CHEK2 is linked to a familial case of multiple primary lung cancer with accompanying cancers in other organs.

In clinical practice, there are a number of cancer patients with clear family histories, but the patients lack mutations in known familial cancer syndrome genes. Recent advances in genomic technologies have enhanced the possibility of identifying causative genes in such cases. Two siblings, an elder sister and a younger brother, were found to have multiple primary lung cancers at the age of 60. The former subsequently developed breast cancer and had a history of uterine myoma. The latter had initially developed prostate cancer at the age of 59 and had a history of colon cancer. Single-nucleotide polymorphism (SNP) genotyping revealed that ∼10% of the genomes were homozygous in both patients. Exome sequencing revealed nonsynonymous mutations in five genes in the runs of homozygosity: CHEK2, FCGRT, INPP5J, MYO18B, and SFI1. Evolutionary conservation of primary protein structures suggested the functional importance of the CHEK2 mutation, p.R474C. This mutation altered the tertiary structure of CHK2 by disrupting the salt bridge between p.R474 and p.E394. No such structural changes were observed with the other mutated genes. Subsequent cell-based transfection analysis revealed that CHK2 p.R474C was unstable and scarcely activated. We concluded that the homozygous CHEK2 variant was contributory in this case of familial cancer. Although homozygous inactivation of CHEK2 in mice led to cancers in multiple organs, accumulation of additional human cases is needed to establish its pathogenic role in humans.

Transient appearance of circulating tumor DNA associated with de novo treatment.

The limitation of circulating tumor DNA (ctDNA) is its inability to detect cancer cell subpopulations with few or no dying cells. Lung cancer patients subjected to the EGFR tyrosine kinase inhibitor (EGFR-TKI) treatment were prospectively collected, and ctDNA levels represented by the activating and T790M mutations were measured. The first data set (21 patients) consisting of samples collected in the period from before initiation of EGFR-TKI to at least 2 weeks after initiation: the ctDNA dynamics generally exhibited a rapid decrease and/or a transient increase. In 4 patients, we detected a transient increase of ctDNA bearing activating mutations not identified in biopsy samples. ctDNA with the same genotypical pattern was identified in 7 out of the 39 patients of the second data set intended to include samples until the onset of disease progression. In 6 of the 7 patients, this unique ctDNA appeared in the early period after treatment initiation, and did not reappear even after disease progression or chemotherapy. In another patient, similar ctDNA appeared upon radiation therapy. The identification of ctDNA with a unique genotype indicates the presence of cancer cell subpopulations that normally contain few or no dying cells, but generate dead cells because of the treatment.

Early responses of EGFR circulating tumor DNA to EGFR tyrosine kinase inhibitors in lung cancer treatment.

OBJECTIVES: Early evaluation of the effect of treatment is helpful in the management of cancer patients. Circulating biomarkers are an ideal tool for this if they are highly specific to tumors and respond rapidly to tumor volume changes. Circulating tumor DNA (ctDNA) is one such candidate. We conducted a prospective study to test the utility of EGFR ctDNA in early evaluation of EGFR-TKI effects. RESULTS: Twenty-one patients with EGFR-mutant lung cancer who were naïve to EGFR-TKI were enrolled. PM scores of EGFR ctDNA with activating mutations decreased rapidly in response to EGFR-TKI. Of the 14 patients with positive pretreatment PM scores, complete disappearance of major EGFR ctDNA was observed in 14.3%, 42.9%, and 57.1% on days 2 - 4, 8, and 15, respectively. These responses of EGFR ctDNA were most prominent among the measures used to evaluate responses, and correlated with early radiologic responses evaluated by chest X-rays. MATERIALS AND METHODS: EGFR ctDNA in serial plasma samples was amplified and 105 copies were sequenced with a next-generation sequencer. Plasma mutation (PM) score was defined as the number of reads containing deletions/substitutions in 105EGFR cell free DNA (cfDNA). When EGFR mutation in ctDNA was the same as that detected in cancer tissue, the ctDNA was defined as major EGFR ctDNA. CONCLUSIONS: The results indicate the usefulness of ctDNA as a highly specific biomarker for prediction of early response to treatment and that it can be applied to various types of cancer.

Numerical indices based on circulating tumor DNA for the evaluation of therapeutic response and disease progression in lung cancer patients.

Monitoring of disease/therapeutic conditions is an important application of circulating tumor DNA (ctDNA). We devised numerical indices, based on ctDNA dynamics, for therapeutic response and disease progression. 52 lung cancer patients subjected to the EGFR-TKI treatment were prospectively collected, and ctDNA levels represented by the activating and T790M mutations were measured using deep sequencing. Typically, ctDNA levels decreased sharply upon initiation of EGFR-TKI, however this did not occur in progressive disease (PD) cases. All 3 PD cases at initiation of EGFR-TKI were separated from other 27 cases in a two-dimensional space generated by the ratio of the ctDNA levels before and after therapy initiation (mutation allele ratio in therapy, MART) and the average ctDNA level. For responses to various agents after disease progression, PD/stable disease cases were separated from partial response cases using MART (accuracy, 94.7%; 95% CI, 73.5-100). For disease progression, the initiation of ctDNA elevation (initial positive point) was compared with the onset of objective disease progression. In 11 out of 28 eligible patients, both occurred within ±100 day range, suggesting a detection of the same change in disease condition. Our numerical indices have potential applicability in clinical practice, pending confirmation with designed prospective studies.

Characterization of the T-cell receptor beta chain repertoire in tumor-infiltrating lymphocytes.

Tumor-infiltrating lymphocytes (TILs) are direct effectors of tumor immunity, and their characterization is important for further development of immunotherapy. Recent advances in high-throughput sequencing technologies have enabled a comprehensive analysis of T-cell receptor (TCR) complementarity-determining region 3 (CDR3) sequences, which may provide information of therapeutic importance. We developed a high-fidelity target sequencing method with the ability for absolute quantitation, and performed large-scale sequencing of TCR beta chain (TCRB) CDR3 regions in TILs and peripheral blood lymphocytes (PBLs). The estimated TCRB repertoire sizes of PBLs from four healthy individuals and TILs from four colorectal cancer tissue samples were 608,664-1,003,098 and 90,228-223,757, respectively. The usage of J- and V-regions was similar in PBLs and TILs. Proportions of CDR3 amino acid (aa) sequences occupying more than 0.01% of the total molecular population were 0.33-0.43% in PBLs and 1.3-3.6% in TILs. Additional low coverage sequencing of 15 samples identified five CDR3 aa sequences that were shared by nine patients, one sequence shared by 10 patients, and one sequence shared by 12 patients. The estimated size of the TCRB repertoire in TILs was significantly smaller than that in PBLs. The proportion of abundant species (>0.01%) in TILs was larger than that in PBLs. Shared CDR3 aa sequences represent a response to common antigens, and the identification of such CDR3 sequences may be beneficial in developing clinical biomarkers.

Monitoring of treatment responses and clonal evolution of tumor cells by circulating tumor DNA of heterogeneous mutant EGFR genes in lung cancer.

OBJECTIVES: Epidermal growth factor receptor tyrosine kinase inhibitors (EGFR-TKIs) have dramatic effects on EGFR-mutant non-small-cell lung cancer (NSCLC). However, most patients experience disease recurrences, approximately half of which are T790M-mediated. Monitoring EGFR status with re-biopsy has spatiotemporal limitations. PATIENTS AND METHODS: EGFR circulating tumor DNA (ctDNA) in serial plasma samples was amplified and 10(5) of them were sequenced with a next-generation sequencer. Plasma mutation (PM) score was defined as the number of reads containing deletions/substitutions in 10(5)EGFR cell free DNA (cfDNA). RESULTS: PM scores of various EGFR mutations showed dynamic, case-specific changes during EGFR-TKI treatments in 52 patients. The effects of the treatment on EGFR ctDNA were evaluated in 38 patients with elevated pre-treatment PM scores. The ctDNA responses correlated well with radiologic responses in radiologic good responders, whereas correlation was poor in non-responders. In addition to the peaks for the most prevalent ctDNA, small peaks of ctDNA with different types of activating EGFR mutations or the T790M mutation (early T790M ctDNA) appeared transiently in 10.5% and 26.3%, respectively. Early T790M ctDNA disappeared in all patients, including 7 who eventually developed acquired resistance accompanied by elevated levels of T790M ctDNA. CONCLUSIONS: Monitoring ctDNA is useful in evaluating treatment responses and monitoring driver oncogene status in NSCLC. ctDNA revealed clonal heterogeneity and genetic processes of cancer evolution in individual patients. The simple presence of the T790M mutation may be insufficient to confer EGFR-TKI resistance to tumor cells.

Dynamics of circulating tumor DNA represented by the activating and resistant mutations in epidermal growth factor receptor tyrosine kinase inhibitor treatment.

Circulating tumor DNA (ctDNA) is an emerging field of cancer research. For lung cancer, non-invasive genotyping of EGFR is the foremost application. The activating mutations represent the ctDNA from all cancer cells, and the T790M-resistant mutation represents that from resistant cells. We examined the ctDNA dynamics of EGFR mutations by using deep sequencing with a massively parallel DNA sequencer. We obtained 190 plasma samples from 57 patients at various times during the treatment course and classified them according to treatment status. The mutation detection rate of exon 19 deletion/L858R in plasma was high at the initiation of treatment with epidermal growth factor receptor tyrosine kinase inhibitor (EGFR-TKI; P = 0.001), suppressed during EGFR-TKI treatment before disease progression, and elevated after the onset of disease progression (P = 0.023). The mutation detection rate of T790M was low until the onset of disease progression and elevated thereafter (P = 0.01). Samples across the development of disease progression were obtained from 10 patients and showed a correlation between increased ctDNA level and disease progression. Decreased ctDNA level in response to the initiation of EGFR-TKI was observed in 4 of 6 eligible patients. In two patients, the ctDNA dynamics suggested the presence of cancer cell populations only with the T790M mutation. In another patient, the T790M ctDNA represented cell subpopulations that respond to cytotoxic agents differently from the major population. Considering the high incidence, ctDNA could be a clinical parameter to complement information from image analyses.

High-fidelity target sequencing of individual molecules identified using barcode sequences: de novo detection and absolute quantitation of mutations in plasma cell-free DNA from cancer patients.

Circulating tumour DNA (ctDNA) is an emerging field of cancer research. However, current ctDNA analysis is usually restricted to one or a few mutation sites due to technical limitations. In the case of massively parallel DNA sequencers, the number of false positives caused by a high read error rate is a major problem. In addition, the final sequence reads do not represent the original DNA population due to the global amplification step during the template preparation. We established a high-fidelity target sequencing system of individual molecules identified in plasma cell-free DNA using barcode sequences; this system consists of the following two steps. (i) A novel target sequencing method that adds barcode sequences by adaptor ligation. This method uses linear amplification to eliminate the errors introduced during the early cycles of polymerase chain reaction. (ii) The monitoring and removal of erroneous barcode tags. This process involves the identification of individual molecules that have been sequenced and for which the number of mutations have been absolute quantitated. Using plasma cell-free DNA from patients with gastric or lung cancer, we demonstrated that the system achieved near complete elimination of false positives and enabled de novo detection and absolute quantitation of mutations in plasma cell-free DNA.

Diagnostic Accuracy of Noninvasive Genotyping of EGFR in Lung Cancer Patients by Deep Sequencing of Plasma Cell-Free DNA.

BACKGROUND: Genotyping of EGFR (epidermal growth factor receptor) mutations is indispensable for making therapeutic decisions regarding whether to use EGFR tyrosine kinase inhibitors (TKIs) for lung cancer. Because some cases might pose challenges for biopsy, noninvasive genotyping of EGFR in circulating tumor DNA (ctDNA) would be beneficial for lung cancer treatment. METHODS: We developed a detection system for EGFR mutations in ctDNA by use of deep sequencing of plasma DNA. Mutations were searched in >100 000 reads obtained from each exon region. Parameters corresponding to the limit of detection and limit of quantification were used as the thresholds for mutation detection. We conducted a multi-institute prospective study to evaluate the detection system, enrolling 288 non-small cell lung cancer (NSCLC) patients. RESULTS: In evaluating the performance of the detection system, we used the genotyping results from biopsy samples as a comparator: diagnostic sensitivity for exon 19 deletions, 50.9% (95% CI 37.9%-63.9%); diagnostic specificity for exon 19 deletions, 98.0% (88.5%-100%); sensitivity for the L858R mutation, 51.9% (38.7%-64.9%); and specificity for L858R, 94.1% (83.5%-98.6%). The overall sensitivities were as follows: all cases, 54.4% (44.8%-63.7%); stages IA-IIIA, 22.2% (11.5%-38.3%); and stages IIIB-IV, 72.7% (60.9%-82.1%). CONCLUSIONS: Deep sequencing of plasma DNA can be used for genotyping of EGFR in lung cancer patients. In particular, the high specificity of the system may enable a direct recommendation for EGFR-TKI on the basis of positive results with plasma DNA. Because sensitivity was low in early-stage NSCLC, the detection system is preferred for stage IIIB-IV NSCLC.

Analysis of ERBB ligand-induced resistance mechanism to crizotinib by primary culture of lung adenocarcinoma with EML4-ALK fusion gene.

INTRODUCTION: Using cell line-based assays, the secretion of erythroblastic leukemia viral oncogene homologue (ERBB) ligands has been reported to contribute to resistance against crizotinib in lung cancer with the echinoderm microtubule-associated protein-like 4 and anaplastic lymphoma kinase fusion gene. However, it is difficult to predict the role of the ligands in each patient. Here, we report an analysis of the mechanism of resistance behind crizotinib resistance using a primary culture of cancer cells from pleural effusion of an anaplastic lymphoma kinase-positive lung cancer patient who was clinically resistant to crizotinib. METHODS: Primary cancer cells were prepared as cancer tissue-originated spheroids (CTOSs) according to previously described methods. CTOSs were maintained in StemPro medium, and a sensitivity assay was performed under growth factor-free conditions, or under stimulation with epidermal growth factor (EGF) or neuregulin-1/heregulin. The effect of treatment with crizotinib alone or a combination of crizotinib and erlotinib was examined. RESULTS: Cancer cells (LB53) were established to be CTOSs from a patient who was clinically resistant to crizotinib. The CTOSs were sensitive to crizotinib under growth factor-free conditions in vitro, whereas resistant under stimulation with EGF or neuregulin-1. These ligands rescued the inhibition of intracellular signaling by crizotinib. Pleural effusion from the patient also activated EGF receptor signaling to the similar extent of EGF stimulation. The resistance to crizotinib by EGF was reversed by blocking EGF receptor signaling by erlotinib in vitro. CONCLUSION: Stimulation by ERBB ligands is suggested to be responsible for resistance to crizotinib in this patient. The CTOS method may enable analysis of resistance mechanism for targeted therapy in individual patients.

Prognostic prediction of glioblastoma by quantitative assessment of the methylation status of the entire MGMT promoter region.

BACKGROUND: O6-methylguanine-DNA methyltransferase (MGMT) promoter methylation is reported to be a prognostic and predictive factor of alkylating chemotherapy for glioblastoma patients. Methylation specific PCR (MSP) has been most commonly used when the methylation status of MGMT is assessed. However, technical obstacles have hampered the implementation of MSP-based diagnostic tests. We quantitatively analyzed the methylation status of the entire MGMT promoter region and applied this information for prognostic prediction using sequencing technology. METHODS: Between 1998 and 2012, the genomic DNA of 85 tumor samples from newly diagnosed glioblastoma patients was subjected to bisulfite treatment and subdivided into a training set, consisting of fifty-three samples, and a test set, consisting of thirty-two samples. The training set was analyzed by deep Sanger sequencing with a sequencing coverage of up to 96 clones per sample. This analysis quantitatively revealed the degree of methylation of each cytidine phosphate guanosine (CpG) site. Based on these data, we constructed a prognostic prediction system for glioblastoma patients using a supervised learning method. We then validated this prediction system by deep sequencing with a next-generation sequencer using a test set of 32 samples. RESULTS: The methylation status of the MGMT promoter was correlated with progression-free survival (PFS) in our patient population in the training set. The degree of correlation differed among the CpG sites. Using the data from the top twenty CpG sites, we constructed a prediction system for overall survival (OS) and PFS. The system successfully classified patients into good and poor prognosis groups in both the training set (OS, p = 0.0381; PFS, p = 0.00122) and the test set (OS, p = 0.0476; PFS, p = 0.0376). Conventional MSP could not predict the prognosis in either of our sets. (training set: OS; p = 0.993 PFS; p = 0.113, test set: OS; p = 0.326 PFS; p = 0.342). CONCLUSIONS: The prognostic ability of our prediction system using sequencing data was better than that of methylation-specific PCR (MSP). Advances in sequencing technologies will make this approach a plausible option for diagnoses based on MGMT promotor methylation.

Quantitative identification of mutant alleles derived from lung cancer in plasma cell-free DNA via anomaly detection using deep sequencing data.

The detection of rare mutants using next generation sequencing has considerable potential for diagnostic applications. Detecting circulating tumor DNA is the foremost application of this approach. The major obstacle to its use is the high read error rate of next-generation sequencers. Rather than increasing the accuracy of final sequences, we detected rare mutations using a semiconductor sequencer and a set of anomaly detection criteria based on a statistical model of the read error rate at each error position. Statistical models were deduced from sequence data from normal samples. We detected epidermal growth factor receptor (EGFR) mutations in the plasma DNA of lung cancer patients. Single-pass deep sequencing (>100,000 reads) was able to detect one activating mutant allele in 10,000 normal alleles. We confirmed the method using 22 prospective and 155 retrospective samples, mostly consisting of DNA purified from plasma. A temporal analysis suggested potential applications for disease management and for therapeutic decision making to select epidermal growth factor receptor tyrosine kinase inhibitors (EGFR-TKI).