Highly plastic genome of Microcystis aeruginosa PCC 7806, a ubiquitous toxic freshwater cyanobacterium

Background

The colonial cyanobacterium Microcystis proliferates in a wide range of freshwater ecosystems and is exposed to changing environmental factors during its life cycle. Microcystis blooms are often toxic, potentially fatal to animals and humans, and may cause environmental problems. There has been little investigation of the genomics of these cyanobacteria.

Results

Deciphering the 5,172,804 bp sequence of Microcystis aeruginosa PCC 7806 has revealed the high plasticity of its genome: 11.7% DNA repeats containing more than 1,000 bases, 6.8% putative transposases and 21 putative restriction enzymes. Compared to the genomes of other cyanobacterial lineages, strain PCC 7806 contains a large number of atypical genes that may have been acquired by lateral transfers. Metabolic pathways, such as fermentation and a methionine salvage pathway, have been identified, as have genes for programmed cell death that may be related to the rapid disappearance of Microcystis blooms in nature. Analysis of the PCC 7806 genome also reveals striking novel biosynthetic features that might help to elucidate the ecological impact of secondary metabolites and lead to the discovery of novel metabolites for new biotechnological applications. M. aeruginosa and other large cyanobacterial genomes exhibit a rapid loss of synteny in contrast to other microbial genomes.

Conclusion

Microcystis aeruginosa PCC 7806 appears to have adopted an evolutionary strategy relying on unusual genome plasticity to adapt to eutrophic freshwater ecosystems, a property shared by another strain of M. aeruginosa (NIES-843). Comparisons of the genomes of PCC 7806 and other cyanobacterial strains indicate that a similar strategy may have also been used by the marine strain Crocosphaera watsonii WH8501 to adapt to other ecological niches, such as oligotrophic open oceans.

The grapevine genome sequence suggests ancestral hexaploidization in major angiosperm phyla

The analysis of the first plant genomes provided unexpected evidence for genome duplication events in species that had previously been considered as true diploids on the basis of their genetics. These polyploidization events may have had important consequences in plant evolution, in particular for species radiation and adaptation and for the modulation of functional capacities. Here we report a high-quality draft of the genome sequence of grapevine (Vitis vinifera) obtained from a highly homozygous genotype. The draft sequence of the grapevine genome is the fourth one produced so far for flowering plants, the second for a woody species and the first for a fruit crop (cultivated for both fruit and beverage). Grapevine was selected because of its important place in the cultural heritage of humanity beginning during the Neolithic period. Several large expansions of gene families with roles in aromatic features are observed. The grapevine genome has not undergone recent genome duplication, thus enabling the discovery of ancestral traits and features of the genetic organization of flowering plants. This analysis reveals the contribution of three ancestral genomes to the grapevine haploid content. This ancestral arrangement is common to many dicotyledonous plants but is absent from the genome of rice, which is a monocotyledon. Furthermore, we explain the chronology of previously described whole-genome duplication events in the evolution of flowering plants.

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.

YAPI, a new Yersinia pseudotuberculosis pathogenicity island

Pathogenicity islands (PAIs) are chromosomal clusters of pathogen-specific virulence genes often found at tRNA loci. In the Yersinia pseudotuberculosis 32777 chromosome, we characterized a 98-kb segment that has all of the characteristic features of a PAI, including insertion in a (phenylalanine) tRNA gene, the presence of a bacteriophage-like integrase-encoding gene, and direct repeats at the integration sites. The G+C content of the segment ranges from 31 to 60%, reflecting a genetic mosaic: this is consistent with the notion that the sequences were horizontally acquired. The PAI, termed YAPI (for Yersinia adhesion pathogenicity island), carries 95 open reading frames and includes (i) the previously described pil operon, encoding a type IV pilus that contributes to pathogenicity (F. Collyn et al., Infect. Immun. 70:6196-6205, 2002); (ii) a block of genes potentially involved in general metabolism; (iii) a gene cluster for a restriction-modification system; and (iv) a large number of mobile genetic elements. Furthermore, the PAI can excise itself from the chromosome at low frequency and in a precise manner, and deletion does not result in a significant decrease of bacterial virulence compared to inactivation of the fimbrial gene cluster alone. The prevalence and size of the PAI vary from one Y. pseudotuberculosis strain to another, and it can be found integrated into either of the two phe tRNA loci present on the species' chromosome. YAPI was not detected in the genome of the genetically closely related species Y. pestis, whereas a homologous PAI is harbored by the Y. enterocolitica chromosome.

