Publications

What is a Publication?
6 Publications visible to you, out of a total of 6

Abstract (Expand)

Most high throughput genomic data analysis pipelines currently rely on over-representation or gene set enrichment analysis (ORA/GSEA) approaches for functional analysis. In contrast, topology-based pathway analysis methods, which offer a more biologically informed perspective by incorporating interaction and topology information, have remained underutilized and inaccessible due to various limiting factors. These methods heavily rely on the quality of pathway topologies and often utilize predefined topologies from databases without assessing their correctness. To address these issues and make topology-aware pathway analysis more accessible and flexible, we introduce the PSF (Pathway Signal Flow) toolkit R package. Our toolkit integrates pathway curation and topology-based analysis, providing interactive and command-line tools that facilitate pathway importation, correction, and modification from diverse sources. This enables users to perform topology-based pathway signal flow analysis in both interactive and command-line modes. To showcase the toolkit’s usability, we curated 36 KEGG signaling pathways and conducted several use-case studies, comparing our method with ORA and the topology-based signaling pathway impact analysis (SPIA) method. The results demonstrate that the algorithm can effectively identify ORA enriched pathways while providing more detailed branch-level information. Moreover, in contrast to the SPIA method, it offers the advantage of being cut-off free and less susceptible to the variability caused by selection thresholds. By combining pathway curation and topology-based analysis, the PSF toolkit enhances the quality, flexibility, and accessibility of topology-aware pathway analysis. Researchers can now easily import pathways from various sources, correct and modify them as needed, and perform detailed topology-based pathway signal flow analysis. In summary, our PSF toolkit offers an integrated solution that addresses the limitations of current topology-based pathway analysis methods. By providing interactive and command-line tools for pathway curation and topology-based analysis, we empower researchers to conduct comprehensive pathway analyses across a wide range of applications.

Authors: Siras Hakobyan, Ani Stepanyan, Lilit Nersisyan, Hans Binder, Arsen Arakelyan

Date Published: 23rd Aug 2023

Publication Type: Journal

Abstract (Expand)

Regulation of messenger RNA stability is pivotal for programmed gene expression in bacteria and is achieved by a myriad of molecular mechanisms. By bulk sequencing of 5' monophosphorylated mRNA decay intermediates (5'P), we show that cotranslational mRNA degradation is conserved among both Gram-positive and -negative bacteria. We demonstrate that, in species with 5'-3' exonucleases, the exoribonuclease RNase J tracks the trailing ribosome to produce an in vivo single-nucleotide toeprint of the 5' position of the ribosome. In other species lacking 5'-3' exonucleases, ribosome positioning alters endonucleolytic cleavage sites. Using our metadegradome (5'P degradome) sequencing approach, we characterize 5'P mRNA decay intermediates in 96 species including Bacillus subtilis, Escherichia coli, Synechocystis spp. and Prevotella copri and identify codon- and gene-level ribosome stalling responses to stress and drug treatment. We also apply 5'P sequencing to complex clinical and environmental microbiomes and demonstrate that metadegradome sequencing provides fast, species-specific posttranscriptional characterization of responses to drug or environmental perturbations. Finally we produce a degradome atlas for 96 species to enable analysis of mechanisms of RNA degradation in bacteria. Our work paves the way for the application of metadegradome sequencing to investigation of posttranscriptional regulation in unculturable species and complex microbial communities.

Authors: S. Huch, L. Nersisyan, M. Ropat, D. Barrett, M. Wu, J. Wang, V. D. Valeriano, N. Vardazaryan, J. Huerta-Cepas, W. Wei, J. Du, L. M. Steinmetz, L. Engstrand, V. Pelechano

Date Published: 5th Jun 2023

Publication Type: Journal

Abstract (Expand)

The sequencing of SARS-CoV-2 provides essential information on viral evolution, transmission, and epidemiology. In this paper, we performed the whole-genome sequencing of SARS-CoV-2 using nanopore and Illumina sequencing to describe the circulation of the virus lineages in Armenia. The analysis of 145 full genomes identified six clades (19A, 20A, 20B, 20I, 21J, and 21K) and considerable intra-clade PANGO lineage diversity. Phylodynamic and transmission analysis allowed to attribute specific clades as well as infer their importation routes. Thus, the first two waves of positive case increase were caused by the 20B clade, the third peak caused by the 20I (Alpha), while the last two peaks were caused by the 21J (Delta) and 21K (Omicron) variants. The functional analyses of mutations in sequences largely affected epitopes associated with protective HLA loci and did not cause the loss of the signal in PCR tests targeting ORF1ab and N genes as confirmed by RT-PCR. We also compared the performance of nanopore and Illumina short-read sequencing and showed the utility of nanopore sequencing as an efficient and affordable alternative for large-scale molecular epidemiology research. Thus, our paper describes new data on the genomic diversity of SARS-CoV-2 variants in Armenia in the global context of the virus molecular genomic surveillance.

