Principles of resilient coding for plant ecophysiologists

Plant ecophysiology is founded on a rich body of physical and chemical theory, but it is challenging to connect theory with data in unambiguous, analytically rigorous and reproducible ways. Custom scripts written in computer programming languages (coding) enable plant ecophysiologists to model plant processes and fit models to data reproducibly using advanced statistical techniques. Since many ecophysiologists lack formal programming education, we have yet to adopt a unified set of coding principles and standards that could make coding easier to learn, use and modify. We identify eight principles to help in plant ecophysiologists without much programming experience to write resilient code: (i) standardized nomenclature, (ii) consistency in style, (iii) increased modularity/extensibility for easier editing and understanding, (iv) code scalability for application to large data sets, (v) documented contingencies for code maintenance, (vi) documentation to facilitate user understanding; (vii) extensive tutorials and (viii) unit testing and benchmarking. We illustrate these principles using a new R package, {photosynthesis}, which provides a set of analytical and simulation tools for plant ecophysiology. Our goal with these principles is to advance scientific discovery in plant ecophysiology by making it easier to use code for simulation and data analysis, reproduce results and rapidly incorporate new biological understanding and analytical tools

Crossed products by endomorphisms and reduction of relations in relative Cuntz-Pimsner algebras

Starting from an arbitrary endomorphism \alpha of a unital C*-algebra A we construct a crossed product. It is shown that the natural construction depends not only on the C*-dynamical system (A,\alpha) but also on the choice of an ideal orthogonal to kernel of \alpha. The article gives an explicit description of the internal structure of this crossed product and, in particular, discusses the interrelation between relative Cuntz-Pimsner algebras and partial isometric crossed products. We present a canonical procedure that reduces any given C*-correspondence to the 'smallest' C*-correspondence yielding the same relative Cuntz-Pimsner algebra as the initial one. In the context of crossed products this reduction procedure corresponds to the reduction of C*-dynamical systems and allow us to establish a coincidence between relative Cuntz-Pimsner algebras and crossed products introduced.Comment: The article is based on papers arXiv:math.OA/0703801 and arXiv:math.OA/0704.3811, and in essence forms their unification, refinement and developmen

Structural response of Caribbean dry forests to hurricane winds: a case study from Guanica Forest, Puerto Rico

Tropical dry forests in the Caribbean have an uniquely short, shrubby structure with a high proportion of multiple-stemmed trees compared to dry forests elsewhere in the Neotropics. Previous studies have shown that this structure can arise without the loss of main stems from cutting, grazing, or other human intervention. The Caribbean has a high frequency of hurricanes, so wind may also influence forest stature. Furthermore, these forests also tend to grow on soils with low amounts of available phosphorus, which may also influence structure. The objective of this study was to assess the role of high winds in structuring dry forest, and to determine whether soil nutrient pools influence forest response following hurricane disturbance. Methods: Over 2000 stems in five plots were sampled for hurricane effects within 1 week after Hurricane Georges impacted field sites in 1998. Sprout initiation, growth, and mortality were analysed for 1407 stems for 2 years after the hurricane. Soil nutrient pools were measured at the base of 456 stems to assess association between nutrients and sprout dynamics. Results: Direct effects of the hurricane were minimal, with stem mortality at \u3c 2% and structural damage to stems at 13%, although damage was biased toward stems of larger diameter. Sprouting response was high . over 10 times as many trees had sprouts after the hurricane as before. The number of sprouts on a stem also increased significantly. Sprouting was common on stems that only suffered defoliation or had no visible effects from the hurricane. Sprout survival after 2 years was also high (\u3e 86%). Soil nutrient pools had little effect on forest response as a whole, but phosphorus supply did influence sprout dynamics on four of the more common tree species. Main Conclusions: Hurricanes are able to influence Caribbean tropical dry forest structure by reducing average stem diameter and basal area and generating significant sprouting responses. New sprouts, with ongoing survival, will maintain the high frequency of multi-stemmed trees found in this region. Sprouting is not limited to damaged stems, indicating that trees are responding to other aspects of high winds, such as short-term gravitational displacement or sway. Soil nutrients play a secondary role in sprouting dynamics of a subset of species. The short, shrubby forest structure common to the Caribbean can arise naturally as a response to hurricane winds

Probing nucleobase photoprotection with soft x-rays

