Statistical inference for the optimal approximating model

In the setting of high-dimensional linear models with Gaussian noise, we investigate the possibility of confidence statements connected to model selection. Although there exist numerous procedures for adaptive (point) estimation, the construction of adaptive confidence regions is severely limited (cf. Li in Ann Stat 17:1001–1008, 1989). The present paper sheds new light on this gap. We develop exact and adaptive confidence regions for the best approximating model in terms of risk. One of our constructions is based on a multiscale procedure and a particular coupling argument. Utilizing exponential inequalities for noncentral χ2-distributions, we show that the risk and quadratic loss of all models within our confidence region are uniformly bounded by the minimal risk times a factor close to one.

IDH/MGMT-driven molecular classification of low-grade glioma is a strong predictor for long-term survival

BACKGROUND Low-grade gliomas (LGGs) are rare brain neoplasms, with survival spanning up to a few decades. Thus, accurate evaluations on how biomarkers impact survival among patients with LGG require long-term studies on samples prospectively collected over a long period. METHODS The 210 adult LGGs collected in our databank were screened for IDH1 and IDH2 mutations (IDHmut), MGMT gene promoter methylation (MGMTmet), 1p/19q loss of heterozygosity (1p19qloh), and nuclear TP53 immunopositivity (TP53pos). Multivariate survival analyses with multiple imputation of missing data were performed using either histopathology or molecular markers. Both models were compared using Akaike's information criterion (AIC). The molecular model was reduced by stepwise model selection to filter out the most critical predictors. A third model was generated to assess for various marker combinations. RESULTS Molecular parameters were better survival predictors than histology (ΔAIC = 12.5, P< .001). Forty-five percent of studied patients died. MGMTmet was positively associated with IDHmut (P< .001). In the molecular model with marker combinations, IDHmut/MGMTmet combined status had a favorable impact on overall survival, compared with IDHwt (hazard ratio [HR] = 0.33, P< .01), and even more so the triple combination, IDHmut/MGMTmet/1p19qloh (HR = 0.18, P< .001). Furthermore, IDHmut/MGMTmet/TP53pos triple combination was a significant risk factor for malignant transformation (HR = 2.75, P< .05). CONCLUSION By integrating networks of activated molecular glioma pathways, the model based on genotype better predicts prognosis than histology and, therefore, provides a more reliable tool for standardizing future treatment strategies.

Seasonal Variation in the Capacity for Plant Trait Measures to Predict Grassland Carbon and Water Fluxes

There is a need for accurate predictions of ecosystem carbon (C) and water fluxes in field conditions. Previous research has shown that ecosystem properties can be predicted from community abundance-weighted means (CWM) of plant functional traits and measures of trait variability within a community (FDvar). The capacity for traits to predict carbon (C) and water fluxes, and the seasonal dependency of these trait-function relationships has not been fully explored. Here we measured daytime C and water fluxes over four seasons in grasslands of a range of successional ages in southern England. In a model selection procedure, we related these fluxes to environmental covariates and plant biomass measures before adding CWM and FDvar plant trait measures that were scaled up from measures of individual plants grown in greenhouse conditions. Models describing fluxes in periods of low biological activity contained few predictors, which were usually abiotic factors. In more biologically active periods, models contained more predictors, including plant trait measures. Field-based plant biomass measures were generally better predictors of fluxes than CWM and FDvar traits. However, when these measures were used in combination traits accounted for additional variation. Where traits were significant predictors their identity often reflected seasonal vegetation dynamics. These results suggest that database derived trait measures can improve the prediction of ecosystem C and water fluxes. Controlled studies and those involving more detailed flux measurements are required to validate and explore these findings, a worthwhile effort given the potential for using simple vegetation measures to help predict landscape-scale fluxes.

Pre-clinical in vivo models for the screening of bone biomaterials for oral/craniofacial indications: focus on small-animal models.

Preclinical in vivo experimental studies are performed for evaluating proof-of-principle concepts, safety and possible unwanted reactions of candidate bone biomaterials before proceeding to clinical testing. Specifically, models involving small animals have been developed for screening bone biomaterials for their potential to enhance bone formation. No single model can completely recreate the anatomic, physiologic, biomechanic and functional environment of the human mouth and jaws. Relevant aspects regarding physiology, anatomy, dimensions and handling are discussed in this paper to elucidate the advantages and disadvantages of small-animal models. Model selection should be based not on the 'expertise' or capacities of the team, but rather on a scientifically solid rationale, and the animal model selected should reflect the question for which an answer is sought. The rationale for using heterotopic or orthotopic testing sites, and intraosseous, periosseous or extraskeletal defect models, is discussed. The paper also discusses the relevance of critical size defect modeling, with focus on calvarial defects in rodents. In addition, the rabbit sinus model and the capsule model in the rat mandible are presented and discussed in detail. All animal experiments should be designed with care and include sample-size and study-power calculations, thus allowing generation of meaningful data. Moreover, animal experiments are subject to ethical approval by the relevant authority. All procedures and the postoperative handling and care, including postoperative analgesics, should follow best practice.

