Whole-transcriptome Sequencing Research Articles

Co-regulated ceRNA network mediated by circRNA and lncRNA in patients with gouty arthritis

Numerous studies have demonstrated the involvement of messenger RNAs (mRNAs) and non-coding RNAs, including long non-coding RNAs (lncRNA), circular RNAs (circRNAs) and microRNA (miRNAs), in gouty arthritis onset; however, the regulatory mechanism has not yet been elucidated. Here, we applied whole-transcriptome sequencing to identify the differentially expressed circRNAs, lncRNAs, miRNAs and mRNAs between the gout patients and normal people, and constructed co-regulated networks of circRNAs and lncRNAs according to the competitive endogenous RNA (ceRNA) theory for gouty arthritis onset to improve our understanding of the pathogenesis of this disease. The most significant finding of this study is the co-regulated ceRNA network of circRNAs and lncRNAs in gouty arthritis. The circRNA novel_circ_0030384 and the lncRNAs AAMP, TRIM16, PKN1, XLOC_184579 and XLOC_189826 were upstream genes in the co-regulated network. These upstream genes upregulated miR550a-5p and miR550a-3-5p, which downregulated PSME1 and FERMT3 expression. These mRNAs participated in proteasome dynamics, antigen processing and presentation, and platelet activation, which are associated with inflammation in gouty arthritis. In addition, the circRNA and lncRNAs upregulated miR550a-5p, which downregulated GRK2 and OS9 expression. Also, it proved that the down-regulated of PSME1, FERMT3, GRK2 and OS9 can aggravate gouty arthritis in vitro. In summary, these genes mediate inflammation in gouty arthritis through chemokine signaling to regulate neutrophil function.

Transcriptional Responses to Priority Effects in Nectar Yeast.

Priority effects, where the order and timing of species arrival influence the assembly of ecological communities, have been observed in a variety of taxa and habitats. However, the genetic and molecular basis of priority effects remains unclear, hindering a better understanding of when priority effects will be strong. We sought to gain such an understanding for the nectar yeast Metschnikowia reukaufii commonly found in the nectar of our study plant, the hummingbird-pollinated Diplacus (Mimulus) aurantiacus. In this plant, M. reukaufii can experience strong priority effects when it reaches flowers after other nectar yeasts, such as M. rancensis. After inoculation into two contrasting types of synthetic nectar simulating early arrival of M. rancensis, we conducted whole-transcriptome sequencing of 108 strains of M. reukaufii. We found that several genes were differentially expressed in M. reukaufii strains when the nectar had been conditioned by growth of M. rancensis. Many of these genes were associated with amino acid metabolism, suggesting that M. reukaufii strains responded molecularly to the reduction in amino acid availability caused by M. rancensis. Furthermore, investigation of expression quantitative trait loci (eQTLs) revealed that genes involved in amino acid transport and resistance to antifungal compounds were enriched in some genetic variants of M. reukaufii. We also found that gene expression was associated with population growth rate, particularly when amino acids were limited. These results suggest that intraspecific genetic variation in the ability of nectar yeasts to respond to nutrient limitation and direct fungal competition underpins priority effects in this microbial system.

Androgen production, uptake, and conversion (APUC) genes define prostate cancer patients with distinct clinical outcomes.

BACKGROUNDProstate cancer (PC) is driven by aberrant signaling of the androgen receptor (AR) or its ligands, and androgen deprivation therapies (ADTs) are a cornerstone of treatment. ADT responsiveness may be associated with germline changes in genes that regulate androgen production, uptake, and conversion (APUC).METHODSWe analyzed whole-exome sequencing (WES) and whole-transcriptome sequencing (WTS) data from prostate tissues (SU2C/PCF, TCGA, GETx). We also interrogated the Caris Precision Oncology Alliance (POA) DNA (592-gene/whole exome) and RNA (whole transcriptome) next-generation sequencing databases. Algorithm for Linking Activity Networks (ALAN) was used to quantify all pairwise gene-to-gene associations. Real-world overall survival was determined from insurance claims data using Kaplan-Meier estimates.RESULTSSix APUC genes (HSD3B1, HSD3B2, CYP3A43, CYP11A1, CYP11B1, CYP17A1) exhibited coalescent gene behavior in a cohort of metastatic tumors (n = 208). In the Caris POA dataset, the 6 APUC genes (APUC-6) exhibited robust clustering in primary prostate (n = 4,490) and metastatic (n = 2,593) biopsies. Surprisingly, tumors with elevated APUC-6 expression had statically lower expression of AR, AR-V7, and AR signaling scores, suggesting ligand-driven disease biology. APUC-6 genes instead associated with the expression of alternative steroid hormone receptors, ESR1/2 and PGR. We used RNA expression of AR or APUC-6 genes to define 2 subgroups of tumors with differential association with hallmark pathways and cell surface targets.CONCLUSIONSThe APUC-6-high/AR-low tumors represented a subgroup of patients with good clinical outcomes, in contrast with the AR-high or neuroendocrine PCs. Altogether, measuring the aggregate expression of APUC-6 genes in current genomic tests identifies PCs that are ligand (rather than AR) driven and require distinct therapeutic strategies.FUNDINGNCI/NIH 1R37CA288972-01, NCI Cancer Center Support P30 CA077598, DOD W81XWH-22-2-0025, R01 CA249279.

