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.     

Adaptive evolution of genes underlying schizophrenia

Schizophrenia poses an evolutionary-genetic paradox because it exhibits strongly negative fitness effects and high heritability, yet it persists at a prevalence of approximately 1% across all human cultures. Recent theory has proposed a resolution: that genetic liability to schizophrenia has evolved as a secondary consequence of selection for human cognitive traits. This hypothesis predicts that genes increasing the risk of this disorder have been subject to positive selection in the evolutionary history of humans and other primates. We evaluated this prediction using tests for recent selective sweeps in human populations and maximum-likelihood tests for selection during primate evolution. Significant evidence for positive selection was evident using one or both methods for 28 of 76 genes demonstrated to mediate liability to schizophrenia, including DISC1, DTNBP1 and NRG1, which exhibit especially strong and well-replicated functional and genetic links to this disorder. Strong evidence of non-neutral, accelerated evolution was found for DISC1, particularly for exon 2, the only coding region within the schizophrenia-associated haplotype. Additionally, genes associated with schizophrenia exhibited a statistically significant enrichment in their signals of positive selection in HapMap and PAML analyses of evolution along the human lineage, when compared with a control set of genes involved in neuronal activities. The selective forces underlying adaptive evolution of these genes remain largely unknown, but these findings provide convergent evidence consistent with the hypothesis that schizophrenia represents, in part, a maladaptive by-product of adaptive changes during human evolution.     

Adaptive evolution in the photosensory domain of phytochrome A in early angiosperms

Flowering plant diversity now far exceeds the combined diversity of all other plant groups. Recently identified extant remnants of the earliest-diverging lines suggest that the first angiosperms may have lived in shady, disturbed, and moist understory habitats, and that the aquatic habit also arose early. This would have required the capacity to begin life in dimly lit environments. If so, evolution in light-sensing mechanisms may have been crucial to their success. The photoreceptor phytochrome A is unique among angiosperm phytochromes in its capacity to serve a transient role under conditions where an extremely high sensitivity is required. We present evidence of altered functional constraints between phytochrome A (PHYA) and its paralog, PHYC. Tests for selection suggest that an elevation in nonsynonymous rates resulted from an episode of selection along the branch leading to all angiosperm PHYA sequences. Most nucleotide sites (95%) are selectively constrained, and the ratio of nonsynonymous to synonymous substitutions on branches within the PHYA clade does not differ from the ratio on the branches in the PHYC clade. Thus, positive selection at a handful of sites, rather than relaxation of selective constraints, apparently has played a major role in the evolution of the photosensory domain of phytochrome A. The episode of selection occurred very early in the history of flowering plants, suggesting that innovation in phyA may have given the first angiosperms some adaptive advantage.     

Cytogenetic mapping of the W chromosome in the genus Leporinus (Teleostei, Anostomidae) using a highly repetitive DNA sequence

The distribution of the Leporinus elongatus LeSpeI repetitive sequence in other Leporinus species was studied in an attempt to elucidate the evolutionary history of sex chromosomes in this genus using chromosome fluorescence in situ hybridization. The presence of fluorescent signals only in species that have differentiated sex chromosomes suggests that this sequence is related to the differentiation of sex chromosomes in this genus. Thus, these data will contribute to a better understanding of chromosome evolution, especially for sex chromosomes, in the Leporinus genus.     

Adaptive evolution of a key gene affecting queen and worker traits in the honey bee, Apis mellifera

The vitellogenin egg yolk precursor protein represents a well-studied case of social pleiotropy in the model organism Apis mellifera. Vitellogenin is associated with fecundity in queens and plays a major role in controlling division of labour in workers, thereby affecting both individual and colony-level fitness. We studied the molecular evolution of vitellogenin and seven other genes sequenced in a large population panel of Apis mellifera and several closely related species to investigate the role of social pleiotropy on adaptive protein evolution. We found a significant excess of nonsynonymous fixed differences between A. mellifera, A. cerana and A. florea relative to synonymous sites indicating high rates of adaptive evolution at vitellogenin. Indeed, 88% of amino acid changes were fixed by selection in some portions of the gene. Further, vitellogenin exhibited hallmark signatures of selective sweeps in A. mellifera, including a significant skew in the allele frequency spectrum, extreme levels of genetic differentiation and linkage disequilibrium. Finally, replacement polymorphisms in vitellogenin were significantly enriched in parts of the protein involved in binding lipid, establishing a link between the gene's structure, function and effects on fitness. Our case study provides unequivocal evidence of historical and ongoing bouts of adaptive evolution acting on a key socially pleiotropic gene in the honey bee.     

