Generation of iPSCs from a Patient with the M694V Mutation in the MEFV Gene Associated with Familial Mediterranean Fever and Their Differentiation into Macrophages.

Familial Mediterranean fever (FMF) is a systemic autoinflammatory disorder caused by inherited mutations in the MEFV (Mediterranean FeVer) gene, located on chromosome 16 (16p13.3) and encoding the pyrin protein. Despite the existing data on MEFV mutations, the exact mechanism of their effect on the development of the pathological processes leading to the spontaneous and recurrent autoinflammatory attacks observed in FMF, remains unclear. Induced pluripotent stem cells (iPSCs) are considered an important tool to study the molecular genetic mechanisms of various diseases due to their ability to differentiate into any cell type, including macrophages, which contribute to the development of FMF. In this study, we developed iPSCs from an Armenian patient with FMF carrying the M694V, p.(Met694Val) (c.2080A>G, rs61752717) pathogenic mutation in exon 10 of the MEFV gene. As a result of direct differentiation, macrophages expressing CD14 and CD45 surface markers were obtained. We found that the morphology of macrophages derived from iPSCs of a patient with the MEFV mutation significantly differed from that of macrophages derived from iPSCs of a healthy donor carrying the wild-type MEFV gene.

Transposon DNA sequences facilitate the tissue-specific gene transfer of circulating tumor DNA between human cells.

The exchange of genes between cells is known to play an important physiological and pathological role in many organisms. We show that circulating tumor DNA (ctDNA) facilitates cell-specific gene transfer between human cancer cells and explain part of the mechanisms behind this phenomenon. As ctDNA migrates into the nucleus, genetic information is transferred. Cell targeting and ctDNA integration require ERVL, SINE or LINE DNA sequences. Chemically manufactured AluSp and MER11C sequences replicated multiple myeloma (MM) ctDNA cell targeting and integration. Additionally, we found that ctDNA may alter the treatment response of MM and pancreatic cancer models. This study shows that retrotransposon DNA sequences promote cancer gene transfer. However, because cell-free DNA has been detected in physiological and other pathological conditions, our findings have a broader impact than just cancer. Furthermore, the discovery that transposon DNA sequences mediate tissue-specific targeting will open up a new avenue for the delivery of genes and therapies.

Pan-cancer analysis of telomere maintenance mechanisms.

Telomeres, protective caps at chromosome ends, maintain genomic stability and control cell lifespan. Dysregulated telomere maintenance mechanisms (TMMs) are cancer hallmarks, enabling unchecked cell proliferation. We conducted a pan-cancer evaluation of TMM using RNA sequencing data from The Cancer Genome Atlas for 33 different cancer types and analyzed the activities of telomerase-dependent (TEL) and alternative lengthening of telomeres (ALT) TMM pathways in detail. To further characterize the TMM profiles, we categorized the tumors based on their ALT and TEL TMM pathway activities into five major phenotypes: ALT high TEL low, ALT low TEL low, ALT middle TEL middle, ALT high TEL high, and ALT low TEL high. These phenotypes refer to variations in telomere maintenance strategies, shedding light on the heterogeneous nature of telomere regulation in cancer. Moreover, we investigated the clinical implications of TMM phenotypes by examining their associations with clinical characteristics and patient outcomes. Specific TMM profiles were linked to specific survival patterns, emphasizing the potential of TMM profiling as a prognostic indicator and aiding in personalized cancer treatment strategies. Gene ontology analysis of the TMM phenotypes unveiled enriched biological processes associated with cell cycle regulation (both TEL and ALT), DNA replication (TEL), and chromosome dynamics (ALT) showing that telomere maintenance is tightly intertwined with cellular processes governing proliferation and genomic stability. Overall, our study provides an overview of the complexity of transcriptional regulation of telomere maintenance mechanisms in cancer.

Temporal changes of gene expression in health, schizophrenia, bipolar disorder, and major depressive disorder.

