Global patterns of human DNA sequence variation in a 10-kb region on chromosome 1

Human DNA variation is currently a subject of intense research because of its importance for studying human origins, evolution, and demographic history and for association studies of complex diseases. A approximately 10-kb region on chromosome 1, which contains only four small exons (each <155 bp), was sequenced for 61 humans (20 Africans, 20 Asians, and 21 Europeans) and for 1 chimpanzee, 1 gorilla, and 1 orangutan. We found 52 polymorphic sites among the 122 human sequences and 382 variant sites among the human, chimpanzee, gorilla, and orangutan sequences. For the introns sequenced (8,991 bp), the nucleotide diversity (pi) was 0.058% among all sequences, 0.076% among the African sequences, 0.047% among the Asian sequences, and 0.045% among the European sequences. A compilation of data revealed that autosomal regions have, on average, the highest pi value (0.091%), X-linked regions have a somewhat lower pi value (0.079%), and Y-linked regions have a very low pi value (0.008%). The lower polymorphism in the present region may be due to a lower mutation rate and/or selection in the gene containing these introns or in genes linked to this region. The present region and two other 10-kb noncoding regions all show a strong excess of low-frequency variants, indicating a relatively recent population expansion. This region has a low mutation rate, which was estimated to be 0.74 x 10 per nucleotide per year. An average estimate of approximately 12,600 for the long-term effective population size was obtained using various methods; the estimate was not far from the commonly used value of 10,000. Fu and Li's tests rejected the assumption of an equilibrium neutral Wright-Fisher population, largely owing to the high proportion of low-frequency variants. The age of the most recent common ancestor of the sequences in our sample was estimated to be more than 1 Myr. Allowing for some unrealistic assumptions in the model, this estimate would still suggest an age of more than 500,000 years, providing further evidence for a genetic history of humans much more ancient than the emergence of modern humans. The fact that many unique variants exist in Europe and Asia also suggests a fairly long genetic history outside of Africa and argues against a complete replacement of all indigenous populations in Europe and Asia by a small Africa stock. Moreover, the ancient genetic history of humans indicates no severe bottleneck during the evolution of humans in the last half million years; otherwise, much of the ancient genetic history would have been lost during a severe bottleneck. We suggest that both the "Out of Africa" and the multiregional models are too simple to explain the evolution of modern humans.     

Nonuniform linkage disequilibrium within a 1,500-kb region of the human immunoglobulin heavy-chain complex.

