Diverse selective modes among orthologs/paralogs of the chalcone synthase (Chs) gene family of Arabidopsis thaliana and its relative A. halleri ssp. gemmifera

As a model system, Arabidopsis thaliana and its wild relatives have played an important role in the study of genomics and evolution in plants. In this study, we examined the genetic diversity of the chalcone synthase (Chs) gene, which encodes a key enzyme of the flavonoid pathway and is located on chromosome five, as well as two Chs-like genes on the first and fourth chromosomes of Arabidopsis. The objectives of the study are to determine if natural selection operates differentially on the paralogs of the Chs gene family in A. thaliana and Arabidopsis halleri ssp. gemmifera. The mode of selection was inferred from Tajima's D values from noncoding and coding regions, as well as from the ratio of nonsynonymous to synonymous substitutions. Both McDonald-Kreitman and HKA tests revealed the effects of selection on the allelic distribution, except for the chromosome 1 paralog in ssp. gemmifera. The Chs gene on chromosome 5 was under purifying selection in both species. Significant, negative Tajima's D values at synonymous sites and positive Fay and Wu's H values within coding region, plus reduced genetic variability in introns, indicated effects of background selection in shaping the evolution of this gene region in A. thaliana. The Chs paralog on chromosome 1 was under positive selection in A. thaliana, while interspecific introgression and balancing selection determined the fates of the paralog and resulted in high heterogeneity in ssp. gemmifera. Local adaptation differentiated populations of Japan and China at the locus. In contrast, the other Chs-paralog of chromosome 4 was shaped by purifying selection in A. thaliana, while under positive selection in ssp. gemmifera, as indicated by dn/ds>1. Moreover, these contrasting patterns of selection have likely resulted in functional divergence in Arabidopsis, as indicated by radical amino acid substitutions at the chalcone synthase/stilbene synthase motif of the Chs genes. Unlike previous studies of the evolutionary history of A. thaliana, the high levels of genetic diversity in most gene regions of Chs paralogs and nonsignificant Tajima's D in the intron sequences of the Chs gene family in A. thaliana did not reflect the effects of a recent demographic expansion.     

Rapid divergence of gene duplicates on the Drosophila melanogaster X chromosome

The recent sequencing of several eukaryotic genomes has generated considerable interest in the study of gene duplication events. The classical model of duplicate gene evolution is that recurrent mutation ultimately results in one copy becoming a pseudogene, and only rarely will a beneficial new function evolve. Here, we study divergence between coding sequence duplications in Drosophila melanogaster as a function of the linkage relationship between paralogs. The mean K(a)/K(s) between all duplicates in the D. melanogaster genome is 0.2803, indicating that purifying selection is maintaining the structure of duplicate coding sequences. However, the mean K(a)/K(s) between duplicates that are both on the X chromosome is 0.4701, significantly higher than the genome average. Further, the distribution of K(a)/K(s) for these X-linked duplicates is significantly shifted toward higher values when compared with the distributions for paralogs in other linkage relationships. Two models of molecular evolution provide qualitative explanations of these observations-relaxation of selective pressure on the duplicate copies and, more likely, positive selection on recessive adaptations. We also show that there is an excess of X-linked duplicates with low K(s), suggesting a larger proportion of relatively young duplicates on the D. melanogaster X chromosome relative to autosomes.     

真核生物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.     

Spatial and temporal divergence of expression in duplicated barley germin-like protein-encoding genes

Subfunctionalization is the process by which a pair of duplicated genes, or paralogs, experiences a reduction of individual expression patterns or function while still reproducing the complete expression pattern and function of the ancestral gene. Two germin-like protein (GLP)-encoding genes, GerB and GerF, are paralogs that belong to a small gene family in barley (Hordeum vulgare). Both genes share high nucleotide sequence similarity in coding and noncoding regions and encode identical apoplastic proteins. The use of RNA gel blots, coupled with single-stranded conformation polymorphism (SSCP) analysis of RT-PCR products, elucidated the developmental and tissue-specific expression patterns of each gene. Individual expression patterns provided evidence of both overlapping redundancy and early subfunctionalization. GerB is predominantly expressed in developing shoots, while GerF is predominantly expressed in seedling roots, developing spikes, and pericarp/testa. GerF promoter deletion studies located a region (-356/-97) responsible for high promoter activity and showed the ability of GerB and GerF upstream regions to drive gfp expression in coleoptiles, epicarps, and lemma/palea of developing spikes. The observed expression patterns are consistent with proposed roles in plant development and defense mechanisms for this gene family. These roles may explain why redundancy has been selectively maintained in this duplicate gene pair.     

