A P2RY12 deficiency results in sex-specific cellular perturbations and sexually dimorphic behavioral anomalies.

Microglia are sexually dimorphic, yet, this critical aspect is often overlooked in neuroscientific studies. Decades of research have revealed the dynamic nature of microglial-neuronal interactions, but seldom consider how this dynamism varies with microglial sex differences, leaving a significant gap in our knowledge. This study focuses on P2RY12, a highly expressed microglial signature gene that mediates microglial-neuronal interactions, we show that adult females have a significantly higher expression of the receptor than adult male microglia. We further demonstrate that a genetic deletion of P2RY12 induces sex-specific cellular perturbations with microglia and neurons in females more significantly affected. Correspondingly, female mice lacking P2RY12 exhibit unique behavioral anomalies not observed in male counterparts. These findings underscore the critical, sex-specific roles of P2RY12 in microglial-neuronal interactions, offering new insights into basal interactions and potential implications for CNS disease mechanisms.

RIPK3 and RIPK1 gene expression in pterygium: unveiling molecular insights into pathogenesis.

BACKGROUND: Pterygium, characterized by the abnormal proliferation of epithelial cells, matrix remodeling, vascularization, and lesion migration, is a prevalent ocular surface disease involving the growth of fibrovascular tissue on the cornea. Despite the unclear underlying causes of pterygium, numerous investigations have indicated the involvement of cell death pathways in the regulation of cell cycle dynamics. Consequently, the objective of this study was to assess the expression levels of necroptosis markers in individuals diagnosed with pterygium, aiming to shed light on the potential role of necroptosis in the pathogenesis of this condition. METHODS: This study aimed to investigate the expression patterns of receptor-interacting serine/threonine kinase 3 (RIPK3) and receptor-interacting serine/threonine kinase 1 (RIPK1) genes in pterygium tissues. 41 patients undergoing pterygium excision surgery were recruited. Resected pterygium samples and normal conjunctival tissues were collected, and RIPK3 and RIPK1 mRNA levels were measured using quantitative real-time PCR. RESULTS: Our findings reveal that the expression of RIPK3 is significantly increased in samples obtained from individuals with pterygium. However, no significant alterations were observed in the expression of RIPK1 in these samples. Results showed significantly higher RIPK3 expression in pterygium tissues compared to controls. Moreover, increased RIPK3 levels correlated negatively with pterygium recurrence rates. CONCLUSIONS: These findings suggest RIPK3 may play a protective role against pterygium recurrence through necroptosis.

Transcription bodies regulate gene expression by sequestering CDK9.

The localization of transcriptional activity in specialized transcription bodies is a hallmark of gene expression in eukaryotic cells. It remains unclear, however, if and how transcription bodies affect gene expression. Here we disrupted the formation of two prominent endogenous transcription bodies that mark the onset of zygotic transcription in zebrafish embryos and analysed the effect on gene expression using enriched SLAM-seq and live-cell imaging. We find that the disruption of transcription bodies results in the misregulation of hundreds of genes. Here we focus on genes that are upregulated. These genes have accessible chromatin and are poised to be transcribed in the presence of the two transcription bodies, but they do not go into elongation. Live-cell imaging shows that disruption of the two large transcription bodies enables these poised genes to be transcribed in ectopic transcription bodies, suggesting that the large transcription bodies sequester a pause release factor. Supporting this hypothesis, we find that CDK9-the kinase that releases paused polymerase II-is highly enriched in the two large transcription bodies. Overexpression of CDK9 in wild-type embryos results in the formation of ectopic transcription bodies and thus phenocopies the removal of the two large transcription bodies. Taken together, our results show that transcription bodies regulate transcription by sequestering machinery, thereby preventing genes elsewhere in the nucleus from being transcribed.

Von Hippel Lindau tumor suppressor controls m6A-dependent gene expression in renal tumorigenesis.

