Whole-transcriptome Sequencing Research Articles

Glycolysis and de novo fatty acid synthesis cooperatively regulate pathological vascular smooth muscle cell phenotypic switching and neointimal hyperplasia.

Switching of vascular smooth muscle cells (VSMCs) from a contractile phenotype to a dedifferentiated (proliferative) phenotype contributes to neointima formation, which has been demonstrated to possess a tumor-like nature. Dysregulated glucose and lipid metabolism is recognized as a hallmark of tumors but has not thoroughly been elucidated in neointima formation. Here, we investigated the cooperative role of glycolysis and fatty acid synthesis in vascular injury-induced VSMC dedifferentiation and neointima formation. We found that the expression of hypoxia-inducible factor-1α (HIF-1α) and its target 6-phosphofructo-2-kinase/fructose-2,6-bisphosphatase (PFKFB3), a critical glycolytic enzyme, were induced in the neointimal VSMCs of human stenotic carotid arteries and wire-injured mouse carotid arteries. HIF-1α overexpression led to elevated glycolysis and resulted in a decreased contractile phenotype while promoting VSMC proliferation and activation of the mechanistic target of rapamycin complex 1 (mTORC1). Conversely, silencing Pfkfb3 had the opposite effects. Mechanistic studies demonstrated that glycolysis generates acetyl coenzyme A to fuel de novo fatty acid synthesis and mTORC1 activation. Whole-transcriptome sequencing analysis confirmed the increased expression of PFKFB3 and fatty acid synthetase (FASN) in dedifferentiated VSMCs. More importantly, FASN upregulation was observed in neointimal VSMCs of human stenotic carotid arteries. Finally, interfering with PFKFB3 or FASN suppressed vascular injury-induced mTORC1 activation, VSMC dedifferentiation, and neointima formation. Together, this study demonstrated that PFKFB3-mediated glycolytic reprogramming and FASN-mediated lipid metabolic reprogramming are distinctive features of VSMC phenotypic switching and could be potential therapeutic targets for treating vascular diseases with neointima formation. © 2023 The Pathological Society of Great Britain and Ireland.

Identification of non-coding RNAs and their functional network associated with optic nerve invasion in retinoblastoma

Optic nerve invasion (ONI) is an important high-risk feature and prognostic indicator of retinoblastoma (RB). Emerging evidence has revealed that non-coding RNAs (ncRNAs) play important roles in tumor perineural invasion (PNI). Nevertheless, the regulatory role of ncRNAs in the ONI of RB is poorly understood. In the current study, whole-transcriptome sequencing was performed to assess the expression profiles of ncRNAs and mRNAs in RB tissues, with or without ONI. Based on Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses, we predicted the biological functions of differentially expressed (DE) mRNAs. We then constructed competing endogenous RNA (ceRNA) regulatory networks based on bioinformatics analysis. The hsa_circ_0015965/lncRNA MEG3-hsa-miR-378a-5p-NOTCH1 pathway was selected and validated by real-time qPCR, western blotting, and dual luciferase reporter assays. Moreover, we demonstrated that NOTCH1 promotes the malignant progression of RB. Taken together, our results provide novel insights into the mechanism underlying optic nerve invasion in RB.

The role of gene expression of CDC37 in colorectal cancer (CRC).

