Center Nacional dAnalisi Genomica

Barcelona, Spain

Center Nacional dAnalisi Genomica

Barcelona, Spain

Time filter

Source Type

PubMed | Center Nacional dAnalisi Genomica, University of Barcelona, CIBER ISCIII, San Sebastián University and 3 more.
Type: Journal Article | Journal: European journal of human genetics : EJHG | Year: 2016

Colorectal cancer (CRC) is one of the most common neoplasms in the world. Fanconi anemia (FA) is a very rare genetic disease causing bone marrow failure, congenital growth abnormalities and cancer predisposition. The comprehensive FA DNA damage repair pathway requires the collaboration of 53 proteins and it is necessary to restore genome integrity by efficiently repairing damaged DNA. A link between FA genes in breast and ovarian cancer germline predisposition has been previously suggested. We selected 74 CRC patients from 40 unrelated Spanish families with strong CRC aggregation compatible with an autosomal dominant pattern of inheritance and without mutations in known hereditary CRC genes and performed germline DNA whole-exome sequencing with the aim of finding new candidate germline predisposition variants. After sequencing and data analysis, variant prioritization selected only those very rare alterations, producing a putative loss of function and located in genes with a role compatible with cancer. We detected an enrichment for variants in FA DNA damage repair pathway genes in our familial CRC cohort as 6 families carried heterozygous, rare, potentially pathogenic variants located in BRCA2/FANCD1, BRIP1/FANCJ, FANCC, FANCE and REV3L/POLZ. In conclusion, the FA DNA damage repair pathway may play an important role in the inherited predisposition to CRC.


Fenouil R.,Aix - Marseille University | Fenouil R.,French Institute of Health and Medical Research | Fenouil R.,French National Center for Scientific Research | Cauchy P.,Aix - Marseille University | And 26 more authors.
Genome Research | Year: 2012

One clear hallmark of mammalian promoters is the presence of CpG islands (CGIs) at more than two-thirds of genes, whereas TATA boxes are only present at a minority of promoters. Using genome-wide approaches, we show that GC content and CGIs are major promoter elements in mammalian cells, able to govern open chromatin conformation and support paused transcription. First, we define three classes of promoters with distinct transcriptional directionality and pausing properties that correlate with their GC content. We further analyze the direct influence of GC content on nucleosome positioning and depletion and show that CpG content and CGI width correlate with nucleosome depletion both in vivo and in vitro. We also show that transcription is not essential for nucleosome exclusion but influences both a weak +1 and a well-positioned nucleosome at CGI borders. Altogether our data support the idea that CGIs have become an essential feature of promoter structure defining novel regulatory properties in mammals. © 2012, Published by Cold Spring Harbor Laboratory Press.


PubMed | CSIC - Biological Research Center, Center Nacional dAnalisi Genomica, University of Santiago de Compostela, University Pompeu Fabra and CSIC - Institute of Marine Research
Type: Journal Article | Journal: DNA research : an international journal for rapid publication of reports on genes and genomes | Year: 2016

The turbot is a flatfish (Pleuronectiformes) with increasing commercial value, which has prompted active genomic research aimed at more efficient selection. Here we present the sequence and annotation of the turbot genome, which represents a milestone for both boosting breeding programmes and ascertaining the origin and diversification of flatfish. We compare the turbot genome with model fish genomes to investigate teleost chromosome evolution. We observe a conserved macrosyntenic pattern within Percomorpha and identify large syntenic blocks within the turbot genome related to the teleost genome duplication. We identify gene family expansions and positive selection of genes associated with vision and metabolism of membrane lipids, which suggests adaptation to demersal lifestyle and to cold temperatures, respectively. Our data indicate a quick evolution and diversification of flatfish to adapt to benthic life and provide clues for understanding their controversial origin. Moreover, we investigate the genomic architecture of growth, sex determination and disease resistance, key traits for understanding local adaptation and boosting turbot production, by mapping candidate genes and previously reported quantitative trait loci. The genomic architecture of these productive traits has allowed the identification of candidate genes and enriched pathways that may represent useful information for future marker-assisted selection in turbot.


