Alterações sistêmicas na assinatura do ciclo celular de pacientes com COVID-19 / Systemic changes in the cell cycle signature of patients with COVID-19

Com base nas perturbações fosfoproteômicas de moléculas associadas ao ciclo celular em células infectadas pelo coronavírus causador da síndrome respiratória aguda grave (SARSCoV)-2, a hipótese de inibidores do ciclo celular como uma terapia potencial para a doença de coronavírus 2019 (COVID-19) foi proposta. No entanto, o cenário das alterações do ciclo celular em COVID-19 permanece inexplorado. Aqui, realizamos uma análise integrativa de sistemas imunológicos de proteoma publicamente disponível (espectrometria de massa) e dados de transcriptoma (sequenciamento de RNA em massa e de célula única [scRNAseq]), com o objetivo de caracterizar mudanças globais na assinatura do ciclo celular de pacientes com COVID-19. Além de módulos de co-expressão de genes significativos enriquecidos associados ao ciclo celular, encontramos uma rede interconectada de proteínas diferencialmente expressas associadas ao ciclo celular (DEPs) e genes (DEGs) integrando dados moleculares de 1.480 indivíduos (974 pacientes infectados por SARS-CoV-2 e 506 controles [controles saudáveis ou indivíduos com outras doenças respiratórias]). Entre esses DEPs e DEGs estão várias ciclinas (CCNs), ciclo de divisão celular (CDCs), quinases dependentes de ciclinas (CDKs) e proteínas de manutenção de minicromossomos (MCMs). Embora os pacientes com COVID-19 compartilhem parcialmente o padrão de expressão de algumas moléculas associadas ao ciclo celular com outras doenças respiratórias, eles exibiram uma expressão significativamente maior de moléculas associadas ao ciclo celular relacionadas à gravidade da doença. Notavelmente, a assinatura do ciclo celular predominou nos leucócitos do sangue dos pacientes, mas não nas vias aéreas superiores. Os dados de scRNAseq de 229 indivíduos (159 pacientes com COVID- 19 e 70 controles) revelaram que as alterações das assinaturas do ciclo celular predominam nas células B, T e NK. Esses resultados fornecem uma compreensão global única das alterações nas moléculas associadas ao ciclo celular em pacientes com COVID-19, sugerindo novas vias putativas para intervenção terapêutica

Based on phosphoproteomics perturbations of cell cycle-associated molecules in severe acute respiratory syndrome coronavirus (SARS-CoV)-2-infected cells, the hypothesis of cell cycle inhibitors as a potential therapy for Coronavirus disease 2019 (COVID-19) has been proposed. However, the landscape of cell cycle alterations in COVID-19 remains mostly unexplored. Here, we performed an integrative systems immunology analysis of publicly available proteome (mass spectrometry) and transcriptome data (bulk and single-cell RNA sequencing [scRNAseq]), aiming to characterize global changes in the cell cycle signature of COVID-19 patients. Beyond significant enriched cell cycle-associated gene co-expression modules, we found an interconnected network of cell cycle-associated differentially expressed proteins (DEPs) and genes (DEGs) by integrating molecular data of 1,480 individuals (974 SARS-CoV- 2 infected patients and 506 controls [either healthy controls or individuals with other respiratory illness]). Among these DEPs and DEGs are several cyclins (CCNs), cell division cycle (CDCs), cyclin-dependent kinases (CDKs), and mini-chromosome maintenance proteins (MCMs). Although COVID-19 patients partially shared the expression pattern of some cell cycleassociated molecules with other respiratory illnesses, they exhibited a significantly higher expression of cell cycle-associated molecules associated with disease severity. Notably, the cell cycle signature predominated in the patients blood leukocytes but not in the upper airways. The scRNAseq data from 229 individuals (159 COVID-19 patients and 70 controls) revealed that the alterations of cell cycle signatures predominate in B, T, and NK cells. These results provide a unique global comprehension of the alterations in cell cycle-associated molecules in COVID-19 patients, suggesting new putative pathways for therapeutic intervention

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.

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.

