Excess of amino acid substitutions relative to polymorphism between X-linked duplications in Drosophila melanogaster

We have obtained sequence polymorphism data from 13 genes belonging to 5 gene families in Drosophila melanogaster where the K(a)/K(s) between copies is greater than 1. Twelve of these 13 loci are X-linked. In general, there is evidence of purifying selection in all families, as inferred both from levels of silent and replacement variation and insertion/deletion variation, suggesting that the loci are likely functional. Shared polymorphisms indicative of gene conversion between paralogs are rare among the X-linked families, in contrast to available data from autosomal duplicates. McDonald-Kreitman tests between duplicates reveal an excess of amino-acid fixations between copies in the X-linked families, suggesting that the divergence between these loci was driven by positive selection. In contrast, available data from autosomal duplicates show a deficit of fixations, consistent with gene conversion being a strong homogenizing force.     

The resolution of sexual antagonism by gene duplication

Disruptive selection between males and females can generate sexual antagonism, where alleles improving fitness in one sex reduce fitness in the other. This type of genetic conflict arises because males and females carry nearly identical sets of genes: opposing selection, followed by genetic mixing during reproduction, generates a population genetic "tug-of-war" that constrains adaptation in either sex. Recent verbal models suggest that gene duplication and sex-specific cooption of paralogs might resolve sexual antagonism and facilitate evolutionary divergence between the sexes. However, this intuitive proximal solution for sexual dimorphism potentially belies a complex interaction between mutation, genetic drift, and positive selection during duplicate fixation and sex-specific paralog differentiation. The interaction of these processes--within the explicit context of duplication and sexual antagonism--has yet to be formally described by population genetics theory. Here, we develop and analyze models of gene duplication and sex-specific differentiation between paralogs. We show that sexual antagonism can favor the fixation and maintenance of gene duplicates, eventually leading to the evolution of sexually dimorphic genetic architectures for male and female traits. The timescale for these evolutionary transitions is sensitive to a suite of genetic and demographic variables, including allelic dominance, recombination, sex linkage, and population size. Interestingly, we find that female-beneficial duplicates preferentially accumulate on the X chromosome, whereas male-beneficial duplicates are biased toward autosomes, independent of the dominance parameters of sexually antagonistic alleles. Although this result differs from previous models of sexual antagonism, it is consistent with several findings from the empirical genomics literature.     

The role of background selection in shaping patterns of molecular evolution and variation: evidence from variability on the Drosophila X chromosome

In the putatively ancestral population of Drosophila melanogaster, the ratio of silent DNA sequence diversity for X-linked loci to that for autosomal loci is approximately one, instead of the expected "null" value of 3/4. One possible explanation is that background selection (the hitchhiking effect of deleterious mutations) is more effective on the autosomes than on the X chromosome, because of the lack of crossing over in male Drosophila. The expected effects of background selection on neutral variability at sites in the middle of an X chromosome or an autosomal arm were calculated for different models of chromosome organization and methods of approximation, using current estimates of the deleterious mutation rate and distributions of the fitness effects of deleterious mutations. The robustness of the results to different distributions of fitness effects, dominance coefficients, mutation rates, mapping functions, and chromosome size was investigated. The predicted ratio of X-linked to autosomal variability is relatively insensitive to these variables, except for the mutation rate and map length. Provided that the deleterious mutation rate per genome is sufficiently large, it seems likely that background selection can account for the observed X to autosome ratio of variability in the ancestral population of D. melanogaster. The fact that this ratio is much less than one in D. pseudoobscura is also consistent with the model's predictions, since this species has a high rate of crossing over. The results suggest that background selection may play a major role in shaping patterns of molecular evolution and variation.     

Intron length evolution in Drosophila

I present data on the evolution of intron lengths among 3 closely related Drosophila species, D. melanogaster, Drosophila simulans, and Drosophila yakuba. Using D. yakuba as an outgroup, I mapped insertion and deletion mutations in 148 introns (spanning approximately 30 kb) to the D. melanogaster and D. simulans lineages. Intron length evolution in the 2 sister species has been different: in D. melanogaster, X-linked introns have increased slightly in size, whereas autosomal ones have decreased slightly in size; in D. simulans, both X-linked and autosomal introns have decreased in size. To understand the possible evolutionary causes of these lineage- and chromosome-specific patterns of intron evolution, I studied insertion-deletion (indel) polymorphism and divergence in D. melanogaster. Small insertion mutations segregate at elevated frequencies and enjoy elevated probabilities of fixation, particularly on the X chromosome. In contrast, there is no detectable X chromosome effect on fixations in D. simulans. These findings suggest X chromosome-specific selection or biased gene conversion-gap repair favoring insertions in D. melanogaster but not in D. simulans. These chromosome- and lineage-specific patterns of indel substitution are not easily explained by existing general population genetic models of intron length evolution. Genomic data from D. melanogaster further suggest that the forces described here affect introns and intergenic regions similarly.     

