Selection-mutation balance models with epistatic selection

We present an application of birth-and-death processes on configuration spaces to a generalized mutation4 selection balance model. The model describes the aging of population as a process of accumulation of mu5 tations in a genotype. A rigorous treatment demands that mutations correspond to points in abstract spaces. 6 Our model describes an infinite-population, infinite-sites model in continuum. The dynamical equation which 7 describes the system, is of Kimura-Maruyama type. The problem can be posed in terms of evolution of states 8 (differential equation) or, equivalently, represented in terms of Feynman-Kac formula. The questions of interest 9 are the existence of a solution, its asymptotic behavior, and properties of the limiting state. In the non-epistatic 10 case the problem was posed and solved in [Steinsaltz D., Evans S.N., Wachter K.W., Adv. Appl. Math., 2005, 11 35(1)]. In our model we consider a topological space X as the space of positions of mutations and the influence of epistatic potentials

Chlorine tolerant, multilayer reverse-somosis membranes with high permeate flux and high salt rejection

A new class of high molecular weight polyethersulfone ionomers is described in which the ionic content can be varied, at will, over a very wide and fully-controllable range. A novel type of coating process enables these materials to be deposited from alcohol-type solvents as cohesive but very thin (50 – 250 nm) films on porous support-membranes, giving high-flux membranes (3.3 – 5.0 L m-2 h-1 bar-1) with very good, though not outstanding salt rejection (typically 92 - 96%). A secondary layer, of formaldehyde-cross-linked polyvinyl alcohol, can be deposited from aqueous solution on the surface of the ionomer membrane, and this layer increases salt rejection to greater than 99% without serious loss of water permeability. The final multi-layer membrane shows excellent chlorine tolerance in reverse-osmosis operation.

Development of a reverse genetics system enabling the rescue of recombinant avian influenza virus A/Turkey/England/50-92/91 (H5N1)

We previously described the use of an established reverse genetics system for the generation of recombinant human influenza A viruses from cloned cDNAs. Here, we have assembled a set of plasmids to allow recovery of the avian H5N1 influenza virus A/Turkey/England/50-92/91 entirely from cDNA. This system enables us to introduce mutations or truncations into the cDNAs to create mutant viruses altered specifically in a chosen gene. These mutant viruses can then be used in future pathogenesis studies in chickens and in studies to understand the host range restrictions of avian influenza viruses in humans.

Inference on microsatellite mutation processes in the invasive mite, Varroa destructor, using reversible jump Markov chain Monte Carlo

Varroa destructor is a parasitic mite of the Eastern honeybee Apis cerana. Fifty years ago, two distinct evolutionary lineages (Korean and Japanese) invaded the Western honeybee Apis mellifera. This haplo-diploid parasite species reproduces mainly through brother sister matings, a system which largely favors the fixation of new mutations. In a worldwide sample of 225 individuals from 21 locations collected on Western honeybees and analyzed at 19 microsatellite loci, a series of de novo mutations was observed. Using historical data concerning the invasion, this original biological system has been exploited to compare three mutation models with allele size constraints for microsatellite markers: stepwise (SMM) and generalized (GSM) mutation models, and a model with mutation rate increasing exponentially with microsatellite length (ESM). Posterior probabilities of the three models have been estimated for each locus individually using reversible jump Markov Chain Monte Carlo. The relative support of each model varies widely among loci, but the GSM is the only model that always receives at least 9% support, whatever the locus. The analysis also provides robust estimates of mutation parameters for each locus and of the divergence time of the two invasive lineages (67,000 generations with a 90% credibility interval of 35,000-174,000). With an average of 10 generations per year, this divergence time fits with the last post-glacial Korea Japan land separation. (c) 2005 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.

Thiamine is synthesized by a salvage pathway in Rhizobium leguminosarum bv. viciae strain 3841

