Construction of a 1-Mb restriction-mapped cosmid contig containing the candidate region for the familial mediterranean fever locus (MEFV) on chromosome 16p13.3

In this paper we describe the assembly and restriction map of a 1.05-Mb cosmid contig spanning the candidate region for familial Mediterranean fever (FMF), a recessively inherited disorder of inflammation localized to 16p13.3. Using a combination of cosmid walking and screening for P1, PAC, BAG, and YAC clones, we have generated a contig of genomic clones spanning similar to 1050 kb that contains the FMF critical region. The map consists of 179 cosmid, 15 P1, 10 PAC, 3 BAG, and 17 YAC clones, anchored by 27 STS markers. Eight additional STSs have been developed from the similar to 700 kb immediately centromeric to this genomic region. Five of the 35 STSs are microsatellites that have not been previously reported. NotI and EcoRI mapping of the overlapping cosmids, hybridization of restriction fragments from cosmids to one another, and STS analyses have been used to validate the assembly of the contig. Our contig totally subsumes the 250-kb interval recently reported, by founder haplotype analysis, to contain the FMF gene. Thus, our high-resolution clone map provides an ideal resource for transcriptional mapping toward the eventual identification of this disease gene. (C) 1997 Academic Press.

Phytophthora sojae avirulence genes Avr4 and Avr6 are located in a 24kb, recombination-rich region of genomic DNA

A cross between two different races (race 7 x race 25) of the soybean root and stem rot pathogen Phytophthora sojae was analyzed to characterize the genomic region flanking two cosegregating avirulence genes, Anur4 and Anur6. Both genes cosegregated in the ratio of 82:17 (avirulent:virulent) in an F-2 population, suggestive of a single locus controlling both phenotypes. A chromosome walk was commenced from RAPD marker OPE7.1C, 2.0 cM distant from the Anur4/6 locus. Three overlapping cosmids were isolated which included genetic markers that flank the Anur4/6 locus. The chromosome walk spanned a physical distance of 67 kb which represented a genetic map distance of 22.3cM, an average recombination frequency of 3.0kb/cM and 11.7-fold greater than the predicted average recombination frequency of 35.3 kb/cM for the entire P. sojae genome. Six genes (cDNA clones) expressed from the Anur4/6 genomic region encompassed by the cosmid contig were identified. Single nucleotide polymorphisms and restriction fragment length polymorphisms showed these six genes were closely linked to the Anur4/6 locus. Physical mapping of the cDNA clones within the cosmid contig made it possible to deduce the precise linkage order of the cDNAs. None of the six cDNA clones appear to be candidates for Anur4/6. We conclude that two of these cDNA clones flank a physical region of approximately 24 kb and 4.3 cM that appears to include the Anur4/6 locus. (C) 2003 Elsevier Inc. All rights reserved.

Estimation of protein coding density in a corpus of DNA sequence data

A number of experimental methods have been reported for estimating the number of genes in a genome, or the closely related coding density of a genome, defined as the fraction of base pairs in codons. Recently, DNA sequence data representative of the genome as a whole have become available for several organisms, making the problem of estimating coding density amenable to sequence analytic methods. Estimates of coding density for a single genome vary widely, so that methods with characterized error bounds have become increasingly desirable. We present a method to estimate the protein coding density in a corpus of DNA sequence data, in which a ‘coding statistic’ is calculated for a large number of windows of the sequence under study, and the distribution of the statistic is decomposed into two normal distributions, assumed to be the distributions of the coding statistic in the coding and noncoding fractions of the sequence windows. The accuracy of the method is evaluated using known data and application is made to the yeast chromosome III sequence and to C.elegans cosmid sequences. It can also be applied to fragmentary data, for example a collection of short sequences determined in the course of STS mapping.

Leishmania major: histone H1 gene expression from the sw3 locus.

Histone H1 in the parasitic protozoan Leishmania is a developmentally regulated protein encoded by the sw3 gene. Here we report that histone H1 variants exist in different Leishmania species and strains of L. major and that they are encoded by polymorphic genes. Amplification of the sw3 gene from the genome of three strains of L. major gave rise to different products in each strain, suggesting the presence of a multicopy gene family. In L. major, these genes were all restricted to a 50-kb Bg/II fragment found on a chromosomal band of 1.3 Mb (chromosome 27). The detection of RFLPs in this locus demonstrated its heterogeneity within several species and strains of Leishmania. Two different copies of sw3 (sw3.0 and sw3.1) were identified after screening a cosmid library containing L. major strain Friedlin genomic DNA. They were identical in their 5' UTRs and open reading frames, but differed in their 3' UTRs. With respect to the originally cloned copy of sw3 from L. major strain LV39, their open reading frames lacked a repeat unit of 9 amino acids. Immunoblots of L. guyanensis parasites transfected with these cosmids revealed that both copies could give rise to the histone H1 protein. The characterization of this locus will now make possible a detailed analysis of the function of histone H1 in Leishmania, as well as permit the dissection of the molecular mechanisms governing the developmental regulation of the sw3 gene.

