New genes originated via multiple recombinational pathways in the beta-globin gene family of rodents

Species differences in the size or membership composition of multigene families can be attributed to lineage-specific additions of new genes via duplication, losses of genes via deletion or inactivation, and the creation of chimeric genes via domain shuffling or gene fusion. In principle, it should be possible to infer the recombinational pathways responsible for each of these different types of genomic change by conducting detailed comparative analyses of genomic sequence data. Here, we report an attempt to unravel the complex evolutionary history of the beta-globin gene family in a taxonomically diverse set of rodent species. The main objectives were: 1) to characterize the genomic structure of the beta-globin gene cluster of rodents; 2) to assign orthologous and paralogous relationships among duplicate copies of beta-like globin genes; and 3) to infer the specific recombinational pathways responsible for gene duplications, gene deletions, and the creation of chimeric fusion genes. Results of our comparative genomic analyses revealed that variation in gene family size among rodent species is mainly attributable to the differential gain and loss of later expressed beta-globin genes via unequal crossing-over. However, two distinct recombinational mechanisms were implicated in the creation of chimeric fusion genes. In muroid rodents, a chimeric gamma/epsilon fusion gene was created by unequal crossing-over between the embryonic epsilon- and gamma-globin genes. Interestingly, this gamma/epsilon fusion gene was generated in the same fashion as the "anti-Lepore" 5'-delta-(beta/delta)-beta-3' duplication mutant in humans (the reciprocal exchange product of the pathological hemoglobin Lepore deletion mutant). By contrast, in the house mouse, Mus musculus, a chimeric beta/delta fusion pseudogene was created by a beta-globin --> delta-globin gene conversion event. Although the gamma/epsilon and beta/delta fusion genes share a similar chimeric gene structure, they originated via completely different recombinational pathways.     

Application of coalescent methods to reveal fine-scale rate variation and recombination hotspots

There has been considerable recent interest in understanding the way in which recombination rates vary over small physical distances, and the extent of recombination hotspots, in various genomes. Here we adapt, apply, and assess the power of recently developed coalescent-based approaches to estimating recombination rates from sequence polymorphism data. We apply full-likelihood estimation to study rate variation in and around a well-characterized recombination hotspot in humans, in the beta-globin gene cluster, and show that it provides similar estimates, consistent with those from sperm studies, from two populations deliberately chosen to have different demographic and selectional histories. We also demonstrate how approximate-likelihood methods can be used to detect local recombination hotspots from genomic-scale SNP data. In a simulation study based on 80 100-kb regions, these methods detect 43 out of 60 hotspots (ranging from 1 to 2 kb in size), with only two false positives out of 2000 subregions that were tested for the presence of a hotspot. Our study suggests that new computational tools for sophisticated analysis of population diversity data are valuable for hotspot detection and fine-scale mapping of local recombination rates.     

Absence of the TAP2 human recombination hotspot in chimpanzees

Recent experiments using sperm typing have demonstrated that, in several regions of the human genome, recombination does not occur uniformly but instead is concentrated in “hotspots” of 1–2 kb. Moreover, the crossover asymmetry observed in a subset of these has led to the suggestion that hotspots may be short-lived on an evolutionary time scale. To test this possibility, we focused on a region known to contain a recombination hotspot in humans, TAP2, and asked whether chimpanzees, the closest living evolutionary relatives of humans, harbor a hotspot in a similar location. Specifically, we used a new statistical approach to estimate recombination rate variation from patterns of linkage disequilibrium in a sample of 24 western chimpanzees (Pan troglodytes verus). This method has been shown to produce reliable results on simulated data and on human data from the TAP2 region. Strikingly, however, it finds very little support for recombination rate variation at TAP2 in the western chimpanzee data. Moreover, simulations suggest that there should be stronger support if there were a hotspot similar to the one characterized in humans. Thus, it appears that the human TAP2 recombination hotspot is not shared by western chimpanzees. These findings demonstrate that fine-scale recombination rates can change between very closely related species and raise the possibility that rates differ among human populations, with important implications for linkage-disequilibrium based association studies.     