Mechanism of transforming growth factor- β1 induce renal fibrosis based on transcriptome sequencing analysis / 浙江大学学报·医学版

OBJECTIVES@#To explore the mechanism of transforming growth factor-β1 (TGF-β1) induce renal fibrosis.@*METHODS@#Renal fibroblast NRK-49F cells treated with and without TGF-β1 were subjected to RNA-seq analysis. DESeq2 was used for analysis. Differentially expressed genes were screened with the criteria of false discovery rate<0.05 and l o g 2 F C >1. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were performed for differentially expressed genes. Genes encoding transcription factors were further screened for differential expression genes. Then, the expression of these genes during renal fibrosis was verified using unilateral ureteral obstruction (UUO)-induced mouse renal fibrosis model and a public gene expression dataset (GSE104954).@*RESULTS@#After TGF-β1 treatment for 6, 12 and 24 h, 552, 1209 and 1028 differentially expressed genes were identified, respectively. GO analysis indicated that these genes were significantly enriched in development, cell death, and cell migration. KEGG pathway analysis showed that in the early stage of TGF-β1 induction (TGF-β1 treatment for 6 h), the changes in Hippo, TGF-β and Wnt signaling pathways were observed, while in the late stage of TGF-β1 induction (TGF-β1 treatment for 24 h), the changes of extracellular matrix-receptor interaction, focal adhesion and adherens junction were mainly enriched. Among the 291 up-regulated differentially expressed genes treated with TGF-β1 for 6 h, 13 genes (Snai1, Irf8, Bhlhe40, Junb, Arid5a, Vdr, Lef1, Ahr, Foxo1, Myc, Tcf7, Foxc2, Glis1) encoded transcription factors. Validation in a cell model showed that TGF-β1 induced expression of 9 transcription factors (encoded by Snai1, Irf8, Bhlhe40, Junb, Arid5a, Vdr, Lef1, Myc, Tcf7), while the expression levels of the other 4 genes did not significantly change after TGF-β1 treatment. Validation results in UUO-induced mouse renal fibrosis model showed that Snai1, Irf8, Bhlhe40, Junb, Arid5a, Myc and Tcf7 were up-regulated after UUO, Vdr was down-regulated and there was no significant change in Lef1. Validation based on the GSE104954 dataset showed that IRF8 was significantly overexpressed in the renal tubulointerstitium of patients with diabetic nephropathy or IgA nephropathy, MYC was highly expressed in diabetic nephropathy, and the expressions of the other 7 genes were not significantly different compared with the control group.@*CONCLUSIONS@#TGF-β1 induces differentially expressed genes in renal fibroblasts, among which Irf8 and Myc were identified as potential targets of chronic kidney disease and renal fibrosis.

Study on the mechanism of haitongpi prescription extract ointment in the treatment of knee osteoarthritis based on transcriptome / 中国骨伤

OBJECTIVE@#To explore the mechanism of Haitongpi Prescription extract in the treatment of knee osteoarthritis based on transcriptome.@*METHODS@#Total of 12 SPF grade rats were divided into control group(group C), model group(group M), and Haitongpi prescription group(group HP). The knee osteoarthritis rat model was established using the Panicker method for group M and group HP, and group HP was intervened by local topical application of Haitongpi Prescription extract for 4 weeks. Total RNA from mouse knee cartilage was extracted and three sets of differential genes were obtained through sequencing.Differential genes were prediction and analysis through GO function and KEGG pathway enrichment analysis.@*RESULTS@#A total of 109 differentially expressed genes were identified in Group C versus Group M, while 118 differentially expressed genes were identified in Group M versus Group HP, resulting in a total of 28 genes. GO functional enrichment analysis showed that the mechanism of HP extract in treating knee osteoarthritis mainly involved immunoglobulin mediated immune response, immunoglobulin complexes, and antigen binding; KEGG pathway enrichment analysis showed correlation with tumor necrosis factor (TNF) signaling pathway, interleukin 17(IL-17) signaling pathway, and estrogen signaling pathway.@*CONCLUSION@#HP extract can exert therapeutic effects on knee osteoarthritis through mechanisms such as immunoglobulin mediated immune response, immunoglobulin complexes, and antigen binding, as well as signaling pathways such as TNF signaling pathway, IL-17 signaling pathway, and estrogen signaling pathway.

