Targeted Sequencing of Large Genomic Regions with CATCH-Seq

Current target enrichment systems for large-scale next-generation sequencing typically require synthetic oligonucleotides used as capture reagents to isolate sequences of interest. The majority of target enrichment reagents are focused on gene coding regions or promoters en masse. Here we introduce development of a customizable targeted capture system using biotinylated RNA probe baits transcribed from sheared bacterial artificial chromosome clone templates that enables capture of large, contiguous blocks of the genome for sequencing applications. This clone adapted template capture hybridization sequencing (CATCH-Seq) procedure can be used to capture both coding and non-coding regions of a gene, and resolve the boundaries of copy number variations within a genomic target site. Furthermore, libraries constructed with methylated adapters prior to solution hybridization also enable targeted bisulfite sequencing. We applied CATCH-Seq to diverse targets ranging in size from 125 kb to 3.5 Mb. Our approach provides a simple and cost effective alternative to other capture platforms because of template-based, enzymatic probe synthesis and the lack of oligonucleotide design costs. Given its similarity in procedure, CATCH-Seq can also be performed in parallel with commercial systems.     

Long-Range Genomic Enrichment,Sequencing, and Assembly to Determine Unknown Sequences Flanking a Known microRNA

Conserved plant microRNAs (miRNAs) modulate important biological processes but little is known about conserved cis-regulatory elements (CREs) surrounding MIRNA genes. We developed a solution-based targeted genomic enrichment methodology to capture, enrich, and sequence flanking genomic regions surrounding conserved MIRNA genes with a locked-nucleic acid (LNA)-modified, biotinylated probe complementary to the mature miRNA sequence. Genomic DNA bound by the probe is captured by streptavidin-coated magnetic beads, amplified, sequenced and assembled de novo to obtain genomic DNA sequences flanking MIRNA locus of interest. We demonstrate the sensitivity and specificity of this enrichment methodology in Arabidopsis thaliana to enrich targeted regions spanning 10–20 kb surrounding known MIR166 and MIR165 loci. Assembly of the sequencing reads successfully recovered all targeted loci. While further optimization for larger, more complex genomes is needed, this method may enable determination of flanking genomic DNA sequence surrounding a known core (like a conserved mature miRNA) from multiple species that currently don''t have a full genome assembly available.     

Supported PCR: an efficient procedure to amplify sequences flanking a known DNA segment

We describe a novel modification of the polymerase chain reaction for efficient in vitro amplification of genomic DNA sequences flanking short stretches of known sequence. The technique utilizes a target enrichment step, based on the selective isolation of biotinylated fragments from the bulk of genomic DNA on streptavidin-containing support. Subsequently, following ligation with a second universal linker primer, the selected fragments can be amplified to amounts suitable for further molecular studies. The procedure has been applied to recover T-DNA flanking sequences in transgenic tomato plants which could subsequently be used to assign the positions of T-DNA to the molecular map of tomato. The method called supported PCR (sPCR) is a simple and efficient alternative to techniques used in the isolation of specific sequences flanking a known DNA segment.     

Pulling out the 1%: Whole-Genome Capture for the Targeted Enrichment of Ancient DNA Sequencing Libraries

Most ancient specimens contain very low levels of endogenous DNA, precluding the shotgun sequencing of many interesting samples because of cost. Ancient DNA (aDNA) libraries often contain <1% endogenous DNA, with the majority of sequencing capacity taken up by environmental DNA. Here we present a capture-based method for enriching the endogenous component of aDNA sequencing libraries. By using biotinylated RNA baits transcribed from genomic DNA libraries, we are able to capture DNA fragments from across the human genome. We demonstrate this method on libraries created from four Iron Age and Bronze Age human teeth from Bulgaria, as well as bone samples from seven Peruvian mummies and a Bronze Age hair sample from Denmark. Prior to capture, shotgun sequencing of these libraries yielded an average of 1.2% of reads mapping to the human genome (including duplicates). After capture, this fraction increased substantially, with up to 59% of reads mapped to human and enrichment ranging from 6- to 159-fold. Furthermore, we maintained coverage of the majority of regions sequenced in the precapture library. Intersection with the 1000 Genomes Project reference panel yielded an average of 50,723 SNPs (range 3,062–147,243) for the postcapture libraries sequenced with 1 million reads, compared with 13,280 SNPs (range 217–73,266) for the precapture libraries, increasing resolution in population genetic analyses. Our whole-genome capture approach makes it less costly to sequence aDNA from specimens containing very low levels of endogenous DNA, enabling the analysis of larger numbers of samples.     

