scholarly journals Testing Proximity of Genomic Regions to Transcription Start Sites and Enhancers Complements Gene Set Enrichment Testing

2020 ◽  
Vol 11 ◽  
Author(s):  
Christopher Lee ◽  
Kai Wang ◽  
Tingting Qin ◽  
Maureen A. Sartor
2016 ◽  
Author(s):  
Francisco Avila Cobos ◽  
Jasper Anckaert ◽  
Pieter-Jan Volders ◽  
Dries Rombaut ◽  
Jo Vandesompele ◽  
...  

AbstractSummaryReconstructing transcript models from RNA-sequencing (RNA-seq) data and establishing these as independent transcriptional units can be a challenging task. The Zipper plot is an application that enables users to interrogate putative transcription start sites (TSSs) in relation to various features that are indicative for transcriptional activity. These features are obtained from publicly available datasets including CAGE-sequencing (CAGE-seq), ChIP-sequencing (ChIP-seq) for histone marks and DNasesequencing (DNase-seq). The Zipper plot application requires three input fields (chromosome, genomic coordinate (hg19) of the TSS and strand) and generates a report that includes a detailed summary table, a Zipper plot and several statistics derived from this plot.Availability and ImplementationThe Zipper plot is implemented using the statistical programming language R and is freely available at http://[email protected]; [email protected]; [email protected] informationSupplementary Methods available online.


2021 ◽  
Vol 22 (1) ◽  
Author(s):  
Simon Bourdareau ◽  
Leila Tirichine ◽  
Bérangère Lombard ◽  
Damarys Loew ◽  
Delphine Scornet ◽  
...  

Abstract Background Brown algae evolved complex multicellularity independently of the animal and land plant lineages and are the third most developmentally complex phylogenetic group on the planet. An understanding of developmental processes in this group is expected to provide important insights into the evolutionary events necessary for the emergence of complex multicellularity. Here, we focus on mechanisms of epigenetic regulation involving post-translational modifications of histone proteins. Results A total of 47 histone post-translational modifications are identified, including a novel mark H2AZR38me1, but Ectocarpus lacks both H3K27me3 and the major polycomb complexes. ChIP-seq identifies modifications associated with transcription start sites and gene bodies of active genes and with transposons. H3K79me2 exhibits an unusual pattern, often marking large genomic regions spanning several genes. Transcription start sites of closely spaced, divergently transcribed gene pairs share a common nucleosome-depleted region and exhibit shared histone modification peaks. Overall, patterns of histone modifications are stable through the life cycle. Analysis of histone modifications at generation-biased genes identifies a correlation between the presence of specific chromatin marks and the level of gene expression. Conclusions The overview of histone post-translational modifications in the brown alga presented here will provide a foundation for future studies aimed at understanding the role of chromatin modifications in the regulation of brown algal genomes.


2021 ◽  
Author(s):  
Roman Hillje ◽  
Lucilla Luzi ◽  
Stefano Amatori ◽  
Mirco Fanelli ◽  
Pier Giuseppe Pelicci ◽  
...  

Abstract To disclose the epigenetic drift of time passing, we determined the genome-wide distributions of mono- and tri-methylated lysine 4 and acetylated and tri-methylated lysine 27 of histone H3 in the livers of healthy 3, 6 and 12 months old C57BL/6 mice. The comparison of different age profiles of histone H3 marks revealed global redistribution of histone H3 modifications with time, in particular in intergenic regions and near transcription start sites, as well as altered correlation between the profiles of different histone modifications. Moreover, feeding mice with caloric restriction diet, a treatment known to retard aging, preserved younger state of histone H3 in these genomic regions.


Blood ◽  
2018 ◽  
Vol 132 (Supplement 1) ◽  
pp. 2717-2717 ◽  
Author(s):  
Martin G. Klatt ◽  
Sung S. Mun ◽  
Nicholas D. Socci ◽  
Tatyana Korontsvit ◽  
Tao Dao ◽  
...  

