Do pseudo-absence selection strategies influence species distribution models and their predictions? An information-theoretic approach based on simulated data

Background Multiple logistic regression is precluded from many practical applications in ecology that aim to predict the geographic distributions of species because it requires absence data, which are rarely available or are unreliable. In order to use multiple logistic regression, many studies have simulated "pseudo-absences" through a number of strategies, but it is unknown how the choice of strategy influences models and their geographic predictions of species. In this paper we evaluate the effect of several prevailing pseudo-absence strategies on the predictions of the geographic distribution of a virtual species whose "true" distribution and relationship to three environmental predictors was predefined. We evaluated the effect of using a) real absences b) pseudo-absences selected randomly from the background and c) two-step approaches: pseudo-absences selected from low suitability areas predicted by either Ecological Niche Factor Analysis: (ENFA) or BIOCLIM. We compared how the choice of pseudo-absence strategy affected model fit, predictive power, and information-theoretic model selection results. Results Models built with true absences had the best predictive power, best discriminatory power, and the "true" model (the one that contained the correct predictors) was supported by the data according to AIC, as expected. Models based on random pseudo-absences had among the lowest fit, but yielded the second highest AUC value (0.97), and the "true" model was also supported by the data. Models based on two-step approaches had intermediate fit, the lowest predictive power, and the "true" model was not supported by the data. Conclusion If ecologists wish to build parsimonious GLM models that will allow them to make robust predictions, a reasonable approach is to use a large number of randomly selected pseudo-absences, and perform model selection based on an information theoretic approach. However, the resulting models can be expected to have limited fit.

Metagenomics, paratransgenesis and the Anopheles microbiome: a portrait of the geographical distribution of the anopheline microbiota based on a meta-analysis of reported taxa

Anophelines harbour a diverse microbial consortium that may represent an extended gene pool for the host. The proposed effects of the insect microbiota span physiological, metabolic and immune processes. Here we synthesise how current metagenomic tools combined with classical culture-dependent techniques provide new insights in the elucidation of the role of the Anopheles-associated microbiota. Many proposed malaria control strategies have been based upon the immunomodulating effects that the bacterial components of the microbiota appear to exert and their ability to express anti-Plasmodium peptides. The number of identified bacterial taxa has increased in the current “omics” era and the available data are mostly scattered or in “tables” that are difficult to exploit. Published microbiota reports for multiple anopheline species were compiled in an Excel® spreadsheet. We then filtered the microbiota data using a continent-oriented criterion and generated a visual correlation showing the exclusive and shared bacterial genera among four continents. The data suggested the existence of a core group of bacteria associated in a stable manner with their anopheline hosts. However, the lack of data from Neotropical vectors may reduce the possibility of defining the core microbiota and understanding the mosquito-bacteria interactive consortium.

Prevalence of human papillomavirus infection, distribution of viral types and risk factors in cervical samples from human immunodeficiency virus-positive women attending three human immunodeficiency virus-acquired immune deficiency syndrome reference centres in northeastern Brazil

Human immunodeficiency virus (HIV)-positive patients have a greater prevalence of coinfection with human papillomavirus (HPV) is of high oncogenic risk. Indeed, the presence of the virus favours intraepithelial squamous cell lesion progression and may induce cancer. The aim of this study was to evaluate the prevalence of HPV infection, distribution of HPV types and risk factors among HIV-positive patients. Cervical samples from 450 HIV-positive patients were analysed with regard to oncotic cytology, colposcopy and HPV presence and type by means of polymerase chain reaction and sequencing. The results were analysed by comparing demographic data and data relating to HPV and HIV infection. The prevalence of HPV was 47.5%. Among the HPV-positive samples, 59% included viral types of high oncogenic risk. Multivariate analysis showed an association between HPV infection and the presence of cytological alterations (p = 0.003), age greater than or equal to 35 years (p = 0.002), number of partners greater than three (p = 0.002), CD4+ lymphocyte count < 200/mm3 (p = 0.041) and alcohol abuse (p = 0.004). Although high-risk HPV was present in the majority of the lesions studied, the low frequency of HPV 16 (3.3%), low occurrence of cervical lesions and preserved immunological state in most of the HIV-positive patients were factors that may explain the low occurrence of precancerous cervical lesions in this population.

Geographic distribution of phlebotomine sandfly species (Diptera: Psychodidae) in Central-West Brazil

