Scalable inference of heterogeneous reaction kinetics from pooled single-cell recordings.

Mathematical methods combined with measurements of single-cell dynamics provide a means to reconstruct intracellular processes that are only partly or indirectly accessible experimentally. To obtain reliable reconstructions, the pooling of measurements from several cells of a clonal population is mandatory. However, cell-to-cell variability originating from diverse sources poses computational challenges for such process reconstruction. We introduce a scalable Bayesian inference framework that properly accounts for population heterogeneity. The method allows inference of inaccessible molecular states and kinetic parameters; computation of Bayes factors for model selection; and dissection of intrinsic, extrinsic and technical noise. We show how additional single-cell readouts such as morphological features can be included in the analysis. We use the method to reconstruct the expression dynamics of a gene under an inducible promoter in yeast from time-lapse microscopy data.

Restoration of connexin 40 (cx40) in Renin-producing cells reduces the hypertension of cx40 null mice.

Connexin 40 (Cx40) is expressed by the renin-producing cells (RSCs) of the kidneys and the endothelial cells of blood vessels. Cx40 null mice (Cx40(-/-)) feature a much increased renin synthesis and secretion, which results in chronic hypertension, and also display an altered endothelium-dependent relaxation of the aorta because of reduced eNOS levels and nitric oxide production. To discriminate the effect of Cx40 in renin secretion and vascular signaling, we targeted Cx40 to either the RSCs or the endothelial cells of Cx40 null mice. When compared with Cx40(-/-) controls, the animals expressing Cx40 in RSCs were less hypertensive and featured reduced renin levels, still numerous RSCs outside the wall of the afferent arterioles. In contrast, mice expressing Cx40 in the endothelial cells were as hypertensive as Cx40(-/-) mice, in spite of control levels of Cx37 and eNOS. Our data show that blood pressure is improved by restoration of Cx40 expression in RSCs but not in endothelial cells, stressing the prominent role of renin in the mouse hypertension linked to loss of Cx40.

Heterogeneous 40Ar* distributions in naturally deformed muscovite: in situ UV-laser ablation evidence for micro structurally controlled intragrain diffusion

