Tracing the origin of HLA-DRB1 alleles by microsatellite polymorphism

Complete nucleotide sequence characterization of DRB5 alleles reveals a homogeneous allele group that is distinct from other DRB genes

Next Generation Sequencing allows for testing and typing of entire genes of the HLA region. A better and comprehensive sequence assessment can be achieved by the inclusion of full gene sequences of all the common alleles at a given locus. The common alleles of DRB5 are under-characterized with the full exon-intron sequence of two alleles available. In the present study the DRB5 genes from 18 subjects alleles were cloned and sequenced; haplotype analysis showed that 17 of them had a single copy of DRB5 and one consanguineous subject was homozygous at all HLA loci. Methodological approaches including robust and efficient long-range PCR amplification, molecular cloning, nucleotide sequencing and de novo sequence assembly were combined to characterize DRB5 alleles. DRB5 sequences covering from 5'UTR to the end of intron 5 were obtained for DRB5*01:01, 01:02 and 02:02; partial coverage including a segment spanning exon 2 to exon 6 was obtained for DRB5*01:03, 01:08N and 02:03. Phylogenetic analysis of the generated sequences showed that the DRB5 alleles group together and have distinctive differences with other DRB loci. Novel intron variants of DRB5*01:01:01, 01:02 and 02:02 were identified. The newly characterized DRB5 intron variants of each DRB5 allele were found in subjects harboring distinct associations with alleles of DRB1, B and/or ethnicity. The new information provided by this study provides reference sequences for HLA typing methodologies. Extending sequence coverage may lead to identify the disease susceptibility factors of DRB5 containing haplotypes while the unexpected intron variations may shed light on understanding of the evolution of the DRB region.

Homology modelling of frequent HLA class-II alleles: A perspective to improve prediction of HLA binding peptide and understand the HLA associated disease susceptibility

Human leukocyte antigen (HLA) plays significant role via the regulation of immune system and contribute in the progression and protection of many diseases. HLA molecules bind and present peptides to T- cell receptors which generate the immune response. HLA peptide interaction and molecular function of HLA molecule is the key to predict peptide binding and understanding its role in different diseases. The availability of accurate three dimensional (3D) structures is the initial step towards this direction. In the present work, homology modelling of important and frequent HLA-DRB1 alleles (07:01, 11:01 and 09:01) was done and acceptable models were generated. These modelled alleles were further refined and cross validated by using several methods including Ramachandran plot, Z-score, ERRAT analysis and root mean square deviation (RMSD) calculations. It is known that numbers of allelic variants are related to the susceptibility or protection of various infectious diseases. Difference in amino acid sequences and structures of alleles were also studied to understand the association of HLA with disease susceptibility and protection. Susceptible alleles showed more amino acid variations than protective alleles in three selected diseases caused by different pathogens. Amino acid variations at binding site were found to be more than other part of alleles. RMSD values were also higher at variable positions within binding site. Higher RMSD values indicate that mutations occurring at peptide binding site alter protein structure more than rest of the protein. Hence, these findings and modelled structures can be used to design HLA-DRB1 binding peptides to overcome low prediction accuracy of HLA class II binding peptides. Furthermore, it may help to understand the allele specific molecular mechanisms involved in susceptibility/resistance against pathogenic diseases.

Microsatellite interruptions stabilize primate genomes and exist as population-specific single nucleotide polymorphisms within individual human genomes

