Simple, multiplexed, PCR-based barcoding of DNA enables sensitive mutation detection in liquid biopsies using sequencing.

Detection of cell-free DNA in liquid biopsies offers great potential for use in non-invasive prenatal testing and as a cancer biomarker. Fetal and tumor DNA fractions however can be extremely low in these samples and ultra-sensitive methods are required for their detection. Here, we report an extremely simple and fast method for introduction of barcodes into DNA libraries made from 5 ng of DNA. Barcoded adapter primers are designed with an oligonucleotide hairpin structure to protect the molecular barcodes during the first rounds of polymerase chain reaction (PCR) and prevent them from participating in mis-priming events. Our approach enables high-level multiplexing and next-generation sequencing library construction with flexible library content. We show that uniform libraries of 1-, 5-, 13- and 31-plex can be generated. Utilizing the barcodes to generate consensus reads for each original DNA molecule reduces background sequencing noise and allows detection of variant alleles below 0.1% frequency in clonal cell line DNA and in cell-free plasma DNA. Thus, our approach bridges the gap between the highly sensitive but specific capabilities of digital PCR, which only allows a limited number of variants to be analyzed, with the broad target capability of next-generation sequencing which traditionally lacks the sensitivity to detect rare variants.

Comparison of cancer-associated genetic abnormalities in columnar-lined esophagus tissues with and without goblet cells.

OBJECTIVE: To determine and compare the frequency of cancer-associated genetic abnormalities in esophageal metaplasia biopsies with and without goblet cells. BACKGROUND: Barrett's esophagus is associated with increased risk of esophageal adenocarcinoma (EAC), but the appropriate histologic definition of Barrett's esophagus is debated. Intestinal metaplasia (IM) is defined by the presence of goblet cells whereas nongoblet cell metaplasia (NGM) lacks goblet cells. Both have been implicated in EAC risk but this is controversial. Although IM is known to harbor genetic changes associated with EAC, little is known about NGM. We hypothesized that if NGM and IM infer similar EAC risk, then they would harbor similar genetic aberrations in genes associated with EAC. METHODS: Ninety frozen NGM, IM, and normal tissues from 45 subjects were studied. DNA copy number abnormalities were identified using microarrays and fluorescence in situ hybridization. Targeted sequencing of all exons from 20 EAC-associated genes was performed on metaplasia biopsies using Ion AmpliSeq DNA sequencing. RESULTS: Frequent copy number abnormalities targeting cancer-associated genes were found in IM whereas no such changes were observed in NGM. In 1 subject, fluorescence in situ hybridization confirmed loss of CDKN2A and amplification of chromosome 8 in IM but not in a nearby NGM biopsy. Targeted sequencing revealed 11 nonsynonymous mutations in 16 IM samples and 2 mutations in 19 NGM samples. CONCLUSIONS: This study reports the largest and most comprehensive comparison of DNA aberrations in IM and NGM genomes. Our results show that IM has a much higher frequency of cancer-associated mutations than NGM.

Cyclin E involved in early stage carcinogenesis of esophageal adenocarcinoma by SNP DNA microarray and immunohistochemical studies.

BACKGROUND: Cyclin E is a cell cycle regulator which is critical for driving G1/S transition. Abnormal levels of cyclin E have been found in many cancers. However, the level changes of cyclin E in esophageal adenocarcinoma and its precancerous lesion have not been well studied. Here, we focus on the gene amplification and expression of cyclin E in these lesions, and aim to ascertain the relationship with clinicopathological characteristics. METHODS: Genomic DNA was analyzed from 116 esophageal adenocarcinoma and 26 precancerous lesion patients using Affymetrix SNP 6.0 arrays. The protein overexpression of cyclin E was also detected using immunohistochemistry from tissue microarrays containing esophageal adenocarcinoma and precancerous lesions. Patient survival and other clinical data were collected and analyzed. The intensity and percentage of the cyclin E expressing cells in tissue microarrays were scored by two pathologists. Fisher exact tests and Kaplan-Meier methods were used to analyze data. RESULTS: By genomic analysis, cyclin E was amplified in 19.0% of the EAC samples. By immunohistochemistry, high expression of cyclin E was observed in 2.3% of squamous mucosa tissues, 3.7% in columnar cell metaplasia, 5.8% in Barrett's esophagus, 19.0% in low grade dysplasia, 35.7% in high grade dysplasia, and 16.7% in esophageal adenocarcinoma. The differences in cyclin E high expression between neoplastic groups and non-dysplasia groups are statistically significant (p < 0.05). The prognosis for patients with high cyclin E expression appeared slightly better than for those with low cyclin E expression although this was not statistically significant (p = 0.13). CONCLUSIONS: The expression of cyclin E significantly increases from non-dysplasia esophageal lesion to low and high grade dysplasia, suggesting that cyclin E plays an important role in the early stage of carcinogenesis. Importantly, cyclin E is also amplified and highly expressed in a subset of esophageal adenocarcinoma patients, but this increase is not associated with worse prognosis.

