Overcoming the rare species modelling paradox: test of a novel hierarchical framework with an Iberian endemic plant

Rare species have restricted geographic ranges, habitat specialization, and/or small population sizes. Datasets on rare species distribution usually have few observations, limited spatial accuracy and lack of valid absences; conversely they provide comprehensive views of species distributions allowing to realistically capture most of their realized environmental niche. Rare species are the most in need of predictive distribution modelling but also the most difficult to model. We refer to this contrast as the "rare species modelling paradox" and propose as a solution developing modelling approaches that deal with a sufficiently large set of predictors, ensuring that statistical models aren't overfitted. Our novel approach fulfils this condition by fitting a large number of bivariate models and averaging them with a weighted ensemble approach. We further propose that this ensemble forecasting is conducted within a hierarchic multi-scale framework. We present two ensemble models for a test species, one at regional and one at local scale, each based on the combination of 630 models. In both cases, we obtained excellent spatial projections, unusual when modelling rare species. Model results highlight, from a statistically sound approach, the effects of multiple drivers in a same modelling framework and at two distinct scales. From this added information, regional models can support accurate forecasts of range dynamics under climate change scenarios, whereas local models allow the assessment of isolated or synergistic impacts of changes in multiple predictors. This novel framework provides a baseline for adaptive conservation, management and monitoring of rare species at distinct spatial and temporal scales.

Hierarchical analyses of genetic differentiation in a hybrid zone of Sorex araneus (Insectivora : Soricidae)

Microsatellites are used to unravel the fine-scale genetic structure of a hybrid zone between chromosome races Valais and Cordon of the common shrew (Sorex araneus) located in the French Alps. A total of 269 individuals collected between 1992 and 1995 was typed for seven microsatellite loci. A modified version of the classical multiple correspondence analysis is carried out. This analysis clearly shows the dichotomy between the two races. Several approaches are used to study genetic structuring. Gene flow is clearly reduced between these chromosome races and is estimated at one migrant every two generations using X-statistics and one migrant per generation using F-statistics. Hierarchical F- and R-statistics are compared and their efficiency to detect inter- and intraracial patterns of divergence is discussed. Within-race genetic structuring is significant, but remains weak. F-ST displays similar values on both sides of the hybrid zone, although no environmental barriers are found on the Cordon side, whereas the Valais side is divided by several mountain rivers. We introduce the exact G-test to microsatellite data which proved to be a powerful test to detect genetic differentiation within as well as among races. The genetic background of karyotypic hybrids was compared with the genetic background of pure parental forms using a CRT-MCA. Our results indicate that, without knowledge of the karyotypes, we would not have been able to distinguish these hybrids from karyotypically pure samples.

A hierarchical approach to model parameter optimization for developmental systems.

In the context of Systems Biology, computer simulations of gene regulatory networks provide a powerful tool to validate hypotheses and to explore possible system behaviors. Nevertheless, modeling a system poses some challenges of its own: especially the step of model calibration is often difficult due to insufficient data. For example when considering developmental systems, mostly qualitative data describing the developmental trajectory is available while common calibration techniques rely on high-resolution quantitative data. Focusing on the calibration of differential equation models for developmental systems, this study investigates different approaches to utilize the available data to overcome these difficulties. More specifically, the fact that developmental processes are hierarchically organized is exploited to increase convergence rates of the calibration process as well as to save computation time. Using a gene regulatory network model for stem cell homeostasis in Arabidopsis thaliana the performance of the different investigated approaches is evaluated, documenting considerable gains provided by the proposed hierarchical approach.

Hierarchical categorisation of tags for delicious

In the scenario of social bookmarking, a user browsing the Web bookmarks web pages and assigns free-text labels (i.e., tags) to them according to their personal preferences. In this technical report, we approach one of the practical aspects when it comes to represent users' interests from their tagging activity, namely the categorization of tags into high-level categories of interest. The reason is that the representation of user profiles on the basis of the myriad of tags available on the Web is certainly unfeasible from various practical perspectives; mainly concerning the unavailability of data to reliably, accurately measure interests across such fine-grained categorisation, and, should the data be available, its overwhelming computational intractability. Motivated by this, our study presents the results of a categorization process whereby a collection of tags posted at Delicious #http://delicious.com# are classified into 200 subcategories of interest.

Transhepatic ipsilateral right portal vein embolization extended to segment IV: improving hypertrophy and resection outcomes with spherical particles and coils.

