Error removal in microchip-synthesized DNA using immobilized MutS

The development of economical de novo gene synthesis methods using microchip-synthesized oligonucleotides has been limited by their high error rates. In this study, a low-cost, effective and improved-throughput (up to 32 oligos per run) error-removal method using an immobilized cellulose column containing the mismatch binding protein MutS was produced to generate high-quality DNA from oligos, particularly microchip-synthesized oligonucleotides. Error-containing DNA in the initial material was specifically retained on the MutS-immobilized cellulose column (MICC), and error-depleted DNA in the eluate was collected for downstream gene assembly. Significantly, this method improved a population of synthetic enhanced green fluorescent protein (720 bp) clones from 0.93% to 83.22%, corresponding to a decrease in the error frequency of synthetic gene from 11.44/kb to 0.46/kb. In addition, a parallel multiplex MICC error-removal strategy was also evaluated in assembling 11 genes encoding ∼21 kb of DNA from 893 oligos. The error frequency was reduced by 21.59-fold (from 14.25/kb to 0.66/kb), resulting in a 24.48-fold increase in the percentage of error-free assembled fragments (from 3.23% to 79.07%). Furthermore, the standard MICC error-removal process could be completed within 1.5 h at a cost as low as $0.374 per MICC.     

3'-modified oligonucleotides by reverse DNA synthesis

Reverse DNA oligonucleotide synthesis (i.e. from 5′→3′) is a strategy that has yet to be exploited fully. While utilized previously for the construction of alternating 3′-3′- and 5′-5′-linked antisense oligonucleotides, the use of nucleoside 5′-phosphoramidites has not generally been used for the elaboration of (modified) oligonucleotides. Presently, the potential of reverse oligonucleotide synthesis for the facile synthesis of 3′-modified DNAs is illustrated using a phosphoramidite derived from tyrosine. The derived oligonucleotide was shown to have chromatographic and electrophoretic properties identical with the modified oligonucleotide resulting from the proteinase K digestion of the vaccinia topoisomerase I–DNA covalent complex. The results confirm the nature of the structure previously assigned to this product, and establish the facility with which proteinase K is able to complete the digestion of the polypeptide backbone of the DNA oligonucleotide-linked topoisomerase I.     

High-quality,chromosome-scale genome assemblies: comparisons of three Diaphorina citri (Asian citrus psyllid) geographic populations

The Asian citrus psyllid, Diaphorina citri, is the insect vector of the causal agent of huanglongbing (HLB), a devastating bacterial disease of commercial citrus. Presently, few genomic resources exist for D. citri. In this study, we utilized PacBio HiFi and chromatin confirmation contact (Hi-C) sequencing to sequence, assemble, and compare three high-quality, chromosome-scale genome assemblies of D. citri collected from California, Taiwan, and Uruguay. Our assemblies had final sizes of 282.67 Mb (California), 282.89 Mb (Taiwan), and 266.67 Mb (Uruguay) assembled into 13 pseudomolecules—a reduction in assembly size of 41–45% compared with previous assemblies which we validated using flow cytometry. We identified the X chromosome in D. citri and annotated each assembly for repetitive elements, protein-coding genes, transfer RNAs, ribosomal RNAs, piwi-interacting RNA clusters, and endogenous viral elements. Between 19,083 and 20,357 protein-coding genes were predicted. Repetitive DNA accounts for 36.87–38.26% of each assembly. Comparative analyses and mitochondrial haplotype networks suggest that Taiwan and Uruguay D. citri are more closely related, while California D. citri are closely related to Florida D. citri. These high-quality, chromosome-scale assemblies provide new genomic resources to researchers to further D. citri and HLB research.     

A modified method for PCR-directed gene synthesis from large number of overlapping oligodeoxyribonucleotides

Here we report an improved, reproducible, simple, rapid, and cost-effective PCR-based DNA synthesis method using short (25–40 bp) overlapping oligodeoxyribonucleotides (oligos). The method involves two steps; (1) assembly of multiple/overlapping oligos by PCR to generate the template DNA and (2) amplification of the template DNA sequence with the two outermost oligos as primers. We have tested this method by synthesizing approximately 35 genes ranging in size between 300 bp and 1700 bp and G + C content from moderate (30%) to high (65%). In addition, we used the method to introduce 29 mutations simultaneously into a single gene. Key to the success of this method is the use of optimized oligo concentrations and the type of DNA polymerase used. This simplified and highly reproducible method is expected to be beneficial for the synthesis of a wide variety of genes.     