Gene transfer of the Ly-3 chain gene of the mouse CD8 molecular complex: co-transfer with the Ly-2 polypeptide gene results in detectable cell surface expression of the Ly-3 antigenic determinants.

The CD8 molecule is a glycoprotein expressed on a subset of mature T lymphocytes. It has been postulated to be a receptor for class I major histocompatibility complex molecules. In the mouse, CD8 is a heterodimer composed of Ly-2 and Ly-3 chains. We have isolated and analyzed cDNA and cosmid clones corresponding to the Ly-3 subunit. One of the isolated, cosmid clones was subsequently transfected, alone or in combination with the Ly-2 gene, into mouse Ltk- cells. Analysis of the Ly-2,3 molecules expressed at the surface of the double transfectants indicated that they are serologically and biochemically indistinguishable from their normal counterparts expressed on lymphoid cells. Ltk- cells transfected with the Ly-2 gene alone were shown to react with a subset of anti-CD8 monoclonal antibodies whereas Ly-3 transfectants did not stain with any of the anti-Ly-3 antibodies employed in this study. Since at least one of these antibodies (53-5.8) has been previously shown to recognize an epitope which is retained on the Ly-3 subunit after dissociation of the heterodimeric Ly-2,3 complex, these observations suggest that the expression of the Ly-2 polypeptide is required to permit the detectable cell surface expression of the antigenic determinants carried by the Ly-3 subunit.

Genomic Analysis of Wild Tomato Introgressions Determining Metabolism- and Yield-Associated Traits

With the aim of determining the genetic basis of metabolic regulation in tomato fruit, we constructed a detailed physical map of genomic regions spanning previously described metabolic quantitative trait loci of a Solanum pennellii introgression line population. Two genomic libraries from S. pennellii were screened with 104 colocated markers from five selected genomic regions, and a total of 614 bacterial artificial chromosome (BAC)/cosmids were identified as seed clones. Integration of sequence data with the genetic and physical maps of Solanum lycopersicum facilitated the anchoring of 374 of these BAC/cosmid clones. The analysis of this information resulted in a genome-wide map of a nondomesticated plant species and covers 10% of the physical distance of the selected regions corresponding to approximately 1% of the wild tomato genome. Comparative analyses revealed that S. pennellii and domesticated tomato genomes can be considered as largely colinear. A total of 1,238,705 bp from both BAC/cosmid ends and nine large insert clones were sequenced, annotated, and functionally categorized. The sequence data allowed the evaluation of the level of polymorphism between the wild and cultivated tomato species. An exhaustive microsynteny analysis allowed us to estimate the divergence date of S. pennellii and S. lycopersicum at 2.7 million years ago. The combined results serve as a reference for comparative studies both at the macrosyntenic and microsyntenic levels. They also provide a valuable tool for fine-mapping of quantitative trait loci in tomato. Furthermore, they will contribute to a deeper understanding of the regulatory factors underpinning metabolism and hence defining crop chemical composition.

Sequence of a 189-kb segment of the chromosome of Rhodobacter capsulatus SB1003

Cosmids from the 1A3–1A10 region of the complete miniset were individually subcloned by using the vector M13 mp18. Sequences of each cosmid were assembled from about 400 DNA fragments generated from the ends of these phage subclones and merged into one 189-kb contig. About 160 ORFs identified by the CodonUse program were subjected to similarity searches. The biological functions of 80 ORFs could be assigned reliably by using the WIT and Magpie genome investigation tools. Eighty percent of these recognizable ORFs were organized in functional clusters, which simplified assignment decisions and increased the strength of the predictions. A set of 26 genes for cobalamin biosynthesis, genes for polyhydroxyalkanoic acid metabolism, DNA replication and recombination, and DNA gyrase were among those identified. Most of the ORFs lacking significant similarity with reference databases also were grouped. There are two large clusters of these ORFs, one located between 45 and 67 kb of the map, and the other between 150 and 183 kb. Nine of the loosely identified ORFs (of 15) of the first of these clusters match ORFs from phages or transposons. The other cluster also has four ORFs of possible phage origin.

