Publications

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

Abstract (Expand)

Spatial transcriptomics (ST) has transformed genomics by mapping gene expression onto intact tissue architecture, uncovering intricate cellular interactions that bulk and single-cell RNA sequencing often overlook. Traditional ST workflows typically involve clustering spots, performing differential expression analyses, and annotating results via gene-set methods such as overrepresentation analysis (ORA) or gene set enrichment analysis (GSEA). More recent spatially-aware techniques extend these approaches by incorporating tissue organization into gene-set scoring. However, because they operate primarily at the level of individual genes, they may overlook the connectivity and topology of biological pathways, limiting their capacity to trace the propagation of signaling events within tissue regions. In this study, we address that gap by translating gene expression into pathway-level activity using the Pathway Signal Flow (PSF) algorithm. PSF integrates expression data with curated interaction networks to compute numeric activity scores for each branch of a biological pathway, producing a functionally annotated feature space that captures downstream signaling effects as branch-specific activity values. We applied PSF to two public 10x Genomics Visium datasets (human melanoma and mouse brain) and compared clustering based on PSF-derived pathway activities from 40 curated Kyoto Encyclopedia of Genes and Genomes (KEGG) signaling pathways and gene expression with standard Seurat Louvain clustering and spatially aware methods (Vesalius, spatialGE). We observed good correspondence between PSF-based and expression-based clustering when spatially aware clustering methods were used. This suggests that branch-level pathway activities can themselves drive clustering and pinpoint spatially deregulated processes. To assess cluster-specific functional annotation, we compared PSF results to conventional ORA (based on marker genes) and GSDensity (based on cluster-specific gene sets). PSF identified a broader set of significant pathways with substantial overlap with both ORA and GSDensity, providing increased sensitivity due to its branch-level resolution. We further demonstrated that PSF-derived activity values can be used to detect spatially deregulated pathway branches, yielding results comparable to those obtained with spatially aware gene set analysis approaches such as GSDensity and spatialGE. The availability of pathway topology and branch-specific information also enabled the identification of potential intercellular communication via ligand-receptor interactions between deregulated pathways in adjacent tumor regions. To support interactive exploration of results, we developed the PSF Spatial Browser, an R Shiny application for visualizing pathway activities, gene expression patterns, and deregulated pathway networks.

Authors: Siras Hakobyan, Maria Schmidt, H. Binder, A. Arakelyan

Date Published: 14th Aug 2025

Publication Type: Journal

Abstract (Expand)

Deep space represents a challenging environment for human exploration and can be accompanied by harmful health-related risks. We aimed to assess the effect of simplified galactic cosmic ray simulated (simGCRsim) and gamma (γ) ionizing radiation (IR) on transcriptome changes in right ventricular (RV) tissue after a single low dose (0.5 Gy, 500 MeV/nucleon) full body exposure in C57BL/6J male and female mice. In females, no differentially expressed genes (DEGs) and only 2 upregulated genes in males exposed to γ-IR were revealed. In contrast, exposure to simGCRsim-IR resulted in 4 DEGs in females and 371 DEGs in males, suggesting longer-lasting and sex-biased DEGs after simGCRsim-IR. Overrepresentation analysis of DEGs in simGCRsim-IR males revealed significant enrichment in pathways related to muscle contraction, hypertrophic cardiomyopathy, oxytocin release, the regulation of cytoskeleton, and genes associated with Alzheimer’s, Huntington’s, and Parkinson’s diseases. Our results suggested the RV transcriptome exhibits distinct responses after exposure based on both the IR and sex.

Authors: Roksana Zakharyan, Siras Hakobyan, Agnieszka Brojakowska, Malik Bisserier, Shihong Zhang, Mary K. Khlgatian, Amit Kumar Rai, Suren Davitavyan, Ani Stepanyan, Tamara Sirunyan, Gisane Khachatryan, Susmita Sahoo, Venkata Naga Srikanth Garikipati, Arsen Arakelyan, David A. Goukassian

Date Published: 21st Jul 2025

Publication Type: Journal

Abstract (Expand)

Telomere maintenance mechanisms (TMMs) play a critical role in cancer biology, particularly in lower-grade gliomas (LGGs), where telomere dynamics and pathway activity remain poorly understood. In this study, we analyzed TCGA-LGG and CGGA datasets, focusing on telomere length variations, pathway activity, and survival data across IDH subtypes. Additional validation was performed using the GEO COPD and GBM datasets, ensuring consistency in data processing and batch effect correction. Our analysis revealed significant differences in TEL pathway activation between Short- and Long-TL groups, emphasizing the central role of TERT in telomere maintenance. In contrast, ALT pathway activation displayed subtype-specific patterns, with IDH-wt tumors exhibiting the highest ALT activity, primarily driven by the RAD51 branch. Validation using CGGA data confirmed these findings, demonstrating consistent TEL and ALT pathway behaviors across datasets. Additionally, genetic subtype analysis revealed substantial telomere length variability associated with ATRX and IDH mutation status. Notably, IDHwt-ATRX WT tumors exhibited the shortest telomere length and the highest ALT pathway activity. These findings highlight distinct telomere regulatory dynamics across genetic subtypes of LGG and provide new insights into potential therapeutic strategies targeting telomere maintenance pathways.

Authors: Meline Hakobyan, Hans Binder, Arsen Arakelyan

Date Published: 28th Apr 2025

Publication Type: Journal

Abstract (Expand)

