Widespread positive selection in synonymous sites of mammalian genes

Evolution of protein sequences is largely governed by purifying selection, with a small fraction of proteins evolving under positive selection. The evolution at synonymous positions in protein-coding genes is not nearly as well understood, with the extent and types of selection remaining, largely, unclear. A statistical test to identify purifying and positive selection at synonymous sites in protein-coding genes was developed. The method compares the rate of evolution at synonymous sites (Ks) to that in intron sequences of the same gene after sampling the aligned intron sequences to mimic the statistical properties of coding sequences. We detected purifying selection at synonymous sites in approximately 28% of the 1,562 analyzed orthologous genes from mouse and rat, and positive selection in approximately 12% of the genes. Thus, the fraction of genes with readily detectable positive selection at synonymous sites is much greater than the fraction of genes with comparable positive selection at nonsynonymous sites, i.e., at the level of the protein sequence. Unlike other genes, the genes with positive selection at synonymous sites showed no correlation between Ks and the rate of evolution in nonsynonymous sites (Ka), indicating that evolution of synonymous sites under positive selection is decoupled from protein evolution. The genes with purifying selection at synonymous sites showed significant anticorrelation between Ks and expression level and breadth, indicating that highly expressed genes evolve slowly. The genes with positive selection at synonymous sites showed the opposite trend, i.e., highly expressed genes had, on average, higher Ks. For the genes with positive selection at synonymous sites, a significantly lower mRNA stability is predicted compared to the genes with negative selection. Thus, mRNA destabilization could be an important factor driving positive selection in nonsynonymous sites, probably, through regulation of expression at the level of mRNA degradation and, possibly, also translation rate. So, unexpectedly, we found that positive selection at synonymous sites of mammalian genes is substantially more common than positive selection at the level of protein sequences. Positive selection at synonymous sites might act through mRNA destabilization affecting mRNA levels and translation.     

Sexual selection drives weak positive selection in protamine genes and high promoter divergence,enhancing sperm competitiveness

Phenotypic adaptations may be the result of changes in gene structure or gene regulation, but little is known about the evolution of gene expression. In addition, it is unclear whether the same selective forces may operate at both levels simultaneously. Reproductive proteins evolve rapidly, but the underlying selective forces promoting such rapid changes are still a matter of debate. In particular, the role of sexual selection in driving positive selection among reproductive proteins remains controversial, whereas its potential influence on changes in promoter regions has not been explored. Protamines are responsible for maintaining DNA in a compacted form in chromosomes in sperm and the available evidence suggests that they evolve rapidly. Because protamines condense DNA within the sperm nucleus, they influence sperm head shape. Here, we examine the influence of sperm competition upon protamine 1 and protamine 2 genes and their promoters, by comparing closely related species of Mus that differ in relative testes size, a reliable indicator of levels of sperm competition. We find evidence of positive selection in the protamine 2 gene in the species with the highest inferred levels of sperm competition. In addition, sperm competition levels across all species are strongly associated with high divergence in protamine 2 promoters that, in turn, are associated with sperm swimming speed. We suggest that changes in protamine 2 promoters are likely to enhance sperm swimming speed by making sperm heads more hydrodynamic. Such phenotypic changes are adaptive because sperm swimming speed may be a major determinant of fertilization success under sperm competition. Thus, when species have diverged recently, few changes in gene-coding sequences are found, while high divergence in promoters seems to be associated with the intensity of sexual selection.     

Codon volatility as an indicator of positive selection: data from eukaryotic genome comparisons

It has been suggested that codon volatility (the proportion of the point-mutation neighbors of a codon that encode different amino acids) can be used as an index of past positive selection. We compared codon volatility with patterns of synonymous and nonsynonymous nucleotide substitution in genome-wide comparisons of orthologous genes between three pairs of related genomes: (1) the protists Plasmodium falciparum and P. yoelii, (2) the fungi Saccharomyces cerevisiae and S. paradoxus, and (3) the mammals mouse and rat. Codon volatility was not consistently associated with an elevated rate of nonsynonymous substitution, as would be expected under positive selection. Rather, the most consistent and powerful correlate of elevated codon volatility was nucleotide content at the second codon position, as expected, given the nature of the genetic code.     