N6-Methyladenosine (m6A) is the most abundant posttranscriptional modification, and its contribution to cancer evolution has recently been appreciated. Renal cancer is the most common adult genitourinary cancer, approximately 85% of which is accounted for by the clear cell renal cell carcinoma (ccRCC) subtype characterized by VHL loss. However, it is unclear whether VHL loss in ccRCC affects m6A patterns. In this study, we demonstrate that VHL binds and promotes METTL3/METTL14 complex formation while VHL depletion suppresses m6A modification, which is distinctive from its canonical E3 ligase role. m6A RNA immunoprecipitation sequencing (RIP-Seq) coupled with RNA-Seq allows us to identify a selection of genes whose expression may be regulated by VHL-m6A signaling. Specifically, PIK3R3 is identified to be a critical gene whose mRNA stability is regulated by VHL in a m6A-dependent but HIF-independent manner. Functionally, PIK3R3 depletion promotes renal cancer cell growth and orthotopic tumor growth while its overexpression leads to decreased tumorigenesis. Mechanistically, the VHL-m6A-regulated PIK3R3 suppresses tumor growth by restraining PI3K/AKT activity. Taken together, we propose a mechanism by which VHL regulates m6A through modulation of METTL3/METTL14 complex formation, thereby promoting PIK3R3 mRNA stability and protein levels that are critical for regulating ccRCC tumorigenesis.

Powerful and accurate detection of temporal gene expression patterns from multi-sample multi-stage single-cell transcriptomics data with TDEseq.

We present a non-parametric statistical method called TDEseq that takes full advantage of smoothing splines basis functions to account for the dependence of multiple time points in scRNA-seq studies, and uses hierarchical structure linear additive mixed models to model the correlated cells within an individual. As a result, TDEseq demonstrates powerful performance in identifying four potential temporal expression patterns within a specific cell type. Extensive simulation studies and the analysis of four published scRNA-seq datasets show that TDEseq can produce well-calibrated p-values and up to 20% power gain over the existing methods for detecting temporal gene expression patterns.

Rhes, a striatal enriched protein, regulates post-translational small-ubiquitin-like-modifier (SUMO) modification of nuclear proteins and alters gene expression.

Rhes (Ras homolog enriched in the striatum), a multifunctional protein that regulates striatal functions associated with motor behaviors and neurological diseases, can shuttle from cell to cell via the formation of tunneling-like nanotubes (TNTs). However, the mechanisms by which Rhes mediates diverse functions remain unclear. Rhes is a small GTPase family member which contains a unique C-terminal Small Ubiquitin-like Modifier (SUMO) E3-like domain that promotes SUMO post-translational modification of proteins (SUMOylation) by promoting "cross-SUMOylation" of the SUMO enzyme SUMO E1 (Aos1/Uba2) and SUMO E2 ligase (Ubc-9). Nevertheless, the identity of the SUMO substrates of Rhes remains largely unknown. Here, by combining high throughput interactome and SUMO proteomics, we report that Rhes regulates the SUMOylation of nuclear proteins that are involved in the regulation of gene expression. Rhes increased the SUMOylation of histone deacetylase 1 (HDAC1) and histone 2B, while decreasing SUMOylation of heterogeneous nuclear ribonucleoprotein M (HNRNPM), protein polybromo-1 (PBRM1) and E3 SUMO-protein ligase (PIASy). We also found that Rhes itself is SUMOylated at 6 different lysine residues (K32, K110, K114, K120, K124, and K245). Furthermore, Rhes regulated the expression of genes involved in cellular morphogenesis and differentiation in the striatum, in a SUMO-dependent manner. Our findings thus provide evidence for a previously undescribed role for Rhes in regulating the SUMOylation of nuclear targets and in orchestrating striatal gene expression via SUMOylation.

Immunostimulation of Asian elephant (Elephas maximus) blood cells by parapoxvirus ovis and CpG motif-containing bacterial plasmid DNA upregulates innate immune gene expression.

