Methylenetetrahydrofolate reductase (MTHFR), an enzyme expressed in mammalian testes exerts direct effect on spermatogenesis; however, its protein characteristics in bovine testes remain unknown. Here, we analysed bovine testicular structure, MTHFR bioinformatics profile, mRNA, and protein expression characteristics in yellow-cattle (y-c) and yak testis using histological procedures, bioinformatics analysis, qRT-PCR, and western blot. Testes from 13 bovines, ≤ 2 years juvenile (y-c, n = 3; yak, n=3) and ≥ 3 years adult (y-c, n = 3; yak, n = 4) were collected and analysed. Anatomical characteristics of testis in y-c and yak were similar except the weight or size for which that of y-c was significantly higher or greater than yak. In y-c, an open reading frame (ORF) for 2600 nucleotides sequence, encoding 655 amino acids showed high homology with zebu cattle (99.51%) and wild yak (98.68%). Secondary and 3D protein structures were similar to that of humans with differences in number of nucleotides, amino acids, and some physico-chemical characteristics. MTHFR mRNA expression in y-c and yak were significantly higher in adult testes compared with juvenile ones. However, its protein expression was higher but not statistically significant in adult y-c and yak compared to the juvenile ones. The highlights and inferences of these and other findings are discussed.

Previous studies on the prognostic impact of survivin expression in gastrointestinal (GI) cancer have yielded inconsistent results. This study was initiated to assess the relationship between survivin expression and overall survival (OS) or disease free survival (DFS) in GI cancer patients. We applied system literature searches on EMBASE, PubMed, Web of science, and the Cochrane library to conduct this up-to-date meta-analysis. Thirty studies with totally 3622 GI cancer patients were collected. The prevalence of high survivin expression in GI cancer was 0.57 (95% CI: 0.51-0.63). High survivin expression was significantly associated with shorter OS (HR 1.57, 95% CI: 1.42-1.74) and DFS (HR 1.38, 95% CI: 1.21-1.58). Subgroup analysis also showed significant association between high survivin expression and poorer OS or DFS in gastric cancer or colorectal cancer. In summary, our study indicated that high survivin expression was related to poor prognosis in GI cancer. Well-designed studies with large sample and more convincing data are needed to confirm our conclusion.

Eimeria tenella rhoptry protein has the properties of a protective antigen. EtROP27 is a pathogenic related gene detected by transcriptome, but its expression pattern, immunogenicity and potency have been unknown. Therefore, a gene segment of EtROP27 was amplified and transplanted to pET28a prokaryotic vector for recombinant protein expression and purified to generate polyclonal antibody. Then, RT-PCR, Western blotting were performed to know the expression pattern of EtROP27. Subsequently, animal experiments were conducted to evaluate the immunoprotective effect of the recombinant protein using different immunizing doses (50, 100, 150 μg). The results showed that the expression of EtROP27 gradually increased with the prolongation of infection time, reaching the highest level after 96 hours and then decreased. And EtROP27 was a natural antigen of coccidia which can stimulate the body to produce high levels of IgY. As recombinant protein vaccines, the immune protection evaluation tests of EtROP27 showed that the average weight gain rate of the immune challenge groups were significantly higher than that of the infection control group. The average lesion scores were significantly lower than that of the control group. The oocyst excretion decreased by 81.25%, 86.21%, 80.01% and the anticoccidial index were 159.45, 171.47, 166.75. EtROP27 is a promising candidate antigen gene for development of coccidiosis vaccine.

One of the most crucial elements in deeply understanding humans on a psychological level is manifested through facial expressions. The analysis of a human behavior can be informed by their facial expressions, making it essential to employ indicators such as expression (Expr), valence-arousal (VA), and action units (AU). In this paper, we introduce the method proposed in the Challenge of the 6th Workshop and Competition on Affective Behavior Analysis in-the-wild (ABAW) at CVPR 2024. Our proposed method utilizes the multi-modal Aff-wild2 dataset, which is splitted into spatial and audio modalities. For the spatial data, we extract features using a SimMiM model that was pre-trained on a diverse set of facial expression data. For the audio data, we extract features using a WAV2VEC model. To fusion the extracted spatial and audio features, we employed the cascaded cross-attention mechanism of a transformer.

Gene expression profiling data contains more information than is routinely extracted with standard approaches. Here we present Fold-change-Specific Enrichment Analysis (FSEA), a new method for functional annotation of differentially expressed genes from transcriptome data with respect to their fold changes. FSEA identifies GO terms, which are shared by the group of genes with a similar magnitude of response, and assesses these changes. GO terms found by FSEA are fold-change-specifically (e.g. weakly, moderately or strongly) affected by a stimulus under investigation. We demonstrate that many responses to abiotic factors, mutations, treatments and diseases occur in a fold-change-specific manner. FSEA analyses suggest that there are two prevailing responses of functionally-related gene groups, either weak or strong. Notably, some of the fold-change-specific GO terms are invisible by classical algorithms for functional gene enrichment, SEA and GSEA. These are GO terms not enriched compared to the genome background but strictly regulated by a factor within specific fold-change intervals. FSEA analysis of a cancer-related transcriptome suggested that the gene groups with a tightly coordinated response can be the valuable source to search for possible regulators, markers and therapeutic targets in oncogenic processes. Availability and Implementation: FSEA is implemented as the FoldGO Bioconductor R package and a web-server https://webfsgor.sysbio.cytogen.ru/ .

Introduction: Depending on the kind of cancer and neurological disorders, there might be a direct or inverse link between neurodegeneration and cancer. The inverse interaction between cancer and Alzheimer's pathology has garnered substantial interest in medical research. While this inverse interaction between Alzheimer's pathology and cancer presents a fascinating scientific puzzle, there is still limited information and studies on the potential mechanisms involved in this inverse rela-tionship. Aim: This study aims to identify potential genes from Alzheimer's disease that may protect against cancer. Method: The Gene Expression Omnibus (GEO) database provided the Alzheimer's pathology dataset, and GEO 2R was utilized to examine genes that were differentially expressed. The expression of upregulated Alzheimer's pathology genes in breast, colorectal, and prostate cancer was compared using VENN 2.1, colorectal, and prostate cancer was compared using VENN 2.1, and colorectal and prostate cancers were compared using VENN 2.1. Results: Among the upregulated genes associated with Alzheimer's disease, 25 genes are shared among breast cancer, colorectal cancer, and prostate cancer, while 28 of these genes were exclusive to Alzheimer's dis-ease. Gene expression profiling interactive analysis among the shared genes revealed that AQP1, C4A, CSRP2, CRYAB, MACF1, HSPB2, and NEAT1 were downregulated in breast, colorectal, and prostate cancers. These genes have varying roles in cancer progression and suppression. Functional analysis of Alzheimer's disease-related upregulated genes revealed that C1orf56 and HIST1H2AC control cell proliferation, with HIST1H2AC playing a role in necroptosis. Additionally, C1orf56 worked as an oncogenic modifier to support DNMT3B's tumor suppressor properties. Conclusion: Given the role of HIST1H2AC and C1orf56 in the negative control of cell proliferation and their upregulation in Alzheimer's disease, there is a need to further investigate HIST1H2AC and C1orf56 as potential tumor suppressors responsible for Alzheimer's disease and cancer inverse association.

Real-time fluorescent quantitative PCR (RT-qPCR) is the most classic and widely used technology for evaluating the expression level of target gene. In order to select the proper internal reference genes for RT-qPCR analysis in Amorphophallus konjac (Araceae), eight candidate internal reference genes, including 25S ribosomal RNA gene (25S rRNA), 18S ribosomal RNA gene (18S rRNA), actin gene (ACT), glyceraldehyde-3-phosphate dehydrogenase gene (GAPDH), Ubiquitin gene (UBQ), β-tubulin gene (β-TUB), eukaryotic elongation factor 1-αgene(eEF-1α), and eukaryotic translation initiation factor 4α-1 gene (eIF-4α) were selected and tested the corresponding expression level in different tissues at different growing stages. The results showed that 25S rRNA, 18S rRNA, and ACT at the reproductive periods, eEF-1α and eIF-4α at the nutritional periods, and eEF-1α, UBQ, and ACT at different leaf developmental periods had a stable level of gene expression, respectively. These results might be useful for the study of gene function in A. konjac.

Caseous lymphadenitis (CL) is a chronic contagious disease that affects small ruminants and is characterized by the formation of pyogranulomas in lymph nodes and other organs. However, the pathogenesis of this disease and the response of the host genome to infection are not yet fully understood. This study aimed to investigate the whole blood transcriptome and evaluate differential gene expression during the later stages of CL in naturally infected ewes. The study included diseased, serologically positive (EP), exposed, serologically negative (EN) ewes from the same infected flock and healthy ewes (CN) from a different flock. RNA sequencing was performed using the Illumina NextSeq system, and differential gene expression was estimated using DESeq2 and Edge R approaches. The analysis identified 191 annotated differentially expressed genes (DEGs) in the EP group (102 upregulated and 89 downregulated) and 256 DEGs in the EN group (106 upregulated and 150 downregulated). Infection influenced numerous immunoregulatory interactions between lymphoid and nonlymphoid cells in both EP and EN ewes. Immune DEGs were preferentially assigned to antigen presentation through the MHC complex, T-lymphocyte mediated immunity, and extracellular matrix interactions. Furthermore, the EP group showed altered regulation of cytokine and chemokine signalling and activation and recombination of B-cell receptors. Conversely, NF-kappa B signalling, apoptosis, and stress response were the main processes influenced in the EN group. In addition, statistically significant enrichment of the essential immune pathways of Binding and Uptake of Ligands by Scavenger Receptors in EP and p53 signalling in the EN group was found. In conclusion, this study provides new insights into the disease course and host‒pathogen interaction in naturally CL-infected sheep by investigating the blood transcriptome.

Background: Triple negative breast cancer (TNBC) is the most aggressive form of breast cancer with poor outcomes. The molecular basis of TNBC remains poorly understood. The objective of this study was to explore the relationship between obesity and TNBC in premenopausal and postmenopausal Caucasian women using whole genome transcription profiling. Methods: We compared gene expression levels of tumor samples drawn from normal weight, overweight and obese in pre and postmenopausal women diagnosed with TNBC. We performed hierarchical clustering to assess similarity in patterns of gene expression profiles, and conducted network and pathway analysis to identify molecular networks and biological pathways. Results: We discovered gene signatures distinguishing normal weight from obese, normal weight from overweight and overweight from obese individuals in both premenopausal and postmenopausal women. The analysis revealed molecular networks and biological pathways dysregulated in response to obesity. Among the discovered pathways included the unfolded protein response, endoplasmatic reticulum stress, B cell receptor and the autophagy signaling pathways in obese premenopausal women and the integrin, axonal guidance, ERK/MAPK and Glutathione biosynthesis signaling pathways obese postmenopausal women. Conclusions: The results suggest that both overweight and obesity are associated with TNBC, highlighting the need for conformation of these results in independent studies.

Deeper understanding of biological processes, disease mechanisms, and prospective therapeutic targets can be attained by Gene Expression Network Analysis (GENA), which offers a potent framework for revealing the intricate regulatory mechanisms controlling gene expression. GENA is a computational method used to understand the intricate interactions and relationships between the genes in a biological system. It entails using network theory, statistical analysis, and gene expression data to pinpoint functional modules, regulatory linkages, and important genes or pathways involved in a particular biological process or illness. This chapter broadly outlines the principles and practice of GENA, to an novice reader and outlines a simple method of performing a GENA using online Rice gene expression datasets available in various websites.

Long COVID, also known as post-acute sequelae of SARS-CoV-2 infection (PASC), has emerged as a significant health concern following the COVID-19 pandemic. Molecular mechanisms underlying the occurrence and progression of long COVID include viral persistence, immune dysregulation, endothelial dysfunction, and neurological involvement, and highlight the need for further research to develop targeted therapies for this condition. While a clearer picture of the clinical symptomatology is shaping, many molecular mechanisms are yet to be unraveled, given their complexity and high level of interaction with other metabolic pathways. This review summarizes some of the most important symptoms and associated molecular mechanisms that occur in long COVID, as well as the most relevant molecular techniques that can be used in understanding the viral pathogen, its affinity towards the host and the possible outcomes of host-pathogen interaction.

In the present study, we aimed to investigate whether an automated molecular diagnostic method based on PCR-reverse blot hybridization assay can discriminate between human Mycobacterium tuberculosis (MTB)-positive and -negative FFPE tissues and to compare the relative mRNA expression levels of various host immune markers between MTB-infected and uninfected human tissues using quantitative reverse transcription (qRT) PCR. A total of 52 human FFPE tissue samples from various regions of the body, including the lungs, lymph nodes, tendons, colon, and appendix, were collected and used for the molecular identification of Mycobacterium species and analysis of cytokine mRNA expression. As a result, IFN-γ, TNF-α, IP-10, CXCL9, CXCL11, and GM-CSF mRNA expression levels in MTB-infected tissues were significantly higher than those in uninfected samples. Additionally, the differences in the mRNA expression levels of IFN-γ, CXCL9, and GM-CSF between MTB-infected and uninfected tissues were statistically significant were statistically significant (p < 0.05). Correlation curve analysis indicated that the mRNA expression of IFN-γ was inversely proportional to that of IP-10 and that the mRNA expression levels of IFN-γ, TNF-α, CXCL9, CXCL11, GM-CSF, and TNFR were proportional and well-correlated. Furthermore, to establish marker profiles for detecting MTB infection, the statistically significant expression levels of three markers were combined. We confirmed that the combined profile of IFN-γ, CXCL9, and GM-CSF expression levels was statistically significant (P < 0.001). Although the mRNA expression patterns of host immune markers may vary according to MTB infection status, these patterns may be highly correlated and can be simultaneously used as an additional indicator for diagnosing TB.

Growth-regulating factor (GRF) is a plant-specific transcription factor family, which is involved in nearly all of the central developmental processes in plants. However, little is known about GRF family genes in cultivated sunflower. In this study, 17 GRF genes were identified and characterized from sunflower genome. Their gene structures, conserved motifs, chromosomal distributions and phylogenetic relationships were analyzed. The expression patterns of these genes were detected in various tissues of sunflower inbred line SK02R, which revealed that the 10 seed-specific GRF genes may play important roles during seed development in sunflower. Additionally, transcripts changes of the GRFs under two major abiotic stresses and phytohormones showed that most of the detected GRFs were reduced significantly by GA3 treatment, and other treatments(ABA, NaCl and PEG6000) differently regulated various sunflower growth-regulating factors at different time points. MiR396 target analysis indicated that there may exist a complicated homeostasis between miR396 and its targets GRF and WRKY transcription factor genes in cultivated sunflower. The phylogenetic and expression analyses of the GRF gene family in sunflower would be useful for further cloning and function exploration of the HaGRF genes.

Leucine-rich repeat receptor-like kinases (LRR-RLKs) have been reported to play important roles in plant growth, development and stress responses. However, no comprehensive analysis of this family has been performed in Gossypium, which are important economic crop and suffer various stresses in growth and development. Here we conducted a comprehensive analysis of LRR-RLK family in four Gossypium species (G. arboreum, G. barbadense, G. hirsutum and G. raimondii). A total of 1641 LRR-RLK genes were identified in the four Gossypium species involved in our study. Maximum-likelihood phylogenetic tree revealed that all the LRR-RLK genes were divided into 21 subgroups. Exon-intron organization structure of LRR-RLK genes kept relative conserved in subfamilies and between Arabidopsis and Gossypium. Subfamilies XI and XII were found dramatically expanded in Gossypium. Tandem duplication acted as an important mechanism in expansion of Gossypium LRR-RLK gene family. Function analysis suggested that plant hormone signaling and plant-pathogen interaction pathway were enriched in Gossypium LRR-RLK genes. Promoters analysis and expression profiles analysis revealed that Gossypium LRR-RLK genes were extensively regulated by TFs, phytohormone and various environmental stimuli, and play key roles in stress defense and diverse development processes. Our study provided valuable information for further function study of Gossypium LRR-RLK genes.