237 Background: CDC37-HSP90 axis is an essential chaperoning system for stabilization of kinases. CDC37 determines selectivity of client kinases recognized by HSP90. We previously showed patients with CDC37-dependent ( CDC37 high expression) colorectal cancer (CRC) derived more benefit from regorafenib and bevacizumab both of which target HSP90 client kinases or signaling pathway, but not from cetuximab which targets HSP90 non-client kinase. However, molecular characteristics and interaction with relevant signaling pathways in CDC37-dependent CRC are largely unknown. Methods: We retrospectively reviewed CRC patient samples submitted to a commercial CLIA-certified laboratory (Caris Life Sciences, Phoenix AZ). Next-generation sequencing of DNA and RNA (whole-transcriptome sequencing) and immunohistochemistry were performed. Molecular profiles between top quartile transcripts per million for CDC37 expression (Q4, CDC37-high: n = 5056) and bottom quartile (Q1, CDC37-low: n = 5056) were compared for gene mutations, microsatellite instability (MSI) status, PD-L1 expression, co-expression of HSP90 genes ( HSP90AA1 and HSP90AB1), MAPK pathway activity score (MPAS), T-cell inflamed signature, and cell infiltration in the tumor microenvironment (TME) assessed by QuantiSEQ. Gene set enrichment analysis (GSEA) was performed between CDC37-high (Q4) and CDC37-low (Q1) tumors. Results: KRAS mutations trended higher in CDC37-high tumors compared to CDC37-low tumors (49.5% vs 47.1%) while APC (76.9% vs 72.3%), TP53 (75.5% vs 71.8%) and TCF7L2 (5.0% vs 3.5%) mutations were significantly higher in CDC37-high tumors compared to CDC37-low tumors (q < 0.05). No significance in expression quartiles were observed for dMMR/MSI-H or PD-L1 expression. HSP90AA1 and HSP90AB1 expression levels were significantly higher in CDC37-high tumors (2.9-fold and 2.6-fold in median, respectively; both q< 0.05). Both MPAS and T-cell inflamed signatures were significantly higher in CDC37-high tumors ( q< 0.05). Infiltration of NK cells, Tregs, neutrophils, M1 macrophages, and B-cells were greater in CDC37-high tumors (q < 0.05). GSEA showed several signaling pathways, such as KRAS, TGF-beta, hypoxia, WNT-beta catenin, and PI3K-AKT-MTOR, were enriched in CDC37-high tumors. Conclusions: Highly CDC37-expressed (CDC37-dependent) CRC was associated with aberrant TME where abundance of immune cell infiltrate and broad activation of kinase-related signaling pathways were evident. Our results support the kinase-stabilizing function of CDC37 in CRC. Further investigations combined with survival data are warranted to address the prognostic and predictive values of CDC37 expression in targeted therapies of CRC.

The role of gene expression of CD47 in colorectal cancer (CRC).

240 Background: CD47 belongs to the immunoglobulin superfamily and is overexpressed in many tumor types. CD47 plays an important role in suppressing phagocytosis through binding to transmembrane protein SIRP-alpha on macrophages. Targeting CD47 is a novel strategy for cancer immunotherapy and is being evaluated in ongoing clinical trials. However, molecular characteristics of CD47-overexpressed colorectal cancer (CRC) are largely unknown. Methods: We retrospectively reviewed CRC patient samples (n = 14786) submitted to a commercial CLIA-certified laboratory (Caris Life Sciences, Phoenix AZ). Next-generation sequencing of DNA and RNA (whole-transcriptome sequencing) and immunohistochemistry were performed. Correlation of CD47 expression with danger-associated molecular pattern (DAMP)-related genes ( HMGB1, CALR, ANXA1, HSP90AA1, HSPA1A, and CXCL10) expressions was tested. DAMP signature calculated as composite z-score of the DAMP related genes was compared between CD47-high and -low patients classified according to the median level of CD47 expression. Distributions of KRAS and BRAF mutations, consensus molecular subtype (CMS), and signatures of oncogenic signaling pathways were compared between CD47-high and -low patients. In addition, overall survival (OS) was compared between CD47-high and -low patients available for survival data. Results: CD47 expression level was significantly higher in metastatic compared to primary lesions (1.07-fold, q < 0.05) and microsatellite instability high tumors compared to microsatellite stable tumors (1.15-fold, q < 0.05). CD47 expression was positively correlated with DAMP genes expression except for HSPA1A, and the DAMP signature (median score [MS]: 2.66 vs -2.29, q < 0.05) was significantly increased in CD47-high patients. KRAS mutations were less prevalent (45.87% vs 50.05%) and CMS1(17.72% vs 14.42%) and CMS4 (40.33% vs 27.28%) were more prevalent in CD47-high patients, while no difference was observed in the prevalence of BRAF mutations between CD47-high and -low patients. Signatures of EMT (MS:2.91 vs -3.19), TGF-beta (MS:3.72 vs -3.52), angiogenesis (MS:2.95 vs -2.78), MAPK (MS:4.99 vs -4.12), PI3K (MS:2.41 vs -2.03), and immune-related signaling pathways (MS:1.20 vs -2.51) were significantly enriched in CD47-high patients (all q < 0.05). CD47-high patients (n = 4873) showed significantly worse OS than CD47-low patients (n = 4898) (median OS, 32.4 vs 37.6 months; hazard ratio = 1.158, p< 0.01). Conclusions: Highly CD47-expressed CRC harbored activation of DAMPs and oncogenic signaling pathways that linked to aberrant tumor microenvironment and worse prognosis. Our results support intensive treatment strategies using CD47 inhibitors combined with cytotoxic agents and molecular targeted agents (such as anti-VEGF agents) in CD47-overexpressed CRC.

An Integrated Analysis Reveals Ciliary Abnormalities in Antrochoanal Polyps.

