Search

The advent of RNA-sequencing techniques has made it possible to generate large, unbiased gene expression datasets of tissues and cell types. Several studies describing gene expression data of microglia from Alzheimer's disease or multiple sclerosis have been published, aiming to generate more insight into the role of microglia in these neurological diseases. Though the raw sequencing data are often deposited in open access databases, the most accessible source of data for scientists is what is reported in published manuscripts. We observed a relatively limited overlap in reported differentially expressed genes between various microglia RNA-sequencing studies from multiple sclerosis or Alzheimer's diseases. It was clear that differences in experimental set up influenced the number of overlapping reported genes. However, even when the experimental set up was very similar, we observed that overlap in reported genes could be low. We identified that papers reporting large numbers of differentially expressed microglial genes generally showed higher overlap with other papers. In addition, though the pathology present within the tissue used for sequencing can greatly influence microglia gene expression, often the pathology present in samples used for sequencing was underreported, leaving it difficult to assess the data. Whereas reanalyzing every raw dataset could reduce the variation that contributes to the observed limited overlap in reported genes, this is not feasible for labs without (access to) bioinformatic expertise. In this study, we thus provide an overview of data present in manuscripts and their supplementary files and how these data can be interpreted.

Aim: Hepatocellular carcinoma (HCC), as one primary liver cancer type, accounts for 75%-85% of liver cancer cases. HCC is the second leading cause of cancer death in East Asia and sub-Saharan Africa and the sixth most common in western countries. Identification of key genes would facilitate the development of therapies and improve the prognosis outcomes of HCC patients. This study was to identify the key biological processes, pathways, and key genes in HCC. Methods: Data were downloaded from Broad GDAC. Differentially expressed genes (DEGs) and weighted gene coexpression network (WGCNA) were analyzed by DESeq2 and WGCNA, respectively. Gene ontology (GO) and KEGG enrichment analyses were performed on all DEGs and the coexpressed genes in two significant modules. Kaplan-Meier plotter online database was used to identify the potential prognostic genes in HCC. Finally, GEO database was used to validate the analysis of gene expression of Broad GDAC data. Results: The authors identified the dark gray and red modules as the significant modules in HCC based on WGCNA. GO and KEGG enrichment of the two significant modules identified the mitochondrion-mediated metabolic processes and pathways, and the cell cycle as the key biological processes and pathways in HCC. Subsequently, 28 hub genes were screened out by constructing protein-protein interaction network using Metascape. Finally, three genes (NDUFAF6, CKAP5, and DSN1 genes) were identified to be potential prognostic and key genes in HCC based on Kaplan-Meier survival analysis, GEO dataset validation, and literature review. Conclusions: The authors found that mitochondrion-mediated metabolic processes and the cell cycle were the key biological processes and pathways in HCC. NDUFAF6, CKAP5, and DSN1 genes were valuable genes with the potential to be prognosis biomarkers and targeted therapies in HCC.

Purpose. Combinations of expressed genes can discriminate radiation-exposed from normal control blood samples by machine learning based signatures (with 8 to 20% misclassification rates). These signatures can quantify therapeutically-relevant as well as accidental radiation exposures. The prodromal symptoms of Acute Radiation Syndrome (ARS) overlap those present in Influenza and Dengue Fever infections. Surprisingly, these human radiation signatures misclassified gene expression profiles of virally infected samples as false positive exposures. The present study investigates these and other confounders, and then mitigates their impact on signature accuracy. Methods. This study investigated recall by previous and novel radiation signatures independently derived from multiple Gene Expression Omnibus datasets on common and rare non-malignant blood disorders and blood-borne infections (thromboembolism, S. aureus bacteremia, malaria, sickle cell disease, polycythemia vera, and aplastic anemia). Normalized expression levels of signature genes are used as input to machine learning-based classifiers to predict radiation exposure in other hematological conditions. Results. Except for aplastic anemia, these blood-borne disorders modify the normal baseline expression values of genes present in radiation signatures, leading to false-positive misclassification of radiation exposures in 8 to 54% of individuals. Shared changes, predominantly in DNA damage response and apoptosis-related gene transcripts in radiation and confounding hematological conditions, compromise the utility of these signatures for radiation assessment. These confounding conditions (sickle cell disease, thromboembolism, S. aureus bacteremia, malaria) induce neutrophil extracellular traps, initiated by chromatin decondensation, DNA damage response and fragmentation followed by programmed cell death. Riboviral infections (for example, Influenza or Dengue fever) have been proposed to bind and deplete host RNA binding proteins, inducing R-loops in chromatin. R-loops that collide with incoming replication forks can result in incompletely repaired DNA damage, inducing apoptosis and releasing mature virus. To mitigate the effects of confounders, we evaluated predicted radiation-positive samples with novel gene expression signatures derived from radiation-responsive transcripts encoding secreted blood plasma proteins whose expression levels are unperturbed by these conditions. Conclusions. This approach identifies and eliminates misclassified samples with underlying hematological or infectious conditions, leaving only samples with true radiation exposures. Diagnostic accuracy is significantly improved by selecting genes that maximize both sensitivity and specificity in the appropriate tissue using combinations of the best signatures for each of these classes of signatures.

MicroRNAs (miRNAs) are single-stranded small non-coding RNA (18–25 nucleotides) that until a few years ago were considered junk RNA. In the last twenty years, they have acquired more importance thanks to the understanding of their influence on gene expression and their role as negative regulators at post-transcriptional level, influencing the stability of messenger RNA (mRNA). Approximately 5% of the genome encodes miRNAs which are responsible for regulating numerous signaling pathways, cellular processes and cell-to-cell communication. In the cardiovascular system, miRNAs control the functions of various cells, such as cardiomyocytes, endothelial cells, smooth muscle cells and fibroblasts, playing a role in physiological and pathological processes and seeming also related to variations in contractility and hereditary cardiomyopathies. They provide a new perspective on the pathophysiology of disorders such as hypertrophy, fibrosis, arrhythmia, inflammation and atherosclerosis. MiRNAs are differentially expressed in diseased tissue and can be released into the circulation and then detected. MiRNAs have become interesting for the development of new diagnostic and therapeutic tools for various diseases, including heart disease. In this review, the concept of miRNAs and their role in cardiomyopathies will be introduced, focusing on their potential as therapeutic and diagnostic targets (as biomarkers).

COVID-19 is a serious infectious disease that has recently swept the world, and research on its causative virus, SARS-CoV-2, remains insufficient. Therefore, this study uses bioinformatics analysis techniques to explore the human digestive tract diseases that may be caused by SARS-CoV-2 infection. The gene expression profile data set, numbered GSE149312, is from the Gene Expression Omnibus (GEO) database and is divided into a 24-h group and a 60-h group. R software is used to analyze and screen out differentially expressed genes (DEGs) and then gene ontology (GO) term and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses are performed. In KEGG, the pathway of non-alcoholic fatty liver disease exists in both the 24-h group and 60-h group. STRING is used to establish a protein–protein interaction (PPI) network, and Cytoscape is then used to visualize the PPI and define the top 12 genes of the node as the hub genes. Through verification, nine statistically significant hub genes are identified: AKT1, TIMP1, NOTCH, CCNA2, RRM2, TTK, BUB1B, KIF20A, and PLK1. In conclusion, the results of this study can provide a certain direction and basis for follow-up studies of SARS-CoV-2 infection of the human digestive tract and provide new insights for the prevention and treatment of diseases caused by SARS-CoV-2. Supplementary Information The online version contains supplementary material available at 10.1007/s10528-021-10144-w.

