gems: An R Package for Simulating from Disease Progression Models

Mathematical models of disease progression predict disease outcomes and are useful epidemiological tools for planners and evaluators of health interventions. The R package gems is a tool that simulates disease progression in patients and predicts the effect of different interventions on patient outcome. Disease progression is represented by a series of events (e.g., diagnosis, treatment and death), displayed in a directed acyclic graph. The vertices correspond to disease states and the directed edges represent events. The package gems allows simulations based on a generalized multistate model that can be described by a directed acyclic graph with continuous transition-specific hazard functions. The user can specify an arbitrary hazard function and its parameters. The model includes parameter uncertainty, does not need to be a Markov model, and may take the history of previous events into account. Applications are not limited to the medical field and extend to other areas where multistate simulation is of interest. We provide a technical explanation of the multistate models used by gems, explain the functions of gems and their arguments, and show a sample application.

Un modèle animal simple pour l'apprentissage des techniques de microanastomoses vasculaires de congruences diefférentes

BACKGROUND Since the pioneering work of Jacobson and Suarez, microsurgery has steadily progressed and is now used in all surgical specialities, particularly in plastic surgery. Before performing clinical procedures it is necessary to learn the basic techniques in the laboratory. OBJECTIVE To assess an animal model, thereby circumventing the following issues: ethical rules, cost, anesthesia and training time. METHODS Between July 2012 and September 2012, 182 earthworms were used for 150 microsurgical trainings to simulate discrepancy microanastomoses. Training was undertaken over 10 weekly periods. Each training session included 15 simulations of microanastomoses performed using the Harashina technique (earthworm diameters >1.5 mm [n=5], between 1.0 mm and 1.5 mm [n=5], and <1.0 mm [n=5]). The technique is presented and documented. A linear model with main variable as the number of the week (as a numeric covariate) and the size of the animal (as a factor) was used to determine the trend in time of anastomosis over subsequent weeks as well as differences between the different size groups. RESULTS The linear model showed a significant trend (P<0.001) in time of anastomosis in the course of the training, as well as significant differences (P<0.001) between the groups of animal of different sizes. For diameter >1.5 mm, mean anastomosis time decreased from 19.6±1.9 min to 12.6±0.7 min between the first and last week of training. For training involving smaller diameters, the results showed a reduction in execution time of 36.1% (P<0.01) (diameter between 1.0 mm and 1.5 mm) and 40.6% (P<0.01) (diameter <1.0 mm) between the first and last weeks. The study demonstrates an improvement in the dexterity and speed of nodes' execution. CONCLUSION The earthworm appears to be a reliable experimental model for microsurgical training of discrepancy microanastomoses. Its numerous advantages, as discussed in the present report, show that this model of training will significantly grow and develop in the near future.

The Initial Slope of the Variogram, Foundation of the Trabecular Bone Score, Is Not or Poorly Associated With Vertebral Strength

Trabecular bone score (TBS) rests on the textural analysis of DXA to reflect the decay in trabecular structure characterising osteoporosis. Yet, its discriminative power in fracture studies remains incomprehensible as prior biomechanical tests found no correlation with vertebral strength. To verify this result possibly due to an unrealistic set-up and to cover a wide range of loading scenarios, the data from three previous biomechanical studies using different experimental settings was used. They involved the compressive failure of 62 human lumbar vertebrae loaded 1) via intervertebral discs to mimic the in vivo situation (“full vertebra”), 2) via the classical endplate embedding (“vertebral body”) or 3) via a ball joint to induce anterior wedge failure (“vertebral section”). HR-pQCT scans acquired prior testing were used to simulate anterior-posterior DXA from which areal bone mineral density (aBMD) and the initial slope of the variogram (ISV), the early definition of TBS, were evaluated. Finally, the relation of aBMD and ISV with failure load (Fexp) and apparent failure stress (σexp) was assessed and their relative contribution to a multi-linear model was quantified via ANOVA. We found that, unlike aBMD, ISV did not significantly correlate with Fexp and σexp, except for the “vertebral body” case (r2 = 0.396, p = 0.028). Aside from the “vertebra section” set-up where it explained only 6.4% of σexp (p = 0.037), it brought no significant improvement to aBMD. These results indicate that ISV, a replica of TBS, is a poor surrogate for vertebral strength no matter the testing set-up, which supports the prior observations and raises a fortiori the question of the deterministic factors underlying the statistical relationship between TBS and vertebral fracture risk.

Serum Testosterone Levels and Symptom-Based Depression Subtypes in Men