Macular atrophy in patients with long-term anti-VEGF treatment for neovascular age-related macular degeneration.

PURPOSE To identify the prevalence and progression of macular atrophy (MA) in neovascular age-related macular degeneration (AMD) patients under long-term anti-vascular endothelial growth factor (VEGF) therapy and to determine risk factors. METHOD This retrospective study included patients with neovascular AMD and ≥30 anti-VEGF injections. Macular atrophy (MA) was measured using near infrared and spectral-domain optical coherence tomography (SD-OCT). Yearly growth rate was estimated using square-root transformation to adjust for baseline area and allow for linearization of growth rate. Multiple regression with Akaike information criterion (AIC) as model selection criterion was used to estimate the influence of various parameters on MA area. RESULTS Forty-nine eyes (47 patients, mean age 77 ± 14) were included with a mean of 48 ± 13 intravitreal anti-VEGF injections (ranibizumab:37 ± 11, aflibercept:11 ± 6, mean number of injections/year 8 ± 2.1) over a mean treatment period of 6.2 ± 1.3 years (range 4-8.5). Mean best-corrected visual acuity improved from 57 ± 17 letters at baseline (= treatment start) to 60 ± 16 letters at last follow-up. The MA prevalence within and outside the choroidal neovascularization (CNV) border at initial measurement was 45% and increased to 74%. Mean MA area increased from 1.8 ± 2.7 mm(2) within and 0.5 ± 0.98 mm(2) outside the CNV boundary to 2.7 ± 3.4 mm(2) and 1.7 ± 1.8 mm(2) , respectively. Multivariate regression determined posterior vitreous detachment (PVD) and presence/development of intraretinal cysts (IRCs) as significant factors for total MA size (R(2) = 0.16, p = 0.02). Macular atrophy (MA) area outside the CNV border was best explained by the presence of reticular pseudodrusen (RPD) and IRC (R(2) = 0.24, p = 0.02). CONCLUSION A majority of patients show MA after long-term anti-VEGF treatment. Reticular pseudodrusen (RPD), IRC and PVD but not number of injections or treatment duration seem to be associated with the MA size.

Mechanisms of species divergence through visual adaptation and sexual selection: Perspectives from a cichlid model system

The theory of ecological speciation suggests that assortative mating evolves most easily when mating preferences are;directly linked to ecological traits that are subject to divergent selection. Sensory adaptation can play a major role in this process,;because selective mating is often mediated by sexual signals: bright colours, complex song, pheromone blends and so on. When;divergent sensory adaptation affects the perception of such signals, mating patterns may change as an immediate consequence.;Alternatively, mating preferences can diverge as a result of indirect effects: assortative mating may be promoted by selection;against intermediate phenotypes that are maladapted to their (sensory) environment. For Lake Victoria cichlids, the visual environment;constitutes an important selective force that is heterogeneous across geographical and water depth gradients. We investigate;the direct and indirect effects of this heterogeneity on the evolution of female preferences for alternative male nuptial colours;(red and blue) in the genus Pundamilia. Here, we review the current evidence for divergent sensory drive in this system, extract;general principles, and discuss future perspectives

Femur Specific Polyaffine Model to Regularize the Log-Domain Demons Registration

Osteoarticular allograft transplantation is a popular treatment method in wide surgical resections with large defects. For this reason hospitals are building bone data banks. Performing the optimal allograft selection on bone banks is crucial to the surgical outcome and patient recovery. However, current approaches are very time consuming hindering an efficient selection. We present an automatic method based on registration of femur bones to overcome this limitation. We introduce a new regularization term for the log-domain demons algorithm. This term replaces the standard Gaussian smoothing with a femur specific polyaffine model. The polyaffine femur model is constructed with two affine (femoral head and condyles) and one rigid (shaft) transformation. Our main contribution in this paper is to show that the demons algorithm can be improved in specific cases with an appropriate model. We are not trying to find the most optimal polyaffine model of the femur, but the simplest model with a minimal number of parameters. There is no need to optimize for different number of regions, boundaries and choice of weights, since this fine tuning will be done automatically by a final demons relaxation step with Gaussian smoothing. The newly developed synthesis approach provides a clear anatomically motivated modeling contribution through the specific three component transformation model, and clearly shows a performance improvement (in terms of anatomical meaningful correspondences) on 146 CT images of femurs compared to a standard multiresolution demons. In addition, this simple model improves the robustness of the demons while preserving its accuracy. The ground truth are manual measurements performed by medical experts.

