Evolution of the primate beta-globin gene region. High rate of variation in CpG dinucleotides and in short repeated sequences between man and chimpanzee

Population dynamics of GC-changing mutations in humans and great apes

The nucleotide composition of the genome is a balance between origin and fixation rates of different mutations. For example, it is well-known that transitions occur more frequently than transversions, particularly at CpG sites. Differences in fixation rates of mutation types are less explored. Specifically, recombination-associated GC-biased gene conversion (gBGC) may differentially impact GC-changing mutations, due to differences in their genomic distributions and efficiency of mismatch repair mechanisms. Given that recombination evolves rapidly across species, we explore gBGC of different mutation types across human populations and great ape species. We report a stronger correlation between segregating GC frequency and recombination for transitions than for transversions. Notably, CpG transitions are most strongly affected by gBGC in humans and chimpanzees. We show that the overall strength of gBGC is generally correlated with effective population sizes in humans, with some notable exceptions, such as a stronger effect of gBGC on non-CpG transitions in populations of European descent. Furthermore, species of the Gorilla and Pongo genus have a greatly reduced gBGC effect on CpG sites. We also study the dependence of gBGC dynamics on flanking nucleotides and show that some mutation types evolve in opposition to the gBGC expectation, likely due to hypermutability of specific nucleotide contexts. Our results highlight the importance of different gBGC dynamics experienced by GC-changing mutations and their impact on nucleotide composition evolution.

"While there is p57, there is hope." The past and the present of diagnosis in first trimester abortions: Diagnostic dilemmas and algorithmic approaches. A review

Distinction of hydatidiform moles (HM) from non-molar (NM) specimens and subclassification of HM as complete hydatidiform mole (CHM) versus partial hydatidiform mole (PHM) are important for clinical practice and investigational studies. The issue of diagnostic reproducibility is still unsolved, the lack of diagnostic accuracy based on morphology is substantial with an important interobserver variability, even between experienced gynecologic pathologists. Many ancillary techniques have been investigated in the last years to refine HM diagnosis. p57 (a paternally imprinted, maternally expressed gene) immunohistochemistry, based on the unique genetics of CHM (purely androgenetic), PHM (diandric triploid), and NM specimens (biparental, with allelic balance) can identify CHMs, which lack p57 expression because of a lack of maternal DNA. However, although its role in HM diagnosis is pivotal, it does not allow the distinction of PHM from NM specimens, both of which express p57 due to the presence of maternal DNA. Molecular genotyping, which compares villous and decidual DNA patterns to determine the parental source and ratios of polymorphic alleles, distinguishes purely androgenetic CHM from diandric triploid PHM, and both of these from NM specimens. Beyond the claim of establishing a "diagnostic truth", exceptions and peculiar genetic scenarios in the origin of rare CHM and PHM should be kept in mind when approaching any ancillary technique. An algorithmic approach, even in settings with limited resources, can help the pathologists in the diagnostic dilemma of diagnosis of first trimester abortions.

On the number of New World founders: a population genetic portrait of the peopling of the Americas

The founding of New World populations by Asian peoples is the focus of considerable archaeological and genetic research, and there persist important questions on when and how these events occurred. Genetic data offer great potential for the study of human population history, but there are significant challenges in discerning distinct demographic processes. A new method for the study of diverging populations was applied to questions on the founding and history of Amerind-speaking Native American populations. The model permits estimation of founding population sizes, changes in population size, time of population formation, and gene flow. Analyses of data from nine loci are consistent with the general portrait that has emerged from archaeological and other kinds of evidence. The estimated effective size of the founding population for the New World is fewer than 80 individuals, approximately 1% of the effective size of the estimated ancestral Asian population. By adding a splitting parameter to population divergence models it becomes possible to develop detailed portraits of human demographic history. Analyses of Asian and New World data support a model of a recent founding of the New World by a population of quite small effective size.

Gene conversion and different population histories may explain the contrast between polymorphism and linkage disequilibrium levels

To characterize linkage disequilibrium (LD) levels in human populations, we have analyzed 10 independent noncoding segments in three population samples from the major ethnic groups--that is, Africans, Asians, and Europeans. Descriptive statistics show that LD decays much faster in the African samples than in the non-African ones. With the assumption of an equilibrium model, we estimated the population crossing-over parameter (4N(e)r(bp), where N(e) is the effective population size and r(bp) is the crossing-over rate per generation between adjacent base pairs) in the presence of gene conversion. In the African sample, LD and polymorphism levels lead to similar estimates of effective population size, as expected under an equilibrium model. Conversely, in both non-African samples, LD levels suggest a smaller effective population size than that implied by polymorphism levels. This observation is paralleled by significant departures from an equilibrium model in the spectrum of allele frequencies of the non-African samples. Besides ruling out the possibility that non-African populations are at equilibrium, these results suggest different demographic history (temporal and spatial) of these groups. Interestingly, the African sample fits the expectations of an equilibrium model based on polymorphism and divergence levels and on frequency spectrum. For this sample, the estimated ratio of gene conversion to crossing-over rates is 7.3 for a mean tract length of 500 bp, suggesting that gene conversion may be more frequent than previously thought. These findings imply that disease-association studies will require a much denser map of polymorphic sites in African than in non-African populations.