Abstract Acute myeloid leukemia (AML) is an aggressive hematological malignancy with a 5-year overall survival rate of less than 30% which causes over 10,000 deaths per year in the United States. Treatment options for this disease increasingly include epigenetic drugs, such as hypomethylating agents (e.g. decitabine) or histone deacetylase (HDAC) inhibitors (e.g. pracinostat) which can function via direct cytotoxic mechanisms and also through altered differentiation of AML blasts; immunomodulatory effects like reactivation and presentation of cancer testis antigens in context of human leukocyte antigen (HLA) complexes have been reported as well, which may result in clearance of cells via the adaptive immune system. However, the landscape of immunogenic T cell epitopes induced by these drugs might be even broader than reported since standard analyses only consider exonic protein sequences and do not take into account typically untranslated genomic regions. Recently, it has been shown that single and combination treatment of decitabine and pracinostat can induce cryptic transcription start sites in generally epigenetically repressed solitary long-terminal repeats (LTRs) of the LTR12C family which give rise to novel mRNAs and resulting protein variants. We hypothesized that the intronic parts of these gene products might provide a source of cryptic T cell epitopes with high immunogenic potential, which are induced through epigenetic drug treatment. To test this hypothesis, we treated 5 different AML cell lines (HL-60, U937, OCI-AML02, MOLM13, AML14) with (1) DMSO, (2) 500 nM decitabine or (3) a combination of 500 nM decitabine and 100 nM pracinostat for 72 hours to induce transcription of non-annotated transcription start sites. Subsequently, HLA class I complexes were immunopurified and peptides presented by these complexes isolated and analyzed by liquid chromatography-tandem mass spectrometry (LC-MS/MS). The activation of silenced genes by epigenetic drug treatment with either decitabine alone or the combination treatment yielded increases of about two-fold in the identified unique HLA ligands. This increase in peptide identifications also led to improved detection of cancer testis antigen-derived epitopes, as has been reported before. Intriguingly, by adding LTR12C derived sequences stretching from the published GATA2 specific binding site until the next genomic exon to the peptide search analyses we were able to identify several cryptic peptides from 4 out of 5 AML cell lines derived from these usually untranscribed genomic regions. The identifications were exclusively dependent on previous treatment with either decitabine alone or in combination with pracinostat. Though the immunogenicity of these HLA ligands has not been determined yet, we assume that due to their genetically repressed state in untreated cells, these new peptide sequences represent a new class of neoepitopes, with potential to be novel targets of existing T cells within patients or after augmentation by other immunotherapies. In summary, we demonstrated for the first time the induced presentation of epitopes from normally untranscribed LTR12C regions through epigenetic drug treatment and therefore provide a previously undescribed source of potential targets for immunotherapy in AML. Disclosures Scheinberg: Eureka: Consultancy; Ensyce: Consultancy.


Epigenomics ◽  
2021 ◽  
Author(s):  
Marika Groleau ◽  
Frédérique White ◽  
Andres Cardenas ◽  
Patrice Perron ◽  
Marie-France Hivert ◽  
...  

Aim: The placenta undergoes DNA methylation (DNAm) programming that is unique compared with all other fetal tissues. We aim to decipher some of the physiologic roles of the placenta by comparing its DNAm profile with that of another fetal tissue. Materials & methods: We performed a comparative analysis of genome-wide DNAm of 444 placentas paired with cord blood samples collected at birth. Gene ontology term analyses were conducted on the resulting differentially methylated regions. Results: Genomic regions upstream of transcription start sites showing lower DNAm in the placenta were enriched with terms related to miRNA functions and genes encoding G protein-coupled receptors. Conclusion: These results highlight genomic regions that are differentially methylated in the placenta in contrast to fetal blood.


2020 ◽  
Vol 2 (1) ◽  
Author(s):  
Christopher T Lee ◽  
Raymond G Cavalcante ◽  
Chee Lee ◽  
Tingting Qin ◽  
Snehal Patil ◽  
...  