The monosaccharide transporter gene family in Arabidopsis and rice: a history of duplications, adaptive evolution, and functional divergence

Current hypotheses of gene duplicate divergence propose that surviving members of a gene duplicate pair may evolve, under conditions of purifying or nearly neutral selection, in one of two ways: with new function arising in one duplicate while the other retains original function (neofunctionalization [NF]) or partitioning of the original function between the 2 paralogs (subfunctionalization [SF]). More recent studies propose that SF followed by NF (subneofunctionalization [SNF]) explains the divergence of many duplicate genes. In this analysis, we evaluate these hypotheses in the context of the large monosaccharide transporter (MST) gene families in Arabidopsis and rice. MSTs have an ancient origin, predating plants, and have evolved in the seed plant lineage to comprise 7 subfamilies. In Arabidopsis, 53 putative MST genes have been identified, with one subfamily greatly expanded by tandem gene duplications. We searched the rice genome for members of the MST gene family and compared them with the MST gene family in Arabidopsis to determine subfamily expansion patterns and estimate gene duplicate divergence times. We tested hypotheses of gene duplicate divergence in 24 paralog pairs by comparing protein sequence divergence rates, estimating positive selection on codon sites, and analyzing tissue expression patterns. Results reveal the MST gene family to be significantly larger (65) in rice with 2 subfamilies greatly expanded by tandem duplications. Gene duplicate divergence time estimates indicate that early diversification of most subfamilies occurred in the Proterozoic (2500-540 Myr) and that expansion of large subfamilies continued through the Cenozoic (65-0 Myr). Two-thirds of paralog pairs show statistically symmetric rates of sequence evolution, most consistent with the SF model, with half of those showing evidence for positive selection in one or both genes. Among 8 paralog pairs showing asymmetric divergence rates, most consistent with the NF model, nearly half show evidence of positive selection. Positive selection does not appear in any duplicate pairs younger than approximately 34 Myr. Our data suggest that the NF, SF, and SNF models describe different outcomes along a continuum of divergence resulting from initial conditions of relaxed constraint after duplication.     

Coalescent processes and relaxation of selective constraints leading to contrasting genetic diversity at paralogs AtHVA22d and AtHVA22e in Arabidopsis thaliana

Duplicate loci offer a very powerful system for understanding the complicated genome structure and adaptive evolution of a gene family. In this study, the genetic variation at paralogs AtHVA22d and AtHVA22e, members of an ABA- and stress-inducible gene family, is examined in the selfing Arabidopsis thaliana. Population genetic analysis indicates contrasting levels of nucleotide diversity at overall exon sequence and nonsynonymous sites between AtHVA22d (pi = 0.00337, pi(rep) = 0.00158) and AtHVA22e (pi = 0.00054, pi(rep) = 0.00023). The fact of Ka/Ks ratios significantly less than 1 in all sequences indicates that both genes are functional and subjected to purifying selection. In addition, rooted at barley HVA22, accelerated evolution is detected at replacement changes in the AtHVA22d locus, indicating relaxation of purifying selection after gene duplication. However, relative rate tests reveal no deviation from the neutrality at synonymous sites between the two paralogs. Based on clock-like evolution, the rate of synonymous substitution is estimated at 1.83 x 10(-9) substitutions per site per year; and the divergence of the two paralogs is traced to 90 MYA, coinciding with a period of the diversification of angiosperms. Given no codon usage bias in both genes, natural selection alone cannot account for the 6.4-fold differences in the nucleotide variation at synonymous sites between the two paralogs. Random processes resulting in different coalescence times, 3.65 MYA at AtHVA22d vs. 1.20 MYA at AtHVA22e, may have predominantly contributed to the evident differences of the genetic diversity. Partially nonoverlapping modes of expression between the two functional paralogs suggest a subfunctionalization hypothesis for explaining the fates of duplicate loci.     

Evolution and Survival on Eutherian Sex Chromosomes