Polymorphism and divergence in the beta-globin replication origin initiation region

DNA sequence polymorphism and divergence was examined in the vicinity of the human beta-globin gene cluster origin of replication initiation region (IR), a 1.3-kb genomic region located immediately 5' of the adult-expressed beta-globin gene. DNA sequence variation in the replication origin IR and 5 kb of flanking DNA was surveyed in samples drawn from two populations, one African (from the Gambia, West Africa) and the other European (from Oxford, England). In these samples, levels of nucleotide and length polymorphism in the IR were found to be more than two times as high as adjacent non-IR-associated regions (estimates of per-nucleotide heterozygosity were 0.30% and 0.12%, respectively). Most polymorphic positions identified in the origin IR fall within or just adjacent to a 52-bp alternating purine-pyrimidine ((RY)n) sequence repeat. Within- and between-populations divergence is highest in this portion of the IR, and interspecific divergence in the same region, determined by comparison with an orthologous sequence from the chimpanzee, is also pronounced. Higher levels of diversity in this subregion are not, however, primarily attributable to slippage-mediated repeat unit changes, as nucleotide substitution contributes disproportionately to allelic heterogeneity. An estimate of helical stability in the sequenced region suggests that the hypervariable (RY)n constitutes the major DNA unwinding element (DUE) of the replication origin IR, the location at which the DNA duplex first unwinds and new strand synthesis begins. These findings suggest that the beta-globin IR experiences a higher underlying rate of neutral mutation than do adjacent genomic regions and that enzyme fidelity associated with the initiation of DNA replication at this origin may be compromised. The significance of these findings for our understanding of eukaryotic replication origin biology is discussed.

Determination of ancestral alleles for human single-nucleotide polymorphisms using high-density oligonucleotide arrays

Here we report the application of high-density oligonucleotide array (DNA chip)-based analysis to determine the distant history of single nucleotide polymorphisms (SNPs) in current human populations. We analysed orthologues for 397 human SNP sites (identified in CEPH pedigrees from Amish, Venezuelan and Utah populations) from 23 common chimpanzee, 19 pygmy chimpanzee and 11 gorilla genomic DNA samples. From this data we determined 214 proposed ancestral alleles (the sequence found in the last common ancestor of humans and chimpanzees). In a diverse human population set, we found that SNP alleles with higher frequencies were more likely to be ancestral than less frequently occurring alleles. There were, however, exceptions. We also found three shared human/pygmy chimpanzee polymorphisms, all involving CpG dinucleotides, and two shared human/gorilla polymorphisms, one involving a CpG dinucleotide. We demonstrate that microarray-based assays allow rapid comparative sequence analysis of intra- and interspecies genetic variation.

X chromosome evidence for ancient human histories

Diverse African and non-African samples of the X-linked PDHA1 (pyruvate dehydrogenase E1 alpha subunit) locus revealed a fixed DNA sequence difference between the two sample groups. The age of onset of population subdivision appears to be about 200 thousand years ago. This predates the earliest modern human fossils, suggesting the transformation to modern humans occurred in a subdivided population. The base of the PDHA1 gene tree is relatively ancient, with an estimated age of 1.86 million years, a late Pliocene time associated with early species of Homo. PDHA1 revealed very low variation among non-Africans, but in other respects the data are consistent with reports from other X-linked and autosomal haplotype data sets. Like these other genes, but in conflict with microsatellite and mitochondrial data, PDHA1 does not show evidence of human population expansion.

Complete nucleotide sequencing of an HPV-1a variant and determination of extant errors in the prototype HPV-1a sequence

The complete nucleotide sequence of an HPV-1a variant has been determined. The variant has over 99% nucleotide identity with the prototype HPV-1a sequence, with the majority of mutations occurring in the long control region. Additionally, remaining errors in the prototype HPV-1a sequence are reported.

Biological basis of germline mutation: comparisons of spontaneous germline mutation rates among drosophila, mouse, and human