Quantitative measures of sexual selection reveal no evidence for sex-role reversal in a sea spider with prolonged paternal care

Taxa in which males alone invest in postzygotic care of offspring are often considered good models for investigating the proffered relationships between sexual selection and mating systems. In the pycnogonid sea spider Pycnogonum stearnsi, males carry large egg masses on their bodies for several weeks, so this species is a plausible candidate for sex-role reversal (greater intensity of sexual selection on females than on males). Here, we couple a microsatellite-based assessment of the mating system in a natural population with formal quantitative measures of genetic fitness to investigate the direction of sexual selection in P. stearnsi. Both sexes proved to be highly polygamous and showed similar standardized variances in reproductive and mating successes. Moreover, the fertility (number of progeny) of males and females appeared to be equally and highly dependent on mate access, as shown by similar Bateman gradients for the two sexes. The absence of sex-role reversal in this population of P. stearnsi is probably attributable to the fact that males are not limited by brooding space but have evolved an ability to carry large numbers of progeny. Body length was not a good predictor of male mating or reproductive success, so the aim of future studies should be to determine what traits are the targets of sexual selection in this species.     

一个新的睾丸特异表达基因的克隆及功能的初步分析

运用“数据库消减杂交”(digital differential display)方法来筛选人类睾丸特异表达新基因,获得了有差异显示的代表新基因的克隆重叠群。挑选其中一个克隆重叠群HS．326528进行多组织RT—PCR,初步获得该重叠群在人睾丸中有高表达。从该重叠群的IMAGE出发,采用生物信息学的方法快速克隆了一个人类新基因的全长cDNA序列,其全长1044bp,开放阅读框为214～529bp,定位于15q26．2,编码由105个氨基酸组成、分子量为11．7kD、等电点为10．09的一个碱性蛋白,该蛋白与已知蛋白无明显的同源性,克隆实验证明该基因的阅读框完全正确,RT—PCR和Northern blot显示该基因在人类睾丸中特异表达,实时PCR结果表明：该基因在成人睾丸中高表达,在精子中有中度表达,在胚胎睾丸中低表达,推测该基因与精子的生成有关,命名为SRG8(homo sapiens spermatogenesis—related gene 8)(GenBank登录号：AY489187),该基因编码的蛋白定位于细胞核。流式结果分析表明,SRG8基因能够促使HeLa细胞由S期向G2期的转变,从而加速细胞的分裂。这些结果表明SRG8基因可能在睾丸的发育及精子的形成过程中起重要的作用。     

Validation of reference genes for RT-qPCR analysis of CYP4T expression in crucian carp

Reference genes are commonly used for normalization of target gene expression during RT-qPCR analysis. However, no housekeeping genes or reference genes have been identified to be stable across different tissue types or under different experimental conditions. To identify the most suitable reference genes for RT-qPCR analysis of target gene expression in the hepatopancreas of crucian carp (Carassius auratus) under various conditions (sex, age, water temperature, and drug treatments), seven reference genes, including beta actin (ACTB), beta-2 microglobulin (B2M), embryonic elongation factor-1 alpha (EEF1A), glyceraldehyde phosphate dehydrogenase (GAPDH), alpha tubulin (TUBA), ribosomal protein l8 (RPL8) and glucose-6-phosphate dehydrogenase (G6PDH), were evaluated in this study. The stability and ranking of gene expression were analyzed using three different statistical programs: GeNorm, Normfinder and Bestkeeper. The expression errors associated with selection of the genes were assessed by the relative quantity of CYP4T. The results indicated that all the seven genes exhibited variability under the experimental conditions of this research, and the combination of ACTB/TUBA/EEF1A or of ACTB/EEF1A was the best candidate that raised the accuracy of quantitative analysis of gene expression. The findings highlighted the importance of validation of housekeeping genes for research on gene expression under different conditions of experiment and species.     