A Method for Multiplex Gene Synthesis Employing Error Correction Based on Expression

Our ability to engineer organisms with new biosynthetic pathways and genetic circuits is limited by the availability of protein characterization data and the cost of synthetic DNA. With new tools for reading and writing DNA, there are opportunities for scalable assays that more efficiently and cost effectively mine for biochemical protein characteristics. To that end, we have developed the Multiplex Library Synthesis and Expression Correction (MuLSEC) method for rapid assembly, error correction, and expression characterization of many genes as a pooled library. This methodology enables gene synthesis from microarray-synthesized oligonucleotide pools with a one-pot technique, eliminating the need for robotic liquid handling. Post assembly, the gene library is subjected to an ampicillin based quality control selection, which serves as both an error correction step and a selection for proteins that are properly expressed and folded in E. coli. Next generation sequencing of post selection DNA enables quantitative analysis of gene expression characteristics. We demonstrate the feasibility of this approach by building and testing over 90 genes for empirical evidence of soluble expression. This technique reduces the problem of part characterization to multiplex oligonucleotide synthesis and deep sequencing, two technologies under extensive development with projected cost reduction.     

Efficient and accurate whole genome assembly and methylome profiling of E. coli

Background

With the price of next generation sequencing steadily decreasing, bacterial genome assembly is now accessible to a wide range of researchers. It is therefore necessary to understand the best methods for generating a genome assembly, specifically, which combination of sequencing and bioinformatics strategies result in the most accurate assemblies. Here, we sequence three E. coli strains on the Illumina MiSeq, Life Technologies Ion Torrent PGM, and Pacific Biosciences RS. We then perform genome assemblies on all three datasets alone or in combination to determine the best methods for the assembly of bacterial genomes.

Results

Three E. coli strains – BL21(DE3), Bal225, and DH5α – were sequenced to a depth of 100× on the MiSeq and Ion Torrent machines and to at least 125× on the PacBio RS. Four assembly methods were examined and compared. The previously published BL21(DE3) genome [GenBank:{"type":"entrez-nucleotide","attrs":{"text":"AM946981.2","term_id":"313848522","term_text":"AM946981.2"}}AM946981.2], allowed us to evaluate the accuracy of each of the BL21(DE3) assemblies. BL21(DE3) PacBio-only assemblies resulted in a 90% reduction in contigs versus short read only assemblies, while N50 numbers increased by over 7-fold. Strikingly, the number of SNPs in PacBio-only assemblies were less than half that seen with short read assemblies (~20 SNPs vs. ~50 SNPs) and indels also saw dramatic reductions (~2 indel >5 bp in PacBio-only assemblies vs. ~12 for short-read only assemblies). Assemblies that used a mixture of PacBio and short read data generally fell in between these two extremes. Use of PacBio sequencing reads also allowed us to call covalent base modifications for the three strains. Each of the strains used here had a known covalent base modification genotype, which was confirmed by PacBio sequencing.

Conclusion

Using data generated solely from the Pacific Biosciences RS, we were able to generate the most complete and accurate de novo assemblies of E. coli strains. We found that the addition of other sequencing technology data offered no improvements over use of PacBio data alone. In addition, the sequencing data from the PacBio RS allowed for sensitive and specific calling of covalent base modifications.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-14-675) contains supplementary material, which is available to authorized users.     

Prion-like domains drive CIZ1 assembly formation at the inactive X chromosome

CIZ1 forms large assemblies at the inactive X chromosome (Xi) in female fibroblasts in an Xist lncRNA-dependent manner and is required for accurate maintenance of polycomb targets genome-wide. Here we address requirements for assembly formation and show that CIZ1 undergoes two direct interactions with Xist, via independent N- and C-terminal domains. Interaction with Xist, assembly at Xi, and complexity of self-assemblies formed in vitro are modulated by two alternatively spliced glutamine-rich prion-like domains (PLD1 and 2). PLD2 is dispensable for accumulation at existing CIZ1–Xi assemblies in wild-type cells but is required in CIZ1-null cells where targeting, assembly, and enrichment for H3K27me3 and H2AK119ub occur de novo. In contrast, PLD1 is required for both de novo assembly and accumulation at preexisting assemblies and, in vitro, drives formation of a stable fibrillar network. Together they impart affinity for RNA and a complex relationship with repeat E of Xist. These data show that alternative splicing of two PLDs modulates CIZ1’s ability to build large RNA–protein assemblies.     

