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.     

高通量测序时代下持家基因定义的发展

持家基因一般指在所有细胞中均稳定表达的一类组成型基因,其产物是维持细胞基本生命活动所必需的。持家基因的定义及其发展,反映了其本质和时代特征。本文从基因表达水平量化技术的角度,结合作者利用转录组数据对猪持家基因的分析,讨论了持家基因的定义及其发展,尤其是技术的进步及数据量的剧增对持家基因的定义所引发的问题。本文总结发现,应该在研究技术的进步中,动态地认识持家基因定义所具有的层次性,从而对持家基因的定义不断赋予新的时代特征。     

基于8种真核生物的整合分析揭示种属特异性小蛋白的功能和进化特征

小蛋白 ( < 100个氨基酸) 广泛存在于三界生命中,具有重要生物功能．早期涉及小蛋白的研究主要集中于少量特殊物种中的蛋白质家族,以及在全基因组尺度预测短小开放读码框(sORFs)的算法开发,但并无跨真核物种的大规模组学分析来揭示小蛋白的功能和进化特征．通过对已知小蛋白和拥有短小开放读码框的基因进行全基因组尺度的计算分析,长度小于100个氨基酸的RefSeq proteins按照其序列保守性被划分为存在于所有8种真核生物、只存在于脊椎动物和只存在于哺乳动物三个进化分类中,此三个进化分类所对应的生物学功能揭示了小蛋白行使种属特异性功能的特征．进一步研究发现,大多数人类特有的小蛋白也是组织表达特异性的,并且绝大多数古老的小蛋白在人体内普遍表达．因此认为,一些真核小蛋白出现并在自然选择压力下富集,行使种属特异性功能,并且以特殊的方式进化和表达．     

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.     

Genomic analysis reveals association of specific SNPs with athletic performance and susceptibility to injuries in professional soccer players

The development of specific and individualized training programs is a possible way to improve athletic performance and minimize injuries in professional athletes. The information regarding the sport's physical demands and the athletes’ physical profile have been, so far, considered as exhaustive for the design of effective training programs. However, it is currently emerging that the genetic profile has to be also taken into consideration. By merging medical and genetic data, it is thus possible to identify the athlete's specific attitude to respond to training, diet, and physical stress. In this context, we performed a study in which 30 professional soccer players, subjected to standard sport medical evaluation and practices, were also screened for genetic polymorphism in five key genes (ACTN3, COL5A1, MCT1, VEGF, and HFE). This genetic analysis represents the central point of a multidisciplinary method that can be adopted by elite soccer teams to obtain an improvement in athletic performance and a concomitant reduction of injuries by tailoring training and nutritional programs. The genetic fingerprinting of single athletes led to the identification of two performance-enhancing polymorphisms (ACTN3 18705C>T, VEGF-634C>G) significantly enriched. Moreover, we derived a genetic model based on the gene set analyzed, which was tentatively used to reduce athletes’ predisposition to injuries, by dictating a personalized nutrition and training program. The potential usefulness of this approach is concordant with data showing that this team has been classified as the healthiest and least injured team in Europe while covering the highest distance/match with the highest number of high-intensity actions/match.     

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

