Study on transcriptome characteristics of respiratory syncytial virus bronchiolitis in children by RNA sequencing / 中华预防医学杂志

To explore the biological characteristics related to the pathogenesis and severity of respiratory syncytial virus (RSV) bronchiolitis by RNA sequencing of white blood cells in children with RSV bronchiolitis. This study is a case-control study. A total of 87 children diagnosed with bronchiolitis and RSV antigen positive and/or RSV nucleic acid positive in the pediatric respiratory department of the Second Affiliated Hospital of Wenzhou Medical University from October 2019 to April 2022 were selected as the case group. The case group was divided into three groups based on the condition: mild, moderate, and severe, and there were two groups according to the presence or absence of atopic symptoms: the atopic group and the non-atopic group, forty healthy children in the same period were selected as the control group. The whole blood leukocyte RNA of the children in the case group and the control group was extracted for RNA sequencing, and the data were analyzed to obtain differentially expressed genes (DEGs). Then, the immunobiological pathways and genes related to the pathogenesis, disease condition, and atopy were screened through Gene Ontology (GO) annotation, Kyoto Gene and Genome Encyclopedia (KEGG) annotation, and protein interaction network (PPI) construction methods. Construct the weighted gene co-expression network analysis (WGCNA) module to identify potential biological indicators related to disease severity.Compared with the control group, the case group had a total of 1 782 DEGs, including 1 586 upregulated genes and 196 downregulated genes. The GO pathway enrichment of DEGs is mainly enriched in molecular functions such as peroxidase activity and oxidoreductase activity. In the cytological components, it is mainly enriched in cytoplasmic vesicle lumen and secretory granule lumen. In biological processes, it is mainly enriched in processes such as neutrophil activation involved in immune responses, neutrophil degranulation, and neutrophil activation. KEGG analysis is mainly concentrated in the signal pathway of the viral protein interaction with cytokine and cytokine receptor. A PPI network was constructed to screen four genes at the core position, including CCL2, IL-10, MMP9 and JUN. The DEGs obtained by comparing different disease groups with the control group are mainly enriched in retrograde endocannabinoid signaling and cell apoptosis pathways. WGCNA analysis showed that the brown module related to oxygen saturation was most closely related to the disease, and its gene was mainly enriched in the RNA helicase retinoic acid inducible gene-I (RIG-I) like receptor signal pathway. There are 230 specific DEGs in the atopic group and 444 in the non-atopic group. KEGG enrichment analysis results show that both groups are enriched to NF-κB signaling pathway, the characteristic does not cause significant changes in immune response and transcriptome characteristics in children with RSV bronchiolitis. In conclusion, neutrophil activation, degranulation pathway and signal pathway of interaction between viral protein and cytokine and cytokine receptor are involved in the immune response of RSV bronchiolitis host. CCL2, IL-10, MMP9 and JUN genes may be associated with the pathogenesis. They might be potential biomarkers related to disease severity in RIG-I like receptors, cell apoptosis, and endogenous cannabinoid related signaling pathways.

Study on transcriptome characteristics of respiratory syncytial virus bronchiolitis in children by RNA sequencing / 中华预防医学杂志

To explore the biological characteristics related to the pathogenesis and severity of respiratory syncytial virus (RSV) bronchiolitis by RNA sequencing of white blood cells in children with RSV bronchiolitis. This study is a case-control study. A total of 87 children diagnosed with bronchiolitis and RSV antigen positive and/or RSV nucleic acid positive in the pediatric respiratory department of the Second Affiliated Hospital of Wenzhou Medical University from October 2019 to April 2022 were selected as the case group. The case group was divided into three groups based on the condition: mild, moderate, and severe, and there were two groups according to the presence or absence of atopic symptoms: the atopic group and the non-atopic group, forty healthy children in the same period were selected as the control group. The whole blood leukocyte RNA of the children in the case group and the control group was extracted for RNA sequencing, and the data were analyzed to obtain differentially expressed genes (DEGs). Then, the immunobiological pathways and genes related to the pathogenesis, disease condition, and atopy were screened through Gene Ontology (GO) annotation, Kyoto Gene and Genome Encyclopedia (KEGG) annotation, and protein interaction network (PPI) construction methods. Construct the weighted gene co-expression network analysis (WGCNA) module to identify potential biological indicators related to disease severity.Compared with the control group, the case group had a total of 1 782 DEGs, including 1 586 upregulated genes and 196 downregulated genes. The GO pathway enrichment of DEGs is mainly enriched in molecular functions such as peroxidase activity and oxidoreductase activity. In the cytological components, it is mainly enriched in cytoplasmic vesicle lumen and secretory granule lumen. In biological processes, it is mainly enriched in processes such as neutrophil activation involved in immune responses, neutrophil degranulation, and neutrophil activation. KEGG analysis is mainly concentrated in the signal pathway of the viral protein interaction with cytokine and cytokine receptor. A PPI network was constructed to screen four genes at the core position, including CCL2, IL-10, MMP9 and JUN. The DEGs obtained by comparing different disease groups with the control group are mainly enriched in retrograde endocannabinoid signaling and cell apoptosis pathways. WGCNA analysis showed that the brown module related to oxygen saturation was most closely related to the disease, and its gene was mainly enriched in the RNA helicase retinoic acid inducible gene-I (RIG-I) like receptor signal pathway. There are 230 specific DEGs in the atopic group and 444 in the non-atopic group. KEGG enrichment analysis results show that both groups are enriched to NF-κB signaling pathway, the characteristic does not cause significant changes in immune response and transcriptome characteristics in children with RSV bronchiolitis. In conclusion, neutrophil activation, degranulation pathway and signal pathway of interaction between viral protein and cytokine and cytokine receptor are involved in the immune response of RSV bronchiolitis host. CCL2, IL-10, MMP9 and JUN genes may be associated with the pathogenesis. They might be potential biomarkers related to disease severity in RIG-I like receptors, cell apoptosis, and endogenous cannabinoid related signaling pathways.

