Construction and partial characterization of a Corynebacterium pseudotuberculosis bacterial artificial chromosome library through genomic survey sequencing.

Corynebacterium pseudotuberculosis is a gram-positive bacterium that causes caseous lymphadenitis in sheep and goats. However, despite the economic losses caused by caseous lymphadenitis, there is little information about the molecular mechanisms of pathogenesis of this bacterium. Genomic libraries constructed in bacterial artificial chromosome (BAC) vectors have become the method of choice for clone development in high-throughput genomic-sequencing projects. Large-insert DNA libraries are useful for isolation and characterization of important genomic regions and genes. In order to identify targets that might be useful for genome sequencing, we constructed a C. pseudotuberculosis BAC library in the vector pBeloBAC11. This library contains about 18,000 BAC clones, with inserts ranging in size from 25 to 120 kb, theoretically representing a 390-fold coverage of the C. pseudotuberculosis genome (estimated to be 2.5-3.1 Mb). Many genomic survey sequences (GSSs) with homology to C. diphtheriae, C. glutamicum, C. efficiens, and C. jeikeium proteins were observed within a sample of 215 sequenced clones, confirming their close phylogenetic relationship. Computer analyses of GSSs did not detect chimeric, deleted, or rearranged BAC clones, showing that this library has low redundancy. This GSSs collection is now available for further genetic and physical analysis of the C. pseudotuberculosis genome. The GSS strategy that we used to develop our library proved to be efficient for the identification of genes and will be an important tool for mapping, assembly, comparative, and functional genomic studies in a C. pseudotuberculosis genome sequencing project that will begin this year.

Construction and partial characterization of a Corynebacterium pseudotuberculosis bacterial artificial chromosome library through genomic survey sequencing

Corynebacterium pseudotuberculosis is a gram-positive bacterium that causes caseous lymphadenitis in sheep and goats. However, despite the economic losses caused by caseous lymphadenitis, there is little information about the molecular mechanisms of pathogenesis of this bacterium. Genomic libraries constructed in bacterial artificial chromosome (BAC) vectors have become the method of choice for clone development in high-throughput genomic-sequencing projects. Large-insert DNA libraries are useful for isolation and characterization of important genomic regions and genes. In order to identify targets that might be useful for genome sequencing, we constructed a C. pseudotuberculosis BAC library in the vector pBeloBAC11. This library contains about 18,000 BAC clones, with inserts ranging in size from 25 to 120 kb, theoretically representing a 390-fold coverage of the C. pseudotuberculosis genome (estimated to be 2.5-3.1 Mb). Many genomic survey sequences (GSSs) with homology to C. diphtheriae, C. glutamicum, C. efficiens, and C. jeikeium proteins were observed within a sample of 215 sequenced clones, confirming their close phylogenetic relationship. Computer analyses of GSSs did not detect chimeric, deleted, or rearranged BAC clones, showing that this library has low redundancy. This GSSs collection is now available for further genetic and physical analysis of the C. pseudotuberculosis genome. The GSS strategy that we used to develop our library proved to be efficient for the identification of genes and will be an important tool for mapping, assembly, comparative, and functional genomic studies in a C. pseudotuberculosis genome sequencing project that will begin this year.

Genome sequence of the plant pathogen Ralstonia solanacearum.

Ralstonia solanacearum is a devastating, soil-borne plant pathogen with a global distribution and an unusually wide host range. It is a model system for the dissection of molecular determinants governing pathogenicity. We present here the complete genome sequence and its analysis of strain GMI1000. The 5.8-megabase (Mb) genome is organized into two replicons: a 3.7-Mb chromosome and a 2.1-Mb megaplasmid. Both replicons have a mosaic structure providing evidence for the acquisition of genes through horizontal gene transfer. Regions containing genetically mobile elements associated with the percentage of G+C bias may have an important function in genome evolution. The genome encodes many proteins potentially associated with a role in pathogenicity. In particular, many putative attachment factors were identified. The complete repertoire of type III secreted effector proteins can be studied. Over 40 candidates were identified. Comparison with other genomes suggests that bacterial plant pathogens and animal pathogens harbour distinct arrays of specialized type III-dependent effectors.

