Curare and GenExVis: a versatile toolkit for analyzing and visualizing RNA-Seq data.

Even though high-throughput transcriptome sequencing is routinely performed in many laboratories, computational analysis of such data remains a cumbersome process often executed manually, hence error-prone and lacking reproducibility. For corresponding data processing, we introduce Curare, an easy-to-use yet versatile workflow builder for analyzing high-throughput RNA-Seq data focusing on differential gene expression experiments. Data analysis with Curare is customizable and subdivided into preprocessing, quality control, mapping, and downstream analysis stages, providing multiple options for each step while ensuring the reproducibility of the workflow. For a fast and straightforward exploration and visualization of differential gene expression results, we provide the gene expression visualizer software GenExVis. GenExVis can create various charts and tables from simple gene expression tables and DESeq2 results without the requirement to upload data or install software packages. In combination, Curare and GenExVis provide a comprehensive software environment that supports the entire data analysis process, from the initial handling of raw RNA-Seq data to the final DGE analyses and result visualizations, thereby significantly easing data processing and subsequent interpretation.

RNA Sequencing Reveals the Inhibitory Effect of High Levels of Arachidonic Acid and Linoleic Acid on C2C12 Differentiation and Myogenic Biomarkers.

Over the past three decades, studies have shown that consuming polyunsaturated fatty acids (PUFAs) can enhance animal and human health and welfare through biological, biochemical, pathological, and pharmacological impacts. Furthermore, omega-6 plays key roles in the cardiopulmonary system, including promoting airway relaxation and inhibiting atherosclerosis and hypertension. However, findings from investigations of the effects of omega-6 fatty acids on molecular and cellular activity and discussions on their influence on biomarkers are still unclear. Therefore, the present study aimed to evaluate omega-6 fatty acids, the arachidonic acid (AA), and linoleic acid (LA) effects on C2C12 proliferation, myogenesis morphology, and relative myogenic biomarker expression through the Wnt pathway. C2C12 cells were cultured with and without 25, 50, 100, and 150 µM of LA and AA and then subjected to CCK8, Giemsa staining, RT qPCR, Western blotting, and RNA Sequencing. The CCK8 Assay results showed that 25, 50, 100, and 150 µM LA significantly decreased the viability after 72 h for 25, 50, 100, and 150 µM concentrations. Also, AA supplementation decreased cell viability after 24 h for 150 µM, 48 h for 150 µM, and 72 h for 50, 100, and 150 µM concentrations. Moreover, the LA and AA inhibitory effects noticed through Gimesa staining were morphological changes during myoblast differentiation. Both LA and AA showed inhibiting IGF1, Cola1, Col6a2, Col6a1, Itga10, Itga11, SFRP2, DAAM2, and NKD2 effects; however, the depressing effect was higher for AA compared to LA. The previous results were confirmed through Western blotting, which showed that 50 µM LA and AA significantly reduced DAAM2 and SFRP2 protein levels compared to the control. Regarding RNA sequencing results, LA and AA increased the number of differentially expressed (DE) Mt-rRNA and snoRNA; however, the numbers of lncRNA detected decreased compared to the control. Our findings demonstrate that high and moderate LA and AA concentrations reduce primary myoblast proliferation and differentiation. Also, they highlight novel biomarkers and regulatory factors to improve our understanding of how the nutrition of fatty acids can control and modulate the myogenesis and differentiation process through different biomarker families.

Integration of single-nucleus RNA sequencing and network disturbance to elucidate crosstalk between multicomponent drugs and trigeminal ganglia cells in migraine.

ETHNOPHARMACOLOGICAL RELEVANCE: Migraine is caused by hyperactivity of the trigeminovascular system, where trigeminal ganglia (TG) plays an important role. TG is composed of multiple neuronal and non-neuronal cell types, which is related to "neuro-inflammation-vascular" disorder in migraine. Tou Tong Ning capsule (TTNC), a CFDA-approved traditional Chinese medicine for treating migraine, has the characteristics of "multicomponents, multitargets, multipathways". AIM OF THE STUDY: To clarify the mechanism of TTNC and elucidate crosstalk between multicomponent drugs and neuronal and non-neuronal functions and cells in migraine. MATERIALS AND METHODS: We integrated single-nucleus RNA sequencing and a quantitative evaluation algorithm of the disturbance of multitarget drugs on the disease network and explored the specific pathology of migraine and corresponding compounds. A cerebrovascular smooth muscle spasmolytic activity experiment was carried out to verify the results of the bioinformatics analysis. RESULTS: TTNC exhibited its regulation activities in neuronal and non-neuronal aspects based on drugs attack to four subnetworks and cell specific networks, which explored the MoA of TTNC in comprehensive and refined perspectives. Compared to neuronal regulation, TTNC showed more significant attack score on non-neuronal biological function (smooth muscle and vessel). And TTNC compound clusters C1, C6 and C7, targeting non-neuronal function and cells, had larger group area than C10, C4 and C6 for neuronal function and cell, which implied that TTNC may mainly regulate the non-neuronal function, e.g., vessel smooth muscle contraction. Contraction of cerebrovascular smooth muscle of mice ex vivo confirmed the vasodilation activity of TTNC and active compounds from C1, C6, C9 (Emodin, Luteolin and Levistilide A). Literature mining confirmed the vasospasmodolytic activity and neuroprotective effect of TTNC. CONCLUSIONS: The study found that TTNC may primarily alleviate non-neuronal functional disorders in migraine by relaxing cerebral vascular smooth muscle cell spasm to alleviate migraine. Integrating single-nucleus RNA sequencing data and network disturbance tools provides a new strategy for the pharmacological mechanism of multicomponent drugs through cell subtyping.

