Population genetic models of duplicated genes

Various population genetic models of duplicated genes are introduced. The problems covered in this review include the fixation process of a duplicated copy, copy number polymorphism, the fates of duplicated genes and single nucleotide polymorphism in duplicated genes. Because of increasing evidence for concerted evolution by gene conversion, this review introduces recently developed gene conversion models. In the first half, models assuming independent evolution of duplicated genes are introduced, and then the effect of gene conversion is considered in the second half.     

Neofunctionalization of duplicated genes under the pressure of gene conversion

Neofunctionalization occurs when a neofunctionalized allele is fixed in one of duplicated genes. This is a simple fixation process if duplicated genes accumulate mutations independently. However, the process is very complicated when duplicated genes undergo concerted evolution by gene conversion. Our simulations demonstrate that the process could be described with three distinct stages. First, a newly arisen neofunctionalized allele increases in frequency by selection, but gene conversion prevents its complete fixation. These two factors (selection and gene conversion) that work in opposite directions create an equilibrium, and the time during which the frequency of the neofunctionalized allele drifts around the equilibrium value is called the temporal equilibrium stage. During this temporal equilibrium stage, it is possible that gene conversion is inactivated by mutations, which allow the complete fixation of the neofunctionalized allele. And then, permanent neofunctionalization is achieved. This article develops basic population genetics theories on the process to permanent neofunctionalization under the pressure of gene conversion. We obtain the probability and time that the frequency of a newly arisen neofunctionalized allele reaches the equilibrium value. It is also found that during the temporal equilibrium stage, selection exhibits strong signature in the divergence in the DNA sequences between the duplicated genes. The spatial distribution of the divergence likely has a peak around the site targeted by selection. We provide an analytical expression of the pattern of divergence and apply it to the human red- and green-opsin genes. The theoretical prediction well fits the data when we assume that selection is operating for the two amino acid differences in exon 5, which are believed to account for the major part of the functional difference between the red and green opsins.     

Preservation of a pseudogene by gene conversion and diversifying selection

Interlocus gene conversion is considered a crucial mechanism for generating novel combinations of polymorphisms in duplicated genes. The importance of gene conversion between duplicated genes has been recognized in the major histocompatibility complex and self-incompatibility genes, which are likely subject to diversifying selection. To theoretically understand the potential role of gene conversion in such situations, forward simulations are performed in various two-locus models. The results show that gene conversion could significantly increase the number of haplotypes when diversifying selection works on both loci. We find that the tract length of gene conversion is an important factor to determine the efficacy of gene conversion: shorter tract lengths can more effectively generate novel haplotypes given the gene conversion rate per site is the same. Similar results are also obtained when one of the duplicated genes is assumed to be a pseudogene. It is suggested that a duplicated gene, even after being silenced, will contribute to increasing the variability in the other locus through gene conversion. Consequently, the fixation probability and longevity of duplicated genes increase under the presence of gene conversion. On the basis of these findings, we propose a new scenario for the preservation of a duplicated gene: when the original donor gene is under diversifying selection, a duplicated copy can be preserved by gene conversion even after it is pseudogenized.     

A three-step model for the rearrangement of the chloroplast trnK-psbA region of the gymnosperm Pinus contorta.

A region of the Pinus contorta chloroplast genome which contains a duplication of the psbA gene was characterized. From previous experiments it was known that the two copies of the psbA gene were located approximately 3.3 kilobase pairs (kbp) apart, that they had the same orientation and that one endpoint of the duplication was 19 base pairs (bp) downstream of the psbA stop codon. In order to determine the size and additional genetic content of the duplicated segment, both copies as well as the intervening DNA were sequenced completely. It was found that the duplicated segment was 1969 bp long, that the two copies were completely identical and were separated by 2431 bp. The duplicated segment carried, in addition to psbA, the 3' exon of the trnK gene, which was partially included in a 124 bp direct repeat. The translocated copy of the duplicated segment was found to be inserted upstream of the trnK(UUU) gene and was immediately followed by a repeated 41 bp stretch from the psbA coding region. The trnK gene was split by a 2509 bp intron which contained an open reading frame of 515 codons. Sequence comparisons of the duplicated segment and its flanking DNA to the corresponding regions of P. sylvestris, a species which lacks the rearrangements found in P. contorta, made it possible to identify 3-9 bp homologies within which recombinations had occurred. A model was derived which would accommodate the conversion of a trnK-psbA locus of the ancestral P. sylvestris-like organization into the rearranged structure found in P. contorta.     