This study updates the geographic distributions of phlebotomine species in Central-West Brazil and analyses the climatic factors associated with their occurrence. The data were obtained from the entomology services of the state departments of health in Central-West Brazil, scientific collections and a literature review of articles from 1962-2014. Ecological niche models were produced for sandfly species with more than 20 occurrences using the Maxent algorithm and eight climate variables. In all, 2,803 phlebotomine records for 127 species were analysed. Nyssomyia whitmani,Evandromyia lenti and Lutzomyia longipalpiswere the species with the greatest number of records and were present in all the biomes in Central-West Brazil. The models, which were produced for 34 species, indicated that the Cerrado areas in the central and western regions of Central-West Brazil were climatically more suitable to sandflies. The variables with the greatest influence on the models were the temperature in the coldest months and the temperature seasonality. The results show that phlebotomine species in Central-West Brazil have different geographical distribution patterns and that climate conditions in essentially the entire region favour the occurrence of at least one Leishmania vector species, highlighting the need to maintain or intensify vector control and surveillance strategies.

Driven fragmentation of granular gases

The dynamics of homogeneously heated granular gases which fragment due to particle collisions is analyzed. We introduce a kinetic model which accounts for correlations induced at the grain collisions and analyze both the kinetics and relevant distribution functions these systems develop. The work combines analytical and numerical studies based on direct simulation Monte Carlo calculations. A broad family of fragmentation probabilities is considered, and its implications for the system kinetics are discussed. We show that generically these driven materials evolve asymptotically into a dynamical scaling regime. If the fragmentation probability tends to a constant, the grain number diverges at a finite time, leading to a shattering singularity. If the fragmentation probability vanishes, then the number of grains grows monotonously as a power law. We consider different homogeneous thermostats and show that the kinetics of these systems depends weakly on both the grain inelasticity and driving. We observe that fragmentation plays a relevant role in the shape of the velocity distribution of the particles. When the fragmentation is driven by local stochastic events, the longvelocity tail is essentially exponential independently of the heating frequency and the breaking rule. However, for a Lowe-Andersen thermostat, numerical evidence strongly supports the conjecture that the scaled velocity distribution follows a generalized exponential behavior f (c)~exp (−cⁿ), with n ≈1.2, regarding less the fragmentation mechanisms

Role of Particle Shape on the Stress Propagation in Granular Packings

We present an experimental and numerical study on the influence that particle aspect ratio has on the mechanical and structural properties of granular packings. For grains with maximal symmetry (squares), the stress propagation in the packing localizes forming chainlike forces analogous to the ones observed for spherical grains. This scenario can be understood in terms of stochastic models of aggregation and random multiplicative processes. As the grains elongate, the stress propagation is strongly affected. The interparticle normal force distribution tends toward a Gaussian, and, correspondingly, the force chains spread leading to a more uniform stress distribution reminiscent of the hydrostatic profiles known for standard liquids

A stratified approach for modeling the distribution of a threatened ant species in the Swiss National Park

We present models predicting the potential distribution of a threatened ant species, Formica exsecta Nyl., in the Swiss National Park ( SNP). Data to fit the models have been collected according to a random-stratified design with an equal number of replicates per stratum. The basic aim of such a sampling strategy is to allow the formal testing of biological hypotheses about those factors most likely to account for the distribution of the modeled species. The stratifying factors used in this study were: vegetation, slope angle and slope aspect, the latter two being used as surrogates of solar radiation, considered one of the basic requirements of F. exsecta. Results show that, although the basic stratifying predictors account for more than 50% of the deviance, the incorporation of additional non-spatially explicit predictors into the model, as measured in the field, allows for an increased model performance (up to nearly 75%). However, this was not corroborated by permutation tests. Implementation on a national scale was made for one model only, due to the difficulty of obtaining similar predictors on this scale. The resulting map on the national scale suggests that the species might once have had a broader distribution in Switzerland. Reasons for its particular abundance within the SNP might possibly be related to habitat fragmentation and vegetation transformation outside the SNP boundaries.

Airline Consolidation and the Distribution of Traffic between Primary and Secondary Hubs

Several airline consolidation events have recently been completed both in Europe and in the United States. The model we develop considers two airlines operating hub-and-spoke networks, using different hubs to connect the same spoke airports. We assume the airlines to be vertically differentiated, which allows us to distinguish between primary and secondary hubs. We conclude that this differentiation in air services becomes more accentuated after consolidation, with an increased number of flights being channeled through the primary hub. However, congestion can act as a brake on the concentration of flight frequency in the primary hub following consolidation. Our empirical application involves an analysis of Delta s network following its merger with Northwest. We find evidence consistent with an increase in the importance of Delta s primary hubs at the expense of its secondary airports. We also find some evidence suggesting that the carrier chooses to divert traffic away from those hub airports that were more prone to delays prior to the merger, in particular New York s JFK airport. Keywords: primary hub; secondary hub; airport congestion; airline consolidation; airline networks JEL Classi fication Numbers: D43; L13; L40; L93; R4

