Ultra-short cell-free DNA fragments enhance cancer early detection in a multi-analyte blood test combining mutation, protein and fragmentomics.

OBJECTIVES: Cancer morbidity and mortality can be reduced if the cancer is detected early. Cell-free DNA (cfDNA) fragmentomics emerged as a novel epigenetic biomarker for early cancer detection, however, it is still at its infancy and requires technical improvement. We sought to apply a single-strand DNA sequencing technology, for measuring genetic and fragmentomic features of cfDNA and evaluate the performance in detecting multiple cancers. METHODS: Blood samples of 364 patients from six cancer types (colorectal, esophageal, gastric, liver, lung, and ovarian cancers) and 675 healthy individuals were included in this study. Circulating tumor DNA mutations, cfDNA fragmentomic features and a set of protein biomarkers were assayed. Sensitivity and specificity were reported by cancer types and stages. RESULTS: Circular Ligation Amplification and sequencing (CLAmp-seq), a single-strand DNA sequencing technology, yielded a population of ultra-short fragments (<100 bp) than double-strand DNA preparation protocols and reveals a more significant size difference between cancer and healthy cfDNA fragments (25.84 bp vs. 16.05 bp). Analysis of the subnucleosomal peaks in ultra-short cfDNA fragments indicates that these peaks are regulatory element "footprints" and correlates with gene expression and cancer stages. At 98 % specificity, a prediction model using ctDNA mutations alone showed an overall sensitivity of 46 %; sensitivity reaches 60 % when protein is added, sensitivity further increases to 66 % when fragmentomics is also integrated. More improvements observed for samples representing earlier cancer stages than later ones. CONCLUSIONS: These results suggest synergistic properties of protein, genetic and fragmentomics features in the identification of early-stage cancers.

Antigen presentation profiling reveals recognition of lymphoma immunoglobulin neoantigens.

Cancer somatic mutations can generate neoantigens that distinguish malignant from normal cells. However, the personalized identification and validation of neoantigens remains a major challenge. Here we discover neoantigens in human mantle-cell lymphomas by using an integrated genomic and proteomic strategy that interrogates tumour antigen peptides presented by major histocompatibility complex (MHC) class I and class II molecules. We applied this approach to systematically characterize MHC ligands from 17 patients. Remarkably, all discovered neoantigenic peptides were exclusively derived from the lymphoma immunoglobulin heavy- or light-chain variable regions. Although we identified MHC presentation of private polymorphic germline alleles, no mutated peptides were recovered from non-immunoglobulin somatically mutated genes. Somatic mutations within the immunoglobulin variable region were almost exclusively presented by MHC class II. We isolated circulating CD4+ T cells specific for immunoglobulin-derived neoantigens and found these cells could mediate killing of autologous lymphoma cells. These results demonstrate that an integrative approach combining MHC isolation, peptide identification, and exome sequencing is an effective platform to uncover tumour neoantigens. Application of this strategy to human lymphoma implicates immunoglobulin neoantigens as targets for lymphoma immunotherapy.

Minimal residual disease negativity using deep sequencing is a major prognostic factor in multiple myeloma.

The introduction of novel agents has led to major improvements in clinical outcomes for patients with multiple myeloma. To shorten evaluation times for new treatments, health agencies are currently examining minimal residual disease (MRD) as a surrogate end point in clinical trials. We assessed the prognostic value of MRD, measured during maintenance therapy by next-generation sequencing (NGS). MRD negativity was defined as the absence of tumor plasma cell within 1 000 000 bone marrow cells (<10-6). Data were analyzed from a recent clinical trial that evaluated the role of transplantation in newly diagnosed myeloma patients treated with lenalidomide, bortezomib, and dexamethasone (RVD). MRD negativity was achieved at least once during maintenance in 127 patients (25%). At the start of maintenance therapy, MRD was a strong prognostic factor for both progression-free survival (adjusted hazard ratio, 0.22; 95% confidence interval, 0.15-0.34; P < .001) and overall survival (adjusted hazard ratio, 0.24; 95% confidence interval, 0.11-0.54; P = .001). Patients who were MRD negative had a higher probability of prolonged progression-free survival than patients with detectable residual disease, regardless of treatment group (RVD vs transplant), cytogenetic risk profile, or International Staging System disease stage at diagnosis. These results were similar after completion of maintenance therapy. Our findings confirm the value of MRD status, as determined by NGS, as a prognostic biomarker in multiple myeloma, and suggest that this approach could be used to adapt treatment strategies in future clinical trials.