Population-genetic models of the fates of duplicate genes

The ultimate fate of a duplicated gene is that it either silenced through inactivating mutations or both copies are maintained by selection. This later fate can occur via neofunctionalization wherein one copy acquires a new function or by subfunctionalization wherein the original function of the gene is partitioned across both copies. The relative probabilities of these three different fates involve often very subtle iterations between of population size, mutation rate, and selection. All three of these fates are critical to the expansion and diversification of gene families.     

Analysis of two genes encoding heat-labile toxins and located on a single Ent plasmid from Escherichia coli.

A clinical isolate of Escherichia coli harbored two copies of the heat-labile toxin (LT)-encoding gene (elt) on a 157-kb plasmid. The arrangement of the gene copies is different from the cholera toxin-encoding gene duplication described for some strains of Vibrio cholerae. The nucleotide sequences of the elt alleles are not identical (differing by 2 bp) and the duplicated region flanking the alleles extends 314 bp on one side and 1122 bp on the other side of each copy. Different partial copies of IS600 were identified 280 bp 3' to the stop codon of each allele. Partial and complete copies of other IS were also identified near the elt alleles. The data suggest that the regions surrounding the genes are hot spots for IS which would account for the observed heterogeneity in DNA flanking elt genes.     

Duplication and Gene Conversion in the Drosophila melanogaster Genome

Using the genomic sequences of Drosophila melanogaster subgroup, the pattern of gene duplications was investigated with special attention to interlocus gene conversion. Our fine-scale analysis with careful visual inspections enabled accurate identification of a number of duplicated blocks (genomic regions). The orthologous parts of those duplicated blocks were also identified in the D. simulans and D. sechellia genomes, by which we were able to clearly classify the duplicated blocks into post- and pre-speciation blocks. We found 31 post-speciation duplicated genes, from which the rate of gene duplication (from one copy to two copies) is estimated to be 1.0×10⁻⁹ per single-copy gene per year. The role of interlocus gene conversion was observed in several respects in our data: (1) synonymous divergence between a duplicated pair is overall very low. Consequently, the gene duplication rate would be seriously overestimated by counting duplicated genes with low divergence; (2) the sizes of young duplicated blocks are generally large. We postulate that the degeneration of gene conversion around the edges could explain the shrinkage of “identifiable” duplicated regions; and (3) elevated paralogous divergence is observed around the edges in many duplicated blocks, supporting our gene conversion–degeneration model. Our analysis demonstrated that gene conversion between duplicated regions is a common and genome-wide phenomenon in the Drosophila genomes, and that its role should be especially significant in the early stages of duplicated genes. Based on a population genetic prediction, we applied a new genome-scan method to test for signatures of selection for neofunctionalization and found a strong signature in a pair of transporter genes.     

Selection for functional uniformity of tuf duplicates in gamma-proteobacteria

Having an extra copy of a gene is thought to provide some functional redundancy, which results in a higher rate of evolution in duplicated genes. In this article, we estimate the impact of gene duplication on the selection of tuf paralogs, and we find that in the absence of gene conversion, tuf paralogs have evolved significantly slower than when gene conversion has been a factor in their evolution. Thus, tuf gene copies evolve under a selective pressure that ensures their functional uniformity, and gene conversion reduces selection against amino acid substitutions that affect the function of the encoded protein, EF-Tu.     

Tandem duplication of the NF1 gene detected by high-resolution FISH in the 17q11.2 region