Novel insights into the mechanism(s) of silicon-induced drought stress tolerance in lentil plants revealed by RNA sequencing analysis.

BACKGROUND: Lentil is an essential cool-season food legume that offers several benefits in human nutrition and cropping systems. Drought stress is the major environmental constraint affecting lentil plants' growth and productivity by altering various morphological, physiological, and biochemical traits. Our previous research provided physiological and biochemical evidence showing the role of silicon (Si) in alleviating drought stress in lentil plants, while the molecular mechanisms are still unidentified. Understanding the molecular mechanisms of Si-mediated drought stress tolerance can provide fundamental information to enhance our knowledge of essential gene functions and pathways modulated by Si during drought stress in plants. Thus, the present study compared the transcriptomic characteristics of two lentil genotypes (drought tolerant-ILL6002; drought sensitive-ILL7537) under drought stress and investigated the gene expression in response to Si supplementation using high-throughput RNA sequencing. RESULTS: This study identified 7164 and 5576 differentially expressed genes (DEGs) from drought-stressed lentil genotypes (ILL 6002 and ILL 7537, respectively), with Si treatment. RNA sequencing results showed that Si supplementation could alter the expression of genes related to photosynthesis, osmoprotection, antioxidant systems and signal transduction in both genotypes under drought stress. Furthermore, these DEGs from both genotypes were found to be associated with the metabolism of carbohydrates, lipids and proteins. The identified DEGs were also linked to cell wall biosynthesis and vasculature development. Results suggested that Si modulated the dynamics of biosynthesis of alkaloids and flavonoids and their metabolism in drought-stressed lentil genotypes. Drought-recovery-related DEGs identified from both genotypes validated the role of Si as a drought stress alleviator. This study identified different possible defense-related responses mediated by Si in response to drought stress in lentil plants including cellular redox homeostasis by reactive oxygen species (ROS), cell wall reinforcement by the deposition of cellulose, lignin, xyloglucan, chitin and xylan, secondary metabolites production, osmotic adjustment and stomatal closure. CONCLUSION: Overall, the results suggested that a coordinated interplay between various metabolic pathways is required for Si to induce drought tolerance. This study identified potential genes and different defence mechanisms involved in Si-induced drought stress tolerance in lentil plants. Si supplementation altered various metabolic functions like photosynthesis, antioxidant defence system, osmotic balance, hormonal biosynthesis, signalling, amino acid biosynthesis and metabolism of carbohydrates and lipids under drought stress. These novel findings validated the role of Si in drought stress mitigation and have also provided an opportunity to enhance our understanding at the genomic level of Si's role in alleviating drought stress in plants.

Sanguisorba officinalis ethyl acetate extract attenuates ulcerative colitis through inhibiting PI3K-AKT/NF-κB/ STAT3 pathway uncovered by single-cell RNA sequencing.

BACKGROUND: Ulcerative colitis (UC) accounts for the untreatable illness nowadays. Bloody stools are the primary symptom of UC, and the first-line drugs used to treat UC are associated with several drawbacks and negative side effects. S. officinalis has long been used as a medicine to treat intestinal infections and bloody stools. However, what the precise molecular mechanism, the exact etiology, and the material basis of the disease remain unclear. PURPOSE: This work aimed to comprehensively explore pharmacological effects as well as molecular mechanisms underlying the active fraction of S. officinalis, and to produce a comprehensive and brand-new guideline map of its chemical base and mechanism of action. METHODS: First, different polarity S. officinalis extracts were orally administered to the DSS-induced UC model mice for the sake of investigating its active constituents. Using the UPLC-orbitrap high-resolution mass spectrometry (UPLC-Q-Orbitrap-HRMS) technique, the most active S. officinalis (S. officinalis ethyl acetate fraction, SOEA) extract was characterized. Subsequently, the effectiveness of its active fraction on UC was evaluated through phenotypic observation (such as weight loss, colon length, and stool characteristics), and histological examination of pathological injuries, mRNA and protein expression. Cell profile, cell-cell interactions and molecular mechanisms of SOEA in different cell types of the colon tissue from UC mice were described using single-cell RNA sequencing (scRNA-seq). As a final step, the molecular mechanisms were validated by appropriate molecular biological methods. RESULTS: For the first time, this study revealed the significant efficacy of SOEA in the treatment of UC. SOEA reduced DAI and body weight loss, recovered the colon length, and mitigated colonic pathological injuries along with mucosal barrier by promoting goblet cell proliferation. Following treatment with SOEA, inflammatory factors showed decreased mRNA and protein expression. SOEA restored the dynamic equilibrium of cell profile and cell-cell interactions in colon tissue. All of these results were attributed to the ability of SOEA to inhibit the PI3K-AKT/NF-κB/STATAT pathway. CONCLUSIONS: By integrating the chemical information of SOEA derived from UPLC-Q-Orbitrap-HRMS with single-cell transcriptomic data extracted from scRNA-seq, this study demonstrates that SOEA exerts the therapeutic effect through suppressing PI3K-AKT/NF-B/STAT3 pathway to improve clinical symptoms, inflammatory response, mucosal barrier, and intercellular interactions in UC, and effectively eliminates the interference of cellular heterogeneity.

