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.

Genetic variants of the FADS gene cluster and ELOVL gene family, colostrums LC-PUFA levels, breastfeeding, and child cognition

Introduction: Breastfeeding effects on cognition are attributed to long-chain polyunsaturated fatty acids (LC-PUFAs), but controversy persists. Genetic variation in fatty acid desaturase (FADS) and elongase (ELOVL) enzymes has been overlooked when studying the effects of LC-PUFAs supply on cognition. We aimed to: 1) to determine whether maternal genetic variants in the FADS cluster and ELOVL genes contribute to differences in LC-PUFA levels in colostrum; 2) to analyze whether these maternal variants are related to child cognition; and 3) to assess whether children's variants modify breastfeeding effects on cognition. Methods: Data come from two population-based birth cohorts (n = 400 mother-child pairs from INMA-Sabadell; and n = 340 children from INMA-Menorca). LC-PUFAs were measured in 270 colostrum samples from INMA-Sabadell. Tag SNPs were genotyped both in mothers and children (13 in the FADS cluster, 6 in ELOVL2, and 7 in ELOVL5). Child cognition was assessed at 14 mo and 4 y using the Bayley Scales of Infant Development and the McCarthy Scales of Children"s Abilities, respectively. Results: Children of mothers carrying genetic variants associated with lower FADS1 activity (regulating AA and EPA synthesis), higher FADS2 activity (regulating DHA synthesis), and with higher EPA/AA and DHA/AA ratios in colostrum showed a significant advantage in cognition at 14 mo (3.5 to 5.3 points). Not being breastfed conferred an 8- to 9-point disadvantage in cognition among children GG homozygote for rs174468 (low FADS1 activity) but not among those with the A allele. Moreover, not being breastfed resulted in a disadvantage in cognition (5 to 8 points) among children CC homozygote for rs2397142 (low ELOVL5 activity), but not among those carrying the G allele. Conclusion: Genetically determined maternal supplies of LC-PUFAs during pregnancy and lactation appear to be crucial for child cognition. Breastfeeding effects on cognition are modified by child genetic variation in fatty acid desaturase and elongase enzymes.

Treatment patterns and health care resource utilization in a 1-year observational cohort study of outpatients with schizophrenia at risk of nonadherence treated with long-acting injectable antipsychotics

Purpose: To describe (1) the clinical profiles and the patterns of use of long-acting injectable (LAI) antipsychotics in patients with schizophrenia at risk of nonadherence with oral antipsychotics, and in those who started treatment with LAI antipsychotics, (2) health care resource utilization and associated costs. Patients and methods: A total of 597 outpatients with schizophrenia at risk of nonadherence, according to the psychiatrist's clinical judgment, were recruited at 59 centers in a noninterventional prospective observational study of 1-year follow-up when their treatment was modified. In a post hoc analysis, the profiles of patients starting LAI or continuing with oral antipsychotics were described, and descriptive analyses of treatments, health resource utilization, and direct costs were performed in those who started an LAI antipsychotic. Results: Therapy modifications involved the antipsychotic medications in 84.8% of patients, mostly because of insufficient efficacy of prior regimen. Ninety-two (15.4%) patients started an LAI antipsychotic at recruitment. Of these, only 13 (14.1%) were prescribed with first-generation antipsychotics. During 1 year, 16.3% of patients who started and 14.9% of patients who did not start an LAI antipsychotic at recruitment relapsed, contrasting with the 20.9% who had been hospitalized only within the prior 6 months. After 1 year, 74.3% of patients who started an LAI antipsychotic continued concomitant treatment with oral antipsychotics. The mean (median) total direct health care cost per patient per month during the study year among the patients starting any LAI antipsychotic at baseline was 1,407 ( 897.7). Medication costs (including oral and LAI antipsychotics and concomitant medication) represented almost 44%, whereas nonmedication costs accounted for more than 55% of the mean total direct health care costs. Conclusion: LAI antipsychotics were infrequently prescribed in spite of a psychiatrist-perceived risk of nonadherence to oral antipsychotics. Mean medication costs were lower than nonmedication costs.