Obesity is a complex disease with multifactorial causes, and its prevalence is becoming a serious health crisis. For this reason, there is a crucial need to identify novel targets and players. With this aim in mind, we analyzed via RNA-sequencing the subcutaneous adipose tissue of normal weight and obesity-affected women, highlighting the differential expression in the two tissues. We specifically focused on long non-coding RNAs, as 6 of these emerged as dysregulated in the diseased-tissue (COL4A2-AS2, RPS21-AS, PELATON, ITGB2-AS1, ACER2-AS and CTEPHA1). For each of them, we performed both a thorough in silico dissection and in vitro validation, to predict their function during adipogenesis. We report the lncRNAs expression during adipose derived stem cells differentiation to adipocytes as model of adipogenesis and their potential modulation by adipogenesis-related transcription factors (C/EBPs and PPARγ). Moreover, inhibiting CTEPHA1 expression we investigated its impact on adipogenesis-related transcription factors, showing its significative dysregulation of C/EBPα expression. Lastly, we dissected the subcellular localization, pathway involvement and disease-correlation for coding differentially expressed genes. Together, these findings highlight a transcriptional deregulation at the basis of obesity, impacted by both coding and long non-coding RNAs.

Glioblastoma is the most aggressive type of brain tumors resistant to a number of antitumor drugs. The problem of therapy and drug treatment course is complicated by extremely high heterogeneity in the benign cell populations, the random arrangement of tumor cells, and polymorphism of their nuclei. The pathogenesis of gliomas needs to be studied using modern cellular technologies, genome- and transcriptome-wide technologies of high-throughput sequencing, analysis of gene expression on microarrays, and methods of modern bioinformatics to find new therapy targets. Functional annotation of genes related to the disease could be retrieved based on genetic databases and cross-validated by integrating complementary experimental data. Gene network reconstruction for a set of genes (proteins) proved to be effective approach to study mechanisms underlying disease progression. We used online bioinformatics tools for annotation of gene list for glioma, reconstruction of gene network and comparative analysis of gene ontology categories. The available tools and the databases for glioblastoma gene analysis are discussed together with the recent progress in this field.

Inhaled chemical/material exposures are a ubiquitous part of daily life around the world. There is a need to evaluate potential adverse effects of both single and repeat exposures for thousands of chemicals and an exponentially larger number of exposure scenarios (eg, repeated exposures). Meeting this challenge will require the development and use of in vitro new approach methodologies (NAMs); however, 2 major challenges face the deployment of NAMs in risk assessment are (1) characterizing what apical outcome(s) acute assays inform regarding the trajectory to long-term events, especially under repeated exposure conditions, and (2) capturing interindividual variability as it informs considerations of potentially susceptible and/or vulnerable populations. To address these questions, we used a primary human bronchial epithelial cell air-liquid interface model exposed to ozone (O3), a model oxidant and ubiquitous environmental chemical. Here we report that O3-induced proinflammatory gene induction is attenuated in repeated exposures thus demonstrating that single acute exposure outcomes do not reliably represent the trajectory of responses after repeated or chronic exposures. Further, we observed 10.1-, 10.3-, 14.2-, and 7-fold ranges of induction of interleukin (IL)-8, IL-6, heme oxygenase 1, and cyclooxygenase 2 transcripts, respectively, within in our population of 25 unique donors. Calculation of sample size estimates that indicated that 27, 24, 299, and 13 donors would be required to significantly power similar in vitro studies to identify a 2-fold change in IL-8, IL-6, HMOX1, and cyclooxygenase 2 transcript induction, respectively, to inform considerations of the uncertainty factors to reflect variability within the human population for in vitro studies.

Objectives: This study aimed to identify specific diagnostic as well as predictive targets of primary myelofibrosis (PMF). Methods: The gene expression profiles of GSE26049 were obtained from Gene Expression Omnibus (GEO) dataset, WGCNA was constructed to identify the most related module of PMF. Subsequently, Gene Ontology (GO), Kyoto Encyclopedia Genes and Genomes (KEGG), Gene Set Enrichment Analysis (GSEA) and Protein-Protein interaction (PPI) network were conducted to fully understand the detailed information of the interested green module. Machine learning, Principal component analysis (PCA), and expression pattern analysis including immunohistochemistry and immunofluorescence of genes and proteins were performed to validate the reliability of these hub genes. Results: Green module was strongly correlated with PMF disease after WGCNA analysis. 20 genes in green module were identified as hub genes responsible for the progression of PMF. GO, KEGG revealed that these hub genes were primarily enriched in erythrocyte differentiation, transcription factor binding, hemoglobin complex, transcription factor complex and cell cycle, etc. Among them, EPB42, CALR, SLC4A1 and MPL had the most correlations with PMF. Machine learning, Principal component analysis (PCA), and expression pattern analysis proved the results in this study. Conclusions: EPB42, CALR, SLC4A1 and MPL were significantly highly expressed in PMF samples. These four genes may be considered as candidate prognostic biomarkers and potential therapeutic targets for early stage of PMF. The effects are worth expected whether in the diagnosis at early stage or as therapeutic target.

Background Liver cancer is a common malignant tumor in China, with high mortality. Its occurrence and development were thoroughly studied by high-throughput expression microarray, which produced abundant data on gene expression, mRNA quantification and the clinical data of liver cancer. However, the hub genes, which can be served as biomarkers for diagnosis and treatment of early liver cancer, are not well screened. Results Here we present a new method for getting 6 key genes, aiming to diagnose and treat the early liver cancer. We firstly analyzed the different expression microarrays based on TCGA database, and a total of 1564 differentially expressed genes were obtained, of which 1400 were up-regulated and 164 were down-regulated. Furthermore, these differentially expressed genes were studied by using GO and KEGG enrichment analysis, a PPI network was constructed based on the STRING database, and 15 hub genes were obtained. Finally, 15 hub genes were verified by applying the survival analysis method on Oncomine database, and 6 key genes were ultimately identified, including PLK1, CDC20, CCNB2, BUB1, MAD2L1 and CCNA2. The robustness analysis of four independent data sets verifies the accuracy of the key gene’s classification of the data set. Conclusions Although there are complicated differences between cancer and normal cells in gene functions, cancer cells could be differentiated in case that a group of special genes expresses abnormally. Here we presented a new method to identify the 6 key genes for diagnosis and treatment of early liver cancer, and these key genes can help us understand the pathogenesis of liver cancer more deeply.

