Structural components in functional data

Analyzing functional data often leads to finding common factors, for which functional principal component analysis proves to be a useful tool to summarize and characterize the random variation in a function space. The representation in terms of eigenfunctions is optimal in the sense of L-2 approximation. However, the eigenfunctions are not always directed towards an interesting and interpretable direction in the context of functional data and thus could obscure the underlying structure. To overcome such difficulty, an alternative to functional principal component analysis is proposed that produces directed components which may be more informative and easier to interpret. These structural components are similar to principal components, but are adapted to situations in which the domain of the function may be decomposed into disjoint intervals such that there is effectively independence between intervals and positive correlation within intervals. The approach is demonstrated with synthetic examples as well as real data. Properties for special cases are also studied.

Increased audiovisual integration in cochlear-implanted deaf patients: independent components analysis of longitudinal positron emission tomography data.

It has been demonstrated in earlier studies that patients with a cochlear implant have increased abilities for audio-visual integration because the crude information transmitted by the cochlear implant requires the persistent use of the complementary speech information from the visual channel. The brain network for these abilities needs to be clarified. We used an independent components analysis (ICA) of the activation (H2 (15) O) positron emission tomography data to explore occipito-temporal brain activity in post-lingually deaf patients with unilaterally implanted cochlear implants at several months post-implantation (T1), shortly after implantation (T0) and in normal hearing controls. In between-group analysis, patients at T1 had greater blood flow in the left middle temporal cortex as compared with T0 and normal hearing controls. In within-group analysis, patients at T0 had a task-related ICA component in the visual cortex, and patients at T1 had one task-related ICA component in the left middle temporal cortex and the other in the visual cortex. The time courses of temporal and visual activities during the positron emission tomography examination at T1 were highly correlated, meaning that synchronized integrative activity occurred. The greater involvement of the visual cortex and its close coupling with the temporal cortex at T1 confirm the importance of audio-visual integration in more experienced cochlear implant subjects at the cortical level.

Differential effect of Charlson index and its components on in-hospital mortality and health costs

Introduction: The Charlson index (Charlson, 1987) is a commonly used comorbidity index in outcome studies. Still, the use of different weights makes its calculation cumbersome, while the sum of its components (comorbidities) is easier to compute. In this study, we assessed the effects of 1) the Charlson index adapted for the Swiss population and 2) the sum of its components (number of comorbidities, maximum 15) on a) in-hospital deaths and b) cost of hospitalization. Methods: Anonymous data was obtained from the administrative database of the department of internal medicine of the Lausanne University Hospital (CHUV). All hospitalizations of adult (>=18 years) patients occurring between 2003 and 2011 were included. For each hospitalization, the Charlson index and the number of comorbidities were calculated. Analyses were conducted using Stata. Results: Data from 32,741 hospitalizations occurring between 2003 and 2011 was analyzed. On bivariate analysis, both the Charlson index and the number of comorbidities were significantly and positively associated with in hospital death. Conversely, multivariate adjustment for age, gender and calendar year using Cox regression showed that the association was no longer significant for the number of comorbidities (table). On bivariate analysis, hospitalization costs increased both with Charlson index and with number of comorbidities, but the increase was much steeper for the number of comorbidities (figure). Robust regression after adjusting for age, gender, calendar year and duration of hospital stay showed that the increase in one comorbidity led to an average increase in hospital costs of 321 CHF (95% CI: 272 to 370), while the increase in one score point of the Charlson index led to a decrease in hospital costs of 49 CHF (95% CI: 31 to 67). Conclusion: Charlson index is better than the number of comorbidities in predicting in-hospital death. Conversely, the number of comorbidities significantly increases hospital costs.

Inferring ultraviolet anatomical exposure patterns while distinguishing the relative contribution of radiation components