Mixed mating in androdioecious Mercurialis annua inferred using progeny arrays and diploid-acting microsatellite loci in a hexaploid background.

Background and Aims The frequency at which males can be maintained with hermaphrodites in androdioecious populations is predicted to depend on the selfing rate, because self-fertilization by hermaphrodites reduces prospective siring opportunities for males. In particular, high selfing rates by hermaphrodites are expected to exclude males from a population. Here, the first estimates are provided of the mating system from two wild hexaploid populations of the androdioecious European wind-pollinated plant M. annua with contrasting male frequencies.Methods Four diploid microsatellite loci were used to genotype 19-20 progeny arrays from two populations of M. annua, one with males and one without. Mating-system parameters were estimated using the program MLTR.Key Results Both populations had similar, intermediate outcrossing rates (t(m) = 0.64 and 0.52 for the population with and without males, respectively). The population without males showed a lower level of correlated paternity and biparental inbreeding and higher allelic richness and gene diversity than the population with males.Conclusions The results demonstrate the utility of new diploid microsatellite loci for mating system analysis in a hexaploid plant. It would appear that androdioecious M. annua has a mixed-mating system in the wild, an uncommon finding for wind-pollinated species. This study sets a foundation for future research to assess the relative importance of the sexual system, plant-density variation and stochastic processes for the regulation of male frequencies in M. annua over space and time.

aIr-1, a newly found locus on mouse chromosome 16 encoding a trans-acting activator factor for MHC class II gene expression.

RJ 2.2.5 is a human B cell line that has lost the capacity to express MHC class II genes. The human class II-positive phenotype is restored in somatic cell hybrids between RJ 2.2.5 and mouse spleen cells. By karyotype and molecular studies of an informative family of hybrids we have now shown that the reexpression of human class II gene products, as well as the maintenance of the mouse class II-positive phenotype, correlates with the presence of mouse chromosome 16. Thus, the existence on this mouse chromosome of a newly found locus, designated by us aIr-1, that determines a trans-acting activator function for class II gene expression, is established. Possible implications of this finding are discussed.

Analysis of the alternative pathways for the beta-oxidation of unsaturated fatty acids using transgenic plants synthesizing polyhydroxyalkanoates in peroxisomes.

Degradation of fatty acids having cis-double bonds on even-numbered carbons requires the presence of auxiliary enzymes in addition to the enzymes of the core beta-oxidation cycle. Two alternative pathways have been described to degrade these fatty acids. One pathway involves the participation of the enzymes 2, 4-dienoyl-coenzyme A (CoA) reductase and Delta(3)-Delta(2)-enoyl-CoA isomerase, whereas the second involves the epimerization of R-3-hydroxyacyl-CoA via a 3-hydroxyacyl-CoA epimerase or the action of two stereo-specific enoyl-CoA hydratases. Although degradation of these fatty acids in bacteria and mammalian peroxisomes was shown to involve mainly the reductase-isomerase pathway, previous analysis of the relative activity of the enoyl-CoA hydratase II (also called R-3-hydroxyacyl-CoA hydro-lyase) and 2,4-dienoyl-CoA reductase in plants indicated that degradation occurred mainly through the epimerase pathway. We have examined the implication of both pathways in transgenic Arabidopsis expressing the polyhydroxyalkanoate synthase from Pseudomonas aeruginosa in peroxisomes and producing polyhydroxyalkanoate from the 3-hydroxyacyl-CoA intermediates of the beta-oxidation cycle. Analysis of the polyhydroxyalkanoate synthesized in plants grown in media containing cis-10-heptadecenoic or cis-10-pentadecenoic acids revealed a significant contribution of both the reductase-isomerase and epimerase pathways to the degradation of these fatty acids.

Identification of Ligands for the Tau Exon 10 Splicing Regulatory Element RNA Using Dynamic Combinatorial Chemistry

