Measuring ecological niche overlap from occurrence and spatial environmental data

1. Aim - Concerns over how global change will influence species distributions, in conjunction with increased emphasis on understanding niche dynamics in evolutionary and community contexts, highlight the growing need for robust methods to quantify niche differences between or within taxa. We propose a statistical framework to describe and compare environmental niches from occurrence and spatial environmental data.¦2. Location - Europe, North America, South America¦3. Methods - The framework applies kernel smoothers to densities of species occurrence in gridded environmental space to calculate metrics of niche overlap and test hypotheses regarding niche conservatism. We use this framework and simulated species with predefined distributions and amounts of niche overlap to evaluate several ordination and species distribution modeling techniques for quantifying niche overlap. We illustrate the approach with data on two well-studied invasive species.¦4. Results - We show that niche overlap can be accurately detected with the framework when variables driving the distributions are known. The method is robust to known and previously undocumented biases related to the dependence of species occurrences on the frequency of environmental conditions that occur across geographic space. The use of a kernel smoother makes the process of moving from geographical space to multivariate environmental space independent of both sampling effort and arbitrary choice of resolution in environmental space. However, the use of ordination and species distribution model techniques for selecting, combining and weighting variables on which niche overlap is calculated provide contrasting results.¦5. Main conclusions - The framework meets the increasing need for robust methods to quantify niche differences. It is appropriate to study niche differences between species, subspecies or intraspecific lineages that differ in their geographical distributions. Alternatively, it can be used to measure the degree to which the environmental niche of a species or intraspecific lineage has changed over time.

Analysis, modelling and classification of geospatial data using machine learning

The research considers the problem of spatial data classification using machine learning algorithms: probabilistic neural networks (PNN) and support vector machines (SVM). As a benchmark model simple k-nearest neighbor algorithm is considered. PNN is a neural network reformulation of well known nonparametric principles of probability density modeling using kernel density estimator and Bayesian optimal or maximum a posteriori decision rules. PNN is well suited to problems where not only predictions but also quantification of accuracy and integration of prior information are necessary. An important property of PNN is that they can be easily used in decision support systems dealing with problems of automatic classification. Support vector machine is an implementation of the principles of statistical learning theory for the classification tasks. Recently they were successfully applied for different environmental topics: classification of soil types and hydro-geological units, optimization of monitoring networks, susceptibility mapping of natural hazards. In the present paper both simulated and real data case studies (low and high dimensional) are considered. The main attention is paid to the detection and learning of spatial patterns by the algorithms applied.

Surface granularity as a discriminating feature of illicit tablets

In this paper we propose an innovative methodology for automated profiling of illicit tablets bytheir surface granularity; a feature previously unexamined for this purpose. We make use of the tinyinconsistencies at the tablet surface, referred to as speckles, to generate a quantitative granularity profileof tablets. Euclidian distance is used as a measurement of (dis)similarity between granularity profiles.The frequency of observed distances is then modelled by kernel density estimation in order to generalizethe observations and to calculate likelihood ratios (LRs). The resulting LRs are used to evaluate thepotential of granularity profiles to differentiate between same-batch and different-batches tablets.Furthermore, we use the LRs as a similarity metric to refine database queries. We are able to derivereliable LRs within a scope that represent the true evidential value of the granularity feature. Thesemetrics are used to refine candidate hit-lists form a database containing physical features of illicittablets. We observe improved or identical ranking of candidate tablets in 87.5% of cases when granularityis considered.

Integration of local-scale hydrological and regional-scale geophysical based on a nonlinear Bayesian sequential simulation approach

Significant progress has been made with regard to the quantitative integration of geophysical and hydrological data at the local scale. However, extending the corresponding approaches to the regional scale represents a major, and as-of-yet largely unresolved, challenge. To address this problem, we have developed a downscaling procedure based on a non-linear Bayesian sequential simulation approach. The basic objective of this algorithm is to estimate the value of the sparsely sampled hydraulic conductivity at non-sampled locations based on its relation to the electrical conductivity, which is available throughout the model space. The in situ relationship between the hydraulic and electrical conductivities is described through a non-parametric multivariate kernel density function. This method is then applied to the stochastic integration of low-resolution, re- gional-scale electrical resistivity tomography (ERT) data in combination with high-resolution, local-scale downhole measurements of the hydraulic and electrical conductivities. Finally, the overall viability of this downscaling approach is tested and verified by performing and comparing flow and transport simulation through the original and the downscaled hydraulic conductivity fields. Our results indicate that the proposed procedure does indeed allow for obtaining remarkably faithful estimates of the regional-scale hydraulic conductivity structure and correspondingly reliable predictions of the transport characteristics over relatively long distances.

