Twin study reveals non-heritable immune perturbations in multiple sclerosis.

Multiple sclerosis (MS) is a chronic inflammatory disorder of the central nervous system underpinned by partially understood genetic risk factors and environmental triggers and their undefined interactions1,2. Here we investigated the peripheral immune signatures of 61 monozygotic twin pairs discordant for MS to dissect the influence of genetic predisposition and environmental factors. Using complementary multimodal high-throughput and high-dimensional single-cell technologies in conjunction with data-driven computational tools, we identified an inflammatory shift in a monocyte cluster of twins with MS, coupled with the emergence of a population of IL-2 hyper-responsive transitional naive helper T cells as MS-related immune alterations. By integrating data on the immune profiles of healthy monozygotic and dizygotic twin pairs, we estimated the variance in CD25 expression by helper T cells displaying a naive phenotype to be largely driven by genetic and shared early environmental influences. Nonetheless, the expanding helper T cells of twins with MS, which were also elevated in non-twin patients with MS, emerged independent of the individual genetic makeup. These cells expressed central nervous system-homing receptors, exhibited a dysregulated CD25-IL-2 axis, and their proliferative capacity positively correlated with MS severity. Together, our matched-pair analysis of the extended twin approach allowed us to discern genetically and environmentally determined features of an MS-associated immune signature.

Ecosystem service change caused by climatological and non-climatological drivers: a Swiss case study.

Understanding the drivers of ecosystem change and their effects on ecosystem services are essential for management decisions and verification of progress towards national and international sustainability policies (e.g., Aichi Biodiversity Targets, Sustainable Development Goals). We aim to disentangle spatially the effect of climatological and non-climatological drivers on ecosystem service supply and trends. Therefore, we explored time series of three ecosystem services in Switzerland between 2004 and 2014: carbon dioxide regulation, soil erosion prevention, and air quality regulation. We applied additive models to describe the spatial variation attributed to climatological (i.e., temperature, precipitation and relative sunshine duration) and non-climatological drivers (i.e., random effects representing other spatially structured processes) that may affect ecosystem service change. Obtained results indicated strong influences of climatological drivers on ecosystem service trends in Switzerland. We identified equal contributions of all three climatological drivers on trends of carbon dioxide regulation and soil erosion prevention, while air quality regulation was more strongly influenced by temperature. Additionally, our results showed that climatological and non-climatological drivers affected ecosystem services both negatively and positively, depending on the regions (in particular lower and higher altitudinal areas), drivers, and services assessed. Our findings highlight stronger effects of climatological compared to non-climatological drivers on ecosystem service change in Switzerland. Furthermore, drivers of ecosystem change display a spatial heterogeneity in their influence on ecosystem service trends. We propose an approach building on an additive model to disentangle the effect of climatological and non-climatological drivers on ecosystem service trends. Such analyses should be extended in the future to ecosystem service flow and demand to complete ecosystem service assessments and to demonstrate and communicate more clearly the benefits of ecosystem services for human well-being.

Additive Bayesian networks for antimicrobial resistance and potential risk factors in non-typhoidal Salmonella isolates from layer hens in Uganda.

BACKGROUND: Multi-drug resistant bacteria are seen increasingly and there are gaps in our understanding of the complexity of antimicrobial resistance, partially due to a lack of appropriate statistical tools. This hampers efficient treatment, precludes determining appropriate intervention points and renders prevention very difficult. METHODS: We re-analysed data from a previous study using additive Bayesian networks. The data contained information on resistances against seven antimicrobials and seven potential risk factors from 86 non-typhoidal Salmonella isolates from laying hens in 46 farms in Uganda. RESULTS: The final graph contained 22 links between risk factors and antimicrobial resistances. Solely ampicillin resistance was linked to the vaccinating person and disposal of dead birds. Systematic associations between ampicillin and sulfamethoxazole/trimethoprim and chloramphenicol, which was also linked to sulfamethoxazole/trimethoprim were detected. Sulfamethoxazole/trimethoprim was also directly linked to ciprofloxacin and trimethoprim. Trimethoprim was linked to sulfonamide and ciprofloxacin, which was also linked to sulfonamide. Tetracycline was solely linked to ciprofloxacin. CONCLUSIONS: Although the results needs to be interpreted with caution due to a small data set, additive Bayesian network analysis allowed a description of a number of associations between the risk factors and antimicrobial resistances investigated.

