Construction of a genomic DNA 'feature map' by sequencing from nested deletions: application to the HLA class I region

Genome sequencing analysis of the 1.8 Mb entire human MHC class I region

The human MHC class I region spans 1.8 Mb from the MICB gene to the HLA-F gene at the telomeric end of the HLA region. There are fewer genes recognized in this region than in the class II or class III region, probably because this region remained uncharacterized for genomic organization. Based on the 1,796,938 bp genomic sequence of the entire class I region determined in our laboratory, the complete gene structure of this region has finally emerged. This region embraces as many as 118 genes (73 known and 45 new genes) with a gene density of one gene every 15.2 kb, which is comparable to that of the gene-rich class III region. The GC content is fairly uniform throughout the class I region, being 45.8% on average, which corresponds to the isochore H1. By investigation of genetic polymorphisms in 26 out of 758 microsatellite repeats identified in the class I region, we could reduce the critical region for Behçet's disease (associated with B51) and psoriasis vulgaris (associated with Cw6) to approximately 50 kb segments, between MICA and HLA-B and between TCF19 and S, respectively. Thus, systematic large-scale genomic sequencing provides an efficient way of identifying genes and of mapping disease-susceptible genes in the genome.

The human major histocompatibility complex: 42,221 bp of genomic sequence, high-density sequence-tagged site map, evolution, and polymorphism for HLA class I

We report the isolation and characterization of newly identified yeast artificial chromosome (YAC) and bacterial artificial chromosome (BAC) clones spanning the HLA class I region between HLA-C and HLA-E and of YACs extending telomeric of HLA-F. When included with previously characterized HLA class I YACs, a contiguous stretch of over 2.4 Mb pairs including the entire class I region has been isolated as a series of overlapping YAC and BAC clones. Evidence that the cloned DNA faithfully represents the source genomic DNA was obtained by extensive characterization of the YACs and by independent isolation of two or more overlapping YACs or BACs spanning the entire region. As a result of this work, over 80 unique sequence probes were identified, the majority of which were sequenced to yield 42,221 bp of new major histocompatibility complex (MHC)-derived sequence. Some of these data were reduced to sequenced tagged site primer sets, facilitating the isolation of all or nearly all of HLA class I from a variety of genomic libraries. The sequence data were analyzed for protein coding capacity and homology to existing expressed tagged sites and tested for conservation of sequences in other mammalian genomes. These results indicated that large portions of the HLA class I region are conserved among mammals. Measurements of polymorphism within non-HLA class I loci generated additional data pointing toward information of potential relevance to MHC-associated diseases. The combined data and clones presented here set the stage for the determination of the complete nucleotide sequence of HLA class I.

Simple and efficient generation in vitro of nested deletions and inversions: Tn5 intramolecular transposition

We have exploited the intramolecular transposition preference of the Tn 5 in vitro transposition system to test its effectiveness as a tool for generation of nested families of deletions and inversions. A synthetic transposon was constructed containing an ori, an ampicillin resistance (Ampr) gene, a multi-cloning site (MCS) and two hyperactive end sequences. The donor DNA that adjoins the transposon contains a kanamycin resistance (Kanr) gene. Any Amprreplicating plasmid that has undergone a transposition event (Kans) will be targeted primarily to any insert in the MCS. Two different size targets were tested in the in vitro system. Synthetic transposon plasmids containing either target were incubated in the presence of purified transposase (Tnp) protein and transformed. Transposition frequencies (Ampr/Kans) for both targets were found to be 30-50%, of which >95% occur within the target sequence, in an apparently random manner. By a conservative estimate 10(5) or more deletions/inversions within a given segment of DNA can be expected from a single one-step 20 microl transposition reaction. These nested deletions can be used for structure-function analysis of proteins and for sequence analysis. The inversions provide nested sequencing templates of the opposite strand from the deletions.

Nucleotide sequencing analysis of the 146-kilobase segment around the IkBL and MICA genes at the centromeric end of the HLA class I region

