Cost-effective experimental design to support modeling of concentration-response functions

Modeling concentration-response function became extremely popular in ecotoxicology during the last decade. Indeed, modeling allows determining the total response pattern of a given substance. However, reliable modeling is consuming in term of data, which is in contradiction with the current trend in ecotoxicology, which aims to reduce, for cost and ethical reasons, the number of data produced during an experiment. It is therefore crucial to determine experimental design in a cost-effective manner. In this paper, we propose to use the theory of locally D-optimal designs to determine the set of concentrations to be tested so that the parameters of the concentration-response function can be estimated with high precision. We illustrated this approach by determining the locally D-optimal designs to estimate the toxicity of the herbicide dinoseb on daphnids and algae. The results show that the number of concentrations to be tested is often equal to the number of parameters and often related to the their meaning, i.e. they are located close to the parameters. Furthermore, the results show that the locally D-optimal design often has the minimal number of support points and is not much sensitive to small changes in nominal values of the parameters. In order to reduce the experimental cost and the use of test organisms, especially in case of long-term studies, reliable nominal values may therefore be fixed based on prior knowledge and literature research instead of on preliminary experiments

Merlin, a "Magic" Linker between Extracellular Cues and Intracellular Signaling Pathways that Regulate Cell Motility, Proliferation, and Survival.

Genetic alterations of neurofibromatosis type 2 (NF2) gene lead to the development of schwannomas, meningiomas, and ependymomas. Mutations of NF2 gene were also found in thyroid cancer, mesothelioma, and melanoma, suggesting that it functions as a tumor suppressor in a wide spectrum of cells. The product of NF2 gene is merlin (moesin-ezrin-radixin-like protein), a member of the Band 4.1 superfamily proteins. Merlin shares significant sequence homology with the ERM (Ezrin-Radixin-Moesin) family proteins and serves as a linker between transmembrane proteins and the actin-cytoskeleton. Merlin is a multifunctional protein and involved in integrating and regulating the extracellular cues and intracellular signaling pathways that control cell fate, shape, proliferation, survival, and motility. Recent studies showed that merlin regulates the cell-cell and cell-matrix adhesions and functions of the cell surface adhesion/extracellular matrix receptors including CD44 and that merlin and CD44 antagonize each other's function and work upstream of the mammalian Hippo signaling pathway. Furthermore, merlin plays important roles in stabilizing the contact inhibition of proliferation and in regulating activities of several receptor tyrosine kinases. Accumulating data also suggested an emerging role of merlin as a negative regulator of growth and progression of several non-NF2 associated cancer types. Together, these recent advances have improved our basic understanding about merlin function, its regulation, and the major signaling pathways regulated by merlin and provided the foundation for future translation of these findings into the clinic for patients bearing the cancers in which merlin function and/or its downstream signaling pathways are impaired or altered.

Prestress in the extracellular matrix sensitizes latent TGF-β1 for activation.

Integrin-mediated force application induces a conformational change in latent TGF-β1 that leads to the release of the active form of the growth factor from the extracellular matrix (ECM). Mechanical activation of TGF-β1 is currently understood as an acute process that depends on the contractile force of cells. However, we show that ECM remodeling, preceding the activation step, mechanically primes latent TGF-β1 akin to loading a mechanical spring. Cell-based assays and unique strain devices were used to produce a cell-derived ECM of controlled organization and prestrain. Mechanically conditioned ECM served as a substrate to measure the efficacy of TGF-β1 activation after cell contraction or direct force application using magnetic microbeads. The release of active TGF-β1 was always higher from prestrained ECM as compared with unorganized and/or relaxed ECM. The finding that ECM prestrain regulates the bioavailability of TGF-β1 is important to understand the context of diseases that involve excessive ECM remodeling, such as fibrosis or cancer.

Loss of mitochondrial functions associated with azole resistance in Candida glabrata results in enhanced virulence in mice.

