Search

Background: The treatment of melanoma, the deadliest form of skin cancer, has greatly benefited from immunotherapy. However, many patients do not show a durable response, which is only partially explained by known resistance mechanisms., Methods: We performed single-cell RNA sequencing of tumor immune infiltrates and matched peripheral blood mononuclear cells of 22 checkpoint inhibitor (CPI)-naive stage III-IV metastatic melanoma patients. After sample collection, the same patients received CPI treatment, and their response was assessed., Findings: CPI responders showed high levels of classical monocytes in peripheral blood, which preferentially transitioned toward CXCL9-expressing macrophages in tumors. Trajectories of tumor-infiltrating CD8 + T cells diverged at the level of effector memory/stem-like T cells, with non-responder cells progressing into a state characterized by cellular stress and apoptosis-related gene expression. Consistently, predicted non-responder-enriched myeloid-T/natural killer cell interactions were primarily immunosuppressive, while responder-enriched interactions were supportive of T cell priming and effector function., Conclusions: Our study illustrates that the tumor immune microenvironment prior to CPI treatment can be indicative of response. In perspective, modulating the myeloid and/or effector cell compartment by altering the described cell interactions and transitions could improve immunotherapy response., Funding: This research was funded by Roche Pharma Research and Early Development., Competing Interests: Declaration of interests R.S., S. Hoves, M.B., T.H., S.D., E.Y., F.L., S. Herter, P.U., A.I., J.S., J.S.-P., T.K.-T., I.G.d.M., and P.C.S. are employed by and hold F. Hoffmann-La Roche Ltd. company stock. P.U. and M.B. disclose ownership of F. Hoffmann-La Roche Ltd. patents. R.D. has intermittent, project-focused consulting and/or advisory relationships with Novartis, Merck Sharp & Dohme (MSD), Bristol-Myers Squibb (BMS), Roche, Amgen, Takeda, Pierre Fabre, Sun Pharma, Sanofi, Catalym, Second Genome, Regeneron, Alligator, T3 Pharma, MaxiVAX SA, Pfizer, and touchIME outside of the submitted work. M.P.L. receives project-specific research support outside of the scope of this work from Roche, Novartis, Molecular Partners, and Oncobit., (Copyright © 2024 The Authors. Published by Elsevier Inc. All rights reserved.)

Background: The immune status of a patient's tumor microenvironment (TME) may guide therapeutic interventions with cancer immunotherapy and help identify potential resistance mechanisms. Currently, patients' immune status is mostly classified based on CD8+tumor-infiltrating lymphocytes. An unmet need exists for comparable and reliable precision immunophenotyping tools that would facilitate clinical treatment-relevant decision-making and the understanding of how to overcome resistance mechanisms., Methods: We systematically analyzed the CD8 immunophenotype of 2023 patients from 14 phase I-III clinical trials using immunohistochemistry (IHC) and additionally profiled gene expression by RNA-sequencing (RNA-seq). CD8 immunophenotypes were classified by pathologists into CD8-desert, CD8-excluded or CD8-inflamed tumors using CD8 IHC staining in epithelial and stromal areas of the tumor. Using regularized logistic regression, we developed an RNA-seq-based classifier as a surrogate to the IHC-based spatial classification of CD8+tumor-infiltrating lymphocytes in the TME., Results: The CD8 immunophenotype and associated gene expression patterns varied across indications as well as across primary and metastatic lesions. Melanoma and kidney cancers were among the strongest inflamed indications, while CD8-desert phenotypes were most abundant in liver metastases across all tumor types. A good correspondence between the transcriptome and the IHC-based evaluation enabled us to develop a 92-gene classifier that accurately predicted the IHC-based CD8 immunophenotype in primary and metastatic samples (area under the curve inflamed=0.846; excluded=0.712; desert=0.855). The newly developed classifier was prognostic in The Cancer Genome Atlas (TCGA) data and predictive in lung cancer: patients with predicted CD8-inflamed tumors showed prolonged overall survival (OS) versus patients with CD8-desert tumors (HR 0.88; 95% CI 0.80 to 0.97) across TCGA, and longer OS on immune checkpoint inhibitor administration (phase III OAK study) in non-small-cell lung cancer (HR 0.75; 95% CI 0.58 to 0.97)., Conclusions: We provide a new precision immunophenotyping tool based on gene expression that reflects the spatial infiltration patterns of CD8+ lymphocytes in tumors. The classifier enables multiplex analyses and is easy to apply for retrospective, reverse translation approaches as well as for prospective patient enrichment to optimize the response to cancer immunotherapy., Competing Interests: Competing interests: AR, IID, PCS, MLS-S, NS, AV and KK are employees of F. Hoffmann-La Roche. AH, CSF, GD and MAC are employees of Roche Diagnostics. IK was an employee of Roche Diagnostics at the time of the study. AH, AR, IID, PCS, MLS-S, NS, CSF, AV, KK and MAC are shareholders of F. Hoffmann-La Roche., (© Author(s) (or their employer(s)) 2024. Re-use permitted under CC BY-NC. No commercial re-use. See rights and permissions. Published by BMJ.)

