Estimating and Modelling Bias of the Hierarchical Partitioning Public-Domain Software: Implications in Environmental Management and Conservation

BACKGROUND: Hierarchical partitioning (HP) is an analytical method of multiple regression that identifies the most likely causal factors while alleviating multicollinearity problems. Its use is increasing in ecology and conservation by its usefulness for complementing multiple regression analysis. A public-domain software "hier.part package" has been developed for running HP in R software. Its authors highlight a "minor rounding error" for hierarchies constructed from >9 variables, however potential bias by using this module has not yet been examined. Knowing this bias is pivotal because, for example, the ranking obtained in HP is being used as a criterion for establishing priorities of conservation. METHODOLOGY/PRINCIPAL FINDINGS: Using numerical simulations and two real examples, we assessed the robustness of this HP module in relation to the order the variables have in the analysis. Results indicated a considerable effect of the variable order on the amount of independent variance explained by predictors for models with >9 explanatory variables. For these models the nominal ranking of importance of the predictors changed with variable order, i.e. predictors declared important by its contribution in explaining the response variable frequently changed to be either most or less important with other variable orders. The probability of changing position of a variable was best explained by the difference in independent explanatory power between that variable and the previous one in the nominal ranking of importance. The lesser is this difference, the more likely is the change of position. CONCLUSIONS/SIGNIFICANCE: HP should be applied with caution when more than 9 explanatory variables are used to know ranking of covariate importance. The explained variance is not a useful parameter to use in models with more than 9 independent variables. The inconsistency in the results obtained by HP should be considered in future studies as well as in those already published. Some recommendations to improve the analysis with this HP module are given

Potent Innate Immune Response to Pathogenic Leptospira in Human Whole Blood

Background: Leptospirosis is caused by pathogenic spirochetes of the genus Leptospira. The bacteria enter the human body via abraded skin or mucous membranes and may disseminate throughout. In general the clinical picture is mild but some patients develop rapidly progressive, severe disease with a high case fatality rate. Not much is known about the innate immune response to leptospires during haematogenous dissemination. Previous work showed that a human THP-1 cell line recognized heat-killed leptospires and leptospiral LPS through TLR2 instead of TLR4. The LPS of virulent leptospires displayed a lower potency to trigger TNF production by THP-1 cells compared to LPS of non-virulent leptospires. Methodology/Principal Findings: We investigated the host response and killing of virulent and non-virulent Leptospira of different serovars by human THP-1 cells, human PBMC's and human whole blood. Virulence of each leptospiral strain was tested in a well accepted standard guinea pig model. Virulent leptospires displayed complement resistance in human serum and whole blood while in-vitro attenuated non-virulent leptospires were rapidly killed in a complement dependent manner. In vitro stimulation of THP-1 and PBMC's with heat-killed and living leptospires showed differential serovar and cell type dependence of cytokine induction. However, at low, physiological, leptospiral dose, living virulent complement resistant strains were consistently more potent in whole blood stimulations than the corresponding non-virulent complement sensitive strains. At higher dose living virulent and non-virulent leptospires were equipotent in whole blood. Inhibition of different TLRs indicated that both TLR2 and TLR4 as well as TLR5 play a role in the whole blood cytokine response to living leptospires. Conclusions/Significance: Thus, in a minimally altered system as human whole blood, highly virulent Leptospira are potent inducers of the cytokine response

A JWST/MIRI and NIRCam Analysis of the Young Stellar Object Population in the Spitzer I region of NGC 6822

We present an imaging survey of the Spitzer~I star-forming region in NGC 6822 conducted with the NIRCam and MIRI instruments onboard JWST. Located at a distance of 490 kpc, NGC 6822 is the nearest non-interacting low-metallicity ($\sim$0.2 $Z_{\odot}$) dwarf galaxy. It hosts some of the brightest known HII regions in the local universe, including recently discovered sites of highly-embedded active star formation. Of these, Spitzer I is the youngest and most active, and houses 90 color-selected candidate young stellar objects (YSOs) identified from Spitzer Space Telescope observations. We revisit the YSO population of Spitzer~I with these new JWST observations. By analyzing color-magnitude diagrams (CMDs) constructed with NIRCam and MIRI data, we establish color selection criteria and construct spectral energy distributions (SEDs) to identify candidate YSOs and characterize the full population of young stars, from the most embedded phase to the more evolved stages. In this way, we have identified 129 YSOs in Spitzer I. Comparing to previous Spitzer studies of the NGC 6822 YSO population, we find that the YSOs we identify are fainter and less massive, indicating that the improved resolution of JWST allows us to resolve previously blended sources into individual stars.Comment: 17 pages, 9 figures, 2 tables, to be submitted to ApJ, comments welcom

JWST MIRI and NIRCam Unveil Previously Unseen Infrared Stellar Populations in NGC 6822