Solanesol is a noncyclic terpene alcohol composed of nine isoprene units and it mainly accumulates in solanaceous plants, especially tobacco (Nicotiana tabacum L.). Here, RNA-seq analyses of tobacco leaves, stems, and roots were used to identify solanesol biosynthesis genes. Six 1-deoxy-d-xylulose 5-phosphate synthase, two 1-deoxy-d-xylulose 5-phosphate reductoisomerase, two 2-C-methyl-d-erythritol 4-phosphate cytidylyltransferase, four 4-diphosphocytidyl-2-C-methyl-d-erythritol kinase, two 2-C-methyl-d-erythritol 2,4-cyclodiphosphate synthase, four 1-hydroxy-2-methyl-2-(E)-butenyl 4-diphosphate synthase, two 1-hydroxy-2-methyl-2-(E)-butenyl 4-diphosphate reductase, six isopentenyl diphosphate isomerase, and two solanesyl diphosphate synthase (SPS) genes were identified to be involved in solanesol biosynthesis. Furthermore, the two N. tabacum SPS (NtSPS1 and NtSPS2), which had two conserved aspartate-rich DDxxD domains, were highly homologous with SPS enzymes from other solanaceous plant species. In addition, the solanesol contents of three organs, and leaves from four growing stages, corresponded with the distribution of chlorophyll. Our findings provide a comprehensive evaluation of the correlation between the expression of different biosynthetic genes and the accumulation of solanesol in tobacco.

Several specialized terpenoids that are used as phytoalexins are produced in the plant of tobacco. They play a key role in plant defense against pests and pathogens. Here, we performed genome-wide identification, classification, and expression profiling of the terpene synthase (TPS) genes in the tobacco genome. A total of 186 TPS genes were identified, which is significantly higher than that in wild Nicotiana, Capsicum and Solanum species. Phylogenetic analysis of these TPS genes in tobacco and Solanaceous relatives as well as Arabidopsis showed that a few tobacco-exclusive TPS clades occurred recently through gene duplication events. These clades are likely to be responsible for a range of compounds that separate N. tabacum cultivars from its wild relatives. Transcriptome sequencing of two cultivars of N. tabacum at several developmental stages revealed that some TPS genes from the TPS-a and TPS-b lineages are constitutively expressed and a few duplicated TPS genes are only expressed post anthesis at the time of harvest. Additionally, some TPS genes underwent alternative splicing, further diversifying the gene products of the TPS genes. Our findings provide an important candidate genes for further understanding of the diversity and evolution as well as functionality of TPS genes in the plant of tobacco.

Primary central nervous system and brain tumors are one of the global burdens that are continuously increasing in cases and requiring more treatment options. Surgery has been the leading treatment option for tumors, however, the localization of the tumor and its infiltrating nature make this option challenging. DNA microarray expression profiles of different CNS tumors provide insight into potential biomarkers for identifying different tumor types and subtypes. Here, we utilized the differentially expressed genes common in four expression profiles, GSE66354, GSE68848, GSE74195, and GSE43290 to reconstruct the gene co-expression network. In this study, we were able to identify preserved cluster genes, hub genes, co-regulating transcription factors, miRNA families, and candidate repurposed drugs. Fourteen identified hub genes, which were, CACNA1A, DNM1, GABRA1, GRIA2, MAPT, SLC17A7, SNAP25, SNAP91, STXBP1, SYT1, COL1A1, COL6A2, FBN2, and FN1 appeared to play a role in tumor progression and may serve as drug targets. We also reported the DEGs in each tumor type, which were ATRT, EPN, PA, MED, PNET, MEN, ACM, ODG, and GBM. Five identified miRNA families, which were, let-7 family, mir-124 family, mir-1 family, mir-103 family, and mir-27 family described in the literature to have tumor-suppressing characteristics. Drug-gene and drug-transcription factor network revealed 32 candidate repurposed drugs and 13 were validated through connectivity map analysis. Here, two main repurposed drugs fit the regulatory requirement for gene interaction, these were, quercetin and vorinostat. Future validation in experimental studies may utilize the future use of the candidate repurposed drugs. Our study provides insights into drug repurposing prospects and understanding tumor expressions.

Calcineurin B-like proteins (CBL)-interacting protein kinases (CIPKs) by interacting with CBLs regulate developmental processes, hormone signalling transduction and mediate stress responses in plants. Although the genome of chickpea is available, information of CIPK gene family has been missing in chickpea. Here, a total of 22 CIPK encoding genes were identified in chickpea and characterized by in silico methods. We found a high structural conservation in chickpea CIPK family. Our analysis showed that chickpea CIPKs have evolved with dicots from common ancestors, and extensive gene duplication events have played an important role in evolution and expansion of CIPK family in chickpea. Most chickpea CIPK proteins localize in cytoplasm and nucleus. Promoter analysis revealed various cis-regulatory elements related to plant development, hormone signaling and abiotic stresses. Expression analysis indicated that CIPKs are significantly expressed in a spectrum of developmental stages, tissue/organs that hinted their important role in plant development. Several CIPK genes had specific and overlapping expressions in different abiotic stresses and seed development stages, suggesting the important role of CIPK family in abiotic stress signaling, and seed development in chickpea. Thus, this study provides the avenue for detailed functional characterization of CIPK family in chickpea and other legume crops.

Cervical cancer is associated with persistent infections by high-risk HPV types that may have nucleotide polymorphisms and, consequently, different oncogenic potentials. Therefore, the objective of this study was to evaluate the genetic variability and structural effects of E7 oncogene of HPV58 in cervical scraping samples from Brazilian women. The study was carried out with patients from hospitals in metropolitan area of Recife, PE, Brazil. The most frequent HPV type was HPV16, 18 and 58, respectively. Phylogenetic analysis showed that the isolates were classified as sublineages A2, C1 and D2. Two positively selected mutations were found in E7: 63G and 64T. The mutations G41R, G63D and T64A in E7 protein were predicted to reduce the stability of the protein structure. Regarding the interaction of the E7 variant of HPV58 with the signaling of the NF-kB pathway, we observed that the variant HPV58/UFPE-54S decreased the activity of the pathway when compared to the prototype and the other variants behaved similarly to the prototype and the prototype increased the activity of the pathway when compared to pcDNA. In this study, it was possible to identify mutations that may interfere in the molecular interaction between the viral oncoproteins and host proteins.

Schizophrenia (SCZ) is a severe mental disorder that may result in hallucinations, delusions, and extremely disordered thinking. How each cell type in the brain contributes to SCZ occurrence is still unclear. Here, we leveraged the human dorsolateral prefrontal cortex bulk RNA-seq data, then used the RNA-seq deconvolution algorithm CIBERSORTx to generate SCZ brain single-cell RNA-seq data for a comprehensive analysis to understand SCZ-associated brain cell types and gene expression changes. Firstly, we observed that the proportions of brain cell types in SCZ differed from normal samples. Among these cell types, astrocyte, pericyte, and PAX6 cells were found to have a higher proportion in SCZ patients (astrocyte: SCZ = 0.163, Control = 0.145, P.adj = 4.9×10-4; pericyte: SCZ = 0.057, Control = 0.066, P.adj = 1.1×10-4; PAX6 : SCZ = 0.014, Control = 0.011, P.adj = 0.014), while the L5/6_IT_CAR3 cells and LAMP5 cells are the exact opposite (L5/6_IT_Car3 : SCZ = 0.102, Control = 0.108, P.adj = 0.016; LAMP5 : SCZ = 0.057, Control = 0.066, P.adj = 2.2×10-6). Next, we investigated gene expression in cell types and functional pathways in SCZ. We observed chemical synaptic transmission dysregulation in two types of GABAergic neurons (PVALB and LAMP5), and immune reaction involvement in GABAergic neurons (SST) and non-neuronal cell types (endothelial and oligodendrocyte). Furthermore, we observed that some differential expression genes from bulk RNA-seq displayed cell-type-specific abnormal in the expression of molecules in SCZ. Finally, the cell types with the SCZ-related transcriptomic changes could be considered to belong to the same module since we observed two major similar coordinated transcriptomic changes across these cell types. Together, our results offer novel insights into cellular heterogeneity and the molecular mechanisms underlying SCZ.

Identification of the differentially-expressed genes is important for clarification of the complex molecular mechanisms under drought conditions. In this experiment, transcriptome profiles of sensitive and tolerant tomato genotypes under drought stress were analyzed. Three up-regulated genes were selected, included CAB3 (Chlorophyll a-b binding protein3), SAMDC (S-adenosylmethionine decarboxylase proenzyme), and ACS9 (1-aminocyclopropane-1-carboxylate synthase 9). After bioinformatics analysis, tomato genotypes were subjected to drought stress and gene expression was determined using Real-Time-PCR. Physiological parameters of genotypes were also measured by spectrophotometer-based methods. According to the results, these three genes play a key role in stress tolerance. Expression of the CAB3 gene in both sensitive and tolerant genotypes was not significantly different compared to the control, but the SAMDC gene decreased in both genotypes and the ACS9 gene decreased in sensitive genotype and increased in tolerant genotype. The physiological analysis also showed that under stress conditions, the photosynthetic system of the plant was disrupted and the chlorophyll content was reduced, but, proline content and antioxidant enzymes activity increased, in which their quantity in the tolerant genotype was significantly higher than sensitive. Under drought stress, due to damage to the lipid membrane, Malondialdehyde content also increased, in which the sensitive genotype was more affected.

Nitrogen (N) is an important mineral nutrient for plant growth, as well as a limiting factor for crop yield, and how to improve the utilization efficiency of N fertilizer by plants is currently a research hotspot. This article uses bioinformatics methods to identify and analyze members of the glutamine synthetase (GS) and glutamate synthase (GOGAT) gene families in pecan. A total of 6 GS genes and 4 GOGAT genes were identified, and their physicochemical properties, gene structures, and homologous evolutionary relationships were analyzed. Analysis of tissue-specific expression of GS and GOGAT genes based on transcriptome data from pecan. The enzyme activities of GS and GOGAT and the gene expression were quantitatively analyzed under different N form ratios in pecan. According to the results, the promoter cis-acting elements of GS and GOGAT genes can be roughly divided into three types: light-responsive elements, hormone-responsive elements, and stress-responsive elements. The results of homologous evolution showed that there was no tandem duplication event for the two gene families, and GS and GOGAT have undergone purification during the evolutionary process. CiGS2s and CiFd-GOGATs were expressed mainly in leaves, and CiNADH-GOGATs were expressed mainly in fruits. The qPCR analysis results showed that T4 treatment significantly increased the expression levels of CiGS and CiGOGAT genes in the leaves. The enzyme activities of GS and GOGAT in pecan were significantly increased under T3, T4, and T5 treatments. In summary, a higher proportion of ammonium nitrogen (NH4+) in the nutrient solution was profit to pecan NH4+ assimilation. This study determined the appropriate nitrogen ratio for pecan, promoting a theoretical basis for reducing environmental pollution caused by nitrogen fertilizer and improving the nitrogen utilization efficiency of pecan. In summary, both CiGSs and CiGOGATs exhibit tissue specificity, and an ammonium-nitrate mixture with a higher proportion of NH4+ is more favorable for NH4+ assimilation in pecan. This study provides a reference basis for further understanding the functions of CiGSs and CiGOGATs in pecan, and offers a theoretical foundation for improving N use efficiency in pecan.

Alfalfa (Medicago sativa L.) is the most cultivated forage legume around the world. Under a variety of growing conditions, forage yield in alfalfa is stymied by biotic and abiotic stresses including heat, salt, drought, and disease. Given the sessile nature of plants, they use strategies such as differential gene expression to respond to environmental cues. Transcription factors control the expression of genes that contribute to or enable tolerance and survival during periods of stress. Basic-leucine zipper (bZIP) transcription factors have been demonstrated to play a critical role in regulating plant growth and development as well as mediate the responses to abiotic stress in several species, including Arabidopsis thaliana, Oryza sativa, Lotus japonicus, and Medicago truncatula. However, there is little information about bZIP transcription factors in cultivated alfalfa. In the present study, 237 bZIP genes were identified in alfalfa from publicly available sequencing data. Multiple sequence alignments showed the presence of intact bZIP motifs in the identified sequences. Based on previous phylogenetic analyses in Arabidopsis thaliana, alfalfa bZIPs were similarly divided and fell into 10 groups. The physicochemical properties, motif analysis, and phylogenetic study of the alfalfa bZIPs revealed high specificity within groups. The differential expression of alfalfa bZIPs in a suite of tissues indicates that particular bZIP genes are specifically expressed at different developmental stages in alfalfa. Similarly, expression analysis in response to ABA, cold, drought, and salt stresses, indicates that a subset of bZIP genes are also differentially expressed and likely play a role in abiotic stress signaling and/or tolerance. These expression patterns were further verified by qRT-PCR. However, further functional characterization of bZIP transcription factors in alfalfa will help illuminate the role they play in stress tolerance mechanisms in legumes and facilitate the molecular breeding of stress tolerance in alfalfa.

Frontal Temporal Dementia (FTD) is a neurological disorder known to have less therapeutic options. So far only a few biomarkers are available for FTD that can be used as potential comorbidity targets. To this end, in the present study we aim to identify potential biomarkers or susceptible genes for FTD that show comorbidities with diseases such as COVID-19 and Breast Cancer. Gene expression Omnibus datasets containing FTD expression profiles from African American and white ethnicity background were included in our study. In FTD samples of GSE193391 dataset, we identified 305 DEGs with 168 genes being up-regulated and 137 genes being down-regulated. We conducted comorbidity analysis for COVID-19 and Breast cancer followed by analysis of potential drug interactions, pathogenicity, analysis of genetic variants and functional enrichment analysis. Resulting genes AKT3, GFAP, ADCYAP1R1, VDAC1, and C4A showed significant transcriptomic alterations in FTD and significant comorbidity status with COVID-19 and Breast cancer. Functional pathway analysis revealed that these comorbid genes were significantly enriched in the pathways such as glioma, JAK/STAT signaling, systematic lupus erythematosus, neurodegeneration-multiple diseases and neuroactive ligand-receptor interaction. Overall, from these results, we concluded that these genes could be recommended as potential therapeutic targets for the treatment of comorbidities in FTD.

The recent outbreak of Coronavirus Disease 2019 (COVID-19) is a major threat to human health and the global economy. In addition to the development of vaccines there is an urgent need for preventive and therapeutic strategies towards severe COVID-19. Key to this would be a better understanding of the molecular mechanismsoiuuuu affected by SARS-COV-2. To address this, we performed a systems biology approach by integrating available RNA-seq datasets from post-mortem lung tissue of COVID-19 patients and cell culture models infected with SARS-COV-2, Respiratory Syncytial virus or influenza virus. We identified two gene-expression modules that are commonly regulated by the three viral diseases and one that is specific to COVID-19. All 3 gene-expression modules represent key inflammatory processes. We identified several proteins within these networks that can be targeted by FDA approved drugs. Key examples are TNF, NFkB, INTERLEUKIN-1 and ALOX5 signaling pathways. Our data also suggest that Vitamin D supplementation and a ketogenic diet should be further analyzed as preventive strategies. In conclusion, our data highlights the potential of transcriptomics to unravel the pathological processes related to COVID-19 and guide the initiation of clinical trials.