Micas are commonly used in Ar-40/Ar-39 thermochronological studies of variably deformed rocks yet the physical basis by which deformation may affect radiogenic argon retention in mica is poorly constrained. This study examines the relationship between deformation and deformation-induced microstructures on radiogenic argon retention in muscovite, A combination of furnace step-heating and high-spatial resolution in situ UV-laser ablation Ar-40/Ar-39 analyses are reported for deformed muscovites sampled from a granitic pegmatite vein within the Siviez-Mischabel Nappe, western Swiss Alps (Penninic domain, Brianconnais unit). The pegmatite forms part of the Variscan (similar to 350 Ma) Alpine basement and exhibits a prominent Alpine S-C fabric including numerous mica `fish' that developed under greenschist facies metamorphic conditions, during the dominant Tertiary Alpine tectonic phase of nappe emplacement. Furnace step-heating of milligram quantities of separated muscovite grains yields an Ar-40/Ar-39 age spectrum with two distinct staircase segments but without any statistical plateau, consistent with a previous study from the same area. A single (3 X 5 mm) muscovite porphyroclast (fish) was investigated by in situ UV-laser ablation. A histogram plot of 170 individual Ar-40/Ar-39 UV-laser ablation ages exhibit a range from 115 to 387 Ma with modes at approximately 340 and 260 Ma. A variogram statistical treatment of the (40)Ad/Ar-39 results reveals ages correlated with two directions; a highly correlated direction at 310 degrees and a lesser correlation at 0 degrees relative to the sense of shearing. Using the highly correlated direction a statistically generated (Kriging method) age contour map of the Ar-40/Ar-39 data reveals a series of elongated contours subparallel to the C-surfaces which where formed during Tertiary nappe emplacement. Similar data distributions and slightly younger apparent ages are recognized in a smaller mica fish. The observed intragrain age variations are interpreted to reflect the partial loss of radiogenic argon during Alpine (similar to 35 Ma) greenschist facies metamorphism. One-dirnensional diffusion modelling results are consistent with the idea that the zones of youngest apparent age represent incipient shear band development within the mica porphyroclasts, thus providing a network of fast diffusion pathways. During Alpine greenschist facies metamorphism the incipient shear bands enhanced the intragrain loss of radiogenic argon. The structurally controlled intragrain age variations observed in this investigation imply that deformation has a direct control on the effective length scale for argon diffusion, which is consistent with the heterogeneous nature of deformation. (C) 2001 Elsevier Science B.V. All rights reserved.

Two functional forms of the Xenopus laevis estrogen receptor translated from a single mRNA species.

Steroid receptors are nuclear proteins that regulate gene transcription in a ligand-dependent manner. Over-expression of the Xenopus estrogen receptor in a vaccinia virus-derived expression system revealed that the receptor localized exclusively in the nucleus of the infected cells, irrespective of the presence or absence of the ligand. Furthermore, two forms of the receptor were produced, a full-length and a N-terminal truncated version, which are translated from a single mRNA species by the use of two AUG within the same reading frame. These 66- and 61-kDa receptors were also observed after in vitro translation of the mRNA as well as in primary Xenopus hepatocytes. Both forms are potent estrogen-dependent transcriptional activators in transient transfection experiments, as well as in in vitro transcription assays.

Beyond Smith's rule: An optimal dynamic index, rule for single machine stochastic scheduling with convex holding costs

Most research on single machine scheduling has assumedthe linearity of job holding costs, which is arguablynot appropriate in some applications. This motivates ourstudy of a model for scheduling $n$ classes of stochasticjobs on a single machine, with the objective of minimizingthe total expected holding cost (discounted or undiscounted). We allow general holding cost rates that are separable,nondecreasing and convex on the number of jobs in eachclass. We formulate the problem as a linear program overa certain greedoid polytope, and establish that it issolved optimally by a dynamic (priority) index rule,whichextends the classical Smith's rule (1956) for the linearcase. Unlike Smith's indices, defined for each class, ournew indices are defined for each extended class, consistingof a class and a number of jobs in that class, and yieldan optimal dynamic index rule: work at each time on a jobwhose current extended class has larger index. We furthershow that the indices possess a decomposition property,as they are computed separately for each class, andinterpret them in economic terms as marginal expected cost rate reductions per unit of expected processing time.We establish the results by deploying a methodology recentlyintroduced by us [J. Niño-Mora (1999). "Restless bandits,partial conservation laws, and indexability. "Forthcomingin Advances in Applied Probability Vol. 33 No. 1, 2001],based on the satisfaction by performance measures of partialconservation laws (PCL) (which extend the generalizedconservation laws of Bertsimas and Niño-Mora (1996)):PCL provide a polyhedral framework for establishing theoptimality of index policies with special structure inscheduling problems under admissible objectives, which weapply to the model of concern.

On the term structure of Interbank interest rates: Jump-diffusion processes and option pricing

In this paper we study the dynamic behavior of the term structureof Interbank interest rates and the pricing of options on interest ratesensitive securities. We posit a generalized single factor model withjumps to take into account external influences in the market. Daily datais used to test for jump effects. Qualitative examination of the linkagebetween Monetary Authorities' interventions and jumps are studied. Pricingresults suggests a systematic underpricing in bonds and call options ifthe jumps component is not included. However, the pricing of put optionson bonds presents indeterminacies.

EWS-FLI-1 modulates miRNA145 and SOX2 expression to initiate mesenchymal stem cell reprogramming toward Ewing sarcoma cancer stem cells.

Cancer stem cells (CSCs) display plasticity and self-renewal properties reminiscent of normal tissue stem cells, but the events responsible for their emergence remain obscure. We recently identified CSCs in Ewing sarcoma family tumors (ESFTs) and showed that they retain mesenchymal stem cell (MSC) plasticity. In the present study, we addressed the mechanisms that underlie ESFT CSC development. We show that the EWS-FLI-1 fusion gene, associated with 85%-90% of ESFTs and believed to initiate their pathogenesis, induces expression of the embryonic stem cell (ESC) genes OCT4, SOX2, and NANOG in human pediatric MSCs (hpMSCs) but not in their adult counterparts. Moreover, under appropriate culture conditions, hpMSCs expressing EWS-FLI-1 generate a cell subpopulation displaying ESFT CSC features in vitro. We further demonstrate that induction of the ESFT CSC phenotype is the result of the combined effect of EWS-FLI-1 on its target gene expression and repression of microRNA-145 (miRNA145) promoter activity. Finally, we provide evidence that EWS-FLI-1 and miRNA-145 function in a mutually repressive feedback loop and identify their common target gene, SOX2, in addition to miRNA145 itself, as key players in ESFT cell differentiation and tumorigenicity. Our observations provide insight for the first time into the mechanisms whereby a single oncogene can reprogram primary cells to display a CSC phenotype.

Increased mobility of major histocompatibility complex I-peptide complexes decreases the sensitivity of antigen recognition.

CD8(+) cytotoxic T lymphocytes (CTL) can recognize and kill target cells expressing only a few cognate major histocompatibility complex (MHC) I-peptide complexes. This high sensitivity requires efficient scanning of a vast number of highly diverse MHC I-peptide complexes by the T cell receptor in the contact site of transient conjugates formed mainly by nonspecific interactions of ICAM-1 and LFA-1. Tracking of single H-2K(d) molecules loaded with fluorescent peptides on target cells and nascent conjugates with CTL showed dynamic transitions between states of free diffusion and immobility. The immobilizations were explained by association of MHC I-peptide complexes with ICAM-1 and strongly increased their local concentration in cell adhesion sites and hence their scanning by T cell receptor. In nascent immunological synapses cognate complexes became immobile, whereas noncognate ones diffused out again. Interfering with this mobility modulation-based concentration and sorting of MHC I-peptide complexes strongly impaired the sensitivity of antigen recognition by CTL, demonstrating that it constitutes a new basic aspect of antigen presentation by MHC I molecules.

Contrasting the functional properties of GABAergic axon terminals with single and multiple synapses in the thalamus

Diverse sources of GABAergic inhibition are a major feature of cortical networks, but distinct inhibitory input systems have not been systematically characterized in the thalamus. Here, we contrasted the properties of two independent GABAergic pathways in the posterior thalamic nucleus of rat, one input from the reticular thalamic nucleus (nRT), and one "extrareticular" input from the anterior pretectal nucleus (APT). The vast majority of nRT-thalamic terminals formed single synapses per postsynaptic target and innervated thin distal dendrites of relay cells. In contrast, single APT-thalamic terminals formed synaptic contacts exclusively via multiple, closely spaced synapses on thick relay cell dendrites. Quantal analysis demonstrated that the two inputs displayed comparable quantal amplitudes, release probabilities, and multiple release sites. The morphological and physiological data together indicated multiple, single-site contacts for nRT and multisite contacts for APT axons. The contrasting synaptic arrangements of the two pathways were paralleled by different short-term plasticities. The multisite APT-thalamic pathway showed larger charge transfer during 50-100 Hz stimulation compared with the nRT pathway and a greater persistent inhibition accruing during stimulation trains. Our results demonstrate that the two inhibitory systems are morpho-functionally distinct and suggest and that multisite GABAergic terminals are tailored for maintained synaptic inhibition even at high presynaptic firing rates. These data explain the efficacy of extrareticular inhibition in timing relay cell activity in sensory and motor thalamic nuclei. Finally, based on the classic nomenclature and the difference between reticular and extrareticular terminals, we define a novel, multisite GABAergic terminal type (F3) in the thalamus.

Optimization of multiclass queueing networks with changeover times via the achievable region approach: Part I, the single-station case

We address the performance optimization problem in a single-stationmulticlass queueing network with changeover times by means of theachievable region approach. This approach seeks to obtainperformance bounds and scheduling policies from the solution of amathematical program over a relaxation of the system's performanceregion. Relaxed formulations (including linear, convex, nonconvexand positive semidefinite constraints) of this region are developedby formulating equilibrium relations satisfied by the system, withthe help of Palm calculus. Our contributions include: (1) newconstraints formulating equilibrium relations on server dynamics;(2) a flow conservation interpretation of the constraintspreviously derived by the potential function method; (3) newpositive semidefinite constraints; (4) new work decomposition lawsfor single-station multiclass queueing networks, which yield newconvex constraints; (5) a unified buffer occupancy method ofperformance analysis obtained from the constraints; (6) heuristicscheduling policies from the solution of the relaxations.

Functional analysis of human POT1 and RPA : insights into the role of single stranded DNA binding proteins at telomeres

Abstract Telomeres, the natural ends of chromosomes, need to be protected from chromosome end fusions, aberrant homologous recombination and degradation. In humans, chromosome ends are specified through arrays of tandemly repeated 5'-TTAGGG-3' hexamers, ending in a 3' overhang. A complex formed by the six proteins TRF1, TRF2, hRap1, TIN2, TPP1 and POT1 specifically assocìates with and protects telomeres. Telomeres are maintained by semiconservative DNA replication and by a specialized reverse transcriptase, telomerase, that carries an RNA subunit which templates new telomeric repeat synthesis. The telomeric single stranded (ss) DNA binding protein POT1 protects the telomeric 3' overhang and modulates telomerase-mediated telomere elongation. It is possible that POT1 also influences DNA synthesis during semiconservative DNA replication, which is initiated by the DNA polymerase alpha-primase complex. The heterotrimeric ss DNA-binding protein RPA plays essential roles during DNA replication. RPA binds to ss DNA with high affinity in order to stabilize ss DNA and facilitate nascent strand synthesis at the replication fork. Here we investigate how the two proteins RPA and POT1 contribute to telomere maintenance by regulating semi-conservative DNA replication and telomerase. Using chromatin immunoprecipitation experiments, we show that RPA associates with telomeres during S-phase. Analysis of telomere structure in cells shRNA-depleted for RPA and POT1 reveals that loss of RPA and POT1 causes exposure of single-stranded DNA at telomeres, suggestive of incomplete DNA replication. Biochemical experiments using purified recombinant POT1 and RPA show that saturating telomeric oligonucleotides with POT1 or RPA reduces the primase activity of the DNA polymerase alpha-primase complex and the overall activity of telomerase. POT1 and RPA also increase the primer extension by DNA polymerase alpha-primase complex and the processivity of telomerase under certain conditions, although POT1 increases the activities to a greater extent than RPA. We propose that POT1 is required for proper replication of the lagging strand of telomeres and that some phenotypes observed in POT1-depleted cells may stern from incomplete DNA replication rather than de-protection of the single-stranded overhang. Résumé Les télomères, les extrémités normales des chromosomes linéaires, doivent être protégés des fusions chromosomiques, d'événements de recombinaison homologue aberrants et de phénomènes de dégradation. Chez l'Homme, les extrémités des chromosomes sont constitués d'ADN double brin répétitif de séquence 5'-TTAGGG-3', d'une extension simple brin 3' sortante et d'un complexe protéique formé des six facteurs TRF1, TRF2, hRap1, TIN2, TPP1 et POT1 qui, s'associant à cette séquence, protègent l'ADN télomèrique. Les télomères sont maintenus par la télomérase, une transcriptase inverse capable d'allonger l'extension 3' sortante télomérique. POT1 lie l'ADN simple brin télomérique et module l'élongation des télomères par la télomérase. POT1 pourrait en théorie également influencer la réplication semi-conservative de l'ADN. L'ADN-polymérase Pal alpha-primase amorce et initie la synthèse d'ADN. Pendant la réplication, l'ADN simple brin est stabilisé par RPA, un complexe hétérotrimèrique qui lie l'ADN simple brin. RPA facilite la synthèse du brin naissant à la fourche de réplication. Ici nous avons étudié comment ces deux protéines qui lient l'ADN simple brin, RPA et POT1, régulent la réplication des télomères par la télomérase et la machinerie classique de réplication de l'ADN. Par immunoprécipitation de chromatine (ChIP), nous montrons que RPA est localisé aux télomères lors de la phase S du cycle cellulaire. De plus, l'analyse de la structure des télomeres indique que !a perte de RPA ou de POT1 conduit à l'apparition d'ADN simple brin télomérique, suggérant une réplication incomplète de l'ADN télomérique in vivo. Par une approche complémentaire biochimique utilisant les protéines POT1 et RPA recombinantes purifiées, nous montrons également que la liaison de POT1 ou de RPA à des oligonucléotides télomériques bloque l'activité primase du complexe polymérase alpha/primase et réduit l'activité télomérase sur ces substrats. En revanche, leur liaison augmente l'activité ADN-polymérase du complexe polymérase alpha/primase, ainsi que fa processivité de la télomérase dans certaines conditions, POT1 étant le plus efficace des deux facteurs. Nous proposons que POT1 est nécessaire à la réplication du brin retardé au niveau des télomères, ce qui suggère que certains phénotypes des cellules déplétés en POT1 puissent résulter d'une réplication incomplète de l'ADN télémétrique plutôt que d'une déprotection de l'extrémité sortante des télomères.

The Bcl-2/Bcl-XL inhibitor ABT-737 promotes death of retinoblastoma cancer cells.

Purpose: Retinoblastoma is a malignant tumor that usually develops in early childhood. During retinoblastoma spreading, RB1 gene inactivation is followed by additional genomic modifications which progressively lead to resistance of tumor cells to death. Drugs that act at downstream levels of death signaling pathways should therefore be interesting in killing retinoblastoma cells. ABT-737, a BH3 mimetic molecule effective at the mitochondrial level, has been shown to induce apoptosis in different human tumoral cell lines as well as in primary patient-derived cells, and in a mouse xenograph model. Methods: In this report, we analyzed the pro-death effect of ABT-737 on two human retinoblastoma cell lines, Y79 and WERI-Rb, as well as on the mouse photoreceptor cell line 661W. Results: We observed that ABT-737 was very effective as a single agent in inducing human WERI-Rb cells apoptosis without affecting the mouse 661W photoreceptor cells. However human Y79 cells were resistant to ABT-737, as a probable consequence of the absence of Bax. The high sensitivity of WERI-Rb to ABT-737 can be increased by downregulating Mcl-1 using the proteasome inhibitor MG-132. Preliminary analysis in primary mouse retinoblastoma tumoral cell lines predicts high sensitivity to ABT-737. Conclusion: Our data suggest that ABT-737 or related compounds could be a highly effective drug in the treatment of some retinoblastomas.

Dissociation of thymic positive and negative selection in transgenic mice expressing major histocompatibility complex class I molecules exclusively on thymic cortical epithelial cells.

Thymic positive and negative selection of developing T lymphocytes confronts us with a paradox: How can a T-cell antigen receptor (TCR)-major histocompatibility complex (MHC)/peptide interaction in the former process lead to transduction of signals allowing for cell survival and in the latter induce programmed cell death or a hyporesponsive state known as anergy? One of the hypotheses put forward states that the outcome of a TCR-MHC/peptide interaction depends on the cell type presenting the selecting ligand to the developing thymocyte. Here we describe the development and lack of self-tolerance of CD8(+) T lymphocytes in transgenic mice expressing MHC class I molecules in the thymus exclusively on cortical epithelial cells. Despite the absence of MHC class I expression on professional antigen-presenting cells, normal numbers of CD8(+) cells were observed in the periphery. Upon specific activation, transgenic CD8(+) T cells efficiently lysed syngeneic MHC class I(+) targets in vitro and in vivo, indicating that thymic cortical epithelium (in contrast to medullary epithelium and antigen-presenting cells of hematopoietic origin) is incapable of tolerance induction. Thus, compartmentalization of the antigen-presenting cells involved in thymic positive selection and tolerance induction can (at least in part) explain the positive/negative selection paradox.

Large TCR diversity of virus-specific CD8 T cells provides the mechanistic basis for massive TCR renewal after antigen exposure.

Ex vivo analysis of virus-specific CD8 T cell populations by anchored PCR has shown that the CD8 TCR repertoire was less oligoclonal (seven to nine clonotypes per individual epitope) than previously thought. In the current study, TCR diversity was investigated by assessing both the overall TCR β-chain variable regions usage as well as the CDR3 regions in ex vivo-isolated CMV- and EBV-specific CD8 T cells from 27 healthy donors. The average number of clonotypes specific to most single viral epitopes comprised between 14 and 77. Changes in the CD8 TCR repertoire were also longitudinally assessed under conditions of HIV-1 chronic infection (i.e., in patients with suppressed virus replication and after treatment interruption and Ag re-exposure). The results showed that a large renewal (≤80%) of the TRB repertoire occurred after Ag re-exposure and was eventually associated with an increased T cell recognition functional avidity. These results demonstrate that the global CD8 TCR repertoire is much more diverse (≤9-fold) than previously estimated and provide the mechanistic basis for supporting massive repertoire renewal during chronic virus infection and Ag re-exposure.

Immunohistochemical and flow cytometric analysis of long-term label-retaining cells in the adult heart.

Cardiac-resident stem/progenitor cells have been identified based on expression of stem cell-associated antigens. However, no single surface marker allows to identify a definite cardiac stem/progenitor cell entity. Hence, functional stem cell markers have been extensively searched for. In homeostatic systems, stem cells divide infrequently and therefore retain DNA labels such as 5-bromo-2'-deoxyuridine, which are diluted with division. We used this method to analyze long-term label-retaining cells in the mouse heart after 14 days of 5-bromo-2'-deoxyuridine administration. Labeled cells were detected using immunohistochemical and flow-cytometric methods after varying chasing periods up to 12 months. Using mathematical models, the observed label dilution could consistently be described in the context of a 2-population model, whereby a population of rapidly dividing cells accounted for an accelerated early decline, and a population of slowly dividing cells accounted for decelerated dilution on longer time scales. Label-retaining cells were preferentially localized in the atria and apical region and stained negative for markers of the major cell lineages present in the heart. Most cells with long-term label-retention expressed stem cell antigen-1 (Sca-1). Sca-1(+)CD31(-) cells formed cell aggregates in culture, out of which lineage-negative (Lin(-))Sca-1(+)CD31(-) cells emerged, which could be cultured for many passages. These cells formed cardiospheres and showed differentiation potential into mesenchymal cell lineages. When cultured in cardiomyogenic differentiation medium, they expressed cardiac-specific genes. In conclusion, recognition of slow-cycling cells provides functional evidence of stem/progenitor cells in the heart. Lin(-)Sca-1(+)CD31(-) cardiac-derived progenitors have a potential for differentiation into cardiomyogenic and mesenchymal cell lineages.