Targeted Genome-Wide Enrichment of Functional Regions

Only a small fraction of large genomes such as that of the human contains the functional regions such as the exons, promoters, and polyA sites. A platform technique for selective enrichment of functional genomic regions will enable several next-generation sequencing applications that include the discovery of causal mutations for disease and drug response. Here, we describe a powerful platform technique, termed “functional genomic fingerprinting” (FGF), for the multiplexed genomewide isolation and analysis of targeted regions such as the exome, promoterome, or exon splice enhancers. The technique employs a fixed part of a uniquely designed Fixed-Randomized primer, while the randomized part contains all the possible sequence permutations. The Fixed-Randomized primers bind with full sequence complementarity at multiple sites where the fixed sequence (such as the splice signals) occurs within the genome, and multiplex amplify many regions bounded by the fixed sequences (e.g., exons). Notably, validation of this technique using cardiac myosin binding protein-C (MYBPC3) gene as an example strongly supports the application and efficacy of this method. Further, assisted by genomewide computational analyses of such sequences, the FGF technique may provide a unique platform for high-throughput sample production and analysis of targeted genomic regions by the next-generation sequencing techniques, with powerful applications in discovering disease and drug response genes.     

Sequence capture using AFLP‐generated baits: A cost‐effective method for high‐throughput phylogenetic and phylogeographic analysis

Target sequence capture is an efficient technique to enrich specific genomic regions for high‐throughput sequencing in ecological and evolutionary studies. In recent years, many sequence capture approaches have been proposed, but most of them rely on commercial synthetic baits which make the experiment expensive. Here, we present a novel sequence capture approach called AFLP‐based genome sequence capture (AFLP Capture). This method uses the AFLP (amplified fragment length polymorphism) technique to generate homemade capture baits without the need for prior genome information, thus is applicable to any organisms. In this approach, biotinylated AFLP fragments representing a random fraction of the genome are used as baits to capture the homologous fragments from genomic shotgun sequencing libraries. In a trial study, by using AFLP Capture, we successfully obtained 511 orthologous loci (>700,000 bp in total length) from 11 Odorrana species and more than 100,000 single nucleotide polymorphisms (SNPs) in four analyzed individuals of an Odorrana species. This result shows that our method can be used to address questions of various evolutionary depths (from interspecies level to intraspecies level). We also discuss the flexibility in bait preparation and how the sequencing data are analyzed. In summary, AFLP Capture is a rapid and flexible tool and can significantly reduce the experimental cost for phylogenetic studies that require analyzing genome‐scale data (hundreds or thousands of loci).     

Segments missing from the draft human genome sequence can be isolated by transformation-associated recombination cloning in yeast

The reported draft human genome sequence includes many contigs that are separated by gaps of unknown sequence. These gaps may be due to chromosomal regions that are not present in the Escherichia coli libraries used for DNA sequencing because they cannot be cloned efficiently, if at all, in bacteria. Using a yeast artificial chromosome (YAC)/ bacterial artificial chromosome (BAC) library generated in yeast, we found that approximately 6% of human DNA sequences tested transformed E. coli cells less efficiently than yeast cells, and were less stable in E. coli than in yeast. When the ends of several YAC/BAC isolates cloned in yeast were sequenced and compared with the reported draft sequence, major inconsistencies were found with the sequences of those YAC/BAC isolates that transformed E. coli cells inefficiently. Two human genomic fragments were re-isolated from human DNA by transformation-associated recombination (TAR) cloning. Re-sequencing of these regions showed that the errors in the draft are the results of both missassembly and loss of specific DNA sequences during cloning in E. coli. These results show that TAR cloning might be a valuable method that could be widely used during the final stages of the Human Genome Project.     

Identification and isolation of Mu-flanking fragments from maize