Aneurysmal subarachnoid hemorrhage (aSAH) is a life-threatening medical condition with a high mortality and disability rate. aSAH has an unclear pathogenesis, and limited treatment options are available. Here, we aimed to identify critical genes involved in aSAH pathogenesis using peripheral blood gene expression data of 43 patients with aSAH due to ruptured intracranial aneurysms and 18 controls with headache, downloaded from Gene Expression Omnibus. These data were used to construct a co-expression network using weighted gene co-expression network analysis (WGCNA). The biological functions of the hub genes were explored, and critical genes were selected by combining with differentially expressed genes analysis. Fourteen modules were identified by WGCNA. Among those modules, red, blue, brown and cyan modules were closely associated with aSAH. Moreover, 364 hub genes in the significant modules were found to play important roles in aSAH. Biological function analysis suggested that protein biosynthesis-related processes and inflammatory responses-related processes were involved in the pathology of aSAH pathology. Combined with differentially expressed genes analysis and validation in 35 clinical samples, seven gene (CD27, ANXA3, ACSL1, PGLYRP1, ALPL, ARG1, and TPST1) were identified as potential biomarkers for aSAH, and three genes (ANXA3, ALPL, and ARG1) were changed with disease development, that may provide new insights into potential molecular mechanisms for aSAH.

The role of B7‐H3 in acute myeloid leukemia (AML) is not fully understood. Two previous studies investigating its expression and significances in AML are partially different. In this study, we aimed to systematically characterize the genomic and immune landscape in AML patients with altered B7‐H3 expression using multi‐omics data in the public domain. We found significantly increased B7‐H3 expression in AML compared to either other hematological malignancies or healthy controls. Clinically, high B7‐H3 expression was associated with old age, TP53 mutations, wild‐type WT1 and CEBPA, and the M3 and M5 FAB subtypes. Moreover, we observed that increased B7‐H3 expression correlated significantly with a poor outcome of AML patients in four independent datasets. Gene set enrichment analysis (GSEA) revealed the enrichment of the “EMT” oncogenic gene signatures in high B7‐H3 expressers. Further investigation suggested that B7‐H3 was more likely to be associated with immune‐suppressive cells (macrophages, neutrophils, dendritic cells, and Th17 cells). B7‐H3 was also positively associated with a number of checkpoint genes, such as VISTA (B7‐H5), CD80 (B7‐1), CD86 (B7‐2), and CD70. In summary, we uncovered distinct genomic and immunologic features associated with B7‐H3 expression in AML. This may lead to a better understanding of the molecular mechanisms underlying B7‐H3 dysregulation in AML and to the development of novel therapeutic strategies., In this paper, we found significantly increased B7‐H3 expression in AML compared to either other hematological malignancies or healthy controls. Importantly, we observed that increased B7‐H3 expression correlated significantly with a poor outcome of AML patients in four independent datasets. We also uncovered distinct genomic and immunologic features associated with B7‐H3 expression in AML.

