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.

Participation of interleukin-5, interleukin-8 and leukotriene B4 in eosinophil accumulation in two different experimental models

There are several experimental models describing in vivo eosinophil (EO) migration, including ip injection of a large volume of saline (SAL) or Sephadex beads (SEP). The aim of this study was to investigate the mechanisms involved in the EO migration in these two models. Two consecutive injections of SAL given 48 hr apart, induced a selective recruitment of EO into peritoneal cavity of rats, which peaked 48 hr after the last injection. SEP, when injected ip, promoted EO accumulation in rats. The phenomenom was dose-related and peaked 48 hr after SEP injection. To investigate the mediators involved in this process we showed that BW A4C, MK 886 and dexamethasone (DXA) inhibited the EO migration induced by SAL and SEP. To investigate the source of the EO chemotactic factor we showed that mast cells, macrophages (MO), but not lymphocytes, incubated in vitro in presence of SAL released a factor which induced EO migration. With SEP, only mast cells release a factor that induced EO migration, which was inhibited by BW A4C, MK 886 and DXA. Furthermore, the chemotactic activity of SAL-stimulated mast cells was inhibited by antisera against IL-5 and IL-8 (interleukin). SAL-stimulated MO were only inhibited by anti-IL-8 antibodies as well SEP-stimulated mast cells. These results suggest that the EO migration induced by SAL may be dependent on resident mast cells and MO and mediated by LTB4, IL-5 and IL-8. SEP-induced EO migration was dependent on mast cells and may be mediated by LTB4 and IL-8. Furthermore, IL-5 and IL-8 induced EO migration, which was also dependent on resident cells and mediated by LTB4 . In conclusion, EO migration induced by SAL is dependent on mast cells and MO, whereas that induced by SEP is dependent on mast cells alone. Stimulated mast cells release LTB4, IL-5 and IL-8 while MO release LTB4 and IL-8. The IL-5 and IL-8 release by the SAL or SEP-stimulated resident cells may act in an autocrine fashion, thus potentiating LTB4 release.

Role of eosinophilic airway inflammation in models of asthma

Eosinophils play a central role in the establishment and outcome of bronchial inflammation in asthma. Animal models of allergy are useful to answer questions related to mechanisms of allergic inflammation. We have used models of sensitized and boosted guinea pigs to investigate the nature of bronchial inflammation in allergic conditions. These animals develop marked bronchial infiltration composed mainly of CD4+ T-lymphocytes and eosinophils. Further provocation with antigen leads to degranulation of eosinophils and ulceration of the bronchial mucosa. Eosinophils are the first cells to increase in numbers in the mucosa after antigen challenge and depend on the expression of alpha 4 integrin to adhere to the vascular endothelium and transmigrate to the mucosa. Blockage of alpha4 integrin expression with specific antibody prevents not only the transmigration of eosinophils but also the development of bronchial hyperresponsiveness (BHR) to agonists in sensitized and challenged animals, clearly suggesting a role for this cell type in this altered functional state. Moreover, introduction of antibody against Major Basic Protein into the airways also prevents the development of BHR in similar model. BHR can also be suppressed by the use of FK506, an immunosuppressor that reduces in almost 100% the infiltration of eosinophils into the bronchi of allergic animals. These data support the concept that eosinophil is the most important pro-inflammatory factor in bronchial inflammation associated with allergy.

Local and global consistency properties for student placement

In the context of resource allocation on the basis of priorities, Ergin (2002) identifies a necessary and sufficient condition on the priority structure such that the student-optimal stable mechanism satisfies a consistency principle. Ergin (2002) formulates consistency as a local property based on a fixed population of agents and fixed resources -- we refer to this condition as local consistency and to his condition on the priority structure as local acyclicity. We identify a related but stronger necessary and sufficient condition (unit acyclicity) on the priority structure such that the student-optimal stable mechanism satisfies a more standard global consistency property. Next, we provide necessary and sufficient conditions for the student-optimal stable mechanism to satisfy converse consistency principles. We identify a necessary and sufficient condition (local shift-freeness) on the priority structure such that the student-optimal stable mechanism satisfies local converse consistency. Interestingly, local acyclicity implies local shift-freeness and hence the student-optimal stable mechanism more frequently satisfies local converse consistency than local consistency. Finally, in order for the student-optimal stable mechanism to be globally conversely consistent, one again has to impose unit acyclicity on the priority structure. Hence, unit acyclicity is a necessary and sufficient condition on the priority structure for the student-optimal stable mechanism to satisfy global consistency or global converse consistency.