Recombination, dominance and selection on amino acid polymorphism in the Drosophila genome: contrasting patterns on the X and fourth chromosomes

Surveys of nucleotide polymorphism and divergence indicate that the average selection coefficient on Drosophila proteins is weakly positive. Similar surveys in mitochondrial genomes and in the selfing plant Arabidopsis show that weak negative selection has operated. These differences have been attributed to the low recombination environment of mtDNA and Arabidopsis that has hindered adaptive evolution through the interference effects of linkage. We test this hypothesis with new sequence surveys of proteins lying in low recombination regions of the Drosophila genome. We surveyed >3800 bp across four proteins at the tip of the X chromosome and >3600 bp across four proteins on the fourth chromosome in 24 strains of D. melanogaster and 5 strains of D. simulans. This design seeks to study the interaction of selection and linkage by comparing silent and replacement variation in semihaploid (X chromosome) and diploid (fourth chromosome) environments lying in regions of low recombination. While the data do indicate very low rates of exchange, all four gametic phases were observed both at the tip of the X and across the fourth chromosome. Silent variation is very low at the tip of the X (thetaS = 0.0015) and on the fourth chromosome (thetaS = 0.0002), but the tip of the X shows a greater proportional loss of variation than the fourth shows relative to normal-recombination regions. In contrast, replacement polymorphism at the tip of the X is not reduced (thetaR = 0.00065, very close to the X chromosome average). MK and HKA tests both indicate a significant excess of amino acid polymorphism at the tip of the X relative to the fourth. Selection is significantly negative at the tip of the X (Nes = -1.53) and nonsignificantly positive on the fourth (Nes approximately 2.9), analogous to the difference between mtDNA (or Arabidopsis) and the Drosophila genome average. Our distal X data are distinct from regions of normal recombination where the X shows a deficiency of amino acid polymorphism relative to the autosomes, suggesting more efficient selection against recessive deleterious replacement mutations. We suggest that the excess amino acid polymorphism on the distal X relative to the fourth chromosome is due to (1) differences in the mutation rate for selected mutations on the distal X or (2) a greater relaxation of selection from stronger linkage-related interference effects on the distal X. This relaxation of selection is presumed to be greater in magnitude than the difference in efficiency of selection between X-linked vs. autosomal selection.     

Reduced X-linked rare polymorphism in males in comparison to females of Drosophila melanogaster

Natural selection is assumed to act more strongly on X-linked loci than on autosomal loci because the fitness effect of a recessive mutation on the X chromosome is fully expressed in hemizygous males. Therefore, selection is expected to fix or remove recessive mutations on the X chromosome more efficiently than those on autosomes. However, the assumption that hemizygosity of the X chromosome selectively accelerates changes in allele frequency has not been confirmed directly. To examine this assumption, we investigated current natural selection on X-linked chemoreceptor genes in a natural population of Drosophila melanogaster by comparing nucleotide diversity, linkage disequilibrium (LD), and departure from the neutrality in 4 chemoreceptor genes on 100 X chromosomes each from female and male flies. The general pattern of nucleotide diversity and LD for the genes investigated was similar in females and males. In contrast, males harbored significantly fewer rare polymorphisms defined as singletons and doubletons. When all the gene sequences were concatenated, Tajima's D showed a significant departure from the neutrality in both females and males, whereas Fu and Li's F* value revealed departure only in males. These results suggest that some rare polymorphisms on the X chromosome from females are recessively deleterious and are removed by stronger purifying selection when transferred to hemizygous males.     

The X chromosome is a hot spot for sexually antagonistic fitness variation

