Molecular evolution of the nontandemly repeated genes of the histone 3 multigene family.

In some species, histone gene clusters consist of tandem arrays of each type of histone gene, whereas in other species the genes may be clustered but not arranged in tandem. In certain species, however, histone genes are found scattered across several different chromosomes. This study examines the evolution of histone 3 (H3) genes that are not arranged in large clusters of tandem repeats. Although H3 amino acid sequences are highly conserved both within and between species, we found that the nucleotide sequence divergence at synonymous sites is high, indicating that purifying selection is the major force for maintaining H3 amino acid sequence homogeneity over long-term evolution. In cases where synonymous-site divergence was low, recent gene duplication appeared to be a better explanation than gene conversion. These results, and other observations on gene inactivation, organization, and phylogeny, indicated that these H3 genes evolve according to a birth-and-death process under strong purifying selection. Thus, we found little evidence to support previous claims that all H3 proteins, regardless of their genome organization, undergo concerted evolution. Further analyses of the structure of H3 proteins revealed that the histones of higher eukaryotes might have evolved from a replication-independent-like H3 gene.     

Early Evolution of Histone Genes: Prevalence of an ‘Orphon’ H1 Lineage in Protostomes and Birth-and-Death Process in the H2A Family

The study of histone evolution has experienced a rebirth, for two main reasons: the identification of new essential histone variants responsible for regulating chromatin dynamics and the subsequent contradictions posed by this variability as it pertains to their long-term evolution process. Although different evolutionary models (e.g., birth-and-death evolution, concerted evolution) may account for the observed divergence of histone genes, conclusive evidence is lacking (e.g., histone H1) or totally nonexistent (e.g., histone H2A). While most of the published work has focused on deuterostomes, very little is known about the diversification and functional differentiation mechanisms followed by histone protein subtypes in protostomes, for which histone variants have only been recently described. In this study, we identify linker and core histone genes in three clam species. Our results demonstrate the prevalence of an 'orphon' H1 lineage in molluscs, a group in which the protostome H1 and sperm nuclear basic proteins are on the verge of diversification. They share an early monophyletic origin with vertebrate-specific variants prior to the differentiation between protostomes and deuterostomes. Given the intringuing evolutionary features of the histone H1 family, we have evaluated the relative importance of gene conversion, point mutation, and selection in maintaining the diversity found among H2A subtypes in eukaryotes. We show evidence for the first time that the long-term evolution of this family is not subject to concerted evolution but, rather, to a gradual evolution following a birth-and-death model under a strong purifying selection at the protein level.     

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.     

Identifying concerted evolution and gene conversion in mammalian gene pairs lasting over 100 million years

Background  

Concerted evolution occurs in multigene families and is characterized by stretches of homogeneity and higher sequence similarity between paralogues than between orthologues. Here we identify human gene pairs that have undergone concerted evolution, caused by ongoing gene conversion, since at least the human-mouse divergence. Our strategy involved the identification of duplicated genes with greater similarity within a species than between species. These genes were required to be present in multiple mammalian genomes, suggesting duplication early in mammalian divergence. To eliminate genes that have been conserved due to strong purifying selection, our analysis also required at least one intron to have retained high sequence similarity between paralogues.     

Nucleotide sequences of Caenorhabditis elegans core histone genes. Genes for different histone classes share common flanking sequence elements

We have determined the nucleotide sequence of core histone genes and flanking regions from two of approximately 11 different genomic histone clusters of the nematode Caenorhabditis elegans. Four histone genes from one cluster (H3, H4, H2B, H2A) and two histone genes from another (H4 and H2A) were analyzed. The predicted amino acid sequences of the two H4 and H2A proteins from the two clusters are identical, whereas the nucleotide sequences of the genes have diverged 9% (H2A) and 12% (H4). Flanking sequences, which are mostly not similar, were compared to identify putative regulatory elements. A conserved sequence of 34 base-pairs is present 19 to 42 nucleotides 3' of the termination codon of all the genes. Within the conserved sequence is a 16-base dyad sequence homologous to the one typically found at the 3' end of histone genes from higher eukaryotes. The C. elegans core histone genes are organized as divergently transcribed pairs of H3-H4 and H2A-H2B and contain 5' conserved sequence elements in the shared spacer regions. One of the sequence elements, 5' CTCCNCCTNCCCACCNCANA 3', is located immediately upstream from the canonical TATA homology of each gene. Another sequence element, 5' CTGCGGGGACACATNT 3', is present in the spacer of each heterotypic pair. These two 5' conserved sequences are not present in the promoter region of histone genes from other organisms, where 5' conserved sequences are usually different for each histone class. They are also not found in non-histone genes of C. elegans. These putative regulatory sequences of C. elegans core histone genes are similar to the regulatory elements of both higher and lower eukaryotes. The coding regions of the genes and the 3' regulatory sequences are similar to those of higher eukaryotes, whereas the presence of common 5' sequence elements upstream from genes of different histone classes is similar to histone promoter elements in yeast.     

