Contribution of gene conversion in the evolution of the human beta-like globin gene family

A 25-year odyssey of genomic technology advances and structural variant discovery

This perspective focuses on advances in genome technology over the last 25 years and their impact on germline variant discovery within the field of human genetics. The field has witnessed tremendous technological advances from microarrays to short-read sequencing and now long-read sequencing. Each technology has provided genome-wide access to different classes of human genetic variation. We are now on the verge of comprehensive variant detection of all forms of variation for the first time with a single assay. We predict that this transition will further transform our understanding of human health and biology and, more importantly, provide novel insights into the dynamic mutational processes shaping our genomes.

A New Mutation, Hb A2-Canakkale [δ10(A7)Ala→Val; HBD: c.32C>T], and Other Well-Known δ Variants Identified in a Selected Cohort with Low Hb A2 Levels

Hemoglobinopathies are the most common single-gene disorders, and β-thalassemia (β-thal) imposes a tremendous health burden on Turkey. Thus, premarital carrier screening is obligatory in Turkey, as it is in some other countries. As a result of this mandatory procedure, at routine clinical checkups, many individuals who had undergone premarital screening but did not have any clinical symptoms and/or hematological findings, have compulsorily been required to undergo further evaluation due to abnormal levels of hemoglobin (Hb) fractions (Hb A, Hb A₂ and Hb F). Many consequences, such as mutations in unrelated gene(s) or someone's nutritional status, have been reported to affect the Hb fractions levels. In the present study, we aimed to determine whether HBD has a molecular causative role in patients with low Hb A₂ levels (below 1.8%). The study was conducted with 20 individuals with low Hb A₂ levels who had applied to our outpatient clinic. All DNA samples were analyzed for the HBD gene. Nineteen of the 20 subjects were diagnosed to carry a mutation with one of four different δ-globin variants. Three of them had been described previously [Hb A₂-Yialousa (HBD: c.82G>T), Hb A₂-Bornova (HBD: c.350G>C) and Hb A₂-Yokoshima (HBD: c.77G>A)]. The novel [δ10(A7)Ala→Val, HBD: c.32C>T] mutation was defined as a new δ variant and reported to the HbVar database as Hb A₂-Canakkale. In conclusion, the molecular characterization of Hb A₂ low levels has been suggested to be significant for a definite diagnosis and counseling.

CRISPR/Cas9 increases mitotic gene conversion in human cells

Gene conversion is a process of transferring genetic material from one homologous sequence to another. Most reported gene conversions are meiotic although mitotic gene conversion is also described. When using CRISPR/Cas9 to target the human hemoglobin subunit beta (HBB) gene, hemoglobin subunit delta (HBD) gene footprints were observed in HBB gene. However, it is unclear whether these were the results of gene conversion or PCR-mediated sequence shuffling between highly homologous sequences. Here we provide evidence that the HBD footprints in HBB were indeed results of gene conversion. We demonstrated that the CRISPR/Cas9 facilitated unidirectional sequence transfer from the homologous gene without double-strand breaks (DSB) to the one with DSBs, and showed that the rates of HBD footprint in HBB were positively correlated to the HBB insertion and deletion rates. We further showed that when targeting HBD gene, HBB footprints could also be observed in HBD gene. The mitotic gene conversion was observed not only in immortalized HEK293T cells, but also in human primary cells. Our work reveals mitotic gene conversion as an often overlooked effect of CRISPR/Cas9-mediated genome editing.

Emergence of a Chimeric Globin Pseudogene and Increased Hemoglobin Oxygen Affinity Underlie the Evolution of Aquatic Specializations in Sirenia

As limits on O₂ availability during submergence impose severe constraints on aerobic respiration, the oxygen binding globin proteins of marine mammals are expected to have evolved under strong evolutionary pressures during their land-to-sea transition. Here, we address this question for the order Sirenia by retrieving, annotating, and performing detailed selection analyses on the globin repertoire of the extinct Steller’s sea cow (Hydrodamalis gigas), dugong (Dugong dugon), and Florida manatee (Trichechus manatus latirostris) in relation to their closest living terrestrial relatives (elephants and hyraxes). These analyses indicate most loci experienced elevated nucleotide substitution rates during their transition to a fully aquatic lifestyle. While most of these genes evolved under neutrality or strong purifying selection, the rate of nonsynonymous/synonymous replacements increased in two genes (Hbz-T1 and Hba-T1) that encode the α-type chains of hemoglobin (Hb) during each stage of life. Notably, the relaxed evolution of Hba-T1 is temporally coupled with the emergence of a chimeric pseudogene (Hba-T2/Hbq-ps) that contributed to the tandemly linked Hba-T1 of stem sirenians via interparalog gene conversion. Functional tests on recombinant Hb proteins from extant and ancestral sirenians further revealed that the molecular remodeling of Hba-T1 coincided with increased Hb–O₂ affinity in early sirenians. Available evidence suggests that this trait evolved to maximize O₂ extraction from finite lung stores and suppress tissue O₂ offloading, thereby facilitating the low metabolic intensities of extant sirenians. In contrast, the derived reduction in Hb–O₂ affinity in (sub)Arctic Steller’s sea cows is consistent with fueling increased thermogenesis by these once colossal marine herbivores.

Whole-Genome-Sequence-Based Haplotypes Reveal Single Origin of the Sickle Allele during the Holocene Wet Phase

Five classical designations of sickle haplotypes are made on the basis of the presence or absence of restriction sites and are named after the ethno-linguistic groups or geographic regions from which the individuals with sickle cell anemia originated. Each haplotype is thought to represent an independent occurrence of the sickle mutation rs334 (c.20A>T [p.Glu7Val] in HBB). We investigated the origins of the sickle mutation by using whole-genome-sequence data. We identified 156 carriers from the 1000 Genomes Project, the African Genome Variation Project, and Qatar. We classified haplotypes by using 27 polymorphisms in linkage disequilibrium with rs334. Network analysis revealed a common haplotype that differed from the ancestral haplotype only by the derived sickle mutation at rs334 and correlated collectively with the Central African Republic (CAR), Cameroon, and Arabian/Indian haplotypes. Other haplotypes were derived from this haplotype and fell into two clusters, one composed of Senegal haplotypes and the other composed of Benin and Senegal haplotypes. The near-exclusive presence of the original sickle haplotype in the CAR, Kenya, Uganda, and South Africa is consistent with this haplotype predating the Bantu expansions. Modeling of balancing selection indicated that the heterozygote advantage was 15.2%, an equilibrium frequency of 12.0% was reached after 87 generations, and the selective environment predated the mutation. The posterior distribution of the ancestral recombination graph yielded a sickle mutation age of 259 generations, corresponding to 7,300 years ago during the Holocene Wet Phase. These results clarify the origin of the sickle allele and improve and simplify the classification of sickle haplotypes.