Robust Estimators of the Generalized Log-Gamma Distribution

We propose robust estimators of the generalized log-gamma distribution and, more generally, of location-shape-scale families of distributions. A (weighted) Q tau estimator minimizes a tau scale of the differences between empirical and theoretical quantiles. It is n(1/2) consistent; unfortunately, it is not asymptotically normal and, therefore, inconvenient for inference. However, it is a convenient starting point for a one-step weighted likelihood estimator, where the weights are based on a disparity measure between the model density and a kernel density estimate. The one-step weighted likelihood estimator is asymptotically normal and fully efficient under the model. It is also highly robust under outlier contamination. Supplementary materials are available online.

Study of the impact of tissue density heterogeneities on 3-dimensional abdominal dosimetry: comparison between dose kernel convolution and direct monte carlo methods.

Dose kernel convolution (DK) methods have been proposed to speed up absorbed dose calculations in molecular radionuclide therapy. Our aim was to evaluate the impact of tissue density heterogeneities (TDH) on dosimetry when using a DK method and to propose a simple density-correction method. METHODS: This study has been conducted on 3 clinical cases: case 1, non-Hodgkin lymphoma treated with (131)I-tositumomab; case 2, a neuroendocrine tumor treatment simulated with (177)Lu-peptides; and case 3, hepatocellular carcinoma treated with (90)Y-microspheres. Absorbed dose calculations were performed using a direct Monte Carlo approach accounting for TDH (3D-RD), and a DK approach (VoxelDose, or VD). For each individual voxel, the VD absorbed dose, D(VD), calculated assuming uniform density, was corrected for density, giving D(VDd). The average 3D-RD absorbed dose values, D(3DRD), were compared with D(VD) and D(VDd), using the relative difference Δ(VD/3DRD). At the voxel level, density-binned Δ(VD/3DRD) and Δ(VDd/3DRD) were plotted against ρ and fitted with a linear regression. RESULTS: The D(VD) calculations showed a good agreement with D(3DRD). Δ(VD/3DRD) was less than 3.5%, except for the tumor of case 1 (5.9%) and the renal cortex of case 2 (5.6%). At the voxel level, the Δ(VD/3DRD) range was 0%-14% for cases 1 and 2, and -3% to 7% for case 3. All 3 cases showed a linear relationship between voxel bin-averaged Δ(VD/3DRD) and density, ρ: case 1 (Δ = -0.56ρ + 0.62, R(2) = 0.93), case 2 (Δ = -0.91ρ + 0.96, R(2) = 0.99), and case 3 (Δ = -0.69ρ + 0.72, R(2) = 0.91). The density correction improved the agreement of the DK method with the Monte Carlo approach (Δ(VDd/3DRD) < 1.1%), but with a lesser extent for the tumor of case 1 (3.1%). At the voxel level, the Δ(VDd/3DRD) range decreased for the 3 clinical cases (case 1, -1% to 4%; case 2, -0.5% to 1.5%, and -1.5% to 2%). No more linear regression existed for cases 2 and 3, contrary to case 1 (Δ = 0.41ρ - 0.38, R(2) = 0.88) although the slope in case 1 was less pronounced. CONCLUSION: This study shows a small influence of TDH in the abdominal region for 3 representative clinical cases. A simple density-correction method was proposed and improved the comparison in the absorbed dose calculations when using our voxel S value implementation.

Skeletal Muscle Mitochondrial and Lipid Droplet Volume Density: Validity of Electron Microscopy Point-counting Measurements

