The development of sequence-tagged sites for human chromosome 4

Association between 10 Microsatellite Markers and Nonsyndromic Cleft Lip Palate in the Chilean Population

Objective

The objective of this case-control study was to evaluate the possible association between nonsyndromic cleft lip/palate (NSCLP) and 10 genetic markers in four chromosomal regions in the admixed Spanish-Amerindian Chilean population.

Setting

Study participants included 56 patients with NSCLP identified and interviewed for positive family history during the course of clinical examinations at different rehabilitation centers in the cities of Santiago and Talca, Chile. A control group of 59 normal individuals without known familial antecedents of clefting was obtained from blood bank donors of the University Hospital, University of Chile. Cases and controls belonged to low- to low-middle socioeconomic strata.

Results

Ten markers from chromosome 4p, 4q, 6p, 17q, and 19q were assessed (MSX1, D4S175, D4S192, F13A1, EDN1, D6S89, D6S105, D6S109, D17S579, BCL3). Four of them showed significant deviations from Hardy-Weinberg expectations in controls, according to the exact test (D4S192, BCL3, F13A1, and D6S89). The case-control comparison by means of the CLUMP program showed significant differences only in BCL3, and D6S109 almost reached statistical significance.

Conclusions

Most of the genetic regions with positive results in Caucasian populations may not be involved in NSCLP in Chile, regardless of the positive evidence for the candidate region on chromosome 19. Similar findings have been reported recently in the Chinese population.

Preimplantation genetic diagnosis for Huntington's disease with exclusion testing

Huntington's disease is an autosomal dominant, late-onset disorder, for which the gene and the causative mutation have been known since 1993. Some at-risk patients choose for presymptomatic testing and can make reproductive choices accordingly. Others however, prefer not to know their carrier status, but may still wish to prevent the birth of a carrier child. For these patients, exclusion testing after prenatal sampling has been an option for many years. A disadvantage of this test is that unaffected pregnancies may be terminated if the parent at risk (50%) has not inherited the grandparental Huntington gene, leading to serious moral and ethical objections. As an alternative, preimplantation genetic diagnosis (PGD) on embryos obtained in vitro may be proposed, after which only embryos free of risk are replaced. Embryos can then be selected, either by the amplification of the CAG repeat in the embryos without communicating results to the patients (ie non-disclosure testing), which brings its own practical and moral problems, or exclusion testing. We describe here the first PGD cycles for exclusion testing for Huntington's disease in five couples. Three couples have had at least one PGD cycle so far. One pregnancy ensued and a healthy female baby was delivered.

Association of the tissue kallikrein gene promoter with ESRD and hypertension

BACKGROUND

Kallikreins have long been implicated in human essential hypertension and associated complications. In particular, low urinary kallikrein excretion has been associated with hypertension and renal disease in African Americans. In an effort to identify the source of differential kallikrein excretion, we investigated the promoter of KLK1, the tissue kallikrein gene. The KLK1 promoter is uniquely polymorphic with a poly-G length polymorphism coupled with multiple single base substitutions. In this report, we genetically evaluated the association of KLK1 gene promoter alleles with end-stage renal disease (ESRD) in African Americans.

METHODS

A total of 15 haplotypes were identified in the KLK1 promoter region through detailed DNA sequence analysis. This polymorphic region was then genetically evaluated for association with ESRD in African Americans with diabetic and non-diabetic etiologies of ESRD.

RESULTS

The complex polymorphic nature of the promoter presents challenges to determining the alleles. We have redefined the region as six separate loci: five substitution loci and one length locus. The length locus was defined as G repeats starting at position -130 and ending at -121 on the gene. Among four relevant substitution loci for this study, one at position -131, just outside the G repeats, is an A-to-G substitution. The other three variant positions are -129, -128, and -127, all G-to-C substitutions within the G repeats. This region was genotyped in African American subjects with and without ESRD using semiautomated sequencing. Four different G repeat alleles ranging from 11.8% for 12 Gs to 52.3% for 10 Gs were observed in 86 control subjects. The C substitution of Gs ranges from 2.9% at position -127 to 8.2% at -129. When affected probands from each of 76 hypertensive ESRD families were genotyped, an association for the 12 G allele, the longest of the length locus alleles, was detected (allele specific P = 0.004 and locus total P = 0.02). When all ESRD affected individuals with hypertension from each family (107 patients in total) were used in the analysis, an even stronger association was observed for this allele (allele specific P = 0.003, locus total P = 0.01). This allele was more frequent in the hypertensive (non-diabetic) patients (0.20 in probands and 0.19 in all ESRD cases) than in the controls (0.12). No evidence of association in diabetic ESRD patients was observed (P = 0.93).