A BAC library and paired-PCR approach to mapping and completing the genome sequence of Sulfolobus solfataricus P2.

The original strategy used in the Sulfolobus solfataricus genome project was to sequence non overlapping, or minimally overlapping, cosmid or lambda inserts without constructing a physical map. However, after only about two thirds of the genome sequence was completed, this approach became counter-productive because there was a high sequence bias in the cosmid and lambda libraries. Therefore, a new approach was devised for linking the sequenced regions which may be generally applicable. BAC libraries were constructed and terminal sequences of the clones were determined and used for both end mapping and PCR screening. The PCR approaches included a novel chromosome walking method termed "paired-PCR". 21 gaps were filled by BAC end sequence analyses and 6 gaps were filled by PCR including three large ones by paired-PCR. The complete map revealed that 0.9 Mb remained to be sequenced and 34 BAC clones were selected for walking over small gaps and preparing template libraries for larger ones. It is concluded that an optimal strategy for sequencing microorganism genomes involves construction of a high-resolution physical map by BAC end analyses, PCR screening and paired-PCR chromosome walking after about half the genome sequence has been accumulated.

Characterization and repeat analysis of the compact genome of the freshwater pufferfish Tetraodon nigroviridis.

Tetraodon nigroviridis is a freshwater pufferfish 20-30 million years distant from Fugu rubripes. The genome of both tetraodontiforms is compact, mostly because intergenic and intronic sequences are reduced in size compared to other vertebrate genomes. The previously uncharacterized Tetraodon genome is described here together with a detailed analysis of its repeat content and organization. We report the sequencing of 46 megabases of bacterial artificial chromosome (BAC) end sequences, which represents a random DNA sample equivalent to 13% of the genome. The sequence and location of rRNA gene clusters, centromeric and subtelocentric satellite sequences have been determined. Minisatellites and microsatellites have been cataloged and notable differences were observed in comparison with microsatellites from Fugu. The genome contains homologies to all known families of transposable elements, including Ty3-gypsy, Ty1-copia, Line retrotransposons, DNA transposons, and retroviruses, although their overall abundance is <1%. This structural analysis is an important prerequisite to sequencing the Tetraodon genome.

Construction of a swine BAC library: application to the characterization and mapping of porcine type C endoviral elements.

A porcine bacterial artificial chromosome (BAC) library was constructed using the pBeloBAC11 vector. It comprised 107,520 clones with an average insert size of 135 kb, representing an almost fivefold coverage of the swine haploid genome. Screening of the library allowed recovery of one to eight clones for 142 unique markers located all over the genome, while it failed for only one marker. About 4% chimeric clones were found. The library was also screened for the protease gene of type C porcine endoviral sequences (PERVs), and 62 clones were recovered, all but two of which contained one protease gene. We found 20 protease sequences (PERV-1 to PERV-20) which, despite differing by point mutations, were all coding sequences. The most frequent sequence, PERV-2, was 100% similar to a protease sequence expressed in the porcine PK-15 cell line. Most of the clones harbored envelope genes. Thirty-three BAC clones were mapped by fluorescence in situ hybridization to 22 distinct locations on 14 chromosomes, including the X and Y chromosomes. These overall results indicate that there is generally one PERV copy per integration site. Although PERV sequences were not tandemly arranged, clusters of integration sites were observed at positions 3p1.5 and 7p1.1. Southern blot experiments revealed 20-30 PERV copies in the Large White pig genome studied here, and variations in PERV content among pigs of different breeds were observed. In conclusion, this BAC collection represents a significant contribution to the swine large genomic DNA cloned insert resources and provides the first detailed map of PERV sequences in the swine genome. This work is the first step toward identification of potential active sites of PERV elements.

