A hierarchical Bayesian model for predicting the functional consequences of amino-acid polymorphisms

Genetic polymorphisms in deoxyribonucleic acid coding regions may have a phenotypic effect on the carrier, e.g. by influencing susceptibility to disease. Detection of deleterious mutations via association studies is hampered by the large number of candidate sites; therefore methods are needed to narrow down the search to the most promising sites. For this, a possible approach is to use structural and sequence-based information of the encoded protein to predict whether a mutation at a particular site is likely to disrupt the functionality of the protein itself. We propose a hierarchical Bayesian multivariate adaptive regression spline (BMARS) model for supervised learning in this context and assess its predictive performance by using data from mutagenesis experiments on lac repressor and lysozyme proteins. In these experiments, about 12 amino-acid substitutions were performed at each native amino-acid position and the effect on protein functionality was assessed. The training data thus consist of repeated observations at each position, which the hierarchical framework is needed to account for. The model is trained on the lac repressor data and tested on the lysozyme mutations and vice versa. In particular, we show that the hierarchical BMARS model, by allowing for the clustered nature of the data, yields lower out-of-sample misclassification rates compared with both a BMARS and a frequen-tist MARS model, a support vector machine classifier and an optimally pruned classification tree.

Mutational meltdown in primary endosymbionts: selection limits Muller's ratchet

Background Primary bacterial endosymbionts of insects (p-endosymbionts) are thought to be undergoing the process of Muller's ratchet where they accrue slightly deleterious mutations due to genetic drift in small populations with negligible recombination rates. If this process were to go unchecked over time, theory predicts mutational meltdown and eventual extinction. Although genome degradation is common among p-endosymbionts, we do not observe widespread p-endosymbiont extinction, suggesting that Muller's ratchet may be slowed or even stopped over time. For example, selection may act to slow the effects of Muller's ratchet by removing slightly deleterious mutations before they go to fixation thereby causing a decrease in nucleotide substitutions rates in older p-endosymbiont lineages. Methodology/Principal Findings To determine whether selection is slowing the effects of Muller's ratchet, we determined the age of the Candidatus Riesia/sucking louse assemblage and analyzed the nucleotide substitution rates of several p-endosymbiont lineages that differ in the length of time that they have been associated with their insect hosts. We find that Riesia is the youngest p-endosymbiont known to date, and has been associated with its louse hosts for only 13–25 My. Further, it is the fastest evolving p-endosymbiont with substitution rates of 19–34% per 50 My. When comparing Riesia to other insect p-endosymbionts, we find that nucleotide substitution rates decrease dramatically as the age of endosymbiosis increases. Conclusions/Significance A decrease in nucleotide substitution rates over time suggests that selection may be limiting the effects of Muller's ratchet by removing individuals with the highest mutational loads and decreasing the rate at which new mutations become fixed. This countering effect of selection could slow the overall rate of endosymbiont extinction.

Non-hemagglutinating flaviviruses: molecular mechanisms for the emergence of new strains via adaptation to European ticks

Tick-borne encephalitis virus (TBEV) causes human epidemics across Eurasia. Clinical manifestations range from inapparent infections and fevers to fatal encephalitis but the factors that determine disease severity are currently undefined. TBEV is characteristically a hemagglutinating (HA) virus; the ability to agglutinate erythrocytes tentatively reflects virion receptor/fusion activity. However, for the past few years many atypical HA-deficient strains have been isolated from patients and also from the natural European host tick, Ixodes persulcatus. By analysing the sequences of HA-deficient strains we have identified 3 unique amino acid substitutions (D67G, E122G or D277A) in the envelope protein, each of which increases the net charge and hydrophobicity of the virion surface. Therefore, we genetically engineered virus mutants each containing one of these 3 substitutions; they all exhibited HA-deficiency. Unexpectedly, each genetically modified non-HA virus demonstrated increased TBEV reproduction in feeding Ixodes ricinus, not the recognised tick host for these strains. Moreover, virus transmission efficiency between infected and uninfected ticks co-feeding on mice was also intensified by each substitution. Retrospectively, the mutation D67G was identified in viruses isolated from patients with encephalitis. We propose that the emergence of atypical Siberian HA-deficient TBEV strains in Europe is linked to their molecular adaptation to local ticks. This process appears to be driven by the selection of single mutations that change the virion surface thus enhancing receptor/fusion function essential for TBEV entry into the unfamiliar tick species. As the consequence of this adaptive mutagenesis, some of these mutations also appear to enhance the ability of TBEV to cross the human blood-brain barrier, a likely explanation for fatal encephalitis. Future research will reveal if these emerging Siberian TBEV strains continue to disperse westwards across Europe by adaptation to the indigenous tick species and if they are associated with severe forms of TBE.