Recurrent antitopographic inhibition mediates competitive stimulus selection in an attention network

Topographically organized neurons represent multiple stimuli within complex visual scenes and compete for subsequent processing in higher visual centers. The underlying neural mechanisms of this process have long been elusive. We investigate an experimentally constrained model of a midbrain structure: the optic tectum and the reciprocally connected nucleus isthmi. We show that a recurrent antitopographic inhibition mediates the competitive stimulus selection between distant sensory inputs in this visual pathway. This recurrent antitopographic inhibition is fundamentally different from surround inhibition in that it projects on all locations of its input layer, except to the locus from which it receives input. At a larger scale, the model shows how a focal top-down input from a forebrain region, the arcopallial gaze field, biases the competitive stimulus selection via the combined activation of a local excitation and the recurrent antitopographic inhibition. Our findings reveal circuit mechanisms of competitive stimulus selection and should motivate a search for anatomical implementations of these mechanisms in a range of vertebrate attentional systems.

Inbreeding and selection on sex ratio in the bark beetle Xylosandrus germanus

Background Local Mate Competition (LMC) theory predicts a female should produce a more female-biased sex ratio if her sons compete with each other for mates. Because it provides quantitative predictions that can be experimentally tested, LMC is a textbook example of the predictive power of evolutionary theory. A limitation of many earlier studies in the field is that the population structure and mating system of the studied species are often estimated only indirectly. Here we use microsatellites to characterize the levels of inbreeding of the bark beetle Xylosandrus germanus, a species where the level of LMC is expected to be high. Results For three populations studied, genetic variation for our genetic markers was very low, indicative of an extremely high level of inbreeding (FIS = 0.88). There was also strong linkage disequilibrium between microsatellite loci and a very strong genetic differentiation between populations. The data suggest that matings among non-siblings are very rare (3%), although sex ratios from X. germanus in both the field and the laboratory have suggested more matings between non-sibs, and so less intense LMC. Conclusions Our results confirm that caution is needed when inferring mating systems from sex ratio data, especially when a lack of biological detail means the use of overly simple forms of the model of interest.

The SURTAVI model: proposal for a pragmatic risk stratification for patients with severe aortic stenosis

Transcatheter aortic valve implantation (TAVI) is a less invasive alternative to surgical aortic valve replacement (SAVR) for patients with symptomatic severe aortic stenosis (AS) and a high operative risk. Risk stratification plays a decisive role in the optimal selection of therapeutic strategies for AS patients. The accuracy of contemporary surgical risk algorithms for AS patients has spurred considerable debate especially in the higher risk patient population. Future trials will explore TAVI in patients at intermediate operative risk. During the design of the SURgical replacement and Transcatheter Aortic Valve Implantation (SURTAVI) trial, a novel concept of risk stratification was proposed based upon age in combination with a fixed number of predefined risk factors, which are relatively prevalent, easy to capture and with a reasonable impact on operative mortality. Retrospective application of this algorithm to a contemporary academic practice dealing with clinically significant AS patients allocates about one-fourth of these patients as being at intermediate operative risk. Further testing is required for validation of this new paradigm in risk stratification. Finally, the Heart Team, consisting of at least an interventional cardiologist and cardiothoracic surgeon, should have the decisive role in determining whether a patient could be treated with TAVI or SAVR.

Effect of Methodological and Ecological Approaches on Heterogeneity of Nest-Site Selection of a Long-Lived Vulture

