The role of gene conversion in determining sequence variation and divergence in the Est-5 gene family in Drosophila pseudoobscura

Purification and characterization of a specific late-larval esterase from two species of the Drosophila repleta group: contributions to understand its evolution

Background

After duplication, one copy of an original gene can become redundant and decay toward a pseudogene status or functionally diverge. Here, we performed the purification and biochemical characterization of EST-4 (a late larval β-esterase) from two Drosophila repleta group species, Drosophila mulleri and Drosophila arizonae, in order to establish comparative parameters between these enzymes in these species and to contribute to better understand their evolution.

Results

In D. mulleri, EST-4 had an optimal activity in temperatures ranging from 40° to 45°C and at pH 7.5, maintaining stability in alkaline pH (8.0 to 10.0). It was classified as serine esterase as its activity was inhibited by PMSF. No ion negatively modulated EST-4 activity, and iron had the most positive modulating effect. In D. arizonae, it showed similar optimum temperature (40°C), pH (8.0), and was also classified as a serine esterase, but the enzymatic stability was maintained in an acidic pH (5.5 to 6.5). Fe+2 had the opposite effect found in D. mulleri, that is, negative modulation. Al+3 almost totally inhibited the EST-4 activity, and Na+ and Cu+2 had a positive modulation effect. Kinetic studies, using ρ-nitrophenyl acetate as substrate, showed that EST-4 from D. mulleri had higher affinity, while in D. arizonae, it showed higher V max and catalytic efficiency in optimal reaction conditions.

Conclusions

EST-4 from D. mulleri and D. arizonae are very closely related and still maintain several similar features; however, they show some degree of differentiation. Considering that EST-4 from D. mulleri has more conspicuous gel mobility difference among all EST-4 studied so far and a lower catalytic efficiency was observed here, we proposed that after duplication, this new copy of the original gene became redundant and started to decay toward a pseudogene status in this species, which probably is not occurring in D. arizonae.

Computational methods of identification of pseudogenes based on functionality: entropy and GC content

Spectral entropy and GC content analyses reveal comprehensive structural features of DNA sequences. To illustrate the significance of these features, we analyze the β-esterase gene cluster, including the Est-6 gene and the ψEst-6 putative pseudogene, in seven species of the Drosophila melanogaster subgroup. The spectral entropies show distinctly lower structural ordering for ψEst-6 than for Est-6 in all species studied. However, entropy accumulation is not a completely random process for either gene and it shows to be nucleotide dependent. Furthermore, GC content in synonymous positions is uniformly higher in Est-6 than in ψEst-6, in agreement with the reduced GC content generally observed in pseudogenes and nonfunctional sequences. The observed differences in entropy and GC content reflect an evolutionary shift associated with the process of pseudogenization and subsequent functional divergence of ψEst-6 and Est-6 after the duplication event. The data obtained show the relevance and significance of entropy and GC content analyses for pseudogene identification and for the comparative study of gene-pseudogene evolution.

Gene duplication and ectopic gene conversion in Drosophila

The evolutionary impact of gene duplication events has been a theme of Drosophila genetics dating back to the Morgan School. While considerable attention has been placed on the genetic novelties that duplicates are capable of introducing, and the role that positive selection plays in their early stages of duplicate evolution, much less attention has been given to the potential consequences of ectopic (non-allelic) gene conversion on these evolutionary processes. In this paper we consider the historical origins of ectopic gene conversion models and present a synthesis of the current Drosophila data in light of several primary questions in the field.

Nonallelic gene conversion in the genus Drosophila

