Direct capture and sequencing reveal ultra-short single-stranded DNA in biofluids.

Cell-free DNA (cfDNA) has become the predominant analyte of liquid biopsy; however, recent studies suggest the presence of subnucleosomal-sized DNA fragments in circulation that are likely single-stranded. Here, we report a method called direct capture and sequencing (DCS) tailored to recover such fragments from biofluids by directly capturing them using short degenerate probes followed by single strand-based library preparation and next-generation sequencing. DCS revealed a new DNA population in biofluids, named ultrashort single-stranded DNA (ussDNA). Evaluation of the size distribution and abundance of ussDNA manifested generality of its presence in humans, animal species, and plants. In humans, red blood cells were found to contain abundant ussDNA; plasma-derived ussDNA exhibited modal size at 50 nt. This work reports the presence of an understudied DNA population in circulation, and yet more work is awaiting to study its generation mechanism, tissue of origin, disease implications, etc.

Low-level parental somatic mosaic SNVs in exomes from a large cohort of trios with diverse suspected Mendelian conditions.

PURPOSE: The goal of this study was to assess the scale of low-level parental mosaicism in exome sequencing (ES) databases. METHODS: We analyzed approximately 2000 family trio ES data sets from the Baylor-Hopkins Center for Mendelian Genomics (BHCMG) and Baylor Genetics (BG). Among apparent de novo single-nucleotide variants identified in the affected probands, we selected rare unique variants with variant allele fraction (VAF) between 30% and 70% in the probands and lower than 10% in one of the parents. RESULTS: Of 102 candidate mosaic variants validated using amplicon-based next-generation sequencing, droplet digital polymerase chain reaction, or blocker displacement amplification, 27 (26.4%) were confirmed to be low- (VAF between 1% and 10%) or very low (VAF <1%) level mosaic. Detection precision in parental samples with two or more alternate reads was 63.6% (BHCMG) and 43.6% (BG). In nine investigated individuals, we observed variability of mosaic ratios among blood, saliva, fibroblast, buccal, hair, and urine samples. CONCLUSION: Our computational pipeline enables robust discrimination between true and false positive candidate mosaic variants and efficient detection of low-level mosaicism in ES samples. We confirm that the presence of two or more alternate reads in the parental sample is a reliable predictor of low-level parental somatic mosaicism.

Nucleic Acid Quantitation with Log-Linear Response Hybridization Probe Sets.

Concentrations of different nucleic acid species in biological samples span many orders of magnitude. A real-time polymerase chain reaction maps the concentration of a target nucleic acid sequence log-linearly into cycle threshold to enable quantitation with a wide dynamic range but suffers from enzymatic biases. Here, we present a general design for constructing hybridization probe sets with highly log-linear response curves to enable accurate enzyme-free quantitation across large ranges (more than 6 logs) of target DNA concentrations. The sensitivity of each component probe is accurately adjusted via formulation stoichiometry to reduce the standard error of target quantitation down to 7%. As a proof of concept, we show multiplexed quantitation of three microRNA species in total RNA of the human brain and liver.

Encoding multiple digital DNA signals in a single analog channel.

For many analytic and biomedical applications, the presence of an analyte above or below a critical concentration is more informative for decision making than the actual concentration value. Straightforward analog-to-digital signal conversion does not take full advantage of the precision and dynamic range of modern sensors. Here, we present and experimentally demonstrate an analog-to-multiple-digital signal conversion, reporting digital signals that indicate whether the concentrations of specific DNA sequences exceed respective threshold values. These threshold values can be individually programmed for each target sequence. Experimentally, we showed representation of four DNA targets' information in a single fluorescence channel.

Highly Sensitive Blocker Displacement Amplification and Droplet Digital PCR Reveal Low-Level Parental FOXF1 Somatic Mosaicism in Families with Alveolar Capillary Dysplasia with Misalignment of Pulmonary Veins.

Detection of low-level somatic mosaicism [alternate allele fraction (AAF) ≤ 10%] in parents of affected individuals with the apparent de novo pathogenic variants enables more accurate estimate of recurrence risk. To date, only a few systematic analyses of low-level parental somatic mosaicism have been performed. Herein, highly sensitive blocker displacement amplification, droplet digital PCR, quantitative PCR, long-range PCR, and array comparative genomic hybridization were applied in families with alveolar capillary dysplasia with misalignment of pulmonary veins. We screened 18 unrelated families with the FOXF1 variant previously determined to be apparent de novo (n = 14), of unknown parental origin (n = 1), or inherited from a parent suspected to be somatic and/or germline mosaic (n = 3). We identified four (22%) families with FOXF1 parental somatic mosaic single-nucleotide variants (n = 3) and copy number variant deletion (n = 1) detected in parental blood samples and an AAF ranging between 0.03% and 19%. In one family, mosaic allele ratio in tissues originating from three germ layers ranged between <0.03% and 0.65%. Because the ratio of parental somatic mosaicism have significant implications for the recurrence risk, this study further implies the importance of a systematic screening of parental samples for low-level and very-low-level (AAF ≤ 1%) somatic mosaicism using methods that are more sensitive than those routinely applied in diagnostics.