Interruptions of microsatellite sequences impact genome evolution and can alter disease manifestation. However, human polymorphism levels at interrupted microsatellites (iMSs) are not known at a genome-wide scale, and the pathways for gaining interruptions are poorly understood. Using the 1000 Genomes Phase-1 variant call set, we interrogated mono-, di-, tri-, and tetranucleotide repeats up to 10 units in length. We detected ∼26,000–40,000 iMSs within each of four human population groups (African, European, East Asian, and American). We identified population-specific iMSs within exonic regions, and discovered that known disease-associated iMSs contain alleles present at differing frequencies among the populations. By analyzing longer microsatellites in primate genomes, we demonstrate that single interruptions result in a genome-wide average two- to six-fold reduction in microsatellite mutability, as compared with perfect microsatellites. Centrally located interruptions lowered mutability dramatically, by two to three orders of magnitude. Using a biochemical approach, we tested directly whether the mutability of a specific iMS is lower because of decreased DNA polymerase strand slippage errors. Modeling the adenomatous polyposis coli tumor suppressor gene sequence, we observed that a single base substitution interruption reduced strand slippage error rates five- to 50-fold, relative to a perfect repeat, during synthesis by DNA polymerases α, β, or η. Computationally, we demonstrate that iMSs arise primarily by base substitution mutations within individual human genomes. Our biochemical survey of human DNA polymerase α, β, δ, κ, and η error rates within certain microsatellites suggests that interruptions are created most frequently by low fidelity polymerases. Our combined computational and biochemical results demonstrate that iMSs are abundant in human genomes and are sources of population-specific genetic variation that may affect genome stability. The genome-wide identification of iMSs in human populations presented here has important implications for current models describing the impact of microsatellite polymorphisms on gene expression.

Microsatellites are short tandem repeat DNA sequences located throughout the human genome that display a high degree of inter-individual variation. This characteristic makes microsatellites an attractive tool for population genetics and forensics research. Some microsatellites affect gene expression, and mutations within such microsatellites can cause disease. Interruption mutations disrupt the perfect repeated array and are frequently associated with altered disease risk, but they have not been thoroughly studied in human genomes. We identified interrupted mono-, di-, tri- and tetranucleotide MSs (iMS) within individual genomes from African, European, Asian and American population groups. We show that many iMSs, including some within disease-associated genes, are unique to a single population group. By measuring the conservation of microsatellites between human and chimpanzee genomes, we demonstrate that interruptions decrease the probability of microsatellite mutations throughout the genome. We demonstrate that iMSs arise in the human genome by single base changes within the DNA, and provide biochemical data suggesting that these stabilizing changes may be created by error-prone DNA polymerases. Our genome-wide study supports the model in which iMSs act to stabilize individual genomes, and suggests that population-specific differences in microsatellite architecture may be an avenue by which genetic ancestry impacts individual disease risk.

Characterising a microsatellite for DRB typing in Aotus vociferans and Aotus nancymaae (Platyrrhini)

Non-human primates belonging to the Aotus genus have been shown to be excellent experimental models for evaluating drugs and vaccine candidates against malaria and other human diseases. The immune system of this animal model must be characterised to assess whether the results obtained here can be extrapolated to humans. Class I and II major histocompatibility complex (MHC) proteins are amongst the most important molecules involved in response to pathogens; in spite of this, the techniques available for genotyping these molecules are usually expensive and/or time-consuming. Previous studies have reported MHC-DRB class II gene typing by microsatellite in Old World primates and humans, showing that such technique provides a fast, reliable and effective alternative to the commonly used ones. Based on this information, a microsatellite present in MHC-DRB intron 2 and its evolutionary patterns were identified in two Aotus species (A. vociferans and A. nancymaae), as well as its potential for genotyping class II MHC-DRB in these primates.

Analytical methods for disease association studies with immunogenetic data

Disease association studies involving highly polymorphic immunogenetic data may involve analyses at one or many units of analysis, including amino acid, allele, genotype and haplotype levels, as well as consideration of gene-gene or gene-environment interactions. The selection of the appropriate statistical tests is critical and will be dependent on the nature of the dataset (e.g., case-control vs. family data) as well as the specific research hypotheses being tested. This paper describes the various study and analysis categories used for such analyses, including the advantages and limitations of such techniques.

The chimpanzee Mhc-DRB region revisited: gene content, polymorphism, pseudogenes, and transcripts

