Multiple-scale hydrological and geophysical data integration through non-linear Bayesian sequential simulation

Significant progress has been made with regard to the quantitative integration of geophysical and hydrological data at the local scale. However, extending the corresponding approaches to the scale of a field site represents a major, and as-of-yet largely unresolved, challenge. To address this problem, we have developed downscaling procedure based on a non-linear Bayesian sequential simulation approach. The main objective of this algorithm is to estimate the value of the sparsely sampled hydraulic conductivity at non-sampled locations based on its relation to the electrical conductivity logged at collocated wells and surface resistivity measurements, which are available throughout the studied site. The in situ relationship between the hydraulic and electrical conductivities is described through a non-parametric multivariatekernel density function. Then a stochastic integration of low-resolution, large-scale electrical resistivity tomography (ERT) data in combination with high-resolution, local-scale downhole measurements of the hydraulic and electrical conductivities is applied. The overall viability of this downscaling approach is tested and validated by comparing flow and transport simulation through the original and the upscaled hydraulic conductivity fields. Our results indicate that the proposed procedure allows obtaining remarkably faithful estimates of the regional-scale hydraulic conductivity structure and correspondingly reliable predictions of the transport characteristics over relatively long distances.

Human invariant V alpha 24-J alpha Q TCR supports the development of CD1d-dependent NK1.1+ and NK1.1- T cells in transgenic mice.

A sizable fraction of T cells expressing the NK cell marker NK1.1 (NKT cells) bear a very conserved TCR, characterized by homologous invariant (inv.) TCR V alpha 24-J alpha Q and V alpha 14-J alpha 18 rearrangements in humans and mice, respectively, and are thus defined as inv. NKT cells. Because human inv. NKT cells recognize mouse CD1d in vitro, we wondered whether a human inv. V alpha 24 TCR could be selected in vivo by mouse ligands presented by CD1d, thereby supporting the development of inv. NKT cells in mice. Therefore, we generated transgenic (Tg) mice expressing the human inv. V alpha 24-J alpha Q TCR chain in all T cells. The expression of the human inv. V alpha 24 TCR in TCR C alpha(-/-) mice indeed rescues the development of inv. NKT cells, which home preferentially to the liver and respond to the CD1d-restricted ligand alpha-galactosylceramide (alpha-GalCer). However, unlike inv. NKT cells from non-Tg mice, the majority of NKT cells in V alpha 24 Tg mice display a double-negative phenotype, as well as a significant increase in TCR V beta 7 and a corresponding decrease in TCR V beta 8.2 use. Despite the forced expression of the human CD1d-restricted TCR in C alpha(-/-) mice, staining with mCD1d-alpha-GalCer tetramers reveals that the absolute numbers of peripheral CD1d-dependent T lymphocytes increase at most by 2-fold. This increase is accounted for mainly by an increased fraction of NK1.1(-) T cells that bind CD1d-alpha-GalCer tetramers. These findings indicate that human inv. V alpha 24 TCR supports the development of CD1d-dependent lymphocytes in mice, and argue for a tight homeostatic control on the total number of inv. NKT cells. Thus, human inv. V alpha 24 TCR-expressing mice are a valuable model to study different aspects of the inv. NKT cell subset.

Cutting edge: influence of the TCR Vbeta domain on the selection of semi-invariant NKT cells by endogenous ligands.

Invariant Valpha14 (Valpha14i) NKT cells are a murine CD1d-dependent regulatory T cell subset characterized by a Valpha14-Jalpha18 rearrangement and expression of mostly Vbeta8.2 and Vbeta7. Whereas the TCR Vbeta domain influences the binding avidity of the Valpha14i TCR for CD1d-alpha-galactosylceramide complexes, with Vbeta8.2 conferring higher avidity binding than Vbeta7, a possible impact of the TCR Vbeta domain on Valpha14i NKT cell selection by endogenous ligands has not been studied. In this study, we show that thymic selection of Vbeta7(+), but not Vbeta8.2(+), Valpha14i NKT cells is favored in situations where endogenous ligand concentration or TCRalpha-chain avidity are suboptimal. Furthermore, thymic Vbeta7(+) Valpha14i NKT cells were preferentially selected in vitro in response to CD1d-dependent presentation of endogenous ligands or exogenously added self ligand isoglobotrihexosylceramide. Collectively, our data demonstrate that the TCR Vbeta domain influences the selection of Valpha14i NKT cells by endogenous ligands, presumably because Vbeta7 confers higher avidity binding.