Changes in temperature sensitivity of spring phenology with recent climate warming in Switzerland are related to shifts of the preseason.

The spring phenology of plants in temperate regions strongly responds to spring temperatures. Climate warming has caused substantial phenological advances in the past, but trends to be expected in the future are uncertain. A simple indicator is temperature sensitivity, the phenological advance statistically associated with a 1°C warmer mean temperature during the "preseason", defined as the most temperature-sensitive period preceding the phenological event. Recent analyses of phenological records have shown a decline in temperature sensitivity of leaf unfolding, but underlying mechanisms were not clear. Here, we propose that climate warming can reduce temperature sensitivity simply by reducing the length of the preseason due to faster bud development during this time period, unless the entire preseason shifts forward so that its temperature does not change. We derive these predictions theoretically from the widely used "thermal time model" for bud development and test them using data for 19 phenological events recorded in 1970-2012 at 108 stations spanning a 1600 m altitudinal range in Switzerland. We consider how temperature sensitivity, preseason start, preseason length and preseason temperature change (i) with altitude, (ii) between the periods 1970-1987 and 1995-2012, which differed mainly in spring temperatures, and (iii) between two non-consecutive sets of 18 years that differed mainly in winter temperatures. On average, temperature sensitivity increased with altitude (colder climate) and was reduced in years with warmer springs, but not in years with warmer winters. These trends also varied among species. Decreasing temperature sensitivity in warmer springs was associated with a limited forward shift of preseason start, higher temperatures during the preseason and reduced preseason length, but not with reduced winter chilling. Our results imply that declining temperature sensitivity can result directly from spring warming and does not necessarily indicate altered physiological responses or stronger constraints such as reduced winter chilling.

Relative contribution of biological variation and technical variables to zone diameter variations of disc diffusion susceptibility testing.

OBJECTIVES: Disc diffusion is still largely based on manual procedures. Technical variations originate from inoculum preparation, variations in materials, individual operator plate streaking and reading accuracy. Resulting measurement imprecision contributes to categorization errors. Biological variation resembles the natural fluctuation of a measured parameter such as antibiotic susceptibility around a mean value. It is deemed to originate from factors such as genetic background or metabolic state. This study analysed the relative contribution of different technical and biological factors to total disc diffusion variation. METHODS: For calculation of relative error factor contribution to disc diffusion variability, five experiments were designed keeping different combinations of error factors constant. A mathematical model was developed to analyse the individual error factor contribution to disc diffusion variation for each of the tested drug-species combinations. RESULTS: The contribution of biological variation to total diameter variance ranged from 10.4% to 98.8% for different drug-species combinations. Highest biological variation was found for Enterococcus faecalis WT and vancomycin (98.8%) and for penicillinase-producing Staphylococcus aureus and penicillin G (96.0%). Average imprecision of automated zone reading revealed that 1.4%-5.3% of total imprecision was due to technical variation, while materials, i.e. antibiotic discs and agar plates, contributed between 2.6% and 3.9%. Inoculum preparation and manual plate streaking contributed 6.8%-24.8% and 6.6%-24.3%, respectively, to total imprecision. CONCLUSIONS: This study illustrates the relative contributions of technical factors that account for a significant part of total variance in disc diffusion. The highest relative contribution originated from the operator, i.e. manual inoculum preparation and plate streaking. Further standardization of inoculum preparation and plate streaking by automation could potentially increase the precision of disc diffusion and improve the correlation of susceptibility reports with clinical outcome.

Protocol for a prospective, randomized study on neurophysiological assessment of lower urinary tract function in a healthy cohort.

