Construction of random sheared fosmid library from Chinese cabbage and its use for Brassica rapa genome sequencing project

As a part of the Multinational Genome Sequencing Project of Brassica rapa, linkage group R9 and R3 were sequenced using a bacterial artificial chromosome (BAC) by BAC strategy. The current physical contigs are expected to cover approximately 90% euchromatins of both chromosomes. As the project progresses, BAC selection for sequence extension becomes more limited because BAC libraries are restriction enzyme-specific. To support the project, a random sheared fosmid library was constructed. The library consists of 97536 clones with average insert size of approximately 40 kb corresponding to seven genome equivalents, assuming a Chinese cabbage genome size of 550 Mb. The library was screened with primers designed at the end of sequences of nine points of scaffold gaps where BAC clones cannot be selected to extend the physical contigs. The selected positive clones were end-sequenced to check the overlap between the fosmid clones and the adjacent BAC clones.Nine fosmid clones were selected and fully sequenced. The sequences revealed two completed gap filling and seven sequence extensions, which can be used for further selection of BAC clones confirming that the fosmid library will facilitate the sequence completion of B. rapa.     

Construction and utility of 10-kb libraries for efficient clone-gap closure for rice genome sequencing

Rice is an important crop and a model system for monocot genomics, and is a target for whole genome sequencing by the International Rice Genome Sequencing Project (IRGSP). The IRGSP is using a clone by clone approach to sequence rice based on minimum tiles of BAC or PAC clones. For chromosomes 10 and 3 we are using an integrated physical map based on two fingerprinted and end-sequenced BAC libraries to identifying a minimum tiling path of clones. In this study we constructed and tested two rice genomic libraries with an average insert size of 10 kb (10-kb library) to support the gap closure and finishing phases of the rice genome sequencing project. The HaeIII library contains 166,752 clones covering approximately 4.6x rice genome equivalents with an average insert size of 10.5 kb. The Sau3AI library contains 138,960 clones covering 4.2x genome equivalents with an average insert size of 11.6 kb. Both libraries were gridded in duplicate onto 11 high-density filters in a 5 x 5 pattern to facilitate screening by hybridization. The libraries contain an unbiased coverage of the rice genome with less than 5% contamination by clones containing organelle DNA or no insert. An efficient method was developed, consisting of pooled overgo hybridization, the selection of 10-kb gap spanning clones using end sequences, transposon sequencing and utilization of in silico draft sequence, to close relatively small gaps between sequenced BAC clones. Using this method we were able to close a majority of the gaps (up to approximately 50 kb) identified during the finishing phase of chromosome-10 sequencing. This method represents a useful way to close clone gaps and thus to complete the entire rice genome.     

The first insight into the <Emphasis Type="Italic">Taxus</Emphasis> genome via fosmid library construction and end sequencing

Taxus mairei is a critically endangered and commercially important cultured medicinal gymnosperm in China and forms an important medicinal resource, but the research of its genome is absent. In this study, we constructed a T. mairei fosmid library and analyzed the fosmid end sequences to provide a preliminary assessment of the genome. The library consists of one million clones with an average insert size of about 39 kb, amounting to 3.9 genome equivalents. Fosmid stability assays indicate that T. mairei DNA was stable during propagation in the fosmid system. End sequencing of both 5′ and 3′ ends of 968 individual clones generated 1,923 sequences after trimming, with an average sequence length of 839 bp. BLASTN searches of the nr and EST databases of GenBank and BLASTX searches of the nr database resulted in 560 (29.1%) significant hits (E < e⁻⁵). Repetitive sequences analysis revealed that 20.8% of end sequences are repetitive elements, which were composed of retroelements, DNA transposons, satellites, simple repeats, and low complexity sequences. The distribution pattern of various repeat types was found to be more similar to the gymnosperm Pinus and Picea than to the monocot and dicot. The satellites of T. mairei were significantly longer than those of P. taeda and P. glauca. The tetra-nucleotide repeats of T. mairei were much longer than those of P. glauca and P. taeda. The fosmid library and the fosmid end sequences, for the first time, will serve as a useful resource for large-scale genome sequencing, physical mapping, SSR marker development and positional cloning, and provide a better understanding of the Taxus genome.     

A fosmid-based genomic map and identification of 474 genes of the hyperthermophilic archaeonPyrobaculum aerophilum

