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.

Improving the chance of cure of follicular lymphoma by combining immunotherapy and radioimmunotherapy based on anti-CD20 antibodies?

Right ventricular and left ventricular systolic time intervals (RVSTIs and LVSTIs) were measured in normal term and preterm infants from 1 hour to 90 days of life. LVSTIs in both term and preterm infants were similar in the first five days of life. The ratio of left pre-ejection period (LPEP) to left ventricular ejection time (LVET) was lower in preterm infants older than age 5 days. Estimated gestational age had no influence on LVSTI. The ratio of right pre-ejection period (RPEP) to right ventricular ejection time (RVET) was lower in preterm infants (0.32) than in term newborns (0.37). The preterm RPEP/RVET ratio decreased with age, but at a slower rate than in term babies. This was consistent with the lower pulmonary vascular resistance present in preterm infants.

A novel method of dynamic culture surface expansion improves mesenchymal stem cell proliferation and phenotype.

Repeated passaging in conventional cell culture reduces pluripotency and proliferation capacity of human mesenchymal stem cells (MSC). We introduce an innovative cell culture method whereby the culture surface is dynamically enlarged during cell proliferation. This approach maintains constantly high cell density while preventing contact inhibition of growth. A highly elastic culture surface was enlarged in steps of 5% over the course of a 20-day culture period to 800% of the initial surface area. Nine weeks of dynamic expansion culture produced 10-fold more MSC compared with conventional culture, with one-third the number of trypsin passages. After 9 weeks, MSC continued to proliferate under dynamic expansion but ceased to grow in conventional culture. Dynamic expansion culture fully retained the multipotent character of MSC, which could be induced to differentiate into adipogenic, chondrogenic, osteogenic, and myogenic lineages. Development of an undesired fibrogenic myofibroblast phenotype was suppressed. Hence, our novel method can rapidly provide the high number of autologous, multipotent, and nonfibrogenic MSC needed for successful regenerative medicine.

A multilevel multiscale finite volume method

The Multiscale Finite Volume (MsFV) method has been developed to efficiently solve reservoir-scale problems while conserving fine-scale details. The method employs two grid levels: a fine grid and a coarse grid. The latter is used to calculate a coarse solution to the original problem, which is interpolated to the fine mesh. The coarse system is constructed from the fine-scale problem using restriction and prolongation operators that are obtained by introducing appropriate localization assumptions. Through a successive reconstruction step, the MsFV method is able to provide an approximate, but fully conservative fine-scale velocity field. For very large problems (e.g. one billion cell model), a two-level algorithm can remain computational expensive. Depending on the upscaling factor, the computational expense comes either from the costs associated with the solution of the coarse problem or from the construction of the local interpolators (basis functions). To ensure numerical efficiency in the former case, the MsFV concept can be reapplied to the coarse problem, leading to a new, coarser level of discretization. One challenge in the use of a multilevel MsFV technique is to find an efficient reconstruction step to obtain a conservative fine-scale velocity field. In this work, we introduce a three-level Multiscale Finite Volume method (MlMsFV) and give a detailed description of the reconstruction step. Complexity analyses of the original MsFV method and the new MlMsFV method are discussed, and their performances in terms of accuracy and efficiency are compared.

Development of a two-step screening ESI-TOF-MS method for rapid determination of significant stress-induced metabolome modifications in plant leaf extracts: the wound response in Arabidopsis thaliana as a case study.

To study the stress-induced effects caused by wounding under a new perspective, a metabolomic strategy based on HPLC-MS has been devised for the model plant Arabidopsis thaliana. To detect induced metabolites and precisely localise these compounds among the numerous constitutive metabolites, HPLC-MS analyses were performed in a two-step strategy. In a first step, rapid direct TOF-MS measurements of the crude leaf extract were performed with a ballistic gradient on a short LC-column. The HPLC-MS data were investigated by multivariate analysis as total mass spectra (TMS). Principal components analysis (PCA) and hierarchical cluster analysis (HCA) on principal coordinates were combined for data treatment. PCA and HCA demonstrated a clear clustering of plant specimens selecting the highest discriminating ions given by the complete data analysis, leading to the specific detection of discrete-induced ions (m/z values). Furthermore, pool constitution with plants of homogeneous behaviour was achieved for confirmatory analysis. In this second step, long high-resolution LC profilings on an UPLC-TOF-MS system were used on pooled samples. This allowed to precisely localise the putative biological marker induced by wounding and by specific extraction of accurate m/z values detected in the screening procedure with the TMS spectra.