RB but not R-HCVAD is a feasible induction regimen prior to auto-HCT in frontline MCL: results of SWOG Study S1106.

Aggressive induction chemotherapy followed by autologous haematopoietic stem cell transplant (auto-HCT) is effective for younger patients with mantle cell lymphoma (MCL). However, the optimal induction regimen is widely debated. The Southwestern Oncology Group S1106 trial was designed to assess rituximab plus hyperCVAD/MTX/ARAC (hyperfractionated cyclophosphamide, vincristine, doxorubicin and dexamethasone, alternating with high dose cytarabine and methotrexate) (RH) versus rituximab plus bendamustine (RB) in a randomized phase II trial to select a pre-transplant induction regimen for future development. Patients had previously untreated stage III, IV, or bulky stage II MCL and received either 4 cycles of RH or 6 cycles of RB, followed by auto-HCT. Fifty-three of a planned 160 patients were accrued; an unacceptably high mobilization failure rate (29%) on the RH arm prompted premature study closure. The estimated 2-year progression-free survival (PFS) was 81% vs. 82% and overall survival (OS) was 87% vs. 88% for RB and RH, respectively. RH is not an ideal platform for future multi-centre transplant trials in MCL. RB achieved a 2-year PFS of 81% and a 78% MRD negative rate. Premature closure of the study limited the sample size and the precision of PFS estimates and MRD rates. However, RB can achieve a deep remission and could be a platform for future trials in MCL.

Noninvasive monitoring of diffuse large B-cell lymphoma by immunoglobulin high-throughput sequencing.

Recent studies have shown limited utility of routine surveillance imaging for diffuse large B-cell lymphoma (DLBCL) patients achieving remission. Detection of molecular disease by immunoglobulin high-throughput sequencing (Ig-HTS) from peripheral blood provides an alternate strategy for surveillance. We prospectively evaluated the utility of Ig-HTS within 311 blood and 105 tumor samples from 75 patients with DLBCL, comparing Ig-HTS from the cellular (circulating leukocytes) and acellular (plasma cell-free DNA) compartments of peripheral blood to clinical outcomes and (18)fluoro-deoxyglucose positron emission tomography combined with computed tomography (PET/CT; n = 173). Clonotypic immunoglobulin rearrangements were detected in 83% of patients with adequate tumor samples to enable subsequent monitoring in peripheral blood. Molecular disease measured from plasma, compared with circulating leukocytes, was more abundant and better correlated with radiographic disease burden. Before treatment, molecular disease was detected in the plasma of 82% of patients compared with 71% in circulating cells (P = .68). However, molecular disease was detected significantly more frequently in the plasma at time of relapse (100% vs 30%; P = .001). Detection of molecular disease in the plasma often preceded PET/CT detection of relapse in patients initially achieving remission. During surveillance time points before relapse, plasma Ig-HTS demonstrated improved specificity (100% vs 56%, P < .0001) and similar sensitivity (31% vs 55%, P = .4) compared with PET/CT. Given its high specificity, Ig-HTS from plasma has potential clinical utility for surveillance after complete remission.

Antigen receptor sequencing of paired bone marrow samples shows homogeneous distribution of acute lymphoblastic leukemia subclones.