Imputation method for dropout in single-cell transcriptome data / 生物医学工程学杂志

Single-cell transcriptome sequencing (scRNA-seq) can resolve the expression characteristics of cells in tissues with single-cell precision, enabling researchers to quantify cellular heterogeneity within populations with higher resolution, revealing potentially heterogeneous cell populations and the dynamics of complex tissues. However, the presence of a large number of technical zeros in scRNA-seq data will have an impact on downstream analysis of cell clustering, differential genes, cell annotation, and pseudotime, hindering the discovery of meaningful biological signals. The main idea to solve this problem is to make use of the potential correlation between cells and genes, and to impute the technical zeros through the observed data. Based on this, this paper reviewed the basic methods of imputing technical zeros in the scRNA-seq data and discussed the advantages and disadvantages of the existing methods. Finally, recommendations and perspectives on the use and development of the method were provided.

Screening and expression analysis of transcription factors involved in genuineness of Codonopsis pilosula in Shanxi / 中国中药杂志

This study aims to mine the transcription factors that affect the genuineness of Codonopsis pilosula in Shanxi based on the transcriptome data of C. pilosula samples collected from Shanxi and Gansu, and then analyze the gene expression patterns, which will provide a theoretical basis for the molecular assisted breeding of C. pilosula. Gene ontology(GO) functional annotation, conserved motif prediction, and gene expression pattern analysis were performed for the differential transcription factors predicted based on the transcriptome data of C. pilosula from different habitats. A total of 61 differentially expressed genes(DEGs) were screened out from the transcriptome data. Most of the DEGs belonged to AP2/ERF-ERF family, with the conserved motif of [2X]-[LG]-[3X]-T-[3X]-[AARAYDRAA]-[3X]-[RG]-[2X]-A-[2X]-[NFP]. Forty-three of the DEGs showed significantly higher gene expression in C. pilosula samples from Shanxi than in the samples from Gansu, including 11 genes in the AP2/ERF-ERF family, 5 genes in the NAC fa-mily, 1 gene in the bHLH family, and 2 genes in the RWP-RK family, while 18 transcription factors showed higher expression levels in the samples from Gansu. GO annotation predicted that most of the DEGs were enriched in GO terms related to transcriptional binding activity(103), metabolic process(26), and stress response(23). The expression of transcription factor genes, CpNAC92, CpNAC100, CpbHLH128, and CpRAP2-7 was higher in the samples from Shanxi and in the roots of C. pilosula. CpNAC92, CpbHLH128, and CpRAP2-7 responded to the low temperature, temperature difference, and iron stresses, while CpNAC100 only responded to low temperature and iron stresses. The screening and expression analysis of the specific transcription factors CpNAC92, CpNAC100, CpbHLH128, and CpRAP2-7 in C. pilosula in Shanxi laid a theoretical foundation for further research on the mechanism of genuineness formation of C. pilosula.

Transcriptional regulation mechanism of differential accumulation of flavonoids in leaves and roots of Sarcandra glabra based on metabonomics and transcriptomics / 中国中药杂志