G protein-coupled receptor 15 (GPR15, also known as BOB) is an extensively studied orphan GPCR involving HIV infection, colonic inflammation and smoking-related diseases. Recently, GPR15 was deorphanized and its corresponding natural ligand demonstrated an ability of inhibiting cancer cell growth. However, no study reported the potential role of GPR15 in a pan-cancer manner. Using large-scale publicly available data from the Cancer Genome Atlas (TCGA) and the Genotype-Tissue Expression (GTEx) databases, we found that GPR15 expression is significantly lower in colon adenocarcinoma (COAD) than in normal tissues. And among 33 cancer types, GPR15 expression is significantly correlated with the prognoses of COAD, neck squamous carcinoma (HNSC), lung adenocarcinoma (LUAD) positively and stomach adenocarcinoma (STAD) negatively. This study also revealed that commonly upregulated gene set in the high GPR15 expression group (stratified via median) of COAD, HNSC, LUAD and STAD are enriched in immune systems, indicating that GPR15 might be considered as a potential target for cancer immunotherapy. Furthermore, we modelled the 3D structure of GPR15 and conducted the structure-based virtual screening. The top 8 hits compounds were screened and then subjected to molecular dynamics (MD) simulation for stability analysis. Our study provided novel insights into the role of GPR15 in a pan-cancer manner and discovered a potential hit compound for GPR15 antagonists.

The Translationally Controlled Tumor Protein (TCTP) is a multifunctional protein found in most eukaryotic taxa. In plants, TCTP regulates cell-cycle progression as well as programmed cell death and is involved in plant regeneration and responses to abiotic and biotic stress, among other functions. Additionally, there is evidence that the transcript and protein are transported over long-distance. Potatoes are a crop of agricultural interest that can reproduce asexually by regeneration from stolons and tubers. Therefore, in this study the effect of transient expression of Solanum tuberosum TCTP (StTCTP) on tuberization was analyzed. StTCTP mRNA was shown to be transported long-distance. Additionally, transient overexpression of StTCTP resulted in tuber bearing shoots with a greater diameter compared to control plants. Furthermore, the early stages of tuberization were induced compared to control plants, in which only mature tubers were observed.

Recent studies suggest that PEBP1/RKIP and YY1, despite having distinct molecular functions, may interact and mutually influence one another's activity. They exhibit reciprocal control over each other's expression through regulatory loops, prompting the hypothesis that their interplay could be pivotal in cancer advancement and resistance to drugs. To delve into this interplay's functional characteristics, we conducted a comprehensive analysis using bioinformatics tools across a range of cancers. Our results confirm the association between elevated YY1 mRNA levels and varying survival outcomes in diverse tumors. Furthermore, we observed differing degrees of inhibitory or stimulatory effects of these two genes in various cancer pathways, along with correlations between their mRNA expression and immune infiltration. Additionally, YY1/PEBP1 expression (or methylation) displayed connections with genomic alterations across cancer types. Notably, we uncovered links between YY1/PEBP1 and different indicators of immunosuppression, such as immune checkpoint blockade response, and T-cell dysfunction/exclusion levels, across different patient groups. Overall, our findings underscore the significant role of the interplay between YY1 and PEBP1 in cancer progression, influencing genomic changes, tumor immunity or the tumor microenvironment. Additionally, these two gene products appear to impact the sensitivity of anticancer drugs, opening new avenues for cancer therapy.

Strigolactones (SLs) are a novel emerging plant hormones, which play important roles in regulating plant organ development and environmental stress tolerance. Even though the SL related genes have been identified and well characterized in some plants. The information of SL related genes in soybean is not fully established yet, especially in response to salt stress. In this study, we identified nine SL biosynthesis genes: two D27, two CCD7, two CCD8, and three MAX1, and seven SL signaling genes: two D14, two MAX2 and three D53 in soybean genome. We found that SL biosynthesis and signaling genes are conserved during evolution in different species. Syntenic analysis of these genes revealed their location on nine chromosomes as well as existence of ten pairs of duplication genes. Moreover, plant hormone and stress-responsive elements were identified in the promoter regions of SL biosynthesis and signaling genes. By using quantitative real-time PCR (qRT-PCR), we confirmed that SL genes have different tissue expression in roots, stems and leaves. Further, we also explored the expression profiles of SL biosynthesis and signaling genes under salt stress. These results suggested that SL signaling genes may play important regulatory roles in response to salt stress. In conclusion, we identified and provided valuable information on the soybean SL biosynthesis and signaling genes, and established a foundation for further functional analysis of soybean SL related genes in response to salt stress.

The basic (region) leucine zippers (bZIPs) are evolutionarily conserved transcription factors widely distributed in eukaryotic organisms. In plants, they are not only involved in growth and development, defense and stress responses, and regulation of physiological processes, but also play a pivotal role in regulating secondary metabolism. To explore the function related to the bZIP gene family in Stevia rebaudiana Bertoni, we identified 105 SrbZIP genes at the genome-wide level and classified into 12 subfamilies using bioinformation methods. Analysis of physicochemical properties showed that the amino acid length readings ranged from 96 aa (SrbZIP13) to 704 aa (SrbZIP17), and the isoelectric point ranged from 4.64 (SrbZIP21) to 10.08 (SrbZIP101). The majority of these proteins (94.29%) were localized in the nucleus. Three main classes of cis-acting elements found in the SrbZIP promoter regions, including development-related elements (CAT-box), defense and stress-responsive elements (LTR/ARE/MBS), and phytohormone respon- sive elements (ABRE/TGACG-mitif/TCA-element). Through protein-protein interac- tion network of 105 SrbZIP proteins, SrbZIP proteins were mainly classified into four major categories: ABF2/ABF4/ABI5 (SrbZIP51/SrbZIP38/SrbZIP7), involved in phy- tohormone signaling, GBF1/GBF3/GBF4 (SrbZIP29/SrbZIP63/SrbZIP60) involved in environmental signaling, AREB3 (SrbZIP88), PAN (SrbZIP12), TGA1 (SrbZIP69), TGA4 (SrbZIP82), TGA7 (SrbZIP31), TGA9 (SrbZIP95), TGA10 (SrbZIP79) and HY5 (SrbZIP96) involved in cryptochrome signaling, and FD (SrbZIP23, SrbZIP24, SrbZIP47, SrbZIP67, SrbZIP68, SrbZIP71, SrbZIP72) promoted flowering. The transcriptomic data showed that SrbZIP genes differentially expressed in six S. rebaudiana cultivars (‘023’, ‘110’, ‘1188’, ‘11-14’, ‘GP’ and ‘GX’). Moreover, the expression levels of selected 15 SrbZIP genes in response to light, abiotic stress (low temperature, salt and drought), phytohormones (methyl jasmonate, gibberellic acid and salicylic acid) treatment and in different tissues were analyzed utilizing qRT-PCR. Six SrbZIP (SrbZIP54, SrbZIP63, SrbZIP32, SrbZIP45, SrbZIP60 and SrbZIP9) genes were further identified to be highly induced by factors affecting glycoside synthesis. Among them, three SrbZIP genes (SrbZIP54, SrbZIP63 and SrbZIP32) were predicted to be related to stress-responsive terpenoid synthesis in S. rebaudiana. The protein-protein interaction network expanded the potential functions of SrbZIP genes. This study firstly provided the comprehensive genome-wide report of the SrbZIP gene family, laying a foundation for further research on evolution, function and regulatory role of bZIP gene family in terpenoid synthesis in S. rebaudiana.

The basic helix-loop-helix (bHLH) transcription factor family is the second-largest transcription factor in plants. Members of this family are involved in the process of growth and development, secondary metabolic biosynthesis, signal transduction, and plant resistance. Loropetalum chinense var. rubrum is a critical woody plant with higher ornamental and economic values, which used as ornamental architecture and traditional Chinese herbal medicine plant. However, the bHLH transcription factor in the L. chinense var. rubrum (Loropetalum chinense var. Rubrum) have yet to be systematically demonstrated, and their role in anthocyanin biosynthesis remains secret. Here, we identified 165 potential LcbHLHs genes by two methods, and they were unequally distributed on chromosomes 1 to 12 of the L. chinense genome. Based on a phylogenetic comparison with proteins from Arabidopsis, these bHLH proteins were categorized into 21 subfamilies. Most LcbHLHs in a given subfamily had similar gene structures and conserved motifs. The gene ontology annotation and cis-elements predicted that the LcbHLHs had many molecular functions and were involved in plant growth processes, including flavonoid and anthocyanin biosynthetic processes, etc. Transcriptomic analysis revealed different expression patterns among different tissues and cultivars in L. chinense. Many LcbHLHs were expressed in leaves, and only a few genes were highly expressed in flowers. Six LcbHLH genes differentially expressed between species and periods in color variation, which may play a key role in anthocyanin synthesis. These further functional predictions of LcbHLHs were established by quantitative real-time PCR analysis and protein-protein interaction network. This study identified the six potential bHLH genes related to anthocyanin in L. chinense and created a solid framework for future research on the function and evolution of bHLHs.

Plant Cysteine Oxidase (PCO) is a plant O2-sensing enzyme of catalyzing oxidation of cysteine to Cys-sulfinic acid at the N-termini of target proteins. To better understand the Brassica napus PCO gene family, PCO genes in B. napus and related species were analyzed. In this study, 20, 7 and 8 PCO genes were identified in Brassica napus, Brassica rapa and Brassica oleracea, respectively. According to phylogenetic analysis, the PCOs were divided into five groups: PCO1, PCO2, PCO3, PCO4 and PCO5. Gene organization and motif distribution analysis suggested that PCO gene family was relatively conserved during evolution. Expression analysis showed that PCO genes were expressed in different tissues at different developmental stages. Interestingly, most of the Bna/Br/BoPCO5 members were expressed in leaves, roots, flowers, buds and siliques, suggesting an important role in both vegetative and reproductive development. Expression of BnaPCO was induced by various abiotic stress, especially waterlogging stress, which was consistent with the result of cis-element analysis. In this study, the PCO gene family of Brassicaceae was analyzed for the first time, which contributes to a comprehensive understanding of the origin and evolution of PCO genes in Brassicaceae and the function of BnaPCO in regulating plant responses to abiotic stresses.

Abstract: The baculovirus expression vector (BEV) system is an efficient, cost effective and scalable method to produce recombinant adeno-associated virus (rAAV) gene therapy vectors. Most BEV designs emulate the wild type AAV transcriptome and translate the AAV capsid proteins, VP1, VP2 and VP3 from a single mRNA transcript with three overlapping open reading frames (ORFs). Non-canonical translation initiation codons for VP1 and VP2 reduce their abundances relative to VP3. Changing capsid ratios to improve rAAV vector efficacy requires theoretical modification of translational context. We have made a Lac repressor inducible system to empirically regulate expression of VP1 and VP2 proteins relative to VP3 in the context of the BEV. We demonstrate the use of this system to tune the abundance, titer and potency of a neurospecific rAAV9 serotype derivative. VP1:VP2:VP3 ratios of 1:1:8 gave optimal potency for this rAAV. It was discovered that ratios of capsid proteins expressed were different than ratios that ultimately were in purified capsids. Over expressed VP1 did not become incorporated into capsids while over expressed VP2 did. Overabundance of VP2 correlated with reduced rAAV titers. This work demonstrates a novel technology for controlling production of rAAV in the BEV system and shows a new perspective on the biology of rAAV capsid assembly.

The CDPK family genes play crucial roles in signal transduction pathways during plant develop-ment and stress response. In this study, we comprehensively analyzed the CDPK family genes in mustard (Brassica juncea L.), resulting in the identification of 101 genes (BjuCDPK1–101) located on chromosomes AA_Chr01 to BB_Chr08. RNA-seq analysis showed that most BjuCDPK genes were predominantly expressed in the root and flower bud, suggesting their organ-specific expression patterns. Furthermore, the expression levels of most BjuCDPK genes were significantly altered by cold stress. The promoter regions of BjuCDPK genes had abundant cis-acting elements related to hormone and stress. The interaction network analysis predicted that BjuCDPK proteins respond to abiotic stress probably via interaction with ABF1, ABF4, and/or DI19. Our findings provide valu-able information for further investigation of the cold stress adaption of mustard via the CDPK signaling pathway.

During the response of plants to adverse stresses, aquaporin (AQP) plays a prominent role in membrane water transport based on received upstream signals. In addition, they have various physical parts for dealing with environmental stresses. Due to the importance of the AQP gene family, studies have been conducted investigating the function and regulatory system of these genes. However, many of their molecular aspects are still unknown. This study aims to carry out a genomic-wide investigation of the AQP gene family in durum wheat using bioinformatics tools and to investigate the expression patterns of some members in response to salt stress. Our results showed that there are 80 TtAQP genes in durum wheat, which are classified into four main groups based on phylogenetic analysis. Many duplications were observed between the members of the TtAQP gene family, and high diversity in response to post-translational modifications was observed between TtAQP family members. The expression pattern of TtAQP genes disclosed that these genes are primarily upregulated in response to salt stress. Besides, qPCR data revealed that TtAQPs are more induced in delayed responses to salinity stress. Overall, our findings illustrate that TtAQP members are diverse in terms of their structure, regulatory systems and expression levels.

Diacylglycerol kinase (DGK) is a lipid kinase that phosphorylates diacylglycerol (DAG) to generate phosphatidic acid (PA). Based on converting one important signaling molecule (DAG) to another (PA), DGK plays an important role in plant responses to abiotic stress including waterlogging stress. However, no studies have been reported on the characterization of the DGK gene family in the waterlogging-tolerant kiwifruit germplasm, Actinidia valvata Dunn. Here, we identified 18 AvDGK genes in the A. valvata genome. The phylogenetic analysis showed that AvDGKs can be classified into three clusters, and members within the same cluster have similar domain distribution, exon-intron structures and conserved motif compositions. Chromosome localization analysis revealed that all the AvDGK genes are located across 18 different chromosomes. There were 29 duplicated gene pairs in kiwifruit and all had undergone purifying selection during evolution. Promoter cis-element analysis revealed that the cis-elements within AvDGK genes are associated with multiple functions, including phytohormone signal transduction, stress response, and plant growth and development. The expression pattern analyses indicated that AvDGK play an important role in the fruit development and plant response to waterlogging stress. AvDGKs gene family in the tetraploid A. valvata genome might promote PA synthesis and subsequent signal transduction both under short-term and long-term waterlogging stress. These results provide information regarding the structural characteristics and potential function of DGK genes within kiwifruit and lay a fundamental basis for further research into breeding for enhancing the kiwifruit’s tolerance to waterlogging stress.

As a major infectious disease in chickens, Newcastle disease causes considerable economic losses in the poultry industry, especially in developing countries where there is limited access to effective vaccination. Therefore, enhancing resistance to the virus in commercial chickens through breeding is a promising way to promote poultry production. In this study, we investigated gene expression changes at 2 and 6 dpi after infection at day21 with a lentogenic Newcastle disease virus in a commercial egg-laying chicken hybrid using RNA sequencing analysis. By comparing NDV challenged and nonchallenged groups, 526 differentially expressed genes (DEGs) (FDR &lt; 0.05) were identified at 2 dpi, and only 36 at 6 dpi. For the DEGs at 2 dpi, IPA analysis predicted inhibition of multiple signaling pathways in response to NDV that regulate immune cell development and activity, neurogenesis and angiogenesis. Upregulation of Interferon Induced Protein with Tetratricopeptide Repeats 5 (IFIT5) in response to NDV was consistent between the current and most previous studies. Sprouty RTK Signaling Antagonist 1 (SPRY1), a DEG in the current study is located in a significant QTL associated with virus load at 6 dpi in the same population. These identified pathways and DEGs provide potential targets to further study breeding strategy to enhance NDV resistance in chickens.