CONCLUSIONS

The KLK1 promoter is uniquely polymorphic. The observed genetic association suggests an etiologic effect of the KLK1 promoter on hypertension and/or hypertension associated ESRD.

Physical and Meiotic Mapping of the Region of Human Chromosome 4q11–q13 Encompassing the Vitamin D Binding Protein DBP/Gc-Globulin and Albumin Multigene Cluster

The vitamin D binding protein/Gc-globulin (DBP) gene is a member of a multigene cluster that includes albumin (ALB), α-fetoprotein (AFP), and α-albumin/afamin (AFM). All four genes have structural and functional similarities and map to the same chromosomal regions in humans (4q11–q13), mice, and rats. An accurate physical map of the region encompassing these genes is a prerequisite for study of their respective transcriptional regulation and identification of potential shared regulatory elements. By refining the physical and meiotic maps of the 4q11–q13 region and creating a local PAC contig, the order and transcriptional orientations of these four genes were determined to be centromere–3′-DBP-5′–5′-ALB-3′–5′-AFP-3′–5′-AFM3′–telomere. The ancestral DBP gene was separated from the ALBgene by >1.5 Mb. This organization and spacing establishes a foundation for ongoing functional studies in this region.

Identification of human plasma kallikrein gene polymorphisms and evaluation of their role in end-stage renal disease

Kallikreins are serine proteases that release kinins from kininogens. Kinins, via their effects on cardiovascular and renal function, may be involved in the pathogenesis of hypertension and renal failure. Two groups of kallikreins exist, glandular or tissue kallikrein and plasma kallikrein. In this study, we examined the human plasma kallikrein gene KLK3 to determine whether it contributed to end-stage renal disease (ESRD) susceptibility. We identified two novel polymorphic sequences closely linked to the KLK3 gene, designated KLK3b and KLK3c (heterozygosities: 0.64 to 0.68 and 0.48 to 0.52, respectively). We mapped the KLK3 gene and the marker KLK3c to the long arm of human chromosome 4 between F11 and D4S426 using a radiation hybrid panel. The study population consisted of 142 sibling pairs concordant for ESRD from 121 African American families. The 142 sibling pairs were stratified into 78 pairs with hypertension- and chronic glomerulonephritis-associated ESRD and 64 with non-insulin-dependent diabetes mellitus-associated ESRD. Linkage analyses, using SIBPAL of SAGE, and exclusion analysis, using MAPMAKERS/SIBS, were performed. Linkage analysis of affected sibling pairs did not reveal any evidence of linkage of KLK3 to ESRD in all 142 sib-pairs or in the two stratified subsets. Exclusion analysis indicated that the KLK3 gene could be excluded from contributing to ESRD at a relative risk of 3 when the maximum log of the odds score of -2 was used as the criterion for exclusion. However, an association analysis using the relative predispositional effect technique showed that alleles 7 and 9 of KLK3b were consistently associated with ESRD. Alleles 7 and 9 were present in 11.2% and 10.8% of the 113 unrelated ESRD probands and in 6.6% and 6.6% of the 204 race-matched control subjects without renal disease (allele P=.0041 and .0016, respectively). Alleles 7 and 9 were also present in 13% and 10.4% of the proband's first siblings (allele P=.00014 and .0087, respectively). The association of KLK3b alleles with ESRD raises the possibility that polymorphisms in KLK3 may play a role in ESRD susceptibility. The lack of linkage might reflect our relatively small family set.

Distribution of olfactory receptor genes in the human genome

We demonstrate that members of the olfactory receptor (OR) gene family are distributed on all but a few human chromosomes. Through FISH analysis, we show that OR sequences reside at more than 25 locations in the human genome. Their distribution is biased for terminal bands. Flow-sorted chromosomes were used to isolate 87 OR sequences derived from 16 chromosomes. Their sequence-relationships are indicative of the inter- and intrachromosomal duplications responsible for OR family expansion. The human genome has accumulated a striking number of dysfunctional copies: 72% of the sequences are pseudogenes. ORF-containing sequences predominate on chromosomes 7, 16 and 17.

