Shotgun sequencing,clone pooling,and comparative strategies for mapping and sequencing

Comparative genomic sequencing and analysis offers new wealth of information for target selection and the development of therapeutics. This article focuses on the following two key innovations in mapping and sequencing: first, shotgun sequencing of clone pools to combine the benefits of whole-genome shotgun and clone-by-clone strategies, and second, the leveraging of newly available assembled genomic sequences to improve the effectiveness of new sequencing projects through comparative mapping and comparative sequence assembly. The following specific sequencing and mapping methods are discussed in detail: clone-array pooled shotgun sequencing (CAPSS); transversal shotgun pooling designs; clone-array pooled shotgun mapping (CAPS-MAP); pooled genomic indexing (PGI); short-tag pooled genomic indexing (ST-PGI); and comparative sequence assembly (the CSA™ method). The methods can be implemented with only modest modifications of current large-scale sequencing pipelines and are highly synergistic with the next generation of sequencing technologies.     

Integration of animal linkage and BAC contig maps using overgo hybridization

The alignment of genome linkage maps, defined primarily by segregation of sequence-tagged site (STS) markers, with BAC contig physical maps and full genome sequences requires high throughput mechanisms to identify BAC clones that contain specific STS. A powerful technique for this purpose is multi-dimensional hybridization of "overgo" probes. The probes are chosen from available STS sequence data by selecting unique probe sequences that have a common melting temperature. We have hybridized sets of 216 overgo probes in subset pools of 36 overgos at a time to filter-spotted chicken BAC clone arrays. A four-dimensional pooling strategy, including one degree of redundancy, has been employed. This requires 24 hybridizations to completely assign BACs for all 216 probes. Results to date are consistent with about a 10% failure rate in overgo probe design and a 15-20% false negative detection rate within a group of 216 markers. Three complete rounds of overgo hybridization, each to sets of about 39,000 BACs (either BAMHI or ECORI partial digest inserts) generated a total of 1853 BAC alignments for 517 mapped chicken genome STS markers. These data are publicly available, and they have been used in the assembly of a first generation BAC contig map of the chicken genome.     

A Positive Detecting Code and Its Decoding Algorithm for DNA Library Screening

The study of gene functions requires high-quality DNA libraries. However, a large number of tests and screenings are necessary for compiling such libraries. We describe an algorithm for extracting as much information as possible from pooling experiments for library screening. Collections of clones are called pools, and a pooling experiment is a group test for detecting all positive clones. The probability of positiveness for each clone is estimated according to the outcomes of the pooling experiments. Clones with high chance of positiveness are subjected to confirmatory testing. In this paper, we introduce a new positive clone detecting algorithm, called the Bayesian network pool result decoder (BNPD). The performance of BNPD is compared, by simulation, with that of the Markov chain pool result decoder (MCPD) proposed by Knill et al. in 1996. Moreover, the combinatorial properties of pooling designs suitable for the proposed algorithm are discussed in conjunction with combinatorial designs and dhbox{-}{rm disjunct} matrices. We also show the advantage of utilizing packing designs or BIB designs for the BNPD algorithm.     

Error-tolerant pooling designs with inhibitors.

Pooling designs are used in clone library screening to efficiently distinguish positive clones from negative clones. Mathematically, a pooling design is just a nonadaptive group testing scheme which has been extensively studied in the literature. In some applications, there is a third category of clones called "inhibitors" whose effect is to neutralize positives. Specifically, the presence of an inhibitor in a pool dictates a negative outcome even though positives are present. Sequential group testing schemes, which can be modified to three-stage schemes, have been proposed for the inhibitor model, but it is unknown whether a pooling design (a one-stage scheme) exists. Another open question raised in the literature is whether the inhibitor model can treat unreliable pool outcomes. In this paper, we answer both open problems by giving a pooling design, as well as a two-stage scheme, for the inhibitor model with unreliable outcomes. The number of pools required by our schemes are quite comparable to the three-stage scheme.     

Comparative BAC‐based physical mapping of Oryza sativa ssp. indica var. 93–11 and evaluation of the two rice reference sequence assemblies