Complete sequencing of the Fugu WAGR region from WT1 to PAX6: Dramatic compaction and conservation of synteny with human chromosome 11p13

The pufferfish Fugu rubripes has a genome ≈7.5 times smaller than that of mammals but with a similar number of genes. Although conserved synteny has been demonstrated between pufferfish and mammals across some regions of the genome, there is some controversy as to what extent Fugu will be a useful model for the human genome, e.g., [Gilley, J., Armes, N. & Fried, M. (1997) Nature (London) 385, 305–306]. We report extensive conservation of synteny between a 1.5-Mb region of human chromosome 11 and <100 kb of the Fugu genome in three overlapping cosmids. Our findings support the idea that the majority of DNA in the region of human chromosome 11p13 is intergenic. Comparative analysis of three unrelated genes with quite different roles, WT1, RCN1, and PAX6, has revealed differences in their structural evolution. Whereas the human WT1 gene can generate 16 protein isoforms via a combination of alternative splicing, RNA editing, and alternative start site usage, our data predict that Fugu WT1 is capable of generating only two isoforms. This raises the question of the extent to which the evolution of WT1 isoforms is related to the evolution of the mammalian genitourinary system. In addition, this region of the Fugu genome shows a much greater overall compaction than usual but with significant noncoding homology observed at the PAX6 locus, implying that comparative genomics has identified regulatory elements associated with this gene.

Multiple-complete-digest restriction fragment mapping: Generating sequence-ready maps for large-scale DNA sequencing

Multiple-complete-digest mapping is a DNA mapping technique based on complete-restriction-digest fingerprints of a set of clones that provides highly redundant coverage of the mapping target. The maps assembled from these fingerprints order both the clones and the restriction fragments. Maps are coordinated across three enzymes in the examples presented. Starting with yeast artificial chromosome contigs from the 7q31.3 and 7p14 regions of the human genome, we have produced cosmid-based maps spanning more than one million base pairs. Each yeast artificial chromosome is first subcloned into cosmids at a redundancy of ×15–30. Complete-digest fragments are electrophoresed on agarose gels, poststained, and imaged on a fluorescent scanner. Aberrant clones that are not representative of the underlying genome are rejected in the map construction process. Almost every restriction fragment is ordered, allowing selection of minimal tiling paths with clone-to-clone overlaps of only a few thousand base pairs. These maps demonstrate the practicality of applying the experimental and software-based steps in multiple-complete-digest mapping to a target of significant size and complexity. We present evidence that the maps are sufficiently accurate to validate both the clones selected for sequencing and the sequence assemblies obtained once these clones have been sequenced by a “shotgun” method.

In vitro reconstruction of the Aspergillus (= Emericella) nidulans genome

A physical map of the 31-megabase Aspergillus nidulans genome is reported, in which 94% of 5,134 cosmids are assigned to 49 contiguous segments. The physical map is the result of a two-way ordering process, in which clones and probes were ordered simultaneously on a binary DNA/DNA hybridization matrix. Compression by elimination of redundant clones resulted in a minimal map, which is a chromosome walk. Repetitive DNA is nonrandomly dispersed in the A. nidulans genome, reminiscent of heterochromatic banding patterns of higher eukaryotes. We hypothesize gene clusters may arise by horizontal transfer and spread by transposition to explain the nonrandom pattern of repeats along chromosomes.

Assembly of large genomic segments in artificial chromosomes by homologous recombination in Escherichia coli

We developed a method for the reconstruction of a 100 kb DNA fragment into a bacterial artificial chromosome (BAC). The procedure makes use of iterative rounds of homologous recombination in Escherichia coli. Smaller, overlapping fragments of cloned DNA, such as cosmid clones, are required. They are transferred first into a temperature-sensitive replicon and then into the BAC of choice. We demonstrated the usefulness of this procedure by assembling a 90 kb genomic segment into an E.coli–Streptomyces artificial chromosome (ESAC). Using this procedure, ESACs are easy to handle and remarkably more stable than the starting cosmids.

Gene discovery and gene function assignment in filamentous fungi

