What do we gain from simplicity versus complexity in species distribution models?

Species distribution models (SDMs) are widely used to explain and predict species ranges and environmental niches. They are most commonly constructed by inferring species' occurrence-environment relationships using statistical and machine-learning methods. The variety of methods that can be used to construct SDMs (e.g. generalized linear/additive models, tree-based models, maximum entropy, etc.), and the variety of ways that such models can be implemented, permits substantial flexibility in SDM complexity. Building models with an appropriate amount of complexity for the study objectives is critical for robust inference. We characterize complexity as the shape of the inferred occurrence-environment relationships and the number of parameters used to describe them, and search for insights into whether additional complexity is informative or superfluous. By building 'under fit' models, having insufficient flexibility to describe observed occurrence-environment relationships, we risk misunderstanding the factors shaping species distributions. By building 'over fit' models, with excessive flexibility, we risk inadvertently ascribing pattern to noise or building opaque models. However, model selection can be challenging, especially when comparing models constructed under different modeling approaches. Here we argue for a more pragmatic approach: researchers should constrain the complexity of their models based on study objective, attributes of the data, and an understanding of how these interact with the underlying biological processes. We discuss guidelines for balancing under fitting with over fitting and consequently how complexity affects decisions made during model building. Although some generalities are possible, our discussion reflects differences in opinions that favor simpler versus more complex models. We conclude that combining insights from both simple and complex SDM building approaches best advances our knowledge of current and future species ranges.

Phenotypic consequences of copy number variation: insights from Smith-Magenis and Potocki-Lupski syndrome mouse models.

A large fraction of genome variation between individuals is comprised of submicroscopic copy number variation of genomic DNA segments. We assessed the relative contribution of structural changes and gene dosage alterations on phenotypic outcomes with mouse models of Smith-Magenis and Potocki-Lupski syndromes. We phenotyped mice with 1n (Deletion/+), 2n (+/+), 3n (Duplication/+), and balanced 2n compound heterozygous (Deletion/Duplication) copies of the same region. Parallel to the observations made in humans, such variation in gene copy number was sufficient to generate phenotypic consequences: in a number of cases diametrically opposing phenotypes were associated with gain versus loss of gene content. Surprisingly, some neurobehavioral traits were not rescued by restoration of the normal gene copy number. Transcriptome profiling showed that a highly significant propensity of transcriptional changes map to the engineered interval in the five assessed tissues. A statistically significant overrepresentation of the genes mapping to the entire length of the engineered chromosome was also found in the top-ranked differentially expressed genes in the mice containing rearranged chromosomes, regardless of the nature of the rearrangement, an observation robust across different cell lineages of the central nervous system. Our data indicate that a structural change at a given position of the human genome may affect not only locus and adjacent gene expression but also "genome regulation." Furthermore, structural change can cause the same perturbation in particular pathways regardless of gene dosage. Thus, the presence of a genomic structural change, as well as gene dosage imbalance, contributes to the ultimate phenotype.

Evaluation of cumulative PCB exposure estimated by a job exposure matrix versus PCB serum concentrations.

Although polychlorinated biphenyls (PCBs) have been banned in many countries for more than three decades, exposures to PCBs continue to be of concern due to their long half-lives and carcinogenic effects. In National Institute for Occupational Safety and Health studies, we are using semiquantitative plant-specific job exposure matrices (JEMs) to estimate historical PCB exposures for workers (n = 24,865) exposed to PCBs from 1938 to 1978 at three capacitor manufacturing plants. A subcohort of these workers (n = 410) employed in two of these plants had serum PCB concentrations measured at up to four times between 1976 and 1989. Our objectives were to evaluate the strength of association between an individual worker's measured serum PCB levels and the same worker's cumulative exposure estimated through 1977 with the (1) JEM and (2) duration of employment, and to calculate the explained variance the JEM provides for serum PCB levels using (3) simple linear regression. Consistent strong and statistically significant associations were observed between the cumulative exposures estimated with the JEM and serum PCB concentrations for all years. The strength of association between duration of employment and serum PCBs was good for highly chlorinated (Aroclor 1254/HPCB) but not less chlorinated (Aroclor 1242/LPCB) PCBs. In the simple regression models, cumulative occupational exposure estimated using the JEMs explained 14-24% of the variance of the Aroclor 1242/LPCB and 22-39% for Aroclor 1254/HPCB serum concentrations. We regard the cumulative exposure estimated with the JEM as a better estimate of PCB body burdens than serum concentrations quantified as Aroclor 1242/LPCB and Aroclor 1254/HPCB.