Bile exposure inhibits expression of squamous differentiation genes in human esophageal epithelial cells.

OBJECTIVE: This study aimed to identify pathways and cellular processes that are modulated by exposure of normal esophageal cells to bile and acid. BACKGROUND: Barrett's esophagus most likely develops as a response of esophageal stem cells to the abnormal reflux environment. Although insights into the underlying molecular mechanisms are slowly emerging, much of the metaplastic process remains unknown. METHODS: We performed a global analysis of gene expression in normal squamous esophageal cells in response to bile or acid exposure. Differentially expressed genes were classified into major biological functions using pathway analysis and interaction network software. Array data were verified by quantitative PCR and western blot both in vitro and in human esophageal biopsies. RESULTS: Bile modulated expression of 202 genes, and acid modulated expression of 103 genes. Genes involved in squamous differentiation formed the largest functional group (n = 45) all of which were downregulated by bile exposure. This included genes such as involucrin (IVL), keratinocyte differentiation-associated protein (KRTDAP), grainyhead-like 1 (GRHL1), and desmoglein1 (DSG1) the downregulation of which was confirmed by quantitative PCR and western blot. Bile also caused expression changes in genes involved in cell adhesion, DNA repair, oxidative stress, cell cycle, Wnt signaling, and lipid metabolism. Analysis of human esophageal biopsies demonstrated greatly reduced expression of IVL, KRTDAP, DSG1, and GRHL1 in metaplastic compared to squamous epithelia. CONCLUSIONS: We report for the first time that bile inhibits the squamous differentiation program of esophageal epithelial cells. This, coordinated with induction of genes driving intestinal differentiation, may be required for the development of Barrett's esophagus.

Clinicopathologic characteristics of high expression of Bmi-1 in esophageal adenocarcinoma and squamous cell carcinoma.

BACKGROUND: High expression of Bmi-1, a key regulatory component of the polycomb repressive complex-1, has been associated with many solid and hematologic malignancies including esophageal squamous cell carcinoma. However, little is known about the role of Bmi-1 in esophageal adenocarcinoma. The aim of this study is to investigate the amplification and high expression of Bmi-1 and the associated clinicopathologic characteristics in esophageal adenocarcinoma and squamous cell carcinoma. METHODS: The protein expression level of Bmi-1 was detected by immunohistochemistry (IHC) from tissue microarrays (TMA) constructed at the University of Rochester from using tissues accrued between 1997 and 2005. Types of tissues included adenocarcinoma, squamous cell carcinoma and precancerous lesions. Patients' survival data, demographics, histologic diagnoses and tumor staging data were collected. The intensity (0-3) and percentage of Bmi-1 expression on TMA slides were scored by two pathologists. Genomic DNA from 116 esophageal adenocarcinoma was analyzed for copy number aberrations using Affymetrix SNP 6.0 arrays. Fisher exact tests and Kaplan-Meier methods were used to analyze data. RESULTS: By IHC, Bmi-1 was focally expressed in the basal layers of almost all esophageal squamous mucosa, which was similar to previous reports in other organs related to stem cells. High Bmi-1 expression significantly increased from squamous epithelium (7%), columnar cell metaplasia (22%), Barrett's esophagus (22%), to low- (45%) and high-grade dysplasia (43%) and adenocarcinoma (37%). The expression level of Bmi-1 was significantly associated with esophageal adenocarcinoma differentiation. In esophageal adenocarcinoma, Bmi-1 amplification was detected by DNA microarray in a low percentage (3%). However, high Bmi-1 expression did not show an association with overall survival in both esophageal adenocarcinoma and squamous cell carcinoma. CONCLUSIONS: This study demonstrates that high expression Bmi-1 is associated with esophageal adenocarcinoma and precancerous lesions, which implies that Bmi-1 plays an important role in early carcinogenesis in esophageal adenocarcinoma.

HER2 amplification, overexpression and score criteria in esophageal adenocarcinoma.