The molecular events underlying the development, manifestation, and course of schizophrenia, bipolar disorder, and major depressive disorder span from embryonic life to advanced age. However, little is known about the early dynamics of gene expression in these disorders due to their relatively late manifestation. To address this, we conducted a secondary analysis of post-mortem prefrontal cortex datasets using bioinformatics and machine learning techniques to identify differentially expressed gene modules associated with aging and the diseases, determine their time-perturbation points, and assess enrichment with expression quantitative trait loci (eQTL) genes. Our findings revealed early, mid, and late deregulation of expression of functional gene modules involved in neurodevelopment, plasticity, homeostasis, and immune response. This supports the hypothesis that multiple hits throughout life contribute to disease manifestation rather than a single early-life event. Moreover, the time-perturbed functional gene modules were associated with genetic loci affecting gene expression, highlighting the role of genetic factors in gene expression dynamics and the development of disease phenotypes. Our findings emphasize the importance of investigating time-dependent perturbations in gene expression before the age of onset in elucidating the molecular mechanisms of psychiatric disorders.

Development of an optimized, non-stem cell line for intranasal delivery of therapeutic cargo to the central nervous system.

Neural stem cells (NSCs) are considered to be valuable candidates for delivering a variety of anti-cancer agents, including oncolytic viruses, to brain tumors. However, owing to the previously reported tumorigenic potential of NSC cell lines after intranasal administration (INA), here we identified the human hepatic stellate cell line LX-2 as a cell type capable of longer resistance to replication of oncolytic adenoviruses (OAVs) as a therapeutic cargo, and that is non-tumorigenic after INA. Our data show that LX-2 cells can longer withstand the OAV XVir-N-31 replication and oncolysis than NSCs. By selecting the highly migratory cell population out of LX-2, an offspring cell line with a higher and more stable capability to migrate was generated. Additionally, as a safety backup, we applied genomic herpes simplex virus thymidine kinase (HSV-TK) integration into LX-2, leading to high vulnerability to ganciclovir (GCV). Histopathological analyses confirmed the absence of neoplasia in the respiratory tracts and brains of immuno-compromised mice 3 months after INA of LX-2 cells. Our data suggest that LX-2 is a novel, robust, and safe cell line for delivering anti-cancer and other therapeutic agents to the brain.

Wild grapes of Armenia: unexplored source of genetic diversity and disease resistance.

The present study is the first in-depth research evaluating the genetic diversity and potential resistance of Armenian wild grapes utilizing DNA-based markers to understand the genetic signature of this unexplored germplasm. In the proposed research, five geographical regions with known viticultural history were explored. A total of 148 unique wild genotypes were collected and included in the study with 48 wild individuals previously collected as seed. A total of 24 nSSR markers were utilized to establish a fingerprint database to infer information on the population genetic diversity and structure. Three nSSR markers linked to the Ren1 locus were analyzed to identify potential resistance against powdery mildew. According to molecular fingerprinting data, the Armenian V. sylvestris gene pool conserves a high genetic diversity, displaying 292 different alleles with 12.167 allele per loci. The clustering analyses and diversity parameters supported eight genetic groups with 5.6% admixed proportion. The study of genetic polymorphism at the Ren1 locus revealed that 28 wild genotypes carried three R-alleles and 34 wild genotypes carried two R-alleles associated with PM resistance among analyzed 107 wild individuals. This gene pool richness represents an immense reservoir of under-explored genetic diversity and breeding potential. Therefore, continued survey and research efforts are crucial for the conservation, sustainable management, and utilization of Armenian wild grape resources in the face of emerging challenges in viticulture.

PSF toolkit: an R package for pathway curation and topology-aware analysis.

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.

Transcriptomic Maps of Colorectal Liver Metastasis: Machine Learning of Gene Activation Patterns and Epigenetic Trajectories in Support of Precision Medicine.

The molecular mechanisms of the liver metastasis of colorectal cancer (CRLM) remain poorly understood. Here, we applied machine learning and bioinformatics trajectory inference to analyze a gene expression dataset of CRLM. We studied the co-regulation patterns at the gene level, the potential paths of tumor development, their functional context, and their prognostic relevance. Our analysis confirmed the subtyping of five liver metastasis subtypes (LMS). We provide gene-marker signatures for each LMS, and a comprehensive functional characterization that considers both the hallmarks of cancer and the tumor microenvironment. The ordering of CRLMs along a pseudotime-tree revealed a continuous shift in expression programs, suggesting a developmental relationship between the subtypes. Notably, trajectory inference and personalized analysis discovered a range of epigenetic states that shape and guide metastasis progression. By constructing prognostic maps that divided the expression landscape into regions associated with favorable and unfavorable prognoses, we derived a prognostic expression score. This was associated with critical processes such as epithelial-mesenchymal transition, treatment resistance, and immune evasion. These factors were associated with responses to neoadjuvant treatment and the formation of an immuno-suppressive, mesenchymal state. Our machine learning-based molecular profiling provides an in-depth characterization of CRLM heterogeneity with possible implications for treatment and personalized diagnostics.