BACKGROUND Intracranial aneurysm (IA) rupture leads to deadly subarachnoid hemorrhages. However, the mechanisms leading to rupture remain poorly understood. Altered gene expression within IA tissue is linked to the pathobiology of aneurysm development and progression. Here, we analyzed expression patterns of control tissue samples and compared them to those of unruptured and ruptured IA tissue samples using data from the Gene Expression Omnibus (GEO). METHODS FASTQ files for 21 ruptured IAs, 21 unruptured IAs, and 16 control tissue samples were accessed from the GEO database. DESeq2 was used for differential expression analysis in three comparisons: unruptured IA versus control, ruptured IA versus control, and ruptured versus unruptured IA. Genes that were differentially expressed in multiple comparisons were evaluated to find those progressively increasing/decreasing from control to unruptured to ruptured. Significance was tested by either analysis of variance/Gabriel or Brown-Forsythe/Games Howell (p

The increasing prevalence of nonalcoholic fatty liver disease (NAFLD) as multifactorial chronic liver disease and the lack of a specific treatment have begun a new era in its treatment using gene expression changes and microRNAs. This study aimed to investigate the potential therapeutic effects of natural compounds in NAFLD by regulating miRNA expression. MicroRNAs play essential roles in regulating the cell’s biological processes, such as apoptosis, migration, lipid metabolism, insulin resistance, and adipocyte differentiation, by controlling the posttranscriptional gene expression level. The impact of current NAFLD pharmacological management, including drug and biological therapies, is uncertain. In this context, various dietary fruits or medicinal herbal sources have received worldwide attention versus NAFLD development. Natural ingredients such as berberine, lychee pulp, grape seed, and rosemary possess protective and therapeutic effects against NAFLD by modifying the gene’s expression and noncoding RNAs, especially miRNAs.

Neuroprotection in acute stroke has not been successfully translated from animals to humans. Animal research on promising agents continues largely in rats and mice which are commonly available to researchers. However, controversies continue on the most suitable species to model the human situation. Generally, putative agents seem less effective in mice as compared with rats. We hypothesized that this may be due to inter-species differences in stroke response and that this might be manifest at a genetic level. Here we used whole-genome microarrays to examine the differential gene regulation in the ischemic penumbra of mice and rats at 2 and 6 h after permanent middle cerebral artery occlusion (pMCAO; Raw microarray CEL data files are available in the GEO database with an accession number GSE163654). Differentially expressed genes (adj. p ≤ 0.05) were organized by hierarchical clustering, correlation plots, Venn diagrams and pathway analyses in each species and at each time-point. Emphasis was placed on genes already known to be associated with stroke, including validation by RT-PCR. Gene expression patterns in the ischemic penumbra differed strikingly between the species at both 2 h and 6 h. Nearly 90% of significantly regulated genes and most pathways modulated by ischemia differed between mice and rats. These differences were evident globally, among stroke-associated genes, immediate early genes, genes implicated in stress response, inflammation, neuroprotection, ion channels, and signal transduction. The findings of this study may have significant implications for the choice of species for screening putative stroke therapies.

Chromatin interactions play important roles in regulating gene expression. However, the availability of genome-wide chromatin interaction data is limited. We develop a computational method, chromatin interaction neural network (ChINN), to predict chromatin interactions between open chromatin regions using only DNA sequences. ChINN predicts CTCF- and RNA polymerase II-associated and Hi-C chromatin interactions. ChINN shows good across-sample performances and captures various sequence features for chromatin interaction prediction. We apply ChINN to 6 chronic lymphocytic leukemia (CLL) patient samples and a published cohort of 84 CLL open chromatin samples. Our results demonstrate extensive heterogeneity in chromatin interactions among CLL patient samples. Ministry of Education (MOE) National Research Foundation (NRF) Published version This research is supported by the National Research Foundation (NRF) Singapore through an NRF Fellowship awarded to M.J.F (NRF-NRFF2012-054) and NTU start-up funds awarded to M.J.F. This research is supported by the RNA Biology Center at the Cancer Science Institute of Singapore, NUS, as part of funding under the Singapore Ministry of Education Academic Research Fund Tier 3 awarded to Daniel Tenen as lead PI with M.J.F as co-investigator (MOE2014-T3-1-006). This research is supported by a National Research Foundation Competitive Research Programme grant awarded to V.T. as lead PI and M.J.F. as co-PI (NRF-CRP17-2017-02). This research is supported by the National Research Foundation Singapore and the Singapore Ministry of Education under its Research Centres of Excellence initiative. This research is supported by a Ministry of Education Tier II grant awarded to M.J.F (T2EP30120-0020).

Background Lung adenocarcinoma (LUAD), as the most common subtype of lung cancer, is the leading cause of cancer deaths in the world. The accumulation of driver gene mutations enables cancer cells to gradually acquire growth advantage. Therefore, it is important to understand the functions and interactions of driver gene mutations in cancer progression. Methods We obtained gene mutation data and gene expression profile of 506 LUAD tumors from The Cancer Genome Atlas (TCGA). The subtypes of tumors with driver gene mutations were identified by consensus cluster analysis. Results We found 21 significantly mutually exclusive pairs consisting of 20 genes among 506 LUAD patients. Because of the increased transcriptomic heterogeneity of mutations, we identified subtypes among tumors with non‐silent mutations in driver genes. There were 494 mutually exclusive pairs found among driver gene mutations within different subtypes. Furthermore, we identified functions of mutually exclusive pairs based on the hypothesis of functional redundancy of mutual exclusivity. These mutually exclusive pairs were significantly enriched in nuclear division and humoral immune response, which played crucial roles in cancer initiation and progression. We also found 79 mutually exclusive triples among subtypes of tumors with driver gene mutations, which were key roles in cell motility and cellular chemical homeostasis. In addition, two mutually exclusive triples and one mutually exclusive triple were associated with the overall survival and disease‐specific survival of LUAD patients, respectively. Conclusions We revealed novel mutual exclusivity and generated a comprehensive functional landscape of driver gene mutations, which could offer a new perspective to understand the mechanisms of cancer development and identify potential biomarkers for LUAD therapy., We revealed novel mutual exclusivity and generated a comprehensive functional landscape of driver gene mutations, which could offer a new perspective to understand the mechanisms of cancer development and identify biomarkers for LUAD therapy.

Background Metazoan cells only utilize a small subset of the potential DNA replication origins to duplicate the whole genome in each cell cycle. Origin choice is linked to cell growth, differentiation, and replication stress. Although various genetic and epigenetic signatures have been linked to the replication efficiency of origins, there is no consensus on how the selection of origins is determined. Results We apply dual-color stochastic optical reconstruction microscopy (STORM) super-resolution imaging to map the spatial distribution of origins within individual topologically associating domains (TADs). We find that multiple replication origins initiate separately at the spatial boundary of a TAD at the beginning of the S phase. Intriguingly, while both high-efficiency and low-efficiency origins are distributed homogeneously in the TAD during the G1 phase, high-efficiency origins relocate to the TAD periphery before the S phase. Origin relocalization is dependent on both transcription and CTCF-mediated chromatin structure. Further, we observe that the replication machinery protein PCNA forms immobile clusters around TADs at the G1/S transition, explaining why origins at the TAD periphery are preferentially fired. Conclusion Our work reveals a new origin selection mechanism that the replication efficiency of origins is determined by their physical distribution in the chromatin domain, which undergoes a transcription-dependent structural re-organization process. Our model explains the complex links between replication origin efficiency and many genetic and epigenetic signatures that mark active transcription. The coordination between DNA replication, transcription, and chromatin organization inside individual TADs also provides new insights into the biological functions of sub-domain chromatin structural dynamics.

The type 2 deiodinase (D2) in the neonatal liver accelerates local thyroid hormone triiodothyronine (T3) production and expression of T3-responsive genes. Here we show that this surge in T3 permanently modifies hepatic gene expression. Liver-specific Dio2 inactivation (Alb-D2KO) transiently increases H3K9me3 levels during post-natal days 1–5 (P1–P5), and results in methylation of 1,508 DNA sites (H-sites) in the adult mouse liver. These sites are associated with 1,551 areas of reduced chromatin accessibility (RCA) within core promoters and 2,426 within intergenic regions, with reduction in the expression of 1,363 genes. There is strong spatial correlation between density of H-sites and RCA sites. Chromosome conformation capture (Hi-C) data reveals a set of 81 repressed genes with a promoter RCA in contact with an intergenic RCA ~300 Kbp apart, within the same topologically associating domain (χ2 = 777; p, Whether thyroid hormones affect gene expression via DNA methylation is not well known. Here the authors show that type 2 deiodinase (D2) converts T4 to produce T3, which prevents DNA methylation of discrete areas in the neonatal liver. In the absence of D2, DNA methylation occurs and is associated with reduced chromatin accessibility in promoters and enhancers and affects gene expression.

Friedreich’s ataxia (FA) is a neurodegenerative disease caused by the epigenetic repression of the Frataxin gene modulating mitochondrial activity in the brain, which has a diffuse phenotypic impact on patients’ motor behavior. Therefore, with current gold-standard clinical scales, it requires 18–24 month-long clinical trials to determine if disease-modifying therapies are at all beneficial. Our high-performance monitoring approach captures the full-movement kinematics from human subjects using wearable body sensor networks from a cohort of FA patients during their regular clinical visits. We then use artificial intelligence to convert these movement data using universal behavior fingerprints into a digital biomarker of disease state. This enables us to predict two different ‘gold-standard’ clinical scores (SCAFI, SARA) that serve as primary clinical endpoints. Crucially, by performing gene expression analysis on each patient their personal Frataxin gene expression levels were poorly, if at all, correlated with their clinical scores – fundamentally failing to establish a link between disease mechanism (dysregulated gene expression) and measures to quantify it in the behavioral phenotype. In contrast, our wearable digital biomarker can accurately predict for each patient their personal FXN gene expression levels, demonstrating the sensitivity of our approach and the importance of FXN levels in FA. Therefore, our data-derived biomarker approach can not only cross-sectionally predict disease and their gene expression levels but also their longitudinal disease trajectory: it is sensitive and accurate enough to detect disease progression with much fewer subjects or shorter time scales than existing primary endpoints. Our work demonstrates that data-derived wearable biomarkers have the potential to substantially reduce clinical trial durations and a first in-human demonstration of reconstructing FXN gene expression levels from behavioral data alone.

To identify novel prognostic and therapeutic targets for osteosarcoma patients, we compared the gene expression profiles of osteosarcoma and control tissues from the GSE42352 dataset in the Gene Expression Omnibus. Differentially expressed genes were subjected to Gene Ontology, Kyoto Encyclopedia of Genes and Genomes, Gene Set Enrichment and protein-protein interaction network analyses. Survival curve analyses indicated that osteosarcoma patients with lower mRNA levels of cyclin-dependent kinase 1 (CDK1) and topoisomerase II alpha had better prognoses. Various computer-aided techniques were used to identify potential CDK1 inhibitors for osteosarcoma patients, and PHA-793887 was predicted to be a safe drug with a high binding affinity for CDK1. In vitro, MTT and colony formation assays demonstrated that PHA-793887 reduced the viability and clonogenicity of osteosarcoma cells, while a scratch assay suggested that PHA-793887 impaired the migration of these cells. Flow cytometry experiments revealed that PHA-793887 dose-dependently induced apoptosis in osteosarcoma cells. Western blotting and enzyme-linked immunosorbent assays indicated that CDK1 expression in osteosarcoma cells declined with increasing PHA-793887 concentrations. These results suggest that PHA-793887 could be a promising new treatment for osteosarcoma.

Objectives. Chronic obstructive pulmonary disease (COPD) is characterized by lung inflammation and remodeling. Macrophage polarization is associated with inflammation and tissue remodeling, as well as immunity. Therefore, this study attempts to investigate the diagnostic value and regulatory mechanism of macrophage polarization-related genes for COPD by bioinformatics analysis and to provide a new theoretical basis for experimental research. Methods. The raw gene expression profile dataset (GSE124180) was collected from the Gene Expression Omnibus (GEO) database. Next, a weighted gene coexpression network analysis (WGCNA) was conducted to screen macrophage polarization-related genes. The differentially expressed genes (DEGs) between the COPD and normal samples were generated using DESeq2 v3.11 and overlapped with the macrophage polarization-related genes. Moreover, functional annotations of overlapped genes were conducted by Database for Annotation, Visualization and Integrated Discovery (DAVID) Bioinformatics Resource. The immune-related genes were selected, and their correlation with the differential immune cells was analyzed by Pearson. Finally, receiver operating characteristic (ROC) curves were used to verify the diagnostic value of genes. Results. A total of 4922 coexpressed genes related to macrophage polarization were overlapped with the 203 DEGs between the COPD and normal samples, obtaining 25 genes related to COPD and macrophage polarization. GEM, S100B, and GZMA of them participated in the immune response, which were considered the candidate biomarkers. GEM and S100B were significantly correlated with marker genes of B cells which had a significant difference between the COPD and normal samples. Moreover, GEM was highly associated with the genes in the PI3K/Akt/GSK3β signaling pathway, regulation of actin cytoskeleton, and calcium signaling pathway based on a Pearson correlation analysis of the candidate genes and the genes in the B cell receptor signaling pathway. PPI network analysis also indicated that GEM might participate in the regulation of the PI3K/Akt/GSK3β signaling pathway. The ROC curve showed that GEM possessed an excellent accuracy in distinguishing COPD from normal samples. Conclusions. The data provide a transcriptome-based evidence that GEM is related to COPD and macrophage polarization likely contributes to COPD diagnosis. At the same time, it is hoped that in-depth functional mining can provide new ideas for exploring the COPD pathogenesis.

Background In the present study, we aimed to explore whether common genetic targets or signaling pathways existed in chemical cystitis. Methods Gene Expression Omnibus (GEO) database was used to search the related gene expression profiles. The differentially expressed genes (DEGs) were identified by using GEO2R. The DAVID 6.8 Beta and R software were used to perform Kyoto Encyclopedia of Genes and Genomes pathway analysis and Gene Ontology function analysis of DEGs. The protein‐protein interaction network was constructed by STRING 11.0 to reveal the potential gene interactions. The expression of cyclin‐dependent kinase 1 (Cdk1) at the messnger RNA (mRNA) and protein levels was examined by real‐time polymerase chain reaction (PCR) and Western blot analysis analysis, respectively. Results The GEO database was searched, and the gene expression profiles of GSE55986 and GSE68539 were downloaded. A total of 262 DEGs and 356 DEGs were identified from GSE55986 and GSE68539, respectively. We found that the p53 signaling pathway might play a key role in the development of chemical cystitis, and Cdk1 acted as a crucial gene in the p53 signaling pathway. Moreover, the experimental results of real‐time PCR and Western blot analysis analysis demonstrated that the expression of Cdk1 at the mRNA and protein levels in cystitis tissues was significantly increased in different animal models of chemical cystitis compared with the control group. Conclusion Cdk1 might be a potential pathogenic genetic target for chemical cystitis., Cdk1 shows to be a potential pathogenic genetic target for chemical cystitis

Objective The aim of this study was to explore expression profiles of long noncoding RNA (lncRNA)‐messenger RNA (mRNA) in abdominal aortic aneurysm (AAA) patients. Further, we explored the mechanisms by which lncRNA SNHG5 modulates the function of vascular smooth muscle cells (VSMC) in AAA. Methods Human gene expression profile GSE57691 dataset, was retrieved from Gene Expression Omnibus database. The dataset included gene expression array data of 49 AAA patients and 10 control aortic specimens from organ donors. To explore the main roles of the biological network, differentially expressed lncRNA and mRNAs in the aortic aneurysm (AAA) and normal aortic specimens were determined. Differentially expressed lncRNA and mRNAs were then used to construct a competing endogenous RNA (ceRNA) network using Cytoscape software, and the five key lncRNA were identified. SNHG5 which was significantly downregulated in the AAA was chosen and analysis showed that it regulates mir‐205‐5p and SMAD4 by binding to mir‐205‐5p. Double luciferase reporter gene assays, RNA immunoprecipitation, and RNA knockdown studies were used to establish the relationship between SNHG5 and mir‐205‐5p. Apoptosis rate was determined using flow cytometry, whereas cell proliferation was evaluated using Edu, and 24 well Transwell assay. Western blot analysis was used to determine protein expression levels. Results The five differentially expressed lncRNAs were significantly correlated with 34 microRNAs and 112 mRNAs. mRNAs in the ceRNA network are implicated in protein binding, signal transduction, DNA and RNA transcription, development, and cell differentiation. SNHG5 was downregulated in the AAA and acts as a molecular sponge for mir‐205. Downregulation of SNHG5 induces expression of mir‐205‐5p. Increased mir‐205‐5p expression level inhibits SMAD4 production, thus inhibiting proliferation and migration and promotes apoptosis of smooth muscle cells. Conclusion Bioinformatics were used to explore molecular mechanism of AAA progression. The findings of this study show that lncRNA SNHG5 regulates proliferation and apoptosis of VSMC cells through modulation of the mir‐205‐5p/SMAD4 axis. Therefore, SNHG5 is a potential therapeutic target for AAA disease., The findings of this study show that downregulation of SNHG5 increases expression levels of mir‐205‐5p and inhibits SMAD4 expression, thus affecting the function of vascular smooth muscle in abdominal aortic aneurysm. Therefore, SNHG5 is an important factor in pathogenesis of abdominal aortic aneurysm and thus is potential therapeutic target for abdominal aortic aneurysm.

Background Malaria is a disease annually causing over 400,000 deaths. Deep understanding of molecular and genetic processes underlying its life cycle and pathogenicity is required to efficiently resist it. RNA interference is a mechanism of the gene expression regulation typical for a wide variety of species. Even though the existence of this phenomenon in Plasmodium falciparum has long been rejected, several recent works pose hypotheses and provide direct and indirect evidence of the existence of mechanisms similar to RNA interference in this organism. In particular, the possibility of regulation of P. falciparum gene expression through human microRNAs is of great importance both for fundamental biology and for medicine. In the present work we address the problem of possibility of the existence in the P. falciparum genome of the nucleotide sequences such that mRNAs transcribed from genes containing these sequences could form duplexes with human microRNAs. Using bioinformatics methods we have analysed genomes of 15 P. falciparum isolates for sequences homological to these microRNAs. Results The analysis has demonstrated the existence of a vast number of genes that could potentially be regulated by the human microRNAs in the plasmodial genome. Conclusions Despite the fact that the numbers of homological intervals vary significantly between isolates, the hsa-miR-451a and hsa-miR-223-3p microRNAs are expected to make the most notable contribution to the pathogenesis of P. falciparum malaria. The majority of homological intervals occur in genes encoding cell adhesion proteins.

A need exists for local (ie, bladder-specific) interventions to treat overactive bladder (OAB) with low risk of unwanted postprocedural outcomes. Gene therapy targeted to leverage endogenous physiology in bladder cells may assist in restoring normal cell and organ function. Herein, we review the potential promise of gene therapy for treating OAB, focusing on gene transfer of URO-902, a non-viral naked plasmid DNA expressing the big potassium (BK) channel. We searched PubMed for articles concerning functional aspects of the BK channel and its potential use for gene transfer as local OAB treatment. Results from preclinical, phase 1, and phase 2 studies of URO-902 for erectile dysfunction and phase 1 studies of URO-902 for OAB are included. The BK channel has been extensively studied; however, URO-902 is the first gene therapy used in clinical trials directed toward treating OAB via the BK channel. In both URO-902 studies, there were no serious adverse events considered treatment related and no adverse events leading to early withdrawal. Both studies included secondary efficacy endpoints with promising results suggesting improvement in OAB symptoms, and quality of life, with use of URO-902 versus placebo. Gene therapy involving the BK channel, such as gene transfer with URO-902, has demonstrated promising safety and efficacy results in women with OAB. Findings warrant further investigation of the use of URO-902 for OAB treatment.

Thyroid hormones (TH; T3 and T4) play a fundamental role in the fetal stage to the adult phase, controlling gene and protein expression in virtually all tissues. The endocrine and CNS systems have relevant interaction, and the TH are pivotal for the proper functioning of the CNS. A slight failure to regulate TH availability during pregnancy and/or childhood can lead to neurological disorders, for example, autism and cognitive impairment, or depression. In this review, we highlight how TH acts in controlling gene expression, its role in the CNS, and what substances widely found in the environment can cause in this tissue. We highlight the role of Endocrine Disruptors used on an everyday basis in the processing of mRNAs responsible for neurodevelopment. We conclude that TH, more precisely T3, acts mainly throughout its nuclear receptors, that the deficiency of this hormone, either due to the lack of its main substrate iodine, or by to incorrect organification of T4 and T3 in the gland, or by a mutation in transporters, receptors and deiodinases may cause mild (dysregulated mood in adulthood) to severe neurological impairment (Allan-Herndon-Dudley syndrome, presented as early as childhood); T3 is responsible for the expression of numerous CNS genes related to oxygen transport, growth factors, myelination, cell maturation. Substances present in the environment and widely used can interfere with the functioning of the thyroid gland, the action of TH, and the functioning of the CNS.

Hepatitis B virus (HBV) infection is a worldwide health problem with high morbidity and mortality rates. The therapeutic vaccine is a promising method of treatment, and HBV polymerase plays a vital role in viral replication. Therefore, a therapeutic vaccine that binds to HBV DNA polymerase may control HBV infection. We predicted and selected epitopes of polymerase using online databases and analysis software. We then performed molecular docking and peptide binding assays to evaluate the binding energies and affinities between polymerase epitopes and the HLA-A0201 molecule. Finally, we induced T cells from the peripheral blood mononuclear cells (PBMCs) of healthy donors using each epitope and quantified the functions of epitope-specific T cells by IFN-γELISPOT assay, T2 cell cytotoxicity assay, HepG2.2.15 cell cytotoxicity assay and HBV gene expression assays. Four epitopes (RVTGGVFLV, GLLGFAAPF, LLDDEAGPL and YMDDVVLGA) had low binding energy and two epitopes (RVTGGVFLV and GLLGFAAPF) had a high binding affinity. The T cells stimulated by two epitopes (GLLGFAAPF and HLYSHPIIL) had a greater ability to induce immune response and suppress HBV. The HBV DNA polymerase epitopes identified in this study are promising targets for designing an epitope-based therapeutic vaccine against HBV.

Background Improvement of radiotherapy efficacy requires better insight in the dynamic responses that occur during irradiation. Here, we aimed to identify the molecular responses that are triggered during clinically applied fractionated irradiation. Methods Gene expression analysis was performed by RNAseq or microarray analysis of cancer cells or xenograft tumors, respectively, subjected to 3–5 weeks of 5 × 2 Gy/week. Validation of altered gene expression was performed by qPCR and/or ELISA in multiple cancer cell lines as well as in pre- and on-treatment biopsies from esophageal cancer patients (NCT02072720). Targeted protein inhibition and CRISPR/Cas-induced gene knockout was used to analyze the role of type I interferons and cGAS/STING signaling pathway in the molecular and cellular response to fractionated irradiation. Results Gene expression analysis identified type I interferon signaling as the most significantly enriched biological process induced during fractionated irradiation. The commonality of this response was confirmed in all irradiated cell lines, the xenograft tumors and in biopsies from esophageal cancer patients. Time-course analyses demonstrated a peak in interferon-stimulated gene (ISG) expression within 2–3 weeks of treatment. The response was accompanied by a variable induction of predominantly interferon-beta and/or -lambda, but blocking these interferons did not affect ISG expression induction. The same was true for targeted inhibition of the upstream regulatory STING protein while knockout of STING expression only delayed the ISG expression induction. Conclusions Collectively, the presented data show that clinically applied fractionated low-dose irradiation can induce a delayed type I interferon response that occurs independently of interferon expression or STING signaling. These findings have implications for current efforts that aim to target the type I interferon response for cancer treatment.

Epigenetic mechanisms play an important part in the regulation of gene expression and these alterations may induce long-term changes in gene function and metabolism. They have received extensive attention in bridging the gap between environmental exposures and disease development via their influence on gene expression. DNA methylation is the earliest discovered epigenetic alteration. In this review, we try to examine the role of DNA methylation and histone modification in Age related macular degeneration (AMD) and Diabetic Retinopathy (DR), its vascular complications and recent progress. Given the complex nature of AMD and DR, it is crucial to improve therapeutics which will greatly enhance the quality of life and reduce the burden for millions of patients living with these potentially blinding conditions.

Brain tumors are associated with adverse outcomes despite improvements in radiation therapy, chemotherapy, and photodynamic therapy. However, treatment approaches are evolving, and new biological phenomena are being explored to identify the appropriate treatment of brain tumors. Long non-coding RNAs (lncRNAs), a type of non-coding RNA longer than 200 nucleotides, regulate gene expression at the transcriptional, post-transcriptional, and epigenetic levels and are involved in a variety of biological functions. Recent studies on lncRNAs have revealed their aberrant expression in various cancers, with distinct expression patterns associated with their instrumental roles in cancer. Abnormal expression of lncRNAs has also been identified in brain tumors. Here, we review the potential roles of lncRNAs and their biological functions in the context of brain tumors. We also summarize the current understanding of the molecular mechanisms and signaling pathways related to lncRNAs that may guide clinical trials for brain tumor therapy.

Background Malignant pleural mesothelioma (MPM) is a rare tumor in the pleura. This study was carried out to identify key genes and pathways that may be involved in MPM. Methods Microarray datasets GSE51024 and GSE2549 were analyzed for differentially expressed genes (DEGs) between normal and MPM tissues. The identified DEGs were subjected to functional analyses using bioinformatics tools. Results A total of 276 DEGs were identified, consisting of 187 downregulated and 79 upregulated genes. Gene ontology and Kyoto encyclopedia of genes and genomes pathway enrichment analysis indicated that the DEGs were enriched in extracellular structure organization, extracellular matrix, and ECM−receptor interaction. Due to high degree of connectivity among 24 hub genes, EZH2 and HMMR are likely to play roles in the carcinogenesis and progression of MPM. The two genes were found over-expressed in MPM tissues. Patients with elevated EZH2 and HMMR expressions had poor overall survival. Conclusions EZH2 and HMMR are identified to be the hub genes for MPM and they may be further characterized to better understand the molecular mechanisms underlying the carcinogenesis of MPM.

In the investigation of molecular mechanisms underlying cell state changes, a crucial analysis is to identify differentially expressed (DE) genes along a continuous cell trajectory, which can be estimated by pseudotime inference from single-cell RNA-sequencing (scRNA-seq) data. However, existing methods that identify DE genes based on inferred pseudotime do not account for the uncertainty in pseudotime inference. Also, they either have ill-posed p-values that hinder the control of false discovery rate (FDR) or have restrictive models that reduce the power of DE gene identification. To overcome these drawbacks, we propose PseudotimeDE, a robust method that accounts for the uncertainty in pseudotime inference and thus identifies DE genes along cell pseudotime with well-calibrated p-values. PseudotimeDE is flexible in allowing users to specify the pseudotime inference method and to choose the appropriate model for scRNA-seq data. Comprehensive simulations and real-data applications verify that PseudotimeDE provides well-calibrated p-values essential for controlling FDR and downstream analysis and that PseudotimeDE is more powerful than existing methods to identify DE genes.

Gene therapy can be used to restore cell function in monogenic disorders or to endow cells with new capabilities, such as improved killing of cancer cells, expression of suicide genes for controlled elimination of cell populations, or protection against chemotherapy or viral infection. While gene therapies were originally most often used to treat monogenic diseases and to improve hematopoietic stem cell transplantation outcome, the advent of genetically modified immune cell therapies, such as chimeric antigen receptor modified T cells, has contributed to the increased numbers of patients treated with gene and cell therapies. The advancement of gene therapy with integrating retroviral vectors continues to depend upon world-wide efforts. As the topic of this special issue is “Spotlight on Germany,” the goal of this review is to provide an overview of contributions to this field made by German clinical and research institutions. Research groups in Germany made, and continue to make, important contributions to the development of gene therapy, including design of vectors and transduction protocols for improved cell modification, methods to assess gene therapy vector efficacy and safety (e.g., clonal imbalance, insertion sites), as well as in the design and conduction of clinical gene therapy trials.

Polycystic ovary syndrome (PCOS) is a prevalent endocrine and metabolic disorder in women of childbearing age. Recent studies suggest important roles for lncRNAs in PCOS development. Based on the hypothesis that lncRNAs are able to regulate mRNA functions by competitive binding to shared miRNAs, the present work sought to construct a PCOS-related lncRNA-mRNA network (PCLMN) to identify key lncRNAs with dysregulated expression and potential prognostic and therapeutic relevance. A global background network was constructed after retrieving lncRNA-miRNA and miRNA-mRNA pairs from the lncRNASNP2, miRTarBase and StarBase databases. Based on gene expression profiles from ovarian granulosa cells from PCOS patients and controls in the GEO’s GSE95728 dataset, the PCLMN was then constructed by applying hypergeometric testing. Using topological analysis, we identified 3 lncRNAs (LINC00667, AC073172.1 and H19) ranking within the top-ten gene lists for all three centrality measures. We then explored their subcellular localization, performed functional module analyses, and identified 4 sex hormone-related transcription factors as potential regulators of their expression. Significant associations with inflammation, oxidative stress, and apoptosis-related processes and pathways were revealed for the key lncRNAs in our PCMLN. Further studies verifying the mRNA/lncRNA relationships identified herein are needed to clarify their clinical significance.

Kinetoplastid parasites are responsible for both human and animal diseases across the globe where they have a great impact on health and economic well-being. Many species and life cycle stages are difficult to study due to limitations in isolation and culture, as well as to their existence as heterogeneous populations in hosts and vectors. Single-cell transcriptomics (scRNA-seq) has the capacity to overcome many of these difficulties, and can be leveraged to disentangle heterogeneous populations, highlight genes crucial for propagation through the life cycle, and enable detailed analysis of host–parasite interactions. Here, we provide a review of studies that have applied scRNA-seq to protozoan parasites so far. In addition, we provide an overview of sample preparation and technology choice considerations when planning scRNA-seq experiments, as well as challenges faced when analysing the large amounts of data generated. Finally, we highlight areas of kinetoplastid research that could benefit from scRNA-seq technologies.

This study aimed to identify effective targets for carcinogenesis of primary myelofibrosis (PMF), as well as to screen ideal lead compounds with potential inhibition effect on Janus kinase 2 to contribute to the medication design and development. Gene expression profiles of GSE26049, GSE53482, GSE61629 were obtained from the Gene Expression Omnibus database. The differentially expressed genes were identified, and functional enrichment analyses such as Gene Ontology, protein-protein interaction network etc., were performed step by step. Subsequently, highly-precise computational techniques were conducted to identify potential inhibitors of JAK2. A series of structural biology methods including virtual screening, ADMET (absorption, distribution, metabolism, excretion, and toxicity) prediction, molecule docking, molecular dynamics simulation etc., were implemented to discover novel natural compounds. Results elucidated that PMF patients had abnormal LCN2, JAK2, MMP8, CAMP, DEFA4, LTF, MPO, HBD, STAT4, EBF1 mRNA expression compared to normal patients. Functional enrichment analysis revealed that these genes were mainly enriched in erythrocyte differentiation, neutrophil degranulation and killing cells of other organisms. Two novel natural compounds, ZINC000013513540 and ZINC000004099068 were found binding to JAK2 with favorable interaction energy together with high binding affinity. They were predicted with non-Ames mutagenicity, low-rodent carcinogenicity, less developmental toxicity potential as well as non-toxicity with liver. Molecular dynamics simulation demonstrated that these two complexes: ZINC000013513540-JAK2 and ZINC000004099068-JAK2 could exist stably under natural circumstances. In conclusion, this study revealed hub genes in the carcinogenesis of PMF. ZINC000013513540 and ZINC000004099068 were promising drugs in dealing with PMF. This study may also accelerate exploration of new drugs.

Background Osteoarthritis (OA) is a common skeletal system disease that has been partially attributed to genetic factors. The hand is frequently affected, which seriously affects the patient’s quality of life. However, the pathogenetic mechanism of hand osteoarthritis (hand OA) is still elusive. Methods A genome-wide association study (GWAS) summary of hand OA was obtained from the UK Biobank dataset, which contains data from a total of 452,264 White British individuals, including 37,782 OA patients. The transcriptome-wide association study (TWAS) of hand OA was performed using FUnctional Summary-based ImputatiON (FUSION) with the skeletal muscle and blood as gene expression references. The significant genes identified by TWAS were further subjected to gene set enrichment analysis (GSEA) with the Database for Annotation, Visualization and Integrated Discovery (DAVID) tool. Furthermore, we compared the genes and gene sets identified by our TWAS with that of a knee OA mRNA expression profile to detect the genes and gene sets shared by TWAS and mRNA expression profiles in OA. The mRNA expression profiles of 18 normal knee cartilages and 20 OA knee cartilages were acquired from the Gene Expression Omnibus database (accession number: GSE114007). Results TWAS identified 177 genes with P P = 0.0001), RIC3 (P = 0.0003), and AC005154.6 (P = 0.0004). TWAS identified 423 genes with P P = 0.0002), ZNF880 (P = 0.0002), and NCKIPSD (P = 0.0003). After comparing the results of the TWAS to those of the mRNA expression profiling of OA, we identified 5 common genes, including DHRS3 (log2fold = − 1.85, P = 3.31 × 10− 9) and SKP2 (log2fold = 1.36, P = 1.62 × 10− 8). GSEA of TWAS identified 51 gene ontology (GO) terms for hand OA, for example, protein binding (P = 0.0003) and cytosol (P = 0.0020). We also detected 6 common GO terms shared by TWAS and mRNA expression profiling, including protein binding (PTWAS = 2.54 × 10− 4, PmRNA = 3.42 × 10− 8), extracellular exosome (PTWAS = 0.02, PmRNA = 1.18 × 10− 4), and cytoplasm (PTWAS = 0.0183, PmRNA = 0.0048). Conclusion In this study, we identified 5 candidate genes and 6 GO terms related to hand OA, which may help to uncover the pathogenesis of hand OA. It should be noted that the possible difference in the gene expression profiles between hand OA and knee OA may affect our study results, which should be interpreted with caution.