Authors: Diana Avetyan, Siras Hakobyan, Maria Nikoghosyan, Lilit Ghukasyan, Gisane Khachatryan, Tamara Sirunyan, Nelli Muradyan, Roksana Zakharyan, Andranik Chavushyan, Varduhi Hayrapetyan, Anahit Hovhannisyan, Shah A. Mohamed Bakhash, Keith R. Jerome, Pavitra Roychoudhury, Alexander L. Greninger, Lyudmila Niazyan, Mher Davidyants, Gayane Melik-Andreasyan, Shushan Sargsyan, Lilit Nersisyan, Arsen Arakelyan

Date Published: 17th May 2022

Publication Type: Journal

Abstract (Expand)

Mutations in the BRCA1 and BRCA2 genes are known risk factors and drivers of breast and ovarian cancers. So far, few studies have been focused on understanding the differences in transcriptome andd functional landscapes associated with the disease (breast vs. ovarian cancers), gene (BRCA1 vs. BRCA2), and mutation type (germline vs. somatic). In this study, we were aimed at systemic evaluation of the association of BRCA1 and BRCA2 germline and somatic mutations with gene expression, disease clinical features, outcome, and treatment. We performed BRCA1/2 mutation centered RNA-seq data analysis of breast and ovarian cancers from the TCGA repository using transcriptome and phenotype “portrayal” with multi-layer self-organizing maps and functional annotation. The results revealed considerable differences in BRCA1- and BRCA2-dependent transcriptome landscapes in the studied cancers. Furthermore, our data indicated that somatic and germline mutations for both genes are characterized by deregulation of different biological functions and differential associations with phenotype characteristics and poly(ADP-ribose) polymerase (PARP)-inhibitor gene signatures. Overall, this study demonstrates considerable variation in transcriptomic landscapes of breast and ovarian cancers associated with the affected gene (BRCA1 vs. BRCA2), as well as the mutation type (somatic vs. germline). These results warrant further investigations with larger groups of mutation carriers aimed at refining the understanding of molecular mechanisms of breast and ovarian cancers.

Authors: Arsen Arakelyan, Ani Melkonyan, Siras Hakobyan, Uljana Boyarskih, Arman Simonyan, Lilit Nersisyan, Maria Nikoghosyan, Maxim Filipenko, Hans Binder

Date Published: 28th Jan 2021

Publication Type: Journal

Abstract (Expand)

Telomere length (TL) regulation is an important factor in ageing, reproduction and cancer development. Genetic, hereditary and environmental factors regulating TL are currently widely investigated, however, their relative contribution to TL variability is still understudied. We have used whole genome sequencing data of 250 family trios from the Genome of the Netherlands project to perform computational measurement of TL and a series of regression and genome-wide association analyses to reveal TL inheritance patterns and associated genetic factors. Our results confirm that TL is a largely heritable trait, primarily with mother's, and, to a lesser extent, with father's TL having the strongest influence on the offspring. In this cohort, mother's, but not father's age at conception was positively linked to offspring TL. Age-related TL attrition of 40 bp/year had relatively small influence on TL variability. Finally, we have identified TL-associated variations in ribonuclease reductase catalytic subunit M1 (RRM1 gene), which is known to regulate telomere maintenance in yeast. We also highlight the importance of multivariate approach and the limitations of existing tools for the analysis of TL as a polygenic heritable quantitative trait.

