1. Comprehensive molecular characterization of craniopharyngiomas using whole transcriptome and spatial transcriptomics approachesŠpela Kert, Alenka Matjašič, Jože Pižem, Jernej Mlakar, Matic Bošnjak, Miha Jerala, Primož Kotnik, Barbara Faganel Kotnik, Lidija Kitanovski, Andrej Zupan, 2025, original scientific article Abstract: Craniopharyngiomas (CPs) are rare benign brain tumors that are classified as WHO grade I, with two subtypes: adamantinomatous craniopharyngioma (ACP) and papillary craniopharyngioma (PCP). ACP is caused by somatic mutations in exon 3 of the CTNNB1 gene activating the Wnt signaling pathway. PCP is associated with somatic BRAF p.V600E mutations activating the MAPK signaling pathway. Understanding their molecular differences is crucial for diagnosis and treatment. This study aimed to analyze common somatic alterations in ACP and PCP using bulk transcriptome sequencing and in situ spatial transcriptomics. RNA sequencing and high-resolution spatial profiling were used to detect mutations and examine gene expression differences among ACP, PCP, and healthy pituitary tissue. Whole transcriptome sequencing was performed on 24 tumor samples, with healthy pituitary data from the GTEx portal. Bioinformatics analysis utilized the CTAT mutation pipeline, with Sanger sequencing for validation. Results confirmed BRAF p.V600E mutations in all PCP samples and CTNNB1 mutations in all ACP samples. Differential gene expression analysis highlighted distinct molecular profiles and reinforced the involvement of Wnt and MAPK signaling. Spatial profiling identified 41 differentially expressed genes between ACP and PCP. This study provides critical insights into CP biology, supporting improved diagnostics and potential therapeutic strategies. Keywords: craniopharyngioma, differential gene expression, in situ spatial profiling, somatic mutation detection, transcriptional analysis Published in DiRROS: 10.04.2026; Views: 58; Downloads: 34
Full text (2,18 MB) This document has many files! More... |
2. Transcriptomic responses of oxidative and genotoxic stress responsive genes after exposure to MxFe3-xO4 (M = Fe, Zn, Mn) in an advanced 3D human hepatic in vitro model : version v1Iza Rozman, Alja Štern, Bojana Žegura, Gerardo F. Goya, Álvaro Gallo-Cordova, María del Puerto Morales, 2025, complete scientific database of research data Abstract: Nanosized spinel-type ferrites have gained recognition as a unique class of engineered nanomaterials with promising applications, but their safety profiles remain insufficiently explored. Although iron (Fe), zinc (Zn), and manganese (Mn) are biologically relevant elements, the use of Zn- and Mn-containing ferrite nanoparticles in biomedical contexts demands careful (geno)toxicity evaluation. In this study, three ferrite nanoparticles – γFe2O3 (FeNPs), Zn0.7Fe2.3O4 (ZnNPs), and Mn0.4Fe2.6O4 (MnNPs) – synthesised through a microwave-assisted polyol route, functionalized with citric acid to improve colloidal stability, were evaluated for their potential (geno)toxic effects in an advanced in vitro 3D cell model, HepG2 spheroids. Cellular stress responses upon exposure to the particle were assessed using toxicogenomic analysis.This approach allows the identification of early molecular events that may precede overt toxicity, supporting a mechanistic understanding of adverse outcomes and facilitating the development of predictive biomarkers for hazard assessment. In the present study, the expression of selected DNA damage-responsive genes (TP53, MDM2, GADD45a, CDKN1A, OGG1, and JUNB), apoptosis-related genes (BCL2 and BAX) and oxidative stress response genes (SOD1, CAT, GPX1, GCLC, and GSR) was evaluated. The expression of the selected genes after exposure to the tested nanoparticles was analysed by qPCR primer assays (Applied Biosystems, USA) and One 48.48 Dynamic Array IFC for Gene Expression (Fluidigm, USA). After 24 and 96 hours of exposure, the spheroids were collected, and total RNA was isolated using the RNeasy Mini Kit from Qiagen (Qiagen, Germany) according to the manufacturer's instructions. 10 µg/mL etoposide served as athe positive control for the toxicogenomic analysis. RNA concentration and purity were assessed using a NanoDrop 1000 spectrophotometer (Thermo Fisher Scientific) by measuring absorbance at 260/280 nm and gele efectrophoresis (Figure 1). Reverse transcription of 1 µg total RNA per sample was performed with the High-Capacity cDNA Reverse Transcription Kit (Applied Biosystems, MA, USA) on a BIO-RAD T100 thermal cycler under conditions listed in Table 3. For preamplification, 4 µL of each of the 24 selected TaqMan assays (SM2) were pooled into a primer mix. The reaction mixture was prepared using TATAA PreAmp GrandMasterMix (Tataa Biocenter, Sweden), the primer pool, and nuclease-free water, following manufacturer instructions. Negative controls (NTC for preamplification and NTCq for qPCR) were included. Each reaction contained 8 µL of mix and 2 µL of 5× diluted cDNA, processed in a 96-deep well plate, sealed, vortexed, and centrifuged (1000 g, 1 min). Preamplification was carried out on a BIO-RAD T100 thermal cycler under conditions in Table 4. Gene expression analysis used TaqMan Universal PCR Master Mix and TaqMan Gene Expression Assays (Table 6). Preamplified samples were diluted 10× with nuclease-free water. Assays were prepared by mixing equal volumes (6 µL) of each assay with Fluidigm Assay Loading Reagent Kit – 10IFCS. The reaction premix combined DNA Sample Loading Reagent and Fast Probe Master Mix (Biotium/Roche) and was added to each diluted cDNA sample. qPCR was performed on 48.48 Dynamic Array™ IFC chips using the Fluidigm BioMark™ HD System under conditions in Table 5. Data were analysed with Fluidigm Gene Expression Analysis Software and quantGenious. Fold changes >1.5 or <0.66 were considered biologically relevant. Statistical significance between NP-exposed cells and solvent controls was assessed using ANOVA and Dunnett’s test in GraphPad Prism v9 (GraphPad Software, CA, USA). Keywords: ferrite-based nanoparticles, HepG2 spheroids, toxicogenomics, changes in gene expression Published in DiRROS: 24.03.2026; Views: 147; Downloads: 142
Research data (172,98 KB) This document has many files! More... |
3. Transcriptome profiles of skeletal muscle in krškopolje and modern hybrid pigs : genotype comparison and effects of dietary protein reductionMilka Vrecl, Gregor Fazarinc, Malan Štrbenc, Klavdija Poklukar Žnidaršič, Marjeta Čandek-Potokar, Martin Škrlep, 2026, complete scientific database of research data Abstract: RNA sequencing was performed to compare two pig breeds (the Slovenian Krškopolje breed and a modern commercial breed) and to evaluate the effect of breed specific dietary protein reduction on the transcriptomic profiles of two skeletal muscles: longissimus thoracis et lumborum (LD) and semispinalis capitis (SSC). Differentially expressed genes (DEGs) were identified using the criteria |log2FC| > 1 and q value < 0.005. The effect of dietary protein reduction was minimal in both breeds and muscles, with only small numbers of DEGs detected. In Krškopolje pigs receiving medium protein (MP) or low protein (LP) diets, no DEGs were identified. In contrast, in the modern breed, 10 DEGs were detected when comparing the high protein (HP) with the MP diet group. A similar pattern was observed in SSC, where 19 DEGs were detected in Krškopolje pigs (MP vs. LP) and 16 DEGs in the modern breed (HP vs. LP). In comparison, when the LD and SSC transcriptome profiles of the modern breed were compared with those of the Krškopolje breed, the number of DEGs was substantially higher:149 in LD and 201 in SSC. The DEGs in LD of modern breed reflect a shift toward faster growing, more glycolytic muscle with distinct immune and neuromuscular regulation. Likewise, the SSC of modern breed shows stronger activation of growth related and metabolic signaling pathways. Keywords: RNA sequences, pigs, adipose tissue, high throughput sequencing, RNA, transcriptome profiles, gene expression, locel breeds, hybrid breeds, data, data set Published in DiRROS: 24.03.2026; Views: 184; Downloads: 33
Full text (361,31 KB) This document has many files! More... |
4. Adipose tissue transcriptome profiles of local Krškopolje pig and modern hybrid pigs receving reduced protein diets: : expression profiling by high throughput sequencingKlavdija Poklukar Žnidaršič, Martin Škrlep, Marjeta Čandek-Potokar, 2025, complete scientific database of research data Abstract: The main objective of this study was to investigate transcriptomic differences in subcutaneous adipose tissue between the local Krškopolje pig and modern hybrids, and to assess the effects of a reduced-protein diet in both genotypes. Comparative analysis between Krškopolje pigs and modern crossbreeds revealed 375 differentially expressed genes, with 189 upregulated and 186 downregulated in Krškopolje pigs. The upregulated genes were enriched in processes related to adipogenesis (SLC7A10, ADIRF, INHBB, SFRP2), extracellular matrix remodeling (COL6A5, COL4A5, COL2A1), calcium signaling (TRPV4, CAMK2A, CALML5), pro-inflammatory cytokines (IL1A, TNFSF9, CXCL8, PTGS2), and cholesterol metabolism (CYP1A1, CYP2A19, CES1). In contrast, the reduced-protein diet induced only minor transcriptional changes at the individual gene level in both Krškopolje pigs and modern crossbreeds. Keywords: RNA sequences, pigs, adipose tissue, high throughput sequencing, RNA, transcriptome profiles, gene expression, locel breeds, hybrid breeds, data, data set Published in DiRROS: 20.03.2026; Views: 220; Downloads: 37
Full text (481,38 KB) This document has many files! More... |
5. Effect of rearing systems and dietary protein levels on the skeletal muscle histomorphology and transcriptome profiles in indigenous Krškopolje Pig : expression profiling by high throughput sequencingMilka Vrecl, Martin Škrlep, Klavdija Poklukar Žnidaršič, Gregor Fazarinc, Malan Štrbenc, Marjeta Čandek-Potokar, 2026, complete scientific database of research data Abstract: RNA-sequencing was performed to compare the effects of production systems (conventional indoor, vs. outdoor) and diets (standard protein, HP, vs. low protein, LP) within each rearing system on the transcriptomic profiles of two skeletal muscles, longissimus thoracis et lumborum (LL) and semispinalis capitis (SSC). Differentially expressed genes (DEGs) were identified with |log2FC| > 1 and q-value < 0.005. The effect of rearing system (indoor vs. outdoor) resulted in 354 DEGs in LL and 334 DEGs in SSC. The effect of diet within individual rearing system was less pronounced. In LL, HP vs. LP resulted in 23 DEGs under indoor and none under outdoor rearing system. In SSC, HP vs. LP resulted in 28 DEGs under indoor and 30 DEGs under outdoor rearing system. Pathway enrichment analysis revealed 42 pathways significantly enriched for HP vs. LP in the outdoor rearing system (across both muscles); additionally, 11 pathways were altered specifically in SSC for HP vs. LP under indoor rearing system, and 27 pathways were altered specifically in LL for HP vs. LP under outdoor rearing system Keywords: RNA sequences, pigs, adipose tissue, high throughput sequencing, RNA, transcriptome profiles, gene expression, locel breeds, hybrid breeds, data, data set Published in DiRROS: 20.03.2026; Views: 217; Downloads: 40
Full text (726,02 KB) This document has many files! More... |
6. Integrating metabolic and gene expression profiling of glucosinolate biosynthesis under drought stress in Brassica oleraceaHajer Ben Ammar, Souhir Kabtni, Donata Arena, Marwen Amari, Nicolas Al Achkar, Ferdinando Branca, Sonia Marghali, 2026, original scientific article Abstract: Drought stress induces pronounced metabolic and transcriptional reprogramming of glucosinolate (GLS) biosynthesis in Brassica oleracea. An integrative approach combining HPLC-based quantification of individual GLSs, quantitative real-time PCR of core biosynthetic and regulatory genes, correlation-based network analysis, and in silico promoter characterization was applied to evaluate drought responses across genetically diverse accessions. Drought triggered strong, accession-specific shifts in GLS composition, with sinigrin content increasing from 35.9% to 55.1% in BR1 and glucoerucin reaching up to 80.2% in CCP1, while indolic GLSs such as glucobrassicin and neoglucobrassicin accounted for >75% of total GLSs in CV2 and CCP3. Hierarchical clustering separated accessions into four distinct drought response clusters independent of morphotype. Correlation analysis revealed drought-induced rewiring of GLS interdependencies, characterized by strengthened positive associations among aliphatic GLSs (r > 0.75). Gene expression profiling identified conserved MYB-centered regulatory modules (MYB28, MYB29, MYB34, MYB122) alongside strong accession-specific induction of CYP79F1 (up to 6.3-fold), FMOGS-OX5 (up to 4.8-fold), and ST5a (up to 5.1-fold). Promoter analysis revealed enrichment of ABA- and stress-responsive cis-regulatory elements. These findings delineate a genotype-dependent regulatory framework underlying GLS plasticity and identify quantitative metabolic and transcriptional markers relevant for breeding drought-resilient Brassica cultivars. Keywords: glucosinolate, gene expression, qRT-PCR, abiotic stress response, drought stress Published in DiRROS: 10.02.2026; Views: 519; Downloads: 180
Full text (2,15 MB) This document has many files! More... |
7. Previous short-term disuse dictates muscle geneexpression and physiological adaptations to subsequentresistance exerciseMartino V. Franchi, Julián Candia, Fabio Sarto, Giuseppe Sirago, Giacomo Valli, Matteo Paganini, Lisa M. Hartnell, Emiliana Giacomello, Luana Toniolo, Elena Monti, Leonardo Nogara, Tatiana Moro, Marco Vincenzo Narici, 2025, original scientific article Abstract: Short-term unloading experienced following injury or hospitalisation induces muscle atrophy and weakness. The effects of exercise following unloading have been scarcely investigated. We investigated the functional and molecular adaptations to a resistance training (RT) programme following short-term unloading. Eleven males (22.09 ± 2.91 years) underwent 10 days of unilateral lower limb suspension (ULLS) followed by 21 days of knee extensor RT (three times/week). Data collection occurred at Baseline (LS0), after ULLS (LS10) and at active recovery (AR21). Knee extensor maximum voluntary contraction (MVC) was evaluated. Quadriceps volume was estimated by ultrasonography. Muscle fibre cross-sectional area, fibre type distribution, glycogen content and succinate dehydrogenase (SDH) activity were measured from vastus lateralis biopsies. Mitochondrial-related proteins were quantified by western blot and transcriptional responses were assessed by RNA sequencing. Following ULLS, quadriceps volume and MVC decreased significantly (3.7%, P < 0.05; 29.3%, P < 0.001). At AR21 (vs. LS10), MVC was fully restored (42%) and quadriceps volume increased markedly (18.6%, P < 0.001). Glycogen content and whole-body water increased at AR21 (14%, P < 0.001; 3.1%, P < 0.05). We observed a marked increase in fibre type I at AR21 (38%, P < 0.05). SDH immunoreactivity increased significantly after exercise (20%, P < 0.001). Mitochondrial fusion (MFN1, MFN2 and OPA1) and fission (DRP1) proteins were markedly increased by RT, and the most differentially expressed genes belonged to oxidative phosphorylation pathways. In contrast with what is usually observed after RT, oxidative metabolism, slow fibre type and mitochondrial dynamics were enhanced beyond expected. We propose that prior exposure to short-term muscle unloading may drive the nature of molecular adaptations to subsequent RT. Keywords: exercise physiology, gene expression, muscle adaptation, muscle atrophy, musle physiology, muscle plasticity, resistance training, unloading responses Published in DiRROS: 13.01.2026; Views: 286; Downloads: 221
Full text (4,75 MB) This document has many files! More... |
8. Deregulations of DNA damage-responsive genes, genes involved in the endocrine system, in an advanced in vitro 3D zebrafish hepatic cell model after exposure to Bisphenol A (BPA) and its emerging alternatives BPAF, BPAP and BPPH : version v1Alja Štern, Špela Rozman, Bojana Žegura, 2025, complete scientific database of research data Abstract: Bisphenol AF (BPAF), Bisphenol AP (BPAP), and Bisphenol PH (BPPH) are being introduced into consumer products to replace BPA and are subsequently detected in ecosystems. This study investigates the genotoxic and endocrine-disruptive potential of these emerging bisphenols using a 3D in vitro liver spheroid model derived from Danio rerio (ZFL cell line), on the transcriptional level. The selected genes targeted DNA damage response pathways (TP53, NER, BER) and endocrine-related signalling (Table 1). ZFL spheroids were prepared by a force floating method as described by Štampar et al. (2019)1. Four-day-old ZFL spheroids were exposed to BPA (50 and 200 µM), BPAF (25 and 100 µM), BPAP (25 and 100 µM), and BPPH (10 and 50 µM) for 24 (Table 2) and 96 (Table 3) hours. Following exposure, total RNA was extracted using the RNeasy® Mini Kit (Qiagen, Germany). RNA quality and quantity were assessed prior to reverse transcription (Applied Biosystems, USA), followed by gene-specific preamplification (TATAA PreAmp GrandMasterMix, Tataa Biocenter, Sweden). Gene expression analysis was performed using TaqMan Gene Expression Assays (Applied Biosystems, USA) on the Fluidigm One 48.48 Dynamic Array IFC microfluidic system as described by Štern et al. (2024)2. The generated data was analysed using the Fluidigm Gene Expression Analysis Software and with a free-accessible web program, quantGenious3. The difference in gene expression greater than 1.5-fold was considered a biologically important up/downregulation (relative expression >1.5 or <0.66, respectively). Statistically significant differences were analysed using ANOVA and Dunnett’s multiple comparison test in GraphPad Prism v9 (GraphPad Software, San Diego, CA, USA). Keywords: bisphenols, genotoxicity, endocrine disruption, ZFL spheroids, gene expression, data Published in DiRROS: 30.09.2025; Views: 582; Downloads: 419
Full text (1,09 MB) This document has many files! More... |
9. The effect of reduced dietary protein on adipose tissue in local krškopolje pigsKlavdija Poklukar Žnidaršič, Marjeta Čandek-Potokar, Milka Vrecl, Jana Brankovič, Matjaž Uršič, Martin Škrlep, 2025, original scientific article Abstract: The Slovenian autochthonous breed, Krškopolje pig, is known for high fatness and better adaptability to different environmental conditions and feed resources. However, the metabolic processes underlying these adaptations, especially in response to different diets, have not yet been studied. A deeper understanding of these mechanisms could provide valuable insights into the breed’s adaptability to different environmental conditions. Therefore, the main objective of this study was to evaluate the effect of a low-protein (LP) diet on adipose tissue in Krškopolje pigs reared in either organic outdoor (n = 2 × 12) or conventional indoor (n = 2 × 14) systems. In the outdoor system, the LP diet had no effect on adipocyte size compared to the control (high-protein) diet, while it increased lipogenic enzyme activities and monounsaturated fatty acid content, and decreased polyunsaturated fatty acid content (p < 0.05). RNA sequencing revealed the upregulation of 28 genes and the downregulation of 37 genes. The upregulated genes were mainly involved in lipid metabolism (ACLY, FASN, ACACA, MOGAT2), oxidative stress, and mitochondrial function. In the indoor system, pigs on the LP diet had smaller adipocytes (p < 0.05), whereas no differences were detected in the lipogenic enzyme activities or fatty acid composition (p > 0.10). RNA sequencing revealed 30 upregulated and 28 downregulated genes. In the indoor system, heat shock proteins (HSP70.2, HSPA6) were upregulated in pigs on the LP diet, while genes involved in the innate immune system (MSR1, TREM2, CSF3R) were downregulated. To conclude, the present study showed that LP diet affected adipose tissue metabolism and gene expression in Krškopolje pigs, with different transcriptomic responses observed in outdoor and indoor rearing conditions. Keywords: pig, local breed, low-protein diet, indoor system, outdoor system, subcutaneous adipose tissue, RNA-sequencing, nCounter gene expression assay Published in DiRROS: 12.05.2025; Views: 986; Downloads: 595
Full text (1,41 MB) This document has many files! More... |
10. An international comparability study on quantification of mRNA gene expression ratios: CCQM-P103.1Alison S. Devonshire, Tina Demšar, Jana Žel, Andrej Blejec, Mojca Milavec, 2016, original scientific article Abstract: Measurement of RNA can be used to study and monitor a range of infectious and non-communicable diseases, with profiling of multiple gene expression mRNA transcripts being increasingly applied to cancer stratification and prognosis. An international comparison study (Consultative Committee for Amount of Substance (CCQM)-P103.1) was performed in order to evaluate the comparability of measurements of RNA copy number ratio for multiple gene targets between two samples. Six exogenous synthetic targets comprising of External RNA Control Consortium (ERCC) standards were measured alongside transcripts for three endogenous gene targets present in the background of human cell line RNA. The study was carried out under the auspices of the Nucleic Acids (formerly Bioanalysis) Working Group of the CCQM. It was coordinated by LGC (United Kingdom) with the support of National Institute of Standards and Technology (USA) and results were submitted from thirteen National Metrology Institutes and Designated Institutes. The majority of laboratories performed RNA measurements using RT-qPCR, with datasets also being submitted by two laboratories based on reverse transcription digital polymerase chain reaction and one laboratory using a next-generation sequencing method. In RT-qPCR analysis, the RNA copy number ratios between the two samples were quantified using either a standard curve or a relative quantification approach. In general, good agreement was observed between the reported results of ERCC RNA copy number ratio measurements. Measurements of the RNA copy number ratios for endogenous genes between the two samples were also consistent between the majority of laboratories. Some differences in the reported values and confidence intervals (‘measurement uncertainties’) were noted which may be attributable to choice of measurement method or quantification approach. This highlights the need for standardised practices for the calculation of fold change ratios and uncertainties in the area of gene expression profiling.
Keywords: RNA copy number ratio, RT-qPCR, gene expression, normalisation, standardisation, molecular diagnostic, transcriptomics, cancer, diagnostics, biomarker identification and validation Published in DiRROS: 04.03.2025; Views: 914; Downloads: 484
Full text (1,78 MB) This document has many files! More... |