Demystifying emerging bulk RNA-Seq applications: the application and utility of bioinformatic methodology

Liquiritin Alleviates Inflammation in Lipopolysaccharide-Induced Human Corneal Epithelial Cells

PURPOSE

This research was designed to elucidate the anti-inflammatory impacts of liquiritin on lipopolysaccharide (LPS)-activated human corneal epithelial cells (HCECs).

METHODS

The Cell Counting kit-8 (CCK-8) assay was adopted to assess cell viability. The enzyme-linked immunosorbent assay (ELISA) was used to detect the secretion levels of the proinflammatory cytokines IL-6, IL-8, and TNF-α. Transcriptome analysis was conducted to identify the genes that exhibited differential expression between different treatment. The model group included cells treated with LPS (10 µg/mL), the treatment group comprised cells treated with liquiritin (80 µM) and LPS (10 µg/mL), and the control group consisted of untreated cells. To further validate the expression levels of the selected genes, including CSF2, CXCL1, CXCL2, CXCL8, IL1A, IL1B, IL24, IL6, and LTB, quantitative real-time PCR was performed. The expression of proteins related to the Akt/NF-κB signaling pathway was assessed through western blot analysis. NF-κB nuclear translocation was evaluated through immunofluorescence staining.

RESULTS

The secretion of IL-6, IL-8, and TNF-α in LPS-induced HCECs was significantly downregulated by liquiritin. Based on the transcriptome analysis, the mRNA expression of pro-inflammatory cytokines, namely IL-6, IL-8, IL-1β, IL-24, TNF-α, and IL-1α was overproduced by LPS stimulation, and suppressed after liquiritin treatment. Furthermore, the Western blot results revealed a remarkable reduction in the phosphorylation degrees of NF-κB p65, IκB, and Akt upon treatment with liquiritin. Additionally, immunofluorescence analysis confirmed liquiritin's inhibition of LPS-induced p65 nuclear translocation.

CONCLUSIONS

Collectively, these findings imply that liquiritin suppresses the expression of proinflammatory cytokines, and the anti-inflammatory impacts of liquiritin may be caused by its repression of the Akt/NF-κB signaling pathway in LPS-induced HCECs. These data indicate that liquiritin could provide a potential therapeutic application for inflammation-associated corneal diseases.

OmicVerse: a framework for bridging and deepening insights across bulk and single-cell sequencing

Single-cell sequencing is frequently affected by "omission" due to limitations in sequencing throughput, yet bulk RNA-seq may contain these ostensibly "omitted" cells. Here, we introduce the single cell trajectory blending from Bulk RNA-seq (BulkTrajBlend) algorithm, a component of the OmicVerse suite that leverages a Beta-Variational AutoEncoder for data deconvolution and graph neural networks for the discovery of overlapping communities. This approach effectively interpolates and restores the continuity of "omitted" cells within single-cell RNA sequencing datasets. Furthermore, OmicVerse provides an extensive toolkit for both bulk and single cell RNA-seq analysis, offering seamless access to diverse methodologies, streamlining computational processes, fostering exquisite data visualization, and facilitating the extraction of significant biological insights to advance scientific research.

Genetic variation mining of the Chinese mitten crab (Eriocheir sinensis) based on transcriptome data from public databases

At present, public databases house an extensive repository of transcriptome data, with the volume continuing to grow at an accelerated pace. Utilizing these data effectively is a shared interest within the scientific community. In this study, we introduced a novel strategy that harnesses SNPs and InDels identified from transcriptome data, combined with sample metadata from databases, to effectively screen for molecular markers correlated with traits. We utilized 228 transcriptome datasets of Eriocheir sinensis from the NCBI database and employed the Genome Analysis Toolkit software to identify 96 388 SNPs and 20 645 InDels. Employing the genome-wide association study analysis, in conjunction with the gender information from databases, we identified 3456 sex-biased SNPs and 639 sex-biased InDels. The KOG and KEGG annotations of the sex-biased SNPs and InDels revealed that these genes were primarily involved in the metabolic processes of E. sinensis. Combined with SnpEff annotation and PCR experimental validation, a highly sex-biased SNP located in the Kelch domain containing 4 (Klhdc4) gene, CHR67-6415071, was found to alter the splicing sites of Klhdc4, generating two splice variants, Klhdc4_a and Klhdc4_b. Additionally, Klhdc4 exhibited robust expression across the ovaries, testes, and accessory glands. The sex-biased SNPs and InDels identified in this study are conducive to the development of unisexual cultivation methods for E. sinensis, and the alternative splicing event caused by the sex-biased SNP in Klhdc4 may serve as a potential mechanism for sex regulation in E. sinensis. The analysis strategy employed in this study represents a new direction for the rational exploitation and utilization of transcriptome data in public databases.

Novel miRNA-inducing drugs enable differentiation of retinoic acid-resistant neuroblastoma cells