Characterization of the genomic organization of the region bordering the centromere of chromosome V of Podospora anserina by direct sequencing

A Podospora anserina BAC library of 4800 clones has been constructed in the vector pBHYG allowing direct selection in fungi. Screening of the BAC collection for centromeric sequences of chromosome V allowed the recovery of clones localized on either sides of the centromere, but no BAC clone was found to contain the centromere. Seven BAC clones containing 322,195 and 156,244bp from either sides of the centromeric region were sequenced and annotated. One 5S rRNA gene, 5 tRNA genes, and 163 putative coding sequences (CDS) were identified. Among these, only six CDS seem specific to P. anserina. The gene density in the centromeric region is approximately one gene every 2.8kb. Extrapolation of this gene density to the whole genome of P. anserina suggests that the genome contains about 11,000 genes. Synteny analyses between P. anserina and Neurospora crassa show that co-linearity extends at the most to a few genes, suggesting rapid genome rearrangements between these two species.

Initiation of the breakage-fusion-bridge mechanism through common fragile site activation in human breast cancer cells: the model of PIP gene duplication from a break at FRA7I

Gene amplification plays a critical role in tumor progression. Hence, understanding the factors triggering this process in human cancers is an important concern. Unfortunately, the structures formed at early stages are usually unavailable for study, hampering the identification of the initiating events in tumors. Here, we show that the region containing the PIP gene, which is overexpressed in 80% of primary and metastatic breast cancers, is duplicated in the breast carcinoma cell line T47D. The two copies are organized as a large palindrome, lying 'in loco' on one chromosome 7. Such features constitute the landmark of the breakage-fusion-bridge (BFB) cycle mechanism. In hamster cells selected in vitro to resist cytotoxic drugs, common fragile site (CFS) activation has been shown to trigger this mechanism. Here, we characterize FRA7I at the molecular level and demonstrate that it lies 2 Mb telomeric to the PIP gene and sets the distal end of the repeated sequence. Moreover, our results suggest that the BFB process was frozen within the first cycle by healing of the broken chromosome. T47D cells thus offer a unique opportunity to observe the earliest products of the BFB cycle mechanism. Our findings constitute the first evidence that this amplification mechanism can be initiated in vivo by fragile site activation.

Bacterial artificial chromosome-based comparative genomic analysis identifies Mycobacterium microti as a natural ESAT-6 deletion mutant

Mycobacterium microti is a member of the Mycobacterium tuberculosis complex that causes tuberculosis in voles. Most strains of M. microti are harmless for humans, and some have been successfully used as live tuberculosis vaccines. In an attempt to identify putative virulence factors of the tubercle bacilli, genes that are absent from the avirulent M. microti but present in human pathogen M. tuberculosis or Mycobacterium bovis were searched for. A minimal set of 50 bacterial artificial chromosome (BAC) clones that covers almost all of the genome of M. microti OV254 was constructed, and individual BACs were compared to the corresponding BACs from M. bovis AF2122/97 and M. tuberculosis H37Rv. Comparison of pulsed-field gel-separated DNA digests of BAC clones led to the identification of 10 regions of difference (RD) between M. microti OV254 and M. tuberculosis. A 14-kb chromosomal region (RD1(mic)) that partly overlaps the RD1 deletion in the BCG vaccine strain was missing from the genomes of all nine tested M. microti strains. This region covers 13 genes, Rv3864 to Rv3876, in M. tuberculosis, including those encoding the potent ESAT-6 and CFP-10 antigens. In contrast, RD5(mic), a region that contains three phospholipase C genes (plcA to -C), was missing from only the vole isolates and was present in M. microti strains isolated from humans. Apart from RD1(mic) and RD5(mic) other M. microti-specific deleted regions have been identified (MiD1 to MiD3). Deletion of MiD1 has removed parts of the direct repeat region in M. microti and thus contributes to the characteristic spoligotype of M. microti strains.

