Physical mapping of human chromosomes by repetitive sequence fingerprinting

Sister chromatid cohesion is mediated by individual cohesin complexes

Eukaryotic genomes are organized by loop extrusion and sister chromatid cohesion, both mediated by the multimeric cohesin protein complex. Understanding how cohesin holds sister DNAs together, and how loss of cohesion causes age-related infertility in females, requires knowledge as to cohesin's stoichiometry in vivo. Using quantitative super-resolution imaging, we identified two discrete populations of chromatin-bound cohesin in postreplicative human cells. Whereas most complexes appear dimeric, cohesin that localized to sites of sister chromatid cohesion and associated with sororin was exclusively monomeric. The monomeric stoichiometry of sororin:cohesin complexes demonstrates that sister chromatid cohesion is conferred by individual cohesin rings, a key prediction of the proposal that cohesion arises from the co-entrapment of sister DNAs.

CRISPR/Cas systems for in situ imaging of intracellular nucleic acids: Concepts and applications

Accurate and precise localization of intracellular nucleic acids is crucial for regulating genetic information transcription and diagnosing diseases. Although intracellular nucleic acid imaging methods are available for various cell types, their widespread utilization is impeded by the intricate nature of the process and its exorbitant cost. Recently, numerous intracellular nucleic acid labeling techniques based on clustered regularly interspaced short palindromic repeats (CRISPR) have been established due to their modularity, flexibility, and specificity. In this work, we present various CRISPR methods that are currently employed for visualizing intracellular genomic sequences and RNA, based on their detection principles and application scenarios. Furthermore, we discuss the advantages and drawbacks of the existing CRISPR imaging methods, as well as future research directions. We anticipate that with continued refinement, more advanced CRISPR-based imaging techniques can be developed to better elucidate the localization and dynamics of intracellular nucleic acids, thereby providing a powerful tool for molecular biology research and clinical molecular pathology diagnosis.

Probing mercury(II)-DNA interactions by nanopore stochastic sensing

In this work, DNA-Hg(II) interactions were investigated by monitoring the translocation of DNA hairpins in a protein ion channel in the absence and presence of metal ions. Our experiments demonstrate that target-specific hairpin structures could be stabilized much more significantly by mercuric ions than by the stem length and the loop size of the hairpin due to the formation of Thymine-Hg(II)-Thymine complexes. In addition, the designed DNA probe allows the development of a highly sensitive nanopore sensor for Hg(2+) with a detection limit of 25 nM. Further, the sensor is specific, and other tested metal ions including Pb(2+), Cu(2+), Cd(2+), and so on with concentrations of up to 2 orders of magnitude greater than that of Hg(2+) would not interfere with the mercury detection.

De novo transcriptome sequencing reveals a considerable bias in the incidence of simple sequence repeats towards the downstream of 'Pre-miRNAs' of black pepper

Next generation sequencing has an advantageon transformational development of species with limited available sequence data as it helps to decode the genome and transcriptome. We carried out the de novo sequencing using illuminaHiSeq™ 2000 to generate the first leaf transcriptome of black pepper (Piper nigrum L.), an important spice variety native to South India and also grown in other tropical regions. Despite the economic and biochemical importance of pepper, a scientifically rigorous study at the molecular level is far from complete due to lack of sufficient sequence information and cytological complexity of its genome. The 55 million raw reads obtained, when assembled using Trinity program generated 2,23,386 contigs and 1,28,157 unigenes. Reports suggest that the repeat-rich genomic regions give rise to small non-coding functional RNAs. MicroRNAs (miRNAs) are the most abundant type of non-coding regulatory RNAs. In spite of the widespread research on miRNAs, little is known about the hair-pin precursors of miRNAs bearing Simple Sequence Repeats (SSRs). We used the array of transcripts generated, for the in silico prediction and detection of ‘43 pre-miRNA candidates bearing different types of SSR motifs’. The analysis identified 3913 different types of SSR motifs with an average of one SSR per 3.04 MB of thetranscriptome. About 0.033% of the transcriptome constituted ‘pre-miRNA candidates bearing SSRs’. The abundance, type and distribution of SSR motifs studied across the hair-pin miRNA precursors, showed a significant bias in the position of SSRs towards the downstream of predicted ‘pre-miRNA candidates’. The catalogue of transcripts identified, together with the demonstration of reliable existence of SSRs in the miRNA precursors, permits future opportunities for understanding the genetic mechanism of black pepper and likely functions of ‘tandem repeats’ in miRNAs.

Chromosomes in the flow to simplify genome analysis

Nuclear genomes of human, animals, and plants are organized into subunits called chromosomes. When isolated into aqueous suspension, mitotic chromosomes can be classified using flow cytometry according to light scatter and fluorescence parameters. Chromosomes of interest can be purified by flow sorting if they can be resolved from other chromosomes in a karyotype. The analysis and sorting are carried out at rates of 10(2)-10(4) chromosomes per second, and for complex genomes such as wheat the flow sorting technology has been ground-breaking in reducing genome complexity for genome sequencing. The high sample rate provides an attractive approach for karyotype analysis (flow karyotyping) and the purification of chromosomes in large numbers. In characterizing the chromosome complement of an organism, the high number that can be studied using flow cytometry allows for a statistically accurate analysis. Chromosome sorting plays a particularly important role in the analysis of nuclear genome structure and the analysis of particular and aberrant chromosomes. Other attractive but not well-explored features include the analysis of chromosomal proteins, chromosome ultrastructure, and high-resolution mapping using FISH. Recent results demonstrate that chromosome flow sorting can be coupled seamlessly with DNA array and next-generation sequencing technologies for high-throughput analyses. The main advantages are targeting the analysis to a genome region of interest and a significant reduction in sample complexity. As flow sorters can also sort single copies of chromosomes, shotgun sequencing DNA amplified from them enables the production of haplotype-resolved genome sequences. This review explains the principles of flow cytometric chromosome analysis and sorting (flow cytogenetics), discusses the major uses of this technology in genome analysis, and outlines future directions.

