Nasopharyngeal carcinoma（NPC） is a distinct type of head and neck cancer closely associated with Epstein-Barr virus（EBV） infection and exhibits significant geographic variations in its incidence. Despite recent advancements in radiotherapy techniques and precision medicine for NPC, the overall survival rate remains unsatisfactory due to tumor metastasis, recurrence, and drug resistance. Single-cell RNA sequencing（scRNA-seq） is an emerging technology that allows for the analysis of gene expression at single-cell resolution, providing a clearer understanding of tumor cell subpopulations, the evolutionary trajectory of tumor cells, and the functional roles and interactions of cells within the tumor microenvironment. This provides new ideas for the development of precision medicine in NPC. Here, we review the applications of scRNA-seq in exploring the mechanisms of NPC pathogenesis, tumor heterogeneity, the tumor microenvironment, drug resistance, and therapeutic response.
Humans ; Nasopharyngeal Neoplasms/genetics* ; Tumor Microenvironment ; Nasopharyngeal Carcinoma ; Single-Cell Analysis ; Sequence Analysis, RNA ; Precision Medicine ; Drug Resistance, Neoplasm ; Epstein-Barr Virus Infections ; Herpesvirus 4, Human ; Single-Cell Gene Expression Analysis

The regulatory processes in developmental biology research are significantly influenced by long non-coding RNAs (lncRNAs). However, the dynamics of lncRNA expression during human tooth development remain poorly understood. In this research, we examined the lncRNAs present in the dental epithelium (DE) and dental mesenchyme (DM) at the late bud, cap, and early bell stages of human fetal tooth development through bulk RNA sequencing. Developmental regulators co-expressed with neighboring lncRNAs were significantly enriched in odontogenesis. Specific lncRNAs expressed in the DE and DM, such as PANCR, MIR205HG, DLX6-AS1, and DNM3OS, were identified through a combination of bulk RNA sequencing and single-cell analysis. Further subcluster analysis revealed lncRNAs specifically expressed in important regions of the tooth germ, such as the inner enamel epithelium and coronal dental papilla (CDP). Functionally, we demonstrated that CDP-specific DLX6-AS1 enhanced odontoblastic differentiation in human tooth germ mesenchymal cells and dental pulp stem cells. These findings suggest that lncRNAs could serve as valuable cell markers for tooth development and potential therapeutic targets for tooth regeneration.
Humans ; RNA, Long Noncoding/metabolism* ; Odontogenesis/genetics* ; Tooth Germ/embryology* ; Cell Differentiation ; Gene Expression Regulation, Developmental ; Mesoderm/metabolism* ; Tooth/embryology* ; Gene Expression Profiling ; Sequence Analysis, RNA ; Dental Pulp/cytology*

Based on whole genome data, the identification and expression pattern analysis of the Atropa belladonna WRKY transcription factor family were conducted to provide a theoretical foundation for studying the biological functions and mechanisms of these transcription factors. In this study, bioinformatics methods were employed to identify members of the A. belladonna WRKY gene family and to predict their physicochemical properties, conserved motifs, promoter cis-acting elements, and chromosomal localization. Additionally, the expression patterns of the A. belladonna WRKY gene family under the regulation of environmental factors such as light quality and temperature were analyzed. The results revealed a total of 28 AbWRKY transcription factors, randomly distributed across 16 chromosomes, encoding 324-707 amino acids. Most AbWRKY proteins were acidic, unstable, and hydrophilic. Based on multiple sequence alignment and phylogenetic analysis, the WRKY gene family members were classified into two subfamilies. Conserved motif and domain analysis indicated that WRKY transcription factors in the same subfamily possessed conserved structural features. Promoter analysis predicted that the A. belladonna WRKY family contained light-responsive elements, hormone-responsive elements, and stress-responsive elements. Collinearity analysis showed that AbWRKY24 plays a crucial role in the expansion of the AbWRKY gene family. Then qRT-PCR results indicated that AbWRKY6, AbWRKY8, AbWRKY14, and AbWRKY24 responded to red light stress, while AbWRKY8, AbWRKY14, and AbWRKY24 responded to yellow light/low-temperature combined stress. AbWRKY6 and AbWRKY8 were significantly expressed in leaves and stems, AbWRKY27 and AbWRKY28 were significantly expressed in fibrous roots, and AbWRKY25 was significantly expressed in flowers. This study is the first to identify and analyze the WRKY gene family in A. belladonna and to examine its expression patterns under light and temperature regulation, laying a foundation for in-depth analysis and functional validation of the molecular mechanisms of A. belladonna WRKY transcription factors in responding to light quality and temperature environmental factors.
Transcription Factors/chemistry* ; Plant Proteins/metabolism* ; Phylogeny ; Gene Expression Regulation, Plant ; Light ; Temperature ; Atropa belladonna/metabolism* ; Multigene Family/genetics* ; Promoter Regions, Genetic/genetics* ; Sequence Alignment ; Amino Acid Sequence ; Genome, Plant/genetics*