The clinical outcome, pathologic spectrum, and genomic landscape for 454 cases of salivary mucoepidermoid carcinoma

Mucoepidermoid carcinoma (MEC) is the most common malignant salivary tumor. A complete understanding of the high heterogeneity of MEC in histology and genetics would help in accurate diagnosis and treatment. Therefore, We evaluated the clinical features, treatment outcomes, and pathological parameters of 454 MECs and analyzed their genomic features using whole-exome sequencing and whole-transcriptome sequencing. We found that MECs predominantly occurred in females and those in their 4th–5th decades. The parotid gland was the most frequently affected site. All patients underwent complete mass resection with lobectomy; 414 patients were alive without relapse at follow-up, after an average period of 62 months (1–116 months). The disease progressed after initial treatment in 40 patients. The lungs were the most common site of distant metastasis. For classical MECs, histologic gradings of the AFIP, modified Healey, and MSK systems were significantly associated with recurrence and lymph nodal metastasis; these gradings were significantly related to lymph nodal metastasis for the subtypes. Older age, minor salivary gland involvement, clinical symptoms, high TNM stage, high-grade tumor, and improper surgical modality were the main prognostic factors. BAP1 was the most frequently mutated gene in MEC. Mutations in CDKN2A, MET, and TP53 were more frequently found in aggressive tumor phenotypes. MAML2 rearrangement was observed in 42% of patients, and EWSR1 rearrangement in 8%. Specific genetic events (in TP53 and FBXW7) with CRTC1::MAML2 fusion superimposed might be associated with unfavorable prognosis. This study provides new insights into precision therapeutic strategies for MEC.

Genome‑wide identification of circular RNAs and MAPKs reveals the regulatory networks in response to green peach aphid infestation in peach (Prunus persica)

The green peach aphid (GPA), Myzus persicae (Sulzer), is a serious agricultural pest with a worldwide distribution and a vector of over 100 plant viruses. Various pathways, such as the mitogen-activated protein kinase (MAPK) cascades, play pivotal roles in signaling plant defense against pest attack, and circular RNAs (circRNAs) regulate the expression of mRNAs in response to pest attack. However, the mechanism underlying peach (Prunus persica) response to GPA attack remains unclear. The present study initially identified and characterized 316 circRNAs and 18 PpMAPKs from healthy and GPA-infested peach leaves by whole-transcriptome sequencing and predicted the differentially expressed circRNAs (DECs) after GPA infestation. PCR and Sanger sequencing confirmed the presence of six DECs in peach samples. Besides, RNA sequencing analysis detected 13 DECs, including 5 upregulated and 8 downregulated ones, in peach in response to the GPA attack. Gene ontology (GO) enrichment analysis indicated that specific DECs play crucial roles in the MAPK signaling pathway, and qRT-PCR revealed that GPA infestation altered the expression patterns of PpMAPKs. Finally, five circRNAs, three microRNA (miRNAs), and two MAPK target genes were identified to interact as a network and perform critical roles in modulating the MAPK pathway in the peach during GPA infestation.

Phase II trial of imatinib mesylate in patients with PDGFRA/B-negative hypereosinophilic syndrome.