Filamentous fungi are a large group of diverse and economically important microorganisms. Large-scale gene disruption strategies developed in budding yeast are not applicable to these organisms because of their larger genomes and lower rate of targeted integration (TI) during transformation. We developed transposon-arrayed gene knockouts (TAGKO) to discover genes and simultaneously create gene disruption cassettes for subsequent transformation and mutant analysis. Transposons carrying a bacterial and fungal drug resistance marker are used to mutagenize individual cosmids or entire libraries in vitro. Cosmids are annotated by DNA sequence analysis at the transposon insertion sites, and cosmid inserts are liberated to direct insertional mutagenesis events in the genome. Based on saturation analysis of a cosmid insert and insertions in a fungal cosmid library, we show that TAGKO can be used to rapidly identify and mutate genes. We further show that insertions can create alterations in gene expression, and we have used this approach to investigate an amino acid oxidation pathway in two important fungal phytopathogens.

Unique chromosomal regions associated with virulence of an avian pathogenic Escherichia coli strain.

The avian pathogenic Escherichia coli strain (chi)7122 (serotype O78:K80:H9) causes airsacculitis and colisepticemia in chickens. To identify genes associated with avian disease, a genomic subtraction technique was performed between strain (chi)7122 and the E. coli K-12 strain (chi)289. The DNA isolated using this method was found only in strain (chi)7122 and was used to identify cosmid clones carrying unique DNA from a library of (chi)7122 that were then used to map the position of unique DNA on the E. coli chromosome. A total of 12 unique regions were found, 5 of which correspond to previously identified positions for unique DNA sequence in E. coli strains. To assess the role each unique region plays in virulence, mutants of (chi)7122 were constructed in which a segment of unique DNA was replaced with E. coli K-12 DNA by cotransduction of linked transposon insertions in DNA flanking the unique sequence. The resulting replacement mutants were assessed for inability to colonize the air sac and cause septicemia in 2-week-old white Leghorn chickens. Two mutants were found to be avirulent when injected into the right caudal air sac of 2-week-old chickens. One avirulent mutant, designated (chi)7145, carries a replacement of the rfb locus at 44 min, generating a rough phenotype. The second mutant is designated (chi)7146, and carries a replacement at position 0.0 min on the genetic map. Both mutants could be complemented to partial virulence by cosmids carrying sequences unique to (chi)7122.

Suppression subtractive hybridization: a method for generating differentially regulated or tissue-specific cDNA probes and libraries.

A new and highly effective method, termed suppression subtractive hybridization (SSH), has been developed for the generation of subtracted cDNA libraries. It is based primarily on a recently described technique called suppression PCR and combines normalization and subtraction in a single procedure. The normalization step equalizes the abundance of cDNAs within the target population and the subtraction step excludes the common sequences between the target and driver populations. In a model system, the SSH technique enriched for rare sequences over 1,000-fold in one round of subtractive hybridization. We demonstrate its usefulness by generating a testis-specific cDNA library and by using the subtracted cDNA mixture as a hybridization probe to identify homologous sequences in a human Y chromosome cosmid library. The human DNA inserts in the isolated cosmids were further confirmed to be expressed in a testis-specific manner. These results suggest that the SSH technique is applicable to many molecular genetic and positional cloning studies for the identification of disease, developmental, tissue-specific, or other differentially expressed genes.

Transposon tools for recombinant DNA manipulation: characterization of transcriptional regulators from yeast, Xenopus, and mouse.

Transposon Tn1000 has been adapted to deliver novel DNA sequences for manipulating recombinant DNA. The transposition procedure for these "tagged" Tn1000s is simple and applicable to most plasmids in current use. For yeast molecular biology, tagged Tn1000s introduce a variety of yeast selective markers and replication origins into plasmids and cosmids. In addition, the beta-globin minimal promoter and lacZ gene of Tn(beta)lac serve as a mobile reporter of eukaryotic enhancer activity. In this paper, Tn(beta)lac was used to localize a mouse HoxB-complex enhancer in transgenic mice. Other tagged transposons create Gal4 DNA-binding-domain fusions, in either Escherichia coli or yeast plasmids, for use in one- and two-hybrid tests of transcriptional activation and protein-protein interaction, respectively. With such fusions, the Saccharomyces cerevisiae Swi6 G1/S-phase transcription factor and the Xenopus laevis Pintallavis developmental regulator are shown to activate transcription. Furthermore, the same transposon insertions also facilitated mapping of the Swi6 and Pintallavis domains responsible for transcriptional activation. Thus, as well as introducing novel sequences, tagged transposons share the numerous other applications of transposition such as producing insertional mutations, creating deletion series, or serving as mobile primer sites for DNA sequencing.