Recurrent positive selection at bgcn, a key determinant of germ line differentiation, does not appear to be driven by simple coevolution with its partner protein bam

Surveys of nucleotide sequence polymorphism in Drosophila melanogaster and Drosophila simulans were performed at 2 interacting loci crucial for gametogenesis: bag-of-marbles (bam) and benign gonial cell neoplasm (bgcn). At the polymorphism level, both loci appear to be evolving under the expectations of the neutral theory. However, ratios of polymorphism and divergence for synonymous and nonsynonymous mutations depart significantly from neutral expectations for both loci consistent with a previous observation of positive selection at bam. The deviations suggest either an excess of synonymous polymorphisms or an excess of nonsynonymous fixations at both loci. Synonymous evolution appears to conform to neutrality at bam. At bgcn, there is evidence of positive selection affecting preferred synonymous mutations along the D. simulans lineage. However, there is also a significantly higher rate of nonsynonymous fixations at bgcn within D. simulans. Thus, the deviation from neutrality detected by the McDonald-Kreitman test at these 2 loci is likely due to the selective acceleration of nonsynonymous fixations. Differences in the pattern of amino acid fixations between these 2 interacting proteins suggest that the detected positive selection is not due to a simple model of coevolution.     

Statistical analysis of the envelope gene and the prM region of Japanese encephalitis virus: Evidence suggestive of positive selection

The variation in nucleotide sequence observed in the envelope (E) gene and the prM (precursor of M protein) region of different strains of Japanese encephalitis virus (JEV) was analysed. Presence of selective forces acting on these regions was investigated by computing the relative rates of synonymous (K _s) and nonsynonymous (K _a) substitutions. The ratioK _s/K _a was used as an indicator of the overall selective constraints on the amino acid sequence of JEV proteins. The possibility that different regions of the gene may be subject to varying selective pressures was tested by dividing the gene into three regions and estimating theK _s/K _a ratio for each region. On the basis of analysis of a limited number (17) of strains of JEV, evidence suggestive of positive selection acting on certain regions of the E gene of the virus, and in some cases on the entire gene, was obtained. Analysis ofK _a diversity in the prM region of 46 JEV strains grouped into three genotypes revealed that strains included in genotype II were more heterogeneous than strains belonging to genotype I, while the differences between meanK _a values for genotypes I and III and genotypes II and III were not statistically significant. Analysis of host-specific heterogeneity in the prM region revealed that pig isolates were more X_a-diverse than human isolates.     

The effect of positive selection on a sexual reproduction gene in Thalassiosira weissflogii (Bacillariophyta): results obtained from maximum-likelihood and parsimony-based methods

Maximum-Likelihood-based and parsimony-based methods were used to test for potential effects of positive selection on the sexually induced gene 1 (Sig1) in Thalassiosira weissflogii. The Sig proteins are thought to play a role in mediating sperm-egg recognition during the sexual reproduction phase. The results obtained from parsimony-based analyses showed that none of the amino acid sites were influenced by positive selection. Maximum-likelihood analyses indicated that positive selection was affecting a maximum of seven and a minimum of four amino acid sites in the polypeptide derived from Sig1. It was concluded that the results obtained from the maximum-likelihood-based method are more reliable than those obtained from the parsimony-based approach. This is apparently the first study that has shown that reproductive proteins in unicellular eukaryotes are influenced by positive selection.     

Molecular evolution of rickettsia surface antigens: evidence of positive selection