The gene for neurofibromatosis type 1 (NF1), mapping to 17q11.2, has one of the highest observed mutation rates, partially because of its large size and gene conversion primed by NF1 pseudogenes. We have previously shown by means of high resolution fluorescence in situ hybridization (FISH) that a number of the loci flanking the NF1 gene are duplicated, in agreement with the reported presence of NF1 repetitive sequences (REPs). We report a direct tandem duplication of the NF1 gene identified in 17q11.2 by high-resolution FISH. FISH on stretched chromosomes with locus-specific probes revealed the duplication of the NF1 gene from the promoter to 3'UTR, but with at least the absence of exon 22. Fiber FISH with P1 artificial and bacterial artifical chromosomes, including the NF1 5'UTR and 3'UTR and flanking regions, visualized the direct tandem duplication with a similar, but not identical, genomic organization of the NF1 duplicon copies. Duplication was probably present in the human-chimpanzee-gorilla common ancestor, as demonstrated here by the finding of the duplicated NF1 gene at orthologous chromosome loci. The NF1 intrachromosomal duplication may contribute to the high whole-gene mutation rate by gene conversion, although the functional activity of the NF1 copy remains to be investigated. Detection of the NF1 duplicon by high-resolution FISH may pave the way to filling the gaps in the human genomic sequence of the pericentromeric 17q11.2 region.     

Recurrent gene amplification and soft selective sweeps during evolution of multidrug resistance in malaria parasites

When selection is strong and beneficial alleles have a single origin, local reductions in genetic diversity are expected. However, when beneficial alleles have multiple origins or were segregating in the population prior to a change in selection regime, the impact on genetic diversity may be less clear. We describe an example of such a "soft" selective sweep in the malaria parasite Plasmodium falciparum that involves adaptive genome rearrangements. Amplification in copy number of genome regions containing the pfmdr1 gene on chromosome 5 confer resistance to mefloquine and spread rapidly in the 1990s. Using flanking microsatellite data and real-time polymerase chain reaction determination of copy number, we show that 5-15 independent amplification events have occurred in parasites on the Thailand/Burma border. The amplified genome regions (amplicons) range in size from 14.7 to 49 kb and contain 2-11 genes, with 2-4 copies arranged in tandem. To examine the impact of drug selection on flanking variation, we genotyped 48 microsatellites on chromosome 5 in 326 parasites from a single Thai location. Diversity was reduced in a 170- to 250-kb (10-15 cM) region of chromosomes containing multiple copies of pfmdr1, consistent with hitchhiking resulting from the rapid recent spread of selected chromosomes. However, diversity immediately flanking pfmdr1 is reduced by only 42% on chromosomes bearing multiple amplicons relative to chromosomes carrying a single copy. We highlight 2 features of these results: 1) All amplicon break points occur in monomeric A/T tracts (9-45 bp). Given the abundance of these tracts in P. falciparum, we expect that duplications will occur frequently at multiple genomic locations and have been underestimated as drivers of phenotypic evolution in this pathogen. 2) The signature left by the spread of amplified genome segments is broad, but results in only limited reduction in diversity. If such "soft" sweeps are common in nature, statistical methods based on diversity reduction may be inefficient at detecting evidence for selection in genome-wide marker screens. This may be particularly likely when mutation rate is high, as appears to be the case for gene duplications, and in pathogen populations where effective population sizes are typically very large.     

Molecular evolution of psbA gene in ferns: unraveling selective pressure and co-evolutionary pattern

ABSTRACT: BACKGROUND: The photosynthetic oxygen-evolving photo system II (PS II) produces almost the entire oxygen in the atmosphere. This unique biochemical system comprises a functional core complex that is encoded by psbA and other genes. Unraveling the evolutionary dynamics of this gene is of particular interest owing to its direct role in oxygen production. psbA underwent gene duplication in leptosporangiates, in which both copies have been preserved since. Because gene duplication is often followed by the non-fictionalization of one of the copies and its subsequent erosion, preservation of both psbA copies pinpoint functional or regulatory specialization events. The aim of this study was to investigate the molecular evolution of psbA among fern lineages. RESULTS: We sequenced psbA, which encodes D1 protein in the core complex of PSII, in 20 species representing 8 orders of extant ferns; then we searched for selection and convolution signatures in psbA across the 11 fern orders. Collectively, our results indicate that: (1) selective constraints among D1 protein relaxed after the duplication in 4 leptosporangiate orders; (2) a handful positively selected codons were detected within species of single copy psbA, but none in duplicated ones; (3) a few sites among D1 protein were involved in co-evolution process which may intimate significant functional/structural communications between them. CONCLUSIONS: The strong competition between ferns and angiosperms for light may have been the main cause for a continuous fixation of adaptive amino acid changes in psbA, in particular after its duplication. Alternatively, a single psbA copy may have undergone bursts of adaptive changes at the molecular level to overcome angiosperms competition. The strong signature of positive Darwinian selection in a major part of D1 protein is testament to this. At the same time, species own two psbA copies hardly have positive selection signals among the D1 protein coding sequences. In this study, eleven co-evolving sites have been detected via different molecules, which may be more important than others.     

