A physical map of human chromosome 7: an integrated YAC contig map with average STS spacing of 79 kb

Reconstructing the demographic history of the human lineage using whole-genome sequences from human and three great apes

The demographic history of human would provide helpful information for identifying the evolutionary events that shaped the humanity but remains controversial even in the genomic era. To settle the controversies, we inferred the speciation times (T) and ancestral population sizes (N) in the lineage leading to human and great apes based on whole-genome alignment. A coalescence simulation determined the sizes of alignment blocks and intervals between them required to obtain recombination-free blocks with a high frequency. This simulation revealed that the size of the block strongly affects the parameter inference, indicating that recombination is an important factor for achieving optimum parameter inference. From the whole genome alignments (1.9 giga-bases) of human (H), chimpanzee (C), gorilla (G), and orangutan, 100-bp alignment blocks separated by ≥5-kb intervals were sampled and subjected to estimate τ = μT and θ = 4μgN using the Markov chain Monte Carlo method, where μ is the mutation rate and g is the generation time. Although the estimated τ_HC differed across chromosomes, τ_HC and τ_HCG were strongly correlated across chromosomes, indicating that variation in τ is subject to variation in μ, rather than T, and thus, all chromosomes share a single speciation time. Subsequently, we estimated Ts of the human lineage from chimpanzee, gorilla, and orangutan to be 6.0–7.6, 7.6–9.7, and 15–19 Ma, respectively, assuming variable μ across lineages and chromosomes. These speciation times were consistent with the fossil records. We conclude that the speciation times in our recombination-free analysis would be conclusive and the speciation between human and chimpanzee was a single event.

Duplication of C7orf58, WNT16 and FAM3C in an obese female with a t(7;22)(q32.1;q11.2) chromosomal translocation and clinical features resembling Coffin-Siris Syndrome

We characterized the t(7;22)(q32;q11.2) chromosomal translocation in an obese female with coarse features, short stature, developmental delay and a hypoplastic fifth digit. While these clinical features suggest Coffin-Siris Syndrome (CSS), we excluded a CSS diagnosis by exome sequencing based on the absence of deleterious mutations in six chromatin-remodeling genes recently shown to cause CSS. Thus, molecular characterization of her translocation could delineate genes that underlie other syndromes resembling CSS. Comparative genomic hybridization microarrays revealed on chromosome 7 the duplication of a 434,682 bp region that included the tail end of an uncharacterized gene termed C7orf58 (also called CPED1) and spanned the entire WNT16 and FAM3C genes. Because the translocation breakpoint on chromosome 22 did not disrupt any apparent gene, her disorder was deemed to result from the rearrangement on chromosome 7. Mapping of yeast and bacterial artificial chromosome clones by fluorescent in situ hybridization on chromosome spreads from this patient showed that the duplicated region and all three genes within it were located on both derivative chromosomes 7 and 22. Furthermore, DNA sequencing of exons and splice junctional regions from C7orf58, WNT16 and FAM3C revealed the presence of potential splice site and promoter mutations, thereby augmenting the detrimental effect of the duplicated genes. Hence, dysregulation and/or disruptions of C7orf58, WNT16 and FAM3C underlie the phenotype of this patient, serve as candidate genes for other individuals with similar clinical features and could provide insights into the physiological role of the novel gene C7orf58.

Isolation of single, intact chromosomes from single, selected ovarian cancer cells for in situ hybridization and sequencing

The first step towards effective therapy of cancer is to reveal molecular profiles of all cell clones propelling tumor growth. The specific aim of this project was to develop a technology helping us to isolate patient's single, living cells based upon their cancer-specific, cell surface biomarkers, to reveal their molecular profiles, and to isolate, from these selected cells, intact chromosomes for in situ hybridization (FISH) and for next generation sequencing (NGS). We attained this aim, while probing the cells from the ovarian cancer patients. Ovarian cancer is the most deadly of all gynecological cancers. In most of the patients with the advanced stages of this cancer, the gene for epidermal growth factors receptor (EGFR) is mutated, as the deletion variant III, resulting in the truncated transcripts and products. From these patients, we collected cells from peritoneal fluid, blood, lymph, and biopsies. We genetically engineered fluorescent and superparamagnetic single chain variable fragments (scFvs) targeting EGFRwt and EGFRvIII. Using these scFvs, we isolated single, living ovarian cancer cells and analyzed their transcripts and products. We further genetically engineered scFv targeting dsDNA. Using these scFvs, we isolated the entire single, intact chromosomes from the selected, single ovarian cancer cells for NGS and for liquid phase FISH. This novel work-flow opens new routes not only for molecular profiling of the entire spectrum of cancer cell clones in the diagnosed patient, one cell clone at a time, but also for manufacturing targeted contrast for in vivo imaging and for designing and guiding targeted delivery of therapeutic genes in cancer therapy.

An autosomal analysis gives no genetic evidence for complex speciation of humans and chimpanzees

There have been conflicting arguments as to what happened in the human-chimpanzee speciation event. Patterson et al. (2006, Genetic evidence for complex speciation of humans and chimpanzees. Nature 441:1103-1108) proposed a hypothesis that the human-chimpanzee speciation event involved a complicated demographic process: that is, the ancestral lineages of humans and chimpanzees experienced temporal isolation followed by a hybridization event. This hypothesis stemmed from two major observations: a wide range of human-chimpanzee nucleotide divergence across the autosomal genome and very low divergence in the X chromosome. In contrast, Innan and Watanabe (2006, The effect of gene flow on the coalescent time in the human-chimpanzee ancestral population. Mol Biol Evol. 23:1040-1047) demonstrated that the null model of instantaneous speciation fits the genome-wide divergence data for the two species better than alternative models involving partial isolation and migration. To reconcile these two conflicting reports, we first reexamined the analysis of autosomal data by Patterson et al. (2006). By providing a theoretical framework for their analysis, we demonstrated that their observation is what is theoretically expected under the null model of instantaneous speciation with a large ancestral population. Our analysis indicated that the observed wide range of autosomal divergence is simply due to the coalescent process in the large ancestral population of the two species. To further verify this, we developed a maximum likelihood function to detect evidence of hybridization in genome-wide divergence data. Again, the null model with no hybridization best fits the data. We conclude that the simplest speciation model with instantaneous split adequately describes the human-chimpanzee speciation event, and there is no strong reason to involve complicated factors in explaining the autosomal data.

Sequencing and Analyzing the t(1;7) Reciprocal Translocation Breakpoints Associated with a Case of Childhood-onset Schizophrenia/Autistic Disorder