Cladistic structure within the human Lipoprotein lipase gene and its implications for phenotypic association studies

Haplotype variation in 9.7 kb of genomic DNA sequence from the human lipoprotein lipase (LPL) gene was scored in three populations: African-Americans from Jackson, Mississippi (24 individuals), Finns from North Karelia, Finland (24), and non-Hispanic whites from Rochester, Minnesota (23). Earlier analyses had indicated that recombination was common but concentrated into a hotspot and that recurrent mutations at multiple sites may have occurred. We show that much evolutionary structure exists in the haplotype variation on either side of the recombinational hotspot. By peeling off significant recombination events from a tree estimated under the null hypothesis of no recombination, we also reveal some cladistic structure not disrupted by recombination during the time to coalescence of this variation. Additional cladistic structure is estimated to have emerged after recombination. Many apparent multiple mutational events at sites still remain after removing the effects of the detected recombination/gene conversion events. These apparent multiple events are found primarily at sites identified as highly mutable by previous studies, strengthening the conclusion that they are true multiple events. This analysis portrays the complexity of the interplay among many recombinational and mutational events that would be needed to explain the patterns of haplotype diversity in this gene. The cladistic structure in this region is used to identify four to six single-nucleotide polymorphisms (SNPs) that would provide disequilibrium coverage over much of this region. These sites may be useful in identifying phenotypic associations with variable sites in this gene. Evolutionary considerations also imply that the SNPs in the 3' region should have general utility in most human populations, but the 5' SNPs may be more population specific. Choosing SNPs at random would generally not provide adequate disequilibrium coverage of the sequenced region.     

Combining sperm typing and linkage disequilibrium analyses reveals differences in selective pressures or recombination rates across human populations

A previous polymorphism survey of the type 2 diabetes gene CAPN10 identified a segment showing an excess of polymorphism levels in all population samples, coinciding with localized breakdown of linkage disequilibrium (LD) in a sample of Hausa from Cameroon, but not in non-African samples. This raised the possibility that a recombination hotspot is present in all populations and we had insufficient power to detect it in the non-African data. To test this possibility, we estimated the crossover rate by sperm typing in five non-African men; these estimates were consistent with the LD decay in the non-African, but not in the Hausa data. Moreover, resequencing the orthologous region in a sample of Western chimpanzees did not show either an excess of polymorphism level or rapid LD decay, suggesting that the processes underlying the patterns observed in humans operated only on the human lineage. These results suggest that a hotspot of recombination has recently arisen in humans and has reached higher frequency in the Hausa than in non-Africans, or that there is no elevation in crossover rate in any human population, and the observed variation results from long-standing balancing selection.     

Evidence of Gene Conversion in the Evolutionary Process of the Codon 41/42 (-CTTT) Mutation Causing β-Thalassemia in Southern China

The 4-bp deletion (-CTTT) at codon 41/42 (CD41/42) of the human beta-globin gene represents one of the most common beta-thalassemia mutations in East Asia and Southeast Asia, which is historically afflicted with endemic malaria, thus hypothetically evolving under natural selection by malaria infection. To understand the evolutionary process of generating the beta(CD41/42) allele and its maintenance, including the effect of natural selection on the pattern of linkage disequilibrium (LD), we sequenced a 15.933-kb region spanning 20.693 kb of the beta-globin cluster surrounding the 4-bp deletion using a sample from a Chinese population consisting of 24 normal individuals and 16 heterozygotes for the deletion. Forty-nine polymorphic sites were found. Analysis of the data, using a variety of methods including formal population genetics analysis and visual approaches, suggests that the spread of the CD41/42 (-CTTT) deletion is most likely mediated by interallelic gene conversion, although independent deletions in different lineages are also possible. The neutrality test resulted in a significant positive Tajima's D for the beta-globin locus, which is consistent with the existence of balancing selection. This suggests that the 4-bp deletion that occurred at this locus may be an event that is subject to natural selection, due to malaria, which leads to the heterozygote advantage, spread widely with further help by conversion and migration. The evolutionary process of this mutant through gene conversion that could conceivably take place between the 4-bp deletion and the normal sequence in the respective region is discussed in detail.     