Gene conversion drives the evolution of HINTW, an ampliconic gene on the female-specific avian W chromosome

The HINTW gene on the female-specific W chromosome of chicken and other birds is amplified and present in numerous copies. Moreover, as HINTW is distinctly different from its homolog on the Z chromosome (HINTZ), is a candidate gene in avian sex determination, and evolves rapidly under positive selection, it shows several common features to ampliconic and testis-specific genes on the mammalian Y chromosome. A phylogenetic analysis within galliform birds (chicken, turkey, quail, and pheasant) shows that individual HINTW copies within each species are more similar to each other than to gene copies of related species. Such convergent evolution is most easily explained by recurrent events of gene conversion, the rate of which we estimated at 10(-6)-10(-5) per site and generation. A significantly higher GC content of HINTW than of other W-linked genes is consistent with biased gene conversion increasing the fixation probability of mutations involving G and C nucleotides. Furthermore, and as a likely consequence, the neutral substitution rate is almost twice as high in HINTW as in other W-linked genes. The region on W encompassing the HINTW gene cluster is not covered in the initial assembly of the chicken genome, but analysis of raw sequence reads indicates that gene copy number is significantly higher than a previous estimate of 40. While sexual selection is one of several factors that potentially affect the evolution of ampliconic, male-specific genes on the mammalian Y chromosome, data from HINTW provide evidence that gene amplification followed by gene conversion can evolve in female-specific chromosomes in the absence of sexual selection. The presence of multiple and highly similar copies of HINTW may be related to protein function, but, more generally, amplification and conversion offers a means to the avoidance of accumulation of deleterious mutations in nonrecombining chromosomes.     

A diffusion model for the fate of tandem gene duplicates in diploids

Suppose one chromosome in one member of a population somehow acquires a duplicate copy of the gene, fully linked to the original gene's locus. Preservation is the event that eventually every chromosome in the population is a descendant of the one which initially carried the duplicate. For a haploid population in which the absence of all copies of the gene is lethal, the probability of preservation has recently been estimated via a diffusion approximation. That approximation is shown to carry over to the case of diploids and arbitrary strong selection against the absence of the gene. The techniques used lead to some new results. In the large population limit, it is shown that the relative probability that descendants of a small number of individuals carrying multiple copies of the gene fix in the population is proportional to the number of copies carried. The probability of preservation is approximated when chromosomes carrying two copies of the gene are subject to additional, fully non-functionalizing mutations, thereby modelling either an additional cost of replicating a longer genome, or a partial duplication of the gene. In the latter case the preservation probability depends only on the mutation rate to null for the duplicated portion of the gene.     

A tandem segmental duplication (TSD) in green revolution gene Rht-D1b region underlies plant height variation

? Rht-D1c (Rht10) carried by Chinese wheat (Triticum aestivum) line Aibian 1 is an allele at the Rht-D1 locus. Among the Rht-1 alleles, little is known about Rht-D1c although it determines an extreme dwarf phenotype in wheat. ? Here, we cloned and functionally characterized Rht-D1c using a combination of Southern blotting, target region sequencing, gene expression analysis and transgenic experiments. ? We found that the Rht-D1c allele was generated through a tandem segmental duplication (TSD) of a >?1?Mb region, resulting in two copies of the Rht-D1b. Two copies of Rht-D1b in the TSD were three-fold more effective in reducing plant height than a single copy, and transformation with a segment containing the tandemly duplicated copy of Rht-D1b resulted in the same level of reduction of plant height as the original copy in Aibian 1. ? Our results suggest that changes in gene copy number are one of the important sources of genetic diversity and some of these changes could be directly associated with important traits in crops.     

How Does the Cell Count the Number of Ectopic Copies of a Gene in the Premeiotic Inactivation Process Acting in Ascobolus Immersus?