Which subgroup of patients with rheumatoid arthritis benefits from switching to rituximab versus alternative anti-tumour necrosis factor (TNF) agents after previous failure of an anti-TNF agent?

BACKGROUND: Patients with rheumatoid arthritis (RA) with an inadequate response to TNF antagonists (aTNFs) may switch to an alternative aTNF or start treatment from a different class of drugs, such as rituximab (RTX). It remains unclear in which clinical settings these therapeutic strategies offer most benefit. OBJECTIVE: To analyse the effectiveness of RTX versus alternative aTNFs on RA disease activity in different subgroups of patients. METHODS: A prospective cohort study of patients with RA who discontinued at least one aTNF and subsequently received either RTX or an alternative aTNF, nested within the Swiss RA registry (SCQM-RA) was carried out. The primary outcome, longitudinal improvement in 28-joint count Disease Activity Score (DAS28), was analysed using multivariate regression models for longitudinal data and adjusted for potential confounders. RESULTS: Of the 318 patients with RA included; 155 received RTX and 163 received an alternative aTNF. The relative benefit of RTX varied with the type of prior aTNF failure: when the motive for switching was ineffectiveness to previous aTNFs, the longitudinal improvement in DAS28 was significantly better with RTX than with an alternative aTNF (p = 0.03; at 6 months, -1.34 (95% CI -1.54 to -1.15) vs -0.93 (95% CI -1.28 to -0.59), respectively). When the motive for switching was other causes, the longitudinal improvement in DAS28 was similar for RTX and alternative aTNFs (p = 0.40). These results were not significantly modified by the number of previous aTNF failures, the type of aTNF switches, or the presence of co-treatment with a disease-modifying antirheumatic drug. CONCLUSION: This observational study suggests that in patients with RA who have stopped a previous aTNF treatment because of ineffectiveness changing to RTX is more effective than switching to an alternative aTNF.

Synchronous versus asynchronous modeling of gene regulatory networks.

MOTIVATION: In silico modeling of gene regulatory networks has gained some momentum recently due to increased interest in analyzing the dynamics of biological systems. This has been further facilitated by the increasing availability of experimental data on gene-gene, protein-protein and gene-protein interactions. The two dynamical properties that are often experimentally testable are perturbations and stable steady states. Although a lot of work has been done on the identification of steady states, not much work has been reported on in silico modeling of cellular differentiation processes. RESULTS: In this manuscript, we provide algorithms based on reduced ordered binary decision diagrams (ROBDDs) for Boolean modeling of gene regulatory networks. Algorithms for synchronous and asynchronous transition models have been proposed and their corresponding computational properties have been analyzed. These algorithms allow users to compute cyclic attractors of large networks that are currently not feasible using existing software. Hereby we provide a framework to analyze the effect of multiple gene perturbation protocols, and their effect on cell differentiation processes. These algorithms were validated on the T-helper model showing the correct steady state identification and Th1-Th2 cellular differentiation process. AVAILABILITY: The software binaries for Windows and Linux platforms can be downloaded from http://si2.epfl.ch/~garg/genysis.html.

Evaluation of tumor vascularisation in two rat sarcoma models for studying isolated lung perfusion : Injection route determines the origin of tumour vessels