Mitochondrial (M) and lipid droplet (L) volume density (vd) are often used in exercise research. Vd is the volume of muscle occupied by M and L. The means of calculating these percents are accomplished by applying a grid to a 2D image taken with transmission electron microscopy; however, it is not known which grid best predicts these values. PURPOSE: To determine the grid with the least variability of Mvd and Lvd in human skeletal muscle. METHODS: Muscle biopsies were taken from vastus lateralis of 10 healthy adults, trained (N=6) and untrained (N=4). Samples of 5-10mg were fixed in 2.5% glutaraldehyde and embedded in EPON. Longitudinal sections of 60 nm were cut and 20 images were taken at random at 33,000x magnification. Vd was calculated as the number of times M or L touched two intersecting grid lines (called a point) divided by the total number of points using 3 different sizes of grids with squares of 1000x1000nm sides (corresponding to 1µm2), 500x500nm (0.25µm2) and 250x250nm (0.0625µm2). Statistics included coefficient of variation (CV), 1 way-BS ANOVA and spearman correlations. RESULTS: Mean age was 67 ± 4 yo, mean VO2peak 2.29 ± 0.70 L/min and mean BMI 25.1 ± 3.7 kg/m2. Mean Mvd was 6.39% ± 0.71 for the 1000nm squares, 6.01% ± 0.70 for the 500nm and 6.37% ± 0.80 for the 250nm. Lvd was 1.28% ± 0.03 for the 1000nm, 1.41% ± 0.02 for the 500nm and 1.38% ± 0.02 for the 250nm. The mean CV of the three grids was 6.65% ±1.15 for Mvd with no significant differences between grids (P>0.05). Mean CV for Lvd was 13.83% ± 3.51, with a significant difference between the 1000nm squares and the two other grids (P<0.05). The 500nm squares grid showed the least variability between subjects. Mvd showed a positive correlation with VO2peak (r = 0.89, p < 0.05) but not with weight, height, or age. No correlations were found with Lvd. CONCLUSION: Different size grids have different variability in assessing skeletal muscle Mvd and Lvd. The grid size of 500x500nm (240 points) was more reliable than 1000x1000nm (56 points). 250x250nm (1023 points) did not show better reliability compared with the 500x500nm, but was more time consuming. Thus, choosing a grid with square size of 500x500nm seems the best option. This is particularly relevant as most grids used in the literature are either 100 points or 400 points without clear information on their square size.

Density-dependent pollen limitation and reproductive assurance in a wind-pollinated herb with contrasting sexual systems

1. Wind pollination is thought to have evolved in response to selection for mechanisms to promote pollination success, when animal pollinators become scarce or unreliable. We might thus expect wind-pollinated plants to be less prone to pollen limitation than their insect-pollinated counterparts. Yet, if pollen loads on stigmas of wind-pollinated species decline with distance from pollen donors, seed set might nevertheless be pollen-limited in populations of plants that cannot self-fertilize their progeny, but not in self-compatible hermaphroditic populations.2. Here, we test this hypothesis by comparing pollen limitation between dioecious and hermaphroditic (monoecious) populations of the wind-pollinated herb Mercurialis annua.3. In natural populations, seed set was pollen-limited in low-density patches of dioecious, but not hermaphroditic, M. annua, a finding consistent with patterns of distance-dependent seed set by females in an experimental array. Nevertheless, seed set was incomplete in both dioecious and hermaphroditic populations, even at high local densities. Further, both factors limited the seed set of females and hermaphrodites, after we manipulated pollen and resource availability in a common garden experiment.4. Synthesis. Our results are consistent with the idea that pollen limitation plays a role in the evolution of combined vs. separate sexes in M. annua. Taken together, they point to the potential importance of pollen transfer between flowers on the same plant (geitonogamy) by wind as a mechanism of reproductive assurance and to the dual roles played by pollen and resource availability in limiting seed set. Thus, seed set can be pollen-limited in sparse populations of a wind-pollinated species, where mates are rare or absent, having potentially important demographic and evolutionary implications.

The effects of nutrition, puberty and dancing on bone density in adolescent ballet dancers.

Ballet dancers have on average a low bone mineral content (BMC), with elevated fracture-risk, low body mass index (BMI) for age (body mass index, kg/m2), low energy intake, and delayed puberty. This study aims at a better understanding of the interactions of these factors, especially with regard to nutrition. During a competition for pre-professional dancers we examined 127 female participants (60 Asians, 67 Caucasians). They averaged 16.7 years of age, started dancing at 5.8 years, and danced 22 hours/week. Assessments were made for BMI, BMC (DXA), and bone mineral apparent density (BMAD) at the lumbar spine and femoral neck, pubertal stage (Tanner score), and nutritional status (EAT-40 questionnaire and a qualitative three-day dietary record). BMI for age was found to be normal in only 42.5% of the dancers, while 15.7% had a more or less severe degree of thinness (12.6% Grade2 and 3.1% Grade 3 thinness). Menarche was late (13.9 years, range 11 to 16.8 years). Food intake, evaluated by number of consumed food portions, was below the recommendations for a normally active population in all food groups except animal proteins, where the intake was more than twice the recommended amount. In this population, with low BMI and intense exercise, BMC was low and associated with nutritional factors; dairy products had a positive and non-dairy proteins a negative influence. A positive correlation between BMAD and years since menarche confirmed the importance of exposure to estrogens and the negative impact of delayed puberty. Because of this and the probable negative influence of a high intake of non-dairy proteins, such as meat, fish, and eggs, and the positive association with a high dairy intake, ballet schools should promote balanced diets and normal weight and should recognize and help dancers avoid eating disorders and delayed puberty caused by extensive dancing and inadequate nutrition.