We describe the use of dynamic combinatorial chemistry (DCC) to identify ligands for the stem-loop structure located at the exon 10-5'-intron junction of Tau pre-mRNA, which is involved in the onset of several tauopathies including frontotemporal dementia with Parkinsonism linked to chromosome 17 (FTDP-17). A series of ligands that combine the small aminoglycoside neamine and heteroaromatic moieties (azaquinolone and two acridines) have been identified by using DCC. These compounds effectively bind the stem-loop RNA target (the concentration required for 50% RNA response (EC(50)): 2-58 μM), as determined by fluorescence titration experiments. Importantly, most of them are able to stabilize both the wild-type and the +3 and +14 mutated sequences associated with the development of FTDP-17 without producing a significant change in the overall structure of the RNA (as analyzed by circular dichroism (CD) spectroscopy), which is a key factor for recognition by the splicing regulatory machinery. A good correlation has been found between the affinity of the ligands for the target and their ability to stabilize the RNA secondary structure.

9-cis-retinoic acid enhances fatty acid-induced expression of the liver fatty acid-binding protein gene.

The role of retinoic acids (RA) on liver fatty acid-binding protein (L-FABP) expression was investigated in the well differentiated FAO rat hepatoma cell line. 9-cis-Retinoic acid (9-cis-RA) specifically enhanced L-FABP mRNA levels in a time- and dose-dependent manner. The higher induction was found 6 h after addition of 10(-6) M 9-cis-RA in the medium. RA also enhanced further both L-FABP mRNA levels and cytosolic L-FABP protein content induced by oleic acid. The retinoid X receptor (RXR) and the peroxisome proliferator-activated receptor (PPAR), which are known to be activated, respectively, by 9-cis-RA and long chain fatty acid (LCFA), co-operated to bind specifically the peroxisome proliferator-responsive element (PPRE) found upstream of the L-FABP gene. Our result suggest that the PPAR-RXR complex is the molecular target by which 9-cis-RA and LCFA regulate the L-FABP gene.

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.

Interactions of Ly49 family receptors with MHC class I ligands in trans and cis.

The Ly49A NK cell receptor interacts with MHC class I (MHC-I) molecules on target cells and negatively regulates NK cell-mediated target cell lysis. We have recently shown that the MHC-I ligand-binding capacity of the Ly49A NK cell receptor is controlled by the NK cells' own MHC-I. To see whether this property was unique to Ly49A, we have investigated the binding of soluble MHC-I multimers to the Ly49 family receptors expressed in MHC-I-deficient and -sufficient C57BL/6 mice. In this study, we confirm the binding of classical MHC-I to the inhibitory Ly49A, C and I receptors, and demonstrate that detectable MHC-I binding to MHC-I-deficient NK cells is exclusively mediated by these three receptors. We did not detect significant multimer binding to stably transfected or NK cell-expressed Ly49D, E, F, G, and H receptors. Yet, we identified the more distantly related Ly49B and Ly49Q, which are not expressed by NK cells, as two novel MHC-I receptors in mice. Furthermore, we show using MHC-I-sufficient mice that the NK cells' own MHC-I significantly masks the Ly49A and Ly49C, but not the Ly49I receptor. Nevertheless, Ly49I was partly masked on transfected tumor cells, suggesting that the structure of Ly49I is compatible in principal with cis binding of MHC-I. Finally, masking of Ly49Q by cis MHC-I was minor, whereas masking of Ly49B was not detected. These data significantly extend the MHC-I specificity of Ly49 family receptors and show that the accessibility of most, but not all, MHC-I-binding Ly49 receptors is modulated by the expression of MHC-I in cis.

Selectivity in the binding of psychotropic drugs to the variants of alpha-1 acid glycoprotein.

The S- and F-forms of alpha-1 acid glycoprotein (AAG) variants have been isolated by isoelectric focusing with immobilines from commercially available AAG. In equilibrium dialysis experiments using a multicompartmental system, a higher affinity for various basic drugs has been found with S- in comparison with F-AAG: Amitriptyline, nortriptyline, imipramine, desipramine, trimipramine, methadone, thioridazine, clomipramine, desmethylclomipramine, and maprotiline. The selectivity (binding to S- vs. F-AAG) is the most pronounced for methadone and the lowest for thioridazine, while it is absent for the acidic drug mephenytoin.