Clinal variation in a brown lemur (Eulemur spp.) hybrid zone: Combining morphological, genetic and climatic data to examine stability.

Studies of hybrid zones can inform our understanding of reproductive isolation and speciation. Two species of brown lemur (Eulemur rufifrons and E. cinereiceps) form an apparently stable hybrid zone in the Andringitra region of south-eastern Madagascar. The aim of this study was to identify factors that contribute to this stability. We sampled animals at 11 sites along a 90-km transect through the hybrid zone and examined variation in 26 microsatellites, the D-loop region of mitochondrial DNA, six pelage and nine morphological traits; we also included samples collected in more distant allopatric sites. Clines in these traits were noncoincident, and there was no increase in either inbreeding coefficients or linkage disequilibrium at the centre of the zone. These results could suggest that the hybrid zone is maintained by weak selection against hybrids, conforming to either the tension zone or geographical selection-gradient model. However, a closer examination of clines in pelage and microsatellites indicates that these clines are not sigmoid or stepped in shape but instead plateau at their centre. Sites within the hybrid zone also occur in a distinct habitat, characterized by greater seasonality in precipitation and lower seasonality in temperature. Together, these findings suggest that the hybrid zone may follow the bounded superiority model, with exogenous selection favouring hybrids within the transitional zone. These findings are noteworthy, as examples supporting the bounded superiority model are rare and may indicate a process of ecologically driven speciation without geographical isolation.

Identification and validation of copy number variants using SNP genotyping arrays from a large clinical cohort.

BACKGROUND: Genotypes obtained with commercial SNP arrays have been extensively used in many large case-control or population-based cohorts for SNP-based genome-wide association studies for a multitude of traits. Yet, these genotypes capture only a small fraction of the variance of the studied traits. Genomic structural variants (GSV) such as Copy Number Variation (CNV) may account for part of the missing heritability, but their comprehensive detection requires either next-generation arrays or sequencing. Sophisticated algorithms that infer CNVs by combining the intensities from SNP-probes for the two alleles can already be used to extract a partial view of such GSV from existing data sets. RESULTS: Here we present several advances to facilitate the latter approach. First, we introduce a novel CNV detection method based on a Gaussian Mixture Model. Second, we propose a new algorithm, PCA merge, for combining copy-number profiles from many individuals into consensus regions. We applied both our new methods as well as existing ones to data from 5612 individuals from the CoLaus study who were genotyped on Affymetrix 500K arrays. We developed a number of procedures in order to evaluate the performance of the different methods. This includes comparison with previously published CNVs as well as using a replication sample of 239 individuals, genotyped with Illumina 550K arrays. We also established a new evaluation procedure that employs the fact that related individuals are expected to share their CNVs more frequently than randomly selected individuals. The ability to detect both rare and common CNVs provides a valuable resource that will facilitate association studies exploring potential phenotypic associations with CNVs. CONCLUSION: Our new methodologies for CNV detection and their evaluation will help in extracting additional information from the large amount of SNP-genotyping data on various cohorts and use this to explore structural variants and their impact on complex traits.

First consumption ever of multiple substances: applying an expert-based taxonomy to a Swiss national sample of adolescents.

The use of multiple legal and illegal substances by adolescents is a growing concern in all countries, but since no consensus about a taxonomy did emerge yet, it is difficult to understand the different patterns of consumption and to implement tailored prevention and treatment programs directed towards specific subgroups of the adolescent population. Using data from a Swiss survey on adolescent health, we analyzed the age at which ten legal and illegal substances were consumed for the first time ever by applying a method combining the strength of both automatic clustering and use of substance experts. Results were then compared to 30 socio-economic factors to establish the usefulness of and to validate our taxonomy. We also analyzed the succession of substance first use for each group. The final taxonomy consists of eight groups ranging from non-consumers to heavy drug addicts. All but four socio-economic factors were significantly associated with the taxonomy, the strongest associations being observed with health, behavior, and sexuality factors. Numerous factors influence adolescents in their decision to first try substances or to use them on a regular basis, and no factor alone can be considered as an absolute marker of problematic behavior regarding substance use. Different processes of experimentation with substances are associated with different behaviors, therefore focusing on only one substance or only one factor is not efficient. Prevention and treatment programs can then be tailored to address specific issues related to different youth subgroups.