The mechanisms underlying the antrochoanal polyps (ACPs) remained unclear. We aimed to identify the differentially expressed genes (DEGs) profile, the cilia-related genes expression levels and the morphological characteristics of ciliated cells in patients with ACPs. We obtained ACPs biopsy samples from 28 patients and uncinate process from 27 healthy controls. Whole-transcriptome RNA sequencing, immunofluorescence staining, quantitative polymerase chain reaction, and scanning electron microscopy were performed. 3739 DEGs were detected between ACPs and controls, and Gene Ontology analysis on these DEGs implicated cilium assembly, cilium motility, cilia component, cilia function, inflammatory response and immune system process were included in ACPs pathogenesis. Gene set enrichment analysis implicated sets of genes regulated in processes associated with cilium organization, cilium morphogenesis, cilium movement, axoneme assembly, axonemal dynein complex assembly and cell projection assembly. The expression levels of cilia-related genes (FOXJ1, DNAI1, DNAH9, RSPH1, RSPH9 and RSPH4A) were validated by quantitative polymerase chain reaction (Fold change >2, P<0.05) and FOXJ1 was positive correlated with DNAI1, DNAH9, RSPH4, RSPH1, RSPH9, DNAH5, DNALI1 in ACPs (all P < 0.05). Based on our semi-quantitative scoring system, median scores of α-Tubulin, DNAI1 and RSPH4A were significantly higher in ACPs than in controls. In addition, loss of ciliated cells and a shorter cilia pattern were further confirmed by immunofluorescence staining and scanning electron microscopy in ACPs. The aberrant expression of cilia-related genes and ciliary structural impairment are an important pathological phenomenon in ACPs, and our findings may provide novel insights into understanding the mysterious mechanisms underlying ACPs.

Identification and validation of coding and non-coding RNAs involved in high-temperature-mediated seed dormancy in common wheat.

Seed dormancy (SD) significantly decreases under high temperature (HT) environment during seed maturation, resulting in pre-harvest sprouting (PHS) damage under prolonged rainfall and wet weather during wheat harvest. However, the molecular mechanism underlying HT-mediated SD remains elusiveSeed dormancy (SD) significantly decreases under high temperature (HT) environment during seed maturation, resulting in pre-harvest sprouting (PHS) damage under prolonged rainfall and wet weather during wheat harvest. However, the molecular mechanism underlying HT-mediated SD remains elusive. Here, the wheat landrace 'Waitoubai' with strong SD and PHS resistance was treated with HT from 21 to 35 days post anthesis (DPA). Then, the seeds under HT and normal temperature (NT) environments were collected at 21 DPA, 28 DPA, and 35 DPA and subjected to whole-transcriptome sequencing. The phenotypic data showed that the seed germination percentage significantly increased, whereas SD decreased after HT treatment compared with NT, consistent with the results of previous studies. In total, 5128 mRNAs, 136 microRNAs (miRNAs), 273 long non-coding RNAs (lncRNAs), and 21 circularRNAs were found to be responsive to HT, and some of them were further verified through qRT-PCR. In particular, the known gibberellin (GA) biosynthesis gene TaGA20ox1 (TraesCS3D02G393900) was proved to be involved in HT-mediated dormancy by using the EMS-mutagenized wheat cultivar Jimai 22. Similarly, a novel gene TaCDPK21 (TraesCS7A02G267000) involved in the calcium signaling pathway was validated to be associated with HT-mediated dormancy by using the EMS mutant. Moreover, TaCDPK21 overexpression in Arabidopsis and functional complementarity tests supported the negative role of TaCDPK21 in SD. We also constructed a co-expression regulatory network based on differentially expressed mRNAs, miRNAs, and lncRNAs and found that a novel miR27319 was located at a key node of this regulatory network. Subsequently, using Arabidopsis and rice lines overexpressing miR27319 precursor or lacking miR27319 expression, we validated the positive role of miR27319 in SD and further preliminarily dissected the molecular mechanism of miR27319 underlying SD regulation through phytohormone abscisic acid and GA biosynthesis, catabolism, and signaling pathways. These findings not only broaden our understanding of the complex regulatory network of HT-mediated dormancy but also provide new gene resources for improving wheat PHS resistance to minimize PHS damage by using the molecular pyramiding approach.

Hypoxia Inhibits Cell Cycle Progression and Cell Proliferation in Brain Microvascular Endothelial Cells via the miR-212-3p/MCM2 Axis.