The consensus molecular subtypes (CMS) of colorectal cancer (CRC) is the most widely-used gene expression-based classification and has contributed to a better understanding of disease heterogeneity and prognosis. Nevertheless, CMS intratumoral heterogeneity restricts its clinical application, stressing the necessity of further characterizing the composition and architecture of CRC. Here, we used Spatial Transcriptomics (ST) in combination with single-cell RNA sequencing (scRNA-seq) to decipher the spatially resolved cellular and molecular composition of CRC. In addition to mapping the intratumoral heterogeneity of CMS and their microenvironment, we identified cell communication events in the tumor-stroma interface of CMS2 carcinomas. This includes tumor growth-inhibiting as well as -activating signals, such as the potential regulation of the ETV4 transcriptional activity by DCN or the PLAU-PLAUR ligand-receptor interaction. Our study illustrates the potential of ST to resolve CRC molecular heterogeneity and thereby help advance personalized therapy., (© 2024. The Author(s).)

Expansion and differentiation of antigen-experienced PD-1 + TCF-1 + stem-like CD8 + T cells into effector cells is critical for the success of immunotherapies based on PD-1 blockade 1-4 . Hashimoto et al. have shown that, in chronic infections, administration of the cytokine interleukin (IL)-2 triggers an alternative differentiation path of stem-like T cells towards a distinct population of 'better effector' CD8 + T cells similar to those generated in an acute infection 5 . IL-2 binding to the IL-2 receptor α-chain (CD25) was essential in triggering this alternative differentiation path and expanding better effectors with distinct transcriptional and epigenetic profiles. However, constitutive expression of CD25 on regulatory T cells and some endothelial cells also contributes to unwanted systemic effects from IL-2 therapy. Therefore, engineered IL-2 receptor β- and γ-chain (IL-2Rβγ)-biased agonists are currently being developed 6-10 . Here we show that IL-2Rβγ-biased agonists are unable to preferentially expand better effector T cells in cancer models and describe PD1-IL2v, a new immunocytokine that overcomes the need for CD25 binding by docking in cis to PD-1. Cis binding of PD1-IL2v to PD-1 and IL-2Rβγ on the same cell recovers the ability to differentiate stem-like CD8 + T cells into better effectors in the absence of CD25 binding in both chronic infection and cancer models and provides superior efficacy. By contrast, PD-1- or PD-L1-blocking antibodies alone, or their combination with clinically relevant doses of non-PD-1-targeted IL2v, cannot expand this unique subset of better effector T cells and instead lead to the accumulation of terminally differentiated, exhausted T cells. These findings provide the basis for the development of a new generation of PD-1 cis-targeted IL-2R agonists with enhanced therapeutic potential for the treatment of cancer and chronic infections., (© 2022. The Author(s).)