In the absence of added thiamine, Rhizobium leguminosarum bv. viciae strain 3841 does not grow in liquid medium and forms only "pin" colonies on agar plates, which contrasts with the good growth of Sinorhizobium meliloti 1021, Mesorhizobium loti 303099, and Rhizobium etli CFN42. These last three organisms have thiCOGE genes, which are essential for de novo thiamine synthesis. While R. leguminosarum bv. viciae 3841 lacks thiCOGE, it does have thiMED. Mutation of thiM prevented formation of pin colonies on agar plates lacking added thiamine, suggesting thiamine intermediates are normally present. The putative functions of ThiM, ThiE, and ThiD are 4-methyl-5-(beta-hydroxyethyl) thiazole (THZ) kinase, thiamine phosphate pyrophosphorylase, and 4-amino-5-hydroxymethyl-2-methyl pyrimidine (HMP) kinase, respectively. This suggests that a salvage pathway operates in R. leguminosarum, and addition of HMP and THZ enabled growth at the same rate as that enabled by thiamine in strain 3841 but elicited no growth in the thiM mutant (RU2459). There is a putative thi box sequence immediately upstream of the thiM, and a gfp-mut3.1 fusion to it revealed the presence of a promoter that is strongly repressed by thiamine. Using fluorescent microscopy and quantitative reverse transcription-PCR, it was shown that thiM is expressed in the rhizosphere of vetch and pea plants, indicating limitation for thiamine. Pea plants infected by RU2459 were not impaired in nodulation or nitrogen fixation. However, colonization of the pea rhizosphere by the thiM mutant was impaired relative to that of the wild type. Overall, the results show that a thiamine salvage pathway operates to enable growth of Rhizobium leguminosarum in the rhizosphere, allowing its survival when thiamine is limiting.

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.

DNA variation in a 13-Mb region including the F9 gene: inferring the genealogical history and causal role of a hemophilia B mutation (IVS 5+13 A -> G)

About 5.5% of all UK hemophilia B patients have the base substitution IVS 5+13 A-->G as the only change in their factor (F)IX gene (F9). This generates a novel donor splice site which fits the consensus better than the normal intron 5 donor splice. Use of the novel splice site should result in a missense mutation followed by the abnormal addition of four amino acids to the patients' FIX. In order to explain the prevalence of this mutation, its genealogical history is examined. Analysis of restriction fragment length polymorphism in the 21 reference UK individuals (from different families) with the above mutation showed identical haplotypes in 19 while two differed from the rest and from each other. In order to investigate the history of the mutation and to verify that it had occurred independently more than once, the sequence variation in 1.5-kb segments scattered over a 13-Mb region including F9 was examined in 18 patients and 15 controls. This variation was then analyzed with a recently developed Bayesian approach that reconstructs the genealogy of the gene investigated while providing evidence of independent mutations that contribute disconnected branches to the genealogical tree. The method also provides minimum estimates of the age of the mutation inherited by the members of coherent trees. This revealed that 17 or 18 mutant genes descend from a founder who probably lived 450 years ago, while one patient carries an independent mutation. The independent recurrence of the IVS5+13 A-->G mutation strongly supports the conclusion that it is the cause of these patients' mild hemophilia.

Likelihood-based estimation of microsatellite mutation rates

Microsatellites are widely used in genetic analyses, many of which require reliable estimates of microsatellite mutation rates, yet the factors determining mutation rates are uncertain. The most straightforward and conclusive method by which to study mutation is direct observation of allele transmissions in parent-child pairs, and studies of this type suggest a positive, possibly exponential, relationship between mutation rate and allele size, together with a bias toward length increase. Except for microsatellites on the Y chromosome, however, previous analyses have not made full use of available data and may have introduced bias: mutations have been identified only where child genotypes could not be generated by transmission from parents' genotypes, so that the probability that a mutation is detected depends on the distribution of allele lengths and varies with allele length. We introduce a likelihood-based approach that has two key advantages over existing methods. First, we can make formal comparisons between competing models of microsatellite evolution; second, we obtain asymptotically unbiased and efficient parameter estimates. Application to data composed of 118,866 parent-offspring transmissions of AC microsatellites supports the hypothesis that mutation rate increases exponentially with microsatellite length, with a suggestion that contractions become more likely than expansions as length increases. This would lead to a stationary distribution for allele length maintained by mutational balance. There is no evidence that contractions and expansions differ in their step size distributions.