The 102-kilobase pgm locus of Yersinia pestis: sequence analysis and comparison of selected regions among different Yersinia pestis and Yersinia pseudotuberculosis strains.

We report the complete 119,443-bp sequence of the pgm locus from Yersinia pestis and its flanking regions. Sequence analysis confirms that the 102-kb unstable pgm locus is composed of two distinct parts: the pigmentation segment and a high-pathogenicity island (HPI) which carries virulence genes involved in iron acquisition (yersiniabactin biosynthetic gene cluster). Within the HPI, three genes coding for proteins related to phage proteins were uncovered. They are located at both extremities indicating that the entire HPI was acquired en bloc by phage-mediated horizontal transfer. We identified, within the pigmentation segment, two novel loci that may be involved in virulence: a fimbriae gene cluster and a locus probably encoding a two component regulatory system similar to the BvgAS regulatory system of Bordetella pertussis. Three genes containing frameshift mutations and two genes interrupted by insertion element insertion were found within this region. To investigate diversity among different Y. pestis and Yersinia pseudotuberculosis strains, the sequence of selected regions of the pgm locus and flanking regions were compared from 20 different Y. pestis and 10 Y. pseudotuberculosis strains. The results showed that the genes interrupted in Y. pestis are intact in Y. pseudotuberculosis. However, one of these mutations, in the bvgS homologue, is only present in Y. pestis strains of biovar Orientalis and not in those of the biovars Antiqua and Medievalis. The results obtained by analysis of variable positions in the sequence are in accordance with historical records, confirming that biovar Orientalis is the most recent lineage. Furthermore, sequence comparisons among 29 Yersinia strains suggest that Y. pestis is a recently emerged pathogen that is probably entering the initial phase of reductive evolution.

Construction and characterization of a sheep BAC library of three genome equivalents.

A sheep BAC library of over three genome equivalents was constructed and arrayed in superpools and row, column, and plate pools. The library contains 90,000 clones distributed in 39 superpools. The average insert size was estimated at 123 kb. The library was screened by PCR with 77 primer pairs corresponding to ovine microsatellites distributed throughout the genome. The probability of finding a random sequence in the library could be estimated at 0.96.

Identification of variable regions in the genomes of tubercle bacilli using bacterial artificial chromosome arrays.

Whole-genome comparisons of the tubercle bacilli were undertaken using ordered bacterial artificial chromosome (BAC) libraries of Mycobacterium tuberculosis and the vaccine strain, Mycobacterium bovis BCG-Pasteur, together with the complete genome sequence of M. tuberculosis H37Rv. Restriction-digested BAC arrays of M. tuberculosis H37Rv were used in hybridization experiments with radiolabelled M. bovis BCG genomic DNA to reveal the presence of 10 deletions (RD1-RD10) relative to M. tuberculosis. Seven of these regions, RD4-RD10, were also found to be deleted from M. bovis, with the three M. bovis BCG-specific deletions being identical to the RD1-RD3 loci described previously. The distribution of RD4-RD10 in Mycobacterium africanum resembles that of M. tuberculosis more closely than that of M. bovis, whereas an intermediate arrangement was found in Mycobacterium microti, suggesting that the corresponding genes may affect host range and virulence of the various tubercle bacilli. Among the known products encoded by these loci are a copy of the proposed mycobacterial invasin Mce, three phospholipases, several PE, PPE and ESAT-6 proteins, epoxide hydrolase and an insertion sequence. In a complementary approach, direct comparison of BACs uncovered a third class of deletions consisting of two M. tuberculosis H37Rv loci, RvD1 and RvD2, deleted from the genome relative to M. bovis BCG and M. bovis. These deletions affect a further seven genes, including a fourth phospholipase, plcD. In summary, the insertions and deletions described here have important implications for our understanding of the evolution of the tubercle complex.