The immune system of Asian elephants (Elephas maximus) is poorly studied, compared to that of livestock, rodents or humans. The innate immune response has become a focus of interest in relation to Elephant endotheliotropic herpesviruses (EEHVs). EEHVs cause a fatal hemorrhagic disease (EEHV-HD) and are a significant threat to captive Asian elephant populations worldwide. Similar to other herpesvirus infections, nearly all animals become infected, but only some develop disease. As progression to EEHV-HD is often acute, a robust innate immune response is crucial to control EEHV infections. This is invariably true of the host in the first instance, but it can also potentially be modulated by intervention strategies. Here, two immunostimulant veterinary medicinal products, authorized for use in domestic species, were tested for their ability to induce innate anti-viral immune responses in Asian elephant blood cells. Sequence data were obtained for a range of previously unidentified Asian elephant immune genes, including C-X-C motif chemokine ligand 10 (CXCL10), interferon stimulated gene 15 (ISG15) and myxovirus GTPase 1 (Mx1), and were employed in the design of species-specific qPCR assays. These assays were subsequently used in analyses to determine fold changes in gene expression over a period of 24 hours. This study demonstrates that both immunostimulant medications are capable of inducing significant innate anti-viral immune responses which suggests that both could be beneficial in controlling EEHV infections in Asian elephants.

Insertion of a neomycin selection cassette in the Amigo1 locus alters gene expression in the olfactory epithelium leading to region-specific defects in olfactory receptor neuron development.

During development of the nervous system, neurons connect to one another in a precisely organized manner. Sensory systems provide a good example of this organization, whereby the composition of the outside world is represented in the brain by neuronal maps. Establishing correct patterns of neural circuitry is crucial, as inaccurate map formation can lead to severe disruptions in sensory processing. In rodents, olfactory stimuli modulate a wide variety of behaviors essential for survival. The formation of the olfactory glomerular map is dependent on molecular cues that guide olfactory receptor neuron axons to broad regions of the olfactory bulb and on cell adhesion molecules that promote axonal sorting into specific synaptic units in this structure. Here, we demonstrate that the cell adhesion molecule Amigo1 is expressed in a subpopulation of olfactory receptor neurons, and we investigate its role in the precise targeting of olfactory receptor neuron axons to the olfactory bulb using a genetic loss-of-function approach in mice. While ablation of Amigo1 did not lead to alterations in olfactory sensory neuron axonal targeting, our experiments revealed that the presence of a neomycin resistance selection cassette in the Amigo1 locus can lead to off-target effects that are not due to loss of Amigo1 expression, including unexpected altered gene expression in olfactory receptor neurons and reduced glomerular size in the ventral region of the olfactory bulb. Our results demonstrate that insertion of a neomycin selection cassette into the mouse genome can have specific deleterious effects on the development of the olfactory system and highlight the importance of removing antibiotic resistance cassettes from genetic loss-of-function mouse models when studying olfactory system development.

A novel insight into the neuroprotective effects of cannabidiol: maintained apelin/dopamine synthesis, NRF2 signaling, and AKT/CREB/BDNF gene expressions.

Neuroinflammation is a process associated with degeneration and loss of neurons in different parts of the brain. The most important damage mechanisms in its formation are oxidative stress and inflammation. This study aimed to investigate the protective effects of cannabidiol (CBD) against neuroinflammation through various mechanisms. Thirty­two female rats were randomly divided into 4 groups as control, lipopolysaccharide (LPS), LPS + CBD and CBD groups. After six hours following LPS administration, rats were sacrificed, brain and cerebellum tissues were obtained. Tissues were stained with hematoxylin­eosin for histopathological analysis. Apelin and tyrosine hydroxylase synthesis were determined immunohistochemically. Total oxidant status and total antioxidant status levels were measured, and an oxidative stress index was calculated. Protein kinase B (AKT), brain-derived neurotrophic factor (BDNF), cyclic­AMP response element­binding protein (CREB) and nuclear factor erythroid 2­related factor 2 (NRF2) mRNA expression levels were also determined. In the LPS group, hyperemia, degeneration, loss of neurons and gliosis were seen in all three tissues. Additionally, Purkinje cell loss in the cerebellum, as well as neuronal loss in the cerebral cortex and hippocampus, were found throughout the LPS group. The expressions of AKT, BDNF, CREB and NRF2, apelin and tyrosine hydroxylase synthesis all decreased significantly. CBD treatment reversed these changes and ameliorated oxidative stress parameters. CBD showed protective effects against neuroinflammation via regulating AKT, CREB, BDNF expressions, NRF2 signaling, apelin and tyrosine hydroxylase synthesis.