Nucleobases absorb strongly in the ultraviolet region, leading to molecular excitation into reactive states. The molecules avoid the photoreactions by funnelling the electronic energy into less reactive states on an ultrafast timescale via non-Born-Oppenheimer dynamics. Current theory on the nucleobase thymine discusses two conflicting pathways for the photoprotective dynamics. We present our first results of our free electron laser based UV-pump soft x-ray-probe study of the photoprotection mechanism of thymine. We use the high spatial sensitivity of the Auger electrons emitted after the soft x-ray pulse induced core ionization. Our transient spetra show two timescales on the order of 200 fs and 5 ps, in agreement with previous (all UV) ultrafast experiments. The timescales appear at different Auger kinetic energies which will help us to decipher the molecular dynamics

Development and external validation of machine learning algorithms for postnatal gestational age estimation using clinical data and metabolomic markers

Background Accurate estimates of gestational age (GA) at birth are important for preterm birth surveillance but can be challenging to obtain in low income countries. Our objective was to develop machine learning models to accurately estimate GA shortly after birth using clinical and metabolomic data. Methods We derived three GA estimation models using ELASTIC NET multivariable linear regression using metabolomic markers from heel-prick blood samples and clinical data from a retrospective cohort of newborns from Ontario, Canada. We conducted internal model validation in an independent cohort of Ontario newborns, and external validation in heel prick and cord blood sample data collected from newborns from prospective birth cohorts in Lusaka, Zambia and Matlab, Bangladesh. Model performance was measured by comparing model-derived estimates of GA to reference estimates from early pregnancy ultrasound. Results Samples were collected from 311 newborns from Zambia and 1176 from Bangladesh. The best-performing model accurately estimated GA within about 6 days of ultrasound estimates in both cohorts when applied to heel prick data (MAE 0.79 weeks (95% CI 0.69, 0.90) for Zambia; 0.81 weeks (0.75, 0.86) for Bangladesh), and within about 7 days when applied to cord blood data (1.02 weeks (0.90, 1.15) for Zambia; 0.95 weeks (0.90, 0.99) for Bangladesh). Conclusions Algorithms developed in Canada provided accurate estimates of GA when applied to external cohorts from Zambia and Bangladesh. Model performance was superior in heel prick data as compared to cord blood data

International consensus on the most useful physical examination tests used by physiotherapists for patients with headache: A Delphi study

Background: A wide range of physical tests have been published for use in the assessment of musculoskeletal dysfunction in patients with headache. Which tests are used depends on a physiotherapist's clinical and scientific background as there is little guidance on the most clinically useful tests. Objectives: To identify which physical examination tests international experts in physiotherapy consider the most clinically useful for the assessment of patients with headache. Design/methods: Delphi survey with pre-specified procedures based on a systematic search of the literature for physical examination tests proposed for the assessment of musculoskeletal dysfunction in patients with headache. Results: Seventeen experts completed all three rounds of the survey. Fifteen tests were included in round one with eleven additional tests suggested by the experts. Finally eleven physical examination tests were considered clinically useful: manual joint palpation, the cranio-cervical flexion test, the cervical flexion-rotation test, active range of cervical movement, head forward position, trigger point palpation, muscle tests of the shoulder girdle, passive physiological intervertebral movements, reproduction and resolution of headache symptoms, screening of the thoracic spine, and combined movement tests. Conclusions: Eleven tests are suggested as a minimum standard for the physical examination of musculoskeletal dysfunctions in patients with headache

Host-pathogen genetic interactions underlie tuberculosis susceptibility in genetically diverse mice

The outcome of an encounter with Mycobacterium tuberculosis (Mtb) depends on the pathogen's ability to adapt to the variable immune pressures exerted by the host. Understanding this interplay has proven difficult, largely because experimentally tractable animal models do not recapitulate the heterogeneity of tuberculosis disease. We leveraged the genetically diverse Collaborative Cross (CC) mouse panel in conjunction with a library of Mtb mutants to create a resource for associating bacterial genetic requirements with host genetics and immunity. We report that CC strains vary dramatically in their susceptibility to infection and produce qualitatively distinct immune states. Global analysis of Mtb transposon mutant fitness (TnSeq) across the CC panel revealed that many virulence pathways are only required in specific host microenvironments, identifying a large fraction of the pathogen's genome that has been maintained to ensure fitness in a diverse population. Both immunological and bacterial traits can be associated with genetic variants distributed across the mouse genome, making the CC a unique population for identifying specific host-pathogen genetic interactions that influence pathogenesis