Nonallelic gene conversion has been proposed as a major force in homogenizing the sequences of paralogous genes. In this work, we investigate the extent and characteristics of gene conversion among gene families in nine species of the genus Drosophila. We carried out a genome-wide study of 2855 gene families (including 17,742 genes) and determined that conversion events involved 2628 genes. The proportion of converted genes ranged across species from 1 to 9% when paralogs of all ages were included. Although higher levels of gene conversion were found among young gene duplicates, at most 1-2% of the coding sequences of these duplicates were affected by conversion. Using a second approach relying on gene family size changes and gene-tree/species-tree reconciliation methods, we estimate that only 1-15% of gene trees are misled by gene conversion, depending on the lineage considered. Several features of paralogous genes correlate with gene conversion, such as intra-/interchromosomal location, level of nucleotide divergence, and GC content, although we found no definitive evidence for biased substitution patterns. After considering species-specific differences in the age and distance between paralogs, we found a highly significant difference in the amount of gene conversion among species. In particular, members of the melanogaster group showed the lowest proportion of converted genes. Our data therefore suggest underlying differences in the mechanistic basis of gene conversion among species.

Gene conversion and DNA sequence polymorphism in the sex-determination gene fog-2 and its paralog ftr-1 in Caenorhabditis elegans

Gene conversion, a form of concerted evolution, bears enormous potential to shape the trajectory of sequence and functional divergence of gene paralogs subsequent to duplication events. fog-2, a sex-determination gene unique to Caenorhabditis elegans and implicated in the origin of hermaphroditism in this species, resulted from the duplication of ftr-1, an upstream gene of unknown function. Synonymous sequence divergence in regions of fog-2 and ftr-1 (excluding recent gene conversion tracts) suggests that the duplication occurred 46 million generations ago. Gene conversion between fog-2 and ftr-1 was previously discovered in experimental fog-2 knockout lines of C. elegans, whereby hermaphroditism was restored in mutant obligately outcrossing male-female populations. We analyzed DNA-sequence variation in fog-2 and ftr-1 within 40 isolates of C. elegans from diverse geographic locations in order to evaluate the contribution of gene conversion to genetic variation in the two gene paralogs. The analysis shows that gene conversion contributes significantly to DNA-sequence diversity in fog-2 and ftr-1 (22% and 34%, respectively) and may have the potential to alter sexual phenotypes in natural populations. A radical amino acid change in a conserved region of the F-box domain of fog-2 was found in natural isolates of C. elegans with significantly lower fecundity. We hypothesize that the lowered fecundity is due to reduced masculinization and less sperm production and that amino acid replacement substitutions and gene conversion in fog-2 may contribute significantly to variation in the degree of inbreeding and outcrossing in natural populations.

Distinguishing among evolutionary models for the maintenance of gene duplicates

Determining the evolutionary forces responsible for the maintenance of gene duplicates is key to understanding the processes leading to evolutionary adaptation and novelty. In his highly prescient book, Susumu Ohno recognized that duplicate genes are fixed and maintained within a population with 3 distinct outcomes: neofunctionalization, subfunctionalization, and conservation of function. Subsequent researchers have proposed a multitude of population genetic models that lead to these outcomes, each differing largely in the role played by adaptive natural selection. In this paper, I present a nonmathematical review of these models, their predictions, and the evidence collected in support of each of them. Though the various outcomes of gene duplication are often strictly associated with the presence or absence of adaptive natural selection, I argue that determining the outcome of duplication is orthogonal to determining whether natural selection has acted. Despite an ever-growing field of research into the fate of gene duplicates, there is not yet clear evidence for the preponderance of one outcome over the others, much less evidence for the importance of adaptive or nonadaptive forces in maintaining these duplicates.

Birth and death of genes and functions in the beta-esterase cluster of Drosophila

Here we analyze the molecular evolution of the β-esterase gene cluster in the Drosophila genus using the recently released genome sequences of 12 Drosophila species. Molecular evolution in this small cluster is noteworthy because it contains contrasting examples of the types and stages of loss of gene function. Specifically, missing orthologs, pseudogenes, and null alleles are all inferred. Phylogenetic analyses also suggest a minimum of 9 gene gain–loss events; however, the exact number and age of these events is confounded by interparalog recombination. A previous enigma, in which allozyme loci were mapped to β-esterase genes that lacked catalytically essential amino acids, was resolved through the identification of neighbouring genes that contain the canonical catalytic residues and thus presumably encode the mapped allozymes. The originally identified genes are evolving with selective constraint, suggesting that they have a “noncatalytic” function. Curiously, 3 of the 4 paralogous β-esterase genes in the D. ananassae genome sequence have single inactivating (frame-shift or nonsense) mutations. To determine whether these putatively inactivating mutations were fixed, we sequenced other D. ananassae alleles of these four loci. We did not find any of the 3 inactivating mutations of the sequenced strain in 12 other strains; however, other inactivating mutations were observed in the same 3 genes. This is reminiscent of the high frequency of null alleles observed in one of the β-esterase genes (Est7/EstP) of D. melanogaster.