The Rickettsia genus is a group of obligate intracellular parasitic alpha-proteobacteria that includes human pathogens responsible for the typhus disease and various types of spotted fevers. rOmpA and rOmpB are two members of the "surface cell antigen" (Sca) autotransporter (AT) protein family that may play key roles in the adhesion of the Rickettsia cells to the host tissue. These molecules are likely determinants for the pathogenicity of the Rickettsia and represent good candidates for vaccine development. We identified the 17 members of this family of outer-membrane proteins in nine fully sequenced Rickettsia genomes. The typical architecture of the Sca proteins is composed of an N-terminal signal peptide and a C-terminal AT domain that promote the export of the central passenger domain to the outside of the bacteria. A characteristic of this family is the frequent degradation of the genes, which results in different subsets of the sca genes being expressed among Rickettsia species. Here, we present a detailed analysis of their phylogenetic relationships and evolution. We provide strong evidence that rOmpA and rOmpB as well as three other members of the Sca protein family--Sca1, Sca2, and Sca4--have evolved under positive selection. The exclusive distribution of the predicted positively selected sites within the passenger domains of these proteins argues that these regions are involved in the interaction with the host and may be locked in "arms race" coevolutionary conflicts.     

Searching for evidence of selection in avian DNA barcodes

The barcode of life project has assembled a tremendous number of mitochondrial cytochrome c oxidase I (COI) sequences. Although these sequences were gathered to develop a DNA-based system for species identification, it has been suggested that further biological inferences may also be derived from this wealth of data. Recurrent selective sweeps have been invoked as an evolutionary mechanism to explain limited intraspecific COI diversity, particularly in birds, but this hypothesis has not been formally tested. In this study, I collated COI sequences from previous barcoding studies on birds and tested them for evidence of selection. Using this expanded data set, I re-examined the relationships between intraspecific diversity and interspecific divergence and sampling effort, respectively. I employed the McDonald-Kreitman test to test for neutrality in sequence evolution between closely related pairs of species. Because amino acid sequences were generally constrained between closely related pairs, I also included broader intra-order comparisons to quantify patterns of protein variation in avian COI sequences. Lastly, using 22 published whole mitochondrial genomes, I compared the evolutionary rate of COI against the other 12 protein-coding mitochondrial genes to assess intragenomic variability. I found no conclusive evidence of selective sweeps. Most evidence pointed to an overall trend of strong purifying selection and functional constraint. The COI protein did vary across the class Aves, but to a very limited extent. COI was the least variable gene in the mitochondrial genome, suggesting that other genes might be more informative for probing factors constraining mitochondrial variation within species.     

人鼻咽组织特异性表达基因NASG编码产物在细胞中的融合表达

采用PCR技术从人鼻咽上皮组织cDNA文库扩增出人鼻咽组织特异性表达基因NASG基因的编码序列,用Xho I和Sal I酶切NASG基因编码序列的PCR产物及pEGFP-C2载体,产生匹配的粘末端。将回收的pEGFP-C2载体分别与NASG基因编码区酶切片段连接,构建了其编码产物的融合蛋白表达载体pEGFP-C2/NASG,通过脂质体介导分别转染非洲绿猴肾永生化细胞系COS7和鼻咽癌细胞系HNEI,瞬时表达后荧光显微镜观察NASG编码蛋白的活细胞定位,结果表明,NASG编码蛋白在COS7细胞和HNEI细胞的胞浆和胞核中均有表达。     

Inference of positive and negative selection on the 5' regulatory regions of Drosophila genes

Both positive selection and negative selection have been shown to drive the evolution of coding regions. It is of interest to know if the corresponding 5' regions of genes may be subjected to selection of comparable intensities. For such a comparison, we chose the Accessory gland protein (Acp) genes as our test case. About 700 bp and 600 bp for the 5' and coding regions, respectively, of eight previously unstudied genes were sequenced from 21 isogenic lines of D. melanogaster and one line from D. simulans. The ratio of divergence at the amino-acid replacement sites (A) over that at the synonymous sites (S) was twice the ratio for common polymorphism. Interestingly, the 5' region shows the same trend, with the 5'/S divergence ratio being 1.8 times higher than the 5'/S ratio for common polymorphism. There are several possible explanations for the 5'/S ratios, including demography, negative selection, and positive selection. Under normal conditions, positive selection is the most likely explanation. If that is true, about 45 to 50 percent of all fixed differences at both the replacement and 5' sites were adaptive, even though the substitution rate in the former is only half that of the latter (K(A)/K(S) approximately 0.3 vs. K(5')/K(S) approximately 0.6). As previous analyses have indicated, the inclusion of slightly deleterious polymorphism confounds the inference of positive selection. The analysis of published polymorphism data covering 97 verified 5' regions of Drosophila suggests more pronounced selective constraint on the 5' untranslated region and the core promoter (together corresponding to approximately 200 bp in this data set) when compared to the more distal portion of the 5' region of genes.     