BACKGROUND: Lower urinary tract symptoms are highly prevalent and a large proportion of these symptoms are known to be associated with a dysfunction of the afferent pathways. Diagnostic tools for an objective and reproducible assessment of afferent nerve function of the lower urinary tract are missing. Previous studies showed first feasibility results of sensory evoked potential recordings following electrical stimulation of the lower urinary tract in healthy subjects and patients. Nevertheless, a refinement of the methodology is necessary. METHODS: This study is a prospective, randomized trial conducted at Balgrist University Hospital, Zürich, Switzerland. Ninety healthy subjects (forty females and fifty males) without lower urinary tract symptoms are planned to be included in the study. All subjects will undergo a screening visit (including standardized questionnaires, 3-day bladder diary, urinalysis, medical history taking, vital signs, physical examination, neuro-urological examination) followed by two measurement visits separated by an interval of 3 to 4 weeks. Electrical stimulations (0.5Hz-5Hz, bipolar, square wave, pulse width 1 ms) will be applied using a custom-made transurethral catheter at different locations of the lower urinary tract including bladder dome, trigone, proximal urethra, membranous urethra and distal urethra. Every subject will be randomly stimulated at one specific site of the lower urinary tract. Sensory evoked potentials (SEP) will be recorded using a 64-channel EEG cap. For an SEP segmental work-up we will place additional electrodes on the scalp (Cpz) and above the spine (C2 and L1). Visit two and three will be conducted identically for reliability assessment. DISCUSSION: The measurement of lower urinary tract SEPs elicited by electrical stimulation at different locations of the lower urinary tract has the potential to serve as a neurophysiological biomarker for lower urinary tract afferent nerve function in patients with lower urinary tract symptoms or disorders. For implementation of such a diagnostic tool into clinical practice, an optimized setup with efficient and reliable measurements and data acquisition is crucial. In addition, normative data from a larger cohort of healthy subjects would provide information on variability, potential confounding factors and cut-off values for investigations in patients with lower urinary tract dysfunction/symptoms. TRIAL REGISTRATION: Clinicaltrials.gov; Identifier: NCT02272309 .

Standardization of Operator-Dependent Variables Affecting Precision and Accuracy of the Disk Diffusion Method for Antibiotic Susceptibility Testing.

Parameters like zone reading, inoculum density, and plate streaking influence the precision and accuracy of disk diffusion antibiotic susceptibility testing (AST). While improved reading precision has been demonstrated using automated imaging systems, standardization of the inoculum and of plate streaking have not been systematically investigated yet. This study analyzed whether photometrically controlled inoculum preparation and/or automated inoculation could further improve the standardization of disk diffusion. Suspensions of Escherichia coli ATCC 25922 and Staphylococcus aureus ATCC 29213 of 0.5 McFarland standard were prepared by 10 operators using both visual comparison to turbidity standards and a Densichek photometer (bioMérieux), and the resulting CFU counts were determined. Furthermore, eight experienced operators each inoculated 10 Mueller-Hinton agar plates using a single 0.5 McFarland standard bacterial suspension of E. coli ATCC 25922 using regular cotton swabs, dry flocked swabs (Copan, Brescia, Italy), or an automated streaking device (BD-Kiestra, Drachten, Netherlands). The mean CFU counts obtained from 0.5 McFarland standard E. coli ATCC 25922 suspensions were significantly different for suspensions prepared by eye and by Densichek (P < 0.001). Preparation by eye resulted in counts that were closer to the CLSI/EUCAST target of 10(8) CFU/ml than those resulting from Densichek preparation. No significant differences in the standard deviations of the CFU counts were observed. The interoperator differences in standard deviations when dry flocked swabs were used decreased significantly compared to the differences when regular cotton swabs were used, whereas the mean of the standard deviations of all operators together was not significantly altered. In contrast, automated streaking significantly reduced both interoperator differences, i.e., the individual standard deviations, compared to the standard deviations for the manual method, and the mean of the standard deviations of all operators together, i.e., total methodological variation.

Prediction of coccidiosis prevalence in extensive backyard chickens in countries and regions of the Horn of Africa.