Non-B DNA structure-induced genetic instability and evolution

Repetitive DNA motifs are abundant in the genomes of various species and have the capacity to adopt non-canonical (i.e., non-B) DNA structures. Several non-B DNA structures, including cruciforms, slipped structures, triplexes, G-quadruplexes, and Z-DNA, have been shown to cause mutations, such as deletions, expansions, and translocations in both prokaryotes and eukaryotes. Their distributions in genomes are not random and often co-localize with sites of chromosomal breakage associated with genetic diseases. Current genome-wide sequence analyses suggest that the genomic instabilities induced by non-B DNA structure-forming sequences not only result in predisposition to disease, but also contribute to rapid evolutionary changes, particularly in genes associated with development and regulatory functions. In this review, we describe the occurrence of non-B DNA-forming sequences in various species, the classes of genes enriched in non-B DNA-forming sequences, and recent mechanistic studies on DNA structure-induced genomic instability to highlight their importance in genomes.

Non-random features of loop-size chromatin fragmentation

Upon isolation of DNA from normal eukaryotic cells by standard methods involving extensive proteolytic treatment, a rather homogeneous population of loop-size, double-stranded DNA fragments is regularly obtained. These DNA molecules can be efficiently end-labeled by the DNA polymerase I Klenow fragment, as well as by a 3'- to -5'-exonuclease-free Klenow enzyme, but not by terminal transferase (TdT) unless the ends have been filled up by Klenow, suggesting that dominantly 5' protruding termini are generated upon fragmentation. The filled-up termini were used for cloning the distal parts of the approximately 50 kb fragments. BLAST analysis of the sequence of several clones allowed us to determine the sequence of the non-cloned side of the breakpoints. Comparison of 25, 600 bp-long breakpoint sequences demonstrated prevalence of repetitive elements. Consensus motives characteristic of the breakpoint sequences have been identified. Several sequences exhibit peculiar computed conformational characteristics, with sharp transition or center of symmetry located exactly at the breakpoint. Our data collectively suggest that chromatin fragmentation involves nucleolytic cleavages at fragile/hypersensitive sites delimiting loop-size fragments in a non-random manner. Interestingly, the sequence characteristics of the breakpoints are reminiscent of certain breakpoint cluster regions frequently subject to gene rearrangements.

An SSLP marker-anchored BAC framework map of the mouse genome

We have constructed a BAC framework map of the mouse genome consisting of 2,808 PCR-confirmed BAC clusters, using a previously described method. Fingerprints of BACs from selected clusters confirm the accuracy of the map. Combined with BAC fingerprint data, the framework map covers 37% of the mouse genome.

ODS2: a multiplatform software application for creating integrated physical and genetic maps

A contig map is a physical map that shows the native order of a library of overlapping genomic clones. One common method for creating such maps involves using hybridization to detect clone overlaps. False- positive and false-negative hybridization errors, the presence of chimeric clones, and gaps in library coverage lead to ambiguity and error in the clone order. Genomes with good genetic maps, such as Neurospora crassa, provide a means for reducing ambiguities and errors when constructing contig maps if clones can be anchored with genetic markers to the genetic map. A software application called ODS2 for creating contig maps based on clone-clone hybridization data is presented. This application is also designed to exploit partial ordering information provided by anchorage of clones to a genetic map. This information, along with clone-clone hybridization data, is used by a clone ordering algorithm and is represented graphically, allowing users to interactively align physical and genetic maps. ODS2 has a graphical user interface and is implemented entirely in Java, so it runs on multiple platforms. Other features include the flexibility of storing data in a local file or relational database and the ability to create full or minimum tiling contig maps.

Molecular characterization and chromosomal distribution of species-specific repetitive DNA sequences from Beta corolliflora, a wild relative of sugar beet

Repetitive DNA sequences have been isolated from a Sau3AI plasmid library of tetraploid Beta corolliflora (2n = 4x = 36), a wild relative of sugar beet (B. vulgaris). The library was screened by differential hybridization with genomic DNA of B. corolliflora and B. vulgaris. When used as probes for Southern hybridization of genomic DNA, six clones were determined to represent highly repetitive DNA families present only in the B. corolliflora genome. Five other sequences were highly repetitive in B. corolliflora and low or single copy in B. vulgaris. The insert size varied between 43 bp and 448 bp. Two sequences pBC1279 and pBC1944 displayed strong homology to a previously cloned satellite DNA from B. nana. With one exception, sequences are tandemly arranged as revealed by a typical ladder pattern after genomic Southern hybridization. The chromosomal distribution of five probes was determined by fluorescence in situ hybridization (FISH) of mitotic metaphases from B. corolliflora and a triploid hybrid between B. vulgaris and B. corolliflora. Three sequences were spread along all chromosome arms of B. corolliflora while one sequence was present on only six chromosomes. The chromosome-specific sequence pBC216 was found in close vicinity to the 5S rDNA located on B. corolliflora chromosome IV. This set of species-specific sequences has the potential to be used as probes for the identification of monosomic alien addition lines and for marker-assisted gene transfer from wild beet to cultivated beet.