Transposable elements have been utilized as mutagens to create mutant libraries for functional genomics.Isolation of genomic seg-ments flanking the insertion Mutator (Mu) is a key step in insertion mutagenesis studies.Herein,we adopted a modified AFLP method to identify and isolate Mu-flanking fragments from maize.The method consists of the following steps: 1) double-digestion of genomic DNA with Bgl ⅡMsp Ⅰ and ligation of digested fragments to the Bgl Ⅱ- and Msp Ⅰ-adaptors; 2) enrichment of a subset of Bgl Ⅱ/Msp Ⅰ fragments followed by selective amplification of the Mu-flanking fragments; 3) simultaneous display of AFLP bands derived from the flanking re-gions for both insert and native Mu transposons; 4) identification and isolation of AFLP bands resulting from Mu insertions by comparing the banding profiles between Mu-induced mutants and their parental lines; and 5) confirmation of flanking fragments related to these Mu insertions.Using this approach,we have isolated flanking fragment(s) resulting from Mu insertion for every Mu-indueed mutant,and one such fragment,M196-FF,is found to contain a partial sequence of the DNA topoisomerase Ⅰ gene Topl.Moreover,the modified AFLP method including all restriction enzymes,adaptors and primers has been optimized in this study.The modified AFLP method has been proved to be simple and efficient in the isolation of Mu-flanking fragments and will find its usefulness in the functional genomics of maize.     

一种改良的启动子序列克隆的染色体步查法

利用染色体步行法,从已知DNA序列克隆侧翼未知序列是非常有效的方法之一,但由于所选用的特定限制性内切酶对目标基因组不能酶解成合适大小的片段,因而受PCR扩增能力的局限,往往扩增不出有效产物. 针对这一点,这里我们介绍一种简单有效的改良方法,它包括以下步骤：首先用不同的限制性内切酶(包括平末端和粘性末端) 酶解目标基因组DNA,接着,选择能将基因组酶切成弥散、分布均匀的限制性内切酶,如DraⅠ和HindⅢ,合成相对应的接头;然后,选择弥散的、分布均匀的限制性内切酶的酶解产物,构建成含相应接头的基因组DNA文库,用作PCR的模板;最后,用接头引物和特异引物,通过巢式PCR扩增目的片段,获得了理想的扩增效果.采用改进后的染色体步查法,有效地从较复杂的棉花核DNA中克隆出6个棉花启动子序列.     

High-throughput cloning of human liver complete open reading frames using homologous recombination in Escherichia coli

In this article, we describe a high-throughput cloning method, seamless enzyme-free cloning (SEFC), which allows one-step assembly of DNA fragments in vivo via homologous recombination in Escherichia coli. In the method, the desired open reading frame (ORF) is amplified by use of ORF-specific primers with flanking sequences identical to the two ends of a linearized vector. The polymerase chain reaction (PCR) product and the linearized vector are then cotransformed into E. coli cells, where the ORF is incorporated into the vector in vivo. SEFC is a simple, reliable, and inexpensive method of cloning in which PCR fragments are fused into expression vectors without unwanted amino acids or extra in vitro manipulations apart from the single PCR amplification step. Using this method, we successfully cloned human liver complete ORFs into the yeast AD and DB vectors and generated a clone resource of 4964 AD-ORFs and 4676 DB-ORFs in 3 months. This approach will be useful for daily DNA cloning and for creating proteome-scale clone resources.     

用于高通量测序的基因组靶序列捕获方法的建立

文章旨在建立一种基因组目标靶序列捕捉文库的方法,并结合第二代测序技术,以实现候选基因区段的深度测序。利用Agilent公司的eArray在线平台,对1250个基因的11824个外显子共2414977bp的基因组序列进行120个碱基长度的捕捉探针(钓饵)设计,并制备成SureSelect液相靶序列捕获试剂。选用2例人基因组DNA,超声打断后末端补平并磷酸化,连接SOLiD接头,回收150bp～200bp的DNA片段,与靶序列探针杂交捕获目标序列,油包水微乳滴PCR扩增后,磁珠分离富集,上SOLiD测序系统通过工作流程分析(WFA)进行文库质量的评价,或正式测序反应。结果显示对所包含的11147个基因外显子片段设计出并合成了46509个捕捉探针,制备成SureSelect试剂盒。探针可有效地捕捉并富集基因组DNA的目标靶片段,定量PCR显示富集效率可达29倍。WFA分析表明文库可以在SOLiD仪器进行正式测序。测序结果显示靶序列区域的测序数占有效总测序数的比例达到70%,覆盖率均在200×以上。结果表明本研究所建立的SureSelect基因组靶序列捕捉、富集建立测序文库的技术路线可行,可直接用于SOLiD测序仪的测序。     

Specific enrichment of nonribosomal peptide synthetase module by an affinity probe for adenylation domains

