Fluctuating selection and its (elusive) evolutionary consequences in a wild rodent population

Temporal fluctuations in the strength and direction of selection are often proposed as a mechanism that slows down evolution, both over geological and contemporary timescales. Both the prevalence of fluctuating selection and its relevance for evolutionary dynamics remain poorly understood however, especially on contemporary timescales: unbiased empirical estimates of variation in selection are scarce, and the question of how much of the variation in selection translates into variation in genetic change has largely been ignored. Using long‐term individual‐based data for a wild rodent population, we quantify the magnitude of fluctuating selection on body size. Subsequently, we estimate the evolutionary dynamics of size and test for a link between fluctuating selection and evolution. We show that, over the past 11 years, phenotypic selection on body size has fluctuated significantly. However, the strength and direction of genetic change have remained largely constant over the study period; that is, the rate of genetic change was similar in years where selection favoured heavier vs. lighter individuals. This result suggests that over shorter timescales, fluctuating selection does not necessarily translate into fluctuating evolution. Importantly however, individual‐based simulations show that the correlation between fluctuating selection and fluctuating evolution can be obscured by the effect of drift, and that substantially more data are required for a precise and accurate estimate of this correlation. We identify new challenges in measuring the coupling between selection and evolution, and provide methods and guidelines to overcome them.     

Among‐year variation in selection during early life stages and the genetic basis of fitness in Arabidopsis thaliana

Incomplete information regarding both selection regimes and the genetic basis of fitness limits our understanding of adaptive evolution. Among‐year variation in the genetic basis of fitness is rarely quantified, and estimates of selection are typically based on single components of fitness, thus potentially missing conflicting selection acting during other life‐history stages. Here, we examined among‐year variation in selection on a key life‐history trait and the genetic basis of fitness covering the whole life cycle in the annual plant Arabidopsis thaliana. We planted freshly matured seeds of >200 recombinant inbred lines (RILs) derived from a cross between two locally adapted populations (Italy and Sweden), and both parental genotypes at the native site of the Swedish population in three consecutive years. We quantified selection against the nonlocal Italian genotype, mapped quantitative trait loci (QTL) for fitness and its components, and quantified selection on timing of germination during different life stages. In all 3 years, the local Swedish genotype outperformed the nonlocal Italian genotype. However, both the contribution of early life stages to relative fitness, and the effects of fitness QTL varied among years. Timing of germination was under conflicting selection through seedling establishment vs. adult survival and fecundity, and both the direction and magnitude of net selection varied among years. Our results demonstrate that selection during early life stages and the genetic basis of fitness can vary markedly among years, emphasizing the need for multiyear studies considering the whole life cycle for a full understanding of natural selection and mechanisms maintaining local adaptation.     

Comparative Cost-Effectiveness of Drugs in Early versus Late Stages of Cancer; Review of the Literature and a Case Study in Breast Cancer

Background

Many oncological drugs that are being used in the adjuvant setting were first submitted for reimbursement in the metastatic stage, with differences in incremental cost-effectiveness ratios (ICERs) in both settings having potential implications for reimbursement and pricing. The aim of this study is to identify a possible trend in the cost-effectiveness for the early/adjuvant and late/metastatic stages of oncological drugs through review and case study.

Methods

We reviewed pairs of cost-effectiveness analyses of the same oncological drug in different stages for Scotland and the Netherlands. The case study in this report was directed at trastuzumab in the Dutch situation. Using a simplified Markov model, the cost-effectiveness in early and late stage of breast cancer was calculated and compared to the findings from the review.

Results

Comparable studies were found for cetuximab, bortezomib and bosutinib. Treatments in the late stage were found to be more expensive per QALY by a factor ranging from 1.5 to 12. The case study provided a similar result; late stage treatment was more expensive by a factor 10. Using, for example, a threshold of €80,000/QALY, the early stage of cetuximab, bosutinib and trastuzumab are deemed cost-effective, while their compared late stage is lifted over the threshold and potentially considered not cost-effective.

Conclusion

ICERs of oncological drugs used in different stages are more unfavourable in the late stage than in the early stage. Applying a reasonable threshold may result in early stage treatment being deemed cost-effective while late stage potentially not. Authorities should be aware of this when assessing oncological drugs and interpreting the corresponding ICERs, in the situation where oncological drugs are generally most submitted for reimbursement in the late stage initially.     

Genetic and environmental sources of covariance among internal reproductive traits in the yellow dung fly