Contigs built with fingerprints, markers, and FPC V4.7

Contigs have been assembled, and over 2800 clones selected for sequencing for human chromosomes 9, 10 and 13. Using the FPC (FingerPrinted Contig) software, the contigs are assembled with markers and complete digest fingerprints, and the contigs are ordered and localised by a global framework. Publicly available resources have been used, such as, the 1998 International Gene Map for the framework and the GSC Human BAC fingerprint database for the majority of the fingerprints. Additional markers and fingerprints are generated in-house to supplement this data. To support the scale up of building maps, FPC V4.7 has been extended to use markers with the fingerprints for assembly of contigs, new clones and markers can be automatically added to existing contigs, and poorly assembled contigs are marked accordingly. To test the automatic assembly, a simulated complete digest of 110 Mb of concatenated human sequence was used to create datasets with varying coverage, length of clones, and types of error. When no error was introduced and a tolerance of 7 was used in assembly, the largest contig with no false positive overlaps has 9534 clones with 37 out-of-order clones, that is, the starting coordinates of adjacent clones are in the wrong order. This paper describes the new features in FPC, the scenario for building the maps of chromosomes 9, 10 and 13, and the results from the simulation.

Contigs built with fingerprints, markers, and FPC V4.7

Contigs have been assembled, and over 2800 clones selected for sequencing for human chromosomes 9, 10 and 13. Using the FPC (FingerPrinted Contig) software, the contigs are assembled with markers and complete digest fingerprints, and the contigs are ordered and localised by a global framework. Publicly available resources have been used, such as, the 1998 International Gene Map for the framework and the GSC Human BAC fingerprint database for the majority of the fingerprints. Additional markers and fingerprints are generated in-house to supplement this data. To support the scale up of building maps, FPC V4.7 has been extended to use markers with the fingerprints for assembly of contigs, new clones and markers can be automatically added to existing contigs, and poorly assembled contigs are marked accordingly. To test the automatic assembly, a simulated complete digest of 110 Mb of concatenated human sequence was used to create datasets with varying coverage, length of clones, and types of error. When no error was introduced and a tolerance of 7 was used in assembly, the largest contig with no false positive overlaps has 9534 clones with 37 out-of-order clones, that is, the starting coordinates of adjacent clones are in the wrong order. This paper describes the new features in FPC, the scenario for building the maps of chromosomes 9, 10 and 13, and the results from the simulation.

High-resolution transcript map of the region spanning D12S1629 and D12S312 at chromosome 12q13: triple A syndrome-linked region

For those searching for human disease-causing genes, information on the position of genes with respect to genetic markers is essential. The physical map composed of ESTs and genetic markers provides the positional information of these markers as well as the starting point of gene identification in the form of genomic clones containing exons. To facilitate the effort of identification of genes in the region spanning D12S1629 and D12S312, we constructed a high-resolution transcript map with PAC/BAC/cosmid clones. The strategy for the construction of such a map involved utilization of STSs for the screening of the large insert bacterial chromosome libraries and a chromosome 12-specific cosmid library by hybridization. The contig was constructed based on the STS contents of the clones. The resulting high-resolution transcript map of the region between P273P14/SP6 and D12S312 spans 4.4 cM from 66.8 to 71.2 cM of the Généthon genetic map and represents approximately 2.4 Mb. It was composed of 81 BAC, 45 PAC, and 91 cosmid clones with a minimal tiling path consisting of 16 BAC and 4 PAC clones. These clones are being used to sequence this part of chromosome 12. We determined the order of 135 STSs including 74 genes and ESTs in the map. Among these, 115 STSs were unambiguously ordered, resulting in one ordered marker per 21 kb. The order of keratin type II locus genes was determined. This map would greatly enhance the positional cloning effort of the responsible genes for those diseases that are linked to this region, including male germ cell tumor as well as palmoplantar keratoderma, Bothnian-type, and triple A syndrome. This transcript map was localized at human chromosome 12q13.

Contig assembly of bacterial artificial chromosome clones through multiplexed fluorescence-labeled fingerprinting

A rapid multiplexed fingerprinting method has been developed for bacterial artificial chromosome (BAC) contig assembly. Defined subsets of BAC DNA fragments that result from digestion by three paired restriction endonucleases are labeled with unique fluorescent F-ddATP for each subset. Lists of the labeled fragment size are generated by an ABI 377 DNA sequencer and the GeneScan analysis software and then processed by an assembly program, FPC (Fingerprinted Contigs), to produce contig maps. Data obtained from the multiplexed labeling permit detection of smaller overlaps than is observed when data from a single double-digest are analyzed. The method has been tested on 98 BACs from chromosome 22 regions where large-scale sequencing is under way and also through simulation, using randomly generated BAC clones derived from existing DNA sequence data. In each case, contig assembly results demonstrated the advantages of multiplexed fingerprinting.

Integrated mapping package--a physical mapping software tool kit