Mixture of bivariate Poisson regression models with an application to insurance

In a recent paper Bermúdez [2009] used bivariate Poisson regression models for ratemaking in car insurance, and included zero-inflated models to account for the excess of zeros and the overdispersion in the data set. In the present paper, we revisit this model in order to consider alternatives. We propose a 2-finite mixture of bivariate Poisson regression models to demonstrate that the overdispersion in the data requires more structure if it is to be taken into account, and that a simple zero-inflated bivariate Poisson model does not suffice. At the same time, we show that a finite mixture of bivariate Poisson regression models embraces zero-inflated bivariate Poisson regression models as a special case. Additionally, we describe a model in which the mixing proportions are dependent on covariates when modelling the way in which each individual belongs to a separate cluster. Finally, an EM algorithm is provided in order to ensure the models’ ease-of-fit. These models are applied to the same automobile insurance claims data set as used in Bermúdez [2009] and it is shown that the modelling of the data set can be improved considerably.

Electrode location and clinical outcome in hippocampal electrical stimulation for mesial temporal lobe epilepsy.

PURPOSE: To study the clinical outcome in hippocampal deep brain stimulation (DBS) for the treatment of patients with refractory mesial temporal lobe epilepsy (MTLE) according to the electrode location. METHODS: Eight MTLE patients implanted in the hippocampus and stimulated with high-frequency DBS were included in this study. Five underwent invasive recordings with depth electrodes to localize ictal onset zone prior to chronic DBS. Position of the active contacts of the electrode was calculated on postoperative imaging. The distances to the ictal onset zone were measured as well as atlas-based hippocampus structures impacted by stimulation were identified. Both were correlated with seizure frequency reduction. RESULTS: The distances between active electrode location and estimated ictal onset zone were 11±4.3 or 9.1±2.3mm for patients with a >50% or <50% reduction in seizure frequency. In patients (N=6) showing a >50% seizure frequency reduction, 100% had the active contacts located <3mm from the subiculum (p<0.05). The 2 non-responders patients were stimulated on contacts located >3mm to the subiculum. CONCLUSION: Decrease of epileptogenic activity induced by hippocampal DBS in refractory MTLE: (1) seems not directly associated with the vicinity of active electrode to the ictal focus determined by invasive recordings; (2) might be obtained through the neuromodulation of the subiculum.

Climate-based empirical models show biased predictions of butterfly communities along environmental gradients

A better understanding of the factors that mould ecological community structure is required to accurately predict community composition and to anticipate threats to ecosystems due to global changes. We tested how well stacked climate-based species distribution models (S-SDMs) could predict butterfly communities in a mountain region. It has been suggested that climate is the main force driving butterfly distribution and community structure in mountain environments, and that, as a consequence, climate-based S-SDMs should yield unbiased predictions. In contrast to this expectation, at lower altitudes, climate-based S-SDMs overpredicted butterfly species richness at sites with low plant species richness and underpredicted species richness at sites with high plant species richness. According to two indices of composition accuracy, the Sorensen index and a matching coefficient considering both absences and presences, S-SDMs were more accurate in plant-rich grasslands. Butterflies display strong and often specialised trophic interactions with plants. At lower altitudes, where land use is more intense, considering climate alone without accounting for land use influences on grassland plant richness leads to erroneous predictions of butterfly presences and absences. In contrast, at higher altitudes, where climate is the main force filtering communities, there were fewer differences between observed and predicted butterfly richness. At high altitudes, even if stochastic processes decrease the accuracy of predictions of presence, climate-based S-SDMs are able to better filter out butterfly species that are unable to cope with severe climatic conditions, providing more accurate predictions of absences. Our results suggest that predictions should account for plants in disturbed habitats at lower altitudes but that stochastic processes and heterogeneity at high altitudes may limit prediction success of climate-based S-SDMs.

Stochastic time-concentration activity models for cytotoxicity in 3D brain cell cultures.