Mutations in the NS2B and NS3 genes affect mouse neuroinvasiveness of a Western European field strain of tick-borne encephalitis virus

An attenuated strain (263) of the tick-borne encephalitis virus, isolated from field ticks, was either serially subcultured, 5 times in mice, or at 40 degrees C in PS cells, producing 2 independent strains, 263-m5 and 263-TR with identical genomes; both strains exhibited increased plaque size, neuroinvasiveness and temperature-resistance. Sequencing revealed two unique amino acid substitutions, one mapping close to the catalytic site of the viral protease. These observations imply that virus adaptation from ticks to mammals occurs by selection of pre-existing virulent variants from the quasispecies population rather than by the emergence of new random mutations. The significance of these observations is discussed. (c) 2008 Elsevier Inc. All rights reserved.

Detection of unusual mutation within the VP1 region of different re-isolates of poliovirus Sabin vaccine

In the present study, a genomic analysis of full VP1 sequence region of 15 clinical re-isolates (14 healthy vaccinees and one bone marrow tumor patient) was conducted, aiming to the identification of mutations and to the assessment of their impact on virus fitness, providing also insights relevant with the natural evolution of Sabin strains. Clinical re-isolates were analyzed by RT-PCR, sequencing and computational analysis. Some re-isolates were characterized by an unusual mutational pattern in which non-synonymous mutations outnumbered the synonymous ones. Furthermore, the majority of amino-acid substitutions were located in the capsid exterior, specifically in N-Ags, near N-Ags and in the north rim of the canyon. Also mutations, which are well-known determinants of attenuation, were identified. The results of this study propose that some re-isolates are characterized by an evolutionary pattern in which non-synonymous mutations with a direct phenotypic impact on viral fitness are fixed in viral genomes, in spite of synonymous ones with no phenotypic impact on viral fitness. Results of the present retrospective characterization of Sabin clinical re-isolates, based on the full VP1 sequence, suggest that vaccine-derived viruses may make their way through narrow breaches and may evolve into transmissible pathogens even in adequately immunized populations. For this reason increased poliovirus laboratory surveillance should be permanent and full VP1 sequence analysis should be conducted even in isolates originating from healthy vaccinees.

Polymorphisms in the apolipoprotein L1 gene and their effects on blood lipid and glucose levels in middle age males

Apolipoprotein L1 in plasma is associated with high- density lipoprotein. Novel APOL1 polymorphisms are investigated along with the association of two common haplotypes (Lys166Glu, Ile244Met, Lys271Arg) with circulating lipid and glucose levels. Although the amino acid substitutions occur in the amphipathic alpha helices region involved in lipid binding, these substitutions were found not to independently account for variability in circulating lipid and glucose levels in 149 middle age males.

Introduction of vinyl and hydroxymethyl functionalities at C-4 of glucose-derived substrates: synthesis of spirocyclic, bicyclic, and tricyclic nucleosides

Installing hydroxymethyl and hydroxyethyl substitutions at C-4 through vinylation and hydroboration-oxidation reactions of the C-4 bis-hydroxymethyl derivative of D-glucose based substrate, and inserting heteroatoms thereafter permitted formation of N-, O-, or S-heterocycles leading to [4,5]or [5,5]-spirocycles and a bicyclo[3.3.0]octane product. Some of the spirocycles were converted to spironucleosides under Vorbruggen glycosidation reaction conditions. Similarly, the bicyclic product was elaborated to the corresponding bicyclic nucleoside as well as an unexpected tricyclic nucleoside.