The HER2 oncogene was recently reported to be amplified and overexpressed in esophageal adenocarcinoma. However, the relationship of HER2 amplification in esophageal adenocarcinoma with prognosis has not been well defined. The scoring systems for clinically evaluating HER2 in esophageal adenocarcinoma are not established. The aims of the study were to establish a HER2 scoring system and comprehensively investigate HER2 amplification and overexpression in esophageal adenocarcinoma and its precursor lesion. Using a tissue microarray, containing 116 cases of esophageal adenocarcinoma, 34 cases of Barrett's esophagus, 18 cases of low-grade dysplasia and 15 cases of high-grade dysplasia, HER2 amplification and overexpression were analyzed by HercepTest and chromogenic in situ hybridization methods. The amplification frequency in an independent series of 116 esophageal adenocarcinoma samples was also analyzed using Affymetrix SNP 6.0 microarrays. In our studies, we have found that HER2 amplification does not associate with poor prognosis in total 232 esophageal adenocarcinoma patients by chromogenic in situ hybridization and high-density microarrays. We further confirm the similar frequency of HER2 amplification by chromogenic in situ hybridization (18%; 21 out of 116) and SNP 6.0 microarrays (16%, 19 out of 116) in esophageal adenocarcinoma. HER2 protein overexpression was observed in 12% (14 out of 116) of esophageal adenocarcinoma and 7% (1 out of 15) of high-grade dysplasia. No HER2 amplification or overexpression was identified in Barrett's esophagus or low-grade dysplasia. All HER2 protein overexpression cases showed HER2 gene amplification. Gene amplification was found to be more frequent by chromogenic in situ hybridization than protein overexpression in esophageal adenocarcinoma (18 vs 12%). A modified two-step model for esophageal adenocarcinoma HER2 testing is recommended for clinical esophageal adenocarcinoma HER2 trial.

Whole genome exon arrays identify differential expression of alternatively spliced, cancer-related genes in lung cancer.

Alternative processing of pre-mRNA transcripts is a major source of protein diversity in eukaryotes and has been implicated in several disease processes including cancer. In this study we have performed a genome wide analysis of alternative splicing events in lung adenocarcinoma. We found that 2369 of the 17 800 core Refseq genes appear to have alternative transcripts that are differentially expressed in lung adenocarcinoma versus normal. According to their known functions the largest subset of these genes (30.8%) is believed to be cancer related. Detailed analysis was performed for several genes using PCR, quantitative RT-PCR and DNA sequencing. We found overexpression of ERG variant 2 but not variant 1 in lung tumors and overexpression of CEACAM1 variant 1 but not variant 2 in lung tumors but not in breast or colon tumors. We also identified a novel, overexpressed variant of CDH3 and verified the existence and overexpression of a novel variant of P16 transcribed from the CDKN2A locus. These findings demonstrate how analysis of alternative pre-mRNA processing can shed additional light on differences between tumors and normal tissues as well as between different tumor types. Such studies may lead to the development of additional tools for tumor diagnosis, prognosis and therapy.

Molecular staging of pathologically negative sentinel lymph nodes from melanoma patients using multimarker, quantitative real-time rt-PCR.

The aim of this study was to evaluate the prognostic potential of quantitative reverse-transcription, polymerase chain reaction (qRT-PCR) in melanoma patients with pathologically negative sentinel lymph nodes (SLN). Our study included 195 node-negative melanoma patients with a Breslow thickness greater than 0.76 mm (n = 158), or less than 0.76 mm but who had Clark level IV-V, microscopic ulceration, or pathological signs of regression (n = 32), and five patients with melanoma of unknown thickness. SLNs were examined by serial-section histopathology. A portion of each SLN was frozen for qRT-PCR analysis using markers Tyrosinase, MART1, SSX2, MAGEA3, PAX3, and GalNAc-T. In addition, two other markers (PLAB and L1CAM) were evaluated for melanoma specificity but not for SLN analysis. Median follow-up was 64 months, during which time there were 15 (7.7%) recurrences. A total of 370 lymph nodes were analyzed by qRT-PCR. No association was found between quantitative expression level of any marker and disease recurrence. Previously published primer designs were tested for PAX3 and GalNAc-T and revealed that alternative PAX3 transcripts are differentially expressed in melanoma and benign lymph nodes. No associations with recurrence were found regardless of the transcripts amplified by different primer sets. PLAB and L1CAM did not appear to differentiate between malignant melanoma and benign melanocytes or lymph nodes in our analysis. We conclude that, in this large cohort of patients, multimarker qRT-PCR analysis of SLNs did not correlate with disease recurrence. Our data support specific PAX3 splice variants but not GalNAc-T, PLAB or L1CAM as possible markers for melanoma metastasis to SLNs.