Cell morphology and intracellular ionic homeostasis explored with a multimodal approach combining epifluorescence and digital holographic microscopy.

The authors have developed a live-cell multimodality microscope combining epifluorescence with digital holographic microscopy; it has been implemented with a decoupling procedure allowing to separately measure from the quantitative phase important cell parameters including absolute volume, shape and integral intracellular refractive index. In combination with the numerous different specific fluorescent cellular probes, this multimodality microscopy can address important issues in cell biology. This is demonstrated by the study of intracellular calcium homeostasis associated with the change in cell volume, which play a critical role in the excitotoxicity-induced neuronal death.

The systematic profiling of false identity documents: method validation and performance evaluation using seizures known to originate from common and different sources

False identity documents constitute a potential powerful source of forensic intelligence because they are essential elements of transnational crime and provide cover for organized crime. In previous work, a systematic profiling method using false documents' visual features has been built within a forensic intelligence model. In the current study, the comparison process and metrics lying at the heart of this profiling method are described and evaluated. This evaluation takes advantage of 347 false identity documents of four different types seized in two countries whose sources were known to be common or different (following police investigations and dismantling of counterfeit factories). Intra-source and inter-sources variations were evaluated through the computation of more than 7500 similarity scores. The profiling method could thus be validated and its performance assessed using two complementary approaches to measuring type I and type II error rates: a binary classification and the computation of likelihood ratios. Very low error rates were measured across the four document types, demonstrating the validity and robustness of the method to link documents to a common source or to differentiate them. These results pave the way for an operational implementation of a systematic profiling process integrated in a developed forensic intelligence model.

An optimized method for establishing high purity murine CD8(+) T cell cultures.

Establishing CD8(+) T cell cultures has been empirical and the published methods have been largely individual laboratory based. In this study, we optimized culturing conditions and show that IL-2 concentration is the most critical factor for the success of establishing CD8(+) T cell cultures. High IL-2 concentration encouraged T cells to non-specifically proliferate, express a B cell marker, B220, and undergo apoptosis. These cells also lose typical irregular T cell morphology and are incapable of sustaining long-term cultures. Using tetramer and intracellular cytokine assessments, we further demonstrated that many antigen-specific T cells have been rendered nonfunctional when expanded under high IL-2 concentration. When IL-2 is used in the correct range, B220-mediated cell depletion greatly enhanced the success rate of such T cell cultures.

Phosphorus saturation of Finnish soils : evaluating an easy oxalate extraction method

Selostus: Alumiini- ja rautaoksidien fosforikyllästysasteen arvioiminen suomalaisista peltomaista

Fecal pellets collection as a method for assessing egesta of the marine cave-dwelling mysid Hemimysis speluncula

Egesta of a cave-dwelling mysid (Hemimysis speluncola Ledoyer, 1963) was studied in a submarine cave of Medes Islands, NW Mediterranean by in situ fecal pellet collecting. Fecal pellet production and gut fullness of mysids during incubation experiments are used to estimate mysid egestion rates. Intrinsic factors related with the natural history of this species such as population structure, density of mysids, daily rhythms and pellet decomposition rates are tested for their influence on the egestion rate. The effects of methodological artifacts, such as the stress induced by both incubation and preservation procedures, are also studied. An average mysid egests about 2.5 pellets per day into the cave. The time of day is the main factor affecting egestion. The highest deposition rate is between 2 to 4 hours after sunrise when about 38 % of the total daily pellet production becomes egested. Fecal pellet morphology changes with mysid demographic classes: immature mysids produce slender and thick pellets, whereas mature mysids produce only thick pellets. Immature classes show higher percentages of full guts than mature ones. Mysid density in the incubators does not affect the results on gut fullness, but it causes a decrease in the number of pellets collected after incubation. Coprorhexia seems to be the only plausible process to explain this paradox. The incubation procedure does not increase deposition rate significantly. Time of incubation is critical because the half-life of fecal pellets is about 2.5 hours. Fixation with liquid nitrogen decreases gut fullness and also deposition rates. Higher values are obtained with 70 % ethanol and 5 % formalin solutions which show very similar results for both gut fullness and pellet deposition rates. Nevertheless, ethanol is not suitable as fixative because it enhances the opacity of the body. Several suggestions are given in order to optimize the reliability of further in situ experiments for evaluation of egesta of Hemimysis speluncola in submarine caves.