Sexually antagonistic alleles are selected discordantly between the sexes. Experimental evidence indicates that sexually antagonistic fitness variation is abundant in the genome of Drosophila melanogaster. Theory predicts that the X chromosome will be enriched with this type of variation. To test this prediction in D. melanogaster, we sampled, and cytogenetically cloned, 20 X chromosomes and compared their fitness variation to genome-wide levels. At the juvenile stage, in which gender roles are most similar, the X chromosome made no detectable contribution to genome-wide fitness variation. At the adult stage, in which gender roles diverge, the X chromosome was estimated to harbour 45% of the genome-wide fitness variation and 97% of the genome-wide sexually antagonistic variation. This genomic structure has important implications for the process of sexual selection because X-linked sexually antagonistic variation contributes to negative intersexual heritability for fitness, i.e. high-fitness males (females) produce, on average, low-fitness daughters (sons).     

Adaptive protein evolution of X-linked and autosomal genes in Drosophila: implications for faster-X hypotheses

Patterns of sex chromosome and autosome evolution can be used to elucidate the underlying genetic basis of adaptative change. Evolutionary theory predicts that X-linked genes will adapt more rapidly than autosomes if adaptation is limited by the availability of beneficial mutations and if such mutations are recessive. In Drosophila, rates of molecular divergence between species appear to be equivalent between autosomes and the X chromosome. However, molecular divergence contrasts are difficult to interpret because they reflect a composite of adaptive and nonadaptive substitutions between species. Predictions based on faster-X theory also assume that selection is equally effective on the X and autosomes; this might not be true because the effective population sizes of X-linked and autosomal genes systematically differ. Here, population genetic and divergence data from Drosophila melanogaster, Drosophila simulans, and Drosophila yakuba are used to estimate the proportion of adaptive amino acid substitutions occurring in the D. melanogaster lineage. After gene composition and effective population size differences between chromosomes are controlled, X-linked and autosomal genes are shown to have equivalent rates of adaptive divergence with approximately 30% of amino acid substitutions driven by positive selection. The results suggest that adaptation is either unconstrained by a lack of beneficial genetic variation or that beneficial mutations are not recessive and are thus highly visible to natural selection whether on sex chromosomes or on autosomes.     

Evolutionary constraints and expression analysis of gene duplications in Rhodobacter Sphaeroides 2.4.1

ABSTRACT: BACKGROUND: Gene duplication is a major force that contributes to the evolution of new metabolic functions in all organisms. Rhodobacter sphaeroides 2.4.1 is a bacterium that displays a wide degree of metabolic versatility and genome complexity and therefore is a fitting model for the study of gene duplications in bacteria. A comprehensive analysis of 234 duplicate gene-pairs in R. sphaeroides was performed using structural constraint and expression analysis. RESULTS: The results revealed that most gene-pairs in in-paralogs are maintained under negative selection (omega [less than or equal to] 0.3), but the strength of selection differed among in-paralog gene-pairs. Although in-paralogs located on different replicons are maintained under purifying selection, the duplicated genes distributed between the primary chromosome (CI) and the second chromosome (CII) are relatively less selectively constrained than the gene-pairs located within each chromosome. The mRNA expression patterns of duplicate gene-pairs were examined through microarray analysis of this organism grown under seven different growth conditions. Results revealed that that ~62% of paralogs have similar expression patterns (cosine [greater than or equal to] 0.90) over all of these growth conditions, while only ~7% of paralogs are very different in their expression patterns (cosine < 0.50). CONCLUSIONS: The overall findings of the study suggest that only a small proportion of paralogs contribute to the metabolic diversity and the evolution of novel metabolic functions in R. sphaeroides. In addition, the lack of relationships between structural constraints and gene-pair expression suggests that patterns of gene-pair expression are likely associated with conservation or divergence of gene-pair promoter regions and other coregulation mechanisms.     

Evolutionary rearrangement of the amylase genomic regions between Drosophila melanogaster and Drosophila pseudoobscura

Two Drosophila pseudoobscura genomic clones have sequence similarity to the Drosophila melanogaster amylase region that maps to the 53CD region on the D. melanogaster cytogenetic map. The two clones with similarity to amylase map to sections 73A and 78C of the D. pseudoobscura third chromosome cytogenetic map. The complete sequences of both the 73A and 78C regions were compared to the D. melanogaster genome to determine if the coding region for amylase is present in both regions and to determine the evolutionary mechanism responsible for the observed distribution of the amylase gene or genes. The D. pseudoobscura 73A and 78C linkage groups are conserved with the D. melanogaster 41E and 53CD regions, respectively. The amylase gene, however, has not maintained its conserved linkage between the two species. These data indicate that amylase has moved via a transposition event in the D. melanogaster or D. pseudoobscura lineage. The predicted genes within the 73A and 78C regions show patterns of molecular evolution in synonymous and nonsynonymous sites that are consistent with previous studies of these two species.     

