Gene Conversion Process Research Articles

Simulation of Y-chromosomal haplotype data

The non-recombining nature of the Y-chromosome determines the non-independence of alleles between loci. The evolution of short tandem repeat (STR) loci in the Y-chromosome is the result of different factors such as differential mutation rates, mutation modes, gene conversion, selection and demographic processes. The degree of correlation between loci is dependent on the magnitude of these processes. The simulation of data is a routine tool used for testing hypotheses in population and evolutionary studies. The most basic parameters hitherto used in lineage haplotype simulations are the allele frequency distributions and mutation rates, assuming either full independence or linkage between loci. In this study we introduce use of the Spearman correlation coefficient to estimate the degree of dependence between non-recombining loci. Then, both the interdependence between loci and the allele frequency distributions at multi-allelic loci are incorporated in an algorithm for simulating haplotypes. We illustrate the method using published and unpublished Y-chromosome STR data.

Heterothallism inSaccharomyces cerevisiaeisolates from nature: effect ofHOlocus on the mode of reproduction

Understanding the evolution of sex and recombination, key factors in the evolution of life, is a major challenge in biology. Studies of reproduction strategies of natural populations are important to complement the theoretical and experimental models. Fungi with both sexual and asexual life cycles are an interesting system for understanding the evolution of sex. In a study of natural populations of yeast Saccharomyces cerevisiae, we found that the isolates are heterothallic, meaning their mating type is stable, while the general belief is that natural S. cerevisiae strains are homothallic (can undergo mating-type switching). Mating-type switching is a gene-conversion process initiated by a site-specific endonuclease HO; this process can be followed by mother-daughter mating. Heterothallic yeast can mate with unrelated haploids (amphimixis), or undergo mating between spores from the same tetrad (intratetrad mating, or automixis), but cannot undergo mother-daughter mating as homothallic yeasts can. Sequence analysis of HO gene in a panel of natural S. cerevisiae isolates revealed multiple mutations. Good correspondence was found in the comparison of population structure characterized using 19 microsatellite markers spread over eight chromosomes and the HO sequence. Experiments that tested whether the mating-type switching pathway upstream and downstream of HO is functional, together with the detected HO mutations, strongly suggest that loss of function of HO is the cause of heterothallism. Furthermore, our results support the hypothesis that clonal reproduction and intratetrad mating may predominate in natural yeast populations, while mother-daughter mating might not be as significant as was considered.

Biased Gene Conversion and the Evolution of Mammalian Genomic Landscapes

Recombination is typically thought of as a symmetrical process resulting in large-scale reciprocal genetic exchanges between homologous chromosomes. Recombination events, however, are also accompanied by short-scale, unidirectional exchanges known as gene conversion in the neighborhood of the initiating double-strand break. A large body of evidence suggests that gene conversion is GC-biased in many eukaryotes, including mammals and human. AT/GC heterozygotes produce more GC- than AT-gametes, thus conferring a population advantage to GC-alleles in high-recombining regions. This apparently unimportant feature of our molecular machinery has major evolutionary consequences. Structurally, GC-biased gene conversion explains the spatial distribution of GC-content in mammalian genomes-the so-called isochore structure. Functionally, GC-biased gene conversion promotes the segregation and fixation of deleterious AT --> GC mutations, thus increasing our genomic mutation load. Here we review the recent evidence for a GC-biased gene conversion process in mammals, and its consequences for genomic landscapes, molecular evolution, and human functional genomics.

Phylogenetic relationships among diploid Aegilops species inferred from 5S rDNA units

