The pattern of alternative splicing and DNA methylation alteration and their interaction in linseed (Linum usitatissimum L) response to repeated drought stresses

BACKGROUND: Drought stress has significantly hampered agricultural productivity worldwide and can also result in modifications to DNA methylation levels. However, the dynamics of DNA methylation and its association with the changes in gene transcription and alternative splicing (AS) under drought stress are unknown in linseed, which is frequently cultivated in arid and semiarid regions. RESULTS: We analysed AS events and DNA methylation patterns in drought-tolerant (Z141) and drought-sensitive (NY-17) linseed under drought stress (DS) and repeated drought stress (RD) treatments. We found that the number of intron-retention (IR) and alternative 3' splice site (Alt3'SS) events were significantly higher in Z141 and NY-17 under drought stress. We found that the linseed response to the DS treatment was mainly regulated by transcription, while the response to the RD treatment was coregulated by transcription and AS. Whole genome-wide DNA methylation analysis revealed that drought stress caused an increase in the overall methylation level of linseed. Although we did not observe any correlation between differentially methylated genes (DMGs) and differentially spliced genes (DSGs) in this study, we found that the DSGs whose gene body region was hypermethylated in Z141 and hypomethylated in NY-17 were enriched in abiotic stress response Gene Ontology (GO) terms. This finding implies that gene body methylation plays an important role in AS regulation in some specific genes. CONCLUSION: Our study is the first comprehensive genome-wide analysis of the relationship between linseed methylation changes and AS under drought and repeated drought stress. Our study revealed different interaction patterns between differentially expressed genes (DEGs) and DSGs under DS and RD treatments and differences between methylation and AS regulation in drought-tolerant and drought-sensitive linseed varieties. The findings will probably be of interest in the future. Our results provide interesting insights into the association between gene expression, AS, and DNA methylation in linseed under drought stress. Differences in these associations may account for the differences in linseed drought tolerance.

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

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.

A single-nucleus transcriptomic atlas of primate testicular aging reveals exhaustion of the spermatogonial stem cell reservoir and loss of Sertoli cell homeostasis

The testis is pivotal for male reproduction, and its progressive functional decline in aging is associated with infertility. However, the regulatory mechanism underlying primate testicular aging remains largely elusive. Here, we resolve the aging-related cellular and molecular alterations of primate testicular aging by establishing a single-nucleus transcriptomic atlas. Gene-expression patterns along the spermatogenesis trajectory revealed molecular programs associated with attrition of spermatogonial stem cell reservoir, disturbed meiosis and impaired spermiogenesis along the sequential continuum. Remarkably, Sertoli cell was identified as the cell type most susceptible to aging, given its deeply perturbed age-associated transcriptional profiles. Concomitantly, downregulation of the transcription factor Wilms' Tumor 1 (WT1), essential for Sertoli cell homeostasis, was associated with accelerated cellular senescence, disrupted tight junctions, and a compromised cell identity signature, which altogether may help create a hostile microenvironment for spermatogenesis. Collectively, our study depicts in-depth transcriptomic traits of non-human primate (NHP) testicular aging at single-cell resolution, providing potential diagnostic biomarkers and targets for therapeutic interventions against testicular aging and age-related male reproductive diseases.

Single-cell transcriptomic analysis uncovers the origin and intratumoral heterogeneity of parotid pleomorphic adenoma / 国际口腔科学杂志·英文版

Pleomorphic adenoma (PA) is the most common benign tumour in the salivary gland and has high morphological complexity. However, the origin and intratumoral heterogeneity of PA are largely unknown. Here, we constructed a comprehensive atlas of PA at single-cell resolution and showed that PA exhibited five tumour subpopulations, three recapitulating the epithelial states of the normal parotid gland, and two PA-specific epithelial cell (PASE) populations unique to tumours. Then, six subgroups of PASE cells were identified, which varied in epithelium, bone, immune, metabolism, stemness and cell cycle signatures. Moreover, we revealed that CD36+ myoepithelial cells were the tumour-initiating cells (TICs) in PA, and were dominated by the PI3K-AKT pathway. Targeting the PI3K-AKT pathway significantly inhibited CD36+ myoepithelial cell-derived tumour spheres and the growth of PA organoids. Our results provide new insights into the diversity and origin of PA, offering an important clinical implication for targeting the PI3K-AKT signalling pathway in PA treatment.

Deciphering the dynamic characteristics of non-neuronal cells in dorsal root ganglion of rat at different developmental stage based on single cell transcriptome data / 生物工程学报

Dorsal root ganglia (DRG) is an essential part of the peripheral nervous system and the hub of the peripheral sensory afferent. The dynamic changes of neuronal cells and their gene expression during the development of dorsal root ganglion have been studied through single-cell RNAseq analysis, while the dynamic changes of non-neuronal cells have not been systematically studied. Using single cell RNA sequencing technology, we conducted a research on the non-neuronal cells in the dorsal root ganglia of rats at different developmental stage. In this study, primary cell suspension was obtained from using the dorsal root ganglions (DRGs, L4-L5) of ten 7-day-old rats and three 3-month-old rats. The 10×Genomics platform was used for single cell dissociation and RNA sequencing. Twenty cell subsets were acquired through cluster dimension reduction analysis, and the marker genes of different types of cells in DRG were identified according to previous researches about DRG single cell transcriptome sequencing. In order to find out the non-neuronal cell subsets with significant differences at different development stage, the cells were classified into different cell types according to markers collected from previous researches. We performed pseudotime analysis of 4 types Schwann cells. It was found that subtype Ⅱ Schwann cells emerged firstly, and then were subtype Ⅲ Schwann cells and subtype Ⅳ Schwann cells, while subtype Ⅰ Schwann cells existed during the whole development procedure. Pseudotime analysis indicated the essential genes influencing cell fate of different subtypes of Schwann cell in DRG, such as Ntrk2 and Pmp2, which affected cell fate of Schwann cells during the development period. GO analysis of differential expressed genes showed that the up-regulated genes, such as Cst3 and Spp1, were closely related to biological process of tissue homeostasis and multi-multicellular organism process. The down regulated key genes, such as Col3a1 and Col4a1, had close relationship with the progress of extracellular structure organization and negative regulation of cell adhesion. This suggested that the expression of genes enhancing cell homestasis increased, while the expression of related genes regulating ECM-receptor interaction pathway decreased during the development. The discovery provided valuable information and brand-new perspectives for the study on the physical and developmental mechanism of Schwann cell as well as the non-neuronal cell changes in DRG at different developmental stage. The differential gene expression results provided crucial references for the mechanism of somatosensory maturation during development.

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.

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.

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.

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.

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.

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.

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.

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.

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.