Establishment of animal models to analyze the kinetics and distribution of human tumor antigen-specific CD8⁺ T cells.

Many patients develop tumor antigen-specific T cell responses detectable in peripheral blood mononuclear cells (PBMCs) following cancer vaccine. However, measurable tumor regression is observed in a limited number of patients receiving cancer vaccines. There is a need to re-evaluate systemically the immune responses induced by cancer vaccines. Here, we established animal models targeting two human cancer/testis antigens, NY-ESO-1 and MAGE-A4. Cytotoxic T lymphocyte (CTL) epitopes of these antigens were investigated by immunizing BALB/c mice with plasmids encoding the entire sequences of NY-ESO-1 or MAGE-A4. CD8(+) T cells specific for NY-ESO-1 or MAGE-A4 were able to be detected by ELISPOT assays using antigen presenting cells pulsed with overlapping peptides covering the whole protein, indicating the high immunogenicity of these antigens in mice. Truncation of these peptides revealed that NY-ESO-1-specific CD8(+) T cells recognized D(d)-restricted 8mer peptides, NY-ESO-181-88. MAGE-A4-specific CD8(+) T cells recognized D(d)-restricted 9mer peptides, MAGE-A4265-273. MHC/peptide tetramers allowed us to analyze the kinetics and distribution of the antigen-specific immune responses, and we found that stronger antigen-specific CD8(+) T cell responses were required for more effective anti-tumor activity. Taken together, these animal models are valuable for evaluation of immune responses and optimization of the efficacy of cancer vaccines.

A new technique to measure micromotion distribution around a cementless femoral stem.

The interfacial micromotion is closely associated to the long-term success of cementless hip prostheses. Various techniques have been proposed to measure them, but only a few number of points over the stem surface can be measured simultaneously. In this paper, we propose a new technique based on micro-Computer Tomography (μCT) to measure locally the relative interfacial micromotions between the metallic stem and the surrounding femoral bone. Tantalum beads were stuck at the stem surface and spread at the endosteal surface. Relative micromotions between the stem and the endosteal bone surfaces were measured at different loading amplitudes. The estimated error was 10μm and the maximal micromotion was 60μm, in the loading direction, at 1400N. This pilot study provided a local measurement of the micromotions in the 3 direction and at 8 locations on the stem surface simultaneously. This technique could be easily extended to higher loads and a much larger number of points, covering the entire stem surface and providing a quasi-continuous distribution of the 3D interfacial micromotions around the stem. The new measurement method would be very useful to compare the induced micromotions of different stem designs and to optimize the primary stability of cementless total hip arthroplasty.

Bilinear models for spatio-temporal point distribution analysis: application to extrapolation of left ventricular, biventricular and whole heart cardiac dynamics

In this work we describe the usage of bilinear statistical models as a means of factoring the shape variability into two components attributed to inter-subject variation and to the intrinsic dynamics of the human heart. We show that it is feasible to reconstruct the shape of the heart at discrete points in the cardiac cycle. Provided we are given a small number of shape instances representing the same heart atdifferent points in the same cycle, we can use the bilinearmodel to establish this. Using a temporal and a spatial alignment step in the preprocessing of the shapes, around half of the reconstruction errors were on the order of the axial image resolution of 2 mm, and over 90% was within 3.5 mm. From this, weconclude that the dynamics were indeed separated from theinter-subject variability in our dataset.

Characterization of myocardial deformation in patients with different etiologies of left ventricular hypertrophy by using strain distribution from Magnetic Resonance Imaging

Introducción y objetivos. Se ha señalado que, en la miocardiopatía hipertrófica (MCH), la desorganización de las fibras regionales da lugar a segmentos en los que la deformación es nula o está gravemente reducida, y que estos segmentos tienen una distribución no uniforme en el ventrículo izquierdo (VI). Esto contrasta con lo observado en otros tipos de hipertrofia como en el corazón de atleta o la hipertrofia ventricular izquierda hipertensiva (HVI-HT), en los que puede haber una deformación cardiaca anormal, pero nunca tan reducida como para que se observe ausencia de deformación. Así pues, proponemos el empleo de la distribución de los valores de strain para estudiar la deformación en la MCH. Métodos. Con el empleo de resonancia magnética marcada (tagged), reconstruimos la deformación sistólica del VI de 12 sujetos de control, 10 atletas, 12 pacientes con MCH y 10 pacientes con HVI-HT. La deformación se cuantificó con un algoritmo de registro no rígido y determinando los valores de strain sistólico máximo radial y circunferencial en 16 segmentos del VI. Resultados. Los pacientes con MCH presentaron unos valores medios de strain significativamente inferiores a los de los demás grupos. Sin embargo, aunque la deformación observada en los individuos sanos y en los pacientes con HVI-HT se concentraba alrededor del valor medio, en la MCH coexistían segmentos con contracción normal y segmentos con una deformación nula o significativamente reducida, con lo que se producía una mayor heterogeneidad de los valores de strain. Se observaron también algunos segmentos sin deformación incluso en ausencia de fibrosis o hipertrofia. Conclusiones. La distribución de strain caracteriza los patrones específicos de deformación miocárdica en pacientes con diferentes etiologías de la HVI. Los pacientes con MCH presentaron un valor medio de strain significativamente inferior, así como una mayor heterogeneidad de strain (en comparación con los controles, los atletas y los pacientes con HVI-HT), y tenían regiones sin deformación.