Genetic and physical analyses of the centromeric and pericentromeric regions of human chromosome 5: recombination across 5cen.

Human centromeres are poorly understood at both the genetic and the physical level. In this paper, we have been able to distinguish the alphoid centromeric sequences of chromosome 5 from those of chromosome 19. This result was obtained by pulsed-field gel electrophoresis after cutting genomic DNA with restriction endonucleases NcoI (chromosome 5) and BamHI (chromosome 19). We could thus define a highly polymorphic marker, representing length variations of the D5Z1 domain located at the q arm boundary of the chromosome 5 centromere. The centromeric region of chromosome 5 was then analyzed in full detail. We established an approximately 4.6-Mb physical map of the whole region with five rare-cutting enzymes by using nonchimeric YACs, two of which were shown to contain the very ends of 5cen on both sides. The p-arm side of 5cen was shown to contain an alphoid subset (D5Z12) different from those described thus far. Two genes and several putative cDNAs could be precisely located close to the centromere. Several L1 elements were shown to be present within alpha satellites at the boundary between alphoid and nonalphoid sequences on both sides of 5cen. They were used to define STSs that could serve as physical anchor points at the junction of 5cen with the p and q arms. Some STSs were placed on a radiation hybrid map. One was polymorphic and could therefore be used as a second centromeric genetic marker at the p arm boundary of 5cen. We could thus estimate recombination rates within and around the centromeric region of chromosome 5. Recombination is highly reduced within 5cen, with zero recombinants in 58 meioses being detected between the two markers located at the two extremities of the centromere. In its immediate vicinity, 5cen indeed exerts a direct negative effect on meiotic recombination within the proximal chromosomal DNA. This effect is, however, less important than expected and is polarized, as different rates are observed on both arms if one compares the 0 cM/Mb of the p proximal first 5.5 Mb and the 0.64 cM/Mb of the q proximal first 5 Mb to the sex-average 1.02 cM/Mb found throughout the entire chromosome 5. Rates then become close to the average when one goes further within the arms. Finally, most recombinants (21/22), irrespective of the arm, are of female origin, thus showing that recombination around 5cen is essentially occurring in the female lineage.

Use of a Mycobacterium tuberculosis H37Rv bacterial artificial chromosome library for genome mapping, sequencing, and comparative genomics.

The bacterial artificial chromosome (BAC) cloning system is capable of stably propagating large, complex DNA inserts in Escherichia coli. As part of the Mycobacterium tuberculosis H37Rv genome sequencing project, a BAC library was constructed in the pBeloBAC11 vector and used for genome mapping, confirmation of sequence assembly, and sequencing. The library contains about 5,000 BAC clones, with inserts ranging in size from 25 to 104 kb, representing theoretically a 70-fold coverage of the M. tuberculosis genome (4.4 Mb). A total of 840 sequences from the T7 and SP6 termini of 420 BACs were determined and compared to those of a partial genomic database. These sequences showed excellent correlation between the estimated sizes and positions of the BAC clones and the sizes and positions of previously sequenced cosmids and the resulting contigs. Many BAC clones represent linking clones between sequenced cosmids, allowing full coverage of the H37Rv chromosome, and they are now being shotgun sequenced in the framework of the H37Rv sequencing project. Also, no chimeric, deleted, or rearranged BAC clones were detected, which was of major importance for the correct mapping and assembly of the H37Rv sequence. The minimal overlapping set contains 68 unique BAC clones and spans the whole H37Rv chromosome with the exception of a single gap of approximately 150 kb. As a postgenomic application, the canonical BAC set was used in a comparative study to reveal chromosomal polymorphisms between M. tuberculosis, M. bovis, and M. bovis BCG Pasteur, and a novel 12.7-kb segment present in M. tuberculosis but absent from M. bovis and M. bovis BCG was characterized. This region contains a set of genes whose products show low similarity to proteins involved in polysaccharide biosynthesis. The H37Rv BAC library therefore provides us with a powerful tool both for the generation and confirmation of sequence data as well as for comparative genomics and other postgenomic applications. It represents a major resource for present and future M. tuberculosis research projects.