This study aims to explore the molecular regulation mechanism of the differential accumulation of flavonoids in the leaves and roots of Sarcandra glabra. Liquid chromatography-mass spectrometry(LC-MS) and high-throughput transcriptome sequencing(RNA-seq) were employed to screen out the flavonoid-related differential metabolites and differentially expressed genes(DEGs) encoding key metabolic enzymes. Eight DEGs were randomly selected for qRT-PCR verification. The results showed that a total of 37 flavonoid-related differential metabolites between the leaves and roots of S. glabra were obtained, including pinocembrin, phlorizin, na-ringenin, kaempferol, leucocyanidin, and 5-O-caffeoylshikimic acid. The transcriptome analysis predicted 36 DEGs associated with flavonoids in the leaves and roots of S. glabra, including 2 genes in the PAL pathway, 3 genes in the 4CL pathway, 2 genes in the CHS pathway, 4 genes in the CHI pathway, 2 genes in the FLS pathway, 1 gene in the DFR pathway, 1 gene in the CYP73A pathway, 1 gene in the CYP75B1 pathway, 3 genes in the PGT1 pathway, 6 genes in the HCT pathway, 2 genes in the C3'H pathway, 1 gene in the CCOAOMT pathway, 1 gene in the ANR pathway, 1 gene in the LAR pathway, 2 genes in the 3AT pathway, 1 gene in the BZ1 pathway, 2 genes in the IFTM7 pathway, and 1 gene in the CYP81E9 pathway. Six transcription factors, including C2H2, bHLH, and bZIP, were involved in regulating the differential accumulation of flavonoids in the leaves and roots of S. glabra. The qRT-PCR results showed that the up-or down-regulated expression of the 8 randomly selected enzyme genes involved in flavonoid synthesis in the leaves and roots of S. glabra was consistent with the transcriptome sequencing results. This study preliminarily analyzed the transcriptional regulation mechanism of differential accumulation of flavonoids in the leaves and roots of S. glabra, laying a foundation for further elucidating the regulatory effects of key enzyme genes and corresponding transcription factors on the accumulation of flavonoids in S. glabra.

Selection and validation of reference genes for quantitative real-time PCR analysis in Paeonia veitchii / 中国中药杂志

Paeonia veitchii and P. lactiflora are both original plants of the famous Chinese medicinal drug Paeoniae Radix Rubra in the Chinese Pharmacopoeia. They have important medicinal value and great potential in the flower market. The selection of stable and reliable reference genes is a necessary prerequisite for molecular research on P. veitchii. In this study, two reference genes, Actin and GAPDH, were selected as candidate genes from the transcriptome data of P. veitchii. The expression levels of the two candidate genes in different tissues(phloem, xylem, stem, leaf, petiole, and ovary) and different growth stages(bud stage, flowering stage, and dormant stage) of P. veitchii were detected using real-time fluorescence quantitative technology(qRT-PCR). Then, the stability of the expression of the two reference genes was comprehensively analyzed using geNorm, NormFinder, BestKeeper, ΔCT, and RefFinder. The results showed that the expression patterns of Actin and GAPDH were stable in different tissues and growth stages of P. veitchii. Furthermore, the expression levels of eight genes(Pv-TPS01, Pv-TPS02, Pv-CYP01, Pv-CYP02, Pv-CYP03, Pv-BAHD01, Pv-UGT01, and Pv-UGT02) in different tissues were further detected based on the transcriptome data of P. veitchii. The results showed that when Actin and GAPDH were used as reference genes, the expression trends of the eight genes in different tissues of P. veitchii were consistent, validating the reliability of Actin and GAPDH as reference genes for P. veitchii. In conclusion, this study finds that Actin and GAPDH can be used as reference genes for studying gene expression levels in different tissues and growth stages of P. veitchii.

Exploration of transcriptome SSR markers and its application in genetic diversity assessment of Asarum sieboldii / 中国中药杂志