We have constructed a physical map of the approximately 1.7-Mb genome of the hyperthermophilic archaeon Pyrobaculum aerophilum. Derived from a 12× coverage genomic fosmid library with an average insert size of 36 Kb, the map consists of a single circular contig of 96 overlapping fosmid clones with 211 markers ordered along them. One hundred of the sequence markers have strong similarities to known genes. Many overlaps were also checked using restriction fingerprint analysis. This map is an important step in the elucidation of the sequence of the entire genome of Pyrobaculum aerophilum. To this end we have determined more than 95% of the genome with 15000 random sequences. Each sequence has been screened against the public sequence databases to identify similarities to known genes. We report here a list of the 474 putative genes we have identified. Received: 18 October 1996 / Accepted: 30 October 1996     

Construction and characterization of a copy number-inducible fosmid library of Xanthomonas oryzae pathovar oryzae MAFF311018

A fosmid library of Xanthomonas oryzae pathovar oryzae MAFF311018 (T7174), the causative agent of bacterial blight on rice, was constructed and characterized. The average fosmid library insert size was > 34 kb, and 967 clones were uniquely positioned on its sequenced genome. The entire Xoo MAFF311018 genome was covered by end-sequenced clones with at least 5 kb of overlap. The fosmid vector contains both the single-copy Escherichia coli fertility factor origin, which enhances fosmid stability, and the multi-copy IncPα origin, allowing amplification of copy number upon induction with l-arabinose. Real-time quantitative PCR on 12 randomly picked fosmid library clones determined that fosmid copy number increased 8- to 58-fold after 5 hour induction. This library provides a new resource for complementation experiments and systematic functional studies in Xoo and related species.     

Construction, characterization, and chromosomal mapping of a fosmid library of the white-cheeked gibbon (Nomascus leucogenys)

Gibbons have experienced extensive karyotype rearrangements during evolution and represent an ideal model for studying the underlying molecular mechanism of evolutionary chromosomal rearrangements. It is anticipated that the cloning and sequence characterization of evolutionary chromosomal breakpoints will provide vital insights into the molecular force that has driven such a radical karyotype reshuffle in gibbons. We constructed and characterized a high-quality fosmid library of the white-cheeked gibbon (Nomascus leucogenys) containing 192,000 non- redundant clones with an average insert size of 38 kb and 2.5-fold genome coverage. By end sequencing of 100 randomly selected fosmid clones, we generated 196 sequence tags for the library. These end-sequenced fosmid clones were then mapped onto the chromosomes of the white-cheeked gibbon by fluorescence in situ hybridization, and no spurious chimeric clone was detected. BLAST search against the human genome showed a good correlation between the number of hit clones and the number of chromosomes, an indication of unbiased chromosomal distribution of the fosmid library. The chromosomal distribution of the mapped clones is also consistent with the BLAST search result against human and white-cheeked gibbon genomes. The fosmid library and the mapped clones will serve as a valuable resource for further studying gibbons' chromosomal rearrangements and the underlying molecular mechanism as well as for comparative genomic study in the lesser apes.     

Using partial genomic fosmid libraries for sequencing complete organellar genomes

Organellar genome sequences provide numerous phylogenetic markers and yield insight into organellar function and molecular evolution. These genomes are much smaller in size than their nuclear counterparts; thus, their complete sequencing is much less expensive than total nuclear genome sequencing, making broader phylogenetic sampling feasible. However; for some organisms, it is challenging to isolate plastid DNA for sequencing using standard methods. To overcome these difficulties, we constructed partial genomic libraries from total DNA preparations of two heterotrophic and two autotrophic angiosperm species using fosmid vectors. We then used macroarray screening to isolate clones containing large fragments of plastid DNA. A minimum tiling path of clones comprising the entire genome sequence of each plastid was selected, and these clones were shotgun-sequenced and assembled into complete genomes. Although this method worked well for both heterotrophic and autotrophic plants, nuclear genome size had a dramatic effect on the proportion of screened clones containing plastid DNA and, consequently, the overall number of clones that must be screened to ensure full plastid genome coverage. This technique makes it possible to determine complete plastid genome sequences for organisms that defy other available organellar genome sequencing methods, especially those for which limited amounts of tissue are available.     

Construction of a gorilla fosmid library and its PCR screening system

A gorilla fosmid library of 261,120 independent clones was constructed and characterized. The fosmid vector is similar to the cosmid in average insert size of ca. 40 kb but contains the F factor for replication, and it is more resistant to recombination. This clone library represents about 3.7 times coverage of the gorilla genome. A simple screening system by PCR was established, and we successfully found 9 clones that cover the entire Hox A gene cluster of the gorilla genome. This gorilla fosmid DNA library is a useful resource for comparative genomics of human and apes.     

Genes and genomes: Towards construction of an overlapping YAC library of the Arabidopsis thaliana genome