Quantitative distribution of parvalbumin, calretinin, and calbindin D-28k immunoreactive neurons in the visual cortex of normal and Alzheimer cases.

The distribution of parvalbumin (PV), calretinin (CR), and calbindin (CB) immunoreactive neurons was studied with the help of an image analysis system (Vidas/Zeiss) in the primary visual area 17 and associative area 18 (Brodmann) of Alzheimer and control brains. In neither of these areas was there a significant difference between Alzheimer and control groups in the mean number of PV, CR, or CB immunoreactive neuronal profiles, counted in a cortical column going from pia to white matter. Significant differences in the mean densities (numbers per square millimeter of cortex) of PV, CR, and CB immunoreactive neuronal profiles were not observed either between groups or areas, but only between superficial, middle, and deep layers within areas 17 and 18. The optical density of the immunoreactive neuropil was also similar in Alzheimer and controls, correlating with the numerical density of immunoreactive profiles in superficial, middle, and deep layers. The frequency distribution of neuronal areas indicated significant differences between PV, CR, and CB immunoreactive neuronal profiles in both areas 17 and 18, with more large PV than CR and CB positive profiles. There were also significantly more small and less large PV and CR immunoreactive neuronal profiles in Alzheimer than in controls. Our data show that, although the brain pathology is moderate to severe, there is no prominent decrease of PV, CR and CB positive neurons in the visual cortex of Alzheimer brains, but only selective changes in neuronal perikarya.

Parameter estimation for the distribution of single cell lag times.

In Quantitative Microbial Risk Assessment, it is vital to understand how lag times of individual cells are distributed over a bacterial population. Such identified distributions can be used to predict the time by which, in a growth-supporting environment, a few pathogenic cells can multiply to a poisoning concentration level. We model the lag time of a single cell, inoculated into a new environment, by the delay of the growth function characterizing the generated subpopulation. We introduce an easy-to-implement procedure, based on the method of moments, to estimate the parameters of the distribution of single cell lag times. The advantage of the method is especially apparent for cases where the initial number of cells is small and random, and the culture is detectable only in the exponential growth phase.

Distribution of hospital admissions in a middle-income country of the African region

Introduction: Few studies have reported the distribution of all hospital admissions at the entire country level in low and middle-income countries (LMICs). We examined this question in Seychelles, a rapidly developing small island state in the Africa region, in which access to health care is provided free of charge to all inhabitants through a national health system and all hospital admissions are routinely registered. Methods: Based on all admissions to all hospitals in Seychelles in 2005-2008, we calculated the distribution of hospital admissions, age at admission, length of stay and bed occupancy (i.e. cumulated number of patients * number of days spent in all hospitals) according to both hospital departments and broad causes of diseases (using codes of the ICD-10 classification of diseases). Results: Bed occupancy was largest in the surgical wards (36.7% of all days spent in all hospitals), followed by the medical wards (24.3%), gynecology/obstetrics wards (18.4%), pediatric wards (11.2%), and psychiatric wards (7.2%). According to broad causes of diseases/conditions, bed occupancy was highest for obstetrics/gynecology conditions (19.9% of all days spent at hospital), mental diseases (8.6%), cardiovascular diseases (8.1%), upper aerodigestive/pulmonary diseases (8%), infectious/parasitic diseases (8%), gastrointestinal diseases (7.2%), and urogenital diseases (6.7%). Adjusted to 100'000 population, 153 hospital beds are needed every day, including 31 for obstetrics/gynecologic conditions, 13 for mental diseases, 12 for cardiovascular diseases, 12 for upper aerodigestive diseases, 12 for infectious/parasitic diseases, and 11 for gastrointestinal diseases. Conclusion: Our findings give a good indication of the overall distribution of admissions according to both hospital departments and broad causes of diseases in a middle-income country. These findings provide important information for health care planning at the national level