Spinal cord injury (SCI) represents a complex pathophysiological process involving the interaction of multiple cell types. Conventional sequencing methods can only detect the average gene expression level of the damaged local cell populations, which is difficult to reflect its heterogeneity. Therefore, new technologies are needed to reveal the intercellular heterogeneity and the complex intercellular interactions of the damaged lesions. The single-cell RNA sequencing (scRNA-seq) technique facilitates high-resolution profiling of gene expression at the single-cell level, providing insights into cellular heterogeneity and function, potential molecular pathways, cell fate transitions, and the intercellular interactions pertinent to disease progression. This technology generates valuable gene expression data that support both basic and translational research efforts aiming at the identification of therapeutic targets for intervention. The scRNA-seq technique and its multifaceted application in the local immune microenvironment of injury after SCI were discussed, which will contribute to a more comprehensive understanding of the pathophysiological processes in the immune microenvironment of SCI.
Spinal Cord Injuries/genetics* ; Humans ; Single-Cell Analysis/methods* ; Sequence Analysis, RNA/methods* ; Animals ; Gene Expression Profiling/methods* ; Cellular Microenvironment/genetics*

Gelsemium elegans, a vine plant of Loganiaceae, has both medicinal and forage values. However, it is susceptible to low temperatures during growth. Exploring low temperature response genes is of great significance for cold resistance breeding of G. elegans. Ethylene response factors (ERFs) are the transcription factors of the AP2/ERF superfamily and play a crucial role in plant stress response. In this study, based on the unigene GeERF involved in the response to low temperature stress in the transcriptome of G. elegans, a full-length cDNA sequence of the transcription factor GeERF4B-1 was cloned from the leaves of G. elegans by reverse transcription-polymerase chain reaction (RT-PCR). Bioinformatics analysis showed that GeERF4B-1 had an open reading frame of 759 bp, encoding a protein composed of 252 amino acid residues and with a relative molecular weight of 27 kDa. The deduced protein was predicted to be an unstable, alkaline, and hydrophilic protein. The phylogenetic tree showed that GeERF4B-1 was in the same clade as the B-4 subfamily of the ERF family. The results of the subcellular localization experiment revealed that GeERF4B-1 was located in the nucleus. Real time quantitative PCR (RT-qPCR) analysis indicated that GeERF4B-1 was expressed in the root, stem, and leaf of G. elegans, with the highest expression level in the root. Compared with the control, the treatments with a low temperature (4 ℃), methyl jasmonate (MeJA), and abscisic acid (ABA) up-regulated the expression level of GeERF4B-1, which reached the peak at 24-48 h. This result revealed that GeERF4B-1 actively responded to low temperature, MeJA, and ABA stresses. However, both sodium chloride (NaCl) and drought treatments down-regulated the expression of GeERF4B-1. In addition, a prokaryotic expression vector of GeERF4B-1 was constructed, and a fusion protein of approximately 52 kDa was yielded after induced expression. The results of the plate stress assay showed that compared with the control, the prokaryotic strain expressing GeERF4B-1 demonstrated enhanced tolerance to low temperatures and sensitivity to salt and mannitol stresses. This study provides theoretical references and potential genetic resources for breeding G. elegans varieties with stress resistance.
Transcription Factors/metabolism* ; Plant Proteins/metabolism* ; Gelsemium/metabolism* ; Acetates/pharmacology* ; Gene Expression Regulation, Plant ; Phylogeny ; Cold Temperature ; Amino Acid Sequence ; Cyclopentanes/metabolism* ; Oxylipins/metabolism* ; Stress, Physiological/genetics* ; Abscisic Acid/metabolism* ; Cloning, Molecular