Bone microarchitecture assessed by TBS predicts osteoporotic fractures independent of bone density: The manitoba study.

The measurement of BMD by dual-energy X-ray absorptiometry (DXA) is the "gold standard" for diagnosing osteoporosis but does not directly reflect deterioration in bone microarchitecture. The trabecular bone score (TBS), a novel gray-level texture measurement that can be extracted from DXA images, correlates with 3D parameters of bone microarchitecture. Our aim was to evaluate the ability of lumbar spine TBS to predict future clinical osteoporotic fractures. A total of 29,407 women 50 years of age or older at the time of baseline hip and spine DXA were identified from a database containing all clinical results for the Province of Manitoba, Canada. Health service records were assessed for the incidence of nontraumatic osteoporotic fracture codes subsequent to BMD testing (mean follow-up 4.7 years). Lumbar spine TBS was derived for each spine DXA examination blinded to clinical parameters and outcomes. Osteoporotic fractures were identified in 1668 (5.7%) women, including 439 (1.5%) spine and 293 (1.0%) hip fractures. Significantly lower spine TBS and BMD were identified in women with major osteoporotic, spine, and hip fractures (all p < 0.0001). Spine TBS and BMD predicted fractures equally well, and the combination was superior to either measurement alone (p < 0.001). Spine TBS predicts osteoporotic fractures and provides information that is independent of spine and hip BMD. Combining the TBS trabecular texture index with BMD incrementally improves fracture prediction in postmenopausal women. © 2011 American Society for Bone and Mineral Research.

Photodynamic induced uptake of liposomal doxorubicin to rat lung tumors parallels tumor vascular density.

BACKGROUND: Visudyne®-mediated photodynamic therapy (PDT) at low drug/light conditions has shown to selectively enhance the uptake of liposomal doxorubicin in subpleural localized sarcoma tumors grown on rodent lungs without causing morphological alterations of the lung. The present experiments explore the impact of low-dose PDT on liposomal doxorubicin (Liporubicin™) uptake to different tumor types grown on rodent lungs. MATERIAL AND METHODS: Three groups of Fischer rats underwent subpleural generation of sarcoma, mesothelioma, or adenocarcinoma tumors on the left lung. At least five animals of each group (sarcoma, n = 5; mesothelioma, n = 7; adenocarcinoma, n = 5) underwent intraoperative low-dose (10 J/cm(2) at 35 mW/cm(2) ) PDT with 0.0625 mg/kg Visudyne® of the tumor and the lower lobe. This was followed by intravenous (IV) administration of 400 µg Liporubicin™. After a circulation time of 60 min, the tumor-bearing lung was processed for HPLC analyses. At least five animals per group underwent the same procedure but without PDT (sarcoma, n = 5; mesothelioma, n = 5; adenocarcinoma, n = 6). Five untreated animals per group underwent CD31 immunostaining of their tumors with histomorphometrical assessment of the tumor vascularization. RESULTS: Low-dose PDT significantly enhanced Liporubicin™ uptake to all tumor types (sarcoma, P = 0.0007; mesothelioma, P = 0.001; adenocarcinoma, P = 0.02) but not to normal lung tissue compared to IV drug administration alone. PDT led to a significantly increased ratio of tumor to lung tissue drug uptake for all three tumor types (P < 0.05). However, the tumor drug uptake varied between tumor types and paralleled tumor vascular density. The vascular density was significantly higher in sarcoma than in adenocarcinoma (P < 0.001) and mesothelioma (P < 0.001), whereas there was no significant difference between adenocarcinoma and mesothelioma. CONCLUSION: Low-dose Visudyne®-mediated PDT selectively enhances the uptake of systemically administered liposomal doxorubicin in tumors without affecting the drug uptake to normal lung. However, drug uptake varied significantly between tumor types and paralleled tumor vascular density.