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.

Bilateral high-frequency stimulation of the subthalamic nucleus on attentional performance: transient deleterious effects and enhanced motivation in both intact and parkinsonian rats.

It is now well established that subthalamic nucleus high-frequency stimulation (STN HFS) alleviates motor problems in Parkinson's disease. However, its efficacy for cognitive function remains a matter of debate. The aim of this study was to assess the effects of STN HFS in rats performing a visual attentional task. Bilateral STN HFS was applied in intact and in bilaterally dopamine (DA)-depleted rats. In all animals, STN HFS had a transient debilitating effect on all the variables measured in the task. In DA-depleted rats, STN HFS did not alleviate the deficits induced by the DA lesion such as omissions and latency to make correct responses, but induced perseverative approaches to the food magazine, an indicator of enhanced motivation. In sham-operated controls, STN HFS significantly reduced accuracy and induced perseverative behaviour, mimicking partially the effects of bilateral STN lesions in the same task. These results are in line with the hypothesis that STN HFS only partially mimics inactivation of STN produced by lesioning and confirm the motivational exacerbation induced by STN inactivation.

Puroindoline-b mutations control the lipid binding interactions in mixed puroindoline-a: puroindoline-b systems

The interactions have been investigated of puroindoline-a (Pin-a) and mixed protein systems of Pin-a and wild-type puroindoline-b (Pin-b+) or puroindoline-b mutants (G46S mutation (Pin bH) or W44R mutation (Pin-bS)) with condensed phase monolayers of an anionic phospholipid (L-α-dipalmitoylphosphatidyl-dl-glycerol (DPPG)) at the air/water interface. The interactions of the mixed systems were studied at three different concentration ratios of Pin-a:Pin-b, namely 3:1, 1:1 and 1:3 in order to establish any synergism in relation to lipid binding properties. Surface pressure measurements revealed that Pin-a interaction with DPPG monolayers led to an equilibrium surface pressure increase of 8.7 ± 0.6 mN m-1. This was less than was measured for Pin-a:Pin-b+ (9.6 to 13.4 mN m-1), but was significantly more than was measured for Pin-a:Pin-bH (4.0 to 6.2 mN m-1) or Pin-a:Pin-bS (3.8 to 6.3 mN m-1) over the complete range of concentration ratio. Consequently, surface pressure increases were shown to correlate to endosperm hardness phenotype, with puroindolines present in hard-textured wheat varieties yielding lower equilibrium surface pressure changes. Integrated amide I peak areas from corresponding external reflectance Fourier-transform infrared (ER-FTIR) spectra, used to indicate levels of protein adsorption to the lipid monolayers, showed that differences in adsorbed amount were less significant. The data therefore suggest that Pin-b mutants having single residue substitutions within their tryptophan-rich loop that are expressed in some hard-textured wheat varieties influence the degree of penetration of Pin-a and Pin-b into anionic phospholipid films. These findings highlight the key role of the tryptophan-rich loop in puroindoline-lipid interactions.

Mutations in H5N1 influenza virus hemagglutinin that confer binding to human tracheal airway epithelium

The emergence in 2009 of a swine-origin H1N1 influenza virus as the first pandemic of the 21st Century is a timely reminder of the international public health impact of influenza viruses, even those associated with mild disease. The widespread distribution of highly pathogenic H5N1 influenza virus in the avian population has spawned concern that it may give rise to a human influenza pandemic. The mortality rate associated with occasional human infection by H5N1 virus approximates 60%, suggesting that an H5N1 pandemic would be devastating to global health and economy. To date, the H5N1 virus has not acquired the propensity to transmit efficiently between humans. The reasons behind this are unclear, especially given the high mutation rate associated with influenza virus replication. Here we used a panel of recombinant H5 hemagglutinin (HA) variants to demonstrate the potential for H5 HA to bind human airway epithelium, the predominant target tissue for influenza virus infection and spread. While parental H5 HA exhibited limited binding to human tracheal epithelium, introduction of selected mutations converted the binding profile to that of a current human influenza strain HA. Strikingly, these amino-acid changes required multiple simultaneous mutations in the genomes of naturally occurring H5 isolates. Moreover, H5 HAs bearing intermediate sequences failed to bind airway tissues and likely represent mutations that are an evolutionary "dead end." We conclude that, although genetic changes that adapt H5 to human airways can be demonstrated, they may not readily arise during natural virus replication. This genetic barrier limits the likelihood that current H5 viruses will originate a human pandemic.

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.

Characterization of recombinant influenza B viruses with key neuraminidase inhibitor resistance mutations