Hypoxia impairs blood-brain barrier (BBB) structure and function, causing pathophysiological changes in the context of stroke and high-altitude brain edema. Brain microvascular endothelial cells (BMECs) are major structural and functional elements of the BBB, and their exact role in hypoxia remains unknown. Here, we first deciphered the molecular events that occur in BMECs under 24 h hypoxia by whole-transcriptome sequencing assay. We found that hypoxia inhibited BMEC cell cycle progression and proliferation and downregulated minichromosome maintenance complex component 2 (Mcm2) expression. Mcm2 overexpression attenuated the inhibition of cell cycle progression and proliferation caused by hypoxia. Then, we predicted the upstream miRNAs of MCM2 through TargetScan and miRanDa and selected miR-212-3p, whose expression was significantly increased under hypoxia. Moreover, the miR-212-3p inhibitor attenuated the inhibition of cell cycle progression and cell proliferation caused by hypoxia by regulating MCM2. Taken together, these results suggest that the miR-212-3p/MCM2 axis plays an important role in BMECs under hypoxia and provide a potential target for the treatment of BBB disorder-related cerebrovascular disease.

Whole-Transcriptome Sequencing Combined with High-Dimensional Proteomic Technologies Reveals the Potential Value of miR-135b-5p as a Biomarker for Hepatocellular Carcinoma.

Hepatocellular carcinoma (HCC) is a disease with great heterogeneity and a high mortality rate. It is crucial to identify reliable biomarkers for diagnosis, prognosis, and treatment to improve clinical outcomes in patients with HCC. Alpha-fetoprotein (AFP) is not only a widely used biomarker in clinical practice but also plays a complicated role in HCC, and it has recently been considered to be related to immunotherapy. MicroRNAs (miRNAs) are regarded as key regulators and promising biomarkers of HCC. We investigated the role of an AFP-related miRNA, miR-135b-5p, in HCC progression. Identification of miR-135b-5p was performed based on a cohort of 65 HCC cases and the liver hepatocellular carcinoma cohort of The Cancer Genome Atlas (Asian people only). A combination of whole-transcriptome sequencing and high-dimensional proteomic technologies was used to study the role of miR-135b-5p in HCC. Upregulation of miR-135b-5p was detected in patients with HCC with high serum AFP levels (AFP > 400 ng/ml). Elevated miR-135b-5p expression was associated with adverse prognosis. We also identified the relevance between high miR-135b-5p expression and tumor-related pathological characteristics, such as Edmondson grade and vascular invasion. We revealed tyrosine kinase nonreceptor 1 as a potential target of miR-135b-5p. Additionally, the transcriptional start site of miR-135b-5p had potential binding sites for SRY-box transcription factor 9, and the stemness properties of tumor cells were more remarkable in HCC with the upregulation of miR-135b-5p. The molecular characterization of the miR-135b-5p-high group was similar to that of the HCC subclasses containing moderately and poorly differentiated tumors. Finally, gene signatures associated with improved clinical outcomes in immune checkpoint inhibitor therapy were upregulated in the miR-135b-5p-high group. miR-135b-5p could be a biomarker for predicting the prognosis and antiprogrammed cell death protein 1 monotherapy response in HCC.

Three-dimensional culture of rat ovarian granulosa cells shows increased SPP1 and FGF7 expression through the PI3K/Akt pathway.

Ovarian granulosa cells (OGCs)play an essential role in the regulation of follicular growth and development. However, previous studies of OGCs have concentrated on traditional 2D cultures. In the present study, we used the hanging drop culture method to culture rat OGCs (rOGCs) and assessed the effects of 3D conditions on their proliferation and gene expression profiles. Compared with those grown in 2D conditions, rOGCs grown in 3D cultures showed a significantly different spatial cell distribution and cell alignment under electron microscopy. In particular, rOGCs in 3D cultures showed abundant rough and microvilli-like structures on their cell surface. Here, we showed that these cells grew slowly following 3D culture; the G0/G1-phase increasedand the S- and G2/M-phases decreased. Using whole-transcriptome sequencing analysis, 501 genes were shown to have been significantly upregulated and 502 were shown to have been downregulated. Differentially expressed genes were most enriched in pathways involved in focal adhesion, MAPK, and PI3K/Akt signaling according to Kyoto Encyclopedia of Genes and Genomes pathway enrichment analysis. Western blotting revealed that SPP1 and FGF7 in the PI3K/Akt pathway were significantly upregulated following 3D culture. These findings improve our understanding of OGCs in real 3D environments in vivo and provide possible avenues for future research on OGCs.

Whole-transcriptome profiling across different developmental stages of Aedes albopictus (Diptera: Culicidae) provides insights into chitin-related non-coding RNA and competing endogenous RNA networks

