Reconstruction of patient-specific 3D bone surface from 2D calibrated fluoroscopic images and point distribution model

Reconstruction of patient-specific 3D bone surface from 2D calibrated fluoroscopic images and a point distribution model is discussed. We present a 2D/3D reconstruction scheme combining statistical extrapolation and regularized shape deformation with an iterative image-to-model correspondence establishing algorithm, and show its application to reconstruct the surface of proximal femur. The image-to-model correspondence is established using a non-rigid 2D point matching process, which iteratively uses a symmetric injective nearest-neighbor mapping operator and 2D thin-plate splines based deformation to find a fraction of best matched 2D point pairs between features detected from the fluoroscopic images and those extracted from the 3D model. The obtained 2D point pairs are then used to set up a set of 3D point pairs such that we turn a 2D/3D reconstruction problem to a 3D/3D one. We designed and conducted experiments on 11 cadaveric femurs to validate the present reconstruction scheme. An average mean reconstruction error of 1.2 mm was found when two fluoroscopic images were used for each bone. It decreased to 1.0 mm when three fluoroscopic images were used.

Accurate and robust reconstruction of proximal femur from sparse intraoperative data and dense point distribution model for surgical navigation

Constructing a 3D surface model from sparse-point data is a nontrivial task. Here, we report an accurate and robust approach for reconstructing a surface model of the proximal femur from sparse-point data and a dense-point distribution model (DPDM). The problem is formulated as a three-stage optimal estimation process. The first stage, affine registration, is to iteratively estimate a scale and a rigid transformation between the mean surface model of the DPDM and the sparse input points. The estimation results of the first stage are used to establish point correspondences for the second stage, statistical instantiation, which stably instantiates a surface model from the DPDM using a statistical approach. This surface model is then fed to the third stage, kernel-based deformation, which further refines the surface model. Handling outliers is achieved by consistently employing the least trimmed squares (LTS) approach with a roughly estimated outlier rate in all three stages. If an optimal value of the outlier rate is preferred, we propose a hypothesis testing procedure to automatically estimate it. We present here our validations using four experiments, which include 1 leave-one-out experiment, 2 experiment on evaluating the present approach for handling pathology, 3 experiment on evaluating the present approach for handling outliers, and 4 experiment on reconstructing surface models of seven dry cadaver femurs using clinically relevant data without noise and with noise added. Our validation results demonstrate the robust performance of the present approach in handling outliers, pathology, and noise. An average 95-percentile error of 1.7-2.3 mm was found when the present approach was used to reconstruct surface models of the cadaver femurs from sparse-point data with noise added.

Candida species distribution and antifungal susceptibility testing according to European Committee on Antimicrobial Susceptibility Testing and new vs. old Clinical and Laboratory Standards Institute clinical breakpoints: a 6-year prospective candidaemia survey from the fungal infection network of Switzerland

We analyzed the species distribution of Candida blood isolates (CBIs), prospectively collected between 2004 and 2009 within FUNGINOS, and compared their antifungal susceptibility according to clinical breakpoints defined by the European Committee on Antimicrobial Susceptibility Testing (EUCAST) in 2013, and the Clinical and Laboratory Standards Institute (CLSI) in 2008 (old CLSI breakpoints) and 2012 (new CLSI breakpoints). CBIs were tested for susceptiblity to fluconazole, voriconazole and caspofungin by microtitre broth dilution (Sensititre® YeastOne™ test panel). Of 1090 CBIs, 675 (61.9%) were C. albicans, 191 (17.5%) C. glabrata, 64 (5.9%) C. tropicalis, 59 (5.4%) C. parapsilosis, 33 (3%) C. dubliniensis, 22 (2%) C. krusei and 46 (4.2%) rare Candida species. Independently of the breakpoints applied, C. albicans was almost uniformly (>98%) susceptible to all three antifungal agents. In contrast, the proportions of fluconazole- and voriconazole-susceptible C. tropicalis and F-susceptible C. parapsilosis were lower according to EUCAST/new CLSI breakpoints than to the old CLSI breakpoints. For caspofungin, non-susceptibility occurred mainly in C. krusei (63.3%) and C. glabrata (9.4%). Nine isolates (five C. tropicalis, three C. albicans and one C. parapsilosis) were cross-resistant to azoles according to EUCAST breakpoints, compared with three isolates (two C. albicans and one C. tropicalis) according to new and two (2 C. albicans) according to old CLSI breakpoints. Four species (C. albicans, C. glabrata, C. tropicalis and C. parapsilosis) represented >90% of all CBIs. In vitro resistance to fluconazole, voriconazole and caspofungin was rare among C. albicans, but an increase of non-susceptibile isolates was observed among C. tropicalis/C. parapsilosis for the azoles and C. glabrata/C. krusei for caspofungin according to EUCAST and new CLSI breakpoints compared with old CLSI breakpoints.

