Restriction site-associated DNA sequencing, genotyping error estimation and de novo assembly optimization for population genetic inference.

Restriction site-associated DNA sequencing (RADseq) provides researchers with the ability to record genetic polymorphism across thousands of loci for nonmodel organisms, potentially revolutionizing the field of molecular ecology. However, as with other genotyping methods, RADseq is prone to a number of sources of error that may have consequential effects for population genetic inferences, and these have received only limited attention in terms of the estimation and reporting of genotyping error rates. Here we use individual sample replicates, under the expectation of identical genotypes, to quantify genotyping error in the absence of a reference genome. We then use sample replicates to (i) optimize de novo assembly parameters within the program Stacks, by minimizing error and maximizing the retrieval of informative loci; and (ii) quantify error rates for loci, alleles and single-nucleotide polymorphisms. As an empirical example, we use a double-digest RAD data set of a nonmodel plant species, Berberis alpina, collected from high-altitude mountains in Mexico.

High-resolution magnetic resonance imaging quantitatively detects individual pancreatic islets.

OBJECTIVE-We studied whether manganese-enhanced high-field magnetic resonance (MR) imaging (MEHFMRI) could quantitatively detect individual islets in situ and in vivo and evaluate changes in a model of experimental diabetes.RESEARCH DESIGN AND METHODS-Whole pancreata from untreated (n = 3), MnCl(2) and glucose-injected mice (n = 6), and mice injected with either streptozotocin (STZ; n = 4) or citrate buffer (n = 4) were imaged ex vivo for unambiguous evaluation of islets. Exteriorized pancreata of MnCl(2) and glucose-injected mice (n = 6) were imaged in vivo to directly visualize the gland and minimize movements. In all cases, MR images were acquired in a 14.1 Testa scanner and correlated with the corresponding (immuno)histological sections.RESULTS-In ex vivo experiments, MEHFMRI distinguished different pancreatic tissues and evaluated the relative abundance of islets in the pancreata of normoglycemic mice. MEHFMRI also detected a significant decrease in the numerical and volume density of islets in STZ-injected mice. However, in the latter measurements the loss of beta-cells was undervalued under the conditions tested. The experiments on the externalized pancreata confirmed that MEHFMRI could visualize native individual islets in living, anesthetized mice.CONCLUSIONS-Data show that MEHFMRI quantitatively visualizes individual islets in the intact mouse pancreas, both ex vivo and in vivo. Diabetes 60:2853-2860, 2011

Gene flow among wild and domesticated almond species: insights from chloroplast and nuclear markers

Hybridization has played a central role in the evolutionary history of domesticated plants. Notably, several breeding programs relying on gene introgression from the wild compartment have been performed in fruit tree species within the genus Prunus but few studies investigated spontaneous gene flow among wild and domesticated Prunus species. Consequently, a comprehensive understanding of genetic relationships and levels of gene flow between domesticated and wild Prunus species is needed. Combining nuclear and chloroplastic microsatellites, we investigated the gene flow and hybridization among two key almond tree species, the cultivated Prunus dulcis and one of the most widespread wild relative Prunus orientalis in the Fertile Crescent. We detected high genetic diversity levels in both species along with substantial and symmetric gene flow between the domesticated P. dulcis and the wild P. orientalis. These results were discussed in light of the cultivated species diversity, by outlining the frequent spontaneous genetic contributions of wild species to the domesticated compartment. In addition, crop-to-wild gene flow suggests that ad hoc transgene containment strategies would be required if genetically modified cultivars were introduced in the northwestern Mediterranean.

In silico prediction of human carboxylesterase-1 (hCES1) metabolism combining docking analyses and MD simulations.