The role of imatinib in PDGFRA/B-negative hypereosinophilic syndromes (HES) is controversial because of the heterogeneity of HES and the scarcity of prospective studies. We conducted a phase II clinical trial to evaluate the efficacy of imatinib in PDGFRA/B-negative HES. Thirty-two patients were treated with imatinib (100-400 mg daily), and the molecular basis of their response was identified using whole-exome sequencing (WES) and whole-transcriptome sequencing (WTS). The haematological response rate was 46.9%, with a complete haematological response (CHR) rate of 18.8%. The median time to response was 1.5 months. Among the six patients who achieved CHR, five maintained it until the 24th cycle of imatinib and one lost response after 20 months. The median progression-free survival was 4.3 months. WES and WTS were conducted for 11 patients. The number of non-silent mutations did not differ between responders and non-responders. Nine differentially expressed genes, including SNORD15A, were downregulated in responders. STAT5B::RARA, PAK2::PIGX, and FIP1L1::CHIC2 fusions were identified in patients with sustained responses, and RNF130::BRAF and WNK1::KDM5A fusions were identified in non-responders. Imatinib, along with an appropriate biomarker, could be a promising option for PDGFRA/B-negative HES.

Inhibition of the FOXO1–ROCK1 axis mitigates cardiomyocyte injury under chronic hypoxia in Tetralogy of Fallot by maintaining mitochondrial quality control

IntroductionPersistent chronic myocardial hypoxia causes disturbances in mitochondrial quality control (MQC), ultimately leading to increased cardiomyocyte injury in patients with Tetralogy of Fallot (TOF). The present study aimed to identify the key effector molecules of cardiomyocyte injury under chronic hypoxia in TOF. MethodsClinical data from TOF patients were collected and whole transcriptome sequencing was performed on myocardial samples. Chronic hypoxia models were established in cardiac-specific knockout mice and cardiomyocytes, and a series of molecular experiments were used to determine the specific mechanisms involved. ResultsClinical cohort data and whole-transcriptome sequencing analysis of myocardial samples from TOF patients revealed that forkhead box O1 (FOXO1) plays an important role in chronic hypoxic cardiomyocyte injury. In a model of chronic hypoxia established in FOXO1 cardiac-specific knockout mice and FOXO1 gene-deficient cardiomyocytes, the AMPK signaling pathway regulates the expression of FOXO1, which in turn disrupts MQC by regulating the transcriptional activation of Rho-associated protein kinase 1 (ROCK1), and increasing the production of mitochondrial ROS, thereby exacerbating damage to cardiomyocytes. Excessive reactive oxygen species (ROS) production during MQC dysfunction further activates Cox7a2L to increase the assembly of the respiratory chain supercomplex. In addition, we found that miR-27b-3p partially binds to the 3′ untranslated region of FOXO1 to exert a protective effect. ConclusionsMaintenance of MQC under chronic hypoxia is achieved through a series of injury-protection mechanisms, suggesting that FOXO1 inhibition may be crucial for future mitigation of chronic hypoxic cardiomyocyte injury in TOF.

The N-terminal fragment of histone deacetylase 4 (1-669aa) promotes chondrocyte apoptosis via the p53-dependent endoplasmic reticulum stress pathway.

Exogenous administration of the histone deacetylation 4 (HDAC4) protein can effectively delay osteoarthritis (OA) progression. However, HDAC4 is unstable and easily degrades into N-terminal (HDAC4-NT) and C-terminal fragments, and the HDAC4-NT can exert biological effects, but little is known about its role in chondrocytes and cartilage. Thus, the roles of HDAC4-NT fragments (1-289aa, 1-326aa and 1-669aa) in chondrocytes and cartilage were evaluated via real-time cell analysis (RTCA), safranin O staining, Sirius Red staining and nanoindentation. Molecular mechanisms were profiled via whole-transcriptome sequencing (RNA-seq) and verified invitro and invivo by a live cell real-time monitoring system, flow cytometry, western blotting and immunohistochemistry. The results showed that 1-669aa induced chondrocyte death and cartilage injury significantly, and the differentially expressed genes (DEGs) were enriched mainly in the apoptotic term and p53 signalling pathway. The validation experiments showed that 1-669aa induced chondrocyte apoptosis via the endoplasmic reticulum stress (ERS) pathway, and up-regulated p53 expression was essential for this process. Thus, we concluded that the HDAC4-NT fragment 1-669aa induces chondrocyte apoptosis via the p53-dependent ERS pathway, suggesting that in addition to overexpressing HDAC4, preventing it from degradation may be a new strategy for the treatment of OA.

Dissecting the Significance of Acid Phosphatase 1 Gene Alterations in Prostate Cancer.

