Detection of Copy-Number Variation in Circulating Cell-Free DNA in Patients With Uveal Melanoma.

PURPOSE: Somatic chromosomal alterations, particularly monosomy 3 and 8q gains, have been associated with metastatic risk in uveal melanoma (UM). Whole genome-scale evaluation of detectable alterations in cell-free DNA (cfDNA) in UM could provide valuable prognostic information. Our pilot study evaluates the correlation between genomic information using ultra-low-pass whole-genome sequencing (ULP-WGS) of cfDNA in UM and associated clinical outcomes. MATERIALS AND METHODS: ULP-WGS of cfDNA was performed on 29 plasma samples from 16 patients, 14 metastatic UM (mUM) and two non-metastatic, including pre- and post-treatment mUM samples from 10 patients treated with immunotherapy and one with liver-directed therapy. We estimated tumor fraction (TFx) and detected copy-number alterations (CNAs) using ichorCNA. Presence of 8q amplification was further analyzed using the likelihood ratio test (LRT). RESULTS: Eleven patients with mUM (17 samples) of 14 had detectable circulating tumor DNA (ctDNA). 8q gain was detected in all 17, whereas monosomy 3 was detectable in 10 of 17 samples. TFx generally correlated with disease status, showing an increase at the time of disease progression (PD). 8q gain detection sensitivity appeared greater with the LRT than with ichorCNA at lower TFxs. The only patient with mUM with partial response on treatment had a high pretreatment TFx and undetectable on-treatment ctDNA, correlating with her profound response and durable survival. CONCLUSION: ctDNA can be detected in mUM using ULP-WGS, and the TFx correlates with DS. 8q gain was consistently detectable in mUM, in line with previous studies indicating 8q gains early in primary UM and higher amplification with PD. Our work suggests that detection of CNAs by ULP-WGS, particularly focusing on 8q gain, could be a valuable blood biomarker to monitor PD in UM.

Liquid biopsy detection of genomic alterations in pediatric brain tumors from cell-free DNA in peripheral blood, CSF, and urine.

BACKGROUND: The ability to identify genetic alterations in cancers is essential for precision medicine; however, surgical approaches to obtain brain tumor tissue are invasive. Profiling circulating tumor DNA (ctDNA) in liquid biopsies has emerged as a promising approach to avoid invasive procedures. Here, we systematically evaluated the feasibility of profiling pediatric brain tumors using ctDNA obtained from plasma, cerebrospinal fluid (CSF), and urine. METHODS: We prospectively collected 564 specimens (257 blood, 240 urine, and 67 CSF samples) from 258 patients across all histopathologies. We performed ultra-low-pass whole-genome sequencing (ULP-WGS) to assess copy number variations and estimate tumor fraction and developed a pediatric CNS tumor hybrid capture panel for deep sequencing of specific mutations and fusions. RESULTS: ULP-WGS detected copy number alterations in 9/46 (20%) CSF, 3/230 (1.3%) plasma, and 0/153 urine samples. Sequencing detected alterations in 3/10 (30%) CSF, 2/74 (2.7%) plasma, and 0/2 urine samples. The only positive results were in high-grade tumors. However, most samples had insufficient somatic mutations (median 1, range 0-39) discoverable by the sequencing panel to provide sufficient power to detect tumor fractions of greater than 0.1%. CONCLUSIONS: Children with brain tumors harbor very low levels of ctDNA in blood, CSF, and urine, with CSF having the most DNA detectable. Molecular profiling is feasible in a small subset of high-grade tumors. The level of clonal aberrations per genome is low in most of the tumors, posing a challenge for detection using whole-genome or even targeted sequencing methods. Substantial challenges therefore remain to genetically characterize pediatric brain tumors from liquid biopsies.

Donor Clonal Hematopoiesis and Recipient Outcomes After Transplantation.