In B-cell precursor acute lymphoblastic leukemia, the initial leukemic cells share the same antigen receptor gene rearrangements. However, due to ongoing rearrangement processes, leukemic cells with different gene rearrangement patterns can develop, resulting in subclone formation. We studied leukemic subclones and their distribution in the bone marrow and peripheral blood at diagnosis. Antigen receptor gene rearrangements (IGH, IGK, TRG, TRD, TRB) were analyzed by next-generation sequencing in seven paired bone marrow samples and five paired bone marrow-peripheral blood samples. Background-thresholds were defined, which enabled identification of leukemic gene rearrangements down to very low levels. Paired bone marrow analysis showed oligoclonality in all 7 patients and up to 34 leukemic clones per patient. Additional analysis of evolutionary-related IGH gene rearrangements revealed up to 171 leukemic clones per patient. Interestingly, overall 86% of all leukemic gene rearrangements, including small subclones, were present in both bone marrow samples (range per patient: 72-100%). Paired bone marrow-peripheral blood analysis showed that 83% of all leukemic gene rearrangements in bone marrow were also found in peripheral blood (range per patient: 81-100%). Remarkably, in the paired bone marrow samples and paired bone marrow-peripheral blood samples the vast majority of leukemic gene rearrangements had a similar frequency (<5-fold frequency difference) (96% and 96% of leukemic rearrangements, respectively). Together, these results indicate that B-cell precursor acute lymphoblastic leukemia is generally highly oligoclonal. Nevertheless, the vast majority of leukemic clones, even the minor antigen receptor-defined subclones, are homogeneously distributed throughout the bone marrow and peripheral blood compartment.

Quantification of Acute Lymphoblastic Leukemia Clonotypes in Leukapheresed Peripheral Blood Progenitor Cells Predicts Relapse Risk after Autologous Hematopoietic Stem Cell Transplantation.

Since the incorporation of tyrosine kinase inhibitors into the treatment of Philadelphia chromosome-positive (Ph+) acute lymphoblastic leukemia (ALL), the notion that all patients with "high-risk" ALL uniformly require allogeneic (allo) hematopoietic cell transplantation (HCT) has received increasing scrutiny. Although multiple studies have shown superiority of alloHCT over autologous (auto) hematopoietic cell transplantation for high-risk patients, these findings may be explained, in part, by contamination of the peripheral blood progenitor cell (PBPC) leukapheresis product by residual leukemic cells in patients undergoing autoHCT. We retrospectively evaluated minimal residual disease (MRD) using next-generation sequencing (NGS) in the PBPC leukapheresis product of 32 ALL patients who underwent autoHCT. Twenty-eight patients (88%) had diagnostic samples with quantifiable immunoreceptor rearrangements to follow for MRD. Twelve (38%) patients had Ph+ B-ALL, 12 (38%) had Philadelphia chromosome-negative (Ph-) B-ALL, and 4 (14%) had T cell ALL. With a median follow-up of 41 months (range, 3 to 217), median relapse-free survival (RFS) and overall survival for the entire cohort were 3.2 and 4.2 years, respectively; at 5 years after transplantation, 42% of patients remain alive and relapse free. Using MRD detection at a threshold of ≥ 1 × 10(-6), median RFS for patients with detectable MRD was 6.5 months and was not reached for patients without detectable disease (P = .0005). In multivariate analysis, the only factor significantly associated with relapse was the presence of MRD ≥1 × 10(-6) (odds ratio, 23.8; confidence interval, 1.8 to 312.9; P = .0158). Our findings suggest that NGS for MRD detection can predict long-term RFS in patients undergoing autoHCT for high-risk ALL.

A phase 2 study of Rituximab-Bendamustine and Rituximab-Cytarabine for transplant-eligible patients with mantle cell lymphoma.

Chemoimmunotherapy followed by autologous stem cell transplantation (ASCT) is a standard therapy for transplant-eligible patients with newly diagnosed mantle cell lymphoma (MCL). The achievement of complete remission (CR) and minimal residual disease (MRD) negativity are associated with better outcomes. We tested an induction regimen of rituximab/bendamustine followed by rituximab/high-dose cytarabine (RB/RC). This phase 2 study (NCT01661881) enrolled 23 transplant-eligible patients aged 42-69, of whom 70% were MCL international prognostic index low-risk. Patients received three cycles of RB followed by three cycles of RC. The primary endpoint of the trial was the rate of CR after six cycles of therapy, with a rate of 75% considered promising. 96% of patients achieved a CR/unconfirmed CR after treatment, meeting the primary objective. One patient progressed on study, one declined ASCT in CR, and the remaining 21 underwent successful stem cell collection and ASCT. After a median follow-up of 13 months, the progression-free survival rate was 96%. Among 15 MRD-evaluable patients who completed treatment, 93% achieved MRD negativity after RB/RC. In conclusion, RB/RC achieves very high CR and MRD negativity rates in transplant-eligible patients, with a favourable safety profile. RB/RC warrants further comparative studies, and may become a useful alternative to RCHOP (rituximab, cyclophosphamide, doxorubicin, vincristine, prednisone)-based induction regimens in this patient population.