We characterized a t(1;7)(p22;q21) reciprocal translocation in a patient with childhood-onset schizophrenia (COS) and autism using genome mapping and sequencing methods. Based on genomic maps of human chromosome 7 and fluorescence in situ hybridization (FISH) studies, we delimited the region of 7q21 harboring the translocation breakpoint to a approximately 16-kb interval. A cosmid containing the translocation-associated 1:7 junction on der(1) was isolated and sequenced, revealing the positions on chromosomes 1 and 7, respectively, where the translocation occurred. PCR-based studies enabled the isolation and sequencing of the reciprocal 7:1 junction on der(7). No currently recognized gene on either chromosome appears to be disrupted by the translocation. We further found no evidence for copy-number differences in the genomic regions flanking the translocation junctions in the patient. Our efforts provide sequence-based information about a schizophrenia/autism-associated translocation, and may facilitate future studies investigating the genetic bases of these disorders.

Estimation of DNA sequence context-dependent mutation rates using primate genomic sequences

It is understood that DNA and amino acid substitution rates are highly sequence context-dependent, e.g., C --> T substitutions in vertebrates may occur much more frequently at CpG sites and that cysteine substitution rates may depend on support of the context for participation in a disulfide bond. Furthermore, many applications rely on quantitative models of nucleotide or amino acid substitution, including phylogenetic inference and identification of amino acid sequence positions involved in functional specificity. We describe quantification of the context dependence of nucleotide substitution rates using baboon, chimpanzee, and human genomic sequence data generated by the NISC Comparative Sequencing Program. Relative mutation rates are reported for the 96 classes of mutations of the form 5' alphabetagamma 3' --> 5' alphadeltagamma 3', where alpha, beta, gamma, and delta are nucleotides and beta not equal delta, based on maximum likelihood calculations. Our results confirm that C --> T substitutions are enhanced at CpG sites compared with other transitions, relatively independent of the identity of the preceding nucleotide. While, as expected, transitions generally occur more frequently than transversions, we find that the most frequent transversions involve the C at CpG sites (CpG transversions) and that their rate is comparable to the rate of transitions at non-CpG sites. A four-class model of the rates of context-dependent evolution of primate DNA sequences, CpG transitions > non-CpG transitions approximately CpG transversions > non-CpG transversions, captures qualitative features of the mutation spectrum. We find that despite qualitative similarity of mutation rates among different genomic regions, there are statistically significant differences.

A composite-likelihood approach for detecting directional selection from DNA sequence data

We present a novel composite-likelihood-ratio test (CLRT) for detecting genes and genomic regions that are subject to recurrent natural selection (either positive or negative). The method uses the likelihood functions of Hartl et al. (1994) for inference in a Wright-Fisher genic selection model and corrects for nonindependence among sites by application of coalescent simulations with recombination. Here, we (1) characterize the distribution of the CLRT statistic (Lambda) as a function of the population recombination rate (R=4Ner); (2) explore the effects of bias in estimation of R on the size (type I error) of the CLRT; (3) explore the robustness of the model to population growth, bottlenecks, and migration; (4) explore the power of the CLRT under varying levels of mutation, selection, and recombination; (5) explore the discriminatory power of the test in distinguishing negative selection from population growth; and (6) evaluate the performance of maximum composite-likelihood estimation (MCLE) of the selection coefficient. We find that the test has excellent power to detect weak negative selection and moderate power to detect positive selection. Moreover, the test is quite robust to bias in the estimate of local recombination rate, but not to certain demographic scenarios such as population growth or a recent bottleneck. Last, we demonstrate that the MCLE of the selection parameter has little bias for weak negative selection and has downward bias for positively selected mutations.

The DNA sequence of human chromosome 7

Human chromosome 7 has historically received prominent attention in the human genetics community, primarily related to the search for the cystic fibrosis gene and the frequent cytogenetic changes associated with various forms of cancer. Here we present more than 153 million base pairs representing 99.4% of the euchromatic sequence of chromosome 7, the first metacentric chromosome completed so far. The sequence has excellent concordance with previously established physical and genetic maps, and it exhibits an unusual amount of segmentally duplicated sequence (8.2%), with marked differences between the two arms. Our initial analyses have identified 1,150 protein-coding genes, 605 of which have been confirmed by complementary DNA sequences, and an additional 941 pseudogenes. Of genes confirmed by transcript sequences, some are polymorphic for mutations that disrupt the reading frame.

Molecular characterization of partial trisomy 16q24.1-qter: clinical report and review of the literature

We describe a 3(1/2)-year-old girl with psychomotor and mental retardation; dysmorphic features, including a high forehead with bitemporal narrowing; a broad nasal bridge and a broadened nose; downslanting palpebral fissures; abnormal ears; vertebral abnormalities; cardiac defect; genital hypoplasia; and anal abnormalities. The karyotype of our patient (550 bands) was normal. Molecular cytogenetic techniques, including comparative genomic hybridization (CGH) and fluorescence in situ hybridization (FISH), revealed that this girl was a carrier of a de novo derivative chromosome 7 arising from a cryptic t(7;16)(p22.3;q24.1) translocation generating a trisomy 16q24.1-qter and a 7p22.3-pter deletion. FISH with a series of specific chromosome 7p and 16q probes allowed us to delineate the chromosome 7 breakpoint between YAC660G6 (WD7S517) and YAC848A12 (D7S521, D7S31, and WI-4829) and the chromosome 16 breakpoint between BAC457K7 (D42053) and BAC44201 (SGC30711). The comparison of the clinical features of our patient with those of 2 cases of pure terminal 7p deletion and 28 cases of trisomy 16q reported in the literature allowed us to establish the following phenotype-genotype correlation for trisomy of the long arm of chromosome 16: distinctive facies (high/prominent forehead, bitemporal narrowing, periorbital edema in the neonatal period); severe mental retardation; vertebral, genital, and anal abnormalities to 16q24; distal joint contractures and camptodactyly to 16q23; cleft palate and renal anomalies to 16q22; beaked nose and gall bladder agenesis to 16q21; gut malrotation; lung and liver anomalies to 16q13; and behavior abnormalities to band 16q11-q13.

Fine mapping of the locus for Shwachman-Diamond syndrome at 7q11, identification of shared disease haplotypes, and exclusion of TPST1 as a candidate gene

Shwachman-Diamond syndrome (SDS) is an autosomal recessive disorder characterised by exocrine pancreatic dysfunction, haematological and skeletal abnormalities. We have previously defined the SDS locus as a 2.7 cM interval spanning the centromere of chromosome 7. To facilitate additional analysis of this complex and poorly characterised region, a framework of ordered genetic markers at 7p11-q11, including six newly identified, has been constructed using somatic cell, radiation hybrid and STS-content mapping. We have identified shared disease haplotypes, that recur in unrelated families of common ethnic origin, and extend across the SDS locus. Detection of ancestral and intrafamilial recombination events in patients refined the SDS locus to a 1.9 cM interval at 7q11, which contains the tyrosylprotein sulfotransferase 1 (TPST1) gene. Patients with SDS were screened for mutations in TPST1 by sequencing of exons and intron-exon junctions. Two single nucleotide polymorphisms, but no disease-causing mutations, were identified. In addition, Southern blot analysis yielded no evidence of large-scale mutations, and RT-PCR analysis failed to detect alterations in expression. These results exclude TPST1 as the causative gene for SDS. The established map of the refined SDS locus will assist in the identification and characterisation of other candidate genes for SDS.

