Fluctuating selection models and Mcdonald-Kreitman type analyses

It is likely that the strength of selection acting upon a mutation varies through time due to changes in the environment. However, most population genetic theory assumes that the strength of selection remains constant. Here we investigate the consequences of fluctuating selection pressures on the quantification of adaptive evolution using McDonald-Kreitman (MK) style approaches. In agreement with previous work, we show that fluctuating selection can generate evidence of adaptive evolution even when the expected strength of selection on a mutation is zero. However, we also find that the mutations, which contribute to both polymorphism and divergence tend, on average, to be positively selected during their lifetime, under fluctuating selection models. This is because mutations that fluctuate, by chance, to positive selected values, tend to reach higher frequencies in the population than those that fluctuate towards negative values. Hence the evidence of positive adaptive evolution detected under a fluctuating selection model by MK type approaches is genuine since fixed mutations tend to be advantageous on average during their lifetime. Never-the-less we show that methods tend to underestimate the rate of adaptive evolution when selection fluctuates

Selective Constraints on Amino Acids Estimated by a Mechanistic Codon Substitution Model with Multiple Nucleotide Changes

Empirical substitution matrices represent the average tendencies of substitutions over various protein families by sacrificing gene-level resolution. We develop a codon-based model, in which mutational tendencies of codon, a genetic code, and the strength of selective constraints against amino acid replacements can be tailored to a given gene. First, selective constraints averaged over proteins are estimated by maximizing the likelihood of each 1-PAM matrix of empirical amino acid (JTT, WAG, and LG) and codon (KHG) substitution matrices. Then, selective constraints specific to given proteins are approximated as a linear function of those estimated from the empirical substitution matrices. Akaike information criterion (AIC) values indicate that a model allowing multiple nucleotide changes fits the empirical substitution matrices significantly better. Also, the ML estimates of transition-transversion bias obtained from these empirical matrices are not so large as previously estimated. The selective constraints are characteristic of proteins rather than species. However, their relative strengths among amino acid pairs can be approximated not to depend very much on protein families but amino acid pairs, because the present model, in which selective constraints are approximated to be a linear function of those estimated from the JTT/WAG/LG/KHG matrices, can provide a good fit to other empirical substitution matrices including cpREV for chloroplast proteins and mtREV for vertebrate mitochondrial proteins. The present codon-based model with the ML estimates of selective constraints and with adjustable mutation rates of nucleotide would be useful as a simple substitution model in ML and Bayesian inferences of molecular phylogenetic trees, and enables us to obtain biologically meaningful information at both nucleotide and amino acid levels from codon and protein sequences.Comment: Table 9 in this article includes corrections for errata in the Table 9 published in 10.1371/journal.pone.0017244. Supporting information is attached at the end of the article, and a computer-readable dataset of the ML estimates of selective constraints is available from 10.1371/journal.pone.001724

The GABBR1 locus and the G1465A variant is not associated with temporal lobe epilepsy preceded by febrile seizures

BACKGROUND: Polymorphism G1465A in the GABBR1 gene has been suggested as a risk factor for non-lesional temporal lobe epilepsy (TLE); however, this genetic association study has not been independently replicated. We attempted to replicate this study in our cohort of patients with TLE. Furthermore, we also analyzed the coding sequence of this gene and searched for disease-causing mutations. METHODS: We included 120 unrelated individuals with TLE that was preceded by febrile seizures (FS) who did not have any evidence of structural lesions suggesting secondary epilepsy. 66 individuals had positive family history of TLE epilepsy and 54 were sporadic. Each patient was genotyped for the presence of G1465A polymorphism. All exons of the GABBR1 gene were screened by single strand confirmation polymorphism method. Genotypes were compared with two independent matched control groups. RESULTS: We detected two A alleles of the G1465A polymorphism in one homozygous control subject (0.87% of all alleles) and one A allele in a patient with TLE (0.45%, not significant). Other detected polymorphisms in coding regions had similar frequencies in epilepsy patients and control groups. No disease causing mutations in the GABBR1 gene were detected in patients with sporadic or familial TLE. CONCLUSION: Our results indicate that TLE preceded by FS is not associated with the polymorphisms or mutations in the GABBR1 gene, including the G1465A polymorphism. The proportion of TLE patients with FS in the original study, reporting this positive association, did not differ between allele A negative and positive cases. Thus, our failure to reproduce this result is likely applicable to all non-lesional TLE epilepsies

Time-resolved thermal lens measurements of the thermo-optical properties of glasses at low temperature down to 20 K

In this work the time resolved thermal lens spectrometry was applied to measure the absolute values of the thermo-optical properties of low silica calcium aluminosilicate and soda lime glasses at low temperatures, in the range between 20 and 300 K. The thermal relaxation calorimetry was used as a complementary technique to determine the specific heat. The results showed a marked decrease of the thermal diffusivity with the temperature rise, with a dependence similar to that of the mean free path (similar to T-1) in the interval between 20 and 70 K, while in the range between 70 and 300 K the dependence was T-(0.33 +/- 0.02). The marked variation of the temperature coefficient of the optical path length change with the temperature rise was attributed to the increase in the coefficient of the electronic polarizability. The results also showed that for the aluminosilicate glass the excess in the specific heat correlated to the so-called boson peak occurred at about 17 K, higher than that of soda lime, which occurs at about 12 K. In conclusion, our results showed the ability of the time resolved thermal lens to determine the thermo-optical properties of glasses at low temperatures, bringing possibilities for experiments in a wide range of optical materials.712

A Strong Deletion Bias in Nonallelic Gene Conversion