Comparative mapping of Z-orthologous genes in vertebrates: implications for the evolution of avian sex chromosomes

Sex chromosomes of birds and mammals are highly differentiated and share several cytological features. However, comparative gene mapping reveals extensive conserved synteny between the chicken Z sex chromosome and human chromosome 9 but not the human X sex chromosome, implying an independent origin of avian and mammalian sex chromosomes. To better understand the evolution of the avian Z chromosome we analysed the synteny of chicken Z-linked genes in zebrafish, which is the best-mapped teleost genome so far. Existing zebrafish maps do not support the existence of an ancestral Z linkage group in the zebrafish genome, whereas mammalian X-linked genes show at least some degree of synteny conservation. This is consistent with in situ hybridisation mapping data in the freshwater pufferfish, Tetraodon nigroviridis where mammalian X-linked genes show a much higher degree of conserved synteny than human chromosome 9 or the avian Z chromosome. Collectively, these data argue in favour of a more recent evolution of the avian Z chromosome, compared with the mammalian X.     

X-linked genes evolve higher codon bias in Drosophila and Caenorhabditis

Comparing patterns of molecular evolution between autosomes and sex chromosomes (such as X and W chromosomes) can provide insight into the forces underlying genome evolution. Here we investigate patterns of codon bias evolution on the X chromosome and autosomes in Drosophila and Caenorhabditis. We demonstrate that X-linked genes have significantly higher codon bias compared to autosomal genes in both Drosophila and Caenorhabditis. Furthermore, genes that become X-linked evolve higher codon bias gradually, over tens of millions of years. We provide several lines of evidence that this elevation in codon bias is due exclusively to their chromosomal location and not to any other property of X-linked genes. We present two possible explanations for these observations. One possibility is that natural selection is more efficient on the X chromosome due to effective haploidy of the X chromosomes in males and persistently low effective numbers of reproducing males compared to that of females. Alternatively, X-linked genes might experience stronger natural selection for higher codon bias as a result of maladaptive reduction of their dosage engendered by the loss of the Y-linked homologs.     

Early evolutionary history and genomic features of gene duplicates in the human genome

Background

Human gene duplicates have been the focus of intense research since the development of array-based and targeted next-generation sequencing approaches in the last decade. These studies have primarily concentrated on determining the extant copy-number variation from a population-genomic perspective but lack a robust evolutionary framework to elucidate the early structural and genomic characteristics of gene duplicates at emergence and their subsequent evolution with increasing age.

Results

We analyzed 184 gene duplicate pairs comprising small gene families in the draft human genome with 10 % or less synonymous sequence divergence. Human gene duplicates primarily originate from DNA-mediated events, taking up genomic residence as intrachromosomal copies in direct or inverse orientation. The distribution of paralogs on autosomes follows random expectations in contrast to their significant enrichment on the sex chromosomes. Furthermore, human gene duplicates exhibit a skewed gradient of distribution along the chromosomal length with significant clustering in pericentromeric regions. Surprisingly, despite the large average length of human genes, the majority of extant duplicates (83 %) are complete duplicates, wherein the entire ORF of the ancestral copy was duplicated. The preponderance of complete duplicates is in accord with an extremely large median duplication span of 36 kb, which enhances the probability of capturing ancestral ORFs in their entirety. With increasing evolutionary age, human paralogs exhibit declines in (i) the frequency of intrachromosomal paralogs, and (ii) the proportion of complete duplicates. These changes may reflect lower survival rates of certain classes of duplicates and/or the role of purifying selection. Duplications arising from RNA-mediated events comprise a small fraction (11.4 %) of all human paralogs and are more numerous in older evolutionary cohorts of duplicates.

Conclusions

The degree of structural resemblance, genomic location and duplication span appear to influence the long-term maintenance of paralogs in the human genome. The median duplication span in the human genome far exceeds that in C. elegans and yeast and likely contributes to the high prevalence of complete duplicates relative to structurally heterogeneous duplicates (partial and chimeric). The relative roles of regulatory sequence versus exon-intron structure changes in the acquisition of novel function by human paralogs remains to be determined.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1827-3) contains supplementary material, which is available to authorized users.     