Vertebrate genome sequencing: building a backbone for comparative genomics

The human genome sequence provides a reference point from which we can compare ourselves with other organisms. Interspecies comparison is a powerful tool for inferring function from genomic sequence and could ultimately lead to the discovery of what makes humans unique. To date, most comparative sequencing has focused on pair-wise comparisons between human and a limited number of other vertebrates, such as mouse. Targeted approaches now exist for mapping and sequencing vertebrate bacterial artificial chromosomes (BACs) from numerous species, allowing rapid and detailed molecular and phylogenetic investigation of multi-megabase loci. Such targeted sequencing is complementary to current whole-genome sequencing projects, and would benefit greatly from the creation of BAC libraries from a diverse range of vertebrates.

A balanced reciprocal translocation t(5;7)(q14;q32) associated with autistic disorder: molecular analysis of the chromosome 7 breakpoint

Autism is a neuropsychiatric disorder characterized by impairments in social interaction, restricted and stereotypic pattern of interest with onset by 3 years of age. The results of genetic linkage studied for autistic disorder (AD) have suggested a susceptibility locus for the disease on the long arm of chromosome 7. We report a girl with AD and a balanced reciprocal translocation t(5;7)(q14;q32). The mother carries the translocation but do not express the disease. Fluorescent in situ hybridization (FISH) analysis with chromosome 7-specific YAC clones showed that the breakpoint coincides with the candidate region for AD. We identified a PAC clone that spans the translocation breakpoint and the breakpoint was mapped to a 2 kb region. Mutation screening of the genes SSBP and T2R3 located just centromeric to the breakpoint was performed in a set of 29 unrelated autistic sibling pairs who shared at least one chromosome 7 haplotype. We found no sequence variations, which predict amino acid alterations. Two single nucleotide polymorphisms were identified in the T2R3 gene, and associations between allele variants and AD in our population were not found. The methylation pattern of different chromosome 7 regions in the patient's genomic DNA appears normal. Here we report the clinical presentation of the patient with AD and the characterization of the genomic organization across the breakpoint at 7q32. The precise localization of the breakpoint on 7q32 may be relevant for further linkage studies and molecular analysis of AD in this region.

Candidate gene isolation and comparative analysis of a commonly deleted segment of 7q22 implicated in myeloid malignancies

Monosomy 7 and deletion of 7q are recurring abnormalities in malignant myeloid diseases. Here we extensively characterize an approximately 2-Mb commonly deleted segment (CDS) of 7q22 bounded by D7S1503 and D7S1841. Approximately 1.8 Mb of sequence have been generated from this interval, facilitating the construction of a transcript map that includes large numbers of genes and ESTs. The intron/exon organization of seven genes and expression patterns of three genes were determined, and leukemia samples were screened for mutations in five genes. We have used polymorphic markers from this region to examine leukemia cells for allelic loss within 7q22. Finally, we isolated mouse genomic clones orthologous to several of the characterized human genes. Fluorescence in situ hybridization studies using these clones indicate that a region of orthologous synteny lies on proximal mouse chromosome 5. These resources should greatly accelerate the pace of candidate gene discovery in this region.

Strategies for the systematic sequencing of complex genomes

Recent spectacular advances in the technologies and strategies for DNA sequencing have profoundly accelerated the detailed analysis of genomes from myriad organisms. The past few years alone have seen the publication of near-complete or draft versions of the genome sequence of several well-studied, multicellular organisms - most notably, the human. As well as providing data of fundamental biological significance, these landmark accomplishments have yielded important strategic insights that are guiding current and future genome-sequencing projects.

Familial chordoma, a tumor of notochordal remnants, is linked to chromosome 7q33

Chordoma is a rare tumor originating from notochordal remnants that is usually diagnosed during midlife. We performed a genomewide analysis for linkage in a family with 10 individuals affected by chordoma. The maximum two-point LOD score based on only the affected individuals was 2.21, at recombination fraction 0, at marker D7S2195 on chromosome 7q. Combined analysis of additional members of this family (11 affected individuals) and of two unrelated families (one with 2 affected individuals and the other with 3 affected individuals), with 20 markers on 7q, showed a maximum two-point LOD score of 4.05 at marker D7S500. Multipoint analysis based on only the affected individuals gave a maximum LOD score of 4.78, with an approximate 2-LOD support interval from marker D7S512 to marker D7S684. Haplotype analysis of the three families showed a minimal disease-gene region from D7S512 to D7S684, a distance of 11.1 cM and approximately 7.1 Mb. No loss of heterozygosity was found at markers D7S1804, D7S1824, and D7S2195 in four tumor samples from affected family members. These results map a locus for familial chordoma to 7q33. Further analysis of this region, to identify this gene, is ongoing.

Construction of a physical map of an autism susceptibility region in 7q32.3-q33

The fast evolving progress of the human genome mapping and sequencing efforts facilitate the detection of genes also for complex traits. We focus on the detection of susceptibility loci for autism, a prototypical pervasive developmental disorder. Five genome screens worldwide have identified several putative locations of susceptibility genes thus far, with the most common region on chromosome 7q. In order to identify new candidate genes for infantile autism we constructed a physical map of bacterial artificial chromosome, P1-derived artificial chromosome and yeast artificial chromosome clones of a 3 Mb region between D7S1575 and D7S500, including a complete contig of the approximately 1.2 Mb region around D7S2533, the marker with the most significant association result. We developed 16 novel sequence tag sites and mapped 23 genes/expressed sequence tags to the contigs. As this map contains a putative autistic disorder locus this integrated physical and transcript map provides a valuable resource for identification of candidate gene(s).

Deletions of PURA, at 5q31, and PURB, at 7p13, in myelodysplastic syndrome and progression to acute myelogenous leukemia

Deletions or monosomy of chromosomes 5 and 7 are frequently observed in myelodysplastic syndromes (MDS) and acute myelogenous leukemia (AML). In this study two genes, PURA and PURB, encoding functionally cooperative proteins in the Pur family, are localized to chromosome bands 5q31.1 and 7p13, respectively. One or both of these loci are shown to be hemizygously deleted in 60 MDS or AML patients using fluorescence in situ hybridization (FISH). High-resolution mapping of PURA localizes it approximately 1.1 Mb telomeric to the EGR-1 gene. Frequency of PURA deletion and segregation with EGR-1 indicate that PURA is within the most commonly deleted segment in myeloid disorders characterized by del(5)(q31). No mutations have been detected within the coding sequence of PURA. Concurrent deletions of PURA and PURB occur in MDS at a rate nearly 1.5-fold higher than statistically expected and in AML at a rate > 5-fold higher. This novel simultaneous deletion of two closely related gene family members may thus have consequences related to progression to AML. Pur alpha, an Rb-binding protein, has been implicated in cell cycle control and differentiation, and Pur alpha and Pur beta are reported to function as heterodimers. Alterations in these genes could affect a delicate balance critical in myeloid development.