De Almeida S.F.,University of Lisbon | Grosso A.R.,University of Lisbon | Koch F.,French Institute of Health and Medical Research | Fenouil R.,French Institute of Health and Medical Research | And 9 more authors.
Nature Structural and Molecular Biology | Year: 2011

Several lines of recent evidence support a role for chromatin in splicing regulation. Here, we show that splicing can also contribute to histone modification, which implies bidirectional communication between epigenetic mechanisms and RNA processing. Genome-wide analysis of histone methylation in human cell lines and mouse primary T cells reveals that intron-containing genes are preferentially marked with histone H3 Lys36 trimethylation (H3K36me3) relative to intronless genes. In intron-containing genes, H3K36me3 marking is proportional to transcriptional activity, whereas in intronless genes, H3K36me3 is always detected at much lower levels. Furthermore, splicing inhibition impairs recruitment of H3K36 methyltransferase HYPB (also known as Setd2) and reduces H3K36me3, whereas splicing activation has the opposite effect. Moreover, the increase of H3K36me3 correlates with the length of the first intron, consistent with the view that splicing enhances H3 methylation. We propose that splicing is mechanistically coupled to recruitment of HYPB/Setd2 to elongating RNA polymerase II. © 2011 Nature America, Inc. All rights reserved. (c) 2011 Nature America, Inc. All rights reserved.


Gonzalez-Porta M.,Center for Genomic Regulation | Gonzalez-Porta M.,European Bioinformatics Institute | Calvo M.,University of Barcelona | Sammeth M.,Center for Genomic Regulation | And 3 more authors.
Genome Research | Year: 2012

DNA arrays have been widely used to perform transcriptome-wide analysis of gene expression, and many methods have been developed to measure gene expression variability and to compare gene expression between conditions. Because RNA-seq is also becoming increasingly popular for transcriptome characterization, the possibility exists for further quantification of individual alternative transcript isoforms, and therefore for estimating the relative ratios of alternative splice forms within a given gene. Changes in splicing ratios, even without changes in overall gene expression, may have important phenotypic effects. Here we have developed statistical methodology to measure variability in splicing ratios within conditions, to compare it between conditions, and to identify genes with condition-specific splicing ratios. Furthermore, we have developed methodology to deconvolute the relative contribution of variability in gene expression versus variability in splicing ratios to the overall variability of transcript abundances. As a proof of concept, we have applied this methodology to estimates of transcript abundances obtained from RNA-seq experiments in lymphoblastoid cells from Caucasian and Yoruban individuals. We have found that protein-coding genes exhibit low splicing variability within populations, with many genes exhibiting constant ratios across individuals. When comparing these two populations, we have found that up to 10% of the studied protein-coding genes exhibit population-specific splicing ratios. We estimate that ∼60% of the total variability observed in the abundance of transcript isoforms can be explained by variability in transcription. A large fraction of the remaining variability can likely result from variability in splicing. Finally, we also detected that variability in splicing is uncommon without variability in transcription. © 2012 by Cold Spring Harbor Laboratory Press.


PubMed | Whitehead Institute For Biomedical Research, Center Nacional dAnalisi Genomica, Barcelona Supercomputing Center, 08908 Lhospitalet Of Llobregat and University of Gottingen
Type: | Journal: Genome biology | Year: 2016