Mitochondrial dysfunction is one of the possible mechanisms by which azole resistance can occur in Candida glabrata. Cells with mitochondrial DNA deficiency (so-called "petite mutants") upregulate ATP binding cassette (ABC) transporter genes and thus display increased resistance to azoles. Isolation of such C. glabrata mutants from patients receiving antifungal therapy or prophylaxis has been rarely reported. In this study, we characterized two sequential and related C. glabrata isolates recovered from the same patient undergoing azole therapy. The first isolate (BPY40) was azole susceptible (fluconazole MIC, 4 μg/ml), and the second (BPY41) was azole resistant (fluconazole MIC, >256 μg/ml). BPY41 exhibited mitochondrial dysfunction and upregulation of the ABC transporter genes C. glabrata CDR1 (CgCDR1), CgCDR2, and CgSNQ2. We next assessed whether mitochondrial dysfunction conferred a selective advantage during host infection by testing the virulence of BPY40 and BPY41 in mice. Surprisingly, even with in vitro growth deficiency compared to BPY40, BPY41 was more virulent (as judged by mortality and fungal tissue burden) than BPY40 in both systemic and vaginal murine infection models. The increased virulence of the petite mutant correlated with a drastic gain of fitness in mice compared to that of its parental isolate. To understand this unexpected feature, genome-wide changes in gene expression driven by the petite mutation were analyzed by use of microarrays during in vitro growth. Enrichment of specific biological processes (oxido-reductive metabolism and the stress response) was observed in BPY41, all of which was consistent with mitochondrial dysfunction. Finally, some genes involved in cell wall remodelling were upregulated in BPY41 compared to BPY40, which may partially explain the enhanced virulence of BPY41. In conclusion, this study shows for the first time that mitochondrial dysfunction selected in vivo under azole therapy, even if strongly affecting in vitro growth characteristics, can confer a selective advantage under host conditions, allowing the C. glabrata mutant to be more virulent than wild-type isolates.

Estimation of jump-diffusion processes via empirical characteristic functions