Bacteria causing bacteremia in pediatric cancer patients presenting with febrile neutropenia – species distribution and susceptibility patterns

PURPOSE Infections are a major cause of morbidity and mortality in pediatric cancer patients. The aim of this study was to establish the microbiological spectrum and the susceptibility patterns of bacteremia-causing bacteria in pediatric cancer patients with febrile neutropenia in relation to the use of prophylactic and empirical antibiotics. METHODS We analyzed positive blood cultures of pediatric cancer patients presenting with febrile neutropenia between 2004 and 2011 in Groningen and Amsterdam (the Netherlands) and in Bern (Switzerland), using different antibiotic prophylactic and empirical regimens. RESULTS A total of 156 patients with 202 bacteremias, due to 248 bacteria species, were enrolled. The majority (73%) of bacteremias were caused by Gram-positive bacteria. Gram-negative bacteria, especially Pseudomonas aeruginosa, were observed significantly more often in Bern, where no fluoroquinolone prophylaxis was used. Ciprofloxacin-resistant bacteria were cultured more often from patients who did receive ciprofloxacin prophylaxis, compared to the patients who did not (57 versus 11%, p = 0.044). CONCLUSIONS Gram-positive bacteria predominated in this study. We showed that the use of prophylactic antibiotics in pediatric cancer patients was associated with increased resistance rates, which needs further study. The strategy for empiric antimicrobial therapy for febrile neutropenia should be adapted to local antibiotic resistance patterns.

Climate refugia: joint inference from fossil records, species distribution models and phylogeography

Climate refugia, locations where taxa survive periods of regionally adverse climate, are thought to be critical for maintaining biodiversity through the glacial–interglacial climate changes of the Quaternary. A critical research need is to better integrate and reconcile the three major lines of evidence used to infer the existence of past refugia – fossil records, species distribution models and phylogeographic surveys – in order to characterize the complex spatiotemporal trajectories of species and populations in and out of refugia. Here we review the complementary strengths, limitations and new advances for these three approaches. We provide case studies to illustrate their combined application, and point the way towards new opportunities for synthesizing these disparate lines of evidence. Case studies with European beech, Qinghai spruce and Douglas-fir illustrate how the combination of these three approaches successfully resolves complex species histories not attainable from any one approach. Promising new statistical techniques can capitalize on the strengths of each method and provide a robust quantitative reconstruction of species history. Studying past refugia can help identify contemporary refugia and clarify their conservation significance, in particular by elucidating the fine-scale processes and the particular geographic locations that buffer species against rapidly changing climate.

A model-data comparison of Holocene timberline changes in the Swiss Alps reveals past and future drivers of mountain forest dynamics

Mountain vegetation is strongly affected by temperature and is expected to shift upwards with climate change. Dynamic vegetation models are often used to assess the impact of climate on vegetation and model output can be compared with paleobotanical data as a reality check. Recent paleoecological studies have revealed regional variation in the upward shift of timberlines in the Northern and Central European Alps in response to rapid warming at the Younger Dryas/Preboreal transition ca. 11700years ago, probably caused by a climatic gradient across the Alps. This contrasts with previous studies that successfully simulated the early Holocene afforestation in the (warmer) Central Alps with a chironomid-inferred temperature reconstruction from the (colder) Northern Alps. We use LandClim, a dynamic landscape vegetation model to simulate mountain forests under different temperature, soil and precipitation scenarios around Iffigsee (2065m a.s.l.) a lake in the Northwestern Swiss Alps, and compare the model output with the paleobotanical records. The model clearly overestimates the upward shift of timberline in a climate scenario that applies chironomid-inferred July-temperature anomalies to all months. However, forest establishment at 9800 cal. BP at Iffigsee is successfully simulated with lower moisture availability and monthly temperatures corrected for stronger seasonality during the early Holocene. The model-data comparison reveals a contraction in the realized niche of Abies alba due to the prominent role of anthropogenic disturbance after ca. 5000 cal. BP, which has important implications for species distribution models (SDMs) that rely on equilibrium with climate and niche stability. Under future climate projections, LandClim indicates a rapid upward shift of mountain vegetation belts by ca. 500m and treeline positions of ca. 2500m a.s.l. by the end of this century. Resulting biodiversity losses in the alpine vegetation belt might be mitigated with low-impact pastoralism to preserve species-rich alpine meadows.