The genetic basis and evolution of red blood cell sickling in deer

Crescent-shaped red blood cells, the hallmark of sickle-cell disease, present a striking departure from the biconcave disc shape normally found in mammals. Characterized by increased mechanical fragility, sickled cells promote haemolytic anaemia and vaso-occlusions and contribute directly to disease in humans. Remarkably, a similar sickle-shaped morphology has been observed in erythrocytes from several deer species, without obvious pathological consequences. The genetic basis of erythrocyte sickling in deer, however, remains unknown. Here, we determine the sequences of human β-globin orthologues in 15 deer species and use protein structural modelling to identify a sickling mechanism distinct from the human disease, coordinated by a derived valine (E22V) that is unique to sickling deer. Evidence for long-term maintenance of a trans-species sickling/non-sickling polymorphism suggests that sickling in deer is adaptive. Our results have implications for understanding the ecological regimes and molecular architectures that have promoted convergent evolution of sickling erythrocytes across vertebrates.

Patterns of Inter-Chromosomal Gene Conversion on the Male-Specific Region of the Human Y Chromosome

The male-specific region of the human Y chromosome (MSY) is characterized by the lack of meiotic recombination and it has long been considered an evolutionary independent region of the human genome. In recent years, however, the idea that human MSY did not have an independent evolutionary history begun to emerge with the discovery that inter-chromosomal gene conversion (ICGC) can modulate the genetic diversity of some portions of this genomic region. Despite the study of the dynamics of this molecular mechanism in humans is still in its infancy, some peculiar features and consequences of it can be summarized. The main effect of ICGC is to increase the allelic diversity of MSY by generating a significant excess of clustered single nucleotide polymorphisms (SNPs) (defined as groups of two or more SNPs occurring in close proximity and on the same branch of the Y phylogeny). On the human MSY, 13 inter-chromosomal gene conversion hotspots (GCHs) have been identified so far, involving donor sequences mainly from the X-chromosome and, to a lesser extent, from autosomes. Most of the GCHs are evolutionary conserved and overlap with regions involved in aberrant X–Y crossing-over. This review mainly focuses on the dynamics and the current knowledge concerning the recombinational landscape of the human MSY in the form of ICGC, on how this molecular mechanism may influence the evolution of the MSY, and on how it could affect the information enclosed within a genomic region which, until recently, appeared to be an evolutionary independent unit.

Y chromosome palindromes and gene conversion

The presence of large and near-identical inverted repeat sequences (called palindromes) is a common feature of the constitutively haploid sex chromosomes of different species. Despite the fact palindromes originated in a non-recombining context, they have evolved a strong recombinational activity in the form of abundant arm-to-arm gene conversion. Their independent appearance in different species suggests they can have a profound biological significance that has yet to be fully clarified. It has been theorized that natural selection may have favored palindromic organization of male-specific genes and that the establishment of intra-palindrome gene conversion has strong adaptive significance. Arm-to-arm gene conversion allows the efficient removal of deleterious mutations, increases the fixation rate of beneficial mutations and has played an important role in modulating the equilibrium between gene loss and acquisition during Y chromosome evolution. Additionally, a palindromic organization of duplicates could favor the formation of unusual chromatin structures and could optimize the use of gene conversion as a mechanism to maintain the structural integrity of male-specific genes. In this review, we describe the structural features of palindromes on mammalian sex chromosomes and summarize different hypotheses regarding palindrome evolution and the functional benefits of arm-to-arm gene conversion on the unique haploid portion of the nuclear genome.

Hemoglobin E in Northeast India: A review on its origin, distribution, migration and health implication

A systematic review of the studies on hemoglobin E in Northeast India has been carried out to understand the magnitude of research undertaken on this aspect during the last seven decades. Owing to the high prevalence of hemoglobin E in this part of India different authors have studied this hemoglobin from different perspectives and found conflicting results. However a systematic review of such studies is lacking from a holistic point of view. Most of the epidemiological, in vitro as well as in vivo studies show signatures of selection with this hemoglobin locus. However, how this polymorphism is maintained at different rates at different geographical region is still a matter of contention. This review will fill the gap from all perspectives starting from the frequency distribution of hemoglobin E and its spread in different parts of Northeast India, its relationship with malaria hypothesis, the population migration, population affinity and most importantly the health implication arising out of it. A probable origin of hemoglobin E among an Austroasiatic population of Northeast India has been postulated with the help of advance molecular anthropological knowledge like the deep rooted markers of mt DNA and Y-chromosome haplotypes.

Understanding the contrasting spatial haplotype patterns of malaria-protective β-globin polymorphisms

The malaria-protective β-globin polymorphisms, sickle-cell (β(S)) and β(0)-thalassaemia, are canonical examples of human adaptation to infectious disease. Occurring on distinct genetic backgrounds, they vary markedly in their patterns of linked genetic variation at the population level, suggesting different evolutionary histories. β(S) is associated with five classical restriction fragment length polymorphism haplotypes that exhibit remarkable specificity in their geographical distributions; by contrast, β(0)-thalassaemia mutations are found on haplotypes whose distributions overlap considerably. Here, we explore why these two polymorphisms display contrasting spatial haplotypic distributions, despite having malaria as a common selective pressure. We present a meta-population genetic model, incorporating individual-based processes, which tracks the evolution of β-globin polymorphisms on different haplotypic backgrounds. Our simulations reveal that, depending on the rate of mutation, a large population size and/or high population growth rate are required for both the β(S)- and the β(0)-thalassaemia-like patterns. However, whilst the β(S)-like pattern is more likely when population subdivision is high, migration low and long-distance migration absent, the opposite is true for β(0)-thalassaemia. Including gene conversion has little effect on the overall probability of each pattern; however, when inter-haplotype fitness variation exists, gene conversion is more likely to have contributed to the diversity of haplotypes actually present in the population. Our findings highlight how the contrasting spatial haplotype patterns exhibited by β(S) and β(0)-thalassaemia may provide important indications as to the evolution of these adaptive alleles and the demographic history of the populations in which they have evolved.