BackgroundThe Asian tiger mosquito, Aedes albopictus, is one of the most invasive species and a vector of numerous arboviruses. The deleterious effects of long-term and inappropriate use of chemical pesticides have stimulated the exploration of new, environmentally friendly control strategies. Non-coding RNAs (ncRNAs) have been proven to participate in almost all biological processes of insects. MethodsIn this study, circular RNAs (circRNAs) and microRNAs (miRNAs) covering five developmental stages [egg, early larvae, late larvae, pupae, adult (female and male)] of A. albopictus were obtained using whole-transcriptome sequencing technology. Combined with long non-coding RNAs (lncRNAs) from previous research, circRNA/lncRNA‒miRNA‒mitochondrial RNA (mRNA) networks were constructed.ResultsA total of 1434 circRNAs and 208 miRNAs were identified. More differentially expressed circRNAs (DE circRNAs) and miRNAs (DE miRNAs) were found in the egg versus early larvae comparison group. Functional enrichment analysis demonstrated that most of the circRNA/lncRNA‒miRNA‒mRNA networks were involved in chitin metabolism. Hub genes of each circRNA/lncRNA‒miRNA‒mRNA network were screened out, which can be used as novel targets to disturb the molting process of A. albopictus.ConclusionsRegulatory relationships obtained from competing endogenous RNA (ceRNA) networks provide more information to manipulate the metamorphosis process and are helpful for developing effective and sustainable methods to control mosquitoes.Graphical

Whole genome sequencing provides comprehensive genetic testing in childhood B-cell acute lymphoblastic leukaemia

Childhood B-cell acute lymphoblastic leukaemia (B-ALL) is characterised by recurrent genetic abnormalities that drive risk-directed treatment strategies. Using current techniques, accurate detection of such aberrations can be challenging, due to the rapidly expanding list of key genetic abnormalities. Whole genome sequencing (WGS) has the potential to improve genetic testing, but requires comprehensive validation. We performed WGS on 210 childhood B-ALL samples annotated with clinical and genetic data. We devised a molecular classification system to subtype these patients based on identification of key genetic changes in tumour-normal and tumour-only analyses. This approach detected 294 subtype-defining genetic abnormalities in 96% (202/210) patients. Novel genetic variants, including fusions involving genes in the MAP kinase pathway, were identified. WGS results were concordant with standard-of-care methods and whole transcriptome sequencing (WTS). We expanded the catalogue of genetic profiles that reliably classify PAX5alt and ETV6::RUNX1-like subtypes. Our novel bioinformatic pipeline improved detection of DUX4 rearrangements (DUX4-r): a good-risk B-ALL subtype with high survival rates. Overall, we have validated that WGS provides a standalone, reliable genetic test to detect all subtype-defining genetic abnormalities in B-ALL, accurately classifying patients for the risk-directed treatment stratification, while simultaneously performing as a research tool to identify novel disease biomarkers.

Midline 1 interacting protein 1 promotes cancer metastasis through FOS-like 1-mediated matrix metalloproteinase 9 signaling in HCC.

Understanding the mechanisms of HCC progression and metastasis is crucial to improve early diagnosis and treatment. This study aimed to identify key molecular targets involved in HCC metastasis. Using whole-transcriptome sequencing of patients' HCCs, we identified and validated midline 1 interacting protein 1 (MID1IP1) as one of the most significantly upregulated genes in metastatic HCCs, suggesting its potential role in HCC metastasis. Clinicopathological correlation demonstrated that MID1IP1 upregulation significantly correlated with more aggressive tumor phenotypes and poorer patient overall survival rates. Functionally, overexpression of MID1IP1 significantly promoted the migratory and invasive abilities and enhanced the sphere-forming ability and expression of cancer stemness-related genes of HCC cells, whereas its stable knockdown abrogated these effects. Perturbation of MID1IP1 led to significant tumor shrinkage and reduced pulmonary metastases in an orthotopic liver injection mouse model and reduced pulmonary metastases in a tail-vein injection model in vivo . Mechanistically, SP1 transcriptional factor was found to be an upstream driver of MID1IP1 transcription. Furthermore, transcriptomic sequencing on MID1IP1-overexpressing HCC cells identified FOS-like 1 (FRA1) as a critical downstream mediator of MID1IP1. MID1IP1 upregulated FRA1 to subsequently promote its transcriptional activity and extracellular matrix degradation activity of matrix metalloproteinase MMP9, while knockdown of FRA1 effectively abolished the MID1IP1-induced migratory and invasive abilities. Our study identified MID1IP1 as a regulator in promoting FRA1-mediated-MMP9 signaling and demonstrated its role in HCC metastasis. Targeting MID1IP1-mediated FRA1 pathway may serve as a potential therapeutic strategy against HCC progression.