BACKGROUND: In vitro aggregating brain cell cultures containing all types of brain cells have been shown to be useful for neurotoxicological investigations. The cultures are used for the detection of nervous system-specific effects of compounds by measuring multiple endpoints, including changes in enzyme activities. Concentration-dependent neurotoxicity is determined at several time points. METHODS: A Markov model was set up to describe the dynamics of brain cell populations exposed to potentially neurotoxic compounds. Brain cells were assumed to be either in a healthy or stressed state, with only stressed cells being susceptible to cell death. Cells may have switched between these states or died with concentration-dependent transition rates. Since cell numbers were not directly measurable, intracellular lactate dehydrogenase (LDH) activity was used as a surrogate. Assuming that changes in cell numbers are proportional to changes in intracellular LDH activity, stochastic enzyme activity models were derived. Maximum likelihood and least squares regression techniques were applied for estimation of the transition rates. Likelihood ratio tests were performed to test hypotheses about the transition rates. Simulation studies were used to investigate the performance of the transition rate estimators and to analyze the error rates of the likelihood ratio tests. The stochastic time-concentration activity model was applied to intracellular LDH activity measurements after 7 and 14 days of continuous exposure to propofol. The model describes transitions from healthy to stressed cells and from stressed cells to death. RESULTS: The model predicted that propofol would affect stressed cells more than healthy cells. Increasing propofol concentration from 10 to 100 μM reduced the mean waiting time for transition to the stressed state by 50%, from 14 to 7 days, whereas the mean duration to cellular death reduced more dramatically from 2.7 days to 6.5 hours. CONCLUSION: The proposed stochastic modeling approach can be used to discriminate between different biological hypotheses regarding the effect of a compound on the transition rates. The effects of different compounds on the transition rate estimates can be quantitatively compared. Data can be extrapolated at late measurement time points to investigate whether costs and time-consuming long-term experiments could possibly be eliminated.

Stroke Aphasia: 1,500 Consecutive Cases.

Background: To study the characteristics of vascular aphasia in a cohort of patients with a first-ever stroke. Methods: All patients admitted to the Lausanne neurology department for a first-ever stroke between 1979 and 2004 were included. Neurological examination including language was performed on admission. Stroke risk factors, stroke origin and location, associated symptoms and Rankin scale scores were recorded for each patient. The influence of these factors on aphasia frequency and subtypes was analyzed using logistic regression models. Results: 1,541 (26%) of patients included in this study had aphasia. The more frequent clinical presentations were expressive-receptive aphasia (38%) and mainly expressive aphasia (37%), whereas mainly receptive aphasia was less frequently observed (25%). In ischemic stroke, the frequency of aphasia increased with age (55% of nonaphasic vs. 61% of aphasic patients were more than 65 years old), female sex (40% of women in the nonaphasia group vs. 44% in the aphasia group) and risk factors for cardioembolic origin (coronary heart disease 20 vs. 26% and atrial fibrillation 15 vs. 24%). Stroke aphasia was more likely associated with superficial middle cerebral artery (MCA) stroke and leads to relevant disability. Clinical subtypes depended on stroke location and associated symptoms. Exceptions to the classic clinical-topographic correlations were not rare (26%). Finally, significant differences were found for patients with crossed aphasia in terms of stroke origin and aphasia subtypes. Conclusions: Risk factors for stroke aphasia are age, cardioembolic origin and superficial MCA stroke. Exceptions to classic clinical-topographic correlations are not rare. Stroke aphasia is associated with relevant disability. Stroke location and associated symptoms strongly influence aphasia subtypes.

A non-possibility theorem for joint-stability in interindustry models

Joint-stability in interindustry models relates to the mutual simultaneous consistency of the demand-driven and supply-driven models of Leontief and Ghosh, respectively. Previous work has claimed joint-stability to be an acceptable assumption from the empirical viewpoint, provided only small changes in exogenous variables are considered. We show in this note, however, that the issue has deeper theoretical roots and offer an analytical demonstration that shows the impossibility of consistency between demand-driven and supply-driven models.

Entrevistant infants pre-escolars víctimes d’abús sexual i/o maltractament familiar : Eficàcia dels models d’entrevista forense