Entropy and GC Content in the beta-esterase gene cluster of the Drosophila melanogaster subgroup

We perform spectral entropy and GC content analyses in the beta-esterase gene cluster, including the Est-6 gene and the psiEst-6 putative pseudogene, in seven species of the Drosophila melanogaster species subgroup. psiEst-6 combines features of functional and nonfunctional genes. The spectral entropies show distinctly lower structural ordering for psiEst-6 than for Est-6 in all species studied. Our observations agree with previous results for D. melanogaster and provide additional support to our hypothesis that after the duplication event Est-6 retained the esterase-coding function and its role during copulation, while psiEst-6 lost that function but now operates in conjunction with Est-6 as an intergene. Entropy accumulation is not a completely random process for either gene. Structural entropy is nucleotide dependent. The relative normalized deviations for structural entropy are higher for G than for C nucleotides. The entropy values are similar for Est-6 and psiEst-6 in the case of A and T but are lower for Est-6 in the case of G and C. The GC content in synonymous positions is uniformly higher in Est-6 than in psiEst-6, which agrees with the reduced GC content generally observed in pseudogenes and nonfunctional sequences. The observed differences in entropy and GC content reflect an evolutionary shift associated with the process of pseudogenization and subsequent functional divergence of psiEst-6 and Est-6 after the duplication event.

Theories for analyzing polymorphism data in duplicated genes

A simple model for the evolutionary process of a pair of duplicated genes under concerted evolution is developed. The model considers mutation, recombination and gene conversion between two genes in a finite population. Based on diffusion theory, the expected amount of DNA variation within and between two genes are obtained. To investigate the pattern of DNA polymorphism, a coalescent tool to simulate patterns of polymorphism is developed. The theoretical results are well in agreement with polymorphism data in duplicated genes. The effect of selection on the pattern of polymorphism is also considered.

The effect of gene conversion on the divergence between duplicated genes

Nonindependent evolution of duplicated genes is called concerted evolution. In this article, we study the evolutionary process of duplicated regions that involves concerted evolution. The model incorporates mutation and gene conversion: the former increases d, the divergence between two duplicated regions, while the latter decreases d. It is demonstrated that the process consists of three phases. Phase I is the time until d reaches its equilibrium value, d(0). In phase II d fluctuates around d(0), and d increases again in phase III. Our simulation results demonstrate that the length of concerted evolution (i.e., phase II) is highly variable, while the lengths of the other two phases are relatively constant. It is also demonstrated that the length of phase II approximately follows an exponential distribution with mean tau, which is a function of many parameters including gene conversion rate and the length of gene conversion tract. On the basis of these findings, we obtain the probability distribution of the level of divergence between a pair of duplicated regions as a function of time, mutation rate, and tau. Finally, we discuss potential problems in genomic data analysis of duplicated genes when it is based on the molecular clock but concerted evolution is common.

Characterization of a mouse recombination hot spot locus encoding a novel non-protein-coding RNA

Our current knowledge of recombination hot spot activity in mammalian systems implicates a role for both the primary DNA sequence and the nature of the chromatin domain around it. In mice, the only recombination hot spots mapped to date have been confined to a cluster within the major histocompatibility complex (MHC) region. We present a high resolution analysis of a new recombination hot spot in the mouse genome which maps to mouse chromosome 8 C-D. Haplotype diversity analysis across 40 different strains of mice has enabled us to map recombination breakpoints to a 1-kb interval. This hot spot has a recombination intensity that is 10- to 100-fold above the genome average and has a mean gene conversion tract length of 371 bp. This meiotically active locus happens to be flanked by a transcribed region encoding a non-protein-coding RNA polymerase II transcript and the previously characterized repair site. Many of the primary DNA sequence features that have been reported for the mouse MHC hot spots are also shared by this hot spot locus and in addition, along with three other MHC hot spot loci, we show a new parallel feature of association of the crossover sites with the nuclear matrix.