Evolutionary constraints in the β-globin cluster: the signature of purifying selection at the δ-globin (HBD) locus and its role in developmental gene regulation

Human hemoglobins, the oxygen carriers in the blood, are composed by two α-like and two β-like globin monomers. The β-globin gene cluster located at 11p15.5 comprises one pseudogene and five genes whose expression undergoes two critical switches: the embryonic-to-fetal and fetal-to-adult transition. HBD encodes the δ-globin chain of the minor adult hemoglobin (HbA2), which is assumed to be physiologically irrelevant. Paradoxically, reduced diversity levels have been reported for this gene. In this study, we sought a detailed portrait of the genetic variation within the β-globin cluster in a large human population panel from different geographic backgrounds. We resequenced the coding and noncoding regions of the two adult β-globin genes (HBD and HBB) in European and African populations, and analyzed the data from the β-globin cluster (HBE, HBG2, HBG1, HBBP1, HBD, and HBB) in 1,092 individuals representing 14 populations sequenced as part of the 1000 Genomes Project. Additionally, we assessed the diversity levels in nonhuman primates using chimpanzee sequence data provided by the PanMap Project. Comprehensive analyses, based on classic neutrality tests, empirical and haplotype-based studies, revealed that HBD and its neighbor pseudogene HBBP1 have mainly evolved under purifying selection, suggesting that their roles are essential and nonredundant. Moreover, in the light of recent studies on the chromatin conformation of the β-globin cluster, we present evidence sustaining that the strong functional constraints underlying the decreased contemporary diversity at these two regions were not driven by protein function but instead are likely due to a regulatory role in ontogenic switches of gene expression.

Assessment of anti-recombination and double-strand break-induced gene conversion in human cells by a chromosomal reporter

Gene conversion is one of the frequent end results of homologous recombination, and it often underlies the inactivation of tumor suppressor genes in cancer cells. Here, we have developed an integrated assay system that allows simultaneous examination of double-strand break (DSB)-induced gene conversion events at the site of a DSB (proximal region) and at a surrounding region ~1 kb away from the break (distal region). Utilizing this assay system, we find that gene conversion events at the proximal and distal regions are relatively independent of one another. The results also indicate that synthesis-dependent strand annealing (SDSA) plays a major role in DSB-induced gene conversion. In addition, our current study has demonstrated that hMLH1 plays an essential role in anti-recombination and gene conversion. Specifically, the anti-recombination activity of hMLH1 is partially dependent on its interaction with hMRE11. Our data suggests that the role of hMLH1 and hMRE11 in the process of gene conversion is complex, and these proteins play different roles in DSB-induced proximal and distal gene conversions. In particular, the involvement of hMLH1 and hMRE11 in the distal gene conversion requires both hMSH2 and heteroduplex formation.

Genetic flux between h1 and h2 haplotypes of the 17q21.31 inversion in European population

The chromosome 17q21.31 inversion is a 900-kb common structural polymorphism found primarily in European population. Although the genetic flux within inversion region was assumed to be considerable suppressed, it is still unclear about the details of genetic exchange between the H1 (non-inverted sequence) and H2 (inverted sequence) haplotypes of this inversion. Here we describe a refined map of genetic exchanges between pairs of gene arrangements within the 17q21.31 region. Using HapMap phase II data of 1,546 single nucleotide polymorphisms, we successfully deduced 96 H1 and 24 H2 haplotypes in European samples by neighbor-joining tree reconstruction. Furthermore, we identified 15 and 26 candidate tracts with reciprocal and non-reciprocal genetic exchanges, respectively. In all 15 regions harboring reciprocal exchange, haplotypes reconstructed by clone sequencing did not support these exchange events, suggesting that such signals of exchange between two sister chromosomes in certain heterozygous individual were caused by phasing error regions. On the other hand, the finished clone sequencing across 4 of 26 tracts with non-reciprocal genetic flux confirmed that this kind of genetic exchange was caused by gene conversion. In summary, as crossover between pairs of gene arrangements had been considerably suppressed, gene conversion might be the most important mechanism for genetic exchange at 17q21.31.

Identical mutations in the paralogous human γ-globin genes leading to hemoglobin variants and nondeletional hereditary persistence of fetal hemoglobin

The human fetal globin genes are highly similar at the DNA sequence level, resulting in a single amino acid difference between the (G)γ- and (A)γ-globin chains. A large proportion of hemoglobin (Hb) variants of the (G)γ- and (A)γ-globin chains result from an identical mutation in the HBG2 and HBG1 genes, respectively, while the same is true for a fraction of mutations leading to nondeletional hereditary persistence of fetal Hb (HPFH). In particular, 11 different Hb variants result from identical mutations on either one of the two human γ-globin paralogous genes, while seven other promoter substitutions result either in nondeletional HPFH or are benign polymorphisms. In the former case, the percentage of the Hb variants due to an HBG2 gene mutation was significantly higher than the percentage of Hb variants due to the same HBG1 gene mutation, following the (G)γ/(A)γ-globin chain ratio seen in wild-type individuals. These γ-globin chain variants have most likely occurred via recurrent mutations, gene conversion events or both and, contrary to the situation observed in the human α-globin genes, these mutations lead to distinct variant Hb molecules.

A de novo unequal cross-over mutation between CYP11B1 and CYP11B2 genes causes familial hyperaldosteronism type I