PURPOSE: Clonal hematopoiesis (CH) can be transmitted from a donor to a recipient during allogeneic hematopoietic cell transplantation. Exclusion of candidate donors with CH is controversial since its impact on recipient outcomes and graft alloimmune function is uncertain. PATIENTS AND METHODS: We performed targeted error-corrected sequencing on samples from 1,727 donors age 40 years or older and assessed the effect of donor CH on recipient clinical outcomes. We measured long-term engraftment of 102 donor clones and cytokine levels in 256 recipients at 3 and 12 months after transplant. RESULTS: CH was present in 22.5% of donors, with DNMT3A (14.6%) and TET2 (5.2%) mutations being most common; 85% of donor clones showed long-term engraftment in recipients after transplantation, including clones with a variant allele fraction < 0.01. DNMT3A-CH with a variant allele fraction ≥ 0.01, but not smaller clones, was associated with improved recipient overall (hazard ratio [HR], 0.79; P = .042) and progression-free survival (HR, 0.72; P = .003) after adjustment for significant clinical variables. In patients who received calcineurin-based graft-versus-host disease prophylaxis, donor DNMT3A-CH was associated with reduced relapse (subdistribution HR, 0.59; P = .014), increased chronic graft-versus-host disease (subdistribution HR, 1.36; P = .042), and higher interleukin-12p70 levels in recipients. No recipient of sole DNMT3A or TET2-CH developed donor cell leukemia (DCL). In seven of eight cases, DCL evolved from donor CH with rare TP53 or splicing factor mutations or from donors carrying germline DDX41 mutations. CONCLUSION: Donor CH is closely associated with clinical outcomes in transplant recipients, with differential impact on graft alloimmune function and potential for leukemic transformation related to mutated gene and somatic clonal abundance. Donor DNMT3A-CH is associated with improved recipient survival because of reduced relapse risk and with an augmented network of inflammatory cytokines in recipients. Risk of DCL in allogeneic hematopoietic cell transplantation is driven by somatic myelodysplastic syndrome-associated mutations or germline predisposition in donors.

Sensitive Detection of Minimal Residual Disease in Patients Treated for Early-Stage Breast Cancer.

PURPOSE: Existing cell-free DNA (cfDNA) methods lack the sensitivity needed for detecting minimal residual disease (MRD) following therapy. We developed a test for tracking hundreds of patient-specific mutations to detect MRD with a 1,000-fold lower error rate than conventional sequencing. EXPERIMENTAL DESIGN: We compared the sensitivity of our approach to digital droplet PCR (ddPCR) in a dilution series, then retrospectively identified two cohorts of patients who had undergone prospective plasma sampling and clinical data collection: 16 patients with ER+/HER2- metastatic breast cancer (MBC) sampled within 6 months following metastatic diagnosis and 142 patients with stage 0 to III breast cancer who received curative-intent treatment with most sampled at surgery and 1 year postoperative. We performed whole-exome sequencing of tumors and designed individualized MRD tests, which we applied to serial cfDNA samples. RESULTS: Our approach was 100-fold more sensitive than ddPCR when tracking 488 mutations, but most patients had fewer identifiable tumor mutations to track in cfDNA (median = 57; range = 2-346). Clinical sensitivity was 81% (n = 13/16) in newly diagnosed MBC, 23% (n = 7/30) at postoperative and 19% (n = 6/32) at 1 year in early-stage disease, and highest in patients with the most tumor mutations available to track. MRD detection at 1 year was strongly associated with distant recurrence [HR = 20.8; 95% confidence interval, 7.3-58.9]. Median lead time from first positive sample to recurrence was 18.9 months (range = 3.4-39.2 months). CONCLUSIONS: Tracking large numbers of individualized tumor mutations in cfDNA can improve MRD detection, but its sensitivity is driven by the number of tumor mutations available to track.

Best practices for benchmarking germline small-variant calls in human genomes.