Dimorphic glioblastoma with glial and epithelioid phenotypes: Clonal evolution and immune selection.

Epithelioid glioblastoma is an unusual histologic variant of malignant glioma. The present study investigates both the genomic and transcriptomic determinants that may promote the development of this tumor. Whole-exome sequencing (WES) and whole-transcriptome sequencing (WTS) were performed on an epithelioid glioblastoma, along with a specific bioinformatic pipeline to generate electronic karyotyping and investigate the tumor immune microenvironment. Microdissected sections containing typical glioblastoma features and epithelioid morphology were analyzed separately using the same methodologies. An epithelioid glioblastoma, with immunopositivity for GFAP, Olig-2, and ATRX but negative for IDH-1 and p53, was identified. The tumor cell content from microdissection was estimated to be 85-90% for both histologic tumor components. WES revealed that both glioma and epithelioid sections contained identical point mutations in PTEN, RB1, TERT promoter, and TP53. Electronic karyotype analysis also revealed similar chromosomal copy number alterations, but the epithelioid component showed additional abnormalities that were not found in the glioblastoma component. The tumor immune microenvironments were strikingly different and WTS revealed high levels of transcripts from myeloid cells as well as M1 and M2 macrophages in the glioma section, while transcripts from CD4+ lymphocytes and NK cells predominated in the epithelioid section. Epithelioid glioblastoma may be genomically more unstable and oncogenically more advanced, harboring an increased number of mutations and karyotype abnormalities, compared to typical glioblastomas. The tumor immune microenvironment is also different.

Comprehensive Analysis of Whole-Transcriptome Profiles in Response to Acute Hypersaline Challenge in Chinese Razor Clam Sinonovacula constricta.

The Chinese razor clam (Sinonovacula constricta) is an important for Chinese aquaculture marine bivalve that naturally occurs across intertidal and estuarine areas subjected to significant changes in salinity level. However, the information on the molecular mechanisms related to high salinity stress in the species remain limited. In this study, nine gill samples of S. constricta treated with 20, 30, and 40 ppt salinity for 24 h were used for whole-transcriptome RNA sequencing, and a regulatory network of competing endogenous RNAs (ceRNAs) was constructed to better understand the mechanisms responsible for adaptation of the species to high salinity. A total of 83,262 lncRNAs, 52,422 mRNAs, 2890 circRNAs, and 498 miRNAs were identified, and 4175 of them displayed differential expression pattern among the three groups examined. The KEGG analyses of differentially expressed RNAs evidenced that amino acid synthesis and membrane transport were the dominant factors involved in the adaptation of the Chinese razor clam to acute salinity increase, while lipid metabolism and signaling played only a supporting role. In addition, lncRNA/circRNA-miRNA-mRNA regulatory networks (ceRNA network) showed clearly regulatory relationships among different RNAs. Moreover, the expression of four candidate genes, including tyrosine aminotransferase (TAT), hyaluronidase 4 (HYAL4), cysteine sulfinic acid decarboxylase (CSAD), and ∆1-pyrroline-5-carboxylate synthase (P5CS) at different challenge time were detected by qRT-PCR. The expression trend of TAT and HYAL4 was consistent with that of the ceRNA network, supporting the reliability of established network. The expression of TAT, CSAD, and P5CS were upregulated in response to increased salinity. This might be associated with increased amino acid synthesis rate, which seems to play an essential role in adaptation of the species to high salinity stress. In contrast, the expression level of HYAL4 gene decreased in response to elevated salinity level, which is associated with reduction Hyaluronan hydrolysis to help maintain water in the cell. Our findings provide a very rich reference for understanding the important role of ncRNAs in the salinity adaptation of shellfish. Moreover, the acquired information may be useful for optimization of the artificial breeding of the Chinese razor clam under aquaculture conditions.

Transcriptome-wide analysis reveals the molecular mechanisms of cannabinoid type II receptor agonists in cardiac injury induced by chronic psychological stress.