Preface The starting point for this work and eventually the subject of the whole thesis was the question: how to estimate parameters of the affine stochastic volatility jump-diffusion models. These models are very important for contingent claim pricing. Their major advantage, availability T of analytical solutions for characteristic functions, made them the models of choice for many theoretical constructions and practical applications. At the same time, estimation of parameters of stochastic volatility jump-diffusion models is not a straightforward task. The problem is coming from the variance process, which is non-observable. There are several estimation methodologies that deal with estimation problems of latent variables. One appeared to be particularly interesting. It proposes the estimator that in contrast to the other methods requires neither discretization nor simulation of the process: the Continuous Empirical Characteristic function estimator (EGF) based on the unconditional characteristic function. However, the procedure was derived only for the stochastic volatility models without jumps. Thus, it has become the subject of my research. This thesis consists of three parts. Each one is written as independent and self contained article. At the same time, questions that are answered by the second and third parts of this Work arise naturally from the issues investigated and results obtained in the first one. The first chapter is the theoretical foundation of the thesis. It proposes an estimation procedure for the stochastic volatility models with jumps both in the asset price and variance processes. The estimation procedure is based on the joint unconditional characteristic function for the stochastic process. The major analytical result of this part as well as of the whole thesis is the closed form expression for the joint unconditional characteristic function for the stochastic volatility jump-diffusion models. The empirical part of the chapter suggests that besides a stochastic volatility, jumps both in the mean and the volatility equation are relevant for modelling returns of the S&P500 index, which has been chosen as a general representative of the stock asset class. Hence, the next question is: what jump process to use to model returns of the S&P500. The decision about the jump process in the framework of the affine jump- diffusion models boils down to defining the intensity of the compound Poisson process, a constant or some function of state variables, and to choosing the distribution of the jump size. While the jump in the variance process is usually assumed to be exponential, there are at least three distributions of the jump size which are currently used for the asset log-prices: normal, exponential and double exponential. The second part of this thesis shows that normal jumps in the asset log-returns should be used if we are to model S&P500 index by a stochastic volatility jump-diffusion model. This is a surprising result. Exponential distribution has fatter tails and for this reason either exponential or double exponential jump size was expected to provide the best it of the stochastic volatility jump-diffusion models to the data. The idea of testing the efficiency of the Continuous ECF estimator on the simulated data has already appeared when the first estimation results of the first chapter were obtained. In the absence of a benchmark or any ground for comparison it is unreasonable to be sure that our parameter estimates and the true parameters of the models coincide. The conclusion of the second chapter provides one more reason to do that kind of test. Thus, the third part of this thesis concentrates on the estimation of parameters of stochastic volatility jump- diffusion models on the basis of the asset price time-series simulated from various "true" parameter sets. The goal is to show that the Continuous ECF estimator based on the joint unconditional characteristic function is capable of finding the true parameters. And, the third chapter proves that our estimator indeed has the ability to do so. Once it is clear that the Continuous ECF estimator based on the unconditional characteristic function is working, the next question does not wait to appear. The question is whether the computation effort can be reduced without affecting the efficiency of the estimator, or whether the efficiency of the estimator can be improved without dramatically increasing the computational burden. The efficiency of the Continuous ECF estimator depends on the number of dimensions of the joint unconditional characteristic function which is used for its construction. Theoretically, the more dimensions there are, the more efficient is the estimation procedure. In practice, however, this relationship is not so straightforward due to the increasing computational difficulties. The second chapter, for example, in addition to the choice of the jump process, discusses the possibility of using the marginal, i.e. one-dimensional, unconditional characteristic function in the estimation instead of the joint, bi-dimensional, unconditional characteristic function. As result, the preference for one or the other depends on the model to be estimated. Thus, the computational effort can be reduced in some cases without affecting the efficiency of the estimator. The improvement of the estimator s efficiency by increasing its dimensionality faces more difficulties. The third chapter of this thesis, in addition to what was discussed above, compares the performance of the estimators with bi- and three-dimensional unconditional characteristic functions on the simulated data. It shows that the theoretical efficiency of the Continuous ECF estimator based on the three-dimensional unconditional characteristic function is not attainable in practice, at least for the moment, due to the limitations on the computer power and optimization toolboxes available to the general public. Thus, the Continuous ECF estimator based on the joint, bi-dimensional, unconditional characteristic function has all the reasons to exist and to be used for the estimation of parameters of the stochastic volatility jump-diffusion models.

High-density collagen gel tubes as a matrix for primary human bladder smooth muscle cells.

Tissue-engineered grafts for the urinary tract are being investigated for the potential treatment of several urologic diseases. These grafts, predominantly tubular-shaped, usually require in vitro culture prior to implantation to allow cell engraftment on initially cell-free scaffolds. We have developed a method to produce tubular-shaped collagen scaffolds based on plastic compression. Our approach produces a ready cell-seeded graft that does not need further in vitro culture prior to implantation. The tubular collagen scaffolds were in particular investigated for their structural, mechanical and biological properties. The resulting construct showed an especially high collagen density, and was characterized by favorable mechanical properties assessed by axial extension and radial dilation. Young modulus in particular was greater than non-compressed collagen tubes. Seeding densities affected proliferation rate of primary human bladder smooth muscle cells. An optimal seeding density of 10(6) cells per construct resulted in a 25-fold increase in Alamar blue-based fluorescence after 2 wk in culture. These high-density collagen gel tubes, ready seeded with smooth muscle cells could be further seeded with urothelial cells, drastically shortening the production time of graft for urinary tract regeneration.

Hierarchy of interactive functions in father-mother-baby three-way games

In developmental research, the family has mainly been studied through dyadic interaction. Three-way interactions have received less attention, partly because of their complexity. This difficulty may be overcome by distinguishing between four hierarchically embedded functions in three-way interactions: (1) participation (inclusion of all participants), (2) organization (partners keeping to their roles), (3) focalization (sharing a common focus) and (4) affective contact (being in tune). We document this hierarchical model on a sample of 80 families observed in the Lausanne Trilogue Play situation across four different sites. Hierarchy between functions was demonstrated by means of Guttman scalability coefficient. Given the importance of the child's development in a threesome, the pertinence of this model for family assessment is discussed.

