High-throughput Transcriptome Sequencing Research Articles

Transcriptome analysis of Vero cells infected with attenuated vaccine strain CDV-QN-1

To better understand the interaction between attenuated vaccines and host antiviral responses, we used bioinformatics and public transcriptomics data to analyze the immune response mechanisms of host cells after canine distemper virus (CDV) infection in Vero cells and screened for potential key effector factors. In this study, CDV-QN-1 infect with Vero cells at an MOI of 0.5, and total RNA was extracted from the cells 24 hours later and reverse transcribed into cDNA. Transcriptome high-throughput sequencing perform using Illumina. The results showed that 438 differentially expressed genes were screened, of which 409 were significantly up-regulated and 29 were significantly down-regulated. Eight differentially expressed genes were randomly selected for RT-qPCR validation, and the change trend was consistent with the transcriptomics data. GO and KEGG analysis of differentially expressed genes revealed that most of the differentially expressed genes in CDV-QN-1 infection in the early stage were related to immune response and antiviral activity. The enriched signaling pathways mainly included the interaction between cytokines and cytokine receptors, the NF-kappa B signaling pathway, the Toll-like receptor signaling pathway, and the NOD-like receptor signaling pathway. This study provides a foundation for further exploring the pathogenesis of CDV and the innate immune response of host cells in the early stage of infection.

The correlation between TRIM28 expression and immune checkpoints in CRPC.

This study delves into the unexplored realm of castration-resistant prostate cancer (CRPC) by investigating the role of TRIM28 and its intricate molecular mechanisms using high-throughput single-cell transcriptome sequencing and advanced bioinformatics analysis. Our comprehensive examination unveiled dynamic TRIM28 expression changes, particularly in immune cells such as macrophages and CD8+ T cells within CRPC. Correlation analyses with TCGA data highlighted the connection between TRIM28 and immune checkpoint expression and emphasized its pivotal influence on the quantity and functionality of immune cells. Using TRIM28 knockout mouse models, we identified differentially expressed genes and enriched pathways, unraveling the potential regulatory involvement of TRIM28 in the cGAS-STING pathway. Invitro, experiments further illuminated that TRIM28 knockout in prostate cancer cells induced a notable anti-tumor immune effect by inhibiting M2 macrophage polarization and enhancing CD8+ T cell activity. This impactful discovery was validated in an insitu transplant tumor model, where TRIM28 knockout exhibited a deceleration in tumor growth, reduced proportions of M2 macrophages, and enhanced infiltration of CD8+ T cells. In summary, this study elucidates the hitherto unknown anti-tumor immune role of TRIM28 in CRPC and unravels its potential regulatory mechanism via the cGAS-STING signaling pathway. These findings provide novel insights into the immune landscape of CRPC, offering promising directions for developing innovative therapeutic strategies.

Uniparental silencing of 5S rRNA genes in plant allopolyploids - insights from Cardamine (Brassicaceae).

While the phenomenon of uniparental silencing of 35S rDNA in interspecific hybrids and allopolyploids is well documented, there is a notable absence of information regarding whether such silencing extends to the 5S RNA component of ribosomes. To address this gap in knowledge, we analyzed the 5S and 35S rDNA expression in Cardamine (Brassicaceae) allopolyploids, namely C. × insueta (2n = 3x = 24, genome composition RRA), C. flexuosa (2n = 4x = 32, AAHH), and C. scutata (2n = 4x = 32, PPAA) which share a common diploid ancestor (AA). We employed high-throughput sequencing of transcriptomes and genomes and phylogenetic analyses of 5S rRNA variants. The genomic organization of rDNA was further scrutinized through clustering and fluorescence in situ hybridization. In the C. × insueta allotriploid, we observed uniparental dominant expression of 5S and 35S rDNA loci. In the C. flexuosa and C. scutata allotetraploids, the expression pattern differed, with the 35S rDNA being expressed from the A subgenome, whereas the 5S rDNA was expressed from the partner subgenome. Both C. flexuosa and C. scutata but not C. × insueta showed copy and locus number changes. We conclude that in stabilized allopolyploids, transcription of ribosomal RNA components occurs from different subgenomes. This phenomenon appears to result in the formation of chimeric ribosomes comprising rRNA molecules derived from distinct parental origins. We speculate that the interplay of epigenetic silencing and rDNA rearrangements introduces an additional layer of variation in multimolecule ribosomal complexes, potentially contributing to the evolutionary success of allopolyploids.

