Ubiquitin genes as a paradigm of concerted evolution of tandem repeats

Summary Ubiquitin is remarkable for its ubiquitous distribution and its extreme protein sequence conservation. Ubiquitin genes comprise direct repeats of the ubiquitin coding unit with no spacers. The nucleotide sequences of several ubiquitin repeats from each of humans, chicken,Xenopus, Drosophila, barley, and yeast have recently been determined. By analysis of these data we show that ubiquitin is evolving more slowly than any other known protein, and that this (together with its gene organization) contributes to an ideal situation for the occurrence of concerted evolution of tandem repeats. By contrast, there is little evidence of between-cluster concerted evolution. We deduce that in ubiquitin genes, concerted evolution involves both unequal crossover and gene conversion, and that the average time since two repeated units within the polyubiquitin locus most recently shared a common ancestor is approximately 38 million years (Myr) in mammals, but perhaps only 11 Myr inDrosophila. The extreme conservatism of ubiquitin evolution also allows the inference that certain synonymous serine codons differing at the first two positions were probably mutated at single steps.     

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.     

On the evolution of multigene families

Multigene families are classified into three groups: small families as exemplified by hemoglobin genes of mammals; middlesize multigene families, by genes of mammalian histocompatibility antigens; and large multigene families, by variable region genes of immunoglobulins. Facts and theories on these evolving multigene families are reviewed, with special reference to the population genetics of their concerted evolution. It is shown that multigene families are evolving under continued occurrence of unequal (but homologous) crossing-over and gene conversion, and that mechanisms for maintaining genetic variability are totally different from the conventional models of population genetics. Thus, in view of widespread occurrence of multigene families in genomes of higher organisms, the evolutionary theory based mainly on change of gene frequency at each locus would appear to need considerable revision.     

Kinetic basis of spontaneous mutation. Misinsertion frequencies, proofreading specificities and cost of proofreading by DNA polymerases of Escherichia coli

Repeating members of multiple-copy sequence families display high levels of sequence homogeneity. In order to examine the rates at which this is achieved, and to compare the rates with those assessed for the ribosomal DNA and histone gene families (Coen et al., 1982, accompanying paper), we have examined the patterns of variation in the Drosophila melanogaster species subgroup for the “complex” noncoding families of high copy-number. Our analysis reveals that the evolution of some of the families has involved the gradual replacement of ancestral repeats by variant repeats, independently within each species. Hybridizations between genomes at different levels of stringency indicate the presence of two basic ancestral families (the “500” and “360” families) within the subgroup. The majority of repeats representative of these families can be characterized by restriction sites and patterns of organization that are uniquely diagnostic for each species, excepting the two most closely related species. Drosophila mauritiana and Drosophila simulans. Another family (the “180” family) is confined to the one species. Drosophila orena, with features suggestive of a more rapid origin. The wide karyotypic distribution of some members of the 500 and 180 families, revealed by hybridization in situ, shows that chromosomes are evolving in concert with respect to gradual and rapidly evolving families. The distribution of sequence and pattern variation within the subgroup shows that the time required for gradual fixation (concerted evolution) of variants within large families, distributed throughout the karyotype, is longer than that required for the smaller and chromosomally restricted families of rDNA and histone genes (Coen et al., 1982). We discuss the forces that might either accelerate or retard the fixation of variants in karyotypically dispersed families.     

Transition stages of molecular drive in multiple-copy DNA families in Drosophila

Multigene and non-genic DNA families are in a state of turnover and hence are continually being replaced throughout a population by new variant repeats. To quantify such molecular processes, in the absence of selection, it is necessary to find and compare stages of transistion during the homogenization of at least two non-genic families evolving in parallel in a closely related group of species. Detailed sequence analysis of patterns of variation, at each nucleotide position considered independently, amongst repeats of two tandem DNA families from seven related Drosophila species, reveals all stages of transition during the spread of randomly produced variant repeats. Variant repeats are found at different stages of homogenization and fixation in a population, irrespective of the loci, chromosomes or individuals from which they were cloned. Differences between the families in the relatively small number of variants at each transition stage and the greater number of fully homogenized and fixed variants between species of greater divergence indicate that the process of spread (molecular drive) is rapid relative to the mutation rate and occurs at seemingly different constant rates for each family. Occasional gene conversions, in addition to unequal exchanges, have contributed to family turnover. The significance of these results to the evolution of functional multigene families and divergence and conservation of sequences is discussed.     

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.     

The Evolution of Multigene Families under Intrachromosomal Gene Conversion

