Integrated methodologies to study human transcriptional regulation from functional genomics data

Abstract : The human body is composed of a huge number of cells acting together in a concerted manner. The current understanding is that proteins perform most of the necessary activities in keeping a cell alive. The DNA, on the other hand, stores the information on how to produce the different proteins in the genome. Regulating gene transcription is the first important step that can thus affect the life of a cell, modify its functions and its responses to the environment. Regulation is a complex operation that involves specialized proteins, the transcription factors. Transcription factors (TFs) can bind to DNA and activate the processes leading to the expression of genes into new proteins. Errors in this process may lead to diseases. In particular, some transcription factors have been associated with a lethal pathological state, commonly known as cancer, associated with uncontrolled cellular proliferation, invasiveness of healthy tissues and abnormal responses to stimuli. Understanding cancer-related regulatory programs is a difficult task, often involving several TFs interacting together and influencing each other's activity. This Thesis presents new computational methodologies to study gene regulation. In addition we present applications of our methods to the understanding of cancer-related regulatory programs. The understanding of transcriptional regulation is a major challenge. We address this difficult question combining computational approaches with large collections of heterogeneous experimental data. In detail, we design signal processing tools to recover transcription factors binding sites on the DNA from genome-wide surveys like chromatin immunoprecipitation assays on tiling arrays (ChIP-chip). We then use the localization about the binding of TFs to explain expression levels of regulated genes. In this way we identify a regulatory synergy between two TFs, the oncogene C-MYC and SP1. C-MYC and SP1 bind preferentially at promoters and when SP1 binds next to C-NIYC on the DNA, the nearby gene is strongly expressed. The association between the two TFs at promoters is reflected by the binding sites conservation across mammals, by the permissive underlying chromatin states 'it represents an important control mechanism involved in cellular proliferation, thereby involved in cancer. Secondly, we identify the characteristics of TF estrogen receptor alpha (hERa) target genes and we study the influence of hERa in regulating transcription. hERa, upon hormone estrogen signaling, binds to DNA to regulate transcription of its targets in concert with its co-factors. To overcome the scarce experimental data about the binding sites of other TFs that may interact with hERa, we conduct in silico analysis of the sequences underlying the ChIP sites using the collection of position weight matrices (PWMs) of hERa partners, TFs FOXA1 and SP1. We combine ChIP-chip and ChIP-paired-end-diTags (ChIP-pet) data about hERa binding on DNA with the sequence information to explain gene expression levels in a large collection of cancer tissue samples and also on studies about the response of cells to estrogen. We confirm that hERa binding sites are distributed anywhere on the genome. However, we distinguish between binding sites near promoters and binding sites along the transcripts. The first group shows weak binding of hERa and high occurrence of SP1 motifs, in particular near estrogen responsive genes. The second group shows strong binding of hERa and significant correlation between the number of binding sites along a gene and the strength of gene induction in presence of estrogen. Some binding sites of the second group also show presence of FOXA1, but the role of this TF still needs to be investigated. Different mechanisms have been proposed to explain hERa-mediated induction of gene expression. Our work supports the model of hERa activating gene expression from distal binding sites by interacting with promoter bound TFs, like SP1. hERa has been associated with survival rates of breast cancer patients, though explanatory models are still incomplete: this result is important to better understand how hERa can control gene expression. Thirdly, we address the difficult question of regulatory network inference. We tackle this problem analyzing time-series of biological measurements such as quantification of mRNA levels or protein concentrations. Our approach uses the well-established penalized linear regression models where we impose sparseness on the connectivity of the regulatory network. We extend this method enforcing the coherence of the regulatory dependencies: a TF must coherently behave as an activator, or a repressor on all its targets. This requirement is implemented as constraints on the signs of the regressed coefficients in the penalized linear regression model. Our approach is better at reconstructing meaningful biological networks than previous methods based on penalized regression. The method is tested on the DREAM2 challenge of reconstructing a five-genes/TFs regulatory network obtaining the best performance in the "undirected signed excitatory" category. Thus, these bioinformatics methods, which are reliable, interpretable and fast enough to cover large biological dataset, have enabled us to better understand gene regulation in humans.

Resting-State Functional Connectivity Emerges from Structurally and Dynamically Shaped Slow Linear Fluctuations.