One of the hallmarks of cancer is the disruption of gene expression patterns. Many molecular lesions contribute to this phenotype, and the importance of aberrant DNA methylation profiles is increasingly recognized. Much of the research effort in this area has examined proximal promoter regions and epigenetic alterations at other loci are not well characterized.Using whole genome bisulfite sequencing to examine uncharted regions of the epigenome, we identify a type of far-reaching DNA methylation alteration in cancer cells of the distal regulatory sequences described as super-enhancers. Human tumors undergo a shift in super-enhancer DNA methylation profiles that is associated with the transcriptional silencing or the overactivation of the corresponding target genes. Intriguingly, we observe locally active fractions of super-enhancers detectable through hypomethylated regions that suggest spatial variability within the large enhancer clusters. Functionally, the DNA methylomes obtained suggest that transcription factors contribute to this local activity of super-enhancers and that trans-acting factors modulate DNA methylation profiles with impact on transforming processes during carcinogenesis.We develop an extensive catalogue of human DNA methylomes at base resolution to better understand the regulatory functions of DNA methylation beyond those of proximal promoter gene regions. CpG methylation status in normal cells points to locally active regulatory sites at super-enhancers, which are targeted by specific aberrant DNA methylation events in cancer, with putative effects on the expression of downstream genes.


Eizirik D.L.,Free University of Colombia | Sammeth M.,Center Nacional dAnalisi Genomica | Bouckenooghe T.,Free University of Colombia | Bottu G.,Free University of Colombia | And 16 more authors.
PLoS Genetics | Year: 2012

Type 1 diabetes (T1D) is an autoimmune disease in which pancreatic beta cells are killed by infiltrating immune cells and by cytokines released by these cells. Signaling events occurring in the pancreatic beta cells are decisive for their survival or death in diabetes. We have used RNA sequencing (RNA-seq) to identify transcripts, including splice variants, expressed in human islets of Langerhans under control conditions or following exposure to the pro-inflammatory cytokines interleukin-1β (IL-1β) and interferon-γ (IFN-γ). Based on this unique dataset, we examined whether putative candidate genes for T1D, previously identified by GWAS, are expressed in human islets. A total of 29,776 transcripts were identified as expressed in human islets. Expression of around 20% of these transcripts was modified by pro-inflammatory cytokines, including apoptosis- and inflammation-related genes. Chemokines were among the transcripts most modified by cytokines, a finding confirmed at the protein level by ELISA. Interestingly, 35% of the genes expressed in human islets undergo alternative splicing as annotated in RefSeq, and cytokines caused substantial changes in spliced transcripts. Nova1, previously considered a brain-specific regulator of mRNA splicing, is expressed in islets and its knockdown modified splicing. 25/41 of the candidate genes for T1D are expressed in islets, and cytokines modified expression of several of these transcripts. The present study doubles the number of known genes expressed in human islets and shows that cytokines modify alternative splicing in human islet cells. Importantly, it indicates that more than half of the known T1D candidate genes are expressed in human islets. This, and the production of a large number of chemokines and cytokines by cytokine-exposed islets, reinforces the concept of a dialog between pancreatic islets and the immune system in T1D. This dialog is modulated by candidate genes for the disease at both the immune system and beta cell level. © 2012 Eizirik et al.


Cnop M.,Free University of Colombia | Cnop M.,Erasmus University College Brussels | Abdulkarim B.,Free University of Colombia | Bottu G.,Free University of Colombia | And 17 more authors.
Diabetes | Year: 2014

Pancreatic β-cell dysfunction and death are central in the pathogenesis of type 2 diabetes (T2D). Saturated fatty acids cause β-cell failure and contribute to diabetes development in genetically predisposed individuals. Here we used RNA sequencing to map transcripts expressed in five palmitate-treated human islet preparations, observing 1,325 modified genes. Palmitate induced fatty acid metabolism and endoplasmic reticulum (ER) stress. Functional studies identified novel mediators of adaptive ER stress signaling. Palmitate modified genes regulating ubiquitin and proteasome function, autophagy, and apoptosis. Inhibition of autophagic flux and lysosome function contributed to lipotoxicity. Palmitate inhibited transcription factors controlling β-cell phenotype, including PAX4 and GATA6. Fifty-nine T2D candidate genes were expressed in human islets, and 11 were modified by palmitate. Palmitate modified expression of 17 splicing factors and shifted alternative splicing of 3,525 transcripts. Ingenuity Pathway Analysis of modified transcripts and genes confirmed that top changed functions related to cell death. Database for Annotation, Visualization and Integrated Discovery (DAVID) analysis of transcription factor binding sites in palmitate-modified transcripts revealed a role for PAX4, GATA, and the ER stress response regulators XBP1 and ATF6. This human islet transcriptome study identified novel mechanisms of palmitate-induced β-cell dysfunction and death. The data point to cross talk between metabolic stress and candidate genes at the β-cell level. © 2014 by the American Diabetes Association.