We have developed an integrated physical mapping computer software package (IMP), originally designed to support the physical mapping of human chromosome 13 and expanded to support several gene-identification projects based on the positional candidate approach. IMP displays map data in a form that provides useful guidelines to the end users. An integrated map with high resolution and confidence is constructed from different types of mapping data, including hybridization experiments, STS-based PCR assays, genetic linkage mapping, cDNA localization, and FISH data. The map is also designed to provide suggestions for specific experiments that are required to obtain maps with even higher resolution and confidence. To this end, the optimization employs multiple constraints that take into account already established STS "scaffold" maps. This software thus serves as an important general tool kit for physical mapping, sequencing, and gene-hunting projects.

Allelotyping demonstrates common and distinct patterns of chromosomal loss in human lung cancer types

Allelic loss is a hallmark of tumor suppressor gene (TSG) inactivation. We have allelotyped 29 paired lymphoblastoid and lung cancer cell lines derived from 11 patients with small cell (SCLC) and 18 patients with non-small cell lung carcinomas (NSCLC). Statistical analysis indicated that a threshold of 30% separated non-random allelic loss from the random genetic deletions of malignancy. We have identified non-random allelic loss at 42 of 54 (78%) specific chromosomal regions examined, with 22 regions (52%) common between the two major lung cancer histologic types. There were 3 regions (7%) with allelic loss specific for SCLC and 17 regions (41%) specific for NSCLC. Furthermore, there were significant differences in loss of heterozygosity (LOH) frequencies between NSCLC and SCLC at 13 regions on eight chromosome arms (3p, 5q, 6q, 9p, 10q, 11p, 13q, and 19p). Eight homozygous deletions were present in seven cell lines at four regions, 3p12, 3p14.2, 9p21, and 10q23-25. We have also identified novel sites of chromosomal deletions. In particular, there was frequent loss at 11p13 in SCLC and loss at 6p21.3 and 13q12.3 in NSCLC. In this study, we demonstrate that a) non-random allelic losses in lung cancer involve multiple regions; b) some losses are common to both NSCLC and SCLC subtypes, whereas others are subtype specific; c) there are genetic deletions at novel chromosomal regions; and d) several homozygous deletions have been noted. Our studies demonstrate the usefulness of continuous cell lines for detailed allelotyping, for comparing genetic abnormalities between SCLC and NSCLC, and for identifying homozygous deletions.

High throughput fingerprint analysis of large-insert clones

As part of the Human Genome Project, the Washington University Genome Sequencing Center has commenced systematic sequencing of human chromsome 7. To organize and supply the effort, we have undertaken the construction of sequence-ready physical maps for defined chromosomal intervals. Map construction is a serial process composed of three main activities. First, candidate STS-positive large-insert PAC and BAC clones are identified. Next, these candidate clones are subjected to fingerprint analysis. Finally, the fingerprint data are used to assemble sequence-ready maps. The fingerprinting method we have devised is key to the success of the overall approach. We present here the details of the method and show that the fingerprints are of sufficient quality to permit the construction of megabase-size contigs in defined regions of the human genome. We anticipate that the high throughput and precision characteristic of our fingerprinting method will make it of general utility.

A sequence-ready physical map of a region of 12q24.1

We developed a sequence-ready map of a part of human chromosome 12q24.1. We utilized a number of sequence-tagged site (STS) markers from 12q24.1 to screen large insert bacterial chromosome libraries and a chromosome 12-specific cosmid library. The clones were assembled into contiguous sets (contigs) by STS-content analysis. Contigs were extended by obtaining end sequences of bacterial clones, generation of additional STSs, rescreening the libraries, and screening the additional clones for the presence of STSs. The resulting contig covers nearly 2 Mb of DNA and provides an average marker resolution of 16 kb. Based on the STS content, we developed fingerprints of a subset of clones. The STS content and fingerprint data allowed us to define a minimal tiling path of clones. These clones are being used to sequence this part of chromosome 12. This contig contains the Ataxin 2 gene, and it covers the interval harboring the gene responsible for Darier disease.

Construction of a 1-Mb restriction-mapped cosmid contig containing the candidate region for the familial Mediterranean fever locus (MEFV) on chromosome 16p 13.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, BAC, and YAC clones, we have generated a contig of genomic clones spanning approximately 1050 kb that contains the FMF critical region. The map consists of 179 cosmid, 15 P1, 10 PAC, 3 BAC, and 17 YAC clones, anchored by 27 STS markers. Eight additional STSs have been developed from the approximately 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.

Genomic organization and the 5' flanking region of the gamma subunit of the human amiloride-sensitive epithelial sodium channel

The amiloride-sensitive epithelial sodium channel (ENaC) complex is made up of at least three different subunits alpha, beta, and gamma, which are developmentally regulated, selectively expressed, and variously up-regulated by steroid hormones. To understand mechanisms involved in regulation of the gamma subunit, we have determined the structure of the human gammaENaC gene. By 5' rapid amplification of cDNA ends, primer extension analysis, and nuclease protection assay, we identified transcription start sites in human brain, kidney, and lung. A human genomic library was screened and overlapping cosmid clones that span approximately 50 kilobases and contain the hgammaENaC gene were identified. The 5'-untranslated region is 141 bases long, and the translation start codon is contained within the second exon. The human gene spans at least 35 kilobases. The 5' end of the gene including portions of 5' flanking genomic DNA and the first intron are G + C rich and contain several CpG dinucleotides, consistent with a CpG island. The 5' flanking region contains no CCAAT or TATA-like elements but does contain two GC boxes as well as several putative transcription factor binding sites including AP-2, Sp1, CRE, PEA-3, and NF-IL6. This is the first description of the structural organization and the 5' flanking region of a member of the epithelial sodium channel complex.

