MethylScreen: DNA methylation density monitoring using quantitative PCR.

Aberrant gene silencing of genes through cytosine methylation has been demonstrated during the development of many types of cancers including prostate cancer Several genes including GSTP1 have been shown to be methylated in prostate cancer leading to the suggestion and demonstration that methylation status of such genes could be used as cancer diagnosis markers alone or in support of histology. We developed a bisulfite-free alternative, MethylScreen technology, an assay for DNA methylation detection utilizing combined restriction from both methylation-sensitive restriction enzymes (MSRE) and methylation-dependent restriction enzymes (MDRE). MethylScreen was used to analyze the 5' region of GSTP1 in cell lines, in vitro methylated DNA populations, and flash-frozen tissue samples in an effort to characterize the output and analytical performance characteristics of the assay. The output from the quantitative PCR assay suggested that it could not only detect fully methylated molecules in a mixed population below the 1% level, but it could also quantify the abundance of intermediately methylated molecules. Interestingly, the interpreted output from the four quantitative PCRs closely resembled the molecular population as described by clone-based bisulfite genomic sequencing.

Comprehensive DNA methylation profiling in a human cancer genome identifies novel epigenetic targets.

Using a unique microarray platform for cytosine methylation profiling, the DNA methylation landscape of the human genome was monitored at more than 21,000 sites, including 79% of the annotated transcriptional start sites (TSS). Analysis of an oligodendroglioma derived cell line LN-18 revealed more than 4000 methylated TSS. The gene-centric analysis indicated a complex pattern of DNA methylation exists along each autosome, with a trend of increasing density approaching the telomeres. Remarkably, 2% of CpG islands (CGI) were densely methylated, and 17% had significant levels of 5 mC, whether or not they corresponded to a TSS. Substantial independent verification, obtained from 95 loci, suggested that this approach is capable of large scale detection of cytosine methylation with an accuracy approaching 90%. In addition, we detected large genomic domains that are also susceptible to DNA methylation reinforced inactivation, such as the HOX cluster on chromosome 7 (CH7). Extrapolation from the data suggests that more than 2000 genomic loci may be susceptible to methylation and associated inactivation, and most have yet to be identified. Finally, we report six new targets of epigenetic inactivation (IRX3, WNT10A, WNT6, RARalpha, BMP7 and ZGPAT). These targets displayed cell line and tumor specific differential methylation when compared with normal brain samples, suggesting they may have utility as biomarkers. Uniquely, hypermethylation of the CGI within an IRX3 exon was correlated with over-expression of IRX3 in tumor tissues and cell lines relative to normal brain samples.

Using a pericentromeric interspersed repeat to recapitulate the phylogeny and expansion of human centromeric segmental duplications.

Despite considerable advances in sequencing of the human genome over the past few years, the organization and evolution of human pericentromeric regions have been difficult to resolve. This is due, in part, to the presence of large, complex blocks of duplicated genomic sequence at the boundary between centromeric satellite and unique euchromatic DNA. Here, we report the identification and characterization of an approximately 49-kb repeat sequence that exists in more than 40 copies within the human genome. This repeat is specific to highly duplicated pericentromeric regions with multiple copies distributed in an interspersed fashion among a subset of human chromosomes. Using this interspersed repeat (termed PIR4) as a marker of pericentromeric DNA, we recovered and sequence-tagged 3 Mb of pericentromeric DNA from a variety of human chromosomes as well as nonhuman primate genomes. A global evolutionary reconstruction of the dispersal of PIR4 sequence and analysis of flanking sequence supports a model in which pericentromeric duplications initiated before the separation of the great ape species (>12 MYA). Further, analyses of this duplication and associated flanking duplications narrow the major burst of pericentromeric duplication activity to a time just before the divergence of the African great ape and human species (5 to 7 MYA). These recent duplication exchange events substantially restructured the pericentromeric regions of hominoid chromosomes and created an architecture where large blocks of sequence are shared among nonhomologous chromosomes. This report provides the first global view of the series of historical events that have reshaped human pericentromeric regions over recent evolutionary time.

A map of human genome sequence variation containing 1.42 million single nucleotide polymorphisms.