Long-term monounsaturated fatty acid diets reduce platelet aggregation in healthy young subjects

The aim of the present study was to compare the response of a range of atherogenic and thrombogenic risk markers to two dietary levels of saturated fatty acid (SFA) substitution with monounsaturated fatty acids (MUFA) in students living in a university hall of residence. Although the benefits of such diets have been reported for plasma lipoproteins in high-risk groups, more needs to be known about effects of more modest SFA-MUFA substitutions over the long term and in young healthy adults. In a parallel design over 16 weeks, fifty-one healthy young subjects were randomised to one of two diets: (1) a moderate-MUFA diet in which 16 g dietary SFA/100 g total fatty acids were substituted with MUFA (n 25); (2) a high-MUFA diet in which 33 g dietary SFA/100 g total fatty acids were substituted with MUFA (n 26). All subjects followed an 8-week run-in diet (reference diet), with a fatty acid composition close to the UK average values. There were no differences in plasma lipid responses between the two diets over 16 weeks of the study with similar reductions in total cholesterol (P<0.001) and LDL-cholesterol (P<0.01) in both groups; a small but significant reduction in HDL-cholesterol was also observed in both groups (P<0.01). Platelet responses to ADP (P<0.01) and arachidonic acid (P<0.05) differed with time on the two diets; at 16 weeks, platelet aggregatory response to ADP was significantly lower on the high-MUFA than the moderate-MUFA (P<0.01) diet; ADP responses were also significantly lower within this group at 8 (P< 0.05) and 16 (P< 0.01) weeks compared with baseline. There were no differences in fasting factor VII activity (factors VIII and VIIag), fibrinogen concentration or tissue-type plasminogen activator activity between the diets. There were no differences in postprandial factor VIII responses to a standard meal (area under the curve) between the diets after 16 weeks, but postprandial factor VIII response was lower than on the high-MUFA diet compared with baseline (P<0.01). In conclusion, a high-MUFA diet sustains potentially beneficial effects on platelet aggregation and postprandial activation of factor VII. Moderate or high substitution of MUFA for SFA achieves similar reductions in fasting blood lipids in young healthy subjects.

Polymorphisms in the apolipoprotein L1 gene and their effects on blood lipid and glucose levels in middle age males

Apolipoprotein L1 in plasma is associated with high- density lipoprotein. Novel APOL1 polymorphisms are investigated along with the association of two common haplotypes (Lys166Glu, Ile244Met, Lys271Arg) with circulating lipid and glucose levels. Although the amino acid substitutions occur in the amphipathic alpha helices region involved in lipid binding, these substitutions were found not to independently account for variability in circulating lipid and glucose levels in 149 middle age males.

Achievement of dietary fatty acid intakes in long-term controlled intervention studies: approach and methodology

Objective: To describe the calculations and approaches used to design experimental diets of differing saturated fatty acid (SFA) and monounsaturated fatty acid (MUFA) compositions for use in a long-term dietary intervention study, and to evaluate the degree to which the dietary targets were met. Design, setting and subjects: Fifty-one students living in a university hall of residence consumed a reference (SFA) diet for 8 weeks followed by either a moderate MUFA (MM) diet or a high MUFA (HM) diet for 16 weeks. The three diets were designed to differ only in their proportions of SFA and MUFA, while keeping total fat, polyunsaturated fatty acids (PUFA), trans-fatty acids, and the ratio of palmitic to stearic acid, and n-6 to n-3 PUFA, unchanged. Results: Using habitual diet records and a standardised database for food fatty acid compositions, a sequential process of theoretical fat substitutions enabled suitable fat sources for use in the three diets to be identified, and experimental margarines for baking, spreading and the manufacture of snack foods to be designed. The dietary intervention was largely successful in achieving the fatty acid targets of the three diets, although unintended differences between the original target and the analysed fatty acid composition of the experimental margarines resulted in a lower than anticipated MUFA intake on the HM diet, and a lower ratio of palmitic to stearic acid compared with the reference or MM diet. Conclusions: This study has revealed important theoretical considerations that should be taken into account when designing diets of specific fatty acid composition, as well as practical issues of implementation.