The acid phosphatase 1 (ACP1) gene encodes low-molecular-weight protein tyrosine phosphatase, which is overexpressed in prostate cancer (PC) and a potential therapeutic target. We analyzed ACP1 expression in primary/metastatic PC and its association with molecular profiles and clinical outcomes. NextGen sequencing of DNA (592-gene/whole-exome sequencing)/RNA(whole-transcriptome sequencing) was performed for 5,028 specimens. ACP1-High/ACP1-Low expression was defined as quartile (Q4/1) of RNA transcripts per million (TPM). DNA mutational profiles were analyzed for ACP1-quartile-stratified samples. Gene set enrichment analysis was used for Hallmark collection of pathways. PD-L1+(≥2+, ≥5%; SP142) was tested by immunohistochemistry. Tumor microenvironment's (TME) immune cell fractions were estimated by RNA deconvolution/quanTIseq. Overall survival (OS) was assessed from initial diagnosis/treatment initiation to death/last follow-up. We included 3,058 (60.8%) samples from the prostate, 634 (12.6%) from lymph node metastases (LNMs), and 1,307 (26.0%) from distant metastases (DMs). ACP1 expression was higher in LNM/DM than prostate (49.8/47.9 v 44.1 TPM; P < .0001). TP53 mutations were enriched in ACP1-Q4 (37.9%[Q4] v 27.0%[Q1]; P < .001) among prostate samples. Pathways associated with cell cycle regulation and oxidative phosphorylation were enriched in ACP1-Q4, whereas epithelial-mesenchymal transition and tumor necrosis factor-alpha signaling via nuclear factor kappa-light-chain-enhancer of activated B-cell pathways were enriched in ACP1-Q1. Neuroendocrine and androgen receptor signaling was increased in ACP1-Q4. M2 macrophages and natural killer cell fractions were increased, whereas T cells and M1 macrophages were decreased in ACP1-Q4. While OS differences between ACP1-Q1/Q4 were not statistically significant, there was a trend for worse OS among ACP1-Q4 prostate samples (Q4 v Q1: hazard ratio [HR], 1.19 [95% CI, 0.99 to 1.42]; P = .06) and DM (HR, 1.12 [95% CI, 0.93 to 1.36]; P = .22) but not LNM (HR, 0.98 [95% CI, 0.74 to 1.29]; P = .87). ACP1-High tumors exhibit a distinct molecular profile and cold TME, highlighting ACP1's potential role in PC pathogenesis and novel therapeutic targeting.

Whole-transcriptome analyses of ovine lung microvascular endothelial cells infected with bluetongue virus

Bluetongue virus (BTV) infection induces profound and intricate changes in the transcriptional profile of the host to facilitate its survival and replication. However, there have been no whole-transcriptome studies on ovine lung microvascular endothelial cells (OLMECs) infected with BTV. In this study, we comprehensively analysed the whole-transcriptome sequences of BTV-1 serotype-infected and mock-infected OLMECs and subsequently performed bioinformatics differential analysis. Our analysis revealed 1215 differentially expressed mRNA transcripts, 82 differentially expressed long noncoding RNAs (lncRNAs) transcripts, 63 differentially expressed microRNAs (miRNAs) transcripts, and 42 differentially expressed circular RNAs (circRNAs) transcripts. Annotation from Gene Ontology, enrichment from the Kyoto Encyclopedia of Genes and Genomes, and construction of endogenous competing RNA network analysis revealed that the differentially expressed RNAs primarily participated in viral sensing and signal transduction pathways, antiviral and immune responses, inflammation, and extracellular matrix (ECM)-related pathways. Furthermore, protein‒protein interaction network analysis revealed that BTV may regulate the conformation of ECM receptor proteins and change their biological activity through a series of complex mechanisms. Finally, on the basis of real-time fluorescence quantitative polymerase chain reaction results, the expression trends of the differentially expressed RNA were consistent with the whole-transcriptome sequencing data, such as downregulation of the expression of COL4A1, ITGA8, ITGB5, and TNC and upregulation of the expression of CXCL10, RNASEL, IRF3, IRF7, and IFIHI. This study provides a novel perspective for further investigations of the mechanism of the ECM in the BTV-host interactome and the pathogenesis of lung microvascular endothelial cells.

Lif-deficiency promote systemic Iron metabolism disorders and increases the susceptibility of osteoblasts to ferroptosis