Arabidopsis thaliana (Thale cress, Arabidopsis) is an ideal model organism for the molecular genetic analysis of many plant processes. The availability of a complete physical map would greatly facilitate the gene cloning steps in these studies. The small genome size of Arabidopsis makes the construction of such a map a feasible goal. One of the approaches to construct an overlapping library of the Arabidopsis genome takes advantage of the many mapped markers and the availability of Arabidopsis yeast artificial chromosome (YAC) libraries. Mapped molecular markers are used to identify corresponding YAC clones and thereby place them on the genetic map. Subsequently, these YAC clones provide the framework for directed walking experiments aimed at closing the gaps between the YAC contigs. Adopting this strategy, YAC clones comprising about 10% of the genome have been assigned to the top halves of Arabidopsis chromosomes 4 and 5. Extensive walking experiments in a 10 cM interval of chromosome 4 have resulted in two contiguous regions in the megabase size range.     

A set of BAC clones spanning the human genome

Using the human bacterial artificial chromosome (BAC) fingerprint-based physical map, genome sequence assembly and BAC end sequences, we have generated a fingerprint-validated set of 32855 BAC clones spanning the human genome. The clone set provides coverage for at least 98% of the human fingerprint map, 99% of the current assembled sequence and has an effective resolving power of 79 kb. We have made the clone set publicly available, anticipating that it will generally facilitate FISH or array-CGH-based identification and characterization of chromosomal alterations relevant to disease.     

Use of the IGF BAC library for physical mapping of the Arabidopsis thaliana genome

In order to generate a physical map of the Arabidopsis thaliana genome based on bacterial artificial chromosome clones (BACs), an iterative high throughput hybridisation strategy was applied and its efficiency was evaluated. Thus, probes generated from both ends of 500 BAC clones selected from the Arabidopsis –IGF–BAC library were hybridised to the entire library gridded on high density filters. The 1000 hybridisation reactions identified 4496 clones (41.8% of the complete library, or 50.3% if organellar, centromeric, and ribosomal DNA carrying clones are excluded) which were assembled into a minimum of 220 contigs. These results demonstrate the viability of the applied ‘double-end clone-limited/sampling without replacement’ hybridisation strategy for the generation of a high resolution physical map, and provide a highly useful resource for map-based gene cloning approaches and further genome analysis.     

Fosmid Library Construction and Initial Analysis of End Sequences in Zhikong Scallop (<Emphasis Type="Italic">Chlamys farreri</Emphasis>)

Zhikong scallop (Chlamys farreri Jones et Preston, 1904) is one of the most commercially important bivalves in China, but research on its genome is underdeveloped. In this study, we constructed the first Zhikong scallop fosmid library, and analyzed the fosmid end sequences to provide a preliminary assessment of the genome. The library consists of 133,851 clones with an average insert size of about 40 kb, amounting to 4.3 genome equivalents. Fosmid stability assays indicate that Zhikong scallop DNA was stable during propagation in the fosmid system. Library screening with two genes and seven microsatellite markers yielded between two and eight positive clones, and none of those tested was absent from the library. End-sequencing of 480 individual clones generated 828 sequences after trimming, with an average sequence length of 624 bp. BLASTN searches of the nr and EST databases of GenBank and BLASTX searches of the nr database resulted in 213 (25.72%) and 44 (5.31%) significant hits (E < e⁻⁵), respectively. Repetitive sequences analysis resulted in 375 repeats, accounting for 15.84% of total length, which were composed of interspersed repetitive sequences, tandem repeats, and low-complexity sequences. The fosmid library, in conjunction with the fosmid end sequences, will serve as a useful resource for physical mapping and positional cloning, and provide a better understanding of the Zhikong scallop genome.     

A fine physical map of the rice chromosome 5

A fine physical map of the rice (Oryza sativa spp. Japonica var. Nipponbare) chromosome 5 with bacterial artificial chromosome (BAC) and PI-derived artificial chromosome (PAC) clones was constructed through integration of 280 sequenced BAC/PAC clones and 232 sequence tagged site/expressed sequence tag markers with the use of fingerprinted contig data of the Nipponbare genome. This map consists of five contigs covering 99% of the estimated chromosome size (30.08 Mb). The four physical gaps were estimated at 30 and 20 kb for gaps 1–3 and gap 4, respectively. We have submitted 42.2-Mb sequences with 29.8 Mb of nonoverlapping sequences to public databases. BAC clones corresponding to telomere and centromere regions were confirmed by BAC-fluorescence in situ hybridization (FISH) on a pachytene chromosome. The genetically centromeric region at 54.6 cM was covered by a minimum tiling path spanning 2.1 Mb with no physical gaps. The precise position of the centromere was revealed by using three overlapping BAC/PACs for ~150 kb. In addition, FISH results revealed uneven chromatin condensation around the centromeric region at the pachytene stage. This map is of use for positional cloning and further characterization of the rice functional genomics. Electronic supplementary material Supplementary material is available in the online version of this article at and is accessible for authorized users. Chia-Hsiung Cheng and Mei-Chu Chung have equal contributions.     