A plant’s early response to pathogen stress is a vital indicator of its disease resistance. In order to study the response mechanism of resistant and susceptible flax cultivars to Fusarium oxysporum f. sp. lini (Foln), we applied RNA-sequencing to analyze transcriptomes of flax with Foln 0.5, 2 and 8 hours post inoculation (hpi). We found a significant difference in the number of differential expression genes (DEGs) between resistant and susceptible flax clutivars. The number of DEGs in the Fusarium-resistant cultivar increased dramatically at 2 hpi, and a large number of DEGs participated in the Fusarium-susceptible cultivar response to Foln infection 0.5 hpi. GO enrichment analysis determined that the up-regulated DEGs of both flax cultivars were enriched such as oxidoreductase activity and oxidation-reduction process. At the same time, the genes involved in diterpenoid synthesis were up-regulated in resistant cultivar, while those involved in extracellular region, cell wall and organophosphate ester transport were down-regulated in susceptible cultivar. KEGG enrichment analysis showed the genes encoded WRKY 22 and WRKY33 which involved in MAPK signaling pathway were up-regulated expressed in S-29 and down-regulated expressed in R-7, negatively regulated the disease resistance of flax; The genes encoded Hsp 90 family which in involved in plant pathogen interaction pathway were up-regulated in R-7 and down-regulated in S-29, which positively regulated the disease resistance of flax; The genes encoded MYC2 transcription factor and TIFY proteins which involved in plant hormone signaling pathway were up-regulated in R-7, and regulated the jasmonic acid metabolism of flax and the signal transduction of plant hormones. Meanwhile seven regulatory genes with the most correlation were screened out, Among Lus10025000.g and Lus10026447.g regulated other genes expressed both in plant hormone signal transduction pathway and MAPK signal pathway. In conclusion, these findings will facilitate further studies on the function of these candidate genes in flax of response to Fusarium stress, and the breeding of disease-resistant flax cultivar.

Deep learning is based on matrix computing with a large amount of hidden parameters that is not visible outside the computing module. Deep learning is a non-linear system and a linear approach is not possible. It is natural for people to visualize the algorithm and to follow some hidden parameters. In this paper, we propose a simple graphical programming of a nonlinear system based on drawing the simplest unit, such as a single neuron model. A more complex scheme of a classical neural network is obtained by commands copy-move-place of one neuron. The number of neurons and layers can be chosen arbitrarily. Once the scheme is complete, the implementation code is evaluated with symbolic parameters and nonlinear activation functions. This cannot be done manually. With the symbolic expression of outputs in terms of inputs and symbolic parameters, including symbolic activation pure functions, some other properties can be derived in closed form. This unique original approach can help scientists and developers and design numerical algorithms for machine learning and to understand how deep learning algorithms work.

To guide analysts to select the right tool and parameters in differential gene expression analysis of single-cell RNA sequencing (scRNA-seq) data, we developed a novel simulator that recapitulates the data characteristics of real scRNA-seq datasets while accounting for all the relevant sources of variation in a multi-subject, multi-condition scRNA-seq experiment: the cell-to-cell variation within a subject, the variation across subjects, the variability across cell types, the mean/variance relationship of gene expression across genes, library size effects, group effects, and covariate effects. By applying it to benchmark 12 differential gene expression analysis methods (including cell-level and pseudo-bulk methods) on simulated multi-condition, multi-subject data of the 10x Genomics platform, we demonstrated that methods originating from the negative binomial mixed model such as glmmTMB and NEBULA-HL outperformed other methods. Utilizing NEBULA-HL in a statistical analysis pipeline (https://github.com/interactivereport/scRNAseq_DE) for single cell analysis will enable scientists to better understand cell-type specific transcriptomic response to disease or treatment effects and to discover new drug targets. Further, application to two real datasets showed the outperformance of our differential expression (DE) pipeline, with unified findings of differentially expressed genes (DEG) and a pseudo-time trajectory transcriptomic result. In the end, we made recommendations of filtering strategies of cells and genes based on simulation results to achieve optimal experimental goals.

The corneal epithelium (CE) is spread between two domains, the outer vascularized limbus and the avascular cornea proper. Epithelial cells undergo constant migration from the limbus to the vision-critical central cornea. Coordinated with this migration the cells undergo differentiation changes where a pool of unique stem/precursor cells at the limbus yields the mature cells that reach the corneal center. Differentiation is heralded by the expression of the corneal-specific Krt12. Processing data acquired by scRNA-Seq we found that the increase of Krt12 expression occurs in four distinct steps within the limbus plus a single continuous increase in the cornea. Differential gene analysis demonstrated that these domains reflect discreet stages of CE differentiation and yielded extensive information of the genes undergoing down or up regulation in the sequential transition from less to more differentiate conditions. The approach allowed the identification of multiple gene cohorts, including a) the genes that have maximal expression in the most primitive, Krt12-negative cell cohort, which is likely to include the stem/precursor cells; b) the sets of genes that undergo continuous increase or decrease along the whole differentiation path, and c) and the genes showing maximal positive or negative correlation with the changes in Krt12.

This research elucidates the dynamic expression of expansin genes during the wheat grain (Triticum aestivum L.) development process using comprehensive meta-analysis and experimental validation. We leveraged RNA-seq data from multiple public databases, applying stringent criteria for selection, and identified 60,852 differentially expressed genes across developmental stages. From this pool, 28,558 DEGs were found to exhibit significant temporal regulation in at least two different datasets, and were enriched for processes integral to grain development such as carbohydrate metabolism and cell wall organization. Notably, 30% of the 241 known expansin genes showed differential expression during grain growth. Hierarchical clustering and expression level analysis revealed temporal regulation and distinct contributions of expansin subfamilies during the early stages of grain development. Further analysis using co-expression networks underscored the significance of expansin genes, revealing their substantial co-expression with genes involved in cell wall modification. Finally, qPCR validation and grain morphological analysis under field conditions indicated a significant negative correlation between expression of select expansin genes and grain size and weight. This study illuminates the potential role of expansin genes in wheat grain development and provides new avenues for targeted genetic improvements in wheat.

The largest gene families in plants were found to be Glutathione transferase (GSTs), which played significant roles in regulating plant growth, development, and stress response. Within the GSTs gene family, members were found to play a crucial role in the low-temperature response process of plants. A comprehensive study identified a total of 70 BraGSTs genes. Cluster analysis results demonstrated that the BraGSTs in Brassica rapa(B.rapa) could be categorized into eight sub-families and were unevenly distributed across ten chromosomes. The 39 BraGSTs genes were found to be organized into 15 tandem gene clusters, with the promoters containing multiple cis-elements associated with low-temperature response. Cold stress was observed to stimulate the expression of 15 genes, with the BraGSTF2 gene exhibiting the highest level of expression, sug-gesting its significant involvement in winter B.rapa's response to low-temperature stress. Sub-cellular localization analysis of the BraGSTF2 protein indicated its potential expression in both the cell membrane and nucleus. The analysis of stress resistance in BraGSTF2 transgenic Arabidopsis thaliana lines emon-strated that the over-expression of this gene resulted in significantly elevated levels of SOD, POD activity, and SP content compared to the wild type following exposure to low temperatures. These levels reached their peak after 24 hours of treatment. Conversely, the MDA content was lower in the transgenic plants compared to the wild type. Additionally, the survival rate of BraGSTF2 transgenic Arabidopsis thaliana was higher than that of the wild type, suggesting that the BraG-STF2 gene may play a crucial role in enhancing the cold stress tolerance of winter B. rapa. This study lay a foundation for further study on the role of BraGSTs gene in the molecular regulation of cold resistance in winter B.rapa.

Barth Syndrome, a rare X-linked disorder affecting 1:300,000 live births, results from defects in Tafazzin, an acyltransferase that remodels cardiolipin and is essential for mitochondrial respiration. Barth Syndrome patients develop cardiomyopathy, muscular hypotonia and cyclic neutropenia during childhood, rarely surviving to middle age. At present, no effective therapy exists and downstream transcriptional effects of Tafazzin dysfunction are incompletely understood. To identify novel, cell-specific, pathological pathways that mediate heart dysfunction, we performed single-nucleus RNA-sequencing (snRNA-seq) on wild-type (WT) and Tafazzin-knockout (Taz-KO) mouse hearts. We determined differentially expressed genes (DEGs) and inferred predicted cell-cell communication networks from these data. Surprisingly, DEGs were distributed heterogeneously across the cell types, with fibroblasts, cardiomyocytes, endothelial cells, macrophages, adipocytes and pericytes exhibiting the greatest number of DEGs between genotypes. One differentially expressed gene was detected for the lymphatic endothelial and mesothelial cell types, while no significant DEGs were found in the lymphocytes. A Gene Ontology (GO) analysis of these DEGs showed cell-specific effects on biological processes such as fatty acid metabolism in adipocytes and cardiomyocytes, increased translation in cardiomyocytes, endothelial cells and fibroblasts, in addition to other cell-specific processes. Analysis of ligand-receptor pair expression, to infer intercellular communication patterns, revealed the strongest dysregulated communication involved adipocytes and cardiomyocytes. For the knockout hearts, there was a strong loss of ligand-receptor pair expression involving adipocytes, and cardiomyocyte expression of ligand-receptor pairs underwent reorganization. These findings suggest that adipocyte and cardiomyocyte mitochondria may be most sensitive to mitochondrial Tafazzin deficiency and that rescuing adipocyte mitochondrial dysfunction, in addition to cardiomyocyte mitochondrial dysfunction, may provide therapeutic benefit in Barth Syndrome patients.

Cotton (Gossypium spp.) is a major source of natural fiber and an important cash crop. The cotton growth habit and architecture determine its productivity and influence management strategies for commercial production. The GATA transcription factors (TFs) control various developmental processes in plants, such as flower, bract and embryo development, and petal differentiation. As stable transformation is still a bottleneck in many plant species, TRV-VIGS was used to manipulate gene expression in different plants, including Gossypium hirsutum L. In this study, we undertook the TRV-based VIGS to functionally characterize two candidate genes, Gohir.D05G103700 and Gohir.D12G153600, identified through the expression QTL analysis for five floral induction and meristem identity genes using the upland cotton mini-core collection. Virus-induced silencing of the Gohir.D05G103700 gene resulted in up to a 1.4-fold reduction in the transcript level in two inoculated plants, G3 and G4, and Gohir.D12G153600 gene resulted in up to a 2.3-fold reduction in transcript level in a single inoculated plant P05 relative to the mock-treated plant. The TRV2-Gohir.D05G103700 inoculated plants precisely G3 and G4 also exhibited loss of the supernumerary (fourth) floral bract in the squares, whereas the TRV2-Gohir.D12G153600 inoculated plants did not show any observable phenotypic change relative to the mock-treated plants.  Altogether, this study suggested that TRV-VIGS can be used to characterize genes in cotton relatively rapidly and the cotton Gohir.D05G103700 gene is a positive regulator of the indeterminate growth habit in cotton, which could be manipulated to obtain a cotton plant with architecture best suited for the cultivation area.

The present study was performed to examine that transforming growth factor beta (TGF-β) in root-surrounding tissues on deciduous teeth during the physiological root resorption regulates the differentiation induction into odontoclast. We prepared root-surrounding tissues with (R) or without (N) physiological root resorption scraped off at three regions (R1-R3 or N1-N3) from the cervical area to the apical area of the tooth and measured both TGF-β and the tartrate-resistant acid phosphatase (TRAP) activities. The TGF-β activity level was increased in N1-N3, whereas the TRAP activity was increased in R2 and R3. In vitro experiments for RANKL-mediated osteoclast differentiation revealed that TGF-β in N1-N3 and R1-R3 enhanced the TRAP activity in RAW264 cells. A genetic study indicated that the mRNA level of TGF-β1 in N1 and N2 was significantly increased, and corresponded with that of osteoprotegerin (OPG). In contrast, the expression level of receptor activator of NF-κB ligand (RANKL) was increased in R2 and R3. Our findings suggest that TGF-β is closely related to the regulation of OPG induction and RANKL-mediated odontoclast differentiation depending on the timing of RANKL and OPG mRNA expression in the root-surrounding tissues of deciduous teeth during physiological root resorption.

Hypertrophic Cardiomyopathy (HCM) is a common inherited disorder that can lead to heart failure and sudden cardiac death, characterized at the histological level by focal areas of myocyte disarray, hypertrophy and fibrosis, and few disease-targeted therapies exist. To identify, focal, spatially restricted alterations in transcriptional pathways and reveal novel therapeutic targets, we performed a spatial transcriptomic analysis of areas of focal myocyte disarray compared to areas of normal tissue, using a commercially available platform (GeoMx, nanoString). We analyzed surgical myectomy tissue from four patients with HCM and control interventricular septum tissue from two unused organ donor hearts that were free of cardiovascular disease. Histological sections were reviewed by an expert pathologist and 72 focal areas with varying degrees of myocyte disarray (normal, mild, moderate, severe) were chosen for analysis. Areas of interest were interrogated with the Human Cancer Transcriptome Atlas designed to profile 1800 transcripts. Differential expression analysis revealed significant changes in gene expression between HCM and Control tissue, and functional enrichment analysis indicated these genes were primarily involved in interferon production and mitochondrial energetics. Within HCM tissue, differentially expressed genes between areas of mild and moderate disarray were enriched for genes related to mitochondrial energetics (moderate disarray) and response to oxygen/cytokine levels (mild disarray). The comparison between areas of moderate and severe disarray were enriched for genes related to the c-Jun N-terminal kinase (JNK) cascade in severe disarray. Analysis of ligand-receptor pair gene expression revealed that HCM tissue exhibited downregulation of platelet-derived growth factor (PDGF), NOTCH, junctional adhesion molecule, and CD46 signaling, while showing upregulation of fibronectin, CD99, cadherin, and amyloid precursor protein signaling. A deconvolution analysis utilizing the matched single nuclei RNA-sequencing (snRNA-seq) data to determine cell type composition in areas of interest revealed significant differences in fibroblast and vascular cell composition in areas of severe disarray when compared to normal areas in HCM samples. Cell composition in normal areas from control tissue was also divergent from normal areas in HCM samples, which was consistent with the differential expression results. Overall, our data identify novel and potential disease-modifying targets for therapy in HCM.

Pea (Pisum sativum L.), like most legumes, forms mutualistic symbioses with nodule bacteria and arbuscular mycorrhizal (AM) fungi. The positive effect of inoculation is partially determined by the plant genotype, thus, pea varieties with high and low symbiotic responsivity have been described, but the molecular genetic basis of this trait remains unknown. Here, we compare the symbiotically responsive breeding line ‘Triumph’ of grain pea with its parental cultivars ‘Vendevil’ (a donor of high symbiotic responsivity) and ‘Classic’ (a donor of agriculturally valuable traits) using genome and transcriptome sequencing. We show that ‘Triumph’ inherited one-fourth of its genome from ‘Vendevil’, including the genes related to AM and nodule formation, and reveal that under combined inoculation with nodule bacteria and AM fungi, ‘Triumph’ and ‘Vendevil’, in contrast to ‘Classic’, demonstrate similar up-regulation of the genes related to solute transport, hormonal regulation and flavonoid biosynthesis in their roots. We also identify the gene PsGLP2, whose expression pattern distinguishing ‘Triumph’ and ‘Vendevil’ from ‘Classic’ correlates with difference within the promoter region sequence, making it a promising marker for the symbiotic responsivity trait. The results of the study may be helpful for future molecular breeding programs aimed at creation of symbiotically responsive cultivars of pea.

Expression analysis of small noncoding RNA (sRNA), including microRNA, piwi-interacting RNA, small rRNA-derived RNA, and tRNA-derived small RNA, is a novel and quickly developing field. Despite a range of proposed approaches, selecting and adapting a particular pipeline for transcriptomic analysis of sRNA remains a challenge. This paper focuses on the identification of the optimal pipeline configurations for each step of human sRNA analysis, including reads trimming, filtering, mapping, transcript abundance quantification and differential expression analysis. Based on our study, we suggest the following parameters for analysis of human sRNA in relation to categorical analyses with two groups of biosamples: (1) trimming with the lower length bound = 15 and the upper length bound = \(Read\ length - 40\% Adapter\ length\); (2) mapping on a reference genome with bowtie aligner with one mismatch allowed (-v 1 parameter); (3) filtering by mean threshold > 5; and (4) analyzing differential expression with DESeq2 with adjusted p-value < 0.05 or limma with p-value < 0.05 if there is very little signal and few transcripts.