Next-generation sequencing-based detection of circulating tumour DNA After allogeneic stem cell transplantation for lymphoma.

Next-generation sequencing (NGS)-based circulating tumour DNA (ctDNA) detection is a promising monitoring tool for lymphoid malignancies. We evaluated whether the presence of ctDNA was associated with outcome after allogeneic haematopoietic stem cell transplantation (HSCT) in lymphoma patients. We studied 88 patients drawn from a phase 3 clinical trial of reduced-intensity conditioning HSCT in lymphoma. Conventional restaging and collection of peripheral blood samples occurred at pre-specified time points before and after HSCT and were assayed for ctDNA by sequencing of the immunoglobulin or T-cell receptor genes. Tumour clonotypes were identified in 87% of patients with adequate tumour samples. Sixteen of 19 (84%) patients with disease progression after HSCT had detectable ctDNA prior to progression at a median of 3·7 months prior to relapse/progression. Patients with detectable ctDNA 3 months after HSCT had inferior progression-free survival (PFS) (2-year PFS 58% vs. 84% in ctDNA-negative patients, P = 0·033). In multivariate models, detectable ctDNA was associated with increased risk of progression/death (Hazard ratio 3·9, P = 0·003) and increased risk of relapse/progression (Hazard ratio 10·8, P = 0·0006). Detectable ctDNA is associated with an increased risk of relapse/progression, but further validation studies are necessary to confirm these findings and determine the clinical utility of NGS-based minimal residual disease monitoring in lymphoma patients after HSCT.

Prognostic value of deep sequencing method for minimal residual disease detection in multiple myeloma.

We assessed the prognostic value of minimal residual disease (MRD) detection in multiple myeloma (MM) patients using a sequencing-based platform in bone marrow samples from 133 MM patients in at least very good partial response (VGPR) after front-line therapy. Deep sequencing was carried out in patients in whom a high-frequency myeloma clone was identified and MRD was assessed using the IGH-VDJH, IGH-DJH, and IGK assays. The results were contrasted with those of multiparametric flow cytometry (MFC) and allele-specific oligonucleotide polymerase chain reaction (ASO-PCR). The applicability of deep sequencing was 91%. Concordance between sequencing and MFC and ASO-PCR was 83% and 85%, respectively. Patients who were MRD(-) by sequencing had a significantly longer time to tumor progression (TTP) (median 80 vs 31 months; P < .0001) and overall survival (median not reached vs 81 months; P = .02), compared with patients who were MRD(+). When stratifying patients by different levels of MRD, the respective TTP medians were: MRD ≥10(-3) 27 months, MRD 10(-3) to 10(-5) 48 months, and MRD <10(-5) 80 months (P = .003 to .0001). Ninety-two percent of VGPR patients were MRD(+). In complete response patients, the TTP remained significantly longer for MRD(-) compared with MRD(+) patients (131 vs 35 months; P = .0009).

Diverse somatic mutation patterns and pathway alterations in human cancers.

The systematic characterization of somatic mutations in cancer genomes is essential for understanding the disease and for developing targeted therapeutics. Here we report the identification of 2,576 somatic mutations across approximately 1,800 megabases of DNA representing 1,507 coding genes from 441 tumours comprising breast, lung, ovarian and prostate cancer types and subtypes. We found that mutation rates and the sets of mutated genes varied substantially across tumour types and subtypes. Statistical analysis identified 77 significantly mutated genes including protein kinases, G-protein-coupled receptors such as GRM8, BAI3, AGTRL1 (also called APLNR) and LPHN3, and other druggable targets. Integrated analysis of somatic mutations and copy number alterations identified another 35 significantly altered genes including GNAS, indicating an expanded role for galpha subunits in multiple cancer types. Furthermore, our experimental analyses demonstrate the functional roles of mutant GNAO1 (a Galpha subunit) and mutant MAP2K4 (a member of the JNK signalling pathway) in oncogenesis. Our study provides an overview of the mutational spectra across major human cancers and identifies several potential therapeutic targets.