Refinement within single yeast artificial chromosome clones of a minimal region commonly deleted on the short arm of chromosome 7 in Wilms tumours

Cytogenetic and molecular data indicate an involvement of genes mapped to the proximal portion of the short arm of chromosome 7 (7p) in Wilms tumours (WTs). We have analysed 38 WTs using a panel of eight microsatellite markers mapped to proximal 7p. Loss of heterozygosity (LOH) in tumour, compared with matched constitutional DNA, was identified in eight cases. To define better the minimal region commonly deleted in these tumours, they were analysed with nine additional markers, mapped within the region of interest. One tumour (case 30) showed LOH for only one marker (D7S510), while maintaining heterozygosity for the two immediately flanking loci (D7S555 and D7S668). This result was confirmed by fluorescence in situ hybridisation analysis, which showed that in the majority (65%) of nuclei from tumour 30 hybridising with a bacterial artificial chromosome clone containing the D7S510 locus, only one signal was visible. Noticeably, both markers defining the limits of the observed deleted region are simultaneously present within two distinct overlapping yeast artificial chromosome (YAC) clones mapped to chromosome bands 7p13-p14. This suggests that the maximum length of the missing DNA fragment was approximately 1.3 Mb, corresponding to the length of the smaller of the two YAC clones. In all other cases that showed LOH, the deletion encompassed the 7p13-p14 region. For this reason, we speculate that the identified interval contains a gene whose inactivation is important for the development of at least a fraction of WTs.

Mutational and functional analyses reveal that ST7 is a highly conserved tumor-suppressor gene on human chromosome 7q31

Loss of heterozygosity (LOH) of markers on human chromosome 7q31 is frequently encountered in a variety of human neoplasias, indicating the presence of a tumor-suppressor gene (TSG). By a combination of microcell-fusion and deletion-mapping studies, we previously established that this TSG resides within a critical region flanked by the genetic markers D7S522 and D7S677. Using a positional cloning strategy and aided by the availability of near-complete sequence of this genomic interval, we have identified a TSG within 7q31, named ST7 (for suppression of tumorigenicity 7; this same gene was recently reported in another context and called RAY1). ST7 is ubiquitously expressed in human tissues. Analysis of a series of cell lines derived from breast tumors and primary colon carcinomas revealed the presence of mutations in ST7. Introduction of the ST7 cDNA into the prostate-cancer-derived cell line PC3 had no effect on the in vitro proliferation of the cells, but abrogated their in vivo tumorigenicity. Our data indicate that ST7 is a TSG within chromosome 7q31 and may have an important role in the development of some types of human cancer.

Comparison of human genetic and sequence-based physical maps

Recombination is the exchange of information between two homologous chromosomes during meiosis. The rate of recombination per nucleotide, which profoundly affects the evolution of chromosomal segments, is calculated by comparing genetic and physical maps. Human physical maps have been constructed using cytogenetics, overlapping DNA clones and radiation hybrids; but the ultimate and by far the most accurate physical map is the actual nucleotide sequence. The completion of the draft human genomic sequence provides us with the best opportunity yet to compare the genetic and physical maps. Here we describe our estimates of female, male and sex-average recombination rates for about 60% of the genome. Recombination rates varied greatly along each chromosome, from 0 to at least 9 centiMorgans per megabase (cM Mb(-1)). Among several sequence and marker parameters tested, only relative marker position along the metacentric chromosomes in males correlated strongly with recombination rate. We identified several chromosomal regions up to 6 Mb in length with particularly low (deserts) or high (jungles) recombination rates. Linkage disequilibrium was much more common and extended for greater distances in the deserts than in the jungles.

A physical map of the human genome

The human genome is by far the largest genome to be sequenced, and its size and complexity present many challenges for sequence assembly. The International Human Genome Sequencing Consortium constructed a map of the whole genome to enable the selection of clones for sequencing and for the accurate assembly of the genome sequence. Here we report the construction of the whole-genome bacterial artificial chromosome (BAC) map and its integration with previous landmark maps and information from mapping efforts focused on specific chromosomal regions. We also describe the integration of sequence data with the map.

Mutational analysis of patients with p47-phox-deficient chronic granulomatous disease: The significance of recombination events between the p47-phox gene (NCF1) and its highly homologous pseudogenes

OBJECTIVE

The aim of this study was to determine the molecular basis of p47-phox-deficient chronic granulomatous disease (CGD), the most common autosomal recessive form of the disease. CGD is an inherited condition characterized by defective oxygen radical production due to defects in the phagocyte nicotinamide adenine dinucleotide phosphate (NADPH) oxidase. Mutational analysis of p47-phox-deficient CGD patients previously demonstrated that the majority of patients have a GT dinucleotide (Delta GT) deletion at the start of exon 2, a signature sequence also observed in the highly homologous pseudogenes of NCF1.

MATERIALS AND METHODS

We performed genetic analysis of NCF1 and its pseudogenes using genomic DNA in 29 p47-phox-deficient CGD patients from 22 separate families. First-strand cDNA analysis was performed in 17 of the 29 patients.

RESULTS

We confirmed the significance of the Delta GT mutation; in 27 of 29 patients, only the Delta GT sequence was detectable. All but one of the 27 had at least one additional signature sequence, specific to the pseudogene, in either intron 1 and/or intron 2. We extended our analysis to look at signature sequence differences in exons 6 and 9 and detected both the wild-type and pseudogene sequences in all patients tested.

CONCLUSIONS

Although detection of only Delta GT sequence accounts for over 85% of affected patients, the molecular basis is most likely due to partial cross-over events between the wild-type and pseudogene(s) of p47-phox at different recombination sites. Our results suggest that complete gene conversion or deletion of the p47-phox gene (NCF1) occurs rarely, if it all.

Computational and experimental analyses reveal previously undetected coding exons of the KRIT1 (CCM1) gene

A notable difficulty in annotating genomic sequence is identifying the correct start codon in a gene. An important such case has been found with KRIT1, the cerebral cavernous malformation type 1 (CCM1) gene. Analysis of human and mouse genomic sequence encompassing the region containing KRIT1/Krit1 using exon/gene-prediction and comparative alignment programs revealed putative exons upstream of the previously described first exon. These additional candidate exons show significant matches to mouse and human ESTs that are contiguous with and extend upstream from the previously designated 5' end of the KRIT1 cDNA sequence. RT-PCR and 5'RACE experiments confirm the presence of four additional upstream coding exons that encode an additional 207 amino acids. Importantly, a novel frameshift mutation in one of these newly identified KRIT1 exons has been found in a CCM1 family. These data establish the authentic KRIT1 amino acid sequence and suggest that the additional KRIT1 exons may harbor mutations in other CCM1 families. In addition, these results provide another example of the utility of rigorous computational and comparative sequence analysis for refining gene structure.