Dynamics of a human interparalog gene conversion hotspot

Gene conversion between paralogs can alter their patterns of sequence identity, thus obscuring their evolutionary relationships and affecting their propensity to sponsor genomic rearrangements. The details of this important process are poorly understood in the human genome because allelic diversity complicates the interpretation of interparalog sequence differences. Here we exploit the haploid nature of the Y chromosome, which obviates complicating interallelic processes, together with its known phylogeny, to understand the dynamics of conversion between two directly repeated HERVs flanking the 780-kb AZFa region on Yq. Sequence analysis of a 787-bp segment of each of the HERVs in 36 Y chromosomes revealed one of the highest nucleotide diversities in the human genome, as well as evidence of a complex patchwork of highly directional gene conversion events. The rate of proximal-to-distal conversion events was estimated as 2.4 x 10(-4) to 1.2 x 10(-3) per generation (3.9 x 10(-7) to 1.9 x 10(-6) per base per generation), and the distal-to-proximal rate as about one-twentieth of this. Minimum observed conversion tract lengths ranged from 1 to 158 bp and maximum lengths from 19 to 1365 bp, with an estimated mean of 31 bp. Analysis of great ape homologs shows that conversion in this hotspot has a deep evolutionary history.

Pseudogenes: are they "junk" or functional DNA?

Pseudogenes have been defined as nonfunctional sequences of genomic DNA originally derived from functional genes. It is therefore assumed that all pseudogene mutations are selectively neutral and have equal probability to become fixed in the population. Rather, pseudogenes that have been suitably investigated often exhibit functional roles, such as gene expression, gene regulation, generation of genetic (antibody, antigenic, and other) diversity. Pseudogenes are involved in gene conversion or recombination with functional genes. Pseudogenes exhibit evolutionary conservation of gene sequence, reduced nucleotide variability, excess synonymous over nonsynonymous nucleotide polymorphism, and other features that are expected in genes or DNA sequences that have functional roles. We first review the Drosophila literature and then extend the discussion to the various functional features identified in the pseudogenes of other organisms. A pseudogene that has arisen by duplication or retroposition may, at first, not be subject to natural selection if the source gene remains functional. Mutant alleles that incorporate new functions may, nevertheless, be favored by natural selection and will have enhanced probability of becoming fixed in the population. We agree with the proposal that pseudogenes be considered as potogenes, i.e., DNA sequences with a potentiality for becoming new genes.

Molecular population genetics of theβ-esterase gene cluster ofDrosophila melanogaster

We have investigated nucleotide polymorphism at the beta-esterase gene cluster including the Est-6 gene and psiEst-6 putative pseudogene in four samples of Drosophila melanogaster derived from natural populations of southern Africa (Zimbabwe), Europe (Spain), North America (USA: California), and South America (Venezuela). A complex haplotype structure is revealed in both Est-6 and psiEst-6. Total nucleotide diversity is twice in psiEst-6 as in Est-6; diversity is higher in the African sample than in the non-African ones. Strong linkage disequilibrium occurs within the beta-esterase gene cluster in non-African samples, but not in the African one. Intragenic gene conversion events are detected within Est-6 and, to a much greater extent, within psiEst-6; intergenic gene conversion events are rare. Tests of neutrality with recombination are significant for the beta-esterase gene cluster in the non-African samples but not significant in the African one. We suggest that the demographic history (bottleneck and admixture of genetically differentiated populations) is the major factor shaping the pattern of nucleotide polymorphism in the beta-esterase gene cluster. However there are some 'footprints' of directional and balancing selection shaping specific distribution of nucleotide polymorphism within the cluster. Intergenic epistatic selection between Est-6 and psiEst-6 may play an important role in the evolution of the beta-esterase gene cluster preserving the putative pseudogene from degenerative destruction and reflecting possible functional interaction between the functional gene and the putative pseudogene. Est-6 and psiEst-6 may represent an indivisible intergenic complex ('intergene') in which each single component (Est-6 or psiEst-6) cannot separately carry out the full functional role.