Duplication of accelerated evolution and growth hormone gene in passerine birds

We report the discovery of a duplication of the growth hormone (GH) gene in a major group of birds, the passerines (Aves: Passeriformes). Phylogenetic analysis of 1.3-kb partial DNA sequences of GH genes for 24 species of passerines and numerous outgroups indicates that the duplication occurred in the ancestral lineage of extant passerines. Both duplicates and their open-reading frames are preserved throughout the passerine clade, and both duplicates are expressed in the zebra finch brain, suggesting that both are likely to be functional. The estimated rates of amino acid evolution are more than 10-fold higher in passerine GH genes than in those of their closest nonpasserine relatives. In addition, although the 84 codons sequenced are generally highly conserved for both passerines and nonpasserines, comparisons of the nonsynonymous/synonymous substitution ratios and the rate of predicted amino acid changes indicate that the 2 gene duplicates are evolving under different selective pressures and may be functionally divergent. The evidence of differential selection, coupled with the preservation of both gene copies in all major lineages since the origin of passerines, suggests that the duplication may be of adaptive significance, with possible implications for the explosive diversification of the passerine clade.     

Adaptive protein evolution and regulatory divergence in Drosophila

Two recent studies demonstrated a positive correlation between divergence in gene expression and protein sequence in Drosophila. This correlation could be driven by positive selection or variation in functional constraint. To distinguish between these alternatives, we compared patterns of molecular evolution for 1,862 genes with two previously reported estimates of expression divergence in Drosophila. We found a slight negative trend (nonsignificant) between positive selection on protein sequence and divergence in expression levels between Drosophila melanogaster and Drosophila simulans. Conversely, shifts in expression patterns during Drosophila development showed a positive association with adaptive protein evolution, though as before the relationship was weak and not significant. Overall, we found no strong evidence for an increase in the incidence of positive selection on protein-coding regions in genes with divergent expression in Drosophila, suggesting that the previously reported positive association between protein and regulatory divergence primarily reflects variation in functional constraint.     

Genomic evidence for a large-Z effect

The 'large-X effect' suggests that sex chromosomes play a disproportionate role in adaptive evolution. Theoretical work indicates that this effect may be most pronounced in genetic systems with female heterogamety under both good-genes and Fisher's runaway models of sexual selection (males ZZ, females ZW). Here, I use a comparative genomic approach (alignments of several thousands of chicken-zebra finch-human-mouse-opossum orthologues) to show that avian Z-linked genes are highly overrepresented among those bird-mammalian orthologues that show evidence of accelerated rate of functional evolution in birds relative to mammals; the data suggest a twofold excess of such genes on the Z chromosome. A reciprocal analysis of genes accelerated in mammals found no evidence for an excess of X-linkage. This would be compatible with theoretical expectations for differential selection on sex-linked genes under male and female heterogamety, although the power in this case was not sufficient to statistically show that 'large-Z' was more pronounced than 'large-X'. Accelerated Z-linked genes include a variety of functional categories and are characterized by higher non-synonymous to synonymous substitution rate ratios than both accelerated autosomal and non-accelerated genes. This points at a genomic 'large-Z effect', which is widespread and of general significance for adaptive divergence in birds.     

Evolution of the Spindlin Gene in Birds: Independent Cessation of the Recombination of Sex Chromosomes at the Spindlin Locus in Neognathous Birds and Tinamous, a Palaeognathous Avian Family