Brain fluctuations at rest are not random but are structured in spatial patterns of correlated activity across different brain areas. The question of how resting-state functional connectivity (FC) emerges from the brain's anatomical connections has motivated several experimental and computational studies to understand structure-function relationships. However, the mechanistic origin of resting state is obscured by large-scale models' complexity, and a close structure-function relation is still an open problem. Thus, a realistic but simple enough description of relevant brain dynamics is needed. Here, we derived a dynamic mean field model that consistently summarizes the realistic dynamics of a detailed spiking and conductance-based synaptic large-scale network, in which connectivity is constrained by diffusion imaging data from human subjects. The dynamic mean field approximates the ensemble dynamics, whose temporal evolution is dominated by the longest time scale of the system. With this reduction, we demonstrated that FC emerges as structured linear fluctuations around a stable low firing activity state close to destabilization. Moreover, the model can be further and crucially simplified into a set of motion equations for statistical moments, providing a direct analytical link between anatomical structure, neural network dynamics, and FC. Our study suggests that FC arises from noise propagation and dynamical slowing down of fluctuations in an anatomically constrained dynamical system. Altogether, the reduction from spiking models to statistical moments presented here provides a new framework to explicitly understand the building up of FC through neuronal dynamics underpinned by anatomical connections and to drive hypotheses in task-evoked studies and for clinical applications.

Ab initio cluster-model study of the on-top chemisorption of F and Cl on Si(111) and Ge(111) surfaces

The interaction of atomic F and Cl with Si4H9 and Ge4H9 cluster models has been studied by using ab initio pseudopotentials and basis sets of increasing complexity. The results show that the effect of d orbitals is important in order to reproduce the experimental findings. However, the use of polarization functions in the atoms which are directly involved in the chemisorption bond leads to results which are very close to those obtained using extended basis sets. The local nature of the chemisorption bond is also interpreted by means of a Mulliken population analysis. For F-Si4H9 and Cl-Si4H9 the present results are in good agreement with previous ab initio all-electron calculations, and for the chemisorption of Cl on Si(111) and Ge(111) surfaces, good agreement is found with respect to the available experimental results as well as with previous slab calculations based on the local-density-functional formalism.

Effective t-J model Hamiltonian parameters of monolayered cuprate superconductors from ab initio electronic structure calculations

The magnetic coupling constant of selected cuprate superconductor parent compounds has been determined by means of embedded cluster model and periodic calculations carried out at the same level of theory. The agreement between both approaches validates the cluster model. This model is subsequently employed in state-of-the-art configuration interaction calculations aimed to obtain accurate values of the magnetic coupling constant and hopping integral for a series of superconducting cuprates. Likewise, a systematic study of the performance of different ab initio explicitly correlated wave function methods and of several density functional approaches is presented. The accurate determination of the parameters of the t-J Hamiltonian has several consequences. First, it suggests that the appearance of high-Tc superconductivity in existing monolayered cuprates occurs with J/t in the 0.20¿0.35 regime. Second, J/t=0.20 is predicted to be the threshold for the existence of superconductivity and, third, a simple and accurate relationship between the critical temperatures at optimum doping and these parameters is found. However, this quantitative electronic structure versus Tc relationship is only found when both J and t are obtained at the most accurate level of theory.

Improving plant functional groups for dynamic models of biodiversity: at the crossroads between functional and community ecology

The pace of on-going climate change calls for reliable plant biodiversity scenarios. Traditional dynamic vegetation models use plant functional types that are summarized to such an extent that they become meaningless for biodiversity scenarios. Hybrid dynamic vegetation models of intermediate complexity (hybrid-DVMs) have recently been developed to address this issue. These models, at the crossroads between phenomenological and process-based models, are able to involve an intermediate number of well-chosen plant functional groups (PFGs). The challenge is to build meaningful PFGs that are representative of plant biodiversity, and consistent with the parameters and processes of hybrid-DVMs. Here, we propose and test a framework based on few selected traits to define a limited number of PFGs, which are both representative of the diversity (functional and taxonomic) of the flora in the Ecrins National Park, and adapted to hybrid-DVMs. This new classification scheme, together with recent advances in vegetation modeling, constitutes a step forward for mechanistic biodiversity modeling.

A dual functional origin of transfer in the ICEclc genomic island of Pseudomonas knackmussii B13.

Genomic islands (GEIs) are large DNA segments, present in most bacterial genomes, that are most likely acquired via horizontal gene transfer. Here, we study the self-transfer system of the integrative and conjugative element ICEclc of Pseudomonas knackmussii B13, which stands model for a larger group of ICE/GEI with syntenic core gene organization. Functional screening revealed that unlike conjugative plasmids and other ICEs ICEclc carries two separate origins of transfer, with different sequence context but containing a similar repeat motif. Conjugation experiments with GFP-labelled ICEclc variants showed that both oriTs are used for transfer and with indistinguishable efficiencies, but that having two oriTs results in an estimated fourfold increase of ICEclc transfer rates in a population compared with having a single oriT. A gene for a relaxase essential for ICEclc transfer was also identified, but in vivo strand exchange assays suggested that the relaxase processes both oriTs in a different manner. This unique dual origin of transfer system might have provided an evolutionary advantage for distribution of ICE, a hypothesis that is supported by the fact that both oriT regions are conserved in several GEIs related to ICEclc.

