Benoit Ballester

Five-Vertebrate ChIP-seq Reveals the Evolutionary Dynamics of Transcription Factor Binding

Subtle Variation Despite vast phenotypic differences, vertebrates have many readily recognizable specific cell types, like liver hepatocytes. The gene expression that defines specific cells depends on evolutionarily conserved orthologous transcription factors. Schmidt et al. (p. 1036, published online 8 April) studied the conservation and divergence in the genome-wide binding of two such transcription factors, CEBPA and HNF4A, in livers from human, dog, mouse, short-tailed opossum, and chicken. Although the sequence bound by orthologous transcription factors was similar, the vast majority of binding events were unique to each species. Binding of two liver-specific transcription factors in several vertebrate species reveals complex regulatory evolution. Transcription factors (TFs) direct gene expression by binding to DNA regulatory regions. To explore the evolution of gene regulation, we used chromatin immunoprecipitation with high-throughput sequencing (ChIP-seq) to determine experimentally the genome-wide occupancy of two TFs, CCAAT/enhancer-binding protein alpha and hepatocyte nuclear factor 4 alpha, in the livers of five vertebrates. Although each TF displays highly conserved DNA binding preferences, most binding is species-specific, and aligned binding events present in all five species are rare. Regions near genes with expression levels that are dependent on a TF are often bound by the TF in multiple species yet show no enhanced DNA sequence constraint. Binding divergence between species can be largely explained by sequence changes to the bound motifs. Among the binding events lost in one lineage, only half are recovered by another binding event within 10 kilobases. Our results reveal large interspecies differences in transcriptional regulation and provide insight into regulatory evolution.

BioMart Central Portal—unified access to biological data

BioMart Central Portal (www.biomart.org) offers a one-stop shop solution to access a wide array of biological databases. These include major biomolecular sequence, pathway and annotation databases such as Ensembl, Uniprot, Reactome, HGNC, Wormbase and PRIDE; for a complete list, visit, http://www.biomart.org/biomart/martview. Moreover, the web server features seamless data federation making cross querying of these data sources in a user friendly and unified way. The web server not only provides access through a web interface (MartView), it also supports programmatic access through a Perl API as well as RESTful and SOAP oriented web services. The website is free and open to all users and there is no login requirement.

Ensembl’s 10th year

