FTIR and UV spectroscopy studies of triplex formation between alpha-oligonucleotides with non-ionic phoshoramidate linkages and DNA targets

The triplexes formed by pyrimidine alpha-oligodeoxynucleotides, 15mers alpha dT(15) or 12mers alpha dCT having dimethoxyethyl (PNHdiME), morpholino (PMOR) or propyl (PNHPr) non-ionic phosphoramidate linkages with DNA duplex targets have been investigated by UV and FTIR spectroscopy. Due to the decrease in the electrostatic repulsion between partner strands of identical lengths all modifications result in triplexes more stable than those formed with unmodified phosphodiester beta-oligodeoxynucleotides (beta-ODNs). Among the alpha-ODN third strands having C and T bases and non-ionic phosphoramidate linkages (alpha dCTPN) the most efficient modification is (PNHdiME). The enhanced third strand stability of the alpha dCTPN obtained as diastereoisomeric mixtures is attenuated by the steric hindrance of the PMOR linkages or by the hydrophobicity of the PNHPr linkages. All alpha dCTPN strands form triplexes even at neutral pH. In the most favorable case (PNHdiME), we show by FTIR spectroscopy that the triplex formed at pH 7 is held by Hoogsteen T*A.T triplets and in addition by an hydrogen bond between O6 of G and C of the third strand (Tm = 30 degrees C). The detection of protonated cytosines is correlated at pH 6 with a high stabilization of the triplex (Tm = 65 degrees C). While unfavorable steric effects are overcome with alpha anomers, the limitation of the pH dependence is not completely suppressed. Different triplexes are evidenced for non pH dependent phosphoramidate alpha-thymidilate strands (alpha dT(15)PN) interacting with a target duplex of identical length. At low ionic strength and DNA concentration we observe the binding to beta dA(15) either of alpha dT(15)PN as duplex strand and beta dT(15) as third strand, or of two hydrophobic alpha dT(15)PNHPr strands. An increase in the DNA and counterion concentration stabilizes the anionic target duplex and then the alpha dT(15)PN binds as Hoogsteen third strand.     

A selective sweep driven by pyrimethamine treatment in southeast asian malaria parasites

Malaria parasites (Plasmodium falciparum) provide an excellent system in which to study the genomic effects of strong selection in a recombining eukaryote because the rapid spread of resistance to multiple drugs during the last the past 50 years has been well documented, the full genome sequence and a microsatellite map are now available, and haplotype data can be easily generated. We examined microsatellite variation around the dihydrofolate reductase (dhfr) gene on chromosome 4 of P. falciparum. Point mutations in dhfr are known to be responsible for resistance to the antimalarial drug pyrimethamine, and resistance to this drug has spread rapidly in Southeast (SE) Asia after its introduction in 1970s. We genotyped 33 microsatellite markers distributed across chromosome 4 in 61 parasites from a location on the Thailand/Myanmar border. We observed minimal microsatellite length variation in a 12-kb (0.7-cM) region flanking the dhfr gene and diminished variation for approximately 100 kb (6 cM), indicative of a single origin of resistant alleles. Furthermore, we found the same or similar microsatellite haplotypes flanked resistant dhfr alleles sampled from 11 parasite populations in five SE Asian countries indicating recent invasion of a single lineage of resistant dhfr alleles in locations 2000 km apart. Three features of these data are of especially interest. (1). Pyrimethamine resistance is generally assumed to have evolved multiple times because the genetic basis is simple and resistance can be selected easily in the laboratory. Yet our data clearly indicate a single origin of resistant dhfr alleles sampled over a large region of SE Asia. (2). The wide valley ( approximately 6 cM) of reduced variation around dhfr provides "proof-of-principle" that genome-wide association may be an effective way to locate genes under strong recent selection. (3). The width of the selective valley is consistent with predictions based on independent measures of recombination, mutation, and selection intensity, suggesting that we have reasonable estimates of these parameters. We conclude that scanning the malaria parasite genome for evidence of recent selection may prove an extremely effective way to locate genes underlying recently evolved traits such as drug resistance, as well as providing an opportunity to study the dynamics of selective events that have occurred recently or are currently in progress.     

Bayes or bootstrap? A simulation study comparing the performance of Bayesian Markov chain Monte Carlo sampling and bootstrapping in assessing phylogenetic confidence