Random-breakage mapping method applied to human DNA sequences

The random-breakage mapping method [Game et al. (1990) Nucleic Acids Res., 18, 4453-4461] was applied to DNA sequences in human fibroblasts. The methodology involves NotI restriction endonuclease digestion of DNA from irradiated calls, followed by pulsed-field gel electrophoresis, Southern blotting and hybridization with DNA probes recognizing the single copy sequences of interest. The Southern blots show a band for the unbroken restriction fragments and a smear below this band due to radiation induced random breaks. This smear pattern contains two discontinuities in intensity at positions that correspond to the distance of the hybridization site to each end of the restriction fragment. By analyzing the positions of those discontinuities we confirmed the previously mapped position of the probe DXS1327 within a NotI fragment on the X chromosome, thus demonstrating the validity of the technique. We were also able to position the probes D21S1 and D21S15 with respect to the ends of their corresponding NotI fragments on chromosome 21. A third chromosome 21 probe, D21S11, has previously been reported to be close to D21S1, although an uncertainty about a second possible location existed. Since both probes D21S1 and D21S11 hybridized to a single NotI fragment and yielded a similar smear pattern, this uncertainty is removed by the random-breakage mapping method.

Mapping the human Y chromosome by fingerprinting cosmid clones

We have used Y-specific cosmid clones in a random fingerprinting approach to build contigs on the human Y chromosome. Clones derived from two libraries have been analyzed. The construction of one library is described here, the second was the Y chromosome-specific library LLOYNCO3 "M" (Lawrence Livermore National Laboratory). To date, we have fingerprinted 4430 cosmids: 377 contigs have been constructed containing from 2 to 39 clones. Along with the singletons, we estimate that we have covered 72.5% of the euchomatic portion of the Y chromosome with fingerprinted clones. Sequence tagged sites are being used to anchor cosmids and contigs onto the YAC framework.

Cloning and characterization of the majority of repetitive DNA in cotton (Gossypium L.)

Repetitive DNA elements representing 60-70% of the total repetitive DNA in tetraploid cotton (Gossypium barbadense L.) and comprising 30-36% of the tetraploid cotton genome were isolated from a genomic library of DNA digested with a mixture of four blunt-end cutting restriction enzymes. A total of 313 clones putatively containing nuclear repetitive sequences were classified into 1103 families, based on cross hybridization and Southern blot analysis. The 103 families were characterized in terms of genome organization, methylation pattern, abundance, and DNA variation. As in many other eukaryotic genomes, interspersed repetitive elements are the most abundant class of repetitive DNA in the cotton genome. Paucity of tandem repeat families with high copy numbers (>10(4)) may be a unique feature of the cotton genome as compared with other higher plant genomes. Interspersed repeats tend to be methylated, while tandem repeats seem to be largely unmethylated in the cotton genome. Minimal variation in repertoire and overall copy number of repetitive DNA elements among different tetraploid cotton species is consistent with the hypothesis of a relatively recent origin of tetraploid cottons.

Repetitive sequence fingerprinting in the long range mapping of mammalian genomes

This review presents some properties of interspersed repeats, particularly human and mouse repeats, and shows how these have been utilized in long-range genome mapping. The link between the distribution of such repeats and their relationship with genome organization is discussed.

Structure of the human CRFB4 gene: comparison with its IFNAR neighbor

The cytokine receptor family consists of a growing number of structurally and evolutionarily related transmembrane receptors. CRFB4 and IFNAR are two of the most similar members of this family. They are encoded by two neighboring genes on both human chromosome 21 and murine chromosome 16. The sequence of the human CRFB4 gene was determined from the first exon to the last intron. The nature of the repetitive sequences present in the introns was analyzed and compared with those present in the human IFNAR gene. This analysis leads to considerations of the antiquity of the duplication that gave rise to both genes from a common ancestor. A pseudogene for USF has been identified in the IFNAR gene and a new definition for the repetitive sequence MER37 is proposed. The polymorphism associated with two CA repeats present in the CRFB4 gene is described.

Phenol sulfotransferases: candidate genes for Batten disease

Batten disease (juvenile-onset neuronal ceroid lipofuscinosis; JNCL) is an autosomal recessive neurodegenerative disorder, characterized by the cytosomal accumulation of autofluorescent proteolipopigments in neurons and other cell types. The Batten disease gene (CLN3) has not yet been identified, but has been mapped to a small region of human chromosome area 16p12.1-p11.2. We recently reported the fortuitous discovery that the cytosolic phenol sulfotransferase gene (STP) is located within this same interval of chromosome 16p. Since phenol sulfotransferase is expressed in neurons, can sulfate lipophilic phenolic compounds, and is mapped near CLN3, STP is considered as a candidate gene for Batten disease. YAC and cosmid cloning results have further substantiated the close proximity of STP and a highly related sulfotransferase (STM), encoding the catecholamine-preferring enzyme, to the CLN3 region of chromosome 16p. In this report, we summarize some of the recent progress in the identification of two phenol sulfotransferase genes (STP and STM) as positional candidate genes for Batten disease.