Objectives and methods: An influenza B virus plasmid-based rescue system was used to introduce site-specific mutations, previously observed in neuraminidase (NA) inhibitor-resistant viruses, into the NA protein of six recombinant viruses. Three mutations observed only among in vitro selected zanamivir-resistant influenza A mutants were introduced into the B/Beijing/1/87 virus NA protein, to change residue E116 to glycine, alanine or aspartic acid. Residue E116 was also mutated to valine, a mutation found in the clinic among oseltamivir-resistant viruses. An arginine to lysine change at position 291 (292 N2 numbering) mimicked that seen frequently in influenza A N2 clinical isolates resistant to oseltamivir. Similarly, an arginine to lysine change at position 149 (152 in N2 numbering) was made to reproduce the change found in the only reported zanamivir-resistant clinical isolate of influenza B virus. In vitro selection and prolonged treatment in the clinic leads to resistance pathways that require compensatory mutations in the haemagglutinin gene, but these appear not to be important for mutants isolated from immunocompetent patients. The reverse genetics system was therefore used to generate mutants containing only the NA mutation. Results and conclusions: With the exception of a virus containing the E116G mutation, mutant viruses were attenuated to different levels in comparison with wild-type virus. This attenuation was a result of altered NA activity or stability depending on the introduced mutation. Mutant viruses displayed increased resistance to zanamivir, oseltamivir and peramivir, with certain viruses displaying cross-resistance to all three drugs.

Identification of mutations conferring insecticide-insensitive AChE in the cotton-melon aphid, Aphis gossypii Glover

We have identified two mutations in the ace1 gene of Aphis gossypii that are associated with insensitivity of acetylcholinesterase (AChE) to carbamate and organophosphate insecticides. The first of these, S431F (equivalent to F331 in Torpedo californica), is associated with insensitivity to the carbamate insecticide pirimicarb in a range of A. gossypii clones. The S431F mutation is also found in the peach-potato aphid, Myzus persicae (Sulzer), and a rapid RFLP diagnostic allows the identification of individuals of both aphid species with a resistant genotype. This diagnostic further revealed the presence of S431 in several other pirimicarb-susceptible aphid species. The serine at this position in the wild-type enzyme has only been reported for aphids and provides a molecular explanation of why pirimicarb has a specific aphicidal action. A less specific insensitivity to a wide range of carbamates and organophosphates is associated with a second mutation, A302S (A201 in T. californica).

Influence of mutations affecting gonadotropin production or responsiveness on expression of inhibin subunit mRNA and protein in the mouse ovary

Measurement of inhibins A and B in the serum of normal cyclic rodents has implicated FSH in the regulation of these peptides within the ovary. To extend these observations we have used a panel of mutant mice carrying mutations which affect either the production of, or the ability to respond to, FSH and LH. As a consequence, the females are infertile and show different degrees of follicular development. The aim of this study was to measure inhibin gene transcription in the ovaries of these mutant females together with inhibin protein levels in ovaries and serum and to relate these to follicular development within the ovary. Comparison was made with a pool of normal/heterozygous females. In hpg females where lack of GnRH production results in the absence of gonadotropin synthesis, in FSHbeta knockout (FSHbetaKO) females where disruption of the gene encoding FSHbeta results in the absence of FSH production, and in FSH receptor knockout (FSHRKO) females which are unable to respond to circulating FSH, follicular development remains at the pre-antral stage in these three mutants. Only in the hpg females were common inhibin alpha subunit mRNA levels significantly lower than normal. In these three mutants, however, mRNA levels for both the betaA and betaB subunits were extremely low compared with normal mice. At the protein level, neither inhibin A nor B was detected in the serum of these three mutants; however inhibin B, albeit at very low levels, was detectable within the ovaries. These observations confirm a major role for FSH in the control of transcription of the RA and betaB genes but suggest that the constitutive transcription of the alpha subunit is less dependent on FSH. In contrast, in LH receptor knockout (LuRKO) female mice inhibin betaA subunit mRNA levels were similar to those measured in normal/heterozygous females but levels of inhibin alpha and betaB subunit mRNAs were significantly higher than in the normal group. This was reflected in significantly higher inhibin B protein levels in ovaries and serum. An inability to respond to LH combined with high circulating levels of FSH leads to a high proportion of antral follicles in LuRKO females, with granulosa cells constituting the major cell type within the ovary. The high percentage of antral granulosa cells is likely to account for the significantly higher levels of inhibin B production in these ovaries.

Attenuating mutations in the influenza virus genome which may increase the safety of vaccine production