Tumor cell heterogeneity in neuroblastoma, a pediatric cancer arising from neural crest-derived progenitor cells, poses a significant clinical challenge. In particular, unlike adrenergic (ADRN) neuroblastoma cells, mesenchymal (MES) cells are resistant to chemotherapy and retinoid therapy and thereby significantly contribute to relapses and treatment failures. Previous research suggested that overexpression or activation of miR-124, a neurogenic microRNA with tumor suppressor activity, can induce the differentiation of retinoic acid-resistant neuroblastoma cells. Leveraging our established screen for miRNA modulatory small molecules, we validated PP121, a dual inhibitor of tyrosine and phosphoinositide kinases, as a robust inducer of miR-124. A combination of PP121 and miR-132-inducing bufalin synergistically arrests proliferation, induces differentiation, and prolongs the survival of differentiated MES SK-N-AS cells for 8 weeks. RNA- seq and deconvolution analyses revealed a collapse of the ADRN core regulatory circuitry (CRC) and the emergence of novel CRCs associated with chromaffin cells and Schwann cell precursors. Using a similar protocol, we differentiated and maintained other MES neuroblastoma, as well as glioblastoma cells, over 16 weeks. In conclusion, our novel protocol suggests a promising treatment for therapy-resistant cancers of the nervous system. Moreover, these long-lived, differentiated cells provide valuable models for studying mechanisms underlying differentiation, maturation, and senescence.

Crosstalk among Oxidative Stress, Autophagy, and Apoptosis in the Protective Effects of Ginsenoside Rb1 on Brain Microvascular Endothelial Cells: A Mixed Computational and Experimental Study

OBJECTIVE

Brain microvascular endothelial cells (BMECs) were found to shift from their usually inactive state to an active state in ischemic stroke (IS) and cause neuronal damage. Ginsenoside Rb1 (GRb1), a component derived from medicinal plants, is known for its pharmacological benefits in IS, but its protective effects on BMECs have yet to be explored. This study aimed to investigate the potential protective effects of GRb1 on BMECs.

METHODS

An in vitro oxygen-glucose deprivation/reperfusion (OGD/R) model was established to mimic ischemia-reperfusion (I/R) injury. Bulk RNA-sequencing data were analyzed by using the Human Autophagy Database and various bioinformatic tools, including gene set enrichment analysis (GSEA), Gene Ontology (GO) classification and enrichment analysis, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis, protein-protein interaction network analysis, and molecular docking. Experimental validation was also performed to ensure the reliability of our findings.

RESULTS

Rb1 had a protective effect on BMECs subjected to OGD/R injury. Specifically, GRb1 was found to modulate the interplay between oxidative stress, apoptosis, and autophagy in BMECs. Key targets such as sequestosome 1 (SQSTM1/p62), autophagy related 5 (ATG5), and hypoxia-inducible factor 1-alpha (HIF-1α) were identified, highlighting their potential roles in mediating the protective effects of GRb1 against IS-induced damage.

CONCLUSION

GRbl protects BMECs against OGD/R injury by influencing oxidative stress, apoptosis, and autophagy. The identification of SQSTM1/p62, ATG5, and HIF-1α as promising targets further supports the potential of GRb1 as a therapeutic agent for IS, providing a foundation for future research into its mechanisms and applications in IS treatment.

Advances in long-read single-cell transcriptomics

Long-read single-cell transcriptomics (scRNA-Seq) is revolutionizing the way we profile heterogeneity in disease. Traditional short-read scRNA-Seq methods are limited in their ability to provide complete transcript coverage, resolve isoforms, and identify novel transcripts. The scRNA-Seq protocols developed for long-read sequencing platforms overcome these limitations by enabling the characterization of full-length transcripts. Long-read scRNA-Seq techniques initially suffered from comparatively poor accuracy compared to short read scRNA-Seq. However, with improvements in accuracy, accessibility, and cost efficiency, long-reads are gaining popularity in the field of scRNA-Seq. This review details the advances in long-read scRNA-Seq, with an emphasis on library preparation protocols and downstream bioinformatics analysis tools.

Construction of an Immunogenic Cell Death-Related Gene Signature and Genetic Subtypes for Predicting Prognosis, Immune Microenvironments, and Drug Sensitivity in Hepatocellular Carcinoma

Purpose

Immunogenic cell death (ICD) is a type of regulated cell death that modifies the immune response by releasing DAMPs or danger signals. Herein, we aimed to develop an ICD-related predictive model for patients with hepatocellular carcinoma (HCC) and investigate its applicability for predicting prognostic outcomes and immunotherapeutic responses.

Methods

Differentially expressed genes of ICD were identified in the HCC and normal liver samples. A prognostic risk model and a nomogram containing clinicopathological features were created. To validate the effectiveness of the model, an external dataset was used. Clinical characteristics, prognosis, tumor mutation burden, immune microenvironments, biological function and chemotherapeutic drug sensitivity were evaluated for different genetic subtypes and risk groups.

Results