Non-word repetition in adolescents with Specific Language Impairment and Autism plus Language Impairments: a qualitative analysis

Non-word repetition (NWR) was investigated in adolescents with typical development, Specific Language Impairment (SLI) and Autism Plus language Impairment (ALI) (n = 17, 13, 16, and mean age 14;4, 15;4, 14;8 respectively). The study evaluated the hypothesis that poor NWR performance in both groups indicates an overlapping language phenotype (Kjelgaard & Tager-Flusberg, 2001). Performance was investigated both quantitatively, e.g. overall error rates, and qualitatively, e.g. effect of length on repetition, proportion of errors affecting phonological structure, and proportion of consonant substitutions involving manner changes. Findings were consistent with previous research (Whitehouse, Barry, & Bishop, 2008) demonstrating a greater effect of length in the SLI group than the ALI group, which may be due to greater short-term memory limitations. In addition, an automated count of phoneme errors identified poorer performance in the SLI group than the ALI group. These findings indicate differences in the language profiles of individuals with SLI and ALI, but do not rule out a partial overlap. Errors affecting phonological structure were relatively frequent, accounting for around 40% of phonemic errors, but less frequent than straight Consonant-for-Consonant or vowel-for-vowel substitutions. It is proposed that these two different types of errors may reflect separate contributory mechanisms. Around 50% of consonant substitutions in the clinical groups involved manner changes, suggesting poor auditory-perceptual encoding. From a clinical perspective algorithms which automatically count phoneme errors may enhance sensitivity of NWR as a diagnostic marker of language impairment. Learning outcomes: Readers will be able to (1) describe and evaluate the hypothesis that there is a phenotypic overlap between SLI and Autism Spectrum Disorders (2) describe differences in the NWR performance of adolescents with SLI and ALI, and discuss whether these differences support or refute the phenotypic overlap hypothesis, and (3) understand how computational algorithms such as the Levenshtein Distance may be used to analyse NWR data.

Detection of somaclonal variation during cocoa somatic embryogenesis characterised using cleaved amplified polymorphic sequence and the new freeware Artbio

The scarcity and stochastic nature of genetic mutations presents a significant challenge for scientists seeking to characterise de novo mutation frequency at specific loci. Such mutations can be particularly numerous during regeneration of plants from in vitro culture and can undermine the value of germplasm conservation efforts. We used cleaved amplified polymorphic sequence (CAPS) analysis to characterise new mutations amongst a clonal population of cocoa plants regenerated via a somatic embryogenesis protocol used previously for cocoa cryopreservation. Efficacy of the CAPS system for mutation detection was greatly improved after an ‘a priori’ in silico screen of reference target sequences for actual and potential restriction enzyme recognition sites using a new freely available software called Artbio. Artbio surveys known sequences for existing restriction enzyme recognition sites but also identifies all single nucleotide polymorphism (SNP) deviations from such motifs. Using this software, we performed an in silico screen of seven loci for restriction sites and their potential mutant SNP variants that were possible from 21 restriction enzymes. The four most informative locus-enzyme combinations were then used to survey the regenerant populations for de novo mutants. We characterised the pattern of point mutations and, using the outputs of Artbio, calculated the ratio of base substitution in 114 somatic embryo-derived cocoa regenerants originating from two explant genotypes. We found 49 polymorphisms, comprising 26.3% of the samples screened, with an inferred rate of 2.8 × 10−3 substitutions/screened base. This elevated rate is of a similar order of magnitude to previous reports of de novo microsatellite length mutations arising in the crop and suggests caution should be exercised when applying somatic embryogenesis for the conservation of plant germplasm.

Replication enhancer elements within the open reading frame of tick-borne encephalitis virus and their evolution within the Flavivirus genus