Screening of housekeeping genes in Gelsemium elegans and expression patterns of genes involved in its alkaloid biosynthesis / 生物工程学报

Gelsemium elegans is a traditional Chinese herb of medicinal importance, with indole terpene alkaloids as its main active components. To study the expression of the most suitable housekeeping reference genes in G. elegans, the root bark, stem segments, leaves and inflorescences of four different parts of G. elegans were used as materials in this study. The expression stability of 10 candidate housekeeping reference genes (18S, GAPDH, Actin, TUA, TUB, SAND, EF-1α, UBC, UBQ, and cdc25) was assessed through real-time fluorescence quantitative PCR, GeNorm, NormFinder, BestKeeper, ΔCT, and RefFinder. The results showed that EF-1α was stably expressed in all four parts of G. elegans and was the most suitable housekeeping gene. Based on the coexpression pattern of genome, full-length transcriptome and metabolome, the key candidate targets of 18 related genes (AS, AnPRT, PRAI, IGPS, TSA, TSB, TDC, GES, G8H, 8-HGO, IS, 7-DLS, 7-DLGT, 7-DLH, LAMT, SLS, STR, and SGD) involved in the Gelsemium alkaloid biosynthesis were obtained. The expression of 18 related enzyme genes were analyzed by qRT-PCR using the housekeeping gene EF-1α as a reference. The results showed that these genes' expression and gelsenicine content trends were correlated and were likely to be involved in the biosynthesis of the Gelsemium alkaloid, gelsenicine.

Preliminary Results of Single-cell Transcriptome Sequencing in Renal Arterial Lesions of Takayasu Arteritis / 中国医学科学院学报

Objective To explore the preliminary application of single-cell RNA sequencing (scRNA-seq) in the renal arterial lesions in Takayasu arteritis (TA) patients. Methods This study included 2 TA patients with renal artery stenosis treated by bypass surgery in the Department of Vascular Surgery,Beijing Hospital.The obtained 2 renal artery samples were digested with two different protocols (GEXSCOPE kit and self-made digestion liquid) before scRNA-seq and bioinformatics analysis. Results A total of 2920 cells were obtained for further analysis.After unbiased cluster analysis,2 endothelial cell subsets,2 smooth muscle cell subsets,1 fibroblast subset,2 mononuclear macrophage subsets,1 T cell subset,and 1 undefined cell subset were identified.Among them,the two subsets of smooth muscle cells were contractile and secretory,respectively.The results of scRNA-seq indicated that enzymatic hydrolysis with GEXSCOPE kit produced a large number of endothelial cells (57.46%) and a small number of immune cells (13.21%).However,immune cells (34.64%) were dominant in the cells obtained by enzymatic hydrolysis with self-made digestive liquid. Conclusion scRNA-seq can be employed to explore the cellular heterogeneity of diseased vessels in TA patients.Different enzymatic digestion protocols may impact the proportion of different cells.

Analysis of WRKY transcription factor family based on full-length transcriptome sequencing in Polygonatum cyrtonema / 中国中药杂志

WRKY transcription factor family plays an important role in plant growth and development, secondary metabolite synthesis, and biotic and abiotic stress responses. The present study performed full-length transcriptome sequencing of Polygonatum cyrtonema by virtue of the PacBio SMRT high-throughput platform, identified the WRKY family by bioinformatics methods, and analyzed the physicochemical properties, subcellular localization, phylogeny, and conserved motifs. The results showed that 30.69 Gb nucleotide bases and 89 564 transcripts were obtained after redundancy removal. These transcripts had a mean length of 2 060 bp and an N50 value of 3 156 bp. Based on the full-length transcriptome sequencing data, 64 candidate proteins were selected from the WRKY transcription factor family, with the protein size of 92-1 027 aa, the relative molecular mass of 10 377.85-115 779.48 kDa, and the isoelectric point of 4.49-9.84. These WRKY family members were mostly located in the nucleus and belonged to the hydrophobic proteins. According to the phylogenetic analysis of WRKY family in P. cyrtonema and Arabidopsis thaliana, all WRKY family members were clustered into seven subfamilies and WRKY proteins from P. cyrtonema were distributed in different numbers in these seven subgroups. Expression pattern analysis confirmed that 40 WRKY family members had distinct expression patterns in the rhizomes of 1-and 3-year-old P. cyrtonema. Except for PcWRKY39, the expression of 39 WRKY family members was down-regulated in 3-year-old samples. In conclusion, this study provides abundant reference data for genetic research on P. cyrtonema and lays a foundation for the in-depth investigation of the biological functions of the WRKY family.