This study aimed to explore the roles of microRNAs (miRNAs) in the post-transcriptional regulation of cocaine- and amphetamine-regulated transcript (CART) peptide in the bovine hypothalamus and to screen key regulatory miRNAs. Targetscan was used to predict the potential miRNAs binding to CART 3' untranslated regions (3'UTR). Bioinformatics analysis predicted 7 miRNA binding sites in the bovine CART 3'UTR, which were bta-miR-377, bta-miR-331-3p, bta-miR-491, bta-miR-493, bta-miR-758, bta-miR-877, and bta-miR-381, respectively. Reverse transcription-PCR (RT-PCR) was carried out to determine the endogenous expression of CART and target miRNAs in the bovine hypothalamus. All the 7 target miRNAs and CART were endogenously expressed in the bovine hypothalamus. The dual-luciferase reporter gene assay was employed to detect the targeted binding relationship between CART 3'UTR and target miRNAs obtained from bioinformatics analysis. The dual-luciferase reporter gene assay confirmed that the 3'UTR of CART had a targeted binding relationship with the 7 target miRNAs. Cell experiments were conducted to examine the effects of target miRNAs on the messenger RNA (mRNA) and protein levels of exogenous CART and screen for key regulatory miRNAs. The results of cell experiments showed that the 7 miRNAs downregulated the mRNA level of CART, with bta-miR-491 demonstrating the strongest downregulating effect. Bta-miR-377, bta-miR-331-3p, bta-miR-491, bta-miR-493, and bta-miR-381 downregulated the protein level of CART, with bta-miR-381 exerting the strongest downregulating effect. Animal experiments were conducted to explore the effects of key regulatory miRNAs on the mRNA and protein levels of CART in the hypothalamus and the CART concentration in the serum. The results from animal experiments showed that miR-491 and miR-381 regulated the endogenous expression of CART in the hypothalamus and the concentration in the serum by binding to the CART 3'UTR. These results suggest that miR-491 and miR-381 are the main miRNAs regulating CART expression in the bovine hypothalamus, which can affect serum CART concentration by modulating endogenous CART expression.
Animals ; MicroRNAs/metabolism* ; Cattle ; Hypothalamus/metabolism* ; 3' Untranslated Regions/genetics* ; Nerve Tissue Proteins/metabolism* ; Gene Expression Regulation ; Binding Sites ; Base Sequence ; Computational Biology/methods* ; Cocaine- and Amphetamine-Regulated Transcript Protein

Dental primary afferent (DPA) neurons and proprioceptive mesencephalic trigeminal nucleus (MTN) neurons, located in the trigeminal ganglion and the brainstem, respectively, are essential for controlling masticatory functions. Despite extensive transcriptomic studies on various somatosensory neurons, there is still a lack of knowledge about the molecular identities of these populations due to technical challenges in their circuit-validated isolation. Here, we employed high-depth single-cell RNA sequencing (scRNA-seq) in combination with retrograde tracing in mice to identify intrinsic transcriptional features of DPA and MTN neurons. Our transcriptome analysis revealed five major types of DPA neurons with cell type-specific gene enrichment, some of which exhibit unique mechano-nociceptive properties capable of transmitting nociception in response to innocuous mechanical stimuli in the teeth. Furthermore, we discovered cellular heterogeneity within MTN neurons that potentially contribute to their responsiveness to mechanical stretch in the masseter muscle spindles. Additionally, DPA and MTN neurons represented sensory compartments with distinct molecular profiles characterized by various ion channels, receptors, neuropeptides, and mechanoreceptors. Together, our study provides new biological insights regarding the highly specialized mechanosensory functions of DPA and MTN neurons in pain and proprioception.
Animals ; Mice ; Neurons ; Proprioception ; Gene Expression Profiling ; Pain ; Sequence Analysis, RNA