Comparative analysis of PacBio and ONT RNA sequencing methods for Nemopilema Nomurai venom identification.

Recent studies on marine organisms have made use of third-generation sequencing technologies such as Pacific Biosciences (PacBio) and Oxford Nanopore Technologies (ONT). While these specialized bioinformatics tools have different algorithmic designs and performance capabilities, they offer scalability and can be applied to various datasets. We investigated the effectiveness of PacBio and ONT RNA sequencing methods in identifying the venom of the jellyfish species Nemopilema nomurai. We conducted a detailed analysis of the sequencing data from both methods, focusing on key characteristics such as CD, alternative splicing, long-chain noncoding RNA, simple sequence repeat, transcription factor, and functional transcript annotation. Our findings indicate that ONT generally produced higher raw data quality in the transcriptome analysis, while PacBio generated longer read lengths. PacBio was found to be superior in identifying CDs and long-chain noncoding RNA, whereas ONT was more cost-effective for predicting alternative splicing events, simple sequence repeats, and transcription factors. Based on these results, we conclude that PacBio is the most specific and sensitive method for identifying venom components, while ONT is the most cost-effective method for studying venogenesis, cnidocyst (venom gland) development, and transcription of virulence genes in jellyfish. Our study has implications for future sequencing technologies in marine jellyfish, and highlights the power of full-length transcriptome analysis in discovering potential therapeutic targets for jellyfish dermatitis.

Extracellular matrix remodelling in dental pulp tissue of carious human teeth through the prism of single-cell RNA sequencing.

Carious lesions are bacteria-caused destructions of the mineralised dental tissues, marked by the simultaneous activation of immune responses and regenerative events within the soft dental pulp tissue. While major molecular players in tooth decay have been uncovered during the past years, a detailed map of the molecular and cellular landscape of the diseased pulp is still missing. In this study we used single-cell RNA sequencing analysis, supplemented with immunostaining, to generate a comprehensive single-cell atlas of the pulp of carious human teeth. Our data demonstrated modifications in the various cell clusters within the pulp of carious teeth, such as immune cells, mesenchymal stem cells (MSC) and fibroblasts, when compared to the pulp of healthy human teeth. Active immune response in the carious pulp tissue is accompanied by specific changes in the fibroblast and MSC clusters. These changes include the upregulation of genes encoding extracellular matrix (ECM) components, including COL1A1 and Fibronectin (FN1), and the enrichment of the fibroblast cluster with myofibroblasts. The incremental changes in the ECM composition of carious pulp tissues were further confirmed by immunostaining analyses. Assessment of the Fibronectin fibres under mechanical strain conditions showed a significant tension reduction in carious pulp tissues, compared to the healthy ones. The present data demonstrate molecular, cellular and biomechanical alterations in the pulp of human carious teeth, indicative of extensive ECM remodelling, reminiscent of fibrosis observed in other organs. This comprehensive atlas of carious human teeth can facilitate future studies of dental pathologies and enable comparative analyses across diseased organs.

The efficacy and mechanism of Angelica sinensis (Oliv.) Diels root aqueous extract based on RNA sequencing and 16S rDNA sequencing in alleviating polycystic ovary syndrome.