Space irradiation (IR) is an important health risk for deep-space missions. We reported heart failure with preserved ejection fraction like cardiac phenotype 660-days following exposure to a single-dose of a simplified galactic cosmic ray simulation (simGCRsim) only in males with functional and structural impairment in left ventricular (LV) function. This sex-based dichotomy prompted us to investigate sex-specific changes in the LV transcriptome in three-month-old male and female mice exposed to 137Cs-γ- or simGCRsim-IR. Non-IR male and female (10 each) mice served as controls. LVs were collected at 440/660- and 440/550-days post-IR, male and female, respectively. RNA sequencing, differential gene expression, and functional annotation were performed on tissues from 5 mice/group. Sex and post-IR time points had the greatest influence on gene expression, surpassing the IR-type effects. SimGCRsim-IR showed more persistent transcriptome changes than γ-IR. We suggest that the single IR effects can persist up to 550-660 days, with overwhelmingly sex-biased responses at individual gene expression level.

Authors: Roksana Zakharyan, Siras Hakobyan, Agnieszka Brojakowska, Suren Davitavyan, Ani Stepanyan, Tamara Sirunyan, Gisane Khachatryan, Mary K. Khlgatian, Malik Bisserier, Shihong Zhang, Susmita Sahoo, Lahouaria Hadri, Venkata Naga Srikanth Garikipati, Arsen Arakelyan, David A. Goukassian

Date Published: 18th Feb 2025

Publication Type: Journal

Abstract (Expand)

The global dissemination of Klebsiella pneumoniae pathotypes with multidrug-resistant (MDR) and hypervirulent traits poses a threat to public health. The situation in Armenia is unclear, and we performed a comprehensive characterisation of 48 clinical isolates of K. pneumoniae, collected from 2018 to 2024. The majority of the isolates (64.58%) were extensively drug-resistant (XDR) and MDR. Genomic analysis of 21 isolates revealed the presence of international high-risk MDR clones (ST395, ST15, and ST307). The ST395 strains were isolated from children and resisted the first-line drugs such as beta-lactams. These isolates harboured a range of virulence determinants, from capsule polysaccharides to siderophores to regulators of the mucoid phenotype. The ST395 strains are enriched by ICEs, plasmids, and prophages, on which antimicrobial resistance (AMR) and virulence genes are located and which may lead to the convergence of MDR and hypervirulent traits. There is a widespread non-specific AMR mechanism among our K. pneumoniae strains. These are mutations in the porin genes, which reduce permeability to antimicrobials, and mutations in the regulators of efflux pumps, which lead to overexpression of drug efflux pumps such as AcrAB. These mechanisms may contribute to the elevated MICs and confer AMR to strains with no specific AMR genes.

Authors: Anahit Sedrakyan, Zaruhi Gevorgyan, Magdalina Zakharyan, Karine Arakelova, Shoghik Hakobyan, Alvard Hovhannisyan, Rustam Aminov

Date Published: 9th Jan 2025

Publication Type: Journal

Abstract (Expand)

Mediterranean fever (FMF) is a genetically determined autoinflammatory disease transmitted mostly by an autosomal recessive mechanism and caused by point mutations of the MEFV (Mediterranean FeVer) gene. The aim of this study was to evaluate the expression of inflammasome genes (p65, Casp1, MEFV, and NLRP3) in patients with FMF compared to controls to understand the changes playing a key role in disease development. We found altered expression levels of the full-length MEFV isoform as well as Casp1 and p65 in FMF patients versus controls. This, once again, highlighted the significance of inflammasome genes in terms of FMF.

Authors: Varduhi Hayrapetyan, Lana Karapetyan, Lilit Ghukasyan, Sofi Atshemyan, Hovsep Ghazaryan, Valentina Vardanyan, Vahan Mukuchyan, Arsen Arakelyan, Roksana Zakharyan

Date Published: 2nd Dec 2024

Publication Type: Journal

Abstract (Expand)

Mediterranean Fever (FMF) is a genetic disorder with complex inheritance patterns and genotype-phenotype associations, and it is highly prevalent in Armenia. FMF typically follows an autosomal recessive inheritance pattern (OMIM: 249100), though it can occasionally display a rare dominant inheritance pattern with variable penetrance (OMIM։134610). The disease is caused by mutations in the MEFV gene, which encodes the pyrin protein. While the 26 most prevalent mutations account for nearly 99% of all FMF cases, more than 60 pathogenic mutations have been identified. In this study, we aimed to develop an affordable nanopore sequencing method for full-length MEFV gene mutation detection to aid in the diagnosis and screening of FMF. We employed a multiplex amplicon sequencing approach, allowing for the processing of up to 12 samples on both Flow cells and Flongle flow cells. The results demonstrated near-complete concordance between nanopore variant calling and qPCR genotypes. Moreover, nanopore sequencing identified additional variants, which were confirmed by whole exome sequencing. Additionally, intronic and UTR variants were detected. Our findings demonstrate the feasibility of full-gene nanopore sequencing for detecting FMF-associated pathogenic variants. The method is cost-effective, with costs comparable to those of the qPCR test, making it particularly suitable for settings with limited laboratory infrastructure. Further clinical validation using larger sample cohorts will be necessary.

Authors: Lilit Ghukasyan, Gisane Khachatryan, Tamara Sirunyan, Arpine Minasyan, Siras Hakobyan, Andranik Chavushyan, Varduhi Hayrapetyan, Hovsep Ghazaryan, Gevorg Martirosyan, Gohar Mkrtchyan, Valentina Vardanyan, Vahan Mukuchyan, Ashot Davidyants, Roksana Zakharyan, Arsen Arakelyan

Date Published: 29th Nov 2024

Publication Type: Journal

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