A total of 35 ICD-related genes (ICDRGs) were identified between HCC and normal samples, 11 of which were significantly associated with overall survival (OS) in HCC patients. Four different genetic subtypes were formed and eight ICDRGs were selected to develop a risk prognostic model. The risk scores were shown to be an independent prognostic factor for HCC and positively correlated with pathological severity. Patients in the high-risk group had a higher frequency of TP53 mutations, increased expression of immune checkpoints and human leukocyte antigen genes. The inhibitory concentrations of chemotherapeutic drugs differed in different populations.

Conclusion

In this study, we developed an ICDRG risk model and demonstrated its applicability in predicting survival outcomes, immune and chemotherapeutic responses in HCC patients. ICDRGs are expected to be used as novel biomarkers in the medical decision-making of HCC.

Integrated analysis of single-cell RNA-seq and bulk RNA-seq reveals immune suppression subtypes and establishes a novel signature for determining the prognosis in lung adenocarcinoma

BACKGROUND

Lung adenocarcinoma (LUAD) is the most common histological type of lung cancer with lower survival rates. Recent advancements in targeted therapies and immunotherapies targeting immune checkpoints have achieved remarkable success, there is still a large percentage of LUAD that lacks available therapeutic options. Due to tumor heterogeneity, the diagnosis and treatment of LUAD are challenging. Exploring the biology of LUAD and identifying new biomarker and therapeutic targets options are essential.

METHOD

We performed single-cell RNA sequencing (scRNA-seq) of 6 paired primary and adjacent LUAD tissues, and integrative omics analysis of the scRNA-seq, bulk RNA-seq and whole-exome sequencing data revealed molecular subtype characteristics. Our experimental results confirm that CDC25C gene can serve as a potential marker for poor prognosis in LUAD.

RESULTS

We investigated aberrant gene expression in diverse cell types in LUAD via the scRNA-seq data. Moreover, multi-omics clustering revealed four subgroups defined by transcriptional profile and molecular subtype 4 (MS4) with poor survival probability, and immune cell infiltration signatures revealed that MS4 tended to be the immunosuppressive subtype. Our study revealed that the CDC25C gene can be a distinct prognostic biomarker that indicates immune infiltration levels and response to immunotherapy in LUAD patients. Our experimental results concluded that CDC25C expression affects lung cancer cell invasion and migration, might play a key role in regulating Epithelial-Mesenchymal Transition (EMT) pathways.

CONCLUSIONS

Our multi-omics result revealed a comprehensive set of molecular attributes associated with prognosis-related genes in LUAD at the cellular and tissue level. Identification of a subtype of immunosuppressive TME and prognostic signature for LUAD. We identified the cell cycle regulation gene CDC25C affects lung cancer cell invasion and migration, which can be used as a potential biomarker for LUAD.

An array of signal-specific MoYpd1 isoforms determines full virulence in the pathogenic fungus Magnaporthe oryzae

Magnaporthe oryzae is placed first on a list of the world's top ten plant pathogens with the highest scientific and economic importance. The locus MGG_07173 occurs only once in the genome of M. oryzae and encodes the phosphotransfer protein MoYpd1p, which plays an important role in the high osmolarity glycerol (HOG) signaling pathway for osmoregulation. Originating from this locus, at least three MoYPD1 isoforms are produced in a signal-specific manner. The transcript levels of these MoYPD1-isoforms were individually affected by external stress. Salt (KCI) stress raised MoYPD1_T0 abundance, whereas osmotic stress by sorbitol elevates MoYPD1_T1 levels. In line with this, signal-specific nuclear translocation of green fluorescent protein-fused MoYpd1p isoforms in response to stress was observed. Mutant strains that produce only one of the MoYpd1p isoforms are less virulent, suggesting a combination thereof is required to invade the host successfully. In summary, we demonstrate signal-specific production of MoYpd1p isoforms that individually increase signal diversity and orchestrate virulence in M. oryzae.

Construction of circ_0071922-miR-15a-5p-mRNA network in intervertebral disc degeneration by RNA-sequencing

Background

Low back pain (LBP) is the main factor of global disease burden. Intervertebral disc degeneration (IVDD) has long been known as the leading reason of LBP. Increasing studies have verified that circular RNAs (circRNAs)-microRNAs (miRNAs)-mRNAs network is widely involved in the pathological processes of IVDD. However, no study was made to demonstrate the circRNAs-mediated ferroptosis, oxidative stress, extracellular matrix metabolism, and immune response in IVDD.

Methods

We collected 3 normal and 3 degenerative nucleus pulposus tissues to conduct RNA-sequencing to identify the key circRNAs and miRNAs in IVDD. Bioinformatics analysis was then conducted to construct circRNAs-miRNAs-mRNAs interaction network associated with ferroptosis, oxidative stress, extracellular matrix metabolism, and immune response. We also performed animal experiments to validate the therapeutic effects of key circRNAs in IVDD.

Results