Cutting edge: influence of the TCR V beta domain on the avidity of CD1d:alpha-galactosylceramide binding by invariant V alpha 14 NKT cells.

CD1d tetramers loaded with alpha-galactosylceramide (alpha-GalCer) bind selectively to mouse invariant Valpha14 (Valpha14i) NKT cells and their human counterparts. Whereas tetramer binding strictly depends on the expression of a Valpha14-Jalpha18 chain in murine NKT cells, the associated beta-chain (typically expressing Vbeta8.2 or Vbeta7) appears not to influence tetramer binding. In this study, we describe novel alpha-GalCer-loaded mouse and human CD1d-IgG1 dimers, which revealed an unexpected influence of the TCR-beta chain on the avidity of CD1d:alpha-GalCer binding. A subset of Valpha14i NKT cells clearly discriminated alpha-GalCer bound to mouse or human CD1d on the basis of avidity differences conferred by the Vbeta domain of the TCR-beta chain, with Vbeta8.2 conferring higher avidity binding than Vbeta7.

Local Distance-Based Generalized Linear Models using the dbstats package for R

This paper introduces local distance-based generalized linear models. These models extend (weighted) distance-based linear models firstly with the generalized linear model concept, then by localizing. Distances between individuals are the only predictor information needed to fit these models. Therefore they are applicable to mixed (qualitative and quantitative) explanatory variables or when the regressor is of functional type. Models can be fitted and analysed with the R package dbstats, which implements several distancebased prediction methods.

Characterization of tumor antigen specific CD8+ T cells and semi-invariant Vα24/Vβ11 NKT cells in tumor invaded liver