Ordered restriction endonuclease maps of yeast artificial chromosomes created by optical mapping on surfaces

We have developed a surface mounting technology for the rapid construction of ordered restriction maps from individual DNA molecules. Optical restriction maps constructed from yeast artificial chromosome DNA molecules mounted on specially derivatized glass surfaces are accurate and reproducible, and the technology is amenable to automation. The mounting procedures described here should also be useful for fluorescence in situ hybridization studies. We believe these improvements to optical mapping will further stimulate the development of nonelectrophoretic approaches to genome analysis.

Primary sequence, copy number, and distribution of mariner transposons in the honey bee

A single honey bee mariner transposon (TnM1a) was sequenced, revealing a transpositionally non-autonomous element of 937 bp delimited by 30 bp perfect inverted terminal repeats. The element is flanked by the TA duplication typical of mariner elements in general. There are approximately 435 copies of TnM1a homologous elements per haploid genome. These elements appear, by Southern blot analysis, to be dispersed throughout the genome. Thirteen individual genomic clones with an average size of 15 kb, were found to contain only a single element each, which also suggests that the elements are not tightly clustered. Finally, mariner elements are neither inactivated by methylation nor sequestered into a methylated fraction of the genome.

High-resolution physical mapping of a 250-kb region of human chromosome 11q24 by genomic sequence sampling (GSS)

A physical map of the region of human chromosome 11q24 containing the FLI1 gene, disrupted by the t(11;22) translocation in Ewing sarcoma and primitive neuroectodermal tumors, was analyzed by genomic sequence sampling. Using a 4- to 5-fold coverage chromosome 11-specific library, 22 region-specific cosmid clones were identified by phenol emulsion reassociation hybridization, with a 245-kb yeast artificial chromosome clone containing the FLI1 gene, and by directed "walking" techniques. Cosmid contigs were constructed by individual clone fingerprinting using restriction enzyme digestion and assembly with the Genome Reconstruction and AsseMbly (GRAM) computer algorithm. The relative orientation and spacing of cosmid contigs with respect to the chromosome was determined by the structural analysis of cosmid clones and by direct visual in situ hybridization mapping. Each cosmid clone in the contig was subjected to "one-pass" end sequencing, and the resulting ordered sequence fragments represent approximately 5% of the complete DNA sequence, making the entire region accessible by PCR amplification. The sequence samples were analyzed for putative exons, repetitive DNAs, and simple sequence repeats using a variety of computer algorithms. Based upon the computer predictions, Southern and Northern blot experiments led to the independent identification and localization of the FLI1 gene as well as a previously unknown gene located in this region of chromosome 11q24. This approach to high-resolution physical analysis of human chromosomes allows the assembly of detailed sequence-based maps and provides a tool for further structural and functional analysis of the genome.

Chromosome 16 microdeletion in a patient with juvenile neuronal ceroid lipofuscinosis (Batten disease)

The gene that is involved in juvenile neuronal ceroid lipofuscinosis (JNCL), or Batten disease--CLN3--has been localized to 16p12, and the mutation shows a strong association with alleles of microsatellite markers D16S298, D16S299, and D16S288. Recently, haplotype analysis of a Batten patient from a consanguineous relationship indicated homozygosity for a D16S298 null allele. PCR analysis with different primers on DNA from the patient and his family suggests the presence of a cytogenetically undetectable deletion, which was confirmed by Southern blot analysis. The microdeletion is embedded in a region containing chromosome 16-specific repeated sequences. However, putative candidates for CLN3, members of the highly homologous sulfotransferase gene family, which are also present in this region in several copies, were not deleted in the patient. If the microdeletion in this patient is responsible for Batten disease, then we conclude that the sulfotransferase genes are probably not involved in JNCL. By use of markers and probes flanking D16S298, the maximum size of the microdeletion was determined to be approximately 29 kb. The microdeletion may affect the CLN3 gene, which is expected to be in close proximity to D16S298.

Fragment collapsing and splitting while assembling high-resolution restriction maps

In the process of constructing high-resolution restriction maps via greedy algorithms, a classical anomaly, known as fragment collapsing, introduces errors into the maps that impedes further map assembly. Fragment collapsing occurs when two different genomic fragments of approximately the same length and occurring in the digestion of two different overlapping clones are incorrectly identified as representative of a single genomic fragment. This introduces a single fragment of commensurate length into an incorrect position in the map. The present work describes techniques for detecting and correcting such anomalies.

Mapping the clonally unstable recombinogenic PRM1-->PRM2-->TNP2 region of human 16p13.2

We have isolated a human cosmid clone encompassing approximately 40 kb of the chromosome 16p13.2 region. This region contains the protamine genes PRM1 and PRM2 and the transition protein gene TNP2. To facilitate the characterization of this under represented region of the chromosome 16 physical map, a detailed map of this cosmid clone was constructed. Analysis revealed that this gene cluster was nestled amid numerous repetitive elements. Under representation of this segment may be attributed to these repetitive sequence elements scattered throughout this region rendering this locus clonally unstable. Verification of the map by cosmid and genomic Southern analysis resolved discrepancies within the literature and showed that this highly repetitive DNA sequence containing cosmid was representative of this region of the human genome. This dynamic and recombinogenic region continues to evolve.