Guang-zhao Huang,1,* Qing-qing Wu,1,* Ze-nan Zheng,1,* Ting-ru Shao,1 Fei Li,1 Xin-yan Lu,1 Heng-yu Ye,1 Gao-xiang Chen,1 Yu-xing Song,1 Wei-sen Zeng,2 Yi-long Ai,3 Xiao-zhi Lv1 1Department of Oral & Maxillofacial Surgery, NanFang Hospital, Southern Medical University, Guangzhou, Guangdong Province, People’s Republic of China; 2Department of Cell Biology, School of Basic Medical Science, Southern Medical University, Guangzhou, Guangdong Province, People’s Republic of China; 3Foshan Stomatological Hospital, School of Stomatology and Medicine, Foshan University, Foshan, Guangdong Province, People’s Republic of China*These authors contributed equally to this workCorrespondence: Xiao-zhi LvDepartment of Oral & Maxillofacial Surgery, NanFang Hospital, Southern Medical University, Baiyun District, Guangzhou, Guangdong Province, 510080, People’s Republic of ChinaTel +86 13724811658Email lxzsurgeon@126.comYi-long AiFoshan Stomatological Hospital, School of Stomatology and Medicine, Foshan University, Chancheng District, Foshan, Guangdong Province, 528000, People’s Republic of ChinaTel +86 13827755774Email aiyilong@126.comPurpose: Plenty of studies showed that the immune system was associated with cancer initiation and progression. This study aimed to explore the prognostic biomarkers from immune-related genes (IRGs) in oral squamous cell carcinoma (OSCC).Materials and Methods: RNA-seq data were downloaded from The Cancer Genome Atlas (TCGA) and IRGs and transcription factors (TFs) were extracted. Then, the co-expression network between IRGs and TFs was constructed using the “WGCNA” package in R software. Furthermore, a gene expression signature according to IRGs was constructed to predict OSCC prognosis and its accuracy was validated by survival analysis. Subsequently, correlation analyses between risk-score and immune cells level and clinical parameters were performed. Finally, immune-related biomarkers were selected and further investigated using gain-of-function assays in vitro.Results: A total of 32 normal cases and 317 OSCC cases were selected in our study. Differentially-expressed analysis indicated that there were 381 differentially-expressed IRGs and 62 TFs in OSCC. Among them, 25 TFs and 21 IRGs were enrolled in the co-expression network. Furthermore, we found that gene expression signature on the basis of 10 IRGs could predict the prognosis accurately and a high-risk score based on gene expression signature meant a high T classification, terminal clinical stage, and low immune cells level in OSCC. Finally, cathepsin G (CTSG) was identified as a potential immune-related biomarker and therapeutic target in OSCC.Conclusion: In conclusion, IRGs were directly involved in the development and progression of OSCC. Furthermore, CTSG was identified as a potential independent biomarker and might be an immunotherapeutic target in OSCC treatment.Keywords: oral squamous cell carcinoma, immune, biomarker, bioinformatics, CTSG

Pediatric therapy-related myeloid neoplasms (tMN) occur in children after exposure to cytotoxic therapy and have a dismal prognosis. The somatic and germline genomic alterations that drive these myeloid neoplasms in children and how they arise have yet to be comprehensively described. We use whole exome, whole genome, and/or RNA sequencing to characterize the genomic profile of 84 pediatric tMN cases (tMDS: n = 28, tAML: n = 56). Our data show that Ras/MAPK pathway mutations, alterations in RUNX1 or TP53, and KMT2A rearrangements are frequent somatic drivers, and we identify cases with aberrant MECOM expression secondary to enhancer hijacking. Unlike adults with tMN, we find no evidence of pre-existing minor tMN clones (including those with TP53 mutations), but rather the majority of cases are unrelated clones arising as a consequence of cytotoxic therapy. These studies also uncover rare cases of lineage switch disease rather than true secondary neoplasms., Paediatric therapy-related myeloid neoplasms (tMN) have a dismal prognosis and have not been comprehensively profiled. Here the authors characterise the molecular landscape of 84 paediatric tMN patients, and find that, unlike adult tMNs, these do not emerge from pre-existing clones and that MECOM dysregulation is frequent.