PURPOSE: To analyze outcomes after right portal vein embolization extended to segment IV (right PVE + IV) before extended right hepatectomy, including liver hypertrophy, resection rates, and complications after embolization and resection, and to assess differences in outcomes with two different particulate embolic agents. MATERIALS AND METHODS: Between 1998 and 2004, transhepatic ipsilateral right PVE + IV with particles and coils was performed in 44 patients with malignant hepatobiliary disease, including metastases (n = 24), biliary cancer (n = 14), and hepatocellular carcinoma (n = 6). Right PVE + IV was considered if the future liver remnant (FLR; segments II/III with or without I) was less than 25% of the total estimated liver volume (TELV). Tris-acryl microspheres (100-700 microm; n = 21) or polyvinyl alcohol (PVA) particles (355-1,000 microm; n = 23) were administered in a stepwise fashion. Smaller particles were used to occlude distal branches, followed by larger particles to occlude proximal branches until near-complete stasis. Coils were then placed in secondary portal branches. Computed tomographic volumetry was performed before and 3-4 weeks after right PVE + IV to assess FLR hypertrophy. Liver volumes and postembolization and postoperative outcomes were measured. RESULTS: After right PVE + IV with PVA particles, FLR volume increased 45.5% +/- 40.9% and FLR/TELV ratio increased 6.9% +/- 5.6%. After right PVE + IV with tris-acryl microspheres, FLR volume increased 69.0% +/- 30.7% and FLR/TELV ratio increased 9.7% +/- 3.3%. Differences in FLR volume (P = .0011), FLR/TELV ratio (P = .027), and resection rates (P = .02) were statistically significant. Seventy-one percent of patients underwent extended right hepatectomy (86% after receiving tris-acryl microspheres, 57% after receiving PVA). Thirteen patients (29%) did not undergo resection (extrahepatic spread [n = 9], inadequate hypertrophy [n = 3], other reasons [n = 1]). No patient developed postembolization syndrome or progressive liver insufficiency after embolization or resection. One death after resection occurred as a result of sepsis and hemorrhage. Median hospital stays were 1 day after right PVE + IV and 7 days after resection. CONCLUSION: Transhepatic ipsilateral right PVE + IV with use of particles and coils is a safe, effective method for inducing contralateral hypertrophy before extended right hepatectomy. Embolization with small spherical particles provides improved hypertrophy and resection rates compared with larger, nonspherical particles.

Hierarchical management of carbon sources is regulated similarly by the CbrA/B systems in Pseudomonas aeruginosa and Pseudomonas putida.

The CbrA/B system in pseudomonads is involved in the utilization of carbon sources and carbon catabolite repression (CCR) through the activation of the small RNAs crcZ in Pseudomonas aeruginosa, and crcZ and crcY in Pseudomonas putida. Interestingly, previous works reported that the CbrA/B system activity in P. aeruginosa PAO1 and P. putida KT2442 responded differently to the presence of different carbon sources, thus raising the question of the exact nature of the signal(s) detected by CbrA. Here, we demonstrated that the CbrA/B/CrcZ(Y) signal transduction pathway is similarly activated in the two Pseudomonas species. We show that the CbrA sensor kinase is fully interchangeable between the two species and, moreover, responds similarly to the presence of different carbon sources. In addition, a metabolomics analysis supported the hypothesis that CCR responds to the internal energy status of the cell, as the internal carbon/nitrogen ratio seems to determine CCR and non-CCR conditions. The strong difference found in the 2-oxoglutarate/glutamine ratio between CCR and non-CCR conditions points to the close relationship between carbon and nitrogen availability, or the relationship between the CbrA/B and NtrB/C systems, suggesting that both regulatory systems sense the same sort or interrelated signal.

Qualitative and Hierarchical Analysis of Protective Factors against Illicit Use of Doping Substances in Athletes Calling a National Anti-Doping Phone-Help Service

Evidence of a sport-specific hierarchy of protective factors against doping would thus be a powerful aid in adapting information and prevention campaigns to target the characteristics of specific athlete groups, and especially those athletes most vulnerable for doping control. The contents of phone calls to a free and anonymous national anti-doping service called 'ecoute dopage' were analysed (192 bodybuilders, 124 cyclists and 44 footballers). The results showed that the protective factors that emerged from analysis could be categorised into two groups. The first comprised 'Health concerns', 'Respect for the law' and 'Doping controls from the environment' and the second comprised 'Doubts about the effectiveness of illicit products, 'Thinking skills' and 'Doubts about doctors'. The ranking of the factors for the cyclists differed from that of the other athletes. The ordering of factors was 1) respect for the law, 2) doping controls from the environment, 3) health concerns 4) doubts about doctors, and 5) doubts about the effectiveness illicit products. The results are analysed in terms of the ranking in each athlete group and the consequences on the athletes' experience and relationship to doping. Specific prevention campaigns are proposed to limit doping behaviour in general and for each sport.

Patterns of calcium-binding proteins support parallel and hierarchical organization of human auditory areas.