Cosmid assembly and anchoring to human chromosome 21

A human chromosome 21-specific cosmid library from the Lawrence Livermore National Laboratory has been analyzed by two complementary methods, fingerprinting and hybridization; 40% coverage of the entire chromosome 21 has been achieved. To prepare a contig pool, approximately 9300 cosmid clones randomly selected from the library were fingerprinted and automatically assembled into 467 overlapping sets by the fluorescence-tagged restriction fragment method. The average size of the overlapping sets was 9.5 cosmids with minimal tiling paths consisting of 5.4 cosmids with a 10-kb extension each. However, as many as 10% of overlaps within members were estimated to be false. For regional localization, we hybridized gridded arrays of cosmids with inter-Alu-PCR probes obtained from YAC clones and somatic cell hybrids and assigned 592 cosmids to 26 subregions of 21q. Of these, 371 clones were incorporated into 139 contigs, anchoring the total 1864 cosmids to the subregion. The remaining 221 clones were mapped as orphans. To correlate the cytogenetic, YAC, and cosmid maps on 21q, the translocation breakpoints of the chromosomes contained in the somatic cell hybrids were mapped with respect to the STS content of the YACs. From the gene cluster regions, 176 ribosomal and 25 alphoid clones were isolated by hybridization. Together, these sets of anchored contigs and cosmids will provide a valuable resource for construction of a high-resolution map and for isolation of genes of interest from chromosome 21.

Exploiting the informativity of 'meaningless' simple repetitive DNA from indirect gene diagnosis to multilocus genome scanning

Most eukaryotic genomes are characterized by excessively large amounts of non-coding DNA sequences among which redundant (repetitive) elements constitute a sizable portion. The functional role of an abundant subclass of repetitive sequences--simple, tandemly arranged repeats--remained mysterious so far. Even the biological meaning of most of these elements appears quite refractory to present-day techniques in molecular genetics. Notwithstanding simple repetitive sequences have been developed into superb tools for various aspects of eukaryotic genome research: Using oligonucleotide probes carrying simple repeat motifs multilocus DNA fingerprinting can be applied for individual identification and genetic relationship analyses in plants, animals and humans. Microsatellite analyses via polymerase chain reaction of simple repeat blocks allow for efficient investigations of such divers subject matters as criminal stains, detailed genome maps and indirect gene diagnoses.

FISH mapping of a human chromosome 16 constitutional pericentric inversion inv(16)(p13q22) found in a large kindred

Fluorescence in situ hybridization analysis (FISH) was used to map the constitutional chromosome 16 pericentric inversion breakpoints inv(16)(p13q22) detected in one individual (II-2) from a large kindred [Bianchi et al., 1992: Am J Med Genet 43:791-795]. The breakpoints found in individual II-2 mapped to distinctly different locations than the chromosome 16 pericentric inversion breakpoints commonly acquired in acute nonlymphocytic leukemia. The constitutional pericentric inversion breakpoints also do not map to regions where low abundance repetitive DNA sequences found in bands 16p13 and q22 are located. The results indicate that low abundance, chromosome 16-specific repetitive DNA sequences in bands p13 and q22 are probably not causally related to the inversion that is found in many members of a large kindred [Bianchi et al., 1992].

Design and analysis of chromosome physical mapping experiments

Mathematical and statistical aspects of constructing ordered-clone physical maps of chromosomes are reviewed. Three broad problems are addressed: analysis of fingerprint data to identify configurations of overlapping clones, prediction of the rate of progress of a mapping strategy and optimal design of pooling schemes for screening large clone libraries.

Evidence of linkage disequilibrium in the Spanish polycystic kidney disease I population

Forty-one Spanish families with polycystic kidney disease 1 (PKD1) were studied for evidence of linkage disequilibrium between the disease locus and six closely linked markers. Four of these loci--three highly polymorphic microsatellites (SM6, CW3, and CW2) and an RFLP marker (BLu24)--are described for the first time in this report. Overall the results reveal many different haplotypes on the disease-carrying chromosome, suggesting a variety of independent PKD1 mutations. However, linkage disequilibrium was found between BLu24 and PKD1, and this was corroborated by haplotype analysis including the microsatellite polymorphisms. From this analysis a group of closely related haplotypes, consisting of four markers, was found on 40% of PKD1 chromosomes, although markers flanking this homogeneous region showed greater variability. This study has highlighted an interesting subpopulation of Spanish PKD1 chromosomes, many of which have a common origin, that may be useful for localizing the PKD1 locus more precisely.

ISWAC: proposed system for the integrated assembly of chromosomes

The generation of a physical map as an integral part of sequence project management is a problem that present computer systems do not address. Primarily, the analysis performed is based solely on the information available from a single knowledge level. Management systems that are currently available do not adequately model the multi-layer top down strategy that is most often utilized to manage large scale sequencing projects. Single layered approaches reflect an algorithmic inadequacy since interacting data sets are required to provide a good solution. The analysis tool that is currently under development termed ISWAC, the Integrated System for Wholistic Assembly of Chromosomes, overcomes these limitations by integrating information available from five layers of knowledge. These knowledge layers utilize information from the linkage map, physical map, restriction map, clone strategy map and the DNA sequence itself. The approach we are implementing, reviews current project status and continually refines the experimental strategy necessary to efficiently complete the sequencing task. To facilitate project completion the system is designed to interactively recommend strategies based on partial information. The utility of this tool is enhanced by implementing knowledge representation techniques that allow reasoning with approximate concepts characteristic of these data-sets. In addition, the raw physical data is maintained within an integrated map database to ease data verification. This paper presents the first discussion of the design specifications for a computer system to assimilate the various forms of data that are being generated as part of the human genome project. It was specifically written to stimulate discussion regarding data standardization, translation, analysis and most important, an understandable user-interphase for the molecular biologist. We would hope that interested readers would respond by assisting in the definition of a set of universal data standards and adopting them in their laboratories.