We found that circ_0015435 was most obviously upregulated and circ_0071922 was most obviously downregulated in IVDD using RNA-sequencing. Then we observed that hsa-miR-15a-5p was the key downstream of circ_0071922, and hsa-miR-15a-5p was the top upregulated miRNA in IVDD. Bioinformatics analysis was conducted to predict that 56 immunity-related genes, 29 ferroptosis-related genes, 23 oxidative stress-related genes and 8 ECM-related genes are the targets mRNAs of hsa-miR-15a-5p. Then we constructed a ceRNA network encompassing 24 circRNAs, 6 miRNAs, and 101 mRNAs. Additionally, we demonstrated that overexpression of circ_0071922 can alleviate IVDD progression in a rat model.

Conclusions

The findings of this study suggested that circ_0071922-miR-15a-5p-mRNA signaling network might affect IVDD by modulating the nucleus pulposus cells ferroptosis, oxidative stress, ECM metabolism, and immune response, which is an effective therapeutic targets of IVDD.

Identification and validation of CCR5 linking keloid with atopic dermatitis through comprehensive bioinformatics analysis and machine learning

There is sufficient evidence indicating that keloid is strongly associated with atopic dermatitis (AD) across ethnic groups. However, the molecular mechanism underlying the association is not fully understood. The aim of this study is to discover the underlying mechanism of the association between keloid and AD by integrating comprehensive bioinformatics techniques and machine learning methods. The gene expression profiles of keloid and AD were downloaded from the Gene Expression Omnibus (GEO) database. A total of 449 differentially expressed genes (DEGs) were found to be shared in keloid and AD using the training datasets of GEO (GSE158395 and GSE121212). The hub genes were identified using the protein-protein interaction network and Cytoscape software. 20 of the most significant hub genes were selected, which were mainly involved in the regulation of the inflammatory and immune response. Through two machine learning algorithms of LASSO and SVM-RFE, CCR5 was identified as the most important key gene. Subsequently, upregulated CCR5 gene expression was confirmed in validation GEO datasets (GSE188952 and GSE32924) and clinical samples of keloid and AD. Immune infiltration analysis showed that T helper (Th) 1, 2 and 17 cells were significantly enriched in the microenvironment of both keloid and AD. Positive correlations were found between CCR5 and Th1, Th2 and Th17 cells. Finally, two TFs of CCR5, NR3C2 and YY1, were identified, both of which were downregulated in keloid and AD tissues. Our study firstly reveals that keloid and AD shared common inflammatory and immune pathways. Moreover, CCR5 plays a key role in the pathogenesis association between keloid and AD. The common pathways and key genes may shed light on further mechanism research and targeted therapy, and may provide therapeutic interventions of keloid with AD.

ClusTrast: a short read de novo transcript isoform assembler guided by clustered contigs

BACKGROUND

Transcriptome assembly from RNA-sequencing data in species without a reliable reference genome has to be performed de novo, but studies have shown that de novo methods often have inadequate ability to reconstruct transcript isoforms. We address this issue by constructing an assembly pipeline whose main purpose is to produce a comprehensive set of transcript isoforms.

RESULTS

We present the de novo transcript isoform assembler ClusTrast, which takes short read RNA-seq data as input, assembles a primary assembly, clusters a set of guiding contigs, aligns the short reads to the guiding contigs, assembles each clustered set of short reads individually, and merges the primary and clusterwise assemblies into the final assembly. We tested ClusTrast on real datasets from six eukaryotic species, and showed that ClusTrast reconstructed more expressed known isoforms than any of the other tested de novo assemblers, at a moderate reduction in precision. For recall, ClusTrast was on top in the lower end of expression levels (<15% percentile) for all tested datasets, and over the entire range for almost all datasets. Reference transcripts were often (35-69% for the six datasets) reconstructed to at least 95% of their length by ClusTrast, and more than half of reference transcripts (58-81%) were reconstructed with contigs that exhibited polymorphism, measuring on a subset of reliably predicted contigs. ClusTrast recall increased when using a union of assembled transcripts from more than one assembly tool as primary assembly.

CONCLUSION

We suggest that ClusTrast can be a useful tool for studying isoforms in species without a reliable reference genome, in particular when the goal is to produce a comprehensive transcriptome set with polymorphic variants.

TP53 to mediate immune escape in tumor microenvironment: an overview of the research progress

Increasing evidence suggests that key cancer-causing driver genes continue to exert a sustained influence on the tumor microenvironment (TME), highlighting the importance of immunotherapeutic targeting of gene mutations in governing tumor progression. TP53 is a prominent tumor suppressor that encodes the p53 protein, which controls the initiation and progression of different tumor types. Wild-type p53 maintains cell homeostasis and genomic instability through complex pathways, and mutant p53 (Mut p53) promotes tumor occurrence and development by regulating the TME. To date, it has been wildly considered that TP53 is able to mediate tumor immune escape. Herein, we summarized the relationship between TP53 gene and tumors, discussed the mechanism of Mut p53 mediated tumor immune escape, and summarized the progress of applying p53 protein in immunotherapy. This study will provide a basic basis for further exploration of therapeutic strategies targeting p53 protein.