Substantial inter- and intraspecific variation is found in reproductive traits, but the evolutionary implications of this variation remain unclear. One hypothesis is that natural selection favours female reproductive morphology that allows females to control mating and fertilization and that diverse male reproductive traits arise as counter adaptations to subvert this control. Such co-evolution predicts the establishment of genetic correlations between male and female reproductive traits that closely interact during mating. Therefore, we measured phenotypic and genetic correlations between male and female reproductive tract characteristics in the yellow dung fly, Scathophaga stercoraria (Diptera: Scathophagidae), using a nested half-sib breeding experiment. We found significant heritabilities for the size of most reproductive tract traits investigated in both females (spermathecae and their ducts, accessory glands and their ducts) and males (testis size but not sperm length). Within the sexes, phenotypic and genetic correlations were mostly nil or positive, suggesting functional integration of or condition-dependent investment in internal reproductive traits. Negative intrasexual genetic correlations, potentially suggestive of resource allocation trade-offs, were not evident. Intersexual genetic correlations were mostly positive, reflecting expected allometries between male and female morphologies. Most interestingly, testis size correlated positively with female accessory gland size and duct length, potentially indicative of a co-evolutionary arms race. We discuss these and alternative explanations for these patterns of genetic covariance.     

Multidrug resistance-associated protein-1 (MRP1) genetic variants,MRP1 protein levels and severity of COPD

Background

Multidrug resistance-associated protein-1 (MRP1) protects against oxidative stress and toxic compounds generated by cigarette smoking, which is the main risk factor for chronic obstructive pulmonary disease (COPD). We have previously shown that single nucleotide polymorphisms (SNPs) in MRP1 significantly associate with level of FEV₁ in two independent population based cohorts. The aim of our study was to assess the associations of MRP1 SNPs with FEV₁ level, MRP1 protein levels and inflammatory markers in bronchial biopsies and sputum of COPD patients.

Methods

Five SNPs (rs212093, rs4148382, rs504348, rs4781699, rs35621) in MRP1 were genotyped in 110 COPD patients. The effects of MRP1 SNPs were analyzed using linear regression models.

Results

One SNP, rs212093 was significantly associated with a higher FEV₁ level and less airway wall inflammation. Another SNP, rs4148382 was significantly associated with a lower FEV₁ level, higher number of inflammatory cells in induced sputum and with a higher MRP1 protein level in bronchial biopsies.

Conclusions

This is the first study linking MRP1 SNPs with lung function and inflammatory markers in COPD patients, suggesting a role of MRP1 SNPs in the severity of COPD in addition to their association with MRP1 protein level in bronchial biopsies.     

Competing interactions in DNA assembly on graphene

We study the patterns that short strands of single-stranded DNA form on the top graphene surface of graphite. We find that the DNA assembles into two distinct patterns, small spherical particles and elongated networks. Known interaction models based on DNA-graphene binding, hydrophobic interactions, or models based on the purine/pyrimidine nature of the bases do not explain our observed crossover in pattern formation. We argue that the observed assembly behavior is caused by a crossover in the competition between base-base pi stacking and base-graphene pi stacking and we infer a critical crossover energy of 0.3-0.5 eV. The experiments therefore provide a projective measurement of the base-base interaction strength.     

Smad gene expression in pulmonary fibroblasts: indications for defective ECM repair in COPD

Background

Endothelial dysfunction is a major complication of pulmonary endarterectomy (PTE) that can lead to pulmonary edema and persistent pulmonary hypertension. We hypothesized that endothelial dysfunction is related to increased endothelial-cell (EC) death.

Methods

In piglets, the left pulmonary artery (PA) was ligated to induce lung ischemia then reimplanted into the main PA to reperfuse the lung. Animals sacrificed 5 weeks after ligation (n = 5), 2 days after reperfusion (n = 5), or 5 weeks after reperfusion (n = 5) were compared to a sham-operated group (n = 5). PA vasoreactivity was studied and eNOS assayed. EC apoptosis was assessed by TUNEL in the proximal and distal PA and by caspase-3 activity assay in the proximal PA. Gene expression of pro-apoptotic factors (thrombospondin-1 (Thsp-1) and plasminogen activator inhibitor 1 (PAI-1)) and anti-apoptotic factors vascular endothelial growth factor (VEGF) and basic fibroblast growth factor (bFGF) was investigated by QRT-PCR.

Results

Endothelium-dependent relaxation was altered 5 weeks after ligation (p = 0.04). The alterations were exacerbated 2 days after reperfusion (p = 0.002) but recovered within 5 weeks after reperfusion. EC apoptosis was increased 5 weeks after PA ligation (p = 0.02), increased further within 2 days after reperfusion (p < 0.0001), and returned to normal within 5 weeks after reperfusion. Whereas VEGF and bFGF expressions remained unchanged, TSP and PAI-1 expressions peaked 5 weeks after ligation (p = 0.001) and returned to normal within 2 days after reperfusion.

Conclusion

