Identification of regions in which positive selection may operate in S-RNase of Rosaceae: Implication for S-allele-specific recognition sites in S-RNase

A stylar S-RNase is associated with gametophytic self-incompatibility in the Rosaceae, Solanaceae, and Scrophulariaceae. This S-RNase is responsible for S-allele-specific recognition in the self-incompatible reaction, but how it functions in specific discrimination is not clear. Window analysis of the numbers of synonymous (d_S) and non-synonymous (d_N) substitutions in rosaceous S-RNases detected four regions with an excess of d_N over d_S in which positive selection may operate (PS regions). The topology of the secondary structure of the S-RNases predicted by the PHD method is very similar to that of fungal RNase Rh whose tertiary structure is known. When the sequences of S-RNases are aligned with the sequence of RNase Rh based on the predicted secondary structures, the four PS regions correspond to two surface sites on the tertiary structure of RNase Rh. These findings suggest that in S-RNases the PS regions also form two sites and are candidates for the recognition sites for S-allele-specific discrimination.     

A systematic search for positive selection in higher plants (Embryophytes)

Background  

Previously, a database characterizing examples of Embryophyte gene family lineages showing evidence of positive selection was reported. Of the gene family trees, 138 Embryophyte branches showed Ka/Ks>>1 and are candidates for functional adaptation. The database and these examples have now been studied in further detail to better understand the molecular basis for plant genome evolution.     

A screen for immunity genes evolving under positive selection in Drosophila

Genes involved in the immune system tend to have higher rates of adaptive evolution than other genes in the genome, probably because they are coevolving with pathogens. We have screened a sample of Drosophila genes to identify those evolving under positive selection. First, we identified rapidly evolving immunity genes by comparing 140 loci in Drosophila erecta and D. yakuba. Secondly, we resequenced 23 of the fastest evolving genes from the independent species pair D. melanogaster and D. simulans, and identified those under positive selection using a McDonald-Kreitman test. There was strong evidence of adaptive evolution in two serine proteases (persephone and spirit) and a homolog of the Anopheles serpin SRPN6, and weaker evidence in another serine protease and the death domain protein dFADD. These results add to mounting evidence that immune signalling pathway molecules often evolve rapidly, possibly because they are sites of host-parasite coevolution.     

Large-scale search of SNPs for type 2 DM susceptibility genes in a Japanese population

The etiology of type 2 diabetes (DM) is polygenic. We investigated here genes and polymorphisms that associate with DM in the Japanese population. Single-nucleotide polymorphisms (SNPs) of 398 derived from 120 candidate genes were examined for association with DM in a population-based case-control study. The study group consisted of 148 cases and 227 controls recruited from Funagata, Japan. No evident subpopulation structure was detected for the tested population. The association tests were conducted with standard allele positivity tables (chi(2) tests) between SNP genotype frequency and case-control status. The independent association of the SNPs from serum triglyceride levels and body mass index was examined by multiple logistic regression analysis. A value of P<0.01 was accepted as statistically significant. Six genes (met proto-oncogene, ATP-binding cassette transporter A1, fatty acid binding protein 2, LDL receptor defect C complementing, aldolase B, and sulfonylurea receptor) were shown to be associated with DM.     

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.     

Homology search for genes

MOTIVATION: Life science researchers often require an exhaustive list of protein coding genes similar to a given query gene. To find such genes, homology search tools, such as BLAST or PatternHunter, return a set of high-scoring pairs (HSPs). These HSPs then need to be correlated with existing sequence annotations, or assembled manually into putative gene structures. This process is error-prone and labor-intensive, especially in genomes without reliable gene annotation. RESULTS: We have developed a homology search solution that automates this process, and instead of HSPs returns complete gene structures. We achieve better sensitivity and specificity by adapting a hidden Markov model for gene finding to reflect features of the query gene. Compared to traditional homology search, our novel approach identifies splice sites much more reliably and can even locate exons that were lost in the query gene. On a testing set of 400 mouse query genes, we report 79% exon sensitivity and 80% exon specificity in the human genome based on orthologous genes annotated in NCBI HomoloGene. In the same set, we also found 50 (12%) gene structures with better protein alignment scores than the ones identified in HomoloGene. AVAILABILITY: The Java implementation is available for download from http://www.bioinformatics.uwaterloo.ca/software.     

Evolutionary models accounting for layers of selection in protein-coding genes and their impact on the inference of positive selection

