Project description:The T cell-mediated immune responses associated with asymptomatic infection (AS) of SARS-CoV-2 remain largely unknown. The diversity of T-cell receptor (TCR) repertoire is essential for generating effective immunity against viral infections in T cell response. Here, we performed the single-cell TCR sequencing of the PBMC samples from five AS subjects, 33 symptomatic COVID-19 patients and eleven healthy controls to investigate the size and the diversity of TCR repertoire. We subsequently analyzed the TCR repertoire diversity, the V and J gene segment deference, and the dominant combination of αβ VJ gene pairing among these three study groups. Notably, we revealed significant TCR preference in the AS group, including the skewed usage of TRAV1-2-J33-TRBV6-4-J2-2 and TRAV1-2-J33-TRBV6-1-J2-3. Our findings may shed new light on understanding the immunopathogenesis of COVID-19 and help identify optimal TCRs for development of novel therapeutic strategies against SARS-CoV-2 infection.
Project description:Our research aimed to identify new therapeutic targets for Lung adenocarcinoma (LUAD), a major subtype of non-small cell lung cancer known for its low 5-year survival rate of 22%. By employing a comprehensive methodological approach, we analyzed bulk RNA sequencing data from 513 LUAD and 59 non-tumorous tissues, identifying 2,688 differentially expressed genes. Using Mendelian randomization (MR), we identified 74 genes with strong evidence for a causal effect on risk of LUAD. Survival analysis on these genes revealed significant differences in survival rates for 13 of them. Our pathway enrichment analysis highlighted their roles in immune response and cell communication, deepening our understanding. We also utilized single-cell RNA sequencing (scRNA-seq) to uncover cell type-specific gene expression patterns within LUAD, emphasizing the tumor microenvironment's heterogeneity. Pseudotime analysis further assisted in assessing the heterogeneity of tumor cell populations. Additionally, protein-protein interaction (PPI) network analysis was conducted to evaluate the potential druggability of these identified genes. The culmination of our efforts led to the identification of five genes (tier 1) with the most compelling evidence, including SECISBP2L, PRCD, SMAD9, C2orf91, and HSD17B13, and eight genes (tier 2) with convincing evidence for their potential as therapeutic targets.
Project description:T-cell receptor (TCR) is crucial in T cell-mediated virus clearance. To date, TCR bias has been observed in various diseases. However, studies on the TCR repertoire of COVID-19 patients are lacking. Here, we used single-cell V(D)J sequencing to conduct comparative analyses of TCR repertoire between 12 COVID-19 patients and 6 healthy controls, as well as other virus-infected samples. We observed distinct T cell clonal expansion in COVID-19. Further analysis of VJ gene combination revealed 6 VJ pairs significantly increased, while 139 pairs significantly decreased in COVID-19 patients. When considering the VJ combination of α and β chains at the same time, the combination with the highest frequency on COVID-19 was TRAV12-2-J27-TRBV7-9-J2-3. Besides, preferential usage of V and J gene segments was also observed in samples infected by different viruses. Our study provides novel insights on TCR in COVID-19, which contribute to our understanding of the immune response induced by SARS-CoV-2.
Project description:Moyamoya disease (MMD) remains a chronic progressive cerebrovascular disease with unknown etiology. A growing number of reports describe the development of MMD relevant to infection or autoimmune diseases. Identifying biomarkers of MMD is to understand the pathogenesis and development of novel targeted therapy and may be the key to improving the patient's outcome. Here, we analyzed gene expression from two GEO databases. To identify the MMD biomarkers, the weighted gene co-expression network analysis (WGCNA) and the differential expression analyses were conducted to identify 266 key genes. The KEGG and GO analyses were then performed to construct the protein interaction (PPI) network. The three machine-learning algorithms of support vector machine-recursive feature elimination (SVM-RFE), random forest and least absolute shrinkage and selection operator (LASSO) were used to analyze the key genes and take intersection to construct MMD diagnosis based on the four core genes found (ACAN, FREM1, TOP2A and UCHL1), with highly accurate AUCs of 0.805, 0.903, 0.815, 0.826. Gene enrichment analysis illustrated that the MMD samples revealed quite a few differences in pathways like one carbon pool by folate, aminoacyl-tRNA biosynthesis, fat digestion and absorption and fructose and mannose metabolism. In addition, the immune infiltration profile demonstrated that ACAN expression was associated with mast cells resting, FREM1 expression was associated with T cells CD4 naive, TOP2A expression was associated with B cells memory, UCHL1 expression was associated with mast cells activated. Ultimately, the four key genes were verified by qPCR. Taken together, our study analyzed the diagnostic biomarkers and immune infiltration characteristics of MMD, which may shed light on the potential intervention targets of moyamoya disease patients.
Project description:Background: Osteosarcoma (OS) is a kind of solid tumor with high heterogeneity at tumor microenvironment (TME), genome and transcriptome level. In view of the regulatory effect of metabolism on TME, this study was based on four metabolic models to explore the intertumoral heterogeneity of OS at the RNA sequencing (RNA-seq) level and the intratumoral heterogeneity of OS at the bulk RNA-seq and single cell RNA-seq (scRNA-seq) level. Methods: The GSVA package was used for single-sample gene set enrichment analysis (ssGSEA) analysis to obtain a glycolysis, pentose phosphate pathway (PPP), fatty acid oxidation (FAO) and glutaminolysis gene sets score. ConsensusClusterPlus was employed to cluster OS samples downloaded from the Target database. The scRNA-seq and bulk RNA-seq data of immune cells from GSE162454 dataset were analyzed to identify the subsets and types of immune cells in OS. Malignant cells and non-malignant cells were distinguished by large-scale chromosomal copy number variation. The correlations of metabolic molecular subtypes and immune cell types with four metabolic patterns, hypoxia and angiogenesis were determined by Pearson correlation analysis. Results: Two metabolism-related molecular subtypes of OS, cluster 1 and cluster 2, were identified. Cluster 2 was associated with poor prognosis of OS, active glycolysis, FAO, glutaminolysis, and bad TME. The identified 28608 immune cells were divided into 15 separate clusters covering 6 types of immune cells. The enrichment scores of 5 kinds of immune cells in cluster-1 and cluster-2 were significantly different. And five kinds of immune cells were significantly correlated with four metabolic modes, hypoxia and angiogenesis. Of the 28,608 immune cells, 7617 were malignant cells. The four metabolic patterns of malignant cells were significantly positively correlated with hypoxia and negatively correlated with angiogenesis. Conclusion: We used RNA-seq to reveal two molecular subtypes of OS with prognosis, metabolic pattern and TME, and determined the composition and metabolic heterogeneity of immune cells in OS tumor by bulk RNA-seq and single-cell RNA-seq.
Project description:Fibrosis is characterized by the excessive production of collagen and other extracellular matrix (ECM) components and represents a leading cause of morbidity and mortality worldwide. Previous studies of nonalcoholic steatohepatitis (NASH) with fibrosis were largely restricted to bulk transcriptome profiles. Thus, our understanding of this disease is limited by an incomplete characterization of liver cell types in general and hepatic stellate cells (HSCs) in particular, given that activated HSCs are the major hepatic fibrogenic cell population. To help fill this gap, we profiled 17,810 non-parenchymal cells derived from six healthy human livers. In conjunction with public single-cell data of fibrotic/cirrhotic human livers, these profiles enable the identification of potential intercellular signaling axes (e.g., ITGAV-LAMC1, TNFRSF11B-VWF and NOTCH2-DLL4) and master regulators (e.g., RUNX1 and CREB3L1) responsible for the activation of HSCs during fibrogenesis. Bulk RNA-seq data of NASH patient livers and rodent models for liver fibrosis of diverse etiologies allowed us to evaluate the translatability of candidate therapeutic targets for NASH-related fibrosis. We identified 61 liver fibrosis-associated genes (e.g., AEBP1, PRRX1 and LARP6) that may serve as a repertoire of translatable drug target candidates. Consistent with the above regulon results, gene regulatory network analysis allowed the identification of CREB3L1 as a master regulator of many of the 61 genes. Together, this study highlights potential cell-cell interactions and master regulators that underlie HSC activation and reveals genes that may represent prospective hallmark signatures for liver fibrosis.
Project description:Coronavirus disease 2019 (COVID-19) is a type of pneumonia caused by the SARS-CoV-2 coronavirus. It can cause acute pulmonary and systemic inflammation, which can lead to death in severely ill patients. This study explores the potential reasons behind severe COVID-19 and its similarities to systemic autoimmune diseases. This study reviewed unbiased high-throughput gene expression datasets, including next-generation and single-cell RNA sequencing. A total of 27 studies and eight meta-analyses were reviewed. The studies indicated that severe COVID-19 is associated with the upregulation of genes involved in pro-inflammatory, interferon, and cytokine/chemokine pathways. Additionally, changes were observed in the proportions of immune cell types in the blood and tissues, along with degenerative alterations in lung epithelial cells. Genomic evidence also supports the association of severe COVID-19 with various inflammatory syndromes, such as neuronal COVID-19, acute respiratory distress syndrome, vascular inflammation, and multisystem inflammatory syndrome. In conclusion, this study suggests that gene expression profiling plays a significant role in elucidating the etiology of severe COVID-19.
Project description:BackgroundPapillary thyroid cancer (PTC) is the most common pathological type of thyroid cancer with a high incidence globally. Increasing evidence reported that fibroblasts infiltration in cancer was correlated with prognostic outcomes. However, fibroblasts related study in thyroid cancer remains deficient.MethodsSingle-cell sequencing data of PTC were analyzed by Seurat R package to explore the ecosystem in PTC and identify fibroblasts cluster. The expression profiles and prognostic values of fibroblast related genes were assessed in TCGA dataset. A fibrosis score model was established for prognosis prediction in thyroid cancer patients. Differentially expressed genes and functional enrichment between high and low fibrosis score groups in TCGA dataset were screened. The correlation of immune cells infiltration and fibrosis score in thyroid cancer patients was explored. Expression levels and prognostic values of key fibroblast related factor were validated in clinical tissues another PTC cohort.ResultsFibroblasts were highly infiltrated in PTC and could interact with other type of cells by single-cell data analysis. 34 fibroblast related terms were differentially expressed in thyroid tumor tissues. COX regression analysis suggested that the constructed fibrosis score model was an independent prognostic predictor for thyroid cancer patients (HR = 5.17, 95%CI 2.31-11.56, P = 6.36E-05). Patients with low fibrosis scores were associated with a significantly better overall survival (OS) than those with high fibrosis scores in TCGA dataset (P = 7.659E-04). Specific immune cells infiltration levels were positively correlated with fibrosis score, including monocytes, M1 macrophages and eosinophils.ConclusionOur research demonstrated a comprehensive horizon of fibroblasts features in thyroid cancer microenvironment, which may provide potential value for thyroid cancer treatment.
Project description:We selected humann intervertebral disc samples to perform proteomics analysis. There were 1 case of grade I , 1 case of grade II, 3 cases of grade Ⅲ and 3 cases of grade Ⅳ according to Pfirrmann classfication. RNA seqencing analysis and single-cell RNA sequencing were integrated with proteomics data to identify the hub genes for intervertebral disc degeneration using bioinformatic method.