Résumé Introduction: La perfusion isolée cytostatique du poumon est une technique attractive qui permet l'administration des doses élevées d'un agent cytostatique tout en épargnant dans la mesure du possible la circulation systémique. Cependant, la perfusion de l'artère pulmonaire risque d'épargner le territoire pulmonaire vascularisé par l'intermédiaire des artères bronchiques, ce qui pourrait diminuer l'efficacité de ce traitement au cas où la lésion ciblée est vascularisée par les artères bronchiques. Ce travail est destiné au développement d'un modèle tumoral au niveau des poumons de rongeur (rat) porteur d'un sarcome pulmonaire afin de déterminer si la voie d'injection des cellules tumorales (intraveineuse, versus intratumorale) influencera la vascularisation des tumeurs (provenant du système artères pulmonaires ou artères bronchiques). Méthod: Des tumeurs de sarcomes pulmonaires ont été générées par injection d'une suspension cellulaire de sarcome, soit par injection intraveineuse, soit directement dans le parenchyme pulmonaire par thoracotomie. Ensuite, une perfusion isolée du poumon porteur de la tumeur à l'aide de l'encre a été effectuée, soit par l'artère pulmonaire, soit par le système des artères bronchiques. La distribution de l'encre dans les vaisseaux tumoraux ainsi que dans les vaisseaux non tumoraux du poumon adjacent a été investiguée à l'aide d'une analyse histologique des poumons perfusés. Résultat: L'administration intraveineuse et intratumorale de la suspension de cellules tumorales résulte en des tumeurs similaires sur le plan histologique. Néanmoins, l'injection intra-parenchymateuse démontre des tumeurs plus homogènes et avec un développement plus prédictible, était associée à une survie plus longue qu'après injection intraveineuse. Les analyses histologiques après perfusion isolée à l'aide de l'encre démontre que les tumeurs résultant de l'injection intraveineuse ont développé une vascularisation se basant sur le système d'artères pulmonaires tandis que les tumeurs émergeant après injection intraparenchymateuse ont développé une vascularisation provenant du système des artères bronchiques. Conclusion: Ce travail démontre pour la première fois l'importance du mode de génération de tumeurs pulmonaires en ce qui concerne leur future vascularisation, ce qui pourrait avoir un impact sur leur traitement par perfusion isolée du poumon. Abstract Isolated cytostatic lung perfusion (ILP) is an attractive technique allowing delivery of a high-dose of cytostatic agents to the lungs while limiting systemic toxicity. In developing a rat model of ILP, we have analysed the effect of the route of tumour cell injection on the source of tumour vessels. Pulmonary sarcomas were estab¬lished by injecting a sarcoma cell suspension either by the intravenous (i.v.) route or directly into the lung paren¬chyma. Ink perfusion through either pulmonary artery (PA) or bronchial arteries (BA) was performed and the characteristics of the tumour deposits defined. i.v. and direct injection methods induced pulmonary sarcoma nodules, with similar histological features. The intraparenchymal injection of tumour cells resulted in more reli¬able and reproducible tumour growth and was associat¬ed with a longer survival of the animals. i.v. injected tumours developed a PA-derived vascular tree whereas directly injected tumours developed a BA-derived vasculature.

Weekend versus weekday admission and mortality after acute pulmonary embolism.

BACKGROUND: Optimal management of acute pulmonary embolism (PE) requires medical expertise, diagnostic testing, and therapies that may not be available consistently throughout the entire week. We sought to assess whether associations exist between weekday or weekend admission and mortality and length of hospital stay for patients hospitalized with PE. METHODS AND RESULTS: We evaluated patients discharged with a primary diagnosis of PE from 186 acute care hospitals in Pennsylvania (January 2000 to November 2002). We used random-effect logistic models to study the association between weekend admission and 30-day mortality and used discrete survival models to study the association between weekend admission and time to hospital discharge, adjusting for hospital (region, size, and teaching status) and patient factors (race, insurance, severity of illness, and use of thrombolytic therapy). Among 15 531 patient discharges with PE, 3286 patients (21.2%) had been admitted on a weekend. Patients admitted on weekends had a higher unadjusted 30-day mortality rate (11.1% versus 8.8%) than patients admitted on weekdays, with no difference in length of stay. Patients admitted on weekends had significantly greater adjusted odds of dying (odds ratio 1.17, 95% confidence interval 1.03 to 1.34) than patients admitted on weekdays. The higher mortality among patients hospitalized on weekends was driven by the increased mortality rate among the most severely ill patients. CONCLUSIONS: Patients with PE who are admitted on weekends have a significantly higher short-term mortality than patients admitted on weekdays. Quality-improvement efforts should aim to ensure a consistent approach to the management of PE 7 days a week.

Predictors of the risk of fibrosis at 10 years after breast conserving therapy for early breast cancer: a study based on the EORTC Trial 22881-10882 'boost versus no boost'.