The main objective of this preliminary study was to further clarify the association between testosterone (T) levels and depression by investigating symptom-based depression subtypes in a sample of 64 men. The data were taken from the ZInEP epidemiology survey. Gonadal hormones of a melancholic (n  = 25) and an atypical (n = 14) depression subtype, derived from latent class analysis, were compared with those of healthy controls (n  = 18). Serum T was assayed using an enzyme-linked immunosorbent assay procedure. Analysis of variance, analysis of covariance, non-parametrical tests, and generalized linear regression models were performed to examine group differences. The atypical depressive subtype showed significantly lower T levels compared with the melancholic depressives. While accumulative evidence indicates that, beyond psychosocial characteristics, the melancholic and atypical depressive subtypes are also distinguishable by biological correlates, the current study expanded this knowledge to include gonadal hormones. Further longitudinal research is warranted to disclose causality by linking the multiple processes in pathogenesis of depression.

CYP17A1 Enzyme Activity Is Linked to Ambulatory Blood Pressure in a Family-Based Population Study

BACKGROUND Genome-wide association studies have linked CYP17A1 coding for the steroid hormone synthesizing enzyme 17α-hydroxylase (CYP17A1) to blood pressure (BP). We hypothesized that the genetic signal may translate into a correlation of ambulatory BP (ABP) with apparent CYP17A1 activity in a family-based population study and estimated the heritability of CYP17A1 activity. METHODS In the Swiss Kidney Project on Genes in Hypertension, day and night urinary excretions of steroid hormone metabolites were measured in 518 participants (220 men, 298 women), randomly selected from the general population. CYP17A1 activity was assessed by 2 ratios of urinary steroid metabolites: one estimating the combined 17α-hydroxylase/17,20-lyase activity (ratio 1) and the other predominantly 17α-hydroxylase activity (ratio 2). A mixed linear model was used to investigate the association of ABP with log-transformed CYP17A1 activities exploring effect modification by urinary sodium excretion. RESULTS Daytime ABP was positively associated with ratio 1 under conditions of high, but not low urinary sodium excretion (P interaction <0.05). Ratio 2 was not associated with ABP. Heritability estimates (SE) for day and night CYP17A1 activities were 0.39 (0.10) and 0.40 (0.09) for ratio 1, and 0.71 (0.09) and 0.55 (0.09) for ratio 2 (P values <0.001). CYP17A1 activities, assessed with ratio 1, were lower in older participants. CONCLUSIONS Low apparent CYP17A1 activity (assessed with ratio 1) is associated with elevated daytime ABP when salt intake is high. CYP17A1 activity is heritable and diminished in the elderly. These observations highlight the modifying effect of salt intake on the association of CYP17A1 with BP.

The use of electromyography interference pattern analysis to determine muscle force of the deep digital flexor muscle in healthy and laminitic horses

BACKGROUND: In equine laminitis, the deep digital flexor muscle (DDFM) appears to have increased muscle force, but evidence-based confirmation is lacking. OBJECTIVES: The purpose of this study was to test if the DDFM of laminitic equines has an increased muscle force detectable by needle electromyography interference pattern analysis (IPA). ANIMALS AND METHODS: The control group included six Royal Dutch Sport horses, three Shetland ponies and one Welsh pony [10 healthy, sound adults weighing 411 ± 217 kg (mean ± SD) and aged 10 ± 5 years]. The laminitic group included three Royal Dutch Sport horses, one Friesian, one Haflinger, one Icelandic horse, one Welsh pony, one miniature Appaloosa and six Shetland ponies (14 adults, weight 310 ± 178 kg, aged 13 ± 6 years) with acute/chronic laminitis. The electromyography IPA measurements included firing rate, turns/second (T), amplitude/turn (M) and M/T ratio. Statistical analysis used a general linear model with outcomes transformed to geometric means. RESULTS: The firing rate of the total laminitic group was higher than the total control group. This difference was smaller for the ponies compared to the horses; in the horses, the geometric mean difference of the laminitic group was 1.73 [geometric 95% confidence interval (CI) 1.29-2.32], and in the ponies this value was 1.09 (geometric 95% CI 0.82-1.45). CONCLUSION AND CLINICAL RELEVANCE: In human medicine, an increased firing rate is characteristic of increased muscle force. Thus, the increased firing rate of the DDFM in the context of laminitis suggests an elevated muscle force. However, this seems to be only a partial effect as in this study, the unchanged turns/second and amplitude/turn failed to prove the recruitment of larger motor units with larger amplitude motor unit potentials in laminitic equids.

Age-specific and sex-specific weight gain norms to monitor antiretroviral therapy in children in low-income and middle-income countries.