Construction and extensive characterization of a goat bacterial artificial chromosome library with threefold genome coverage.

A goat Bacterial Artificial Chromosome (BAC) library of 61,440 independent clones was constructed and characterized. The average size of the inserts was estimated at 153 kilobases by analyzing almost 500 clones using Not1 digestion followed by FIGE (Field Inverted Gel Electrophoresis) analysis. The library represents about three genome equivalents, which yields a theoretical probability of 0.95 of isolating a particular DNA sequence. After individual growth, the clones were arrayed in 40 superpools, which were organized in three dimension pools. A rapid technique for pool DNA preparation by microwave treatment was set up. This technique was compatible with PCR analysis. Primer pairs from 166 sequences (microsatellites, coding sequences from goat, and conserved Expressed Sequence Tags (ESTs) from humans) enabled the library to be successfully searched in 165 cases, with an average of 3.52 positive superpools. Only one sequence could not be found. The degree of chimerism was evaluated by FISH analysis with DNA from over 110 clones and was estimated at 4%. This BAC library will constitute an invaluable tool for positional cloning in ruminants, as well as for more general comparative mapping studies in mammals.

The CIC library: a large insert YAC library for genome mapping in Arabidopsis thaliana.

A new Arabidopsis thaliana (ecotype Columbia) genomic library has been constructed in Yeast Artificial Chromosomes: the CIC library (for CEPH, INRA and CNRS). Optimization of plant culture conditions and protoplast preparation allowed the recovery of large amounts of viable protoplasts. Mechanical shearing of DNA was minimized by isolation of DNA from protoplasts embedded in agarose. Cloning of large inserts was favored by including two successive size fractionation steps (after partial EcoRI digestion and after ligation with the vector arms), which selected DNA fragments larger than 350 kb. The library consists of 1152 clones with an average insert size of 420 kb. Clones carrying chloroplast DNA and various nuclear repeated sequences have been identified. Twenty-one per cent of the clones are found to contain chloroplast DNA. Therefore, the library represents around four nuclear genome equivalents. The clones containing 5S rDNA genes, 18S-25S rDNA sequences and the 180 bp paracentromeric repeated element account for 3.6%, 8.9% and 5.8%, respectively. Only one clone was found to carry the 160 bp paracentromeric repeated element. Given the smaller size of clones carrying Arabidopsis repeated DNA, the average size of remaining clones is around 480 kb. The library was screened by PCR amplification using pairs of primers corresponding to sequences dispersed in the genome. Seventy out of 76 pairs of primers identified from one to seven YAC clones. Thus at least 92% of the genome is represented in the CIC library. The survey of the library for clones containing unlinked DNA sequences indicates that the proportion of chimeric clones is lower than 10%.

A YAC contig map of the human genome.

A yeast artificial chromosome library containing 33,000 clones with an average insert size of one megabase of human genomic DNA was extensively analysed by several different procedures for detecting overlaps and positional information. We developed an analysis strategy that resulted, after confirmatory tests, in a YAC contig map reliably covering about 75% of the human genome in 225 contigs having an average size of about ten megabases.

Cloning the human major histocompatibility complex in YACs.

To clone the human major histocompatibility complex (MHC), 53 YACs, with an average size of 490 kb, were isolated and characterized from the CEPH YAC library. These YACs were organized in a single large contig covering more than 4000 kb. Furthermore, a complete physical map of the previously uncloned HLA class I region was established from partial and/or total digestions of 15 YACs spanning 2000 kb. This resulted in the establishment of the first YAC contig that spans the entire MHC region and constitutes an essential step in the isolation of all of the genes present in the region.

Recombinations in individuals homozygous by descent localize the Friedreich ataxia locus in a cloned 450-kb interval.