Identification of a novel isoform of mouse dentin matrix protein 1: spatial expression in mineralized tissues

Dentin matrix protein 1 (Dmp1) is an acidic phosphoprotein first identified by cDNA cloning from a rat tooth library. Northern blot hybridization of a variety of tissues detected Dmp1 mRNAs only in odontoblasts, suggesting that this protein was odontoblast specific. In situ hybridization studies showed expression of Dmp1 in odontoblasts with transient expression in secretory ameloblasts. The purpose of this study was to isolate and characterize a mouse Dmp1 cDNA and determine its spatial expression pattern related to other mineralizing tissues. A mouse molar cDNA library was screened with a 32P-labeled Dmp1 polymerase chain reaction amplification product in order to isolate a full-length clone. DNA sequence analysis of the largest mouse Dmp1 cDNA (2802 base pairs [bp]) revealed an open reading frame of 1509 nucleotides encoding a 503 amino acid protein with a single polyadenylation signal. Comparison with rat and bovine Dmp1 sequence showed high homology and the identification of a 45 bp (15 amino acid) insert, representing an alternative spliced mRNA. This 45 bp segment was shown to represent a small exon by DNA analysis of a mouse genomic Dmp1 clone. In situ hybridization studies revealed a much broader Dmp1 tissue expression pattern than previously reported. Dmp1 transcripts were detected in the odontoblast and ameloblasts, osteoblasts, and cementoblasts. Our data indicate that Dmp1 is alternatively spliced, and the primary full-length transcript contains a 45 bp insert which is encoded by a small exon. Therefore, Dmp1 is not a tooth-specific protein but rather is expressed in a number of mineralizing tissues including enamel, bone, and cementum.

Members of the olfactory receptor gene family are contained in large blocks of DNA duplicated polymorphically near the ends of human chromosomes

We have identified three new members of the olfactory receptor (OR) gene family within a large segment of DNA that is duplicated with high similarity near many human telomeres. This segment is present at 3q, 15q, and 19p in each of 45 unrelated humans sampled from various populations. Additional copies are present polymorphically at 11 other subtelomeric locations. The frequency with which the block is present at some locations varies among populations. While humans carry seven to 11 copies of the OR-containing block, it is located in chimpanzee and gorilla predominantly at a single site, which is not orthologous to any of the locations in the human genome. The observation that sequences flanking the OR-containing segment are duplicated on larger and different sets of chromosomes than the OR block itself demonstrates that the segment is part of a much larger, complex patchwork of subtelomeric duplications. The population analyses and structural results suggest the types of processes that have shaped these regions during evolution. From its sequence, one of the OR genes in this duplicated block appears to be potentially functional. Our findings raise the possibility that functional diversity in the OR family is generated in part through duplications and inter-chromosomal rearrangements of the DNA near human telomeres.

Molecular evolution of an ancient mariner transposon, Hsmar1, in the human genome

A confident consensus sequence for Hsmar1, the first mariner transposon recognized in the human genome, was generated using three genomic and 15 cDNA sequences. It is thought to represent the ancestrally active copy that invaded an early primate genome. The consensus is 1287 base pairs (bp) long, has 30 bp perfect inverted terminal repeats (ITRs), and encodes a 343 amino acid (aa) mariner transposase. Each copy has diverged from the consensus largely independently of the others and mostly neutrally, and most are now defective. They differ from the consensus by an average of 7.8% in DNA sequence and 7.5 indels per kilobase, both of which values indicate that the copies were formed about 50 Myr ago. On average, only 20% of the 73 surmised CpG hypermutable sites in the consensus remain. A remarkable exception to this loss of functionality is revealed by a set of ten cDNA clones derived from a particular genomic copy that has diverged only 2.4% from the consensus, retained 54% of its hypermutable CpG pairs, and which has a full-length transposase open reading frame. The complete sequence of one of these cDNAs (NIB1543) indicates that the transposase gene of this copy may have been conserved because it is spliced to a human cellular gene encoding a SET domain protein. A specific PCR assay was used to reveal the presence of Hsmar1 copies in all primates examined representing all major lineages, but not in close relatives of primates. PCR fragments cloned and sequenced from a representative sample of primates confirmed that Hsmar1 copies are present in all major lineages, and also revealed another cecropia subfamily mariner in prosimians only, and a third highly divergent mariner present in the greater slow loris Nycticebus coucang. There are about 200 copies of Hsmar1 in the human genome, as well as +/-2400 copies of a derived 80 bp paired ITR structure and +/-4600 copies of solo ITRs. Thus, this transposon had a considerable insertional mutagenic effect on past primate genomes.