Cytogenetic localization of 136 genes in the horse: comparative mapping with the human genome

The aim of this study was to increase the number of type I markers on the horse cytogenetic map and to improve comparison with maps of other species, thus facilitating positional candidate cloning studies. BAC clones from two different sources were FISH mapped: homologous horse BAC clones selected from our newly extended BAC library using consensus primer sequences and heterologous goat BAC clones. We report the localization of 136 genes on the horse cytogenetic map, almost doubling the number of cytogenetically mapped genes with 48 localizations from horse BAC clones and 88 from goat BAC clones. For the first time, genes were mapped to ECA13p, ECA29, and probably ECA30. A total of 284 genes are now FISH mapped on the horse chromosomes. Comparison with the human map defines 113 conserved segments that include new homologous segments not identified by Zoo-FISH on ECA7 and ECA13p.

Analysis of molecular markers genetically linked to the Leptosphaeria maculans avirulence gene AvrLm1 in field populations indicates a highly conserved event leading to virulence on Rlm1 genotypes

Map-based cloning of the avirulence gene AvrLm1 of Leptosphaeria maculans was initiated utilizing a genetic map of the fungus and a BAC library constructed from an AvrLm1 isolate. Seven polymorphic DNA markers closely linked to AvrLm1 were identified. Of these, two were shown to border the locus on its 5' end and were present, with size polymorphism, in both the virulent and the avirulent isolates. In contrast, three markers, J19-1.1, J53-1.3 (in coupling phase with avirulence), and Vir1 (in repulsion phase with avirulence), cosegregated with AvrLm1 in 312 progeny from five in vitro crosses. J19-1.1 and J53-1.3 were never amplified in the virulent parents or progeny, whereas Vir1 was never amplified in the avirulent parents or progeny. J19-1.1 and J53-1.3 were shown to be separated by 40 kb within a 184-kb BAC contig. In addition, the 1.6-cM genetic distance between J53-1.3 and the nearest recombinant marker corresponded to a 121-kb physical distance. When analyzing a European Union-wide collection of 192 isolates, J53-1.3, J19-1.1, and Vir1 were found to be closely associated with the AvrLm1 locus. The results of polymerase chain reaction amplification with primers for the three markers were in accordance with the interaction phenotype for 92.2% (J53-1.3), 90.6% (J19-1.1), and 88.0% (Vir1) of the isolates. In addition, genome organization of the AvrLm1 region was highly conserved in field isolates, because 89.1% of the avirulent isolates and 79.0% of the virulent isolates showed the same association of markers as that of the parents of in vitro crosses. The large-scale analysis of field isolates with markers originating from the genetic map therefore confirms (i) the physical proximity between the markers and the target locus and (ii) that AvrLm1 is located in (or close to) a recombination-deficient genome region. As a consequence, map-based markers provided us with high-quality markers for an overview of the occurrence of race "AvrLm1" at the field scale. These data were used to propose hypotheses on evolution towards virulence in field isolates.

Strategy to sequence the genome of Corynebacterium glutamicum ATCC 13032: use of a cosmid and a bacterial artificial chromosome library

The initial strategy of the Corynebacterium glutamicum genome project was to sequence overlapping inserts of an ordered cosmid library. High-density colony grids of approximately 28 genome equivalents were used for the identification of overlapping clones by Southern hybridization. Altogether 18 contiguous genomic segments comprising 95 overlapping cosmids were assembled. Systematic shotgun sequencing of the assembled cosmid set revealed that only 2.84 Mb (86.6%) of the C. glutamicum genome were represented by the cosmid library. To obtain a complete genome coverage, a bacterial artificial chromosome (BAC) library of the C. glutamicum chromosome was constructed in pBeloBAC11 and used for genome mapping. The BAC library consists of 3168 BACs and represents a theoretical 63-fold coverage of the C. glutamicum genome (3.28 Mb). Southern screening of 2304 BAC clones with PCR-amplified chromosomal markers and subsequent insert terminal sequencing allowed the identification of 119 BACs covering the entire chromosome of C. glutamicum. The minimal set representing a 100% genome coverage contains 44 unique BAC clones with an average overlap of 22 kb. A total of 21 BACs represented linking clones between previously sequenced cosmid contigs and provided a valuable tool for completing the genome sequence of C. glutamicum.

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.

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 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.

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.