Mechanism of Huangjing Qianshi Decoction in treatment of prediabetic mice based on transcriptome sequencing / 中国中药杂志

Based on transcriptome sequencing technology, the mouse model of prediabetes treated with Huangjing Qianshi Decoction was sequenced to explore the possible mechanism of treating prediabetes. First of all, transcriptome sequencing was performed on the normal BKS-DB mouse group, the prediabetic model group, and the Huangjing Qianshi Decoction treatment group(treatment group) to obtain differentially expressed genes in the skeletal muscle samples of mice. The serum biochemical indexes were detected in each group to screen out the core genes of Huangjing Qianshi Decoction in prediabetes. Gene Ontology(GO) database and Kyoto Encyclopedia of Genes and Genomes(KEGG) database were used to conduct signaling pathway enrichment analysis of differentially expressed genes, and real-time quantitative polymerase chain reaction(RT-qPCR) was used to verify them. The results showed that the levels of fasting blood glucose(FBG), fasting insulin(FINS), insulin resistance index(HOMA-IR), total cholesterol(TC), triglycerides(TG), and low-density lipoprotein cholesterol(LDL-C) in the mouse model were significantly decreased after treatment with Huangjing Qianshi Decoction. In the results of differential gene screening, there were 1 666 differentially expressed genes in the model group as compared with the normal group, and there were 971 differentially expressed genes in the treatment group as compared with the model group. Among them, interleukin-6(IL-6) and NR3C2 genes, which were closely related to the regulation of insulin resis-tance function, were significantly up-regulated between the model group and the normal group, and vascular endothelial growth factor A(VEGFA) genes were significantly down-regulated between the model group and the normal group. However, the expression results of IL-6, NR3C2, and VEGFA genes were adverse between the treatment group and the model group. GO functional enrichment analysis found that the biological process annotation mainly focused on cell synthesis, cycle, and metabolism; cell component annotation mainly focused on organelles and internal components; and molecular function annotation mainly focused on binding molecular functions. KEGG pathway enrichment analysis found that it involved the protein tyrosine kinase 6(PTK6) pathway, CD28-dependent phosphoinositide 3-kinase/protein kinase B(PI3K/AKT) pathway, p53 pathway, etc. Therefore, Huangjing Qianshi Decoction can improve the state of prediabetes, and the mechanism may be related to cell cycle and apoptosis, PI3K/AKT pathway, p53 pathway, and other biological pathways regulated by IL-6, NR3C2, and VEGFA.

Mechanism of artesunate on bone destruction in experimental rheumatoid arthritis based on transcriptomics and network pharmacology / 中国中药杂志

The present study investigated the mechanism of artesunate in the treatment of bone destruction in experimental rheumatoid arthritis(RA) based on transcriptomics and network pharmacology. The transcriptome sequencing data of artesunate in the inhibition of osteoclast differentiation were analyzed to obtain differentially expressed genes(DEGs). GraphPad Prism 8 software was used to plot volcano maps and heat maps were plotted through the website of bioinformatics. GeneCards and OMIM were used to collect information on key targets of bone destruction in RA. The DEGs of artesunate in inhibiting osteoclast differentiation and key target genes of bone destruction in RA were intersected by the Venny 2.1.0 platform, and the intersection target genes were analyzed by Gene Ontology(GO)/Kyoto Encyclopedia of Genes and Genomes(KEGG) enrichment. Finally, the receptor activator of nuclear factor-κB(RANKL)-induced osteoclast differentiation model and collagen-induced arthritis(CIA) model were established. Quantitative real time polymerase chain reaction(q-PCR), immunofluorescence, and immunohistochemistry were used to verify the pharmacological effect and molecular mechanism of artesunate in the treatment of bone destruction in RA. In this study, the RANKL-induced osteoclast differentiation model in vitro was established and intervened with artesunate, and transcriptome sequencing data were analyzed to obtain 744 DEGs of artesunate in inhibiting osteoclast differentiation. A total of 1 291 major target genes of bone destruction in RA were obtained from GeneCards and OMIM. The target genes of artesunate in inhibiting osteoclast differentiation and the target genes of bone destruction in RA were intersected to obtain 61 target genes of artesunate against bone destruction in RA. The intersected target genes were analyzed by GO/KEGG enrichment. According to the results previously reported, the cytokine-cytokine receptor interaction signaling pathway was selected for experimental verification. Artesunate intervention in the RANKL-induced osteoclast differentiation model showed that artesunate inhibited CC chemokine receptor 3(CCR3), CC chemokine receptor 1(CCR1) and leukemia inhibitory factor(LIF) mRNA expression in osteoclasts in a dose-dependent manner compared with the RANKL-induced group. Meanwhile, the results of immunofluorescence and immunohistochemistry showed that artesunate could dose-dependently reduce the expression of CCR3 in osteoclasts and joint tissues of the CIA rat model in vitro. This study indicated that artesunate regulated the CCR3 in the cytokine-cytokine receptor interaction signaling pathway in the treatment of bone destruction in RA and provided a new target gene for the treatment of bone destruction in RA.