Spontaneous mutation rates per generation are similar among the three species considered here--Drosophila, mouse, and human--and are not related to time, as is often assumed. Spontaneous germline mutation rates per generation averaged among loci are less variable among species than they are among loci and tests and between gender. Mutation rates are highly variable over time in diverse lineages. Recent estimates of the number of germ cell divisions per generation are: for humans, 401 (30-year generation) in males and 31 in females; for mice, 62 (9-month generation) in males and 25 in females; and for Drosophila melanogaster, 35.5 (18-day generation) in males and 36.5 (25-day generation) in females. The relationships between germ cell division estimates of the two sexes in the three species closely reflect those between mutation rates in the sexes, although mutation rates per cell division vary among species. Whereas the overall rate per generation is constant among species, this consistency must be achieved by diverse mechanisms. Modifiers of mutation rates, on which selection might act, include germline characteristics that contribute disproportionately to the total mutation rates. The germline mutation rates between the sexes within a species are largely influenced by germ cell divisions per generation. Also, a large portion of the total mutations occur during the interval between the beginning of meiosis and differentiation of the soma from the germline. Significant genetic events contributing to mutations during this time may include meiosis, lack of DNA repair in sperm cells, methylation of CpG dinucleotides in mammalian sperm and early embryo, gonomeric fertilization, and rapid cleavage divisions.

Molecular and population genetic analysis of allelic sequence diversity at the human beta-globin locus

Allelic sequence polymorphism at the beta-globin locus was investigated in a group of 36 Melanesians. A 3-kilobase fragment containing the gene and its flanking regions was sequenced in 60 normal (beta A) and 12 thalassemic (intron 1, position 5, G-->C) chromosomes. Haplotype relationships between linked polymorphisms were derived by allele-specific PCR amplification and sequencing. Seventeen nucleotide polymorphisms and 2 length variants were identified, and these sites segregated as 17 sequence haplotypes in the normal chromosomes. This haplotype diversity is higher than that expected on the basis of the nucleotide polymorphism observed and is probably due to recombination and gene conversion. Nucleotide diversity at synonymous sites in the sample is 0.14%, suggesting an average age of sequence divergence of approximately 450,000 years, consistent with that expected for a neutrally evolving human nuclear locus.

Possible implications of CpG avoidance in the flatworm Schistosoma mansoni

We report the analysis of the biases of CpG, TpG, and CpA of all the DNA sequences data from the Trematode Schistosoma mansoni. Our results show CpG avoidance whereas TpG and CpA frequencies are over the expected values. These characteristics are similar to the biases displayed by methylated genomes, but in platyhelminths 5mC has not been detected by biochemical methods. The possible implications of this CpG shortage are discussed.

Evolution of repeated sequences in non-coding regions of the genome

Repeated sequences are found ubiquitously in the eukaryotic genome. Population genetic studies on the evolution of such repeated sequences are reviewed while paying special attention to those sequences found in the non-coding regions of the genome. Specifically, the evolution of dispersed repeated sequences by the transposition as well as the evolution of short tandemly repeated sequences due to either replication slippage or unequal sister chromatid exchange are considered. The approach of combining both model and data analyses which has been successfully employed in the development of the neutral theory is also considered to be useful in better understanding the evolution and biological meaning of these sequences.

The beta globin gene cluster of the prosimian primate Galago crassicaudatus: nucleotide sequence determination of the 41-kb cluster and comparative sequence analyses

The nucleotide sequence of the beta globin gene cluster of the prosimian Galago crassicaudatus has been determined. A total sequence spanning 41,101 bp contains and links together previously published sequences of the five galago beta-like globin genes (5'-epsilon-gamma-psi eta-delta-beta-3'). A computer-aided search for middle interspersed repetitive sequences identified 10 LINE (L1) elements, including a 5' truncated repeat that is orthologous to the full-length L1 element found in the human epsilon-gamma intergenic region. SINE elements that were identified included one Alu type I repeat, four Alu type II repeats, and two methionine tRNA-derived Monomer (type III) elements. Alu type II and Monomer sequences are unique to the galago genome. Structural analyses of the cluster sequence reveals that it is relatively A+T rich (about 62%) and regions with high G+C content are associated primarily with globin coding regions. Comparative analyses with the beta globin cluster sequences of human, rabbit, and mouse reveal extensive sequence homologies in their genic regions, but only human, galago, and rabbit sequences share extensive intergenic sequence homologies. Divergence analyses of aligned intergenic and flanking sequences from orthologous human, galago, and rabbit sequences show a gradation in the rate of nucleotide sequence evolution along the cluster where sequences 5' of the epsilon globin gene region show the least sequence divergence and sequences just 5' of the beta globin gene region show the greatest sequence divergence.

Evolution of the primate beta-globin gene region: nucleotide sequence of the delta-beta-globin intergenic region of gorilla and phylogenetic relationships between African apes and man