Nucleotide Variation and Divergence in the Histone Multigene Family in Drosophila Melanogaster

Nucleotide differences in the histone H3 gene family in Drosophila melanogaster were studied on three levels: (1) within a chromosome, (2) within a population and (3) between species (D. melanogaster and Drosophila simulans). The average difference within the H3 gene within a chromosome was 0.0040 per nucleotide site, about 52% of that within a population (0.0077). The proportion of divergent sites between the two species was 0.0575, which is about 8.5 times the difference within a species. The distribution of divergence between species was similar to that of variation within a species. Divergence and variation were noted to be greatest in the 3' noncoding region and least in the coding region. Values intermediate between these were found for the 5' noncoding region. Divergence and variation in silent sites exceeded those in the total coding region, thus indicating possible purifying selection for amino-acid-altering change. Phylogenetic relations among H3 genes and genetic differences on these three levels are evidence for the concerted evolution of the histone gene family. The molecular mechanism by which variation is produced and maintained is discussed.     

Histone genes inPhysarum polycephalum: Transcription and analysis of the flanking regions of the two H4 genes

The histone H4 multigene family of Physarum polycephalum consists of two genes, H41 and H42. Both genes have an unusual structure in that they are interrupted by a small intron. The structure of the P. polycephalum H4 genes is discussed and compared to the structure of histone genes of other organisms. S1 nuclease analysis was used to map the 5' and 3' ends of the histone H4 messengers. We show that the histone H4 genes have a hybrid structure; they are interrupted by an intervening sequence, as in replacement variant histone genes of higher eukaryotes, but their 5' and 3' noncoding regions have the properties of replication-dependent histone genes: the 5' and 3' leader and trailer sequences are short, possess a 3'-hyphenated dyad symmetry element, and a CAGA sequence is found 3' to the hyphenated hairpin structure. This report also provides evidence that both genes are expressed in late G2 phase as well as in S phase and that their expression is temporally coordinated and quantitatively similar during the cell cycle.     

An unusual evolutionary behaviour of a sea urchin histone gene cluster

DNA sequences of cloned histone coding sequences and spacers of sea urchin species that diverged long ago in evolution were compared. The highly repeated H4 and H3 genes active during early embryogenesis had evolved (in their silent sites) at a rate (0.5-0.6% base changes/Myr) similar to single-copy protein-coding genes and nearly as fast as spacer DNA (0.7% base changes/Myr) and unique DNA. Thus, evolution in the major histone genes conforms to a universal evolutionary clock based on the rate of base sequence change. By contrast, the H4 and H3 coding sequences and a non-transcribed spacer of the DNA clone h19 of Psammechinus miliaris show an exceptionally low rate of sequence evolution only 1/100 to 1/200 that predicted from the clock hypothesis. According to the classical model of gene inheritance, the h19 DNA sequences in the Psammechinus genome require unusual conservation mechanisms by selection at the level of the gene and spacer sequences. An alternative explanation could be recent horizontal gene transfer of a histone gene cluster from the very distantly related Strongylocentrotus dröbachiensis to the P. miliaris genome.     

Generalized linear mixed model for segregation distortion analysis

Background

Concerted evolution refers to the pattern in which copies of multigene families show high intraspecific sequence homogeneity but high interspecific sequence diversity. Sequence homogeneity of these copies depends on relative rates of mutation and recombination, including gene conversion and unequal crossing over, between misaligned copies. The internally repetitive intergenic spacer (IGS) is located between the genes for the 28S and 18S ribosomal RNAs. To identify patterns of recombination and/or homogenization within IGS repeat arrays, and to identify regions of the IGS that are under functional constraint, we analyzed 13 complete IGS sequences from 10 individuals representing four species in the Daphnia pulex complex.

Results