To explore the genetic diversity of Asarum sieboldii this study developed SSR markers based on transcriptome sequencing results and five populations of A.sieboldii from different regions were used as samples for genetic diversity assessment using software such as GenALEx 6.5, NTSYS 2.1, and Structure 2.3.4. The results showed that 16 SSR markers with high polymorphism and good repeatability were selected from the A.sieboldii transcriptome. Primers designed based on the flanking sequences of these markers successfully amplified 56 polymorphic fragments from 150 individual samples of the five A.sieboldii populations. On average, each primer amplified 3.5 polymorphic fragments, ranging from 2 to 8. The mean values of expected heterozygosity(H_e), Shannon's diversity index(I), Nei's gene diversity index(H), and the polymorphic information content(PIC) were 0.172, 0.281, 0.429, and 0.382, respectively. The mean population differentiation coefficient(F_(ST)) was 0.588, consistent with the analysis of molecular variance(AMOVA) results, which indicated greater genetic variation among A.sieboldii populations(69%) than that within populations(31%). The percentage of polymorphic loci(PPL) ranged from highest to lowest as SNJ>LN>SY>SZ>TB. Principal coordinate analysis(PCoA) and UPGMA clustering analysis further revealed genetic clustering of A.sieboldii individuals based on their geographical distribution, consistent with the results of the structure clustering analysis. In summary, the SSR markers developed from the transcriptome effectively assessed the genetic differentiation and population structure of natural A.sieboldii populations, revealing a relatively low genetic diversity in A.sieboldii, with genetic variation primarily observed at the population level and a correlation between population differentiation and geographic distance.

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.

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.

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.

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.

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.

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.

Single-nucleus transcriptomics reveals a gatekeeper role for FOXP1 in primate cardiac aging

Aging poses a major risk factor for cardiovascular diseases, the leading cause of death in the aged population. However, the cell type-specific changes underlying cardiac aging are far from being clear. Here, we performed single-nucleus RNA-sequencing analysis of left ventricles from young and aged cynomolgus monkeys to define cell composition changes and transcriptomic alterations across different cell types associated with age. We found that aged cardiomyocytes underwent a dramatic loss in cell numbers and profound fluctuations in transcriptional profiles. Via transcription regulatory network analysis, we identified FOXP1, a core transcription factor in organ development, as a key downregulated factor in aged cardiomyocytes, concomitant with the dysregulation of FOXP1 target genes associated with heart function and cardiac diseases. Consistently, the deficiency of FOXP1 led to hypertrophic and senescent phenotypes in human embryonic stem cell-derived cardiomyocytes. Altogether, our findings depict the cellular and molecular landscape of ventricular aging at the single-cell resolution, and identify drivers for primate cardiac aging and potential targets for intervention against cardiac aging and associated diseases.

Single-cell transcriptomic atlas of mouse cochlear aging

Progressive functional deterioration in the cochlea is associated with age-related hearing loss (ARHL). However, the cellular and molecular basis underlying cochlear aging remains largely unknown. Here, we established a dynamic single-cell transcriptomic landscape of mouse cochlear aging, in which we characterized aging-associated transcriptomic changes in 27 different cochlear cell types across five different time points. Overall, our analysis pinpoints loss of proteostasis and elevated apoptosis as the hallmark features of cochlear aging, highlights unexpected age-related transcriptional fluctuations in intermediate cells localized in the stria vascularis (SV) and demonstrates that upregulation of endoplasmic reticulum (ER) chaperon protein HSP90AA1 mitigates ER stress-induced damages associated with aging. Our work suggests that targeting unfolded protein response pathways may help alleviate aging-related SV atrophy and hence delay the progression of ARHL.

Multiplexed single-cell transcriptome analysis reveals molecular characteristics of monkey pluripotent stem cell lines / 浙江大学学报（英文版）（B辑：生物医学和生物技术）

Efforts have been made to establish various human pluripotent stem cell lines. However, such methods have not yet been duplicated in non-human primate cells. Here, we introduce a multiplexed single-cell sequencing technique to profile the molecular features of monkey pluripotent stem cells in published culture conditions. The results demonstrate suboptimized maintenance of pluripotency and show that the selected signaling pathways for resetting human stem cells can also be interpreted for establishing monkey cell lines. Overall, this work legitimates the translation of novel human cell line culture conditions to monkey cells and provides guidance for exploring chemical cocktails for monkey stem cell line derivation.

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.