Repeated genes, artificially introduced in Ascobolus immersus by integrative transformation, are frequently inactivated during the sexual phase. Inactivation is observed in about 50% of meioses if duplicated genes are at ectopic chromosomal locations, and in 90% of meioses if genes are tandemly repeated. Inactivation is associated with extensive methylation of the cytosine residues of the duplicated sequences and is induced in the still haploid nuclei of the dikaryotic cell which will undergo karyogamy and subsequent meiosis. Only repeated sequences become methylated. This raises the intriguing question of how the premeiotic inactivation machinery is informed that a nucleus contains multiple copies of a gene. By using in crosses recombinant strains of A. immersus in which either one, two or three genetically independent copies of the exogenous amdS gene from Aspergillus nidulans had been introduced, we could follow the premeiotic inactivation of each one of the ectopic amdS copies. This led us to propose that a prerequisite for inactivation is a premeiotic pairing of repeated sequences and that each copy can undergo successive cycles of pairing. In fact, once methylated, a copy can pair with a still unmethylated copy, so that an uneven number of copies can be subject to inactivation.     

Evidence for gene conversion in the amylase multigene family of Drosophila pseudoobscura

The alpha-amylase (Amy) multigene family in Drosophila pseudoobscura is located on the third chromosome, which is polymorphic for more than 40 inverted gene arrangements. The number of copies in this family ranges from one to three, depending on the arrangement in question. A previous study of the three Amy genes from the Standard (ST) arrangement suggested either that duplicated copies (Amy2 and Amy3) are functionally constrained or that they are undergoing gene conversion with Amy1. In order to elucidate further the pattern of molecular evolution in this family, we cloned and sequenced four additional Amy genes, two from the Santa Cruz (SC) and two from the Chiricahua (CH) gene arrangement. Of the two alternatives, only the hypothesis of gene conversion is supported by the sequence analysis. The homogenization effect of gene conversion has been strongest in SC, whose copies differ by only two nucleotides, less noticeable in ST, and negligible in the CH. Furthermore, the action of gene conversion is apparently localized, occurring only in the coding region. Interestingly, these results concur with the findings of other workers for the duplicated Amy genes in the Drosophila melanogaster group. Thus, the occurrence of gene conversion in the Amy multigene family seems to be a common feature in the Drosophila species studied so far.      

The evolution of MHC diversity: evidence of intralocus gene conversion and recombination in a single-locus system

Gene conversion, the unidirectional exchange of genetic material between homologous sequences, is thought to strongly influence patterns of genetic diversity. The high diversity of major histocompatibility complex (MHC) genes in many species is thought to reflect a long history of gene conversion events both within and among loci. Theoretical work suggests that intra- and interlocus gene conversion leave characteristic signatures of nucleotide diversity, but empirical studies of MHC variation have rarely been able to analyze the effects of conversion events in isolation, due to the presence of multiple gene copies in most species. The potbellied seahorse (Hippocampus abdominalis), a species with a single copy of the MH class II beta-chain gene (MHIIb), provides an ideal system in which to explore predictions on the effects of intralocus gene conversion on patterns of genetic diversity. The genetic diversity of the MHIIb peptide binding region (PBR) is high in the seahorse, similar to other vertebrate species. In contrast, the remainder of the gene shows a total absence of synonymous variation and low levels of intronic sequence diversity, concentrated in 3 short repetitive regions and 1-12 SNPs per intron. The distribution of substitutions across the gene results in a patchwork pattern of shared polymorphism between otherwise divergent sequences. The pattern of nucleotide diversity observed in the seahorse MHIIb gene is congruent with theoretical expectations for intralocus gene conversion, indicating that this evolutionary mechanism has played an important role in MHC gene evolution, contributing to both the high diversity in the PBR and the low diversity outside this region. Neutral variation at this locus may be further reduced due to biases in nucleotide composition and functional constraints.     

The impact of natural selection on an ABCC11 SNP determining earwax type