A 6.0-kb DNA fragment from Gorilla gorilla including the 5' part of the beta-globin gene and about 4.5 kb of its upstream flanking region was cloned and sequenced. The sequence was compared to the human, chimpanzee, and macaque delta-beta intergenic region. This analysis reveals four tandemly repeated sequences (RS), at the same location in the four species, showing a variable number of repeats generating both intraspecific (polymorphism) and interspecific variability. These tandem arrays delimit five regions of unique sequence called IG for intergenic. The divergence for these IG sequences is 1.85 +/- 0.22% between human and gorilla, which is not significantly different from the value estimated in the same region between chimpanzee and human (1.62 +/- 0.21%). The CpG and TpA dinucleotides are avoided. CpGs evolve faster than other sequence sites but do not confuse phylogenetic inferences by producing parallel mutations in different lineages. About 75% of CpG doublets have become TpG or CpA since the common ancestor, in agreement with the methylation/deamination pattern. Comparison of this intergenic region gives information on branching order within Hominoidea. Parsimony and distance-based methods when applied to the delta-beta intergenic region provide evidence (although not statistically significant) that human and chimpanzee are more closely related to each other than to gorilla. CpG sites are indeed rich in information by carrying substitutions along the short internal branch. Combining these results with those on the psi eta-delta intergenic region, shows in a statistically significant way that chimpanzee is the closest relative of human.

Nonrandom CpG mutations affect the synonymous codon usage of moderately GC-rich single copy actin genes

In species where actin genes exist as single copies, analysis of their synonymous codon usage and of the substitutions occurring between the genes of closely related species shows that there is a positive selection for codons that do not have highly mutable CpG dinucleotides in codon positions 2 and 3 when the GC content of these genes is less than 57%.

Ionizing radiation and genetic risks. I. Epidemiological, population genetic, biochemical and molecular aspects of Mendelian diseases