The EORTC 22881-10882 trial in 5178 conservatively treated early breast cancer patients showed that a 16 Gy boost dose significantly improved local control, but increased the risk of breast fibrosis. To investigate predictors for the long-term risk of fibrosis, Cox regression models of the time to moderate or severe fibrosis were developed on a random set of 1797 patients with and 1827 patients without a boost, and validated in the remaining set. The median follow-up was 10.7 years. The risk of fibrosis significantly increased (P<0.01) with increasing maximum whole breast irradiation (WBI) dose and with concomitant chemotherapy, but was independent of age. In the boost arm, the risk further increased (P<0.01) if patients had post-operative breast oedema or haematoma, but it decreased (P<0.01) if WBI was given with >6 MV photons. The c-index was around 0.62. Nomograms with these factors are proposed to forecast the long-term risk of moderate or severe fibrosis.

An evaluation of new parsimony-based versus parametric inference methods in biogeography: a case study using the globally distributed plant family Sapindaceae

Aim  Recently developed parametric methods in historical biogeography allow researchers to integrate temporal and palaeogeographical information into the reconstruction of biogeographical scenarios, thus overcoming a known bias of parsimony-based approaches. Here, we compare a parametric method, dispersal-extinction-cladogenesis (DEC), against a parsimony-based method, dispersal-vicariance analysis (DIVA), which does not incorporate branch lengths but accounts for phylogenetic uncertainty through a Bayesian empirical approach (Bayes-DIVA). We analyse the benefits and limitations of each method using the cosmopolitan plant family Sapindaceae as a case study.Location  World-wide.Methods  Phylogenetic relationships were estimated by Bayesian inference on a large dataset representing generic diversity within Sapindaceae. Lineage divergence times were estimated by penalized likelihood over a sample of trees from the posterior distribution of the phylogeny to account for dating uncertainty in biogeographical reconstructions. We compared biogeographical scenarios between Bayes-DIVA and two different DEC models: one with no geological constraints and another that employed a stratified palaeogeographical model in which dispersal rates were scaled according to area connectivity across four time slices, reflecting the changing continental configuration over the last 110 million years.Results  Despite differences in the underlying biogeographical model, Bayes-DIVA and DEC inferred similar biogeographical scenarios. The main differences were: (1) in the timing of dispersal events - which in Bayes-DIVA sometimes conflicts with palaeogeographical information, and (2) in the lower frequency of terminal dispersal events inferred by DEC. Uncertainty in divergence time estimations influenced both the inference of ancestral ranges and the decisiveness with which an area can be assigned to a node.Main conclusions  By considering lineage divergence times, the DEC method gives more accurate reconstructions that are in agreement with palaeogeographical evidence. In contrast, Bayes-DIVA showed the highest decisiveness in unequivocally reconstructing ancestral ranges, probably reflecting its ability to integrate phylogenetic uncertainty. Care should be taken in defining the palaeogeographical model in DEC because of the possibility of overestimating the frequency of extinction events, or of inferring ancestral ranges that are outside the extant species ranges, owing to dispersal constraints enforced by the model. The wide-spanning spatial and temporal model proposed here could prove useful for testing large-scale biogeographical patterns in plants.

Cost-effectiveness of opportunistic versus organised mammography screening in Switzerland

BACKGROUND: Various centralised mammography screening programmes have shown to reduce breast cancer mortality at reasonable costs. However, mammography screening is not necessarily cost-effective in every situation. Opportunistic screening, the predominant screening modality in several European countries, may under certain circumstances be a cost-effective alternative. In this study, we compared the cost-effectiveness of both screening modalities in Switzerland. METHODS: Using micro-simulation modelling, we predicted the effects and costs of biennial mammography screening for 50-69 years old women between 1999 and 2020, in the Swiss female population aged 30-70 in 1999. A sensitivity analysis on the test sensitivity of opportunistic screening was performed. RESULTS: Organised mammography screening with an 80% participation rate yielded a breast cancer mortality reduction of 13%. Twenty years after the start of screening, the predicted annual breast cancer mortality was 25% lower than in a situation without screening. The 3% discounted cost-effectiveness ratio of organised mammography screening was euro11,512 per life year gained. Opportunistic screening with a similar participation rate was comparably effective, but at twice the costs: euro22,671-24,707 per life year gained. This was mainly related to the high costs of opportunistic mammography and frequent use of imaging diagnostics in combination with an opportunistic mammogram. CONCLUSION: Although data on the performance of opportunistic screening are limited, both opportunistic and organised mammography screening seem effective in reducing breast cancer mortality in Switzerland. However, for opportunistic screening to become equally cost-effective as organised screening, costs and use of additional diagnostics should be reduced.