Metabolic problems lead to numerous failures during clinical trials, and much effort is now devoted in developing in silico models predicting metabolic stability and metabolites. Such models are well known for cytochromes P450 and some transferases, whereas little has been done to predict the hydrolytic activity of human hydrolases. The present study was undertaken to develop a computational approach able to predict the hydrolysis of novel esters by human carboxylesterase hCES1. The study involves both docking analyses of known substrates to develop predictive models, and molecular dynamics (MD) simulations to reveal the in situ behavior of substrates and products, with particular attention being paid to the influence of their ionization state. The results emphasize some crucial properties of the hCES1 catalytic cavity, confirming that as a trend with several exceptions, hCES1 prefers substrates with relatively smaller and somewhat polar alkyl/aryl groups and larger hydrophobic acyl moieties. The docking results underline the usefulness of the hydrophobic interaction score proposed here, which allows a robust prediction of hCES1 catalysis, while the MD simulations show the different behavior of substrates and products in the enzyme cavity, suggesting in particular that basic substrates interact with the enzyme in their unprotonated form.

Development of a method to investigate the hydrolysis of xenobiotic esters by a Mycobacterium smegmatis homogenate.

One of the main problems in combating tuberculosis is caused by a poor penetration of drugs into the mycobacterial cells. A prodrug approach via activation inside mycobacterial cells is a possible strategy to overcome this hurdle and achieve efficient drug uptake. Esters are attractive candidates for such a strategy and we and others communicated previously the activity of esters of weak organic acids against mycobacteria. However very little is known about ester hydrolysis by mycobacteria and no biological model is available to study the activation of prodrugs by these microorganisms. To begin filling this gap, we have embarked in a project to develop an in vitro method to study prodrug activation by mycobacteria using Mycobacterium smegmatis homogenates. Model ester substrates were ethyl nicotinate and ethyl benzoate whose hydrolysis was monitored and characterized kinetically. Our studies showed that in M. smegmatis most esterase activity is associated with the soluble fraction (cytosol) and is preserved by storage at 5°C or at room temperature for one hour, or by storage at -80°C up to one year. In the range of homogenate concentrations studied (5-80% in buffer), k(obs) varied linearly with homogenate concentration for both substrates. We also found that the homogenates showed Michaelis-Menten kinetics behavior with both prodrugs. Since ethyl benzoate is a good substrate for the mycobacterial esterases, this compound can be used to standardize the esterasic activity of homogenates, allowing results of incubations of prodrugs with homogenates from different batches to be readily compared.

The complex history of the olive tree: from Late Quaternary diversification of Mediterranean lineages to primary domestication in the northern Levant.

The location and timing of domestication of the olive tree, a key crop in Early Mediterranean societies, remain hotly debated. Here, we unravel the history of wild olives (oleasters), and then infer the primary origins of the domesticated olive. Phylogeography and Bayesian molecular dating analyses based on plastid genome profiling of 1263 oleasters and 534 cultivated genotypes reveal three main lineages of pre-Quaternary origin. Regional hotspots of plastid diversity, species distribution modelling and macrofossils support the existence of three long-term refugia; namely the Near East (including Cyprus), the Aegean area and the Strait of Gibraltar. These ancestral wild gene pools have provided the essential foundations for cultivated olive breeding. Comparison of the geographical pattern of plastid diversity between wild and cultivated olives indicates the cradle of first domestication in the northern Levant followed by dispersals across the Mediterranean basin in parallel with the expansion of civilizations and human exchanges in this part of the world.

Restriction site-associated DNA sequencing, genotyping error estimation and de novo assembly optimization for population genetic inference.

Restriction site-associated DNA sequencing (RADseq) provides researchers with the ability to record genetic polymorphism across thousands of loci for nonmodel organisms, potentially revolutionizing the field of molecular ecology. However, as with other genotyping methods, RADseq is prone to a number of sources of error that may have consequential effects for population genetic inferences, and these have received only limited attention in terms of the estimation and reporting of genotyping error rates. Here we use individual sample replicates, under the expectation of identical genotypes, to quantify genotyping error in the absence of a reference genome. We then use sample replicates to (i) optimize de novo assembly parameters within the program Stacks, by minimizing error and maximizing the retrieval of informative loci; and (ii) quantify error rates for loci, alleles and single-nucleotide polymorphisms. As an empirical example, we use a double-digest RAD data set of a nonmodel plant species, Berberis alpina, collected from high-altitude mountains in Mexico.

Decoupled post-glacial history in mutualistic plant-insect interactions: insights from the yellow loosestrife (Lysimachia vulgaris) and its associated oil-collecting bees (Macropis europaea and M-fulvipes)

