The evolution of hexapod engrailed-family genes: evidence for conservation and concerted evolution

Phylogenetic analyses imply that multiple engrailed-family gene duplications occurred during hexapod evolution, a view supported by previous reports of only a single engrailed-family gene in members of the grasshopper genus Schistocerca and in the beetle Tribolium castaneum. Here, we report the cloning of a second engrailed-family gene from Schistocerca gregaria and present evidence for two engrailed-family genes from four additional hexapod species. We also report the existence of a second engrailed-family gene in the Tribolium genome. We suggest that the engrailed and invected genes of Drosophila melanogaster have existed as a conserved gene cassette throughout holometabolous insect evolution. In total 11 phylogenetically diverse hexapod orders are now known to contain species that possess two engrailed-family paralogues, with in each case only one paralogue encoding the RS-motif, a characteristic feature of holometabolous insect invected proteins. We propose that the homeoboxes of hexapod engrailed-family paralogues are evolving in a concerted fashion, resulting in gene trees that overestimate the frequency of gene duplication. We present new phylogenetic analyses using non-homeodomain amino acid sequence that support this view. The S. gregaria engrailed-family paralogues provide strong evidence that concerted evolution might in part be explained by recurrent gene conversion. Finally, we hypothesize that the RS-motif is part of a serine-rich domain targeted for phosphorylation.     

Duplication,concerted evolution and purifying selection drive the evolution of mosquito vitellogenin genes

Background  

Mosquito vitellogenin (Vtg) genes belong to a small multiple gene family that encodes the major yolk protein precursors required for egg production. Multiple Vtg genes have been cloned and characterized from several mosquito species, but their origin and molecular evolution are poorly understood.     

The evolutionary rate of duplicated genes under concerted evolution

The effect of directional selection on the fixation process of a single mutation that spreads in a multigene family by gene conversion is investigated. A simple two-locus model with two alleles, A and a, is first considered in a random-mating diploid population with size N. There are four haplotypes, AA, Aa, aA, and aa, and selection works on the number of alleles A in a diplod (i = 0, 1, 2, 3, 4). Because gene conversion is allowed between the two loci, when the mutation rate is very low, either AA or aa will fix in the population eventually. We consider a situation where a single mutant, A, arises in one locus when a is fixed in both loci. Then, we derive the fixation probability analytically, and the fixation time is investigated by simulations. It is found that gene conversion has an effect to increase the "effective" population size, so that weak selection works more efficiently in a multigene family. With these results, we discuss the effect of gene conversion on the rate of molecular evolution in a multigene family undergoing concerted evolution. We also argue about the applicability of the theoretical results to models of multigene families with more than two loci.     

Organization and concerted evolution of the ampliconic Y-chromosomal TSPY genes from cattle

The Y-chromosomal gene TSPY (testis-specific protein Y-encoded) is probably involved in early spermatogenesis and has a variable copy number in different mammalian species. Analysis of bovine BAC clones leads to an estimate of 90 TSPY loci on the bovine Y chromosome. Half of these loci (TSPY-M1 and TSPY-M2) contain a single copy, while the other loci (TSPY-C) contain a cluster of three, possibly four, truncated pseudogenes. Fluorescence in situ hybridization indicated that the TSPY loci are located mainly on the short arm (Yp). The TSPY genes appear to account for about 2.5% of the Y chromosome and contain several published bovine Y-chromosomal microsatellites. The homology of TSPY and the major Y-chromosomal repetitive elements BRY.2 from cattle and OY.1 from sheep (80-85% similarity) further illustrates how the Y chromosome is shaped by rearrangements and horizontal spreading of the most abundant sequences. A comparison of TSPY-M1 sequences from different BAC clones and from related bovine species suggests concerted evolution as one of the mechanisms of the rapid evolution of the mammalian Y chromosome.     

Characterization of the cow stomach lysozyme genes: Repetitive DNA and concerted evolution

Cow stomach lysozyme genes have evolved in a mosaic pattern. The majority of the intronic and flanking sequences show an amount of sequence difference consistent with divergent evolution since duplication of the genes 40–50 million years ago. In contrast, exons 1, 2, and 4 and immediately adjacent intronic sequences differ little between genes and show evidence of recent concerted evolution. Exon 3 appears to be evolving divergently. The three characterized genes vary from 5.6 to 7.9 kilobases in length. Different distributions of repetitive DNA are found in each gene, which accounts for the majority of length differences between genes. The different distributions of repetitive DNA in each gene suggest the repetitive elements were inserted into each gene after the duplications that give rise to these three genes and provide additional support for divergent evolution for the majority of each gene. The observation that intronic and flanking sequences are evolving divergently suggests that the concerted evolution events involved in homogenizing the coding regions of lysozyme genes involve only one exon at a time. This model of concerted evolution would allow the shuffling of exon-sized pieces of information between genes, a phenomenon that may have aided in the early adaptive evolution of stomach lysozyme.Deceased July 21, 1991 Correspondence to: D.M. Irwin     