Gene conversion is the unidirectional transfer of genetic information between orthologous (allelic) or paralogous (nonallelic) genomic segments. Though a number of studies have examined nucleotide replacements, little is known about length difference mutations produced by gene conversion. Here, we investigate insertions and deletions produced by nonallelic gene conversion in 338 Drosophila and 10,149 primate paralogs. Using a direct phylogenetic approach, we identify 179 insertions and 614 deletions in Drosophila paralogs, and 132 insertions and 455 deletions in primate paralogs. Thus, nonallelic gene conversion is strongly deletion-biased in both lineages, with almost 3.5 times as many conversion-induced deletions as insertions. In primates, the deletion bias is considerably stronger for long indels and, in both lineages, the per-site rate of gene conversion is orders of magnitudes higher than that of ordinary mutation. Due to this high rate, deletion-biased nonallelic gene conversion plays a key role in genome size evolution, leading to the cooperative shrinkage and eventual disappearance of selectively neutral paralogs

Detecting Clusters of Mutations

Positive selection for protein function can lead to multiple mutations within a small stretch of DNA, i.e., to a cluster of mutations. Recently, Wagner proposed a method to detect such mutation clusters. His method, however, did not take into account that residues with high solvent accessibility are inherently more variable than residues with low solvent accessibility. Here, we propose a new algorithm to detect clustered evolution. Our algorithm controls for different substitution probabilities at buried and exposed sites in the tertiary protein structure, and uses random permutations to calculate accurate P values for inferred clusters. We apply the algorithm to genomes of bacteria, fly, and mammals, and find several clusters of mutations in functionally important regions of proteins. Surprisingly, clustered evolution is a relatively rare phenomenon. Only between 2% and 10% of the genes we analyze contain a statistically significant mutation cluster. We also find that not controlling for solvent accessibility leads to an excess of clusters in terminal and solvent-exposed regions of proteins. Our algorithm provides a novel method to identify functionally relevant divergence between groups of species. Moreover, it could also be useful to detect artifacts in automatically assembled genomes

DNA sequence diversity and the efficiency of natural selection in animal mitochondrial DNA

Selection is expected to be more efficient in species that are more diverse because both the efficiency of natural selection and DNA sequence diversity are expected to depend upon the effective population size. We explore this relationship across a data set of 751 mammal species for which we have mitochondrial polymorphism data. We introduce a method by which we can examine the relationship between our measure of the efficiency of natural selection, the nonsynonymous relative to the synonymous nucleotide site diversity (πN/πS), and synonymous nucleotide diversity (πS), avoiding the statistical non-independence between the two quantities. We show that these two variables are strongly negatively and linearly correlated on a log scale. The slope is such that as πS doubles, πN/πS is reduced by 34%. We show that the slope of this relationship differs between the two phylogenetic groups for which we have the most data, rodents and bats, and that it also differs between species with high and low body mass, and between those with high and low mass-specific metabolic rate

A systematic review and meta-analysis of low intensity CBT for psychosis

Sixteen sessions of individual cognitive behavior therapy for people with psychosis (CBTp) is recommended. However, access to CBTp is poor, so the potential of low intensity CBTp (fewer than 16 sessions of face-to-face contact) is being explored. A systematic review and meta-analysis was conducted of 10 controlled trials evaluating low intensity CBTp. Significant between-group effects were found on the primary outcome, symptoms of psychosis, at post-intervention (d = − 0.46, 95% CI: − 0.06, − 0.86) and follow-up (d = − 0.40, 95% CI: − 0.06, − 0.74). Study quality did not moderate post-intervention psychosis outcomes, nor did contact time/number of sessions or therapy format (individual versus group). Between-group effects on secondary outcomes (depression, anxiety and functioning) were not significant at post-intervention, but became significant at follow-up for depression and functioning outcomes (but not for anxiety). Overall, findings suggest that low intensity CBTp shows promise with effect sizes comparable to those found in meta-analyses of CBTp more broadly. We suggest that low intensity CBTp could help widen access. Future research is called for to identify mechanisms of change and to ascertain moderators of outcome so that low intensity CBTp targets key mechanisms (so that scarce therapy time is used effectively) and so that interventions offered are matched to patient need

The Embedding Problem for Markov Models of Nucleotide Substitution

10.1371/journal.pone.0069187PLoS ONE87-POLN

Correlated Mutations: A Hallmark of Phenotypic Amino Acid Substitutions

Point mutations resulting in the substitution of a single amino acid can cause severe functional consequences, but can also be completely harmless. Understanding what determines the phenotypical impact is important both for planning targeted mutation experiments in the laboratory and for analyzing naturally occurring mutations found in patients. Common wisdom suggests using the extent of evolutionary conservation of a residue or a sequence motif as an indicator of its functional importance and thus vulnerability in case of mutation. In this work, we put forward the hypothesis that in addition to conservation, co-evolution of residues in a protein influences the likelihood of a residue to be functionally important and thus associated with disease. While the basic idea of a relation between co-evolution and functional sites has been explored before, we have conducted the first systematic and comprehensive analysis of point mutations causing disease in humans with respect to correlated mutations. We included 14,211 distinct positions with known disease-causing point mutations in 1,153 human proteins in our analysis. Our data show that (1) correlated positions are significantly more likely to be disease-associated than expected by chance, and that (2) this signal cannot be explained by conservation patterns of individual sequence positions. Although correlated residues have primarily been used to predict contact sites, our data are in agreement with previous observations that (3) many such correlations do not relate to physical contacts between amino acid residues. Access to our analysis results are provided at http://webclu.bio.wzw.tum.de/~pagel/supplements/correlated-positions/