Summary: Detailed knowledge on tumor antigen expression and specific immune cells is required for a rational design of immunotherapy for patients with tumor invaded liver. In this study, we confirmed that Cancer/Testis (CT) tumor-associated antigens are frequently expressed in hepatocellular carcinoma (HCC) and searched for the presence of CD8+ T cells specific for these antigens. In 2/10 HLA-A2+ patients with HCC, we found that MAGE-A10 and/or SSX-2 specific CD8+ T cells naturally responded to the disease, since they were enriched in tumor lesions but not in non-tumoral liver. Isolated T cells specifically and strongly killed tumor cells in vitro, suggesting that these CTL were selected in vivo for high avidity antigen recognition, providing the rational for specific immunotherapy of HCC, based on immunization with CT antigens such as MAGE-Al 0 and SSX-2. Type 1 NKT cells express an invariant TCR α chain (Vα24.1α18, paired with Vβ11 in human) and share a specific reactivity to αGalactosylceramide (αGC) presented by CD1d. These cells can display paradoxical immuno-regulatory properties including strong anti-tumor effects upon αGC administration in murine models. To understand why NKT cells were not sufficiently protective against tumor development in patients with tumor invaded liver, we characterized the diversity of Vα24/Vβ11 NKT cells in healthy donors (HD) and cancer patients: NKT cells from HD and patients were generally diverse in terms of TCR β chain (Vβ11) variability and NKT cells from HD showed a variable recognition of αGC loaded CD 1 d multimers. Vα24/ Vβ11 NKT cells can be divided in 3 populations, the CD4, DN (CD4-/CD8-) and CD8 NKT cell subsets that show distinct ability of cytokine production. In addition, our functional analysis revealed that DN and CD8 subsets displayed a higher cytolytic potential and a weaker IFNγ release than the CD4 NKT cell subset. NKT cell subsets were variably represented in the blood of HD and cancer patients. However, HD with high NKT cell frequencies displayed an enrichment of the DN and CD8 subsets, and few of them were suggestive of an oligoclonal expansion in vivo. Comparable NKT cell frequencies were found between blood, non-tumoral liver and tumor of patients. In contrast, we identified a gradual enrichment of CD4 NKT cells from blood to the liver and to the tumor, together with a decrease of DN and CD8 NKT cell subsets. Most patient derived NKT cells were unresponsive upon αGalactosylceramide stimulation ex vivo; NKT cells from few patients displayed a weak responsiveness with different cytokine polarization. The NKT cell repertoire was thus different in tumor tissue, suggesting that CD4 NKT cells infiltrating tumors may be detrimental for protection against tumors and instead may favour the tumor growth/recurrence as recently reported in mice. Résumé en français scientifique : Afin de développer le traitement des patients porteurs d'une tumeur dans le foie par immunothérapie, de nouvelles connaissances sont requises concernant l'expression d'antigènes par les tumeurs et les cellules immunitaires spécifiques de ces antigènes. Nous avons vérifié que des antigènes associés aux tumeurs, tels que les antigènes « Cancer-Testis » (CT), sont fréquemment exprimés par le carcinome hepatocéllulaire (CHC). La recherche de lymphocytes T CD8+ spécifiques (CTL) de ces antigènes a révélé que des CTL spécifiques de MAGE-A10 et/ou SSX-2 ont répondu naturellement à la tumeur chez 2/10 patients étudiés. Ces cellules étaient présentes dans les lésions tumorales mais pas dans le foie adjacent. De plus, ces CTL ont démontré une activité cytolytique forte et spécifique contre les cellules tumorales in vitro, ce qui suggère que ces CTL ont été sélectionnés pour une haute avidité de reconnaissance de l'antigène in vivo. Ces données fournissent une base pour l'immunothérapie spécifique du CHC, en proposant de cibler les antigènes CT tels que MAGE-A10 ou SSX-2. Les cellules NKT de type 1 ont une chaîne α de TCR qui est invariante (chez l'homme, Vα24Jα18, apparié avec Vβ11) et reconnaissent spécifiquement l'αGalactosylceramide (αGC) présenté par CD1d. Ces cellules ont des propriétés immuno¬régulatrices qui peuvent être parfois contradictoires et leur activation par l'αGC induit une forte protection anti-tumorale chez la souris: Afin de comprendre pourquoi ces cellules ne sont pas assez protectrices contre le développement des tumeurs dans le foie chez l'homme, nous avons étudié la diversité des cellules NKT Vα24/Vβ11 d'individus sains (IS) et de patients cancéreux. Les cellules NKT peuvent être sous-divisées en 3 populations : Les CD4, DN (CD4- /CD8-) ou CDS, qui ont la capacité de produire des cytokines différentes. Nos analyses fonctionnelles ont aussi révélé que les sous-populations DN et CD8 ont un potentiel cytolytique plus élevé et une production d'IFNγ plus faible que la sous-population CD4. Ces sous-populations sont représentées de manière variable dans le sang des IS ou des patients. Cependant, les IS avec un taux élevé de cellules NKT ont un enrichissement des sous- populations DN ou CDS, et certains suggèrent qu'il s'agit d'une expansion oligo-clonale in vivo. Les patients avaient des fréquences comparables de cellules NKT entre le sang, le foie et la tumeur. Par contre, la sous-population CD4 était progressivement enrichie du sang vers le foie et la tumeur, tandis que les sous-populations DN ou CD8 était perdues. La plupart des cellules NKT des patients ne réagissaient pas lors de stimulation avec l'αGC ex vivo et les cellules NKT de quelques patients répondaient faiblement et avec des polarisations de cytokines différentes. Ces données suggèrent que les cellules NKT CD4, prédominantes dans les tumeurs, sont inefficaces pour la lutte anti-tumorale et pourraient même favoriser la croissance ou la récurrence tumorale. Donc, une mobilisation spécifique des cellules NKT CD4 négatives par immunothérapie pourrait favoriser l'immunité contre des tumeurs chez l'homme. Résumé en français pour un large public Au sein des globules blancs, les lymphocytes T expriment un récepteur (le TCR), qui est propre à chacun d'entre eux et leur permet d'accrocher de manière très spécifique une molécule appelée antigène. Ce TCR est employé par les lymphocytes pour inspecter les antigènes associés avec des molécules présentatrices à la surface des autres cellules. Les lymphocytes T CD8 reconnaissent un fragment de protéine (ou peptide), qui est présenté par une des molécules du Complexe Majeur d'Histocompatibilité de classe I et tuent la cellule qui présente ce peptide. Ils sont ainsi bien adaptés pour éliminer les cellules qui présentent un peptide issu d'un virus quand la cellule est infectée. D'autres cellules T CD8 reconnaissent des peptides comme les antigènes CT, qui sont produits anormalement par les cellules cancéreuses. Nous avons confirmé que les antigènes CT sont fréquemment exprimés par le cancer du foie. Nous avons également identifié des cellules T CD8 spécifiques d'antigènes CT dans la tumeur, mais pas dans le foie normal de 2 patients sur 10. Cela signifie que ces lymphocytes peuvent être naturellement activés contre la tumeur et sont capables de la trouver. De plus les lymphocytes issus d'un patient ont démontré une forte sensibilité pour reconnaître l'antigène et tuent spécifiquement les cellules tumorales. Les antigènes CT représentent donc des cibles intéressantes qui pourront être intégrés dans des vaccins thérapeutiques du cancer du foie. De cette manière, les cellules T CD8 du patient lui-même pourront être induites à détruire de manière spécifique les cellules cancéreuses. Un nouveau type de lymphocytes T a été récemment découvert: les lymphocytes NKT. Quand ils reconnaissent un glycolipide présenté par la molécule CD1d, ils sont capables, de manière encore incomprise, d'initier, d'augmenter, ou à l'inverse d'inhiber la défense immunitaire. Ces cellules NKT ont démontré qu'elles jouent un rôle important dans la défense contre les tumeurs et particulièrement dans le foie des souris. Nous avons étudié les cellules NKT de patients atteints d'une tumeur dans le foie, afin de comprendre pourquoi elles ne sont pas assez protectrice chez l'homme. Les lymphocytes NKT peuvent être sous-divisés en 3 populations: Les CD4, les DN (CD4-/CD8-) et les CD8. Ces 3 classes de NKT peuvent produire différents signaux chimiques appelés cytokines. Contrairement aux cellules NKT DN ou CDS, seules les cellules NKT CD4 sont capables de produire des cytokines qui sont défavorables pour la défense anti-tumorale. Par ailleurs nous avons trouvé que les cellules NKT CD4 tuent moins bien les cellules cancéreuses que les cellules NKT DN ou CD8. L'analyse des cellules NKT, fraîchement extraites du sang, du foie et de la tumeur de patients a révélé que les cellules NKT CD4 sont progressivement enrichies du sang vers le foie et la tumeur. La large prédominance des NKT CD4 à l'intérieur des tumeurs suggère que, chez l'homme, ces cellules sont inappropriées pour la lutte anti-tumorale. Par ailleurs, la plupart des cellules NKT de patients n'étaient pas capables de produire des cytokines après stimulation avec un antigène. Cela explique également pourquoi ces cellules ne protègent pas contre les tumeurs dans le foie.