The evolution of concerted evolution

Concerted evolution is a consequence of processes that convert copies of a gene in a multigene family into the same copy. Here we ask whether this homogenization may be adaptive. Analysis of a modifier of homogenization reveals (1) that the trait is most likely to spread if interactions between deleterious mutations are not strongly synergistic; (2) that selection on the modifier is of the order of the mutation rate, hence the modifier is most likely to be favoured by selection when the species has a large effective population size and/or if the modifier affects many genes simultaneously; and (3) that linkage between the genes in the family, and between these genes and the modifier, makes invasion of the modifier easier, suggesting that selection may favour multigene families being in clustered arrays. It follows from the first conclusion that genes for which mutations may often be dominant or semi-dominant should undergo concerted evolution more commonly than others. By analysis of the mouse knockout database, we show that mutations affecting growth-related genes are more commonly associated with dominant lethality than expected by chance. We predict then that selection will favour homogenization of such genes, and possibly others that are significantly dosage dependent, more often than it favours homogenization in other genes. The first condition is almost the opposite of that required for the maintenance of sexual reproduction according to the mutation-deterministic theory. The analysis here therefore suggests that sexual organisms can simultaneously minimize both the effects of deleterious, strongly synergistically, interacting mutations and those that interact either weakly synergistically, multiplicatively, or antagonistically, assuming the latter class belong to a multicopy gene family. Recombination and an absence of homogenization are efficient in purging deleterious mutations in the former class, homogenization and an absence of recombination are efficient at minimizing the costs imposed by the latter classes.     

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.     

Functional conservation of the yeast and Arabidopsis RAD54-like genes

The Saccharomyces cerevisiae RAD54 gene has critical roles in DNA double-strand break repair, homologous recombination, and gene targeting. Previous results show that the yeast gene enhances gene targeting when expressed in Arabidopsis thaliana. In this work we address the trans-species compatibility of Rad54 functions. We show that overexpression of yeast RAD54 in Arabidopsis enhances DNA damage resistance severalfold. Thus, the yeast gene is active in the Arabidopsis homologous-recombination repair system. Moreover, we have identified an A. thaliana ortholog of yeast RAD54, named AtRAD54. This gene, with close sequence similarity to RAD54, complements methylmethane sulfonate (MMS) sensitivity but not UV sensitivity or gene targeting defects of rad54Delta mutant yeast cells. Overexpression of AtRAD54 in Arabidopsis leads to enhanced resistance to DNA damage. This gene's assignment as a RAD54 ortholog is further supported by the interaction of AtRad54 with AtRad51 and the interactions between alien proteins (i.e., yeast Rad54 with AtRAD51 and yeast Rad51 with AtRad54) in a yeast two-hybrid experiment. These interactions hint at the molecular nature of this interkingdom complementation, although the stronger effect of the yeast Rad54 in plants than AtRad54 in yeast might be explained by an ability of the Rad54 protein to act alone, independently of its interaction with Rad51.     

Estimating the time to the whole-genome duplication and the duration of concerted evolution via gene conversion in yeast

A maximum-likelihood (ML) method is developed to estimate the duration of concerted evolution and the time to the whole-genome duplication (WGD) event in baker's yeast (Saccharomyces cerevisiae). The models with concerted evolution fit the data significantly better than the molecular clock model, indicating a crucial role of concerted evolution via gene conversion after gene duplication in yeast. Our ML estimate of the time to the WGD is nearly identical to the time to the speciation event between S. cerevisiae and Kluyveromyces waltii, suggesting that the WGD occurred in very early stages after speciation or the WGD might have been involved in the speciation event.     

Molecular characterization and concerted evolution of two genes encoding RING-C2 type proteins in rice

RING (Really Interesting New Gene) finger proteins are believed to play a critical role in mediating the transfer of ubiquitin to heterogeneous substrate(s). While the two canonical types, RING-H2 and RING-HC, have been well-characterized, the molecular functions of the modified types, particularly the RING-C2 types, remain elusive. We isolated two rice genes harboring the RING-C2 domain on the distal parts of rice chromosomes 11 and 12, termed OsRINGC2-1 and OsRINGC2-2, respectively. A comparison of sequence divergences between 10 duplicate pairs on the distal parts of rice chromosomes 11 and 12 and randomly selected duplicate pairs suggested that OsRINGC2-1 and OsRINGC2-2 have evolved in concert via gene conversion. An in vitro ubiquitination assay revealed that both proteins possess E3 ligase activity, suggesting that the innate functions of these RING domains have not been affected by their modifications during evolution. Subcellular localizations were strikingly different; OsRINGC2-1 was found only in the cytoplasm with many punctate complexes, whereas OsRINGC2-2 was observed in both the nucleus and cytoplasm. The expression patterns of both genes showed striking differences in response to salt stress, whereas plants heterogeneous for both genes mediated salt tolerance in Arabidopsis, supporting the notion of concerted evolution. These results shed light on the molecular functions of OsRINGC2-1 and OsRINGC2-2 and provide insight into their molecular evolution.     

