Conflicting mitochondrial and nuclear phylogeographic signals and evolution of host-plant shifts in the boreo-montane leaf beetle Chrysomela lapponica

We conducted a phylogeographic study on the cold-adapted leaf beetle Chrysomela lapponica, that feeds on willow or birch, by sampling several populations throughout most of the geographic distribution of the species, and by sequencing for each individual one mitochondrial and two nuclear DNA fragments. Patterns of DNA sequence variation from the mitochondrial and nuclear loci, as displayed in the median-joining networks, appear to display contradicting historical signal: a deep genealogical divergence is observed with the mitochondrial genome between the Alpine population and all other populations found in the Euro-Siberian distribution of the species, that is completely absent with both nuclear loci. We use coalescence simulations of DNA sequence evolution to test the hypothesis that this apparent conflict is compatible with a neutral model of sequence evolution (i.e., to check whether the stochastic nature of the coalescence process can explain these patterns). Because the simulations show that this is highly unlikely, we consider two alternative hypotheses: (1) introgression of the mitochondrial genome of another species and (2) the effect of natural selection. Although introgression is the most plausible explanation, we fail to identify the source species of the introgressed mitochondrial genome among all known species closely related to C. lapponica. We therefore suggest that the putative introgression event is ancient and the source species is either extinct or currently outside the geographic range of C. lapponica explored in this study. The observed DNA sequence variation also suggests that a host-plant shift from willow to birch has occurred recently and independently in each of the three birch-feeding populations. This emphasizes further the relative ease with which these beetles can escape their ancestral host-plant specialization on willow, but shows at the same time that host-plant shifts are highly constrained, as they only occur between willow and birch.     

Molecular genetics of the transcription factor GLIS3 identifies its dual function in beta cells and neurons

The GLIS family zinc finger 3 isoform (GLIS3) is a risk gene for Type 1 and Type 2 diabetes, glaucoma and Alzheimer's disease endophenotype. We identified GLIS3 binding sites in insulin secreting cells (INS1) (FDR q < 0.05; enrichment range 1.40–9.11 fold) sharing the motif wrGTTCCCArTAGs, which were enriched in genes involved in neuronal function and autophagy and in risk genes for metabolic and neuro-behavioural diseases. We confirmed experimentally Glis3-mediated regulation of the expression of genes involved in autophagy and neuron function in INS1 and neuronal PC12 cells. Naturally-occurring coding polymorphisms in Glis3 in the Goto-Kakizaki rat model of type 2 diabetes were associated with increased insulin production in vitro and in vivo, suggestive alteration of autophagy in PC12 and INS1 and abnormal neurogenesis in hippocampus neurons. Our results support biological pleiotropy of GLIS3 in pathologies affecting β-cells and neurons and underline the existence of trans?nosology pathways in diabetes and its co-morbidities.     

Colitis induced by proteinase-activated receptor-2 agonists is mediated by a neurogenic mechanism

Proteinase-activated receptor-2 (PAR2) activation induces colonic inflammation by an unknown mechanism. We hypothesized that PAR2 agonists administered intracolonically in mice induce inflammation via a neurogenic mechanism. Pretreatment of mice with neurokinin-1 and calcitonin-gene-related peptide (CGRP) receptor antagonists or with capsaicin showed attenuated PAR2-agonist-induced colitis. Immunohistochemistry demonstrated a differential expression of a marker for the type-1 CGRP receptor during the time course of PAR2-agonist-induced colitis, further suggesting a role for CGRP. We conclude that PAR2-agonist-induced intestinal inflammation involves the release of neuropeptides, which by acting on their receptors cause inflammation. These results implicate PAR2 as an important mediator of intestinal neurogenic inflammation.     

The relationship between postpartum vaginal discharge symptoms and progesterone profile characteristics in lactating dairy cows in Denmark

In this paper, the effect of clinical symptoms of uterine inflammation on progesterone profile characteristics was quantified in dairy cows. A continuous scale based on visual observation of vaginal discharge (the previously developed D-index) was used to describe the clinical symptoms. Progesterone profiles in milk were used to describe the ovarian cycles, and to determine the distinguishing features of these profiles, a multivariate statistical procedure (principal component analysis) was performed.Significant negative effects of the D-index were seen during the first and second postpartum ovarian cycles. The D-index had a significant effect on the shape of progesterone profiles and the length of the ovarian cycles but it only accounted for a small proportion of the variation in these ovarian cycle features. The D-index was not a significant risk factor for the length of postpartum anovulatory period in the present study.     

Distinguishing Antimicrobial Models with Different Resistance Mechanisms via Population Pharmacodynamic Modeling