Prediction of the Coding Sequences of Unidentified Human Genes. IV. The Coding Sequences of 40 New Genes (KIAA0121-KIAA0160) Deduced by Analysis of cDNA Clones from Human Cell Line KG-1

In this series of projects regarding the accumulation of sequenceinformation of unidentified human genes, we newly deduced thesequences of 40 full-length cDNA clones of human cell line KG-1,and predicted the coding sequences of the corresponding genes,named KIAA0121 to 0160. The results of a computer search ofpublic databases indicated that the sequences of 13 genes wereunrelated to any reported genes, while the remaining 27 genescarried sequences which showed some similarities to known genes.Obvious unique sequences noted were as follows. A stretch oftriplet repeats was contained in each of three genes: Thesewere GAG(Glu) in KIAA0122 and KIAA0147, and TCC(Ser) in KIAA0150.A stretch of 10 amino acidresidues was repeated 21 times inKIAA0139, and a homologous sequence of 76–78 nucleotideswas found repeated 6 times in the untranslated region of KIAA0125.northern hybridization analysis demonstrated that 13 genes wereexpressed in a cell- or tissue-specific manner. Although a vastnumber of expressed sequence tags (ESTs) have been registeredfor comprehensive analysis of cDNA clones, our sequence dataindicated that their distribution is very unbalanced: e.g. whileno EST hit 7 genes, 85 ESTs fell in a single gene.     

Identification of positive selection in disease response genes within members of the Poaceae

Millions of years of coevolution between plants and pathogens can leave footprints on their genomes and genes involved on this interaction are expected to show patterns of positive selection in which novel, beneficial alleles are rapidly fixed within the population. Using information about upregulated genes in maize during Colletotrichum graminicola infection and resources available in the Phytozome database, we looked for evidence of positive selection in the Poaceae lineage, acting on protein coding sequences related with plant defense. We found six genes with evidence of positive selection and another eight with sites showing episodic selection. Some of them have already been described as evolving under positive selection, but others are reported here for the first time including genes encoding isocitrate lyase, dehydrogenases, a multidrug transporter, a protein containing a putative leucine-rich repeat and other proteins with unknown functions. Mapping positively selected residues onto the predicted 3-D structure of proteins showed that most of them are located on the surface, where proteins are in contact with other molecules. We present here a set of Poaceae genes that are likely to be involved in plant defense mechanisms and have evidence of positive selection. These genes are excellent candidates for future functional validation.     

Population genetics of Plasmodium resistance genes in Anopheles gambiae: no evidence for strong selection

Anopheles mosquitoes are the primary vectors for malaria in Africa, transmitting the disease to more than 100 million people annually. Recent functional studies have revealed mosquito genes that are crucial for Plasmodium development, but there is presently little understanding of which genes mediate vector competence in the wild, or evolve in response to parasite-mediated selection. Here, we use population genetic approaches to study the strength and mode of natural selection on a suite of mosquito immune system genes, CTL4, CTLMA2, LRIM1, and APL2 (LRRD7), which have been shown to affect Plasmodium development in functional studies. We sampled these genes from two African populations of An. gambiae s.s., along with several closely related species, and conclude that there is no evidence for either strong directional or balancing selection on these genes. We highlight a number of challenges that need to be met in order to apply population genetic tests for selection in Anopheles mosquitoes; in particular the dearth of suitable outgroup species and the potential difficulties that arise when working within a closely-related species complex.