Childhood-onset schizophrenia/autistic disorder and t(1;7) reciprocal translocation: identification of a BAC contig spanning the translocation breakpoint at 7q21

Childhood-onset schizophrenia (COS) is defined by the development of first psychotic symptoms by age 12. While recruiting patients with COS refractory to conventional treatments for a trial of atypical antipsychotic drugs, we discovered a unique case who has a familial t(1;7)(p22;q21) reciprocal translocation and onset of psychosis at age 9. The patient also has symptoms of autistic disorder, which are usually transient before the first psychotic episode among 40-50% of the childhood schizophrenics but has persisted in him even after the remission of psychosis. Cosegregating with the translocation, among the carriers in the family available for the study, are other significant psychopathologies, including alcohol/drug abuse, severe impulsivity, and paranoid personality and language delay. This case may provide a model for understanding the genetic basis of schizophrenia or autism. Here we report the progress toward characterization of genomic organization across the translocation breakpoint at 7q21. The polymorphic markers, D7S630/D7S492 and D7S2410/D7S646, immediately flanking the breakpoint, may be useful for further confirming the genetic linkage for schizophrenia or autism in this region. Am. J. Med. Genet. (Neuropsychiatr. Genet.) 96:749-753, 2000. Published 2000 Wiley-Liss, Inc.

A novel genetic locus for low renin hypertension: familial hyperaldosteronism type II maps to chromosome 7 (7p22)

Familial hyperaldosteronism type II (FH-II) is caused by adrenocortical hyperplasia or aldosteronoma or both and is frequently transmitted in an autosomal dominant fashion. Unlike FH type I (FH-I), which results from fusion of the CYP11B1 and CYP11B2 genes, hyperaldosteronism in FH-II is not glucocorticoid remediable. A large family with FH-II was used for a genome wide search and its members were evaluated by measuring the aldosterone:renin ratio. In those with an increased ratio, FH-II was confirmed by fludrocortisone suppression testing. After excluding most of the genome, genetic linkage was identified with a maximum two point lod score of 3.26 at theta=0, between FH-II in this family and the polymorphic markers D7S511, D7S517, and GATA24F03 on chromosome 7, a region that corresponds to cytogenetic band 7p22. This is the first identified locus for FH-II; its molecular elucidation may provide further insight into the aetiology of primary aldosteronism.

Integration of cytogenetic data with genome maps and available probes: present status and future promise