The selective forces acting on a protein-coding gene are commonly inferred using evolutionary codon models by contrasting the rate of nonsynonymous substitutions to the rate of synonymous substitutions. These models usually assume that the synonymous substitution rate, Ks, is homogenous across all sites, which is justified if synonymous sites are free from selection. However, a growing body of evidence indicates that the DNA and RNA levels of protein-coding genes are subject to varying degrees of selective constraints due to various biological functions encoded at these levels. In this paper, we develop evolutionary models that account for these layers of selection by allowing for both among-site variability of substitution rates at the DNA/RNA level (which leads to Ks variability among protein-coding sites) and among-site variability of substitution rates at the protein level (Ka variability). These models are constructed so that positive selection is either allowed or not. This enables statistical testing of positive selection when variability at the DNA/RNA substitution rate is accounted for. Using this methodology, we show that variability of the baseline DNA/RNA substitution rate is a widespread phenomenon in coding sequence data of mammalian genomes, most likely reflecting varying degrees of selection at the DNA and RNA levels. Additionally, we use simulations to examine the impact that accounting for the variability of the baseline DNA/RNA substitution rate has on the inference of positive selection. Our results show that ignoring this variability results in a high rate of erroneous positive-selection inference. Our newly developed model, which accounts for this variability, does not suffer from this problem and hence provides a likelihood framework for the inference of positive selection on a background of variability in the baseline DNA/RNA substitution rate.     

Identifying candidate positive selection genes in Korean imported pig breeds

Domestication and artificial selection have modified the genome landscape of the pig. The identification of selection signatures in the genome can help to elucidate the selection mechanisms and uncover the causal genes related to the phenotypic variations between domestic pig breeds. We scanned the genomes of Korean imported pig breeds against native breeds using Z-transformed Fst (ZF_ST) and Heterozygosity (ZH_p) statistics to search for the signatures of positive selection. We identified 411 (ZF_ST?=?175; ZH_p?=?256) putatively selected genes in commercial breeds. The gene regions identified were harboring those related to immunity, coat color, reproduction function and other traits. Several genes (e.g., PLSCR4, AGTR1 and CORIN) were related to reproduction traits such as fertility, ovulation rate, and uterine function. This study revealed genes which improve our understanding of the biological mechanisms of higher litter sizes, the phenotype of interest, in higher litter pig breeds.     

Genome-wide scans provide evidence for positive selection of genes implicated in Lassa fever

Rapidly evolving viruses and other pathogens can have an immense impact on human evolution as natural selection acts to increase the prevalence of genetic variants providing resistance to disease. With the emergence of large datasets of human genetic variation, we can search for signatures of natural selection in the human genome driven by such disease-causing microorganisms. Based on this approach, we have previously hypothesized that Lassa virus (LASV) may have been a driver of natural selection in West African populations where Lassa haemorrhagic fever is endemic. In this study, we provide further evidence for this notion. By applying tests for selection to genome-wide data from the International Haplotype Map Consortium and the 1000 Genomes Consortium, we demonstrate evidence for positive selection in LARGE and interleukin 21 (IL21), two genes implicated in LASV infectivity and immunity. We further localized the signals of selection, using the recently developed composite of multiple signals method, to introns and putative regulatory regions of those genes. Our results suggest that natural selection may have targeted variants giving rise to alternative splicing or differential gene expression of LARGE and IL21. Overall, our study supports the hypothesis that selective pressures imposed by LASV may have led to the emergence of particular alleles conferring resistance to Lassa fever, and opens up new avenues of research pursuit.     

Background selection in single genes may explain patterns of codon bias

Background selection involves the reduction in effective population size caused by the removal of recurrent deleterious mutations from a population. Previous work has examined this process for large genomic regions. Here we focus on the level of a single gene or small group of genes and investigate how the effects of background selection caused by nonsynonymous mutations are influenced by the lengths of coding sequences, the number and length of introns, intergenic distances, neighboring genes, mutation rate, and recombination rate. We generate our predictions from estimates of the distribution of the fitness effects of nonsynonymous mutations, obtained from DNA sequence diversity data in Drosophila. Results for genes in regions with typical frequencies of crossing over in Drosophila melanogaster suggest that background selection may influence the effective population sizes of different regions of the same gene, consistent with observed differences in codon usage bias along genes. It may also help to cause the observed effects of gene length and introns on codon usage. Gene conversion plays a crucial role in determining the sizes of these effects. The model overpredicts the effects of background selection with large groups of nonrecombining genes, because it ignores Hill-Robertson interference among the mutations involved.     

Persistence of costly novel genes in the absence of positive selection

Many genetic changes that ultimately lead to adaptive evolution come with a short‐term cost expressed in terms of reduced survival and reproduction. In the absence of genetic drift, it is unclear how such costly mutations may persist. Here we experimentally demonstrate that parasites can promote the persistence of costly genetic variants. We employed a genetically engineered strain (GMMO) of the bacterium Pseudomonas fluorescens as a model of the acquisition of a new gene either through a major mutation or through horizontal transfer, and examined its persistence in different evolving communities comprising an ancestral strain and a lytic bacteriophage. Whereas competition resulted in the elimination of the GMMO, inclusion of the phage promoted GMMO persistence. We provide evidence for why this effect is due to the differential susceptibility of GMMO and ancestral bacteria to phage.     

Two mitochondrial genes under episodic positive selection in subterranean octodontoid rodents