Since the two eutherian sex chromosomes diverged from an ancestral autosomal pair, the X has remained relatively gene-rich, while the Y has lost most of its genes through the accumulation of deleterious mutations in nonrecombining regions. Presently, it is unclear what is distinctive about genes that remain on the Y chromosome, when the sex chromosomes acquired their unique evolutionary rates, and whether X-Y gene divergence paralleled that of paralogs located on autosomes. To tackle these questions, here we juxtaposed the evolution of X and Y homologous genes (gametologs) in eutherian mammals with their autosomal orthologs in marsupial and monotreme mammals. We discovered that genes on the X and Y acquired distinct evolutionary rates immediately following the suppression of recombination between the two sex chromosomes. The Y-linked genes evolved at higher rates, while the X-linked genes maintained the lower evolutionary rates of the ancestral autosomal genes. These distinct rates have been maintained throughout the evolution of X and Y. Specifically, in humans, most X gametologs and, curiously, also most Y gametologs evolved under stronger purifying selection than similarly aged autosomal paralogs. Finally, after evaluating the current experimental data from the literature, we concluded that unique mRNA/protein expression patterns and functions acquired by Y (versus X) gametologs likely contributed to their retention. Our results also suggest that either the boundary between sex chromosome strata 3 and 4 should be shifted or that stratum 3 should be divided into two strata.     

Divergence of Gene Body DNA Methylation and Evolution of Plant Duplicate Genes

It has been shown that gene body DNA methylation is associated with gene expression. However, whether and how deviation of gene body DNA methylation between duplicate genes can influence their divergence remains largely unexplored. Here, we aim to elucidate the potential role of gene body DNA methylation in the fate of duplicate genes. We identified paralogous gene pairs from Arabidopsis and rice (Oryza sativa ssp. japonica) genomes and reprocessed their single-base resolution methylome data. We show that methylation in paralogous genes nonlinearly correlates with several gene properties including exon number/gene length, expression level and mutation rate. Further, we demonstrated that divergence of methylation level and pattern in paralogs indeed positively correlate with their sequence and expression divergences. This result held even after controlling for other confounding factors known to influence the divergence of paralogs. We observed that methylation level divergence might be more relevant to the expression divergence of paralogs than methylation pattern divergence. Finally, we explored the mechanisms that might give rise to the divergence of gene body methylation in paralogs. We found that exonic methylation divergence more closely correlates with expression divergence than intronic methylation divergence. We show that genomic environments (e.g., flanked by transposable elements and repetitive sequences) of paralogs generated by various duplication mechanisms are associated with the methylation divergence of paralogs. Overall, our results suggest that the changes in gene body DNA methylation could provide another avenue for duplicate genes to develop differential expression patterns and undergo different evolutionary fates in plant genomes.     

Evolution on the X chromosome: unusual patterns and processes

Although the X chromosome is usually similar to the autosomes in size and cytogenetic appearance, theoretical models predict that its hemizygosity in males may cause unusual patterns of evolution. The sequencing of several genomes has indeed revealed differences between the X chromosome and the autosomes in the rates of gene divergence, patterns of gene expression and rates of gene movement between chromosomes. A better understanding of these patterns should provide valuable information on the evolution of genes located on the X chromosome. It could also suggest solutions to more general problems in molecular evolution, such as detecting selection and estimating mutational effects on fitness.     

Faster evolution of Z-linked duplicate genes in chicken

It has been shown that duplicate genes on the X chromosome evolve much faster than duplicate genes on autosomes in Drosophila melanogaster.However,whether this phenomenon is general and can be applied to other species is not known.Here we examined this issue in chicken that have heterogametic females(females have ZW sex chromosome).We compared sequence divergence of duplicate genes on the Z chromosome with those on autosomes.We found that duplications on the Z chromosome indeed evolved faster than those on autosomes and show distinct patterns of molecular evolution from autosomal duplications.Examination of the expression of duplicate genes revealed an enrichment of duplications on the Z chromosome having male-biased expression and an enrichment of duplications on the autosomes having female-biased expression.These results suggest an evolutionary trend of the recruitment of duplicate genes towards reproduction-specific function.The faster evolution of duplications on Z than on the autosomes is most likely contributed by the selective forces driving the fixation of adaptive mutations on Z.Therefore,the common phenomena observed in both flies and chicken suggest that duplicate genes on sex chromosomes have distinct dynamics and are more influenced by natural selection than antosomal duplications,regardless of the kind of sex determination systems.     