Relationships among the currently recognized 11 diploid species within the genus Aegilops have been investigated. Sequence similarity analysis, based upon 363 sequenced 5S rDNA clones from 44 accessions plus 15 sequences retrieved from GenBank, depicted two unit classes labeled the long AE1 and short AE1. Several different analytical methods were applied to infer relationships within haplomes, between haplomes and among the species, including maximum parsimony and maximum likelihood analyses of consensus sequences, “total evidence” phylogeny analysis and “matrix representation with parsimony” analysis. None were able to depict suites of markers or unit classes that could discern among the seven haplomes as is observed among established haplomes in other genera within the tribe Triticeae; however, most species could be separated when displayed on gene trees. These results suggest that the haplomes currently recognized are so refined that they may be relegated as sub-haplomes or haplome variants. Amblyopyrum shares the same 5S rDNA unit classes with the diploid Aegilops species suggesting that it belongs within the latter. Comparisons of the Aegilops sequences with those of Triticum showed that the long AE1 unit class of Ae. tauschii shared the clade with the equivalent long D1 unit class, i.e., the putative D haplome donor, but the short AE1 unit class did not. The long AE1 unit class but not the short, of Ae. speltoides and Ae. searsii both share the clade with the previously identified long {S1 and long G1 unit classes meaning that both Aegilops species can be equally considered putative B haplome donors to tetraploid Triticum species. The semiconserved nature of the nontranscribed spacer in Aegilops and in Triticeae in general is discussed in view that it may have originated by processes of incomplete gene conversion or biased gene conversion or birth-and-death evolution.

Comment on "Human-Specific Gain of Function in a Developmental Enhancer"

Prabhakar et al. (Reports, 5 September 2008, p. 1346) argued that the conserved noncoding sequence HACNS1 has undergone positive selection and contributed to human adaptation. However, the pattern of substitution in HACNS1 is more consistent with the neutral process of biased gene conversion (BGC). The reported human-specific gain of function is likely due to the accumulation of deleterious mutations driven by BGC, not positive selection.

A Likely Pathway for Formation of Mobile Group I Introns

Mobile group I introns are RNA splicing elements that have been invaded by endonuclease genes. These endonucleases facilitate intron mobility by a unidirectional, duplicative gene-conversion process known as homing [1]. Survival of the invading endonuclease depends upon its ability to promote intron mobility. Therefore, the endonuclease must either quickly change its cleavage specificity to match the site of intron insertion, or it must already be preadapted to cleave this sequence. Here we show that the group I intron in the DNA polymerase gene of T7-like bacteriophage PhiI is mobile, dependent upon its intronic HNH homing endonuclease gene, I-TslI. We also show that gene 5.3 of phage T3, located adjacent to its intronless DNA polymerase gene, is a homologous homing endonuclease gene whose protein product initiates efficient spread of gene 5.3 into empty sites in related phages. Both of these endonucleases cleave intronless DNA polymerase genes at identical positions. This shared feature between an intronic and free-standing endonuclease is unprecedented. Based on this evidence, we propose that introns and their homing endonucleases evolve separately to target the same highly conserved sequences, uniting afterwards to create a composite mobile element.

Genetic diversity of canine olfactory receptors

BackgroundEvolution has resulted in large repertoires of olfactory receptor (OR) genes, forming the largest gene families in mammalian genomes. Knowledge of the genetic diversity of olfactory receptors is essential if we are to understand the differences in olfactory sensory capability between individuals. Canine breeds constitute an attractive model system for such investigations.ResultsWe sequenced 109 OR genes considered representative of the whole OR canine repertoire, which consists of more than 800 genes, in a cohort of 48 dogs of six different breeds. SNP frequency showed the overall level of polymorphism to be high. However, the distribution of SNP was highly heterogeneous among OR genes. More than 50% of OR genes were found to harbour a large number of SNP, whereas the rest were devoid of SNP or only slightly polymorphic. Heterogeneity was also observed across breeds, with 25% of the SNP breed-specific. Linkage disequilibrium within OR genes and OR clusters suggested a gene conversion process, consistent with a mean level of polymorphism higher than that observed for introns and intergenic sequences. A large proportion (47%) of SNP induced amino-acid changes and the Ka/Ks ratio calculated for all alleles with a complete ORF indicated a low selective constraint with respect to the high level of redundancy of the olfactory combinatory code and an ongoing pseudogenisation process, which affects dog breeds differently.ConclusionOur demonstration of a high overall level of polymorphism, likely to modify the ligand-binding capacity of receptors distributed differently within the six breeds tested, is the first step towards understanding why Labrador Retrievers and German Shepherd Dogs have a much greater potential for use as sniffer dogs than Pekingese dogs or Greyhounds. Furthermore, the heterogeneity in OR polymorphism observed raises questions as to why, in a context in which most OR genes are highly polymorphic, a subset of these genes is not? This phenomenon may be related to the nature of their ligands and their importance in everyday life.

The Impact of Recombination on Nucleotide Substitutions in the Human Genome

Unraveling the evolutionary forces responsible for variations of neutral substitution patterns among taxa or along genomes is a major issue for detecting selection within sequences. Mammalian genomes show large-scale regional variations of GC-content (the isochores), but the substitution processes at the origin of this structure are poorly understood. We analyzed the pattern of neutral substitutions in 1 Gb of primate non-coding regions. We show that the GC-content toward which sequences are evolving is strongly negatively correlated to the distance to telomeres and positively correlated to the rate of crossovers (R2 = 47%). This demonstrates that recombination has a major impact on substitution patterns in human, driving the evolution of GC-content. The evolution of GC-content correlates much more strongly with male than with female crossover rate, which rules out selectionist models for the evolution of isochores. This effect of recombination is most probably a consequence of the neutral process of biased gene conversion (BGC) occurring within recombination hotspots. We show that the predictions of this model fit very well with the observed substitution patterns in the human genome. This model notably explains the positive correlation between substitution rate and recombination rate. Theoretical calculations indicate that variations in population size or density in recombination hotspots can have a very strong impact on the evolution of base composition. Furthermore, recombination hotspots can create strong substitution hotspots. This molecular drive affects both coding and non-coding regions. We therefore conclude that along with mutation, selection and drift, BGC is one of the major factors driving genome evolution. Our results also shed light on variations in the rate of crossover relative to non-crossover events, along chromosomes and according to sex, and also on the conservation of hotspot density between human and chimp.

N-acetyltransferase 8, a positional candidate for blood pressure and renal regulation: resequencing, association and in silico study

BackgroundKidneys have an important function in blood pressure (BP) regulation and elevated BP may lead to kidney failure. Chr2p12-p13 region linked to BP traits in multiple studies harbours a potential candidate for BP and renal function, N-acetyltransferase 8 (NAT8) expressed in embryonic and adult kidney and associated with nephrotoxicity response.Methods/ResultsWe report the first study exploring NAT8 as a potential candidate gene for blood pressure and kidney function. The resequencing (n = 42, random Estonian samples) identified 15 NAT8 polymorphisms, including 6 novel variants. The diversity of NAT8 5' upstream region (π/bp = 0.00320) exceeded up to 10 times the variation in the NAT8 genic region (π/bp = 0.00037) as well as the average variation (π/bp = 0.00040) for the promoters of 29 reference genes associated with hypertension. We suggest that a potential source for such high variation could be an active gene conversion process from NAT8B duplicate gene to NAT8. Similarly to NAT8, several reference genes with the most variable upstream regions have also duplicate copies. The NAT8 promoter SNPs were targeted with pilot quantitative association studies for blood pressure (n = 137, healthy unrelated individuals) and for the index of kidney function – estimated glomerular filtration rate (eGFR; n = 157 hypertensives with and without nephropathy). Minor alleles of these polymorphisms revealed a significant protective effect against elevated systolic BP as well as kidney failure in hypertension patients (p < 0.05; linear regression model, addictive effect).ConclusionThe full resequencing and pilot association study of a novel positional candidate gene for blood pressure and renal function, human N-acetyltransferase 8, suggested a contribution of highly variable NAT8 promoter polymorphisms in determination of systolic blood pressure and eGFR. Based on in silico analysis, we raise the hypothesis that the alternative SNP alleles of the NAT8 upstream region may have differential effect on gene expression.

Analysis of Spontaneous Gene Conversion Tracts within and between Mammalian Chromosomes

In the present study, we report the first characterization of gene conversion tract length, continuity and fidelity for pathways of gene targeting, ectopic and intrachromosomal homologous recombination using the same locus and mammalian somatic cell type. In this isogenic cell system, the vast majority of recombinants (>97%) are generated by homologous recombination and display a high degree of fidelity in the gene conversion process. Individual gene conversion tracts are highly likely to involve single, independent recombination events and proceed through a heteroduplex DNA intermediate. In all recombination pathways, gene conversion tracts are long, extending up to ∼2 kb. Most gene conversion tracts are continuous in favor of donor region sequences, but in a small fraction of recombinants (15%), discontinuous gene conversion tracts are observed. In most cases, the recombination donor sequence is unaltered, although in two cases of intrachromosomal recombination, both recombination donor and recipient sequences bear gene conversion tracts. Overall, gene conversion events are similar, both qualitatively and quantitatively, for homologous recombination within and between mammalian chromosomes.

Structure and Evolution of the Atypical Mitochondrial Genome of Armadillidium vulgare (Isopoda, Crustacea)

The crustacean isopod Armadillidium vulgare is characterized by an unusual approximately 42-kb-long mitochondrial genome consisting of two molecules co-occurring in mitochondria: a circular approximately 28-kb dimer formed by two approximately 14-kb monomers fused in opposite polarities and a linear approximately 14-kb monomer. Here we determined the nucleotide sequence of the fundamental monomeric unit of A. vulgare mitochondrial genome, to gain new insight into its structure and evolution. Our results suggest that the junction zone between monomers of the dimer structure is located in or near the control region. Direct sequencing indicated that the nucleotide sequences of the different monomer units are virtually identical. This suggests that gene conversion and/or replication processes play an important role in shaping nucleotide sequence variation in this mitochondrial genome. The only heteroplasmic site we identified predicts an alloacceptor tRNA change from tRNA(Ala) to tRNA(Val). Therefore, in A. vulgare, tRNA(Ala) and tRNA(Val) are found at the same locus in different monomers, ensuring that both tRNAs are present in mitochondria. The presence of this heteroplasmic site in all sequenced individuals suggests that the polymorphism is selectively maintained, probably because of the necessity of both tRNAs for maintaining proper mitochondrial functions. Thus, our results provide empirical evidence for the tRNA gene recruitment model of tRNA evolution. Moreover, interspecific comparisons showed that the A. vulgare mitochondrial gene order is highly derived compared to the putative ancestral arthropod type. By contrast, an overall high conservation of mitochondrial gene order is observed within crustacean isopods.

The neoselectionist theory of genome evolution.

The vertebrate genome is a mosaic of GC-poor and GC-rich isochores, megabase-sized DNA regions of fairly homogeneous base composition that differ in relative amount, gene density, gene expression, replication timing, and recombination frequency. At the emergence of warm-blooded vertebrates, the gene-rich, moderately GC-rich isochores of the cold-blooded ancestors underwent a GC increase. This increase was similar in mammals and birds and was maintained during the evolution of mammalian and avian orders. Neither the GC increase nor its conservation can be accounted for by the random fixation of neutral or nearly neutral single-nucleotide changes (i.e., the vast majority of nucleotide substitutions) or by a biased gene conversion process occurring at random genome locations. Both phenomena can be explained, however, by the neoselectionist theory of genome evolution that is presented here. This theory fully accepts Ohta's nearly neutral view of point mutations but proposes in addition (i) that the AT-biased mutational input present in vertebrates pushes some DNA regions below a certain GC threshold; (ii) that these lower GC levels cause regional changes in chromatin structure that lead to deleterious effects on replication and transcription; and (iii) that the carriers of these changes undergo negative (purifying) selection, the final result being a compositional conservation of the original isochore pattern in the surviving population. Negative selection may also largely explain the GC increase accompanying the emergence of warm-blooded vertebrates. In conclusion, the neoselectionist theory not only provides a solution to the neutralist/selectionist debate but also introduces an epigenomic component in genome evolution.

Sequencing and comparative analysis of fugu protocadherin clusters reveal diversity of protocadherin genes among teleosts

BackgroundThe synaptic cell adhesion molecules, protocadherins, are a vertebrate innovation that accompanied the emergence of the neural tube and the elaborate central nervous system. In mammals, the protocadherins are encoded by three closely-linked clusters (α, β and γ) of tandem genes and are hypothesized to provide a molecular code for specifying the remarkably-diverse neural connections in the central nervous system. Like mammals, the coelacanth, a lobe-finned fish, contains a single protocadherin locus, also arranged into α, β and γ clusters. Zebrafish, however, possesses two protocadherin loci that contain more than twice the number of genes as the coelacanth, but arranged only into α and γ clusters. To gain further insight into the evolutionary history of protocadherin clusters, we have sequenced and analyzed protocadherin clusters from the compact genome of the pufferfish, Fugu rubripes.ResultsFugu contains two unlinked protocadherin loci, Pcdh1 and Pcdh2, that collectively consist of at least 77 genes. The fugu Pcdh1 locus has been subject to extensive degeneration, resulting in the complete loss of Pcdh1γ cluster. The fugu Pcdh genes have undergone lineage-specific regional gene conversion processes that have resulted in a remarkable regional sequence homogenization among paralogs in the same subcluster. Phylogenetic analyses show that most protocadherin genes are orthologous between fugu and zebrafish either individually or as paralog groups. Based on the inferred phylogenetic relationships of fugu and zebrafish genes, we have reconstructed the evolutionary history of protocadherin clusters in the teleost fish lineage.ConclusionOur results demonstrate the exceptional evolutionary dynamism of protocadherin genes in vertebrates in general, and in teleost fishes in particular. Besides the 'fish-specific' whole genome duplication, the evolution of protocadherin genes in teleost fishes is influenced by lineage-specific gene losses, tandem gene duplications and regional sequence homogenization. The dynamic protocadherin clusters might have led to the diversification of neural circuitry among teleosts, and contributed to the behavioral and physiological diversity of teleosts.

Targeted engineering of the Caenorhabditis elegans genome following Mos1-triggered chromosomal breaks

The Drosophila element Mos1 is a class II transposon, which moves by a 'cut-and-paste' mechanism and can be experimentally mobilized in the Caenorhabditis elegans germ line. Here, we triggered the excision of identified Mos1 insertions to create chromosomal breaks at given sites and further manipulate the broken loci. Double-strand break (DSB) repair could be achieved by gene conversion using a transgene containing sequences homologous to the broken chromosomal region as a repair template. Consequently, mutations engineered in the transgene could be copied to a specific locus at high frequency. This pathway was further characterized to develop an efficient tool--called MosTIC--to manipulate the C. elegans genome. Analysis of DSB repair during MosTIC experiments demonstrated that DSBs could also be sealed by end-joining in the germ line, independently from the evolutionarily conserved Ku80 and ligase IV factors. In conjunction with a publicly available Mos1 insertion library currently being generated, MosTIC will provide a general tool to customize the C. elegans genome.

Insights into Mechanisms of Bacterial Antigenic Variation Derived from the Complete Genome Sequence of Anaplasma marginale

Persistence of Anaplasma spp. in the animal reservoir host is required for efficient tick-borne transmission of these pathogens to animals and humans. Using A. marginale infection of its natural reservoir host as a model, persistent infection has been shown to reflect sequential cycles in which antigenic variants emerge, replicate, and are controlled by the immune system. Variation in the immunodominant outer-membrane protein MSP2 is generated by a process of gene conversion, in which unique hypervariable region sequences (HVRs) located in pseudogenes are recombined into a single operon-linked msp2 expression site. Although organisms expressing whole HVRs derived from pseudogenes emerge early in infection, long-term persistent infection is dependent on the generation of complex mosaics in which segments from different HVRs recombine into the expression site. The resulting combinatorial diversity generates the number of variants both predicted and shown to emerge during persistence.

Unbiased Interpretation of Haplotypes at Duplicated Microsatellites

The Y-chromosome is a powerful tool for population geneticists to study human evolutionary history. Haploid and largely non-recombining, it should contain a simple record of past mutational events. However, this apparent simplicity is compromised by Y-linked duplicons, which make up approximately 35% of this chromosome; 25% of these duplicons are large inverted repeats (palindromes). For microsatellites lying in these palindromes, two loci cannot be easily distinguished due to PCR co-amplification, and this order misspecification of alleles generates an additional variance component. Due to this ambiguity, population geneticists have traditionally used an arbitrary method to assign the alleles (shorter allele to locus 1, larger allele to locus 2). Here, we simulate these posterior estimate distributions under three different novel allele assignment priors and compare this with the original method. We use a sample of 33 human populations, typed for duplicated microsatellites lying within palindrome P8, to illustrate our approach. We show that both intra- and inter-population statistics can be dramatically affected by order misspecification. Surprisingly, matrices of pairwise F-statistics or distance estimates appear far less sensitive to order misspecification and remain relatively unchanged under the priors considered, suggesting that these microsatellites can be considered as useful markers for population genetic studies using an appropriate data treatment. Duplicated microsatellites represent an attractive source of information to investigate the extensive structural polymorphism observed among human Y chromosomes, as well as processes of intra-chromosomal gene conversion acting between duplicons.

High levels of diversity characterize mandrill (Mandrillus sphinx) Mhc-DRB sequences

The major histocompatibility complex (MHC) is highly polymorphic in most primate species studied thus far. The rhesus macaque (Macaca mulatta) has been studied extensively and the Mhc-DRB region demonstrates variability similar to humans. The extent of MHC diversity is relatively unknown for other Old World monkeys (OWM), especially among genera other than Macaca. A molecular survey of the Mhc-DRB region in mandrills (Mandrillus sphinx) revealed extensive variability, suggesting that other OWMs may also possess high levels of Mhc-DRB polymorphism. In the present study, 33 Mhc-DRB loci were identified from only 13 animals. Eleven were wild-born and presumed to be unrelated and two were captive-born twins. Two to seven different sequences were identified for each individual, suggesting that some mandrills may have as many as four Mhc-DRB loci on a single haplotype. From these sequences, representatives of at least six Mhc-DRB loci or lineages were identified. As observed in other primates, some new lineages may have arisen through the process of gene conversion. These findings indicate that mandrills have Mhc-DRB diversity not unlike rhesus macaques and humans.

Diversified bursal medullary B cells survive and expand independently after depletion following neonatal infectious bursal disease virus infection

The primary immunoglobulin repertoire of chickens is generated not by gene rearrangement but by a subsequent process of gene conversion in proliferating immature B cells within the follicles of a specialized gut-associated lymphoid organ, the bursa of Fabricius. Neonatal infection with infectious bursal disease virus can eliminate almost the entire bursal B-cell compartment. Thereafter, two types of follicle reappear. Larger follicles, with rapidly proliferating B cells and normal structure, are correlated with partial recovery of antibody response. Smaller follicles, lacking distinct cortex and medulla, appear unable to produce antigen-responsive B cells. To understand the genesis of the two types of follicle, we analysed their VL sequences and activation-induced deaminase mRNA levels. The results provide a model of bursal repopulation in which surviving bursal stem cells generate new follicles with normal morphology and function, while surviving medullary B cells continue to proliferate slowly, under the influence of stromal cells, giving rise to the smaller follicles. The latter remain fixed in a stage of development incapable of further gene diversification.

MDM2 can interact with the C-terminus of AID but it is inessential for antibody diversification in DT40 B cells

Activation-induced deaminase (AID) is essential for immunoglobulin gene diversification by the distinct processes of class switch recombination, somatic hypermutation and gene conversion. Most evidence indicates that AID triggers these reactions through the direct deamination of cytosine residues in the DNA. However, AID is predominantly cytoplasmic and the mechanism that directs it to the immunoglobulin loci remains elusive. Like its homolog APOBEC1, which requires at least one additional factor to efficiently edit APOB RNA, other proteins are likely to be required for the proper targeting of AID to the immunoglobulin loci. Here, we show that AID can interact with MDM2, an oncoprotein that shuttles between the nucleus and the cytoplasm and targets p53 for nuclear export and degradation. This interaction mapped to the carboxy-terminal region of AID that harbors a nuclear export sequence, suggesting that MDM2 may be involved in the nucleo-cytoplasmic trafficking of AID. We therefore assessed the role of MDM2 in immunoglobulin gene diversification by disrupting MDM2 in DT40, an avian B cell line that constitutively undergoes AID-dependent immunoglobulin gene diversification. The subcellular localization of AID was unaffected in MDM2-deficient DT40 cells. However, slight hyper-and hypo-conversion phenotypes were caused by MDM2-abrogation and overexpression, respectively. These observations suggested that MDM2 has the capacity to negatively regulate AID. Intriguingly, the same carboxy-terminal residues of AID were recently shown to be inessential for somatic hypermutation and immunoglobulin gene conversion but they were strictly required for class switch recombination.

Birth-and-death evolution with strong purifying selection in the histone H1 multigene family and the origin of orphon H1 genes.

Histones are small basic nuclear proteins with critical structural and functional roles in eukaryotic genomes. The H1 multigene family constitutes a very interesting histone class gathering the greatest number of isoforms, with many different arrangements in the genome, including clustered and solitary genes, and showing replication-dependent (RD) or replication-independent (RI) expression patterns. The evolution of H1 histones has been classically explained by concerted evolution through a rapid process of interlocus recombination or gene conversion. Given such intriguing features, we have analyzed the long-term evolutionary pattern of the H1 multigene family through the evaluation of the relative importance of gene conversion, point mutation, and selection in generating and maintaining the different H1 subtypes. We have found the presence of an extensive silent nucleotide divergence, both within and between species, which is always significantly greater than the nonsilent variation, indicating that purifying selection is the major factor maintaining H1 protein homogeneity. The results obtained from phylogenetic analysis reveal that different H1 subtypes are no more closely related within than between species, as they cluster by type in the topologies, and that both RD and RI H1 variants follow the same evolutionary pattern. These findings suggest that H1 histones have not been subject to any significant effect of interlocus recombination or concerted evolution. However, the diversification of the H1 isoforms seems to be enhanced primarily by mutation and selection, where genes are subject to birth-and-death evolution with strong purifying selection at the protein level. This model is able to explain not only the generation and diversification of RD H1 isoforms but also the origin and long-term persistence of orphon RI H1 subtypes in the genome, something that is still unclear, assuming concerted evolution.