We describe a map of 1.42 million single nucleotide polymorphisms (SNPs) distributed throughout the human genome, providing an average density on available sequence of one SNP every 1.9 kilobases. These SNPs were primarily discovered by two projects: The SNP Consortium and the analysis of clone overlaps by the International Human Genome Sequencing Consortium. The map integrates all publicly available SNPs with described genes and other genomic features. We estimate that 60,000 SNPs fall within exon (coding and untranslated regions), and 85% of exons are within 5 kb of the nearest SNP. Nucleotide diversity varies greatly across the genome, in a manner broadly consistent with a standard population genetic model of human history. This high-density SNP map provides a public resource for defining haplotype variation across the genome, and should help to identify biomedically important genes for diagnosis and therapy.

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.

A physical map of the human Y chromosome.

The non-recombining region of the human Y chromosome (NRY), which comprises 95% of the chromosome, does not undergo sexual recombination and is present only in males. An understanding of its biological functions has begun to emerge from DNA studies of individuals with partial Y chromosomes, coupled with molecular characterization of genes implicated in gonadal sex reversal, Turner syndrome, graft rejection and spermatogenic failure. But mapping strategies applied successfully elsewhere in the genome have faltered in the NRY, where there is no meiotic recombination map and intrachromosomal repetitive sequences are abundant. Here we report a high-resolution physical map of the euchromatic, centromeric and heterochromatic regions of the NRY and its construction by unusual methods, including genomic clone subtraction and dissection of sequence family variants. Of the map's 758 DNA markers, 136 have multiple locations in the NRY, reflecting its unusually repetitive sequence composition. The markers anchor 1,038 bacterial artificial chromosome clones, 199 of which form a tiling path for sequencing.

The syntenic relationship of the zebrafish and human genomes.

The zebrafish is an important vertebrate model for the mutational analysis of genes effecting developmental processes. Understanding the relationship between zebrafish genes and mutations with those of humans will require understanding the syntenic correspondence between the zebrafish and human genomes. High throughput gene and EST mapping projects in zebrafish are now facilitating this goal. Map positions for 523 zebrafish genes and ESTs with predicted human orthologs reveal extensive contiguous blocks of synteny between the zebrafish and human genomes. Eighty percent of genes and ESTs analyzed belong to conserved synteny groups (two or more genes linked in both zebrafish and human) and 56% of all genes analyzed fall in 118 homology segments (uninterrupted segments containing two or more contiguous genes or ESTs with conserved map order between the zebrafish and human genomes). This work now provides a syntenic relationship to the human genome for the majority of the zebrafish genome.

Radiation hybrid mapping of the zebrafish genome.

The zebrafish is an excellent genetic system for the study of vertebrate development and disease. In an effort to provide a rapid and robust tool for zebrafish gene mapping, a panel of radiation hybrids (RH) was produced by fusion of irradiated zebrafish AB9 cells with mouse B78 cells. The overall retention of zebrafish sequences in the 93 RH cell lines that constitute the LN54 panel is 22%. Characterization of the LN54 panel with 849 simple sequence length polymorphism markers, 84 cloned genes and 122 expressed sequence tags allowed the production of an RH map whose total size was 11,501 centiRays. From this value, we estimated the average breakpoint frequency of the LN54 RH panel to correspond to 1 centiRay = 148 kilobase. Placement of a group of 235 unbiased markers on the RH map suggests that the map generated for the LN54 panel, at present, covers 88% of the zebrafish genome. Comparison of marker positions in RH and meiotic maps indicated a 96% concordance. Mapping expressed sequence tags and cloned genes by using the LN54 panel should prove to be a valuable method for the identification of candidate genes for specific mutations in zebrafish.

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 (designated CMT2D). 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 approximately 3.4 Mb have been assembled across the combined CMT2D 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.

Promoter analysis of the human centromeric and telomeric survival motor neuron genes (SMNC and SMNT).