Isolation of chromosome 21-specific yeast artificial chromosomes from a total human genome library

A new approach for the isolation of chromosome-specific subsets from a human genomic yeast artificial chromosome (YAC) library is described. It is based on the hybridization with an Alu polymerase chain reaction (PCR) probe. We screened a 1.5 genome equivalent YAC library of megabase insert size with Alu PCR products amplified from hybrid cell lines containing human chromosome 21, and identified a subset of 63 clones representative of this chromosome. The majority of clones were assigned to chromosome 21 by the presence of specific STSs and in situ hybridization. Twenty-nine of 36 STSs that we tested were detected in the subset, and a contig spanning 20 centimorgans in the genetic map and containing 8 STSs in 4 YACs was identified. The proposed approach can greatly speed efforts to construct physical maps of the human genome.

The CEPH YAC library

Because of their large size, YACs provide is a powerful tool for physical mapping studies of complex genomes. As it will be advantageous to have genomic libraries of clones with large inserts for analyzing megabase sized regions of the human genome, we have investigated a number of parameters in order to increase the insert size of the YACs. We constructed a genomic library currently containing more than 85,000 YAC clones. Mean sizes of YACs produced at several stages of construction of the library range from 430 kb to 1,200 kb, representing 13 haploid equivalents of the human genome. This library was organized in order to allow rapid screen of YACs for large scale physical mapping of the human genome.

Physical linkage of a guanine nucleotide-binding protein-related gene to the chicken major histocompatibility complex

Several genes were found closely associated with major histocompatibility class I and class II beta-chain genes in chicken genomic DNA clusters by hybridizing tissue-specific cDNA probes to cosmid clones. A cDNA probe for one of these genes, probe C12.3 isolated from a chicken liver cDNA library, was used to clone the homologous sequence H12.3 from a human B-lymphoblastoid cell line cDNA library. C12.3 and H12.3 encode exactly the same 317-residue-long protein. The sequence of 12.3 shows significant homology with the two known guanine nucleotide-binding protein beta subunits (GP beta 1 and GP beta 2) and other proteins that all share the same segmented structure with seven internal homologous repeats about 45 residues in length. Unlike the chicken gene, the human H12.3 gene and its mouse counterpart are not located on the same chromosome as the major histocompatibility complex. A possible involvement of the C12.3 gene product in major histocompatibility complex-linked control of lymphocyte proliferation in chickens is discussed.

A molecular map of the chicken major histocompatibility complex: the class II beta genes are closely linked to the class I genes and the nucleolar organizer

We have cloned in a cosmid vector four DNA clusters covering 320 kb of the chicken MHC (B complex), including five class II (B-L) beta genes defining two related isotypic families. Additional B complex genes have been revealed using tissue-specific cDNA probes. A cosmid fragment has been used to isolate a cDNA for a class I (B-F) transcript. This transcript, that is by far the most divergent known member of the class I gene family, hybridized to six B-F genes present in the cosmids. One of the clusters was shown to contain two rRNA transcriptional units from the nucleolar organizer region (NOR), marking the telomeric boundary of the B complex. None of the other B complex genes hybridizes to, or has the transcriptional characteristics of mammalian MHC class II alpha or class III genes. The map we have obtained shows that the B complex does not contain well defined class I and class II regions since B-F and B-L beta genes are closely associated with unrelated genes. Moreover, class II beta genes are very closely linked to class I genes in two clusters, and to the NOR in a third one.

Isolation of chicken major histocompatibility complex class II (B-L) beta chain sequences: comparison with mammalian beta chains and expression in lymphoid organs