Evolution in heterogeneous populations: from migration models to fixation probabilities.

Although dispersal is recognized as a key issue in several fields of population biology (such as behavioral ecology, population genetics, metapopulation dynamics or evolutionary modeling), these disciplines focus on different aspects of the concept and often make different implicit assumptions regarding migration models. Using simulations, we investigate how such assumptions translate into effective gene flow and fixation probability of selected alleles. Assumptions regarding migration type (e.g. source-sink, resident pre-emption, or balanced dispersal) and patterns (e.g. stepping-stone versus island dispersal) have large impacts when demes differ in sizes or selective pressures. The effects of fragmentation, as well as the spatial localization of newly arising mutations, also strongly depend on migration type and patterns. Migration rate also matters: depending on the migration type, fixation probabilities at an intermediate migration rate may lie outside the range defined by the low- and high-migration limits when demes differ in sizes. Given the extreme sensitivity of fixation probability to characteristics of dispersal, we underline the importance of making explicit (and documenting empirically) the crucial ecological/ behavioral assumptions underlying migration models.

Incidents occuring during anesthesia for total hip arthroplasty: a comparison of general versus regional anaesthesia

Background: The type of anesthesia to be used for total hip arthroplasty (THA) is still a matter of debate. We compared the occurrence of per- and post-anesthesia incidents in patients receiving either general (GA) or regional anesthesia (RA). Methods: We used data from 29 hospitals, routinely collected in the Anaesthesia Databank Switzerland register between January 2001 and December 2003. We used multi-level logistic regression models. Results: There were more per- and post-anesthesia incidents under GA compared to RA (35.1% vs 32.7 %, n = 3191, and 23.1% vs 19.4%, n = 3258, respectively). In multi-level logistic regression analysis, RA was significantly associated with a lower incidence of per-anesthetic problems, especially hypertension, compared with GA. During the post-anesthetic period, RA was also less associated with pain. Conversely, RA was more associated with post-anesthetic hypotension, especially for epidural technique. In addition, age and ASA were more associated with incidents under GA compared to RA. Men were more associated with per-anesthetic problems under RA compared to GA. Whereas increased age (>67), gender (male), and ASA were linked with the choice of RA, we noticed that this choice depended also on hospital practices after we adjusted for the other variables. Conclusions: Compared to RA, GA was associated with an increased proportion of per- and post-anesthesia incidents. Although this study is only observational, it is rooted in daily practice. Whereas RA might be routinely proposed, GA might be indicated because of contraindications to RA, patients' preferences or other surgical or anaesthesiology related reasons. Finally, the choice of a type of anesthesia seems to depend on local practices that may differ between hospitals.

Incidents occurring during anesthesia for total hip arthroplasty: a comparison of general versus regional anaesthesia

Background: The type of anesthesia to be used for total hip arthroplasty (THA) is still a matter of debate. We compared the occurrence of per- and post-anesthesia incidents in patients receiving either general (GA) or regional anesthesia (RA). Methods: We used data from 29 hospitals, routinely collected in the Anaesthesia Databank Switzerland register between January 2001 and December 2003. We used multi-level logistic regression models. Results: There were more per- and post-anesthesia incidents under GA compared to RA (35.1% vs 32.7 %, n = 3191, and 23.1% vs 19.4%, n = 3258, respectively). In multi-level logistic regression analysis, RA was significantly associated with a lower incidence of per-anesthetic problems, especially hypertension, compared with GA. During the post-anesthetic period, RA was also less associated with pain. Conversely, RA was more associated with post-anesthetic hypotension, especially for epidural technique. In addition, age and ASA were more associated with incidents under GA compared to RA. Men were more associated with per-anesthetic problems under RA compared to GA. Whereas increased age (>67), gender (male), and ASA were linked with the choice of RA, we noticed that this choice depended also on hospital practices after we adjusted for the other variables. Conclusions: Compared to RA, GA was associated with an increased proportion of per- and post-anesthesia incidents. Although this study is only observational, it is rooted in daily practice. Whereas RA might be routinely proposed, GA might be indicated because of contraindications to RA, patients' preferences or other surgical or anaesthesiology related reasons. Finally, the choice of a type of anesthesia seems to depend on local practices that may differ between hospitals.