Fosmid Library Construction and Initial Analysis of End Sequences in Female Half-smooth Tongue Sole (Cynoglossus semilaevis)

Half-smooth tongue sole (Cynoglossus semilaevis: Pleuronectiformes) is a commercially important cultured marine flatfish in China and forms an important fishery resource, but the research of its genome is underdeveloped. In this study, we constructed a female C. semilaevis fosmid library and analyzed the fosmid end sequences to provide a preliminary assessment of the genome. The library consists of 49,920 clones with an average insert size of about 39 kb, amounting to 3.23 genome equivalents. Fosmid stability assays indicate that female C. semilaevis DNA was stable during propagation in the fosmid system. Library screening with eight microsatellite markers yielded between two and five positive clones, and none of those tested was absent from the library. End-sequencing of both 5′ and 3′ ends of 1,152 individual clones generated 2,247 sequences after trimming, with an average sequence length of 855 bp. BLASTN searches of the nr and EST databases of GenBank and BLASTX searches of the nr database resulted in 259 (11.53%) and 287 (12.77%) significant hits (E < e ⁻⁵), respectively. Repetitive sequences analysis resulted in 5.23% of base pairs masked using both the Fugu and Danio databases, repetitive elements were composed of retroelements, DNA transposons, satellites, simple repeats, and low-complexity sequences. The fosmid library, in conjunction with the fosmid end sequences, will serve as a useful resource for large-scale genome sequencing, physical mapping, and positional cloning, and provide a better understanding of female C. semilaevis genome.     

Finishing a whole-genome shotgun: release 3 of the Drosophila melanogaster euchromatic genome sequence

Background

The Drosophila melanogaster genome was the first metazoan genome to have been sequenced by the whole-genome shotgun (WGS) method. Two issues relating to this achievement were widely debated in the genomics community: how correct is the sequence with respect to base-pair (bp) accuracy and frequency of assembly errors? And, how difficult is it to bring a WGS sequence to the accepted standard for finished sequence? We are now in a position to answer these questions.

Results

Our finishing process was designed to close gaps, improve sequence quality and validate the assembly. Sequence traces derived from the WGS and draft sequencing of individual bacterial artificial chromosomes (BACs) were assembled into BAC-sized segments. These segments were brought to high quality, and then joined to constitute the sequence of each chromosome arm. Overall assembly was verified by comparison to a physical map of fingerprinted BAC clones. In the current version of the 116.9 Mb euchromatic genome, called Release 3, the six euchromatic chromosome arms are represented by 13 scaffolds with a total of 37 sequence gaps. We compared Release 3 to Release 2; in autosomal regions of unique sequence, the error rate of Release 2 was one in 20,000 bp.

Conclusions

The WGS strategy can efficiently produce a high-quality sequence of a metazoan genome while generating the reagents required for sequence finishing. However, the initial method of repeat assembly was flawed. The sequence we report here, Release 3, is a reliable resource for molecular genetic experimentation and computational analysis.     

An integrated physical and genetic map of the rice genome

Rice was chosen as a model organism for genome sequencing because of its economic importance, small genome size, and syntenic relationship with other cereal species. We have constructed a bacterial artificial chromosome fingerprint–based physical map of the rice genome to facilitate the whole-genome sequencing of rice. Most of the rice genome (~90.6%) was anchored genetically by overgo hybridization, DNA gel blot hybridization, and in silico anchoring. Genome sequencing data also were integrated into the rice physical map. Comparison of the genetic and physical maps reveals that recombination is suppressed severely in centromeric regions as well as on the short arms of chromosomes 4 and 10. This integrated high-resolution physical map of the rice genome will greatly facilitate whole-genome sequencing by helping to identify a minimum tiling path of clones to sequence. Furthermore, the physical map will aid map-based cloning of agronomically important genes and will provide an important tool for the comparative analysis of grass genomes.     

Bacterial artificial chromosome-based physical map of Gibberella zeae (Fusarium graminearum).