To elucidate the complete gene structure and to identify new genes involved in the development of HLA class I antigen-associated diseases in the class I region of the human major histocompatibility complex on chromosome 6, a YAC clone (745D12) covering the 146-kb segment around the IkBL and MICA loci was isolated from a YAC library constructed from the B-cell line, BOLETH. A physical map of this region was constructed by isolation of overlapping cosmid clones derived from 745D12. Of these, five contiguous cosmids were chosen for DNA sequencing by the shotgun strategy to give a single contig of 146,601 bp from 2.8 kb telomeric of the IkBL gene to exon 6 of MICA. This region was confirmed to contain five known genes, IkBL, BAT1, MICB, P5-1, and HLA-X (class I fragment), from centromere to telomere, and their exon-intron organizations were determined. The 3.8-1 homologue gene (3.8-1-hom) showing 99.7% identity with the 3.8-1 cDNA clone, which was originally isolated using the 3.8-kb EcoRI fragment between the HLA-54/H and the HLA-G genes, was detected between MICA and MICB and was suggested to represent the cognate 3.8-1 genomic sequence from which the cDNA clone was derived. No evidence for the presence of expressed new genes could be obtained in this region by homology and EST searches or coding and exon prediction analyses. One TA microsatellite repeat spanning 2545 bases with as many as 913 repetitions was found on the centromeric side of the MICA gene and was indicated to be a potential hot spot for genetic recombination. The two segments of approximately 35 kb upstream of the MICA and MICB genes showed high sequence homology (about 85%) to each other, suggesting that segmental genome duplication including the MICA and MICB genes must have occurred during the evolution of the human MHC.

Two pathogenicity islands in uropathogenic Escherichia coli J96: cosmid cloning and sample sequencing

Many of the virulence genes of pathogenic strains of Escherichia coli are carried in large multigene chromosomal segments called pathogenicity islands (PAIs) that are absent from normal fecal and laboratory K-12 strains of this bacterium. We are studying PAIs in order to better understand factors that govern virulence and to assess how such DNA segments are gained or lost during evolution. The isolation and sample sequencing of a set of 11 cosmid clones that cover all of one and much of a second large PAI in the uropathogenic E. coli J96 are described. These PAIs were mapped to the 64- and 94-min regions of the E. coli K-12 chromosome, which differ from the locations of three PAIs identified in other pathogenic E. coli strains. Analysis of the junction sequences with E. coli K-12-like DNAs showed that the insert at 94 min is within the 3' end of a phenylalanine tRNA gene, pheR, and is flanked by a 135-bp imperfect direct repeat. Analysis of the one junction recovered from the insert at 64 min indicated that it lies near another tRNA gene, pheV. To identify possible genes unique to these PAIs, 100 independent subclones of the cosmids were made by PstI digestion and ligation into a pBS+ plasmid and used in one-pass sample DNA sequencing from primer binding sites at the cloning site in the vector DNA. Database searches of the J96 PAI-specific sequences identified numerous instances in which the cloned DNAs shared significant sequence similarities to adhesins, toxins, and other virulence determinants of diverse pathogens. Several likely insertion sequence elements (IS100, IS630, and IS911) and conjugative R1 plasmid and P4 phage genes were also found. We propose that such mobile genetic elements may have facilitated the spread of virulence determinants within PAIs among bacteria.

Restriction map of a 35-kb HLA fragment constructed by nested deletion 'drop-out' mapping

An efficient method for generating detailed restriction maps of large cloned DNA segments is demonstrated. The mapping strategy entails comparing restriction fragments from a parent clone and from nested deletion derivatives of that clone. In a set of deletion plasmids of decreasing size, an individual fragment will be lost, or 'drop-out', according to its position in the cloned fragment. In this demonstration, nested deletions were generated in both directions in a 35-kb DNA segment from the human leukocyte antigen (HLA) region by intramolecular transposition of an engineered gamma delta (Tn1000) element present in a special 'deletion factory' cloning vector [Wang et al., Proc. Natl. Acad. Sci. USA 90 (1993) 7874-7878]. Fifteen plasmids with deletions extending in one direction and eleven plasmids with deletions extending in the opposite direction were digested singly by each of four restriction enzymes. A total of 36 cleavage sites were mapped in the 35-kb HLA fragment. This drop-out approach using nested deletions provides a simple and efficient means of mapping restriction sites, genes and other features of interest in cosmid-sized cloned DNA segments or DNAs.