In humans, great apes, and different monkey species, the major histocompatibility complex (MHC) class II DRB region is known to display considerable copy number variation. The microsatellite D6S2878 has been shown to be a valuable marker for haplotyping the DR region in humans and macaque species. The present report illustrates that chimpanzee haplotypes also can be discriminated with this marker. The analyses resulted in the description of nine different region configurations, of which seven are present within the West African chimpanzee population studied. The region configurations vary in gene content from two up to five DRB genes. Subsequent cDNA sequencing increased the number of known full-length Patr-DRB sequences from 3 to 32, and shows that one to three Patr-DRB genes per haplotype apparently produce functional transcripts. This is more or less comparable to humans and rhesus macaques. Moreover, microsatellite analysis in concert with full-length DRB gene sequencing showed that the Patr-DRB*W9 and -DRB3*01/02 lineages most likely arose from a common ancestral lineage: hence, the Patr-DRB*W9 lineage was renamed to Patr-DRB3*07. Overall, the data demonstrate that the D6S2878 microsatellite marker allows fast and accurate haplotyping of the Patr-DRB region. In addition, the limited amount of allelic variation observed at the various Patr-DRB genes is in agreement with the fact that chimpanzees experienced a selective sweep that may have been caused by an ancient retroviral infection.

High resolution definition of HLA-DRB haplotypes by a simplified microsatellite typing technique

In humans, the region configurations DR1, DR8, DR51, DR52 and DR53 are known to display copy number as well as allelic variation, rendering high resolution typing of HLA-DRB haplotypes cumbersome. Advantage was taken of microsatellite D6S2878, present in all DRB genes/pseudogenes with an intact exon 2-intron 2 segment. This DRB-STR is highly polymorphic in composition and length. Recently, it was proven that all exon 2 sequences could be linked to a certain DRB-STR that segregates with the respective DRB allele. Because haplotypes show differential copy numbers and compositions of exon 2-positive DRB genes/pseudogenes, unique DRB-STR patterns could be described that appear to be specific for a particular DRB haplotype. The aim of this workshop project was to approve and to qualify this simple typing protocol in a larger panel covering different European populations. All participants succeeded in correctly defining the DRB-STR amplicons varying from 135 to 222 base pair (bp) lengths. The panel of 101 samples covered 50 DRB alleles distributed over 37 different haplotypes as defined by exon 2 sequence-based typing. These haplotypes could be refined into 105 haplotypes by DRB-STR typing. Thus, discrimination of exon 2-identical DRB alleles was feasible, as well as the exact description of three different crossing-over events that resulted in the generation of hybrid DR region configurations. This typing procedure appears to be a quick and highly robust technique that can easily be performed by different laboratories, even without experience in microsatellite typing; thus, it is suitable for a variety of researchers in diverse research areas.

Comparative genetics of a highly divergent DRB microsatellite in different macaque species

The DRB region of the major histocompatibility complex (MHC) of cynomolgus and rhesus macaques is highly plastic, and extensive copy number variation together with allelic polymorphism makes it a challenging enterprise to design a typing protocol. All intact DRB genes in cynomolgus monkeys (Mafa) appear to possess a compound microsatellite, DRB-STR, in intron 2, which displays extensive length polymorphism. Therefore, this STR was studied in a large panel of animals, comprising pedigreed families as well. Sequencing analysis resulted in the detection of 60 Mafa-DRB exon 2 sequences that were unambiguously linked to the corresponding microsatellite. Its length is often allele specific and follows Mendelian segregation. In cynomolgus and rhesus macaques, the nucleotide composition of the DRB-STR is in concordance with the phylogeny of exon 2 sequences. As in humans and rhesus monkeys, this protocol detects specific combinations of different DRB-STR lengths that are unique for each haplotype. In the present panel, 22 Mafa-DRB region configurations could be defined, which exceeds the number detected in a comparable cohort of Indian rhesus macaques. The results suggest that, in cynomolgus monkeys, even more frequently than in rhesus macaques, new haplotypes are generated by recombination-like events. Although both macaque species are known to share several identical DRB exon 2 sequences, the lengths of the corresponding microsatellites often differ. Thus, this method allows not only fast and accurate DRB haplotyping but may also permit discrimination between highly related macaque species.

Pinpointing a selective sweep to the chimpanzee MHC class I region by comparative genomics