Human leukocyte antigen class II variants and adult-onset asthma: does occupational allergen exposure play a role?

Recently, a locus centred on rs9273349 in the HLA-DQ region emerged from genome-wide association studies of adult-onset asthma. We aimed to further investigate the role of human leukocyte antigen (HLA) class II in adult-onset asthma and a possible interaction with occupational exposures. We imputed classical HLA-II alleles from 7579 single-nucleotide polymorphisms in 6025 subjects (1202 with adult-onset asthma) from European cohorts: ECRHS, SAPALDIA, EGEA and B58C, and from surveys of bakers and agricultural workers. Based on an asthma-specific job-exposure matrix, 2629 subjects had ever been exposed to high molecular weight (HMW) allergens. We explored associations between 23 common HLA-II alleles and adult-onset asthma, and tested for gene-environment interaction with occupational exposure to HMW allergens. Interaction was also tested for rs9273349. Marginal associations of classical HLA-II alleles and adult-onset asthma were not statistically significant. Interaction was detected between the DPB1*03:01 allele and exposure to HMW allergens (p = 0.009), in particular to latex (p = 0.01). In the unexposed group, the DPB1*03:01 allele was associated with adult-onset asthma (OR 0.67, 95%CI 0.53-0.86). HMW allergen exposures did not modify the association of rs9273349 with adult-onset asthma. Common classical HLA-II alleles were not marginally associated with adult-onset asthma. The association of latex exposure and adult-onset asthma may be modified by DPB1*03:01.

Variants : The Journal of the European Society for Textual Scholarship

This volume is the 10th issue of Variants . In keeping with the mission of the European Society for Textual Scholarship, the articles are richly interdisciplinary and transnational. They bring to bear a wide range of topics and disciplines on the field of textual scholarship: historical linguistics, digital scholarly editing, classical philology, Dutch, English, Finnish and Swedish Literature, publishing traditions in Japan, book history, cultural history and folklore. The questions that are explored - what texts are worth editing? what is the nature of the relationship between text, work, document and book? what is a critical digital edition? - all return to fundamental issues that have been at the heart of the editorial discipline for decades. With refreshing insight they assess the increasingly hybrid nature of the theoretical considerations and practical methodologies employed by textual scholars, while reasserting the relevance and need for producing scholarly editions, whether in print or digital, and continuing advanced research in bibliographical codes, textual transmissions, genetic dossiers, the fluidity of texts and other such subjects that connect textual scholarship with broader investigations into our nations' literary culture and written heritage.

Temperature-responsive sensing regulates biocontrol factor expression in Pseudomonas fluorescens CHA0.

In the plant-beneficial, root-colonizing strain Pseudomonas fluorescens CHA0, the Gac/Rsm signal transduction pathway positively regulates the synthesis of biocontrol factors (mostly antifungal secondary metabolites) and contributes to oxidative stress response via the stress sigma factor RpoS. The backbone of this pathway consists of the GacS/GacA two-component system, which activates the expression of three small regulatory RNAs (RsmX, RsmY, RsmZ) and thereby counters translational repression exerted by the RsmA and RsmE proteins on target mRNAs encoding biocontrol factors. We found that the expression of typical biocontrol factors, that is, antibiotic compounds and hydrogen cyanide (involving the phlA and hcnA genes), was significantly lower at 35 degrees C than at 30 degrees C. The expression of the rpoS gene was affected in parallel. This temperature control depended on RetS, a sensor kinase acting as an antagonist of the GacS/GacA system. An additional sensor kinase, LadS, which activated the GacS/GacA system, apparently did not contribute to thermosensitivity. Mutations in gacS or gacA were epistatic to (that is, they overruled) mutations in retS or ladS for expression of the small RNAs RsmXYZ. These data are consistent with a model according to which RetS-GacS and LadS-GacS interactions shape the output of the Gac/Rsm pathway and the environmental temperature influences the RetS-GacS interaction in P. fluorescens CHA0.