Origins of new male germ-line functions from X-derived autosomal retrogenes in the mouse

Recent literature demonstrates that retrogenes tend to leave the X chromosome and integrate onto the autosomes and evolve male-biased expression patterns. Several selection-based evolutionary mechanisms have been proposed to explain this observation. Testing these selection-based models requires examining the evolutionary history and functional properties of new retrogenes, particularly those that show evidence of directional movement between the X and the autosomes (X-related retrogenes). This includes autosomal retrogenes with parental paralogs on the X chromosome (X-derived autosomal retrogenes) and those retrogenes integrated onto the X chromosomes (X-linked retrogenes). In order to understand why retrogenes tend to move nonrandomly in genomes, we examined the expression patterns and evolutionary mechanisms concerning gene pairs having young retrogenes--originating less than 20 MYA (after mouse-rat split). We demonstrate that these X-derived autosomal retrogenes evolved a more restricted male-biased expression pattern: they are expressed exclusively or predominantly in the testis, in particular, during the late stages of spermatogenesis. In contrast, the parental counterparts have relatively broad expression patterns in various tissues and spermatogenetic stages. We further observed that positive selection is targeting these X-derived autosomal retrogenes with novel male-biased expression patterns. This suggests that such retrogenes evolved new male germ-line functions that may be complementary to the functions of the parental paralogs, which themselves contribute little during spermatogenesis. Such evolutionary changes may be beneficial to the populations. Furthermore, most identified X-related retrogenes have recruited novel adjacent sequences as their untranslated regions (UTRs), suggesting that these UTRs, acquired de novo, may play an important role in establishing new regulatory mechanisms to carry out the new male germ-line functions.     

Differential selection after duplication in mammalian developmental genes

Gene duplication provides the opportunity for subsequent refinement of distinct functions of the duplicated copies. Either through changes in coding sequence or changes in regulatory regions, duplicate copies appear to obtain new or tissue-specific functions. If this divergence were driven by natural selection, we would expect duplicated copies to have differentiated patterns of substitutions. We tested this hypothesis using genes that duplicated before the human/mouse split and whose orthologous relations were clear. The null hypothesis is that the number of amino acid changes between humans and mice was distributed similarly across different paralogs. We used a method modified from Tang and Lewontin to detect heterogeneity in the amino acid substitution pattern between those different paralogs. Our results show that many of the paralogous gene pairs appear to be under differential selection in the human/mouse comparison. The properties that led to diversification appear to have arisen before the split of the human and mouse lineages. Further study of the diverged genes revealed insights regarding the patterns of amino acid substitution that resulted in differences in function and/or expression of these genes. This approach has utility in the study of newly identified members of gene families in genomewide data mining and for contrasting the merits of alternative hypotheses for the evolutionary divergence of function of duplicated genes.     

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.     

Complex signatures of selection and gene conversion in the duplicated globin genes of house mice

Results of electrophoretic surveys have suggested that hemoglobin polymorphism may be maintained by balancing selection in natural populations of house mice, Mus musculus. Here we report a survey of nucleotide variation in the adult globin genes of house mice from South America. We surveyed nucleotide polymorphism in two closely linked alpha-globin paralogs and two closely linked beta-globin paralogs to test whether patterns of variation are consistent with a model of long-term balancing selection. Surprisingly high levels of nucleotide polymorphism at the two beta-globin paralogs were attributable to the segregation of two highly divergent haplotypes, Hbbs (which carries two identical beta-globin paralogs) and Hbbd (which carries two functionally divergent beta-globin paralogs). Interparalog gene conversion on the Hbbs haplotype has produced a highly unusual situation in which the two paralogs are more similar to one another than either one is to its allelic counterpart on the Hbbd haplotype. Levels of nucleotide polymorphism and linkage disequilibrium at the two beta-globin paralogs suggest a complex history of diversity-enhancing selection that may be responsible for long-term maintenance of alternative protein alleles. The alternative two-locus beta-globin haplotypes are associated with pronounced differences in intraerythrocyte glutathione and nitric oxide metabolism, suggesting a possible mechanism for selection on hemoglobin function.     