Ensembl (http://www.ensembl.org) integrates genomic information for a comprehensive set of chordate genomes with a particular focus on resources for human, mouse, rat, zebrafish and other high-value sequenced genomes. We provide complete gene annotations for all supported species in addition to specific resources that target genome variation, function and evolution. Ensembl data is accessible in a variety of formats including via our genome browser, API and BioMart. This year marks the tenth anniversary of Ensembl and in that time the project has grown with advances in genome technology. As of release 56 (September 2009), Ensembl supports 51 species including marmoset, pig, zebra finch, lizard, gorilla and wallaby, which were added in the past year. Major additions and improvements to Ensembl since our previous report include the incorporation of the human GRCh37 assembly, enhanced visualisation and data-mining options for the Ensembl regulatory features and continued development of our software infrastructure.

Waves of Retrotransposon Expansion Remodel Genome Organization and CTCF Binding in Multiple Mammalian Lineages

Summary CTCF-binding locations represent regulatory sequences that are highly constrained over the course of evolution. To gain insight into how these DNA elements are conserved and spread through the genome, we defined the full spectrum of CTCF-binding sites, including a 33/34-mer motif, and identified over five thousand highly conserved, robust, and tissue-independent CTCF-binding locations by comparing ChIP-seq data from six mammals. Our data indicate that activation of retroelements has produced species-specific expansions of CTCF binding in rodents, dogs, and opossum, which often functionally serve as chromatin and transcriptional insulators. We discovered fossilized repeat elements flanking deeply conserved CTCF-binding regions, indicating that similar retrotransposon expansions occurred hundreds of millions of years ago. Repeat-driven dispersal of CTCF binding is a fundamental, ancient, and still highly active mechanism of genome evolution in mammalian lineages. PaperClip

Gene expression profiling identifies molecular subgroups among nodal peripheral T-cell lymphomas.

The classification of peripheral T-cell lymphomas (PTCL) is still a matter of debate. To establish a molecular classification of PTCL, we analysed 59 primary nodal T-cell lymphomas using cDNA microarrays, including 56 PTCL and three T-lymphoblastic lymphoma (T-LBL). The expression profiles could discriminate angioimmunoblastic lymphoma, anaplastic large-cell lymphoma and T-LBL. In contrast, cases belonging to the broad category of ‘PTCL, unspecified’ (PTCL-U) did not share a single molecular profile. Using a multiclass predictor, we could separate PTCL-U into three molecular subgroups called U1, U2 and U3. The U1 gene expression signature included genes known to be associated with poor outcome in other tumors, such as CCND2. The U2 subgroup was associated with overexpression of genes involved in T-cell activation and apoptosis, including NFKB1 and BCL-2. The U3 subgroup was mainly defined by overexpression of genes involved in the IFN/JAK/STAT pathway. It comprised a majority of histiocyte-rich PTCL samples. Gene Ontology annotations revealed different functional profile for each subgroup. These results suggest the existence of distinct subtypes of PTCL-U with specific molecular profiles, and thus provide a basis to improve their classification and to develop new therapeutic targets.

Integrative analysis of public ChIP-seq experiments reveals a complex multi-cell regulatory landscape

The large collections of ChIP-seq data rapidly accumulating in public data warehouses provide genome-wide binding site maps for hundreds of transcription factors (TFs). However, the extent of the regulatory occupancy space in the human genome has not yet been fully apprehended by integrating public ChIP-seq data sets and combining it with ENCODE TFs map. To enable genome-wide identification of regulatory elements we have collected, analysed and retained 395 available ChIP-seq data sets merged with ENCODE peaks covering a total of 237 TFs. This enhanced repertoire complements and refines current genome-wide occupancy maps by increasing the human genome regulatory search space by 14% compared to ENCODE alone, and also increases the complexity of the regulatory dictionary. As a direct application we used this unified binding repertoire to annotate variant enhancer loci (VELs) from H3K4me1 mark in two cancer cell lines (MCF-7, CRC) and observed enrichments of specific TFs involved in biological key functions to cancer development and proliferation. Those enrichments of TFs within VELs provide a direct annotation of non-coding regions detected in cancer genomes. Finally, full access to this catalogue is available online together with the TFs enrichment analysis tool (http://tagc.univ-mrs.fr/remap/).

JASPAR 2018: update of the open-access database of transcription factor binding profiles and its web framework

Abstract JASPAR (http://jaspar.genereg.net) is an open-access database of curated, non-redundant transcription factor (TF)-binding profiles stored as position frequency matrices (PFMs) and TF flexible models (TFFMs) for TFs across multiple species in six taxonomic groups. In the 2018 release of JASPAR, the CORE collection has been expanded with 322 new PFMs (60 for vertebrates and 262 for plants) and 33 PFMs were updated (24 for vertebrates, 8 for plants and 1 for insects). These new profiles represent a 30% expansion compared to the 2016 release. In addition, we have introduced 316 TFFMs (95 for vertebrates, 218 for plants and 3 for insects). This release incorporates clusters of similar PFMs in each taxon and each TF class per taxon. The JASPAR 2018 CORE vertebrate collection of PFMs was used to predict TF-binding sites in the human genome. The predictions are made available to the scientific community through a UCSC Genome Browser track data hub. Finally, this update comes with a new web framework with an interactive and responsive user-interface, along with new features. All the underlying data can be retrieved programmatically using a RESTful API and through the JASPAR 2018 R/Bioconductor package.

Multi-species, multi-transcription factor binding highlights conserved control of tissue-specific biological pathways

As exome sequencing gives way to genome sequencing, the need to interpret the function of regulatory DNA becomes increasingly important. To test whether evolutionary conservation of cis-regulatory modules (CRMs) gives insight into human gene regulation, we determined transcription factor (TF) binding locations of four liver-essential TFs in liver tissue from human, macaque, mouse, rat, and dog. Approximately, two thirds of the TF-bound regions fell into CRMs. Less than half of the human CRMs were found as a CRM in the orthologous region of a second species. Shared CRMs were associated with liver pathways and disease loci identified by genome-wide association studies. Recurrent rare human disease causing mutations at the promoters of several blood coagulation and lipid metabolism genes were also identified within CRMs shared in multiple species. This suggests that multi-species analyses of experimentally determined combinatorial TF binding will help identify genomic regions critical for tissue-specific gene control. DOI: http://dx.doi.org/10.7554/eLife.02626.001

The chromatin environment shapes DNA replication origin organization and defines origin classes

To unveil the still-elusive nature of metazoan replication origins, we identified them genome-wide and at unprecedented high-resolution in mouse ES cells. This allowed initiation sites (IS) and initiation zones (IZ) to be differentiated. We then characterized their genetic signatures and organization and integrated these data with 43 chromatin marks and factors. Our results reveal that replication origins can be grouped into three main classes with distinct organization, chromatin environment, and sequence motifs. Class 1 contains relatively isolated, low-efficiency origins that are poor in epigenetic marks and are enriched in an asymmetric AC repeat at the initiation site. Late origins are mainly found in this class. Class 2 origins are particularly rich in enhancer elements. Class 3 origins are the most efficient and are associated with open chromatin and polycomb protein-enriched regions. The presence of Origin G-rich Repeated elements (OGRE) potentially forming G-quadruplexes (G4) was confirmed at most origins. These coincide with nucleosome-depleted regions located upstream of the initiation sites, which are associated with a labile nucleosome containing H3K64ac. These data demonstrate that specific chromatin landscapes and combinations of specific signatures regulate origin localization. They explain the frequently observed links between DNA replication and transcription. They also emphasize the plasticity of metazoan replication origins and suggest that in multicellular eukaryotes, the combination of distinct genetic features and chromatin configurations act in synergy to define and adapt the origin profile.

Gene profiling reveals specific oncogenic mechanisms and signaling pathways in oncocytic and papillary thyroid carcinoma

The oncogenic pathways in mitochondrial-rich thyroid carcinomas are not clearly understood. To investigate the possible implication of mitochondrial abundance in the genesis of thyroid tumors, we have explored the gene expression profile of six oncocytic carcinomas and six mitochondrial-rich papillary carcinomas using cDNA-microarray technology. A supervised approach allowed us to identify 83 genes differentially expressed in the two types of carcinoma. These genes were classified according to their ontologic profiles. Three genes, NOS3, alpha-actinin-2 and alpha-catenin, suspected of playing a role in tumor genesis, were explored by quantitative RT–PCR analysis and immunohistochemistry. Of the 59 genes overexpressed in papillary carcinomas, 51% were involved in cell communication. Of the 24 genes overexpressed in oncocytic carcinomas, 84% were involved in mitochondrial and cellular metabolism. Our results suggest that mitochondrial respiratory chain complexes III and IV play a significant role in the regulation of reactive oxygen species production by oncocytic tumors.