Chimpanzees experienced a reduction of the allelic repertoire at the major histocompatibility complex (MHC) class I A and B loci, which may have been caused by a retrovirus belonging to the simian immunodeficiency virus (SIV) family. Extended MHC haplotypes were defined in a pedigreed chimpanzee colony. Comparison of genetic variation at microsatellite markers mapping inside and outside the Mhc region was carried out in humans and chimpanzees to investigate the genomic extent of the repertoire reduction. Multilocus demographic analyses underscored that chimpanzees indeed experienced a selective sweep that mainly targeted the chromosomal segment carrying the Mhc class I region. Probably due to genetic linkage, the sweep also affected other polymorphic loci, mapping in the close vicinity of the Mhc class I region genes. Nevertheless, although the allelic repertoire at particular Mhc class I and II loci appears to be limited, naturally occurring recombination events allowed the establishment of haplotype diversity after the sweep. However, recombination did not have sufficient time to erase the signal of the selective sweep.

Sequence-dependent effect of interruptions on microsatellite mutation rate in mismatch repair-deficient human cells

Although microsatellite mutation rates generally increase with increasing length of the repeat tract, interruptions in a microsatellite may stabilize it. We have performed a direct analysis of the effect of microsatellite interruptions on mutation rate and spectrum in cultured mammalian cells. Two mononucleotide sequences (G(17) and A(17)) and a dinucleotide [(CA)(17)] were compared with interrupted repeats of the same size and with sequences of 8 repeat units. MMR-deficient (MMR(-)) cells were used for these studies to eliminate effects of this repair process. Mutation rates were determined by fluctuation analysis on cells containing a microsatellite sequence at the 5' end of an antibiotic-resistance gene; the vector carrying this sequence was integrated in the genome of the cells. In general, interrupted sequences had lower mutation rates than perfect ones of the same size, but the magnitude of the difference was dependent upon the sequence of the interrupting base(s). Some interrupted repeats had mutation rates that were lower than those of perfect sequences of the same length but similar to those of half the length. This suggests that interrupting bases effectively divide microsatellites into smaller repeat runs with mutational characteristics different from those of the corresponding full-length microsatellite. We conclude that interruptions decrease microsatellite mutation rate and influence the spectrum of frameshift mutations. The sequence of the interrupting base(s) determines the magnitude of the effect on mutation rate.

A highly divergent microsatellite facilitating fast and accurate DRB haplotyping in humans and rhesus macaques

The DRB region of the MHC in primate species is known to display abundant region configuration polymorphism with regard to the number and content of genes present per haplotype. Furthermore, depending on the species studied, the different DRB genes themselves may display varying degrees of allelic polymorphism. Because of this combination of diversity (differential gene number) and polymorphism (allelic variation), molecular typing methods for the primate DRB region are cumbersome. All intact DRB genes present in humans and rhesus macaques appear to possess, however, a complex and highly divergent microsatellite. Microsatellite analysis of a sizeable panel of outbred rhesus macaques, covering most of the known Mamu-DRB haplotypes, resulted in the definition of unique genotyping patterns that appear to be specific for a given haplotype. Subsequent examination of a representative panel of human cells illustrated that this approach also facilitates high-resolution HLA-DRB typing in an easy, quick, and reproducible fashion. The genetic composition of this complex microsatellite is shown to be in concordance with the phylogenetic relationships of various HLA-DRB and Mamu-DRB exon 2 gene/lineage sequences. Moreover, its length variability segregates with allelic variation of the respective gene. This simple protocol may find application in a variety of research avenues such as transplantation biology, disease association studies, molecular ecology, paternity testing, and forensic medicine.

Estimation of the species-specific mutation rates at the DRB1 locus in humans and chimpanzee