Proximal spinal muscular atrophy (SMA) is caused by mutations in the telomeric (SMNT), but not centromeric (SMNC), survival motor neuron gene. Here we have identified and analyzed the two SMN promoters. We show that a 750-bp 5'-flanking fragment from each is capable of driving expression from a reporter construct. Within this fragment, we define a approximately 200-bp element that results in high expression in a motor neuron cell line. Sequence comparison of a 3. 4-kb upstream fragment from each gene shows minimal differences. Although these differences produce a 2-fold difference in reporter activity between the two promoters, this is not sufficiently high to explain why SMNT, but not SMNC, is the disease determining gene. Our data thus demonstrate, for the first time, almost complete equivalence between the SMN promoters and rule out the important possibility that differences in them might explain why mutations in only the telomeric SMN gene cause SMA.

A single nucleotide difference that alters splicing patterns distinguishes the SMA gene SMN1 from the copy gene SMN2.

Spinal muscular atrophy (SMA) is a recessive disorder characterized by loss of motor neurons in the spinal cord. It is caused by mutations in the telomeric survival motor neuron 1 ( SMN1 ) gene. Alterations within an almost identical copy gene, the centromeric survival motor neuron 2 ( SMN2 ) gene produce no known phenotypic effect. The exons of the two genes differ by just two nucleotides, neither of which alters the encoded amino acids. At the genomic level, only five nucleotides that differentiate the two genes from one another have been reported. The entire genomic sequence of the two genes has not been determined. Thus, differences which might explain why SMN1 is the SMA gene are not readily apparent. In this study, we have completely sequenced and compared genomic clones containing the SMN genes. The two genes show striking similarity, with the homology being unprecedented between two different yet functional genes. The only critical difference in an approximately 32 kb region between the two SMN genes is the C->T base change 6 bp inside exon 7. This alteration but not other variations in the SMN genes affects the splicing pattern of the genes. The majority of the transcript from the SMN1 locus is full length, whereas the majority of the transcript produced by the SMN2 locus lacks exon 7. We suggest that the exon 7 nucleotide change affects the activity of an exon splice enhancer. In SMA patients, the loss of SMN1 but the presence of SMN2 results in low levels of full-length SMN transcript and therefore low SMN protein levels which causes SMA.

Molecular refinement of the 1p36 deletion syndrome reveals size diversity and a preponderance of maternally derived deletions.

The deletion of chromosome 1p36 is a newly recognized, relatively common contiguous gene deletion syndrome with a variable phenotype. The clinical features have recently been delineated and molecular analysis indicates that the prevalence of certain phenotypic features appears to correlate with deletion size. Phenotype/genotype comparisons have allowed the assignment of certain clinical features to specific deletion intervals, significantly narrowing the regions within which to search for candidate genes. We have extensively characterized the deletion regions in 30 cases using microsatellite markers and fluorescence in situ hybridization analyses. The map order of 28 microsatellite markers spanning the deletion region was obtained by a combination of genotypic analysis and physical mapping. The deletion region was divided into six intervals and breakpoints were found to cluster in mainly two regions. Molecular analysis of the deletions showed that two patients had complex re-arrangements; these cases shared their distal and proximal breakpoints in the two common breakpoint regions. Of the de novo deletions ( n = 28) in whichparental samples were available and the analysis was informative ( n = 27), there were significantly morematernally derived deletions ( n = 21) than paternally derived deletions ( n = 6) (chi1(2) = 8.35, P < 0.0001). Phenotype/genotype correlations and refinements of critical regions in our naturally occurring deletion panel have delineated specific areas in which to focus the search for the causative genes for the features of this syndrome.

Molecular cytogenetic delineation of the critical deleted region in the 5q- syndrome.

