13 research outputs found
Single-Cell Epigenomics and Functional Fine-Mapping of Atherosclerosis GWAS Loci
Rationale: Genome-wide association studies have identified hundreds of loci associated with coronary artery disease (CAD). Many of these loci are enriched in cisregulatory elements but not linked to cardiometabolic risk factors nor to candidate causal genes, complicating their functional interpretation. Objective: Single-nucleus chromatin accessibility profiling of the human atherosclerotic lesions was used to investigate cell type-specific patterns of cisregulatory elements, to understand transcription factors establishing cell identity, and to interpret CAD-relevant, noncoding genetic variation. Methods and Results: We used single-nucleus ATAC-seq (assay for transposase-accessible chromatin with sequencing) to generate DNA accessibility maps in >7000 cells derived from human atherosclerotic lesions. We identified 5 major lesional cell types including endothelial cells, smooth muscle cells, monocyte/macrophages, natural killer/T cells, and B cells and further investigated subtype characteristics of macrophages and smooth muscle cells transitioning into fibromyocytes. We demonstrated that CAD-associated genetic variants are particularly enriched in endothelial and smooth muscle cell-specific open chromatin. Using single-cell coaccessibility and cis-expression quantitative trait loci information, we prioritized putative target genes and candidate regulatory elements for approximate to 30% of all known CAD loci. Finally, we performed genome-wide experimental fine-mapping of the CAD variants identified in genome-wide association studies using epigenetic quantitative trait loci analysis in primary human aortic endothelial cells and self-transcribing active regulatory region sequencing (STARR-Seq) massively parallel reporter assay in smooth muscle cells. This analysis identified potential causal single-nucleotide polymorphisms (SNPs) and the associated target gene for over 30 CAD loci. We present several examples where the chromatin accessibility and gene expression could be assigned to one cell type predicting the cell type of action for CAD loci. Conclusions: These findings highlight the potential of applying single-nucleus ATAC-seq to human tissues in revealing relative contributions of distinct cell types to diseases and in identifying genes likely to be influenced by noncoding genome-wide association study variants.</p
MGEx-Udb: A Mammalian Uterus Database for Expression-Based Cataloguing of Genes across Conditions, Including Endometriosis and Cervical Cancer
Gene expression profiling of uterus tissue has been performed in various contexts, but a significant amount of the data remains underutilized as it is not covered by the existing general resources.). The database can be queried with gene names/IDs, sub-tissue locations, as well as various conditions such as the cervical cancer, endometrial cycles and disorders, and experimental treatments. Accordingly, the output would be a) transcribed and dormant genes listed for the queried condition/location, or b) expression profile of the gene of interest in various uterine conditions. The results also include the reliability score for the expression status of each gene. MGEx-Udb also provides information related to Gene Ontology annotations, protein-protein interactions, transcripts, promoters, and expression status by other sequencing techniques, and facilitates various other types of analysis of the individual genes or co-expressed gene clusters.In brief, MGEx-Udb enables easy cataloguing of co-expressed genes and also facilitates bio-marker discovery for various uterine conditions
Human Ezrin: Prediction of the 3D structure and interactions with mTOR
The protein Ezrin, is a member of the ERM family (Ezrin, Radixin and Moesin) that links the F-actin to the plasma membrane. The protein is made of three domains namely the FERM domain, a central α-helical domain and the CERMAD domain. The residues in Ezrin such as Ser66, Tyr145, Tyr353 and Tyr477 regulate the function of the protein through phosphorylation. The protein is found in two distinct conformations of active and dormant (inactive) state. The initial step during the conformation change is the breakage of intramolecular interaction in dormant Ezrin by phosphorylation of residue Thr567. The dormant structure of human Ezrin was predicted computationally since only partial active form structure was available. The validation analysis showed that 99.7% residues were positioned in favored, allowed and generously allowed regions of the Ramachandran plot. The Z-score of Ezrin was â7.36, G-factor was 0.1, and the QMEAN score of the model was 0.61 indicating a good model for human Ezrin. The comparison of the conformations of the activated and dormant Ezrin showed a major shift in the F2 lobe (residues 142-149 and 161-177) while changes in the conformation induced mobility shifts in lobe F3 (residues 261 to 267). The 3D positions of the phosphorylation sites Tyr145, Tyr353, Tyr477, Tyr482 and Thr567 were also located. Using targeted molecular dynamic simulation, the molecular movements during conformational change from active to dormant were visualized. The dormant Ezrin auto-inhibits itself by a head-to-tail interaction of the N-terminal and C-terminal residues. The trajectory shows the breakage of the interactions and mobility of the CERMAD domain away from the FERM domain. Protein docking and clustering analysis were used to predict the residues involved in the interaction between dormant Ezrin and mTOR. Residues Tyr477 and Tyr482 were found to be involved in interaction with mTOR.Siirretty Doriast
Schematic representation of MGEx-Udb.
<p>The figure represents the data collection (top portion), architecture (central portion) and operation (bottom portion) of the database.</p
Datasets
<p>(<b>with gene count</b>) <b>collected from various sources.</b> In case of âPubMed & GEOâ and âPubMed & ArrayExpressâ, smaller gene lists came from validation experiments and were collected from PubMed, while raw/processed data were always collected from the repositories (GEO/ArrayExpress).</p
Example hierarchy of the conditions and sub-conditions.
<p>An example <i>(âstage IIA non-keratinizing squamous cell cervical carcinomaâ)</i> hierarchy of the conditions and sub-conditions, for which data have been collected, and drop-down options provided in the query and upload pages of MGEx-Udb. Currently the database allows up to four levels of the hierarchy to query.</p
Number of datasets
<p>(<b>and studies</b>) <b>in MGEx-Udb corresponding to various physiological and pathological uterine conditions.</b> âOthersâ represent post-parturition, genetic-ablation, artificial insemination and embryo implantation. Studies considering tissues that are used as controls but may not be absolutely ânormalâ have been grouped in <i>âmay be normalâ</i> category (examples: ânormal tissue adjacent to tumor/cancer tissueâ, âvehicle-treatedâ).</p
Source of data across various mammalian species in MGEx-Udb.
<p>Other species include cow and pig. Among the data collected from GEO or âPubMed & GEOâ, 85% of the studies were also present in ArrayExpress, even though this is not indicated in the figure.</p
Oncogenic KEAP1 mutations activate TRAF2-NFÎșB signaling to prevent apoptosis in lung cancer cells
The Kelch-like ECH-associated protein 1 (KEAP1) - Nuclear factor erythroid 2 -related factor 2 (NRF2) pathway is the major transcriptional stress response system in cells against oxidative and electrophilic stress. NRF2 is frequently constitutively active in many cancers, rendering the cells resistant to chemo- and radiotherapy. Lossof-function (LOF) mutations in the repressor protein KEAP1 are common in non-small cell lung cancer, particularly adenocarcinoma. While the mutations can occur throughout the gene, they are enriched in certain areas, indicating that these may have unique functional importance. In this study, we show that in the GSEA analysis of TCGA lung adenocarcinoma RNA-seq data, the KEAP1 mutations in R320 and R470 were associated with enhanced Tumor Necrosis Factor alpha (TNF alpha) - Nuclear Factor kappa subunit B (NF kappa B) signaling as well as MYC and MTORC1 pathways. To address the functional role of these hotspot mutations, affinity purification and mass spectrometry (AP-MS) analysis of wild type (wt) KEAP1 and its mutation forms, R320Q and R470C were employed to interrogate differences in the protein interactome. We identified TNF receptor associated factor 2 (TRAF2) as a putative protein interaction partner. Both mutant KEAP1 forms showed increased interaction with TRAF2 and other anti-apoptotic proteins, suggesting that apoptosis signalling could be affected by the protein interactions. A549 lung adenocarcinoma cells overexpressing mutant KEAP1 showed high TRAF2-mediated NF kappa B activity and increased protection against apoptosis, XIAP being one of the key proteins involved in anti-apoptotic signalling. To conclude, KEAP1 R320Q and R470C and its interaction with TRAF2 leads to activation of NF kappa B pathway, thereby protecting against apoptosis.Peer reviewe
A pan-cancer analysis shows immunoevasive characteristics in NRF2 hyperactive squamous malignancies
The NRF2 pathway is frequently activated in various cancer types, yet a comprehensive analysis of its effects across different malignancies is currently lacking. We developed a NRF2 activity metric and utilized it to conduct a pan-cancer analysis of oncogenic NRF2 signaling. We identified an immunoevasive phenotype where high NRF2 activity is associated with low interferon-gamma (IFNÎł), HLA-I expression and T cell and macrophage infiltration in squamous malignancies of the lung, head and neck area, cervix and esophagus. Squamous NRF2 overactive tumors comprise a molecular phenotype with SOX2/TP63 amplification, TP53 mutation and CDKN2A loss. These immune cold NRF2 hyperactive diseases are associated with upregulation of immunomodulatory NAMPT, WNT5A, SPP1, SLC7A11, SLC2A1 and PD-L1. Based on our functional genomics analyses, these genes represent candidate NRF2 targets, suggesting direct modulation of the tumor immune milieu. Single-cell mRNA data shows that cancer cells of this subtype exhibit decreased expression of IFNÎł responsive ligands, and increased expression of immunosuppressive ligands NAMPT, SPP1 and WNT5A that mediate signaling in intercellular crosstalk. In addition, we discovered that the negative relationship of NRF2 and immune cells are explained by stromal populations of lung squamous cell carcinoma, and this effect spans multiple squamous malignancies based on our molecular subtyping and deconvolution data