Long-term environmental metal exposure is associated with hypomethylation of CpG sites in NFKB1 and other genes related to oncogenesis.

BACKGROUND: Long-term environmental exposure to metals leads to epigenetic changes and may increase risks to human health. The relationship between the type and level of metal exposure and epigenetic changes in subjects exposed to high concentrations of metals in the environment is not yet clear. The aim of our study is to find the possible association of environmental long-term exposure to metals with DNA methylation changes of genes related to immune response and carcinogenesis. We investigated the association of plasma levels of 21 essential and non-essential metals detected by ICP-MS and the methylation level of 654 CpG sites located on NFKB1, CDKN2A, ESR1, APOA5, IGF2 and H19 genes assessed by targeted bisulfite sequencing in a cohort of 40 subjects living near metal mining area and 40 unexposed subjects. Linear regression was conducted to find differentially methylated positions with adjustment for gender, age, BMI class, smoking and metal concentration. RESULTS: In the metal-exposed group, five CpGs in the NFKB1 promoter region were hypomethylated compared to unexposed group. Four differentially methylated positions (DMPs) were associated with multiple metals, two of them are located on NFKB1 gene, and one each on CDKN2A gene and ESR1 gene. Two DMPs located on NFKB1 (chr4:102500951, associated with Be) and IGF2 (chr11:2134198, associated with U) are associated with specific metal levels. The methylation status of the seven CpGs located on NFKB1 (3), ESR1 (2) and CDKN2A (2) positively correlated with plasma levels of seven metals (As, Sb, Zn, Ni, U, I and Mn). CONCLUSIONS: Our study revealed methylation changes in NFKB1, CDKN2A, IGF2 and ESR1 genes in individuals with long-term human exposure to metals. Further studies are needed to clarify the effect of environmental metal exposure on epigenetic mechanisms and pathways involved.

Generation of three induced pluripotent stem cell lines (RAUi001-A, RAUi001-B and RAUi001-C) from peripheral blood mononuclear cells of a healthy Armenian individual.

The study of pathological processes in cells carrying mutations should be carried out in comparison with a healthy control group. Familial Mediterranean fever (FMF), which is caused by a mutation in the MEFV gene, is predominantly found in people of Armenian nationality with the prevalence of 14-100 per 10000. We have obtained induced pluripotent stem cells (iPSCs) from Armenian healthy patient, which will be included as a control group in the study of this disease. iPSCs rapidly proliferate in colonies of cells with a typical pluripotent-like morphology, have a normal karyotype (46,XX). iPSCs express pluripotency markers (OCT4, SOX2, TRA-1-60, NANOG) and are able to give derivatives of three germ layers.

Dual domestications and origin of traits in grapevine evolution.

We elucidate grapevine evolution and domestication histories with 3525 cultivated and wild accessions worldwide. In the Pleistocene, harsh climate drove the separation of wild grape ecotypes caused by continuous habitat fragmentation. Then, domestication occurred concurrently about 11,000 years ago in Western Asia and the Caucasus to yield table and wine grapevines. The Western Asia domesticates dispersed into Europe with early farmers, introgressed with ancient wild western ecotypes, and subsequently diversified along human migration trails into muscat and unique western wine grape ancestries by the late Neolithic. Analyses of domestication traits also reveal new insights into selection for berry palatability, hermaphroditism, muscat flavor, and berry skin color. These data demonstrate the role of the grapevines in the early inception of agriculture across Eurasia.

Transcriptional states of CAR-T infusion relate to neurotoxicity - lessons from high-resolution single-cell SOM expression portraying.