Fusarium graminearum is the primary causal pathogen of Fusarium head blight of wheat and barley. To accelerate genomic analysis of F. graminearum, we developed a bacterial artificial chromosome (BAC)-based physical map and integrated it with the genome sequence and genetic map. One BAC library, developed in the HindIII restriction enzyme site, consists of 4608 clones with an insert size of approximately 107 kb and covers about 13.5 genome equivalents. The other library, developed in the BamHI restriction enzyme site, consists of 3072 clones with an insert size of approximately 95 kb and covers about 8.0 genome equivalents. We fingerprinted 2688 clones from the HindIII library and 1536 clones from the BamHI library and developed a physical map of F. graminearum consisting of 26 contigs covering 39.2 Mb. Comparison of our map with the F. graminearum genome sequence showed that the size of our physical map is equivalent to the 36.1 Mb of the genome sequence. We used 31 sequence-based genetic markers, randomly spaced throughout the genome, to integrate the physical map with the genetic map. We also end-sequenced 17 BamHI BAC clones and identified 4 clones that spanned gaps in the genome sequence. Our new integrated map is highly reliable and useful for a variety of genomics studies.     

Physical mapping of the rice genome with YAC clones

Construction of a rice physical map covered by YAC clones which have been arranged over half of the genome length is presented here. A total of 1285 RFLP and RAPD markers almost evenly distributed on the rice genetic map could select 2974 YAC clones and 2443 clones of them were located on their original positions. Rice YACs carrying 350 kb average insert fragments of 2443 clones could cover 222 megabase length of the rice genome, corresponding to 52% of the whole genome size (4.3 Mb). Chromosome landing with many YAC clones on the high-density genetic map loci efficiently integrated the genetic map with a physical map. This is the first step to generate a comprehensive genome map of rice. An integrated genome map should be an indispensable tool to figure out genome structure as well as to clone trait genes by map-based cloning.     

Comparative physical mapping of targeted regions of the rat genome

The comparative mapping and sequencing of vertebrate genomes is now a key priority for the Human Genome Project. In addition to finishing the human genome sequence and generating a `working draft' of the mouse genome sequence, significant attention is rapidly turning to the analysis of other model organisms, such as the laboratory rat (Rattus norvegicus). As a complement to genome-wide mapping and sequencing efforts, it is often important to generate detailed maps and sequence data for specific regions of interest. Using an adaptation of our previously described approach for constructing mouse comparative and physical maps, we have established a general strategy for targeted mapping of the rat genome. Specifically, we constructed a framework comparative map of human Chromosome (Chr) 7 and the orthologous regions of the rat genome, as well as two large (>1-Mb) P1-derived artificial chromosome (PAC)-based physical maps. Generation of these physical maps involved the use of mouse-derived probes that cross-hybridized with rat PAC clones. The first PAC map encompasses the cystic fibrosis transmembrane conductance regulator gene (Cftr), while the second map allows a three-species comparison of a genomic region containing intra- and inter-chromosomal evolutionary rearrangements. The studies reported here further demonstrate that cross-species hybridization between related animals, such as rat and mouse, can be readily used for the targeted construction of clone-based physical maps, thereby accelerating the analysis of biologically interesting regions of vertebrate genomes. Received: 5 December 2000 / Accepted: 27 February 2001     

One large-insert plant-transformation-competent BIBAC library and three BAC libraries of Japonica rice for genome research in rice and other grasses

We report one large-insert BIBAC library and three BAC libraries for japonica rice cv Nipponbare. The BIBAC library was constructed in the HindIII site of a plant-transformation-competent binary vector (pCLD04541) and the three BAC libraries were constructed in the BamHI, HindIII and EcoRI sites of a BAC vector (pECBAC1), respectively. Each library contains 23,040 clones, has an average insert size of 130 kb, 170 kb, 150 kb and 156 kb, and covers 6.7x, 8.7x, 7.7x and 8.0 x rice haploid genomes, respectively. The combined libraries contain 92,160 clones in total, covering 31.1 x rice haploid genomes. To demonstrate their utility, we screened the libraries with 55 DNA markers mapped to chromosome 8 of the rice genetic maps and analyzed a number of clones by the restriction fingerprinting and contig assembly method. The results indicate that the libraries completely cover the rice genome and, thus, are well-suited for genome research in rice and other gramineous crops. The BIBAC library represents the first plant-transformation-competent large-insert DNA library for rice, which will streamline map-based cloning, functional analysis of the rice genome sequence and molecular breeding in rice and other grass species. These libraries are being used in the development of a whole-genome, BAC/BIBAC-based, integrated physical, genetic and sequence map of rice and in the research of genome-wide comparative genomics of grass species.