Background: Dendritic cells (DCs) are professional antigen presenting cells that initiate immune defense to pathogens and tumor cells. Human tumors contain only few DCs that mostly display a non-activated phenotype. Hence, activation of tumor-associated DCs may improve efficacy of cancer immunotherapies. Toll-like receptor (TLR) agonists and interferons are known to promote DC maturation. However, it is unclear if DCs in human tumors respond to activation signals and which stimuli induce the optimal activation of human tumor DCs., Methods: We first screened combinations of TLR agonists, a STING agonist and interferons (IFNs) for their ability to activate human conventional DCs (cDCs). Two combinations: TL8-506 (a TLR8 agonist)+IFN-γ and TL8-506+Poly(I:C) (a TLR3 agonist) were studied in more detail. cDC1s and cDC2s derived from cord blood stem cells, blood or patient tumor samples were stimulated with either TL8-506+IFN-γ or TL8-506+Poly(I:C). Different activation markers were analyzed by ELISA, flow cytometry, NanoString nCounter Technology or single-cell RNA-sequencing. T cell activation and migration assays were performed to assess functional consequences of cDC activation., Results: We show that TL8-506 synergized with IFN-γ or Poly(I:C) to induce high expression of different chemokines and cytokines including interleukin (IL)-12p70 in human cord blood and blood cDC subsets in a combination-specific manner. Importantly, both combinations induced the activation of cDC subsets in patient tumor samples ex vivo. The expression of immunostimulatory genes important for anticancer responses including CD40 , IFNB1 , IFNL1 , IL12A and IL12B were upregulated on stimulation. Furthermore, chemokines associated with CD8 + T cell recruitment were induced in tumor-derived cDCs in response to TL8-506 combinations. In vitro activation and migration assays confirmed that stimulated cDCs induce T cell activation and migration., Conclusions: Our data suggest that cord blood-derived and blood-derived cDCs are a good surrogate to study treatment responses in human tumor cDCs. While most cDCs in human tumors display a non-activated phenotype, TL8-506 combinations drive human tumor cDCs towards an immunostimulatory phenotype associated with Th1 responses on stimulation. Hence, TL8-506-based combinations may be promising candidates to initiate or boost antitumor responses in patients with cancer., Competing Interests: Competing interests: All authors, except MPL, RD and MvdB are employees of Roche., (© Author(s) (or their employer(s)) 2022. Re-use permitted under CC BY-NC. No commercial re-use. See rights and permissions. Published by BMJ.)

Single-cell RNA sequencing (scRNA-seq) revolutionized our understanding of disease biology. The promise it presents to also transform translational research requires highly standardized and robust software workflows. Here, we present the toolkit Besca , which streamlines scRNA-seq analyses and their use to deconvolute bulk RNA-seq data according to current best practices. Beyond a standard workflow covering quality control, filtering, and clustering, two complementary Besca modules, utilizing hierarchical cell signatures and supervised machine learning, automate cell annotation and provide harmonized nomenclatures. Subsequently, the gene expression profiles can be employed to estimate cell type proportions in bulk transcriptomics data. Using multiple, diverse scRNA-seq datasets, some stemming from highly heterogeneous tumor tissue, we show how Besca aids acceleration, interoperability, reusability and interpretability of scRNA-seq data analyses, meeting crucial demands in translational research and beyond., (© The Author(s) 2021. Published by Oxford University Press on behalf of NAR Genomics and Bioinformatics.)

The recruitment of thermogenic brite adipocytes within white adipose tissue attenuates obesity and metabolic comorbidities, arousing interest in understanding the underlying regulatory mechanisms. The molecular network of brite adipogenesis, however, remains largely unresolved. In this light, long noncoding RNAs (lncRNAs) emerged as a versatile class of modulators that control many steps within the differentiation machinery. Leveraging the naturally varying propensities of different inbred mouse strains for white adipose tissue browning, we identify the nuclear lncRNA Ctcflos as a pivotal orchestrator of thermogenic gene expression during brite adipocyte differentiation. Mechanistically, Ctcflos acts as a pleiotropic regulator, being essential for the transcriptional recruitment of the early core thermogenic regulatory program and the modulation of alternative splicing to drive brite adipogenesis. This is showcased by Ctcflos-regulated gene transcription and splicing of the key browning factor Prdm16 toward the isoform that is specific for the thermogenic gene program. Conclusively, our findings emphasize the mechanistic versatility of lncRNAs acting at several independent levels of gene expression for effective regulation of key differentiation factors to direct cell fate and function., (© 2021 The Authors. Published under the terms of the CC BY 4.0 license.)

Recruitment of brite/beige cells, known as browning of white adipose tissue (WAT), is an efficient way to turn an energy-storing organ into an energy-dissipating one and may therefore be of therapeutic value in combating obesity. However, a comprehensive understanding of the regulatory mechanisms mediating WAT browning is still lacking. Here, we exploit the large natural variation in WAT browning propensity between inbred mouse strains to gain an inclusive view of the core regulatory network coordinating this cellular process. Combining comparative transcriptomics, perturbation-based validations, and gene network analyses, we present a comprehensive gene regulatory network of inguinal WAT browning, revealing up to four distinct regulatory modules with key roles for uncovered transcriptional factors, while also providing deep insights into the genetic architecture of brite adipogenesis. The presented findings therefore greatly increase our understanding of the molecular drivers mediating the intriguing cellular heterogeneity and plasticity of adipose tissue., (Copyright © 2019 The Author(s). Published by Elsevier Inc. All rights reserved.)