Bayesian Markov chain Monte Carlo sampling has become increasingly popular in phylogenetics as a method for both estimating the maximum likelihood topology and for assessing nodal confidence. Despite the growing use of posterior probabilities, the relationship between the Bayesian measure of confidence and the most commonly used confidence measure in phylogenetics, the nonparametric bootstrap proportion, is poorly understood. We used computer simulation to investigate the behavior of three phylogenetic confidence methods: Bayesian posterior probabilities calculated via Markov chain Monte Carlo sampling (BMCMC-PP), maximum likelihood bootstrap proportion (ML-BP), and maximum parsimony bootstrap proportion (MP-BP). We simulated the evolution of DNA sequence on 17-taxon topologies under 18 evolutionary scenarios and examined the performance of these methods in assigning confidence to correct monophyletic and incorrect monophyletic groups, and we examined the effects of increasing character number on support value. BMCMC-PP and ML-BP were often strongly correlated with one another but could provide substantially different estimates of support on short internodes. In contrast, BMCMC-PP correlated poorly with MP-BP across most of the simulation conditions that we examined. For a given threshold value, more correct monophyletic groups were supported by BMCMC-PP than by either ML-BP or MP-BP. When threshold values were chosen that fixed the rate of accepting incorrect monophyletic relationship as true at 5%, all three methods recovered most of the correct relationships on the simulated topologies, although BMCMC-PP and ML-BP performed better than MP-BP. BMCMC-PP was usually a less biased predictor of phylogenetic accuracy than either bootstrapping method. BMCMC-PP provided high support values for correct topological bipartitions with fewer characters than was needed for nonparametric bootstrap.     

Extreme divergence of a homeotic gene: the bicoid case

Evolutionary developmental genetics (evo-devo) reveals that the plasticity of development is so important that every developmental biology project should carefully take this point into consideration. The example of bicoid, the first discovered morphogen, illustrates how an essential gene can change its function during evolution. The search for bicoid homologues showed that this gene is surprisingly specific to flies (cyclorraphan diptera) and absent in other insects. In fact, recent studies demonstrate that bicoid is a very derived Hox3 homeotic gene. During insect evolution, the ancestral Hox3 gene lost its homeotic function and acquired new roles in oocytes and embryonic annexes. Then, in the lineage leading to modern flies, a duplication of this new gene, followed by functional divergence, led to the formation of bicoid and zerknüllt. Both genes are located within the Drosophila Hox complex; however, they have no homeotic function. Thanks to the power of Drosophila genetics, it is possible to suggest that torso and hunchback may constitute the insect primitive anterior organizer. The bicoid evolutionary history reveals several fundamental mechanisms of the evolution of developmental genes, such as changes of gene regulation, modifications of protein sequences and gene duplication. It also shows the need for studying a wider range of model organisms before generalisations can be made from data obtained with one particular species.     

Influence of mutational and sampling factors on the estimation of demographic parameters in a "continuous" population under isolation by distance

In numerous species, individual dispersal is restricted in space so that "continuous" populations evolve under isolation by distance. A method based on individual genotypes assuming a lattice population model was recently developed to estimate the product Dsigma2, where D is the population density and sigma2 is the average squared parent-offspring distance. We evaluated the influence on this method of both mutation rate and mutation model, with a particular reference to microsatellite markers, as well as that of the spatial scale of sampling. Moreover, we developed and tested a nonparametric bootstrap procedure allowing the construction of confidence intervals for the estimation of Dsigma2. These two objectives prompted us to develop a computer simulation algorithm based on the coalescent theory giving individual genotypes for a continuous population under isolation by distance. Our results show that the characteristics of mutational processes at microsatellite loci, namely the allele size homoplasy generated by stepwise mutations, constraints on allele size, and change of slippage rate with repeat number, have little influence on the estimation of Dsigma2. In contrast, a high genetic diversity (approximately 0.7-0.8), as is commonly observed for microsatellite markers, substantially increases the precision of the estimation. However, very high levels of genetic diversity (>0.85) were found to bias the estimation. We also show that statistics taking into account allele size differences give unreliable estimations (i.e., high variance of Dsigma2 estimation) even under a strict stepwise mutation model. Finally, although we show that this method is reasonably robust with respect to the sampling scale, sampling individuals at a local geographical scale gives more precise estimations of Dsigma2.     

Memory in the grip of emotions and stress: a necessarily harmful impact?