To estimate the species-specific mutation rates at the DRB1 locus in humans and chimpanzee, we analyzed the nucleotide sequence of a 37.6-kb chimpanzee chromosomal segment containing the entire Patr-DRB1*0701 allele and the flanking nongenic region and we compared it with two corresponding human sequences containing the HLA-DRB1*070101 allele using the sequence of HLA-DRB1*04011 as an outgroup. Because the allelic pair of HLA-DRB1*070101 and Patr-DRB1*0701 shows the lowest number of substitutions between the two species, it appears that these sequences diverged close to the time of the humans-chimpanzee divergence (6 million years ago). Alignment of the nucleotide sequences for HLA-DRB1*070101 and Patr-DRB1*0701 alleles showed that they share a high degree of similarity, suggesting that the studied chromosomal segments with these sequences have not been subjected to recombination since the humans-chimpanzee divergence. Comparison of the flanking 10.6 kb of nongenic sequences revealed an average of 41.5 and 83 single nucleotide substitutions in humans and chimpanzee, respectively. Thus, the species-specific nucleotide substitution rates in the flanking nongenic region were estimated to be 6.53 x 10(-10) and 1.31 x 10(-9) per site per year in humans and chimpanzee, respectively. Unexpectedly, the estimated rate in humans was twofold lower than in chimpanzee (P < 10(-3), Tajima's relative rate test) and lower than the average substitution rate in the human genome. Because the nucleotide substitution rate in nongenic regions free from selection is expected to be equal to the mutation rate, the estimated substitution rate should correspond to the species-specific mutation rate at the DRB1 locus. Our results strongly suggest that the mutation rate at DRB1 locus differs among species.

Identification and analysis of MHC class II DRB1 (Patr-DRB1) alleles in chimpanzees

The MHC-DRB1 gene is known to display the most extensive allelic polymorphisms among MHC class II genes. We attempted the selective identification of chimpanzee (Pan troglodytes) DRB1 (Patr-DRB1) alleles using the polymerase chain reaction (PCR) technique in three steps: first, we performed Patr-DRB1*02 lineage-specific 8-kb PCR for *02 lineage detection in each chimpanzee; second, we performed 620-bp PCR for amplification of full-length exon 2; and finally, we carried out an insert check using the pattern of microsatellite repeat length variability. In the genomic DNA of 23 chimpanzees, nine Patr-DRB1 alleles containing two new alleles were detected. Our approach provides a relatively effective method of identifying Patr-DRB1 alleles in individual chimpanzees and should also contribute to our understanding of the features of MHC molecules in non-human primates.

Polymorphism scan for differences between transmitted and nontransmitted DRB1*030101 alleles outside of exon 2 for type 1 diabetes: The frequency of polymorphisms is similar

DRB1*030101 is a major genetic risk factor for type 1 diabetes mellitus (T1DM) and is the only DRB1*03 allele usually seen in T1DM probands. Approximately 16% of parental DRB1*030101 alleles were not transmitted to T1DM probands in our Genetics of Kidneys and Diabetes study trio families. We performed a polymorphism screen to determine whether variations exist in DRB1*030101 alleles outside of exon 2 that may modify risk for developing T1DM. A combination of long-range and sequence-specific priming polymerase chain reaction was used to amplify a hemizygous template from both transmitted and nontransmitted parental DRB1*030101 chromosomes. Exon 2 DRB1*030101-specific and flanking DRB1-specific primers amplified the entire genomic locus as a 10.6-kb 5' fragment and a 5.3-kb 3' fragment, respectively. All exons and intron/exon borders of introns 1 and 2, all of introns 3-5, and flanking regulatory regions of 32 transmitted and 31 nontransmitted alleles (99% power to detect a 5% minimal allele frequency) were analyzed through fluorescent DNA sequencing. The only polymorphic sites detected, a previously described intron 2 complex dinucleotide repeat and an additional complex repeat approximately 1.8 kb downstream of exon 6, do not significantly differ between T1DM patients and controls in this small data set.

Recent divergence of the HLA-DRB1*04 allelic lineage from the DRB1*0701 lineage after the separation of the human and chimpanzee species

Conventional phylogenetic trees for the human leukocyte antigen (HLA)-DRB1 alleles constructed by the neighbor-joining (Saitou and Nei 1987) and UPGMA (Sneath and Sokal 1973) methods using nucleotide sequences of the DRB1 alleles suggest that DRB1*0701 may have diverged from other DRB1 alleles before the separation of the human and chimpanzee species, because of a large number of nucleotide changes in DRB1*0701 compared with any of the other DRB1 alleles. Here we show new evidence that the haplotypes centering on DRB1*0701 and DRB1*04 alleles are the most homologous. This suggests that these haplotypes have derived from the common ancestral haplotype, and that they have likely retained complete linkage disequilibrium even after the divergence of the DRB1*0701 and DRB1*04 allelic lineages. Together with the corresponding haplotype carrying chimpanzee DRB1*0701, which has a high sequence homology to HLA-DRB1*0701, these haplotypes reveal that: (1) the DRB1*04 allelic lineage may have been generated from the DRB1*0701 lineage after the separation of the human and chimpanzee species; (2) the DRB1*04 allelic lineage possibly has a higher substitution rate of DRB1 compared with pseudogene and neutral region; (3) there could be a significant difference in the substitution rate of DRB1 between the DRB1*0701 and DRB1*04 allelic lineages. Based on the difference between the present and previous results, we would like to propose that phylogenetic studies using not only nucleotide sequences of the DRB1 alleles but also haplotypes centering on the alleles should be conducted for understanding detailed phylogenetic relationships of the DRB1 alleles.

Recombinant DRB sequences produced by mismatch repair of heteroduplexes during cloning in Escherichia coli

Recombinant chimeric sequences originating from a mixture of the sequences of two different alleles are frequently found after amplification and cloning in Escherichia coli of exon 2 of the major histocompatibility complex (MHC) DRB genes. Several authors have suggested that the recombinant molecules result from in vitro recombination during PCR; nevertheless, a clear experimental demonstration of this hypothesis is lacking. In order to understand the mechanism producing the chimeric sequences, we set up a simple experiment based on the different restriction patterns of parental and recombinant sequences. Our data demonstrate that in the analysed case most of the recombinant variants were not produced by in vitro recombination during PCR, but were the result of the mismatch repair of heteroduplex molecules during cloning in E. coli. The high mutation rate in the alpha-helix region of DRB expressed genes, both after cloning in E. coli and after the germ-line differentiation process in vertebrates, suggests that the observed mutations are the result of similar gene conversion processes, probably favoured by chi-dependent microrecombination events.

Genotyping for human whole-genome scans: past, present, and future

Efficient and effective whole-genome 10-cM short tandem repeat polymorphism (STRP) scans are now available. Doubling or tripling STRP density to an average spacing of 3-5 cM is readily achievable. However, if typing costs for diallelic polymorphisms can be brought close to, or preferably less than, one-third those of STRPs, then diallelics may gradually supplement or supplant STRPs in whole-genome scans. The power of higher density genome scans for gene map ping by association and for many other research and clinical applications is great. It would be wise to continue investing heavily for many years in genotyping technology.

Microsatellites in the HLA region: 1999 update

In this third update of a series of reviews on microsatellites in the HLA region or close to it we report 155 microsatellites, corresponding to 51 newly described markers, in addition to the 103 reported in the 1997 and 1998 reviews. This work is based both on a literature review and on data publicly available in molecular databases on the internet (http://www.gdb.org; http://bioinfo.weizmann.ac.il/cards/; http://cedar.genetics.soton.ac.uk/) up to September 1999. Thanks to numerous studies involving major histocompatibility complex (MHC) microsatellites, documentation on HLA region is proposed, including information on microsatellites described through MHC sequence projects and presenting documented location, polymorphism and amplification condition, together with additional information on previously described microsatellites when available and information on data in the literature regarding gametic associations between HLA region loci and alleles and microsatellite alleles. As basic information are presented various documents: i) a table showing the following characteristics of the 155 microsatellites: name, localisation, polymorphism, primer sequences, reference; ii) an integrated map of some HLA region genes and the 155 microsatellites considered; and iii) a summary table on HLA and microsatellites association patterns. In addition, an overview on HLA microsatellite analysis application is presented, with a special focus on disease genetics studies in the form of recent references where the use of microsatellites of the HLA region was a key tool. This review aims at providing the human immunogenetics community with a tool for helping optimal choice of microsatellites to be used in various studies.

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.