Spatial transcriptome analysis of long non-coding RNAs reveals tissue specificity and functional roles in cancer / 浙江大学学报（英文版）（B辑：生物医学和生物技术）

Long non-coding RNAs (lncRNAs) play a significant role in maintaining tissue morphology and functions, and their precise regulatory effectiveness is closely related to expression patterns. However, the spatial expression patterns of lncRNAs in humans are poorly characterized. Here, we constructed five comprehensive transcriptomic atlases of human lncRNAs covering thousands of major tissue samples in normal and disease states. The lncRNA transcriptomes exhibited high consistency within the same tissues across resources, and even higher complexity in specialized tissues. Tissue-elevated (TE) lncRNAs were identified in each resource and robust TE lncRNAs were refined by integrative analysis. We detected 1 to 4684 robust TE lncRNAs across tissues; the highest number was in testis tissue, followed by brain tissue. Functional analyses of TE lncRNAs indicated important roles in corresponding tissue-related pathways. Moreover, we found that the expression features of robust TE lncRNAs made them be effective biomarkers to distinguish tissues; TE lncRNAs also tended to be associated with cancer, and exhibited differential expression or were correlated with patient survival. In summary, spatial classification of lncRNAs is the starting point for elucidating the function of lncRNAs in both maintenance of tissue morphology and progress of tissue-constricted diseases.

Identification and expression pattern analysis of RcACA gene family in castor under abiotic stresses / 生物工程学报

Auto-inhibited Ca2+-ATPase (ACA) is one of the Ca2+-ATPase subfamilies that plays an important role in maintaining Ca2+ concentration balance in plant cells. To explore the function and gene expression pattern of the RcACA gene family in castor, bioinformatics analysis was used to identify the members of the RcACA gene family in castor. The basic physical and chemical properties, subcellular location, protein secondary and tertiary structure, conserved domain, conserved motif, gene structure, chromosome location and collinear relationship, as well as the evolutionary characteristics and promoter cis-acting elements were predicted and analyzed. The expression pattern of the RcACA gene under abiotic stress was analyzed by expression (fragments per kilobase of exon model per million mapped fragments, FPKM) in castor transcriptome data. The results showed that 8 RcACA gene family members were identified in castor, acidic proteins located in the plasma membrane. In the secondary structure of all proteins, the α-helix and random coil is more; the RcACA genes were clustered into three categories, and the design of the genes in the same category was similar to the conserved motif. Both of them had four typical domains, RcACA3-RcACA8 had a Ca2+-ATPase N-terminal autoinhibitory domain. The RcACA gene is mostly located on the long arm of the chromosome and has 2 pairs of collinear relationships. There are more light response elements but fewer hormone-induced elements located upstream of the RcACA coding region. Interspecific clustering showed that the evolution of ACA genes among species was conservative. Tissue expression pattern analysis showed that RcACA genes showed apparent tissue expression specificity, and most of the genes showed the highest expression level in male flowers. Expression analysis under abiotic stress showed that RcACA2-RcACA8 were up-regulated under high salt and drought stress, and RcACA1 was up-regulated at 0-24 h under low-temperature stress, indicating that RcACA genes positively responded to abiotic stresses. The above results provide a theoretical basis for exploring the role of the RcACA gene in castor growth, development and stress response.

Comparative analysis of differentially expressed genes for biosynthesis of active ingredients in fruits of different cultivars of Lycium barbarum L. based on transcriptome sequencing / 生物工程学报

To explore the differentially expressed genes (DEGs) related to biosynthesis of active ingredients in wolfberry fruits of different varieties of Lycium barbarum L. and reveal the molecular mechanism of the differences of active ingredients, we utilized Illumina NovaSeq 6000 high-throughput sequencing technology to conduct transcriptome sequencing on the fruits of 'Ningqi No.1' and 'Ningqi No.7' during the green fruit stage, color turning stage and maturity stage. Subsequently, we compared the profiles of related gene expression in the fruits of the two varieties at different development stages. The results showed that a total of 811 818 178 clean reads were obtained, resulting in 121.76 Gb of valid data. There were 2 827, 2 552 and 2 311 DEGs obtained during the green fruit stage, color turning stage and maturity stage of 'Ningqi No. 1' and 'Ningqi No. 7', respectively, among which 2 153, 2 050 and 1 825 genes were annotated in six databases, including gene ontology (GO), Kyoto encyclopedia of genes and genomes (KEGG) and clusters of orthologous groups of proteins (KOG). In GO database, 1 307, 865 and 624 DEGs of green fruit stage, color turning stage and maturity stage were found to be enriched in biological processes, cell components and molecular functions, respectively. In the KEGG database, the DEGs at three developmental stages were mainly concentrated in metabolic pathways, biosynthesis of secondary metabolites and plant-pathogen interaction. In KOG database, 1 775, 1 751 and 1 541 DEGs were annotated at three developmental stages, respectively. Searching the annotated genes against the PubMed database revealed 18, 26 and 24 DEGs related to the synthesis of active ingredients were mined at the green fruit stage, color turning stage and maturity stage, respectively. These genes are involved in carotenoid, flavonoid, terpenoid, alkaloid, vitamin metabolic pathways, etc. Seven DEGs were verified by RT-qPCR, which showed consistent results with transcriptome sequencing. This study provides preliminary evidences for the differences in the content of active ingredients in different Lycium barbarum L. varieties from the transcriptional level. These evidences may facilitate further exploring the key genes for active ingredients biosynthesis in Lycium barbarum L. and analyzing their expression regulation mechanism.