Profiling of long non-coding RNAs in hippocampal-entorhinal system subfields: impact of RN7SL1 on neuroimmune response modulation in Alzheimer's disease.

Alzheimer's disease (AD) is recognized as the predominant cause of dementia, and neuroimmune processes play a pivotal role in its pathological progression. The involvement of long non-coding RNAs (lncRNAs) in AD has attracted widespread attention. Herein, transcriptomic analysis of 262 unique samples extracted from five hippocampal-entorhinal system subfields of individuals with AD pathology and without AD pathology revealed distinctive lncRNA expression profiles. Through differential expression and coexpression analyses, we identified 16 pivotal lncRNAs. Notably, RN7SL1 knockdown significantly modulated microglial responses upon oligomeric amyloid-ß stimulation, resulting in a considerable decrease in proinflammatory cytokine production and subsequent neuronal damage. These findings highlight RN7SL1 as an essential neuroimmune-related lncRNA that could serve as a prospective target for AD diagnosis and treatment.

Circadian disruption dysregulates lung gene expression associated with inflammatory lung injury.

Rationale: Circadian systems drive the expression of multiple genes in nearly all cells and coordinate cellular-, tissue-, and system-level processes that are critical to innate immunity regulation. Objective: We examined the effects of circadian rhythm disorganization, produced by light shift exposure, on innate immunity-mediated inflammatory lung responses including vascular permeability and gene expression in a C57BL/6J murine model of inflammatory lung injury. Methods: A total of 32 C57BL/6J mice were assigned to circadian phase shifting (CPS) with intratracheal phosphate-buffered saline (PBS), CPS with intratracheal lipopolysaccharide (LPS), control (normal lighting) condition with intratracheal PBS, and control condition with intratracheal LPS. Bronchoalveolar lavage (BAL) protein, cell counts, tissue immunostaining, and differentially expressed genes (DEGs) were measured in lung tissues at 2 and 10 weeks. Measurements and results: In mice exposed to both CPS and intratracheal LPS, both BAL protein and cell counts were increased at both 2 and 10 weeks compared to mice exposed to LPS alone. Multiple DEGs were identified in CPS-LPS-exposed lung tissues compared to LPS alone and were involved in transcriptional pathways associated with circadian rhythm disruption, regulation of lung permeability, inflammation with Rap1 signaling, and regulation of actin cytoskeleton. The most dysregulated pathways included myosin light chain kinase, MAP kinase, profilin 2, fibroblast growth factor receptor, integrin b4, and p21-activated kinase. Conclusion: Circadian rhythm disruption results in exacerbated immune response and dysregulated expression of cytoskeletal genes involved in the regulation of epithelial and vascular barrier integrity-the mechanistic underpinnings of acute lung injury. Further studies need to explore circadian disorganization as a druggable target.

Penicillium-Infected Apples Benefit Larval Development of Conogethes punctiferalis via Alterations of Their Gut Bacteria Community and Gene Expression.

Pathogenic microorganisms can impact the behavior and physiology of herbivores by direct or indirect means. This study demonstrated that yellow peach moth Conogethes punctiferalis larvae feeding on Penicillium-infected apples exhibited significantly longer body length and weight parameters compared to the control group. The sequencing of gut 16S rRNA showed a significant increase in the diversity and abundance of bacteria in the larvae feeding on Penicillium-infected apples. Additionally, transcriptomic sequencing of the larval gut indicated significant upregulation of genes related to digestion and cuticle formation after consuming Penicillium-infected apples. Furthermore, enzyme activity assays revealed notable changes in the trypsin and lipase activity. Consequently, these alterations in gut microbiota structure, diversity, and gene expression levels may underlie the observed growth and developmental variations in C. punctiferalis larvae mediated by pathogenic microorganisms. This study holds theoretical significance for a deeper understanding of the tripartite interaction among microorganisms, insects, and plants as well as for the development of novel pest control measures based on gut microbiota.