The eyes are the sense organs through we view the world around us, and the cornea is the transparent layer which covers the outer, visible part of the eye. It is known that the gene expression of corneal endothelium depends on age. And the expressed genes of the endothelium of cornea for pediatric samples is different than that of the adult samples. The purpose or objective of this study was to characterize human corneal endothelial cell (HCEnC) gene expression and differential gene expression and to detect the expressed genes mapped to chromosomal loci associated with some corneal endothelial dystrophies. Both upregulated and downregulated genes were analyzed. For this purpose, total RNA was isolated from ex vivo corneal endothelium taken from six pediatric and five adult corneas. The complementary DNA was hybridized to the Affymetrix GeneChip 1.1 ST array. The data analysis was performed using the Enrichr software for both upregulated and downregulated genes. These are described in more detail in the results and discussion section. This study uses bioinformatics tools to identify and analyse gene sets present in the transcriptome of the corneal endothelium, and tries to find out and observe the relation of the aging effect on the corneal endothelium gene expression. The human subjects had participated voluntarily and informed consent was obtained from all before carrying out any testing procedure. Proper guidelines from the hospital ethical committee were also followed and no harmful chemicals were used on the participants. This study simply aims to raise some awareness of the given topic among the local people so that they are better able to take informed choices about their health in the near future, and also so that they seek medical help when necessary and have no inhibitions in doing so.

Facial expression has lots of applications in human-computer interaction. Although feature extraction and selection have been well studied, the specificity of each expression variation is not fully explored in state-of-the-art works. In this work, the problem of multiclass expression recognition is converted into triplet-wise expression recognition. For each expression triplet, a new feature optimization model based on Action Unit (AU) weighting and patch weight optimization is proposed to represent the specificity of the expression triplet. Sparse representation based approach is then proposed to detect the active AUs of testing sample for better generalization. The algorithm achieved competitive accuracies of 89.67% and 94.09% for Jaffe and CK+ databases, respectively. Better cross-database performance has also been observed.

The COVID-19 pandemic has led to a significant and enduring influence on global health, in-cluding maternal and fetal well-being. Evidence suggests that placental dysfunction is a potential consequence of SARS-CoV-2 infection during pregnancy, which may result in adverse outcomes such as preeclampsia and preterm birth. However, the molecular mechanisms underlying this association remain unclear, and it is uncertain whether a mature placenta can protect the fetus from SARS-CoV-2 infection. To address the above hiatus, we conducted a transcriptome-based study of the placenta in both maternal and fetal compartments. We collected placental samples from 16 women, immediately after term delivery, of which seven had confirmed SARS-CoV-2 in-fection by PCR before parturition. Notably, we did not detect any viral load in either the maternal or fetal compartments of the placenta, regardless of symptomatic status. We extracted total RNA from placental tissues, separately from maternal and fetal compartments, constructed cDNA li-braries, and sequenced them to assess mRNA and small RNA expression. Our analysis revealed 727 differentially expressed genes (DEG) in the maternal placental tissue, with 608 upregulated and 109 downregulated in SARS-CoV-2-positive women compared with healthy, negative wom-en. In contrast, the fetal compartment did not exhibit any significant changes in gene expression with SARS-CoV-2 infection. Specifically, we observed significant downregulation of seven genes belonging to the pregnancy-specific glycoprotein (PSG), related to the immunoglobulin superfam-ily in the maternal compartment with active SARS-CoV-2 infection (fold change range from -13.70 to -5.28, FDR ≤ 0.05). Additionally, comparing symptomatic women with healthy, we identified 14,223 DEGs, with high expression of the inflammatory cytokine IL6 in the maternal placenta of the symptomatic women. Furthermore, KEGG analysis revealed that pathways relat-ed to viral infection, vascular smooth muscle contraction, and oxytocin signaling were altered significantly in symptomatic women. Overall, our study sheds light on the molecular mechanisms underlying the reported clinical risk of preeclampsia and preterm delivery in women with SARS-CoV-2 infection. Nonetheless, studies with larger sample sizes are warranted to further deepen our understanding of the molecular mechanisms of the placenta’s anti-viral effects in ma-ternal SARS-CoV-2 infection.

Recently intraductal carcinoma of the prostate (IDCP) has attracted increasing interest owing to its unfavorable prognoses. To effectively identify the IDCP-specific gene expression profile, we took a novel approach of characterizing a typical IDCP case using spatial gene expression analysis. A formalin-fixed, paraffin-embedded sample was subjected to Visium CytAssist Spatial Gene Expression analysis. IDCP within invasive prostate cancer sites was recognized as a distinct cluster separate from other invasive cancer clusters. Highly expressed genes defining the IDCP cluster, such as MUC6, MYO16, NPY, and KLK12, reflected the aggressive nature of high-grade prostate cancer. IDCP sites also showed increased hypoxia markers HIF1A, BNIP3L, PDK1, and POGLUT1; decreased fibroblast markers COL1A2, DCN, and LUM; and decreased immune cell markers CCR5 and FCGR3A. Overall, these findings indicate that the hypoxic tumor microenvironment and reduced accessibility to immune cells, which reflect the pathological characteristics of IDCP, may influence the aggressiveness of high-grade prostate cancer.

Abstract: Kiwifruit, a perennial horticultural crop known for its nutritional value, faces significant challenges due to biotic and abiotic stresses, affecting its yield and quality. Leucine-rich repeat receptor-like proteins (LRR-RLPs), a major group of receptor-like proteins in plants, have diverse functions in plant physiology, including growth, development, signal transduction, and stress responses. Despite their importance, the specific roles of kiwifruit LRR-RLPs in response to biotic and abiotic stresses remain poorly understood. In this study, we investigated the expression patterns of kiwifruit LRR-RLPs under various stress conditions, including pathogen stress, heat stress, waterlogging stress, and copper stress. Through genome and transcriptome analyses, we identified a total of 101, 164, and 105 LRR-RLPs in Actinidia Chinensis 'Hongyang', Actinidia eriantha 'Huate', and Actinidia chinensis 'Red5' kiwifruit, respectively, classified into 9 subgroups based on phylogenetic relationships. Synteny analysis revealed that the expansion of kiwifruit LRR-RLPs was primarily attributed to segmental duplication events. Using RNA-seq data from pathogen-infected kiwifruits, we identified specific LRR-RLP genes potentially involved in different stages of pathogen infection. Additionally, we observed the potential involvement of kiwifruit LRR-RLPs in abiotic stress responses, with upstream transcription factors possibly regulating their expression. Furthermore, protein interaction network analysis unveiled the participation of kiwifruit LRR-RLP genes in the regulatory network of abiotic stress responses. These findings highlight the crucial roles of LRR-RLPs in mediating both biotic and abiotic stress responses in kiwifruit, offering valuable insights for the breeding of stress-resistant kiwifruit varieties.

Mulberry (Morus alba L), a significant fruit tree crop, requires magnesium (Mg) for its optimal growth and productivity. Nonetheless, our understanding of the molecular basis underlying magnesium stress tolerance in mulberry plants remains unexplored. In our previous study, we identified several differential candidate genes associated with Mg homeostasis via transcriptome analysis including xyloglucan endotransglucosylase/hydrolase (XTH) gene family. The XTH gene family is crucial for plant cell wall reconstruction and stress responses. These genes have been identified and thoroughly investigated in various plant species. However, there is no research pertaining to XTH genes within the M. alba plant. This research systematically examined the M. alba XTH (MaXTH) gene family at the genomic level using a bioinformatic approach. In total, 22 MaXTH genes were discovered and contained the Glyco_hydro_16 and XET_C conserved domains. The MaXTH’s were categorized into five distinct groups by their phylogenetic relationships. The gene structure possesses four exons and three introns. Furthermore, the MaXTH gene promoter analysis reveals a plethora of cis-regulatory elements; mainly stress responsiveness, growth and development, light responsiveness, and phytohormone responsiveness. GO analysis indicated that MaXTHs encode proteins that exhibit xyloglucan xyloglucosyl transferase and hydrolase activities in addition cell wall biogenesis as well as xyloglucan and carbohydrate metabolic processes. Moreover, a synteny analysis unveiled an evolutionary relationship between the XTH genes in M. alba and those in three other species: A. thaliana, P. trichocarpa, and Zea mays. Expression profiles from RNA-Seq data displayed distinct expression patterns of XTH genes in M. alba leaf tissue during Mg treatments. Real-time quantitative PCR analysis confirmed the expression of the MaXTH genes in Mg stress response. Overall, this research enhances our understanding of the characteristics of MaXTH gene family members and lays the foundation for future functional genomic study in M. alba.

Transporter associated with antigen processing 1 (TAP1) gene codes for a transporter protein, which is responsible for tumor antigen presentation in the MHC I or HLA complex. A defect in the gene results in an inadequate tumor tracking. TAP1 may also influence multi drug resistance, which is an extreme threat in providing treatment by drugs which are chemotherapeutic. The gene of TAP1 was analyzed bioinformatically. It gave us prognostic data as a confirmation of whether it should be used as a biomarker. The expression level and pattern analysis were conducted using ONCOMINE, GENT2 and GEPIA2 online platforms. Samples with different clinical outcomes were investigated for expression and promoter methylation analysis was done in cancer vs normal tissues using UALCAN. The copy number alteration and mutation frequency and expression in different cancer studies were analyzed using cBioPortal. The PrognoScan and KM plotter survival analysis of significant data (p-value<0.05) was representing graphically. Pathway and Gene ontology analysis of gene correlated to TAP1 gene was presented using bar charts. After arranging the data in a single panel and correlating expression to prognosis, understanding mutational and alterations and comparing pathways, TAP1 may be a potential novel target to evade a threat against chemotherapy and the study gives new aspects to consider for immunotherapy in human breast, lung, liver and ovarian cancer.

Genome-wide association studies (GWAS) have identified and reproduced thousands of diseases associated loci but many of them are not directly interpretable due to the strong linkage disequilibrium among variants. Transcriptome-wide association studies (TWAS) incorporated expression quantitative trait loci (eQTL) cohorts as reference panel to detect associations with the phenotype at the gene level and were gaining popularity in recent years. For nicotine addiction, several important susceptible genetic variants were identified by GWAS, but TWAS that detected genes associated with nicotine addiction and unveiled the underlying molecular mechanism were still lacking. In this study, we used eQTL data from the Genotype-Tissue Expression (GTEx) consortium as reference panel to conduct tissue specific TWAS on cigarettes per day (CPD) over 13 brain tissues in two large cohorts: UK Biobank (UKBB; N=142,202) and the GWAS &amp; Sequencing Consortium of Alcohol and Nicotine use (GSCAN; N=143,210), and then meta-analyzed the results across tissues while considering the heterogeneity across tissues. We identified three major clusters of genes with different meta-patterns across tissues consistent in both cohorts, including homogenous genes associated with CPD in all brain tissues, partially homogeneous genes associated with CPD in cortex, cerebellum and hippocampus tissues, and lastly the tissue-specific genes associated with CPD in only few specific brain tissues. Downstream enrichment analyses on each gene cluster identified unique biological pathways associated with CPD and provided important biological insights into the regulatory mechanism of nicotine dependence in the brain.

Background: Immunoglobulin A nephropathy (IgAN) is a complex autoimmune disease, and the exact pathogenesis remains to be elucidated. Methods: We conducted summary data-based Mendelian randomization (SMR) analysis and performed functional mapping and annotation using FUMA to explore genetic loci that are po-tentially involved in the pathogenies of IgAN. Both analyses used summarized data of a recent genome-wide association study (GWAS) on IgANs, which included 477,784 Europeans (15,587 cases and 462,197 controls) and 175,359 East Asians (71 cases and 175,288 controls). We performed separate SMR analysis using CAGE and GTEx eQTL data. Results: Using the CAGE eQTL data, our SMR analysis identified 32 probes tagging 25 unique genes that were pleiotropically/potentially causally associated with IgAN, with the top three probes being ILMN_2150787 (tagging HLA-C, PSMR=2.10×10-18), ILMN_1682717 (tagging IER3, PSMR=1.07×10-16) and ILMN_1661439 (tagging FLOT1, PSMR=1.16×10-14). Using GTEx eQTL data, our SMR analysis identified 24 probes tagging 24 unique genes, with the top three probes being ENSG00000271581.1 (tagging XXbac-BPG248L24.12, PSMR=1.44×10-10), ENSG00000186470.9 (tagging BTN3A2, PSMR=2.28×10-10), and ENSG00000224389.4 (tagging C4B, PSMR=1.23×10-9). FUMA analysis identified 3 independent, significant and lead SNPs, 2 genomic risk loci and 39 genes. Conclusion: We identified many genetic variants/loci that are potentially involved in the patho-genesis of IgAN.

Expansins (EXPs) which are the essential components of plant cell wall function as relaxation factors to directly promote turgor-driven expansion of cell wall and control plant growth and development and diverse environmental stress responses. A lot of EXPs genes have been widely identified and characterized in numerous plant species, but not in sweetpotato. In the present study, a total of 59 EXPs genes that unevenly distributed onto 14 of 15 chromosomes were identified from sweetpotato genome, and segmental and tandem duplications were exhibited to act as dominant functions in IbEXPs diversity. Phylogenetic analysis exhibited that IbEXPs members were clustered into four subfamilies according to EXPs in Arabidopsis and rice and the regularity of motif, domain and gene structures was in accordance with the subfamily classification. The collinearity analysis between IbEXPs and related homologous sequences in 9 plants provided the phylogenetic insights of EXPs gene family. Cis-element analysis futhrer revealed the potential roles of IbEXPs in sweetpotato development and stress responses. Evidence from RNA-seq and qRT-PCR analysis of eight selected IbEXPs genes displayed their various tissue specificity in different tissues, and their transcripts were varyingly induced or suppressed under multiple hormone treatments (ABA, SA, JA and ACC) and abiotic stresses (low and high temperature). These results supply worthy foundations to further comprehensively investigate the exact functions of IbEXPs genes, and multiple members may function as promising regulators to control plant development and enhance the stress resistance of transgenic plants.

Mutations that provide environment dependent selective advantages drive adaptive divergence among species. Many phenotypic differences among related species are more likely to result from gene expression divergence rather than from non-synonymous mutations. In this regard, cis-regulatory mutations play an important part in generating functionally significant variation. Some proposed mechanisms that explore the role of cis-regulatory mutations in gene expression divergence involve microsatellites. Microsatellites exhibit high mutation rates and are abundant in both coding and non-coding regions and could influence gene function and products. Here we tested the hypothesis that microsatellites contribute to gene expression divergence among species with 50 individuals from nine closely related Helianthus species using an RNA-seq approach. Differential expression analyses of the transcriptomes revealed that genes containing microsatellites in non-coding regions (UTRs and introns) are more likely to be differentially expressed among species when compared to genes with microsatellites in the coding regions and transcripts lacking microsatellites. We detected a greater proportion of shared microsatellites in 5’UTRs and coding regions compared to 3’UTRs and non-coding transcripts among Helianthus spp. Further, allele frequency differences measured by pairwise FST at single nucleotide polymorphisms (SNPs), indicate greater genetic divergence in transcripts containing microsatellites compared to those lacking microsatellites. A gene ontology (GO) analysis revealed that microsatellite-containing differentially expressed genes are significantly enriched for GO terms associated with regulation of transcription and transcription factor activity. Collectively, our study provides compelling evidence to support the role of microsatellites in gene expression divergence.