BACKGROUND: Polycystic ovary syndrome (PCOS) leads to persistent anovulation, hyperandrogenism, insulin resistance, and polycystic ovary, and is mainly characterized by menstrual disorders, and reproductive dysfunction. Angelica sinensis (Oliv.) Diels root has been used in many classical formulas of traditional Chinese medicine, and is commonly used to treat various gynecological diseases. PURPOSE: To investigate the protective effect of water extract of A. sinensis root (WEA) on PCOS rats, and the mechanism by RNA sequencing, and 16S rDNA sequencing. METHODS: The PCOS rat model was established by letrozole combined with high-fat diet (gavage; 2 months), and treated with WEA (gavage; 2 g/kg, 4 g/kg or 8 g/kg; 1 month). To evaluate the therapeutic effect of WEA on PCOS rats, vaginal smear, hematoxylin-eosin staining, and biochemical indicators detection were performed. The rat ovarian tissue was analyzed by RNA sequencing, and the results were verified by qRT-PCR, and Western blot. 16S rDNA sequencing was used to analyze the gut microbiota of rats. RESULTS: The results of the vaginal smear, and hematoxylin-eosin staining showed that WEA improved estrous cycle disorder, and ovarian tissue lesions. WEA (4 g/kg or 8 g/kg; 1 months) alleviated hormone disorders, insulin resistance, and dyslipidemia. RNA sequencing showed that WEA intervention significantly changed the expressions of 2756 genes, which were enriched in phosphatidylinositol3-kinase/phosphorylated protein kinase B (PI3K/AKT), peroxisome proliferator-activated receptor (PPAR), mitogen-activated protein kinase (MAPK), AMP-activated protein kinase (AMPK), and insulin signaling pathways. 16S rDNA sequencing found that WEA increased the species diversity of gut microbiota, and regulated the abundance of some microbiota (genus level: Dubosiella, Bifidobacterium, Coriobacteriaceae (UCG-002), and Treponema; species level: Bifidobacterium animalis, Lactobacillus murinus, and Lactobacillus johnsonii). CONCLUSION: WEA regulated hormone, and glycolipid metabolism disorders, thereby relieving the PCOS induced by letrozole combined with high-fat diet. The mechanism was related to the regulation of PI3K/AKT, PPAR, MAPK, AMPK, and insulin signaling pathways in ovarian tissues, and the maintenance of gut microbiota homeostasis. Clarifying the efficacy and mechanism of WEA in alleviating PCOS based on RNA sequencing and 16S rDNA sequencing will guide the more reasonable clinical use of WEA.

Selenium metabolism heterogeneity in pan-cancer: insights from bulk and single-cell RNA sequencing.

BACKGROUND: Selenium, a natural microelement with both nutritional and toxicological properties, is intertwined with tumorigenesis and progression. However, it is not fully understood how selenium metabolism affects immune response and cancer biology. METHODS: We estimated selenium metabolism by Gene Set Enrichment Analysis (GSEA) to delineate the selenium metabolism landscape using The Cancer Genome Atlas (TCGA), Genotype-Tissue Expression (GTEx), Cancer Cell Line Encyclopedia (CCLE) and a integrated pan-cancer single-cell dataset. We systematically explored the prognostic implications of selenium metabolism and selenium-related regulatory patterns. The therapeutic value of selenium metabolism was explored through machine learning and examined in several immunotherapy cohorts. The heterogeneity and underlying mechanism of selenium metabolism were investigated by cell‒cell communication analysis at the single-cell level. RESULTS: A GSEA analysis based on 86 genes was used to evaluate the selenium metabolism landscape. The selenium metabolism score exhibited prognostic value in predicting the lower risk of mortality, possibly due to its correlation with multiple cancer hallmarks, including a positive correlation with complement (R = 0.761, P < 0.001), inflammatory response (R = 0.663, P < 0.001), apoptosis (R = 0.626, P < 0.001), hypoxia (R = 0.587, P < 0.001), reactive oxygen species (ROS) (R = 0.558, P < 0.001), and interferon gamma response (R = 0.539, P < 0.001). We also observed heterogeneity in the relationship between selenium metabolism and immunity across different cancers. Based on selenium-related genes, we constructed a machine learning model with area under the ROC curve (AUC) of 0.82 in predicting immune checkpoint inhibitor (ICI)-based immunotherapy response. Single-cell selenium metabolism quantification revealed that adjacent and tumor tissues had higher selenium metabolism compared with normal tissues, especially in epithelial cells, fibroblasts and macrophages. The communication between high-selenium epithelium and high-selenium fibroblast was significantly higher than other cells, especially in cytokines, chemokines, collagen, Wnt, VEGF, IGF and FGF pathways. CONCLUSION: Our study provides a comprehensive landscape of selenium metabolism levels and diverse regulatory patterns in different cancers, deepening the understanding of selenium's roles in tumorigenesis and immunity.

Single-cell RNA sequencing reveals enhanced antitumor immunity after combined application of PD-1 inhibitor and Shenmai injection in non-small cell lung cancer.