Tinamous (Aves, Palaeognathae, Tinamiformes) are primitive birds, generally considered to be the sister group to the ratites. Tinamous possess a W sex-chromosome, intermediate in heterochromatization between the largely euchromatic W chromosome of the ratites and the highly condensed W chromosome of the neognathous birds. Of the four genes which are known to have diverged copies on the neognathous avian W and Z chromosome (ATP5A1, CHD1, PKC and SPIN) only the spindlin gene has W- and Z-chromosomal forms in the tinamiformes. This paper describes experiments which show that the sequences of these forms are more similar to each other and to the homologous undifferentiated spindlin gene sequences in the ratite genome than to the W or Z forms of the spindlin gene in other, neognathous species. This suggests that cessation of recombination at the spindlin locus of the ancestral W and Z chromosomes of the paleognathous tinamiformes and the neognathous avian species were independent events.     

Male-biased mutation, sex linkage, and the rate of adaptive evolution

An interaction between sex-linked inheritance and sex-biased mutation rates may affect the rate of adaptive evolution. Males have much higher mutation rates than females in several vertebrate and plant taxa. When evolutionary rates are limited by the supply of favorable new mutations, then genes will evolve faster when located on sex chromosomes that spend more time in males. For mutations with additive effects, Y-linked genes evolve fastest, followed by Z-linked genes, autosomal genes, X-linked genes, and finally W-linked and cytoplasmic genes. This ordering can change when mutations show dominance. The predicted differences in substitution rates may be detectable at the molecular level. Male-biased mutation could cause adaptive changes to accumulate more readily on certain kinds of chromosomes and favor animals with Z-W sex determination to have rapidly evolving male sexual displays.     

HINTW,a W‐chromosome HINT gene in chick,is expressed ubiquitously and is a robust female cell marker applicable in intraspecific chimera studies

Grafting and transplantation experiments in embryology require proper distinction between host and donor tissues. For the avian model this has traditionally been achieved by using two closely related species (e.g., chick and quail) followed by species‐specific antibody staining. Here, we show that an in situ hybridization probe against the HINTW gene is a robust and reliable marker for female‐derived chicken cells. At all pre‐circulation stages tested, all cells in female embryos, independently confirmed by PCR analysis, were strongly positive for HINTW, whereas all male embryos were negative. This probe is broadly applicable in intra‐specific chick/chick chimera studies, and as a proof of principle, we utilized this probe to detect female cells in three experimental settings: (1) to mark female donor cells in a node transplantation assay; (2) to distinguish female cells in male/female twins generated by the Cornish pasty culture; and (3) to detect female half of the embryo in artificially generated bilateral gynandromorphs. A rapid, PCR based pre‐screening step increases the efficiency of obtaining desired donor/host sex combination from 25% to 100%. For most avian chimera studies, this female‐specific in situ probe is a low cost alternative to the commonly used QCPN antibody and to ubiquitous‐GFP chicken strains which are not widely available to the research community. genesis 52:424–430, 2014. © 2014 Wiley Periodicals, Inc.     

Molecular evolution of the avian growth hormone gene and comparison with its mammalian counterpart

The molecular evolution of all available avian growth hormone (GH) gene sequences was investigated using both maximum-likelihood and parsimony methods, and the patterns compared to those found in mammals. In contrast to the rapid bursts of evolution observed for mammalian GH, the evolutionary rate of the avian GH mature peptide appears to have been more constant. However several positively selected sites were identified at functionally important positions in the avian signal peptide by the site-specific likelihood method. This implies that sequence variation in the avian GH signal peptide may be adaptive, although more conservative parsimony methods failed to confirm this. Nevertheless, the differing patterns of avian and mammalian GH signal peptide molecular evolution are consistent with the apparently differing roles of GH in controlling growth in these taxonomic groups and support the hypothesis that signal peptide sequence variation may in fact be the basis for increased functional complexity.     

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.     

Origin of a Giant Sex Chromosome