运用“数据库消减杂交”(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.     

The signature of positive selection on standing genetic variation

Considerable interest is focused on the use of polymorphism data to identify regions of the genome that underlie recent adaptations. These searches are guided by a simple model of positive selection, in which a mutation is favored as soon as it arises. This assumption may not be realistic, as environmental changes and range expansions may lead previously neutral or deleterious alleles to become beneficial. We examine what effect this mode of selection has on patterns of variation at linked neutral sites by implementing a new coalescent model of positive directional selection on standing variation. In this model, a neutral allele arises and drifts in the population, then at frequency f becomes beneficial, and eventually reaches fixation. Depending on the value of f, this scenario can lead to a large variance in allele frequency spectra and in levels of linkage disequilibrium at linked, neutral sites. In particular, for intermediate f, the beneficial substitution often leads to a loss of rare alleles--a pattern that differs markedly from the signature of directional selection currently relied on by researchers. These findings highlight the importance of an accurate characterization of the effects of positive selection, if we are to reliably identify recent adaptations from polymorphism data.     

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.     

人管家基因启动子序列中的组合转录调控元件分析

研究表明,第一内含子可能参与基因转录调控.利用统计方法提取人管家基因上游至第一内含子序列中潜在的组合转录调控模体,分析模体间的距离、区域分布等特征,探讨内含子参与基因转录调控的可能性及其参与方式.在管家基因中共获得960对潜在转录调控模体对,其中57%与实验已知的具有转录相互作用的因子对吻合,共涉及12组因子对.分析发现,绝大多数模体对(80%)偏向于上游区域及"上游-内含子"区域,进一步支持了内含子参与基因转录调控的假设,并据此推测内含子与上游序列之间具有转录协同作用,模体在基因转录起始位点(TSS)附近较为集中,模体对的两个模体之间距离较近,60%左右距离在200 bp以内,特别地,65%的模体对特征距离在100 bp以内,短距离间隔有利于转录因子间的协同作用.这些结果将有助于对人基因转录调控机制及内含子功能的深入认识.     

New SNPs for population genetic analysis reveal possible cryptic speciation of eastern Australian sea mullet (Mugil cephalus)

Sustainable management of sea mullet (Mugil cephalus) fisheries needs to account for recent observations of regional‐scale differentiation. Population genetic analysis is sought to assess the situation of this ecologically and economically important fish species in eastern Australian waters. Here, we report (i) new population genetic markers [single nucleotide polymorphisms (SNPs) and potential microsatellites], (ii) first estimates of spatial genetic differentiation and (iii) prospective power tests for designing more comprehensive studies. Six DNA samples from three sampling regions (North Queensland, South Queensland and central New South Wales) on the eastern coast of Australia were used to prepare restriction site associated DNA (RAD) tag libraries from genomic DNA digested with EcoRI and MseI. A pooled sample of regional RAD tag libraries was sequenced using the Roche GS‐FLX Titanium platform. A total of 172 837 raw reads (17.4 Mbp) were retrieved, 95 500 of which were used to discover 1267 SNPs and 1417 microsatellites. A subset of 161 SNPs was validated based on 63 additional DNA samples genotyped using the Sequenom MassArray (iPLEX Gold chemistry). Altogether 92 SNPs (57%) were confirmed, with 40% of these marking fixed variants between northern and southern sampling regions. Our preliminary findings indicate a multispecies fishery stock of M. cephalus in eastern Australian waters, but suggest that strong genetic differentiation occurs north of major fishing grounds. Low potential differentiation within major fishing grounds (e.g. F_ST = 0.0025) can be resolved with a likely power ≥67% by using standard sample sizes of 50 and validated subsets of available markers.     

优秀长跑运动员NADPH氧化酶p22phox亚基基因多态性分子标记的探讨*

目的:解析NADPH 氧化酶p22^phox(NADPH Oxidase p22phox, CYBA)基因的C242T、A930G和A675T多态位点,探讨用于我国北方汉族优秀长跑运动员分子选材标记的可行性。方法:选取优秀长跑运动员123人作为运动员组,127名中国北方汉族普通大学生为对照组,采用基质辅助激光解吸附电离飞行时间质谱检测技术,对CYBA基因的C242T、A930G和A675T位点进行解析并分析比较。结果:与对照组相比,运动员组A930G和A675T位点的基因型和等位基因频率均无显著性差异(P＞0.05),而C242T位点的TC基因型在5 km/10 km运动员组的分布频率存在显著性差异(P＜0.05)。结论:C242T位点的TC基因型可作为5 km/10 km优秀运动员选材的分子标记。