Anti-CD19 CAR-T cell immunotherapy is a hopeful treatment option for patients with B cell lymphomas, however it copes with partly severe adverse effects like neurotoxicity. Single-cell resolved molecular data sets in combination with clinical parametrization allow for comprehensive characterization of cellular subpopulations, their transcriptomic states, and their relation to the adverse effects. We here present a re-analysis of single-cell RNA sequencing data of 24 patients comprising more than 130,000 cells with focus on cellular states and their association to immune cell related neurotoxicity. For this, we developed a single-cell data portraying workflow to disentangle the transcriptional state space with single-cell resolution and its analysis in terms of modularly-composed cellular programs. We demonstrated capabilities of single-cell data portraying to disentangle transcriptional states using intuitive visualization, functional mining, molecular cell stratification, and variability analyses. Our analysis revealed that the T cell composition of the patient's infusion product as well as the spectrum of their transcriptional states of cells derived from patients with low ICANS grade do not markedly differ from those of cells from high ICANS patients, while the relative abundancies, particularly that of cycling cells, of LAG3-mediated exhaustion and of CAR positive cells, vary. Our study provides molecular details of the transcriptomic landscape with possible impact to overcome neurotoxicity.

Molecular Epidemiology and Virulence of Non-Typhoidal Salmonella in Armenia.

In this work, we analysed human isolates of nontyphoidal Salmonella enterica subsp. enterica (NTS), which were collected from salmonellosis cases in Armenia from 1996 to 2019. This disease became a leading food-borne bacterial infection in the region, with the younger age groups especially affected. The isolates were characterised by serotyping, Enterobacterial Repetitive Intergenic Consensus (ERIC-PCR) typing, and whole genome sequencing (WGS). The main serotypes were S. Typhimurium, S. Enteritidis, and S. Arizonae. ERIC-PCR indicated a high degree of clonality among S. Typhimurium strains, which were also multidrug-resistant and produced extended spectrum beta-lactamases. During the study period, the frequency of S. Typhimurium and S. Arizonae isolations decreased, but with the increase in S. Enteritidis and other NTS. A total of 42 NTS isolates were subjected to WGS and explored for virulence-related traits and the corresponding genetic elements. Some virulence and genetic factors were shared by all NTS serotypes, while the main differences were attributed to the serotype-specific diversity of virulence genes, SPIs, virulence plasmids, and phages. The results indicated the variability and dynamics in the epidemiology of salmonellosis and a high virulence potential of human NTS isolates circulating in the region.

Spaceflight-Associated Changes of snoRNAs in Peripheral Blood Mononuclear Cells and Plasma Exosomes-A Pilot Study.

During spaceflight, astronauts are exposed to various physiological and psychological stressors that have been associated with adverse health effects. Therefore, there is an unmet need to develop novel diagnostic tools to predict early alterations in astronauts' health. Small nucleolar RNA (snoRNA) is a type of short non-coding RNA (60-300 nucleotides) known to guide 2'-O-methylation (Nm) or pseudouridine (ψ) of ribosomal RNA (rRNA), small nuclear RNA (snRNA), or messenger RNA (mRNA). Emerging evidence suggests that dysregulated snoRNAs may be key players in regulating fundamental cellular mechanisms and in the pathogenesis of cancer, heart, and neurological disease. Therefore, we sought to determine whether the spaceflight-induced snoRNA changes in astronaut's peripheral blood (PB) plasma extracellular vesicles (PB-EV) and peripheral blood mononuclear cells (PBMCs). Using unbiased small RNA sequencing (sRNAseq), we evaluated changes in PB-EV snoRNA content isolated from astronauts (n = 5/group) who underwent median 12-day long Shuttle missions between 1998 and 2001. Using stringent cutoff (fold change > 2 or log2-fold change >1, FDR < 0.05), we detected 21 down-and 9-up-regulated snoRNAs in PB-EVs 3 days after return (R + 3) compared to 10 days before launch (L-10). qPCR validation revealed that SNORA74A was significantly down-regulated at R + 3 compared to L-10. We next determined snoRNA expression levels in astronauts' PBMCs at R + 3 and L-10 (n = 6/group). qPCR analysis further confirmed a significant increase in SNORA19 and SNORA47 in astronauts' PBMCs at R + 3 compared to L-10. Notably, many downregulated snoRNA-guided rRNA modifications, including four Nms and five ψs. Our findings revealed that spaceflight induced changes in PB-EV and PBMCs snoRNA expression, thus suggesting snoRNAs may serve as potential novel biomarkers for monitoring astronauts' health.

Astronauts Plasma-Derived Exosomes Induced Aberrant EZH2-Mediated H3K27me3 Epigenetic Regulation of the Vitamin D Receptor.