Chromosome size and morphology vary within and among species, but little is known about the proximate or ultimate causes of these differences. Cichlid fish species in the tribe Oreochromini share an unusual giant chromosome that is ∼3 times longer than the other chromosomes. This giant chromosome functions as a sex chromosome in some of these species. We test two hypotheses of how this giant sex chromosome may have evolved. The first hypothesis proposes that it evolved by accumulating repetitive elements as recombination was reduced around a dominant sex determination locus, as suggested by canonical models of sex chromosome evolution. An alternative hypothesis is that the giant sex chromosome originated via the fusion of an autosome with a highly repetitive B chromosome, one of which carried a sex determination locus. We test these hypotheses using comparative analysis of chromosome-scale cichlid and teleost genomes. We find that the giant sex chromosome consists of three distinct regions based on patterns of recombination, gene and transposable element content, and synteny to the ancestral autosome. The WZ sex determination locus encompasses the last ∼105 Mb of the 134-Mb giant chromosome. The last 47 Mb of the giant chromosome shares no obvious homology to any ancestral chromosome. Comparisons across 69 teleost genomes reveal that the giant sex chromosome contains unparalleled amounts of endogenous retroviral elements, immunoglobulin genes, and long noncoding RNAs. The results favor the B chromosome fusion hypothesis for the origin of the giant chromosome.     

Adaptive evolution of metabolic pathways in Drosophila

The adaptive significance of enzyme variation has been of central interest in population genetics. Yet, how natural selection operates on enzymes in the larger context of biochemical pathways has not been broadly explored. A basic expectation is that natural selection on metabolic phenotypes will target enzymes that control metabolic flux, but how adaptive variation is distributed among enzymes in metabolic networks is poorly understood. Here, we use population genetic methods to identify enzymes responding to adaptive selection in the pathways of central metabolism in Drosophila melanogaster and Drosophila simulans. We report polymorphism and divergence data for 17 genes that encode enzymes of 5 metabolic pathways that converge at glucose-6-phosphate (G6P). Deviations from neutral expectations were observed at five loci. Of the 10 genes that encode the enzymes of glycolysis, only aldolase (Ald) deviated from neutrality. The other 4 genes that were inconsistent with neutral evolution (glucose-6-phosphate dehydrogenase [G6pd]), phosphoglucomutase [Pgm], trehalose-6-phosphate synthetase [Tps1], and glucose-6phosphatase [G6pase] encode G6P branch point enzymes that catalyze reactions at the entry point to the pentose-phosphate, glycogenic, trehalose synthesis, and gluconeogenic pathways. We reconcile these results with population genetics theory and existing arguments on metabolic regulation and propose that the incidence of adaptive selection in this system is related to the distribution of flux control. The data suggest that adaptive evolution of G6P branch point enzymes may have special significance in metabolic adaptation.     

Adaptive molecular evolution of MC1R gene reveals the evidence for positive diversifying selection in indigenous goat populations

Detecting signatures of selection can provide a new insight into the mechanism of contemporary breeding and artificial selection and further reveal the causal genes associated to the phenotypic variation. However, the signatures of selection on genes entailing for profitable traits between Chinese commercial and indigenous goats have been poorly interpreted. We noticed footprints of positive selection at MC 1R gene containing SNP s genotyped in five Chinese native goat breeds. An experimental distribution of F _ST was built based on approximations of F _ST for each SNP across five breeds. We identified selection using the high F _ST outlier method and found that MC 1R candidate gene show evidence of positive selection. Furthermore, adaptive selection pressure on specific codons was determined using different codon based on maximum‐likelihood methods; signature of positive selection in mammalian MC 1R was explored in individual codons. Evolutionary analyses were inferred under maximum likelihood models, the HyPhy package implemented in the DATAMONKEY Web Server. The results of codon selection displayed positive diversifying selection at the sites were mainly involved in development of genetic variations in coat color in various mammalian species. Positive diversifying selection inferred with recent evolutionary changes in domesticated goat MC 1R provides new insights that the gene evolution may have been modulated by domestication events in goats.     

Evaluation of whether accelerated protein evolution in chordates has occurred before, after, or simultaneously with gene duplication