Counterbalancing effects of competition for resources and facilitation against grazing in alpine snowbed communities

Alpine snowbeds are habitats where the major limiting factors for plant growth are herbivory and a small time window for growth due to late snowmelt. Despite these limitations, snowbed vegetation usually forms a dense carpet of palatable plants due to favourable abiotic conditions for plant growth within the short growing season. These environmental characteristics make snowbeds particularly interesting to study the interplay of facilitation and competition. We hypothesised an interplay between resource competition and facilitation against herbivory. Further, we investigated whether these predicted neighbour effects were species-specific and/or dependent on ontogeny, and whether the balance of positive and negative plant–plant interactions shifted along a snowmelt gradient. We determined the neighbour effects by means of neighbour removal experiments along the snowmelt gradient, and linear mixed model analyses. The results showed that the effects of neighbour removal were weak but generally consistent among species and snowmelt dates, and depended on whether biomass production or survival was considered. Higher total biomass and increased fruiting in removal plots indicated that plants competed for nutrients, water, and light, thereby supporting the hypothesis of prevailing competition for resources in snowbeds. However, the presence of neighbours reduced herbivory and thereby also facilitated survival. For plant growth the facilitative effects against herbivores in snowbeds counterbalanced competition for resources, leading to a weak negative net effect. Overall the neighbour effects were not species-specific and did not change with snowmelt date. Our finding of counterbalancing effects of competition and facilitation within a plant community is of special theoretical value for species distribution models and can explain the success of models that give primary importance to abiotic factors and tend to overlook interrelations between biotic and abiotic effects on plants.

Species-level selection reduces selfishness through competitive exclusion

Adaptation does not necessarily lead to traits which are optimal for the population. This is because selection is often the strongest at the individual or gene level. The evolution of selfishness can lead to a 'tragedy of the commons', where traits such as aggression or social cheating reduce population size and may lead to extinction. This suggests that species-level selection will result whenever species differ in the incentive to be selfish. We explore this idea in a simple model that combines individual-level selection with ecology in two interacting species. Our model is not influenced by kin or trait-group selection. We find that individual selection in combination with competitive exclusion greatly increases the likelihood that selfish species go extinct. A simple example of this would be a vertebrate species that invests heavily into squabbles over breeding sites, which is then excluded by a species that invests more into direct reproduction. A multispecies simulation shows that these extinctions result in communities containing species that are much less selfish. Our results suggest that species-level selection and community dynamics play an important role in regulating the intensity of conflicts in natural populations.

Conservation of threatened habitat types under future climate change – Lessons from plant-distribution models and current extinction trends in southern Germany

A higher risk of future range losses as a result of climate change is expected to be one of the main drivers of extinction trends in vascular plants occurring in habitat types of high conservation value. Nevertheless, the impact of the climate changes of the last 60 years on the current distribution and extinction patterns of plants is still largely unclear. We applied species distribution models to study the impact of environmental variables (climate, soil conditions, land cover, topography), on the current distribution of 18 vascular plant species characteristic of three threatened habitat types in southern Germany: (i) xero-thermophilous vegetation, (ii) mesophilous mountain grasslands (mountain hay meadows and matgrass communities), and (iii) wetland habitats (bogs, fens, and wet meadows). Climate and soil variables were the most important variables affecting plant distributions at a spatial level of 10 × 10 km. Extinction trends in our study area revealed that plant species which occur in wetland habitats faced higher extinction risks than those in xero-thermophilous vegetation, with the risk for species in mesophilous mountain grasslands being intermediary. For three plant species characteristic either of mesophilous mountain grasslands or wetland habitats we showed exemplarily that extinctions from 1950 to the present day have occurred at the edge of the species’ current climatic niche, indicating that climate change has likely been the main driver of extinction. This is largely consistent with current extinction trends reported in other studies. Our study indicates that the analysis of past extinctions is an appropriate means to assess the impact of climate change on species and that vulnerability to climate change is both species- and habitat-specific.

