Improving generalized regression analysis for the spatial prediction of forest communities

Aim This study used data from temperate forest communities to assess: (1) five different stepwise selection methods with generalized additive models, (2) the effect of weighting absences to ensure a prevalence of 0.5, (3) the effect of limiting absences beyond the environmental envelope defined by presences, (4) four different methods for incorporating spatial autocorrelation, and (5) the effect of integrating an interaction factor defined by a regression tree on the residuals of an initial environmental model. Location State of Vaud, western Switzerland. Methods Generalized additive models (GAMs) were fitted using the grasp package (generalized regression analysis and spatial predictions, http://www.cscf.ch/grasp). Results Model selection based on cross-validation appeared to be the best compromise between model stability and performance (parsimony) among the five methods tested. Weighting absences returned models that perform better than models fitted with the original sample prevalence. This appeared to be mainly due to the impact of very low prevalence values on evaluation statistics. Removing zeroes beyond the range of presences on main environmental gradients changed the set of selected predictors, and potentially their response curve shape. Moreover, removing zeroes slightly improved model performance and stability when compared with the baseline model on the same data set. Incorporating a spatial trend predictor improved model performance and stability significantly. Even better models were obtained when including local spatial autocorrelation. A novel approach to include interactions proved to be an efficient way to account for interactions between all predictors at once. Main conclusions Models and spatial predictions of 18 forest communities were significantly improved by using either: (1) cross-validation as a model selection method, (2) weighted absences, (3) limited absences, (4) predictors accounting for spatial autocorrelation, or (5) a factor variable accounting for interactions between all predictors. The final choice of model strategy should depend on the nature of the available data and the specific study aims. Statistical evaluation is useful in searching for the best modelling practice. However, one should not neglect to consider the shapes and interpretability of response curves, as well as the resulting spatial predictions in the final assessment.

Prediction of cross-recognition of peptide-HLA A2 by Melan-A-specific cytotoxic T lymphocytes using three-dimensional quantitative structure-activity relationships.

The cross-recognition of peptides by cytotoxic T lymphocytes is a key element in immunology and in particular in peptide based immunotherapy. Here we develop three-dimensional (3D) quantitative structure-activity relationships (QSARs) to predict cross-recognition by Melan-A-specific cytotoxic T lymphocytes of peptides bound to HLA A*0201 (hereafter referred to as HLA A2). First, we predict the structure of a set of self- and pathogen-derived peptides bound to HLA A2 using a previously developed ab initio structure prediction approach [Fagerberg et al., J. Mol. Biol., 521-46 (2006)]. Second, shape and electrostatic energy calculations are performed on a 3D grid to produce similarity matrices which are combined with a genetic neural network method [So et al., J. Med. Chem., 4347-59 (1997)] to generate 3D-QSAR models. The models are extensively validated using several different approaches. During the model generation, the leave-one-out cross-validated correlation coefficient (q (2)) is used as the fitness criterion and all obtained models are evaluated based on their q (2) values. Moreover, the best model obtained for a partitioned data set is evaluated by its correlation coefficient (r = 0.92 for the external test set). The physical relevance of all models is tested using a functional dependence analysis and the robustness of the models obtained for the entire data set is confirmed using y-randomization. Finally, the validated models are tested for their utility in the setting of rational peptide design: their ability to discriminate between peptides that only contain side chain substitutions in a single secondary anchor position is evaluated. In addition, the predicted cross-recognition of the mono-substituted peptides is confirmed experimentally in chromium-release assays. These results underline the utility of 3D-QSARs in peptide mimetic design and suggest that the properties of the unbound epitope are sufficient to capture most of the information to determine the cross-recognition.

Risk prediction of prevalent diabetes in a Swiss population using a weighted genetic score--the CoLaus Study.

AIMS/HYPOTHESIS: Several susceptibility genes for type 2 diabetes have been discovered recently. Individually, these genes increase the disease risk only minimally. The goals of the present study were to determine, at the population level, the risk of diabetes in individuals who carry risk alleles within several susceptibility genes for the disease and the added value of this genetic information over the clinical predictors. METHODS: We constructed an additive genetic score using the most replicated single-nucleotide polymorphisms (SNPs) within 15 type 2 diabetes-susceptibility genes, weighting each SNP with its reported effect. We tested this score in the extensively phenotyped population-based cross-sectional CoLaus Study in Lausanne, Switzerland (n = 5,360), involving 356 diabetic individuals. RESULTS: The clinical predictors of prevalent diabetes were age, BMI, family history of diabetes, WHR, and triacylglycerol/HDL-cholesterol ratio. After adjustment for these variables, the risk of diabetes was 2.7 (95% CI 1.8-4.0, p = 0.000006) for individuals with a genetic score within the top quintile, compared with the bottom quintile. Adding the genetic score to the clinical covariates improved the area under the receiver operating characteristic curve slightly (from 0.86 to 0.87), yet significantly (p = 0.002). BMI was similar in these two extreme quintiles. CONCLUSIONS/INTERPRETATION: In this population, a simple weighted 15 SNP-based genetic score provides additional information over clinical predictors of prevalent diabetes. At this stage, however, the clinical benefit of this genetic information is limited.

Setting boundaries: brain dynamics of modal and amodal illusory shape completion in humans.

Normal visual perception requires differentiating foreground from background objects. Differences in physical attributes sometimes determine this relationship. Often such differences must instead be inferred, as when two objects or their parts have the same luminance. Modal completion refers to such perceptual "filling-in" of object borders that are accompanied by concurrent brightness enhancement, in turn termed illusory contours (ICs). Amodal completion is filling-in without concurrent brightness enhancement. Presently there are controversies regarding whether both completion processes use a common neural mechanism and whether perceptual filling-in is a bottom-up, feedforward process initiating at the lowest levels of the cortical visual pathway or commences at higher-tier regions. We previously examined modal completion (Murray et al., 2002) and provided evidence that the earliest modal IC sensitivity occurs within higher-tier object recognition areas of the lateral occipital complex (LOC). We further proposed that previous observations of IC sensitivity in lower-tier regions likely reflect feedback modulation from the LOC. The present study tested these proposals, examining the commonality between modal and amodal completion mechanisms with high-density electrical mapping, spatiotemporal topographic analyses, and the local autoregressive average distributed linear inverse source estimation. A common initial mechanism for both types of completion processes (140 msec) that manifested as a modulation in response strength within higher-tier visual areas, including the LOC and parietal structures, is demonstrated, whereas differential mechanisms were evident only at a subsequent time period (240 msec), with amodal completion relying on continued strong responses in these structures.

The twist, writhe and overall shape of supercoiled DNA change during counterion-induced transition from a loosely to a tightly interwound superhelix. Possible implications for DNA structure in vivo.

A cryo-electron microscopy study of supercoiled DNA molecules freely suspended in cryo-vitrified buffer was combined with Monte Carlo simulations and gel electrophoretic analysis to investigate the role of intersegmental electrostatic repulsion in determining the shape of supercoiled DNA molecules. It is demonstrated here that a decrease of DNA-DNA repulsion by increasing concentrations of counterions causes a higher fraction of the linking number deficit to be partitioned into writhe. When counterions reach concentrations likely to be present under in vivo conditions, naturally supercoiled plasmids adopt a tightly interwound conformation. In these tightly supercoiled DNA molecules the opposing segments of interwound superhelix seem to directly contact each other. This form of supercoiling, where two DNA helices interact laterally, may represent an important functional state of DNA. In the particular case of supercoiled minicircles (178 bp) the delta Lk = -2 topoisomers undergo a sharp structural transition from almost planar circles in low salt buffers to strongly writhed "figure-eight" conformations in buffers containing neutralizing concentrations of counterions. Possible implications of this observed structural transition in DNA are discussed.

Coloration signals the ability to cope with elevated stress hormones: effects of corticosterone on growth of barn owls are associated with melanism.

Stressful situations during development can shape the phenotype for life by provoking a trade-off between development and survival. Stress hormones, mainly glucocorticoids, play an important orchestrating role in this trade-off. Hence, how stress sensitive an animal is critically determines the phenotype and ultimately fitness. In several species, darker eumelanic individuals are less sensitive to stressful conditions than less eumelanic conspecifics, which may be due to the pleiotropic effects of genes affecting both coloration and physiological traits. We experimentally tested whether the degree of melanin-based coloration is associated with the sensitivity to an endocrine response to stressful situations in the barn owl. We artificially administered the mediator of a hormonal stress response, corticosterone, to nestlings to examine the prediction that corticosterone-induced reduction in growth rate is more pronounced in light eumelanic nestlings than in darker nest mates. To examine whether such an effect may be genetically determined, we swapped hatchlings between randomly chosen pairs of nests. We first showed that corticosterone affects growth and, thus, shapes the phenotype. Second, we found that under corticosterone administration, nestlings with large black spots grew better than nestlings with small black spots. As in the barn owl the expression of eumelanin-based coloration is heritable and not sensitive to environmental conditions, it is therefore a reliable, genetically based sign of the ability to cope with an increase in blood corticosterone level.

Validity of impedance-based equations for the prediction of total body water as measured by deuterium dilution in African women.

BACKGROUND: Little information is available on the validity of simple and indirect body-composition methods in non-Western populations. Equations for predicting body composition are population-specific, and body composition differs between blacks and whites. OBJECTIVE: We tested the hypothesis that the validity of equations for predicting total body water (TBW) from bioelectrical impedance analysis measurements is likely to depend on the racial background of the group from which the equations were derived. DESIGN: The hypothesis was tested by comparing, in 36 African women, TBW values measured by deuterium dilution with those predicted by 23 equations developed in white, African American, or African subjects. These cross-validations in our African sample were also compared, whenever possible, with results from other studies in black subjects. RESULTS: Errors in predicting TBW showed acceptable values (1.3-1.9 kg) in all cases, whereas a large range of bias (0.2-6.1 kg) was observed independently of the ethnic origin of the sample from which the equations were derived. Three equations (2 from whites and 1 from blacks) showed nonsignificant bias and could be used in Africans. In all other cases, we observed either an overestimation or underestimation of TBW with variable bias values, regardless of racial background, yielding no clear trend for validity as a function of ethnic origin. CONCLUSIONS: The findings of this cross-validation study emphasize the need for further fundamental research to explore the causes of the poor validity of TBW prediction equations across populations rather than the need to develop new prediction equations for use in Africa.

Human Primary Auditory Cortex Follows the Shape of Heschl's Gyrus.

The primary auditory cortex (PAC) is central to human auditory abilities, yet its location in the brain remains unclear. We measured the two largest tonotopic subfields of PAC (hA1 and hR) using high-resolution functional MRI at 7 T relative to the underlying anatomy of Heschl's gyrus (HG) in 10 individual human subjects. The data reveals a clear anatomical-functional relationship that, for the first time, indicates the location of PAC across the range of common morphological variants of HG (single gyri, partial duplications, and complete duplications). In 20/20 individual hemispheres, two primary mirror-symmetric tonotopic maps were clearly observed with gradients perpendicular to HG. PAC spanned both divisions of HG in cases of partial and complete duplications (11/20 hemispheres), not only the anterior division as commonly assumed. Specifically, the central union of the two primary maps (the hA1-R border) was consistently centered on the full Heschl's structure: on the gyral crown of single HGs and within the sulcal divide of duplicated HGs. The anatomical-functional variants of PAC appear to be part of a continuum, rather than distinct subtypes. These findings significantly revise HG as a marker for human PAC and suggest that tonotopic maps may have shaped HG during human evolution. Tonotopic mappings were based on only 16 min of fMRI data acquisition, so these methods can be used as an initial mapping step in future experiments designed to probe the function of specific auditory fields.

Patients' prediction of extubation success.

The spontaneous breathing trial (SBT)-relying on objective criteria assessed by the clinician-is the major diagnostic tool to determine if patients can be successfully extubated. However, little is known regarding the patient's subjective perception of autonomous breathing. We performed a prospective observational study in 211 mechanically ventilated adult patients successfully completing a SBT. Patients were randomly assigned to be interviewed during this trial regarding their prediction of extubation success. We compared post-extubation outcomes in three patient groups: patients confident (confidents; n = 115) or not (non-confidents; n = 38) of their extubation success and patients not subjected to interview (control group; n = 58). Extubation success was more frequent in confidents than in non-confidents (90 vs. 45%; p < 0.001/positive likelihood ratio = 2.00) or in the control group (90 vs. 78%; p = 0.04). On the contrary, extubation failure was more common in non-confidents than in confidents (55 vs. 10%; p < 0.001/negative likelihood ratio = 0.19). Logistic regression analysis showed that extubation success was associated with patient's prediction [OR (95% CI): 9.2 (3.74-22.42) for confidents vs.non-confidents] as well as to age [0.72 (0.66-0.78) for age 75 vs. 65 and 1.31 (1.28-1.51) for age 55 vs. 65]. Our data suggest that at the end of a sustained SBT, extubation success might be correlated to the patients' subjective perception of autonomous breathing. The results of this study should be confirmed by a large multicenter trial.

Geometry and physics of knots

KNOTS are usually categorized in terms of topological properties that are invariant under changes in a knot's spatial configuration(1-4). Here we approach knot identification from a different angle, by considering the properties of particular geometrical forms which we define as 'ideal'. For a knot with a given topology and assembled from a tube of uniform diameter, the ideal form is the geometrical configuration having the highest ratio of volume to surface area. Practically, this is equivalent to determining the shortest piece of tube that can be closed to form the knot. Because the notion of an ideal form is independent of absolute spatial scale, the length-to-diameter ratio of a tube providing an ideal representation is constant, irrespective of the tube's actual dimensions. We report the results of computer simulations which show that these ideal representations of knots have surprisingly simple geometrical properties. In particular, there is a simple linear relationship between the length-to-diameter ratio and the crossing number-the number of intersections in a two-dimensional projection of the knot averaged over all directions. We have also found that the average shape of knotted polymeric chains in thermal equilibrium is closely related to the ideal representation of the corresponding knot type. Our observations provide a link between ideal geometrical objects and the behaviour of seemingly disordered systems, and allow the prediction of properties of knotted polymers such as their electrophoretic mobility(5).

Psychoses a l'adolescence, les neurosciences ameliorent-elles la prediction? [Adolescent psychosis, can neuroscience improve prediction?]

Developments in the field of neuroscience have created a high level of interest in the subject of adolescent psychosis, particularly in relation to prediction and prevention. As the medical practice of adolescent psychosis and its treatment is characterised by a heterogeneity which is both symptomatic and evolutive, the somewhat poor prognosis of chronic development justifies the research performed: apparent indicators of schizophrenic disorders on the one hand and specific endophenotypes on the other are becoming increasingly important. The significant progresses made on the human genome show that the genetic predetermination in current psychiatric pathologies is complex and subject to moderating effects and there is therefore significant potential for nature-nurture interactions (between the environment and the genes). The road to be followed in researching the phenotypic expression of a psychosis gene is long and winding and is susceptible to many external influences at various levels with different effects. Neurobiological, neurophysiological, neuropsychological and neuroanatomical studies help to identify endophenotypes, which allow researchers to create identifying "markers" along this winding road. The endophenotypes could make it possible to redefine the nosological categories and enhance understanding of the physiopathology of schizophrenia. In a predictive approach, large-scale retrospective and prospective studies make it possible to identify risk factors, which are compatible with the neurodevelopmental hypothesis of schizophrenia. However, the predictive value of such markers or risk indicators is not yet sufficiently developed to offer a reliable early-detection method or possible schizophrenia prevention measures. Nonetheless, new developments show promise against the background of a possible future nosographic revolution, based on a paradigm shift. It is perhaps on the basis of homogeneous endophenotypes in particular that we will be able to understand what protects against, or indeed can trigger, psychosis irrespective of the clinical expression or attempts to isolate the common genetic and biological bases according to homogeneous clinical characteristics, which have to date, proved unsuccessful