Nuclear DNA analyses in genetic studies of populations: practice, problems and prospects

Population-genetic studies have been remarkably productive and successful in the last decade following the invention of PCR technology and the introduction of mitochondrial and microsatellite DNA markers. While mitochondrial DNA has proven powerful for genealogical and evolutionary studies of animal populations, and microsatellite sequences are the most revealing DNA markers available so far for inferring population structure and dynamics, they both have important and unavoidable limitations. To obtain a fuller picture of the history and evolutionary potential of populations, genealogical data from nuclear loci are essential, and the inclusion of other nuclear markers, i.e. single copy nuclear polymorphic (scnp) sequences, is clearly needed. Four major uncertainties for nuclear DNA analyses of populations have been facing us, i.e. the availability of scnp markers for carrying out such analysis, technical laboratory hurdles for resolving haplotypes, difficulty in data analysis because of recombination, low divergence levels and intraspecific multifurcation evolution, and the utility of scnp markers for addressing population-genetic questions. In this review, we discuss the availability of highly polymorphic single copy DNA in the nuclear genome, describe patterns and rate of evolution of nuclear sequences, summarize past empirical and theoretical efforts to recover and analyse data from scnp markers, and examine the difficulties, challenges and opportunities faced in such studies. We show that although challenges still exist, the above-mentioned obstacles are now being removed. Recent advances in technology and increases in statistical power provide the prospect of nuclear DNA analyses becoming routine practice, allowing allele-discriminating characterization of scnp loci and microsatellite loci. This certainly will increase our ability to address more complex questions, and thereby the sophistication of genetic analyses of populations.

The coalescent and infinite-site model of a small multigene family

The infinite-site model of a small multigene family with two duplicated genes is studied. The expectations of the amounts of nucleotide variation within and between two genes and linkage disequilibrium are obtained, and a coalescent-based method for simulating patterns of polymorphism in a small multigene family is developed. The pattern of DNA variation is much more complicated than that in a single-copy gene, which can be simulated by the standard coalescent. Using the coalescent simulation of duplicated genes, the applicability of statistical tests of neutrality to multigene families is considered.

Molecular population genetics of X-linked genes in Drosophila pseudoobscura

This article presents a nucleotide sequence analysis of 500 bp determined in each of five X-linked genes, runt, sisterlessA, period, esterase 5, and Heat-shock protein 83, in 40 Drosophila pseudoobscura strains collected from two populations. Estimates of the neutral migration parameter for the five loci show that gene flow among D. pseudoobscura populations is sufficient to homogenize inversion frequencies across the range of the species. Nucleotide diversity at each locus fails to reject a neutral model of molecular evolution. The sample of 40 chromosomes included six Sex-ratio inversions, a series of three nonoverlapping inversions that are associated with a strong meiotic drive phenotype. The selection driven by the Sex-ratio meiotic drive element has not fixed variation across the X chromosome of D. pseudoobscura because, while significant linkage disequilibrium was observed within the sisterlessA, period, and esterase 5 genes, we did not find evidence for nonrandom association among loci. The Sex-ratio chromosome was estimated to be 25,000 years old based on the decomposition of linkage disequilibrium between esterase 5 and Heat-shock protein 83 or 1 million years old based on the net divergence of esterase 5 between Standard and Sex-ratio chromosomes. Genetic diversity was depressed within esterase 5 within Sex-ratio chromosomes, while the four other genes failed to show a reduction in heterozygosity in the Sex-ratio background. The reduced heterogeneity in esterase 5 is due either to its location near one of the Sex-ratio inversion breakpoints or that it is closely linked to a gene or genes responsible for the Sex-ratio meiotic drive system.

Contrasting patterns of nonneutral evolution in proteins encoded in nuclear and mitochondrial genomes