This paper reviews the currently available information on naturally occurring Mendelian diseases in man; it is aimed at providing a background and framework for discussion of experimental data on radiation-induced mutations (papers II and III) and for the estimation of the risk of Mendelian disease in human populations exposed to ionizing radiation (paper IV). Current consensus estimates indicate that a total of about 125 per 10(4) livebirths are directly affected by one or another naturally occurring Mendelian disease (autosomal dominants, 95/10(4); X-linked ones, 5/10(4); and autosomal recessives, 25/10(4). These estimates are conservative and take into account conditions which are very rare and for which prevalence estimates are unavailable. Most, although not all, of the recognized "common" dominants have onset in adult ages while most sex-linked and autosomal recessives have onset at birth or in childhood. Autosomal dominant and X-linked diseases (i.e., the responsible mutant alleles) presumed to be maintained in the population due to a balance between mutation and selection are the ones which may be expected to increase in frequency as a result of radiation exposures. Viewed from this standpoint, the above assumption seems safe only for a small proportion of such diseases; for the remainder, there is no easy way to discriminate between different mechanisms that may be responsible or to rigorously exclude some in favor of some others. Mutations in genes that code for enzymic proteins are more often recessive in contrast to those that code for non-enzymic proteins, which are more often dominant. At the molecular level, with recessives, a wide variety of changes is possible and these include specific types of point mutations, small and large intragenic deletions, multilocus deletions and rearrangements. In the case of dominants, however, the kinds of recoverable point mutations and deletion-type changes are less extensive because of functional constraints. The mutational potential of genes varies, depending on the gene, its size, sequence content and arrangement, location and its normal functions, and can be grouped into three groups: those in which only point mutations have been found to occur, those in which only deletions or other gross changes have been recovered and those in which both kinds of changes are known. Molecular data are available for about 75 Mendelian conditions and these suggest that in approximately 50% of them, the changes categorized to date are point mutations and in the remainder, intragenic deletions or other gross changes; there does not seem to be any fundamental difference between dominants and recessives with respect to the underlying molecular defect.(ABSTRACT TRUNCATED AT 400 WORDS)

Cloning and nucleotide sequence of the chimpanzee c-myc gene

A DNA fragment covering the chimpanzee c-myc locus was cloned from the DNA of peripheral blood lymphocytes, sequenced, and compared to its human c-myc counterpart. The two nucleotide sequences were found to be highly homologous (99%). The divergence rate between the two species was 0.4% in exons and 1.7% in introns. The different TATA-boxes described in the human myc gene were also identified in the chimpanzee sequence and an open reading frame (ORF) was observed which overlaps the chimpanzee c-myc first exon. This latter ORF contained three silent mutations with regard to the human region, whereas the chimpanzee Myc oncoprotein coded by exons 2 and 3 differed by two amino acids from the human one.

The mutational spectrum of single base-pair substitutions causing human genetic disease: patterns and predictions

Reports of single base-pair substitutions that cause human genetic disease and that have been located and characterized in an unbiased fashion were collated; 32% of point mutations were CG----TG or CG----CA transitions consistent with a chemical model of mutation via methylation-mediated deamination. This represents a 12-fold higher frequency than that predicted from random expectation, confirming that CG dinucleotides are indeed hotspots of mutation causing human genetic disease. However, since CG also appears hypermutable irrespective of methylation-mediated deamination, a second mechanism may also be involved in generating CG mutations. The spectrum of point mutations occurring outwith CG dinucleotides is also non-random, at both the mono- and dinucleotide, levels. An intrinsic bias in clinical detection was excluded since frequencies of specific amino acid substitutions did not correlate with the 'chemical difference' between the amino acids exchanged. Instead, a strong correlation was observed with the mutational spectrum predicted from the experimentally measured mispairing frequencies of vertebrate DNA polymerases alpha and beta in vitro. This correlation appears to be independent of any difference in the efficiency of enzymatic proofreading/mismatch-repair mechanisms but is consistent with a physical model of mutation through nucleotide misincorporation as a result of transient misalignment of bases at the replication fork. This model is further supported by an observed correlation between dinucleotide mutability and stability, possibly because transient misalignment must be stabilized long enough for misincorporation to occur. Since point mutations in human genes causing genetic disease neither arise by random error nor are independent of their local sequence environment, predictive models may be considered. We present a computer model (MUTPRED) based upon empirical data; it is designed to predict the location of point mutations within gene coding regions causing human genetic disease. The mutational spectrum predicted for the human factor IX gene was shown to resemble closely the observed spectrum of point mutations causing haemophilia B. Further, the model was able to predict successfully the rank order of disease prevalence and/or mutation rates associated with various human autosomal dominant and sex-linked recessive conditions. Although still imperfect, this model nevertheless represents an initial attempt to relate the variable prevalence of human genetic disease to the mutability inherent in the nucleotide sequences of the underlying genes.

A directly repeated sequence in the beta-globin promoter regulates transcription in murine erythroleukemia cells

We have identified a previously undetected cis-acting element in the mouse beta-major globin promoter region that is necessary for maximal transcription levels of the gene in the inducible preerythroid murine erythroleukemia (MEL) cell line. This element, termed the beta-globin direct-repeat element (beta DRE), consists of a directly repeated 10-base-pair sequence, 5'-AGGGCAG(G)AGC-3', that lies just upstream from the TATA box of the promoter. The beta DRE motif is highly conserved in all adult mammalian beta-globin promoter sequences known. Mutation of either single repeat alone caused less than a twofold decrease in transcript levels. However, simultaneous mutation of both repeated regions resulted in a ninefold decrease in accumulated transcripts when the gene was transiently transfected into MEL cells. Attachment of the beta DRE to a heterologous promoter had little effect on levels of accumulated transcripts initiated from the promoter in undifferentiated MEL cells but resulted in a threefold increase in transcript levels in induced (differentiated) MEL cells. Similarly, a comparison of the relative effects of mutations in the beta DRE in uninduced and induced MEL cells indicated that the element was more active in induced cells. The increase in beta DRE activity upon MEL cell differentiation and the more pronounced effects of mutations in both repeats of the beta DRE have implications for the mechanism of action of the element in regulating beta-globin transcription and for mutational studies of other repetitive or redundant transcription elements.

A directly repeated sequence in the beta-globin promoter regulates transcription in murine erythroleukemia cells

We have identified a previously undetected cis-acting element in the mouse beta-major globin promoter region that is necessary for maximal transcription levels of the gene in the inducible preerythroid murine erythroleukemia (MEL) cell line. This element, termed the beta-globin direct-repeat element (beta DRE), consists of a directly repeated 10-base-pair sequence, 5'-AGGGCAG(G)AGC-3', that lies just upstream from the TATA box of the promoter. The beta DRE motif is highly conserved in all adult mammalian beta-globin promoter sequences known. Mutation of either single repeat alone caused less than a twofold decrease in transcript levels. However, simultaneous mutation of both repeated regions resulted in a ninefold decrease in accumulated transcripts when the gene was transiently transfected into MEL cells. Attachment of the beta DRE to a heterologous promoter had little effect on levels of accumulated transcripts initiated from the promoter in undifferentiated MEL cells but resulted in a threefold increase in transcript levels in induced (differentiated) MEL cells. Similarly, a comparison of the relative effects of mutations in the beta DRE in uninduced and induced MEL cells indicated that the element was more active in induced cells. The increase in beta DRE activity upon MEL cell differentiation and the more pronounced effects of mutations in both repeats of the beta DRE have implications for the mechanism of action of the element in regulating beta-globin transcription and for mutational studies of other repetitive or redundant transcription elements.

Cytosine methylation and the fate of CpG dinucleotides in vertebrate genomes

The dinucleotide CpG is a "hotspot" for mutation in the human genome as a result of (1) the modification of the 5' cytosine by cellular DNA methyltransferases and (2) the consequent high frequency of spontaneous deamination of 5-methyl cytosine (5mC) to thymidine. DNA methylation thus contributes significantly, albeit indirectly, to the incidence of human genetic disease. We have attempted to estimate for the first time the in vivo rate of deamination of 5mC from the measured rate of 5mC deamination in vitro and the known error frequency of the cellular G/T mismatch-repair system. The accuracy and utility of this estimate (md) was then assessed by comparison with clinical data, and an improved estimate of md (1.66 X 10(-16) s-1) was derived. Comparison of the CpG mutation rates exhibited by globin gene and pseudogene sequences from human, chimpanzee and macaque provided further estimates of md, all of which were consistent with the first. Use of this value in a mathematical model then permitted the estimation of the length of time required to produce the level of "CpG suppression" currently found in the "bulk DNA" of vertebrate genomes. This time span, approximately 450 million years, corresponds closely to the estimated time since the emergence and adaptive radiation of the vertebrates and thus coincides with the probable advent of heavily methylated genomes. An accurate estimate of the 5mC deamination rate is important not only for clinical medicine but also for studies of gene evolution. Our data suggest both that patterns of vertebrate gene methylation may be comparatively stable over relatively long periods of evolutionary time, and that the rate of CpG deamination can, under certain limited conditions, serve as a "molecular clock".

Nucleotide sequence of the rat gamma-crystallin gene region and comparison with an orthologous human region

The sequences of a 51-kb region containing the cluster of five rat gamma-crystallin-coding genes (CRYG) and of a 7-kb region surrounding the sixth rat CRYG gene were determined. Approximately 78% of the total sequence represents intergenic DNA. We also sequenced 22 kb of DNA from the human CRYG gene cluster. All CRYG genes are associated with CpG-rich regions. The sequence similarity between the human and rat gene regions drops sharply (to 65%) in intronic and 3'-flanking regions but decreases only gradually in the 5'-flanking region. Highly conserved regions (greater than 80%) are found as far upstream as 1.5 kb. Overall intergenic distances are conserved. The human region contains much more repetitive DNA (24% vs. 10%) but less simple-sequence (sps) DNA (0.7% vs. 4%) than the rat region. Almost all repeats and spsDNA elements are located in the intergenic region. The location of repetitive and spsDNA differs between the orthologous regions and these elements were probably inserted after the evolutionary separation of rat and man. The Alu repeats in man and the B3 repeats in the rat are close copies of their respective consensus sequences and bordered by virtually perfect repeats. In contrast, the B1 and B2 repeats in the rat have diverged considerably from the consensus sequence and the surrounding direct repeats are usually imperfect. Thus the dispersion of the B1 and B2 repeats in the rat probably preceded that of the B3 repeats. Within the rat genomic region the spacing of Z-DNA elements is surprisingly regular, they are located about 12 kb apart. A search for putative matrix-associated regions suggests that the rat CRYG gene cluster is organized into two chromosomal domains.

Structural analysis of the 5' flanking region of the beta-globin gene in African sickle cell anemia patients: further evidence for three origins of the sickle cell mutation in Africa

Haplotype analysis of the beta-globin gene cluster shows two regions of DNA characterized by nonrandom association of restriction site polymorphisms. These regions are separated by a variable segment containing the repeated sequences (ATTTT)n and (AT)xTy, which might be involved in recombinational events. Studies of haplotypes linked to the sickle cell gene in Africa provide strong argument for three origins of the mutation: Benin, Senegal, and the Central African Republic. Nevertheless, the haplotype determination does not give any information about the variable segment and does not totally exclude the possibility of recombination leading to different haplotypes linked to the mutation. The structure of the variable segment in the three African populations was studied by S1 nuclease mapping of genomic DNA, which allows a comparison of several samples. A 1080-base-pair DNA segment was sequenced for one sample from each population. S1 nuclease mapping confirmed the homogeneity of each population with regard to both (ATTTT)n and (AT)xTy repeats. We found three additional structures for (AT)xTy correlating with the geographic origin of the patients. Ten other nucleotide positions, 5' and 3' to the (AT)xTy copies, were found to be variable when compared to homologous sequences from human and monkey DNAs. These results allow us to propose an evolutionary scheme for the polymorphisms in the 5' flanking region of the beta-globin gene. The results strongly support the hypothesis of three origins for the sickle mutation in Africa.

The CpG dinucleotide and human genetic disease

Reports of single base-pair mutations within gene coding regions causing human genetic disease were collated. Thirty-five per cent of mutations were found to have occurred within CpG dinucleotides. Over 90% of these mutations were C----T or G----A transitions, which thus occur within coding regions at a frequency 42-fold higher than that predicted from random mutations. These findings are consistent with methylation-induced deamination of 5-methyl cytosine and suggest that methylation of DNA within coding regions may contribute significantly to the incidence of human genetic disease.

Phylogenetic relations of humans and African apes from DNA sequences in the psi eta-globin region

Sequences from the upstream and downstream flanking DNA regions of the psi eta-globin locus in Pan troglodytes (common chimpanzee), Gorilla gorilla (gorilla), and Pongo pygmaeus (orangutan, the closest living relative to Homo, Pan, and Gorilla) provided further data for evaluating the phylogenetic relations of humans and African apes. These newly sequenced orthologs [an additional 4.9 kilobase pairs (kbp) for each species] were combined with published psi eta-gene sequences and then compared to the same orthologous stretch (a continuous 7.1-kbp region) available for humans. Phylogenetic analysis of these nucleotide sequences by the parsimony method indicated (i) that human and chimpanzee are more closely related to each other than either is to gorilla and (ii) that the slowdown in the rate of sequence evolution evident in higher primates is especially pronounced in humans. These results indicate that features (for example, knuckle-walking) unique to African apes (but not to humans) are primitive and that even local molecular clocks should be applied with caution.

Recent insertion of an Alu sequence in the beta-globin gene cluster of the gorilla

We present the nucleotide sequence of a new Alu family member that lies between the delta- and beta-globin genes in gorilla DNA. The sequence exhibits 91% similarity with a consensus sequence of the Alu family. It is flanked by a perfect repetition of a 16-nucleotide target sequence and terminates with 24 adenylic residues. As this sequence is absent at this locus in other primate DNAs, its insertion occurred less than 8 million years ago, thus supporting the idea that Alu sequences are still mobile elements in the hominoid genome.

The distribution of the dinucleotide CpG and cytosine methylation in the vitellogenin gene family

Sequence data from regions of five vertebrate vitellogenin genes were used to examine the frequency, distribution, and mutability of the dinucleotide CpG, the preferred modification site for eukaryotic DNA methyltransferases. The observed level of the CpG dinucleotide in all five genes was markedly lower than that expected from the known mononucleotide frequencies. CpG suppression was greater in introns than in exons. CpG-containing codons were found to be avoided in the vitellogenin genes, but not completely despite the redundancy of the genetic code. Frequency and distribution patterns of this dinucleotide varied dramatically among these otherwise closely related genes. Dense clusters of CpG dinucleotides tended to appear in regions of either functional or structural interest (e.g., in the transposon-like Vi-element of Xenopus) and these clusters contained 5-methylcytosine (5 mC). 5 mC is known to undergo deamination to form thymidine, but the extent to which this transition occurs in the heavily methylated genomes of vertebrates and its contribution to CpG suppression are still unclear. Sequence comparison of the methylated vitellogenin gene regions identified C----T and G----A substitutions that were found to occur at relatively high frequencies. The predicted products of CpG deamination, TpG and CpA, were elevated. These findings are consistent with the view that CpG distribution and methylation are interdependent and that deamination of 5 mC plays an important role in promoting evolutionary change at the nucleotide sequence level.

Nucleotide sequence of the beta-globin genes in gorilla and macaque: the origin of nucleotide polymorphisms in human

Part of the beta-globin genes of Macaca cynomolgus and Gorilla gorilla has been cloned and sequenced. Ten putatively neutral nucleotide polymorphisms have been described at the beta-globin locus in humans. They are associated in seven combinations, which define seven different haplotypes of the beta-globin gene: four major frameworks--1, 2, 3, and 3--and three minor frameworks, which we term KI1, KA1, and OR1. The nucleotide sequences of these frameworks are compared with those of homologous sequences in chimpanzee, colobus, macaque, and gorilla. This comparison provides strong evidence that framework 2 was the earliest framework in the human lineage. From framework 2, a rooted parsimonious tree for the six other frameworks is constructed. This phylogenetic tree is discussed in terms of the evolution of nucleotide polymorphisms as well as in terms of genetic affinities between human populations. For each position at which there is base difference in comparing human, gorilla, and chimpanzee beta-globin genes, the phyletic lineage where the corresponding substitution occurred has been identified using the maximum parsimony procedure. The data provide evidence that polymorphisms may represent a significant component of differences between closely related species. If so, nucleotide polymorphisms may strongly bias estimates of small evolutionary distances.

Nucleotide sequence of the delta-beta-globin intergenic segment in the macaque: structure and evolutionary rates in higher primates

A 5600-base-pair (bp) fragment including the beta-globin gene and about 4000 bp of its 5' flanking sequence was cloned from the DNA of Macaca cynomolgus (an Old World monkey), and the 5' flanking region was sequenced. Comparison with human, chimpanzee, mouse, rabbit, and Xenopus orthologous sequences reveals a tandemly repeated sequence called RS4 at the same position (about 500 bp 5' from the transcription start of the adult beta-globin gene) in all six species. We suggest that a tandemly repeated sequence has been maintained by functional constraints since the divergence between amphibians and reptiles. Excluding tandemly repeated sequences as well as about 400 nucleotides upstream from the cap site, the average base substitution frequencies among human, chimpanzee, and macaque intergenic sequences were calculated. They appear to be strongly correlated with the delta T50 values measured between the corresponding nuclear DNAs. They are also similar to base substitution frequencies calculated by Chang and Slightom (1984) at the pseudo-eta-globin locus. Thus, exclusion of sequences involved in specific modes of variation might allow the use of intergenic sequences for the accurate calculation of genetic distances. Using a time scale based on the dating of the Atlantic split, we estimate the base substitution rate of primate noncoding DNA to be 1.0 X 10(-9) substitution/site/year.

(A-T)n tracts embedded in random sequence DNA--formation of a structure which is chemically reactive and torsionally deformable

Alternating d(A-T)n sequences which are contiguous with DNA of effectively random sequence have an abnormal conformation in linear DNA molecules. These regions are strongly reactive towards chemical modification by osmium tetroxide, and are preferentially cleaved by micrococcal nuclease. Both the chemical modification and the enzymic cutting occur uniformly through the alternating tract, and there is no evidence for enzyme or chemical sensitivity in the interfaces between the tract and DNA of normal conformation. These reactivities have a requirement for an alternating sequence. In addition to chemical reactivity, alternating (A-T)n sequences exhibit anomalously small twist changes on cruciform formation, suggesting that the pre-extruded DNA is underwound. We propose that the alternating sequences adopt an altered conformation which is subject to easy torsional deformation.

DNA cytosine methylation and heat-induced deamination

The heat-induced conversion of 5-methylcytosine (m5C) residues to thymine residues and of cytosine to uracil residues in single-stranded DNA was studied. The calculated rates for deamination at 37 degrees C and pH 7.4 were approximately 9.5 X 10(-10) and 2.1 X 10(-10) sec-1, respectively. N4-Methyldeoxycytidine, which is in the DNA of certain thermophilic bacteria, was more heat-resistant than was deoxycytidine and much more than was 5-methyldeoxycytidine. Thermophilic bacteria which contain N4-methylcytosine rather than m5C in their genomes may thereby largely avoid heat-induced mutation due to deamination, which is incurred by the many organisms that contain m5C in their DNA.

Rates of nuclear DNA evolution in pheasant-like birds: evidence from restriction maps

To examine the tempo of genomic evolution in birds, we mapped 161 restriction sites in the nuclear DNA of seven species of birds belonging to the pheasant superfamily Phasianoidea. The three regions mapped lie on different chromosomes and bear eight genes, coding for lysozyme c, three "alpha-like" globins, and four "beta-like" globins. Together, the three regions span about 56 kilobases, most of which is presumably noncoding. The maps differed from one another at a minimum of 77 sites and by 9 length mutations. The extent of sequence divergence due to base substitutions was inferred to be similar for all three regions, even though the three coding regions differ by 5-fold from one another in mean rate of evolution at the amino acid level. A tree relating the maps differs in branching order from that implied by the traditional classification of phasianoid birds and is supported by published protein comparisons. Five of the nodes in the tree were associated with fossil evidence and historical biogeographic information, allowing us to estimate the mean rate of DNA divergence to be 0.34-0.40% per million years. This rate is similar to that estimated for the globin gene regions of higher primates, which validates the concept of an evolutionary clock at the DNA level. Our fossil-based calibration of DNA evolution differs by a factor of almost 2 from that proposed by others on the basis of biogeography. In consequence, published estimates of divergence times for birds and primates that are based on a biogeographically calibrated DNA clock may be too long.

Characterization of a set of X-linked sequences and of a panel of somatic cell hybrids useful for the regional mapping of the human X chromosome

We have characterized 19 DNA fragments originating from the human X chromosome. Most of them have been isolated from an X chromosome genomic library (Davies et al. 1981) using a systematic screening procedure. These DNA probes have been used to search for restriction fragment length polymorphisms (RFLP). The frequency of restriction polymorphisms (1 per 350 bp analysed) was lower than expected from data obtained with autosomal fragments. The various probes have been mapped within 12 subchromosomal regions using a panel of human-rodent hybrid cell lines. The validity of the panel was established by hybridization experiments performed with 27 X-specific DNA probes, which yielded information on the relative position of translocation breakpoints on the X chromosome. The DNAs from the various hybrid lines are blotted onto a reusable support which allows one to quickly map any new X-specific DNA fragment. The probes already isolated should be of use to map unbalanced X chromosome aberrations or to characterize new somatic cell hybrid lines. The probes which detect RFLPs define new genetic markers which will help to construct a detailed linkage map of the human X chromosome, and might also serve for the diagnosis of carriers or prenatal diagnosis.

Nucleotide sequence, evolution, and expression of the fetal globin gene of the spider monkey Ateles geoffroyi

The single gamma-globin gene of the New World spider monkey Ateles geoffroyi is similar to other primate genes of the beta-globin gene cluster ("beta-like" globin genes). The number of nonsynonymous nucleotide substitutions between the coding regions of Ateles and other primate gamma-globin genes suggests that the Platyrrhine and Catarrhine evolutionary lines diverged approximately equal to 40 million years ago, an estimate reasonably consistent with the fossil record. However, the number of synonymous coding region and noncoding base differences is much smaller than predicted by various molecular "clocks." This suggests that the rate of synonymous coding and noncoding base substitution has not been constant per absolute time in primate lineages. Expression of the cloned Ateles gamma-globin gene in cultured monkey cells showed that the sequence AAUAAA near the mRNA 3' terminus is not sufficient to define the site of transcript polyadenylylation.