Expanding the Drosophila toolkit for dual control of gene expression.

The ability to independently control gene expression in two different tissues in the same animal is emerging as a major need, especially in the context of inter-organ communication studies. This type of study is made possible by technologies combining the GAL4/UAS and a second binary expression system such as the LexA system or QF system. Here, we describe a resource of reagents that facilitate combined use of the GAL4/UAS and a second binary system in various Drosophila tissues. Focusing on genes with well-characterized GAL4 expression patterns, we generated a set of more than 40 LexA-GAD and QF2 insertions by CRISPR knock-in and verified their tissue specificity in larvae. We also built constructs that encode QF2 and LexA-GAD transcription factors in a single vector. Following successful integration of this construct into the fly genome, FLP/FRT recombination is used to isolate fly lines that express only QF2 or LexA-GAD. Finally, using new compatible shRNA vectors, we evaluated both LexA and QF systems for in vivo gene knockdown and are generating a library of such RNAi fly lines as a community resource. Together, these LexA and QF system vectors and fly lines will provide a new set of tools for researchers who need to activate or repress two different genes in an orthogonal manner in the same animal.

In order for researchers to understand how organisms develop and function, they often switch specific genes on or off in certain tissues or at selected times. This can be achieved using genetic tools called binary expression systems. In the fruit fly ­ a popular organism for studying biological processes ­ the most common is the GAL4/UAS system. In this system, a protein called GAL4 is expressed in a specific organ or tissue where it activates a UAS element ­ a genetic sequence that is inserted in front of the gene that is to be switched on. This can also include genes inserted into the fruit fly encoding fluorescent proteins or stretches of DNA coding for factors that can silence specific genes. For example, fruit flies expressing GAL4 protein specifically in nerve cells and a UAS element in front of a gene for a fluorescent protein will display fluorescent nerve cells, which can then be examined using fluorescence microscopy. Studying how organs communicate with one other can require controlled expression of multiple genes at the same time. In fruit flies, other binary expression systems that are analogous to the GAL4/UAS system (known as LexA/LexAop and QF/QUAS) can be used in tandem. For example, to study gut-brain communication, the GAL4/UAS system might be used to switch on the gene for an insulin-like protein in the gut, with one of the other systems controlling the expression of its corresponding receptor in the brain. However, these experiments are currently difficult because, while there are thousands of GAL4/UAS genetic lines, there are only a few LexA/LexAop and QF/QUAS genetic lines. To address this lack of resources, Zirin et al. produced a range of genetically engineered fruit flies containing the LexA/LexAop and QF/QUAS binary expression systems. The flies expressed LexA or QF in each of the major fly organs, including the brain, heart, muscles, and gut. A fluorescent reporter gene linked to the LexAop or QUAS elements, respectively, was then used to test the specificity to single organs and compare the different systems. In some organs the LexA/LexAop system was more reliable than the QF/QUAS system. However, both systems could be successfully combined with genetic elements to switch on a fluorescent reporter gene or switch off a gene of interest in the intended organ. The resources developed by Zirin et al. expand the toolkit for studying fruit fly biology. In future, it will be important to understand the differences between GAL4, LexA and QF systems, and to increase the number of fruit fly lines containing the newer binary expression systems.

Gene Expression Patterns Associated with Survival in Glioblastoma.