BACKGROUND: Immune checkpoint inhibitors (ICIs) have altered the clinical management of non-small cell lung cancer (NSCLC). However, the low response rate, severe immune-related adverse events (irAEs), and hyperprogressive disease following ICIs monotherapy require attention. Combination therapy may overcome these limitations and traditional Chinese medicine with immunomodulatory effects provides a promising approach. Shenmai injection (SMI) is a clinically effective adjuvant treatment for cancer with chemotherapy and radiotherapy. Therefore, the combined effects and mechanisms of SMI and programmed death-1 (PD-1) inhibitor against NSCLC was focused on this study. METHODS: A Lewis lung carcinoma mouse model and a lung squamous cell carcinoma humanized mouse model were used to investigate the combined efficacy and safety of SMI and PD-1 inhibitor. The synergistic mechanisms of the combination therapy against NSCLC were explored using single-cell RNA sequencing. Validation experiments were performed using immunofluorescence analysis, in vitro experiment, and bulk transcriptomic datasets. RESULTS: In both models, combination therapy alleviated tumor growth and prolonged survival without increasing irAEs. The GZMAhigh and XCL1high natural killer (NK) cell subclusters with cytotoxic and chemokine signatures increased in the combination therapy, while malignant cells from combination therapy were mainly in the apoptotic state, suggesting that mediating tumor cell apoptosis through NK cells is the main synergistic mechanisms of combination therapy. In vitro experiment confirmed that combination therapy increased secretion of Granzyme A by NK cells. Moreover, we discovered that PD-1 inhibitor and SMI combination blocked inhibitory receptors on NK and T cells and restores their antitumoral activity in NSCLC better than PD-1 inhibitor monotherapy, and immune and stromal cells exhibited a decrease of angiogenic features and attenuated cancer metabolism reprogramming in microenvironment of combination therapy. CONCLUSIONS: This study demonstrated that SMI reprograms tumor immune microenvironment mainly by inducing NK cells infiltration and synergizes with PD-1 inhibitor against NSCLC, suggested that targeting NK cells may be an important strategy for combining with ICIs. Video Abstract.

Exploring the effect and mechanism of Haizao Yuhu decoction containing three variants of glycyrrhiza on goiter using an integrated strategy of network pharmacology and RNA sequencing.

ETHNOPHARMACOLOGICAL RELEVANCE: Haizao Yuhu decoction (HYD) is a classic Chinese herbal formula described in the surgical monographs of the Ming Dynasty "Waikezhengzong." It has been widely used to treat goiter for approximately 500 years and found to be particularly effective. HYD contains glycyrrhiza and sargassum. This pair of herbs belongs to "18 incompatible medicaments" of traditional Chinese medicine theory. Although these two herbs are opposite, our preliminary study proved that they have superior effect when added into HYD at 2 times the dose of Chinese Pharmacopoeia. However, the species of glycyrrhiza in HYD that are the most effective have not been recorded in ancient Chinese medical texts. According to the Chinese Pharmacopoeia, glycyrrhiza is divided into the following three species: Glycyrrhiza uralensis Fish., G. glabra L., and G. inflata Bat. The effect of HYD containing different species of glycyrrhiza and their mechanisms remain to be further explored. AIM OF THE STUDY: To investigate the effect of HYD containing three species of glycyrrhiza on goiter, and to elucidate the molecular mechanism using network pharmacology combined with RNA sequencing (RNA-seq). MATERIALS AND METHODS: A rat model of goiter was established by 14 days of intragastric gavage of propylthiouracil (PTU), and the rats were treated for 4 weeks with HYD containing three different species of glycyrrhiza. The body weight and rectal temperature of rats were tested weekly. At the end of the experiment, the serum and thyroid tissues of rats were collected. The effect of the three HYDs was assessed based on general observations (including body weight, rectal temperature, and living status of rats), absolute/relative thyroid weight, thyroid function (including triiodothyronine, thyroxine, free triiodothyronine, free thyroxine, and thyroid-stimulating hormone levels), and thyroid tissue pathology. Next, we explored their pharmacological mechanisms using network pharmacology combined with RNA-seq and validated key targets using real-time quantitative reverse-transcription polymerase chain reaction (RT-qPCR), western blotting (WB), and immunofluorescence (IF) assays. RESULTS: The three HYDs reduced the absolute/relative weights of thyroid tissues and improved the pathological structure, thyroid function, and general findings of rats with goiter. Overall, the effect of HYD-G. uralensis Fish. (HYD-U) was better. Results from network pharmacology and RNA-seq jointly suggested that both the pathogenesis of goiter and the mechanism of action of HYD for goiter were related to the phosphatidylinositol 3-kinase-protein kinase B (PI3K-Akt) pathway. We validated the key targets in the pathway, namely, vascular endothelial growth factor (VEGF) A, VEGF receptor 2, phosphoinositide-3-kinase regulatory subunit 1 (PIK3R1) and its encoded protein PI3K (p85), AKT serine/threonine kinase 1 (AKT1), phospho-AKT and cyclin D1 using RT-qPCR, WB, and IF assays. The PI3K-Akt pathway was hyperactivated in rats with PTU-induced goiter, whereas the three HYDs could inhibit the pathway. CONCLUSION: This study confirmed the definite effect of the three HYDs in the treatment of goiter, and HYD-U was found to be more effective. The three HYDs inhibited angiogenesis and cell proliferation in goiter tissue by inhibiting the PI3K-Akt signaling pathway.