Exposure to solar ultraviolet (UV) radiation is the main causative factor for skin cancer. UV exposure depends on environmental and individual factors, but individual exposure data remain scarce. While ground UV irradiance is monitored via different techniques, it is difficult to translate such observations into human UV exposure or dose because of confounding factors. A multi-disciplinary collaboration developed a model predicting the dose and distribution of UV exposure on the basis of ground irradiation and morphological data. Standard 3D computer graphics techniques were adapted to develop a simulation tool that estimates solar exposure of a virtual manikin depicted as a triangle mesh surface. The amount of solar energy received by various body locations is computed for direct, diffuse and reflected radiation separately. Dosimetric measurements obtained in field conditions were used to assess the model performance. The model predicted exposure to solar UV adequately with a symmetric mean absolute percentage error of 13% and half of the predictions within 17% range of the measurements. Using this tool, solar UV exposure patterns were investigated with respect to the relative contribution of the direct, diffuse and reflected radiation. Exposure doses for various body parts and exposure scenarios of a standing individual were assessed using erythemally-weighted UV ground irradiance data measured in 2009 at Payerne, Switzerland as input. For most anatomical sites, mean daily doses were high (typically 6.2-14.6 Standard Erythemal Dose, SED) and exceeded recommended exposure values. Direct exposure was important during specific periods (e. g. midday during summer), but contributed moderately to the annual dose, ranging from 15 to 24% for vertical and horizontal body parts, respectively. Diffuse irradiation explained about 80% of the cumulative annual exposure dose.

Genetic and phenotypic architecture of metabolic syndrome-associated components in dyslipidemic and normolipidemic subjects: the GEMS Study.

Atherogenic dyslipidemia, manifest by low HDL-cholesterol and high TG levels, is an important component of ATP-III defined metabolic syndrome. Here, we dissected the phenotypic and genetic architecture of these traits by assessing their relationships with other metabolically relevant measures, including plasma adipo-cytokines, highly sensitive C-reactive protein (hsCRP) and LDL particle size, in a large family data set (n=2800) and in an independent set of dyslipidemic cases (n=716) and normolipidemic controls (n=1073). We explored the relationships among these phenotypes using variable clustering and then estimated their genetic heritabilities and cross-trait correlations. In families, four clusters explained 61% of the total variance, with one adiposity-related cluster (including hsCRP), one BP-related cluster, and two lipid-related clusters (HDL-C, TG, adiponectin and LDL particle size; apoB and non-HDL-C). A similar structure was observed in dyslipidemic cases and normolipidemic controls. The genetic correlations in the families largely paralleled the phenotype clustering results, suggesting that common genes having pleiotropic effects contributed to the correlations observed. In summary, our analyses support a model of metabolic syndrome with two major components, body fat and lipids, each with two subcomponents, and quantifies their degree of overlap with each other and with metabolic-syndrome related measures (adipokines, LDL particle size and hsCRP).

Inactivation des pathogènes des produits sanguins labiles: entre enjeux financiers et conséquences possibles à long terme [Pathogen reduction of blood components: from financial issues to possible long-term consequences].

Pathogen inactivation of blood products represents a global and major paradigm shift in transfusion medicine. In the next near future, it is likely that most blood products will be inactivated by various physicochemical approaches. The concept of blood safety will be challenged as well as transfusion medicine practice, notably for donor selection or biological qualification. In this context, it seems mandatory to develop analytical economic approaches by assessing costs-benefits ratio of blood transfusion as well as to set up cohorts of patients based on hemovigilance networks allowing rigorous scientific analysis of the benefits and the risks of blood transfusion at short- and long-term.

Distinct neuronal circuits underlying the bahavioral and physiological components of the fear response