Molecular studies regarding regulatory elements such as small ncRNAs and their mechanisms are poorly understood in infectious diseases. Tuberculosis is one of the oldest infectious diseases of humanity, and it is still a challenge to prevent and treat it. The control of the infection as well as its diagnosis are still complex, and treatments used are linked to several side effects. This study aimed to investigate miRNA’s expression profile to identify possible biomarkers for tuberculosis. We applied NGS techniques to investigate miRNA’s global expression profile from blood samples of infected patients with tuberculosis, their respective healthy physicians, and external healthy individuals as controls. Samples from 22 individuals run through a differential expression, target genes, gene set enrichment, and miRNA-gene network analysis. We observed 153 altered miRNAs, among which, only three DEmiRNAs (hsa-let-7g-5p, hsa-miR-486-3p and hsa-miR-4732-5p) were found between the investigated patients and their respective physicians. These DEmiRNAs are suggested to play an important role in granuloma regulation and their immune physiopathology. Our results propose that miRNAs may be involved in immune modulation, regulating the repertoire of genes expressed in the immune system’s cells. Our findings encourage the application of miRNAs as potential biomarkers for tuberculosis.

In the paper, the author presents explicit and unified expressions for a sequence of improper integrals in terms of the beta functions and the Wallis ratios. Hereafter, the author derives integral representations for the Catalan numbers originating from combinatorics.

This paper establishes methods that quantify the structure of statistical interactions within a given data set using the characterization of information theory in cohomology by finite methods, and provides their expression in terms of statistical physic and machine learning. Following [1–3], we show directly that k multivariate mutual-informations (I_k) are k-coboundaries. The k-cocycles are given by I_k = 0, which generalize statistical independence to arbitrary dimension k. The topological approach allows to investigate Shannon’s information in the multivariate case without the assumptions of independent identically distributed variables. We develop the computationally tractable subcase of simplicial information cohomology represented by entropy H_k and information I_k landscapes. The I₁ component defines a self-internal energy functional U_k, and (−1)^k I_k,k≥2 components define the contribution to a free energy functional G_k of the k-body interactions. The set of information paths in simplicial structures is in bijection with the symmetric group and random processes, provides a topological expression of the 2nd law and points toward a discrete Noether theorem (1st law). The local minima of free-energy, related to conditional information negativity and the non-Shannonian cone of Yeung [4], characterize a minimum free energy complex. This complex formalizes the minimum free-energy principle in topology, provides a definition of a complex system, and characterizes a multiplicity of local minima that quantifies the diversity observed in biology. Finite data size effects and estimation bias severely constrain the effective computation of the information topology on data, and we provide simple statistical tests for the undersampling bias and for the k-dependences following [5]. We give an example of application of these methods to genetic expression and cell-type classification. The maximal positive I_k identifies the variables that co-vary the most in the population, whereas the minimal negative I_k identifies clusters and the variables that differentiate-segregate the most. The methods unravel biologically relevant I₁₀ with a sample size of 41. It establishes generic methods to quantify the epigenetic information storage and a unified epigenetic unsupervised learning formalism.

The IgM and IgT classes were previously identified and characterized in the Antarctic teleost Trematomus bernacchii, a species belonging to the Perciform suborder Notothenoidei. Herein we characterized the gene encoding the polymeric immunoglobulin receptor (pIgR) in the same species and compared it to pIgR of multiple teleost species belonging to five perciform sub-orders, including 11 Antarctic and one non-Antarctic (Cottoperca gobio) notothenioid species, the latter living in less cold periantarctic sea. Antarctic pIgR genes displayed particularly long in-trons marked by sites of transposable elements and transcription factors. Furthermore, analysis of T. bernacchii pIgR cDNA unveiled multiple amino acid substitutions unique to Antarctic species, all introducing adaptive features, including N-glycosylation sequons. Interestingly, C. gobio shared most features with the other perciforms rather than with the cold adapted relatives. T. bernacchii pIgR transcripts were predominantly expressed in mucosal tissues, as indicated by q-PCR and in situ hybridization analysis. These results suggest that in cold adapted species pIgR preserved its fundamental role in mucosal immune defense, although remarkable gene structure modifications occurred.

(1) Background: Intracortical microelectrodes (IMEs) are an important part of interfacing with the central nervous system (CNS) and recording neural signals. However, recording electrodes have shown a characteristic steady decline in recording performance owing to chronic neuroinflammation. The topography of implanted devices has been explored to mimic the nanoscale three-dimensional architecture of the extracellular matrix. Our previous work used histology to study the implant sites of non-recording probes and showed that a nanoscale topography at the probe surface mitigated the neuroinflammatory response compared to probes with smooth surfaces. Here, we hypothesized that the improvement in the neuroinflammatory response for probes with nanoscale surface topography would extend to improved recording performance. (2) Methods: A novel design modification was implemented on planar silicon-based neural probes by etching nanopatterned grooves (with a 500 nm pitch) into the probe shank. To assess the hypothesis, two groups of rats were implanted with either nanopatterned (n=6) or smooth control (n=6) probes, and their recording performance was assessed over 4 weeks. Postmortem gene expression analysis was performed to compare the neuroinflammatory response from the two groups. (3) Results: Nanopatterned probes demonstrated increased impedance and noise floor compared to controls. However, the recording performance of the nanopatterned and smooth probes was similar, with active electrode yield for control probes and nanopatterned probes being approximately 50% and 45%, respectively by 4 weeks post-implantation. Gene expression analysis showed 1 gene, Sirt1, differentially expressed out of 152 in the panel. (4) Conclusions: This study provides a foundation for investigating novel nanoscale topographies on neural probes.

In this paper, an analytical form of Rump Function is presented. This seminal function constitutes a fundamental concept of digital signal processing theory and is also involved in many other areas of applied mathematics. In particular, Rump Function is performed in a simple manner as the limit of a sequence of real functions letting ? tend to infinity. This limit is zero for strictly negative values of the real variable ? whereas it coincides with the independent variable ? for strictly positive values of the variable ?. The novelty of this work compared to other research studies concerning analytical expressions of the Ramp Function, is that the proposed formula is not exhibited in terms of miscellaneous special functions, e.g. Gamma Function, Biexponential Function or any other special functions such as Error Function, Hyperbolic Function, Orthogonal polynomials etc. Hence, this formula may be much more practical, flexible and useful in the computational procedures which are inserted into digital signal processing techniques and other engineering practices.

The aim of present study was to link the gene expression profile of selected candidate genes with blood profile, growth performance and carcass traits of Barki lambs. Thirty-eight Barki lambs were divided into 3 groups (fast, intermediate and slow growing) according to growth perfor-mance. Body tissues (muscle, liver and fat) were taken from for RNA isolation and Real-time PCR. The results indicated that, the final body weight hot carcass weight were heavier (P ≤ 0.05) in fast (49.9 Kg and 24.57) than intermediate (40.7 and 19.07 Kg) and slow (30.8 and 15.10 Kg) growing animals. The blood profile of total protein, total lipids, calcium, T3 and T4 hormones did not differ among sheep groups. Genes involved in protein biosynthesis (RPL7), fatty acid oxidation (CPT1) and lipolysis (FABP4) were up regulated in fast and intermediate growing lambs in all studied tissues. While, gene-regulating lipogenesis (ADIPOQ) was expressed simi-larly in fat and liver tissues, but increased its expression in muscle of fast and intermediate growing lambs. Expression of CAPN3 was increased in fast and intermediate growing compared to slow growing lambs. In conclusion, the current study providing an evidence for the im-portance of co-expression of these genes in main body tissues linked with growth performance of Barki lambs.

Cyanobacteria are highly diverse, widely distributed photosynthetic bacteria inhabiting various environments ranging from deserts to the cryosphere. Throughout this range of niches, they have to cope with various stresses and kinds of deprivation which threaten their growth and viability. In order to adapt to these stresses and survive, they have developed several global adaptive responses which modulate the patterns of gene expression and the cellular functions at work. Sigma factors, two-component systems, transcriptional regulators and small regulatory RNAs acting either separately or collectively, for example, induce appropriate cyanobacterial stress responses. The aim of this review is to summarize our current knowledge about the diversity of the sensors and regulators involved in the perception and transduction of light, oxidative and thermal stresses and nutrient starvation responses. The studies discussed here point to the fact that various stresses affecting the photosynthetic capacity are transduced by common mechanisms.

MicroRNA (miRNA) is a typical class of small RNAs that could modulate gene expression in trans at the post-transcriptional level. miRNAs bind to the miRNA binding sites (MBSs) in target mRNAs by sequence complementarity. Alternative splicing (AS) is another commonly occurred process in pre-mRNAs that changes the isoforms of a gene. It is hypothesized that there should be an interaction for gene regulation that involves both AS and miRNA targeting. Studies have verified this hypothesis in the model organism Arabidopsis thaliana. High-throughput sequencing data suggested that in Arabidopsis a considerably large fraction of MBSs are affected by AS events. The overlapping between MBS and AS exceeds the randomly simulated number. Functional experiments have indicated that the AS events are required for the gene expression changes of miRNA targets. Therefore, AS and MBS are mutually favored. The observed expression changes caused by miRNAs could also be contributed by AS events. In the present perspective article, we propose that the AS analysis should be incorporated in the differential-expression analysis of miRNA studies. When defining a differentially-expressed gene, it should be clarified whether the change in gene expression is caused by AS events or solely by miRNA targeting.

Research in various fields of evolutionary biology has shown that divergence in gene expression is a key driver for phenotypic variation. An exceptional contribution of cis-regulatory evolution has for instance been found to contribute to morphological diversification. In the light of these findings, the analysis of genome-wide expression data has become one of the central tools to link genotype and phenotype information on a more mechanistic level. However, in many studies, especially if general conclusions are drawn from such data, a key feature of gene regulation is often neglected. With our article, we want to raise awareness that gene regulation and thus gene expression is highly context dependent. Genes show tissue- and developmental stage-specific expression. We argue that the regulatory context must be considered when studying evolution of gene expression.

Hippocampal neuronal activity generates dendritic and somatic Ca2+ signals, which depending on stimulus intensity, rapidly propagate to the nucleus and induce the expression of transcription factors and genes with crucial roles in cognitive functions. Soluble Amyloid-beta Oligomers (AβOs), the main synaptotoxins engaged in the pathogenesis of Alzheimer's disease, generate aberrant Ca2+ signals in primary hippocampal neurons, increase their oxidative tone and disrupt structural plasticity. Here, we explored the effects of sub-lethal AβOs concentrations on activity-generated nuclear Ca2+ signals and on the Ca2+-dependent expression of neuroprotective genes. To induce neuronal activity, neuron-enriched primary hippocampal cultures were treated with the GABAA receptor blocker gabazine (GBZ), and nuclear Ca2+ signals were measured in AβOs-treated or control neurons transfected with a genetically encoded nuclear Ca2+ sensor. Incubation (6 h) with AβOs significantly reduced the nuclear Ca2+ signals and the enhanced phosphorylation of cyclic AMP response element binding protein (CREB) induced by GBZ. Likewise, incubation (6 h) with AβOs significantly reduced the GBZ-induced increases in the mRNA levels of Neuronal Per Arnt Sim domain protein 4 (Npas4), Brain-derived Neurotrophic Factor (Bdnf), Ryanodine Receptor type-2 (RyR2), and the antioxidant enzyme NADPH-Quinone-Oxidoreductase (Nqo1). Based on these findings we propose that AβOs, by inhibiting the generation of activity induced nuclear Ca2+ signals, disrupt key neuroprotective gene expression pathways required for hippocampal-dependent learning and memory processes.

Purpose: This study aims to analyze the factors that affect the behavior of user-defined avatars in different virtual environments, and compare the differences in public self-awareness, self-expression, and emotional expression among customized avatars in multiple virtual contexts. Methods: Using a between-subjects experimental design, two random groups of participants were asked to customize avatars for themselves in two contexts, a multiplayer online social game (MOSG) and a virtual meeting (VM). i.e. a relaxed and a serious social environment. Results: When subjects perceived a more relaxed environment, the customized avatars had less self-similarity, and the subjects exhibited a stronger self-disclosure willingness and enhanced avatar wishful identification; nevertheless, public self-consciousness was not increased. When subjects perceived a more serious environment, the customized avatars exhibited a higher degree of self-similarity, and the subjects exhibited a greater self-presentation willingness, along with enhanced identification of avatar similarity, and increased public self-consciousness. Conclusions: Participants expressed positive emotions, suggesting that avatars play a positive role in various virtual contexts. The virtual context affects the self-similarity of user-customized avatars, and avatar self-similarity affects self-presentation and self-disclosure willingness, and these factors will affect the behavior of the user-customized avatar. This study contributes suggestions to the Metaverse avatar customization platform design.

Decades of research identified numerous gene biomarkers of cardiac diseases whose restored sequence or/and expression level was hoped to recover the normal cardiac function. However, each human has unique and dynamic pathophysiological characteristics resulting from the unrepeatable combination of favoring factors such are: race, sex, age, medical history, diet, stress, exposure to toxins, habits etc. As such, no treatment fits everybody and finding personalized solutions is a top priority for medicine of 21st century. The Genomic Fabric Paradigm (GFP) provides the most theoretically possible comprehensive characterization of the transcriptome, its alterations in disease and recovery following a treatment. By attaching to each gene the independent average expression level, expression variation and expression coordination with each other gene, GFP delivers thousands times more information than the traditional analysis. This report presents the theoretical bases of the GFP and some applications to our microarray data from mouse models of post ischemic, and constant and intermittent hypoxia-induced heart failure. The GFP analyses revealed novel transcriptomic aspects of the gene expression control and networking under ischemic conditions. Through all-inclusive characterization of the transcriptome and the unrepeatable gene hierarchy in each condition, GFP is an essential avenue towards development of a truly personalized cardiogenomic therapy.

The AREB/ABF (ABA response element binding) proteins in plants are vital for plant stress responses, but the understanding of AREB/ABFs in orchid species, an important traditional medicinal and ornamental plants, is limited. Here, twelve AREB/ABF genes were identified within three orchids’ complete genome and classified into three groups (groups I, II, and III) via a phylogenetic analysis, which was further supported by the analysis of their conserved motifs and gene structures. The cis-elements analysis showed that hormone response elements followed by light and stress responses were widely rich in AREB/ABFs. A prediction analysis of orchid ABRE/ABF-mediated regulatory network was further constructed through cis-regulatory elements (CREs) analysis of their promoter regions and revealed that several transcriptional factor (TF) gene family were abundant as the potential regulators for orchid AREB/ABFs. Expression profile analysis based on public transcriptomic data showed that most AREB/ABF genes have distinct tissue-specific expression patterns in orchid plants. Additionally, DcaABI5 as the homolog of ABA INSENSITIVE 5 (ABI5) from Arabidopsis was selected for further analysis. The results showed that Arabidopsis transgenic overexpressing DcaABI5 could rescue the ABA insensitive phenotype in the mutant abi5. Collectively, these findings will provide valuable information on AREB/ABF genes in orchids.

Zanthoxylum armatum is an forest tree with apomixis characteristic, and RWP-RKs are plant-specific transcription factors known to regulate nitrogen response and embryonic development. However, the genome-wide analysis of RWP-RK family genes in Z.armatum is unclear. In this study, 36 ZaRWP-RK transcription factors were identified in the genome of Z.armatum, among which 15 genes were belonged to the RKD subfamily and 21 belonged to the NLP subfamily. Duplication events of ZaRWP-RK genes were mainly segmental duplication, and synteny analysis revealed close phylogenetic relationship between Z.armatum and Arabidopsis. The analysis of cis-elements indicated that ZaRWP-RK genes may be involved in the regulation of embryonic development of Z.armatum by responding to plant hormones such as abscisic acid, auxin, and gibberellin. Real-time PCR showed that the expression levels of most ZaRWP-RK genes were significantly increased from flowers to young fruits. Protein-protein interaction network analysis further revealed the potential roles of the ZaRWP-RK proteins in apomixis. Collectively, this study provide a theoretical basis for future investigations into the ZaRWP-RK genes and their regulatory mechanisms in the apomixis process of Z.armatum.