The National Cancer Institute has established an initiative, called the Cancer Chromosome Aberration Project (Ccap), in order to link and integrate the physical and genetic maps of the human genome with cytogenetic data and the location of chromosomal rearrangements in human diseases. This goal will be achieved by high-resolution fluorescence in situ hybridization (FISH) mapping of colony-purified bacterial artificial chromosome (BAC) clones spaced at 1-to 2-Mb intervals across the entire genome. All BAC clones will be anchored on the physical map by the presence of a mapped sequence tagged site (STS). The generation of a publicly accessible clone repository will allow convenient distribution of these BACs. Ccap data can be correlated with other cancer-associated and genomic databases, such as the catalog of chromosomal aberrations in cancer and the emerging full genomic sequence. We anticipate that the use of Ccap clones will expedite and refine the mapping of chromosomal breakpoints. The eventual set of approximately 3,000 Ccap BACs should facilitate the production of BAC-containing DNA chips for assessing copy number of genomic segments by matrix comparative genomic hybridization. In addition, the repository will provide genome-wide tools for defining chromosomal aberrations in cytological specimens by interphase cytogenetics. The Ccap Web site illustrates goals and progress of this initiative (http://www.ncbi.nlm.nih.gov/CCAP/).

The SPCH1 region on human 7q31: genomic characterization of the critical interval and localization of translocations associated with speech and language disorder

The KE family is a large three-generation pedigree in which half the members are affected with a severe speech and language disorder that is transmitted as an autosomal dominant monogenic trait. In previously published work, we localized the gene responsible (SPCH1) to a 5.6-cM region of 7q31 between D7S2459 and D7S643. In the present study, we have employed bioinformatic analyses to assemble a detailed BAC-/PAC-based sequence map of this interval, containing 152 sequence tagged sites (STSs), 20 known genes, and >7.75 Mb of completed genomic sequence. We screened the affected chromosome 7 from the KE family with 120 of these STSs (average spacing <100 kb), but we did not detect any evidence of a microdeletion. Novel polymorphic markers were generated from the sequence and were used to further localize critical recombination breakpoints in the KE family. This allowed refinement of the SPCH1 interval to a region between new markers 013A and 330B, containing approximately 6.1 Mb of completed sequence. In addition, we have studied two unrelated patients with a similar speech and language disorder, who have de novo translocations involving 7q31. Fluorescence in situ hybridization analyses with BACs/PACs from the sequence map localized the t(5;7)(q22;q31.2) breakpoint in the first patient (CS) to a single clone within the newly refined SPCH1 interval. This clone contains the CAGH44 gene, which encodes a brain-expressed protein containing a large polyglutamine stretch. However, we found that the t(2;7)(p23;q31.3) breakpoint in the second patient (BRD) resides within a BAC clone mapping >3.7 Mb distal to this, outside the current SPCH1 critical interval. Finally, we investigated the CAGH44 gene in affected individuals of the KE family, but we found no mutations in the currently known coding sequence. These studies represent further steps toward the isolation of the first gene to be implicated in the development of speech and language.

The familial mediterranean fever protein interacts and colocalizes with a putative Golgi transporter

The biological function of pyrin, the protein mutated in Familial Mediterranean Fever (FMF), has not been elucidated. Based on sequence homology, a transcription factor activity was proposed for this neutrophil-specific protein. In a yeast two-hybrid assay, neither transcription activation activity nor any self interaction was detected for pyrin. Screening of an expression cDNA library of peripheral blood leukocytes using as bait the carboxyl portion of pyrin (amino acids 557-781), which contains most of the FMF mutations, led to the identification of P/M-IP1 (pyrin/marenostrin interacting protein 1). A splice variant of P/M-IP1, GTC-90, had previously been described as a component of the 13S hetero-oligomeric protein complex that stimulates in vitro Golgi transport. We have now shown that P/M-IP1 colocalizes with pyrin in the perinuclear cytoplasm of Cos-7 cells and that the interaction between these two proteins is impaired by FMF causing mutations in pyrin. These data suggest that, at some stage of its functional pathway, pyrin resides in the cytoplasm and might be involved in, or impacted by, cellular protein sorting by the Golgi apparatus. The data also imply that P/M-IP1 may be involved in the abnormal inflammatory response that occurs in patients with FMF.

Comparative genome mapping in the sequence-based era: early experience with human chromosome 7

The success of the ongoing Human Genome Project has resulted in accelerated plans for completing the human genome sequence and the earlier-than-anticipated initiation of efforts to sequence the mouse genome. As a complement to these efforts, we are utilizing the available human sequence to refine human-mouse comparative maps and to assemble sequence-ready mouse physical maps. Here we describe how the first glimpses of genomic sequence from human chromosome 7 are directly facilitating these activities. Specifically, we are actively enhancing the available human-mouse comparative map by analyzing human chromosome 7 sequence for the presence of orthologs of mapped mouse genes. Such orthologs can then be precisely positioned relative to mapped human STSs and other genes. The chromosome 7 sequence generated to date has allowed us to more than double the number of genes that can be placed on the comparative map. The latter effort reveals that human chromosome 7 is represented by at least 20 orthologous segments of DNA in the mouse genome. A second component of our program involves systematically analyzing the evolving human chromosome 7 sequence for the presence of matching mouse genes and expressed-sequence tags (ESTs). Mouse-specific hybridization probes are designed from such sequences and used to screen a mouse bacterial artificial chromosome (BAC) library, with the resulting data used to assemble BAC contigs based on probe-content data. Nascent contigs are then expanded using probes derived from newly generated BAC-end sequences. This approach produces BAC-based sequence-ready maps that are known to contain a gene(s) and are homologous to segments of the human genome for which sequence is already available. Our ongoing efforts have thus far resulted in the isolation and mapping of >3,800 mouse BACs, which have been assembled into >100 contigs. These contigs include >250 genes and represent approximately 40% of the mouse genome that is homologous to human chromosome 7. Together, these approaches illustrate how the availability of genomic sequence directly facilitates studies in comparative genomics and genome evolution.

Support for linkage of autism and specific language impairment to 7q3 from two chromosome rearrangements involving band 7q31

Childhood autism is characterised by impairments in communication and reciprocal social interaction together with restricted/stereotyped interests, which are evident before 3 years of age. Specific developmental disorders of speech and language (SDDSL) are characterised by impairment in the development of expressive and/or receptive language skills which is not associated with intellectual, sensory, physical, or neurological impairment. Family and twin studies indicate a substantial genetic component in the aetiology of both disorders. They also reveal increased rates of SDDSL in relatives of autistic individuals, suggesting that this phenotype can represent one manifestation of the genetic liability for autism. Modelling of the recurrence risk for autism and milder phenotypes, such as SDDSL, suggest that three or four epistatic loci may be aetiologically involved. A recently published linkage study of an exceptional family with an apparently dominantly inherited SDDSL implicated chromosome band 7q31 as the site of the putative susceptibility locus (SPCH1). This region of chromosome 7 also shows strong linkage in multiplex families with autism. We present two individuals (one has autism, the other SDDSL) with different, apparently balanced chromosome rearrangements involving a breakpoint at 7q31.3. Fluorescence in situ hybridisation was used to localise the breakpoints to an approximately 1 cM interval between CFTR and D7S643. Our findings may be of interest and relevance to the genetic aetiology of autism, and helpful in the search for susceptibility loci for SDDSL and autism. Am. J. Med. Genet. (Neuropsychiatr. Genet. ) 96:228-234, 2000.

The reliability of the Aspergillus nidulans physical map

Here we report an evaluation of the Aspergillus nidulans physical map (a cosmid contig map) emphasizing quantification and description of obvious mapping errors. Classification and appraisal of mapping errors should be helpful to researchers working on particular regions of the map. We estimate between 47 (4.1%) and 63 (5.4%) probe/clone-linking errors. The majority of identified false links (38) permit reciprocal exchanges among linking clones located on disconnected mapping regions. The order of adjacent clones or probes on the affected contigs remains unchanged. In addition we describe an Internet-accessible resource in which genetic and physical maps were integrated through a graphic interface. A simple search engine allows retrieval of cosmids from redundant clone lists and provides links to the minimal clone order. Integration of genetic and physical maps provides an additional level of accountability in which mapping discrepancies are visually located.

Comparative mapping of the human homologue of the rat diabetes susceptibility gene lyp to a 1.3-Mb segment on HSA7

The rat diabetes susceptibility gene, Lyp or Lymphopenia, has been localized to RNO4. Proximal to Lyp are the genes caspase-2 (Casp2) and pancreatic trypsin 1 (Prss1), while neuropeptide Y (Npy) is the closest distally positioned gene. In human, the three genes are syntenic on HSA7, but they are not on a conserved segment: CASP2 and PRSS1 are localized to 7q35, while NPY is localized to 7p15.1. This raises the question whether the human homologue of Lyp is linked to CASP2/PRSS1 or to NPY. We present a comparative map of the Lyp region in rat and human, assigning the gene to a 1.3-Mb segment between RNY3 and ABP1 at 7q35.

Detecting ancient admixture in humans using sequence polymorphism data

A debate of long-standing interest in human evolution centers around whether archaic human populations (such as the Neanderthals) have contributed to the modern gene pool. A model of ancient population structure with recent mixing is introduced, and it is determined how much information (i.e., sequence data from how many unlinked nuclear loci) would be necessary to distinguish between different demographic scenarios. It is found that approximately 50-100 loci are necessary if plausible parameter estimates are used. There are not enough data available at the present to support either the "single origin" or the "multiregional" model of modern human evolution. However, this information should be available in a few years.

Comparative genomic sequence analysis of the human and mouse cystic fibrosis transmembrane conductance regulator genes

The identification of the cystic fibrosis transmembrane conductance regulator gene (CFTR) in 1989 represents a landmark accomplishment in human genetics. Since that time, there have been numerous advances in elucidating the function of the encoded protein and the physiological basis of cystic fibrosis. However, numerous areas of cystic fibrosis biology require additional investigation, some of which would be facilitated by information about the long-range sequence context of the CFTR gene. For example, the latter might provide clues about the sequence elements responsible for the temporal and spatial regulation of CFTR expression. We thus sought to establish the sequence of the chromosomal segments encompassing the human CFTR and mouse Cftr genes, with the hope of identifying conserved regions of biologic interest by sequence comparison. Bacterial clone-based physical maps of the relevant human and mouse genomic regions were constructed, and minimally overlapping sets of clones were selected and sequenced, eventually yielding approximately 1.6 Mb and approximately 358 kb of contiguous human and mouse sequence, respectively. These efforts have produced the complete sequence of the approximately 189-kb and approximately 152-kb segments containing the human CFTR and mouse Cftr genes, respectively, as well as significant amounts of flanking DNA. Analyses of the resulting data provide insights about the organization of the CFTR/Cftr genes and potential sequence elements regulating their expression. Furthermore, the generated sequence reveals the precise architecture of genes residing near CFTR/Cftr, including one known gene (WNT2/Wnt2) and two previously unknown genes that immediately flank CFTR/Cftr.

Genomic structure of the human p47-phox (NCF1) gene

The cytosolic factor p47-phox, encoded by the NCF1 gene, is an essential component of the phagocyte NADPH-oxidase system. Upon activation of this multicomponent system, p47-phox translocates to the membrane and participates in the electron transfer from NADPH to molecular oxygen. A deficiency or absence of p47-phox is the most common autosomal form of chronic granulomatous disease (CGD). We have cloned and characterized the NCF1 gene from four bacteriophage clones, a P1 clone and genomic DNA from normal individuals. The gene is 15,236 base pairs long and includes 11 exons. It is 98.6% homologous in sequence to at least one pseudogene that maps to the same region of chromosome 7q11.23. Slightly more than half (50.37%) of the wild-type NCF1 gene consists of repetitive elements. In particular, the density of Alu sequences is high (1.4 Alu/kb); there are 21 Alu repeats interspersed through 10 introns. These findings are consistent with the observation that recombination events between the wild-type gene and its highly homologous pseudogenes account for the majority of potentially lethal mutations in p47-phox-deficient chronic granulomatous disease. Analysis of 1.96 kb of sequence 5' of the start of translation revealed a high homology (99.6%) between wild-type and pseudogene clones. Characterization of NCF1 establishes a foundation for detailed molecular analysis of p47-phox-deficient CGD patients as well as for the study of the regulation of the NCF1 gene and pseudogenes, both of which are present as full-length transcripts in normal individuals.

A comparison of estimators of the population recombination rate

Three new estimators of the population recombination rate C = 4Nr are introduced. These estimators summarize the data using the number of distinct haplotypes and the estimated minimum number of recombination events, then calculate the value of C that maximizes the likelihood of obtaining the summarized data. They are compared with a number of previously proposed estimators of the recombination rate. One of the newly proposed estimators is generally better than the others for the parameter values considered here, while the three programs that calculate maximum-likelihood estimates give conflicting results.

A physical map, including a BAC/PAC clone contig, of the Williams-Beuren syndrome--deletion region at 7q11.23

Williams-Beuren syndrome (WBS) is a developmental disorder caused by haploinsufficiency for genes in a 2-cM region of chromosome band 7q11.23. With the exception of vascular stenoses due to deletion of the elastin gene, the various features of WBS have not yet been attributed to specific genes. Although >/=16 genes have been identified within the WBS deletion, completion of a physical map of the region has been difficult because of the large duplicated regions flanking the deletion. We present a physical map of the WBS deletion and flanking regions, based on assembly of a bacterial artificial chromosome/P1-derived artificial chromosome contig, analysis of high-throughput genome-sequence data, and long-range restriction mapping of genomic and cloned DNA by pulsed-field gel electrophoresis. Our map encompasses 3 Mb, including 1.6 Mb within the deletion. Two large duplicons, flanking the deletion, of >/=320 kb contain unique sequence elements from the internal border regions of the deletion, such as sequences from GTF2I (telomeric) and FKBP6 (centromeric). A third copy of this duplicon exists in inverted orientation distal to the telomeric flanking one. These duplicons show stronger sequence conservation with regard to each other than to the presumptive ancestral loci within the common deletion region. Sequence elements originating from beyond 7q11.23 are also present in these duplicons. Although the duplicons are not present in mice, the order of the single-copy genes in the conserved syntenic region of mouse chromosome 5 is inverted relative to the human map. A model is presented for a mechanism of WBS-deletion formation, based on the orientation of duplicons' components relative to each other and to the ancestral elements within the deletion region.

An integrated physical map for the short arm of human chromosome 5

The short arm of human chromosome 5 contains approximately 48 Mb of DNA and comprises 1.5% of the genome. We have constructed a mega-YAC/ STS map of this region that includes 436 YACs anchored by 216 STSs. By combining and integrating our map with the 5p maps of other groups using the same recombinant DNA library, a comprehensive map was constructed that includes 552 YACs and 504 markers. The YAC map covers >94% of 5p in four YAC contigs, bridges the centromere, and includes an additional 5 Mb of 5q DNA. The average marker density is 95 kb. This integrated 5p map will serve as a resource for the continuing localization of genes on the short arm of human chromosome 5 and as a framework for both generating and aligning the DNA sequence of this region.

Mutations in the gene encoding KRIT1, a Krev-1/rap1a binding protein, cause cerebral cavernous malformations (CCM1)

Cerebral cavernous malformations (CCM) are congenital vascular anomalies of the brain that can cause significant neurological disabilities, including intractable seizures and hemorrhagic stroke. One locus for autosomal dominant CCM ( CCM1 ) maps to chromosome 7q21-q22. Recombination events in linked family members define a critical region of approximately 2 Mb and a shared disease haplotype associated with a presumed founder effect in families of Mexican-American descent points to a potentially smaller region of interest. Using a genomic sequence-based positional cloning strategy, we have identified KRIT1, encoding a protein that interacts with the Krev-1/rap1a tumor suppressor, as the CCM1 gene. Seven different KRIT1 mutations have been identified in 23 distinct CCM1 families. The identical mutation is present in 16 of 21 Mexican-American families analyzed, substantiating a founder effect in this population. Other Mexican-American and non-Hispanic Caucasian CCM1 kindreds harbor other KRIT1 mutations. Identification of a common Mexican-American mutation has potential clinical significance for presymptomatic diagnosis of CCM in this population. In addition, these data point to a key role for the Krev-1/rap1a signaling pathway in angiogenesis and cerebrovascular disease.

Loss of heterozygosity of the long arm of chromosome 7 in follicular and anaplastic thyroid cancer, but not in papillary thyroid cancer

Papillary thyroid cancer (PTC), but neither the follicular nor the anaplastic histotype [follicular thyroid cancer (FTC), anaplastic thyroid cancer (ATC)], overexpresses simultaneously the protooncogene HGF (hepatocyte growth factor) and its receptor HGF-R (or c-met). Because 1) HGF and c-met map to chromosome 7q21 and 7q31, respectively, 2) FTC loses genetic material at multiple loci with a frequency much higher than PTC, and 3) loss of heterozygosity (LOH) on 7q has been previously found in various tumors, we tested the hypothesis that both FTC and ATC, but not PTC, could harbor LOH in segments of 7q encompassing the loci for HGF and c-met. We screened 6 normal thyroids, 10 colloid nodules, 10 follicular hyperplasias, 10 oncocytic adenomas, 10 follicular adenomas (FA), 10 FTC, 6 ATC, 12 PTC using two microsatellite markers for HGF, and two for c-met. LOH for all 4 markers was found in 100% of FTC, 100% of ATC, and (for only 1 or 2 markers) in 10-29% of FA. This is the first demonstration of an LOH that separates both FTC and ATC from PTC, in the best possible manner: 100% vs. 0%. Clearly, each of the two segments we have probed contains at least one tumor suppressor gene, whose inactivation is crucial for the establishment of the FTC (and ATC) phenotype. This loss of genetic material explains why FTC and ATC, but not PTC, fail to express both HGF and c-met. Our findings may also have immediate diagnostic application, in the context of assisting pathologists in the often difficult task of distinguishing FA from FTC.