Nonuniform recombination within the human beta-globin gene cluster.

Population genetic analysis of 15 restriction site polymorphisms demonstrates nonuniform recombination within the human beta-globin gene cluster. These DNA polymorphisms show two clusters of high nonrandom associations, one 5' and another 3' to the beta-globin structural gene, with no significant linkage disequilibrium between the two clusters. The 5'- and 3'-association clusters are 34.6 kilobases (kb) and 19.4 kb long, respectively, and are separated by 9.1 kb of DNA immediately 5' to the beta-globin gene. For each of these three DNA regions, we have observed a relationship between nonrandom associations and physical distance between the polymorphisms. However, this relationship differed for each of these regions. On the assumption that the effective population size (Ne) is 5,000-50,000, we estimate the total recombination rate to be 0.0017%-0.0002% in the 5' cluster, 0.0931%-0.0093% in the 3' cluster, and 0.2912%-0.0219% in the 9.1-kb region between them. The beta cluster thus shows nonuniformity in recombination. Moreover, the recombination rate in the 9.1-kb DNA segment is 3-30 times greater than expected and is thus a hot spot for meiotic recombination.     

Molecular characterization of a meiotic recombinational hotspot enhancing homologous equal crossing-over.

We have cloned and sequenced a meiotic recombinational hotspot between the A beta 3 and A beta 2 genes in the major histocompatibility complex (MHC) of the mouse. This recombinational hotspot in the Mus musculus castaneus cas3 haplotype was previously localized to a region of 9.5 kb of DNA in which five independent crossing-over events occurred at the unusually high frequency of 0.6%. Aside from cas3, the hotspot appears to be absent in many other MHC haplotypes. We have now confined the five recombinational breakpoints to a stretch of 3.5 kb of DNA. From the nucleotide sequence around the recombinational breakpoints, determined in the parental cas3 and b haplotypes as well as for two recombinant haplotypes, we show that the two recombinant haplotypes were generated by homologous equal crossing-over and place the breakpoints within two non-overlapping stretches of 10 and 36 bp, respectively. Comparison of the DNA sequences of the hotspot-positive cas3 and the hotspot-negative b haplotypes reveals a number of differences, in particular, a CAGA-repeat sequence which is present in CAS3 in six, but only four copies in C57BL/6 DNA. This repeat sequence is reminiscent of one in a previously characterized hotspot in the E beta gene.     

Genetic control of sex-dependent meiotic recombination in the major histocompatibility complex of the mouse.

Meiotic recombination within the proximal region of the major histocompatibility complex (MHC) of the mouse is not random but occurs in clusters at certain restricted sites, so-called recombinational hotspots. The wm7 haplotype of the MHC, derived from the wild mouse, enhances recombination specifically during female meiosis within a fragment of 1.3 kb of DNA located between the A beta 3 and A beta 2 genes in genetic crosses with laboratory haplotypes. Previous studies revealed no significant strain differences in nucleotide sequences around the hotspot, irrespective of the ability of the strain to enhance the recombination. It appeared that a distant genetic element might, therefore, control the rate of recombination. In the present study, original recombinants whose breakpoints were defined by direct sequencing of PCR-amplified DNAs were tested for the rate of secondary recombination in the crosses with laboratory strains in order to determine the location of such a genetic element. The results clearly demonstrated that the chromosomal segment proximal to the hotspot is essential for enhancement of recombination. Moreover, the male recombination is suppressed by a segment distal to the hotspot.     