A model for the evolution of the probabilities of genetic identity within and between loci of a multigene family in a finite population is formulated and investigated. Unbiased intrachromosomal gene conversion, equal crossing over between tandemly repeated genes, random genetic drift and mutation to new alleles are incorporated. Generations are discrete and nonoverlapping; the diploid, monoecious population mates at random. Formulas for the equilibrium values of the probabilities of identity and a cubic equation for the rate of convergence are deduced. Numerical examples indicate the following. The amount of homology at equilibrium generally decreases as the mutation rate, the population size and the number of repeats increase; it may increase or decrease with increasing crossover rate. The intralocus homology has an intermediate minimum, whereas the interlocus homology increases, as the rate of gene conversion increases. The intralocus homology decreases, whereas the interlocus homology increases, as the proportion of symmetric heteroduplexes increases. The characteristic convergence time can be sufficiently short to imply that intrachromosomal gene conversion may be an important mechanism for maintaining sequence homogeneity among repeated genes. The convergence time decreases as the conversion rate and the proportion of symmetric heteroduplexes increase; although exceptions occur, it generally increases as the population size and the number of repeats increase; it may increase or decrease with increasing crossover rate.     

A Concertedly Evolving Region in Chironomus, Unique Within the Telomere

Chromosome terminal, complex repeats in the dipteran Chironomus pallidivittatus show rapid concerted evolution during which there is remarkably efficient homogenization of the repeat units within and between chromosome ends. It has been shown previously that gene conversion is likely to be an important component during these changes. The sequence evolution could be a result of different processes—exchanges between repeats in the tandem array as well as information transfer between units in different chromosomes—and is therefore difficult to analyze in detail. In this study the concerted evolution of a region present only once per chromosome, at the junction between the telomeric complex repeats and the subtelomeric DNA was therefore investigated in the two sibling species C. pallidivittatus and C. tentans. Material from individual microdissected chromosome ends was used, as well as clones from bulk genomic DNA. On the telomeric side of the border pronounced species-specific sequence differences were observed, the patterns being similar for clones of different origin within each species. Mutations had been transmitted efficiently between chromosomes also when adjoining, more distally localized DNA showed great differences in sequence, suggesting that gene conversion had taken place. The evolving telomeric region bordered proximally to subtelomeric DNA with high evolutionary constancy. More proximally localized, subtelomeric DNA evolved more rapidly and showed heterogeneity between species and chromosomes. Received: 24 September 1997 / Accepted: 24 November 1997     

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.     

Concerted evolution of exons and introns in the MHC-linked tenascin-X gene of mammals.

Two modes of evolution of repeated domains in proteins have been described: (1) a conservative mode, whereby individual domains are conserved across gene duplication and speciation events, and (2) a concerted mode, whereby repeat domains become homogenized within a gene, presumably by intragenic partial duplication and/or gene conversion. The evolution of repeated EGF-like and fibronection-type-III-like (Fn-III) domains in the vertebrate extracellular matrix proteins tenascin-X (TNX) and tenascin-C (TNC) was studied by comparisons between human and mouse orthologs and between the paralogous TNC and TNX genes. The EGF-like repeats have largely been homogenized within each gene by concerted evolution since the duplication of the two genes but have been conserved since the divergence of rodents and primates. The Fn-III domains of TNC have likewise mainly evolved in a conservative fashion since the divergence of rodents and primates. In contrast, the Fn-III repeats of TNX fall into three distinct categories with regard to mode of evolution: (1) The three C-terminal repeats have been conserved since before duplication of the TNX and TNC genes. (2) Certain other repeats have been homogenized within each gene since gene duplication but have been conserved since the divergence of rodents and primates. (3) Still other repeats have evolved in a concerted fashion in rodent and primate lineages since their divergence. Remarkably, certain introns adjacent to the exons encoding these concertedly evolving Fn-III repeats have themselves evolved in a concerted fashion. This is the first known example of concerted evolution of repeated introns within a protein-coding gene.     

Adaptive evolution at the molecular level of the duplicatedAmy gene system inDrosophila

Analyses of the nucleotide sequences of the duplicatedAmy genes in the eight species of theDrosophila melanogaster species subgroup have revealed concerted evolution of the coding regions and divergent evolution between the duplicated genes of the 5’-flanking regions. Homogenization between the duplicated genes in the coding region is maintained by frequent genetic exchange in various portions of the coding region. On the other hand, such genetic exchange seems to produce a large amount of DNA sequence variation and protein polymorphism at the two loci within a species. The puzzling observation that concerted evolution is restricted to the coding regions seems to be explained by not only adaptive evolution of the AMY proteins in speciation but also adaptive fixation of selectively advantageous mutations in the intergenic region that differentiate expression of the twoAmy genes. We review molecular work on theAmy gene system inDrosophila, including evidence from biochemical characterization of the AMY proteins and molecular characterization of the cis regulatory elements.     

Complete Chloroplast Genome of Sedum sarmentosum and Chloroplast Genome Evolution in Saxifragales