A Rapid and Simple Method for DNA Engineering Using Cycled Ligation Assembly

DNA assembly techniques have developed rapidly, enabling efficient construction of complex constructs that would be prohibitively difficult using traditional restriction-digest based methods. Most of the recent methods for assembling multiple DNA fragments in vitro suffer from high costs, complex set-ups, and diminishing efficiency when used for more than a few DNA segments. Here we present a cycled ligation-based DNA assembly protocol that is simple, cheap, efficient, and powerful. The method employs a thermostable ligase and short Scaffold Oligonucleotide Connectors (SOCs) that are homologous to the ends and beginnings of two adjacent DNA sequences. These SOCs direct an exponential increase in the amount of correctly assembled product during a reaction that cycles between denaturing and annealing/ligating temperatures. Products of early cycles serve as templates for later cycles, allowing the assembly of many sequences in a single reaction. To demonstrate the method’s utility, we directed the assembly of twelve inserts, in one reaction, into a transformable plasmid. All the joints were precise, and assembly was scarless in the sense that no nucleotides were added or missing at junctions. Simple, efficient, and low-cost cycled ligation assemblies will facilitate wider use of complex genetic constructs in biomedical research.     

Directed evolution of green fluorescent protein by a new versatile PCR strategy for site-directed and semi-random mutagenesis

To develop a simple, speedy, economical and widely applicable method for multiple-site mutagenesis, we have substantially modified the Quik-Change™ Site-Directed Mutagenesis Kit protocol (Stratagene, La Jolla, CA). Our new protocol consists of (i) a PCR reaction using an in vitro technique, LDA (ligation-during-amplification), (ii) a DpnI treatment to digest parental DNA and to make megaprimers and (iii) a synthesis of double-stranded plasmid DNA for bacterial transformation. While the Quik Change™ Kit protocol introduces mutations at a single site, requiring two complementary mutagenic oligonucleotides, our new protocol requires only one mutagenic oligonucleotide for a mutation site, and can introduce mutations in a plasmid at multiple sites simultaneously. A targeting efficiency >70% was consistently achieved for multiple-site mutagenesis. Furthermore, the new protocol allows random mutagenesis with degenerative primers, because it does not use two complementary primers. Our mutagenesis strategy was successfully used to alter the fluorescence properties of green fluorescent protein (GFP), creating a new-color GFP mutant, cyan-green fluorescent protein (CGFP). An eminent feature of CGFP is its remarkable stability in a wide pH range (pH 4–12). The use of CGFP would allow us to monitor protein localization quantitatively in acidic organelles in secretory pathways.     

Whole-genome sequencing and assembly with high-throughput, short-read technologies

While recently developed short-read sequencing technologies may dramatically reduce the sequencing cost and eventually achieve the $1000 goal for re-sequencing, their limitations prevent the de novo sequencing of eukaryotic genomes with the standard shotgun sequencing protocol. We present SHRAP (SHort Read Assembly Protocol), a sequencing protocol and assembly methodology that utilizes high-throughput short-read technologies. We describe a variation on hierarchical sequencing with two crucial differences: (1) we select a clone library from the genome randomly rather than as a tiling path and (2) we sample clones from the genome at high coverage and reads from the clones at low coverage. We assume that 200 bp read lengths with a 1% error rate and inexpensive random fragment cloning on whole mammalian genomes is feasible. Our assembly methodology is based on first ordering the clones and subsequently performing read assembly in three stages: (1) local assemblies of regions significantly smaller than a clone size, (2) clone-sized assemblies of the results of stage 1, and (3) chromosome-sized assemblies. By aggressively localizing the assembly problem during the first stage, our method succeeds in assembling short, unpaired reads sampled from repetitive genomes. We tested our assembler using simulated reads from D. melanogaster and human chromosomes 1, 11, and 21, and produced assemblies with large sets of contiguous sequence and a misassembly rate comparable to other draft assemblies. Tested on D. melanogaster and the entire human genome, our clone-ordering method produces accurate maps, thereby localizing fragment assembly and enabling the parallelization of the subsequent steps of our pipeline. Thus, we have demonstrated that truly inexpensive de novo sequencing of mammalian genomes will soon be possible with high-throughput, short-read technologies using our methodology.     