MicroRNAs from Holarrhena pubescens stems: Identification by small RNA Sequencing and their Potential Contribution to Human Gene Targets.

Holarrhena pubescens is an effective medicinal plant from the Apocynaceae family, widely distributed over the Indian subcontinent and extensively used by Ayurveda and ethno-medicine systems without apparent side effects. We postulated that miRNAs, endogenous non-coding small RNAs that regulate gene expression at the post-transcriptional level, may, after ingestion into the human body, contribute to the medicinal properties of plants of this species by inducing regulated human gene expression to modulate. However, knowledge is scarce about miRNA in Holarrhena. In addition, to test the hypothesis on the potential pharmacological properties of miRNA, we performed a high-throughput sequencing analysis using the Next Generation Sequencing Illumina platform; 42,755,236 raw reads have been generated from H. pubescens stems from a library of small RNA isolated, identifying 687 known and 50 new miRNAs led. The novel H. pubescens miRNAs were predicted to regulate specific human genes, and subsequent annotations of gene functions suggested a possible role in various biological processes and signaling pathways, such as Wnt, MAPK, PI3K-Akt, and AMPK signaling pathways and endocytosis. The association of these putative targets with many diseases, including cancer, congenital malformations, nervous system disorders, and cystic fibrosis, has been demonstrated. The top hub proteins STAT3, MDM2, GSK3B, NANOG, IGF1, PRKCA, SNAP25, SRSF1, HTT, and SNCA show their interaction with human diseases, including cancer and cystic fibrosis. To our knowledge, this is the first report of uncovering H. pubescens miRNAs based on high-throughput sequencing and bioinformatics analysis. This study has provided new insight into a potential cross-species control of human gene expression. The potential for miRNA transfer should be evaluated as one possible mechanism of action to account for the beneficial properties of this valuable species.

A comparative analysis of small RNA sequencing data in tubers of purple potato and its red mutant reveals small RNA regulation in anthocyanin biosynthesis.

Anthocyanins are a group of natural pigments acting as stress protectants induced by biotic/abiotic stress in plants. Although the metabolic pathway of anthocyanin has been studied in potato, the roles of miRNAs on the metabolic pathway remain unclear. In this study, a purple tetraploid potato of SD92 and its red mutant of SD140 were selected to explore the regulation mechanism of miRNA in anthocyanin biosynthesis. A comparative analysis of small RNAs between SD92 and SD140 revealed that there were 179 differentially expressed miRNAs, including 65 up- and 114 down-regulated miRNAs. Furthermore, 31 differentially expressed miRNAs were predicted to potentially regulate 305 target genes. KEGG pathway enrichment analysis for these target genes showed that plant hormone signal transduction pathway and plant-pathogen interaction pathway were significantly enriched. The correlation analysis of miRNA sequencing data and transcriptome data showed that there were 140 negative regulatory miRNA-mRNA pairs. The miRNAs included miR171 family, miR172 family, miR530b_4 and novel_mir170. The mRNAs encoded transcription factors, hormone response factors and protein kinases. All these results indicated that miRNAs might regulate anthocyanin biosynthesis through transcription factors, hormone response factors and protein kinase.

Detection of multiple types of cancer driver mutations using targeted RNA sequencing in non-small cell lung cancer.

BACKGROUND: DNA-based next-generation sequencing has been widely used in the selection of target therapies for patients with nonsmall cell lung cancer (NSCLC). RNA-based next-generation sequencing has been proven to be valuable in detecting fusion and exon-skipping mutations and is recommended by National Comprehensive Cancer Network guidelines for these mutation types. METHODS: The authors developed an RNA-based hybridization panel targeting actionable driver oncogenes in solid tumors. Experimental and bioinformatics pipelines were optimized for the detection of fusions, single-nucleotide variants (SNVs), and insertion/deletion (indels). In total, 1253 formalin-fixed, paraffin-embedded samples from patients with NSCLC were analyzed by DNA and RNA panel sequencing in parallel to assess the performance of the RNA panel in detecting multiple types of mutations. RESULTS: In analytical validation, the RNA panel achieved a limit of detection of 1.45-3.15 copies per nanogram for SNVs and 0.21-6.48 copies per nanogram for fusions. In 1253 formalin-fixed, paraffin-embedded NSCLC samples, the RNA panel identified a total of 124 fusion events and 26 MET exon 14-skipping events, in which 14 fusions and six MET exon 14-skipping mutations were missed by DNA panel sequencing. By using the DNA panel as the reference, the positive percent agreement and the positive predictive value of the RNA panel were 98.08% and 98.62%, respectively, for detecting targetable SNVs and 98.15% and 99.38%, respectively, for detecting targetable indels. CONCLUSIONS: Parallel DNA and RNA sequencing analyses demonstrated the accuracy and robustness of the RNA sequencing panel in detecting multiple types of clinically actionable mutations. The simplified experimental workflow and low sample consumption will make RNA panel sequencing a potentially effective method in clinical testing.