Influenza virus epidemics occur on an annual basis and cause severe disease in the very young and old. The vaccine administered to high-risk groups is generated by amplifying reassortant viruses, with chronologically relevant viral surface antigens, in eggs. Every 20 years or so, influenza pandemics occur causing widespread fatality in all age groups. These viruses display novel viral surface antigens acquired from a zoonotic source, and vaccination against them poses new issues since production of large amounts of a respiratory virus containing novel surface antigens could be dangerous for those involved in manufacture. To minimise risks, it is advisable to use a virus whose genetic backbone is highly attenuated in man. Traditionally, the A/PR/8/34 strain of virus is used, however, the genetic basis of its attenuation is unclear. Cold-adapted (CA) strains of the influenza virus are all based on the H2N2 subtype, itself a virus with pandemic potential, and again the genetic basis of temperature sensitivity is not yet established. Reverse genetics technology allows us to engineer designer influenza viruses to order. Using this technology, we have been investigating mutations in several different gene segments to effectively attenuate potential vaccine strains allowing the safe production of vaccine to protect against the next pandemic. (C) 2003 Elsevier B.V. All rights reserved.

Role of oxygen radicals in DNA damage and cancer incidence

The development of cancer in humans and animals is a multistep process. The complex series of cellular and molecular changes participating in cancer development are mediated by a diversity of endogenous and exogenous stimuli. One type of endogenous damage is that arising from intermediates of oxygen (dioxygen) reduction - oxygen-free radicals (OFR), which attacks not only the bases but also the deoxyribosyl backbone of DNA. Thanks to improvements in analytical techniques, a major achievement in the understanding of carcinogenesis in the past two decades has been the identification and quantification of various adducts of OFR with DNA. OFR are also known to attack other cellular components such as lipids, leaving behind reactive species that in turn can couple to DNA bases. Endogenous DNA lesions are genotoxic and induce mutations. The most extensively studied lesion is the formation of 8-OH-dG. This lesion is important because it is relatively easily formed and is mutagenic and therefore is a potential biomarker of carcinogenesis. Mutations that may arise from formation of 8-OH-dG involve GC. TA transversions. In view of these findings, OFR are considered as an important class of carcinogens. The effect of OFR is balanced by the antioxidant action of non-enzymatic antioxidants as well as antioxidant enzymes. Non-enzymatic antioxidants involve vitamin C, vitamin E, carotenoids (CAR), selenium and others. However, under certain conditions, some antioxidants can also exhibit a pro-oxidant mechanism of action. For example, beta-carotene at high concentration and with increased partial pressure of dioxygen is known to behave as a pro-oxidant. Some concerns have also been raised over the potentially deleterious transition metal ion-mediated (iron, copper) pro-oxidant effect of vitamin C. Clinical studies mapping the effect of preventive antioxidants have shown surprisingly little or no effect on cancer incidence. The epidemiological trials together with in vitro experiments suggest that the optimal approach is to reduce endogenous and exogenous sources of oxidative stress, rather than increase intake of anti-oxidants. In this review, we highlight some major achievements in the study of DNA damage caused by OFR and the role in carcinogenesis played by oxidatively damaged DNA. The protective effect of antioxidants against free radicals is also discussed.

The impact of genome defense on mobile elements in Microbotryum

Repeat induced point mutation (RIP), a mechanism causing hypermutation of repetitive DNA sequences in fungi, has been described as a ‘genome defense’ which functions to inactivate mobile elements and inhibit their deleterious effects on genome stability. Here we address the interactions between RIP and transposable elements in the Microbotryum violaceum species complex. Ten strains of M. violaceum, most of which belong to different species of the fungus, were all found to contain intragenomic populations of copia-like retrotransposons. Intragenomic DNA sequence variation among the copia-like elements was analyzed for evidence of RIP. Among species with RIP, there was no significant correlation between the frequency of RIP-induced mutations and inferred transposition rate based on diversity. Two strains of M. violaceum, from two different plant species but belonging to the same fungal lineage, contained copia-like elements with very low diversity, as would result from a high transposition rate, and these were also unique in showing no evidence of the hypermutation patterns indicative of the RIP genome defense. In this species, evidence of RIP was also absent from a Class II helitron-like transposable element. However, unexpectedly the absolute repetitive element load was lower than in other strains.

Dietary, lifestyle and clinico-pathological factors associated with APC mutations and promoter methylation in colorectal cancers from the EPIC-Norfolk Study