The 5q- syndrome is a distinct type of myelodysplastic syndrome (MDS) characterised by refractory anaemia, morphological abnormalities of megakaryocytes, and del(5q) as the sole cytogenetic abnormality. In contrast to patients with therapy-related MDS with 5q deletions, 5q- syndrome patients have a favourable prognosis and a low rate of transformation to acute leukaemia. We have previously delineated a common deleted region of 5.6 Mb between the gene for fibroblast growth factor acidic (FGF1) and the subunit of interleukin 12 (IL12B) in two patients with 5q- syndrome and small deletions, del(5)(q31q33). The present study used fluorescence in situ hybridisation (FISH) analysis of these and a third 5q- syndrome patient with a small deletion, del(5)(q33q34), to refine further the critical deleted region. This resulted in the narrowing of the common deleted region within 5q31.3-5q33 to approximately 3 Mb, flanked by the adrenergic receptor beta 2 (ADRB2) and IL/2B genes. The common region of loss in these three 5q- syndrome patients includes the macrophage colony-stimulating factor-1 receptor (CSF1R), secreted protein, acidic, cysteine-rich (SPARC), and glutamate receptor (GR1A1) genes. This 5q- syndrome critical region is telomeric to and distinct from the other critical regions on 5q associated with MDS and acute myeloid leukaemia.

Molecular cloning and mapping of human semaphorin F from the Cri-du-chat candidate interval.

Cri-du-chat is a human contiguous gene deletion syndrome resulting from hemizygous deletions of chromosome 5p. Here we describe the isolation from within this interval of the human Semaphorin F (SEMAF) gene, a member of a family of proteins that has been implicated in axonal pathfinding. The human SEMAF gene covers at least 10% of the deleted region and defines a new class within this large gene family characterized by the presence of seven type 1 thrombospondin repeats. Prominent expression of murine semaphorin F (Semaf) was observed in the mouse brain, consistent with a role for semaphorin F as a signaling molecule that guides axons or migrating neuronal precursors during development. The known functions of semaphorins and the interesting pattern of expression for Semaf suggest that haploinsufficiency for SEMAF may disrupt normal brain development and might lead to some of the features of Cri-du-chat.

A radiation hybrid map of human chromosome 5 with integration of cytogenetic, genetic, and transcript maps.

One of the major goals of the human genome project is to establish a physical map of each human chromosome with a density of sequence-tagged site (STS) markers exceeding one every 100 kb. We report here the generation of a human chromosome 5-specific radiation hybrid (RH) map that includes 556 markers. Of these markers, 132 loci are ordered with a maximum likelihood ratio of >1000:1 compared with the next most likely order. An additional 113 loci were ordered relative to these backbone markers with a maximum likelihood ratio of >10:1 but <1000:1. Together, these 245 loci form an ordered framework map for the chromosome. Using this framework, >300 more markers were localized based on two-point analysis with the ordered set. On average, there are 50 markers in common with the RH map presented here and other chromosome 5 maps included in the current whole genome cytogenetic, genetic, and physical maps. The accuracy of all the maps is evident in that there are no more than two discrepancies between any one of them and these data. All of the maps encompassing chromosome 5 complement each other providing excellent STS coverage with >2200 loci combined. The chromosome 5-specific RH map contains 20% of these independent loci. In addition, our RH map contains STSs derived from clones suitable for fluorescent in situ hybridization, allowing alignment to the cytogenetic map. Together, these maps will assist in the assembly of sequence-ready contigs and will aid in the identification of disease loci on chromosome 5 by positional cloning and positional candidate approaches.

High resolution mapping of the renal sodium-phosphate cotransporter gene (NPT2) confirms its localization to human chromosome 5q35.

The precise chromosomal localization of the type II renal-specific Na+-phosphate (Pi) cotransporter (NPT2) gene (gene symbol SLC17A2) is necessary for the identification of closely linked polymorphic markers to determine whether NPT2 is a candidate gene for inherited disorders of renal Pi reabsorption. Recent studies by two different groups localized NPT2 to human chromosome 5q35 and 5q13, respectively. To resolve this discrepancy, we used three independent methods. The results using a human chromosome 5/rodent somatic cell hybrid deletion panel, fluorescence in situ hybridization with a PAC clone containing the NPT2 locus, and analysis of a chromosome 5-specific radiation hybrid panel were all consistent with the 5q35 assignment of the NPT2 gene. The radiation hybrid results placed NPT2 between polymorphic microsatellite markers D5S498 and D5S469. These findings will allow the initiation of linkage analysis to determine if NPT2 has a causative role in Mendelian disorders of renal Pi wasting.

High throughput fingerprint analysis of large-insert clones.

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