AimWe take a comparative phylogeographical approach to assess whether three species involved in a specialized oil-rewarding pollination system (i.e. Lysimachia vulgaris and two oil-collecting bees within the genus Macropis) show congruent phylogeographical trajectories during post-glacial colonization processes. Our working hypothesis is that within specialized mutualistic interactions, where each species relies on the co-occurrence of the other for survival and/or reproduction, partners are expected to show congruent evolutionary trajectories, because they are likely to have followed parallel migration routes and to have shared glacial refugia. LocationWestern Palaearctic. MethodsOur analysis relies on the extensive sampling of 104 Western Palaearctic populations (totalling 434, 159 and 74 specimens of Lysimachiavulgaris, Macropiseuropaea and Macropisfulvipes, respectively), genotyped with amplified fragment length polymorphism. Based on this, we evaluated the regional genetic diversity (Shannon diversity and allele rarity index) and genetic structure (assessed using structure, population networks, isolation-by-distance and spatial autocorrelation metrics) of each species. Finally, we compared the general phylogeographical patterns obtained. ResultsContrary to our expectations, the analyses revealed phylogeographical signals suggesting that the investigated organisms demonstrate independent post-glacial trajectories as well as distinct contemporaneous demographic parameters, despite their mutualistic interaction. Main conclusionsThe mutualistic partners investigated here are likely to be experiencing distinct and independent evolutionary dynamics because of their contrasting life-history traits (e.g. dispersal abilities), as well as distinct hubs and migration routes. Such conditions would prevent and/or erase any signature of co-structuring of lineages in space and time. As a result, the lack of phylogeographical congruence driven by differences in life-history traits might have arisen irrespective of the three species having shared similar Pleistocene glacial refugia.

Does a shift in host plants trigger speciation in the Alpine leaf beetle Oreina speciosissima (Coleoptera, Chrysomelidae)?

Background: Within the Coleoptera, the largest order in the animal kingdom, the exclusively herbivorous Chrysomelidae are recognized as one of the most species rich beetle families. The evolutionary processes that have fueled radiation into the more than thirty-five thousand currently recognized leaf beetle species remain partly unresolved. The prominent role of leaf beetles in the insect world, their omnipresence across all terrestrial biomes and their economic importance as common agricultural pest organisms make this family particularly interesting for studying the mechanisms that drive diversification. Here we specifically focus on two ecotypes of the alpine leaf beetle Oreina speciosissima (Scop.), which have been shown to exhibit morphological differences in male genitalia roughly corresponding to the subspecies Oreina speciosissima sensu stricto and Oreina speciosissima troglodytes. In general the two ecotypes segregate along an elevation gradient and by host plants: Oreina speciosissima sensu stricto colonizes high forb vegetation at low altitude and Oreina speciosissima troglodytes is found in stone run vegetation at higher elevations. Both host plants and leaf beetles have a patchy geographical distribution. Through use of gene sequencing and genome fingerprinting (AFLP) we analyzed the genetic structure and habitat use of Oreina speciosissima populations from the Swiss Alps to examine whether the two ecotypes have a genetic basis. By investigating a wide range of altitudes and focusing on the structuring effect of habitat types, we aim to provide answers regarding the factors that drive adaptive radiation in this phytophagous leaf beetle.Results: While little phylogenetic resolution was observed based on the sequencing of four DNA regions, the topology and clustering resulting from AFLP genotyping grouped specimens according to their habitat, mostly defined by plant associations. A few specimens with intermediate morphologies clustered with one of the two ecotypes or formed separate clusters consistent with habitat differences. These results were discussed in an ecological speciation framework.Conclusions: The question of whether this case of ecological differentiation occurred in sympatry or allopatry remains open. Still, the observed pattern points towards ongoing divergence between the two ecotypes which is likely driven by a recent shift in host plant use.

Homology modeling and metabolism prediction of human carboxylesterase-2 using docking analyses by GriDock: a parallelized tool based on AutoDock 4.0.