Impact of reduced cerebral perfusion pressure on outcome after severe traumatic brain injury is dependent on brain tissue oxygen pressure

INTRODUCTION. Reduced cerebral perfusion pressure (CPP) may worsen secondary damage and outcome after severe traumatic brain injury (TBI), however the optimal management of CPP is still debated. STUDY HYPOTHESIS: We hypothesized that the impact of CPP on outcome is related to brain tissue oxygen tension (PbtO2) level and that reduced CPP may worsen TBI prognosis when it is associated with brain hypoxia. DESIGN. Retrospective analysis of prospective database. METHODS. We analyzed 103 patients with severe TBI who underwent continuous PbtO2 and CPP monitoring for an average of 5 days. For each patient, duration of reduced CPP (\60 mm Hg) and brain hypoxia (PbtO2\15 mm Hg for[30 min [1]) was calculated with linear interpolation method and the relationship between CPP and PbtO2 was analyzed with Pearson's linear correlation coefficient. Outcome at 30 days was assessed with the Glasgow Outcome Score (GOS), dichotomized as good (GOS 4-5) versus poor (GOS 1-3). Multivariable associations with outcome were analyzed with stepwise forward logistic regression. RESULTS. Reduced CPP (n=790 episodes; mean duration 10.2 ± 12.3 h) was observed in 75 (74%) patients and was frequently associated with brain hypoxia (46/75; 61%). Episodes where reduced CPP were associated with normal brain oxygen did not differ significantly between patients with poor versus those with good outcome (8.2 ± 8.3 vs. 6.5 ± 9.7 h; P=0.35). In contrast, time where reduced CPP occurred simultaneously with brain hypoxia was longer in patients with poor than in those with good outcome (3.3±7.4 vs. 0.8±2.3 h; P=0.02). Outcome was significantly worse in patients who had both reduced CPP and brain hypoxia (61% had GOS 1-3 vs. 17% in those with reduced CPP but no brain hypoxia; P\0.01). Patients in whom a positive CPP-PbtO2 correlation (r[0.3) was found also were more likely to have poor outcome (69 vs. 31% in patients with no CPP-PbtO2 correlation; P\0.01). Brain hypoxia was an independent risk factor of poor prognosis (odds ratio for favorable outcome of 0.89 [95% CI 0.79-1.00] per hour spent with a PbtO2\15 mm Hg; P=0.05, adjusted for CPP, age, GCS, Marshall CT and APACHE II). CONCLUSIONS. Low CPP may significantly worsen outcome after severe TBI when it is associated with brain tissue hypoxia. PbtO2-targeted management of CPP may optimize TBI therapy and improve outcome of head-injured patients.

Impact of normalization on miRNA microarray expression profiling.

Profiling miRNA levels in cells with miRNA microarrays is becoming a widely used technique. Although normalization methods for mRNA gene expression arrays are well established, miRNA array normalization has so far not been investigated in detail. In this study we investigate the impact of normalization on data generated with the Agilent miRNA array platform. We have developed a method to select nonchanging miRNAs (invariants) and use them to compute linear regression normalization coefficients or variance stabilizing normalization (VSN) parameters. We compared the invariants normalization to normalization by scaling, quantile, and VSN with default parameters as well as to no normalization using samples with strong differential expression of miRNAs (heart-brain comparison) and samples where only a few miRNAs are affected (by p53 overexpression in squamous carcinoma cells versus control). All normalization methods performed better than no normalization. Normalization procedures based on the set of invariants and quantile were the most robust over all experimental conditions tested. Our method of invariant selection and normalization is not limited to Agilent miRNA arrays and can be applied to other data sets including those from one color miRNA microarray platforms, focused gene expression arrays, and gene expression analysis using quantitative PCR.

Laboratory estimation of the effects of increasing temperatures on the duration of gonotrophic cycle of Anopheles albimanus (Diptera: Culicidae)

The increase of malaria transmission in the Pacific Coast of Colombia during the occurrence of El Niño warm event has been found not to be linked to increases in the density of the vector Anopheles albimanus, but to other temperature-sensitive variables such as longevity, duration of the gonotrophic cycle or the sporogonic period of Plasmodium. The present study estimated the effects of temperature on duration of the gonotrophic cycle and on maturation of the ovaries of An. albimanus. Blood fed adult mosquitoes were exposed to temperatures of 24, 27, and 30°C, held individually in oviposition cages and assessed at 12 h intervals. At 24, 27, and 30°C the mean development time of the oocytes was 91.2 h (95% C.I.: 86.5-96), 66.2 h (61.5-70.8), and 73.1 h (64-82.3), respectively. The mean duration of the gonotrophic cycle for these three temperatures was 88.4 h (81.88-94.9), 75 h (71.4-78.7), and 69.1 h (64.6-73.6) respectively. These findings indicate that both parameters in An. albimanus are reduced when temperatures rose from 24 to 30°C, in a nonlinear manner. According to these results the increase in malaria transmission during El Niño in Colombia could be associated with a shortening of the gonotrophic cycle in malaria vectors, which could enhance the frequency of man-vector contact, affecting the incidence of the disease.

Binary classification of dyslipidemia from the waist-to-hip ratio and body mass index: a comparison of linear, logistic, and CART models.

BACKGROUND: We sought to improve upon previously published statistical modeling strategies for binary classification of dyslipidemia for general population screening purposes based on the waist-to-hip circumference ratio and body mass index anthropometric measurements. METHODS: Study subjects were participants in WHO-MONICA population-based surveys conducted in two Swiss regions. Outcome variables were based on the total serum cholesterol to high density lipoprotein cholesterol ratio. The other potential predictor variables were gender, age, current cigarette smoking, and hypertension. The models investigated were: (i) linear regression; (ii) logistic classification; (iii) regression trees; (iv) classification trees (iii and iv are collectively known as "CART"). Binary classification performance of the region-specific models was externally validated by classifying the subjects from the other region. RESULTS: Waist-to-hip circumference ratio and body mass index remained modest predictors of dyslipidemia. Correct classification rates for all models were 60-80%, with marked gender differences. Gender-specific models provided only small gains in classification. The external validations provided assurance about the stability of the models. CONCLUSIONS: There were no striking differences between either the algebraic (i, ii) vs. non-algebraic (iii, iv), or the regression (i, iii) vs. classification (ii, iv) modeling approaches. Anticipated advantages of the CART vs. simple additive linear and logistic models were less than expected in this particular application with a relatively small set of predictor variables. CART models may be more useful when considering main effects and interactions between larger sets of predictor variables.

A topological invariant to predict the three-dimensional writhe of ideal configurations of knots and links.

We present herein a topological invariant of oriented alternating knots and links that predicts the three-dimensional (3D) writhe of the ideal geometrical configuration of the considered knot/link. The fact that we can correlate a geometrical property of a given configuration with a topological invariant supports the notion that the ideal configuration contains important information about knots and links. The importance of the concept of ideal configuration was already suggested by the good correlation between the 3D writhe of ideal knot configurations and the ensemble average of the 3D writhe of random configurations of the considered knots. The values of the new invariant are quantized: multiples of 4/7 for links with an odd number of components (including knots) and 2/7 plus multiples of 4/7 for links with an even number of components.

Convex linear combination processes for compositions

Aitchison and Bacon-Shone (1999) considered convex linear combinations ofcompositions. In other words, they investigated compositions of compositions, wherethe mixing composition follows a logistic Normal distribution (or a perturbationprocess) and the compositions being mixed follow a logistic Normal distribution. Inthis paper, I investigate the extension to situations where the mixing compositionvaries with a number of dimensions. Examples would be where the mixingproportions vary with time or distance or a combination of the two. Practicalsituations include a river where the mixing proportions vary along the river, or acrossa lake and possibly with a time trend. This is illustrated with a dataset similar to thatused in the Aitchison and Bacon-Shone paper, which looked at how pollution in aloch depended on the pollution in the three rivers that feed the loch. Here, I explicitlymodel the variation in the linear combination across the loch, assuming that the meanof the logistic Normal distribution depends on the river flows and relative distancefrom the source origins