Background: Growing evidence has supported that chronic psychological stress would cause heart damage, However the mechanisms involved are not clear and effective interventions are insufficient. Cannabinoid type 2 receptor (CB2R) can be a potential treatment for cardiac injury. This study is aimed to investigate the protective mechanism of CB2R agonist against chronic psychological stress-induced cardiac injury. Methods: A mouse chronic psychological stress model was constructed based on a chronic unpredictable stress pattern. Mice were performed a three-week psychological stress procedure, and cardiac tissues of them were collected for whole-transcriptome sequencing. Overlap analysis was performed on differentially expressed mRNAs (DE-mRNAs) and ER stress-related genes (ERSRGs), and bioinformatic methods were used to predict the ceRNA networks and conduct pathway analysis. The expressions of the DE-ERSRGs were validated by RT-qPCR. Results: In the comparison of DE mRNA in Case group, Control group and Treatment group, three groups of ceRNA networks and ceRNA (circ) networks were constructed. The DE-mRNAs were mainly enriched in chromatid-relevant terms and Hematopoietic cell lineage pathway. Additionally, 13 DE-ERSRGs were obtained by the overlap analysis, which were utilized to establish a ceRNA network with 15 nodes and 14 edges and a ceRNA (circ) network with 23 nodes and 28 edges. Furthermore, four DE-ERSRGs (Cdkn1a, Atf3, Fkbp5, Gabarapl1) in the networks were key, which were mainly enriched in response to extracellular stimulus, response to nutrient levels, cellular response to external stimulus, and FoxO signaling pathway. Finally, the RT-qPCR results showed almost consistent expression patterns of 13 DE-ERSRGs between the transcriptome and tissue samples. Conclusion: The findings of this study provide novel insights into the molecular mechanisms of chronic psychological stress-induced cardiac diseases and reveal novel targets for the cardioprotective effects of CB2R agonists.

Circular RNA circEMB promotes osteosarcoma progression and metastasis by sponging miR-3184-5p and regulating EGFR expression

BackgroundOsteosarcoma (OSA) is the most prevalent type of bone cancer with a high rate of metastasis. Circular RNAs (CircRNAs) play an essential role in multiple aspects of tumour biology. This study aimed to elucidate the role of circEMB in OSA.MethodscircRNAs related to OSA invasion were identified via RNA sequencing and qRT-PCR. The relationship between circEMB levels and clinicopathological features of OSA was examined using the clinical specimens and data of 53 patients with OSA. Several in vivo and in vitro experiments, including intravital imaging, whole-transcriptome sequencing, transwell assay, flow cytometry, dual-luciferase reporter assay, RIP assay, RNA pull-down assay and RNA-FISH, were performed to examine the effects of circEMB on the malignant behaviour of OSA.ResultsA novel circRNA, named circEMB (hsa_circ_001310), was identified in this study. circEMB can promote the malignant behaviour of OSA. In vitro experiments revealed that circEMB knockdown decreased cell proliferation, inhibited tumour invasion and metastasis; increased apoptosis and resulted in G1/S phase arrest. In vivo experiments revealed that circEMB knockdown inhibited tumour growth and metastasis in xenograft-bearing mice. Mechanistically, circEMB affects the malignant behaviour of OSA by mediating EGFR as an miR-3184-5p sponge. In addition, the circEMB/miR-3184-5p/EGFR axis modulates methotrexate (MTX) resistance in OSA.ConclusionsCircEMB plays a critical role in promoting cancer via the miR-3184-5p/EGFR pathway, indicating that circEMB may serve as a therapeutic target for OSA.

Genomic Characterization of Rare Primary Cardiac Sarcoma Entities.

Primary cardiac sarcomas are considered rare malignant entities associated with poor prognosis. In fact, knowledge regarding their gene signature and possible treatments is still limited. In our study, whole-transcriptome sequencing on formalin-fixed paraffin-embedded (FFPE) samples from one cardiac osteosarcoma and one cardiac leiomyosarcoma was performed, to investigate their mutational profiles and to highlight differences and/or similarities to other cardiac histotypes. Both cases have been deeply detailed from a pathological point of view. The osteosarcoma sample presented mutations involving ATRX, ERCC5, and COL1A1, while the leiomyosarcoma case showed EXT2, DNM2, and PSIP1 alterations. Altered genes, along with the most differentially expressed genes in the leiomyosarcoma or osteosarcoma sample versus the cardiac angiosarcomas and intimal sarcomas (e.g., YAF2, PAK5, and CRABP1), appeared to be associated with cell growth, proliferation, apoptosis, and the repair of DNA damage, which are key mechanisms involved in tumorigenesis. Moreover, a distinct gene expression profile was detected in the osteosarcoma sample when compared to other cardiac sarcomas. For instance, WIF1, a marker of osteoblastic differentiation, was upregulated in our bone tumor. These findings pave the way for further studies on these entities, in order to identify targeted therapies and, therefore, improve patients' prognoses.