NGC 6822 is a nearby (\sim490 kpc) non-interacting low-metallicity (0.2 Z_\odot) dwarf galaxy which hosts several prominent Hii regions, including sites of highly embedded active star formation. In this work, we present an imaging survey of NGC 6822 conducted with the NIRCam and MIRI instruments onboard JWST. We provide a description of the data reduction, source extraction, and stellar population identifications from combined near- and mid-infrared (IR) photometry. Our NIRCam observations reach seven magnitudes deeper than previous JHKs surveys of this galaxy, which were sensitive to just below the tip of the red giant branch (TRGB). These JWST observations thus reveal for the first time in the near-IR the red clump stellar population and extend nearly three magnitudes deeper. In the mid-IR, we observe roughly two magnitudes below the TRGB with the MIRI F770W and F1000W filters. With these improvements in sensitivity, we produce a catalogue of \sim900,000 point sources over an area of \sim 6.0 x 4.3 arcmin2. We present several NIRCam and MIRI colour-magnitude diagrams and discuss which colour combinations provide useful separations of various stellar populations to aid in future JWST observation planning. Finally, we find populations of carbon- and oxygen-rich asymptotic giant branch stars which will assist in improving our understanding of dust production in low-metallicity, early Universe analogue galaxies

The identification of informative genes from multiple datasets with increasing complexity

Background In microarray data analysis, factors such as data quality, biological variation, and the increasingly multi-layered nature of more complex biological systems complicates the modelling of regulatory networks that can represent and capture the interactions among genes. We believe that the use of multiple datasets derived from related biological systems leads to more robust models. Therefore, we developed a novel framework for modelling regulatory networks that involves training and evaluation on independent datasets. Our approach includes the following steps: (1) ordering the datasets based on their level of noise and informativeness; (2) selection of a Bayesian classifier with an appropriate level of complexity by evaluation of predictive performance on independent data sets; (3) comparing the different gene selections and the influence of increasing the model complexity; (4) functional analysis of the informative genes. Results In this paper, we identify the most appropriate model complexity using cross-validation and independent test set validation for predicting gene expression in three published datasets related to myogenesis and muscle differentiation. Furthermore, we demonstrate that models trained on simpler datasets can be used to identify interactions among genes and select the most informative. We also show that these models can explain the myogenesis-related genes (genes of interest) significantly better than others (P < 0.004) since the improvement in their rankings is much more pronounced. Finally, after further evaluating our results on synthetic datasets, we show that our approach outperforms a concordance method by Lai et al. in identifying informative genes from multiple datasets with increasing complexity whilst additionally modelling the interaction between genes. Conclusions We show that Bayesian networks derived from simpler controlled systems have better performance than those trained on datasets from more complex biological systems. Further, we present that highly predictive and consistent genes, from the pool of differentially expressed genes, across independent datasets are more likely to be fundamentally involved in the biological process under study. We conclude that networks trained on simpler controlled systems, such as in vitro experiments, can be used to model and capture interactions among genes in more complex datasets, such as in vivo experiments, where these interactions would otherwise be concealed by a multitude of other ongoing events

Surgical experience and identification of errors in laparoscopic cholecystectomy.

BACKGROUND: Surgical errors are acts or omissions resulting in negative consequences and/or increased operating time. This study describes surgeon-reported errors in laparoscopic cholecystectomy. METHODS: Intraoperative videos were uploaded and annotated on Touch SurgeryTM Enterprise. Participants evaluated videos for severity using a 10-point intraoperative cholecystitis grading score, and errors using Observational Clinical Human Reliability Assessment, which includes skill, consequence, and mechanism classifications. RESULTS: Nine videos were assessed by 8 participants (3 junior (specialist trainee (ST) 3-5), 2 senior trainees (ST6-8), and 3 consultants). Participants identified 550 errors. Positive relationships were seen between total operating time and error count (r2 = 0.284, P < 0.001), intraoperative grade score and error count (r2 = 0.578, P = 0.001), and intraoperative grade score and total operating time (r2 = 0.157, P < 0.001). Error counts differed significantly across intraoperative phases (H(6) = 47.06, P < 0.001), most frequently at dissection of the hepatocystic triangle (total 282; median 33.5 (i.q.r. 23.5-47.8, range 15-63)), ligation/division of cystic structures (total 124; median 13.5 (i.q.r. 12-19.3, range 10-26)), and gallbladder dissection (total 117; median 14.5 (i.q.r. 10.3-18.8, range 6-26)). There were no significant differences in error counts between juniors, seniors, and consultants (H(2) = 0.03, P = 0.987). Errors were classified differently. For dissection of the hepatocystic triangle, thermal injuries (50 in total) were frequently classified as executional, consequential errors; trainees classified thermal injuries as step done with excessive force, speed, depth, distance, time or rotation (29 out of 50), whereas consultants classified them as incorrect orientation (6 out of 50). For ligation/division of cystic structures, inappropriate clipping (60 errors in total), procedural errors were reported by junior trainees (6 out of 60), but not consultants. For gallbladder dissection, inappropriate dissection (20 errors in total) was reported in incorrect planes by consultants and seniors (6 out of 20), but not by juniors. Poor economy of movement (11 errors in total) was reported more by consultants (8 out of 11) than trainees (3 out of 11). CONCLUSION: This study suggests that surgical experience influences error interpretation, but the benefits for surgical training are currently unclear