Chronic lung ischemia induces over-expression of pro-apoptotic factors. Lung reperfusion is followed by a dramatic transient increase in EC death that may explain the development of endothelial dysfunction after PE. Anti-apoptotic agents may hold considerable potential for preventing postoperative complications.     

Attenuated AMPA Receptor Expression Allows Glioblastoma Cell Survival in Glutamate-Rich Environment

Background

Glioblastoma multiforme (GBM) cells secrete large amounts of glutamate that can trigger AMPA-type glutamate receptors (AMPARs). This commonly results in Na⁺ and Ca²⁺-permeability and thereby in excitotoxic cell death of the surrounding neurons. Here we investigated how the GBM cells themselves survive in a glutamate-rich environment.

Methods and Findings

In silico analysis of published reports shows down-regulation of all ionotropic glutamate receptors in GBM as compared to normal brain. In vitro, in all GBM samples tested, mRNA expression of AMPAR subunit GluR1, 2 and 4 was relatively low compared to adult and fetal total brain mRNA and adult cerebellum mRNA. These findings were in line with primary GBM samples, in which protein expression patterns were down-regulated as compared to the normal tissue. Furthermore, mislocalized expression of these receptors was found. Sequence analysis of GluR2 RNA in primary and established GBM cell lines showed that the GluR2 subunit was found to be partly unedited.

Conclusions

Together with the lack of functional effect of AMPAR inhibition by NBQX our results suggest that down-regulation and afunctionality of AMPARs, enable GBM cells to survive in a high glutamate environment without going into excitotoxic cell death themselves. It can be speculated that specific AMPA receptor inhibitors may protect normal neurons against the high glutamate microenvironment of GBM tumors.     

Down-regulation of E-cadherin in human bronchial epithelial cells leads to epidermal growth factor receptor-dependent Th2 cell-promoting activity

Airway epithelial cells are well-known producers of thymus- and activation-regulated chemokine (TARC), a Th2 cell-attracting chemokine that may play an important role in the development of allergic airway inflammation. However, the mechanism responsible for up-regulation of TARC in allergy is still unknown. In the asthmatic airways, loss of expression of the cell-cell contact molecule E-cadherin and reduced epithelial barrier function has been observed, which may be the result of an inadequate repair response. Because E-cadherin also suppressed multiple signaling pathways, we studied whether disruption of E-cadherin-mediated cell contact may contribute to increased proallergic activity of epithelial cells, e.g., production of the chemokine TARC. We down-regulated E-cadherin in bronchial epithelial cells by small interference RNA and studied effects on electrical resistance, signaling pathways, and TARC expression (by electric cell-substrate impedance sensing, immunodetection, immunofluorescent staining, and real-time PCR). Small interference RNA silencing of E-cadherin resulted in loss of E-cadherin-mediated junctions, enhanced phosphorylation of epidermal growth factor receptor (EGFR), and the downstream targets MEK/ERK-1/2 and p38 MAPK, finally resulting in up-regulation of TARC as well as thymic stromal lymphopoietin expression. The use of specific inhibitors revealed that the effect on TARC is mediated by EGFR-dependent activation of the MAPK pathways. In contrast to TARC, expression of the Th1/Treg cell-attracting chemokine RANTES was unaffected by E-cadherin down-regulation. In summary, we show that loss of E-cadherin-mediated epithelial cell-cell contact by damaging stimuli, e.g., allergens, may result in reduced suppression of EGFR-dependent signaling pathways and subsequent induction of Th2 cell-attracting molecule TARC. Thus, disruption of intercellular epithelial contacts may specifically promote Th2 cell recruitment in allergic asthma.     

Dynamic and static calcium gradients inside large snail (Helix aspersa) neurones detected with calcium-sensitive microelectrodes

We have used quartz Ca2+-sensitive microelectrodes (CASMs) in large voltage-clamped snail neurones to investigate the inward spread of Ca2+ after a brief depolarisation. Both steady state and [Ca2+]i transients changed with depth of penetration. When the CASM tip was within 20 microm of the far side of the cell the [Ca2+]i transient time to peak was 4.4+/-0.5s, rising to 14.7+/-0.7s at a distance of 80 microm. We estimate that the Ca2+ transients travelled centripetally at an average speed of 6 microm2 s(-1) and decreased in size by half over a distance of about 45 microm. Cyclopiazonic acid had little effect on the size and time to peak of Ca2+ transients but slowed their recovery significantly. This suggests that the endoplasmic reticulum curtails rather than reinforces the transients. Injecting the calcium buffer BAPTA made the Ca2+ transients more uniform in size and increased their times to peak and rates of recovery near the membrane. We have developed a computational model for the transients, which includes diffusion, uptake and Ca2+ extrusion. Good fits were obtained with a rather large apparent diffusion coefficient of about 90+/-20 microm2 s(-1).This may assist fast recovery by extrusion.