Association between CCND1 G/A870 polymorphism, allele-specific amplification, cyclin D1 expression, and survival in esophageal and lung carcinoma.

PURPOSE: Cyclin D1 is found on 11q13, which is a region frequently amplified in several tumor types. The CCND1 locus gives rise to at least two protein isoforms of D1 (D1a and D1b). A common G/A polymorphism (G/A870) is thought to influence the expression levels of D1a and D1b. D1b has been suggested to be increased in the presence of the A allele and more oncogenic than D1a. Furthermore, the A allele has been reported to correlate with increased risk of carcinoma in several tumor types, suggesting that this polymorphism and D1b are important in tumor progression. However, contradictory data about the polymorphism, D1 variant expression, and correlation with survival have been reported. We explored the relationship between gene amplification, G/A870 genotype, D1a and D1b expression, and overall survival in esophageal adenocarcinoma and non-small cell lung cancer. EXPERIMENTAL DESIGN: DNA and RNA were isolated from 54 esophageal adenocarcinoma samples and 89 non-small cell lung cancer samples and were analyzed for gene amplification, genotype at the polymorphism, gene expression, and association with overall survival. RESULTS: The D1 variant expression did not correlate with amplification, genotype, or overall survival in either tumor type. The total D1 expression correlated with decreased patient survival. Several other genes on 11q13 also seem to be overexpressed and correlated with decreased survival. CONCLUSIONS: We report that the G/A870 polymorphism does not correlate with patient survival, or with D1a or D1b expression. However, the total D1 expression and the expression of several other genes on 11q13 seem to be associated with esophageal adenocarcinoma patient survival.

The Molecular Biologic Basis of Esophageal and Gastric Cancers.

Esophageal cancer and gastric cancer are leading causes of cancer-related mortality worldwide. In this article, the authors discuss the molecular biology of esophageal and gastric cancer with a focus on esophageal adenocarcinoma. They review data from The Cancer Genome Atlas project and advances in the molecular stratification and classification of esophageal carcinoma and gastric cancer. They also summarize advances in microRNA, molecular staging, gene expression profiling, tumor microenvironment, and detection of circulating tumor DNA. Finally, the authors summarize some of the implications of understanding the molecular basis of esophageal cancer and future directions in the management of esophageal cancer.

Impact of Polymerase Fidelity on Background Error Rates in Next-Generation Sequencing with Unique Molecular Identifiers/Barcodes.

Liquid biopsy and detection of tumor-associated mutations in cell-free circulating DNA often requires the ability to identify single nucleotide variants at allele frequencies below 0.1%. Standard sequencing protocols cannot achieve this level of sensitivity due to background noise from DNA damage and polymerase induced errors. Addition of unique molecular identifiers allows identification and removal of errors responsible for this background noise. Theoretically, high fidelity enzymes will also reduce error rates in barcoded NGS but this has not been thoroughly explored. We evaluated the impact of polymerase fidelity on the magnitude of error reduction at different steps of barcoded NGS library construction. We find that barcoding itself displays largest impact on error reduction, even with low fidelity polymerases. Use of high fidelity polymerases in the barcoding step of library construction further suppresses error in barcoded NGS, and allows detection of variant alleles below 0.1% allele frequency. However, the improvement in error correction is modest and is not directly proportional to polymerase fidelity. Depending on the specific application, other polymerase characteristics such as multiplexing capacity, PCR efficiency, buffer requirements and ability to amplify targets with high GC content may outweigh the relatively small additional decrease in error afforded by ultra-high fidelity polymerases.

Considerations and quality controls when analyzing cell-free tumor DNA.

Circulating cell-free tumor DNA (ctDNA) is a promising biomarker in cancer. Ultrasensitive technologies enable detection of low (< 0.1%) mutant allele frequencies, a pre-requisite to fully utilize the potential of ctDNA in cancer diagnostics. In addition, the entire liquid biopsy workflow needs to be carefully optimized to enable reliable ctDNA analysis. Here, we discuss important considerations for ctDNA detection in plasma. We show how each experimental step can easily be evaluated using simple quantitative PCR assays, including detection of cellular DNA contamination and PCR inhibition. Furthermore, ctDNA assay performance is also demonstrated to be affected by both DNA fragmentation and target sequence. Finally, we show that quantitative PCR is useful to estimate the required sequencing depth and to monitor DNA losses throughout the workflow. The use of quality control assays enables the development of robust and standardized workflows that facilitate the implementation of ctDNA analysis into clinical routine.