Circulating tumour DNA and CT monitoring in patients with untreated diffuse large B-cell lymphoma: a correlative biomarker study.

BACKGROUND: Diffuse large-B-cell lymphoma is curable, but when treatment fails, outcome is poor. Although imaging can help to identify patients at risk of treatment failure, they are often imprecise, and radiation exposure is a potential health risk. We aimed to assess whether circulating tumour DNA encoding the clonal immunoglobulin gene sequence could be detected in the serum of patients with diffuse large-B-cell lymphoma and used to predict clinical disease recurrence after frontline treatment. METHODS: We used next-generation DNA sequencing to retrospectively analyse cell-free circulating tumour DNA in patients assigned to one of three treatment protocols between May 8, 1993, and June 6, 2013. Eligible patients had diffuse large-B-cell lymphoma, no evidence of indolent lymphoma, and were previously untreated. We obtained serial serum samples and concurrent CT scans at specified times during most treatment cycles and up to 5 years of follow-up. VDJ gene segments of the rearranged immunoglobulin receptor genes were amplified and sequenced from pretreatment specimens and serum circulating tumour DNA encoding the VDJ rearrangements was quantitated. FINDINGS: Tumour clonotypes were identified in pretreatment specimens from 126 patients who were followed up for a median of 11 years (IQR 6·8-14·2). Interim monitoring of circulating tumour DNA at the end of two treatment cycles in 108 patients showed a 5-year time to progression of 41·7% (95% CI 22·2-60·1) in patients with detectable circulating tumour DNA and 80·2% (69·6-87·3) in those without detectable circulating tumour DNA (p<0·0001). Detectable interim circulating tumour DNA had a positive predictive value of 62·5% (95% CI 40·6-81·2) and a negative predictive value of 79·8% (69·6-87·8). Surveillance monitoring of circulating tumour DNA was done in 107 patients who achieved complete remission. A Cox proportional hazards model showed that the hazard ratio for clinical disease progression was 228 (95% CI 51-1022) for patients who developed detectable circulating tumour DNA during surveillance compared with patients with undetectable circulating tumour DNA (p<0·0001). Surveillance circulating tumour DNA had a positive predictive value of 88·2% (95% CI 63·6-98·5) and a negative predictive value of 97·8% (92·2-99·7) and identified risk of recurrence at a median of 3·5 months (range 0-200) before evidence of clinical disease. INTERPRETATION: Surveillance circulating tumour DNA identifies patients at risk of recurrence before clinical evidence of disease in most patients and results in a reduced disease burden at relapse. Interim circulating tumour DNA is a promising biomarker to identify patients at high risk of treatment failure. FUNDING: National Cancer Institute and Adaptive Biotechnologies.

Detection of classical Hodgkin lymphoma specific sequence in peripheral blood using a next-generation sequencing approach.

We applied a highly sensitive next-generation sequencing method to identify lymphoma-specific immunoglobulin gene segments in classical Hodgkin lymphoma (CHL) at initial diagnosis or recurrence, and assessed the ability of detecting such lymphoma-specific sequences in peripheral blood (PB). Seventeen CHL cases were tested and lymphoma-specific sequences were identified in 12 of the primary tumour biopsies. In 11 of these patients whose paired PB samples were available, tumour-specific clonotypes were detected in PB in eight patients. This data demonstrates the feasibility of detecting circulating tumour-specific sequences, creating an unprecedented opportunity to optimize the future treatment and monitoring strategies for patients with CHL.

Immunoglobulin and T cell receptor gene high-throughput sequencing quantifies minimal residual disease in acute lymphoblastic leukemia and predicts post-transplantation relapse and survival.