Epithelial sodium channel/degenerin family of ion channels: a variety of functions for a shared structure.

The recently discovered epithelial sodium channel (ENaC)/degenerin (DEG) gene family encodes sodium channels involved in various cell functions in metazoans. Subfamilies found in invertebrates or mammals are functionally distinct. The degenerins in Caenorhabditis elegans participate in mechanotransduction in neuronal cells, FaNaC in snails is a ligand-gated channel activated by neuropeptides, and the Drosophila subfamily is expressed in gonads and neurons. In mammals, ENaC mediates Na+ transport in epithelia and is essential for sodium homeostasis. The ASIC genes encode proton-gated cation channels in both the central and peripheral nervous system that could be involved in pain transduction. This review summarizes the physiological roles of the different channels belonging to this family, their biophysical and pharmacological characteristics, and the emerging knowledge of their molecular structure. Although functionally different, the ENaC/DEG family members share functional domains that are involved in the control of channel activity and in the formation of the pore. The functional heterogeneity among the members of the ENaC/DEG channel family provides a unique opportunity to address the molecular basis of basic channel functions such as activation by ligands, mechanotransduction, ionic selectivity, or block by pharmacological ligands.

Natural killer cell genes and functions after kidney transplantation

Natural Killer (NK) cells are of special interest in solid organ transplantation (SOT) because classical immunosuppressive drugs could enhance NK cells activity.We studied NK cells after kidney transplantation in three different situations. First, we analysed the peripheral repertoire reconstitution and function of NK cells after a polyclonal rabbit anti-thymocytes globulin (rATG) induction therapy, in 20 patients transplanted with living donor and with a low immunological risk. Second, we analysed the influence of KIR genes on the risk of CMV primo-infection or reactivation in 224 transplanted patients during the first year. Finally, we studied the risk of rejection and graft function during the first 5 years according to the KIR genes. Our study demonstrates that after an intial drop, NK cell reconstitution is fast with a ratio of CD56+/CD3− cells versus CD3+ cells that remains identical. The fraction of NK cells expressing the inhibitory receptor NKG2A significantly increases and the activating receptor NKG2D decreases after transplantation to retrieve the pretransplantation value after one year. The secretion of INF-f × and the cytotoxicity is maintained over time after transplantation. Then, we demonstrated that the presence of 2 KIR missing ligands and a large number of activating KIR gene protected against CMV primo-infection or reactivation during the first year post transplantation. Finally, the KIR genes and their HLA ligands do not influence the long term graft function after univariate and multivariate analysis. Our data suggest that despite the modification of the receptor repertoire, NK cell activity is preserved. NK cells are an important player of the immune response in the first year after transplantation mainly thanks to their anti-infectious activity.

Brain Connectivity Analysis in Children. Effects of Prematurity. and Prenatal Growth Restriction