Abstract : Expression of fear involves changes in a number of behavioral and physiological parameters that are triggered by the central amygdala (CeA). The fear circuit also includes a series of brain stem nuclei that are the final effectors of the changes induced by the fear reaction. The CeA expresses many different neuropeptide receptors that can modulate fear responses. Today, the precise organization and the modulation of projections from the amygdala to the brain stem are still poorly understood. The aim of this project was to better understand the organization and the modulation of the fear circuit. To investigate this we first determined whether the CeA is composed of separate neuronal populations, where each one projects to specific brain stem nuclei, or whether single CeA neurons project to several nuclei. For this purpose, we first selected two brain stem nuclei implicated in the modulation of different components of the fear reactions, the periaqueductal gray (implicated in freezing) and the nucleus of solitary tract (implicated in heart rate modulation). We then performed double injections of two different retrograde tracers in these two nuclei and we quantified the subsequent presence of co-labelling in the CeA. We found that neurons projecting to the PAG and to the NTS are organized in separate populations. Subsequent electrophysiological recordings of the two populations revealed that PAG and NTS projecting neurons also have different electrophysiological characteristics. We then verified in vitro whether the neurons projecting to different brain stem nuclei express specific combinations of neuropeptide receptors, and whether a neuropeptide acting pre-synaptically (oxytocin) specifically modulates one of these two projections. We did not find differences at the level of expression of neurópeptide receptors, but we observed that oxytocin, a neuropeptide with anxiolytic properties, modulates PAG projecting neurons without affecting NTS projecting neurons. As oxytocin appeared to specifically modulate projections to the PAG, involved in the modulation of the freezing reaction, but did not affect the projections to the NTS, implicated in the modulation of cardiovascular parameters, we verified how this modulation translates in living animals. We investigated the effects of infra-amygdala injection of oxytocin on cardiovascular and behavioral changes induced by contextual fear conditioning. We found that oxytocin decreased the freezing response without affecting the cardiovascular system. Finally, as neuropeptides are considered potential future anxiolytics, we investigated whether diazepam and oxytocin, acting on the same circuit, had additive effects. This question was addressed exclusively with an in vitro electrophysiological approach. We obtained that oxytocin and diazepam, when co-applied, had an additive effect on both synaptic transmission and neuronal activity. These results open new perspectives for the possible clinical applications of oxytocin. Résumé : L'expression de la peur est accompagnée par de nombreux changements physiologiques et comportementaux qui sont déclenchés par l'amygdale centrale (CeA). Le circuit inclue aussi une série de noyaux du tronc cérébrale qui sont les effecteurs des différentes composantes de la réaction de peur. On sait que CeA envoie des projections aux noyaux du tronc cérébral et que ces neurones expriment une grande variété de récepteurs aux neuropeptides. Par contre, l'organisation des projections, ainsi que la modulation de ces projections par les neuropeptides reste encore peu connue. Avec ce projet, on premièrement voulu déterminer si CeA est composée de populations neuronales séparées qui projettent vers un noyau spécifique, ou bien si chaque neurones envoie des projections vers plusieurs noyaux. A ce propos, on a effectué des doubles injections de deux traceurs rétrogrades différentes dans deux noyaux du tronc cérébral impliqués dans des différentes composantes des réactions de peur. On a injecté la substance grise périaqueducale (PAG), qui est impliquée dans la réponse d'immobilisation, ainsi que le noyau du tractus solitaire (NTS) qui est responsable des changements cardiovasculaires. On a ensuite quantifié la présence de neurones contenant les deux traceurs dans CeA. On a trouvé que la plupart des neurones de l'amygdale centrale projettent vers un noyau spécifique, et on peut donc dire que l'amygdale semble être composée de populations neuronales séparées. On a ensuite mesuré les caractéristiques électrophysiologiques de ces deux projections et on a trouvé des différences substantielles concernant la résistance membranaire, la capacitance, le potentiel membranaire de repos ainsi que la fréquence des potentiels d'action spontanés. Puis, comme beaucoup de neuropéptides dans l'amygdale exercent un effet modulatoire sûr les réactions de peur et sur l'anxiété, on a étudié les effets directs et indirects d'une série de neuropeptides sur les différentes projections pour évaluer s'il y a des neuropeptides qui agissent spécifiquement sur une. On n'a pas trouvé de différences entre neurones qui projettent vers le PAG et neurones qui projettent vers le NTS concernant les effets de neuropeptides qui agissent directement sur ces cellules. Par contre, on a trouvé que l'ocytocine, un neuropeptide qui se lie à des récepteurs dans la partie latérale de l'amygdale centrale et inhibe de façon indirecte les neurones de l'amygdala centrale médiale, module les projections vers le PAG sans affecter celles qui vont vers le NTS. Comme le PAG est impliqué dans la réponse d'immobilisation, alors que le NTS est impliqué dans la modulation cardiovasculaire, on a ensuite étudié les effets de l'ocytocine injectée dans l'amygdale de rat vivants sur les réactions de peur conditionnées. On a trouvé que l'ocytocine diminue la réponse d'immobilisation sans par contre affecter la réponse cardiovasculaire. Pour terminer, on a vérifié si l'ocytocine potentialise les effets d'un médicament anxiolytique, le diazeparn. Avec une étude in vitro on a trouvé qu'une co-application d'ocytocine et diazeparn résulte en un effet additionnel à la fois sur la transmission synaptique ainsi que sur l'activité neuronale des neurones de l'amygdale centrale médiale. Ces résultats ouvrent des nouvelles perspectives pour une potentielle utilisation clinique de l'ocytocine.

