Efficient and unique cobarcoding of second-generation sequencing reads from long DNA molecules enabling cost-effective and accurate sequencing, haplotyping, and de novo assembly.

Here, we describe single-tube long fragment read (stLFR), a technology that enables sequencing of data from long DNA molecules using economical second-generation sequencing technology. It is based on adding the same barcode sequence to subfragments of the original long DNA molecule (DNA cobarcoding). To achieve this efficiently, stLFR uses the surface of microbeads to create millions of miniaturized barcoding reactions in a single tube. Using a combinatorial process, up to 3.6 billion unique barcode sequences were generated on beads, enabling practically nonredundant cobarcoding with 50 million barcodes per sample. Using stLFR, we demonstrate efficient unique cobarcoding of more than 8 million 20- to 300-kb genomic DNA fragments. Analysis of the human genome NA12878 with stLFR demonstrated high-quality variant calling and phase block lengths up to N50 34 Mb. We also demonstrate detection of complex structural variants and complete diploid de novo assembly of NA12878. These analyses were all performed using single stLFR libraries, and their construction did not significantly add to the time or cost of whole-genome sequencing (WGS) library preparation. stLFR represents an easily automatable solution that enables high-quality sequencing, phasing, SV detection, scaffolding, cost-effective diploid de novo genome assembly, and other long DNA sequencing applications.

Accurate whole-genome sequencing and haplotyping from 10 to 20 human cells.

Recent advances in whole-genome sequencing have brought the vision of personal genomics and genomic medicine closer to reality. However, current methods lack clinical accuracy and the ability to describe the context (haplotypes) in which genome variants co-occur in a cost-effective manner. Here we describe a low-cost DNA sequencing and haplotyping process, long fragment read (LFR) technology, which is similar to sequencing long single DNA molecules without cloning or separation of metaphase chromosomes. In this study, ten LFR libraries were made using only ∼100 picograms of human DNA per sample. Up to 97% of the heterozygous single nucleotide variants were assembled into long haplotype contigs. Removal of false positive single nucleotide variants not phased by multiple LFR haplotypes resulted in a final genome error rate of 1 in 10 megabases. Cost-effective and accurate genome sequencing and haplotyping from 10-20 human cells, as demonstrated here, will enable comprehensive genetic studies and diverse clinical applications.

Structure of the SWI2/SNF2 chromatin-remodeling domain of eukaryotic Rad54.

SWI2/SNF2 chromatin-remodeling proteins mediate the mobilization of nucleosomes and other DNA-associated proteins. SWI2/SNF2 proteins contain sequence motifs characteristic of SF2 helicases but do not have helicase activity. Instead, they couple ATP hydrolysis with the generation of superhelical torsion in DNA. The structure of the nucleosome-remodeling domain of zebrafish Rad54, a protein involved in Rad51-mediated homologous recombination, reveals that the core of the SWI2/SNF2 enzymes consist of two alpha/beta-lobes similar to SF2 helicases. The Rad54 helicase lobes contain insertions that form two helical domains, one within each lobe. These insertions contain SWI2/SNF2-specific sequence motifs likely to be central to SWI2/SNF2 function. A broad cleft formed by the two lobes and flanked by the helical insertions contains residues conserved in SWI2/SNF2 proteins and motifs implicated in DNA-binding by SF2 helicases. The Rad54 structure suggests that SWI2/SNF2 proteins use a mechanism analogous to helicases to translocate on dsDNA.

3' Branch ligation: a novel method to ligate non-complementary DNA to recessed or internal 3'OH ends in DNA or RNA.