By cross-hybridization in low stringency conditions, using a probe derived from an HLA-DQ beta cDNA clone, we have isolated several chicken genomic DNA clones. These clones were mapped to the major histocompatibility complex (MHC) of the chick (B complex) by virtue of their ability to detect restriction enzyme length polymorphisms between congenic lines of chicken. Evidence was obtained for the presence of at least three B-L beta genes in the chicken genome. The B-L beta genes are transcribed specifically in tissues containing cells of the B lymphocyte and myeloid lineages and expressing the B-L antigens. Exons encoding the beta 1, beta 2 and transmembrane domains of a B-L beta chain have been identified with 63, 66 and 62% similarity with the HLA-DQ beta sequence. This first isolation of an MHC class II gene outside of the mammalian class provides insight into the evolution of MHC genes based on the comparison of avian and mammalian class II beta chain amino acid and nucleotide sequences.

A complete and a truncated U1 snRNA gene of Drosophila melanogaster are found as inverted repeats at region 82E of the polytene chromosomes

A phage containing two sequences homologous to U1 snRNA was isolated from a Drosophila melanogaster genomic library, and identified with a previously cloned D. melanogaster U1 snRNA gene. DNA sequence analysis showed that complete and truncated U1 snRNA genes are present, both of which have base substitutions relative to U1 snRNA. These genes show conservation of 5' and 3' flanking regions relative to other U1 and U2 snRNA genes of Drosophila. Intramolecular renaturation experiments and electron microscope mapping demonstrates that the two U1 snRNA sequences are present as inverted repeats about 2.7kb apart, separated by a smaller pair of inverted repeats of an unrelated sequence. These U1 snRNA sequences were located by in situ hybridization at 82E, and related sequences were found at 21D and 95C on the polytene chromosome map. The results are discussed with reference to the origin and function of snRNAs.

Construction of a colony bank of E. coli containing hybrid plasmids representative of the Bacillus subtilis 168 genome. Expression of functions harbored by the recombinant plasmids in B. subtilis

A collection of about 2500 clones containing hybrid plasmids representative of nearly the entire genome of B. subtilis 168 was established in E. coli SK1592 by using the poly(dA).poly(dT) joining method with randomly sheared DNA fragments and plasmid pHV33, a bifunctional vector which can replicate in both E. coli and B. subtilis. Detection of cloned recombinant DNA molecules was based on the insertional inactivation of the Tc gene occurring at the unique BamHI cleavage site present in the vector plasmid. Thirty individual clones of the collection were shown to hybridize specifically with a B. subtilis rRNA probe. CCC-recombinant plasmids extracted from E. coli were pooled in lots of 100 and used to transform auxotrophic mutants of B. subtilis 168. Complementation of these auxotrophic mutations was observed for several markers such at thr, leuA, hisA, glyB and purB. In several cases, markers carried by the recombinant plasmids were lost from the plasmid and integrated into the chromosomal DNA. Loss of genetic markers from the hybrid plasmids did not occur when a rec- recipient strain of B. subtilis was used.

Construction of a kit of reference strains for rapid genetic mapping in Bacillus subtilis 168

A set of nine reference strains bringing convenient markers in the genetic background of Bacillus subtilis Marburg 168 has been prepared to allow rapid mapping of new markers.

Revision of the linkage map of Bacillus subtilis 168: indications for circularity of the chromosome

A revision of the linkage map of the Bacillus subtilis 168 chromosome has been undertaken with the use of the generalized transducing phage PBS1. The mapping of four new markers (narB1, mtlB1, aroI906, and tre-12) has allowed a determination of the relative orientation of the purB-dal segment and its linkage with the lin markers. The chromosomal segment comprised between the sacQ36 and gtaA12 markers has been linked with the narA1, ctrA1, and sacA321 markers. The recA1 marker has been mapped relative to the thyA and citB17 markers. Indications of linkage have been found between the tre-12 and catA markers and the aroG932 and sacQ36 markers. According to these results, a circular genetic map of the chromosome of B. subtilis 168 is presented. Taken together, the transduction data and the order of marker replication determined by Harford in the accompanying paper support strongly the hypothesis of a symmetrical and fully bidirectional mode of replication for the B. subtilis 168 chromosome.