Pancreatic beta cells synthesize and secrete the neurotransmitter γ-aminobutyric acid (GABA) as a paracrine and autocrine signal to help regulate hormone secretion and islet homeostasis. Islet GABA release has classically been described as a secretory vesicle-mediated event. Yet, a limitation of the hypothesized vesicular GABA release from islets is the lack of expression of a vesicular GABA transporter in beta cells. Consequentially, GABA accumulates in the cytosol. Here we provide evidence that the human beta cell effluxes GABA from a cytosolic pool in a pulsatile manner, imposing a synchronizing rhythm on pulsatile insulin secretion. The volume regulatory anion channel (VRAC), functionally encoded by LRRC8A or Swell1, is critical for pulsatile GABA secretion. GABA content in beta cells is depleted and secretion is disrupted in islets from type 1 and type 2 diabetic patients, suggesting that loss of GABA as a synchronizing signal for hormone output may correlate with diabetes pathogenesis., Competing Interests: COMPETING INTERESTS STATEMENT The authors declare no competing interests.

Krüppel-associated box zinc finger proteins (KZFPs) constitute the largest family of mammalian transcription factors, but most remain completely uncharacterized. While initially proposed to primarily repress transposable elements, recent reports have revealed that KFZPs contribute to a wide variety of other biological processes. Using murine and human in vitro and in vivo models, we demonstrate here that one poorly studied KZFP, ZFP30, promotes adipogenesis by directly targeting and activating a retrotransposon-derived Pparg2 enhancer. Through mechanistic studies, we further show that ZFP30 recruits the co-regulator KRAB-associated protein 1 (KAP1), which, surprisingly, acts as a ZFP30 co-activator in this adipogenic context. Our findings provide an understanding of both adipogenic and KZFP-KAP1 complex-mediated gene regulation, showing that the KZFP-KAP1 axis can also function in a non-repressive manner.

Adipose tissue in the mammary gland undergoes dramatic remodeling during reproduction. Adipocytes are replaced by mammary alveolar structures during pregnancy and lactation, then reappear upon weaning. The fate of the original adipocytes during lactation and the developmental origin of the re-appearing adipocyte post involution are unclear. Here, we reveal that adipocytes in the mammary gland de-differentiate into Pdgfrα + preadipocyte- and fibroblast-like cells during pregnancy and remain de-differentiated during lactation. Upon weaning, de-differentiated fibroblasts proliferate and re-differentiate into adipocytes. This cycle occurs over multiple pregnancies. These observations reveal the potential of terminally differentiated adipocytes to undergo repeated cycles of de-differentiation and re-differentiation in a physiological setting., (Copyright © 2018 Elsevier Inc. All rights reserved.)

Adipocyte development and differentiation have an important role in the aetiology of obesity and its co-morbidities 1,2 . Although multiple studies have investigated the adipogenic stem and precursor cells that give rise to mature adipocytes 3-14 , our understanding of their in vivo origin and properties is incomplete 2,15,16 . This is partially due to the highly heterogeneous and unstructured nature of adipose tissue depots 17 , which has proven difficult to molecularly dissect using classical approaches such as fluorescence-activated cell sorting and Cre-lox lines based on candidate marker genes 16,18 . Here, using the resolving power of single-cell transcriptomics 19 in a mouse model, we reveal distinct subpopulations of adipose stem and precursor cells in the stromal vascular fraction of subcutaneous adipose tissue. We identify one of these subpopulations as CD142 + adipogenesis-regulatory cells, which can suppress adipocyte formation in vivo and in vitro in a paracrine manner. We show that adipogenesis-regulatory cells are refractory to adipogenesis and that they are functionally conserved in humans. Our findings point to a potentially critical role for adipogenesis-regulatory cells in modulating adipose tissue plasticity, which is linked to metabolic control, differential insulin sensitivity and type 2 diabetes.