Nucleic acid ligases are crucial enzymes that repair breaks in DNA or RNA during synthesis, repair and recombination. Various genomic tools have been developed using the diverse activities of DNA/RNA ligases. Herein, we demonstrate a non-conventional ability of T4 DNA ligase to insert 5' phosphorylated blunt-end double-stranded DNA to DNA breaks at 3'-recessive ends, gaps, or nicks to form a Y-shaped 3'-branch structure. Therefore, this base pairing-independent ligation is termed 3'-branch ligation (3'BL). In an extensive study of optimal ligation conditions, the presence of 10% PEG-8000 in the ligation buffer significantly increased ligation efficiency to more than 80%. Ligation efficiency was slightly varied between different donor and acceptor sequences. More interestingly, we discovered that T4 DNA ligase efficiently ligated DNA to the 3'-recessed end of RNA, not to that of DNA, in a DNA/RNA hybrid, suggesting a ternary complex formation preference of T4 DNA ligase. These novel properties of T4 DNA ligase can be utilized as a broad molecular technique in many important genomic applications, such as 3'-end labelling by adding a universal sequence; directional tagmentation for NGS library construction that achieve theoretical 100% template usage; and targeted RNA NGS libraries with mitigated structure-based bias and adapter dimer problems.

Development of coupling controlled polymerizations by adapter-ligation in mate-pair sequencing for detection of various genomic variants in one single assay.

The diversity of disease presentations warrants one single assay for detection and delineation of various genomic disorders. Herein, we describe a gel-free and biotin-capture-free mate-pair method through coupling Controlled Polymerizations by Adapter-Ligation (CP-AL). We first demonstrated the feasibility and ease-of-use in monitoring DNA nick translation and primer extension by limiting the nucleotide input. By coupling these two controlled polymerizations by a reported non-conventional adapter-ligation reaction 3' branch ligation, we evidenced that CP-AL significantly increased DNA circularization efficiency (by 4-fold) and was applicable for different sequencing methods but at a faction of current cost. Its advantages were further demonstrated by fully elimination of small-insert-contaminated (by 39.3-fold) with a ∼50% increment of physical coverage, and producing uniform genome/exome coverage and the lowest chimeric rate. It achieved single-nucleotide variants detection with sensitivity and specificity up to 97.3 and 99.7%, respectively, compared with data from small-insert libraries. In addition, this method can provide a comprehensive delineation of structural rearrangements, evidenced by a potential diagnosis in a patient with oligo-atheno-terato-spermia. Moreover, it enables accurate mutation identification by integration of genomic variants from different aberration types. Overall, it provides a potential single-integrated solution for detecting various genomic variants, facilitating a genetic diagnosis in human diseases.

cPAS-based sequencing on the BGISEQ-500 to explore small non-coding RNAs.

BACKGROUND: We present the first sequencing data using the combinatorial probe-anchor synthesis (cPAS)-based BGISEQ-500 sequencer. Applying cPAS, we investigated the repertoire of human small non-coding RNAs and compared it to other techniques. RESULTS: Starting with repeated measurements of different specimens including solid tissues (brain and heart) and blood, we generated a median of 30.1 million reads per sample. 24.1 million mapped to the human genome and 23.3 million to the miRBase. Among six technical replicates of brain samples, we observed a median correlation of 0.98. Comparing BGISEQ-500 to HiSeq, we calculated a correlation of 0.75. The comparability to microarrays was similar for both BGISEQ-500 and HiSeq with the first one showing a correlation of 0.58 and the latter one correlation of 0.6. As for a potential bias in the detected expression distribution in blood cells, 98.6% of HiSeq reads versus 93.1% of BGISEQ-500 reads match to the 10 miRNAs with highest read count. After using miRDeep2 and employing stringent selection criteria for predicting new miRNAs, we detected 74 high-likely candidates in the cPAS sequencing reads prevalent in solid tissues and 36 candidates prevalent in blood. CONCLUSIONS: While there is apparently no ideal platform for all challenges of miRNome analyses, cPAS shows high technical reproducibility and supplements the hitherto available platforms.

Rad54 protein possesses chromatin-remodeling activity stimulated by the Rad51-ssDNA nucleoprotein filament.

In Saccharomyces cerevisiae, the Rad54 protein participates in the recombinational repair of double-strand DNA breaks together with the Rad51, Rad52, Rad55 and Rad57 proteins. In vitro, Rad54 interacts with Rad51 and stimulates DNA strand exchange promoted by Rad51 protein. Rad54 is a SWI2/SNF2-related protein that possesses double-stranded DNA-dependent ATPase activity and changes DNA topology in an ATP hydrolysis-dependent manner. Here we show that Rad54 catalyzes bidirectional nucleosome redistribution by sliding nucleosomes along DNA. Nucleosome redistribution is greatly stimulated by the Rad51 nucleoprotein filament but does not require the presence of homologous single-stranded DNA within the filament. On the basis of these data, we propose that Rad54 facilitates chromatin remodeling and, perhaps more generally, protein clearing at the homology search step of genetic recombination.

A novel function of Rad54 protein. Stabilization of the Rad51 nucleoprotein filament.

Homologous recombination is important for the repair of double-stranded DNA breaks in all organisms. Rad51 and Rad54 proteins are two key components of the homologous recombination machinery in eukaryotes. In vitro, Rad51 protein assembles with single-stranded DNA to form the helical nucleoprotein filament that promotes DNA strand exchange, a basic step of homologous recombination. Rad54 protein interacts with this Rad51 nucleoprotein filament and stimulates its DNA pairing activity, suggesting that Rad54 protein is a component of the nucleoprotein complex involved in the DNA homology search. Here, using physical criteria, we demonstrate directly the formation of Rad54-Rad51-DNA nucleoprotein co-complexes that contain equimolar amounts of each protein. The binding of Rad54 protein significantly stabilizes the Rad51 nucleoprotein filament formed on either single-stranded DNA or double-stranded DNA. The Rad54-stabilized nucleoprotein filament is more competent in DNA strand exchange and acts over a broader range of solution conditions. Thus, the co-assembly of an interacting partner with the Rad51 nucleoprotein filament represents a novel means of stabilizing the biochemical entity central to homologous recombination, and reveals a new function of Rad54 protein.