Introduction: Survival of children born prematurely or with very low birth weight has increased dramatically, but the long term developmental outcome remains unknown. Many children have deficits in cognitive capacities, in particular involving executive domains and those disabilities are likely to involve a central nervous system deficit. To understand their neurostructural origin, we use DTI. Structurally segregated and functionally regions of the cerebral cortex are interconnected by a dense network of axonal pathways. We noninvasively map these pathways across cortical hemispheres and construct normalized structural connection matrices derived from DTI MR tractography. Group comparisons of brain connectivity reveal significant changes in fiber density in case of children with poor intrauterine grown and extremely premature children (gestational age<28 weeks at birth) compared to control subjects. This changes suggest a link between cortico-axonal pathways and the central nervous system deficit. Methods: Sixty premature born infants (5-6 years old) were scanned on clinical 3T scanner (Magnetom Trio, Siemens Medical Solutions, Erlangen, Germany) at two hospitals (HUG, Geneva and CHUV, Lausanne). For each subject, T1-weighted MPRAGE images (TR/TE=2500/2.91,TI=1100, resolution=1x1x1mm, matrix=256x154) and DTI images (30 directions, TR/TE=10200/107, in-plane resolution=1.8x1.8x2mm, 64 axial, matrix=112x112) were acquired. Parent(s) provided written consent on prior ethical board approval. The extraction of the Whole Brain Structural Connectivity Matrix was performed following (Cammoun, 2009 and Hagmann, 2008). The MPARGE images were registered using an affine registration to the non-weighted-DTI and WM-GM segmentation performed on it. In order to have equal anatomical localization among subjects, 66 cortical regions with anatomical landmarks were created using the curvature information, i.e. sulcus and gyrus (Cammoun et al, 2007; Fischl et al, 2004; Desikan et al, 2006) with freesurfer software (http://surfer.nmr.mgh.harvard.edu/). Tractography was performed in WM using an algorithm especially designed for DTI/DSI data (Hagmann et al., 2007) and both information were then combined in a matrix. Each row and column of the matrix corresponds to a particular ROI. Each cell of index (i,j) represents the fiber density of the bundle connecting the ROIs i and j. Subdividing each cortical region, we obtained 4 Connectivity Matrices of different resolution (33, 66, 125 and 250 ROI/hemisphere) for each subject . Subjects were sorted in 3 different groups, namely (1) control, (2) Intrauterine Growth Restriction (IUGR), (3) Extreme Prematurity (EP), depending on their gestational age, weight and percentile-weight score at birth. Group-to-group comparisons were performed between groups (1)-(2) and (1)-(3). The mean age at examination of the three groups were similar. Results: Quantitative analysis were performed between groups to determine fibers density differences. For each group, a mean connectivity matrix with 33ROI/hemisphere resolution was computed. On the other hand, for all matrix resolutions (33,66,125,250 ROI/hemisphere), the number of bundles were computed and averaged. As seen in figure 1, EP and IUGR subjects present an overall reduction of fibers density in both interhemispherical and intrahemispherical connections. This is given quantitatively in table 1. IUGR subjects presents a higher percentage of missing fiber bundles than EP when compared to control subjects (~16% against 11%). When comparing both groups to control subjects, for the EP subjects, the occipito-parietal regions seem less interhemispherically connected whilst the intrahemispherical networks present lack of fiber density in the lymbic system. Children born with IUGR, have similar reductions in interhemispherical connections than the EP. However, the cuneus and precuneus connections with the precentral and paracentral lobe are even lower than in the case of the EP. For the intrahemispherical connections the IUGR group preset a loss of fiber density between the deep gray matter structures (striatum) and the frontal and middlefrontal poles, connections typically involved in the control of executive functions. For the qualitative analysis, a t-test comparing number of bundles (p-value<0.05) gave some preliminary significant results (figure 2). Again, even if both IUGR and EP appear to have significantly less connections comparing to the control subjects, the IUGR cohort seems to present a higher lack of fiber density specially relying the cuneus, precuneus and parietal areas. In terms of fiber density, preliminary Wilcoxon tests seem to validate the hypothesis set by the previous analysis. Conclusions: The goal of this study was to determine the effect of extreme prematurity and poor intrauterine growth on neurostructural development at the age of 6 years-old. This data indicates that differences in connectivity may well be the basis for the neurostructural and neuropsychological deficit described in these populations in the absence of overt brain lesions (Inder TE, 2005; Borradori-Tolsa, 2004; Dubois, 2008). Indeed, we suggest that IUGR and prematurity leads to alteration of connectivity between brain structures, especially in occipito-parietal and frontal lobes for EP and frontal and middletemporal poles for IUGR. Overall, IUGR children have a higher loss of connectivity in the overall connectivity matrix than EP children. In both cases, the localized alteration of connectivity suggests a direct link between cortico-axonal pathways and the central nervous system deficit. Our next step is to link these connectivity alterations to the performance in executive function tests.

RAE-1, a novel PHR binding protein, is required for axon termination and synapse formation in Caenorhabditis elegans.

Previous studies in Caenorhabditis elegans showed that RPM-1 (Regulator of Presynaptic Morphology-1) regulates axon termination and synapse formation. To understand the mechanism of how rpm-1 functions, we have used mass spectrometry to identify RPM-1 binding proteins, and have identified RAE-1 (RNA Export protein-1) as an evolutionarily conserved binding partner. We define a RAE-1 binding region in RPM-1, and show that this binding interaction is conserved and also occurs between Rae1 and the human ortholog of RPM-1 called Pam (protein associated with Myc). rae-1 loss of function causes similar axon and synapse defects, and synergizes genetically with two other RPM-1 binding proteins, GLO-4 and FSN-1. Further, we show that RAE-1 colocalizes with RPM-1 in neurons, and that rae-1 functions downstream of rpm-1. These studies establish a novel postmitotic function for rae-1 in neuronal development.

Behavior of nuclear matrix proteins during camptothecin-induced apoptosis in HL-60 human leukemia cells.

In this study we focused our attention on the behavior of four nuclear matrix proteins during the various stages of apoptosis in the HL-60 cell line exposed to the DNA topoisomerase I inhibitor, camptothecin. We have examined the following antigens by immunocytochemical techniques: (i) the 180-kDa nucleolar isoform of DNA topoisomerase II; (ii) a 126-kDa polypeptide of nuclear bodies; (iii) a 125-kDa protein; and (iv) a 160-kDa polypeptide which are known to be components of the matrix inner network. Indirect immunofluorescence experiments were performed to follow these nuclear matrix antigens during apoptosis. Moreover, the ultrastructural localization of both 125- and 160-kDa proteins was investigated by electron microscope immunocytochemistry with gold-conjugated secondary antibodies. While the antibody to the nucleolar isoform of DNA topoisomerase II gave a fluorescent pattern that was well-maintained until the late phases of apoptosis, the other three nuclear antigens showed marked modifications in their distribution. A common feature, particularly evident for 125- and 160-kDa proteins, was their absence from cap-shaped chromatin marginations, whereas they were present in the areas of remaining decondensed chromatin. The 126-kDa polypeptide concentrated progressively in an irregular mass at the opposite side of the crescentic caps and then broke up in fine spots. The 125- and 160-kDa proteins localized in the nucleolus and precisely within certain granules which are known to appear in the nucleolar area after camptothecin administration. These results show that, in addition to the well-known chromatin changes, nuclear organization undergoes other rearrangements during the apoptotic process.

Characterization and classification of matrix effects in biological samples analyses.

An exhaustive classification of matrix effects occurring when a sample preparation is performed prior to liquid-chromatography coupled to mass spectrometry (LC-MS) analyses was proposed. A total of eight different situations were identified allowing the recognition of the matrix effect typology via the calculation of four recovery values. A set of 198 compounds was used to evaluate matrix effects after solid phase extraction (SPE) from plasma or urine samples prior to LC-ESI-MS analysis. Matrix effect identification was achieved for all compounds and classified through an organization chart. Only 17% of the tested compounds did not present significant matrix effects.

GAP-independent functions of DLC1 in metastasis.

Metastases are responsible for most cancer-related deaths. One of the hallmarks of metastatic cells is increased motility and migration through extracellular matrixes. These processes rely on specific small GTPases, in particular those of the Rho family. Deleted in liver cancer-1 (DLC1) is a tumor suppressor that bears a RhoGAP activity. This protein is lost in most cancers, allowing malignant cells to proliferate and disseminate in a Rho-dependent manner. However, DLC1 is also a scaffold protein involved in alternative pathways leading to tumor and metastasis suppressor activities. Recently, substantial information has been gathered on these mechanisms and this review is aiming at describing the potential and known alternative GAP-independent mechanisms allowing DLC1 to impair migration, invasion, and metastasis formation.