We targeted the development of an affinity probe for adenylation (A) domains that can facilitate enrichment, identification, and quantification of A domain-containing modules in nonribosomal peptide synthetase (NRPS)-polyketide synthase (PKS) hybrids and NRPSs. A 5′-O-sulfamoyladenosine (AMS) non-hydrolyzable analogue of adenosine monophosphate (AMP) has been reported as a scaffold for the design of inhibitors exhibiting tight binding of adenylation enzymes. Here we describe the application of an affinity probe for A domains. Our synthetic probe, a biotinylated l-Phe-AMS (l-Phe-AMS-biotin) specifically targets the A domains in NRPS modules that activates l-Phe to an aminoacyladenylate intermediate in both recombinant NRPS enzyme systems and whole proteomes.     

Targeted multiplex next‐generation sequencing: advances in techniques of mitochondrial and nuclear DNA sequencing for population genomics

Next‐generation sequencing (NGS) is emerging as an efficient and cost‐effective tool in population genomic analyses of nonmodel organisms, allowing simultaneous resequencing of many regions of multi‐genomic DNA from multiplexed samples. Here, we detail our synthesis of protocols for targeted resequencing of mitochondrial and nuclear loci by generating indexed genomic libraries for multiplexing up to 100 individuals in a single sequencing pool, and then enriching the pooled library using custom DNA capture arrays. Our use of DNA sequence from one species to capture and enrich the sequencing libraries of another species (i.e. cross‐species DNA capture) indicates that efficient enrichment occurs when sequences are up to about 12% divergent, allowing us to take advantage of genomic information in one species to sequence orthologous regions in related species. In addition to a complete mitochondrial genome on each array, we have included between 43 and 118 nuclear loci for low‐coverage sequencing of between 18 kb and 87 kb of DNA sequence per individual for single nucleotide polymorphisms discovery from 50 to 100 individuals in a single sequencing lane. Using this method, we have generated a total of over 500 whole mitochondrial genomes from seven cetacean species and green sea turtles. The greater variation detected in mitogenomes relative to short mtDNA sequences is helping to resolve genetic structure ranging from geographic to species‐level differences. These NGS and analysis techniques have allowed for simultaneous population genomic studies of mtDNA and nDNA with greater genomic coverage and phylogeographic resolution than has previously been possible in marine mammals and turtles.     

SLiCE: a novel bacterial cell extract-based DNA cloning method

We describe a novel cloning method termed SLiCE (Seamless Ligation Cloning Extract) that utilizes easy to generate bacterial cell extracts to assemble multiple DNA fragments into recombinant DNA molecules in a single in vitro recombination reaction. SLiCE overcomes the sequence limitations of traditional cloning methods, facilitates seamless cloning by recombining short end homologies (≥15 bp) with or without flanking heterologous sequences and provides an effective strategy for directional subcloning of DNA fragments from Bacteria Artificial Chromosomes (BACs) or other sources. SLiCE is highly cost effective as a number of standard laboratory bacterial strains can serve as sources for SLiCE extract. In addition, the cloning efficiencies and capabilities of these strains can be greatly improved by simple genetic modifications. As an example, we modified the DH10B Escherichia coli strain to express an optimized λ prophage Red recombination system. This strain, termed PPY, facilitates SLiCE with very high efficiencies and demonstrates the versatility of the method.     

Characterizing seamless ligation cloning extract for synthetic biological applications

Synthetic biology aims at designing and engineering organisms. The engineering process typically requires the establishment of suitable DNA constructs generated through fusion of multiple protein coding and regulatory sequences. Conventional cloning techniques, including those involving restriction enzymes and ligases, are often of limited scope, in particular when many DNA fragments must be joined or scar-free fusions are mandatory. Overlap-based-cloning methods have the potential to overcome such limitations. One such method uses seamless ligation cloning extract (SLiCE) prepared from Escherichia coli cells for straightforward and efficient in vitro fusion of DNA fragments. Here, we systematically characterized extracts prepared from the unmodified E. coli strain DH10B for SLiCE-mediated cloning and determined DNA sequence-associated parameters that affect cloning efficiency. Our data revealed the virtual absence of length restrictions for vector backbone (up to 13.5 kbp) and insert (90 bp to 1.6 kbp). Furthermore, differences in GC content in homology regions are easily tolerated and the deletion of unwanted vector sequences concomitant with targeted fragment insertion is straightforward. Thus, SLiCE represents a highly versatile DNA fusion method suitable for cloning projects in virtually all molecular and synthetic biology projects.     