Two Low Coverage Bird Genomes and a Comparison of Reference-Guided versus De Novo Genome Assemblies

As a greater number and diversity of high-quality vertebrate reference genomes become available, it is increasingly feasible to use these references to guide new draft assemblies for related species. Reference-guided assembly approaches may substantially increase the contiguity and completeness of a new genome using only low levels of genome coverage that might otherwise be insufficient for de novo genome assembly. We used low-coverage (∼3.5–5.5x) Illumina paired-end sequencing to assemble draft genomes of two bird species (the Gunnison Sage-Grouse, Centrocercus minimus, and the Clark''s Nutcracker, Nucifraga columbiana). We used these data to estimate de novo genome assemblies and reference-guided assemblies, and compared the information content and completeness of these assemblies by comparing CEGMA gene set representation, repeat element content, simple sequence repeat content, and GC isochore structure among assemblies. Our results demonstrate that even lower-coverage genome sequencing projects are capable of producing informative and useful genomic resources, particularly through the use of reference-guided assemblies.     

Parallel gene synthesis in a microfluidic device

The ability to synthesize custom de novo DNA constructs rapidly, accurately and inexpensively is highly desired by researchers, as synthetic genes and longer DNA constructs are enabling to numerous powerful applications in both traditional molecular biology and the emerging field of synthetic biology. However, the current cost of de novo synthesis—driven largely by reagent and handling costs—is a significant barrier to the widespread availability of such technology. In this work, we demonstrate, to our knowledge, the first gene synthesis in a microfluidic environment. The use of microfluidic technology greatly reduces reaction volumes and the corresponding reagent and handling costs. Additionally, microfluidic technology enables large numbers of complex reactions to be performed in parallel. Here, we report the fabrication of a multi-chamber microfluidic device and its use in carrying out the syntheses of several DNA constructs. Genes up to 1kb in length were synthesized in parallel at minute starting oligonucleotide concentrations (10–25nM) in four 500nl reactors. Such volumes are one to two orders of magnitude lower than those utilized in conventional gene synthesis. The identity of all target genes was verified by sequencing, and the resultant error rate was determined to be 1 per 560 bases.     

Experimental analysis of gene assembly with TopDown one-step real-time gene synthesis

Herein we present a simple, cost-effective TopDown (TD) gene synthesis method that eliminates the interference between the polymerase chain reactions (PCR) assembly and amplification in one-step gene synthesis. The method involves two key steps: (i) design of outer primers and assembly oligonucleotide set with a melting temperature difference of >10°C and (ii) utilization of annealing temperatures to selectively control the efficiencies of oligonucleotide assembly and full-length template amplification. In addition, we have combined the proposed method with real-time PCR to analyze the step-wise efficiency and the kinetics of the gene synthesis process. Gel electrophoresis results are compared with real-time fluorescence signals to investigate the effects of oligonucleotide concentration, outer primer concentration, stringency of annealing temperature, and number of PCR cycles. Analysis of the experimental results has led to insights into the gene synthesis process. We further discuss the conditions for preventing the formation of spurious DNA products. The TD real-time gene synthesis method provides a simple and efficient method for assembling fairly long DNA sequence, and aids in optimizing gene synthesis conditions. To our knowledge, this is the first report that utilizes real-time PCR for gene synthesis.     

Combinatorial Pooling Enables Selective Sequencing of the Barley Gene Space