Subclinical inflammation in protocol biopsies relates to tacrolimus exposure and human leukocyte antigen (HLA) matching. We aimed to characterize transcripts associated with rejection and tacrolimus exposure and its association with transplant outcome. We tested whether gene expression is associated with rejection using strictly normal protocol biopsies (n = 17) and biopsies with TCMR or ABMR according to Banff criteria (n = 12). Subsequently, we analyzed these transcripts in a set of 4-month protocol biopsies (n = 137) to assess their association with donor and recipient characteristics, the intensity of immunosuppression, and graft outcome. Differential expression (FDR&gt;3) between normal and rejection biopsies yielded a set of 111 genes. Nineteen out of these 111 genes were correlated with tacrolimus trough levels (TAC-C0) at the time of biopsy in the protocol biopsy cohort (n=137) and allowed us to define 2 clusters of biopsies. The 2 clusters differed in donor age and tacrolimus levels. Subclinical rejection, including borderline lesions, tended to occur in the same cluster. Logistic regression analysis indicated that TAC-C0 at the time of biopsy (OR: 0.83, 95%CI:0.72-0.06, p=0.0117) was associated with cluster 2. In a follow-up averaging 70±30 months, this patient group displayed a significant decline in renal function (p=0.0135). The expression of rejection-associated transcripts in early protocol biopsies is associated with tacrolimus exposure and a faster decline in renal function.

Background: We screened key genes involved in the dedifferentiation process of thyroid cancer by analyzing differentially expressed genes in anaplastic thyroid carcinoma (ATC) using Gene Expression Omnibus (GEO) database. Methods: We analyzed differentially expressed genes in ATC from GEO database and constructed a co-expression network using Weighted Gene Co-expression Network Analysis (WGCNA) . Cytoscape was used to calculate top ten key genes based on their gene scores. DEPDC1 was selected as interest gene for further validation, including examination of its prognostic relationship with thyroid cancer, correlation analysis with thyroid differentiation markers, immune cell infiltration, and tumor microenvironment. The impact of DEPDC1 on cancer cell migration, invasion, and proliferation was validated through downregulation and overexpression of DEPDC1 gene in thyroid cancer cell lines c643 and BCPAP. Results: We identified top 10 key genes in ATC: AURKA, TPX2, RACGAP1, MELK, DLGAP5, DEPDC1, KIF2C, NCAPG, CCNB1, and PBK. DEPDC1 showed a strong correlation with the prognosis of thyroid cancer. The odds ratio (OR) for overall survival was 3.5 (P=0.016) in individuals with high DEPDC1 expression, and the OR for disease-free survival was 2.6 (P=0.003). Additionally, the analysis revealed a significant negative correlation between the expression of DEPDC1 and thyroid differentiation markers. Gene Ontology enrichment analysis demonstrated that DEPDC1 co-expressed genes were enriched in cellular processes, cell cycle pathways, and dynamic protein groups. Matrix and immune scoring indicated tumor microenvironment and immune scores were associated with DEPDC1 expression. TIMER2.0 online analysis revealed a strong correlation between DEPDC1 expression and immune cell infiltration in papillary thyroid carcinoma. After downregulation of DEPDC1 expression, the invasion, migration, and proliferation abilities of tumor cells were significantly reduced. Conversely, overexpression of DEPDC1 exhibited the opposite effect. Conclusion: DEPDC1, a newly identified dedifferentiation gene in thyroid cancer, inversely correlated with the degree of differentiation. High expression of DEPDC1 promotes cancer progression and leads to poor prognosis in thyroid cancer.

Coordination of activities between nuclei and organelles in plant cells involves information exchange, in which phytohormones may play an essential role. Therefore, dissection of the mechanisms of hormone-related integration between phytohormones and mitochondria is an important and challenging task. Here, we found that inputs from multiple hormones may cause changes in transcript accumulation of mitochondrial-encoded genes and nuclear genes encoding mitochondrial (mt) proteins. In particular, treatments with exogenous hormones induced changes in GUS expression in the reporter line possessing a 5'-deletion fragment of the RPOTmp promoter. These changes corresponded in part to up- or downregulation of RPOTmp in wild-type plants, which affected the transcription of mt-encoded genes, implying that promoter fragments of the RPOTmp gene are functionally involved in responses to IAA (indole-3-acetic acid), ACC (1-aminocyclopropane-1-carboxylic acid), and ABA (abscisic acid). Hormone-dependent modulations in the expression of mt-encoded genes can also be mediated through mitochondrial transcription termination factors 15, 17, and 18 of the mTERF family and genes for tetratricopeptide repeat proteins that are coexpressed with mTERF genes, in addition to SWIB5 encoding a mitochondrial SWI/SNF (nucleosome remodelling) complex B protein. These genes specifically responded to hormone treatment, displaying both negative and positive regulation in a context-dependent manner. According to bioinformatic resources, their promoter regions possess putative cis-acting elements involved in responses to phytohormones. Alternatively, hormone-related transcriptional activity of these genes may be modulated indirectly, which is especially relevant for brassinosteroids (BS). In general, the results of the study indicate that hormones are essential mediators that are able to cause alterations in the transcript accumulation of mt-related nuclear genes, which in turn trigger the expression of mt genes.

Autophagy is a highly conserved self-degradation process involving the degradation and recycling of cellular components and organelles. Although the involvement of autophagy in metabolic changes during fruit ripening has been preliminarily demonstrated, the variations in autophagic flux and specific functional roles in fruit ripening remain to be elucidated. In this study, we analyzed the occurrence of autophagy during tomato fruit ripening. The results revealed differential expression of SlATG8 family member during tomato fruit ripening. Transmission electron microscopy observations and dansylcadaverine (MDC) staining confirmed the presence of autophagy at the cellular level in tomato fruits. Furthermore, the overexpression of SlATG8f induced the formation of autophagosomes, increased autophagic flux within tomato fruits, effectively enhanced the expression of ATG8 proteins during the color transition phase of fruit ripening, thereby promoting tomato fruit maturation. This overexpression also led to the accumulation of vitamin C (VC) and soluble solids while reducing acidity in the fruit. Collectively, our findings highlight the pivotal role of SlATG8f in enhancing tomato fruit ripening, providing insights into the mechanistic involvement of autophagy in this process. This research contributes to a better understanding of key factors regulating tomato fruit quality and offers a theoretical basis for tomato variety improvement.

Water deficit greatly affects global crop growth and productivity, particularly in water-limited environments like upland rice cultivation, leading to reduced grain yield. Plants activate various defense mechanisms during water deficit, involving numerous genes and complex metabolic pathways. Exploring orthologous genes, linked to enhanced drought tolerance, can aid functional validation in target species using genomic data from model organisms. Here, we studied the OsCPK5 gene in upland rice, an AtCPK6 ortholog from A. thaliana, by overexpressing it in the BRSMG Curinga cultivar. Transformants were assessed under two conditions: water deficit applied 79 days after seeding, lasting 14 days, followed by 7 days of irrigation at 80% field capacity. Physiological data and samples were collected at R3, R6, and R8 stages. GM plants consistently exhibited higher OsCPK5 gene expression across stages, peaking during grain filling. During this phase, GM plants displayed reduced stomatal conductance, photosynthetic rate, and increased water use efficiency compared to non-GM plants under drought. GM plants also exhibited higher filled grain percentage in both irrigation conditions. Their drought susceptibility index was 0.9 times lower than NGM plants, and they maintained a higher chlorophyll a/b index, indicating sustained photosynthesis. NGM plants under water deficit exhibited more leaf senescence, while OsCPK5-overexpressing plants retained green leaves. Overall, OsCPK5 overexpression induced diverse drought tolerance mechanisms, hinting at potential for future development of more drought-tolerant rice cultivars.

This paper presents ConF, a novel deep learning model designed for accurate and efficient prediction of non-coding RNA families. NcRNAs are essential functional RNA molecules involved in various cellular processes, including replication, transcription, and gene expression. Identifying ncRNA families is crucial for comprehensive RNA research, as ncRNAs within the same family often exhibit similar functionalities. Traditional experimental methods for identifying ncRNA fam-ilies are time-consuming and labor-intensive. Computational approaches relying on annotated secondary structure data face limitations in handling complex structures like pseudoknots and have restricted applicability, resulting in suboptimal prediction performance. To overcome these chal-lenges, ConF integrates mainstream techniques such as residual networks with dilated convolutions and cross multi-head attention mechanisms. By employing a combination of dual-layer convolu-tional networks and BiLSTM, ConF effectively captures intricate features embedded within RNA sequences. This feature extraction process leads to significantly improved prediction accuracy compared to existing methods. Experimental evaluations conducted on a ten-fold publicly available dataset demonstrate the superiority of ConF in terms of accuracy, sensitivity, and other perfor-mance metrics. Overall, ConF represents a promising solution for accurate and efficient ncRNA family prediction, addressing the limitations of traditional experimental and computational methods.

Prostate cancer (PCa) Lung metastases are rarely resected, therefore PCa lung metastases are insufficiently molecularly characterized. We recently identified carcinoembryonic antigen cell adhesion molecule 6 (CEACAM6) as a potential driver of pulmonary metastatic spread. Here, we show the biological significance of CEACAM6 in PCa cell proliferation, apoptosis and migration. CEACAM6 was silenced by siRNA in PC3 cells. Functional assessment of apoptosis, cell viability, proliferation and migration were performed in siRNA-CEACAM6 PC3 cells. Non-treated and scrambled (scr) RNA PC3 cells were used as control. Following a specific knockdown of CEACAM6 in PC3 cells, the expression of CEACAM6 protein was significantly decreased in comparison to controls. Cell viability and cell counts decreased in CEACAM6 silent PC3 cells. In contrast, caspase-3 activity was highly elevated in siRNA-CEACAM6 PC3. Furthermore, by performing a cell scratch assay, the migration ratio in siRNA-CEACAM6 PC3 cells were significantly diminished compared with the control group after 48 hours of post transfection incubation. CEACAM6 as a cell adhesion molecule has been implicated in promoting metastatic disease in several solid tumours such as colorectal or gastric cancer. We could show that silencing of CEACAM6 has a significant functional effect on PCa cells. CEACAM6 might play an important role in fostering metastatic spread to the lung of PCa patients via enhancing proliferation and suppressing apoptosis. CEACAM6 might therefore pose an attractive therapeutic target to prevent metastatic disease.

At 2.3 megabases in length, the dystrophin gene is enormous: transcription of a single mRNA requires approximately 16 hours. Principally expressed in skeletal muscle, the dystrophin protein product protects the muscle sarcolemma against contraction-induced injury, and dystrophin deficiency results in the fatal muscle-wasting disease, Duchenne muscular dystrophy. This gene is thus of key clinical interest, and therapeutic strategies aimed at eliciting dystrophin restoration require quantitative analysis of its expression. Approaches for quantifying dystrophin at the protein level are well established, however study at the mRNA level warrants closer scrutiny: measured expression values differ in a sequence-dependent fashion, with significant consequences for data interpretation. In this manuscript we discuss these nuances of expression and present evidence to support a transcriptional model whereby the long transcription time is coupled to a short mature mRNA half-life, with dystrophin transcripts being predominantly nascent as consequence. We explore the effects of such a model on cellular transcriptional dynamics, and then discuss key implications for the study of dystrophin gene expression, focussing both on conventional (qPCR) and next-gen (RNAseq) approaches.

Diabetes is a chronic condition that worsens living conditions and causes significant problems, particularly in the vascular system. Many different treatment options are used to prevent these negative effects of diabetes. The most important of these is changing living conditions and exercising regularly.This study aimed todetermining the impacts of a 6-week anaerobic exercise protocol on the blood glucose levels and Glut4 gene expression in the muscle tissues of streptozocin (STZ)-induced diabetic rats. A total of 45 rats were separate control, sedentary, and exercise groups (n=15 each). The STZ used to induce diabetes in rats was applied once with a single intraperitoneal injection. No diabetes was created in the control group, and no exercise was performed, whereas in the sedentary group, diabetes was induced, and no exercise was applied. It was determined that the Glut4 gene expression in the muscle tissue of the exercise group was importantly increased when compared with the sedentary diabetes group.Glut4 gene is an important gene involved in blood glucose regulation and anaerobic exercise significantly increased the expression level of this gene in our study.

This concept paper gives a narrative about intelligence from insects to the human brain, showing where evolution may have been influenced by the structures in these simpler organisms. The ideas also come from the author's own cognitive model, where a number of algorithms have been developed over time and the precursor structures should be codable to some level. Through developing and trying to implement the design, ideas like separating the data from the function have become architecturally appropriate and there have been several opportunities to make the system more orthogonal. Similarly for the human brain, neural structures may work in-sync with the neural functions, or may be slightly separate from them. Each section discusses one of the neural assemblies with a potential functional result, that cover ideas such as timing or scheduling, structural intelligence and neural binding. Another aspect of self-representation or expression is interesting and may help the brain to realise higher-level functionality based on these lower-level processes.

For developmental processes we know most of the gene networks controlling specific cell responses. We still have to determine how these networks cooperate and how signals become integrated. The JNK pathway is one of the key elements modulating cellular responses during development. Yet, we still know little on how the core components of the pathway interact with additional regulators or how this network modulates cellular responses in the whole organism in homeostasis or during tissue morphogenesis. We have performed a promoter analysis searching for potential regulatory sequences of puc and identified different specific enhancers directing gene expression in different tissues and at different developmental times. Remarkably, some of these domains respond to the JNK activity, but not all. Altogether, these analyses show that puc expression regulation is very complex and that JNK activities participate in non-previously known processes during the development of Drosophila.

Background: Hypoxia-inducible factors (HIFs) are transcription factors that get activated and stabilized in the heterodimerized form under hypoxic conditions. The three members of the HIF alpha factors share high structural similarity but have tissue-specific expression patterns. A majority of studies have reported the importance of the HIF1A and HIF2A activity in the survival, proliferation, metastatic potential, and metabolic regulation of hypoxic cancer cells. However, the importance of the expression pattern and activity of HIF3A in a variety of cancers remains unknown. Method and materials: The expression profile of 13 different types of The Cancer Genome Atlas (TCGA) cancer samples were downloaded, normalized and differential gene expression analysis (DGE) was performed to compare the expression pattern of HIF alpha family members in cancer and adjacent normal tissues, as well as at different stages and tumor-sizes. Receiver operating characteristic (ROC) test and survival analysis were carried out to estimate the diagnostic potential of HIF alpha isomers in different cancers, as well as the survival rate of patients with the varying expression levels of HIF alpha factors. Results: The expression status of HIF3A was notably less in all cancer samples in contrast to their adjacent normal tissues. The expression degree of HIF1A varied among distinct types of cancer and the expression degree of HIF2A was lower in nearly all types of cancers. The expression level of HIF alpha isomers did not significantly correlate with different sizes of tumor samples and stages of different tumor tissue samples. HIF3A had very weak diagnostic potential, while HIF2A had better diagnostic potential in most types of cancers compared to HIF1A. Patients who had a higher level of HIF3A had better survival, while the higher expression levels of HIF1A and HIF2A were associated with worse survival in many types of cancers. Conclusion: Our study shows the heterogenous expression pattern of HIF alpha subunits in distinctive kinds of cancers and the influence of HIF3A expression level in the survival of patients with varying types of cancers.