The application of scientific-based conservation measures requires that sampling methodologies in studies modelling similar ecological aspects produce comparable results making easier their interpretation. We aimed to show how the choice of different methodological and ecological approaches can affect conclusions in nest-site selection studies along different Palearctic meta-populations of an indicator species. First, a multivariate analysis of the variables affecting nest-site selection in a breeding colony of cinereous vulture (Aegypius monachus) in central Spain was performed. Then, a meta-analysis was applied to establish how methodological and habitat-type factors determine differences and similarities in the results obtained by previous studies that have modelled the forest breeding habitat of the species. Our results revealed patterns in nesting-habitat modelling by the cinereous vulture throughout its whole range: steep and south-facing slopes, great cover of large trees and distance to human activities were generally selected. The ratio and situation of the studied plots (nests/random), the use of plots vs. polygons as sampling units and the number of years of data set determined the variability explained by the model. Moreover, a greater size of the breeding colony implied that ecological and geomorphological variables at landscape level were more influential. Additionally, human activities affected in greater proportion to colonies situated in Mediterranean forests. For the first time, a meta-analysis regarding the factors determining nest-site selection heterogeneity for a single species at broad scale was achieved. It is essential to homogenize and coordinate experimental design in modelling the selection of species' ecological requirements in order to avoid that differences in results among studies would be due to methodological heterogeneity. This would optimize best conservation and management practices for habitats and species in a global context.

A genetically explicit model of speciation by sensory drive within a continuous population in aquatic environments

BACKGROUND: The sensory drive hypothesis predicts that divergent sensory adaptation in different habitats may lead to premating isolation upon secondary contact of populations. Speciation by sensory drive has traditionally been treated as a special case of speciation as a byproduct of adaptation to divergent environments in geographically isolated populations. However, if habitats are heterogeneous, local adaptation in the sensory systems may cause the emergence of reproductively isolated species from a single unstructured population. In polychromatic fishes, visual sensitivity might become adapted to local ambient light regimes and the sensitivity might influence female preferences for male nuptial color. In this paper, we investigate the possibility of speciation by sensory drive as a byproduct of divergent visual adaptation within a single initially unstructured population. We use models based on explicit genetic mechanisms for color vision and nuptial coloration. RESULTS: We show that in simulations in which the adaptive evolution of visual pigments and color perception are explicitly modeled, sensory drive can promote speciation along a short selection gradient within a continuous habitat and population. We assumed that color perception evolves to adapt to the modal light environment that individuals experience and that females prefer to mate with males whose nuptial color they are most sensitive to. In our simulations color perception depends on the absorption spectra of an individual's visual pigments. Speciation occurred most frequently when the steepness of the environmental light gradient was intermediate and dispersal distance of offspring was relatively small. In addition, our results predict that mutations that cause large shifts in the wavelength of peak absorption promote speciation, whereas we did not observe speciation when peak absorption evolved by stepwise mutations with small effect. CONCLUSION: The results suggest that speciation can occur where environmental gradients create divergent selection on sensory modalities that are used in mate choice. Evidence for such gradients exists from several animal groups, and from freshwater and marine fishes in particular. The probability of speciation in a continuous population under such conditions may then critically depend on the genetic architecture of perceptual adaptation and female mate choice.

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.

Evolutionary demography of long-lived monocarpic perennials: a time-lagged integral projection model

1. The evolution of flowering strategies (when and at what size to flower) in monocarpic perennials is determined by balancing current reproduction with expected future reproduction, and these are largely determined by size-specific patterns of growth and survival. However, because of the difficulty in following long-lived individuals throughout their lives, this theory has largely been tested using short-lived species (< 5 years). 2. Here, we tested this theory using the long-lived monocarpic perennial Campanula thyrsoides which can live up to 16 years. We used a novel approach that combined permanent plot and herb chronology data from a 3-year field study to parameterize and validate integral projection models (IPMs). 3. Similar to other monocarpic species, the rosette leaves of C. thyrsoides wither over winter and so size cannot be measured in the year of flowering. We therefore extended the existing IPM framework to incorporate an additional time delay that arises because flowering demography must be predicted from rosette size in the year before flowering. 4. We found that all main demographic functions (growth, survival probability, flowering probability and fecundity) were strongly size-dependent and there was a pronounced threshold size of flowering. There was good agreement between the predicted distribution of flowering ages obtained from the IPMs and that estimated in the field. Mostly, there was good agreement between the IPM predictions and the direct quantitative field measurements regarding the demographic parameters lambda, R-0 and T. We therefore conclude that the model captures the main demographic features of the field populations. 5. Elasticity analysis indicated that changes in the survival and growth function had the largest effect (c. 80%) on lambda and this was considerably larger than in short-lived monocarps. We found only weak selection pressure operating on the observed flowering strategy which was close to the predicted evolutionary stable strategy. 6. Synthesis. The extended IPM accurately described the demography of a long-lived monocarpic perennial using data collected over a relatively short period. We could show that the evolution of flowering strategies in short- and long-lived monocarps seem to follow the same general rules but with a longevity-related emphasis on survival over fecundity.