Corrigendum to "Considerations and quality controls when analyzing cell-free tumor DNA" [Biomol. Detect. Quantif. 17 (2019) 100078].

[This corrects the article DOI: 10.1016/j.bdq.2018.12.003.].

Plasma circulating tumor DNA as a potential tool for disease monitoring in head and neck cancer.

BACKGROUND: Recommendations for perioperative therapy in head and neck cancer are not explicit and recurrence occurs frequently. Circulating tumor DNA is an emerging cancer biomarker, but has not been extensively explored for detection of recurrence in head and neck cancer. METHODS: Patients diagnosed with head and neck squamous cell carcinoma were recruited into the study protocol. Tumors were sequenced to identify patient-specific mutations. Mutations were then identified in plasma circulating tumor DNA from pre-treatment blood samples and longitudinally during standard follow-up. Circulating tumor DNA status during follow-up was correlated to disease recurrence. RESULTS: Samples were taken from eight patients. Tumor mutations were verified in seven patients. Baseline circulating tumor DNA was positive in six patients. Recurrence occurred in four patients, two of whom had detectable circulating tumor DNA prior to recurrence. CONCLUSION: Circulating tumor DNA is a potential tool for disease and recurrence monitoring following curative therapy in head and neck cancer, allowing for better prognostication, and/or modification of treatment strategies.

Combined EsophaCap cytology and MUC2 immunohistochemistry for screening of intestinal metaplasia, dysplasia and carcinoma.

Purpose: The incidence of esophageal adenocarcinoma (EAC) has increased by 700% in Western countries over the last 30 years. Although clinical guidelines call for endoscopic surveillance for EAC among high-risk populations, fewer than 5% of new EAC patients are under surveillance at the time of diagnosis. We studied the accuracy of combined cytopathology and MUC2 immunohistochemistry (IHC) for screening of Intestinal Metaplasia (IM), dysplasia and EAC, using specimens collected from the EsophaCap swallowable encapsulated cytology sponge from Canada and United States. Patients and methods: By comparing the EsophaCap cytological diagnosis with concurrent endoscopic biopsies performed on the same patients in 28 cases, we first built up the cytology diagnostic categories and criteria. Based on these criteria, 136 cases were evaluated by both cytology and MUC2 IHC with blinded to patient biopsy diagnosis. Results: We first set up categories and criteria for cytological diagnosis of EscophaCap samples. Based on these, we divided our evaluated cytological samples into two groups: non-IM group and IM or dysplasia or adenocarcinoma group. Using the biopsy as our gold standard to screen IM, dysplasia and EAC by combined cytology and MUC2 IHC, the sensitivity and specificity were 68% and 91%, respectively, which is in the range of clinically useful cytological screening tests such as the cervical Pap smear. Conclusions: Combined EsophaCap cytology and MUC2 IHC could be a good screening test for IM and Beyond.

Is occult lymph node disease in colorectal cancer patients clinically significant? A review of the relevant literature.

The clinical significance of micrometastasis of colorectal cancer (CRC) to regional lymph nodes remains controversial. In this review, we analyze publications that have evaluated the clinical significance of occult lymph node metastasis in CRC. An extensive literature search identified 19 publications that evaluated the clinical significance of micrometastatic CRC by various methods, including immunohistochemistry (IHC; n = 13) and reverse transcription-polymerase chain reaction (RT-PCR, n = 6). These studies were reviewed for methodology and findings. Significant limitations in methodology were identified, including inconsistent histological definitions of micrometastatic disease, poor sampling because of an inadequate number of lymph nodes or number of sections per lymph node analyzed, lack of conformity with respect to IHC antibody or RT-PCR marker, and inadequate power because of small sample size. Micrometastatic lymph node metastasis identified by RT-PCR was consistently found to be prognostically significant, but this was not true of micrometastatic disease identified by IHC. RT-PCR analysis of lymph nodes with specific markers can help identify pN0 (pathological-negative lymph node) CRC patients at increased risk for recurrence. The identification of occult disease by IHC techniques may also ultimately prove to be associated with worse outcome, but a number of inadequately powered studies have concluded conversely.