Familial hyperaldosteronism type I (FH-I) is an autosomal dominant disorder caused by an unequal cross-over of the gene encoding steroid 11β-hydroxylase (CYP11B1) and aldosterone synthase (CYP11B2), giving rise to a chimeric CYP11B1/CYP11B2 gene that displays aldosterone synthase activity regulated by ACTH instead of angiotensin II.

AIM

To report an unprecedented case of a de novo unequal crossover mutation between CYP11B1 and CYP11B2 genes causing FH-I.

PATIENTS AND METHODS

The index case is a 45-yr-old Chilean male diagnosed with primary aldosteronism (PA). All family members were also studied: his biological parents, 1 brother, 6 sisters, 2 daughters, and 1 son. Plasma renin activity, serum aldosterone, and its ratio were measured in all patients. Genetic analyses were performed using long-extension PCR (XL-PCR), DNA sequencing and Southern blot methods.

RESULTS

PA was diagnosed for the index case, 1 of his daughters, his son but not for his parents or siblings. XLPCR and Southern blotting demonstrated the presence of the chimeric CYP11B1/CYP11B2 gene solely in PA-affected subjects, suggesting a case of a de novo mutation. Sequence analysis showed the unequal cross-over CYP11B1/CYP11B2 at intron 2 (c.2600-273 CYP11B2). We also identified a polymorphism at the same intron (c.2600-145C>A CYP11B2) in the genome of the index case's father.

CONCLUSION

We describe an unprecedented case of unequal cross-over mutation for the chimeric CYP11B1/CYP11B2 gene causing FH-I, which may be linked to a polymorphism in the index case's father germ line.

Recombination-associated sequence homogenization of neighboring Alu elements: signature of nonallelic gene conversion

Recently, researchers have begun to recognize that, in order to establish neutral models for disease association and evolutionary genomics studies, it is crucial to have a clear understanding of the genomic impact of nonallelic gene conversion. Drawing on previous successes in characterizing this phenomenon over protein-coding gene families, we undertook a computational analysis of neighboring Alu sequences in the genome scale. For this purpose, we developed adjusted comutation rate (aCMR), a novel statistical method measuring the excess number of identical point mutations shared by adjacent Alu sequences, vis-à-vis random pairs. Using aCMR, we uncovered a remarkable genome-wide sequence homogenization of neighboring Alus, with the strongest signal observed in the pseudoautosomal regions of the X and Y chromosomes. The magnitude of sequence homogenization between Alu pairs is greater with shorter interlocus distance, higher sequence identity, and parallel orientation. Moreover, shared substitutions show a strong directionality toward GC nucleotides, with multiple substitutions tending to cluster within the Alu sequence. Taken together, these observed recombination-associated sequence homogenization patterns are best explained by frequent ubiquitous gene conversion events between neighboring Alus. We believe that these observations help to illuminate the nature and impact of the enigmatic phenomenon of gene conversion.

Evolution of duplicated beta-globin genes and the structural basis of hemoglobin isoform differentiation in Mus

The functional diversification of multigene families may be strongly influenced by mechanisms of concerted evolution such as interparalog gene conversion. The beta-globin gene family of house mice (genus Mus) represents an especially promising system for evaluating the effects of gene conversion on the functional divergence of duplicated genes. Whereas the majority of mammalian species possess tandemly duplicated copies of the adult beta-globin gene that are identical in sequence, natural populations of house mice are often polymorphic for distinct two-locus haplotypes that differ in levels of functional divergence between duplicated beta-globin genes, HBB-T1 and HBB-T2. Here, we use a phylogenetic approach to unravel the complex evolutionary history of the HBB-T1 and HBB-T2 paralogs in a taxonomically diverse set of species in the genus Mus. The main objectives of this study were 1) to reconstruct the evolutionary history of the different HBB haplotypes of house mice, 2) to assess the role of recombinational exchange between HBB-T1 and HBB-T2 in promoting concerted evolution, 3) to assess the role of recombinational exchange between HBB-T1 and HBB-T2 in creating chimeric genes, and 4) to assess the structural basis of hemoglobin isoform differentiation in species that possess distinct HBB paralogs. Results of our phylogenetic survey revealed that the HBB-T1 and HBB-T2 genes in different species of Mus exhibit the full range of evolutionary outcomes with respect to levels of interparalog divergence. At one end of the spectrum, the two identical HBB paralogs on the Hbb(s) haplotype (shared by Mus domesticus, Mus musculus, and Mus spretus) represent a classic example of concerted evolution. At the other end of the spectrum, the two distinct HBB paralogs on the Hbb(d), Hbb(p), Hbb(w1), and Hbb(w2) haplotypes (shared by multiple species in the subgenus Mus) show no trace of gene conversion and are distinguished by a number of functionally important amino acid substitutions. Because the possession of distinct HBB paralogs expands the repertoire of functionally distinct hemoglobin isoforms that can be synthesized during fetal development and postnatal life, variation in the level of functional divergence between HBB-T1 and HBB-T2 may underlie important physiological variation within and among species.

Gene conversion transfers the GAF-A domain of phosphodiesterase TbrPDEB1 to one allele of TbrPDEB2 of Trypanosoma brucei

Background

Chromosome 9 of Trypanosoma brucei contains two closely spaced, very similar open reading frames for cyclic nucleotide specific phosphodiesterases TbrPDEB1 and TbrPDEB2. They are separated by 2379 bp, and both code for phosphodiesterases with two GAF domains in their N-terminal moieties and a catalytic domain at the C-terminus.

Methods and Findings

The current study reveals that in the Lister427 strain of T. brucei, these two genes have undergone gene conversion, replacing the coding region for the GAF-A domain of TbrPDEB2 by the corresponding region of the upstream gene TbrPDEB1. As a consequence, these strains express two slightly different versions of TbrPDEB2. TbrPDEB2a represents the wild-type phosphodiesterase, while TbrPDEB2b represents the product of the converted gene. Earlier work on the subcellular localization of TbrPDEB1 and TbrPDEB2 had demonstrated that TbrPDEB1 is predominantly located in the flagellum, whereas TbrPDEB2 partially locates to the flagellum but largely remains in the cell body. The current findings raised the question of whether this dual localization of TbrPDEB2 may reflect the two alleles. To resolve this, TbrPDEB2 of strain STIB247 that is homozygous for TbrPDEB2a was tagged in situ, and its intracellular localization was analyzed.