A recurrent 8 bp frameshifting indel in FOXF1 defines a novel mutation hotspot associated with alveolar capillary dysplasia with misalignment of pulmonary veins.

Alveolar capillary dysplasia with misalignment of pulmonary veins (ACDMPV) is a rare lethal lung developmental disease. Affected infants manifest with severe respiratory distress and refractory pulmonary hypertension and uniformly die in the first month of life. Heterozygous point mutations or copy-number variant deletions involving FOXF1 and/or its upstream lung-specific enhancer on 16q24.1 have been identified in the vast majority of ACDMPV patients. We have previously described two unrelated families with a de novo pathogenic frameshift variant c.691_698del (p.Ala231Argfs*61) in the exon 1 of FOXF1. Here, we present a third unrelated ACDMPV family with the same de novo variant and propose that a direct tandem repeat of eight consecutive nucleotides GCGGCGGC within the ~4 kb CpG island in FOXF1 exon 1 is a novel mutation hotspot causative for ACDMPV.

Predicting DNA hybridization kinetics from sequence.

Hybridization is a key molecular process in biology and biotechnology, but so far there is no predictive model for accurately determining hybridization rate constants based on sequence information. Here, we report a weighted neighbour voting (WNV) prediction algorithm, in which the hybridization rate constant of an unknown sequence is predicted based on similarity reactions with known rate constants. To construct this algorithm we first performed 210 fluorescence kinetics experiments to observe the hybridization kinetics of 100 different DNA target and probe pairs (36 nt sub-sequences of the CYCS and VEGF genes) at temperatures ranging from 28 to 55 °C. Automated feature selection and weighting optimization resulted in a final six-feature WNV model, which can predict hybridization rate constants of new sequences to within a factor of 3 with ∼91% accuracy, based on leave-one-out cross-validation. Accurate prediction of hybridization kinetics allows the design of efficient probe sequences for genomics research.

Multiplexed enrichment of rare DNA variants via sequence-selective and temperature-robust amplification.

Rare DNA-sequence variants hold important clinical and biological information, but existing detection techniques are expensive, complex, allele-specific, or don't allow for significant multiplexing. Here, we report a temperature-robust polymerase-chain-reaction method, which we term blocker displacement amplification (BDA), that selectively amplifies all sequence variants, including single-nucleotide variants (SNVs), within a roughly 20-nucleotide window by 1,000-fold over wild-type sequences. This allows for easy detection and quantitation of hundreds of potential variants originally at ≤0.1% in allele frequency. BDA is compatible with inexpensive thermocycler instrumentation and employs a rationally designed competitive hybridization reaction to achieve comparable enrichment performance across annealing temperatures ranging from 56 °C to 64 °C. To show the sequence generality of BDA, we demonstrate enrichment of 156 SNVs and the reliable detection of single-digit copies. We also show that the BDA detection of rare driver mutations in cell-free DNA samples extracted from the blood plasma of lung-cancer patients is highly consistent with deep sequencing using molecular lineage tags, with a receiver operator characteristic accuracy of 95%.

Publisher Correction: Multiplexed enrichment of rare DNA variants via sequence-selective and temperature-robust amplification.

In the version of this Article originally published, owing to a technical error, the Life Sciences Reporting Summary was not included; this summary is now available.

Continuously tunable nucleic acid hybridization probes.

In silico-designed nucleic acid probes and primers often do not achieve favorable specificity and sensitivity tradeoffs on the first try, and iterative empirical sequence-based optimization is needed, particularly in multiplexed assays. We present a novel, on-the-fly method of tuning probe affinity and selectivity by adjusting the stoichiometry of auxiliary species, which allows for independent and decoupled adjustment of the hybridization yield for different probes in multiplexed assays. Using this method, we achieved near-continuous tuning of probe effective free energy. To demonstrate our approach, we enforced uniform capture efficiency of 31 DNA molecules (GC content, 0-100%), maximized the signal difference for 11 pairs of single-nucleotide variants and performed tunable hybrid capture of mRNA from total RNA. Using the Nanostring nCounter platform, we applied stoichiometric tuning to simultaneously adjust yields for a 24-plex assay, and we show multiplexed quantitation of RNA sequences and variants from formalin-fixed, paraffin-embedded samples.