Minimal residual disease (MRD) quantification is an important predictor of outcome after treatment for acute lymphoblastic leukemia (ALL). Bone marrow ALL burden ≥ 10(-4) after induction predicts subsequent relapse. Likewise, MRD ≥ 10(-4) in bone marrow before initiation of conditioning for allogeneic (allo) hematopoietic cell transplantation (HCT) predicts transplantation failure. Current methods for MRD quantification in ALL are not sufficiently sensitive for use with peripheral blood specimens and have not been broadly implemented in the management of adults with ALL. Consensus-primed immunoglobulin (Ig), T cell receptor (TCR) amplification and high-throughput sequencing (HTS) permit use of a standardized algorithm for all patients and can detect leukemia at 10(-6) or lower. We applied the LymphoSIGHT HTS platform (Sequenta Inc., South San Francisco, CA) to quantification of MRD in 237 samples from 29 adult B cell ALL patients before and after allo-HCT. Using primers for the IGH-VDJ, IGH-DJ, IGK, TCRB, TCRD, and TCRG loci, MRD could be quantified in 93% of patients. Leukemia-associated clonotypes at these loci were identified in 52%, 28%, 10%, 35%, 28%, and 41% of patients, respectively. MRD ≥ 10(-4) before HCT conditioning predicted post-HCT relapse (hazard ratio [HR], 7.7; 95% confidence interval [CI], 2.0 to 30; P = .003). In post-HCT blood samples, MRD ≥10(-6) had 100% positive predictive value for relapse with median lead time of 89 days (HR, 14; 95% CI, 4.7 to 44, P < .0001). The use of HTS-based MRD quantification in adults with ALL offers a standardized approach with sufficient sensitivity to quantify leukemia MRD in peripheral blood. Use of this approach may identify a window for clinical intervention before overt relapse.

Deep-sequencing approach for minimal residual disease detection in acute lymphoblastic leukemia.

The persistence of minimal residual disease (MRD) during therapy is the strongest adverse prognostic factor in acute lymphoblastic leukemia (ALL). We developed a high-throughput sequencing method that universally amplifies antigen-receptor gene segments and identifies all clonal gene rearrangements (ie, leukemia-specific sequences) at diagnosis, allowing monitoring of disease progression and clonal evolution during therapy. In the present study, the assay specifically detected 1 leukemic cell among greater than 1 million leukocytes in spike-in experiments. We compared this method with the gold-standard MRD assays multiparameter flow cytometry and allele-specific oligonucleotide polymerase chain reaction (ASO-PCR) using diagnostic and follow-up samples from 106 patients with ALL. Sequencing detected MRD in all 28 samples shown to be positive by flow cytometry and in 35 of the 36 shown to be positive by ASO-PCR and revealed MRD in 10 and 3 additional samples that were negative by flow cytometry and ASO-PCR, respectively. We conclude that this new method allows monitoring of treatment response in ALL and other lymphoid malignancies with great sensitivity and precision. The www.clinicaltrials.gov identifier number for the Total XV study is NCT00137111.

Massive evolution of the immunoglobulin heavy chain locus in children with B precursor acute lymphoblastic leukemia.

The ability to distinguish clonal B-cell populations based on the sequence of their rearranged immunoglobulin heavy chain (IgH) locus is an important tool for diagnosing B-cell neoplasms and monitoring treatment response. Leukemic precursor B cells may continue to undergo recombination of the IgH gene after malignant transformation; however, the magnitude of evolution at the IgH locus is currently unknown. We used next-generation sequencing to characterize the repertoire of IgH sequences in diagnostic samples of 51 children with B precursor acute lymphoblastic leukemia (B-ALL). We identified clonal IgH rearrangements in 43 of 51 (84%) cases and found that the number of evolved IgH sequences per patient ranged dramatically from 0 to 4024. We demonstrate that the evolved IgH sequences are not the result of amplification artifacts and are unique to leukemic precursor B cells. In addition, the evolution often follows an allelic exclusion pattern, where only 1 of 2 rearranged IgH loci exhibit ongoing recombination. Thus, precursor B-cell leukemias maintain evolution at the IgH locus at levels that were previously underappreciated. This finding sheds light on the mechanisms associated with leukemic clonal evolution and may fundamentally change approaches for monitoring minimal residual disease burden.

Human transcriptome array for high-throughput clinical studies.