Abstract Gene set enrichment (GSE) testing enhances the biological interpretation of ChIP-seq data and other large sets of genomic regions. Our group has previously introduced two GSE methods for genomic regions: ChIP-Enrich for narrow regions and Broad-Enrich for broad regions. Here, we introduce Poly-Enrich, which has wider applicability, additional capabilities and models the number of peaks assigned to a gene using a generalized additive model with a negative binomial family to determine gene set enrichment, while adjusting for gene locus length. As opposed to ChIP-Enrich, Poly-Enrich works well even when nearly all genes have a peak, illustrated by using Poly-Enrich to characterize pathways and types of genic regions enriched with different families of repetitive elements. By comparing Poly-Enrich and ChIP-Enrich results with ENCODE ChIP-seq data, we found that the optimal test depends more on the pathway being regulated than on properties of the transcription factors. Using known transcription factor functions, we discovered clusters of related biological processes consistently better modeled with Poly-Enrich. This suggests that the regulation of certain processes may be modified by multiple binding events, better modeled by a count-based method. Our new hybrid method automatically uses the optimal method for each gene set, with correct FDR-adjustment.


2018 ◽  
Author(s):  
Christopher T Lee ◽  
Raymond G Cavalcante ◽  
Chee Lee ◽  
Tingting Qin ◽  
Snehal Patil ◽  
...  