Reference sequences are sequences that are used for public consultation, and therefore must be of high quality. Using the whole‐genome shotgun/next‐generation sequencing approach, many genome sequences of complex higher plants have been generated in recent years, and are generally considered reference sequences. However, none of these sequences has been experimentally evaluated at the whole‐genome sequence assembly level. Rice has a relatively simple plant genome, and the genome sequences for its two sub‐species obtained using different sequencing approaches were published approximately 10 years ago. This provides a unique system for a case study to evaluate the qualities and utilities of published plant genome sequences. We constructed a robust BAC physical map embedding a large number of BAC end sequences forrice variety 93–11. Through BAC end sequence alignments and tri‐assembly comparisons of the 93–11 physical map and the two reference sequences, we found that the Nipponbare reference sequence generated using the clone‐by‐clone approach has a high quality but still contains small artifact inversions and missing sequences. In contrast, the 93–11 reference sequence generated using the whole‐genome shotgun approach contains many large and varied assembly errors, such as inversions, duplications and translocations, as well as missing sequences. The 93–11 physical map provides an invaluable resource for evaluation and improvements toward completion of both Nipponbare and 93–11 reference sequences.     

Comparative analysis of BAC and whole genome shotgun sequences from an Anopheles gambiae region related to Plasmodium encapsulation

The only natural mechanism of malaria transmission in sub-Saharan Africa is the mosquito, generally Anopheles gambiae. Blocking malaria parasite transmission by stopping the development of Plasmodium in the insect vector would provide a useful alternative to the current methods of malaria control. Toward this end, it is important to understand the molecular basis of the malaria parasite refractory phenotype in An. gambiae mosquito strains. We have selected and sequenced six bacterial artificial chromosome (BAC) clones from the Pen-1 region that is the major quantitative trait locus involved in Plasmodium encapsulation. The sequence and the annotation of five overlapping BAC clones plus one adjacent, but not contiguous clone, totaling 585kb of genomic sequence from the centromeric end of the Pen-1 region of the PEST strain were compared to that of the genome sequence of the same strain produced by the whole genome shotgun technique. This project identified 23 putative mosquito genes plus putative copies of the retrotransposable elements BEL12 and TRANSIBN1_AG in the six BAC clones. Nineteen of the predicted genes are most similar to their Drosophila melanogaster homologs while one is more closely related to vertebrate genes. Comparison of these new BAC sequences plus previously published BAC sequences to the cognate region of the assembled genome sequence identified three retrotransposons present in one sequence version but not the other. One of these elements, Indy, has not been previously described. These observations provide evidence for the recent active transposition of these elements and demonstrate the plasticity of the Anopheles genome. The BAC sequences strongly support the public whole genome shotgun assembly and automatic annotation while also demonstrating the benefit of complementary genome sequences and of human curation. Importantly, the data demonstrate the differences in the genome sequence of an individual mosquito compared to that of a hypothetical, average genome sequence generated by whole genome shotgun assembly.     

The Design of Pooling Experiments for Screening a Clone Map

We consider nonadaptive pooling designs for unique-sequence screening of a 1530-clone map ofAspergillus nidulans.The map has the properties that the clones are, with possibly a few exceptions, ordered and no more than 2 of them cover any point on the genome. We propose two subdesigns of the Steiner systemS(3, 5, 65), one with 65 pools and approximately 118 clones per pool, the other with 54 pools and about 142 clones per pool. Each design allows 1 or 2 positive clones to be detected, even in the presence of substantial experimental error rates. More efficient designs are possible if the overlap information in the map is exploited, if there is no constraint on the number of clones in a pool, and if no error tolerance is required. An information theory lower bound requires at least 12 pools to satisfy these minimal criteria, and an “interleaved binary” design can be constructed on 20 pools, with about 380 clones per pool. However, the designs with more pools have important properties of robustness to various possible errors and general applicability to a wider class of pooling experiments.     

Multiplex sequencing of bacterial artificial chromosomes for assembling complex plant genomes

Hierarchical shotgun sequencing remains the method of choice for assembling high‐quality reference sequences of complex plant genomes. The efficient exploitation of current high‐throughput technologies and powerful computational facilities for large‐insert clone sequencing necessitates the sequencing and assembly of a large number of clones in parallel. We developed a multiplexed pipeline for shotgun sequencing and assembling individual bacterial artificial chromosomes (BACs) using the Illumina sequencing platform. We illustrate our approach by sequencing 668 barley BACs (Hordeum vulgare L.) in a single Illumina HiSeq 2000 lane. Using a newly designed parallelized computational pipeline, we obtained sequence assemblies of individual BACs that consist, on average, of eight sequence scaffolds and represent >98% of the genomic inserts. Our BAC assemblies are clearly superior to a whole‐genome shotgun assembly regarding contiguity, completeness and the representation of the gene space. Our methods may be employed to rapidly obtain high‐quality assemblies of a large number of clones to assemble map‐based reference sequences of plant and animal species with complex genomes by sequencing along a minimum tiling path.     