Leukemia inhibitory factor (LIF) is a multifunctional cytokine that plays a crucial role in various biological processes. However, LIF involvement in iron metabolism remains almost unexplored. This study aimed to explore the impact of LIF on systemic iron transportation and its potential role in ferroptosis in osteoblasts. We observed that the Lif-deficient (Lif−/−) mice is characterized by a reduction in bone mass and a decrease in bone mineral density compared with wild-type (WT) mice. Energy-dispersive X-ray spectroscopy revealed a marked increase in iron content on the surface of femurs from Lif−/− mice. Meanwhile, iron stores test lower iron levels in the spleens and higher levels in the femurs of Lif−/− mice. Besides, Lif−/− mice display increased levels of serum iron, total iron-binding capacity, unsaturated iron-binding capacity, and transferrin saturation and serum ferritin relative to WT mice. Hepcidin mRNA expression reduction in the liver of Lif−/− mice. It also holds true in the AML-12 hepatocyte cell line after Lif-knockdown. Immunohistochemistry and RT-PCR revealed elevated ferroportin (FPN) in duodenal cells of Lif−/− mice. Lif-deficiency decreases SLC7A11 levels in osteoblasts. In addition, overexpression of LIF downregulates CD71, DCYTB, and DMT1, thereby reducing iron uptake in iron-overloaded cells. Femur immunohistochemistry (IHC) revealed increased ACSL4 and decreased GPX4 and SLC7A11, indicating an increase in ferroptosis of osteoblasts in Lif−/− mice. Whole-transcriptome sequencing showed gene expression changes after Lif-knockdown, exhibiting a negative correlation with genes involved in long-chain fatty acid transport, mitochondrial organization, and the p38 MAPK signaling pathway. These results demonstrate that Lif-deficiency alter systemic iron metabolism and increases the susceptibility of osteoblasts to ferroptosis.

RNA Sequencing of Sperm from Healthy Cattle and Horses Reveals the Presence of a Large Bacterial Population.

RNA molecules within ejaculated sperm can be characterized through whole-transcriptome sequencing, enabling the identification of pivotal transcripts that may influence reproductive success. However, the profiling of sperm transcriptomes through next-generation sequencing has several limitations impairing the identification of functional transcripts. In this study, we explored the nature of the RNA sequences present in the sperm transcriptome of two livestock species, cattle and horses, using RNA sequencing (RNA-seq) technology. Through processing of transcriptomic data derived from bovine and equine sperm cell preparations, low mapping rates to the reference genomes were observed, mainly attributed to the presence of ribosomal RNA and bacteria in sperm samples, which led to a reduced sequencing depth of RNAs of interest. To explore the presence of bacteria, we aligned the unmapped reads to a complete database of bacterial genomes and identified bacteria-associated transcripts which were characterized. This analysis examines the limitations associated with sperm transcriptome profiling by reporting the nature of the RNA sequences among which bacterial RNA was found. These findings can aid researchers in understanding spermatozoal RNA-seq data and pave the way for the identification of molecular markers of sperm performance.

Multiomics Profiling Distinguishes Sebaceous Carcinoma from Benign Sebaceous Neoplasms and Provides Insight into the Genetic Evolution of Sebaceous Carcinogenesis.

Sebaceous carcinoma is the third most common nonkeratinocyte skin cancer in the United States with 1,000 cases per year. The clinicopathologic features of sebaceous carcinoma and benign sebaceous neoplasms (adenomas, sebaceomas) can overlap, highlighting the need for molecular biomarkers to improve classification. This study describes the genomic and transcriptomic landscape of sebaceous neoplasms in order to understand tumor etiology and biomarkers relevant for diagnosis and treatment. We performed whole-genome sequencing (WGS) and whole-transcriptome sequencing (WTS) of sebaceous neoplasms from six academic and two federal healthcare facilities in the United States diagnosed between January 1, 1999, and December 31, 2021. We evaluated 98 sebaceous neoplasms: 64 tumors (32 adenomas, 2 sebaceomas, 5 atypical sebaceous neoplasms, 25 carcinomas) had sufficient material for WGS, 96 tumors (42 adenomas, 11 sebaceomas, 8 atypical sebaceous neoplasms, 35 carcinomas) had sufficient material for WTS, and 62 tumors (31 adenomas, 2 sebaceomas, 5 atypical sebaceous neoplasms, 24 carcinomas) had sufficient material for combined WGS and WTS. Overall, we found decreased cholesterol biosynthesis and increased TP53 mutations, copy number gains (chromosome 6, 8q, and/or 18), and tumor mutation burden-high (>10 mutations/MB) in carcinomas compared to adenomas. Although diminished compared to adenomas, most carcinomas still had higher cholesterol biosynthesis than nonmalignant skin. Multiomics profiling also supported a precancerous model of tumor evolution with sebaceomas and atypical sebaceous neoplasms being likely intermediate lesions. The study findings highlight key diagnostic biomarkers for sebaceous carcinoma and suggest that immunotherapy and modulation of cholesterol biosynthesis could be effective treatment strategies.