The human primary auditory cortex (AI) is surrounded by several other auditory areas, which can be identified by cyto-, myelo- and chemoarchitectonic criteria. We report here on the pattern of calcium-binding protein immunoreactivity within these areas. The supratemporal regions of four normal human brains (eight hemispheres) were processed histologically, and serial sections were stained for parvalbumin, calretinin or calbindin. Each calcium-binding protein yielded a specific pattern of labelling, which differed between auditory areas. In AI, defined as area TC [see C. von Economo and L. Horn (1930) Z. Ges. Neurol. Psychiatr.,130, 678-757], parvalbumin labelling was dark in layer IV; several parvalbumin-positive multipolar neurons were distributed in layers III and IV. Calbindin yielded dark labelling in layers I-III and V; it revealed numerous multipolar and pyramidal neurons in layers II and III. Calretinin labelling was lighter than that of parvalbumin or calbindin in AI; calretinin-positive bipolar and bitufted neurons were present in supragranular layers. In non-primary auditory areas, the intensity of labelling tended to become progressively lighter while moving away from AI, with qualitative differences between the cytoarchitectonically defined areas. In analogy to non-human primates, our results suggest differences in intrinsic organization between auditory areas that are compatible with parallel and hierarchical processing of auditory information.

Spatial distribution of a spherical gall (Hymenoptera, Eulophidae) on Caryocar brasiliense (Caryocaraceae)

We obtained the first data on spatial distribution of a spherical galling insect (Hymenoptera, Eulophidae) at the Caryocar brasiliense Camb. (Caryocaraceae) tree level. This work was developed in two pastures in Montes Claros, Minas Gerais State, Brazil. The areas studied were: pasture 1 (in activity) and pasture 2 (abandoned pasture = savanna in recovery). We evaluated the distribution of spherical galls in: foliage orientation (slope), among leaves (border and interior of the tree crown), among leaflets (right, central, left), distal, median, and proximal as well as border, central area, and adjacent to the mid leaf vein of the leaflet, and difference between areas in 10 infested trees per area. The smaller number of spherical gall/leaflet was observed in pasture 1 than in pasture 2. More spherical galls were found on the northern in pasture 1, but in the pasture 2, the lower spherical galls were observed on the northeast than other slopes. The average number of spherical galls did not differ statistically among the three leaflets of C. brasiliense in pasture 2. However, in pasture 1, we observed highest number of spherical galls in the central leaflet. More spherical galls were found in the border than interior of the tree crown. The average number of spherical galls did not differ statistically among the longitudinal region on leaflet of C. brasiliense. The spherical gall insect preferred to colonize the leaf margin than the central portion or near mid vein on transversal regions on a leaflet.

Hierarchical location-allocation models for congested systems

In this paper we address the issue of locating hierarchical facilities in the presence of congestion. Two hierarchical models are presented, where lower level servers attend requests first, and then, some of the served customers are referred to higher level servers. In the first model, the objective is to find the minimum number of servers and theirlocations that will cover a given region with a distance or time standard. The second model is cast as a Maximal Covering Location formulation. A heuristic procedure is then presented together with computational experience. Finally, some extensions of these models that address other types of spatial configurations are offered.

HIERFSTAT, a package for R to compute and test hierarchical F-statistics

The package HIERFSTAT for the statistical software R, created by the R Development Core Team, allows the estimate of hierarchical F-statistics from a hierarchy with any numbers of levels. In addition, it allows testing the statistical significance of population differentiation for these different levels, using a generalized likelihood-ratio test. The package HIERFSTAT is available at http://www.unil.ch/popgen/softwares/hierfstat.htm.

Sample size and statistical power in the hierarchical analysis of variance: applications in morphometry of the nervous system.

Analysis of variance is commonly used in morphometry in order to ascertain differences in parameters between several populations. Failure to detect significant differences between populations (type II error) may be due to suboptimal sampling and lead to erroneous conclusions; the concept of statistical power allows one to avoid such failures by means of an adequate sampling. Several examples are given in the morphometry of the nervous system, showing the use of the power of a hierarchical analysis of variance test for the choice of appropriate sample and subsample sizes. In the first case chosen, neuronal densities in the human visual cortex, we find the number of observations to be of little effect. For dendritic spine densities in the visual cortex of mice and humans, the effect is somewhat larger. A substantial effect is shown in our last example, dendritic segmental lengths in monkey lateral geniculate nucleus. It is in the nature of the hierarchical model that sample size is always more important than subsample size. The relative weight to be attributed to subsample size thus depends on the relative magnitude of the between observations variance compared to the between individuals variance.

A program for the computation of power and determination of sample size in hierarchical experimental designs.

We have devised a program that allows computation of the power of F-test, and hence determination of appropriate sample and subsample sizes, in the context of the one-way hierarchical analysis of variance with fixed effects. The power at a fixed alternative is an increasing function of the sample size and of the subsample size. The program makes it easy to obtain the power of F-test for a range of values of sample and subsample sizes, and therefore the appropriate sizes based on a desired power. The program can be used for the 'ordinary' case of the one-way analysis of variance, as well as for hierarchical analysis of variance with two stages of sampling. Examples are given of the practical use of the program.