中国美利奴细毛羊BAC文库的三维PCR筛选

本研究利用中国美利奴细毛羊全基因组BAC文库,构建了可供快速筛选的两级水平的混合池,一级混合池和二级混合池(Primary pools and secondary pools).一级混合池基于每一384-well盘而构建,由盘、行,列三维混合池组成,二级混合池基于整个BAC文库而构建.设计了一种基于PCR技术的快速筛选方法,先筛选二级混合池m再根据结果筛选相应的一级混合池.利用此方法只需一步共66个PCR反应即可从BAC丈库中7.4万个克隆中筛选出1个阳性克隆,或三步100个以内的PCR反应筛选出多个阳性克隆.以绵羊基因组多态性分子标记BF94-1为引物,用一步共66个PCR反应成功筛选到1个阳性克隆373D13.     

Generation of a soybean BAC library, and identification of DNA sequences tightly linked to the Rps1-k disease resistance gene

 A soybean bacterial artificial chromosome (BAC) library, comprising approximately 45 000 clones, was constructed from high-molecular-weight nuclear DNA of cultivar Williams 82, which carries the Rps1-k gene for resistance against Phytophthora sojae. The library is stored in 130 pools with about 350 clones per pool. Completeness of the library was evaluated for 21 random sequences including four markers linked to the Rps1 locus and 16 cDNAs. We identified pools containing BACs for all sequences except for one cDNA. Additionally, when screened for possible contaminating BAC clones carrying chloroplast genes, no sequences homologous to two barley chloroplast genes were found. The estimated average insert size of the BAC clones was about 105 kb. The library comprises about four genome equivalents of soybean DNA. Therefore, this gives a probability of 0.98 of finding a specific sequence from this library. This library should be a useful resource for the positional cloning of Rps1-k, and other soybean genes. We have also evaluated the feasibility of an RFLP-based screening procedure for the isolation of BAC clones specific for markers that are members of repetitive sequence families, and are linked to the Rps1-k gene. We show that BAC clones isolated for two genetically linked marker loci, Tgmr and TC1-2, are physically linked. Application of this method in expediting the map-based cloning of a gene, especially from an organism, such as soybean, maize and wheat, with a complex genome is discussed. Received: 12 May 1998/Accepted: 24 August 1998     

Construction and characterization of a gridded cattle BAC library