Conclusions

The results obtained were very similar to those found earlier with Lister427, indicating that the dual localization of TbrPDEB2 reflects its true function and is not simply due to the presence of the two different alleles. Notably, the gene conversion event is unique for the Lister427 strain and all its derivatives. Based on this finding, a convenient PCR test has been developed that allows the stringent discrimination between Lister-derived strains that are common in many laboratories and other isolates. The technique is likely very useful to resolve questions about potential mix-ups of precious field isolates with the ubiquitous Lister strain.

Cyclic nucleotide specific phosphodiesterases are important regulators of cyclic nucleotide signalling in eukaryotes. In many organisms, including humans and trypanosomes, some of these enzymes contain specific domains (GAF domains) that bind cyclic nucleotides, and that are involved in the regulation of the catalytic domain. In the parasitic protozoon that causes human sleeping sickness, Trypanosoma brucei, two closely related phosphodiesterases each contain two such GAF domains, GAF-A and GAF-B. Their genes are tandemly located on chromosome 9, spaced by only a few thousand nucleotides. We here show that a gene conversion event has exchanged the region that codes for the GAF-A domain of the downstream gene by the closely similar corresponding sequence of the upstream gene. This domain exchange has no effect on intracellular localization of the two enzymes. The gene conversion event has occurred in one particular strain of trypanosomes (Lister427) and is found in all its derivatives, but not in any other strain or isolate. The presence or absence of this gene conversion represents a useful analytical marker for the stringent discrimination of Lister427 derivatives from other trypanosome strains.

Mutations in the paralogous human alpha-globin genes yielding identical hemoglobin variants

The human alpha-globin genes are paralogues, sharing a high degree of DNA sequence similarity and producing an identical alpha-globin chain. Over half of the alpha-globin structural variants reported to date are only characterized at the amino acid level. It is likely that a fraction of these variants, with phenotypes differing from one observation to another, may be due to the same mutation but on a different alpha-globin gene. There have been very few previous examples of hemoglobin variants that can be found at both HBA1 and HBA2 genes. Here, we report the results of a systematic multicenter study in a large multiethnic population to identify such variants and to analyze their differences from a functional and evolutionary perspective. We identified 14 different Hb variants resulting from identical mutations on either one of the two human alpha-globin paralogue genes. We also showed that the average percentage of hemoglobin variants due to a HBA2 gene mutation (alpha2) is higher than the percentage of hemoglobin variants due to the same HBA1 gene mutation (alpha1) and that the alpha2/alpha1 ratio varied between variants. These alpha-globin chain variants have most likely occurred via recurrent mutations, gene conversion events, or both. Based on these data, we propose a nomenclature for hemoglobin variants that fall into this category.

Hb Gerland [alpha 55(E4)Val-->Ala]: a mutation found on the alpha1-globin gene

The Hb Gerland [alpha 55(E4)Val-->Ala] mutation has been described in the alpha2-globin gene. We report here the same mutation in the paralogous alpha1-globin gene. This variant was found in a healthy 18-month-old boy of Chinese origin. Abnormal hemoglobin (Hb) fractions were visible on isoelectric focusing (IEF) and cation exchange high performance liquid chromatography (HPLC), with elution patterns differing from one system to another. Direct sequencing of the alpha-globin genes revealed a GTT>GCT (Val-->Ala) transversion at codon 55 of the alpha1-globin gene. We propose to name the variant encoded by the alpha1-globin gene Hb Gerland [A1], and the variant that is encoded by the alpha2-globin gene Hb Gerland [A2].

Ectopic gene conversions in the human genome

We used the GENCONV method to characterize the gene conversions that occurred amongst the 1434 protein coding human gene families with three or more genes. Conversions occur at a frequency of 0.88% (483 conversion events/55,050 gene pairs compared) and have an average length of 371+/-752 bp (+/-standard deviation). Both the size and the frequency of conversions are positively correlated with the similarity of the sequences involved in these conversions. The frequency of conversions and the local recombination rate are also positively correlated. Intrachromosomal conversions are almost 5 times more frequent than interchromosomal conversions and the frequency of intrachromosomal conversions increases as the distance between genes decreases. However, the higher frequency of conversions between nearby genes with the same transcriptional orientation is due to the fact that most functional duplicated genes are found next to one another and in the same transcriptional orientation. The average length of a conversion spanning only an intron region is significantly smaller than conversions spanning both exons and introns or only exons. This suggests that the smaller degree of sequence similarity of introns limits the size of conversions between duplicated human genes. The significant excess of conversions at the 3'-end of human genes suggests that incomplete cDNA molecules are often involved in conversions with chromosomal gene copies.

Gene conversion: mechanisms, evolution and human disease

Gene conversion, one of the two mechanisms of homologous recombination, involves the unidirectional transfer of genetic material from a 'donor' sequence to a highly homologous 'acceptor'. Considerable progress has been made in understanding the molecular mechanisms that underlie gene conversion, its formative role in human genome evolution and its implications for human inherited disease. Here we assess current thinking about how gene conversion occurs, explore the key part it has played in fashioning extant human genes, and carry out a meta-analysis of gene-conversion events that are known to have caused human genetic disease.

Characterization of the −16C>G sequence variation in the promoters of both HBG1 and HBG2: Convergent evolution of the human γ-globin genes

We encountered a homozygous -16C>G mutation in cis at identical positions in the promoters of both human gamma-globin genes in a subject who was also homozygous for Hemoglobin C (HbC). Subsequent analysis of normal control individuals of African American ancestry revealed that both mutations were always present in cis with an allelic frequency of 3%. Furthermore, 10 out of 11 HbC subjects carried the -16C>G sequence variations, suggesting an association with HbC. The -16C>G mutation disrupts a putative CACCC box positioned between the TATA box and the transcriptional start site. However, the absence of high levels of HbF in HbC subjects homozygous and heterozygous for the -16C>G sequence variation suggested no effect of this mutation on gamma-globin gene expression in the adult stage of development. Further functional characterization by means of transient transfections in human erythroleukemic K562 cells showed that the -16C>G promoter sequence variation did not have an effect on gamma-globin expression in the fetal stage of development either. We therefore conclude that the -16C>G gamma-globin sequence variations are not beneficial to the clinical phenotype of HbC. The unique concurrent presence of this non-functional sequence variation is likely the result of a gene conversion event, and supports the concept of sequence homogenization between the two human gamma-globin genes.