RNA sequencing in Artemisia annua L explored the genetic and metabolic responses to hardly soluble aluminum phosphate treatment.

Artemisia annua L. is a medicinal plant valued for its ability to produce artemisinin, a molecule used to treat malaria. Plant nutrients, especially phosphorus (P), can potentially influence plant biomass and secondary metabolite production. Our work aimed to explore the genetic and metabolic response of A. annua to hardly soluble aluminum phosphate (AlPO4, AlP), using soluble monopotassium phosphate (KH2PO4, KP) as a control. Liquid chromatography-mass spectrometry (LC-MS) was used to analyze artemisinin. RNA sequencing, gene ontology (GO), and the Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses were applied to analyze the differentially expressed genes (DEGs) under poor P conditions. Results showed a significant reduction in plant growth parameters, such as plant height, stem diameter, number of leaves, leaf areas, and total biomass of A. annua. Conversely, LC-MS analysis revealed a significant increase in artemisinin concentration under the AlP compared to the KP. Transcriptome analysis revealed 762 differentially expressed genes (DEGs) between the AlP and the KP. GH3, SAUR, CRE1, and PYL, all involved in plant hormone signal transduction, showed differential expression. Furthermore, despite the downregulation of HMGR in the artemisinin biosynthesis pathway, the majority of genes (ACAT, FPS, CYP71AV1, and ALDH1) were upregulated, resulting in increased artemisinin accumulation in the AlP. In addition, 12 transcription factors, including GATA and MYB, were upregulated in response to AlP, confirming their importance in regulating artemisinin biosynthesis. Overall, our findings could contribute to a better understanding the parallel transcriptional regulation of plant hormone transduction and artemisinin biosynthesis in A. annua L. in response to hardly soluble phosphorus fertilizer.

Transcriptomic profiling and differential analysis reveal the renal toxicity mechanisms of mice under cantharidin exposure.

Cantharidin (CTD), extracted from the traditional Chinese medicine mylabris, has shown significant curative effects against a variety of tumors, but its clinical application is limited by its high toxicity. Studies have revealed that CTD can cause toxicity in the kidneys; however, the underlying molecular mechanisms remain unclear. In this study, we investigated the toxic effects in mouse kidneys following CTD treatment by pathological and ultrastructure observations, biochemical index detection, and transcriptomics, and explored the underlying molecular mechanisms by RNA sequencing (RNA-seq). The results showed that after CTD exposure, the kidneys had different degrees of pathological damage, altered uric acid and creatinine levels in serum, and the antioxidant indexes in tissues were significantly increased. These changes were more pronounced at medium and high doses of CTD. RNA-seq analysis revealed 674 differentially expressed genes compared with the control group, of which 131 were upregulated and 543 were downregulated. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses showed that many differentially expressed genes were closely related to the stress response, the CIDE protein family, and the transporter superfamily, as well as the MAPK, AMPK, and HIF-1 pathways. The reliability of the RNA-seq results was verified by qRT-PCR of the six target genes. These findings offer insight into the molecular mechanisms of renal toxicity caused by CTD and provide an important theoretical basis for the clinical treatment of CTD-induced nephrotoxicity.

Zinc Supplementation Induced Transcriptional Changes in Primary Human Retinal Pigment Epithelium: A Single-Cell RNA Sequencing Study to Understand Age-Related Macular Degeneration.

Zinc supplementation has been shown to be beneficial to slow the progression of age-related macular degeneration (AMD). However, the molecular mechanism underpinning this benefit is not well understood. This study used single-cell RNA sequencing to identify transcriptomic changes induced by zinc supplementation. Human primary retinal pigment epithelial (RPE) cells could mature for up to 19 weeks. After 1 or 18 weeks in culture, we supplemented the culture medium with 125 µM added zinc for one week. RPE cells developed high transepithelial electrical resistance, extensive, but variable pigmentation, and deposited sub-RPE material similar to the hallmark lesions of AMD. Unsupervised cluster analysis of the combined transcriptome of the cells isolated after 2, 9, and 19 weeks in culture showed considerable heterogeneity. Clustering based on 234 pre-selected RPE-specific genes divided the cells into two distinct clusters, we defined as more and less differentiated cells. The proportion of more differentiated cells increased with time in culture, but appreciable numbers of cells remained less differentiated even at 19 weeks. Pseudotemporal ordering identified 537 genes that could be implicated in the dynamics of RPE cell differentiation (FDR < 0.05). Zinc treatment resulted in the differential expression of 281 of these genes (FDR < 0.05). These genes were associated with several biological pathways with modulation of ID1/ID3 transcriptional regulation. Overall, zinc had a multitude of effects on the RPE transcriptome, including several genes involved in pigmentation, complement regulation, mineralization, and cholesterol metabolism processes associated with AMD.