Gene conversion and unequal crossing over between misaligned IGS repeats generates variation in copy number between arrays, as has been observed in previous studies. Moreover, terminal repeats are rarely involved in these events. Despite the occurrence of recombination, orthologous repeats in different species are more similar to one another than are paralogous repeats within species that diverged less than 4 million years ago. Patterns consistent with concerted evolution of these repeats were observed between species that diverged 8-10 million years ago. Sequence homogeneity varies along the IGS; the most homogeneous regions are downstream of the 28S rRNA gene and in the region containing the core promoter. The inadvertent inclusion of interspecific hybrids in our analysis uncovered evidence of both inter- and intrachromosomal recombination in the nonrepetitive regions of the IGS.

Conclusions

Our analysis of variation in ribosomal IGS from Daphnia shows that levels of homogeneity within and between species result from the interaction between rates of recombination and selective constraint. Consequently, different regions of the IGS are on substantially different evolutionary trajectories.     

Concerted evolution of sea anemone neurotoxin genes is revealed through analysis of the Nematostella vectensis genome

Gene families, which encode toxins, are found in many poisonous animals, yet there is limited understanding of their evolution at the nucleotide level. The release of the genome draft sequence for the sea anemone Nematostella vectensis enabled a comprehensive study of a gene family whose neurotoxin products affect voltage-gated sodium channels. All gene family members are clustered in a highly repetitive approximately 30-kb genomic region and encode a single toxin, Nv1. These genes exhibit extreme conservation at the nucleotide level which cannot be explained by purifying selection. This conservation greatly differs from the toxin gene families of other animals (e.g., snakes, scorpions, and cone snails), whose evolution was driven by diversifying selection, thereby generating a high degree of genetic diversity. The low nucleotide diversity at the Nv1 genes is reminiscent of that reported for DNA encoding ribosomal RNA (rDNA) and 2 hsp70 genes from Drosophila, which have evolved via concerted evolution. This evolutionary pattern was experimentally demonstrated in yeast rDNA and was shown to involve unequal crossing-over. Through sequence analysis of toxin genes from multiple N. vectensis populations and 2 other anemone species, Anemonia viridis and Actinia equina, we observed that the toxin genes for each sea anemone species are more similar to one another than to those of other species, suggesting they evolved by manner of concerted evolution. Furthermore, in 2 of the species (A. viridis and A. equina) we found genes that evolved under diversifying selection, suggesting that concerted evolution and accelerated evolution may occur simultaneously.     

Concerted evolution of the immunoglobulin VH gene family

With the aim of understanding the concerted evolution of the immunoglobulin VH multigene family, a phylogenetic tree for the DNA sequences of 16 mouse and five human germ line genes was constructed. This tree indicates that all genes in this family have undergone substantial evolutionary divergence. The most closely related genes so far identified in the mouse genome seem to have diverged about 6 million years (MY) ago, whereas the most distantly related genes diverged about 300 MY ago. This suggests that gene duplication caused by unequal crossing-over or gene conversion occurs very slowly in this gene family. The rate of occurrence of gene duplication in the VH gene family has been estimated to be 5 x 10(-7) per gene per year, which seems to be at least about 100 times lower than that for the rRNA gene family. This low rate of concerted evolution in the VH gene family helps retain intergenic genetic variability that in turn contributes to antibody diversity. Because of accumulation of destructive mutations, however, about one-third of the mouse and human VH genes seem to have become nonfunctional. Many of these pseudogenes have apparently originated recently, but some of them seem to have existed in the genome for more than 10 MY. The rate of nucleotide substitution for the complementarity-determining regions (CDRs) is as high as that of pseudogenes. This suggests that there is virtually no purifying selection operating in the CDRs and that germ line mutations are effectively used for generating antibody diversity.      

Either of the major H2A genes but not an evolutionarily conserved H2A.F/Z variant of Tetrahymena thermophila can function as the sole H2A gene in the yeast Saccharomyces cerevisiae.