The Brassica family contains several economically important crops, including rapeseed (Brassica napus, 2n = 38, AACC), the second largest source of seed oil and protein meal worldwide. However, research in rapeseed is hampered because it is complicated and time-consuming for researchers to access different types of expression data. We therefore developed the Brassica Expression Database (BrassicaEDB, https://biodb.swu.edu.cn/brassica/) for the research community. We conducted RNA sequencing (RNA-Seq) of 103 tissues from rapeseed cultivar ZhongShuang11 (ZS11) at seven developmental stages (seed germination, seedling, bolting, initial flowering, full-bloom, podding, and maturation). We determined the expression patterns of 101,040 genes via FPKM analysis and displayed the results using the eFP browser. We also analyzed transcriptome data for rapeseed from 70 BioProjects in the SRA database and obtained three types of expression level data (FPKM, TPM, and read counts). We used this information to develop the BrassicaEDB, including eFP, Treatment, Coexpression, and SRA Project modules based on gene expression profiles and Gene Feature, qPCR Primer, and BLAST modules based on gene sequences. The BrassicaEDB provides comprehensive gene expression profile information and a user-friendly visualization interface for Brassica crop researchers. Using this database, researchers can quickly retrieve the expression level data for target genes in different tissues and in response to different treatments to elucidate gene functions and explore the biology of rapeseed at the transcriptome level.

MicroRNAs are small non-coding RNAs that are implicated in several physiological processes such as cell development, proliferation, differentiation, apoptosis, immune response and angiogenesis. In the last couple of decades, several studies on miRNA profiling in OSCC have tried to associate miRNAs with clinical characteristics such as relapse, metastasis, and survival, however, the results have been found to vary considerably, sometimes even being contradictory. The main objective of our study was to analyse and verify miRNA expression in oral squamous cell carcinoma in a Spanish population. Second, we attempted to associate the identified deregulated miRNAs with molecular pathways. 8 Oral Squamous Cell Carcinoma patients and 8 healthy control samples were analysed by a microarray Affymetrix® miRNA 4.1 array plate and validated with 8 more cases using RT-qPCR. Deregulated miRNAs were studied using Diana Tools miRPath 3.0 to associate miRNA targets with molecular pathways. Microarray analysis identified 80 deregulated miRNAs (35 over-expressed and 45 under-expressed). Only miR-497-5p and miR-4417 maintained its deregulated expression after validation with qPCR. Among the molecular pathways in which deregulated miRNAs could be implicated, the most statistically significant pathway was ‘proteoglycans in cancer’. No relationship was found between miR-497-5p or miR-4417 expression and clinical or pathological parameters except of nodular affectation (high miR-4417 expression in patients with nodular affectation, p = 0.035) and radiotherapy (diminished miR-497-5p expression in patients who needed radiotherapy, p = 0.05). We have verified the altered expression of miR-497-5p and miR-4417 in Oral Squamous Cell Carcinoma samples and related the deregulated miRNAs with the ‘proteoglycans in cancer’ pathway.

Alfalfa, an important legume forage, is an ideal crop for sustainable agriculture and a potential bioenergy plant. Drought, one of the most common environmental stresses, substantially affects plants’ growth, development and productivity. MicroRNAs (miRNAs) are newly discovered gene expression regulators that have been linked to several plant stress responses. To elucidate the role of miRNAs in drought stress regulation of alfalfa, a high-throughput sequencing approach was used to analyze 12 small RNA libraries comprising of 4 samples, each with 3 biological replicates. We identified 348 known miRNAs, belonging to 80 miRNA families, from the 12 libraries and 281 novel miRNAs using Mireap software. 18 known miRNAs in roots and 12 known miRNAs in leaves were screened out as drought-responsive miRNAs. Except for miR319d and miR157a which were upregulated under drought stress, the expression pattern of drought-responsive miRNAs were different between roots and leaves in alfalfa. This is the first study discovering miR157a, miR1507, miR3512, miR3630, miR5213, miR5294, miR5368 and miR6173 are drought-responsive miRNAs. Target transcripts of drought-responsive miRNAs were computationally predicted. All 447 target genes for the known miRNAs were predicted using an online tool. This study provides a significant insight on understanding drought-responsive mechanisms of alfalfa.

: F-box proteins are a large gene family in plants, and play crucial roles in plant growth, development, and stress response. However, a comprehensive investigation of F-box family genes in cultivated peanut (Arachis hypogaea L.), has not been undertaken. In this study, we used a bioinformatics approach to identify and analyze F-box protein family genes in the cultivated peanut genome, and examined their expression profiles. Herein, a total of 95 F-box protein family members on 18 chromosomes were identified in cultivated peanut, which were sequentially designed as AhFBX1-AhFBX95. Phylogenetic analysis classified these AhFBXs into four different groups (Group I - IV), and eight subfamilies based on the structural domains of the C-terminal. Structural analysis revealed the presence of ten unique motifs within these genes. Additionally, the investigation of non-synonymous and synonymous substitution rates (Ka/Ks), and synteny in the F-box proteins, provided valuable insights into the evolutionary and domestication processes of AhFBX genes in peanut. Furthermore, in silico expression profiling of 22 different tissues demonstrated that most genes were predominantly expressed in leaves, pistils, and the aerial gynophore tip, while some members showed tissue-specific expression at the apex of both vegetative and reproductive shoots. Subsequently, we validated the expression of 9 differentially expressed genes (DEGs) using transcriptomeic data of lateral branch development and qRT-PCR analysis. These DEGs may be involved in regulating the development of lateral branches and subsequently affect the architecture of peanut plants.

Abstract: There are a large number of cross-linked chemical bonds such as tight disulfide bonds, hydrogen bonds,and salt bonds in keratin, which makes it difficult to break the structural stability of the protein using conventional methods. The degradation of feather waste using keratinases produced by microorganisms has the advantage of being efficient and environmentally friendly. In this study, 23 strains capable of producing transparent circles were obtained via preliminary screening with skim milk as the sole carbon and nitrogen source, and 9 strains capable of degrading feathers were obtained via rescreening with feathers as the sole carbon and nitrogen source. Finally, through the determination of enzyme activity, the strain CY-A with a high ability to degrade feathers was obtained. The strain was identified as Bacillus tequilensis via morphological observation, physiological and biochemical identification, and 16 S rRNA molecular biology. The keratinase gene (bta) of Bacillus tequila was cloned, and its full-length gene sequence was 1110 bp, containing an open reading frame of 1089, encoding 369 amino acids (GeneBank accession number OR733336). The physicochemical properties, secondary structure, and tertiary structure of the recombinant protein were analyzed using bioinformatics tools such as ProtParm, SOPMA, MEGA, and SWISS-MODELS. It was found that the relative molecular mass was 37.953 kD and that the isoelectric point was 8.48. The protein was a hydrophilic protein. The keratinase bta gene was ligated to the expression vector pET28a(+) using T4 ligase, and the recombination expression plasmid pET28-bta was constructed and transformed into E. coli BL21(DE3) to obtain recombinant engineered bacteria. The expressed keratinase was purified using a Ni2+ affinity chromatography column to obtain recombinant bta keratinase bta with a molecular weight of 37.953 kD. The identification and characterization of this keratinase gene provides some theoretical basis for the further enhancement of keratinase activity by means of genetic engineering.

Epitranscriptomics refers to post-transcriptional regulation of gene expression via RNA modifications and editing that affect RNA functions. Many kinds of modifications of mRNA have been described, among which N6-Methyladenosine (m6A), N1-methyladenosine (m1A), 7-methylguanosine (m7G), pseudouridine (Ψ), 5-methylcytidine (m5C). They alter mRNA structure and consequently stability, localization and translation efficiency. Perturbation of the epitranscriptome is associated with human diseases, thus opening the opportunity for potential manipulations as therapeutic approach. In this review, we overview the functional roles of epitranscriptomic marks in the skeletal muscle system, in particular in embryonic myogenesis, muscle cell differentiation and muscle homeostasis processes. Further, we explored high throughput epitranscriptome sequencing data to identify RNA chemical modifications in muscle-specific genes and we discuss the possible functional role and the potential therapeutic applications.

Background: TP53 alterations have a significant prognostic effect in myeloid neoplasms. Our objective was to investigate the TP53 gene mutation status, p53 protein expression, and their re-lationship in dysplasia-related myeloid neoplasms with varying levels of myeloblast counts. Methods: 76 bone marrow biopsy samples with different blast counts were analyzed. Total and strong (3+) p53 expression was determined. Dual immunohistochemical staining was performed to determine the cell population associated with p53 expression. NGS analysis was performed using the Accel-Amplicon Comprehensive TP53 panel. Results: Both p53 expression and TP53 VAF showed a significant correlation with the myeloblast ratio (p<0.0001), however, p53 expres-sion was present in other cell lineages as well. The VAF value exhibited a significant correlation with p53 expression. A high specificity (0.9800) was observed for TP53 mutation using the ≥10% strong (3+) p53 cut-off value, although the sensitivity (0.4231) was low. Conclusion: Strong (3+) p53 expression using a ≥10% cut-off value accurately predicts TP53 mutation but doesn't reveal the allelic state. The p53 expression is significantly influenced by myeloblast count, and histo-logical interpretation should consider the presence of intermixed non-neoplastic marrow cells with varying physiological p53 expression.

Cell identity is determined by chromatin structure and the profiles of gene expression, which are dependent on chromatin accessibility and DNA methylation of regions critical for gene expression, such as enhancers and promoters. These epigenetic modifications are required for mammalian development and are essential for the establishment and the maintenance of the cellular identity. DNA methylation was once thought to be a permanent repressive epigenetic mark, but systematic analysis in various genomic contexts reveals a more dynamic regulation than previously thought. In fact, both active DNA methylation and demethylation occur during cell fate commitment and terminal differentiation. To link methylation signatures of specific genes to their expression profiles, we have determined the methyl-CpG configurations of the promoters of five genes switched on and off during murine postnatal brain differentiation by bisulfite-targeted sequencing. We report here the structure of significant, dynamic and stable methyl-CpG profiles associated with silencing or activation of expression of genes during neural stem cell and brain postnatal differentiation. Strikingly, these methylation cores mark different mouse brain areas and cell types derived from the same areas during differentiation.

Matrix metalloproteinases (MMPs) are often considered biomarkers of skin fibrosis. At the early stages of the pathological process, an elevation of their enzymatic activity causes significant changes in the composition of the extracellular matrix. MMPs secreted by immune cells facilitate their migration to the site of damage. Then, the immune cells eliminate the affected cells and biomolecules. Moreover, bidirectional changes in the activity of proteolytic enzymes, including MMPs, accompany wound healing. This study aimed to assess changes in the expression of Mmp2, Mmp3, and Mmp9 after treating the mice with laser therapy using the experimental model of bleomycin-induced skin fibrosis. Using immunohistochemistry, we characterized the histological features of scarred skin. We also analyzed changes in the expression of MMPs using real-time polymerase chain reaction before and after the irradiation with laser. We showed that treatment of the mice with CO2 laser partially normalized the histological features of scarred skin. We also noticed a decrease in the expression of Mmp2, Mmp3 (both p &lt; 0.05), and Mmp9 (p = 0.065) during scar healing. The obtained results suggest that normalization of skin homeostasis requires a control of MMPs activities via induction of their genes.

Drought stress on soybean is a research-demanding matter for negative influence that agricultural drought brings about. This study was designated to evaluate the effect of drought stress on some gene expression in flowering and pod elongation stages in soybean. This experiment was carried out in split-plot format with RCBD design with four replicates. Drought stress as the main factor included three levels (irrigation after 50, 100, and 150 millimeters evaporation from the A-class evaporation pot) of which 50 millimeters evaporation is considered as control. The sub-factor included a factorial combination of 3 varieties (DPX, Sari and WE6) and two sampling stages (flowering and pod elongation). The gene expression analysis was carried out by using the QRT-PCR technique. According to our results, all genes have shown overexpression in drought stress despite this result was not the same for all genotypes and stress levels. Some genes have up-regulated in mediate stress (treatment 100) level (like as Gmdreb 2, Gmdreb 5, GmRD20A, GmaxACD2) and other genes up-regulated in serve stress (treatment 150) level. Between genotypes, DPX cultivar and WE6 line were better than of the sari cultivar for all genes up-regulated.

Background Polycystic ovary syndrome (PCOS) is a common hormonal disorder worldwide among women of reproductive age. It is characterized by endocrine, reproductive, and metabolic abnormalities. Insulin resistance (IR) is one of its most important clinical features, associated with metabolic disorders and increased risk of type 2 diabetes (T2D). This study aimed to explore the whole blood gene expression profiling related to IR in PCOS patients compared to controls. Methods Blood RNA was extracted from 5 PCOS and 5 non-PCOS women with matched age and BMI. Homeostasis model assessment (HOMA-IR) was used to estimate the IR. The expression of IR genes was analyzed by Profiler PCR array. Results Both groups have similar levels of HOMA-IR (p>0.05). However, differential expression levels were observed between them. Fourteen genes were upregulated and 26 genes were downregulated in PCOS samples. Among the upregulated genes (>2 fold-change, p-value<0.05) are ADIPOQ, ADIPOR1, OLR, IGF-1, and APOE. Downregulated genes (>-2 fold-change, p-value<0.05) include HK-2, IRS1, and SERPINE1. These genes are involved in insulin and adipokines signaling, commonly dysregulated in T2D. They are also involved in innate immunity and inflammatory processes and are essential for lipid and carbohydrate metabolism. Conclusion Our finding suggests that despite both groups having no difference in IR level, there are differentially expressed genes involved in the IR pathway.

Given the importance of menstrual blood in the pathogenesis of endometriosis and the multifunctional roles of menstrual mesenchymal stem cells (MenSCs) in regenerative medicine, this issue has gained prominence in the scientific community. Moreover, recent reviews highlight how robust the integrated assessment of omics data is for endometriosis. To our knowledge, no study has applied the multi-omics approaches to endometriosis MenSCs. It is a case-control study at a university-affiliated hospital. MenSCs transcriptome and proteome data were obtained by RNA-seq and UHPLC-MS/MS detection. Among the differentially expressed proteins and genes, we emphasize ATF3, ID1, ID3, FOSB, SNAI1, NR4A1, EGR1, LAMC3, and ZFP36 genes and MT2A, TYMP, COL1A1, COL6A2, and NID2 proteins that were already reported in the endometriosis. Our functional enrichment analysis reveals integrated modulating signaling pathways such as epithelial-mesenchymal transition (↑) and PI3K signaling via AKT to mTORC1 (↓in proteome), mTORC1 signaling, TGF beta signaling, TNFA signaling via NFkB, and response to hypoxia via HIF1A targets (↑in transcriptome). Our findings highlight primary changes in the endometriosis MenSCs, suggesting that the chronic inflammatory endometrial microenvironment can modulate these cells, providing opportunities for endometriosis etiopathogenesis. Moreover, they identify challenges for future research leveraging knowledge for regenerative and precision medicine in endometriosis.

Middle East Respiratory Syndrome coronavirus (MERS-CoV) infects dromedary camels and zoonotically infects humans, causing a respiratory disease with severe pneumonia and death. With no approved antiviral or vaccine interventions for MERS, vaccines are being developed for camels to prevent virus transmission into humans. We have previously developed a chimpanzee adenoviral vector-based vaccine for MERS-CoV (ChAdOx1 MERS) and reported its strong humoral immunogenicity in dromedary camels. Here, we looked back at total RNA isolated from three immunised dromedaries pre and post-vaccination during the first day; and performed RNA sequencing and bioinformatic analysis in order to shed light on the molecular immune responses following a ChAdOx1 MERS vaccination. Our finding shows that a number of transcripts were differentially regulated as an effect of the vaccination, including genes that are involved in innate and adaptive immunity, such as type I and II interferon responses. The camel Bcl-3 and Bcl-6 transcripts were significantly upregulated, indicating a strong activation of Tfh cells, B cell, and NF-kB pathways. In conclusion, this study gives an overall view of the first changes in the immune transcriptome of dromedaries after vaccination; it supports the potency of ChAdOx1 MERS as a potential camel vaccine to block transmission and prevent new human cases and outbreaks.