A population genetics-phylogenetics approach to inferring natural selection in coding sequences

Through an analysis of polymorphism within and divergence between species, we can hope to learn about the distribution of selective effects of mutations in the genome, changes in the fitness landscape that occur over time, and the location of sites involved in key adaptations that distinguish modern-day species. We introduce a novel method for the analysis of variation in selection pressures within and between species, spatially along the genome and temporally between lineages. We model codon evolution explicitly using a joint population genetics-phylogenetics approach that we developed for the construction of multiallelic models with mutation, selection, and drift. Our approach has the advantage of performing direct inference on coding sequences, inferring ancestral states probabilistically, utilizing allele frequency information, and generalizing to multiple species. We use a Bayesian sliding window model for intragenic variation in selection coefficients that efficiently combines information across sites and captures spatial clustering within the genome. To demonstrate the utility of the method, we infer selective pressures acting in Drosophila melanogaster and D. simulans from polymorphism and divergence data for 100 X-linked coding regions.     

Selection, recombination and demographic history in Drosophila miranda

Selection, recombination, and the demographic history of a species can all have profound effects on genomewide patterns of variability. To assess the impact of these forces in the genome of Drosophila miranda, we examine polymorphism and divergence patterns at 62 loci scattered across the genome. In accordance with recent findings in D. melanogaster, we find that noncoding DNA generally evolves more slowly than synonymous sites, that the distribution of polymorphism frequencies in noncoding DNA is significantly skewed toward rare variants relative to synonymous sites, and that long introns evolve significantly slower than short introns or synonymous sites. These observations suggest that most noncoding DNA is functionally constrained and evolving under purifying selection. However, in contrast to findings in the D. melanogaster species group, we find little evidence of adaptive evolution acting on either coding or noncoding sequences in D. miranda. Levels of linkage disequilibrium (LD) in D. miranda are comparable to those observed in D. melanogaster, but vary considerably among chromosomes. These patterns suggest a significantly lower rate of recombination on autosomes, possibly due to the presence of polymorphic autosomal inversions and/or differences in chromosome sizes. All chromosomes show significant departures from the standard neutral model, including too much heterogeneity in synonymous site polymorphism relative to divergence among loci and a general excess of rare synonymous polymorphisms. These departures from neutral equilibrium expectations are discussed in the context of nonequilibrium models of demography and selection.     

Evolution of cis-regulatory elements in duplicated genes of yeast

An increasing number of studies report that functional divergence in duplicated genes is accompanied by gene expression changes, although the evolutionary mechanism behind this process remains unclear. Our genomic analysis on the yeast Saccharomyces cerevisiae shows that the number of shared regulatory motifs in the duplicates decreases with evolutionary time, whereas the total number of regulatory motifs remains unchanged. Moreover, genes with numerous paralogs in the yeast genome do not have especially low number of regulatory motifs. These findings indicate that degenerative complementation is not the sole mechanism behind expression divergence in yeast. Moreover, we found some evidence for the action of positive selection on cis-regulatory motifs after gene duplication. These results suggest that the evolution of functional novelty has a substantial role in yeast duplicate gene evolution.     

真核生物SMC基因家族中拷贝数目的长期稳定进化

Members of the Structural Maintenance of Chromosome (SMC) family have long been of interest to molecular and evolutionary biologists for their role in chromosome structural dynamics, particularly sister chromatid cohesion, condensation, and DNA repair. SMC and related proteins are found in all major groups of living organisms and share a common structure of conserved N and C globular domains separated from the conserved hinge domain by long coiled-coil regions. In eukaryotes there are six paralogous proteins that form three heterodimeric pairs, whereas in prokaryotes there is only one SMC protein that homodimerizes. From recently completed genome sequences, we have identified SMC genes from 34 eukaryotes that have not been described in previous reports. Our phylogenetic analysis of these and previously identified SMC genes supports an origin for the vertebrate meiotic SMC1 in the most recent common ancestor since the divergence from invertebrate animals. Additionally, we have identified duplicate copies due to segmental duplications for some of the SMC paralogs in plants and yeast, mainly SMC2 and SMC6, and detected evidence that duplicates of other paralogs were lost, suggesting differential evolution for these genes. Our analysis indicates that the SMC paralogs have been stably maintained at very low copy numbers, even after segmental (genome-wide) duplications. It is possible that such low copy numbers might be selected during eukaryotic evolution, although other possibilities are not ruled out.