There are unique stressors in the spaceflight environment. Exposure to such stressors may be associated with adverse effects on astronauts' health, including increased cancer and cardiovascular disease risks. Small extracellular vesicles (sEVs, i.e., exosomes) play a vital role in intercellular communication and regulate various biological processes contributing to their role in disease pathogenesis. To assess whether spaceflight alters sEVs transcriptome profile, sEVs were isolated from the blood plasma of 3 astronauts at two different time points: 10 days before launch (L-10) and 3 days after return (R+3) from the Shuttle mission. AC16 cells (human cardiomyocyte cell line) were treated with L-10 and R+3 astronauts-derived exosomes for 24 h. Total RNA was isolated and analyzed for gene expression profiling using Affymetrix microarrays. Enrichment analysis was performed using Enrichr. Transcription factor (TF) enrichment analysis using the ENCODE/ChEA Consensus TF database identified gene sets related to the polycomb repressive complex 2 (PRC2) and Vitamin D receptor (VDR) in AC16 cells treated with R+3 compared to cells treated with L-10 astronauts-derived exosomes. Further analysis of the histone modifications using datasets from the Roadmap Epigenomics Project confirmed enrichment in gene sets related to the H3K27me3 repressive mark. Interestingly, analysis of previously published H3K27me3-chromatin immunoprecipitation sequencing (ChIP-Seq) ENCODE datasets showed enrichment of H3K27me3 in the VDR promoter. Collectively, our results suggest that astronaut-derived sEVs may epigenetically repress the expression of the VDR in human adult cardiomyocytes by promoting the activation of the PRC2 complex and H3K27me3 levels.

Integrated Multi-Omics Maps of Lower-Grade Gliomas.

Multi-omics high-throughput technologies produce data sets which are not restricted to only one but consist of multiple omics modalities, often as patient-matched tumour specimens. The integrative analysis of these omics modalities is essential to obtain a holistic view on the otherwise fragmented information hidden in this data. We present an intuitive method enabling the combined analysis of multi-omics data based on self-organizing maps machine learning. It "portrays" the expression, methylation and copy number variations (CNV) landscapes of each tumour using the same gene-centred coordinate system. It enables the visual evaluation and direct comparison of the different omics layers on a personalized basis. We applied this combined molecular portrayal to lower grade gliomas, a heterogeneous brain tumour entity. It classifies into a series of molecular subtypes defined by genetic key lesions, which associate with large-scale effects on DNA methylation and gene expression, and in final consequence, drive with cell fate decisions towards oligodendroglioma-, astrocytoma- and glioblastoma-like cancer cell lineages with different prognoses. Consensus modes of concerted changes of expression, methylation and CNV are governed by the degree of co-regulation within and between the omics layers. The method is not restricted to the triple-omics data used here. The similarity landscapes reflect partly independent effects of genetic lesions and DNA methylation with consequences for cancer hallmark characteristics such as proliferation, inflammation and blocked differentiation in a subtype specific fashion. It can be extended to integrate other omics features such as genetic mutation, protein expression data as well as extracting prognostic markers.

The Effect of Low-Energy Laser-Driven Ultrashort Pulsed Electron Beam Irradiation on Erythropoiesis and Oxidative Stress in Rats.

Anemia is a commonly observed consequence of whole-body exposure to a dose of X-ray or gamma irradiation of the order of the mean lethal dose in mammals, and it is an important factor for the determination of the survival of animals. The aim of this study was to unravel the effect of laser-driven ultrashort pulsed electron beam (UPEB) irradiation on the process of erythropoiesis and the redox state in the organism. Wistar rats were exposed to laser-driven UPEB irradiation, after which the level of oxidative stress and the activities of different antioxidant enzymes, as well as blood smears, bone marrow imprints and sections, erythroblastic islets, hemoglobin and hematocrit, hepatic iron, DNA, and erythropoietin levels, were assessed on the 1st, 3rd, 7th, 14th, and 28th days after irradiation. Despite the fact that laser-driven UPEB irradiation requires quite low doses and repetition rates to achieve the LD50 in rats, our findings suggest that whole-body exposure with this new type of irradiation causes relatively mild anemia in rats, with subsequent fast recovery up to the 28th day. Moreover, this novel type of irradiation causes highly intense processes of oxidative stress, which, despite being relatively extinguished, did not reach the physiologically stable level even at the 28th day after irradiation due to the violations in the antioxidant system of the organism.

Molecular Analysis of SARS-CoV-2 Lineages in Armenia.

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.