Modeling the structure of amorphous MoS3: a neutron diffraction and reverse Monte Carlo study

A model for the structure of amorphous molybdenum trisulfide, a-MoS3, has been created using reverse Monte Carlo methods. This model, which consists of chains Of MoS6 units sharing three sulfurs with each of its two neighbors and forming alternate long, nonbonded, and short, bonded, Mo-Mo separations, is a good fit to the neutron diffraction data and is chemically and physically realistic. The paper identifies the limitations of previous models based on Mo-3 triangular clusters in accounting for the available experimental data.

A facile one step synthesis of N-2 substituted 3-phenyliminoisoindolinones from N-(2-carboxybenzoyl)-anthranilic acid and the design of reverse-turn mimetics

Stirring of N-(2-carboxybenzoyl) anthranilic acid with anilines and amines such as p-toluidine, benzylamine, methyl esters of Leu, Phe, Ile and Val in presence of DCC produces N- 2 substituted 3-phenyliminoisoindolinones in very good yields. Single crystal X-ray diffraction studies and solution phase NMR and CD studies reveal that the 3-phenyliminoisoindolinone moiety is a turn-inducing scaffold which should be useful for reverse-turn mimetics.

Enzyme glycation influences product yields during oligosaccharide synthesis by reverse hydrolysis

Possible evidence is presented for Maillard glycation of enzymes during oligosaccharide synthesis by reverse hydrolysis. In 70% (w/v) mannose solutions, 1,2-alpha-mannosidase from Penicillium citrinum lost 40% and alpha-mannosidase from almonds lost 60% activity at 55 degreesC over 2 weeks. Oligosaccharide yields were 15 and 45% respectively. Higher molecular weight glycation adducts were formed in a time-dependent manner as seen by MALDI-TOF. Inhibitors of the Maillard. reaction were able to partially alleviate these effects resulting in reduced loss of enzyme activity and oligosaccharide yield increases of 27-53% relative to the control. (C) 2004 Elsevier B.V. All rights reserved.

Synthesis of FimH receptor-active manno-oligosaccharides by reverse hydrolysis using alpha-mannosidases from Penicillium citrinum, Aspergillus phoenicis and almond

Recombinant Penicillium citrinum alpha-1,2-mannosidase, expressed in Aspergillus oryzae, was employed to carry out regioselective synthesis of alpha-D-mannopyranosyl-(1-->2)-D-mannose. Yields (w/w) of 16.68% disaccharide, 3.07% trisaccharide and 0.48% tetrasaccharide were obtained, with alpha1-->2 linkages present at 98.5% of the total linkages formed. Non-specific alpha-mannosidase from almond was highly efficient in reverse hydrolysis and oligosaccharide yields of 45-50% were achieved. The products of the almond mannosidase were a mixture of disaccharides (30.75%, w/w), trisaccharides (12.26%, w/w) and tetrasaccharides (1.89%, w/w) with 1-->2, 1-->3 and 1-->6 isomers. alpha-1,2-linkage specific mannosidase from P. citrinum and alpha-1,6-linkage-specific mannosidase from Aspergillus phoenicis were used in combination to hydrolyse the respective linkages from the mixture of isomers, resulting in alpha-D-mannopyranosyl-(1-->3)-D-mannose in 86.4% purity. The synthesised oligosaccharides can potentially inhibit the adhesion of pathogens by acting as 'decoys' of receptors of type-1 fimbriae carried by enterobacteria.