CircAGK regulates high dihydrotestosterone-induced apoptosis in DPCs through the miR-3180-5p/BAX axis.

The incidence of androgen alopecia (AGA), also known as seborrheic alopecia, has surged in recent years, and onset is occurring at younger ages. Dermal papilla cells (DPCs) are key to maintaining hair cycling, and apoptosis-driven processes in DPCs are closely related to hair follicle regeneration. Circular RNAs (circRNAs) are widely present in the human body and are closely related to the occurrence and development of many diseases. Currently, the biological functions of circRNAs in AGA are largely unknown. Whole-transcriptome sequencing was used to screen differential circRNA expression profiles between AGA patients and non-AGA patients. We found that hsa_circ_0002980 (circAGK) was significantly highly expressed in the AGA group. CircAGK promoted DPC apoptosis in the presence of high dihydrotestosterone (DHT) (15 nmol/L). By regulating the miR-3180-5p/BAX axis, circAGK promotes DPC apoptosis in a high DHT environment in vitro and inhibits hair growth in AGA mice in vivo, indicating that circAGK is a potential target for the clinical treatment of AGA.

Integrated whole transcriptome analysis for the crucial regulators and functional pathways related to cardiac fibrosis in rats.

Cardiac fibrosis has gradually gained significance in the field of cardiovascular disease; however, its specific pathogenesis remains unclear. This study aims to establish the regulatory networks based on whole-transcriptome RNA sequencing analyses and reveal the underlying mechanisms of cardiac fibrosis. An experimental model of myocardial fibrosis was induced using the chronic intermittent hypoxia (CIH) method. Expression profiles of long non-coding RNA (lncRNA), microRNA (miRNA), and messenger RNA (mRNA) were acquired from right atrial tissue samples of rats. Differentially expressed RNAs (DERs) were identified, and functional enrichment analysis was performed. Moreover, a protein-protein interaction (PPI) network and competitive endogenous RNA (ceRNA) regulatory network that are related to cardiac fibrosis were constructed, and the relevant regulatory factors and functional pathways were identified. Finally, the crucial regulators were validated using qRT-PCR. DERs, including 268 lncRNAs, 20 miRNAs, and 436 mRNAs, were screened. Further, 18 relevant biological processes, such as "chromosome segregation, " and 6 KEGG signaling pathways, such as "cell cycle, " were significantly enriched. The regulatory relationship of miRNA-mRNA-KEGG pathways showed eight overlapping disease pathways, including "pathways in cancer." In addition, crucial regulatory factors, such as Arnt2, WNT2B, GNG7, LOC100909750, Cyp1a1, E2F1, BIRC5, and LPAR4, were identified and verified to be closely related to cardiac fibrosis. This study identified the crucial regulators and related functional pathways in cardiac fibrosis by integrating the whole transcriptome analysis in rats, which might provide novel insights into the pathogenesis of cardiac fibrosis.

Transformer 2β regulates the alternative splicing of cell cycle regulatory genes to promote the malignant phenotype of ovarian cancer.

Late-stage ovarian cancer (OC) has a poor prognosis and a high metastasis rate, but the underlying molecular mechanism is unclear. RNA binding proteins (RBPs) play important roles in posttranscriptional regulation in the contexts of neoplasia and tumor metastasis. In this study, we explored the molecular functions of a canonical RBP, Transformer 2β homolog (TRA2B), in cancer cells. TRA2B knockdown in HeLa cells and subsequent whole-transcriptome RNA sequencing (RNA-seq) analysis revealed the TRA2B-regulated alternative splicing (AS) profile. We disrupted TRA2B expression in epithelial OC cells and performed a series of experiments to confirm the resulting effects on OC cell proliferation, apoptosis and invasion. TRA2B-regulated AS was tightly associated with the mitotic cell cycle, apoptosis and several cancer pathways. Moreover, the expression of hundreds of genes was regulated by TRA2B, and these genes were enriched in the functions of cell proliferation, cell adhesion and angiogenesis, which are related to the malignant phenotype of OC. By integrating the alternatively spliced and differentially expressed genes, we found that AS events and gene expression were regulated independently. We then explored and validated the oncogenic functions of TRA2B by knocking down its expression in OC cells. The high TRA2B expression was associated with poor prognosis in patients with OC. In ovarian tissue, TRA2B expression showed a gradual increasing trend with increasing malignancy. We demonstrated the important roles of TRA2B in ovarian neoplasia and aggressive OC behaviors and identified the underlying molecular mechanisms, facilitating the targeted treatment of OC.