A bovine genomic large-insert bacterial artificial chromosome (BAC) library has been constructed from leukocytes of a Holstein-Friesian male. Size fractionated DpnII-digested genomic DNA was ligated to the dephosphorylated BamH1 ends of a pBACe3.6 vector. Approximately 8.3 x 10(4) individual BAC clones were picked into 384-well plates. Two-hundred and sixty-seven randomly chosen clones were characterized by pulsed-field gel electrophoresis (PFGE). The average insert size was 104 kb with a frequency of clones without inserts of 5.5%. Thirty-four BAC clones were mapped by fluorescence in situ hybridization (FISH) to cattle chromosomes. Three showed signals at more than one location, one of them on the centromeric regions of all autosomes, indicating that the clone contains centromeric repeats. A subset of these BAC clones was used for the development of sequence tagged sites. Both subcloning and direct sequencing of the BACs were used for generating sequence tagged site information. The clones from the library were gridded onto high-density membranes, and PCR superpools were produced from the same set of clones. Membranes and superpools are available through the Resource Centre of the German Human Genome Project in Berlin (http:// www.rzpd.de).     

A second-generation integrated map of the silkworm reveals synteny and conserved gene order between lepidopteran insects

A second-generation linkage map was constructed for the silkworm, Bombyx mori, focusing on mapping Bombyx sequences appearing in public nucleotide databases and bacterial artificial chromosome (BAC) contigs. A total of 874 BAC contigs containing 5067 clones (22% of the library) were constructed by PCR-based screening with sequence-tagged sites (STSs) derived from whole-genome shotgun (WGS) sequences. A total of 523 BAC contigs, including 342 independent genes registered in public databases and 85 expressed sequence tags (ESTs), were placed onto the linkage map. We found significant synteny and conserved gene order between B. mori and a nymphalid butterfly, Heliconius melpomene, in four linkage groups (LGs), strongly suggesting that using B. mori as a reference for comparative genomics in Lepidotera is highly feasible.     

Construction of a bacterial artificial chromosome library for the Rongchang pig breed and its use for the identification of genes involved in intramuscular fat deposition

In a search for genes affecting intramuscular fat deposition, we constructed a bacterial artificial chromosome (BAC) library for the whole genome of Rongchang pig, a domestic Chinese swine breed. The library consisted of approximately 192,000 clones, with an averaged insert size of 116 kb. Frequency of non-insert clone of the BAC library was no higher than 1.8%, based on estimation of 220 BAC clones randomly selected. We estimated the coverage of the library to be more than seven porcine genome equivalents. Subsequent screening of the BAC library with a three-step PCR procedure resulted in identification of seven candidate genes that were potentially involved in intramuscular fat deposition. The number of positive BAC clones ranged from 2 to 4 for each of the seven genes. One positive clone, containing the lipin1 gene, was fully sequenced by shotgun method to generate 118,041 bp porcine genomic sequences. The BAC clone contained complete DNA sequence of porcine lipin1 gene including all the exons and introns. Our results indicate that this BAC library is a useful tool for gene identification and help to serve as an important resource for future porcine genomic study.     

Construction of Japanese BAC library Yamato-2 (JY2): a set of 330K clone resources of damage-minimized DNA taken from a genetically established Japanese individual

A bacterial artificial chromosome (BAC) library referred to as Yamato-2 (JY2), was constructed from a Japanese individual and contained 330,000 clones. Library construction was based on 2 concepts: Japanese pedigree and non-immortalization. Genomic DNA was extracted from white blood cells from umbilical cord blood of a Japanese male individual. Four traits of the sample, (1) amelogenin DNA, (2) short tandem repeat (STR), (3) mitochondrial DNA (mtDNA), and (4) HLA-allele typing, were investigated to verify attribution of the donor. One of the samples with quite good Japanese characteristics was named JY2 and used as a resource for construction of a BAC library. Amelogenin DNA indicated male. STR indicated Mongoloid. MtDNA suggested haplogroup B, which is different from any other diploid whose sequence has been reported. The HLA gene was classified into east-Asian specific haplotype. These results revealed that JY2 was obtained from a Japanese male. We sequenced both ends of 185,012 BAC clones. By using the BLAST search, BAC end sequences (BESs) were mapped on the human reference sequence provided by NCBI. Inserts of individual BAC clones were mapped with both ends properly placed. As a result, 103,647 BAC clones were successfully mapped. The average insert size of BAC calculated from the mapping information was 130?kb. Coverage and redundancy of the reference sequence by successfully mapped BAC clones were 96.4% and 3.9-fold, respectively. This library will be especially suitable as a Japanese standard genome resource. The availability of an accurate library is indispensable for diagnostics or drug-design based on genome information, and JY2 will provide an accurate sequence of the Japanese genome as an important addition to the human genome.     

PCR筛选BAC文库和直接BAC末端测序方法的建立

建立了一种用PCR方法筛选富含高度重复序列的大麦BAC DNA 文库和直接对 BAC DNA进行末端测序的方法.用PCR技术进行大麦BAC DNA 文库(含816个平板,每个平板含384个克隆)的筛选分4步进行.在实验中,建立了两个水平的BAC DNA池(一级池和二级池).一个二级池由一个平板(含有384个克隆)的DNA 组成,一个一级池由连续10个稀释100倍的二级池的DNA混合而成(如1～10,11～20等),共82个一级池.BAC DNA 文库筛选的第一步是对82个一级池的筛选.得到阳性一级池后(如2号一级池),对其所含的10个二级池(从11～20)进行第二步筛选.得到阳性二级池后,培养相应的阳性平板的所有克隆(384个),从头开始(左上侧),每相邻的4个克隆为一组,在96孔板上(4 X 96=384) 进行第三轮PCR反应;之后对筛选结果为阳性的4个克隆分别进行菌落 PCR(第四轮)得到单一阳性克隆.根据BAC DNA Hind III 酶切指纹图谱,对同一引物筛选的BACs进行重叠群作图(Contig).对代表contig 的两端的BAC DNA直接进行末端测序并对测序结果Blast,以检测其在大麦中是否属于单拷贝序列.根据测序和Blast结果设计引物,用中国春附加系(附加大麦染色体)对来自BAC克隆的引物进行染色体定位并用分离群体进行遗传学作图,以确定是否可以用作下一步的染色体步行.     

A BAC-based integrated linkage map of the silkworm Bombyx mori

Background

In 2004, draft sequences of the model lepidopteran Bombyx mori were reported using whole-genome shotgun sequencing. Because of relatively shallow genome coverage, the silkworm genome remains fragmented, hampering annotation and comparative genome studies. For a more complete genome analysis, we developed extended scaffolds combining physical maps with improved genetic maps.

Results

We mapped 1,755 single nucleotide polymorphism (SNP) markers from bacterial artificial chromosome (BAC) end sequences onto 28 linkage groups using a recombining male backcross population, yielding an average inter-SNP distance of 0.81 cM (about 270 kilobases). We constructed 6,221 contigs by fingerprinting clones from three BAC libraries digested with different restriction enzymes, and assigned a total of 724 single copy genes to them by BLAST (basic local alignment search tool) search of the BAC end sequences and high-density BAC filter hybridization using expressed sequence tags as probes. We assigned 964 additional expressed sequence tags to linkage groups by restriction fragment length polymorphism analysis of a nonrecombining female backcross population. Altogether, 361.1 megabases of BAC contigs and singletons were integrated with a map containing 1,688 independent genes. A test of synteny using Oxford grid analysis with more than 500 silkworm genes revealed six versus 20 silkworm linkage groups containing eight or more orthologs of Apis versus Tribolium, respectively.

Conclusion

The integrated map contains approximately 10% of predicted silkworm genes and has an estimated 76% genome coverage by BACs. This provides a new resource for improved assembly of whole-genome shotgun data, gene annotation and positional cloning, and will serve as a platform for comparative genomics and gene discovery in Lepidoptera and other insects.     

Construction of a physical and transcript map for a 1-Mb genomic region containing the urofacial (Ochoa) syndrome gene on 10q23-q24 and localization of the disease gene within two overlapping BAC clones (<360 kb).

Urofacial (Ochoa) syndrome is an autosomal recessive disease characterized by distorted facial expression and urinary abnormalities. Previously, we mapped the UFS gene to chromosome 10q23-q24 and narrowed the interval to one YAC clone of 1410 kb. Here, we have constructed a BAC/PAC contig of the 1-Mb region using STS content mapping with 42 BAC/PAC-end sequences, 9 previously reported and 16 newly identified microsatellite markers, and 14 EST markers. A total of 26 polymorphic microsatellite markers were genotyped for 31 UFS patients from Colombia and 2 patients from the United States. Haplotype analyses suggest that the UFS gene is located within two overlapping BAC clones, a region of <360 kb of DNA sequence. We tested 42 EST markers previously mapped to the D10S1709-D10S603 interval against the BAC/PAC contig and identified 11 ESTs located in the 1-Mb region. Four of the 11 ESTs mapped to the 360-kb UFS critical region. Shotgun sequencing of the two BAC clones and BLASTN search of the EST databases revealed 3 other ESTs contained in the UFS critical region. These results will facilitate the cloning and identification of the UFS gene.     

水稻中一个NBS-LRR抗病同源基因家族的克隆和分析

利用克隆的抗病基因同源序列RS13作为探针,从水稻IR64的BAC文库中筛选到4个阳性克隆,其中一个克隆14E19能够覆盖其余3个克隆。对14E19进行全序列测定和分析,获得了73kb的全长DNA序列,基因预测显示其上有4个编码NBS-LRR结构域的基因(NL),分别命名为NL-A,B,C和D。对具有相同基因组背景的IRBB56同一染色体位置上跨度更大的BAC克隆106P13进行分析,发现其上有10个NL同源拷贝,其中4个同14E19上的NL一样。搜索日本晴、93—11、广陆矮4号的序列,发现三者有类似的同源序列。但与已知的NBS-LRR抗病基因同源性较低,说明NL是一个至少由10个成员(分别命名为NL-A至J)组成的新基因家族。对NL家族进行RT-PCR和cDNA库筛选分析,发现NL-B基因能够在抗白叶枯病品系IRBB4中表达,暗示该基因参与了抗病反应。