Spinal cord T-cell infiltration in the spared nerve injury model of neuropathic pain: a time course study

Background : Numerous studies have shown that immune cells infiltrate the spinal cord after peripheral nerve injury and that they play a major contribution to sensory hypersensitivity in rodents. In particular, the role of monocyte-derived cells and T lymphocytes seems to be prominent in this process. This exciting new perspective in research on neuropathic pain opens many different areas of work, including the understanding of the function of these cells and how they impact on neural function. However, no systematic description of the time course or cell types that characterize this infiltration has been published yet, although this seems to be the rational first step of an overall understanding of the phenomenon. Objective : Describe the time course and cell characteristics of T lymphocyte infiltration in the spinal cord in the Spared Nerve Injury (SNI) model of neuropathic pain in rats. Methods : Collect of lumbar spinal cords of rats at days 2, 7, 21 and 40 after SNI or sham operation (n=4). Immunofluorescence detecting different proteins of T cell subgroups (CD2+CD4+, CD2+CD8+, Th1 markers, Th2 markers, Th17 markers). Quantification of the infiltration rate of the different subgroups. Expected results : First, we expect to see an infiltration of T cells in the spinal cord ipsilateral to nerve injury, higher in SNI rats than in sham animals. Second, we anticipate that different subtypes of T cells penetrate at different time points. Finally, the number of T lymphocytes are expected to decrease at the latest time point, showing a resolution of the process underlying their infiltrating the spinal cord in the first place. Impact : A systematic description of the infiltration of T cells in the spinal cord after peripheral nerve injury is needed to have a better understanding of the role of immune cells in neuropathic pain. The time course that we want to establish will provide the scientific community with new perspectives. First, it will confirm that T cells do indeed infiltrate the spinal cord after SNI in rats. Second, the type of T cells infiltrating at different time points will give clues about their function, in particular their inflammatory or anti-inflammatory profile. From there on, other studies could be lead, investigating the functional side of the specific subtypes put to light by us. Ultimately, this could lead to the discovery of new drugs targeting T cells or their infiltration, in the hope of improving neuropathic pain.

Why is the cladoceran Simocephalus vetulus (Müller) not a 'bang-bang strategist '? A critique of the optimal-body-size model

Lynch's (1980a) optimal-body-size model is designed to explain some major trends in cladoceran life histories; in particular the fact that large and littoral species seem to be bang-bang strategists (they grow first and the reproduce) whereas smaller planktonic species seem to be intermediate strategists (they grow and reproduce simultaneously). Predation is assumed to be an important selective pressure for these trends. Simocephalus vetulus (Müller) does not fit this pattern; being a littoral and relatively large species but an intermediate strategist. As shown by computer simulations, this species would reduce its per capita rate of increase by adopting the strategy predicted by the optimal-body-size model. Two aspects of the model are criticized: (1) the optimization criterion is shown to be incorrect and (2) the prediction of an intermediate strategy is not justified. Structural constraints are suggested to be responsible for the intermediate strategy of S.vetulus. Biotic interactions seem to have little effect on the observed life-history patterns of this species.

Functional homogenization of bumblebee communities in alpine landscapes under projected climate change

Background: Bumblebees represent an active pollinator group in mountain regions and assure the pollination of many different plant species from low to high elevations. Plant-pollinator interactions are mediated by functional traits. Shift in bumblebee functional structure under climate change may impact plant-pollinator interactions in mountains. Here, we estimated bumblebee upward shift in elevation, community turnover, and change in functional structure under climate change. Method: We sampled bumblebee species at 149 sites along the elevation gradient. We used stacked species distribution models (S-SDMs) forecasted under three climate change scenarios (A2, A1B, RCP3PD) to model the potential distribution of the Bombus species. Furthermore, we used species proboscis length measurements to assess the functional change in bumblebee assemblages along the elevation gradient. Results: We found species-specific response of bumblebee species to climate change. Species differed in their predicted rate of range contraction and expansion. Losers were mainly species currently restricted to high elevation. Under the most severe climate change scenarios (A2), we found a homogenization of proboscis length structure in bumblebee communities along the elevation gradient through the upward colonization of high elevation by species with longer proboscides. Conclusions: Here, we show that in addition to causing the shift in the distribution of bumblebee species, climate change may impact the functional structure of communities. The colonization of high elevation areas by bumblebee species with long proboscides may modify the structure of plant-pollination interaction networks by increasing the diversity of pollination services at high elevation.