Molecular evolution of intergenic DNA in higher primates: pattern of DNA changes, molecular clock, and evolution of repetitive sequences

A 3.1-kb intergenic DNA fragment located between the psi beta-globin and delta-globin genes in the beta-globin gene cluster was cloned from gorilla, orangutan, rhesus monkey, and spider monkey, and the nucleotide sequence of each fragment was determined. The phylogeny of these four sequences, together with two previously published allelic sequences from humans and one from chimpanzee, was constructed, and the accumulation of mutations in the region was analyzed. The sites of base substitutions are not evenly distributed within the region: two Alu repeats have accumulated 0.21 + 0.02 substitutions/site with 0.15 + 0.008 substitutions/site in the remainder of the fragment. The occurrence of substitutions at neighboring sites is more frequent than would be expected if they were independent. The observed excesses disappear when ancestral -CG- dinucleotide sites are excluded. The phylogenetic relationships of the sequences indicate that the human sequence shares a most recent coancestor with the chimpanzee sequence. The data also show that great apes have accumulated fewer mutations in this part of the genome than has the rhesus monkey. The relative rates of accumulation of 12 kinds of nucleotide substitution in the region during primate evolution are asymmetric in the DNA strands. From these rates of accumulation, the origin of a simple stretch of sequence near the 3' end of the 3.1-kb fragment was deduced to be a sequence comprising 50% T and 50% C on one strand. The two oppositely oriented Alu sequences in the 3.1-kb region were inserted at their present positions before the divergence of the New-World monkeys from other lineages. Our analysis shows that the nucleotide sequences of the two Alu repeats in spider monkey are unexpectedly similar both to each other and to the deduced ancestral sequence of Alu repeats. The data suggest that there has been some type of recombinational event between the spider monkey Alu repeats but that it was not a simple gene conversion.      

Archaic African and Asian lineages in the genetic ancestry of modern humans.

A 3-kb region encompassing the beta-globin gene has been analyzed for allelic sequence polymorphism in nine populations from Africa, Asia, and Europe. A unique gene tree was constructed from 326 sequences of 349 in the total sample. New maximum-likelihood methods for analyzing gene trees on the basis of coalescence theory have been used. The most recent common ancestor of the beta-globin gene tree is a sequence found only in Africa and estimated to have arisen approximately 800,000 years ago. There is no evidence for an exponential expansion out of a bottlenecked founding population, and an effective population size of approximately 10,000 has been maintained. Modest differences in levels of beta-globin diversity between Africa and Asia are better explained by greater African effective population size than by greater time depth. There may have been a reduction of Asian effective population size in recent evolutionary history. Characteristically Asian ancestry is estimated to be older than 200,000 years, suggesting that the ancestral hominid population at this time was widely dispersed across Africa and Asia. Patterns of beta-globin diversity suggest extensive worldwide late Pleistocene gene flow and are not easily reconciled with a unidirectional migration out of Africa 100,000 years ago and total replacement of archaic populations in Asia.     

Comparison of population genetic patterns in two widespread freshwater mussels with contrasting life histories in western North America

We investigate population genetic structuring in Margaritifera falcata, a freshwater mussel native to western North America, across the majority of its geographical range. We find shallow rangewide genetic structure, strong population‐level structuring and very low population diversity in this species, using both mitochondrial sequence and nuclear microsatellite data. We contrast these patterns with previous findings in another freshwater mussel species group (Anodonta californiensis/A. nuttalliana) occupying the same continental region and many of the same watersheds. We conclude that differences are likely caused by contrasting life history attributes between genera, particularly host fish requirements and hermaphroditism. Further, we demonstrate the occurrence of a ‘hotspot’ for genetic diversity in both groups of mussels, occurring in the vicinity of the lower Columbia River drainage. We suggest that stream hierarchy may be responsible for this pattern and may produce similar patterns in other widespread freshwater species.     