Principal components of functional connectivity: A new approach to study dynamic brain connectivity during rest.

Functional connectivity (FC) as measured by correlation between fMRI BOLD time courses of distinct brain regions has revealed meaningful organization of spontaneous fluctuations in the resting brain. However, an increasing amount of evidence points to non-stationarity of FC; i.e., FC dynamically changes over time reflecting additional and rich information about brain organization, but representing new challenges for analysis and interpretation. Here, we propose a data-driven approach based on principal component analysis (PCA) to reveal hidden patterns of coherent FC dynamics across multiple subjects. We demonstrate the feasibility and relevance of this new approach by examining the differences in dynamic FC between 13 healthy control subjects and 15 minimally disabled relapse-remitting multiple sclerosis patients. We estimated whole-brain dynamic FC of regionally-averaged BOLD activity using sliding time windows. We then used PCA to identify FC patterns, termed "eigenconnectivities", that reflect meaningful patterns in FC fluctuations. We then assessed the contributions of these patterns to the dynamic FC at any given time point and identified a network of connections centered on the default-mode network with altered contribution in patients. Our results complement traditional stationary analyses, and reveal novel insights into brain connectivity dynamics and their modulation in a neurodegenerative disease.

Swain: Stata module to correct the SEM chi-square overidentification test in small sample sizes or complex models. Statistical Software Components S457617, Boston College Department of Economics.

Swain corrects the chi-square overidentification test (i.e., likelihood ratio test of fit) for structural equation models whethr with or without latent variables. The chi-square statistic is asymptotically correct; however, it does not behave as expected in small samples and/or when the model is complex (cf. Herzog, Boomsma, & Reinecke, 2007). Thus, particularly in situations where the ratio of sample size (n) to the number of parameters estimated (p) is relatively small (i.e., the p to n ratio is large), the chi-square test will tend to overreject correctly specified models. To obtain a closer approximation to the distribution of the chi-square statistic, Swain (1975) developed a correction; this scaling factor, which converges to 1 asymptotically, is multiplied with the chi-square statistic. The correction better approximates the chi-square distribution resulting in more appropriate Type 1 reject error rates (see Herzog & Boomsma, 2009; Herzog, et al., 2007).

Age-specific fitness components and their temporal variation in the barn owl.