The Evolution of Single-Cell RNA Sequencing Technology and Application: Progress and Perspectives.

As an emerging sequencing technology, single-cell RNA sequencing (scRNA-Seq) has become a powerful tool for describing cell subpopulation classification and cell heterogeneity by achieving high-throughput and multidimensional analysis of individual cells and circumventing the shortcomings of traditional sequencing for detecting the average transcript level of cell populations. It has been applied to life science and medicine research fields such as tracking dynamic cell differentiation, revealing sensitive effector cells, and key molecular events of diseases. This review focuses on the recent technological innovations in scRNA-Seq, highlighting the latest research results with scRNA-Seq as the core technology in frontier research areas such as embryology, histology, oncology, and immunology. In addition, this review outlines the prospects for its innovative application in traditional Chinese medicine (TCM) research and discusses the key issues currently being addressed by scRNA-Seq and its great potential for exploring disease diagnostic targets and uncovering drug therapeutic targets in combination with multiomics technologies.

Single-cell RNA sequencing reveals unique monocyte-derived interstitial macrophage subsets during lipopolysaccharide-induced acute lung inflammation.

Interstitial macrophages (IMs) reside in the lung tissue surrounding key structures including airways, vessels, and alveoli. Recent work has described IM heterogeneity during homeostasis, however, there are limited data on IMs during inflammation. We sought to characterize IM origin, subsets, and transcriptomic profiles during homeostasis and lipopolysaccharide (LPS) induced acute lung inflammation. During homeostasis, we used three complementary methods, spectral flow cytometry, single-cell RNA-sequencing, and gene regulatory network enrichment, to demonstrate that IMs can be divided into two core subsets distinguished by surface and transcriptional expression of folate receptor ß (Folr2/FRß). These subsets inhabited distinct niches within the lung interstitium. Within FRß+ IMs we identified a subpopulation marked by coexpression of LYVE1. During acute LPS-induced inflammation, lung IM numbers expand. Lineage tracing revealed IM expansion was due to recruitment of monocyte-derived IMs. At the peak of inflammation, recruited IMs were comprised two unique subsets defined by expression of genes associated with interferon signaling and glycolytic pathways. As recruited IMs matured, they adopted the overall transcriptional state of FRß- resident IMs but retained expression in several origin-specific genes, such as IL-1ß. FRß+ IMs were of near-pure resident origin. Taken together our data show that during LPS-induced inflammation, there are distinct populations of IMs that likely have unique functions. FRΒ+ IMs comprise a stable, resident population, whereas FRß- ΙΜs represent a mixed population of resident and recruited IMs.

Small RNA sequencing and identification of Andrographis paniculata miRNAs with potential cross­kingdom human gene targets.

Cross-species post-transcriptional regulatory potential of plant derived small non-coding microRNAs (miRNAs) has been well documented by plenteous studies. MicroRNAs are transferred to host cells via oral ingestion wherein they play a decisive role in regulation of host genes; thus, miRNAs have evolved as the nascent bioactive molecules imparting pharmacological values to traditionally used medicinal plants. The present study aims to investigate small RNA profiling in order to uncover the potential regulatory role of miRNAs derived from Andrographis paniculata, one of the most widely used herb by tribal communities for liver disorders and document the pharmacological properties of A. paniculata miRNAs. In this study, high-throughput sequencing method was used to generate raw data, ~ 60 million sequences were generated from A. paniculata leaves. Using computational tools and bioinformatics approach, analyses of 3,480,097 clean reads resulted in identification of 3440 known and 51 putative novel miRNAs regulating 1365 and 192 human genes respectively. Remarkably, the identified plausible novel miRNAs apa-miR-5, apa-miR-1, apa-miR-26, and apa-miR-30 are projected to target significant host genes including CDK6, IKBKB, TRAF3, CHD4, MECP2, and ADIPOQ. Subsequent annotations revealed probable involvement of the target genes in various pathways for instance p38-MAPK, AKT, AMPK, NF-Kß, ERK, WNT signalling, MYD88 dependant cascade, and pathways in cancer. Various diseases such as human papilloma virus infection, Alzheimer's, Non-alcoholic Fatty Liver, Alcoholic liver diseases, HepatoCellular Carcinoma (HCC), and numerous other cancers were predominantly found to be linked with target genes. Our findings postulate novel interpretations regarding modulation of human transcripts by A. paniculata miRNAs and exhibit the regulation of human diseases by plant-derived miRNAs. Though our study elucidates miRNAs as novel therapeutic agents, however, experimental validations for assessment of therapeutic potential of these miRNAs are still warranted.