H2A.F/Z histones are conserved variants that diverged from major H2A proteins early in evolution, suggesting they perform an important function distinct from major H2A proteins. Antisera specific for hv1, the H2A.F/Z variant of the ciliated protozoan Tetrahymena thermophila, cross-react with proteins from Saccharomyces cerevisiae. However, no H2A.F/Z variant has been reported in this budding yeast species. We sought to distinguish among three explanations for these observations: (i) that S. cerevisiae has an undiscovered H2A.F/Z variant, (ii) that the major S. cerevisiae H2A proteins are functionally equivalent to H2A.F/Z variants, or (iii) that the conserved epitope is found on a non-H2A molecule. Repeated attempts to clone an S. cerevisiae hv1 homolog only resulted in the cloning of the known H2A genes yHTA1 and yHTA2. To test for functional relatedness, we attempted to rescue strains lacking the yeast H2A genes with either the Tetrahymena major H2A genes (tHTA1 or tHTA2) or the gene (tHTA3) encoding hv1. Although they differ considerably in sequence from the yeast H2A genes, the major Tetrahymena H2A genes can provide the essential functions of H2A in yeast cells, the first such case of trans-species complementation of histone function. The Tetrahymena H2A genes confer a cold-sensitive phenotype. Although expressed at high levels and transported to the nucleus, hv1 cannot replace yeast H2A proteins. Proteins from S. cerevisiae strains lacking yeast H2A genes fail to cross-react with anti-hv1 antibodies. These studies make it likely that S. cerevisiae differs from most other eukaryotes in that it does not have an H2A.F/Z homolog. A hypothesis is presented relating the absence of H2A.F/Z in S. cerevisiae to its function in other organisms.     

Evolution of histone H3: emergence of variants and conservation of post-translational modification sites

Histone H3 proteins are highly conserved across all eukaryotes and are dynamically modified by many post-translational modifications (PTMs). Here we describe a method that defines the evolution of the family of histone H3 proteins, including the emergence of functionally distinct variants. It combines information from histone H3 protein sequences in eukaryotic species with the evolution of these species as described by the tree of life (TOL) project. This so-called TOL analysis identified the time when the few observed protein sequence changes occurred and when distinct, co-existing H3 protein variants arose. Four distinct ancient duplication events were identified where replication-coupled (RC) H3 variants diverged from replication-independent (RI) forms, like histone H3.3 in animals. These independent events occurred in ancestral lineages leading to the clades of metazoa, viridiplantae, basidiomycota, and alveolata. The proto-H3 sequence in the last eukaryotic common ancestor (LECA) was expanded to at least 133 of its 135 residues. Extreme conservation of known acetylation and methylation sites of lysines and arginines predicts that these PTMs will exist across the eukaryotic crown phyla and in protists with canonical chromatin structures. Less complete conservation was found for most serine and threonine phosphorylation sites. This study demonstrates that TOL analysis can determine the evolution of slowly evolving proteins in sequence-saturated datasets.     

Low codon bias and high rates of synonymous substitution in Drosophila hydei and D. melanogaster histone genes

We have evaluated codon usage bias in Drosophila histone genes and have obtained the nucleotide sequence of a 5,161-bp D. hydei histone gene repeat unit. This repeat contains genes for all five histone proteins (H1, H2a, H2b, H3, and H4) and differs from the previously reported one by a second EcoRI site. These D. hydei repeats have been aligned to each other and to the 5.0-kb (i.e., long) and 4.8-kb (i.e., short) histone repeat types from D. melanogaster. In each species, base composition at synonymous sites is similar to the average genomic composition and approaches that in the small intergenic spacers of the histone gene repeats. Accumulation of synonymous changes at synonymous sites after the species diverged is quite high. Both of these features are consistent with the relatively low codon usage bias observed in these genes when compared with other Drosophila genes. Thus, the generalization that abundantly expressed genes in Drosophila have high codon bias and low rates of silent substitution does not hold for the histone genes.      

Dramatic diversity of ciliate histone H4 genes revealed by comparisons of patterns of substitutions and paralog divergences among eukaryotes

The accumulation of divergent histone H4 amino acid sequences within and between ciliate lineages challenges traditional views of the evolution of this essential eukaryotic protein. We analyzed histone H4 sequences from 13 species of ciliates and compared these data with sequences from well-sampled eukaryotic clades. Ciliate histone H4s differ from one another at as many as 46% of their amino acids, in contrast with the highly conserved character of this protein in most other eukaryotes. Equally striking, we find paralogs of histone H4 within ciliate genomes that differ by up to 25% of their amino acids, whereas paralogs in other eukaryotes share identical or nearly identical amino acid sequences. Moreover, the most divergent H4 proteins within ciliates are found in the lineages with highly processed macronuclear genomes. Our analyses demonstrate that the dual nature of ciliate genomes-the presence of a "germline" micronucleus and a "somatic" macronucleus within each cell-allowed the dramatic variation in ciliate histone genes by altering functional constraints or enabling adaptive evolution of the histone H4 protein, or both.