An algorithm for the assembly of robust physical maps based on a combination of multi-level hybridization data and fingerprinting data

We have developed an algorithm which combines data obtained from restriction digestion experiments and hybridization experiments to construct robust physical maps of whole chromosomes. The algorithm has been incorporated into a program which accepts hybridization data consisting of an unordered hybridization matrix and fingerprinting data containing band coordinates for each clone. The combined data is used to produce a non-redundant, ordered matrix which can be further reduced to represent a minimum tile coverage of the chromosome. In addition, the method also takes into account multi-level hybridization events which allows for an improved treatment of the hybridization data. The program is evaluated against several other contig building programs using simulated and real data sets. Finally, it is applied to construct a physical map of the 4.1 Mb genome of Ochrobactrum anthropi based on 1387 clones and 70 probes, as well as 624 fingerprints.

The CMT2D Locus: Refined Genetic Position and Construction of a Bacterial Clone-Based Physical Map

Charcot-Marie-Tooth (CMT) disease is a progressive neuropathy of the peripheral nervous system, typically characterized by muscle weakness of the distal limbs. CMT is noted for its genetic heterogeneity, with four distinct loci already identified for the axonal form of the disease (CMT2). In 1996, linkage analysis of a single large family revealed the presence of a CMT2 locus on chromosome 7p14 (designatedCMT2D). Additional families have been linked subsequently to the same genomic region, including one with distal spinal muscular atrophy (dSMA) and one with mixed features of dSMA and CMT2; symptoms in both of these latter families closely resemble those seen in the original CMT2D family. There is thus a distinct possibility that CMT2 and dSMA encountered in these families reflect allelic heterogeneity at a single chromosome 7 locus. In the study reported here, we have performed more detailed linkage analysis of the original CMT2D family based on new knowledge of the physical locations of various genetic markers. The region containing the CMT2D gene, as defined by the original family, overlaps with those defined by at least two other families with CMT2 and/or dSMA symptoms. Both yeast artificial chromosome (YAC) and bacterial clone-based [bacterial artificial chromosome (BAC) and P1-derived artificial chromosome (PAC)] contig maps spanning ∼3.4 Mb have been assembled across the combinedCMT2D critical region, with the latter providing suitable clones for systematic sequencing of the interval. Preliminary analyses have already revealed at least 28 candidate genes and expressed-sequence tags (ESTs). The mapping information reported here in conjunction with the evolving sequence data should expedite the identification of the CMT2D/dSMA gene or genes.