A new bacteriophage P1-derived vector for the propagation of large human DNA fragments

We have designed a P1 vector (pCYPAC-1) for the introduction of recombinant DNA into E. coli using electroporation procedures. The new cloning system, P1-derived artificial chromosomes (PACs), was used to establish an initial 15,000 clone library with an average insert size of 130-150 kilobase pairs (kb). No chimaerism has been observed in 34 clones, by fluorescence in situ hybridization. Similarly, no insert instability has been observed after extended culturing, for 20 clones. We conclude that the PAC cloning system will be useful in the mapping and detailed analysis of complex genomes.

Detection of cryptic chromosomal abnormalities in unexplained mental retardation: a general strategy using hypervariable subtelomeric DNA polymorphisms

Given the availability of DNA from both parents, unusual segregation of hypervariable DNA polymorphisms (HVPs) in the offspring may be attributable to deletion, unbalanced chromosomal translocation, or uniparental disomy. The telomeric regions of chromosomes are rich in both genes and hypervariable minisatellite sequences and may also be particularly prone to cryptic breakage events. Here I describe and analyze a general approach to the detection of subtelomeric abnormalities and uniparental disomy in patients with unexplained mental retardation. With 29 available polymorphic systems, approximately 50%-70% of these abnormalities could currently be detected. Development of subtelomeric HVPs physically localized with respect to their telomeres should provide a valuable resource in routine diagnostics.

Use of reference libraries and hybridisation fingerprinting for relational genome analysis

The concept of relational genome analysis by hybridisation has been developed into a working system. Various genomic and cDNA libraries have been generated and are distributed via a reference system. Analysis procedures have been tested successfully in the mapping of the entire Schizosaccharomyces pombe genome. In another test-case for their refinement, analyses on the Drosophila genome are well under way. Human and mouse libraries are being studied on all levels, from generating YAC maps to partially sequencing representative cDNA libraries. The automation of the involved processes and the development of improved image detection and analysis are well advanced.

The human genome project

The Human Genome Project in the United States is now well underway. Its programmatic direction was largely set by a National Research Council report issued in 1988. The broad framework supplied by this report has survived almost unchanged despite an upheaval in the technology of genome analysis. This upheaval has primarily affected physical and genetic mapping, the two dominant activities in the present phase of the project. Advances in mapping techniques have allowed good progress toward the specific goals of the project and are also providing strong corollary benefits throughout biomedical research. Actual DNA sequencing of the genomes of the human and model organisms is still at an early stage. There has been little progress in the intrinsic efficiency of DNA-sequence determination. However, refinements in experimental protocols, instrumentation, and project management have made it practical to acquire sequence data on an enlarged scale. It is also increasingly apparent that DNA-sequence data provide a potent means of relating knowledge gained from the study of model organisms to human biology. There is as yet little indication that the infusion of technology from outside biology into the Human Genome Project has been effectively stimulated. Opportunities in this area remain large, posing substantial technical and policy challenges.

Low-frequency chimeric yeast artificial chromosome libraries from flow-sorted human chromosomes 16 and 21

Construction of chromosome-specific yeast artificial chromosome (YAC) libraries from sorted chromosomes was undertaken (i) to eliminate drawbacks associated with first-generation total genomic YAC libraries, such as the high frequency of chimeric YACs, and (ii) to provide an alternative method for generating chromosome-specific YAC libraries in addition to isolating such collections from a total genomic library. Chromosome-specific YAC libraries highly enriched for human chromosomes 16 and 21 were constructed. By maximizing the percentage of fragments with two ligatable ends and performing yeast transformations with less than saturating amounts of DNA in the presence of carrier DNA, YAC libraries with a low percentage of chimeric clones were obtained. The smaller number of YAC clones in these chromosome-specific libraries reduces the effort involved in PCR-based screening and allows hybridization methods to be a manageable screening approach.

Genome mapping by nonrandom anchoring: a discrete theoretical analysis

As part of our effort to construct a physical map of the genome of the fission yeast Schizosaccharomyces pombe, we have made theoretical predictions for the progress expected, as measured by the expected length fraction of island coverage and by the expected properties of the anchored islands such as the number and the size of islands. Our experimental strategy is to construct a random clone library and screen the library for clones having unique sequence at both ends. This scheme is essentially the same as the clone-limited double sequence-tagged-site selection scheme which was used in a computer simulation by Palazzolo et al. [Palazzolo, M. J., Sawyer, S. A., Martin, C. H., Smoller, D. A. & Hartl, D. L. (1991) Proc. Natl. Acad. Sci. USA 88, 8034-8038]. Both simulation and ongoing experiments in our laboratory have shown that the nonrandom anchoring method is far superior to random anchoring. In this paper, we propose a theoretical model to explain the simulated data and the experimental data.