Semi-mechanistic pharmacokinetic-pharmacodynamic (PK-PD) modeling is increasingly used for antimicrobial drug development and optimization of dosage regimens, but systematic simulation-estimation studies to distinguish between competing PD models are lacking. This study compared the ability of static and dynamic in vitro infection models to distinguish between models with different resistance mechanisms and support accurate and precise parameter estimation. Monte Carlo simulations (MCS) were performed for models with one susceptible bacterial population without (M1) or with a resting stage (M2), a one population model with adaptive resistance (M5), models with pre-existing susceptible and resistant populations without (M3) or with (M4) inter-conversion, and a model with two pre-existing populations with adaptive resistance (M6). For each model, 200 datasets of the total bacterial population were simulated over 24h using static antibiotic concentrations (256-fold concentration range) or over 48h under dynamic conditions (dosing every 12h; elimination half-life: 1h). Twelve-hundred random datasets (each containing 20 curves for static or four curves for dynamic conditions) were generated by bootstrapping. Each dataset was estimated by all six models via population PD modeling to compare bias and precision. For M1 and M3, most parameter estimates were unbiased (<10%) and had good imprecision (<30%). However, parameters for adaptive resistance and inter-conversion for M2, M4, M5 and M6 had poor bias and large imprecision under static and dynamic conditions. For datasets that only contained viable counts of the total population, common statistical criteria and diagnostic plots did not support sound identification of the true resistance mechanism. Therefore, it seems advisable to quantify resistant bacteria and characterize their MICs and resistance mechanisms to support extended simulations and translate from in vitro experiments to animal infection models and ultimately patients.     

Signature of selective sweep associated with the evolution of sex-ratio drive in Drosophila simulans

In several Drosophila species, the XY Mendelian ratio is disturbed by X-linked segregation distorters (sex-ratio drive). We used a collection of recombinants between a nondistorting chromosome and a distorting X chromosome originating from the Seychelles to map a candidate sex-ratio region in Drosophila simulans using molecular biallelic markers. Our data were compatible with the presence of a sex-ratio locus in the 7F cytological region. Using sequence polymorphism at the Nrg locus, we showed that sex-ratio has induced a strong selective sweep in populations from Madagascar and Réunion, where distorting chromosomes are close to a 50% frequency. The complete association between the marker and the sex-ratio phenotype and the near absence of mutations and recombination in the studied fragment after the sweep event indicate that this event is recent. Examples of selective sweeps are increasingly reported in a number of genomes. This case identifies the causal selective force. It illustrates that all selective sweeps are not necessarily indicative of an increase in the average fitness of populations.     

Impaired integrin‐mediated adhesion contributes to reduced barrier properties in VASP‐deficient microvascular endothelium

Recent studies point to a significant role of vasodilator‐stimulated phosphoprotein (VASP) in the maintenance of endothelial barrier functions in vivo and in vitro. Moreover, it has been reported that VASP is required for activation of the small GTPase Rac 1. However, little is known whether VASP is involved in the regulation of cell adhesion molecules that are critical for maintenance of the endothelial barrier. Here we demonstrate that impaired barrier properties in VASP‐deficient (VASP?/?) microvascular myocardial endothelial cells (MyEnd) correlated with both impaired integrin‐mediated adhesion as revealed by laser tweezer trapping and reduced integrin‐dependent cell migration. This was paralleled by reduction of focal adhesions at the cell periphery as well as of β₁‐integrin and VE‐cadherin cytoskeletal anchorage. Incubation of MyEnd VASP wt with RGD peptide to block interaction of integrins with extracellular matrix (ECM) reduced barrier properties and Rac 1 activity in wt endothelial monolayers mimicking the situation in VASP (?/?) cells under resting conditions. Moreover, cAMP‐mediated Rac 1 activation was reduced under conditions of impaired integrin‐mediated adhesion in wt cells and cAMP‐induced increase in VE‐cadherin cytoskeletal anchorage was abolished in VASP (?/?) endothelium. In summary, these data indicate that VASP is required for integrin‐mediated adhesion which stabilizes endothelial barrier properties at least in part by facilitating Rac 1 activation. J. Cell. Physiol. 220: 357–366, 2009. © 2009 Wiley‐Liss, Inc.     

Maximizing the Power of Principal-Component Analysis of Correlated Phenotypes in Genome-wide Association Studies

Many human traits are highly correlated. This correlation can be leveraged to improve the power of genetic association tests to identify markers associated with one or more of the traits. Principal component analysis (PCA) is a useful tool that has been widely used for the multivariate analysis of correlated variables. PCA is usually applied as a dimension reduction method: the few top principal components (PCs) explaining most of total trait variance are tested for association with a predictor of interest, and the remaining components are not analyzed. In this study we review the theoretical basis of PCA and describe the behavior of PCA when testing for association between a SNP and correlated traits. We then use simulation to compare the power of various PCA-based strategies when analyzing up to 100 correlated traits. We show that contrary to widespread practice, testing only the top PCs often has low power, whereas combining signal across all PCs can have greater power. This power gain is primarily due to increased power to detect genetic variants with opposite effects on positively correlated traits and variants that are exclusively associated with a single trait. Relative to other methods, the combined-PC approach has close to optimal power in all scenarios considered while offering more flexibility and more robustness to potential confounders. Finally, we apply the proposed PCA strategy to the genome-wide association study of five correlated coagulation traits where we identify two candidate SNPs that were not found by the standard approach.