Metabolic problems lead to numerous failures during clinical trials, and much effort is now devoted to developing in silico models predicting metabolic stability and metabolites. Such models are well known for cytochromes P450 and some transferases, whereas less has been done to predict the activity of human hydrolases. The present study was undertaken to develop a computational approach able to predict the hydrolysis of novel esters by human carboxylesterase hCES2. The study involved first a homology modeling of the hCES2 protein based on the model of hCES1 since the two proteins share a high degree of homology (congruent with 73%). A set of 40 known substrates of hCES2 was taken from the literature; the ligands were docked in both their neutral and ionized forms using GriDock, a parallel tool based on the AutoDock4.0 engine which can perform efficient and easy virtual screening analyses of large molecular databases exploiting multi-core architectures. Useful statistical models (e.g., r (2) = 0.91 for substrates in their unprotonated state) were calculated by correlating experimental pK(m) values with distance between the carbon atom of the substrate's ester group and the hydroxy function of Ser228. Additional parameters in the equations accounted for hydrophobic and electrostatic interactions between substrates and contributing residues. The negatively charged residues in the hCES2 cavity explained the preference of the enzyme for neutral substrates and, more generally, suggested that ligands which interact too strongly by ionic bonds (e.g., ACE inhibitors) cannot be good CES2 substrates because they are trapped in the cavity in unproductive modes and behave as inhibitors. The effects of protonation on substrate recognition and the contrasting behavior of substrates and products were finally investigated by MD simulations of some CES2 complexes.

Passive lipoidal diffusion and carrier-mediated cell uptake are both important mechanisms of membrane permeation in drug disposition.

Recently, it has been proposed that drug permeation is essentially carrier-mediated only and that passive lipoidal diffusion is negligible. This opposes the prevailing hypothesis of drug permeation through biological membranes, which integrates the contribution of multiple permeation mechanisms, including both carrier-mediated and passive lipoidal diffusion, depending on the compound's properties, membrane properties, and solution properties. The prevailing hypothesis of drug permeation continues to be successful for application and prediction in drug development. Proponents of the carrier-mediated only concept argue against passive lipoidal diffusion. However, the arguments are not supported by broad pharmaceutics literature. The carrier-mediated only concept lacks substantial supporting evidence and successful applications in drug development.

The bioavailability of bromazepam, omeprazole and paracetamol given by nasogastric feeding tube.

AIMS: To characterize and compare the pharmacokinetic profiles of bromazepam, omeprazole and paracetamol when administered by the oral and nasogastric routes to the same healthy cohort of volunteers. METHODS: In a prospective, monocentric, randomized crossover study, eight healthy volunteers received the three drugs by the oral (OR) and nasogastric routes (NT). Sequential plasma samples were analyzed by high-performance liquid chromatography-UV, pharmacokinetic parameters (Cmax, AUC(0-infinity), t(1/2), k(e), tmax) were compared statistically, and Cmax, AUC(0-infinity) and t(max) were analyzed for bioequivalence. RESULTS: A statistically significant difference was seen in the AUC(0-infinity) of bromazepam, with nasogastric administration decreasing availability by about 25%: AUC(OR) = 2501 ng mL(-1) h; AUC(NT) = 1855 ng mL(-1) h (p < 0.05); ratio (geometric mean) = 0.74 [90% confidence interval (CI) 0.64-0.87]. However, this does not appear to be clinically relevant given the usual dosage range and the drug's half-life (approx. 30 h). A large interindividual variability in omeprazole parameters prevented any statistical conclusion from being drawn in terms of both modes of administration despite their similar average profile: AUC(OR) = 579 ng mL(-1) h; AUC(NT) = 587 ng mL(-1) h (p > 0.05); ratio (geometric mean) = 1.01 (90% CI 0.64-1.61). An extended study with a larger number of subjects may possibly provide clearer answers. The narrow 90% confidence limits of paracetamol indicate bioequivalence: AUC(OR) = 37 microg mL(-1) h; AUC(NT) = 41 microg mL(-1) h(p > 0.05); ratio (geometric mean) = 1.12 (90% CI 0.98-1.28). CONCLUSION: The results of this study show that the nasogastric route of administration does not appear to cause marked, clinically unsuitable alterations in the bioavailability of the tested drugs.