Emerging applications of single-cell profiling in precision medicine of atherosclerosis

Atherosclerosis is a chronic, progressive, inflammatory disease that occurs in the arterial wall. Despite recent advancements in treatment aimed at improving efficacy and prolonging survival, atherosclerosis remains largely incurable. In this review, we discuss emerging single-cell sequencing techniques and their novel insights into atherosclerosis. We provide examples of single-cell profiling studies that reveal phenotypic characteristics of atherosclerosis plaques, blood, liver, and the intestinal tract. Additionally, we highlight the potential clinical applications of single-cell analysis and propose that combining this approach with other techniques can facilitate early diagnosis and treatment, leading to more accurate medical interventions.

Nanostrategy of Targeting at Embryonic Trophoblast Cells Using CuO Nanoparticles for Female Contraception

Effective contraceptives have been comprehensively adopted by women to prevent the negative consequences of unintended pregnancy for women, families, and societies. With great contributions of traditional hormonal drugs and intrauterine devices (IUDs) to effective female contraception by inhibiting ovulation and deactivating sperm, their long-standing side effects on hormonal homeostasis and reproductive organs for females remain concerns. Herein, we proposed a nanostrategy for female contraceptives, inducing embryonic trophoblast cell death using nanoparticles to prevent embryo implantation. Cupric oxide nanoparticles (CuO NPs) were adopted in this work to verify the feasibility of the nanostrategy and its contraceptive efficacy. We carried out the in vitro assessment on the interaction of CuO NPs with trophoblast cells using the HTR8/SVneo cell line. The results showed that the CuO NPs were able to be preferably uptaken into cells and induced cell damage via a variety of pathways including oxidative stress, mitochondrial damage, DNA damage, and cell cycle arrest to induce cell death of apoptosis, ferroptosis, and cuproptosis. Moreover, the key regulatory processes and the key genes for cell damage and cell death caused by CuO NPs were revealed by RNA-Seq. We also conducted in vivo experiments using a rat model to examine the contraceptive efficacy of both the bare CuO NPs and the CuO/thermosensitive hydrogel nanocomposite. The results demonstrated that the CuO NPs were highly effective for contraception. There was no sign of disrupting the homeostasis of copper and hormone, or causing inflammation and organ damage in vivo. In all, this nanostrategy exhibited huge potential for contraceptive development with high biosafety, efficacy, clinical translation, nonhormonal style, and on-demand for women.

Identification of hub genes based on integrated analysis of single-cell and microarray transcriptome in patients with pulmonary arterial hypertension

BACKGROUND

Pulmonary arterial hypertension (PAH) is a devastating chronic cardiopulmonary disease without an effective therapeutic approach. The underlying molecular mechanism of PAH remains largely unexplored at single-cell resolution.

METHODS

Single-cell RNA sequencing (scRNA-seq) data from the Gene Expression Omnibus (GEO) database (GSE210248) was included and analyzed comprehensively. Additionally, microarray transcriptome data including 15 lung tissue from PAH patients and 11 normal samples (GSE113439) was also obtained. Seurat R package was applied to process scRNA-seq data. Uniform manifold approximation and projection (UMAP) was utilized for dimensionality reduction and cluster identification, and the SingleR package was performed for cell annotation. FindAllMarkers analysis and ClusterProfiler package were applied to identify differentially expressed genes (DEGs) for each cluster in GSE210248 and GSE113439, respectively. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genome (KEGG) were used for functional enrichment analysis of DEGs. Microenvironment Cell Populations counter (MCP counter) was applied to evaluate the immune cell infiltration. STRING was used to construct a protein-protein interaction (PPI) network of DEGs, followed by hub genes selection through Cytoscape software and Veen Diagram.

RESULTS

Nineteen thousand five hundred seventy-six cells from 3 donors and 21,896 cells from 3 PAH patients remained for subsequent analysis after filtration. A total of 42 cell clusters were identified through UMAP and annotated by the SingleR package. 10 cell clusters with the top 10 cell amounts were selected for consequent analysis. Compared with the control group, the proportion of adipocytes and fibroblasts was significantly reduced, while CD8+ T cells and macrophages were notably increased in the PAH group. MCP counter revealed decreased distribution of CD8+ T cells, cytotoxic lymphocytes, and NK cells, as well as increased infiltration of monocytic lineage in PAH lung samples. Among 997 DEGs in GSE113439, module 1 with 68 critical genes was screened out through the MCODE plug-in in Cytoscape software. The top 20 DEGs in each cluster of GSE210248 were filtered out by the Cytohubba plug-in using the MCC method. Eventually, WDR43 and GNL2 were found significantly increased in PAH and identified as the hub genes after overlapping these DEGs from GSE210248 and GSE113439.

CONCLUSION