Large-scale population genetic analysis for hemoglobinopathies reveals different mutation spectra in central Greece compared to the rest of the country

We have undertaken a large population screening study to identify the molecular basis of hemoglobinopathies in the central Greece region. A total of 845 unrelated beta-thalassemia patients and alpha-, beta-, and deltabeta-thalassemia carriers have been recruited and screened for mutations in the alpha- and beta-globin gene clusters. The alpha(-MED) deletion and the Turkish inversion/deletion are the most frequent genetic rearrangements leading to alpha- and deltabeta-thalassemia respectively, contrary to the situation in the rest of the country, while the beta -101 (C>T) promoter mutation is surprisingly frequent in the central part of Greece. Our data indicate that determination of mutation frequencies in different regions is vital for accurate provision of genetic services and counseling and for precise estimation of genetic diversity.

Comparative genomics and concerted evolution of beta-tubulin paralogs in Leishmania spp

Background

Tubulin isotypes and expression patterns are highly regulated in diverse organisms. The genome sequence of the protozoan parasite Leishmania major contains three distinct β-tubulin loci. To investigate the diversity of β-tubulin genes, we have compared the published genome sequence to draft genome sequences of two further species L. infantum and L. braziliensis. Untranscribed regions and coding sequences for each isoform were compared within and between species in relation to the known diversity of β-tubulin transcripts in Leishmania spp.

Results

All three β-tubulin loci were present in L. infantum and L. braziliensis, showing conserved synteny with the L. major sequence, hence confirming that these loci are paralogous. Flanking regions suggested that the chromosome 21 locus is an amastigote-specific isoform and more closely related (either structurally or functionally) to the chromosome 33 'array' locus than the chromosome 8 locus. A phylogenetic network of all isoforms indicated that paralogs from L. braziliensis and L. mexicana were monophyletic, rather than clustering by locus.

Conclusion

L. braziliensis and L. mexicana sequences appeared more similar to each other than each did to its closest relative in another species; this indicates that these sequences have evolved convergently in each species, perhaps through ectopic gene conversion; a process not yet evident among the more recently derived L. major and L. infantum isoforms. The distinctive non-coding regions of each β-tubulin locus showed that it is the regulatory regions of these loci that have evolved most during the diversification of these genes in Leishmania, while the coding regions have been conserved and concerted. The various loci in Leishmania satisfy a need for innovative expression of β-tubulin, rather than elaboration of its structural role.

Segmental duplications and gene conversion: Human luteinizing hormone/chorionic gonadotropin beta gene cluster

Segmental duplicons (>1 kb) of high sequence similarity (>90%) covering >5% of the human genome are characterized by complex sequence variation. Apart from a few well-characterized regions (MHC, beta-globin), the diversity and linkage disequilibrium (LD) patterns of duplicons and the role of gene conversion in shaping them have been poorly studied. To shed light on these issues, we have re-sequenced the human Luteinizing Hormone/Chorionic Gonadotropin beta (LHB/CGB) cluster (19q13.32) of three population samples (Estonians, Mandenka, and Han). The LHB/CGB cluster consists of seven duplicated genes critical in human reproduction. In the LHB/CGB region, high sequence diversity, concentration of gene-conversion acceptor sites, and strong LD colocalize with peripheral genes, whereas central loci are characterized by lower variation, gene-conversion donor activity, and breakdown of LD between close markers. The data highlight an important role of gene conversion in spreading polymorphisms among duplicon copies and generating LD around them. The directionality of gene-conversion events seems to be determined by the localization of a predicted recombination "hotspot" and "warm spot" in the vicinity of the most active acceptor genes at the periphery of the cluster. The data suggest that enriched crossover activity in direct and inverted segmental repeats is in accordance with the formation of palindromic secondary structures promoting double-strand breaks rather than fixed DNA sequence motifs. Also, this first detailed coverage of sequence diversity and structure of the LHB/CGB gene cluster will pave the way for studying the identified polymorphisms as well as potential genomic rearrangements in association with an individual's reproductive success.

Evolution of the ATP-binding cassette (ABC) transporter superfamily in vertebrates

The ATP-binding cassette (ABC) superfamily of genes encode membrane proteins that transport a diverse set of substrates across membranes. Mutations in ABC transporters cause or contribute to many different Mendelian and complex disorders including adrenoleukodystrophy, cystic fibrosis, retinal degeneration, hypercholesterolemia, and cholestasis. The genes play important roles in protecting organisms from xenobiotics and transport compounds across the intestine, blood-brain barrier, and the placenta. There are 48 ABC genes in the human genome divided into seven subfamilies based on amino acid sequence similarities and phylogeny. These seven subfamilies are represented in all eukaryotic genomes and are therefore of ancient origin. Sequencing the genomes of numerous vertebrate organisms has allowed the complement of ABC transporters to be characterized and the evolution of the genes to be assessed. Most ABC transporters are conserved in all vertebrates, but there are also several examples of recent duplication and gene loss. For genes with a conserved ortholog, animal models have been identified or developed that can be used to probe the function and regulation of selected genes. Genes that are restricted to a specific group of animals may represent specialized functions that could provide insight into unique biological properties of that organism. Further characterization of all ABC transporters from the human genome and from model organisms will lead to additional insights into normal physiology and human disease.

Phylogenetic analysis reveals a novel protein family closely related to adenosine deaminase