Critical single proximal left arterial descending coronary artery stenosis to mimic chronic myocardial ischemia: a new model induced by minimal invasive technology.

BACKGROUND/AIMS: The present report examines a new pig model for progressive induction of high-grade stenosis, for the study of chronic myocardial ischemia and the dynamics of collateral vessel growth. METHODS: Thirty-nine Landrace pigs were instrumented with a novel experimental stent (GVD stent) in the left anterior descending coronary artery. Eight animals underwent transthoracic echocardiography at rest and under low-dose dobutamine. Seven animals were examined by nuclear PET and SPECT analysis. Epi-, mid- and endocardial fibrosis and the numbers of arterial vessels were examined by histology. RESULTS: Functional analysis showed a significant decrease in global left ventricular ejection fraction (24.5 +/- 1.6%) 3 weeks after implantation. There was a trend to increased left ventricular ejection fraction after low-dose dobutamine stress (36.0 +/- 6.6%) and a significant improvement of the impaired regional anterior wall motion. PET and SPECT imaging documented chronic hibernation. Myocardial fibrosis increased significantly in the ischemic area with a gradient from epi- to endocardial. The number of arterial vessels in the ischemic area increased and coronary angiography showed abundant collateral vessels of Rentrop class 1. CONCLUSION: The presented experimental model mimics the clinical situation of chronic myocardial ischemia secondary to 1-vessel coronary disease.

Deleterious Actions Of Vegf On The Blood Retinal Barrier And Photoreceptor Survival In The Light-damage Model

Purpose: The retinal balance between pro- and anti-angiogenic factors is critical for angiogenesis control, but is also involved in cell survival. We previously reported upregulation of VEGF and photoreceptor (PR) cell death in the Light-damage (LD) model. Preliminary results showed that anti-VEGF can rescue PR from cell death. Thus, we investigated the role of VEGF on the retina and we herein described the effect of anti-VEGF antibody delivered by lentiviral gene transfer in this model.Methods: To characterize the action of VEGF during the LD, we exposed Balb/c mice subretinally injected with LV-anti-VEGF, or not, to 5'000 lux for 1h. We next evaluated the retinal function, PR survival and protein expression (VEGF, VEGFR1/2, Src, PEDF, p38MAPK, Akt, Peripherin, SWL-opsin) after LD. We analyzed Blood retinal barrier (BRB) integrity on flat-mounted RPE and cryosections stained with β-catenin, ZO-1, N-cadherin and albumin.Results: Results indicate that the VEGF pathway is modulated after LD. LD leads to extravascular albumin leakage and BRB breakdown: β-catenin, ZO-1 and N-cadherin translocate to the cytoplasm of RPE cells showing loss of cell cohesion. This phenomenon is in adequacy with the VEGF time-course expression. Assessment of the retinal function reveals that PR rescue correlates with the level of LV-anti-VEGF expression. Rhodopsin content was higher in the LV-anti-VEGF group than in controls and measures of the ONL thickness indicate that LV-anti-VEGF preserves by 82% the outer nuclear layer from degeneration. Outer segments (OS) appeared well organized with an appropriate length in the LV-anti-VEGF group compared to controls, and the expression of SWL-opsin is maintained in the OS without being mislocalized as in the LV-GFP group. Finally, LV-anti-VEGF treatment prevents BRB breakdown and maintained RPE cell integrity.Conclusions: This study involves VEGF in LD and highlights the prime importance of the BRB integrity for PR survival. Taken together, these results show that anti-VEGF is neuroprotective in this model and maintains functional PR layer in LD-treated mice.

Evaluating fear of falling among community-dwelling older people: A qualitative pattern model based on the experience of falls