Generation and characterization of biotinylated recombinant Fab antibody fragment against cortisol

This paper describes the in vivo generation method of biotinylated recombinant Fab antibody fragments. The original molecular vectors for Escherichia coli Fab fragments expression were designed. In vivo biotinylated recombinant antibody Fab fragment against cortisol was generated. The kinetic parameters of interaction of these antibody fragments with cortisol-BSA complex were measured via the biolayer interferometry method. An equilibrium dissociation constant of this interaction K_D is 5.48 × 10^–10 M. The interaction is reversible and it could be competitively inhibited by free cortisol addition. These results could be used in generating of immunoassays for quantitative cortisol determinations. The in vivo biotinylation system under review is universal and suitable for expression of any biotinylated Fab fragments in the E. coli system.     

Performance of microarray and liquid based capture methods for target enrichment for massively parallel sequencing and SNP discovery

Targeted sequencing is a cost-efficient way to obtain answers to biological questions in many projects, but the choice of the enrichment method to use can be difficult. In this study we compared two hybridization methods for target enrichment for massively parallel sequencing and single nucleotide polymorphism (SNP) discovery, namely Nimblegen sequence capture arrays and the SureSelect liquid-based hybrid capture system. We prepared sequencing libraries from three HapMap samples using both methods, sequenced the libraries on the Illumina Genome Analyzer, mapped the sequencing reads back to the genome, and called variants in the sequences. 74-75% of the sequence reads originated from the targeted region in the SureSelect libraries and 41-67% in the Nimblegen libraries. We could sequence up to 99.9% and 99.5% of the regions targeted by capture probes from the SureSelect libraries and from the Nimblegen libraries, respectively. The Nimblegen probes covered 0.6 Mb more of the original 3.1 Mb target region than the SureSelect probes. In each sample, we called more SNPs and detected more novel SNPs from the libraries that were prepared using the Nimblegen method. Thus the Nimblegen method gave better results when judged by the number of SNPs called, but this came at the cost of more over-sampling.     

Rapid generation of incremental truncation libraries for protein engineering using α-phosphothioate nucleotides

Incremental truncation for the creation of hybrid enzymes (ITCHY) is a novel tool for the generation of combinatorial libraries of hybrid proteins independent of DNA sequence homology. We herein report a fundamentally different methodology for creating incremental truncation libraries using nucleotide triphosphate analogs. Central to the method is the polymerase catalyzed, low frequency, random incorporation of α-phosphothioate dNTPs into the region of DNA targeted for truncation. The resulting phosphothioate internucleotide linkages are resistant to 3′→5′ exonuclease hydrolysis, rendering the target DNA resistant to degradation in a subsequent exonuclease III treatment. From an experimental perspective the protocol reported here to create incremental truncation libraries is simpler and less time consuming than previous approaches by combining the two gene fragments in a single vector and eliminating additional purification steps. As proof of principle, an incremental truncation library of fusions between the N-terminal fragment of Escherichia coli glycinamide ribonucleotide formyltransferase (PurN) and the C-terminal fragment of human glycinamide ribonucleotide formyltransferase (hGART) was prepared and successfully tested for functional hybrids in an auxotrophic E.coli host strain. Multiple active hybrid enzymes were identified, including ones fused in regions of low sequence homology.     

Use of competitive DNA hybridization to identify differences in the genomes of bacteria

Although recent technological advances in DNA sequencing and computational biology now allow scientists to compare entire microbial genomes, comparisons of closely related bacterial species and individual isolates by whole-genome sequencing approaches remains prohibitively expensive for most laboratories. Here we report the development and testing of a biochemical approach for targeted sequencing of only those chromosomal regions that differ between two DNA preparations. The method, designated GFE (genome fragment enrichment) uses competitive solution hybridization and positive selection to obtain genomic DNA fragments that are present in one pool of fragments but not another. Repeated comparisons of the genomes of Enterococcus faecalis and E. faecium led to the identification of 225 putative genome-specific DNA fragments. Species and strain variations within these fragments were confirmed by both experimental and bioinformatic analyses. The E. faecalis genome-specific sequences identified included both a preponderance of those predicted to encode surface-exposed proteins, as well as several previously described unique marker regions embedded within highly conserved rrn operons. The GFE strategy we describe efficiently identified genomic differences between two enterococcal genomes, and will be widely applicable for studying genetic variation among closely related bacterial species.