Amino acid metabolism characteristics of Banxia Baizhu Tianma Decoction in realizing drug withdrawal based on transcriptomic analysis / 中国中药杂志

This study aimed to demonstrate the effect of Banxia Baizhu Tianma Decoction(BBTD) on realizing withdrawal of anti-epileptic drugs and explore the relationship between BBTD and the amino acid metabolism by transcriptomic analysis in the rat model of epilepsy induced by lithium chloride-pilocarpine. The rats with epilepsy were divided into a control group(Ctrl), an epilepsy group(Ep), a BBTD & antiepileptic drug integrative group(BADIG), and an antiepileptic drug withdrawal group(ADWG). The Ctrl and Ep were given ultrapure water by gavage for 12 weeks. The BADIG was given BBTD extract and carbamazepine solution by gavage for 12 weeks. The ADWG was given carbamazepine solution and BBTD extract by gavage for the former 6 weeks, and then only given BBTD extract for the latter 6 weeks. The therapeutic effect was evaluated by behavioral observation, electroencephalogram(EEG), and hippocampal neuronal morphological changes. High-throughput sequencing was used to obtain amino acid metabolism-related differen-tial genes in the hippocampus, and the mRNA expression in the hippocampus of each group was verified by real-time quantitative polymerase chain reaction(RT-qPCR). The hub genes were screened out through protein-protein interaction(PPI) network, and Gene Ontology(GO) functional enrichment analysis and Kyoto Encyclopedia of Genes and Genomes(KEGG) pathway enrichment analysis were performed. Two ceRNA networks, namely circRNA-miRNA-mRNA and lncRNA-miRNA-mRNA, were constructed for ADWG vs BADIG. The experimental results showed that compared with those in Ep, rats in ADWG were significantly improved in the behavioral observation, EEG, and hippocampal neuronal impairment. Thirty-four amino acid metabolism-related differential genes were obtained by transcriptomic analysis, and the sequencing results were confirmed by RT-qPCR. Eight hub genes were obtained through PPI network, involving several biological processes, molecular functions, and signal pathways related to amino acid metabolism. Finally, the circRNA-miRNA-mRNA ternary transcription network of 17 circRNA, 5 miRNA, and 2 mRNA, and a lncRNA-miRNA-mRNA ternary network of 10 lncRNA, 5 miRNA, and 2 mRNA were constructed in ADWG vs BADIG. In conclusion, BBTD can effectively achieve the withdrawal of antiepileptic drugs, which may be related to the transcriptomic regulation of amino acid metabolism.

Immune regulation mechanism of Saposhnikoviae Radix polysaccharide based on zebrafish model / 中国中药杂志

The immunomodulatory effect of Saposhnikoviae Radix polysaccharide(SRP) was evaluated based on the zebrafish mo-del, and its mechanism was explored by transcriptome sequencing and real-time fluorescence-based quantitative PCR(RT-qPCR). The immune-compromised model was induced by navelbine in the immunofluorescence-labeled transgenic zebrafish Tg(lyz: DsRed), and the effect of SRP on the density and distribution of macrophages in zebrafish was evaluated. The effect of SRP on the numbers of macrophages and neutrophils in wild-type AB zebrafish was detected by neutral red and Sudan black B staining. The content of NO in zebrafish was detected by DAF-FM DA fluorescence probe. The content of IL-1β and IL-6 in zebrafish was detected by ELISA. The differentially expressed genes(DEGs) of zebrafish in the blank control group, the model group, and the SRP treatment group were analyzed by transcriptome sequencing. The immune regulation mechanism was analyzed by Gene Ontology(GO)and Kyoto Encyclopedia of Genes and Genomes(KEGG)enrichment, and the expression levels of key genes were verified by RT-qPCR. The results showed that SRP could significantly increase the density of immune cells in zebrafish, increase the number of macrophages and neutrophils, and reduce the content of NO, IL-1β, and IL-6 in immune-compromised zebrafish. The results of transcriptome sequencing analysis showed that SRP could affect the expression level of immune-related genes on Toll-like receptor pathway and herpes simplex infection pathway to affect the release of downstream cytokines and interferon, thereby completing the activation process of T cells and playing a role in regulating the immune activity of the body.