The locus for Friedreich ataxia (FRDA), a severe neurodegenerative disease, is tightly linked to markers D9S5 and D9S15, and analysis of rare recombination events has suggested the order cen-FRDA-D9S5-D9S15-qter. We report here the construction of a YAC contig extending 800 kb centromeric to D9S5 and the isolation of five new microsatellite markers from this region. In order to map these markers with respect to the FRDA locus, all within a 1-cM confidence interval, we sought to increase the genetic information of available FRDA families by considering homozygosity by descent and association with founder haplotypes in isolated populations. This approach allowed us to identify one phase-known recombination and one probable historic recombination on haplotypes from Réunion Island patients, both of which place three of the five markers proximal to FRDA. This represents the first identification of close FRDA flanking markers on the centromeric side. The two other markers allowed us to narrow the breakpoint of a previously identified distal recombination that is > 180 kb from D9S5 (26P). Taken together, the results place the FRDA locus in a 450-kb interval, which is small enough for direct search of candidate genes. A detailed rare cutter restriction map and a cosmid contig covering this interval were constructed and should facilitate the search of genes in this region.

Relationship between Charcot-Marie-Tooth 1A and Smith-Magenis regions. snU3 may be a candidate gene for the Smith-Magenis syndrome.

The juxtacentromeric region of the human chromosome 17 short arm (17p11.2-p12) contains genes involved in the Charcot-Marie-Tooth type 1A disease (CMT1A) and the Smith-Magenis syndrome (SMS). CMT1A is associated with a duplication of a short segment whereas SMS is linked to microdeletions, extending toward the centromere. We describe the construction and analysis of a 5 Mb YAC contig spanning the CMT1A duplicated segment and the distal part of four SMS microdeletions. We concluded that the YAC contig contains about 1Mb of genomic DNA which is deleted in the four SMS patients analysed. Moreover two YACs contain both STS deleted in SMS (U3) and STS duplicated in CMT1A (5H5), but the proximal breakpoint associated with the CMT1A duplication is not the same as the distal SMS breakpoint we studied. Finally we located five new STS in SMS deletion. Two of them, a microsatellite (D17S805(23)) and the gene coding for small nuclear RNA U3, have been localized in the contig we described. We may also note that snU3 is the first expressed sequence localized in an SMS deletion so far. The possible participation of this gene in the SMS phenotype is discussed.

Construction of a yeast artificial chromosome contig spanning the pseudoautosomal region and isolation of 25 new sequence-tagged sites.

Thirty-one yeast artificial chromosomes (YACs) from the human pseudoautosomal region were identified by a combination of sequence-tagged site (STS) screenings and colony hybridizations, using a subtelomeric interspersed repetitive element mapping predominantly to the pseudoautosomal region. Twenty-five new pseudoautosomal STSs were generated, of which 4 detected restriction fragment length polymorphisms. A total of 33 STSs were used to assemble the 31 YACs into a single contiguous set of overlapping DNA fragments spanning at least 2.3 megabases of the pseudoautosomal region. In addition, four pseudoautosomal genes including hydroxyindole O-methyltransferase have been positioned on this set of fragments.

Mapping the whole human genome by fingerprinting yeast artificial chromosomes.

Physical mapping of the human genome has until now been envisioned through single chromosome strategies. We demonstrate that by using large insert yeast artificial chromosomes (YACs) a whole genome approach becomes feasible. YACs (22,000) of 810 kb mean size (5 genome equivalents) have been fingerprinted to obtain individual patterns of restriction fragments detected by a LINE-1 (L1) probe. More than 1000 contigs were assembled. Ten randomly chosen contigs were validated by metaphase chromosome fluorescence in situ hybridization, as well as by analyzing the inter-Alu PCR patterns of their constituent YACs. We estimate that 15% to 20% of the human genome, mainly the L1-rich regions, is already covered with contigs larger than 3 Mb.