Theory predicts that temporal variability plays an important role in the evolution of life histories, but empirical studies evaluating this prediction are rare. In constant environments, fitness can be measured by the population growth rate lambda, and the sensitivity of lambda to changes in fitness components estimates selection on these traits. In variable environments, fitness is measured by the stochastic growth rate lambda(S), and stochastic sensitivities estimate selection pressure. Here we examine age-specific schedules for reproduction and survival in a barn owl population (Tyto alba). We estimated how temporal variability affected fitness and selection, accounting for sampling variance. Despite large sample sizes of old individuals, we found no strong evidence for senescence. The most variable fitness components were associated with reproduction. Survival was less variable. Stochastic simulations showed that the observed variation decreased fitness by about 30%, but the sensitivities of lambda and lambda(S) to changes in all fitness components were almost equal, suggesting that temporal variation had negligible effects on selection. We obtained these results despite high observed variability in the fitness components and relatively short generation time of the study organism, a situation in which temporal variability should be particularly important for natural selection and early senescence is expected.

Geographical and temporal body size variation in a reptile: roles of sex, ecology, phylogeny and ecology structured in phylogeny.

Geographical body size variation has long interested evolutionary biologists, and a range of mechanisms have been proposed to explain the observed patterns. It is considered to be more puzzling in ectotherms than in endotherms, and integrative approaches are necessary for testing non-exclusive alternative mechanisms. Using lacertid lizards as a model, we adopted an integrative approach, testing different hypotheses for both sexes while incorporating temporal, spatial, and phylogenetic autocorrelation at the individual level. We used data on the Spanish Sand Racer species group from a field survey to disentangle different sources of body size variation through environmental and individual genetic data, while accounting for temporal and spatial autocorrelation. A variation partitioning method was applied to separate independent and shared components of ecology and phylogeny, and estimated their significance. Then, we fed-back our models by controlling for relevant independent components. The pattern was consistent with the geographical Bergmann's cline and the experimental temperature-size rule: adults were larger at lower temperatures (and/or higher elevations). This result was confirmed with additional multi-year independent data-set derived from the literature. Variation partitioning showed no sex differences in phylogenetic inertia but showed sex differences in the independent component of ecology; primarily due to growth differences. Interestingly, only after controlling for independent components did primary productivity also emerge as an important predictor explaining size variation in both sexes. This study highlights the importance of integrating individual-based genetic information, relevant ecological parameters, and temporal and spatial autocorrelation in sex-specific models to detect potentially important hidden effects. Our individual-based approach devoted to extract and control for independent components was useful to reveal hidden effects linked with alternative non-exclusive hypothesis, such as those of primary productivity. Also, including measurement date allowed disentangling and controlling for short-term temporal autocorrelation reflecting sex-specific growth plasticity.

Comparison of the polynomial model against explicit measurements of noise components for different mammography systems.

Given the adverse impact of image noise on the perception of important clinical details in digital mammography, routine quality control measurements should include an evaluation of noise. The European Guidelines, for example, employ a second-order polynomial fit of pixel variance as a function of detector air kerma (DAK) to decompose noise into quantum, electronic and fixed pattern (FP) components and assess the DAK range where quantum noise dominates. This work examines the robustness of the polynomial method against an explicit noise decomposition method. The two methods were applied to variance and noise power spectrum (NPS) data from six digital mammography units. Twenty homogeneously exposed images were acquired with PMMA blocks for target DAKs ranging from 6.25 to 1600 µGy. Both methods were explored for the effects of data weighting and squared fit coefficients during the curve fitting, the influence of the additional filter material (2 mm Al versus 40 mm PMMA) and noise de-trending. Finally, spatial stationarity of noise was assessed.Data weighting improved noise model fitting over large DAK ranges, especially at low detector exposures. The polynomial and explicit decompositions generally agreed for quantum and electronic noise but FP noise fraction was consistently underestimated by the polynomial method. Noise decomposition as a function of position in the image showed limited noise stationarity, especially for FP noise; thus the position of the region of interest (ROI) used for noise decomposition may influence fractional noise composition. The ROI area and position used in the Guidelines offer an acceptable estimation of noise components. While there are limitations to the polynomial model, when used with care and with appropriate data weighting, the method offers a simple and robust means of examining the detector noise components as a function of detector exposure.