A comparative study of vicilin genes in Lens: negative evidence of concerted evolution

Genes for vicilin, a component of legume seed storage proteins, have been identified in the cultivated lentil (Lens culinaris ssp. culinaris) and in wild species of the genus Lens. Five different types of vicilin sequences (designated A-E) have been identified in each lentil individual. The different types of sequences, and some possible variants of them (also present in each individual) are part of the vicilin family of genes. Type D sequences have the characteristics of nonprocessed pseudogenes. Comparison of nucleotide sequences indicates that lentil vicilin sequences are similar to vicilin sequences of other legume species, in particular to those of the tribe Vicieae, in which the genes Lens is included. Sequence comparison and distance and parsimony trees indicated that two groups or subfamilies of sequences, including, respectively, types A, B, and E (47 kDa vicilins) and types C and D (50 kDa), can be distinguished in lentil and other Vicieae species, and that in the Vicieae species there is no evidence of concerted evolution among the vicilin sequences of different gene subfamilies or sequences groups, as has been suggested for other legume species.      

Gene conversion and concerted evolution in bacterial genomes

Gene conversion is defined as the non-reciprocal transfer of information between homologous sequences. Despite methodological problems to establish non-reciprocity, gene conversion has been demonstrated in a wide variety of bacteria. Besides examples of high-frequency reversion of mutations in repeated genes, gene conversion in bacterial genomes has been implicated in concerted evolution of multigene families. Gene conversion also has a prime importance in the generation of antigenic variation, an interesting mechanism whereby some bacterial pathogens are able to avoid the host immune system. In this review, we analyze examples of bacterial gene conversion (some of them spawned from the current genomic revolution), as well as the molecular models that explain gene conversion and its association with crossovers.     

Chromosomal inversions promote genomic islands of concerted evolution of Hsp70 genes in the Drosophila subobscura species subgroup

Heat‐shock (HS) assays to understand the connection between standing inversion variation and evolutionary response to climate change in Drosophila subobscura found that “warm‐climate” inversion O₃₊₄ exhibits non‐HS levels of Hsp70 protein like those of “cold‐climate” O_ST after HS induction. This was unexpected, as overexpression of Hsp70 can incur multiple fitness costs. To understand the genetic basis of this finding, we have determined the genomic sequence organization of the Hsp70 family in four different inversions, including O_ST, O₃₊₄ , O₃₊₄₊₈ and O₃₊₄₊₁₆ , using as outgroups the remainder of the subobscura species subgroup, namely Drosophila madeirensis and Drosophila guanche. We found (i) in all the assayed lines, the Hsp70 family resides in cytological locus 94A and consists of only two genes, each with four HS elements (HSEs) and three GAGA sites on its promoter. Yet, in O_ST, the family is comparatively more compact; (ii) the two Hsp70 copies evolve in concert through gene conversion, except in D. guanche; (iii) within D. subobscura, the rate of concerted evolution is strongly structured by inversion, being higher in O_ST than in O₃₊₄ ; and (iv) in D. guanche, the two copies accumulated multiple differences, including a newly evolved “gap‐type” HSE2. The absence of concerted evolution in this species may be related to a long‐gone‐unnoticed observation that it lacks Hsp70 HS response, perhaps because it has evolved within a narrow thermal range in an oceanic island. Our results point to a previously unrealized link between inversions and concerted evolution, with potentially major implications for understanding genome evolution.     

Rapid concerted evolution in animal mitochondrial DNA

Recombinational genetic processes are thought to be rare in the uniparentally inherited mitochondrial (mt) DNA molecules of vertebrates and other animals. Here, however, we document extremely rapid concerted microevolution, probably mediated by frequent gene conversion events, of duplicated sequences in the mtDNA control region of mangrove killifishes (Kryptolebias marmoratus). In local populations, genetic distances between paralogous loci within an individual were typically smaller (and often zero) than those between orthologous loci in different specimens. These findings call for the recognition of concerted evolution as a microevolutionary process and gene conversion as a likely recombinational force in animal mtDNA. The previously unsuspected power of these molecular phenomena could greatly impact mtDNA dynamics within germ cell lineages and in local animal populations.     