This research was conducted in the context of the project IRIS 8A Health and Society (2002-2008) and financially supported by the University of Lausanne. It was aomed at developping a model based on the elder people's experience and allowed us to develop a "Portrait evaluation" of fear of falling using their examples and words. It is a very simple evaluation, which can be used by professionals, but by the elder people themselves. The "Portrait evaluation" and the user's guide are on free access, but we would very much approciate to know whether other people or scientists have used it and collect their comments. (contact: Chantal.Piot-Ziegler@unil.ch)The purpose of this study is to create a model grounded in the elderly people's experience allowing the development of an original instrument to evaluate FOF.In a previous study, 58 semi-structured interviews were conducted with community-dwelling elderly people. The qualitative thematic analysis showed that fear of falling was defined through the functional, social and psychological long-term consequences of falls (Piot-Ziegler et al., 2007).In order to reveal patterns in the expression of fear of falling, an original qualitative thematic pattern analysis (QUAlitative Pattern Analysis - QUAPA) is developed and applied on these interviews.The results of this analysis show an internal coherence across the three dimensions (functional, social and psychological). Four different patterns are found, corresponding to four degrees of fear of falling. They are formalized in a fear of falling intensity model.This model leads to a portrait-evaluation for fallers and non-fallers. The evaluation must be confronted to large samples of elderly people, living in different environments. It presents an original alternative to the concept of self-efficacy to evaluate fear of falling in older people.The model of FOF presented in this article is grounded on elderly people's experience. It gives an experiential description of the three dimensions constitutive of FOF and of their evolution as fear increases, and defines an evaluation tool using situations and wordings based on the elderly people's discourse.

An integer-based score to predict functional outcome in acute ischemic stroke: The ASTRAL score.

OBJECTIVE: To develop and validate a simple, integer-based score to predict functional outcome in acute ischemic stroke (AIS) using variables readily available after emergency room admission. METHODS: Logistic regression was performed in the derivation cohort of previously independent patients with AIS (Acute Stroke Registry and Analysis of Lausanne [ASTRAL]) to identify predictors of unfavorable outcome (3-month modified Rankin Scale score >2). An integer-based point-scoring system for each covariate of the fitted multivariate model was generated by their β-coefficients; the overall score was calculated as the sum of the weighted scores. The model was validated internally using a 2-fold cross-validation technique and externally in 2 independent cohorts (Athens and Vienna Stroke Registries). RESULTS: Age (A), severity of stroke (S) measured by admission NIH Stroke Scale score, stroke onset to admission time (T), range of visual fields (R), acute glucose (A), and level of consciousness (L) were identified as independent predictors of unfavorable outcome in 1,645 patients in ASTRAL. Their β-coefficients were multiplied by 4 and rounded to the closest integer to generate the score. The area under the receiver operating characteristic curve (AUC) of the score in the ASTRAL cohort was 0.850. The score was well calibrated in the derivation (p = 0.43) and validation cohorts (0.22 [Athens, n = 1,659] and 0.49 [Vienna, n = 653]). AUCs were 0.937 (Athens), 0.771 (Vienna), and 0.902 (when pooled). An ASTRAL score of 31 indicates a 50% likelihood of unfavorable outcome. CONCLUSIONS: The ASTRAL score is a simple integer-based score to predict functional outcome using 6 readily available items at hospital admission. It performed well in double external validation and may be a useful tool for clinical practice and stroke research.

Electric dipole polarizability in 208Pb: insights from the droplet model

We study the electric dipole polarizability α D in 208 Pb based on the predictions of a large and representative set of relativistic and nonrelativistic nuclear mean-field models. We adopt the droplet model as a guide to better understand the correlations between α D and other isovector observables. Insights from the droplet model suggest that the product of α D and the nuclear symmetry energy at saturation density J is much better correlated with the neutron skin thickness r np of 208 Pb than the polarizability alone. Correlations of α D J with r np and with the symmetry energy slope parameter L suggest that α D J is a strong isovector indicator. Hence, we explore the possibility of constraining the isovector sector of the nuclear energy density functional by comparing our theoretical predictions against measurements of both α D and the parity-violating asymmetry in 208 Pb. We find that the recent experimental determination of α D in 208 Pb in combination with the range for the symmetry energy at saturation density J = [31 ± (2) est] MeV suggests r np (208 Pb) = 0 . 165 ± (0 . 009) expt ± (0 . 013) theor ± (0.021) est fm and L = 43 ± (6) expt ± (8) theor ± (12) est MeV

A dynamic model of keratinocyte stem cell renewal and differentiation: role of the p21WAF1/Cip1 and Notch1 signaling pathways.

We present here a dynamic model of functional equilibrium between keratinocyte stem cells, transit amplifying populations and cells that are reversibly versus irreversibly committed to differentiation. According to this model, the size of keratinocyte stem cell populations can be controlled at multiple levels, including relative late steps in the sequence of events leading to terminal differentiation and by the influences of a heterogeneous extra-cellular environment. We discuss how work in our laboratory, on the interconnection between the cyclin/CDK inhibitor p21WAF1/Cip1 and the Notch1 signaling pathways, provides strong support to this dynamic model of stem cell versus committed and/or differentiated keratinocyte populations.