Standardized benchmarking approaches are required to assess the accuracy of variants called from sequence data. Although variant-calling tools and the metrics used to assess their performance continue to improve, important challenges remain. Here, as part of the Global Alliance for Genomics and Health (GA4GH), we present a benchmarking framework for variant calling. We provide guidance on how to match variant calls with different representations, define standard performance metrics, and stratify performance by variant type and genome context. We describe limitations of high-confidence calls and regions that can be used as truth sets (for example, single-nucleotide variant concordance of two methods is 99.7% inside versus 76.5% outside high-confidence regions). Our web-based app enables comparison of variant calls against truth sets to obtain a standardized performance report. Our approach has been piloted in the PrecisionFDA variant-calling challenges to identify the best-in-class variant-calling methods within high-confidence regions. Finally, we recommend a set of best practices for using our tools and evaluating the results.

Author Correction: Best practices for benchmarking germline small-variant calls in human genomes.

In the version of this article initially published online, two pairs of headings were switched with each other in Table 4: "Recall (PCR free)" was switched with "Recall (with PCR)," and "Precision (PCR free)" was switched with "Precision (with PCR)." The error has been corrected in the print, PDF and HTML versions of this article.

A synthetic-diploid benchmark for accurate variant-calling evaluation.

Existing benchmark datasets for use in evaluating variant-calling accuracy are constructed from a consensus of known short-variant callers, and they are thus biased toward easy regions that are accessible by these algorithms. We derived a new benchmark dataset from the de novo PacBio assemblies of two fully homozygous human cell lines, which provides a relatively more accurate and less biased estimate of small-variant-calling error rates in a realistic context.

Characterization and remediation of sample index swaps by non-redundant dual indexing on massively parallel sequencing platforms.

BACKGROUND: Here we present an in-depth characterization of the mechanism of sequencer-induced sample contamination due to the phenomenon of index swapping that impacts Illumina sequencers employing patterned flow cells with Exclusion Amplification (ExAmp) chemistry (HiSeqX, HiSeq4000, and NovaSeq). We also present a remediation method that minimizes the impact of such swaps. RESULTS: Leveraging data collected over a two-year period, we demonstrate the widespread prevalence of index swapping in patterned flow cell data. We calculate mean swap rates across multiple sample preparation methods and sequencer models, demonstrating that different library methods can have vastly different swapping rates and that even non-ExAmp chemistry instruments display trace levels of index swapping. We provide methods for eliminating sample data cross contamination by utilizing non-redundant dual indexing for complete filtering of index swapped reads, and share the sequences for 96 non-combinatorial dual indexes we have validated across various library preparation methods and sequencer models. Finally, using computational methods we provide a greater insight into the mechanism of index swapping. CONCLUSIONS: Index swapping in pooled libraries is a prevalent phenomenon that we observe at a rate of 0.2 to 6% in all sequencing runs on HiSeqX, HiSeq 4000/3000, and NovaSeq. Utilizing non-redundant dual indexing allows for the removal (flagging/filtering) of these swapped reads and eliminates swapping induced sample contamination, which is critical for sensitive applications such as RNA-seq, single cell, blood biopsy using circulating tumor DNA, or clinical sequencing.

Corrigendum: Landscape of X chromosome inactivation across human tissues.

This corrects the article DOI: 10.1038/nature24265.

Landscape of X chromosome inactivation across human tissues.

X chromosome inactivation (XCI) silences transcription from one of the two X chromosomes in female mammalian cells to balance expression dosage between XX females and XY males. XCI is, however, incomplete in humans: up to one-third of X-chromosomal genes are expressed from both the active and inactive X chromosomes (Xa and Xi, respectively) in female cells, with the degree of 'escape' from inactivation varying between genes and individuals. The extent to which XCI is shared between cells and tissues remains poorly characterized, as does the degree to which incomplete XCI manifests as detectable sex differences in gene expression and phenotypic traits. Here we describe a systematic survey of XCI, integrating over 5,500 transcriptomes from 449 individuals spanning 29 tissues from GTEx (v6p release) and 940 single-cell transcriptomes, combined with genomic sequence data. We show that XCI at 683 X-chromosomal genes is generally uniform across human tissues, but identify examples of heterogeneity between tissues, individuals and cells. We show that incomplete XCI affects at least 23% of X-chromosomal genes, identify seven genes that escape XCI with support from multiple lines of evidence and demonstrate that escape from XCI results in sex biases in gene expression, establishing incomplete XCI as a mechanism that is likely to introduce phenotypic diversity. Overall, this updated catalogue of XCI across human tissues helps to increase our understanding of the extent and impact of the incompleteness in the maintenance of XCI.