We have characterized 10 VH polymorphic loci of the VH2, VH3, VH4, and VH5 families. Eight of 10 VH polymorphisms were found to be insertion/deletion polymorphisms, probably the result of nonhomologous recombination over the course of evolution of the current human VH repertoire. The 10 VH polymorphic loci were analyzed in 10 three-generation and 10 two-generation Canadian caucasoid families. Linkage disequilibrium (allelic association) was measured between pairs of VH polymorphic loci, and 12 significant associations were found. The degree of linkage disequilibrium measured between IGH polymorphic loci was then compared with the physical distance separating the loci. The physical distance between IGH polymorphic loci does not entirely determine the degree of linkage disequilibrium between polymorphic loci. Two regions, one in the VH region (between VH3f-2 and VH5-2 and one in the CH region (between C delta and C gamma 3), were found to have linkage disequilibrium values approximately 1/3,000 of that observed in other portions of the IGH region. The previous identification of recombinants in the C delta-to C gamma 3 region indicates that these areas of low linkage disequilibrium are consistent with the presence of recombination hot spots. The observed high amount of recombination in the subtelomeric portion of chromosome 14 therefore appears to be the result of specific hot spots for recombination, rather than a general increase in recombination in this region.     

Closing in on the Rieger syndrome gene on 4q25: mapping translocation breakpoints within a 50-kb region.

Rieger syndrome (RGS) is an autosomal dominant disorder of morphogenesis affecting mainly the formation of the anterior eye chamber and of the teeth. RGS has been localized to human chromosome 4q25 by linkage to epidermal growth factor (EGF). We have constructed a detailed physical map and a YAC contig of the genomic region encompassing the EGF locus. Using FISH, several YACs could be shown to cross the breakpoint in two independent RGS patients with balanced 4q translocations. Alu- and LINE-fragmentation of a 2.4-Mb YAC generated a panel of shorter YACs ranging in size from 2.4 Mb to 75 kb. Several fragmentation YACs were subcloned in cosmids, which were mapped to specific subregions of the original YAC by hybridization to the fragmentation panel to further refine the localization of the translocation breakpoints, allowing mapping of the breakpoints to within the most-telomeric 200 kb of the original 2.4-Mb YAC. FiberFISH of cosmids located in this 200-kb region mapped the two translocation breakpoints within a 50-kb region approximately 100-150 kb centromeric to D4S193, significantly narrowing down the candidate region for RGS. The mapping data and resources reported here should facilitate the identification of a gene implicated in Rieger syndrome.     

Conservation of pericentromeric duplications of a 200-kb part of the human 21q22.1 region in primates

We analyzed the conservation of large paralogous regions (more than 200 kb) on human chromosome regions 21q22.1 and 21q11.2 and on pericentromeric regions of chromosomes 2, 13, and 18 in three nonhuman primate species. Orthologous regions were found by FISH analysis of metaphase chromosomes from Gorilla gorilla, Pan troglodytes, and Pongo pygmaeus. Only one orthologous region was detected in chromosomes of P. pygmaeus, showing that the original locus was at 21q22.1 and that the duplication arose after the separation of Asian orangutans from the other hominoids. Surprisingly, the paralogous regions were more highly conserved in gorilla than in chimpanzee. PCR amplification of STSs derived from sequences of the chromosome 21 loci and low-stringency FISH analysis showed that this duplication occurred recently in the evolution of the genome. Different rates of sequence evolution through substitutions or deletions, after the duplication, may have resulted in diversity between closely related primates.     

An 800-kb region of deletion at 13q14 in human prostate and other carcinomas

Deletions of regions at 13q14 have been detected by various genetic approaches in human cancers including prostate cancer. Several studies have defined one region of loss of heterozygosity (LOH) at 13q14 that seems to reside in a DNA segment of 7.1 cM between genetic markers D13S263 and D13S153. To define the smallest region of overlap (SRO) for deletion at 13q14, we first applied tissue microdissection and multiplex PCR to detect homozygous deletion and/or hemizygous deletion at 13q14 in 134 prostate cancer specimens from 114 patients. We detected deletions at markers D13S1227, D13S1272, and A005O48 in 13 (10%) of these tumor specimens. Of the 13 tumors with deletions, 12 were either poorly differentiated primary tumors or metastases of prostate cancer. To fine-map the deletion region, we then constructed a high-resolution YAC/BAC/STS/EST physical map based on experimental and database analyses. Several markers encompassing the deletion region were analyzed for homozygous deletion and/or hemizygous deletion in 61 cell lines/xenografts derived from human cancers of the prostate, breast, ovary, endometrium, cervix, and bladder, and a region of deletion was defined by duplex PCR assay between markers A005X38 and WI-7773. We also analyzed LOH at 13q14 in the 61 cell lines/xenografts using the homozygosity mapping of deletion approach and 26 microsatellite markers. We found 24 (39%) of the cell lines/xenografts to show LOH at 13q14 and defined a region of LOH by markers M1 and M5. Combination of homozygous or hemizygous deletion and LOH results defined the SRO for deletion to be an 800-kb DNA interval between A005X38 and M5. There are six known genes located in or close to the SRO for deletion. This region of deletion is at least 2 Mb centromeric to the RB1 tumor-suppressor gene and the leukemia-associated genes 1 and 2, each of which is located at 13q14. These data suggest that the 800-kb DNA segment with deletion contains a gene whose deletion may be important for the development of prostate and other cancers. This study also provides a framework for the fine-mapping, cloning, and identification of a novel tumor-suppressor gene at 13q14.     

Physical mapping by PFGE localizes the COL3A1 and COL5A2 genes to a 35-kb region on human chromosome 2

The genes encoding the alpha 1 chain of Type III collagen (COL3A1) and the alpha 2 chain of Type V (COL5A2) collagen have been mapped to the long arm of human chromosome 2. Linkage analysis in CEPH families indicated that these two genes are close to each other, with no recombination in 37 informative meioses. In the present study, DNA probes from the 3' ends of each gene have been physically mapped by pulsed-field gel electrophoresis. The probes recognized 11 macrorestriction fragments in common, ranging from greater than 1000 kb MluI and NotI fragments to a 35-kb SfiI fragment. Therefore, the COL3A1 and COL5A2 genes appear to exist as a gene cluster on chromosome 2. This is the third example of a collagen gene cluster. Other examples include the COL4A1-COL4A2 genes on chromosome 13q and the COL6A1-COL6A2 genes on chromosome 21q. The physical proximity of these genes may indicate common evolution and/or regulation.     

Clustering of two fragile sites and seven homeobox genes in human chromosome region 2q31-->q32.1

In this study we have used FISH to examine the relationship between a group of homeobox genes, namely DLX1/DLX2, EVX2 and four HOXD genes (10, 11, 12, 13), that map to region q31 on chromosome 2, and the FRA2G and FRA2H fragile sites located at 2q31 and 2q32.1 respectively. Our results indicate that these homeobox genes lie between the two fragile regions.     

Variable arrangement of 5S ribosomal genes within the ribosomal DNA repeats of arthropods.

The polymerase chain reaction and hybridization to genomic blots were used to investigate whether the previously observed inclusion of 5S ribosomal RNA genes in the 28S-18S ribosomal DNA intergenic regions of some crustacean species (copepods) could also be detected in other arthropods. Such an arrangement was found not only in other calanoid copepod species but also in a cirriped, an euphausid, and a spider. It is interesting that species from two different calanoid copepod genera do not have this type of arrangement. We conclude that the inclusion of 5S ribosomal RNA genes within the ribosomal DNA repeats has probably occurred repeatedly during the evolution of arthropod species and that the mechanism(s) responsible for these insertions could also be responsible for their loss.     

Co-localization of the ketohexokinase and glucokinase regulator genes to a 500-kb region of Chromosome 2p23

The glucokinase regulator (GCKR) is a 65-kDa protein that inhibits glucokinase (hexokinase IV) in liver and pancreatic islet. The role of glucokinase (GCK) as pancreatic β cell glucose sensor and the finding of GCK mutations in maturity onset diabetes of the young (MODY) suggest GCKR as a further candidate gene for type 2 diabetes. The inhibition of GCK by GCKR is relieved by the binding of fructose-1-phosphate (F-1-P) to GCKR. F-1-P is the end product of ketohexokinase (KHK, fructokinase), which, like GCK and GCKR, is present in both liver and pancreatic islet. KHK is the first enzyme of the specialized pathway that catabolizes dietary fructose. We have isolated genomic clones containing the human GCKR and KHK genes. By fluorescent in situ hybridization (FISH), KHK maps to Chromosome (Chr) 2p23.2-23.3, a new assignment corroborated by somatic cell hybrid analysis. The localization of GCKR, originally reported by others as 2p22.3, has been reassessed by high-resolution FISH, indicating that, like KHK, GCKR maps to 2p23.2-23.3. The proximity of GCKR and KHK was further demonstrated both by two-color interphase FISH, which suggests that the two genes lie within 500 kb of each other, and by analysis of overlapping YAC and P1 clones spanning the interval between GCKR and KHK. A new microsatellite polymorphism was used to place the GCKR-KHK locus between D2S305 and D2S165 on the genetic map. The co-localization of these two metabolically connected genes has implications for the interpretation of linkage or allele association studies in type 2 diabetes. It also raises the possibility of coordinate regulation of GCKR and KHK by common cis-acting regulatory elements. Received: 8 December 1995 / Accepted: 27 January 1996     

DNA sequences near a meiotic recombinational breakpoint within the human HLA-DQ region

The molecular organization of HLA-DQ regions derived from DR7, DQw2, and DR4, DQw3 parental haplotypes and DR7, DQw3, a presumed recombinant haplotype, have been studied to define the sequences between DQA1 and DQB1 which may have been involved in this recombinational event. The breakpoint was localized in the intergenic region near the 3 end of the DQB1 gene by restriction mapping. DNA sequences in the immediate vicinity of the breakpoint in DR7, DQw2 (parental), and DR7, DQw3 (recombinant) haplotypes revealed the presence of (CA)₂₂ repeats, minisatellite-related sequences and GC-rich sequences. The intergenic regions varied considerably depending on the haplotype and contained several additional types of repetitive sequences including Alu and LINE repeats. Some of these sequences are related to sequences previously suggested to be involved in meiotic or somatic recombination. In particular, (CA)_n repeats, which can adopt the Z-DNA conformation, have previously been shown to promote recombination in several systems. Address correspondence and offprint requests to: J. Strominger.     

Both hypomethylation and hypermethylation in a 0.2-kb region of a DNA repeat in cancer

NBL2 is a tandem 1.4-kb DNA repeat, whose hypomethylation in hepatocellular carcinomas was shown previously to be an independent predictor of disease progression. Here, we examined methylation of all cytosine residues in a 0.2-kb subregion of NBL2 in ovarian carcinomas, Wilms' tumors, and diverse control tissues by hairpin-bisulfite PCR. This new genomic sequencing method detects 5-methylcytosine on covalently linked complementary strands of a DNA fragment. All DNA clones from normal somatic tissues displayed symmetrical methylation at seven CpG positions and no methylation or only hemimethylation at two others. Unexpectedly, 56% of cancer DNA clones had decreased methylation at some normally methylated CpG sites as well as increased methylation at one or both of the normally unmethylated sites. All 146 DNA clones from 10 cancers could be distinguished from all 91 somatic control clones by assessing methylation changes at three of these CpG sites. The special involvement of DNA methyltransferase 3B in NBL2 methylation was indicated by analysis of cells from immunodeficiency, centromeric region instability, and facial anomalies syndrome patients who have mutations in the gene encoding DNA methyltransferase 3B. Blot hybridization of 33 cancer DNAs digested with CpG methylation-sensitive enzymes confirmed that NBL2 arrays are unusually susceptible to cancer-linked hypermethylation and hypomethylation, consistent with our novel genomic sequencing findings. The combined Southern blot and genomic sequencing data indicate that some of the cancer-linked alterations in CpG methylation are occurring with considerable sequence specificity. NBL2 is an attractive candidate for an epigenetic cancer marker and for elucidating the nature of epigenetic changes in cancer.     

Physical mapping of four serpin genes: alpha 1-antitrypsin, alpha 1-antichymotrypsin, corticosteroid-binding globulin, and protein C inhibitor, within a 280-kb region on chromosome I4q32.1.

Alpha 1-antitrypsin (alpha 1AT; protease inhibitor [PI] locus), alpha 1-antichymotrypsin (alpha 1ACT; AACT locus), corticosteroid-binding globulin (CBG; CBG locus), and protein C inhibitor (PCI; PCI locus) are members of the serine protease inhibitor (serpin) superfamily. A noncoding PI-like (PIL) gene has been located 12 kb 3' of the PI gene. The PI, PIL, and AACT loci have been localized to 14q32.1, the CBG locus has been localized to 14q31-14q32.1, and PCI has been mapped to chromosome 14. Genetic linkage analysis suggests tight linkage between PI and AACT. We have used pulsed-field gel electrophoresis to generate a physical map linking these five serpin genes. The order of the genetic loci is AACT/PCI-PI-PIL-CBG, with a maximum distance of about 220 kb between the AACT/PCI and PI genes. These genes form a PI cluster at 14q32.1, similar to that of the homologous genes on murine chromosome l2. The close proximity of these genes has implications for disease-association studies.     

Comparative human-mouse-rat sequence analysis of the ICAM gene cluster on HSA 19p13.2 and a 185-kb porcine region from SSC 2q

The human intercellular adhesion molecule gene (ICAM) cluster is located in a GC-rich and gene-rich region on HSA 19p13.2. We determined the complete DNA sequence of a 185-kb porcine bacterial artificial chromosome (BAC) clone containing parts of the ICAM gene cluster. We used the porcine sequence for a detailed comparative analysis between human, pig, mouse and rat. The 185 kb of porcine sequence covered 220 kb of homologous sequence in the human genome, which adds to the growing evidence that the porcine genome is somewhat smaller than the human genome. The genomic sequences of the four species showed a high level of conserved synteny and no rearrangements in gene order were observed. During evolution, the ICAM3 gene was inactivated by mutation in the mouse and rat genome, whereas it is still present in the human and pig genome. The loss of Icam3 in rodent genomes might be relevant for rodent-specific properties of the T-cell-mediated immune response. All the other investigated genes are conserved across all four investigated sequences.     

A transcript map of an 800-kb region on human chromosome 11q13, part of the candidate region for SCA5 and BBS1

The exon-amplification method was used to identify putative transcribed sequences from an 800-kb region that includes the genes for phospholipase Cβ3 and PYGM on human chromosome 11q13. The clone contig consisted of ten cosmids, three bacterial artificial chromosomes, and one P1 artificial chromosome. A total of 83 exons were generated of which 23 were derived from known genes and expressed sequence tags (ESTs). Five different EST cDNA clones were identified and mapped on the contig. One is a homolog of the human p70S6 kinase (p70s6 k) gene whose function involves the translational regulation of ribosomal protein synthesis and thereby impacts on ribosomal biogenesis. The gene for p70s6 k is expressed universally, including within adipose cells and retina, and it could play a role in Bardet-Biedl syndrome type 1, which has been mapped to 11q13. Received: 22 July 1998 / Accepted: 24 August 1998     

FISH-mapping of a 100-kb terminal 22q13 deletion

Both cytogenetically visible and cryptic deletions of the terminal region of chromosome 22q are associated with a clinical phenotype including mental retardation, delay in expressive speech development, hypotonia, normal to accelerated growth and minor facial dysmorphic features. The genes responsible for the development of the phenotype have not yet been identified, but a distal localization is probable, since the cytogenetically visible and the cryptic deletions show a similar pattern of symptoms. We report a 33-year-old woman with a submicroscopic 22q13 deletion, mild mental retardation, speech delay, autistic symptoms and mild facial dysmorphic features. The deletion was mapped by FISH using cosmid probes from terminal 22q13, and the size of the deletion was estimated to be 100 kb. Three genes are affected by the deletion in this patient. ACR and RABL2B are deleted and proSAP2 is disrupted. This observation, together with recently published data, supports the notion that proSAP2 is the most important contributor to the 22q13 deletion phenotype.     

Characterization of five novel human genes in the 11q13-q22 region

The redundancy of sequences in dbEST has approached a level where contiguous cDNA sequences of genes can be assembled, without the need to physically handle the clones from which the ESTs are derived. This is termed EST based in silico gene cloning. With the availability of sequence chromatogram files for a subset of ESTs, the quality of EST sequences can be ascertained accurately and used in contig assembly. In this report, we performed a study using this approach and isolated five novel human genes, C11orf1-C11orf5, in the 11q13-q22 region. The full open reading frames of these genes were determined by comparison with their orthologs, of which four mouse orthologs were isolated (c11orf1, c11orf2, c11orf3 and c11orf5). These genes were then analyzed using several proteomics tools. Both C11orf1 and C11orf2 are nuclear proteins with no other distinguishing features. C11orf3 is a cytoplasmic protein containing an ATP/GTP binding site, a signal peptide located in the N-terminus and a similarity to the C. elegans protein "Probable ARP 2/3 complex 20kD subunit." C11orf4 is a peptide which displays four putative transmembrane domains and is predicted to have a cytoplasmic localization. It contains signal peptides at the N- and C-termini. C11orf5 is a putative nuclear protein displaying a central coiled coil domain. Here, we propose that this purely EST-based cloning approach can be used by modestly sized laboratories to rapidly and accurately characterize and map a significant number of human genes without the need of further sequencing.