A long-standing question in biology is whether multipotent somatic stem and progenitor cells (SSPCs) feature molecular properties that could guide their system-independent identification. Population-based transcriptomic studies have so far not been able to provide a definite answer, given the rarity and heterogeneous nature of these cells. Here, we exploited the resolving power of single-cell RNA-sequencing to develop a computational model that is able to accurately distinguish SSPCs from differentiated cells across tissues. The resulting classifier is based on the combined expression of 23 genes including known players in multipotency, proliferation, and tumorigenesis, as well as novel ones, such as Lcp1 and Vgll4 that we functionally validate in intestinal organoids. We show how this approach enables the identification of stem-like cells in still ambiguous systems such as the pancreas and the epidermis as well as the exploration of lineage commitment hierarchies, thus facilitating the study of biological processes such as cellular differentiation, tissue regeneration, and cancer. Stem Cells 2017;35:2390-2402., (© 2017 The Authors STEM CELLS published by Wiley Periodicals, Inc. on behalf of AlphaMed Press.)

Motivation: Single-cell RNA-sequencing (scRNA-seq) allows whole transcriptome profiling of thousands of individual cells, enabling the molecular exploration of tissues at the cellular level. Such analytical capacity is of great interest to many research groups in the world, yet these groups often lack the expertise to handle complex scRNA-seq datasets., Results: We developed a fully integrated, web-based platform aimed at the complete analysis of scRNA-seq data post genome alignment: from the parsing, filtering and normalization of the input count data files, to the visual representation of the data, identification of cell clusters, differentially expressed genes (including cluster-specific marker genes), and functional gene set enrichment. This Automated Single-cell Analysis Pipeline (ASAP) combines a wide range of commonly used algorithms with sophisticated visualization tools. Compared with existing scRNA-seq analysis platforms, researchers (including those lacking computational expertise) are able to interact with the data in a straightforward fashion and in real time. Furthermore, given the overlap between scRNA-seq and bulk RNA-seq analysis workflows, ASAP should conceptually be broadly applicable to any RNA-seq dataset. As a validation, we demonstrate how we can use ASAP to simply reproduce the results from a single-cell study of 91 mouse cells involving five distinct cell types., Availability and Implementation: The tool is freely available at asap.epfl.ch and R/Python scripts are available at github.com/DeplanckeLab/ASAP., Contact: bart.deplancke@epfl.ch., Supplementary Information: Supplementary data are available at Bioinformatics online., (© The Author(s) 2017. Published by Oxford University Press.)

Brown adipocytes regulate energy expenditure via mitochondrial uncoupling, which makes them attractive therapeutic targets to tackle obesity. However, the regulatory mechanisms underlying brown adipogenesis are still poorly understood. To address this, we profiled the transcriptome and chromatin state during mouse brown fat cell differentiation, revealing extensive gene expression changes and chromatin remodeling, especially during the first day post-differentiation. To identify putatively causal regulators, we performed transcription factor binding site overrepresentation analyses in active chromatin regions and prioritized factors based on their expression correlation with the bona-fide brown adipogenic marker Ucp1 across multiple mouse and human datasets. Using loss-of-function assays, we evaluated both the phenotypic effect as well as the transcriptomic impact of several putative regulators on the differentiation process, uncovering ZFP467, HOXA4 and Nuclear Factor I A (NFIA) as novel transcriptional regulators. Of these, NFIA emerged as the regulator yielding the strongest molecular and cellular phenotypes. To examine its regulatory function, we profiled the genomic localization of NFIA, identifying it as a key early regulator of terminal brown fat cell differentiation., Competing Interests: The authors declare no competing financial interests.

Adipose tissue is a key determinant of whole body metabolism and energy homeostasis. Unraveling the regulatory mechanisms underlying adipogenesis is therefore highly relevant from a biomedical perspective. Our current understanding of fat cell differentiation is centered on the transcriptional cascades driven by the C/EBP protein family and the master regulator PPARγ. To elucidate further components of the adipogenic gene regulatory network, we performed a large-scale transcription factor (TF) screen overexpressing 734 TFs in mouse pre-adipocytes and probed their effect on differentiation. We identified 22 novel pro-adipogenic TFs and characterized the top ranking TF, ZEB1, as being essential for adipogenesis both in vitro and in vivo. Moreover, its expression levels correlate with fat cell differentiation potential in humans. Genomic profiling further revealed that this TF directly targets and controls the expression of most early and late adipogenic regulators, identifying ZEB1 as a central transcriptional component of fat cell differentiation.