Tuco-tucos (Ctenomys) and related coruros (Spalacopus) are South American subterranean rodents. An energetically demanding lifestyle within the hypoxic, underground atmosphere may change the selective regime on oxidative phosphorylation. We examined whether weak and/or episodic positive directional selection affected the evolution of two mitochondrial genes (COX2, CytB), in a background of purifying selection in these lineages. We estimated rates of synonymous (dS) and non-synonymous (dN) substitutions and found: 1) significantly higher dN/dS ratio in subterranean groups relative to non-subterranean related species, and 2) two codons in each gene under episodic selection: 94 and 277 of COX2 and 269 and 307 of CytB.     

Tabu search model selection for SVM

A model selection method based on tabu search is proposed to build support vector machines (binary decision functions) of reduced complexity and efficient generalization. The aim is to build a fast and efficient support vector machines classifier. A criterion is defined to evaluate the decision function quality which blends recognition rate and the complexity of a binary decision functions together. The selection of the simplification level by vector quantization, of a feature subset and of support vector machines hyperparameters are performed by tabu search method to optimize the defined decision function quality criterion in order to find a good sub-optimal model on tractable times.     

The search for human osteoporosis genes

Osteoporosis is the most prevalent metabolic bone disease and a major clinical and public health problem. Heredity plays an important and well-established role in determining the lifetime risk of this disease. Major efforts are currently underway to identify the specific genes and their allelic variations that contribute to the heritable component to osteoporosis. A number of laboratories are using quantitative trait locus (QTL) methods of genome scanning in families and animal models to identify candidate genomic regions and, ultimately, the genes and genetic variations that lead to osteoporosis. Several chromosomal regions of the human genome have now been linked to osteoporosis-related phenotypes. Although the specific genes contributing to the majority of these linkage signals have not been identified, two positional candidate genes have now been identified: low density lipoprotein receptor-related protein 5 (LRP5) and bone morphogenetic protein 2 (BMP2). A number of QTL has also been identified by cross-breeding strains of mice with variable bone density and several of these QTL have been fine mapped, providing a rich new base for understanding osteoporosis. Genetic association analyses have also provided evidence for a modest relationship between allelic variants in several biological candidate genes and bone mass and the risk of fracture. These ongoing animal and human studies will provide a continuing source of new insight into the genetic regulation of bone and mineral metabolism and the molecular etiology of osteoporosis. The new insight that will emerge from this ongoing research should lead to new ways of diagnosing, preventing and treating the growing clinical and public health problem of osteoporosis.     

Genome re-sequencing analysis uncovers pathogenecity-related genes undergoing positive selection in Magnaporthe oryzae

Rice blast caused by Magnaporthe oryzae(M.oryzae) is one of the most destructive diseases,which causes significant rice yield losses and affects global food security.To better understand genetic variations among different isolates of M.oryzae in nature,we re-sequenced the genomes of two field isolates,CH43 and Zhong-10-8-14,which showed distinct pathogenecity on most of the rice cultivars.Genome-wide genetic variation analysis reveals that ZHONG-10-8-14 exhibits higher sequence variations than CH43.Structural variations(SVs) detection shows that the sequence variations primarily occur in exons and intergenic regions.Bioinformatics analysis for gene variations reveals that many pathogenecity-related pathways are enriched.In addition,193 candidate effectors with various DNA polymorphisms were identified,including two known effectors AVR-Pik and AVR-Pital.Comparative polymorphism analysis of thirteen randomly selected effectors suggests that the genetic variations of effectors are under positive selection.The expression pattern analysis of several pathogenecity-related variant genes indicates that these genes are differentially regulated in two isolates,with much higher expression levels in Zhong-10-8-14 than CH43.Our data demonstrate that the genetic variations of effectors and pathogenecity-related genes are under positive selection,resulting in the distinct pathogenicities of CH43 and Zhong-10-8-14 on rice.     

Rapid detection of positive selection in genes and genomes through variation clusters

Positive selection in genes and genomes can point to the evolutionary basis for differences among species and among races within a species. The detection of positive selection can also help identify functionally important protein regions and thus guide protein engineering. Many existing tests for positive selection are excessively conservative, vulnerable to artifacts caused by demographic population history, or computationally very intensive. I here propose a simple and rapid test that is complementary to existing tests and that overcomes some of these problems. It relies on the null hypothesis that neutrally evolving DNA regions should show a Poisson distribution of nucleotide substitutions. The test detects significant deviations from this expectation in the form of variation clusters, highly localized groups of amino acid changes in a coding region. In applying this test to several thousand human-chimpanzee gene orthologs, I show that such variation clusters are not generally caused by relaxed selection. They occur in well-defined domains of a protein's tertiary structure and show a large excess of amino acid replacement over silent substitutions. I also identify multiple new human-chimpanzee orthologs subject to positive selection, among them genes that are involved in reproductive functions, immune defense, and the nervous system.     

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.     

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.