Identification of human immunodeficiency virus type 1 envelope genes recombinant between subtypes B and F in two epidemiologically linked individuals from Brazil.

Sequence analysis of a human immunodeficiency virus type 1 env gene PCR amplified from a Brazilian woman's peripheral blood mononuclear cell DNA (sample RJIO1) showed that it was likely to have been derived from a double recombination event between human immunodeficiency virus type 1 subtypes B and F. The major portion of the gp120 coding sequence belonged to the B lineage, but a segment of the C2 to V3 region (approximately 135 nucleotides) clearly associated with sequences of the F lineage. The subtype F-like segment had 15 noncontiguous signature nucleotides in common with Brazilian subtype F sequences that were not found, or were rare, in subtype B sequences. In contrast, this same segment had only 3 signature nucleotides shared with subtype B sequences and not present in the Brazilian subtype F sequences. Phylogenetic analysis, amino acid signature pattern analysis, and the pattern of synonymous mutations all supported the hypothesis of a recombinational origin of the RJIO1 sequence. Related recombinant genes were also detected in peripheral blood mononuclear cell DNA obtained from the woman's recent sexual partner, indicating that the recombination event probably occurred at some previous time in the chain of virus transmission. Divergent viral sequences in the V3 region were found in the male sexual partner, while a relatively homogeneous viral population was detected in the woman, consistent with her recent infection.     

Genomic and phylogenic comparisons of the alpha-globin and beta-globin intergenic sequences between zebra fish (Danio rerio) and six closely related Cyprinindae species

The alpha-globin and beta-globin genes of the zebrafish are tightly linked on the same chromosome in a 3'-5' and 5'-3' configuration, respectively. Although the location of the controlling sequences has been mapped to the intergenic region, analysis determining the uniqueness of this unusual arrangement to zebrafish has not been undertaken. To explore this, we isolated, sequenced, and phylogenetically analyzed the intergenic region between globin gene families of seven Cyprinindae species including zebrafish. These species were grouped into an in group (immediate relatives, not so distant relatives), and an out group (distant relative). Cellulose acetate electrophoresis of hemoglobin (Hb) detected multiple variants in each species, but a band with electrophoretic mobility (EM) of 6.7 x 10(-5) cm(2).volt(-1).sec(-1) was shared between species. Polymerase chain reaction (PCR) amplification of the intergenic globin gene region also detected a 1.0-kb fragment that was repeated in the in group and a 1.2-kb fragment in the out group. Sequence comparison confirmed that the genetic orientation and controlling sequences location were conserved throughout this region in all seven species. This phylogenic footprinting indicated that the configuration was not exclusive to zebrafish. To confirm sequence alignment, maximum parsimony phylogenic analysis, was performed. Only one member of that group the giant danio, was not closely clustered, being located almost equidistance between the immediate relative and the species of the other clusters. This may represent an ancestral configuration prior to transposition of the alpha globin and beta-globin genes families to nonsynteny.     

Position effects in mice carrying a lacZ transgene in cis with the beta-globin LCR can be explained by a graded model.

We studied transgenic mice carrying the lacZ reporter gene linked to the erythroid-specific beta-globin promoter and beta-globin locus control region (LCR). Previously, we had demonstrated that the total level of expression of beta-galactosidase enzyme, which is the product of the lacZ gene, varies widely between different transgenic mice due to position effects at the sites of transgene integration. Here, using the X-gal based in situ assay for beta-galactosidase activity, we found that the percent erythroid cells that expressed the transgene also varied widely between the mice. Moreover, a kinetic analysis showed that the average beta-galactosidase content per expressing cell varied both between samples of different transgenic descent and between erythroid cells within each sample, demonstrating that the variable expression of this lacZ transgene was being controlled in a graded manner. These results suggest that the beta-globin LCR enhancers function through a graded model, which is described, rather than the binary mechanism that has been proposed previously for other enhancers.     