While intense negative events are vividly recalled, information learned during stressful situations is poorly remembered. These differential effects of emotions and stress on memory have been attributed to the physiological manifestations generated during those affective states. Intense emotional and stressful events trigger the secretion of catecholamines and of glucocorticoids, in particular. These hormones would be modulatory agents of memory functions. In the first part of this paper, we review the specific effects emotions and stress have on memory. We then summarize the psychological and biological determinants responsible for these effects. Finally, we discuss different methodological issues that could explain the discrepancy found between the impact of emotions and stress on memory. Defining more precisely the effects emotion and stress have on memory will lead to a better comprehension of the cognitive problems that characterize patients dealing with emotional turmoil, such as patients suffering from depression or post-traumatic stress disorder.     

Role of vascularization in determining the time of hypoxic-ischemic encephalopathy in the neonate

OBJECTIVE: To examine the role of vascularization in determining the time of hypoxic-ischemic encephalopathy (HIEP). STUDY DESIGN: Brain sections of 126 neonatal autopsy cases were examined for edema, gliosis, congestion, inflammation and ischemia. Capillary vessels were examined with both reticulum stains and antibody against CD34. Vascular surface density (VSD) and number of vessels per stroma (NVES) were calculated by stereologic methods. RESULTS: Among 126 cases, 64 were male (50.8%) and 62 female (49.2%). In 25 cases HIEP was observed; 14 had a pregnancy history of hypertension, eclampsia or diabetes mellitus in the mother, with fetal distress or underdeveloped features. Statistically, NVES was strongly related to primary HIEP. However, the HIEP and non-HIEP cases revealed no differences in NVES and VSD means. CONCLUSION: Vascularization, especially NVES, helps in determining whether an HIEP case is pregnancy related or due to end-stage changes of dying, but is not an indicator of HIEP.     

Characterization of carbonic anhydrases from Riftia pachyptila,a symbiotic invertebrate from deep-sea hydrothermal vents

The symbiotic hydrothermal vent tubeworm Riftia pachyptila needs to supply its internal bacterial symbionts with carbon dioxide, their inorganic carbon source. Our aim in this study was to characterize the carbonic anhydrase (CA) involved in CO(2) transport and conversion at various steps in the plume and the symbiotic tissue, the trophosome. A complete 1209 kb cDNA has been sequenced from the trophosome and identified as a putative alpha-CA based on BLAST analysis and the similarities of total deduced amino-acid sequence with those from the GenBank database. In the plume, the putative CA sequence obtained from cDNA library screening was 90% identical to the trophosome CA, except in the first 77 nucleotides downstream from the initiation site identified on trophosome CA. A phylogenetic analysis showed that the annelidan Riftia CA (CARp) emerges clustered with invertebrate CAs, the arthropodan Drosophila CA and the cnidarian Anthopleura CA. This invertebrate cluster appeared as a sister group of the cluster comprising mitochondrial and cytosolic isoforms in vertebrates: CAV, CAI II and III, and CAVII. However, amino acid sequence alignment showed that Riftia CA was closer to cytosolic CA than to mitochondrial CA. Combined biochemical approaches revealed two cytosolic CAs with different molecular weights and pI's in the plume and the trophosome, and the occurrence of a membrane-bound CA isoform in addition to the cytosolic one in the trophosome. The physiologic roles of cytosolic CA in both tissues and supplementary membrane-bound CA isoform in the trophosome in the optimization of CO(2) transport and conversion are discussed.     

Electrostatic properties in the catalytic site of papain: A possible regulatory mechanism for the reactivity of the ion pair

We present an analysis of the electrostatic properties in the catalytic site of papain (EC 3.4.22.2), an archetype enzyme of the C1 cysteine proteinase family, and we investigate their possible role in the formation, stabilization and regulation of the Cys25((-))...His159((+)) catalytic ion pair. The electrostatic properties were computed using a reassociation method based in multicentered multipolar expansions obtained from ab initio quantum calculations of overlapping protein fragments. Solvent effects were introduced by coupling the use of multicentered multipolar expansions to two continuum boundary element methods to solve the Poisson and the linearized Poisson-Boltzmann equations. The electrostatic profile found in the proton transfer region of papain showed that this enzyme has a well-defined electrostatic environment to favor the formation and stabilization of the catalytic ion pair. The papain catalytic site electrostatic profile can be considered as an electrostatic fingerprint of the papain family with the following characteristics: (i) the presence of a net electric field highly aligned in the (Cys25)-SG-->(His159)-ND1 direction; (ii) the electrostatic profile has a saddle-point character; (iii) it is basically a local environmental effect. Furthermore, our analysis describes a possible regulatory mechanism (the E(SG-->ND1) attenuation effect) controlling the ion pair reactivity and permits to infer the Asp57 acidic residue as the most probable candidate to act as the electrostatic modulator.