Flavanone 3-hydroxylase (F3H) is a key enzyme in the synthesis of phycocyanidins. In this experiment, the petals of red Rhododendron hybridum Hort. at different developmental stages were used as experimental materials. The R. hybridum flavanone 3-hydroxylase (RhF3H) gene was cloned using reverse transcription PCR (RT-PCR) and rapid-amplification of cDNA ends (RACE) techniques, and bioinformatics analyses were performed. Petal RhF3H gene expression at different developmental stages were analyzed by using quantitative real-time polymerase chain reaction (qRT-PCR). A pET-28a-RhF3H prokaryotic expression vector was constructed for the preparation and purification of RhF3H protein. A pCAMBIA1302-RhF3H overexpression vector was constructed for genetic transformation in Arabidopsis thaliana by Agrobacterium-mediated method. The results showed that the R. hybridum Hort. RhF3H gene is 1 245 bp long, with an open reading frame of 1 092 bp, encoding 363 amino acids. It contains a Fe²⁺ binding motif and a 2-ketoglutarate binding motif of the dioxygenase superfamily. Phylogenetic analysis showed that the R. hybridum RhF3H protein is most closely related to the Vaccinium corymbosum F3H protein. qRT-PCR analysis showed that the expression level of the red R. hybridum RhF3H gene tended to increase and then decrease in the petals at different developmental stages, with the highest expression at middle opening stage. The results of the prokaryotic expression showed that the size of the induced protein of the constructed prokaryotic expression vector pET-28a-RhF3H was about 40 kDa, which was similar to the theoretical value. Transgenic RhF3H Arabidopsis thaliana plants were successfully obtained, and PCR identification and β-glucuronidase (GUS) staining demonstrated that the RhF3H gene was integrated into the genome of A. thaliana plants. qRT-PCR, total flavonoid and anthocyanin contentanalysis showed that RhF3H was significantly higher expressed in the transgenic A. thaliana relative to that of the wild type, and its total flavonoid and anthocyanin content were significantly increased. This study provides a theoretical basis for investigating the function of RhF3H gene, as well as for studying the molecular mechanism of flower color in R. simsiib Planch.
Arabidopsis/metabolism* ; Rhododendron/metabolism* ; Amino Acid Sequence ; Anthocyanins/metabolism* ; Phylogeny ; Flavonoids/metabolism* ; Cloning, Molecular ; Gene Expression Regulation, Plant ; Plant Proteins/metabolism*

Transcriptome analysis based on high-throughput sequencing of a cDNA library has been widely applied to functional genomic studies. However, the cDNA dependence of most RNA sequencing techniques constrains their ability to detect base modifications on RNA, which is an important element for the post-transcriptional regulation of gene expression. To comprehensively profile the N⁶-methyladenosine (m⁶A) and N⁵-methylcytosine (m⁵C) modifications on RNA, direct RNA sequencing (DRS) using the latest Oxford Nanopore Technology was applied to analyze the transcriptome of six tissues in rice. Approximately 94 million reads were generated, with an average length ranging from 619 nt to 1013 nt, and a total of 45,707 transcripts across 34,763 genes were detected. Expression profiles of transcripts at the isoform level were quantified among tissues. Transcriptome-wide mapping of m⁶A and m⁵C demonstrated that both modifications exhibited tissue-specific characteristics. The transcripts with m⁶A modifications tended to be modified by m⁵C, and the transcripts with modifications presented higher expression levels along with shorter poly(A) tails than transcripts without modifications, suggesting the complexity of gene expression regulation. Gene Ontology analysis demonstrated that m⁶A- and m⁵C-modified transcripts were involved in central metabolic pathways related to the life cycle, with modifications on the target genes selected in a tissue-specific manner. Furthermore, most modified sites were located within quantitative trait loci that control important agronomic traits, highlighting the value of cloning functional loci. The results provide new insights into the expression regulation complexity and data resource of the transcriptome and epitranscriptome, improving our understanding of the rice genome.
RNA ; Oryza/genetics* ; RNA, Messenger ; Gene Expression Profiling ; Transcriptome ; Sequence Analysis, RNA ; High-Throughput Nucleotide Sequencing/methods* ; RNA Processing, Post-Transcriptional

Tumors are complex ecosystems in which heterogeneous cancer cells interact with their microenvironment composed of diverse immune, endothelial, and stromal cells. Cancer biology had been studied using bulk genomic and gene expression profiling, which however mask the cellular diversity and average the variability among individual molecular programs. Recent advances in single-cell transcriptomic sequencing have enabled a detailed dissection of tumor ecosystems and promoted our understanding of tumorigenesis at single-cell resolution. In the present review, we discuss the main topics of recent cancer studies that have implemented single-cell RNA sequencing (scRNA-seq). To study cancer cells, scRNA-seq has provided novel insights into the cancer stem-cell model, treatment resistance, and cancer metastasis. To study the tumor microenvironment, scRNA-seq has portrayed the diverse cell types and complex cellular states of both immune and non-immune cells interacting with cancer cells, with the promise to discover novel targets for future immunotherapy.
Ecosystem ; Gene Expression Profiling ; Genomics ; Humans ; Neoplasms/pathology* ; Sequence Analysis, RNA ; Single-Cell Analysis ; Transcriptome ; Tumor Microenvironment/genetics*