吸虱亚目trL1(tag)和trnL2(taa)基因间独立进化和协同进化的研究

[目的]吸虱亚目是真兽类动物的永久专性体表寄生虫.吸虱亚目的线粒体基因组发生了剧烈的裂化,形成了不同于以往典型单一大环的多个裂化微环.本文对17种吸虱以及外群尖叫虱Bothriometopus macrocnemis的trnL1(tag)和trnL2(taa)基因序列进行比较分析,探讨trnL1(tag)和trnL2(taa)基因间的独立进化和协同进化.[方法]对云南采集到的3科3属4种吸虱(弯多板虱Polyplax reclinata、锯多板虱Polyplax serrata、太平洋甲胁虱Hoplopleura pacifica和麝鼩钩板虱Ancistroplax crocidurae),用Illumina MiSeq PE250平台高通量测序后与GenBank中查找的其它13种吸虱及尖叫虱的trnL1(tag)和trnL2(taa)基因序列进行比较,采用最大简约法(Maximum parsimony,MP)分析17种吸虱trnL1(tag)和trnL2(taa)基因的进化关系.[结果]17种吸虱的trnL1(tag)和trnL2(taa)基因均形成典型的三叶草结构,trnL1(tag)和trnL2(taa)基因在虱属、阴虱属和猴虱属中有较长的等同序列;在甲胁虱属和钩板虱属中有很短的等同序列;在多板虱属、微胸虱属和血虱属中的等同序列介于二者之间.常见典型单一环状线粒体基因组物种的trnL1(tag)和trnL2(taa)基因等同序列的长度在6-10 bp之间.对17种吸虱的trnL1(tag)和trnL2(taa)基因序列进行同源性比较分析,发现不同属吸虱的trnL1(tag)和trnL2(taa)基因序列相似度差异较大,同属内吸虱的trnL1(tag)和trnL2(taa)基因序列相似度差异较小.基于17种吸虱及尖叫虱的trnL1(tag)和trnL taa)基因构建系统进化树,结果表明trnL1(tag)、trnL2(taa)及trnL1(tag)和trnL2(taa)基因间既有协同进化又有独立进化.[结论]trnL 1(tag)和trnL 2(taa)基因等同序列较长的吸虱易发生协同进化,等同序列较短的吸虱易发生独立进化.吸虱亚目的trnL tag)和trnL2(taa)基因协同进化是长期的或发生在远缘物种间,而独立进化是短期(两次重组事件间)的或发生在近缘物种间.吸虱亚目线粒体基因组的裂化模式可能影响trnL1(tag)和trnL2(taa)基因等同序列的长短.     

Selection for more of the same product as a force to enhance concerted evolution of duplicated genes

The duration of concerted evolution after gene duplication is highly variable across genes. To identify the cause of the variation, we analyzed of duplicated genes in yeast that originate from a whole genome duplication event. There appears to be a strong positive correlation between the duration of concerted evolution and the gene expression level. This observation can be explained by selection favoring more of the same product, which could enhance concerted evolution in dosage-sensitive genes.     

Copy number variants in pharmacogenetic genes

Variation in drug efficacy and toxicity remains an important clinical concern. Presently, single nucleotide polymorphisms (SNPs) only explain a portion of this problem, even in situations where the pharmacological trait is clearly heritable. The Human CNV Project identified copy number variations (CNVs) across approximately 12% of the human genome, and these CNVs were considered causes of diseases. Although the contribution of CNVs to the pathogenesis of many common diseases is questionable, CNVs play a clear role in drug-related genes by altering drug metabolizing and drug response. In this review, we provide a comprehensive evaluation of the clinical relevance of CNVs to drug efficacy, toxicity, and disease prevalence in world populations, and discuss the implication of using CNVs as a diagnostic tool in clinical intervention.     

The 1.688 repetitive DNA of Drosophila: concerted evolution at different genomic scales and association with genes

Concerted evolution leading to homogenization of tandemly repeated DNA arrays is widespread and important for genome evolution. We investigated the range and nature of the process at chromosomal and array levels using the 1.688 tandem repeats of Drosophila melanogaster where large arrays are present in the heterochromatin of chromosomes 2, 3, and X, and short arrays are found in the euchromatin of the same chromosomes. Analysis of 326 euchromatic and heterochromatic repeats from 52 arrays showed that the homogenization of 1.688 repeats occurred differentially for distinct genomic regions, from euchromatin to heterochromatin and from local arrays to chromosomes. We further found that most euchromatic arrays are either close to, or are within introns of, genes. The short size of euchromatic arrays (one to five repeats) could be selectively constrained by their role as gene regulators, a situation similar to the so-called "tuning knobs."