WDR43 and GNL2 might provide novel insight into revealing the new molecular mechanisms and potential therapeutic targets for PAH.

Distinct regulatory machineries underlying divergent chromatin landscapes distinguish innate lymphoid cells from T helper cells

Innate lymphoid cells (ILCs), as the innate counterpart of CD4+ T helper (Th) cells, play crucial roles in maintaining tissue homeostasis. While the ILC subsets and their corresponding Th subsets demonstrate significant similarities in core programming related to effector function and regulatory mechanisms, their principal distinctions, given their innate and adaptive lymphocyte nature, remain largely unknown. In this study, we have employed an integrative analysis of 294 bulk RNA-sequencing results across all ILC and Th subsets, using scRNA-seq algorithms. Consequently, we identify two genesets that predominantly differentiate ILCs from Th cells, as well as three genesets that distinguish various immune responses. Furthermore, through chromatin accessibility analysis, we find that the ILC geneset tends to rely on specific transcriptional regulation at promoter regions compared with the Th geneset. Additionally, we observe that ILCs and Th cells are under differential transcriptional regulation. For example, ILCs are under stronger regulation by multiple transcription factors, including RORα, GATA3, and NF-κB. Otherwise, Th cells are under stronger regulation by AP-1. Thus, our findings suggest that, despite the acknowledged similarities in effector functions between ILC subsets and corresponding Th subsets, the underlying regulatory machineries still exhibit substantial distinctions. These insights provide a comprehensive understanding of the unique roles played by each cell type during immune responses.

Non-Coding RNAs in Human Cancer and Other Diseases: Overview of the Diagnostic Potential

Non-coding RNAs (ncRNAs) are abundant single-stranded RNA molecules in human cells, involved in various cellular processes ranging from DNA replication and mRNA translation regulation to genome stability defense. MicroRNAs are multifunctional ncRNA molecules of 18-24 nt in length, involved in gene silencing through base-pair complementary binding to target mRNA transcripts. piwi-interacting RNAs are an animal-specific class of small ncRNAs sized 26-31 nt, responsible for the defense of genome stability via the epigenetic and post-transcriptional silencing of transposable elements. Long non-coding RNAs are ncRNA molecules defined as transcripts of more than 200 nucleotides, their function depending on localization, and varying from the regulation of cell differentiation and development to the regulation of telomere-specific heterochromatin modifications. The current review provides recent data on the several forms of small and long non-coding RNA's potential to act as diagnostic, prognostic or therapeutic target for various human diseases.

Beneficial effects of ginkgetin on improving nonalcoholic steatohepatitis characterized by bulk and single-cell RNA sequencing analysis

Background and aims: Nonalcoholic steatohepatitis (NASH) has become one of the major causes of cirrhosis and liver failure. However, there are currently no approved medications for managing NASH. Our study was designed to assess the effects of ginkgetin on NASH and the involved mechanisms. Methods: We constructed a mouse model of NASH by high-fat diet for 24 weeks. The effects of ginkgetin on NASH were evaluated by histological study, Western blot, and biochemical analysis. RNA Sequencing (RNA-Seq) analysis was used to investigate the alteration in gene expression and signaling pathways at bulk and single-cell levels. Results: Administration of ginkgetin resulted in a marked improvement in hepatic lipid accumulation, inflammation, and fibrosis in the NASH model. And these results were supported by bulk RNA-Seq analysis, in which the related signaling pathways and gene expression were markedly downregulated. Furthermore, single-cell RNA-Seq (scRNA-Seq) analysis revealed that the effects of ginkgetin on NASH were associated with the reprogramming of macrophages, hepatic stellate cells, and endothelial cells. Especially, ginkgetin induced a marked decrease in macrophages and a shift from pro-inflammatory to anti-inflammatory phenotype in NASH mice. And the NASH-associated macrophages (NAMs), which emerge during NASH, were also significantly downregulated by ginkgetin. Conclusion: Ginkgetin exhibits beneficial effects on improving NASH, supported by bulk and single-cell RNA-Seq. Our study may promote pharmacological therapy for NASH and raise the existent understanding of NASH.

Deciphering the heterogeneity and immunosuppressive function of regulatory T cells in osteosarcoma using single-cell RNA transcriptome