BACKGROUND Viral load and CD4% are often not available in resource-limited settings for monitoring children's responses to antiretroviral therapy (ART). We aimed to construct normative curves for weight gain at 6, 12, 18, and 24 months following initiation of ART in children, and to assess the association between poor weight gain and subsequent responses to ART. DESIGN Analysis of data from HIV-infected children younger than 10 years old from African and Asian clinics participating in the International epidemiologic Databases to Evaluate AIDS. METHODS The generalized additive model for location, scale, and shape was used to construct normative percentile curves for weight gain at 6, 12, 18, and 24 months following ART initiation. Cox proportional models were used to assess the association between lower percentiles (< 50th) of weight gain distribution at the different time points and subsequent death, virological suppression, and virological failure. RESULTS Among 7173 children from five regions of the world, 45% were underweight at baseline. Weight gain below the 50th percentile at 6, 12, 18, and 24 months of ART was associated with increased risk of death, independent of baseline characteristics. Poor weight gain was not associated with increased hazards of virological suppression or virological failure. CONCLUSION Monitoring weight gain on ART using age-specific and sex-specific normative curves specifically developed for HIV-infected children on ART is a simple, rapid, sustainable tool that can aid in the identification of children who are at increased risk of death in the first year of ART.

A legacy of human-induced ecosystem changes: spatial processes drive the taxonomic and functional diversities of testate amoebae in Sphagnum peatlands of the Galápagos

Aim Our aims were to compare the composition of testate amoeba (TA) communities from Santa Cruz Island, Galápagos Archipelago, which are likely in existence only as a result of anthropogenic habitat transformation, with similar naturally occurring communities from northern and southern continental peatlands. Additionally, we aimed at assessing the importance of niche-based and dispersal-based processes in determining community composition and taxonomic and functional diversity. Location The humid highlands of the central island of Santa Cruz, Galápagos Archipelago. Methods We survey the alpha, beta and gamma taxonomic and functional diversities of TA, and the changes in functional traits along a gradient of wet to dry habitats. We compare the TA community composition, abundance and frequency recorded in the insular peatlands with that recorded in continental peatlands of Northern and Southern Hemispheres. We use generalized linear models to determine how environmental conditions influence taxonomic and functional diversity as well as the mean values of functional traits within communities. We finally apply variance partitioning to assess the relative importance of niche- and dispersal-based processes in determining community composition. Results TA communities in Santa Cruz Island were different from their Northern Hemisphere and South American counterparts with most genera considered as characteristic for Northern Hemisphere and South American Sphagnum peatlands missing or very rare in the Galápagos. Functional traits were most correlated with elevation and site topography and alpha functional diversity to the type of material sampled and site topography. Community composition was more strongly correlated with spatial variables than with environmental ones. Main conclusions TA communities of the Sphagnum peatlands of Santa Cruz Island and the mechanisms shaping these communities contrast with Northern Hemisphere and South American peatlands. Soil moisture was not a strong predictor of community composition most likely because rainfall and clouds provide sufficient moisture. Dispersal limitation was more important than environmental filtering because of the isolation of the insular peatlands from continental ones and the young ecological history of these ecosystems.

Disentangling mental imagery and perceptual expectation

Mental imagery and perception are thought to rely on similar neural circuits, and many recent behavioral studies have attempted to demonstrate interactions between actual physical stimulation and sensory imagery in the corresponding sensory modality. However, there has been a lack of theoretical understanding of the nature of these interactions, and both interferential and facilitatory effects have been found. Facilitatory effects appear strikingly similar to those that arise due to experimental manipulations of expectation. Using a self-motion discrimination task, we try to disentangle the effects of mental imagery from those of expectation by using a hierarchical drift diffusion model to investigate both choice data and response times. Manipulations of expectation are reasonably well understood in terms of their selective influence on parameters of the drift diffusion model, and in this study, we make the first attempt to similarly characterize the effects of mental imagery. We investigate mental imagery within the computational framework of control theory and state estimation. • Mental imagery and perception are thought to rely on similar neural circuits; however, on more theoretical grounds, imagery seems to be closely related to the output of forward models (sensory predictions). • We reanalyzed data from a study of imagined self-motion. • Bayesian modeling of response times may allow us to disentangle the effects of mental imagery on behavior from other cognitive (top-down) effects, such as expectation.

The nasal microbiota in infants with cystic fibrosis in the first year of life: a prospective cohort study.