Coccidiosis is one of the leading morbidity causes in chickens, causing a reduction of body weight and egg production. Backyard chickens are at risk of developing clinical and subclinical coccidiosis due to outdoor housing and scavenging behaviour, jeopardizing food security in households. The objectives of this study were to estimate clinical prevalence of coccidiosis at country and regional levels in the Horn of Africa in extensive backyard chickens. A binomial random effects model was developed to impute prevalence of coccidiosis. Previously gathered prevalence data (n = 40) in backyard chickens was used to define the model. Precipitation (OR: 1.09 (95% CI: 1.05-1.13) and the presence of seasonal rainfall (OR: 1.85, 95% CI: 1.27-2.70) significantly increase prevalence. Results showed an overall prevalence of coccidiosis in the Horn of Africa of 0.21 (95% CI: 0.15-0.29). Ethiopia, the Republic of South Sudan and Kenya showed the highest prevalence and Djibouti the lowest. Significant differences between Djibouti and the countries with highest prevalence were found. However, no evidence of a significant difference between the rest of the countries. Kenya and Ethiopia showed larger prevalence differences between regions. Results could assist with the targeting of testing for coccidiosis, the observation for clinical disease of chickens living in specific regions and as a baseline for the evaluation of future control measures.

Multiorder correction algorithms to remove image distortions from mass spectrometry imaging data sets.

Time-of-flight secondary ion mass spectrometry imaging is a rapidly evolving technology. Its main application is the study of the distribution of small molecules on biological tissues. The sequential image acquisition process remains susceptible to measurement distortions that can render imaging data less analytically useful. Most of these artifacts show a repetitive nature from tile to tile. Here we statistically describe these distortions and derive two different algorithms to correct them. Both a generalized linear model approach and the linear discriminant analysis approach are able to increase image quality for negative and positive ion mode data sets. Additionally, performing simulation studies with repetitive and nonrepetitive tiling error we show that both algorithms are only removing repetitive distortions. It is further shown that the spectral component of the data set is not altered by the use of these correction methods. Both algorithms presented in this work greatly increase the image quality and improve the analytical usefulness of distorted images dramatically.

Simultaneously estimating food web connectance and structure with uncertainty.

Food web models explain and predict the trophic interactions in a food web, and they can infer missing interactions among the organisms. The allometric diet breadth model (ADBM) is a food web model based on the foraging theory. In the ADBM, the foraging parameters are allometrically scaled to body sizes of predators and prey. In Petchey et al. (Proceedings of the National Academy of Sciences, 2008; 105: 4191), the parameterization of the ADBM had two limitations: (a) the model parameters were point estimates and (b) food web connectance was not estimated.The novelty of our current approach is: (a) We consider multiple predictions from the ADBM by parameterizing it with approximate Bayesian computation, to estimate parameter distributions and not point estimates. (b) Connectance emerges from the parameterization, by measuring model fit using the true skill statistic, which takes into account prediction of both the presences and absences of links.We fit the ADBM using approximate Bayesian computation to 12 observed food webs from a wide variety of ecosystems. Estimated connectance was consistently greater than previously found. In some of the food webs, considerable variation in estimated parameter distributions occurred and resulted in considerable variation (i.e., uncertainty) in predicted food web structure.These results lend weight to the possibility that the observed food web data is missing some trophic links that do actually occur. It also seems likely that the ADBM likely predicts some links that do not exist. The latter could be addressed by accounting in the ADBM for additional traits other than body size. Further work could also address the significance of uncertainty in parameter estimates for predicted food web responses to environmental change.

Bayesian network analysis reveals the interplay of intracranial aneurysm rupture risk factors.

Clinical decision making regarding the treatment of unruptured intracranial aneurysms (IA) benefits from a better understanding of the interplay of IA rupture risk factors. Probabilistic graphical models can capture and graphically display potentially causal relationships in a mechanistic model. In this study, Bayesian networks (BN) were used to estimate IA rupture risk factors influences. From 1248 IA patient records, a retrospective, single-cohort, patient-level data set with 9 phenotypic rupture risk factors (n=790 complete entries) was extracted. Prior knowledge together with score-based structure learning algorithms estimated rupture risk factor interactions. Two approaches, discrete and mixed-data additive BN, were implemented and compared. The corresponding graphs were learned using non-parametric bootstrapping and Markov chain Monte Carlo, respectively. The BN models were compared to standard descriptive and regression analysis methods. Correlation and regression analyses showed significant associations between IA rupture status and patient's sex, familial history of IA, age at IA diagnosis, IA location, IA size and IA multiplicity. BN models confirmed the findings from standard analysis methods. More precisely, they directly associated IA rupture with familial history of IA, IA size and IA location in a discrete framework. Additive model formulation, enabling mixed-data, found that IA rupture was directly influenced by patient age at diagnosis besides additional mutual influences of the risk factors. This study establishes a data-driven methodology for mechanistic disease modelling of IA rupture and shows the potential to direct clinical decision-making in IA treatment, allowing personalised prediction.

Discussion on Competition for Spatial Statistics for Large Datasets.

We discuss the experiences and results of the AppStatUZH team's participation in the comprehensive and unbiased comparison of different spatial approximations conducted in the Competition for Spatial Statistics for Large Datasets. In each of the different sub-competitions, we estimated parameters of the covariance model based on a likelihood function and predicted missing observations with simple kriging. We approximated the covariance model either with covariance tapering or a compactly supported Wendland covariance function.

Remotely sensed between-individual functional trait variation in a temperate forest.

Trait-based ecology holds the promise to explain how plant communities work, for example, how functional diversity may support community productivity. However, so far it has been difficult to combine field-based approaches assessing traits at the level of plant individuals with limited spatial coverage and approaches using remote sensing (RS) with complete spatial coverage but assessing traits at the level of vegetation pixels rather than individuals. By delineating all individual-tree crowns within a temperate forest site and then assigning RS-derived trait measures to these trees, we combine the two approaches, allowing us to use general linear models to estimate the influence of taxonomic or environmental variation on between- and within-species variation across contiguous space.We used airborne imaging spectroscopy and laser scanning to collect individual-tree RS data from a mixed conifer-angiosperm forest on a mountain slope extending over 5.5 ha and covering large environmental gradients in elevation as well as light and soil conditions. We derived three biochemical (leaf chlorophyll, carotenoids, and water content) and three architectural traits (plant area index, foliage-height diversity, and canopy height), which had previously been used to characterize plant function, from the RS data. We then quantified the contributions of taxonomic and environmental variation and their interaction to trait variation and partitioned the remaining within-species trait variation into smaller-scale spatial and residual variation. We also investigated the correlation between functional trait and phylogenetic distances at the between-species level. The forest consisted of 13 tree species of which eight occurred in sufficient abundance for quantitative analysis.On average, taxonomic variation between species accounted for more than 15% of trait variation in biochemical traits but only around 5% (still highly significant) in architectural traits. Biochemical trait distances among species also showed a stronger correlation with phylogenetic distances than did architectural trait distances. Light and soil conditions together with elevation explained slightly more variation than taxonomy across all traits, but in particular increased plant area index (light) and reduced canopy height (elevation). Except for foliage-height diversity, all traits were affected by significant interactions between taxonomic and environmental variation, the different responses of the eight species to the within-site environmental gradients potentially contributing to the coexistence of the eight abundant species.We conclude that with high-resolution RS data it is possible to delineate individual-tree crowns within a forest and thus assess functional traits derived from RS data at individual level. With this precondition fulfilled, it is then possible to apply tools commonly used in field-based trait ecology to partition trait variation among individuals into taxonomic and potentially even genetic variation, environmental variation, and interactions between the two. The method proposed here presents a promising way of assessing individual-based trait information with complete spatial coverage and thus allowing analysis of functional diversity at different scales. This information can help to better understand processes shaping community structure, productivity, and stability of forests.

Bayesian spatial modelling of terrestrial radiation in Switzerland.

The geographic variation of terrestrial radiation can be exploited in epidemiological studies of the health effects of protracted low-dose exposure. Various methods have been applied to derive maps of this variation. We aimed to construct a map of terrestrial radiation for Switzerland. We used airborne γ-spectrometry measurements to model the ambient dose rates from terrestrial radiation through a Bayesian mixed-effects model and conducted inference using Integrated Nested Laplace Approximation (INLA). We predicted higher levels of ambient dose rates in the alpine regions and Ticino compared with the western and northern parts of Switzerland. We provide a map that can be used for exposure assessment in epidemiological studies and as a baseline map for assessing potential contamination.

Association between environmental and climatic risk factors and the spatial distribution of cystic and alveolar echinococcosis in Kyrgyzstan.

BACKGROUND: Cystic and alveolar echinococcosis (CE and AE) are neglected tropical diseases caused by Echinococcus granulosus sensu lato and E. multilocularis, and are emerging zoonoses in Kyrgyzstan. In this country, the spatial distribution of CE and AE surgical incidence in 2014-2016 showed marked heterogeneity across communities, suggesting the presence of ecological determinants underlying CE and AE distributions. METHODOLOGY/PRINCIPAL FINDINGS: For this reason, in this study we assessed potential associations between community-level confirmed primary CE (no.=2359) or AE (no.=546) cases in 2014-2016 in Kyrgyzstan and environmental and climatic variables derived from satellite-remote sensing datasets using conditional autoregressive models. We also mapped CE and AE relative risk. The number of AE cases was negatively associated with 10-year lag mean annual temperature. Although this time lag should not be considered as an exact measurement but with associated uncertainty, it is consistent with the estimated 10-15-year latency following AE infection. No associations were detected for CE. We also identified several communities at risk for CE or AE where no disease cases were reported in the study period. CONCLUSIONS/SIGNIFICANCE: Our findings support the hypothesis that CE is linked to an anthropogenic cycle and is less affected by environmental risk factors compared to AE, which is believed to result from spillover from a wild life cycle. As CE was not affected by factors we investigated, hence control should not have a geographical focus. In contrast, AE risk areas identified in this study without reported AE cases should be targeted for active disease surveillance in humans. This active surveillance would confirm or exclude AE transmission which might not be reported with the present passive surveillance system. These areas should also be targeted for ecological investigations in the animal hosts.

rasterdiv-An Information Theory tailored R package for measuring ecosystem heterogeneity from space: To the origin and back.

Ecosystem heterogeneity has been widely recognized as a key ecological indicator of several ecological functions, diversity patterns and change, metapopulation dynamics, population connectivity or gene flow.In this paper, we present a new R package-rasterdiv-to calculate heterogeneity indices based on remotely sensed data. We also provide an ecological application at the landscape scale and demonstrate its power in revealing potentially hidden heterogeneity patterns.The rasterdiv package allows calculating multiple indices, robustly rooted in Information Theory, and based on reproducible open-source algorithms.

Generation and application of river network analogues for use in ecology and evolution.

Several key processes in freshwater ecology are governed by the connectivity inherent to dendritic river networks. These have extensively been analyzed from a geomorphological and hydrological viewpoint, yet structures classically used in ecological modeling have been poorly representative of the structure of real river basins, often failing to capture well-known scaling features of natural rivers. Pioneering work identified optimal channel networks (OCNs) as spanning trees reproducing all scaling features characteristic of natural stream networks worldwide. While OCNs have been used to create landscapes for studies on metapopulations, biodiversity, and epidemiology, their generation has not been generally accessible.Given the increasing interest in dendritic riverine networks by ecologists and evolutionary biologists, we here present a method to generate OCNs and, to facilitate its application, we provide the R-package OCNet. Owing to the stochastic process generating OCNs, multiple network replicas spanning the same surface can be built; this allows performing computational experiments whose results are irrespective of the particular shape of a single river network. The OCN construct also enables the generation of elevational gradients derived from the optimal network configuration, which can constitute three-dimensional landscapes for spatial studies in both terrestrial and freshwater realms. Moreover, the package provides functions that aggregate OCNs into an arbitrary number of nodes, calculate several descriptors of river networks, and draw relevant network features.We describe the main functionalities of the package and its integration with other R-packages commonly used in spatial ecology. Moreover, we exemplify the generation of OCNs and discuss an application to a metapopulation model for an invasive riverine species.In conclusion, OCNet provides a powerful tool to generate realistic river network analogues for various applications. It thereby allows the design of spatially realistic studies in increasingly impacted ecosystems and enhances our knowledge on spatial processes in freshwater ecology in general.

Epidemic cystic and alveolar echinococcosis in Kyrgyzstan: an analysis of national surveillance data.

BACKGROUND: Human cystic and alveolar echinococcosis are among the priority neglected zoonotic diseases for which WHO advocates control. The incidence of both cystic echinococcosis and alveolar echinococcosis has increased substantially in the past 30 years in Kyrgyzstan. Given the scarcity of adequate data on the local geographical variation of these focal diseases, we aimed to investigate within-country incidence and geographical variation of cystic echinococcosis and alveolar echinococcosis at a high spatial resolution in Kyrgyzstan. METHODS: We mapped all confirmed surgical cases of cystic echinococcosis and alveolar echinococcosis reported through the national echinococcosis surveillance system in Kyrgyzstan between Jan 1, 2014, and Dec 31, 2016, from nine regional databases. We then estimated crude surgical incidence, standardised incidence, and standardised incidence ratios (SIRs) of primary cases (ie, excluding relapses) based on age and sex at country, region, district, and local community levels. Finally, we tested the SIRs for global and local spatial autocorrelation to identify disease hotspots at the local community level. All incidence estimates were calculated per 100 000 population and averaged across the 3-year study period to obtain annual estimates. FINDINGS: The surveillance system reported 2359 primary surgical cases of cystic echinococcosis and 546 primary surgical cases of alveolar echinococcosis. Country-level crude surgical incidence was 13·1 per 100 000 population per year for cystic echinococcosis and 3·02 per 100 000 population per year for alveolar echinococcosis. At the local community level, we found annual crude surgical incidences up to 176 per 100 000 population in Sary-Kamysh (Jalal-Abad region) for cystic echinococcosis and 246 per 100 000 population in Uch-Dobo (Alay district, Osh region) for alveolar echinococcosis. Significant hotspots of cystic echinococcosis were found in four regions: Osh (five local communities in Uzgen district and four in Alay district), Naryn (three local communities in Jumgal district and one in Naryn district), Talas (three local communities in Talas district), and Chuy (one local community in Jayyl district). Significant alveolar echinococcosis hotspots were detected in the Osh region (11 communities in Alay district, including the local community of Sary Mogol, and one in Chong-Alay district) and in the Naryn region (five communities in Jumgal district and three in At-Bashy district), in the southwest and centre of the country. INTERPRETATION: Our analyses reveal remarkable within-country variation in the surgical incidence of cystic echinococcosis and alveolar echinococcosis in Kyrgyzstan. These high-resolution maps identify precise locations where interventions and epidemiological research should be targeted to reduce the burden of human cystic echinococcosis and alveolar echinococcosis. FUNDING: Swiss National Science Foundation.

Bayesian Network Modeling Applied to Feline Calicivirus Infection Among Cats in Switzerland.

Bayesian network (BN) modeling is a rich and flexible analytical framework capable of elucidating complex veterinary epidemiological data. It is a graphical modeling technique that enables the visual presentation of multi-dimensional results while retaining statistical rigor in population-level inference. Using previously published case study data about feline calicivirus (FCV) and other respiratory pathogens in cats in Switzerland, a full BN modeling analysis is presented. The analysis shows that reducing the group size and vaccinating animals are the two actionable factors directly associated with FCV status and are primary targets to control FCV infection. The presence of gingivostomatitis and Mycoplasma felis is also associated with FCV status, but signs of upper respiratory tract disease (URTD) are not. FCV data is particularly well-suited to a network modeling approach, as both multiple pathogens and multiple clinical signs per pathogen are involved, along with multiple potentially interrelated risk factors. BN modeling is a holistic approach-all variables of interest may be mutually interdependent-which may help to address issues, such as confounding and collinear factors, as well as to disentangle directly vs. indirectly related variables. We introduce the BN methodology as an alternative to the classical uni- and multivariable regression approaches commonly used for risk factor analyses. We advise and guide researchers about how to use BNs as an exploratory data tool and demonstrate the limitations and practical issues. We present a step-by-step case study using FCV data along with all code necessary to reproduce our analyses in the open-source R environment. We compare and contrast the findings of the current case study using BN modeling with previous results that used classical regression techniques, and we highlight new potential insights. Finally, we discuss advanced methods, such as Bayesian model averaging, a common way of accounting for model uncertainty in a Bayesian network context.

Progression and risk factors of pododermatitis in part-time group housed rabbit does in Switzerland.

In rabbits (Oryctolagus cuniculus L.), pododermatitis is a chronic multifactorial skin disease that appears mainly on the plantar surface of the hind legs. This presumably progressive disease can cause pain leading to poor welfare, yet the progression of this disease has not been thoroughly assessed on the level of individual animals. The aim of this longitudinal study thus was to investigate the possible risk factors and the progression of pododermatitis in group housed breeding does in Switzerland on litter and plastic slats. Three commercial rabbit farms with part-time group housing on litter and plastic slats were visited every four weeks throughout one year. During every visit, the same 201 adult female breeding rabbits (67 does per farm) were evaluated for the presence and severity of pododermatitis. Additionally, the does' age, parity, body weight, reproductive state, hybrid, claw length, cleanliness and moisture of the paws and the temperature and humidity inside the barns were recorded as potential risk factors. The risk factors were analysed through general linear models and additive Bayesian network (ABN) modelling using a directed acyclic graph (DAG) for visualising associations between potential risk factors. The progression of pododermatitis was analysed with a transition matrix. Relative humidity inside the barns, body weight, number of kindlings, age, and claw length were the most important risk factors, all being positively associated with pododermatitis. In contrast to expectations, the cleanliness of the left hind paw was negatively associated with the occurrence of pododermatitis, but the effect was small. In young does, the severity of pododermatitis quickly increased and in some rabbits proceeded to ulcerated spots. It was shown that 60.00%, 14.17% and 3.33% of ulcerated lesions recovered to a state without ulceration within 4, 8 or >12 weeks, respectively.