Transcriptome analysis and development of EST-SSR markers in the mushroom Auricularia heimuer

Auricularia heimuer, the third most frequently cultivated edible mushroom species worldwide, has high medicinal value. However, a shortage of molecular marker hinders the efficiency and accuracy of genetic breeding efforts for A. heimuer. High-throughput transcriptome sequencing data are essential for gene discovery and molecular markers development. This study aimed to clarify the distribution of SSR loci across the A. heimuer transcriptome and to develop highly informative EST-SSR markers. These tools can be used for phylogenetic analysis, functional gene mining, and molecular marker-assisted breeding of A. heimuer. This study used Illumina high-throughput sequencing technology to obtain A. heimuer transcriptome data. The results revealed 37,538 unigenes in the A. heimuer transcriptome. Of these unigenes, 24,777 (66.01%) were annotated via comparison with the COG, Pfam, and NR databases. Overall, 2510 SSRs were identified from the unigenes, including 6 types of SSRs. The most abundant type of repeats were trinucleotides (1425, 56.77%), followed by mononucleotides (391, 15.58%) and dinucleotides (456, 18.17%). Primer pairs for 102 SSR loci were randomly designed for validity confirmation and polymorphism identification; this process yielded 53 polymorphic EST-SSR markers. Finally, 13 pairs of highly polymorphic EST-SSR primers were used to analyze the genetic diversity and population structure of 52 wild A. heimuer germplasms, revealing that the 52 germplasms could be divided into three categories. These results indicated that SSR loci were abundant in types, numbers, and frequencies, providing a potential basis for germplasm resource identification, genetic diversity analysis, and molecular marker-assisted breeding of A. heimuer.

Panoramic RNA expression of fibroblast growth factors in human glioblastoma tissues and the impact on the survival of patients.

Fibroblast growth factors (FGFs) have a key role in various critical steps of tumor growth and progression through effects on angiogenesis, inflammation and the growth and invasion of malignant cells. Nevertheless, the role of the FGF family in human glioblastoma (GBM) has been rarely studied. The objective of the present study was to assess the RNA expression of all FGF family members in tissues obtained from patients with GBM and to analyze the association between FGF expression and the survival of these patients. For this, the RNA expression of FGF family members in the malignant and proximal tissues of 12 patients with GBM was determined by analyzing high-throughput RNA transcriptome sequencing data uploaded to the National Center for Biotechnology Information database. The relationship between FGF genes and the survival of patients with GBM and glioma was also respectively studied by analyzing data from The Cancer Genome Atlas database using the Gene Expression Profiling Interactive Analysis tool. The results showed that the expression of FGF1, FGF17, FGF20 and FGF22 in GBM tissues was lower than that in adjacent tissues, with a difference of >2 times. Analysis of the overall survival of patients with GBM indicated there were no significant relationships between the expression of FGF1, FGF17, FGF20, FGF22 and overall survival. Analysis of the overall survival of patients with glioma showed that glioma patients with low FGF1 expression achieved a longer survival time than patients with high FGF1 expression; however, high expression of FGF17 and FGF22 indicated a longer survival time. In summary, the results of the present study demonstrated the panoramic expression of FGF family members in patients with GBM, and indicated that FGF1, FGF17 and FGF22 did not affect the survival of patients with GBM, but had a notable influence on the survival of patients with glioma.

Cytochrome P450 CYP736A12 is crucial for Trichoderma asperellum-induced alleviation of phoxim phytotoxicity and reduction of pesticide residue in tomato roots

Trichoderma can enhance the metabolism of organophosphate pesticides in plants, but the mechanism is unclear. Here, we performed high-throughput transcriptome sequencing of roots upon Trichoderma asperellum (TM) inoculation and phoxim (P) application in tomato (Solanum lycopersicum L.). A total of 4059 differentially expressed genes (DEGs) were obtained, including 2110 up-regulated and 1949 down-regulated DEGs in P vs TM+P. COG and KOG analysis indicated that DEGs were mainly enriched in signal transduction mechanisms. We then focused on the pesticide detoxification pathway and screened out cytochrome P450 CYP736A12 as a putative gene for functional analysis. We suppressed the expression of CYP736A12 in tomato plants by virus-induced gene silencing and analyzed tissue-specific phoxim residues, oxidative stress markers, glutathione pool, GST activity and related gene expression. Silencing CYP736A12 significantly increased phoxim residue and induced oxidative stress in tomato plants, by attenuating the TM-induced increased activity of antioxidant and detoxification enzymes, redox homeostasis and transcripts of detoxification genes including CYP724B2, GSH1, GSH2, GR, GPX, GST1, GST2, GST3, and ABC. The study revealed a critical mechanism by which TM promotes the metabolism of phoxim in tomato roots, which can be useful for further understanding the Trichoderma-induced xenobiotic detoxification and improving food safety.