Developmental effects of TCIPP and TnBP on zebrafish (Danio rerio) embryos / 中华预防医学杂志

Objective: To investigate the toxicity of tris (2-chloropropyl) phosphate (TCIPP) and tributyl phosphate (TnBP) on the growth and development of zebrafish embryos, as well as to explore the underlying mechanisms at the transcriptional level. Methods: With zebrafish as a model, two hpf zebrafish embryos were exposed to TCIPP and TnBP (0.1, 1, 10, 100, 500, and 1 000 μmol/L) using the semi-static method, and their rates of lethality and hatchability were determined. The transcriptome changes of 120 hpf juvenile zebrafish exposed to environmentally relevant concentrations of 0.1 and 1 μmol/L were measured. Results: The 50% lethal concentrations (LC50) of TCIPP and TnBP for zebrafish embryos were 155.30 and 27.62 μmol/L (96 hpf), 156.5 and 26.05 μmol/L (120 hpf), respectively. The 72 hpf hatching rates of TCIPP (100 μmol/L) and TnBP (10 μmol/L) were (23.33±7.72)% and (91.67±2.97)%, which were significantly decreased compared with the control group (P<0.05). Transcriptome analysis showed that TnBP had more differential genes (DEGs) than TCIPP, with a dose-response relationship. These DEGs were enriched in 32 pathways in total, including those involved in oxidative stress, energy metabolism, lipid metabolism, and nuclear receptor-related pathways, using the IPA pathway analysis. Among them, three enriched pathways overlapped between TCIPP and TnBP, including TR/RXR activation and CAR/RXR activation. Additionally, DEGs were also mapped onto pathways of LXR/RXR activation and oxidative stress for TnBP exposure only. Conclusion: Both TCIPP and TnBP have growth and developmental toxicities in zebrafish embryos, with distinct biomolecular mechanisms, and TnBP has a stronger effect than TCIPP.

Research progress of transcriptomics and proteomics in schizophrenia / 中华预防医学杂志

Schizophrenia is a severe psychiatric disorder with an unclear etiology and various clinical manifestations. The diagnosis and consequent treatment of schizophrenia mainly rely on clinical symptoms. Multiple risk sites associated with schizophrenia have been identified, yet objective indicators have not been found to facilitate clinical diagnosis and treatment of schizophrenia. The development of omics technology provides different perspectives on the etiology of schizophrenia and make the early identification, diagnosis and treatment of the disorder possible. This article summarizes the prevalence of schizophrenia, reviews the research results and shortcomings of transcriptomics and proteomics, as well as the latest achievements and prospects of multi-omics, aiming to reveal the use of omics in SZ, provide more comprehensive biological evidence to reveal the complex pathogenesis of schizophrenia and provide a theoretical basis for the early identification, accurate diagnosis, disease progression control, and prognosis improvement of schizophrenia.

Research progress of transcriptomics and proteomics in schizophrenia / 中华预防医学杂志

Schizophrenia is a severe psychiatric disorder with an unclear etiology and various clinical manifestations. The diagnosis and consequent treatment of schizophrenia mainly rely on clinical symptoms. Multiple risk sites associated with schizophrenia have been identified, yet objective indicators have not been found to facilitate clinical diagnosis and treatment of schizophrenia. The development of omics technology provides different perspectives on the etiology of schizophrenia and make the early identification, diagnosis and treatment of the disorder possible. This article summarizes the prevalence of schizophrenia, reviews the research results and shortcomings of transcriptomics and proteomics, as well as the latest achievements and prospects of multi-omics, aiming to reveal the use of omics in SZ, provide more comprehensive biological evidence to reveal the complex pathogenesis of schizophrenia and provide a theoretical basis for the early identification, accurate diagnosis, disease progression control, and prognosis improvement of schizophrenia.

Identification of core pathogenic genes and pathways in elderly osteoporosis based on bioinformatics analysis / 中华预防医学杂志

Objective: Using bioinformatics methods to analyze the core pathogenic genes and related pathways in elderly osteoporosis. Methods: From November 2020 and August 2021, eight elderly osteoporosis patients who received treatment and five healthy participants who underwent physical examinations in Beijing Jishuitan Hospital were selected as subjects. The expression level of RNA in the peripheral blood of eight elderly osteoporosis patients and five healthy participants was collected for high-throughput transcriptome sequencing and analysis. The gene ontology (GO) analysis Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis was performed for the differentially expressed genes (DEGs). The protein-protein interaction (PPI) network was constructed using the STRING website and Cytoscape software, and the most significant modules and hub genes were screened out. Results: Among the eight elderly osteoporosis patients, there were seven females and one male, with an average age of 72.4 years (SD=4.2). Among the five healthy participants, there were four females and one male, with an average age of 68.2 years (SD=5.7). A total of 1 635 DEGs (847 up-regulated and 788 down-regulated) were identified. GO analysis revealed that the molecular functions of DEGs were mainly enriched in structural constituents of the ribosome, protein dimerization activity, and cellular components were mainly enriched in the nucleosome, DNA packaging complex, cytosolic part, protein-DNA complex and the cytosolic ribosome. KEGG pathway analysis showed that DEGs were mainly enriched in systemic lupus erythematosus and ribosome. Gene UBA52, UBB, RPS27A, RPS15, RPS12, RPL13A, RPL23A, RPL10A, RPS25 and RPS6 were selected and seven of them could encode ribosome proteins. Conclusion: The pathogenesis of elderly osteoporosis may be associated with ribosome-related genes and pathways.