The tumour suppressor APC is the most commonly altered gene in colorectal cancer (CRC). Genetic and epigenetic alterations of APC may therefore be associated with dietary and lifestyle risk factors for CRC. Analysis of APC mutations in the extended mutation cluster region (codons 1276-1556) and APC promoter 1A methylation was performed on 185 archival CRC samples collected from participants of the European Prospective Investigation into Cancer (EPIC)-Norfolk Study, with the aim of relating these to high quality seven-day dietary and lifestyle data collected prospectively. Truncating APC mutations (APC+) and promoter 1A methylation (PM+) were identified in 43% and 23% of CRCs analysed, respectively. Distal CRCs were more likely than proximal CRCs to be APC+ or PM+ (P = 0.04). APC+ CRCs were more likely to be moderately/well differentiated and microsatellite stable than APC- CRCs (P = 0.05 and 0.03). APC+ CRC cases consumed more alcohol than their counterparts (P = 0.01) and PM+ CRC cases consumed lower levels of folate and fibre (P = 0.01 and 0.004). APC+ or PM+ CRC cases consumedhigher levels of processed meat and iron from red meat and red meat products (P=0.007 and 0.006). Specifically, CRC cases harbouring GC to AT transition mutations consumed higher levels of processed meat (35 versus 24 g/day, P = 0.04) and iron from red meat and red meat products (0.8 versus 0.6 mg/day, P = 0.05). In a logistic regression model adjusted for age, sex and cigarette smoking status, each 19g/day (1SD) increment increase in processed meat consumption was associated with cases with GC to AT mutations (OR 1.68, 95% CI 1.03-2.75). In conclusion, APC+ and PM+ CRCs may be influenced by diet and GC to AT mutations in APC are associated with processed meat consumption, suggesting a mechanistic link with dietary alkylating agents, such as N-nitroso compounds.

Pathogenic LRRK2 mutations do not alter gene expression in cell model systems or human brain tissue

Point mutations in LRRK2 cause autosomal dominant Parkinson's disease. Despite extensive efforts to determine the mechanism of cell death in patients with LRRK2 mutations, the aetiology of LRRK2 PD is not well understood. To examine possible alterations in gene expression linked to the presence of LRRK2 mutations, we carried out a case versus control analysis of global gene expression in three systems: fibroblasts isolated from LRRK2 mutation carriers and healthy, non-mutation carrying controls; brain tissue from G2019S mutation carriers and controls; and HEK293 inducible LRRK2 wild type and mutant cell lines. No significant alteration in gene expression was found in these systems following correction for multiple testing. These data suggest that any alterations in basal gene expression in fibroblasts or cell lines containing mutations in LRRK2 are likely to be quantitatively small. This work suggests that LRRK2 is unlikely to play a direct role in modulation of gene expression, although it remains possible that this protein can influence mRNA expression under pathogenic cicumstances.

Use of a LightCycler gyrA mutation assay for rapid identification of mutations conferring decreased susceptibility to ciprofloxacin in multiresistant Salmonella enterica serotype typhimurium DT104 isolates

A LightCycler-based PCR-hybridization gyrA mutation assay (GAMA) was developed to rapidly detect gyrA point mutations in multiresistant (MR) Salmonella enterica serotype Typhimurium DT104 with decreased susceptibility to ciprofloxacin (MIC, 0.25 to 1.0 mg/liter). Ninety-two isolates (49 human, 43 animal) were tested with three individual oligonucleotide probes directed against an Asp-87-to-Asn (GAC --> AAC) mutation, an Asp-87-to-Gly (GAC --> GGC) mutation, and a Ser-83-to-Phe (TCC --> TTC) mutation. Strains homologous to the probes could be distinguished from strains that had different mutations by their probe-target melting temperatures. Thirty-seven human and 30 animal isolates had an Asp-87-to-Asn substitution, 6 human and 6 animal isolates had a Ser-83-to-Phe substitution, and 5 human and 2 animal isolates had an Asp-87-to-Gly substitution. The remaining six strains all had mismatches with the three probes and therefore different gyrA mutations. The sequencing of gyrA from these six isolates showed that one human strain and two animal strains had an Asp-87-to-Tyr (GAC --> TAC) substitution and two animal strains had a Ser-83-to-Tyr (TCC --> TAC) substitution. One animal strain had no gyrA mutation, suggesting that this isolate had a different mechanism of resistance. Fifty-eight of the strains tested were indistinguishable by several different typing methods including antibiograms, pulsed-field gel gel electrophoresis, and plasmid profiling, although they could be further subdivided according to gyrA mutation. This study confirmed that MR DT104 with decreased susceptibility to ciprofloxacin from humans and food animals in England and Wales may have arisen independently against a background of clonal spread of MR DT104.