Gene duplication and loss are predicted to be at least of the order of the substitution rate and are key contributors to the development of novel gene function and overall genome evolution. Although it has been established that proteins evolve more rapidly after gene duplication, we were interested in testing to what extent this reflects causation or association. Therefore, we investigated the rate of evolution prior to gene duplication in chordates. Two patterns emerged; firstly, branches, which are both preceded by a duplication and followed by a duplication, display an elevated rate of amino acid replacement. This is reflected in the ratio of nonsynonymous to synonymous substitution (mean nonsynonymous to synonymous nucleotide substitution rate ratio [Ka:Ks]) of 0.44 compared with branches preceded by and followed by a speciation (mean Ka:Ks of 0.23). The observed patterns suggest that there can be simultaneous alteration in the selection pressures on both gene duplication and amino acid replacement, which may be consistent with co-occurring increases in positive selection, or alternatively with concurrent relaxation of purifying selection. The pattern is largely, but perhaps not completely, explained by the existence of certain families that have elevated rates of both gene duplication and amino acid replacement. Secondly, we observed accelerated amino acid replacement prior to duplication (mean Ka:Ks for postspeciation preduplication branches was 0.27). In some cases, this could reflect adaptive changes in protein function precipitating a gene duplication event. In conclusion, the circumstances surrounding the birth of new proteins may frequently involve a simultaneous change in selection pressures on both gene-copy number and amino acid replacement. More precise modeling of the relative importance of preduplication, postduplication, and simultaneous amino acid replacement will require larger and denser genomic data sets from multiple species, allowing simultaneous estimation of lineage-specific fluctuations in mutation rates and adaptive constraints.     

Sexually dimorphic gene expression and transcriptome evolution provide mixed evidence for a fast‐Z effect in Heliconius

Sex chromosomes have different evolutionary properties compared to autosomes due to their hemizygous nature. In particular, recessive mutations are more readily exposed to selection, which can lead to faster rates of molecular evolution. Here, we report patterns of gene expression and molecular evolution for a group of butterflies. First, we improve the completeness of the Heliconius melpomene reference annotation, a neotropical butterfly with a ZW sex determination system. Then, we analyse RNA from male and female whole abdomens and sequence female ovary and gut tissue to identify sex‐ and tissue‐specific gene expression profiles in H. melpomene. Using these expression profiles, we compare (a) sequence divergence and polymorphism; (b) the strength of positive and negative selection; and (c) rates of adaptive evolution, for Z and autosomal genes between two species of Heliconius butterflies, H. melpomene and H. erato. We show that the rate of adaptive substitutions is higher for Z than autosomal genes, but contrary to expectation, it is also higher for male‐biased than female‐biased genes. Additionally, we find no significant increase in the rate of adaptive evolution or purifying selection on genes expressed in ovary tissue, a heterogametic‐specific tissue. Our results contribute to a growing body of literature from other ZW systems that also provide mixed evidence for a fast‐Z effect where hemizygosity influences the rate of adaptive substitutions.     

Adaptive evolution of the insulin gene in caviomorph rodents

Insulin is a conservative molecule among mammals, maintaining both its structure and function. Rodents that belong to the Suborder Hystricognathi represent an exception, having a very divergent molecule with unusual physiological properties. In this work, we analyzed the evolutionary pattern of the insulin gene in caviomorph rodents (South American hystricomorph rodents). We found that these rodents have higher rates of nonsynonymous:synonymous substitutions (d(N)/d(S)) than nonhystricomorph rodents and that values are heterogeneous inside the group. We estimated codons under positive selection, specifically the second binding site (A13 and B17) and others related with hexamerization (B18, B20, and B22). In the monomer structure, all selected sites formed a single patch around the second binding site. In the hexamer structure, these amino acids were grouped into three major patches. In this structure, contacts between B chains involved all selected sites (except B18), and between faces in the center of the molecule, all contacts were among selected sites. While there is no clear hypothesis regarding the cause of this drastic change, experimental evidence does show that this group of rodents has some peculiarities in growth function, and, whether coincidental or not, these changes appeared together with important changes in life-history traits.