Comparative chloroplast genome analyses are mostly carried out at lower taxonomic levels, such as the family and genus levels. At higher taxonomic levels, chloroplast genomes are generally used to reconstruct phylogenies. However, little attention has been paid to chloroplast genome evolution within orders. Here, we present the chloroplast genome of Sedum sarmentosum and take advantage of several available (or elucidated) chloroplast genomes to examine the evolution of chloroplast genomes in Saxifragales. The chloroplast genome of S. sarmentosum is 150,448 bp long and includes 82,212 bp of a large single-copy (LSC) region, 16.670 bp of a small single-copy (SSC) region, and a pair of 25,783 bp sequences of inverted repeats (IRs).The genome contains 131 unique genes, 18 of which are duplicated within the IRs. Based on a comparative analysis of chloroplast genomes from four representative Saxifragales families, we observed two gene losses and two pseudogenes in Paeonia obovata, and the loss of an intron was detected in the rps16 gene of Penthorum chinense. Comparisons among the 72 common protein-coding genes confirmed that the chloroplast genomes of S. sarmentosum and Paeonia obovata exhibit accelerated sequence evolution. Furthermore, a strong correlation was observed between the rates of genome evolution and genome size. The detected genome size variations are predominantly caused by the length of intergenic spacers, rather than losses of genes and introns, gene pseudogenization or IR expansion or contraction. The genome sizes of these species are negatively correlated with nucleotide substitution rates. Species with shorter duration of the life cycle tend to exhibit shorter chloroplast genomes than those with longer life cycles.     

Population genetics theory of concerted evolution and its application to the immunoglobulin V gene tree

Summary The previous simple model for treating concerted evolution of multigene families has been revised to be compatible with various new observations on the immunoglobulin variable region family and other families. In the previous model, gene conversion and unequal crossing-over were considered, and it was assumed that genes are randomly arranged on the chromosome; neither subdivision nor correlation of gene identity and chromosomal distance were considered. Although this model satisfactorily explains the observed amino acid diversity within and between species, it fails to predict the very ancient branching of the mouse immunoglobulin heavy chain V-gene family. By incorporating subdivided structure and genetic correlation with chromosomal distance into the simple model, the data of divergence may be satisfactorily explained, as well as the rate of nucleotide substitution and the amino acid diversity. The rate at which a V-gene is duplicated or deleted by conversion or by unequal crossing-over is estimated by the new model to be on the order of 10^–6 per year. The model may be applicable to other multigene families, such as those coding for silkmoth chorion or mammalian kallikrein.Contribution no. 1560 from the National Institute of Genetics, Mishima, 411 Japan     

Concerted Evolution in the Ribosomal RNA Cistron

Background

Gene conversion is the mechanism proposed to be responsible for the homogenization of multigene families such as the nuclear ribosomal gene clusters. This concerted evolutionary process prevents individual genes in gene clusters from accumulating mutations. The mechanism responsible for concerted evolution is not well understood but recombination during meiosis has been hypothesized to play a significant role in this homogenization. In this study we tested the hypothesis of unequal crossing over playing a significant role in gene conversion events within the ribosomal RNA cistron during meiosis, mitosis or both life stages in the fungal tree pathogen Ceratocystis manginecans.

Methods

Ceratocystis manginecans, a haploid ascomycete, reproduces homothallically and was found to have two distinct sequences within the internally transcribed spacer (ITS) region of the ribosomal RNA cistron. The different ITS types were scored using PCR-RFLP assays and chi-square analyses to determine the level of significance of the changes in the ratios of the ITS types.

Results

The relative ratios of the two ITS sequence types changed when the fungal isolates were cultured vegetatively or allowed to produced sexual structures and spores. These active changes were shown to occur more frequently during meiosis than mitosis.

Conclusion

The evidence presented provides concrete support for homogenization in the rRNA gene clusters found in this fungus and that the most reasonable explanation for this process is unequal crossing over.     

Complete Chloroplast Genome Sequence of Omani Lime (Citrus aurantiifolia) and Comparative Analysis within the Rosids

The genus Citrus contains many economically important fruits that are grown worldwide for their high nutritional and medicinal value. Due to frequent hybridizations among species and cultivars, the exact number of natural species and the taxonomic relationships within this genus are unclear. To compare the differences between the Citrus chloroplast genomes and to develop useful genetic markers, we used a reference-assisted approach to assemble the complete chloroplast genome of Omani lime (C. aurantiifolia). The complete C. aurantiifolia chloroplast genome is 159,893 bp in length; the organization and gene content are similar to most of the rosids lineages characterized to date. Through comparison with the sweet orange (C. sinensis) chloroplast genome, we identified three intergenic regions and 94 simple sequence repeats (SSRs) that are potentially informative markers with resolution for interspecific relationships. These markers can be utilized to better understand the origin of cultivated Citrus. A comparison among 72 species belonging to 10 families of representative rosids lineages also provides new insights into their chloroplast genome evolution.     