We provide experimental evidence of a replication enhancer element (REE) within the capsid gene of tick-borne encephalitis virus (TBEV, genus Flavivirus). Thermodynamic and phylogenetic analyses predicted that the REE folds as a long stable stem–loop (designated SL6), conserved among all tick-borne flaviviruses (TBFV). Homologous sequences and potential base pairing were found in the corresponding regions of mosquito-borne flaviviruses, but not in more genetically distant flaviviruses. To investigate the role of SL6, nucleotide substitutions were introduced which changed a conserved hexanucleotide motif, the conformation of the terminal loop and the base-paired dsRNA stacking. Substitutions were made within a TBEV reverse genetic system and recovered mutants were compared for plaque morphology, single-step replication kinetics and cytopathic effect. The greatest phenotypic changes were observed in mutants with a destabilized stem. Point mutations in the conserved hexanucleotide motif of the terminal loop caused moderate virus attenuation. However, all mutants eventually reached the titre of wild-type virus late post-infection. Thus, although not essential for growth in tissue culture, the SL6 REE acts to up-regulate virus replication. We hypothesize that this modulatory role may be important for TBEV survival in nature, where the virus circulates by non-viraemic transmission between infected and non-infected ticks, during co-feeding on local rodents.

Patterns of nucleotide substitution in angiosperm cpDNA trnL (UAA)-trnF (GAA) regions

Patterns of substitution in chloroplast encoded trnL_F regions were compared between species of Actaea (Ranunculales), Digitalis (Scrophulariales), Drosera (Caryophyllales), Panicoideae (Poales), the small chromosome species clade of Pelargonium (Geraniales), each representing a different order of flowering plants, and Huperzia (Lycopodiales). In total, the study included 265 taxa, each with > 900-bp sequences, totaling 0.24 Mb. Both pairwise and phylogeny-based comparisons were used to assess nucleotide substitution patterns. In all six groups, we found that transition/transversion ratios, as estimated by maximum likelihood on most-parsimonious trees, ranged between 0.8 and 1.0 for ingroups. These values occurred both at low sequence divergences, where substitutional saturation, i.e., multiple substitutions having occurred at the same (homologous) nucleotide position, was not expected, and at higher levels of divergence. This suggests that the angiosperm trnL-F regions evolve in a pattern different from that generally observed for nuclear and animal mtDNA (transitional/transversion ratio > or = 2). Transition/transversion ratios in the intron and the spacer region differed in all alignments compared, yet base compositions between the regions were highly similar in all six groups. A>-C transversions were significantly less frequent than the other four substitution types. This correlates with results from studies on fidelity mechanisms in DNA replication that predict A<->T and G<->C transversions to be least likely to occur. It therefore strengthens confidence in the link between mutation bias at the polymerase level and the actual fixation of substitutions as recorded on evolutionary trees, and concomitantly, in the neutrality of nucleotide substitutions as phylogenetic markers.

Evolutionary analysis of novel serine proteases in the venom gland transcriptome of Bitis gabonica rhinoceros

Background: Serine proteases are major components of viper venom and target various stages of the blood coagulation system in victims and prey. A better understanding of the diversity of serine proteases and other enzymes present in snake venom will help to understand how the complexity of snake venom has evolved and will aid the development of novel therapeutics for treating snake bites. Methodology and Principal Findings: Four serine protease-encoding genes from the venom gland transcriptome of Bitis gabonica rhinoceros were amplified and sequenced. Mass spectrometry suggests the four enzymes corresponding to these genes are present in the venom of B. g. rhinoceros. Two of the enzymes, rhinocerases 2 and 3 have substitutions to two of the serine protease catalytic triad residues and are thus unlikely to be catalytically active, though they may have evolved other toxic functions. The other two enzymes, rhinocerases 4 and 5, have classical serine protease catalytic triad residues and thus are likely to be catalytically active, however they have glycine rather than the more typical aspartic acid at the base of the primary specificity pocket (position 189). Based on a detailed analysis of these sequences we suggest that alternative splicing together with individual amino acid mutations may have been involved in their evolution. Changes within amino acid segments which were previously proposed to undergo accelerated change in venom serine proteases have also been observed. Conclusions and Significance: Our study provides further insight into the diversity of serine protease isoforms present within snake venom and discusses their possible functions and how they may have evolved. These multiple serine protease isoforms with different substrate specificities may enhance the envenomation effects and help the snake to adapt to new habitats and diets. Our findings have potential for helping the future development of improved therapeutics for snake bites.