A 6.9 million-feature oligonucleotide array of the human transcriptome [Glue Grant human transcriptome (GG-H array)] has been developed for high-throughput and cost-effective analyses in clinical studies. This array allows comprehensive examination of gene expression and genome-wide identification of alternative splicing as well as detection of coding SNPs and noncoding transcripts. The performance of the array was examined and compared with mRNA sequencing (RNA-Seq) results over multiple independent replicates of liver and muscle samples. Compared with RNA-Seq of 46 million uniquely mappable reads per replicate, the GG-H array is highly reproducible in estimating gene and exon abundance. Although both platforms detect similar expression changes at the gene level, the GG-H array is more sensitive at the exon level. Deeper sequencing is required to adequately cover low-abundance transcripts. The array has been implemented in a multicenter clinical program and has generated high-quality, reproducible data. Considering the clinical trial requirements of cost, sample availability, and throughput, the GG-H array has a wide range of applications. An emerging approach for large-scale clinical genomic studies is to first use RNA-Seq to the sufficient depth for the discovery of transcriptome elements relevant to the disease process followed by high-throughput and reliable screening of these elements on thousands of patient samples using custom-designed arrays.

Population structure, differential bias and genomic control in a large-scale, case-control association study.

The main problems in drawing causal inferences from epidemiological case-control studies are confounding by unmeasured extraneous factors, selection bias and differential misclassification of exposure. In genetics the first of these, in the form of population structure, has dominated recent debate. Population structure explained part of the significant +11.2% inflation of test statistics we observed in an analysis of 6,322 nonsynonymous SNPs in 816 cases of type 1 diabetes and 877 population-based controls from Great Britain. The remainder of the inflation resulted from differential bias in genotype scoring between case and control DNA samples, which originated from two laboratories, causing false-positive associations. To avoid excluding SNPs and losing valuable information, we extended the genomic control method by applying a variable downweighting to each SNP.

Ultra-sensitive molecular residual disease detection through whole genome sequencing with single-read error correction.

While whole genome sequencing (WGS) of cell-free DNA (cfDNA) holds enormous promise for detection of molecular residual disease (MRD), its performance is limited by WGS error rate. Here we introduce AccuScan, an efficient cfDNA WGS technology that enables genome-wide error correction at single read-level, achieving an error rate of 4.2 × 10-7, which is about two orders of magnitude lower than a read-centric de-noising method. The application of AccuScan to MRD demonstrated analytical sensitivity down to 10-6 circulating variant allele frequency at 99% sample-level specificity. AccuScan showed 90% landmark sensitivity (within 6 weeks after surgery) and 100% specificity for predicting relapse in colorectal cancer. It also showed 67% sensitivity and 100% specificity in esophageal cancer using samples collected within one week after surgery. When AccuScan was applied to monitor immunotherapy in melanoma patients, the circulating tumor DNA (ctDNA) levels and dynamic profiles were consistent with clinical outcomes. Overall, AccuScan provides a highly accurate WGS solution for MRD detection, empowering ctDNA detection at parts per million range without requiring high sample input or personalized reagents.

Ultra-sensitive molecular residual disease detection through whole genome sequencing with single-read error correction.

While whole genome sequencing (WGS) of cell-free DNA (cfDNA) holds enormous promise for molecular residual disease (MRD) detection, its performance is limited by WGS error rate. Here we introduce AccuScan, an efficient cfDNA WGS technology that enables genome-wide error correction at single read level, achieving an error rate of 4.2×10 -7 , which is about two orders of magnitude lower than a read-centric de-noising method. When applied to MRD detection, AccuScan demonstrated analytical sensitivity down to 10 -6 circulating tumor allele fraction at 99% sample level specificity. In colorectal cancer, AccuScan showed 90% landmark sensitivity for predicting relapse. It also showed robust MRD performance with esophageal cancer using samples collected as early as 1 week after surgery, and predictive value for immunotherapy monitoring with melanoma patients. Overall, AccuScan provides a highly accurate WGS solution for MRD, empowering circulating tumor DNA detection at parts per million range without high sample input nor personalized reagents. One Sentence Summary: AccuScan showed remarkable ultra-low limit of detection with a short turnaround time, low sample requirement and a simple workflow for MRD detection.