High-resolution genetic map of the human glutamyl aminopeptidase gene (ENPEP)

The murine B-lymphocyte differentiation antigen BP-1/6C3 has been identified as glutamyl aminopeptidase (EAP), the gene symbol for which is ENPEP. Using genomic DNA encoding for human EAP as a probe, we identified the ENPEP gene location on human chromosome 4q25 by polymerase chain reaction analysis of a human/rodent somatic cell hybrid mapping panel and by fluorescence in situ hybridization. Using a radiation hybrid panel, the gene order around ENPEP was determined to be centromere-D4S1236-(570 kb)-ENPEP-(210 kb)-D4S262-(270 kb)-D4S953-(270 kb)-D4S474-(570 kb)-IF. The linkage of ENPEP to complement factor I (IF) confirms the human chromosome band 4q25 localization predicted from the chromosomal location of murine ENPEP. Human ENPEP thus provides an additional marker for the long arm of chromosome 4 that should facilitate studies of this genomic region.

A novel tandem repeat sequence located on human chromosome 4p: isolation and characterization

In an effort to analyze the genomic region of the distal half of human chromosome 4p, to where Huntington disease and other diseases have been mapped, we have isolated the cosmid clone (CRS447) that was likely to contain a region with specific repeat sequences. Clone CRS447 was subjected to detailed analysis, including chromosome mapping, restriction mapping, and DNA sequencing. Chromosome mapping by both a human-CHO hybrid cell panel and FISH revealed that CRS447 was predominantly located in the 4p15.1-15.3 region. CRS447 was shown to consist of tandem repeats of 4.7-kb units present on chromosome 4p. A single EcoRI unit was subcloned (pRS447), and the complete sequence was determined as 4752 nucleotides. When pRS447 was used as a probe, the number of copies of this repeat per haploid genome was estimated to be 50-70. Sequence analysis revealed that it contained two internal CA repeats and one putative ORF. Database search established that this sequence was unreported. However, two homologous STS markers were found in the database. We concluded that CRS447/pRS447 is a novel tandem repeat sequence that is mainly specific to human chromosome 4p.

An STS-based radiation hybrid map of the human genome

We have constructed a physical map of the human genome by using a panel of 83 whole genome radiation hybrids (the Stanford G3 panel) in conjunction with 10,478 sequence-tagged sites (STSs) derived from random genomic DNA sequences, previously mapped genetic markers, and expressed sequences. Of these STSs, 5049 are framework markers that fall into 1766 high-confidence bins. An additional 945 STSs are indistinguishable in their map location from one or more of the framework markers. These 5994 mapped STSs have an average spacing of 500 kb. An additional 4484 STSs are positioned with respect to the framework markers. Comparison of the orders of markers on this map with orders derived from independent meiotic and YAC STS-content maps indicates that the error rate in defining high-confidence bins is < 5%. Analysis of 322 random cDNAs indicates that the map covers the vast majority of the human genome. This STS-based radiation hybrid map of the human genome brings us one step closer to the goal of a physical map containing 30,000 unique ordered landmarks with an average marker spacing of 100 kb.

A collection of 1814 human chromosome 7-specific STSs

An established goal of the ongoing Human Genome Project is the development and mapping of sequence-tagged sites (STSs) every 100 kb, on average, across all human chromosomes. En route to constructing such a physical map of human chromosome 7, we have generated 1814 chromosome 7-specific STSs. The corresponding PCR assays were designed by the use of DNA sequence determined in our laboratory (79%) or generated elsewhere (21%) and were demonstrated to be suitable for screening yeast artificial chromosome (YAC) libraries. This collection provides the requisite landmarks for constructing a physical map of chromosome 7 at < 100-kb average spacing of STSs.