Solvent Role in the Formation of Electric Double Layers with Surface Charge Regulation: A Bystander or a Key Participant?

The charge formation at interfaces involving electrolyte solutions is due to the chemical equilibrium between the surface reactive groups and the potential determining ions in the solution (i.e., charge regulation). In this Letter we report our findings that this equilibrium is strongly coupled to the precise molecular structure of the solution near the charged interface. The neutral solvent molecules dominate this structure due to their overwhelmingly large number. Treating the solvent as a structureless continuum leads to a fundamentally inadequate physical picture of charged interfaces. We show that a proper account of the solvent effect leads to an unexpected and complex system behavior that is affected by the molecular and ionic excluded volumes and van der Waals interactions.

Charge regulation at semiconductor-electrolyte interfaces.

The interface between a semiconductor material and an electrolyte solution has interesting and complex electrostatic properties. Its behavior will depend on the density of mobile charge carriers that are present in both phases as well as on the surface chemistry at the interface through local charge regulation. The latter is driven by chemical equilibria involving the immobile surface groups and the potential determining ions in the electrolyte solution. All these lead to an electrostatic potential distribution that propagate such that the electrolyte and the semiconductor are dependent on each other. Hence, any variation in the charge density in one phase will lead to a response in the other. This has significant implications on the physical properties of single semiconductor-electrolyte interfaces and on the electrostatic interactions between semiconductor particles suspended in electrolyte solutions. The present paper expands on our previous publication (Fleharty et al., 2014) and offers new results on the electrostatics of single semiconductor interfaces as well as on the interaction of charged semiconductor colloids suspended in electrolyte solution.

Manipulating semiconductor colloidal stability through doping.

The interface between a doped semiconductor material and electrolyte solution is of considerable fundamental interest, and is relevant to systems of practical importance. Both adjacent domains contain mobile charges, which respond to potential variations. This is exploited to design electronic and optoelectronic sensors, and other enabling semiconductor colloidal materials. We show that the charge mobility in both phases leads to a new type of interaction between semiconductor colloids suspended in aqueous electrolyte solutions. This interaction is due to the electrostatic response of the semiconductor interior to disturbances in the external field upon the approach of two particles. The electrostatic repulsion between two charged colloids is reduced from the one governed by the charged groups present at the particles surfaces. This type of interaction is unique to semiconductor particles and may have a substantial effect on the suspension dynamics and stability.

The effect of surface charge regulation on conductivity in fluidic nanochannels.

The precise electrostatic potential distribution is very important for the electrokinetic transport in fluidic channels. This is especially valid for small nanochannels where the electric double layers formed at the walls are comparable to the channel width. It can be expected that due to the large surface to volume ratio in such systems, they will exhibit properties that are not detectable in larger channels, capillaries and pores. We present a detailed numerical analysis of the current transport in fluidic nanochannels. It is based on solving the Poisson-Boltzmann equation with charge regulation boundary conditions that account for the surface-aqueous solution chemical equilibria. The focus is on studying the effect of the pH on the current transport. The pH is varied by adding either HCl or KOH. The analysis predicts non-monotonous and sometimes counterintuitive dependence of the conductivity on the pH. The channel conductivity exhibits practically no change over a range of pH values due to a buffering exerted by the chemical groups at the walls. An unexpected drop of the conductivity is observed around the wall isoelectric point and also in the vicinity of pH=7 even though the concentration of ions in the channel increases. These observations are explained in the framework of charge regulation theory.