EHHADH is a key gene in fatty acid metabolism pathways in hepatocellular carcinoma: a transcriptomic analysis / 南方医科大学学报

OBJECTIVE@#To explore the driving gene of hepatocellular carcinoma (HCC) occurrence and progression and its potential as new therapeutic target of HCC.@*METHODS@#The transcriptome and genomic data of 858 HCC tissues and 493 adjacent tissues were obtained from TCGA, GEO, and ICGC databases. Gene Set Enrichment Analysis (GSEA) identified EHHADH (encoding enoyl-CoA hydratase/L-3-hydroxyacyl-CoA dehydrogenase) as the hub gene in the significantly enriched differential pathways in HCC. The downregulation of EHHADH expression at the transcriptome level was found to correlate with TP53 mutation based on analysis of the TCGA- HCC dataset, and the mechanism by which TP53 mutation caused EHHADH downregulation was explored through correlation analysis. Analysis of the data from the Metascape database suggested that EHHADH was strongly correlated with the ferroptosis signaling pathway in HCC progression, and to verify this result, immunohistochemical staining was used to examine EHHADH expression in 30 HCC tissues and paired adjacent tissues.@*RESULTS@#All the 3 HCC datasets showed signficnatly lowered EHHADH expression in HCC tissues as compared with the adjacent tissues (P < 0.05) with a close correlation with the degree of hepatocyte de-differentiation (P < 0.01). The somatic landscape of HCC cohort in TCGA dataset showed that HCC patients had the highest genomic TP53 mutation rate. The transcriptomic level of PPARGC1A, the upstream gene of EHHADH, was significantly downregulated in HCC patients with TP53 mutation as compared with those without the mutation (P < 0.05), and was significantly correlated with EHHADH expression level. GO and KEGG enrichment analyses showed that EHHADH expression was significantly correlated with abnormal fatty acid metabolism in HCC. The immunohistochemical results showd that the expression level of EHHADH in HCC tissues was down-regulated, and its expression level was related to the degree of hepatocytes de-differentiation and the process of ferroptosis.@*CONCLUSION@#TP53 mutations may induce abnormal expression of PPARGC1A to cause downregulation of EHHADH expression in HCC. The low expression of EHHADH is closely associated with aggravation of de-differentiation and ferroptosis escape in HCC tissues, suggesting the potential of EHHADH as a therapeutic target for HCC.

Single-cell RNA sequencing reveals the transcriptomic landscape of kidneys in patients with ischemic acute kidney injury / 中华医学杂志(英文版)

BACKGROUND@#Ischemic acute kidney injury (AKI) is a common syndrome associated with considerable mortality and healthcare costs. Up to now, the underlying pathogenesis of ischemic AKI remains incompletely understood, and specific strategies for early diagnosis and treatment of ischemic AKI are still lacking. Here, this study aimed to define the transcriptomic landscape of AKI patients through single-cell RNA sequencing (scRNA-seq) analysis in kidneys.@*METHODS@#In this study, scRNA-seq technology was applied to kidneys from two ischemic AKI patients, and three human public scRNA-seq datasets were collected as controls. Differentially expressed genes (DEGs) and cell clusters of kidneys were determined. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis, as well as the ligand-receptor interaction between cells, were performed. We also validated several DEGs expression in kidneys from human ischemic AKI and ischemia/reperfusion (I/R) injury induced AKI mice through immunohistochemistry staining.@*RESULTS@#15 distinct cell clusters were determined in kidney from subjects of ischemic AKI and control. The injured proximal tubules (PT) displayed a proapoptotic and proinflammatory phenotype. PT cells of ischemic AKI had up-regulation of novel pro-apoptotic genes including USP47 , RASSF4 , EBAG9 , IER3 , SASH1 , SEPTIN7 , and NUB1 , which have not been reported in ischemic AKI previously. Several hub genes were validated in kidneys from human AKI and renal I/R injury mice, respectively. Furthermore, PT highly expressed DEGs enriched in endoplasmic reticulum stress, autophagy, and retinoic acid-inducible gene I (RIG-I) signaling. DEGs overexpressed in other tubular cells were primarily enriched in nucleotide-binding and oligomerization domain (NOD)-like receptor signaling, estrogen signaling, interleukin (IL)-12 signaling, and IL-17 signaling. Overexpressed genes in kidney-resident immune cells including macrophages, natural killer T (NKT) cells, monocytes, and dendritic cells were associated with leukocyte activation, chemotaxis, cell adhesion, and complement activation. In addition, the ligand-receptor interactions analysis revealed prominent communications between macrophages and monocytes with other cells in the process of ischemic AKI.@*CONCLUSION@#Together, this study reveals distinct cell-specific transcriptomic atlas of kidney in ischemic AKI patients, altered signaling pathways, and potential cell-cell crosstalk in the development of AKI. These data reveal new insights into the pathogenesis and potential therapeutic strategies in ischemic AKI.