The aim of this study was to investigate gene expression alterations associated with overall survival (OS) in glioblastoma (GBM). Using the Nanostring nCounter platform, we identified four genes (COL1A2, IGFBP3, NGFR, and WIF1) that achieved statistical significance when comparing GBM with non-neoplastic brain tissue. The four genes were included in a multivariate Cox Proportional Hazard model, along with age, extent of resection, and O6-methylguanine-DNA methyltransferase (MGMT) promotor methylation, to create a unique glioblastoma prognostic index (GPI). The GPI score inversely correlated with survival: patient with a high GPI had a median OS of 7.5 months (18-month OS = 9.7%) whereas patients with a low GPI had a median OS of 20.1 months (18-month OS = 54.5%; log rank p-value = 0.004). The GPI score was then validated in 188 GBM patients from The Cancer Genome Atlas (TCGA) from a national data base; similarly, patients with a high GPI had a median OS of 10.5 months (18-month OS = 12.4%) versus 16.9 months (18-month OS = 41.5%) for low GPI (log rank p-value = 0.0003). We conclude that this novel mRNA-based prognostic index could be useful in classifying GBM patients into risk groups and refine prognosis estimates to better inform treatment decisions or stratification into clinical trials.

Claudin-10 Expression and the Gene Expression Pattern of Thick Ascending Limb Cells.

Many genomic, anatomical and functional differences exist between the medullary (MTAL) and the cortical thick ascending limb of the loop of Henle (CTAL), including a higher expression of claudin-10 (CLDN10) in the MTAL than in the CTAL. Therefore, we assessed to what extent the Cldn10 gene expression is a determinant of differential gene expression between MTAL and CTAL. RNAs extracted from CTAL and MTAL microdissected from wild type (WT) and Cldn10 knock out mice (cKO) were analyzed by RNAseq. Differential and enrichment analyses (GSEA) were performed with interactive R Shiny software. Between WT and cKO MTAL, 637 genes were differentially expressed, whereas only 76 were differentially expressed between WT and cKO CTAL. Gene expression patterns and GSEA analyses in all replicates showed that WT MTAL did not cluster with the other replicates; no hierarchical clustering could be found between WT CTAL, cKO CTAL and cKO MTAL. Compared to WT replicates, cKO replicates were enriched in Cldn16, Cldn19, Pth1r, (parathyroid hormone receptor type 1), Casr (calcium sensing receptor) and Vdr (Vitamin D Receptor) mRNA in both the cortex and medulla. Cldn10 is associated with gene expression patterns, including genes specifically involved in divalent cations reabsorption in the TAL.

High NaCl Concentrations in Water Are Associated with Developmental Abnormalities and Altered Gene Expression in Zebrafish.

Salt is frequently introduced in ecosystems, where it acts as a pollutant. This study examined how changes in salinity affect the survival and development of zebrafish from the two-cell to the blastocyst stage and from the blastocyst to the larval stage. Control zebrafish embryos were cultured in E3 medium containing 5 mM Sodium Chloride (NaCl), 0.17 mM Potassium Chloride (KCL), 0.33 mM Calcium Chloride (CaCl2), and 0.33 mM Magnesium Sulfade (MgSO4). Experiments were conducted using increasing concentrations of each individual salt at 5×, 10×, 50×, and 100× the concentration found in E3 medium. KCL, CaCl2, and MgSO4 did not result in lethal abnormalities and did not affect early embryo growth at any of the concentrations tested. Concentrations of 50× and 100× NaCl caused embryonic death in both stages of development. Concentrations of 5× and 10× NaCl resulted in uninflated swim bladders in 12% and 65% of larvae, compared to 4.2% of controls, and caused 1654 and 2628 genes to be differentially expressed in blastocysts, respectively. The ATM signaling pathway was affected, and the Sonic Hedgehog pathway genes Shh and Ptc1 implicated in swim bladder development were downregulated. Our findings suggest that increased NaCl concentrations may alter gene expression and cause developmental abnormalities in animals found in affected ecosystems.

Estrogen α and ß Receptor Expression in the Various Regions of Resected Glioblastoma Multiforme Tumors and in an In Vitro Model.