Authors: Lilit Nersisyan, Maria Nikoghosyan, Laurent C. Francioli, Androniki Menelaou, Sara L. Pulit, Clara C. Elbers, Wigard P. Kloosterman, Jessica van Setten, Isaäc J. Nijman, Ivo Renkens, Paul I. W. de Bakker, Freerk van Dijk, Pieter B. T. Neerincx, Patrick Deelen, Alexandros Kanterakis, Martijn Dijkstra, Heorhiy Byelas, K. Joeri van der Velde, Mathieu Platteel, Morris A. Swertz, Cisca Wijmenga, Pier Francesco Palamara, Itsik Pe’er, Kai Ye, Eric-Wubbo Lameijer, Matthijs H. Moed, Marian Beekman, Anton J. M. de Craen, H. Eka D. Suchiman, P. Eline Slagboom, Victor Guryev, Abdel Abdellaoui, Jouke Jan Hottenga, Mathijs Kattenberg, Gonneke Willemsen, Dorret I. Boomsma, Elisabeth M. van Leeuwen, Lennart C. Karssen, Najaf Amin, Fernando Rivadeneira, Aaron Isaacs, Albert Hofman, André G. Uitterlinden, Cornelia M. van Duijn, Mannis van Oven, Manfred Kayser, Martijn Vermaat, Jeroen F. J. Laros, Johan T. den Dunnen, David van Enckevort, Hailiang Mei, Mingkun Li, Mark Stoneking, Barbera D. C. van Schaik, Jan Bot, Tobias Marschall, Alexander Schönhuth, Jayne Y. Hehir-Kwa, Robert E. Handsaker, Paz Polak, Mashaal Sohail, Dana Vuzman, Karol Estrada, Steven A. McCarroll, Shamil R. Sunyaev, Fereydoun Hormozdiari, Vyacheslav Koval, Carolina Medina-Gomez, Ben Oostra, Jan H. Veldink, Leonard H. van den Berg, Steven J. Pitts, Shobha Potluri, Purnima Sundar, David R. Cox, Peter de Knijff, Qibin Li, Yingrui Li, Yuanping Du, Ruoyan Chen, Hongzhi Cao, Jun Wang, Ning Li, Sujie Cao, Jasper A. Bovenberg, Gert-Jan B. van Ommen, Arsen Arakelyan

Date Published: 10th Dec 2019

Publication Type: Journal

Abstract (Expand)

Background: Activation of telomere maintenance mechanisms (TMMs) is a hallmark of most cancers, and is required to prevent genome instability and to establish cellular immortality through reconstitution of capping of chromosome ends. TMM depends on the cancer type. Comparative studies linking tumor biology and TMM have potential impact for evaluating cancer onset and development. Methods: We have studied alterations of telomere length, their sequence composition and transcriptional regulation in mismatch repair deficient colorectal cancers arising in Lynch syndrome (LS-CRC) and microsatellite instable (MSI) sporadic CRC (MSI s-CRC), and for comparison, in microsatellite stable (MSS) s-CRC and in benign colon mucosa. Our study applied bioinformatics analysis of whole genome DNA and RNA sequencing data and a pathway model to study telomere length alterations and the potential effect of the "classical" telomerase (TEL-) and alternative (ALT-) TMM using transcriptomic signatures. Results: We have found progressive decrease of mean telomere length in all cancer subtypes compared with reference systems. Our results support the view that telomere attrition is an early event in tumorigenesis. TMM gets activated in all tumors studied due to concerted overexpression of a large fraction of genes with direct relation to telomere function, where only a very small fraction of them showed recurrent mutations. TEL-related transcriptional state was dominating in all CRC subtypes, showing, however, subtype-specific activation patterns; while contribution of the ALT-TMM was slightly more prominent in the hypermutated MSI s-CRC and LS-CRC. TEL-TMM is mainly activated by over-expression of DKC1 and/or TERT genes and their interaction partners, where DKC1 is more prominent in MSS than in MSI s-CRC and can serve as a transcriptomic marker of TMM activity. Conclusions: Our results suggest that transcriptional patterns are indicative for TMM pathway activation with subtle differences between TEL and ALT mechanisms in a CRC subtype-specific fashion. Sequencing data potentially provide a suited measure to study alterations of telomere length and of underlying transcriptional regulation. Further studies are needed to improve this method. Keywords: DNAseq and RNAseq data analysis; colorectal cancer; mismatch repair; pathway models; telomerase and alternative telomere maintenance; telomere attrition; telomere length; telomere repeat variants.

Authors: Lilit Nersisyan, Lydia Hopp, Henry Loeffler-Wirth, Jörg Galle, Markus Loeffler, Arsen Arakelyan, Hans Binder

Date Published: 5th Nov 2019

Publication Type: Journal

Powered by
(v.1.15.0-main)
Copyright © 2008 - 2024 The University of Manchester and HITS gGmbH