Spatial and temporal variation in population genetic structure of wild Nile tilapia (Oreochromis niloticus) across Africa

Background: Reconstructing the evolutionary history of a species is challenging. It often depends not only on the past biogeographic and climatic events but also the contemporary and ecological factors, such as current connectivity and habitat heterogeneity. In fact, these factors might interact with each other and shape the current species distribution. However, to what extent the current population genetic structure reflects the past and the contemporary factors is largely unknown. Here we investigated spatio-temporal genetic structures of Nile tilapia (Oreochromis niloticus) populations, across their natural distribution in Africa. While its large biogeographic distribution can cause genetic differentiation at the paleo-biogeographic scales, its restricted dispersal capacity might induce a strong genetic structure at micro-geographic scales. Results: Using nine microsatellite loci and 350 samples from ten natural populations, we found the highest genetic differentiation among the three ichthyofaunal provinces and regions (Ethiopian, Nilotic and Sudano-Sahelian) (R(ST) = 0.38 - 0.69). This result suggests the predominant effect of paleo-geographic events at macro-geographic scale. In addition, intermediate divergences were found between rivers and lakes within the regions, presumably reflecting relatively recent interruptions of gene flow between hydrographic basins (R(ST) = 0.24 - 0.32). The lowest differentiations were observed among connected populations within a basin (R(ST) = 0.015 in the Volta basin). Comparison of temporal sample series revealed subtle changes in the gene pools in a few generations (F = 0 - 0.053). The estimated effective population sizes were 23 - 143 and the estimated migration rate was moderate (m similar to 0.094 - 0.097) in the Volta populations. Conclusions: This study revealed clear hierarchical patterns of the population genetic structuring of O. niloticus in Africa. The effects of paleo-geographic and climatic events were predominant at macro-geographic scale, and the significant effect of geographic connectivity was detected at micro-geographic scale. The estimated effective population size, the moderate level of dispersal and the rapid temporal change in genetic composition might reflect a potential effect of life history strategy on population dynamics. This hypothesis deserves further investigation. The dynamic pattern revealed at micro-geographic and temporal scales appears important from a genetic resource management as well as from a biodiversity conservation point of view.

An integrated approach for reconstructing surface models of the proximal femur from sparse input data for surgical navigation

A patient-specific surface model of the proximal femur plays an important role in planning and supporting various computer-assisted surgical procedures including total hip replacement, hip resurfacing, and osteotomy of the proximal femur. The common approach to derive 3D models of the proximal femur is to use imaging techniques such as computed tomography (CT) or magnetic resonance imaging (MRI). However, the high logistic effort, the extra radiation (CT-imaging), and the large quantity of data to be acquired and processed make them less functional. In this paper, we present an integrated approach using a multi-level point distribution model (ML-PDM) to reconstruct a patient-specific model of the proximal femur from intra-operatively available sparse data. Results of experiments performed on dry cadaveric bones using dozens of 3D points are presented, as well as experiments using a limited number of 2D X-ray images, which demonstrate promising accuracy of the present approach.

Growth response to climate and drought change along an aridity gradient in the southernmost Pinus nigra relict forests

Tree populations at the rear edge of species distribution are sensitive to climate stress and drought. However, growth responses of these tree populations to those stressors may vary along climatic gradients. To analyze growth responses to climate and drought using dendrochronology in rear-edge Pinus nigra populations located along an aridity gradient. Tree-ring width chronologies were built for the twentieth century and related to monthly climatic variables, a drought index (Standardized Precipitation-Evapotranspiration Index), and two atmospheric circulation patterns (North Atlantic and Western Mediterranean Oscillations). Growth was enhanced by wet and cold previous autumns and warm late winters before tree-ring formation. The influence of the previous year conditions on growth increased during the past century. Growth was significantly related to North Atlantic and Western Mediterranean Oscillations in two out of five sites. The strongest responses of growth to the drought index were observed in the most xeric sites. Dry conditions before tree-ring formation constrain growth in rear-edge P. nigra populations. The comparisons of climate-growth responses along aridity gradients allow characterizing the sensitivity of relict stands to climate warming.