Osteosarcoma (OS) is a highly invasive malignant neoplasm with poor prognosis. The tumor microenvironment (TME) plays an essential role in the occurrence and development of OS. Regulatory T cells (Tregs) are known to facilitate immunosuppression, tumor progression, invasion, and metastasis. However, the effect of Tregs in the TME of OS remains unclear. In this study, single-cell RNA sequencing (scRNA-seq) data was used to identify Tregs and various other cell clusters in the TME of OS. Gene set variation analysis (GSVA) was used to investigate the signaling pathways in Tregs from OS and adjacent tissues. The CellChat and iTALK packages were used to analyze cellular communication. In addition, a prognostic model was established based on the Tregs-specific genes using bulk RNA-seq from the TARGET database, and it was verified using a Gene Expression Omnibus dataset. The pRRophetic package was used to predict drug sensitivity. Immunohistochemistry was used to verify the expression of candidate genes in OS. Based on the above methods, we showed that the OS samples were highly infiltrated with Tregs. GSVA revealed that oxidative phosphorylation, angiogenesis and mammalian target of rapamycin complex 1 (mTORC1) were highly activated in Tregs from OS compared with those from adjacent tissues. Using cellular communication analysis, we found that Tregs interacted with osteoblastic, endothelial, and myeloid cells via C-X-C motif chemokine ligand (CXCL) signaling; particularly, they strongly affected the expression of C-X-C motif chemokine receptor 4 (CXCR4) and interacted with other cell clusters through CXCL12/transforming growth factor β1 (TGFB1) to collectively enable tumor growth and progression. Subsequently, two Tregs-specific genes-CD320 and MAF-were screened through univariate, least absolute shrinkage and selection operator regression (LASSO) and multivariate analysis to construct a prognostic model, which showed excellent prognostic accuracy in two independent cohorts. In addition, drug sensitivity analysis demonstrated that OS patients at high Tregs risk were sensitive to sunitinib, sorafenib, and axitinib. We also used immunohistochemistry to validate that CD320 and MAF were significantly upregulated in OS tissues compared with adjacent tissues. Overall, this study reveals the heterogeneity of Tregs in the OS TME, providing new insights into the invasion and treatment of this cancer.

Divergence of cochlear transcriptomics between reference‑based and reference‑free transcriptome analyses among Rhinolophus ferrumequinum populations

Differences in gene expression within tissues can lead to differences in tissue function. Understanding the transcriptome of a species helps elucidate the molecular mechanisms underlying phenotypic divergence. According to the presence or absence of a reference genome of for a studied species, transcriptome analyses can be divided into reference‑based and reference‑free methods, respectively. Presently, comparisons of complete transcriptome analysis results between those two methods are still rare. In this study, we compared the cochlear transcriptome analysis results of greater horseshoe bats (Rhinolophus ferrumequinum) from three lineages in China with different acoustic phenotypes using reference‑based and reference‑free methods to explore their differences in subsequent analysis. The results gained by reference-based results had lower false-positive rates and were more accurate because differentially expressed genes among the three populations obtained by this method had greater reliability and a higher annotation rate. Some phenotype-related enrichment terms, including those related to inorganic molecules and proton transmembrane channels, were also obtained only by the reference-based method. However, the reference‑based method might have the limitation of incomplete information acquisition. Thus, we believe that a combination of reference‑free and reference‑based methods is ideal for transcriptome analyses. The results of our study provided a reference for the selection of transcriptome analysis methods in the future.

Neoantigens: promising targets for cancer therapy

Recent advances in neoantigen research have accelerated the development and regulatory approval of tumor immunotherapies, including cancer vaccines, adoptive cell therapy and antibody-based therapies, especially for solid tumors. Neoantigens are newly formed antigens generated by tumor cells as a result of various tumor-specific alterations, such as genomic mutation, dysregulated RNA splicing, disordered post-translational modification, and integrated viral open reading frames. Neoantigens are recognized as non-self and trigger an immune response that is not subject to central and peripheral tolerance. The quick identification and prediction of tumor-specific neoantigens have been made possible by the advanced development of next-generation sequencing and bioinformatic technologies. Compared to tumor-associated antigens, the highly immunogenic and tumor-specific neoantigens provide emerging targets for personalized cancer immunotherapies, and serve as prospective predictors for tumor survival prognosis and immune checkpoint blockade responses. The development of cancer therapies will be aided by understanding the mechanism underlying neoantigen-induced anti-tumor immune response and by streamlining the process of neoantigen-based immunotherapies. This review provides an overview on the identification and characterization of neoantigens and outlines the clinical applications of prospective immunotherapeutic strategies based on neoantigens. We also explore their current status, inherent challenges, and clinical translation potential.

Anti-Inflammatory Effects of Polyphenols from Plum (Prunus salicina Lindl) on RAW264.7 Macrophages Induced by Monosodium Urate and Potential Mechanisms

Acute gouty arthritis is an acute inflammatory reaction caused by the deposition of monosodium urate (MSU) crystals in joints and surrounding soft tissues. Controlling inflammation is the key to preventing acute gouty arthritis. Anti-inflammatory activities and the possible molecular mechanisms of plum (Prunus salicina Lindl cv. "furong") polyphenols (PSLP) on RAW264.7 macrophage cells induced by monosodium urate were investigated. PPSF significantly inhibited the activity of inflammatory factors such as tumor necrosis factor-α (TNF-α), interleukin-1β (IL-1β), and interleukin-18 (IL-18). In addition, PPSF exhibited excellent activation of superoxide dismutase (SOD) activity and reduction of intracellular reactive oxygen species (ROS) and malondialdehyde (MDA) levels in RAW264.7 macrophages. The results of global screening of all transcripts by RNA-seq revealed 8585 differentially expressed genes between the PSLP-treated group and the MUS group. From GO analysis, PSLP could affect the occurrence and development of RAW264.7 macrophage inflammation through biological processes, such as organic substance metabolism, intracellular organelles, and binding function. The regulation mechanism of PSLP on MSU-induced RAW264.7 macrophage inflammation may be achieved through the HIF-1 signaling pathway, renal cell carcinoma, the ErbB signaling pathway, and the FoxO signaling pathway. Therefore, PSLP has great prospects in the prevention of gout and similar inflammatory diseases.