Adenosine deaminase (ADA) is a well-characterized enzyme involved in the depletion of adenosine levels. A group of proteins with similarity to ADA, the adenosine deaminase-related growth factors (ADGF; known as CECR1 in vertebrates), has been described recently in various organisms. We have determined the phylogenetic relationships of various gene products with significant amino acid similarity to ADA using parsimony and Bayesian methods, and discovered a novel paralogue, termed ADA-like (ADAL). The ADGF proteins share a novel amino acid motif, "MPKG," within which the proline and lysine residues are also conserved in the ADAL and ADA subfamilies. The significance of this new domain is unknown, but it is located just upstream of two ADA catalytic residues, of which all eight are conserved among the ADGF and ADAL proteins. This conservation suggests that ADGF and ADAL may share the same catalytic function as ADA, which has been proven for some ADGF members. These analyses also revealed that some genes previously thought to be classic ADAs are instead ADAL or ADGFs. We here define the ADGF, ADAL, ADA, adenine deaminase (ADE), and AMP deaminase (AMPD) groups as subfamilies of the adenyl-deaminase family. The availability of genomic data for the members of this family allowed us to reconstruct the intron evolution within the phylogeny and strengthen the introns-late hypothesis of the synthetic introns theory. This study shows that ADA activity is clearly more complex than once thought, perhaps involving a delicately balanced pattern of temporal and spatial expression of a number of paralogous proteins.

Molecular diagnosis of inherited disorders: lessons from hemoglobinopathies

Hemoglobinopathies constitute a major health problem worldwide, with a high carrier frequency, particularly in certain regions where malaria has been endemic. These disorders are characterized by a vast clinical and hematological phenotypic heterogeneity. Over 1,200 different genetic alterations that affect the DNA sequence of the human alpha-like (HBZ, HBA2, HBA1, and HBQ1) and beta-like (HBE1, HBG2, HBG1, HBD, and HBB) globin genes are mainly responsible for the observed clinical heterogeneity. These mutations, together with detailed information about the resulting phenotype, are documented in the globin locus-specific HbVar database. Family studies and comprehensive hematological analyses provide useful insights for accurately diagnosing thalassemia at the DNA level. For this purpose, numerous techniques can provide accurate, rapid, and cost-effective identification of the underlying genetic defect in affected individuals. The aim of this article is to review the diverse methodological and technical platforms available for the molecular diagnosis of inherited disorders, using thalassemia and hemoglobinopathies as a model. This article also attempts to shed light on issues closely related to thalassemia diagnostics, such as prenatal and preimplantation genetic diagnoses and genetic counseling, for better-quality disease management.

Origins of chromosomal rearrangement hotspots in the human genome: evidence from the AZFa deletion hotspots

BACKGROUND

The origins of the recombination hotspots that are a common feature of both allelic and non-allelic homologous recombination in the human genome are poorly understood. We have investigated, by comparative sequencing, the evolution of two hotspots of non-allelic homologous recombination on the Y chromosome that lie within paralogous sequences known to sponsor deletions resulting in male infertility.

RESULTS

These recombination hotspots are characterized by signatures of concerted evolution, which indicate that gene conversion between paralogs has been predominant in shaping their recent evolution. By contrast, the paralogous sequences that surround the hotspots exhibit little evidence of gene conversion. A second feature of these rearrangement hotspots is the extreme interspecific sequence divergence (around 2.5%) that places them among the most divergent orthologous sequences between humans and chimpanzees.

CONCLUSIONS

Several hominid-specific gene conversion events have rendered these hotspots better substrates for chromosomal rearrangements in humans than in chimpanzees or gorillas. Monte Carlo simulations of sequence evolution suggest that extreme sequence divergence is a direct consequence of gene conversion between paralogs. We propose that the coincidence of signatures of concerted evolution and recurrent breakpoints of chromosomal rearrangement (mapped at the sequence level) may enable the identification of putative rearrangement hotspots from analysis of comparative sequences from great apes.

Evidence against reciprocal recombination as the basis for tuf gene conversion in Salmonella enterica serovar Typhimurium

The duplicate tuf genes on the Salmonella enterica serovar Typhimurium chromosome co-evolve by a RecA-, RecB-dependent gene conversion mechanism. Gene conversion is defined as a non-reciprocal transfer of genetic information. However, in a replicating bacterial chromosome there is a possibility that a reciprocal genetic exchange between different tuf genes sitting on sister chromosomes could result in "apparent" gene conversion. We asked whether the major mechanism of tuf gene conversion was classical or apparent. We devised a genetic selection that allowed us to isolate and examine both expected products from a reciprocal recombination event between the tuf genes. Using this selection we tested within individual cultures for a correlation in the frequency of jackpots as expected if recombination were reciprocal. We found no correlation, either in the frequency of each type of recombinant product, or in the DNA sequences of the products resulting from each recombination event. We conclude that the evidence argues in favor of a non-reciprocal gene conversion mechanism as the basis for tuf gene co-evolution.

Dynamics of a human interparalog gene conversion hotspot

Gene conversion between paralogs can alter their patterns of sequence identity, thus obscuring their evolutionary relationships and affecting their propensity to sponsor genomic rearrangements. The details of this important process are poorly understood in the human genome because allelic diversity complicates the interpretation of interparalog sequence differences. Here we exploit the haploid nature of the Y chromosome, which obviates complicating interallelic processes, together with its known phylogeny, to understand the dynamics of conversion between two directly repeated HERVs flanking the 780-kb AZFa region on Yq. Sequence analysis of a 787-bp segment of each of the HERVs in 36 Y chromosomes revealed one of the highest nucleotide diversities in the human genome, as well as evidence of a complex patchwork of highly directional gene conversion events. The rate of proximal-to-distal conversion events was estimated as 2.4 x 10(-4) to 1.2 x 10(-3) per generation (3.9 x 10(-7) to 1.9 x 10(-6) per base per generation), and the distal-to-proximal rate as about one-twentieth of this. Minimum observed conversion tract lengths ranged from 1 to 158 bp and maximum lengths from 19 to 1365 bp, with an estimated mean of 31 bp. Analysis of great ape homologs shows that conversion in this hotspot has a deep evolutionary history.

Reciprocal and nonreciprocal recombination at the glucocerebrosidase gene region: implications for complexity in Gaucher disease