Background: Alterations in glucose metabolism and epithelial-mesenchymal transition (EMT) constitute two important characteristics of carcinoma progression toward invasive cancer. Despite an extensive characterization of each of them separately, the links between EMT and glucose metabolism of tumor cells remain elusive. Here we show that the neuronal glucose transporter GLUT3 contributes to glucose uptake and proliferation of lung tumor cells that have undergone an EMT., Results: Using a panel of human non-small cell lung cancer (NSCLC) cell lines, we demonstrate that GLUT3 is strongly expressed in mesenchymal, but not epithelial cells, a finding corroborated in hepatoma cells. Furthermore, we identify that ZEB1 binds to the GLUT3 gene to activate transcription. Importantly, inhibiting GLUT3 expression reduces glucose import and the proliferation of mesenchymal lung tumor cells, whereas ectopic expression in epithelial cells sustains proliferation in low glucose. Using a large microarray data collection of human NSCLCs, we determine that GLUT3 expression correlates with EMT markers and is prognostic of poor overall survival., Conclusions: Altogether, our results reveal that GLUT3 is a transcriptional target of ZEB1 and that this glucose transporter plays an important role in lung cancer, when tumor cells loose their epithelial characteristics to become more invasive. Moreover, these findings emphasize the development of GLUT3 inhibitory drugs as a targeted therapy for the treatment of patients with poorly differentiated tumors.

Background: The genomic binding of CTCF is highly conserved across mammals, but the mechanisms that underlie its stability are poorly understood. One transcription factor known to functionally interact with CTCF in the context of X-chromosome inactivation is the ubiquitously expressed YY1. Because combinatorial transcription factor binding can contribute to the evolutionary stabilization of regulatory regions, we tested whether YY1 and CTCF co-binding could in part account for conservation of CTCF binding., Results: Combined analysis of CTCF and YY1 binding in lymphoblastoid cell lines from seven primates, as well as in mouse and human livers, reveals extensive genome-wide co-localization specifically at evolutionarily stable CTCF-bound regions. CTCF-YY1 co-bound regions resemble regions bound by YY1 alone, as they enrich for active histone marks, RNA polymerase II and transcription factor binding. Although these highly conserved, transcriptionally active CTCF-YY1 co-bound regions are often promoter-proximal, gene-distal regions show similar molecular features., Conclusions: Our results reveal that these two ubiquitously expressed, multi-functional zinc-finger proteins collaborate in functionally active regions to stabilize one another's genome-wide binding across primate evolution.

Chromatin organization and gene activity are responsive to developmental and environmental cues. Although many genes are transcribed throughout development and across cell types, much of gene regulation is highly cell-type specific. To readily track chromatin features at the resolution of cell types within complex tissues, we developed and validated chromatin affinity purification from specific cell types by chromatin immunoprecipitation (CAST-ChIP), a broadly applicable biochemical procedure. RNA polymerase II (Pol II) CAST-ChIP identifies ~1,500 neuronal and glia-specific genes in differentiated cells within the adult Drosophila brain. In contrast, the histone H2A.Z is distributed similarly across cell types and throughout development, marking cell-type-invariant Pol II-bound regions. Our study identifies H2A.Z as an active chromatin signature that is refractory to changes across cell fates. Thus, CAST-ChIP powerfully identifies cell-type-specific as well as cell-type-invariant chromatin states, enabling the systematic dissection of chromatin structure and gene regulation within complex tissues such as the brain., (Copyright © 2013 The Authors. Published by Elsevier Inc. All rights reserved.)

At least half of the human genome is derived from repetitive elements, which are often lineage specific and silenced by a variety of genetic and epigenetic mechanisms. Using a transchromosomic mouse strain that transmits an almost complete single copy of human chromosome 21 via the female germline, we show that a heterologous regulatory environment can transcriptionally activate transposon-derived human regulatory regions. In the mouse nucleus, hundreds of locations on human chromosome 21 newly associate with activating histone modifications in both somatic and germline tissues, and influence the gene expression of nearby transcripts. These regions are enriched with primate and human lineage-specific transposable elements, and their activation corresponds to changes in DNA methylation at CpG dinucleotides. This study reveals the latent regulatory potential of the repetitive human genome and illustrates the species specificity of mechanisms that control it., (Copyright © 2013 Elsevier Inc. All rights reserved.)