Abstract LB316: Endocrine pancreas alteration in pancreatic ductal adenocarcinoma

Abstract Background: The pancreas is central to energy homeostasis, playing essential roles in digestion and regulation of blood glucose. Exocrine pancreas contains large amounts of of hydrolytic enzyme which is responsible for the digestion. Acinar cells, as well as ductal cells, represent the prime components of exocrine pancreas. They are both supposed to be the origin of pancreatic ductal adenocarcinoma(PDAC) and the researches focused on these exocrine cells, as well as other components in pancreatic cancer, are numerous. Islet cells are crucial for glucose homeostasis and their identity are proved to be altered significantly during diabetes by single-cell RNA-sequencing(scRNA-seq). However, the presence or extent of variation of islet cells in PDAC is not established yet. Moreover, due to the low proportion of endocrine cells in pancreas, it is hard to analyze the data of endocrine cells by using published scRNA-seq data of PDAC. As a consequence of this underrepresentation, the function of endocrine cells in PDAC remains unclear. Methods: We used scRNA-seq, multiplex immunohistochemistry (mIHC), high-throughput spatial transcriptomic sequencing (Stereo-seq) to detect the alteration of endocrine pancreas in PDAC of human and mouse model. KPC mice (KrasLSL-G12Dp53LoxPPdx1-CreER) and pancreatic orthotopic transplantation mice were applied to verify the islet cell alteration in multiple stages of PDAC. We also performed islet isolation and purification to enhance the proportion of endocrine cells when preparing single-cell suspension. Cell co-culture experiments were used to determine the interaction between exocrine-endocrine compartments. Results: Through mIHC of islet cells in PDAC tissues and normal pancreas tissues of mice (KPC vs. C57BL/6J) as well as human (PDAC patients vs. healthy donors), we found the decreased proportions of inulin positive (INS+) cells and increased proportions of glucagon positive (GCG+) cells in PDAC tissues. Furthermore, the proportion of pancreatic polypeptide positive (PPY+) cells was significantly elevated in PDAC tissue in both human and murine. Then, we quantified the double-positive cells and found the proportions of INS+GCG+, INS+PPY+ and GCG+PPY+ were markedly increased in PDAC tissue. Additionally, through the experiments in mouse model, we found the islet cells altered gradually within the multiple stage of PDAC(normal-PanIN-PDAC). Finally, it was found that PDAC cells could induce the transdifferentiation of pancreatic α and β cells into GCG+PPY+ and INS+PPY+ double-positive cells, which further promoted PDAC proliferation in a paracrine-dependent manner and formed a reciprocal interaction. Conclusions: Our study systematically maps the alteration of pancreatic endocrine cells in PDAC and elucidates the potential interaction mechanisms between endocrine and exocrine cells during PDAC carcinogenesis. Meanwhile, it is suggested that pancreatic endocrine cells should be considered as an important component of the PDAC microenvironment. Citation Format: Yuan Chen, Xinpeng Yin, Ruiyuan Xu, Rexiati Ruze, Jianlu Song, Chenglin Hu, Chengcheng Wang, Yupei Zhao. Endocrine pancreas alteration in pancreatic ductal adenocarcinoma [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2024; Part 2 (Late-Breaking, Clinical Trial, and Invited Abstracts); 2024 Apr 5-10; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2024;84(7_Suppl):Abstract nr LB316.

Transcriptomic analysis of the response of Avena sativa to Bacillus amyloliquefaciens DGL1.

Bacillus amyloliquefaciens DGL1, isolated from the arid sandy areas in Dagler, Qinghai Province, China, promotes the growth of Avena sativa variety "Qing Yan 1". To elucidate the transcriptomic changes in the oat root system following interaction with DGL1 and to reveal the molecular mechanism by which DGL1 promotes oat growth, treatment and control groups of oat roots at 2, 4, 8, and 12 h after inoculation with a suspension of strain DGL1 were analyzed using Illumina high-throughput transcriptome sequencing technology. The differentially expressed genes were determined through Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses, and the metabolic pathways and key genes were analyzed. The results showed that 7874, 13,392, 13,169, and 19,026 differentially expressed genes were significantly enriched in the glycolysis/gluconeogenesis pathway, amino acid metabolism, nitrogen metabolism, plant hormone signal transduction, and other related metabolic pathways in the oat roots at 2, 4, 8, and 12 h after inoculation with a DGL1 suspension. The GO and KEGG enrichment analyses revealed that the genes encoding plasma membrane ATPase, phosphoglycerate kinase gene PGK, ammonium transporter protein gene AMT, cellulose synthase gene CSLF6, and growth hormone response family gene IAA18 were significantly upregulated. It is hypothesized that the pro-growth mechanism of strain DGL1 in oats is the result of the coordination of multiple pathways through the promotion of oat energy metabolism, phytohormone signaling, secondary metabolite synthesis, and amino acid metabolism.

Prognostic Value of IGF2BP3 Gene Expression Levels in Patients with Acute Myeloid Leukemia

To investigate the relationship between IGF2BP3 gene expression and prognosis in patients with acute myeloid leukemia (AML). High throughput transcriptome sequencing was performed on bone marrow primary leukemia cells from 27 patients with AML in our center, the relationship between IGF2BP3 expression levels and clinical characteristics were analyzed and verify the samples from patients with newly treated AML and refractory AML. The expression level of IGF2BP3 gene were analyzed in 20 healthy subjects and 26 patients with AML. The expression of IGF2BP3 in two anthracycline-resistant cell lines (HL60/ADR, K562/ADR) was detected by RT-qPCR and Western blot, and the expression difference of IGF2BP3 was compared with that in sensitive cells (HL60, K562). The relationship between the expression level of IGF2BP3 in patients with AML and prognostic were analyzed through data analysis of 746 patients with AML, and the prognostic value of IGF2BP3 in AML was analyzed by multivariate Cox regression analysis. In the bone marrow primary leukemia cells of 27 AML patients in our center, the expression level of IGF2BP3 in patients with refractory AML was significantly higher than that in chemotherapy sensitive patients (P =0.0343). The expression of IGF2BP3 in leukemia patients with extramedullary infiltration (EMI) was significantly higher than that in AML patients without extramedullary infiltration (P =0.0049). Compared with healthy subjects (n=20), IGF2BP3 expression in AML patients (n=26) was higher (P =0.0009). The expression of IGF2BP3 mRNA in the anthracycline resistant cell lines (HL60/ADR, K562/ADR) was significantly higher than that in the sensitive cell lines (K562/ADR vs K562,P =0.0430; HL60/ADR vs HL60, P =0.7369). Western blot results showed that the expression of IGF2BP3 protein in mycin resistant cells was significantly higher than that in sensitive cells (P < 0.001). qPCR results showed that the expression level of IGF2BP3 mRNA in refractory AML patients was significantly higher than that in patients with chemotherapy sensitive (P =0.002). High expression of IGF2BP3 was associated with poor prognosis in AML (P < 0.05) in 3 large sample cohorts of AML patients. Univariate and multivariate prognostic analyses demonstrated that high expression of IGF2BP3 was significantly associated with shorter event-free survival (EFS, HR=1.887, P =0.024) and overall survival (OS, HR=1.619, P =0.016). The high expression of IGF2BP3 gene may be an important factor in the poor prognosis of AML, suggesting that IGF2BP3 gene may be a new molecular marker for the clinical prognosis evaluation and treatment strategy of AML.

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.

Alterations in HCN1 expression and distribution during epileptogenesis in rats

BackgroundThe hyperpolarization-activated cyclic nucleotide-gated cation channel (HCN1) is predominantly located in key regions associated with epilepsy, such as the neocortex and hippocampus. Under normal physiological conditions, HCN1 plays a crucial role in the excitatory and inhibitory regulation of neuronal networks. In temporal lobe epilepsy, the expression of HCN1 is decreased in the hippocampi of both animal models and patients. However, whether HCN1 expression changes during epileptogenesis preceding spontaneous seizures remains unclear. ObjectiveThe aim of this study was to determine whether the expression of HCN1 is altered during the epileptic prodromal phase, thereby providing evidence for its role in epileptogenesis. MethodsWe utilized a cobalt wire-induced rat epilepsy model to observe changes in HCN1 during epileptogenesis and epilepsy. Additionally, we also compared HCN1 alterations in epileptogenic tissues between cobalt wire- and pilocarpine-induced epilepsy rat models. Long-term video EEG recordings were used to confirm seizures development. Transcriptional changes, translation, and distribution of HCN1 were assessed using high-throughput transcriptome sequencing, total protein extraction, membrane and cytoplasmic protein fractionation, western blotting, immunohistochemistry, and immunofluorescence techniques. ResultsIn the cobalt wire-induced rat epilepsy model during the epileptogenesis phase, total HCN1 mRNA and protein levels were downregulated. Specifically, the membrane expression of HCN1 was decreased, whereas cytoplasmic HCN1 expression showed no significant change. The distribution of HCN1 in the distal dendrites of neurons decreased. During the epilepsy period, similar HCN1 alterations were observed in the neocortex of rats with cobalt wire-induced epilepsy and hippocampus of rats with lithium pilocarpine-induced epilepsy, including downregulation of mRNA levels, decreased total protein expression, decreased membrane expression, and decreased distal dendrite expression. ConclusionsAlterations in HCN1 expression and distribution are involved in epileptogenesis beyond their association with seizure occurrence. Similarities in HCN1 alterations observed in epileptogenesis-related tissues from different models suggest a shared pathophysiological pathway in epileptogenesis involving HCN1 dysregulation. Therefore, the upregulation of HCN1 expression in neurons, maintenance of the HCN1 membrane, and distal dendrite distribution in neurons may represent promising disease-modifying strategies in epilepsy.

Transcriptome profiling and analysis of patients with esophageal squamous cell carcinoma from Kazakhstan.

Esophageal squamous cell carcinoma (ESCC) is the predominant subtype of esophageal cancer in Central Asia, often diagnosed at advanced stages. Understanding population-specific patterns of ESCC is crucial for tailored treatments. This study aimed to unravel ESCC's genetic basis in Kazakhstani patients and identify potential biomarkers for early diagnosis and targeted therapies. ESCC patients from Kazakhstan were studied. We analyzed histological subtypes and conducted in-depth transcriptome sequencing. Differential gene expression analysis was performed, and significantly dysregulated pathways were identified using KEGG pathway analysis (p-value < 0.05). Protein-protein interaction networks were constructed to elucidate key modules and their functions. Among Kazakhstani patients, ESCC with moderate dysplasia was the most prevalent subtype. We identified 42 significantly upregulated and two significantly downregulated KEGG pathways, highlighting molecular mechanisms driving ESCC pathogenesis. Immune-related pathways, such as viral protein interaction with cytokines, rheumatoid arthritis, and oxidative phosphorylation, were elevated, suggesting immune system involvement. Conversely, downregulated pathways were associated with extracellular matrix degradation, crucial in cancer invasion and metastasis. Protein-protein interaction network analysis revealed four distinct modules with specific functions, implicating pathways in esophageal cancer development. High-throughput transcriptome sequencing elucidated critical molecular pathways underlying esophageal carcinogenesis in Kazakhstani patients. Insights into dysregulated pathways offer potential for early diagnosis and precision treatment strategies for ESCC. Understanding population-specific patterns is essential for personalized approaches to ESCC management.

Pan-histone deacetylase inhibitor vorinostat suppresses osteoclastic bone resorption through modulation of RANKL-evoked signaling and ameliorates ovariectomy-induced bone loss

BackgroundEstrogen deficiency-mediated hyperactive osteoclast represents the leading role during the onset of postmenopausal osteoporosis. The activation of a series of signaling cascades triggered by RANKL-RANK interaction is crucial mechanism underlying osteoclastogenesis. Vorinostat (SAHA) is a broad-spectrum pan-histone deacetylase inhibitor (HDACi) and its effect on osteoporosis remains elusive.MethodsThe effects of SAHA on osteoclast maturation and bone resorptive activity were evaluated using in vitro osteoclastogenesis assay. To investigate the effect of SAHA on the osteoclast gene networks during osteoclast differentiation, we performed high-throughput transcriptome sequencing. Molecular docking and the assessment of RANKL-induced signaling cascades were conducted to confirm the underlying regulatory mechanism of SAHA on the action of RANKL-activated osteoclasts. Finally, we took advantage of a mouse model of estrogen-deficient osteoporosis to explore the clinical potential of SAHA.ResultsWe showed here that SAHA suppressed RANKL-induced osteoclast differentiation concentration-dependently and disrupted osteoclastic bone resorption in vitro. Mechanistically, SAHA specifically bound to the predicted binding site of RANKL and blunt the interaction between RANKL and RANK. Then, by interfering with downstream NF-κB and MAPK signaling pathway activation, SAHA negatively regulated the activity of NFATc1, thus resulting in a significant reduction of osteoclast-specific gene transcripts and functional osteoclast-related protein expression. Moreover, we found a significant anti-osteoporotic role of SAHA in ovariectomized mice, which was probably realized through the inhibition of osteoclast formation and hyperactivation.ConclusionThese data reveal a high affinity between SAHA and RANKL, which results in blockade of RANKL-RANK interaction and thereby interferes with RANKL-induced signaling cascades and osteoclastic bone resorption, supporting a novel strategy for SAHA application as a promising therapeutic agent for osteoporosis.

Whole transcriptome scanning and validation of negatively related genes in UC-MSCs

BackgroundHuman umbilical cord mesenchymal stem cells (UC-MSCs) are one of the most extensively researched stem cell types due to their potential for multi-lineage differentiation, secretion of regenerative factors, modulations of immunological activities, and the release of regenerative substances and influence immunological processes. Since UC-MSCs must be cultivated on a large scale for clinical use, selecting the appropriate storing passage, such as the usage-based passage of UC-MSCs, is critical for long-term autologous or allogeneic usage. Long-term cultivation of stem cells, on the other hand, causes them to lose their pluripotent differentiation capacity. As a result, distinguishing between high and low passages of UC-MSCs and identifying the particular variations associated with stem cells and their modes of action is essential for regenerative medicine. Therefore, we investigated the biological features and transcriptional changes of UC-MSCs over passages. MethodsUC-MSCs were isolated from the tissues of the human umbilical cord, and UC-MSCs from five passages (P1, P3, P5, P10 and P15) with three repetitions were compared and identified based on morphology, cell markers, differentiation capacity, and aging-related characteristics. It was previously assumed that the phenotype of cells before the P10 passage was stable, defined as early passage, and that culture could be continued until the 15th passage, defined as late passage. Next, the five passages of UC-MSCs were sequenced using high-throughput complete transcriptome sequencing. Fuzzy C-Means Clustering (FCM) and Weighted Gene Co-expression Network Analysis (WGCNA) were used to find hub genes, and gene silencing was performed to investigate the impact of missing genes on the stemness of UC-MSC cells. ResultsUC-MSCs of different passages displayed similar surface markers, including CD73, CD105, CD90, CD34, CD45 and HLA-DR. However, the proliferation time of late-phase UC-MSCs was longer than that of early-phase UC-MSCs, and the expression of the senescence-associated (SA)-β-gal staining marker was higher. At the same time, pluripotency markers (NANOG, OCT4, SOX2 and KIF4A) were down-regulated, and the multi-differentiation potential was reduced. Meanwhile, KIFC1 and UBE2C were down-regulated in late-phase UC-MSCs, which were involved in the maintenance of stemness. ConclusionsKIFC1 and UBE2C were highly expressed in early-UC-MSCs and showed a downward gradient trend with cell expansion in vitro. They regulated UC-MSC proliferation, colony sphere formation, multiple differentiation, stemness maintenance, and other biological manifestations. Therefore, they are anticipated to be new biomarkers for UC-MSCs quality identification in regenerative medicine applications.

Color components determination and full-length comparative transcriptomic analyses reveal the potential mechanism of carotenoid synthesis during Paphiopedilum armeniacum flowering.

Paphiopedilum armeniacum (P. armeniacum), an ornamental plant native to China, is known for its distinctive yellow blossoms. However, the mechanisms underlying P. armeniacum flower coloration remain unclear. We selected P. armeniacum samples from different flowering stages and conducted rigorous physicochemical analyses. The specimens were differentiated based on their chemical properties, specifically their solubilities in polar solvents. This key step enabled us to identify the main metabolite of flower color development of P. armeniacum, and to complete the identification by High-performance liquid chromatography (HPLC) based on the results. Additionally, we employed a combined approach, integrating both third-generation full-length transcriptome sequencing and second-generation high-throughput transcriptome sequencing, to comprehensively explore the molecular components involved. We combined physical and chemical analysis with transcriptome sequencing to reveal that carotenoid is the main pigment of P. armeniacum flower color. Extraction colorimetric method and HPLC were used to explore the characteristics of carotenoid accumulation during flowering. We identified 28 differentially expressed carotenoid biosynthesis genes throughout the flowering process, validated their expression through fluorescence quantification, and discovered 19 potential positive regulators involved in carotenoid synthesis. Among these candidates, three RCP2 genes showed a strong potential for governing the PDS and ZDS gene families. In summary, our study elucidates the fundamental mechanisms governing carotenoid synthesis during P. armeniacum flowering, enhancing our understanding of this process and providing a foundation for future research on the molecular mechanisms driving P. armeniacum flowering.

Uncovering lupus nephritis-specific genes and the potential of TNFRSF17-targeted immunotherapy: a high-throughput sequencing study.

Lupus nephritis (LN) is a severe manifestation of systemic lupus erythematosus (SLE). This study aimed to identify LN specific-genes and potential therapeutic targets. We performed high-throughput transcriptome sequencing on peripheral blood mononuclear cells (PBMCs) from LN patients. Healthy individuals and SLE patients without LN were used as controls. To validate the sequencing results, qRT-PCR was performed for 5 upregulated and 5 downregulated genes. Furthermore, the effect of the TNFRSF17-targeting drug IBI379 on patient plasma cells and B cells was evaluated by flow cytometry. Our analysis identified 1493 and 205 differential genes in the LN group compared to the control and SLE without LN groups respectively, with 70 genes common to both sets, marking them as LN-specific. These LN-specific genes were significantly enriched in the 'regulation of biological quality' GO term and the cell cycle pathway. Notably, several genes including TNFRSF17 were significantly overexpressed in the kidneys of both LN patients and NZB/W mice. TNFRSF17 levels correlated positively with urinary protein levels, and negatively with complement C3 and C4 levels in LN patients. The TNFRSF17-targeting drug IBI379 effectively induced apoptosis in patient plasma cells without significantly affecting B cells. Our findings suggest that TNFRSF17 could serve as a potential therapeutic target for LN. Moreover, IBI379 is presented as a promising treatment option for LN.

Uncovering the potential molecular mechanism of liraglutide to alleviate the effects of high glucose on myoblasts based on high-throughput transcriptome sequencing technique

BackgroundMyoblasts play an important role in muscle growth and repair, but the high glucose environment severely affects their function. The purpose of this study is to explore the potential molecular mechanism of liraglutide in alleviating the effects of high glucose environments on myoblasts.MethodsMTT, western blot, and ELISA methods were used to investigate the role of liraglutide on C2C12 myoblasts induced by high glucose. The high-throughput transcriptome sequencing technique was used to sequence C2C12 myoblasts from different treated groups. The DESeq2 package was used to identify differentially expressed-mRNAs (DE-mRNAs). Then, functional annotations and alternative splicing (AS) were performed. The Cytoscape-CytoHubba plug-in was used to identify multicentric DE-mRNAs.ResultsThe MTT assay results showed that liraglutide can alleviate the decrease of myoblasts viability caused by high glucose. Western blot and ELISA tests showed that liraglutide can promote the expression of AMPKα and inhibit the expression of MAFbx, MuRF1 and 3-MH in myoblasts. A total of 15 multicentric DE-mRNAs were identified based on the Cytoscape-CytoHubba plug-in. Among them, Top2a had A3SS type AS. Functional annotation identifies multiple signaling pathways such as metabolic pathways, cytokine-cytokine receptor interaction, cAMP signaling pathway and cell cycle.ConclusionLiraglutide can alleviate the decrease of cell viability and degradation of muscle protein caused by high glucose, and improves cell metabolism and mitochondrial activity. The molecular mechanism of liraglutide to alleviate the effect of high glucose on myoblasts is complex. This study provides a theoretical basis for the clinical effectiveness of liraglutide in the treatment of skeletal muscle lesions in diabetes.

Identification of Key Genes Associated with 1,2,6-Tri-O-galloyl-β-D-glucopyranose Accumulation in Camellia sinensis Based on Transcriptome Sequencing

Hydrolyzed tannin 1,2,6-tri-O-galloyl-β-D-glucopyranose (1,2,6-TGGP) possesses significant medicinal properties. However, little is known about its underlying molecular mechanisms. In this study, the levels of 1,2,6-TGGP in tea materials from different cultivars and leaf positions were compared. Additionally, one leaf and one bud sample from six tea cultivars with significant variations in 1,2,6-TGGP levels were analyzed using transcriptome high-throughput sequencing to identify the genes that are responsible for 1,2,6-TGGP accumulation. The sequencing results were mapped to the reference tea genome, revealing a total of 2735 differentially expressed genes (DEGs). This set included four UDP glycosyltransferase (UGTs) and six serine carboxypeptidases-like (SCPLs) genes. Among them, the upregulated SCPLs (CSS0032817) may directly participate in the acylation reaction of 1,2,6-TGGP. In addition, several classes of DEGs, including cytochrome P450, were significantly associated with the 1,2,6-TGGP content, which is potentially involved in their regulation. Overall, these results provide new insights into the molecular mechanism of 1,2,6-TGGP accumulation.

Genome-wide identification of the B3 transcription factor family in pepper (Capsicum annuum) and expression patterns during fruit ripening

In plants, B3 transcription factors play important roles in a variety of aspects of their growth and development. While the B3 transcription factor has been extensively identified and studied in numerous species, there is limited knowledge regarding its B3 superfamily in pepper. Through the utilization of genome-wide sequence analysis, we identified a total of 106 B3 genes from pepper (Capsicum annuum), they are categorized into four subfamilies: RAV, ARF, LAV, and REM. Chromosome distribution, genetic structure, motif, and cis-acting element of the pepper B3 protein were analyzed. Conserved gene structure and motifs outside the B3 domain provided strong evidence for phylogenetic relationships, allowing potential functions to be deduced by comparison with homologous genes from Arabidopsis. According to the high-throughput transcriptome sequencing analysis, expression patterns differ during different phases of fruit development in the majority of the 106 B3 pepper genes. By using qRT-PCR analysis, similar expression patterns in fruits from various time periods were discovered. In addition, further analysis of the CaRAV4 gene showed that its expression level decreased with fruit ripening and located in the nucleus. B3 transcription factors have been genome-wide characterized in a variety of crops, but the present study is the first genome-wide analysis of the B3 superfamily in pepper. More importantly, although B3 transcription factors play key regulatory roles in fruit development, it is uncertain whether B3 transcription factors are involved in the regulation of the fruit development and ripening process in pepper and their specific regulatory mechanisms because the molecular mechanisms of the process have not been fully explained. The results of the study provide a foundation and new insights into the potential regulatory functions and molecular mechanisms of B3 genes in the development and ripening process of pepper fruits, and provide a solid theoretical foundation for the enhancement of the quality of peppers and their selection and breeding of high-yield varieties.

Transcriptome Analysis and Identification of Cadmium-Induced Oxidative Stress Response Genes in Different Meretrix meretrix Developmental Stages.

Cadmium (Cd) is one of the major pollutants in the aquatic environment, and it can easily accumulate in aquatic animals and result in toxic effects by changing the metabolism of the body, causing a serious impact on the immune system, reproductive system, and the development of offspring. The clam Meretrix meretrix is one of the commercially important species that is cultivated in large-scale aquaculture in China. To elucidate the underlying molecular mechanisms of Cd2+ in the developmental processes, fertilized eggs and larvae of M. meretrix at different developmental stages were exposed to Cd2+ (27.2 mg L-1 in natural seawater) or just natural seawater without Cd2+ (control), and high-throughput transcriptome sequencing and immunohistochemistry techniques were used to analyze the toxic effects of Cd on larvae at different early developmental stages. The results revealed 31,914 genes were differentially expressed in the different stages of M. meretrix development upon treatment with Cd2+. Ten of these genes were differentially expressed in all stages of development examined, but they comprised only six unigenes (CCO, Ndh, HPX, A2M, STF, and pro-C3), all of which were related to the oxidative stress response. Under Cd exposure, the expression levels of CCO and Ndh were significantly upregulated in D-shaped and pediveliger larvae, while pro-C3 expression was significantly upregulated in the fertilized egg, D-shaped larva, and pediveliger. Moreover, HPX, A2M, and STF expression levels in the fertilized egg and pediveliger larvae were also significantly upregulated. In contrast, CCO, Ndh, HPX, A2M, STF, and pro-C3 expression levels in the postlarva were all downregulated under Cd exposure. Besides the genes with changes in expression identified by the transcriptome, the expression of two other oxidative stress-related genes (MT and Nfr2) was also found to change significantly in the different developmental stages of M. meretrix upon Cd exposure, confirming their roles in combating oxidative stress. Overall, the findings of this study indicated that Cd would interfere with cellular respiration, ion transport, and immune response through inducing oxidative stress, and changes in the expression of oxidative stress-related genes might be an important step for M. meretrix to deal with the adverse effects of Cd at different stages of its development.