A combination of molecular markers accurately detects lymph node metastasis in non-small cell lung cancer patients.

Occult lymph node metastasis (micrometastasis) is a good prognostic indicator in non-small cell lung cancer (NSCLC) and could be used to direct adjuvant chemotherapy in stage I patients. This study was designed to evaluate molecular markers for detection of occult lymph node metastasis in NSCLC, define the best marker or marker combination to distinguish positive from benign lymph nodes, and evaluate these markers in lymph nodes from pathologically node-negative (pN(0)) NSCLC patients. Potential markers were identified through literature and database searches and all markers were analyzed by quantitative reverse transcription-PCR in a primary screen of six NSCLC specimens and 10 benign nodes. Selected markers were further evaluated on 21 primary NSCLC specimens, 21 positive nodes, and 21 benign nodes, and the best individual markers and combinations were identified. A combination of three markers was further validated on an independent set of 32 benign lymph nodes, 38 histologically positive lymph nodes, and 462 lymph nodes from 68 pN(0) NSCLC patients. Forty-two markers were evaluated in the primary screen and eight promising markers were selected for further analysis. A combination of three markers (SFTPB, TACSTD1, and PVA) was identified that provided perfect classification of benign and positive nodes in all sample sets. PVA and SFTPB are particularly powerful in tumors of squamous and adenocarcinoma histologies, respectively, whereas TACSTD1 is a good general marker for NSCLC metastasis. The combination of these genes identified 32 of 462 (7%) lymph nodes from 20 of 68 (29%) patients as potentially positive for occult metastasis. Long-term follow-up will determine the clinical relevance of these findings.

Molecular staging of cervical lymph nodes in squamous cell carcinoma of the head and neck.

Clinical staging of cervical lymph nodes from patients with squamous cell carcinoma of the head and neck (SCCHN) has only 50% accuracy compared with definitive pathologic assessment. Consequently, both clinically positive and clinically negative patients frequently undergo neck dissections that may not be necessary. To address this potential overtreatment, sentinel lymph node (SLN) biopsy is currently being evaluated to provide better staging of the neck. However, to fully realize the potential improvement in patient care afforded by the SLN procedure, a rapid and accurate SLN analysis is necessary. We used quantitative reverse transcription-PCR (QRT-PCR) to screen 40 potential markers for their ability to detect SCCHN metastases to cervical lymph nodes. Seven markers were identified with good characteristics for identifying metastatic disease, and these were validated using a set of 26 primary tumors, 19 histologically positive lymph nodes, and 21 benign nodes from patients without cancer. Four markers discriminated between positive and benign nodes with accuracy >97% but only one marker, pemphigus vulgaris antigen (PVA), discriminated with 100% accuracy in both the observed data and a statistical bootstrap analysis. A rapid QRT-PCR assay for PVA was then developed and incorporated into a prototype instrument capable of performing fully automated RNA isolation and QRT-PCR. The automated analysis with PVA provided perfect discrimination between histologically positive and benign lymph nodes and correctly identified two lymph nodes with micrometastatic tumor deposits. These assays were completed (from tissue to result) in approximately 30 minutes, thus demonstrating the feasibility of intraoperative staging of SCCHN SLNs by QRT-PCR.

Simple multiplexed PCR-based barcoding of DNA for ultrasensitive mutation detection by next-generation sequencing.

Detection of extremely rare variant alleles within a complex mixture of DNA molecules is becoming increasingly relevant in many areas of clinical and basic research, such as the detection of circulating tumor DNA in the plasma of cancer patients. Barcoding of DNA template molecules early in next-generation sequencing (NGS) library construction provides a way to identify and bioinformatically remove polymerase errors that otherwise make detection of these rare variants very difficult. Several barcoding strategies have been reported, but all require long and complex library preparation protocols. Simple, multiplexed, PCR-based barcoding of DNA for sensitive mutation detection using sequencing (SiMSen-seq) was developed to generate targeted barcoded libraries with minimal DNA input, flexible target selection and a very simple, short (∼4 h) library construction protocol. The protocol comprises a three-cycle barcoding PCR step followed directly by adaptor PCR to generate the library and then bead purification before sequencing. Thus, SiMSen-seq allows detection of variant alleles at <0.1% frequency with easy customization of library content (from 1 to 40+ PCR amplicons) and a protocol that can be implemented in any molecular biology laboratory. Here, we provide a detailed protocol for assay development and describe software to process the barcoded sequence reads.