For the vast majority of species – including many economically or ecologically important organisms, progress in biological research is hampered due to the lack of a reference genome sequence. Despite recent advances in sequencing technologies, several factors still limit the availability of such a critical resource. At the same time, many research groups and international consortia have already produced BAC libraries and physical maps and now are in a position to proceed with the development of whole-genome sequences organized around a physical map anchored to a genetic map. We propose a BAC-by-BAC sequencing protocol that combines combinatorial pooling design and second-generation sequencing technology to efficiently approach denovo selective genome sequencing. We show that combinatorial pooling is a cost-effective and practical alternative to exhaustive DNA barcoding when preparing sequencing libraries for hundreds or thousands of DNA samples, such as in this case gene-bearing minimum-tiling-path BAC clones. The novelty of the protocol hinges on the computational ability to efficiently compare hundred millions of short reads and assign them to the correct BAC clones (deconvolution) so that the assembly can be carried out clone-by-clone. Experimental results on simulated data for the rice genome show that the deconvolution is very accurate, and the resulting BAC assemblies have high quality. Results on real data for a gene-rich subset of the barley genome confirm that the deconvolution is accurate and the BAC assemblies have good quality. While our method cannot provide the level of completeness that one would achieve with a comprehensive whole-genome sequencing project, we show that it is quite successful in reconstructing the gene sequences within BACs. In the case of plants such as barley, this level of sequence knowledge is sufficient to support critical end-point objectives such as map-based cloning and marker-assisted breeding.     

Chromosome-Specific Painting in Cucumis Species Using Bulked Oligonucleotides

Chromosome-specific painting is a powerful technique in molecular cytogenetic and genome research. We developed an oligonucleotide (oligo)-based chromosome painting technique in cucumber (Cucumis sativus) that will be applicable in any plant species with a sequenced genome. Oligos specific to a single chromosome of cucumber were identified using a newly developed bioinformatic pipeline and then massively synthesized de novo in parallel. The synthesized oligos were amplified and labeled with biotin or digoxigenin for use in fluorescence in situ hybridization (FISH). We developed three different probes with each containing 23,000–27,000 oligos. These probes spanned 8.3–17 Mb of DNA on targeted cucumber chromosomes and had the densities of 1.5–3.2 oligos per kilobases. These probes produced FISH signals on a single cucumber chromosome and were used to paint homeologous chromosomes in other Cucumis species diverged from cucumber for up to 12 million years. The bulked oligo probes allowed us to track a single chromosome in early stages during meiosis. We were able to precisely map the pairing between cucumber chromosome 7 and chromosome 1 of Cucumis hystrix in a F₁ hybrid. These two homeologous chromosomes paired in 71% of prophase I cells but only 25% of metaphase I cells, which may provide an explanation of the higher recombination rates compared to the chiasma frequencies between homeologous chromosomes reported in plant hybrids.     

Paired-End Sequencing of Long-Range DNA Fragments for De Novo Assembly of Large,Complex Mammalian Genomes by Direct Intra-Molecule Ligation

Background

The relatively short read lengths from next generation sequencing (NGS) technologies still pose a challenge for de novo assembly of complex mammal genomes. One important solution is to use paired-end (PE) sequence information experimentally obtained from long-range DNA fragments (>1 kb). Here, we characterize and extend a long-range PE library construction method based on direct intra-molecule ligation (or molecular linker-free circularization) for NGS.

Results

We found that the method performs stably for PE sequencing of 2- to 5- kb DNA fragments, and can be extended to 10–20 kb (and even in extremes, up to ∼35 kb). We also characterized the impact of low quality input DNA on the method, and develop a whole-genome amplification (WGA) based protocol using limited input DNA (<1 µg). Using this PE dataset, we accurately assembled the YanHuang (YH) genome, the first sequenced Asian genome, into a scaffold N50 size of >2 Mb, which is over100-times greater than the initial size produced with only small insert PE reads(17 kb). In addition, we mapped two 7- to 8- kb insertions in the YH genome using the larger insert sizes of the long-range PE data.

Conclusions

In conclusion, we demonstrate here the effectiveness of this long-range PE sequencing method and its use for the de novo assembly of a large, complex genome using NGS short reads.     

Unbiased Deep Sequencing of RNA Viruses from Clinical Samples

Here we outline a next-generation RNA sequencing protocol that enables de novo assemblies and intra-host variant calls of viral genomes collected from clinical and biological sources. The method is unbiased and universal; it uses random primers for cDNA synthesis and requires no prior knowledge of the viral sequence content. Before library construction, selective RNase H-based digestion is used to deplete unwanted RNA — including poly(rA) carrier and ribosomal RNA — from the viral RNA sample. Selective depletion improves both the data quality and the number of unique reads in viral RNA sequencing libraries. Moreover, a transposase-based ''tagmentation'' step is used in the protocol as it reduces overall library construction time. The protocol has enabled rapid deep sequencing of over 600 Lassa and Ebola virus samples-including collections from both blood and tissue isolates-and is broadly applicable to other microbial genomics studies.