NOTCH1 can initiate NF-κB activation via cytosolic interactions with components of the T cell Signalosome.

T cell stimulation requires the input and integration of external signals. Signaling through the T cell receptor (TCR) is known to induce formation of the membrane-tethered CBM complex, comprising CARMA1, BCL10, and MALT1, which is required for TCR-mediated NF-κB activation. TCR signaling has been shown to activate NOTCH proteins, transmembrane receptors also implicated in NF-κB activation. However, the link between TCR-mediated NOTCH signaling and early events leading to induction of NF-κB activity remains unclear. In this report, we demonstrate a novel cytosolic function for NOTCH1 and show that it is essential to CBM complex formation. Using a model of skin allograft rejection, we show in vivo that NOTCH1 acts in the same functional pathway as PKCθ, a T cell-specific kinase important for CBM assembly and classical NF-κB activation. We further demonstrate in vitro NOTCH1 associates physically with PKCθ and CARMA1 in the cytosol. Unexpectedly, when NOTCH1 expression was abrogated using RNAi approaches, interactions between CARMA1, BCL10, and MALT1 were lost. This failure in CBM assembly reduced inhibitor of kappa B alpha phosphorylation and diminished NF-κB-DNA binding. Finally, using a luciferase gene reporter assay, we show the intracellular domain of NOTCH1 can initiate robust NF-κB activity in stimulated T cells, even when NOTCH1 is excluded from the nucleus through modifications that restrict it to the cytoplasm or hold it tethered to the membrane. Collectively, these observations provide evidence that NOTCH1 may facilitate early events during T cell activation by nucleating the CBM complex and initiating NF-κB signaling.

Cenomanian (early Late Cretaceous) ammonoid faunas of Western Europe - Part II: Diversity patterns and the end-Cenomanian anoxic event

Diversity patterns of ammonoids are analyzed and compared with the timing of anoxic deposits around the Cenomanian/Turonian (C/T) boundary in the Vocontian, Anglo-Paris, and Monster basins of Western Europe. Differing from most previous studies, which concentrate on a narrow time span bracketing the C/T boundary, the present analysis covers the latest Albian to Early Turonian interval for which a high resolution, ammonoid-based biochronology, including 34 Unitary Associations zones, is now available. During the latest Albian-Middle Cenomanian interval, species richness of ammonoids reveals a dynamical equilibrium oscillating around an average of 20 species, whereas the Late Cenomanian-Early Turonian interval displays an equilibrium centered on an average value of 6 species. The abrupt transition between these two successive equilibria lasted no longer than two Unitary Associations. The onset of the decline of species richness thus largely predates the spread of oxygen-poor water masses onto the shelves, while minimal values of species richness coincide with the Cenomanian-Turonian boundary only. The decline of species richness during the entire Late Cenomanian seems to result from lower origination percentages rather than from higher extinction percentages. This result is also supported by the absence of statistically significant changes in the extinction probabilities of the poly-cohorts. Separate analyses of species richness for acanthoceratids and heteromorphs, the two essential components of the Cenomanian ammonoid community, reveal that heteromorphs declined sooner than acanthoceratids. Moreover, acanthoceratids showed a later decline at the genus level than at the species level. Such a decoupling is accompanied by a significant increase in morphological disparity of acanthoceratids, which is expressed by the appearance of new genera. Last, during the Late Cenomanian, paedomorphic processes, juvenile innovations and reductions of adult size dominated the evolutionary radiation of acanthoceratids. Hence, the decrease in ammonoid species richness and their major evolutionary changes significantly predates the spread of anoxic deposits. Other environmental constraints such as global flooding of platforms, warmer and more equable climate, as well as productivity changes better correlate with the timing of diversity changes and evolutionary patterns of ammonoids and therefore, provide more likely causative mechanisms than anoxia alone.