TRANSCRIPTOME PROFILING AND EXAMINATION OF INDIVIDUALS DIAGNOSED WITH ESOPHAGEAL SQUAMOUS CELL CARCINOMA ORIGINATING FROM KAZAKHSTAN

Esophageal cancer stands as the eighth most prevalent cancer globally and the sixth most prevalent in Kazakhstan, with esophageal squamous cell carcinoma (ESCC) as its primary histological subtype, often diagnosed in advanced stages. This project aimed to unveil the genetic underpinnings of ESCC by examining differentially expressed genes (DEGs) via whole-transcriptome sequencing of Kazakhstani patients. Tissue samples from 22 individuals undergoing Ivor-Lewis esophagectomy at the Oncology Center in Astana were collected. Utilizing STAR software and DESeq2 package, DEGs were mapped and defined, followed by functional analysis using various R packages. The study encompassed 13 males and 9 females, averaging 65.72 ± 8.26 years old, with 86% diagnosed at advanced stages T3-T4. Analysis of tumor and normal esophageal tissues unveiled 6689 DEGs, including 2056 upregulated and 4633 downregulated genes (adjusted p-value &lt;0.05). Additionally, 42 up-regulated and 2 down-regulated KEGG pathways were identified (p-value &lt;0.05). The top 300 DEGs were mapped to a protein-protein interaction (PPI) network, and functional enrichment analysis revealed three modules comprising closely connected genes. Within our patient cohort, ESCC with moderate dysplasia emerged as the most prevalent histologic subtype, associated with poor prognosis. High-throughput sequencing techniques facilitate the identification of molecular pathways implicated in esophageal carcinogenesis, potentially revolutionizing diagnosis and treatment strategies. This research has been funded by the Science Committee of the Ministry of Science and Higher Education of the Republic of Kazakhstan (Program targeted funding BR18574184 and grant AP09058660), and Nazarbayev University funding CRP grant 021220CRP2222.

Identification of competing endogenous RNA networks associated with circRNA and lncRNA in TCDD-induced cleft palate development

2,3,7,8 -tetrachlorodibenzo-p-dioxin (TCDD) is a teratogen that can induce cleft palate formation, a common birth defect. Competing endogenous RNAs (ceRNAs), including circular RNAs (circRNAs) and long non-coding RNAs (lncRNAs), indirectly regulate gene expression via sharing microRNAs (miRNAs). Nevertheless, the mechanism by which they act as ceRNAs to regulate palatal development remains to be explored in greater detail. Here, the cleft palate model of C57BL/6 N pregnant mice was constructed by gavage of TCDD (64 ug/kg) on gestation day (GD) 10.5, and the palatal shelves were taken on gestation day (GD) 14.5 for whole-transcriptome sequencing to investigate the underlying mechanisms of the roles of circRNAs and lncRNAs as ceRNAs in cleft palate. Sequencing results revealed that 293 lncRNA, 589 circRNA, 47 miRNA, and 138 messenger RNA (mRNA) were significantly dysregulated, and the cytochrome P450 (CYP) enzymes and the aryl hydrocarbon receptor (AhR) pathway play key roles in the induction of cleft palate upon exposure to TCDD. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis revealed the function of TCDD function was mainly related to the metabolic processes of intracellular compounds, including the metabolic processes of cellular aromatic compounds and the metabolism of exogenous drugs by cytochrome P450, etc. Furthermore, quantitative reverse transcription polymerase chain reaction (qRT-PCR) indicated that the circRNA_1781/miR-30c-1–3p/PKIB and XR_380026.2/miR-1249–3p/DNAH10 ceRNA networks were hypothesized to be a hub involved in palatal development suggesting that the circRNA_1781/miR-30c-1–3p/PKIB and XR_380026.2/miR-1249–3p/DNAH10 ceRNA networks may be critical for palatogenesis, setting the foundation for the investigation of cleft palate.

Long non-coding RNA MAGEA4-AS1 binding to p53 enhances MK2 signaling pathway and promotes the proliferation and metastasis of oral squamous cell carcinoma