Glioblastoma multiforme (GBM) is a malignant tumor with a higher prevalence in men and a higher survival rate in transmenopausal women. It exhibits distinct areas influenced by changing environmental conditions. This study examines how these areas differ in the levels of estrogen receptors (ERs) which play an important role in the development and progression of many cancers, and whose expression levels are often correlated with patient survival. This study utilized two research models: an in vitro model employing the U87 cell line and a second model involving tumors resected from patients (including tumor core, enhancing tumor region, and peritumoral area). ER expression was assessed at both gene and protein levels, with the results validated using confocal microscopy and immunohistochemistry. Under hypoxic conditions, the U87 line displayed a decrease in ERß mRNA expression and an increase in ERα mRNA expression. In patient samples, ERß mRNA expression was lower in the tumor core compared to the enhancing tumor region (only in males when the study group was divided by sex). In addition, ERß protein expression was lower in the tumor core than in the peritumoral area (only in women when the study group was divided by sex). Immunohistochemical analysis indicated the highest ERß protein expression in the enhancing tumor area, followed by the peritumoral area, and the lowest in the tumor core. The findings suggest that ER expression may significantly influence the development of GBM, exhibiting variability under the influence of conditions present in different tumor areas.

Investigating the level of vitamin D receptor gene expression in two tumoral and healthy breast tissues in breast cancer patients and its association with prognostic factors.

BACKGROUND: Breast cancer is one of the most common cancers known among women. This study aimed to investigate the level of vitamin D receptor gene expression in two tumoral and healthy breast tissues in breast cancer patients and its association with prognostic factors. METHODS: This descriptive cross-sectional study was conducted in 2022 on 50 patients with high suspicion of breast cancer who were candidates for mastectomy and lumpectomy in a learning hospital. From the patients, two tissue samples were prepared, and there was a total of 100 samples. The samples were subjected to H/E staining and evaluated by a pathologist. The presence or absence of malignancy in each sample was confirmed by two pathologists, and HER2/ER/PR indices were determined. Descriptive and analytical statistical methods and SPSS version 22 software were used. RESULTS: The average age of the patients was 51.60 ± 11.22 years old, and the average tumor size was 3.17 ± 1.28. Most tumors were grade 2 (48%). The expression of HER2, ER, and PR was positive in 24, 64, and 54%, respectively. The largest number of cases were in stage 2A. The expression level of vitamin D receptor (VDR) gene in healthy tissue (2.08 ± 1.01) was higher than tumoral tissue (0.25 ± 1.38) (P = 0.001). In tumoral and healthy tissue, VDR expression was not significant according to tumor grade, HER2, ER, PR, LVI, LN, disease stage, age, and tumor size. CONCLUSIONS: The expression level of VDR in healthy tissue was significantly higher than tumoral tissue. However, there was no significant relationship between VDR and tumor grade, HER2, ER, PR, LVI, LN, disease stage, age, and tumor size.

SEGCECO: Subgraph Embedding of Gene expression matrix for prediction of CEll-cell COmmunication.

Recent advances in single-cell RNA sequencing technology have eased analyses of signaling networks of cells. Recently, cell-cell interaction has been studied based on various link prediction approaches on graph-structured data. These approaches have assumptions about the likelihood of node interaction, thus showing high performance for only some specific networks. Subgraph-based methods have solved this problem and outperformed other approaches by extracting local subgraphs from a given network. In this work, we present a novel method, called Subgraph Embedding of Gene expression matrix for prediction of CEll-cell COmmunication (SEGCECO), which uses an attributed graph convolutional neural network to predict cell-cell communication from single-cell RNA-seq data. SEGCECO captures the latent and explicit attributes of undirected, attributed graphs constructed from the gene expression profile of individual cells. High-dimensional and sparse single-cell RNA-seq data make converting the data into a graphical format a daunting task. We successfully overcome this limitation by applying SoptSC, a similarity-based optimization method in which the cell-cell communication network is built using a cell-cell similarity matrix which is learned from gene expression data. We performed experiments on six datasets extracted from the human and mouse pancreas tissue. Our comparative analysis shows that SEGCECO outperforms latent feature-based approaches, and the state-of-the-art method for link prediction, WLNM, with 0.99 ROC and 99% prediction accuracy. The datasets can be found at https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE84133 and the code is publicly available at Github https://github.com/sheenahora/SEGCECO and Code Ocean https://codeocean.com/capsule/8244724/tree.