BACKGROUND Respiratory tract infections and subsequent airway inflammation occur early in the life of infants with cystic fibrosis. However, detailed information about the microbial composition of the respiratory tract in infants with this disorder is scarce. We aimed to undertake longitudinal in-depth characterisation of the upper respiratory tract microbiota in infants with cystic fibrosis during the first year of life. METHODS We did this prospective cohort study at seven cystic fibrosis centres in Switzerland. Between Feb 1, 2011, and May 31, 2014, we enrolled 30 infants with a diagnosis of cystic fibrosis. Microbiota characterisation was done with 16S rRNA gene pyrosequencing and oligotyping of nasal swabs collected every 2 weeks from the infants with cystic fibrosis. We compared these data with data for an age-matched cohort of 47 healthy infants. We additionally investigated the effect of antibiotic treatment on the microbiota of infants with cystic fibrosis. Statistical methods included regression analyses with a multivariable multilevel linear model with random effects to correct for clustering on the individual level. FINDINGS We analysed 461 nasal swabs taken from the infants with cystic fibrosis; the cohort of healthy infants comprised 872 samples. The microbiota of infants with cystic fibrosis differed compositionally from that of healthy infants (p=0·001). This difference was also found in exclusively antibiotic-naive samples (p=0·001). The disordering was mainly, but not solely, due to an overall increase in the mean relative abundance of Staphylococcaceae in infants with cystic fibrosis compared with healthy infants (multivariable linear regression model stratified by age and adjusted for season; second month: coefficient 16·2 [95% CI 0·6-31·9]; p=0·04; third month: 17·9 [3·3-32·5]; p=0·02; fourth month: 21·1 [7·8-34·3]; p=0·002). Oligotyping analysis enabled differentiation between Staphylococcus aureus and coagulase-negative Staphylococci. Whereas the analysis showed a decrease in S aureus at and after antibiotic treatment, coagulase-negative Staphylococci increased. INTERPRETATION Our study describes compositional differences in the microbiota of infants with cystic fibrosis compared with healthy controls, and disordering of the microbiota on antibiotic administration. Besides S aureus, coagulase-negative Staphylococci also contributed to the disordering identified in these infants. These findings are clinically important in view of the crucial role that bacterial pathogens have in the disease progression of cystic fibrosis in early life. Our findings could be used to inform future studies of the effect of antibiotic treatment on the microbiota in infants with cystic fibrosis, and could assist in the prevention of early disease progression in infants with this disorder. FUNDING Swiss National Science Foundation, Fondation Botnar, the Swiss Society for Cystic Fibrosis, and the Swiss Lung Association Bern.

A Single-blinded, Randomized Clinical Trial of How to Implement an Evidence-based Treatment for Generalized Anxiety Disorder [IMPLEMENT] — Effects of Three Different Strategies of Implementation

BACKGROUND: Despite long-standing calls to disseminate evidence-based treatments for generalized anxiety (GAD), modest progress has been made in the study of how such treatments should be implemented. The primary objective of this study was to test three competing strategies on how to implement a cognitive behavioral treatment (CBT) for out-patients with GAD (i.e., comparison of one compensation vs. two capitalization models). METHODS: For our three-arm, single-blinded, randomized controlled trial (implementation of CBT for GAD [IMPLEMENT]), we recruited adults with GAD using advertisements in high-circulation newspapers to participate in a 14-session cognitive behavioral treatment (Mastery of your Anxiety and Worry, MAW-packet). We randomly assigned eligible patients using a full randomization procedure (1:1:1) to three different conditions of implementation: adherence priming (compensation model), which had a systematized focus on patients' individual GAD symptoms and how to compensate for these symptoms within the MAW-packet, and resource priming and supportive resource priming (capitalization model), which had systematized focuses on patients' strengths and abilities and how these strengths can be capitalized within the same packet. In the intention-to-treat population an outcome composite of primary and secondary symptoms-related self-report questionnaires was analyzed based on a hierarchical linear growth model from intake to 6-month follow-up assessment. This trial is registered at ClinicalTrials.gov (identifier: NCT02039193) and is closed to new participants. FINDINGS: From June 2012 to Nov. 2014, from 411 participants that were screened, 57 eligible participants were recruited and randomly assigned to three conditions. Forty-nine patients (86%) provided outcome data at post-assessment (14% dropout rate). All three conditions showed a highly significant reduction of symptoms over time. However, compared with the adherence priming condition, both resource priming conditions indicated faster symptom reduction. The observer ratings of a sub-sample of recorded videos (n = 100) showed that the therapists in the resource priming conditions conducted more strength-oriented interventions in comparison with the adherence priming condition. No patients died or attempted suicide. INTERPRETATION: To our knowledge, this is the first trial that focuses on capitalization and compensation models during the implementation of one prescriptive treatment packet for GAD. We have shown that GAD related symptoms were significantly faster reduced by the resource priming conditions, although the limitations of our study included a well-educated population. If replicated, our results suggest that therapists who implement a mental health treatment for GAD might profit from a systematized focus on capitalization models. FUNDING: Swiss Science National Foundation (SNSF-Nr. PZ00P1_136937/1) awarded to CF. KEYWORDS: Cognitive behavioral therapy; Evidence-based treatment; Implementation strategies; Randomized controlled trial