How gene families evolve

Theories and facts of gene family evolution are reviewed. Concerted evolution is commonly observed for gene families which originated a long time ago, however there are many different types of multigene families, from uniform to diverse. The rate of homogenization by unequal crossing-over, gene conversion, etc. has been evolutionarily adjusted for each gene family. When new functions are needed by organisms, gene families may evolve into superfamilies, in which no further concerted evolution takes place, and each member of the family may acquire an indispensable function. The homeobox-containing gene family is a most exciting example of such superfamily.     

The Distribution of GYR- and YLP-Like Motifs in Drosophila Suggests a General Role in Cuticle Assembly and Other Protein-Protein Interactions

Background

Arthropod cuticle is composed predominantly of a self-assembling matrix of chitin and protein. Genes encoding structural cuticular proteins are remarkably abundant in arthropod genomes, yet there has been no systematic survey of conserved motifs across cuticular protein families.

Methodology/Principal Findings

Two short sequence motifs with conserved tyrosines were identified in Drosophila cuticular proteins that were similar to the GYR and YLP Interpro domains. These motifs were found in members of the CPR, Tweedle, CPF/CPFL, and (in Anopheles gambiae) CPLCG cuticular protein families, and the Dusky/Miniature family of cuticle-associated proteins. Tweedle proteins have a characteristic motif architecture that is shared with the Drosophila protein GCR1 and its orthologs in other species, suggesting that GCR1 is also cuticular. A resilin repeat, which has been shown to confer elasticity, matched one of the motifs; a number of other Drosophila proteins of unknown function exhibit a motif architecture similar to that of resilin. The motifs were also present in some proteins of the peritrophic matrix and the eggshell, suggesting molecular convergence among distinct extracellular matrices. More surprisingly, gene regulation, development, and proteolysis were statistically over-represented ontology terms for all non-cuticular matches in Drosophila. Searches against other arthropod genomes indicate that the motifs are taxonomically widespread.

Conclusions

This survey suggests a more general definition for GYR and YLP motifs and reveals their contribution to several types of extracellular matrix. They may define sites of protein interaction with DNA or other proteins, based on ontology analysis. These results can help guide experimental studies on the biochemistry of cuticle assembly.     

Gene conversions within the skeletal myosin multigene family.

Comparisons of the nucleotide sequences of the light meromyosin (LMM) region of developmentally regulated fast chicken myosin heavy chain (MHC) isoforms indicates that chicken MHC isoforms are more similar to each other than to MHC isoforms in other species. The sequence data provide evidence that gene conversion events have occurred recently among the isoforms. An embryonic (Cemb1) isoform and neonatal isoform have the most extensive regions of sequence identity. Similar gene conversion events are present in the rat alpha- and beta-cardiac MHCs, but were not obvious in the LMM of developmentally regulated fast human MHC isoforms. The data suggest that gene conversion events can play a significant role in the evolution of the MHC multigene families and that concerted evolution of the chicken multigene family occurred after the divergence of mammals and avians.     

The Chloroplast Genome of Hyoscyamus niger and a Phylogenetic Study of the Tribe Hyoscyameae (Solanaceae)

The tribe Hyoscyameae (Solanaceae) is restricted to Eurasia and includes the genera Archihyoscyamus, Anisodus, Atropa, Atropanthe, Hyoscyamus, Physochlaina, Przewalskia and Scopolia. Even though the monophyly of Hyoscyameae is strongly supported, the relationships of the taxa within the tribe remain unclear. Chloroplast markers have been widely used to elucidate plant relationships at low taxonomic levels. Identification of variable chloroplast intergenic regions has been developed based on comparative genomics of chloroplast genomes, but these regions have a narrow phylogenetic utility. In this study, we present the chloroplast genome sequence of Hyoscyamus niger and make comparisons to other solanaceous plastid genomes in terms of gene order, gene and intron content, editing sites, origins of replication, repeats, and hypothetical open reading frames. We developed and sequenced three variable plastid markers from eight species to elucidate relationships within the tribe Hyoscyameae. The presence of a horizontally transferred intron in the mitochondrial cox1 gene of some species of the tribe is considered here a likely synapomorphy uniting five genera of the Hyoscyameae. Alternatively, the cox1 intron could be a homoplasious character acquired twice within the tribe. A homoplasious inversion in the intergenic plastid spacer trnC-psbM was recognized as a source of bias and removed from the data set used in the phylogenetic analyses. Almost 12 kb of plastid sequence data were not sufficient to completely resolve relationships among genera of Hyoscyameae but some clades were identified. Two alternative hypotheses of the evolution of the genera within the tribe are proposed.