Construction of a high-resolution physical map of the approximate 1-Mb region of human chromosome 7q31.1-q31.2 harboring a putative tumor suppressor gene

Reports of frequent loss of heterozygosity (LOH) of markers on human chromosome 7q in malignant myeloid disorders as well as breast, prostate, ovarian, colon, head and neck, gastric, pancreatic, and renal cell carcinomas suggest the presence of a tumor suppressor gene (TSG). Functional assays have demonstrated that the introduction of an intact copy of human chromosome 7 (hchr7) can restore senescence to immortalized human fibroblast cell lines having LOH of markers within 7q31-q32 and can inhibit the tumorigenic phenotype of a murine squamous cell carcinoma cell line. To facilitate the cloning of the putative TSG, we have constructed a high-resolution physical map of this region of hchr7, specifically that encompassing the markers D7S522 and D7S677 within 7q31.1-q31.2. By using a lower resolution yeast artificial chromosome-based map as a starting framework, we established complete clone coverage of the implicated critical region in bacterial-artificial chromosomes (BACs) and P1-derived artificial chromosomes (PACs). The resulting BAC/PAC-based contig map has provided suitable clones for the systematic sequencing of the entire interval. In addition, we have already identified 29 clusters of overlapping expressed-sequence tags (ESTs) and 4 known genes contained within these clones. Together, the physical map reported here coupled with the evolving sequence and gene maps should hasten the identification of the putative TSG residing within this region of hchr7.

Leishmania major Friedlin chromosome 1 has an unusual distribution of protein-coding genes

Leishmania are evolutionarily ancient protozoans (Kinetoplastidae) and important human pathogens that cause a spectrum of diseases ranging from the asymptomatic to the lethal. The Leishmania genome is relatively small [ approximately 34 megabases (Mb)], lacks substantial repetitive DNA, and is distributed among 36 chromosomes pairs ranging in size from 0.3 Mb to 2.5 Mb, making it a useful candidate for complete genome sequence determination. We report here the nucleotide sequence of the smallest chromosome, chr1. The sequence of chr1 has a 257-kilobase region that is densely packed with 79 protein-coding genes. This region is flanked by telomeric and subtelomeric repetitive elements that vary in number and content among the chr1 homologs, resulting in an approximately 27.5-kilobase size difference. Strikingly, the first 29 genes are all encoded on one DNA strand, whereas the remaining 50 genes are encoded on the opposite strand. Based on the gene density of chr1, we predict a total of approximately 9,800 genes in Leishmania, of which 40% may encode unknown proteins.

“KARIBIN,” an Information Resource for Obtaining Genomic Information in a Cytogenetic Band

KARIBIN (http://sgiweb.ncbi.nlm.nih.gov:80/Zjing/yac.html) is a karyotypicregion-based integrated information resource that provides a comprehensive view of the integrated mapping and sequencing data for the human genome. A cytogenetic band is linked to a genetic or physical location using fluorescence in situ hybridization (FISH) mapping data. The genetic, physical mapping data and the sequencing data are integrated using STS markers positioned on multiple maps. For each cytogenetic band, the user can obtain the most up-to-date information that includes genetic and physical maps, human transcript gene map, YAC and PAC/BAC clone coverage, disease gene phenotype, and high throughput genomic sequences from the major human genome sequencing centers. This information provides a framework for future experiments and may accelerate the process of disease gene hunting. It is envisioned that other cytogenetic-based information such as chromosome aberrations can be linked to this framework.

Toward a complete human genome sequence

We have begun a joint program as part of a coordinated international effort to determine a complete human genome sequence. Our strategy is to map large-insert bacterial clones and to sequence each clone by a random shotgun approach followed by directed finishing. As of September 1998, we have identified the map positions of bacterial clones covering approximately 860 Mb for sequencing and completed >98 Mb ( approximately 3.3%) of the human genome sequence. Our progress and sequencing data can be accessed via the World Wide Web (http://webace.sanger.ac.uk/HGP/ or http://genome.wustl.edu/gsc/).