AbstractGene set enrichment (GSE) testing enhances the biological interpretation of ChIP-seq data and other large sets of genomic regions. Our group has previously introduced two GSE methods for genomic regions: ChIP-Enrich for narrow regions and Broad-Enrich for broad genomic regions, such as histone modifications. Here, we introduce new methods and extensions that more appropriately analyze sets of genomic regions with vastly different properties. First, we introduce Poly-Enrich, which models the number of peaks assigned to a gene using a generalized additive model with a negative binomial family to determine gene set enrichment, while adjusting for gene locus length (#bps associated with each gene). This is the first method that controls for locus length while accounting for the number of peaks per gene and variability among genes. We also introduce a flexible weighting approach to incorporate region scores, a hybrid enrichment approach, and support for new gene set databases and reference genomes/species.As opposed to ChIP-Enrich, Poly-Enrich works well even when nearly all genes have a peak. To illustrate this, we used Poly-Enrich to characterize the pathways and types of genic regions (introns, promoters, etc) enriched with different families of repetitive elements. By comparing ChIP-Enrich and Poly-Enrich results from ENCODE ChIP-seq data, we found that the optimal test depends more on the pathway being regulated than on the transcription factor or other properties of the dataset. Using known transcription factor functions, we discovered clusters of related biological processes consistently better modeled with either the binary score method (ChIP-Enrich) or count based method (Poly-Enrich). This suggests that the regulation of certain processes is more often modified by multiple binding events (count-based), while others tend to require only one (binary). Our new hybrid method handles this by automatically choosing the optimal method, with correct FDR-adjustment.Author SummaryAlthough every cell in our body contains the same DNA, our cells perform vastly different functions due to differences in how our genes are regulated. Certain regions of the genome are bound by DNA binding proteins (transcription factors), which regulate the expression of nearby genes. After an experiment to identify a large set of these regions, we can then model the association of these regions with various cellular pathways and biological processes. This analysis helps understand the overall biological effect that the binding events have on the cells. For example, if genes relating to apoptosis tend to have the transcription factor, Bcl-2, bind more often nearby, then Bcl-2 is likely to have a vital role in regulating apoptosis. The specifics of how to best perform this analysis is still being researched and depends on properties of the set of genomic regions. Here, we introduce a new, more flexible method that counts the number of occurrences per gene and models that in a sophisticated statistical test, and compare it to a previous method. We show that the optimal method depends on multiple factors, and the new method, Poly-Enrich, allows interesting findings in scenarios where the previous method failed.


2020 ◽  
Vol 8 (Suppl 3) ◽  
pp. A646-A647
Author(s):  
Max Meneveau ◽  
Pankaj Kumar ◽  
Kevin Lynch ◽  
Karlyn Pollack ◽  
Craig Slingluff

BackgroundVaccines are a promising therapeutic for patients with advanced cancer, but achieving robust T-cell responses remains a challenge. Melanoma-associated antigen-A3 (MAGE-A3) in combination with adjuvant AS15 (a formulation of Toll-Like-Receptor (TLR)-4 and 9 agonists and a saponin), induced systemic CD4+ T-cell responses in 50% of patients when given subcutaneously/intradermally. Little is known about the transcriptional landscape of the vaccine-site microenvironment (VSME) of patients with systemic T-cell responses versus those without. We hypothesized that patients with systemic T-cell responses to vaccination would exhibit increased immune activation in the VSME, higher dendritic cell (DC) activation/maturation, TLR-pathway activation, and enhanced Th1 signatures.MethodsBiopsies of the VSME were obtained from participants on the Mel55 clinical trial (NCT01425749) who were immunized with MAGE-A3/AS15. Biopsies were taken 8 days after immunization. T-cell response to MAGE-A3 was assessed in PBMC after in-vitro stimulation with recMAGE-A3, by IFNγ ELISPOT assay. Gene expression was assessed by RNAseq using DESeq2. Comparisons were made between immune-responders (IR), non-responders (NR), and normal skin controls. FDR p<0.01 was considered significant.ResultsFour IR, four NR, and three controls were evaluated. The 500 most variable genes were used for principal component analysis (PCA). Two IR samples were identified as outliers on PCA and excluded from further analysis. There were 882 differentially expressed genes (DEGs) in the IR group vs the NR group (figure 1A). Unsupervised clustering of the top 500 DEGs revealed clustering according to the experimental groups (figure 1B). Of the 10 most highly upregulated DEGs, 9 were immune-related (figure 1C). Gene-set enrichment analysis revealed that immune-related pathways were highly enriched in IRs vs NRs (figure 1D). CD4 and CD8 expression did not differ between IR and NR (figure 2A), though both were higher in IR compared to control. Markers of DC activation/maturation were higher in IR vs NR (figure 2B), as were several Th1 associated genes (figure 2C). Interestingly, markers of exhaustion were higher in IR v NR (figure 2D). Expression of numerous TLR-pathway genes was higher in IR vs NR, including MYD88, but not TICAM1 (figure 2E).Abstract 611 Figure 1Gene expression profiling of vaccine site samples from patients immunized with MAGE-A3/AS15. (A) Volcano plots showing the distribution of differentially expressed genes (DEGs) between immune responders (IR) and non-responders (NR), IR and control, and NR and control. (B) Heatmap of the top 500 most differentially expressed genes demonstrating hierarchical clustering of sequenced samples according to IR, NR, and control. (C) Table showing the 10 most highly up and down-regulated genes in IR compared to NR. 9 of the top 10 most highly up-regulated genes are related to the immune response. (D) Enrichment plots from a gene set enrichment analysis highlighting the upregulation of immune related pathways in IR compared to NR. Gene set enrichment data was generated from the Reactome gene set database and included all expressed genes. Significance was set at FDR p <0.01Abstract 611 Figure 2Expression of T-cell markers in IR vs NR vs Control samples in the vaccine site microenvironment (VSME). (A) T-cell markers showing similar expression in IR vs NR but higher expression in IR vs control. (B) Markers of dendritic cell activation and maturation in the VSME which are higher in IR vs control but not IR vs NR. (B) Transcription factors and genes associated with Th1/Th2 responses within the VSME. (D) Genes associated with T-cell exhaustion at the VSME. (E) Expression of TLR pathway genes in the VSME. Expression data is provided in terms of normalized counts. Bars demonstrate median and interquartile range. N=9. IR = immune responder, NR = non-responder, TLR = Toll-like Receptor. * = <0.01, ** < 0.001, *** <0.0001, **** < 0.00001ConclusionsThese findings suggest a unique immune-transcriptional landscape in the VSME is associated with circulating T-cell responses to immunization, with differences in DC activation/maturation, Th1 response, and TLR signaling. Thus, immunologic changes in the VSME are useful predictors of systemic immune response, and host factors that modulate immune-related signaling at the vaccine site may have concordant systemic effects on promoting or limiting immune responses to vaccines.Trial RegistrationSamples for this work were collected from patients enrolled on the Mel55 clinical trial NCT01425749.Ethics ApprovalThis work was completed after approval from the UVA institutional review board IRB-HSR# 15398.


Sign in / Sign up

Export Citation Format

Share Document