Mundry M.,Westfaelische Wilhelms University | Bornberg-Bauer E.,Westfaelische Wilhelms University | Sammeth M.,Center Nacional dAnalisi Genomica | Feulner P.G.D.,Westfaelische Wilhelms University
PLoS ONE | Year: 2012

Background: The quantity of transcriptome data is rapidly increasing for non-model organisms. As sequencing technology advances, focus shifts towards solving bioinformatic challenges, of which sequence read assembly is the first task. Recent studies have compared the performance of different software to establish a best practice for transcriptome assembly. Here, we adapted a simulation approach to evaluate specific features of assembly programs on 454 data. The novelty of our study is that the simulation allows us to calculate a model assembly as reference point for comparison. Findings: The simulation approach allows us to compare basic metrics of assemblies computed by different software applications (CAP3, MIRA, Newbler, and Oases) to a known optimal solution. We found MIRA and CAP3 are conservative in merging reads. This resulted in comparably high number of short contigs. In contrast, Newbler more readily merged reads into longer contigs, while Oases produced the overall shortest assembly. Due to the simulation approach, reads could be traced back to their correct placement within the transcriptome. Together with mapping reads onto the assembled contigs, we were able to evaluate ambiguity in the assemblies. This analysis further supported the conservative nature of MIRA and CAP3, which resulted in low proportions of chimeric contigs, but high redundancy. Newbler produced less redundancy, but the proportion of chimeric contigs was higher. Conclusion: Our evaluation of four assemblers suggested that MIRA and Newbler slightly outperformed the other programs, while showing contrasting characteristics. Oases did not perform very well on the 454 reads. Our evaluation indicated that the software was either conservative (MIRA) or liberal (Newbler) about merging reads into contigs. This suggested that in choosing an assembly program researchers should carefully consider their follow up analysis and consequences of the chosen approach to gain an assembly. © 2012 Mundry et al.


Hintermair C.,Helmholtz Center Munich | Heidemann M.,Helmholtz Center Munich | Koch F.,Aix - Marseille University | Descostes N.,Aix - Marseille University | And 9 more authors.
EMBO Journal | Year: 2012

Eukaryotic RNA polymerase II (Pol II) has evolved an array of heptad repeats with the consensus sequence Tyr1-Ser2-Pro3-Thr4-Ser5-Pro6-Ser7 at the carboxy-terminal domain (CTD) of the large subunit (Rpb1). Differential phosphorylation of Ser2, Ser5, and Ser7 in the 5′ and 3′ regions of genes coordinates the binding of transcription and RNA processing factors to the initiating and elongating polymerase complexes. Here, we report phosphorylation of Thr4 by Polo-like kinase 3 in mammalian cells. ChIPseq analyses indicate an increase of Thr4-P levels in the 3′ region of genes occurring subsequently to an increase of Ser2-P levels. A Thr4/Ala mutant of Pol II displays a lethal phenotype. This mutant reveals a global defect in RNA elongation, while initiation is largely unaffected. Since Thr4 replacement mutants are viable in yeast we conclude that this amino acid has evolved an essential function(s) in the CTD of Pol II for gene transcription in mammalian cells. © 2012 European Molecular Biology Organization.

Loading Center Nacional dAnalisi Genomica collaborators
Loading Center Nacional dAnalisi Genomica collaborators