Characterization of candidate factors associated with the metastasis and progression of high-grade serous ovarian cancer / 中华医学杂志(英文版)

BACKGROUND@#High-grade serous ovarian cancer (HGSOC) is the biggest cause of gynecological cancer-related mortality because of its extremely metastatic nature. This study aimed to explore and evaluate the characteristics of candidate factors associated with the metastasis and progression of HGSOC.@*METHODS@#Transcriptomic data of HGSOC patients' samples collected from primary tumors and matched omental metastatic tumors were obtained from three independent studies in the National Center for Biotechnology Information (NCBI) Gene Expression Omnibus (GEO) database. Differentially expressed genes (DEGs) were selected to evaluate the effects on the prognosis and progression of ovarian cancer using data from The Cancer Genome Atlas (TCGA) database. Hub genes' immune landscapes were estimated by the Tumor Immune Estimation Resource (TIMER) database. Finally, using 25 HGSOC patients' cancer tissues and 10 normal fallopian tube tissues, immunohistochemistry (IHC) was performed to quantify the expression levels of hub genes associated with International Federation of Gynecology and Obstetrics (FIGO) stages.@*RESULTS@#Fourteen DEGs, ADIPOQ , ALPK2 , BARX1 , CD37 , CNR2 , COL5A3 , FABP4 , FAP , GPR68 , ITGBL1 , MOXD1 , PODNL1 , SFRP2 , and TRAF3IP3 , were upregulated in metastatic tumors in every database while CADPS , GATA4 , STAR , and TSPAN8 were downregulated. ALPK2 , FAP , SFRP2 , GATA4 , STAR , and TSPAN8 were selected as hub genes significantly associated with survival and recurrence. All hub genes were correlated with tumor microenvironment infiltration, especially cancer-associated fibroblasts and natural killer (NK) cells. Furthermore, the expression of FAP and SFRP2 was positively correlated with the International Federation of Gynecology and Obstetrics (FIGO) stage, and their increased protein expression levels in metastatic samples compared with primary tumor samples and normal tissues were confirmed by IHC ( P = 0.0002 and P = 0.0001, respectively).@*CONCLUSIONS@#This study describes screening for DEGs in HGSOC primary tumors and matched metastasis tumors using integrated bioinformatics analyses. We identified six hub genes that were correlated with the progression of HGSOC, particularly FAP and SFRP2 , which might provide effective targets to predict prognosis and provide novel insights into individual therapeutic strategies for HGSOC.

An atlas of immune cell transcriptomes in human immunodeficiency virus-infected immunological non-responders identified marker genes that control viral replication / 中华医学杂志(英文版)

BACKGROUND@#Previous studies have examined the bulk transcriptome of peripheral blood immune cells in acquired immunodeficiency syndrome patients experiencing immunological non-responsiveness. This study aimed to investigate the characteristics of specific immune cell subtypes in acquired immunodeficiency syndrome patients who exhibit immunological non-responsiveness.@*METHODS@#A single-cell transcriptome sequencing of peripheral blood mononuclear cells obtained from both immunological responders (IRs) (CD4 + T-cell count >500) and immunological non-responders (INRs) (CD4 + T-cell count <300) was conducted. The transcriptomic profiles were used to identify distinct cell subpopulations, marker genes, and differentially expressed genes aiming to uncover potential genetic factors associated with immunological non-responsiveness.@*RESULTS@#Among the cellular subpopulations analyzed, the ratios of monocytes, CD16 + monocytes, and exhausted B cells demonstrated the most substantial differences between INRs and IRs, with fold changes of 39.79, 11.08, and 2.71, respectively. In contrast, the CD4 + T cell ratio was significantly decreased (0.39-fold change) in INRs compared with that in IRs. Similarly, the ratios of natural killer cells and terminal effector CD8 + T cells were also lower (0.37-fold and 0.27-fold, respectively) in the INRs group. In addition to several well-characterized immune cell-specific markers, we identified a set of 181 marker genes that were enriched in biological pathways associated with human immunodeficiency virus (HIV) replication. Notably, ISG15 , IFITM3 , PLSCR1 , HLA-DQB1 , CCL3L1 , and DDX5 , which have been demonstrated to influence HIV replication through their interaction with viral proteins, emerged as significant monocyte marker genes. Furthermore, the differentially expressed genes in natural killer cells were also enriched in biological pathways associated with HIV replication.@*CONCLUSIONS@#We generated an atlas of immune cell transcriptomes in HIV-infected IRs and INRs. Host genes associated with HIV replication were identified as markers of, and were found to be differentially expressed in, different types of immune cells.