We report that patterns of nonneutral DNA sequence evolution among published nuclear and mitochondrially encoded protein-coding loci differ significantly in animals. Whereas an apparent excess of amino acid polymorphism is seen in most (25/31) mitochondrial genes, this pattern is seen in fewer than half (15/36) of the nuclear data sets. This differentiation is even greater among data sets with significant departures from neutrality (14/15 vs. 1/6). Using forward simulations, we examined patterns of nonneutral evolution using parameters chosen to mimic the differences between mitochondrial and nuclear genetics (we varied recombination rate, population size, mutation rate, selective dominance, and intensity of germ line bottleneck). Patterns of evolution were correlated only with effective population size and strength of selection, and no single genetic factor explains the empirical contrast in patterns. We further report that in Arabidopsis thaliana, a highly self-fertilizing plant with effectively low recombination, five of six published nuclear data sets also exhibit an excess of amino acid polymorphism. We suggest that the contrast between nuclear and mitochondrial nonneutrality in animals stems from differences in rates of recombination in conjunction with a distribution of selective effects. If the majority of mutations segregating in populations are deleterious, high linkage may hinder the spread of the occasional beneficial mutation.

The expected number of alleles in a gene conversion model with mutation

Gene conversions and mutations affect the number of different alleles in a population as it evolves. Here these effects are studied for a sample of r individuals represented by one chromosome, each with n loci. The long-term evolution of the sample is modelled with the coalescent. Gene conversions can only occur within a chromosome between a pair of loci. Recombination is not included in the model. Mutations occur according to the infinitely many alleles model. A new expression for the expected number of alleles is derived using the Backward time approach. It is then possible to simulate the expectation for various sample paths. The new estimator is computationally inexpensive and has a lower variance than the conventional, forward time, approximation by simulation.

A note on the correlation between heterozygosity and recombination rate

Mechanisms responsible for the correlation between heterozygosity and recombination rate have been thought to be background selection and selective sweep. In addition to these two, diversity enhancing selection and/or gene conversion is suggested to be contributing to the correlation, by surveying some examples of highly polymorphic loci that tend to locate in high-recombination regions.

Within- and between-species DNA sequence variation and the 'footprint' of natural selection

Extensive DNA data emerging from genome-sequencing projects have revitalized interest in the mechanisms of molecular evolution. Although the contribution of natural selection at the molecular level has been debated for over 30 years, the relevant data and appropriate statistical methods to address this issue have only begun to emerge. This paper will first present the predominant models of neutral, nearly neutral, and adaptive molecular evolution. Then, a method to identify the role of natural selection in molecular evolution by comparing within- and between-species DNA sequence variation will be presented. Computer simulations show that such methods are powerful for detecting even very weak selection. Examination of DNA variation data within and between Drosophila species suggests that 'silent' sites evolve under a balance between weak selection and genetic drift. Simulated data also show that sequence comparisons are a powerful method to detect adaptive protein evolution, even when selection is weak or affects a small fraction of nucleotide sites. In the Drosophila data examined, positive selection appears to be a predominant force in protein evolution.

Inferring the fitness effects of DNA mutations from polymorphism and divergence data: statistical power to detect directional selection under stationarity and free recombination

The fitness effects of classes of DNA mutations can be inferred from patterns of nucleotide variation. A number of studies have attributed differences in levels of polymorphism and divergence between silent and replacement mutations to the action of natural selection. Here, I investigate the statistical power to detect directional selection through contrasts of DNA variation among functional categories of mutations. A variety of statistical approaches are applied to DNA data simulated under Sawyer and Hartl's Poisson random field model. Under assumptions of free recombination and stationarity, comparisons that include both the frequency distributions of mutations segregating within populations and the numbers of mutations fixed between populations have substantial power to detect even very weak selection. Frequency distribution and divergence tests are applied to silent and replacement mutations among five alleles of each of eight Drosophila simulans genes. Putatively "preferred" silent mutations segregate at higher frequencies and are more often fixed between species than "unpreferred" silent changes, suggesting fitness differences among synonymous codons. Amino acid changes tend to be either rare polymorphisms or fixed differences, consistent with a combination of deleterious and adaptive protein evolution. In these data, a substantial fraction of both silent and replacement DNA mutations appear to affect fitness.