Adaptive evolution of recently duplicated accessory gland protein genes in desert Drosophila

The relationship between animal mating system variation and patterns of protein polymorphism and divergence is poorly understood. Drosophila provides an excellent system for addressing this issue, as there is abundant interspecific mating system variation. For example, compared to D. melanogaster subgroup species, repleta group species have higher remating rates, delayed sexual maturity, and several other interesting differences. We previously showed that accessory gland protein genes (Acp's) of Drosophila mojavensis and D. arizonae evolve more rapidly than Acp's in the D. melanogaster subgroup and that adaptive Acp protein evolution is likely more common in D. mojavensis/D. arizonae than in D. melanogaster/D. simulans. These findings are consistent with the idea that greater postcopulatory selection results in more adaptive evolution of seminal fluid proteins in the repleta group flies. Here we report another interesting evolutionary difference between the repleta group and the D. melanogaster subgroup Acp's. Acp gene duplications are present in D. melanogaster, but their high sequence divergence indicates that the fixation rate of duplicated Acp's has been low in this lineage. Here we report that D. mojavensis and D. arizonae genomes contain several very young duplicated Acp's and that these Acp's have experienced very rapid, adaptive protein divergence. We propose that rapid remating of female desert Drosophila generates selection for continuous diversification of the male Acp complement to improve male fertilization potential. Thus, mating system variation may be associated with adaptive protein divergence as well as with duplication of Acp's in Drosophila.     

Evolutionary genomics of the recently duplicated amphioxus Hairy genes

Amphioxus Hairy genes have gone through a number of lineage-specific duplications, resulting in eight members, some of which are differentially expressed in the embryo. In order to gain insights into the evolution and function of this gene family we have compared their genomic structure and searched for conserved non-coding sequence elements. We have found that introns have been lost independently from these genes at least twice and after the duplication events. By carrying out phylogenetic footprinting between paralogues expressed in the embryo, we have found a differential distribution of conserved elements that could explain the limited overlap in expression patterns of Hairy genes in the amphioxus embryo. Furthermore, clustering of RBP-Jk binding sites in these conserved elements suggests that amphioxus Hairy genes are downstream targets of the Notch signaling pathway, as occurs in vertebrates. All of this evidence suggests that amphioxus Hairy genes have gone through a process of subfunctionalization shortly after their duplication, representing an extreme and rapid case of the duplication-degeneration-complementation model.     

Molecular and functional evolution of human DHRS2 and DHRS4 duplicated genes

Human DHRS2 and DHRS4 genes code for similar NADP-dependent short-chain carbonyl-reductase enzymes having different substrate specificity. Human DHRS2 and DHRS4 enzymes share several common sequence motives including residues responsible for coenzyme binding as well as for the intimate catalytic oxido-reductase mechanism, while their substrate-binding sequences have very low similarity. We found that DHRS2 and DHRS4 genes are syntenic outparalogues originated from a duplication of the DHRS4 gene that took place before the formation of the mammalian clade. DHRS2 gene evolved more rapidly and underwent positive selection on more sites than the DHRS4 gene. DHRS2 sites under positive selection were mainly located on the enzyme active site thus showing that substrate specificity drove the divergence from the DHRS4 enzyme. Rapid divergent evolution brought the human DHRS2 enzyme to have subcellular localization, synthesis regulation and specialized cellular functions very different from those of the human DHRS4 enzyme.     

The major shifts of human duplicated genes

Since many gene duplications in the human genome are ancient duplications going back to the origin of vertebrates, the question may be asked about the fate of such duplicated genes at the compositional genome transitions that occurred between cold- and warm-blooded vertebrates. Indeed, at that transition, about half of the (GC-poor) genes of cold-blooded vertebrates (the genes of the gene-dense "ancestral genome core") underwent a GC enrichment to become the genes of the "genome core" of warm-blooded vertebrates. Since the compositional distribution of the human duplicated genes investigated (1111 pairs) mimics the general distribution of human genes (about 50% GC(3)-poor and 50% GC(3)-rich genes, the border being at 60% GC(3)), we considered two possibilities, namely that the compositional transition affected either (i) about half of the copies on a random basis, or (ii) preferentially only one copy of the duplicated genes. The two possibilities could be distinguished if each copy is put into one of two subsets according to its GC(3) level. Indeed, in the first case, the two distributions would be similar, whereas in the second case, the two distributions would be different, one copy having maintained the ancestral GC-poor composition, and one copy having undergone the compositional change. Using this approach, we could show that, by far and large, one copy of the duplicated genes preferentially underwent the GC enrichment. This result implies that this copy, which had possibly acquired a different function and/or regulation, was preferentially translocated into the gene-dense compartment of the genome, the "ancestral genome core", namely the "gene space" which underwent the compositional transition at the emergence of warm-blooded vertebrates.     

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.     

Selective maintenance of Drosophila tandemly arranged duplicated genes during evolution

Background

The physical organization and chromosomal localization of genes within genomes is known to play an important role in their function. Most genes arise by duplication and move along the genome by random shuffling of DNA segments. Higher order structuring of the genome occurs in eukaryotes, where groups of physically linked genes are co-expressed. However, the contribution of gene duplication to gene order has not been analyzed in detail, as it is believed that co-expression due to recent duplicates would obscure other domains of co-expression.

Results

We have catalogued ordered duplicated genes in Drosophila melanogaster, and found that one in five of all genes is organized as tandem arrays. Furthermore, among arrays that have been spatially conserved over longer periods than would be expected on the basis of random shuffling, a disproportionate number contain genes encoding developmental regulators. Using in situ gene expression data for more than half of the Drosophila genome, we find that genes in these conserved clusters are co-expressed to a much higher extent than other duplicated genes.

Conclusions

These results reveal the existence of functional constraints in insects that retain copies of genes encoding developmental and regulatory proteins as neighbors, allowing their co-expression. This co-expression may be the result of shared cis-regulatory elements or a shared need for a specific chromatin structure. Our results highlight the association between genome architecture and the gene regulatory networks involved in the construction of the body plan.     

Molecular evolution and population genetics of duplicated accessory gland protein genes in Drosophila

To investigate the potential importance of gene duplication in D. melanogaster accessory gland protein (Acp) gene evolution we carried out a computational analysis comparing annotated D. melanogaster Acp genes to the entire D. melanogaster genome. We found that two known Acp genes are actually members of small multigene families. Polymorphism and divergence data from these duplicated genes suggest that in at least four cases, protein divergence between D. melanogaster and D. simulans is a result of directional selection. One putative Acp revealed by our computational analysis shows evidence of a recent selective sweep in a non-African population (but not in an African population). These data support the idea that selection on reproduction-related genes may drive divergence of populations within species, and strengthen the conclusion that Acps may often be under directional selection in Drosophila.     

Independent recombination events between the duplicated human alpha globin genes; implications for their concerted evolution.

We have examined the molecular structure of the human alpha globin gene complex from individuals with a common form of alpha thalassaemia in which one of the duplicated pair of alpha genes (alpha alpha) has been deleted (-alpha 3-7). Restriction mapping and DNA sequence analysis of the mutants indicate that different -alpha 3.7 chromosomes are the result of at least three independent events. In each case the genetic crossover has occurred within a region of complete homology between the alpha 1 and alpha 2 genes. Since the -alpha chromosomes may reflect the processes of crossover fixation and gene conversion between the two genes, their structures may provide some insight into the mechanism by which the concerted evolution of the human alpha globin genes occurs.     

Nonuniform evolution of duplicated,developmentally controlled chorion genes in a silkmoth

Summary We report the sequence of A/B.L1, a tightly linked pair of genes from the A and B chorion families inBombyx mori. Comparison with the previously characterized A/B.L11 and A/B.L12 pairs from the same species reveals moderate sequence divergence, which is clearly nonuniform. Although the average divergence of A/B.L12 from the other two pairs is more than double that between A/B.L11 and A/B.L1, the ratio differs by more than 30-fold in different DNA regions. One domain of the A gene is highly divergent between A/B.L12 and A/B.L1 or A/B.L11, but essentially invariable in the latter two. In well-aligned DNA segments, nearly all mutated sites (111/112) show variants shared by two of the three sequences, in 42% of the cases between the more distantly related pairs (A/B.L12 and either A/B.L1 or A/B.L11). Eight of the variants shared by distantly related pairs are clustered within 51 bp, suggesting the possibility that they arose through sequence transfers between gene pairs, rather than being primitive or resulting from independent mutations. The short intergenic, putatively regulatory DNa of the developmentally middle A/B.L1 and A/B.L11 pairs resembles that of the late HcA/HcB pairs, particularly in patches that may correspond tocis-regulatory elements.     

Conservative evolution in duplicated genes of the primate Class I ADH cluster

Humans have seven alcohol dehydrogenase genes (ADH) falling into five classes. Three out of the seven genes (ADH1A, ADH1B and ADH1C) belonging to Class I are expressed primarily in liver and code the main enzymes catalyzing ethanol oxidization. The three genes are tandemly arrayed within the ADH cluster on chromosome 4 and have very high nucleotide similarity to each other (exons: >90%; introns: >70%), suggesting the genes have been generated by duplication event(s). One explanation for maintaining similarity of such clustered genes is homogenization via gene conversion(s). Alternatively, recency of the duplications or some other functional constraints might explain the high similarities among the genes. To test for gene conversion, we sequenced introns 2, 3, and 8 of all three Class I genes (total>15.0 kb) for five non-human primates--four great apes and one Old World Monkey (OWM)--and compared them with those of humans. The phylogenetic analysis shows each intron sequence clusters strongly within each gene, giving no evidence for gene conversion(s). Several lines of evidence indicate that the first split was between ADH1C and the gene that gave rise to ADH1A and ADH1B. We also analyzed cDNA sequences of the three genes that have been previously reported in mouse and Catarrhines (OWMs, chimpanzee, and humans) and found that the synonymous and non-synonymous substitution (dN/dS) ratios in all pairs are less than 1 representing purifying selection. This suggests that purifying selection is more important than gene conversion(s) in maintaining the overall sequence similarity among the Class I genes. We speculate that the highly conserved sequences on the three duplicated genes in primates have been achieved essentially by maintaining stability of the hetero-dimer formation that might have been related to dietary adaptation in primate evolution.     

Selective maintenance of Drosophila tandemly arranged duplicated genes during evolution

Background  

The physical organization and chromosomal localization of genes within genomes is known to play an important role in their function. Most genes arise by duplication and move along the genome by random shuffling of DNA segments. Higher order structuring of the genome occurs in eukaryotes, where groups of physically linked genes are co-expressed. However, the contribution of gene duplication to gene order has not been analyzed in detail, as it is believed that co-expression due to recent duplicates would obscure other domains of co-expression.     

Asymmetric evolution in two fish-specifically duplicated receptor tyrosine kinase paralogons involved in teleost coloration

The occurrence of a fish-specific genome duplication (FSGD) in the lineage leading to teleost fishes is widely accepted, but the consequences of this event remain elusive. Teleosts, and the cichlid fishes from the species flocks in the East African Great Lakes in particular, evolved a unique complexity and diversity of body coloration and color patterning. Several genes involved in pigment cell development have been retained in duplicate copies in the teleost genome after the FSGD. Here we investigate the evolutionary fate of one of these genes, the type III receptor tyrosine kinase (RTK) colony-stimulating factor 1 receptor (csf1r). We isolated and shotgun sequenced two paralogous csf1r genes from a bacterial artificial chromosome library of the cichlid fish Astatotilapia burtoni that are both linked to paralogs of the pdgfr beta gene, another type III RTK. Two pdgfr beta-csf1r paralogons were also identified in the genomes of pufferfishes and medaka, and our phylogenetic analyses suggest that the pdgfr beta-csf1r locus was duplicated during the course of the FSGD. Comparisons of teleosts and tetrapods suggest asymmetrical divergence at different levels of genomic organization between the teleost-specific pdgfr beta-csf1r paralogons, which seem to have evolved as coevolutionary units. The high-evolutionary rate in the teleost B-paralogon, consisting of csf1rb and pdgfr betab, further suggests neofunctionalization by functional divergence of the extracellular, ligand-binding region of these cell-surface receptors. Finally, we hypothesize that genome duplications and the associated expansion of the RTK family might be causally linked to the evolution of coloration in vertebrates and teleost fishes in particular.     

The evolutionary fate of recently duplicated retrogenes in mice

Inferences about the evolutionary impact of gene duplications often rely on the analysis of their long-term outcome. The fate of the majority of them must, however, be decided shortly after duplication. Here we analysed the evolutionary pattern of 10 mouse genes very recently duplicated by retrotransposition, by sequencing the retroposed copy in five to 10 closely related mouse species. In all cases the retroposed copy experienced accelerated nonsynonymous evolution whereas the divergence pattern of the source copy appeared unaffected by the duplication, consistent with the neofunctionalization model. The analysis further revealed that most retrogenes, including pseudogenes, did not experience a period of relaxed neutral evolution, but have been submitted to purifying selection ever since their retroposition. We propose that these duplicates play a biochemical role but are not indispensable. Purifying selection prevents them from acquiring a negative role until they are lost or silenced. This period of unnecessary redundancy could in rare cases give the time for new functions to evolve.     

The evolution and sequence comparison of two recently diverged mouse chromosomal beta--globin genes

We have determined the entire nucleotide sequence of a cloned mouse beta--globinminor gene and compared it to the closely related sequence of the betamajor gene. These two genes differ by nine amino acids and presumably evolved from a common ancestral gene as recently as 50 million years ago. Since these genes are closely linked and coordinately expressed, they provide an especially favorable opportunity to assess selection and mutation as these processes affect genes under similar constraints. We find that evolution has preserved these two genes in two short segments of DNA which include their immediately adjacent flanking regions. These regions presumably encode functions that are necessary for proper globin gene expression. In contrast, the more distal flanking sequences and major segments of the long intervening sequences have diverged much more sharply. The homology pattern in these genes also provides considerable insight into the mechanisms by which less constrained nucleotide sequences diverge rapidly. Change in such regions apparently occurs less by point mutation than by insertion, deletion and duplication of relatively short segments of the genome.     

The ovalbumin gene family: Structure of the X gene and evolution of duplicated split genes

The X, Y and ovalbumin genes, which are found within a 40 kb region of the chicken genome, are all expressed in oviduct under steroid hormone control, and share some sequence homologies. We have now cloned the complete X gene and have analyzed its structure. It codes for two RNA species, X and X′; both are coded by eight exons and appear to differ only by the size of their 3′ untranslated region, X′ RNA being 1400 nucleotides longer than X RNA. The striking similarity in the number and length of the exons which constitute the X, Y or ovalbumin genes establishes that they have evolved from a common ancestor gene by duplication events. Comparison of selected regions of the X and ovalbumin genes indicates that the exon sequences coding for protein and the location of the splice junctions have been well-conserved. The introns and the 3′ untranslated exonic sequences have diverged much more rapidly. Four regions of apparently unrelated repetitive sequences are found both outside the X gene and within it (in two introns and in the sequence coding for the 3′ untranslated part of X′RNA). The intragenic repetitive sequences have no counterpart in the ovalbumin and Y genes.