Relaxed Selection Among Duplicate Floral Regulatory Genes in Lamiales

Polyploidization is a prevalent mode of genome diversification within plants. Most gene duplicates arising from polyploidization (paralogs) are typically lost, although a subset may be maintained under selection due to dosage, partitioning of gene function, or acquisition of novel functions. Because they experience selection in the presence of other duplicate loci across the genome, interactions among genes may also play a significant role in the maintenance of paralogs resulting from polyploidization. Previously, we identified duplicates of the genes LFY/FLO and AP3/DEF that directly interact in a floral regulatory pathway and are thought to be the result of ancient polyploidization in the Lamiales (> 50 mya). Although duplicates of MADS box genes including AP3/DEF are common throughout the angiosperm lineage, LFY/FLO duplicates in Lamiales are the first reported outside of tetraploid taxa. In order to explore hypotheses for the joint preservation of these interacting floral regulatory genes including novel LFY/FLO paralogs, here we clone FLO and DEF duplicates from additional Lamiales taxa and apply codon substitution models to test how selection acts on both genes following duplication. We find acceleration in the ratio of nonsynonymous-to-synonymous nucleotide substitutions for one (FLO) or both (DEF) paralogs that appears to be due to relaxed purifying selection as opposed to positive selection and shows a different pattern among functional domains of these genes. Several mechanisms are discussed that might be responsible for preservation of co-orthologs of FLO and DEF in Lamiales, including interactions among the genes of this regulatory pathway. Electronic Supplementary Material Electronic Supplementary material is available for this article at and accessible for authorised users. [Reviewing Editor: Dr. Yves Van de Peer]     

Gene Conversion Facilitates the Adaptive Evolution of Self-Resistance in Highly Toxic Newts

Reconstructing the histories of complex adaptations and identifying the evolutionary mechanisms underlying their origins are two of the primary goals of evolutionary biology. Taricha newts, which contain high concentrations of the deadly toxin tetrodotoxin (TTX) as an antipredator defense, have evolved resistance to self-intoxication, which is a complex adaptation requiring changes in six paralogs of the voltage-gated sodium channel (Na_v) gene family, the physiological target of TTX. Here, we reconstruct the origins of TTX self-resistance by sequencing the entire Na_v gene family in newts and related salamanders. We show that moderate TTX resistance evolved early in the salamander lineage in three of the six Na_v paralogs, preceding the proposed appearance of tetrodotoxic newts by ∼100 My. TTX-bearing newts possess additional unique substitutions across the entire Na_v gene family that provide physiological TTX resistance. These substitutions coincide with signatures of positive selection and relaxed purifying selection, as well as gene conversion events, that together likely facilitated their evolution. We also identify a novel exon duplication within Na_v1.4 encoding an expressed TTX-binding site. Two resistance-conferring changes within newts appear to have spread via nonallelic gene conversion: in one case, one codon was copied between paralogs, and in the second, multiple substitutions were homogenized between the duplicate exons of Na_v1.4. Our results demonstrate that gene conversion can accelerate the coordinated evolution of gene families in response to a common selection pressure.     

Probabilistic cross-species inference of orthologous genomic regions created by whole-genome duplication in yeast

Identification of orthologous genes across species becomes challenging in the presence of a whole-genome duplication (WGD). We present a probabilistic method for identifying orthologs that considers all possible orthology/paralogy assignments for a set of genomes with a shared WGD (here five yeast species). This approach allows us to estimate how confident we can be in the orthology assignments in each genomic region. Two inferences produced by this model are indicative of purifying selection acting to prevent duplicate gene loss. First, our model suggests that there are significant differences (up to a factor of seven) in duplicate gene half-life. Second, we observe differences between the genes that the model infers to have been lost soon after WGD and those lost more recently. Gene losses soon after WGD appear uncorrelated with gene expression level and knockout fitness defect. However, later losses are biased toward genes whose paralogs have high expression and large knockout fitness defects, as well as showing biases toward certain functional groups such as ribosomal proteins. We suggest that while duplicate copies of some genes may be lost neutrally after WGD, another set of genes may be initially preserved in duplicate by natural selection for reasons including dosage.