Entrevistant infants pre-escolars víctimes d’abús sexual i/o maltractament familiar: eficàcia dels models d’entrevista forense Entrevistar infants en edat preescolar que han viscut una situació traumàtica és una tasca complexa que dins l’avaluació psicològica forense necessita d’un protocol perfectament delimitat, clar i temporalitzat. Per això, s’han seleccionat 3 protocols d’entrevista: el Protocol de Menors (PM) de Bull i Birch, el model del National Institute for Children Development (NICHD) de Michel Lamb, a partir del qual es va desenvolupar l’EASI (Evaluación del Abuso Sexual Infantojuvenil) i l’Entrevista Cognitiva (EC) de Fisher i Geiselman. La hipòtesi de partida vol comprovar si els anteriors models permeten obtenir volums informatius diferents en infants preescolars. Conseqüentment, els objectius han estat determinar quin dels models d’entrevista permet obtenir un volum informatiu amb més precisions i menys errors, dissenyar un model d’entrevista propi i consensuar aquest model. En el treball s’afegeixen esquemes pràctics que facilitin l’obertura, desenvolupament i tancament de l’entrevista forense. La metodologia ha reproduït el binomi infant - esdeveniment traumàtic, mitjançant la visualització i l’explicació d’un fet emocionalment significatiu amb facilitat per identificar-se: l’accident en bicicleta d’un infant que cau, es fa mal, sagna i el seu pare el cura. A partir d’aquí, hem entrevistat 135 infants de P3, P4 i P5, mitjançant els 3 models d’entrevista referits, enfrontant-los a una demanda específica: recordar i narrar aquest esdeveniment. S’ha conclòs que el nivell de record correcte, quan s’utilitza un model d’entrevista adequat amb els infants en edat preescolar, oscil•la entre el 70-90%, fet que permet defensar la confiança en els records dels infants. Es constata que el percentatge d’emissions incorrectes dels infants en edat preescolar és mínim, al voltant d’un 5-6%. L’estudi remarca la necessitat d’establir perfectament les regles de l’entrevista i, per últim, en destaca la ineficàcia de les tècniques de memòria de l’entrevista cognitiva en els infants de P3 i P4. En els de P5 es comencen a veure beneficis gràcies a la tècnica de la reinstauració contextual (RC), estant les altres tècniques fora de la comprensió i utilització dels infants d’aquestes edats. Interviewing preschoolers victims of sexual abuse and/or domestic abuse: Effectiveness of forensic interviews models 135 preschool children were interviewed with 3 different interview models in order to remember a significant emotional event. Authors conclude that the correct recall of children ranging from 70-90% and the percentage of error messages is 5-6%. It is necessary to fully establish the rules of the interview. The present research highlights the effectiveness of the cognitive interview techniques in children from P3 and P4. Entrevistando niños preescolares víctimas de abuso sexual y/o maltrato familiar: eficacia de los modelos de entrevista forense Se han entrevistado 135 niños preescolares con 3 modelos de entrevista diferentes para recordar un hecho emocionalmente significativo. Se concluye que el recuerdo correcto de los niños oscila entre el 70-90% y el porcentaje de errores de mensajes es del 5-6%. El estudio remarca la necesidad de establecer perfectamente las reglas de la entrevista y se destaca la ineficacia de las técnicas de la entrevista cognitiva en los niños de P3 y P4.

Three-dimensional models of 14α-sterol demethylase (Cyp51A) from Aspergillus lentulus and Aspergillus fumigatus: an insight into differences in voriconazole interaction.

Aspergillus lentulus, an Aspergillus fumigatus sibling species, is increasingly reported in corticosteroid-treated patients. Its clinical significance is unknown, but the fact that A. lentulus shows reduced antifungal susceptibility, mainly to voriconazole, is of serious concern. Heterologous expression of cyp51A from A. fumigatus and A. lentulus was performed in Saccharomyces cerevisiae to assess differences in the interaction of Cyp51A with the azole drugs. The absence of endogenous ERG11 was efficiently complemented in S. cerevisiae by the expression of either Aspergillus cyp51A allele. There was a marked difference between azole minimum inhibitory concentration (MIC) values of the clones expressing each Aspergillus spp. cyp51A. Saccharomyces cerevisiae clones expressing A. lentulus alleles showed higher MICs to all of the azoles tested, supporting the hypothesis that the intrinsic azole resistance of A. lentulus could be associated with Cyp51A. Homology models of A. fumigatus and A. lentulus Cyp51A protein based on the crystal structure of Cyp51p from Mycobacterium tuberculosis in complex with fluconazole were almost identical owing to their mutual high sequence identity. Molecular dynamics (MD) was applied to both three-dimensional protein models to refine the homology modelling and to explore possible differences in the Cyp51A-voriconazole interaction. After 20ns of MD modelling, some critical differences were observed in the putative closed form adopted by the protein upon voriconazole binding. A closer study of the A. fumigatus and A. lentulus voriconazole putative binding site in Cyp51A suggested that some major differences in the protein's BC loop could differentially affect the lock-up of voriconazole, which in turn could correlate with their different azole susceptibility profiles.