Thyroid hormone action controls multiple components of cell junctions at the ventricular zone in the newborn rat brain

Thyroid hormone (TH) action controls brain development in a spatiotemporal manner. Previously, we demonstrated that perinatal hypothyroidism led to formation of a periventricular heterotopia in developing rats. This heterotopia occurs in the posterior telencephalon, and its formation was preceded by loss of radial glia cell polarity. As radial glia mediate cell migration and originate in a progenitor cell niche called the ventricular zone (VZ), we hypothesized that TH action may control cell signaling in this region. Here we addressed this hypothesis by employing laser capture microdissection and RNA-Seq to evaluate the VZ during a known period of TH sensitivity. Pregnant rats were exposed to a low dose of propylthiouracil (PTU, 0.0003%) through the drinking water during pregnancy and lactation. Dam and pup THs were quantified postnatally and RNA-Seq of the VZ performed in neonates. The PTU exposure resulted in a modest increase in maternal thyroid stimulating hormone and reduced thyroxine (T4). Exposed neonates exhibited hypothyroidism and T4 and triiodothyronine (T3) were also reduced in the telencephalon. RNA-Seq identified 358 differentially expressed genes in microdissected VZ cells of hypothyroid neonates as compared to controls (q-values ≤0.05). Pathway analyses showed processes like maintenance of the extracellular matrix and cytoskeleton, cell adhesion, and cell migration were significantly affected by hypothyroidism. Immunofluorescence also demonstrated that collagen IV, F-actin, radial glia, and adhesion proteins were reduced in the VZ. Immunohistochemistry of integrin αvβ3 and isoforms of both thyroid receptors (TRα/TRβ) showed highly overlapping expression patterns, including enrichment in the VZ. Taken together, our results show that TH action targets multiple components of cell junctions in the VZ, and this may be mediated by both genomic and nongenomic mechanisms. Surprisingly, this work also suggests that the blood-brain and blood-cerebrospinal fluid barriers may also be affected in hypothyroid newborns.

A Framework for Comparison and Assessment of Synthetic RNA-Seq Data

The ever-growing number of methods for the generation of synthetic bulk and single cell RNA-seq data have multiple and diverse applications. They are often aimed at benchmarking bioinformatics algorithms for purposes such as sample classification, differential expression analysis, correlation and network studies and the optimization of data integration and normalization techniques. Here, we propose a general framework to compare synthetically generated RNA-seq data and select a data-generating tool that is suitable for a set of specific study goals. As there are multiple methods for synthetic RNA-seq data generation, researchers can use the proposed framework to make an informed choice of an RNA-seq data simulation algorithm and software that are best suited for their specific scientific questions of interest.

CLE peptides join the plant longevity club

Leaf senescence, the final step of leaf development, is an essential adaptive process that involves intricate regulatory networks mediated by various developmental and environmental clues. Two recent reports, by Zhang Z. et al. and Zhang Y. et al., shed light on how CLE peptides recruit reactive oxygen species (ROS) and ethylene signaling to promote plant leaf longevity.

Transcriptomic Harmonization as the Way for Suppressing Cross-Platform Bias and Batch Effect

(1) Background: Emergence of methods interrogating gene expression at high throughput gave birth to quantitative transcriptomics, but also posed a question of inter-comparison of expression profiles obtained using different equipment and protocols and/or in different series of experiments. Addressing this issue is challenging, because all of the above variables can dramatically influence gene expression signals and, therefore, cause a plethora of peculiar features in the transcriptomic profiles. Millions of transcriptomic profiles were obtained and deposited in public databases of which the usefulness is however strongly limited due to the inter-comparison issues; (2) Methods: Dozens of methods and software packages that can be generally classified as either flexible or predefined format harmonizers have been proposed, but none has become to the date the gold standard for unification of this type of Big Data; (3) Results: However, recent developments evidence that platform/protocol/batch bias can be efficiently reduced not only for the comparisons of limited transcriptomic datasets. Instead, instruments were proposed for transforming gene expression profiles into the universal, uniformly shaped format that can support multiple inter-comparisons for reasonable calculation costs. This forms a basement for universal indexing of all or most of all types of RNA sequencing and microarray hybridization profiles; (4) Conclusions: In this paper, we attempted to overview the landscape of modern approaches and methods in transcriptomic harmonization and focused on the practical aspects of their application.