Long non-coding RNAs (lncRNAs) regulate the occurrence, development and progression of oral squamous cell carcinoma (OSCC). We elucidated the expression features of MAGEA4-AS1 in patients with OSCC and its activity as an OSCC biomarker. Furthermore, the impact of up-regulation of MAGEA4-AS1 on the cellular behaviors (proliferation, migration and invasion) of OSCC cells and intrinsic signal mechanisms were evaluated. Firstly, we analyzed MAGEA4-AS1 expression data in The Cancer Genome Atlas (TCGA) OSCC using a bioinformatics approach and in 45 pairs of OSCC tissues using qPCR. Then CCK-8, ethynyl deoxyuridine, colony formation, transwell and wound healing assays were conducted to assess changes in the cell proliferation, migration and invasion protential of shMAGEA4-AS1 HSC3 and CAL27 cells. The RNA sequence of MAGEA4-AS1 was identified using the rapid amplification of cDNA ends (RACE) assay. And whole-transcriptome sequencing was used to identify MAGEA4-AS1 affected genes. Additionally, dual-luciferase reporter system, RNA-binding protein immunoprecipitation (RIP), and rescue experiments were performed to clarify the role of the MAGEA4-AS1-p53-MK2 signaling pathway. As results, we found MAGEA4-AS1 was up-regulated in OSCC tissues. We identified a 418 nucleotides length of the MAGEA4-AS1 transcript and it primarily located in the cell nucleus. MAGEA4-AS1 stable knockdown weakened the proliferation, migration and invasion abilities of OSCC cells. Mechanistically, p53 protein was capable to activate MK2 gene transcription. RIP assay revealed an interaction between p53 and MAGEA4-AS1. MK2 up-regulation in MAGEA4-AS1 down-regulated OSCC cells restored MK2 and epithelial-to-mesenchymal transition related proteins’ expression levels. In conclusion, MAGEA4-AS1-p53 complexes bind to MK2 promoter, enhancing the transcription of MK2 and activating the downstream signaling pathways, consequently promoting the proliferation and metastasis of OSCC cells. MAGEA4-AS1 may serve as a diagnostic marker and therapeutic target for OSCC patients.

Consensus molecular subtyping of metastatic colorectal cancer expands biomarker-directed therapeutic benefit for patients with CMS1 and CMS2 tumors.

We developed a whole transcriptome sequencing (WTS)-based Consensus Molecular Subtypes (CMS) classifier using FFPE tissue and investigated its prognostic and predictive utility in a large clinico-genomic database of CRC patients (n = 24,939). The classifier was trained against the original CMS datasets using an SVM model and validated in an independent blinded TCGA dataset (88.0% accuracy). Kaplan-Meier estimates of overall survival (OS) and time-on-treatment (TOT) were calculated for each CMS (p < 0.05 considered significant). CMS2 tumors were enriched on left-side of colon and conferred the longest median OS. In RAS-wildtype mCRC, left-sided tumors and CMS2 classification were associated with longer TOT with anti-EGFR antibodies (cetuximab and panitumumab). When restricting to only CMS2, there was no significant difference in TOT between right- versus left-sided tumors. CMS1 tumors were associated with a longer median TOT with pembrolizumab relative to other CMS groups, even when analyzing only microsatellite stable (MSS) tumors. A WTS-based CMS classifier allowed investigation of a large multi-institutional clinico-genomic mCRC cohort, suggesting anti-EGFR therapy benefit for right-sided RAS-WT CMS2 tumors and immune checkpoint inhibitor benefit for MSS CMS1. Routine CMS classification of CRC provides important treatment associations that should be further investigated.

Exploration and validation of ceRNA regulatory networks in colorectal cancer based on associations whole transcriptome sequencing

Colorectal cancer (CRC) is a wide-spread gastrointestinal cancer that is associated with augmented morbidity and mortality, and we do not yet have a deep understanding of its epidemiology and carcinogenicity. The transcriptome can reveal the complexity and heterogeneity of tumors and uncover new biomarkers or treatment options. In this study, we identified messenger RNAs (mRNAs), long non-coding RNAs (lncRNAs), round RNAs (circRNAs), and microRNAs (miRNAs) using whole-transcriptome sequencing and generated competing endogenous RNA (ceRNA) modulatory axes. We conducted whole transcriptome sequencing on 10 CRC and para-cancer (CRCP) samples and discovered 2465 differentially expressed (DE) mRNAs (DEmRNAs), 77 DE miRNAs (DEmiRNAs). 2852 DE lncRNAs (DElncRNAs) and 1477 DE circRNAs (DEcircRNAs). In addition, utilizing co-DE analysis, we generated the ceRNA axis. Subsequently, we employed the ceRNA axis to identify essential genes and corresponding associations with lncRNAs, circRNAs, and miRNAs in CRC. ceRNA regulatory network including mRNA-miRNA-lncRNA and mRNA-miRNA-circRNA. These modulatory axes potentially modulate the positive regulation of smooth muscle contraction, melanosome, plasma membrane, integral plasma membrane component and so on. Finally, the results of RNA sequencing (RNA-SEQ) were combined with the TCGA and GEO databases, and the DEGs strongly correlated with the TCGA-COAD overall survival (OS) as estimated by univariate cox and logarithmic rank analyses were cross-analyzed, and the co-upregulated DEGs were screened. Among the many DEs, KPNA2 was chosen for additional analysis. Using invitro experimentations, western blot, CCK8, EdU and other experiments were performed to verify the results. We found siRNA-based KPNA2 depletion reduces bladder cancer cells’ viability, migratory, and proliferative activities, which showed that the DEmRNA profiles were comparable to the sequencing information, confirming that the sequencing data were very reliable. These evidences highlight the ceRNA regulatory mechanisms in CRC and will aid future research into the molecular mechanisms behind colorectal cancer prevention and treatment.