The carboxyl-terminal domain of phosphophoryn contains unique extended triplet amino acid repeat sequences forming ordered carboxyl-phosphate interaction ridges that may be essential in the biomineralization process

Phosphophoryns (PPs), a family of Asp and Ser(P)-rich dentin proteins, are considered to be archetypal regulators of several aspects of extracellular matrix (ECM) biomineralization. We have cloned a rat incisor PP gene, Dmp2, from our odontoblast cDNA library and localized it to mouse chromosome 5q21 within 2 centimorgans of Dmp1, another tooth-specific ECM protein. The carboxyl-terminal region of Dmp2 protein (60 residue % Ser, 31 residue % Asp) is divided into two domains, one with unique repetitive blocks of [DSS]n,3</=14, the other with [SD]m = 2,3. Conformational analysis shows the phosphorylated form of the [DS*S*]n repeats to have a unique structure with well defined ridges of phosphates and carboxyls available for counter ion binding. The [S*D]m domains have different phosphate and carboxylate interaction edges and thus different calcium ion and apatite surface binding properties. These two domains and the colocalization of Dmp1 and Dmp2 genes at a position equivalent to the dentinogenesis imperfecta type II location on human 4q21 all suggest that the PPs are indeed involved in some aspect of ECM mineralization.

Simultaneous assessment of loss of heterozygosity at multiple microsatellite loci using semi-automated fluorescence-based detection: subregional mapping of chromosome 4 in cervical carcinoma

Detection of loss of heterozygosity (LOH) by comparison of normal and tumor genotypes using PCR-based microsatellite loci provides considerable advantages over traditional Southern blotting-based approaches. However, current methodologies are limited by several factors, including the numbers of loci that can be evaluated for LOH in a single experiment, the discrimination of true alleles versus "stutter bands," and the use of radionucleotides in detecting PCR products. Here we describe methods for high throughput simultaneous assessment of LOH at multiple loci in human tumors; these methods rely on the detection of amplified microsatellite loci by fluorescence-based DNA sequencing technology. Data generated by this approach are processed by several computer software programs that enable the automated linear quantitation and calculation of allelic ratios, allowing rapid ascertainment of LOH. As a test of this approach, genotypes at a series of loci on chromosome 4 were determined for 58 carcinomas of the uterine cervix. The results underscore the efficacy, sensitivity, and remarkable reproducibility of this approach to LOH detection and provide subchromosomal localization of two regions of chromosome 4 commonly altered in cervical tumors.

Gene-based sequence-tagged-sites (STSs) as the basis for a human gene map

Using our data set of 3,143 single pass sequences from human brain cDNA libraries, we have developed a strategy in which gene-based sequence-tagged-sites (STSs), derived from 3'untranslated regions of human cDNAs, are rapidly assigned to megabase-insert yeast artificial chromosomes and somatic cell hybrids to generate regional gene mapping data. Employing this approach, we have mapped 318 cDNAs, representing 308 human genes. Ninety-two of these mapped to regions implicated in human genetic diseases, identifying them as candidate genes. Extension of this strategy has the potential to result in virtually every human gene having, at its 3' end, its own associated STS, with each STS in turn specifying both a corresponding genomic clone and a specific regional location in the genome.

Summary of the 1993 ASHG ancillary meeting "recent research on chromosome 4p syndromes and genes"

The following is a summary of presentations given during an ancillary meeting to the 1993 American Society of Human Genetics Meeting in New Orleans, LA. This ancillary meeting, entitled "Recent Research on Chromosome 4p Syndromes and Genes," reviewed the history of the Wolf-Hirschhorn syndrome (WHS), the natural history of patients with WHS, and the smallest region of deletion associated with the WHS. The proximal 4p deletion syndrome and the duplication 4p syndrome were also described and advice was offered regarding detection of chromosome 4p deletions, duplications, and rearrangements. The current status of the physical map of chromosome 4p with emphasis on the genes that map to the 4p16 region was presented along with a preliminary phenotypic map of 4p16. The goal of this format was to provide a comprehensive review of the clinical presentations, diagnostic capabilities, and genetic mapping advances involving chromosome 4p.

Use of porcine interspersed repeat sequences in PCR-mediated genotyping