Common 5' beta-globin RFLP haplotypes harbour a surprising level of ancestral sequence mosaicism

Blocks of linkage disequilibrium (LD) in the human genome represent segments of ancestral chromosomes. To investigate the relationship between LD and genealogy, we analysed diversity associated with restriction fragment length polymorphism (RFLP) haplotypes of the 5' beta-globin gene complex. Genealogical analyses were based on sequence alleles that spanned a 12.2-kb interval, covering 3.1 kb around the psibeta gene and 6.2 kb of the delta-globin gene and its 5' flanking sequence known as the R/T region. Diversity was sampled from a Kenyan Luo population where recent malarial selection has contributed to substantial LD. A single common sequence allele spanning the 12.2-kb interval exclusively identified the ancestral chromosome bearing the "Bantu" beta(s) (sickle-cell) RFLP haplotype. Other common 5' RFLP haplotypes comprised interspersed segments from multiple ancestral chromosomes. Nucleotide diversity was similar between psibeta and R/T-delta-globin but was non-uniformly distributed within the R/T-delta-globin region. High diversity associated with the 5' R/T identified two ancestral lineages that probably date back more than 2 million years. Within this genealogy, variation has been introduced into the 3' R/T by gene conversion from other ancestral chromosomes. Diversity in delta-globin was found to lead through parts of the main genealogy but to coalesce in a more recent ancestor. The well-known recombination hotspot is clearly restricted to the region 3' of delta-globin. Our analyses show that, whereas one common haplotype in a block of high LD represents a long segment from a single ancestral chromosome, others are mosaics of short segments from multiple ancestors related in genealogies of unsuspected complexity.     

Linkage disequilibrium in the human insulin/insulin-like growth factor II region of human chromosome II.

Caucasian (N = 128) and Chinese (N = 84) subjects were typed for RFLPs in the insulin (INS)/insulin-like growth factor II (IGF2) region of chromosome 11. Both the analysis of extended haplotypes and the pairwise measures of linkage disequilibrium among the RFLPs indicate that there is extensive linkage disequilibrium in the INS/IGF2 region. The disequilibrium extends across the hypervariable region (HVR) located just 5' to the INS gene and encompasses a region of at least 40 kbp. Previous studies had suggested that linkage disequilibrium in the INS region was negligible and that this region may therefore contain a "recombinational hotspot" (Chakravarti et al. 1986). However, results of this and another recent study (Thompson et al. 1988) highlight the importance of the frequencies of associated alleles in the ability to detect linkage disequilibrium. Thus, the previous failure to detect disequilibrium in the INS region may have reflected a lack of power, rather than a true absence of disequilibrium in this region.     

Comparative molecular population genetics of the Xdh locus in the cactophilic sibling species Drosophila buzzatii and D. koepferae

The Xdh (rosy) gene is one of the best studied in the Drosophila genus from an evolutionary viewpoint. Here we analyze nucleotide variation in a 1875-bp fragment of the second exon of Xdh in Argentinian populations of the cactophilic D. buzzatii and its sibling D. koepferae. The major electrophoretic alleles of D. buzzatii not only lack diagnostic amino acids in the region studied but also differ on average from each other by four to 13 amino acid changes. Our data also suggest that D. buzzatii populations belonging to different phytogeographic regions are not genetically differentiated, whereas D. koepferae exhibits a significant pattern of population structure. The Xdh region studied is twice as polymorphic in D. buzzatii as in D. koepferae. Differences in historical population size or in recombinational environment between species could account for the differences in the level of nucleotide variation. In both species, the Xdh region exhibits a great number of singletons, which significantly departs from the frequency spectrum expected under neutrality for nonsynonymous sites and also for synonymous sites in D. buzzatii. These excesses of singletons could be the signature of a recent population expansion in D. buzzatii, whereas they may be simply explained as the result of negative selection in D. koepferae.