A nonsynonymous single nucleotide polymorphism (SNP), rs17822931-G/A (538G>A; Gly180Arg), in the ABCC11 gene determines human earwax type (i.e., wet or dry) and is one of most differentiated nonsynonymous SNPs between East Asian and African populations. A recent genome-wide scan for positive selection revealed that a genomic region spanning ABCC11, LONP2, and SIAH1 genes has been subjected to a selective sweep in East Asians. Considering the potential functional significance as well as the population differentiation of SNPs located in that region, rs17822931 is the most plausible candidate polymorphism to have undergone geographically restricted positive selection. In this study, we estimated the selection intensity or selection coefficient of rs17822931-A in East Asians by analyzing two microsatellite loci flanking rs17822931 in the African (HapMap-YRI) and East Asian (HapMap-JPT and HapMap-CHB) populations. Assuming a recessive selection model, a coalescent-based simulation approach suggested that the selection coefficient of rs17822931-A had been approximately 0.01 in the East Asian population, and a simulation experiment using a pseudo-sampling variable revealed that the mutation of rs17822931-A occurred 2006 generations (95% credible interval, 1,023-3,901 generations) ago. In addition, we show that absolute latitude is significantly associated with the allele frequency of rs17822931-A in Asian, Native American, and European populations, implying that the selective advantage of rs17822931-A is related to an adaptation to a cold climate. Our results provide a striking example of how local adaptation has played a significant role in the diversification of human traits.     

Variable SNP density in aspartyl-protease genes of the malaria parasite Plasmodium falciparum

An analysis of the diversity of the aspartyl proteases of Plasmodium falciparum, known as plasmepsins (PMs), was completed in view of their possible role as drug targets. DNA sequence polymorphisms were identified in nine pm genes including their non-coding (introns and 5' flanking) sequences. All genes contained at least one single nucleotide polymorphism (SNP). Extensive microsatellite diversity was observed predominantly in non-coding sequences. All but one non-synonymous polymorphism (a conservative substitution) were mapped to the surface of the predicted protein, contradicting a possible role in enzymatic activity. The distribution of SNPs was found to be non-random among pm genes, with pm6 and pm10 having significantly higher SNP densities, suggesting they were under selection. For pm6 the majority of the SNPs were in introns and some of these may contribute to splice site variation. SNPs were found at a high density in both the coding and non-coding sequences of pm10. Recombination was important in generating additional diversity at this locus. Although direct selection for pm10 mutations could not be ruled out, the presence of balancing selection and a high density of SNPs in non-coding sequence led us to propose that another gene under selection may be influencing the diversity in the region. By sequencing short DNA tags in a 200 kb region flanking pm10 we show that a cluster of antigen genes, known to be under diversifying selection, may contribute to the observed diversity. We discuss the importance of diversity and local selection effects when choosing drug targets for intervention strategies.     

Delineation of the viral products of recombination in vaccinia virus-infected cells.

Plasmids containing the vaccinia virus thymidine kinase gene, its flanking DNA sequences, and the Escherichia coli beta-galactosidase gene were used in conjunction with a thymidine kinase-deficient virus to examine the viral products of recombination. Progeny derived from single-crossover events could be distinguished from those generated by gene conversion or double-crossover events when the beta-galactosidase gene was separated from the thymidine kinase gene by the flanking sequences. Using methotrexate to select for recombinant virus and a chromogenic indicator to detect beta-galactosidase, the generation of viral recombinants was measured over a 48-h period. Recombinant progeny were first observed at 12 h and increased to a maximum of 2.5% at 48 h. Single-crossover products, as determined by beta-galactosidase expression, reached a maximum of 57% of the recombinant population at 24 h and thereafter declined. DNA hybridization analysis was used to examine genomic structures of the progeny of the initial viral plaques, plaques purified three times, and those subject to a 10(4)-fold amplification. These analyses confirmed that single-crossover events within either the 5'- or 3'-homologous flanking sequences generated unstable recombinant structures. These structures were shown to contain a single copy of the intact thymidine kinase gene within the corresponding copy of the duplicated thymidine kinase flanking sequences, separated by the beta-galactosidase gene and plasmid DNA. Significantly, these duplicated structures could undergo further recombination to produce repeats of either the intact or the deleted thymidine kinase sequences. These intermediate structures ultimately degenerated to produce either the parental thymidine kinase-deleted or the wild-type genome. The wild-type genome was also shown to be generated directly by gene conversion or double-crossover events.