Whole-transcriptome sequencing analysis to identify key circRNAs, miRNAs, and mRNAs in the development of yak testes

BackgroundThe Testis is an important reproductive organ in male mammals and the site for spermatogenesis, androgen synthesis, and secretion. Non-coding RNAs (ncRNAs) play an important regulatory role in various biological processes. However, the regulatory role of ncRNAs in the development of yak testes and spermatogenesis remains largely unclear.ResultIn this study, we compared the expression profiles of circular RNAs (circRNAs), microRNAs (miRNAs), and messenger RNAs (mRNAs) in yak testicular tissue samples collected at 6 months (Y6M), 18 months (Y18M), and 4 years (Y4Y). Using RNA sequencing (RNA-Seq), we observed a significant difference in the expression patterns of ncRNAs in the samples collected at different testicular development stages. Twenty-two differentially expressed (DE) circRNAs, 69 DE miRNAs, and 64 DE mRNAs were detected in Y6M, Y18M, and Y4Y testicular samples, respectively. The results of gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses showed that the source genes of DE circRNAs, predicted target genes of DE miRNAs, and DE mRNAs were specifically associated with signaling pathways and GO terms that were related to sperm synthesis, sperm vitality, and testicular development, such as cell cycle, Wnt signaling pathway, MAPK signaling pathway, GnRH signaling pathway, and spermatogenesis. The analysis of the circRNA-miRNA-mRNA network revealed that some DE ncRNAs, including miR-574, miR-449a, CDC42, and CYP11A1, among others, may be involved in testicular spermatogenesis. Concurrently, various circRNA-miRNA interaction pairs were observed.ConclusionOur findings provide a database of circRNAs, miRNAs, and mRNAs expression profiles in testicular tissue of yaks at different developmental stages and a detailed understanding of the regulatory network of ncRNAs in yak testicular development and provide data that can help elucidate the molecular mechanisms underlying yak testicular development.

Unraveling the impact of lysosomal dysfunction on myeloproliferative neoplasm.

Lysosomal dysfunction (LD) impacts cytokine regulation, inflammation, and immune responses, influencing the development and progression of cancer. Inflammation is implicated in the pathogenesis of myeloproliferative neoplasm (MPN). With a hypothesis that LD significantly contributes to MPN carcinogenesis by inducing abnormal inflammation, our objective was to elucidate the pathophysiological mechanisms of MPN arising from an LD background. Genotyping of the LD background was performed in a cohort of MPN patients (n = 190) and healthy controls (n = 461). Logistic regression modeling, utilizing genotype data, was employed to estimate the correlation between LD and MPN. Whole transcriptome sequencing (WTS) (LD carriers = 8, non-carriers = 6) and single-cell RNA sequencing data (LD carriers = 2, non-carriers = 2, healthy controls = 2) were generated and analyzed. A higher variant frequency of LD was observed in MPN compared to healthy controls (healthy, 4.9%; MPN, 7.8%), with the highest frequency seen in polycythemia vera (PV) (odds ratio = 2.33, p = 0.03). WTS revealed that LD carriers exhibited upregulated inflammatory cytokine ligand-receptor genes, pathways, and network modules in MPNs compared to non-carriers. At the single-cell level, there was monocyte expansion and elevation of cytokine ligand-receptor interactions, inflammatory transcription factors, and network modules centered on monocytes. Notably, Oncostatin-M (OSM) consistently emerged as a candidate molecule involved in the pathogenesis of LD-related PV. In summary, an LD background is prevalent in MPN patients and leads to increased cytokine dysregulation and inflammation. OSM, as one of the potential molecules, plays a crucial role in PV pathogenesis by impairing lysosomal function.