PCR primers derived from porcine short and long interspersed repeat sequences were used to amplify DNA samples isolated from individual members of three-generation pig reference pedigrees. Subsequent high-resolution gel electrophoresis of both SINE and LINE-PCR products allowed direct visualisation of polymorphisms that segregated in a Mendelian manner. Additional polymorphisms were detected by Southern blotting of the gels described above followed by hybridization with simple sequence DNA. Genotyping by interspersed repeat-PCR exploits the natural architecture of the pig genome and allows the typing of polymorphisms by utilizing pre-existing sequence information.

Toward the complete genomic map and molecular pathology of human chromosome 4

The identification of disease genes via molecular DNA cloning has revolutionized human genetics and medicine. Both the candidate gene approach and positional cloning have been used successfully. The defects causing Huntington's disease, facioscapulohumeral muscular dystrophy, piebaldism, Hurler/Scheie syndrome, one form of autosomal recessive retinitis pigmentosa, and a second locus for autosomal dominant polycystic kidney disease have recently been localized to chromosome 4. In addition to the rapid progress in the cloning of the 203-megabase chromosome, the presence of more than 60 closely spaced microsatellites on this chromosome will undoubtedly lead to the localization of additional disease genes. In order to consider cloned genes as potential candidates for disorders assigned to chromosome 4, it is important to collect and order all genes with respect to their chromosomal localization. Analysis of cytogenetically visible interstitial and terminal deletions should also be helpful in defining new disease gene loci and in mapping novel genes. These data represent the status quo of the integrated molecular map for chromosome 4.

Seven loci on human chromosome 4 map onto sheep chromosome 6: a proposal to restore the original nomenclature of this sheep chromosome

Seven new loci, casein alpha-S1 (CSN1S1), casein alpha-S2 (CSN1S2), casein beta (CSN2), the Hardy-Zuckerman 4 feline sarcoma viral (v-kit) oncogene homolog (KIT), albumin (ALB), phosphodiesterase cyclic GMP (rod receptor) beta polypeptide (PDEB), and complement component 1 (IF), were assigned to sheep Chromosome (Chr) 6 by Southern hybridization to a panel of chromosomally characterized sheep x hamster cell hybrids. By isotopic in situ hybridization, CSN2 was regionally localized to sheep Chr (OOV) 6q22-q31, anchoring this syntenic group of markers on to OOV6 and confirming its homology at a molecular and cytological level with cattle Chr 6. The assignment of these loci, from PDEB (located on human Chr 4p16.3) to IF (on HSA4q24-q25), and the observation that interleukin 2 (IL2, on HSA4q26-q27) and tryptophan 2,3-dioxygenase (TDO2, on HSA4q31) are not located on OOV6, is further evidence of the close evolutionary relationship of sheep and cattle and the conserved synteny in these species of this extensive region of human Chr 4. On the basis of this conserved synteny, and the similar G- and Q-banding patterns of this chromosome in cattle and sheep, we propose that this sheep chromosome be numbered as 6, not 4 as recommended by ISCNDA (1990).

The 1993-94 Généthon human genetic linkage map

In 1992, we described a second-generation genetic linkage map of the human genome. Using 1,267 new microsatellite markers, we now present a new genetic linkage map containing a total of 2,066 (AC)n short tandem repeats, 60% of which show a heterozygosity of over 0.7. Statistical linkage analysis based on the genotyping of eight large CEPH families placed these markers in the 23 linkage groups. The map includes 1,266 intervals and spans a total distance of 3690 centiMorgans (cM). A total of 1,041 markers could be ordered with odds ratios greater than 1000:1. About 56% of this map is at a distance of 1 cM or less from one of its markers.

Chromosome 4 localization of a second gene for autosomal dominant polycystic kidney disease

Autosomal dominant polycystic kidney disease (ADPKD) is a genetically heterogeneous disorder. A gene defect located on the short arm of chromosome 16 is responsible for the disease in roughly 86% of affected European families. Using highly polymorphic microsatellite DNA markers, we have assigned a second gene for ADPKD to chromosome 4. In eight families with clear evidence against linkage to chromosome 16 markers, linkage analysis with the markers D4S231 and D4S423, demonstrated a multipoint lod score of 22.42.