Gaucher disease results from an autosomal recessive deficiency of the lysosomal enzyme glucocerebrosidase. The glucocerebrosidase gene is located in a gene-rich region of 1q21 that contains six genes and two pseudogenes within 75 kb. The presence of contiguous, highly homologous pseudogenes for both glucocerebrosidase and metaxin at the locus increases the likelihood of DNA rearrangements in this region. These recombinations can complicate genotyping in patients with Gaucher disease and contribute to the difficulty in interpreting genotype-phenotype correlations in this disorder. In the present study, DNA samples from 240 patients with Gaucher disease were examined using several complementary approaches to identify and characterize recombinant alleles, including direct sequencing, long-template polymerase chain reaction, polymorphic microsatellite repeats, and Southern blots. Among the 480 alleles studied, 59 recombinant alleles were identified, including 34 gene conversions, 18 fusions, and 7 downstream duplications. Twenty-two percent of the patients evaluated had at least one recombinant allele. Twenty-six recombinant alleles were found among 310 alleles from patients with type 1 disease, 18 among 74 alleles from patients with type 2 disease, and 15 among 96 alleles from patients with type 3 disease. Several patients carried two recombinations or mutations on the same allele. Generally, alleles resulting from nonreciprocal recombination (gene conversion) could be distinguished from those arising by reciprocal recombination (crossover and exchange), and the length of the converted sequence was determined. Homozygosity for a recombinant allele was associated with early lethality. Ten different sites of crossover and a shared pentamer motif sequence (CACCA) that could be a hotspot for recombination were identified. These findings contribute to a better understanding of genotype-phenotype relationships in Gaucher disease and may provide insights into the mechanisms of DNA rearrangement in other disorders.

A review on the origin and spread of deleterious mutants of the beta-globin gene in Indian populations

Deleterious mutations of the human beta-globin gene are responsible for beta-thalassaemia and other haemoglobinopathies, which are the most common genetic diseases in Indian populations. A highly heterogeneous distribution of those mutations is observed in India and certain mutations are restricted to some extent to particular groups only. The reasons behind the geographical clustering and origin of the mutations in India is a highly debated issue and the evidence is conflicting. Our present article aims at tracing the origin of the deleterious beta-globin mutation and evaluates the role of different evolutionary forces responsible for the spread and present distribution of those mutations in Indian populations, using data from molecular biology and statistical methods. Mutations are generated essentially randomly, but "hot-spot" sites for mutation are reported for the beta-globin gene cluster, indicating sequence dependency of mutation. A single origin of a deleterious beta-globin mutation, followed by recombination (in a hot spot region) and/or interallelic gene conversion (within beta-globin gene) through time is the most plausible hypothesis to explain the association of those mutations with multiple haplotype backgrounds and frameworks. It is suggested that India is the place of origin of HbE and HbD mutations and that they dispersed to other parts of the would by migration. HbS mutants present in Indian populations are not of Middle East origin but rather a fresh mutation is the probable explanation for the prevalence among tribal groups. beta-thalassaemia represents a heterogeneous group of mutant alleles in India. Five common and twelve rare mutations have been reported in variable frequencies among different Indian populations. Gene flow of those mutant alleles from different populations of the world by political, military and commercial interactions possibly accounts for the heterogenous nature of beta-thalassaemia among Indians. A multiple allelic polymorphic system of the beta-globin gene exists in different populations. Dynamic interaction of the mutant alleles in the presence of different selective forces including falciparum malaria and biosocial patterns of Indian populations is discussed in order to explain the variable distribution and maintenance of those mutant alleles.

Evolutionary implications of three novel members of the human sarcomeric myosin heavy chain gene family

Sarcomeric myosin heavy chain (MyHC) is the major contractile protein of striated muscle. Six tandemly linked skeletal MyHC genes on chromosome 17 and two cardiac MyHC genes on chromosome 14 have been previously described in the human genome. We report the identification of three novel human sarcomeric MyHC genes on chromosomes 3, 7, and 20, which are notable for their atypical size and intron-exon structure. Two of the encoded proteins are structurally most like the slow-beta MyHC, whereas the third one is closest to the adult fast IIb isoform. Data from pairwise comparisons of aligned coding sequences imply the existence of ancestral genomes with four sarcomeric genes before the emergence of a dedicated smooth muscle MyHC gene. To further address the evolutionary relationships of the distinct sarcomeric and nonsarcomeric rod sequences, we have identified and further annotated human genomic DNA sequences corresponding to 14 class-II MyHCs. An extensive analysis provides a timeline for intron gain and loss, gene contraction and expansion, and gene conversion among genes encoding class-II myosins. One of the novel human genes is found to have introns at positions shared only with the molluscan catchin/MyHC gene, providing evidence for the structure of a pre-Cambrian ancestral gene.

The evolution of the vertebrate beta-globin gene promoter

Complexity analysis is capable of highlighting those gross evolutionary changes in gene promoter regions (loosely termed "promoter shuffling") that are undetectable by conventional DNA sequence alignment. Complexity analysis was therefore used here to identify the modular components (blocks) of the orthologous beta-globin gene promoter sequences of 22 vertebrate species, from zebrafish to humans. Considerable variation between the beta-globin gene promoters was apparent in terms of block presence/absence, copy number, and relative location. Some sequence blocks appear to be ubiquitous, whereas others are restricted to a specific taxon. Block similarities were also evident between the promoters of the paralogous human beta-like globin genes. It may be inferred that a wide variety of different mutational mechanisms have operated upon the beta-globin gene promoter over evolutionary time. Because these include gross changes such as deletion, duplication, amplification, elongation, contraction, and fusion, as well as the steady accumulation of single base-pair substitutions, it is clear that some redefinition of the term "promoter shuffling" is required. This notwithstanding, and as previously described for the vertebrate growth hormone gene promoter, the modular structure of the beta-globin promoter region and those of its paralogous counterparts have continually been rearranged into new combinations through the alteration, or shuffling, of preexisting blocks. Some of these changes may have had no influence on promoter function, but others could have altered either the level of gene expression or the responsiveness of the promoter to external stimuli. The comparative study of vertebrate beta-globin gene promoter regions described here confirms the generality of the phenomenon of sequence block shuffling and thus supports the view that it could have played an important role in the evolution of differential gene expression.