Genetic differentiation for size at first reproduction through male versus female functions in the widespread Mediterranean tree Pinus pinaster.

BACKGROUND AND AIMS: The study of local adaptation in plant reproductive traits has received substantial attention in short-lived species, but studies conducted on forest trees are scarce. This lack of research on long-lived species represents an important gap in our knowledge, because inferences about selection on the reproduction and life history of short-lived species cannot necessarily be extrapolated to trees. This study considers whether the size for first reproduction is locally adapted across a broad geographical range of the Mediterranean conifer species Pinus pinaster. In particular, the study investigates whether this monoecious species varies genetically among populations in terms of whether individuals start to reproduce through their male function, their female function or both sexual functions simultaneously. Whether differences among populations could be attributed to local adaptation across a climatic gradient is then considered. METHODS: Male and female reproduction and growth were measured during early stages of sexual maturity of a P. pinaster common garden comprising 23 populations sampled across the species range. Generalized linear mixed models were used to assess genetic variability of early reproductive life-history traits. Environmental correlations with reproductive life-history traits were tested after controlling for neutral genetic structure provided by 12 nuclear simple sequence repeat markers. KEY RESULTS: Trees tended to reproduce first through their male function, at a size (height) that varied little among source populations. The transition to female reproduction was slower, showed higher levels of variability and was negatively correlated with vegetative growth traits. Several female reproductive traits were correlated with a gradient of growth conditions, even after accounting for neutral genetic structure, with populations from more unfavourable sites tending to commence female reproduction at a lower individual size. CONCLUSIONS: The study represents the first report of genetic variability among populations for differences in the threshold size for first reproduction between male and female sexual functions in a tree species. The relatively uniform size at which individuals begin reproducing through their male function probably represents the fact that pollen dispersal is also relatively invariant among sites. However, the genetic variability in the timing of female reproduction probably reflects environment-dependent costs of cone production. The results also suggest that early sex allocation in this species might evolve under constraints that do not apply to other conifers.

Home versus in-hospital treatment of outpatients with acute deep venous thrombosis of the lower limbs.

BACKGROUND: Some physicians are still concerned about the safety of treatment at home of patients with acute deep venous thrombosis (DVT). METHODS: We used data from the RIETE (Registro Informatizado de la Enfermedad TromboEmbólica) registry to compare the outcomes in consecutive outpatients with acute lower limb DVT according to initial treatment at home or in the hospital. A propensity score-matching analysis was carried out with a logistic regression model. RESULTS: As of December 2012, 13,493 patients had been enrolled. Of these, 4456 (31%) were treated at home. Patients treated at home were more likely to be male and younger and to weigh more; they were less likely than those treated in the hospital to have chronic heart failure, lung disease, renal insufficiency, anemia, recent bleeding, immobilization, or cancer. During the first week of anticoagulation, 27 patients (0.20%) suffered pulmonary embolism (PE), 12 (0.09%) recurrent DVT, and 51 (0.38%) major bleeding; 80 (0.59%) died. When only patients treated at home were considered, 12 (0.27%) had PE, 4 (0.09%) had recurrent DVT, 6 (0.13%) bled, and 4 (0.09%) died (no fatal PE, 3 fatal bleeds). After propensity analysis, patients treated at home had a similar rate of venous thromboembolism recurrences and a lower rate of major bleeding (odds ratio, 0.4; 95% confidence interval, 0.1-1.0) or death (odds ratio, 0.2; 95% confidence interval, 0.1-0.7) within the first week compared with those treated in the hospital. CONCLUSIONS: In outpatients with DVT, home treatment was associated with a better outcome than treatment in the hospital. These data may help safely treat more DVT patients at home.