The cohesin protein complex contributes to transcriptional regulation in a CTCF-independent manner by colocalizing with master regulators at tissue-specific loci. The regulation of transcription involves the concerted action of multiple transcription factors (TFs) and cohesin's role in this context of combinatorial TF binding remains unexplored. To investigate cohesin-non-CTCF (CNC) binding events in vivo we mapped cohesin and CTCF, as well as a collection of tissue-specific and ubiquitous transcriptional regulators using ChIP-seq in primary mouse liver. We observe a positive correlation between the number of distinct TFs bound and the presence of CNC sites. In contrast to regions of the genome where cohesin and CTCF colocalize, CNC sites coincide with the binding of master regulators and enhancer-markers and are significantly associated with liver-specific expressed genes. We also show that cohesin presence partially explains the commonly observed discrepancy between TF motif score and ChIP signal. Evidence from these statistical analyses in wild-type cells, and comparisons to maps of TF binding in Rad21-cohesin haploinsufficient mouse liver, suggests that cohesin helps to stabilize large protein-DNA complexes. Finally, we observe that the presence of mirrored CTCF binding events at promoters and their nearby cohesin-bound enhancers is associated with elevated expression levels.

Gorillas are humans' closest living relatives after chimpanzees, and are of comparable importance for the study of human origins and evolution. Here we present the assembly and analysis of a genome sequence for the western lowland gorilla, and compare the whole genomes of all extant great ape genera. We propose a synthesis of genetic and fossil evidence consistent with placing the human-chimpanzee and human-chimpanzee-gorilla speciation events at approximately 6 and 10 million years ago. In 30% of the genome, gorilla is closer to human or chimpanzee than the latter are to each other; this is rarer around coding genes, indicating pervasive selection throughout great ape evolution, and has functional consequences in gene expression. A comparison of protein coding genes reveals approximately 500 genes showing accelerated evolution on each of the gorilla, human and chimpanzee lineages, and evidence for parallel acceleration, particularly of genes involved in hearing. We also compare the western and eastern gorilla species, estimating an average sequence divergence time 1.75 million years ago, but with evidence for more recent genetic exchange and a population bottleneck in the eastern species. The use of the genome sequence in these and future analyses will promote a deeper understanding of great ape biology and evolution.

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., (Copyright © 2012 Elsevier Inc. All rights reserved.)

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.

The cohesin protein complex holds sister chromatids in dividing cells together and is essential for chromosome segregation. Recently, cohesin has been implicated in mediating transcriptional insulation, via its interactions with CTCF. Here, we show in different cell types that cohesin functionally behaves as a tissue-specific transcriptional regulator, independent of CTCF binding. By performing matched genome-wide binding assays (ChIP-seq) in human breast cancer cells (MCF-7), we discovered thousands of genomic sites that share cohesin and estrogen receptor alpha (ER) yet lack CTCF binding. By use of human hepatocellular carcinoma cells (HepG2), we found that liver-specific transcription factors colocalize with cohesin independently of CTCF at liver-specific targets that are distinct from those found in breast cancer cells. Furthermore, estrogen-regulated genes are preferentially bound by both ER and cohesin, and functionally, the silencing of cohesin caused aberrant re-entry of breast cancer cells into cell cycle after hormone treatment. We combined chromosomal interaction data in MCF-7 cells with our cohesin binding data to show that cohesin is highly enriched at ER-bound regions that capture inter-chromosomal loop anchors. Together, our data show that cohesin cobinds across the genome with transcription factors independently of CTCF, plays a functional role in estrogen-regulated transcription, and may help to mediate tissue-specific transcriptional responses via long-range chromosomal interactions.

When taking their blood meal on the mammalian host, ticks transfer a multitude of different proteins from their saliva into the host. Some of these proteins are hijacked by pathogens for their own purposes. Borrelia burgdorferi, the Lyme disease agent, is critically dependent on the presence of the tick protein Salp15 when infecting the host. Similarly, Anaplasma phagocytophilum, which causes anaplasmosis, needs Salp16, a homologue of Salp15, to get transferred from the host into the tick. Here we analyzed whether adaptive evolution has shaped the Salp15 protein family. Using site-specific estimates of K(A)/K(S) ratios, we identified different positions within the Salp15 protein family which have undergone a phase of positive selection. Additionally, we analyzed the B. burgdorferi protein interacting with Salp15, OspC. Again, sites showing signs of positive selection were identified, although they are more likely a result of the antigenic features of OspC than of the influence of Salp15. The identification of probably functionally relevant sites in the Salp15 family might direct the detailed experimental analysis of their interaction with human and bacterial proteins.