Computational deconvolution: extracting cell type-specific information from heterogeneous samples

The cGAS-STING-interferon regulatory factor 7 pathway regulates neuroinflammation in Parkinson's disease

JOURNAL/nrgr/04.03/01300535-202508000-00026/figure1/v/2024-09-30T120553Z/r/image-tiff Interferon regulatory factor 7 plays a crucial role in the innate immune response. However, whether interferon regulatory factor 7-mediated signaling contributes to Parkinson's disease remains unknown. Here we report that interferon regulatory factor 7 is markedly up-regulated in a 1-methyl-4-phenyl-1,2,3,6-tetrahydropyridine-induced mouse model of Parkinson's disease and co-localizes with microglial cells. Both the selective cyclic guanosine monophosphate adenosine monophosphate synthase inhibitor RU.521 and the stimulator of interferon genes inhibitor H151 effectively suppressed interferon regulatory factor 7 activation in BV2 microglia exposed to 1-methyl-4-phenylpyridinium and inhibited transformation of mouse BV2 microglia into the neurotoxic M1 phenotype. In addition, siRNA-mediated knockdown of interferon regulatory factor 7 expression in BV2 microglia reduced the expression of inducible nitric oxide synthase, tumor necrosis factor α, CD16, CD32, and CD86 and increased the expression of the anti-inflammatory markers ARG1 and YM1. Taken together, our findings indicate that the cyclic guanosine monophosphate adenosine monophosphate synthase-stimulator of interferon genes-interferon regulatory factor 7 pathway plays a crucial role in the pathogenesis of Parkinson's disease.

Identifying effective immune biomarkers in alopecia areata diagnosis based on machine learning methods

BACKGROUND

Alopecia areata (AA) is a common non-scarring hair loss disorder associated with autoimmune conditions. However, the pathobiology of AA is not well understood, and there is no targeted therapy available for AA.  METHODS: In this study, differential gene expression analysis, immune status assessment, weighted correlation network analysis (WGCNA), and functional enrichment analysis were performed to identify shared genes associated with both immunological response and AA. Machine learning methods were then used to identify three hub genes as potential diagnostic markers for AA. External validation was performed, and the correlation of hub genes with immune infiltration, immune checkpoint genes, and key marker genes and pathways were evaluated.

RESULTS

Three hub genes were identified, which accurately predicted the progression of AA and the immune status. The hub genes were found to be diagnostic markers for AA with high predictive accuracy. External validation confirmed the efficacy of these markers in identifying AA patients.

CONCLUSION

Overall, the study provides a novel approach for the diagnosis, prevention, and treatment of AA. The findings could potentially lead to the development of targeted therapies for AA based on the identified hub genes. The study also highlights the potential of machine learning and bioinformatics analysis in identifying new biomarkers for autoimmune diseases.

Using Callus as an Ex Vivo System for Chromatin Analysis

Next-generation sequencing has revolutionized epigenetics research, enabling a comprehensive analysis of DNA methylation and histone modification profiles to explore complex biological systems at unprecedented depth. Deciphering the intricate epigenetic mechanisms that regulate gene activity presents significant challenges, including the issue of analyzing heterogeneous cell populations in bulk. Bulk analysis introduces bias and can obscure crucial information by averaging readouts from distinct cells. Various approaches have been developed to address this issue, such as cell-type-specific enrichment or single-cell sequencing techniques. However, the need for transgenic lines with fluorescent markers, along with technical challenges such as efficient protoplast isolation and low yield, limits their widespread adoption and use in multi-omic studies. This review discusses the pros and cons of these approaches, providing a valuable basis for selecting the most suitable strategy to minimize heterogeneity. We will also highlight the use of cotyledon-derived callus as an ex vivo system as a simple, accessible, and robust platform for enabling high-throughput multi-omic analyses.

A twin analysis to estimate genetic and environmental factors contributing to variation in weighted gene co-expression network module eigengenes

Multivariate network-based analytic methods such as weighted gene co-expression network analysis are frequently applied to human and animal gene-expression data to estimate the first principal component of a module, or module eigengene (ME). MEs are interpreted as multivariate summaries of correlated gene-expression patterns and network connectivity across genes within a module. As such, they have the potential to elucidate the mechanisms by which molecular genomic variation contributes to individual differences in complex traits. Although increasingly used to test for associations between modules and complex traits, the genetic and environmental etiology of MEs has not been empirically established. It is unclear if, and to what degree, individual differences in blood-derived MEs reflect random variation versus familial aggregation arising from heritable or shared environmental influences. We used biometrical genetic analyses to estimate the contribution of genetic and environmental influences on MEs derived from blood lymphocytes collected on a sample of N = 661 older male twins from the Vietnam Era Twin Study of Aging (VETSA) whose mean age at assessment was 67.7 years (SD = 2.6 years, range = 62-74 years). Of the 26 detected MEs, 14 (56%) had statistically significant additive genetic variation with an average heritability of 44% (SD = 0.08, range = 35%-64%). Despite the relatively small sample size, this demonstration of significant family aggregation including estimates of heritability in 14 of the 26 MEs suggests that blood-based MEs are reliable and merit further exploration in terms of their associations with complex traits and diseases.

VarNMF: non-negative probabilistic factorization with source variation

MOTIVATION

Non-negative matrix factorization (NMF) is a powerful tool often applied to genomic data to identify non-negative latent components that constitute linearly mixed samples. It is useful when the observed signal combines contributions from multiple sources, such as cell types in bulk measurements of heterogeneous tissue. NMF accounts for two types of variation between samples - disparities in the proportions of sources and observation noise. However, in many settings, there is also a non-trivial variation between samples in the contribution of each source to the mixed data. This variation cannot be accurately modeled using the NMF framework.

RESULTS

We present VarNMF, a probabilistic extension of NMF that explicitly models this variation in source values. We show that by modeling sources as non-negative distributions, we can recover source variation directly from mixed samples without observing any of the sources directly. We apply VarNMF to a cell-free ChIP-seq dataset of two cancer cohorts and a healthy cohort, demonstrating that VarNMF provides a better estimation of the data distribution. Moreover, VarNMF extracts cancer-associated source distributions that decouple the tumor characteristics from the amount of tumor contribution, and identify patient-specific disease behaviors. This decomposition highlights the inter-tumor variability that is obscured in the mixed samples.

AVAILABILITY AND IMPLEMENTATION

Code is available at https://github.com/Nir-Friedman-Lab/VarNMF.

Luminal epithelial cells integrate variable responses to aging into stereotypical changes that underlie breast cancer susceptibility

Effects from aging in single cells are heterogenous, whereas at the organ- and tissue-levels aging phenotypes tend to appear as stereotypical changes. The mammary epithelium is a bilayer of two major phenotypically and functionally distinct cell lineages: luminal epithelial and myoepithelial cells. Mammary luminal epithelia exhibit substantial stereotypical changes with age that merit attention because these cells are the putative cells-of-origin for breast cancers. We hypothesize that effects from aging that impinge upon maintenance of lineage fidelity increase susceptibility to cancer initiation. We generated and analyzed transcriptomes from primary luminal epithelial and myoepithelial cells from younger <30 (y)ears old and older >55 y women. In addition to age-dependent directional changes in gene expression, we observed increased transcriptional variance with age that contributed to genome-wide loss of lineage fidelity. Age-dependent variant responses were common to both lineages, whereas directional changes were almost exclusively detected in luminal epithelia and involved altered regulation of chromatin and genome organizers such as SATB1. Epithelial expression variance of gap junction protein GJB6 increased with age, and modulation of GJB6 expression in heterochronous co-cultures revealed that it provided a communication conduit from myoepithelial cells that drove directional change in luminal cells. Age-dependent luminal transcriptomes comprised a prominent signal that could be detected in bulk tissue during aging and transition into cancers. A machine learning classifier based on luminal-specific aging distinguished normal from cancer tissue and was highly predictive of breast cancer subtype. We speculate that luminal epithelia are the ultimate site of integration of the variant responses to aging in their surrounding tissue, and that their emergent phenotype both endows cells with the ability to become cancer-cells-of-origin and represents a biosensor that presages cancer susceptibility.

TNF-α and RPLP0 drive the apoptosis of endothelial cells and increase susceptibility to high-altitude pulmonary edema

High-altitude pulmonary edema (HAPE) is a fatal threat for sojourners who ascend rapidly without sufficient acclimatization. Acclimatized sojourners and adapted natives are both insensitive to HAPE but have different physiological traits and molecular bases. In this study, based on GSE52209, the gene expression profiles of HAPE patients were compared with those of acclimatized sojourners and adapted natives, with the common and divergent differentially expressed genes (DEGs) and their hub genes identified, respectively. Bioinformatic methodologies for functional enrichment analysis, immune infiltration, diagnostic model construction, competing endogenous RNA (ceRNA) analysis and drug prediction were performed to detect potential biological functions and molecular mechanisms. Next, an array of in vivo experiments in a HAPE rat model and in vitro experiments in HUVECs were conducted to verify the results of the bioinformatic analysis. The enriched pathways of DEGs and immune landscapes for HAPE were significantly different between sojourners and natives, and the common DEGs were enriched mainly in the pathways of development and immunity. Nomograms revealed that the upregulation of TNF-α and downregulation of RPLP0 exhibited high diagnostic efficiency for HAPE in both sojourners and natives, which was further validated in the HAPE rat model. The addition of TNF-α and RPLP0 knockdown activated apoptosis signaling in endothelial cells (ECs) and enhanced endothelial permeability. In conclusion, TNF-α and RPLP0 are shared biomarkers and molecular bases for HAPE susceptibility during the acclimatization/adaptation/maladaptation processes in sojourners and natives, inspiring new ideas for predicting and treating HAPE.

Multi-modal transcriptomic analysis reveals metabolic dysregulation and immune responses in chronic obstructive pulmonary disease

Chronic obstructive pulmonary disease (COPD), a progressive inflammatory condition of the airways, emerges from the complex interplay between genetic predisposition and environmental factors. Notably, its incidence is on the rise, particularly among the elderly demographic. Current research increasingly highlights cellular senescence as a key driver in chronic lung pathologies. Despite this, the detailed mechanisms linking COPD with senescent genomic alterations remain elusive. To address this gap, there is a pressing need for comprehensive bioinformatics methodologies that can elucidate the molecular intricacies of this link. This approach is crucial for advancing our understanding of COPD and its association with cellular aging processes. Utilizing a spectrum of advanced bioinformatics techniques, this research delved into the potential mechanisms linking COPD with aging-related genes, identifying four key genes (EP300, MTOR, NFE2L1, TXN) through machine learning and weighted gene co-expression network analysis (WGCNA) analyses. Subsequently, a precise diagnostic model leveraging an artificial neural network was developed. The study further employed single-cell analysis and molecular docking to investigate senescence-related cell types in COPD tissues, particularly focusing on the interactions between COPD and NFE2L1, thereby enhancing the understanding of COPD's molecular underpinnings. Leveraging artificial neural networks, we developed a robust classification model centered on four genes-EP300, MTOR, NFE2L1, TXN-exhibiting significant predictive capability for COPD and offering novel avenues for its early diagnosis. Furthermore, employing various single-cell analysis techniques, the study intricately unraveled the characteristics of senescence-related cell types in COPD tissues, enriching our understanding of the disease's cellular landscape. This research anticipates offering novel biomarkers and therapeutic targets for early COPD intervention, potentially alleviating the disease's impact on individuals and healthcare systems, and contributing to a reduction in global COPD-related mortality. These findings carry significant clinical and public health ramifications, bolstering the foundation for future research and clinical strategies in managing and understanding COPD.

Deconvolution from bulk gene expression by leveraging sample-wise and gene-wise similarities and single-cell RNA-Seq data

BACKGROUND

The widely adopted bulk RNA-seq measures the gene expression average of cells, masking cell type heterogeneity, which confounds downstream analyses. Therefore, identifying the cellular composition and cell type-specific gene expression profiles (GEPs) facilitates the study of the underlying mechanisms of various biological processes. Although single-cell RNA-seq focuses on cell type heterogeneity in gene expression, it requires specialized and expensive resources and currently is not practical for a large number of samples or a routine clinical setting. Recently, computational deconvolution methodologies have been developed, while many of them only estimate cell type composition or cell type-specific GEPs by requiring the other as input. The development of more accurate deconvolution methods to infer cell type abundance and cell type-specific GEPs is still essential.

RESULTS

We propose a new deconvolution algorithm, DSSC, which infers cell type-specific gene expression and cell type proportions of heterogeneous samples simultaneously by leveraging gene-gene and sample-sample similarities in bulk expression and single-cell RNA-seq data. Through comparisons with the other existing methods, we demonstrate that DSSC is effective in inferring both cell type proportions and cell type-specific GEPs across simulated pseudo-bulk data (including intra-dataset and inter-dataset simulations) and experimental bulk data (including mixture data and real experimental data). DSSC shows robustness to the change of marker gene number and sample size and also has cost and time efficiencies.

CONCLUSIONS

DSSC provides a practical and promising alternative to the experimental techniques to characterize cellular composition and heterogeneity in the gene expression of heterogeneous samples.

Novel Insights into Post-Myocardial Infarction Cardiac Remodeling through Algorithmic Detection of Cell-Type Composition Shifts

Background

Recent advances in single cell sequencing have led to an increased focus on the role of cell-type composition in phenotypic presentation and disease progression. Cell-type composition research in the heart is challenging due to large, frequently multinucleated cardiomyocytes that preclude most single cell approaches from obtaining accurate measurements of cell composition. Our in silico studies reveal that ignoring cell type composition when calculating differentially expressed genes (DEGs) can have significant consequences. For example, a relatively small change in cell abundance of only 10% can result in over 25% of DEGs being false positives.

Methods

We have implemented an algorithmic approach that uses snRNAseq datasets as a reference to accurately calculate cell type compositions from bulk RNAseq datasets through robust data cleaning, gene selection, and multi-sample cross-subject and cross-cell-type deconvolution. We applied our approach to cardiomyocyte-specific α1A adrenergic receptor (CM-α1A-AR) knockout mice. 8-12 week-old mice (either WT or CM-α1A-KO) were subjected to permanent left coronary artery (LCA) ligation or sham surgery (n=4 per group). Transcriptomes from the infarct border zones were collected 3 days later and analyzed using our algorithm to determine cell-type abundances, corrected differential expression calculations using DESeq2, and validated these findings using RNAscope.

Results

Uncorrected DEGs for the CM-α1A-KO X LCA interaction term featured many cell-type specific genes such as Timp4 (fibroblasts) and Aplnr (cardiomyocytes) and overall GO enrichment for terms pertaining to cardiomyocyte differentiation (P=3.1E-4). Using our algorithm, we observe a striking loss of cardiomyocytes and gain in fibroblasts in the α1A-KO + LCA mice that was not recapitulated in WT + LCA animals, although we did observe a similar increase in macrophage abundance in both conditions. This recapitulates prior results that showed a much more severe heart failure phenotype in CM-α1A-KO + LCA mice. Following correction for cell-type, our DEGs now highlight a novel set of genes enriched for GO terms such as cardiac contraction (P=3.7E-5) and actin filament organization (P=6.3E-5).

Conclusions

Our algorithm identifies and corrects for cell-type abundance in bulk RNAseq datasets opening new avenues for research on novel genes and pathways as well as an improved understanding of the role of cardiac cell types in cardiovascular disease.

The Mechanism of "Treating Different Diseases with the Same Treatment" by Qiangji Jianpi Decoction in Ankylosing Spondylitis Combined with Inflammatory Bowel Disease: A Comprehensive Analysis of Multiple Methods

Background

Ankylosing spondylitis (AS) and inflammatory bowel disease (IBD) are prevalent autoimmune disorders that often co-occur, posing significant treatment challenges. This investigation adopts a multidisciplinary strategy, integrating bioinformatics, network pharmacology, molecular docking, and Mendelian randomization, to elucidate the relationship between AS and IBD and to investigate the potential mechanisms of traditional Chinese medicine formulations, represented by Qiangji Jianpi (QJJP) decoction, in treating these comorbid conditions.

Methods

We utilized databases to pinpoint common targets among AS, IBD, and QJJP decoction's active compounds through intersection analysis. Through Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses, we mapped a network in Cytoscape, isolating critical targets. Molecular docking with AutoDock validated the affinity between targets and compounds. ROC analysis and dataset validation assessed diagnostic performance, while Gene Set Enrichment Analysis (GSEA) offered pathway insights. Mendelian randomization explored the AS-IBD causal relationship.

Results

Screening identified 105 targets for QJJP decoction, 414 for AS, and 2420 for IBD, with 85 overlapping. These targets predominantly participate in organismal responses and DNA transcription factor binding, with a significant cellular presence in the endoplasmic reticulum and vesicle lumen. Molecular docking, facilitated by Cytoscape, confirmed IL1A, IFNG, TGFB1, and EDN1 as critical targets, with IFNG demonstrating diagnostic potential through GEO dataset validation. The integration of GSEA with network pharmacology highlighted the therapeutic significance of the relaxin, osteoclast differentiation, HIF-1, and AGE-RAGE signaling pathways in QJJP decoction's action. Mendelian randomization analysis indicated a positive causal relationship between IBD and AS, pinpointing rs2193041 as a key SNP influencing IFNG.

Conclusion

Based on the principle of "treating different diseases with the same method" in traditional Chinese medicine theory, we explored the intricate mechanisms through which QJJP decoction addresses AS and IBD comorbidity. Our research spotlighted the pivotal role of the IFNG gene. IFNG emerges not only as a key therapeutic target but also assumes significance as a potential diagnostic biomarker through its genetic underpinnings. This investigation establishes a solid base for subsequent experimental inquiries. Our findings introduce novel approaches for incorporating traditional Chinese medicine into the treatment of AS-IBD comorbidity, setting the stage for groundbreaking research directions.

Benchmarking of methods for DNA methylome deconvolution

Defining the number and abundance of different cell types in tissues is important for understanding disease mechanisms as well as for diagnostic and prognostic purposes. Typically, this is achieved by immunohistological analyses, cell sorting, or single-cell RNA-sequencing. Alternatively, cell-specific DNA methylome information can be leveraged to deconvolve cell fractions from a bulk DNA mixture. However, comprehensive benchmarking of deconvolution methods and modalities was not yet performed. Here we evaluate 16 deconvolution algorithms, developed either specifically for DNA methylome data or more generically. We assess the performance of these algorithms, and the effect of normalization methods, while modeling variables that impact deconvolution performance, including cell abundance, cell type similarity, reference panel size, method for methylome profiling (array or sequencing), and technical variation. We observe differences in algorithm performance depending on each these variables, emphasizing the need for tailoring deconvolution analyses. The complexity of the reference, marker selection method, number of marker loci and, for sequencing-based assays, sequencing depth have a marked influence on performance. By developing handles to select the optimal analysis configuration, we provide a valuable source of information for studies aiming to deconvolve array- or sequencing-based methylation data.

Predicting drug response through tumor deconvolution by cancer cell lines

Large-scale cancer drug sensitivity data have become available for a collection of cancer cell lines, but only limited drug response data from patients are available. Bridging the gap in pharmacogenomics knowledge between in vitro and in vivo datasets remains challenging. In this study, we trained a deep learning model, Scaden-CA, for deconvoluting tumor data into proportions of cancer-type-specific cell lines. Then, we developed a drug response prediction method using the deconvoluted proportions and the drug sensitivity data from cell lines. The Scaden-CA model showed excellent performance in terms of concordance correlation coefficients (>0.9 for model testing) and the correctly deconvoluted rate (>70% across most cancers) for model validation using Cancer Cell Line Encyclopedia (CCLE) bulk RNA data. We applied the model to tumors in The Cancer Genome Atlas (TCGA) dataset and examined associations between predicted cell viability and mutation status or gene expression levels to understand underlying mechanisms of potential value for drug repurposing.

Interaction molecular QTL mapping discovers cellular and environmental modifiers of genetic regulatory effects

Bulk-tissue molecular quantitative trait loci (QTLs) have been the starting point for interpreting disease-associated variants, and context-specific QTLs show particular relevance for disease. Here, we present the results of mapping interaction QTLs (iQTLs) for cell type, age, and other phenotypic variables in multi-omic, longitudinal data from the blood of individuals of diverse ancestries. By modeling the interaction between genotype and estimated cell-type proportions, we demonstrate that cell-type iQTLs could be considered as proxies for cell-type-specific QTL effects, particularly for the most abundant cell type in the tissue. The interpretation of age iQTLs, however, warrants caution because the moderation effect of age on the genotype and molecular phenotype association could be mediated by changes in cell-type composition. Finally, we show that cell-type iQTLs contribute to cell-type-specific enrichment of diseases that, in combination with additional functional data, could guide future functional studies. Overall, this study highlights the use of iQTLs to gain insights into the context specificity of regulatory effects.

scMD facilitates cell type deconvolution using single-cell DNA methylation references

The proliferation of single-cell RNA-sequencing data has led to the widespread use of cellular deconvolution, aiding the extraction of cell-type-specific information from extensive bulk data. However, those advances have been mostly limited to transcriptomic data. With recent developments in single-cell DNA methylation (scDNAm), there are emerging opportunities for deconvolving bulk DNAm data, particularly for solid tissues like brain that lack cell-type references. Due to technical limitations, current scDNAm sequences represent a small proportion of the whole genome for each single cell, and those detected regions differ across cells. This makes scDNAm data ultra-high dimensional and ultra-sparse. To deal with these challenges, we introduce scMD (single cell Methylation Deconvolution), a cellular deconvolution framework to reliably estimate cell type fractions from tissue-level DNAm data. To analyze large-scale complex scDNAm data, scMD employs a statistical approach to aggregate scDNAm data at the cell cluster level, identify cell-type marker DNAm sites, and create precise cell-type signature matrixes that surpass state-of-the-art sorted-cell or RNA-derived references. Through thorough benchmarking in several datasets, we demonstrate scMD's superior performance in estimating cellular fractions from bulk DNAm data. With scMD-estimated cellular fractions, we identify cell type fractions and cell type-specific differentially methylated cytosines associated with Alzheimer's disease.

Baseline Peripheral Blood Mononuclear Cell Transcriptomics Before Ustekinumab Treatment Is Linked With Crohn's Disease Clinical Response at 1 Year

INTRODUCTION

Ustekinumab, a monoclonal antibody to the p40 subunit of interleukin (IL)-12 and IL-23, is used for Crohn's disease (CD), and the documented clinical remission rate after 1 year was observed in approximately 50% of patients. We aimed to identify predictors for a clinical response using peripheral blood obtained from patients with CD just before ustekinumab treatment initiation.

METHODS

RNA extraction from peripheral blood mononuclear cells was followed by mRNA paired-end sequencing. Differential gene expression was performed using DESeq2.

RESULTS

We processed samples from 36 adults with CD (13 men, 36%) obtained at baseline before starting ustekinumab treatment. Twenty-two of 36 (61%) were defined as responders and 14/36 (39%) as nonresponders after 1 year based on Physician Global Assessment. Differential gene expression between responders (n = 22) and nonresponders (n = 14) did not show a gene expression signature that passed false discovery rate (FDR) correction. However, the analyses identified 68 genes, including CXCL1/2/3, which were induced in nonresponders vs responders with P < 0.05 and fold change above 1.5. Functional annotation enrichments of these 68 genes using ToppGene indicated enrichment for cytokine activity (FDR = 1.98E-05), CXCR chemokine receptor binding (FDR = 2.11E-05), IL-10 signaling (FDR = 5.03E-07), genes encoding secreted soluble factors (FDR = 1.73E-05), and myeloid dendritic cells (FDR = 1.80E-08).

DISCUSSION

No substantial differences were found in peripheral blood mononuclear cell transcriptomics between responders and nonresponders. However, among the nonresponders, we noted an increased inflammatory response enriched for pathways linked with cytokine activity and chemokine receptor binding and innate myeloid signature. A larger cohort is required to validate and further explore these findings.

Widespread sex dimorphism across single-cell transcriptomes of adult African turquoise killifish tissues

The African turquoise killifish (Nothobranchius furzeri), the shortest-lived vertebrate that can be bred in captivity, is an emerging model organism for aging research. Here, we describe a multitissue, single-cell gene expression atlas of female and male blood, kidney, liver, and spleen. We annotate 22 cell types, define marker genes, and infer differentiation trajectories. We find pervasive sex-dimorphic gene expression across cell types. Sex-dimorphic genes tend to be linked to lipid metabolism, consistent with clear differences in lipid storage in female vs. male turquoise killifish livers. We use machine learning to predict sex using single-cell gene expression and identify potential markers for molecular sex identity. As a proof of principle, we show that our atlas can be used to deconvolute existing bulk RNA sequencing (RNA-seq) data to obtain accurate estimates of cell type proportions. This atlas can be a resource to the community that could be leveraged to develop cell-type-specific expression in transgenic animals.

LncRNA RARA-AS1 could serve as a novel prognostic biomarker in pan-cancer and promote proliferation and migration in glioblastoma

Long non-coding RNAs (lncRNAs) have emerged as crucial regulators of cancer progression and are potential biomarkers for diagnosis and treatment. This study investigates the role of RARA Antisense RNA 1 (RARA-AS1) in cancer and its implications for diagnosis and treatment. Various bioinformatics tools were conducted to analyze the expression patterns, immune-related functions, methylation, and gene expression correlations of RARA-AS1, mainly including the comparisons of different subgroups and correlation analyses between RARA-AS1 expression and other factors. Furthermore, we used short hairpin RNA to perform knockdown experiments, investigating the effects of RARA-AS1 on cell proliferation, invasion, and migration in glioblastoma. Our results revealed that RARA-AS1 has distinct expression patterns in different cancers and exhibits notable correlation with prognosis. Additionally, RARA-AS1 is highly correlated with certain immune checkpoints and mismatch repair genes, indicating its potential role in immune infiltration and related immunotherapy. Further analysis identified potential effective drugs for RARA-AS1 and demonstrated its potential RNA binding protein (RBP) mechanism in glioblastoma. Besides, a series of functional experiments indicated inhibiting RARA-AS1 could decrease cell proliferation, invasion, and migration of glioblastoma cell lines. Finally, RARA-AS1 could act as an independent prognostic factor for glioblastoma patients and may serve as a promising therapeutic target. All in all, Our study provides a comprehensive understanding of the functions and implications of RARA-AS1 in pan-cancer, highlighting it as a promising biomarker for survival. It is also an independent risk factor affecting prognosis in glioblastoma and an important factor affecting proliferation and migration in glioblastoma, setting the stage for further mechanistic investigations.

The T Cell Immunoscore as a Reference for Biomarker Development Utilizing Real-World Data from Patients with Advanced Malignancies Treated with Immune Checkpoint Inhibitors

BACKGROUND

We aimed to determine the prognostic value of an immunoscore reflecting CD3+ and CD8+ T cell density estimated from real-world transcriptomic data of a patient cohort with advanced malignancies treated with immune checkpoint inhibitors (ICIs) in an effort to validate a reference for future machine learning-based biomarker development.

METHODS

Transcriptomic data was collected under the Total Cancer Care Protocol (NCT03977402) Avatar® project. The real-world immunoscore for each patient was calculated based on the estimated densities of tumor CD3+ and CD8+ T cells utilizing CIBERSORTx and the LM22 gene signature matrix. Then, the immunoscore association with overall survival (OS) was estimated using Cox regression and analyzed using Kaplan-Meier curves. The OS predictions were assessed using Harrell's concordance index (C-index). The Youden index was used to identify the optimal cut-off point. Statistical significance was assessed using the log-rank test.

RESULTS

Our study encompassed 522 patients with four cancer types. The median duration to death was 10.5 months for the 275 participants who encountered an event. For the entire cohort, the results demonstrated that transcriptomics-based immunoscore could significantly predict patients at risk of death (p-value < 0.001). Notably, patients with an intermediate-high immunoscore achieved better OS than those with a low immunoscore. In subgroup analysis, the prediction of OS was significant for melanoma and head and neck cancer patients but did not reach significance in the non-small cell lung cancer or renal cell carcinoma cohorts.

CONCLUSIONS

Calculating CD3+ and CD8+ T cell immunoscore using real-world transcriptomic data represents a promising signature for estimating OS with ICIs and can be used as a reference for future machine learning-based biomarker development.

GTM-decon: guided-topic modeling of single-cell transcriptomes enables sub-cell-type and disease-subtype deconvolution of bulk transcriptomes

Cell-type composition is an important indicator of health. We present Guided Topic Model for deconvolution (GTM-decon) to automatically infer cell-type-specific gene topic distributions from single-cell RNA-seq data for deconvolving bulk transcriptomes. GTM-decon performs competitively on deconvolving simulated and real bulk data compared with the state-of-the-art methods. Moreover, as demonstrated in deconvolving disease transcriptomes, GTM-decon can infer multiple cell-type-specific gene topic distributions per cell type, which captures sub-cell-type variations. GTM-decon can also use phenotype labels from single-cell or bulk data to infer phenotype-specific gene distributions. In a nested-guided design, GTM-decon identified cell-type-specific differentially expressed genes from bulk breast cancer transcriptomes.

Fine-grained cell-type specific association studies with human bulk brain data using a large single-nucleus RNA sequencing based reference panel

Brain disorders are leading causes of disability worldwide. Gene expression studies provide promising opportunities to better understand their etiology but it is critical that expression is studied on a cell-type level. Cell-type specific association studies can be performed with bulk expression data using statistical methods that capitalize on cell-type proportions estimated with the help of a reference panel. To create a fine-grained reference panel for the human prefrontal cortex, we performed an integrated analysis of the seven largest single nucleus RNA-seq studies. Our panel included 17 cell-types that were robustly detected across all studies, subregions of the prefrontal cortex, and sex and age groups. To estimate the cell-type proportions, we used an empirical Bayes estimator that substantially outperformed three estimators recommended previously after a comprehensive evaluation of methods to estimate cell-type proportions from brain transcriptome data. This is important as being able to precisely estimate the cell-type proportions may avoid unreliable results in downstream analyses particularly for the multiple cell-types that had low abundances. Transcriptome-wide association studies performed with permuted bulk expression data showed that it is possible to perform transcriptome-wide association studies for even the rarest cell-types without an increased risk of false positives.

Identification and analysis of key hypoxia- and immune-related genes in hypertrophic cardiomyopathy

BACKGROUND

Hypertrophic cardiomyopathy (HCM), an autosomal dominant genetic disease, is the main cause of sudden death in adolescents and athletes globally. Hypoxia and immune factors have been revealed to be related to the pathology of HCM. There is growing evidence of a role for hypoxia and inflammation as triggers and enhancers in the pathology in HCM. However, the role of hypoxia- and immune-related genes in HCM have not been reported.

METHODS

Firstly, we obtained four HCM-related datasets from the Gene Expression Omnibus (GEO) database for differential expression analysis. Immune cells significantly expressed in normal samples and HCM were then screened by a microenvironmental cell population counter (MCP-counter) algorithm. Next, hypoxia- and immune-related genes were screened by the LASSO + support vector machine recursive feature elimination (SVM-RFE) and weighted gene co-expression network analysis (WGCNA). Single-gene enrichment analysis and expression validation of key genes were then performed. Finally, we constructed a competing endogenous RNA (ceRNA) network of key genes.

RESULTS

In this study, 35 differentially expressed hypoxia genes were found. By using LASSO + SVM-RFE analysis, 10 more targets with differentially expressed hypoxia genes were identified. The MCP-count algorithm yielded five differentially expressed immune cells, and after assessing them for WGCNA characteristics, 612 immune genes were discovered. When hypoxia and immune genes were combined for cross-tabulation analysis, three hypoxia- and immune-related genes (ATP2A2, DDAH1, and OMA1) were identified.

CONCLUSION

Based on hypoxia characteristic genes, three key genes were identified. These were also significantly related to immune activation, which proves a theoretical basis and reference value for studying the relationship between HCM and hypoxia and immunity.

scMD: cell type deconvolution using single-cell DNA methylation references

The proliferation of single-cell RNA sequencing data has led to the widespread use of cellular deconvolution, aiding the extraction of cell type-specific information from extensive bulk data. However, those advances have been mostly limited to transcriptomic data. With recent development in single-cell DNA methylation (scDNAm), new avenues have been opened for deconvolving bulk DNAm data, particularly for solid tissues like the brain that lack cell-type references. Due to technical limitations, current scDNAm sequences represent a small proportion of the whole genome for each single cell, and those detected regions differ across cells. This makes scDNAm data ultra-high dimensional and ultra-sparse. To deal with these challenges, we introduce scMD (single cell Methylation Deconvolution), a cellular deconvolution framework to reliably estimate cell type fractions from tissue-level DNAm data. To analyze large-scale complex scDNAm data, scMD employs a statistical approach to aggregate scDNAm data at the cell cluster level, identify cell-type marker DNAm sites, and create a precise cell-type signature matrix that surpasses state-of-the-art sorted-cell or RNA-derived references. Through thorough benchmarking in several datasets, we demonstrate scMD's superior performance in estimating cellular fractions from bulk DNAm data. With scMD-estimated cellular fractions, we identify cell type fractions and cell type-specific differentially methylated cytosines associated with Alzheimer's disease.

Methylome-wide association study of anxiety disorders

Anxiety Disorders (ANX) such as panic disorder, generalized anxiety disorder, and phobias, are highly prevalent conditions that are moderately heritable. Evidence suggests that DNA methylation may play a role, as it is involved in critical adaptations to changing environments. Applying an enrichment-based sequencing approach covering nearly 28 million autosomal CpG sites, we conducted a methylome-wide association study (MWAS) of lifetime ANX in 1132 participants (618 cases/514 controls) from the Netherlands Study of Depression and Anxiety. Using epigenomic deconvolution, we performed MWAS for the main cell types in blood: granulocytes, T-cells, B-cells and monocytes. Cell-type specific analyses identified 280 and 82 methylome-wide significant associations (q-value < 0.1) in monocytes and granulocytes, respectively. Our top finding in monocytes was located in ZNF823 on chromosome 19 (p = 1.38 × 10-10) previously associated with schizophrenia. We observed significant overlap (p < 1 × 10-06) with the same direction of effect in monocytes (210 sites), T-cells (135 sites), and B-cells (727 sites) between this Discovery MWAS signal and a comparable replication dataset from the Great Smoky Mountains Study (N = 433). Overlapping Discovery-Replication MWAS signal was enriched for findings from published GWAS of ANX, major depression, and post-traumatic stress disorder. In monocytes, two specific sites in the FZR1 gene showed significant replication after Bonferroni correction with an additional 15 nominally replicated sites in monocytes and 4 in T-cells. FZR1 regulates neurogenesis in the hippocampus, and its knockout leads to impairments in associative fear memory and long-term potentiation in mice. In the largest and most extensive methylome-wide study of ANX, we identified replicable methylation sites located in genes of potential relevance for brain mechanisms of psychiatric conditions.

Identification of the Key Immune Cells and Genes for the Diagnostics and Therapeutics of Meningioma

BACKGROUND

Dysregulation of immune infiltration critically contributes to the tumorigenesis and progression of meningiomas. However, the landscape of immune microenvironment and key genes correlated with immune cell infiltration remains unclear.

METHODS

Four Gene Expression Omnibus data sets were included. CIBERSORT algorithm was utilized to analyze the immune cell infiltration in samples. Wilcoxon test, Random Forest algorithm, and Least Absolute Shrinkage and Selection Operator regression were adopted in identifying significantly different infiltrating immune cells and differentially expressed genes (DEGs). Functional enrichment analysis was performed by Kyoto Encyclopedia of Genes and Genomes and Gene Ontology. The correlation between genes and immune cells was evaluated via Spearman's correlation analysis. Receiver Operator Characteristic curve analysis evaluated the markers' diagnostic effectiveness. The mRNA-miRNA and Drug-Gene-Immune cell interaction networks were constructed to identify potential diagnostic and therapeutic targets.

RESULTS

Plasma cells, M1 macrophages, M2 macrophages, neutrophils, eosinophils, and activated NK cells were the significantly different infiltrating immune cells in meningioma. A total of 951 DEGs, associated with synaptic function and structure, ion transport regulation, brain function, and immune-related pathways, were identified. Among 11 hub DEGs, RYR2 and TTR were correlated with plasma cells; SNCG was associated with NK cells; ADCY1 exhibited excellent diagnostic effectiveness; and ADCY1, BMX, KCNA5, SLCO4A1, and TTR could be considered as therapeutic targets.

CONCLUSIONS

ADCY1 can be identified as a diagnostic marker; ADCY1, BMX, KCNA5, SLCO4A1, and TTR are potential therapeutic targets, and their associations with macrophages, neutrophils, NK cells, and plasma cells might impact the tumorigenesis of meningiomas.

Therapeutic implications for localized prostate cancer by multiomics analyses of the ageing microenvironment landscape

Background: Numerous studies have substantiated the association between aging and the progression of malignant tumors in humans, notably prostate cancer (PCa). Nevertheless, to the best of our knowledge, no studies have comprehensively elucidated the intricate characteristics of the aging microenvironment (AME) in PCa. Methods: AME regulatory patterns were determined using the NMF algorithm. Then an ageing microenvironment index (AMI) was constructed, with excellent prognostic and immunotherapy prediction ability, and its' clinical relevance was surveyed through spatial transcriptomics. Further, the drug response was analysed using the Genomics of Drug Sensitivity in Cancer (GDSC), the Connectivity Map (CMap) and CellMiner database for patients with PCa. Finally, the AME was studied using in vitro and vivo experiments. Results: Three different AME regulatory patterns were identified across 813 PCa patients, associated with distinct clinical prognosis and physiological pathways. Based on the AMI, patients with PCa were divided into the high-score and low-score subsets. Higher AMI score was significantly infiltrated with more immune cells, higher rate of biochemical recurrence (BCR) and worse response to immunotherapy, antiandrogen therapy and chemotherapy in PCa. In addition, we found that the combination of bicalutamide and embelin was capable of suppressing tumor growth of PCa. Besides, as the main components of AMI, COL1A1 and BGLAP act as oncogenes and were verified via in vivo and in vitro experiments. Conclusions: AME regulation is significantly associated with the diversity and complexity of TME. Quantitative evaluation of the AME regulatory patterns may provide promising novel molecular markers for individualised therapy in PCa.

Identification of immune infiltration and cuproptosis-related molecular clusters in tuberculosis

Background

Tuberculosis (TB) is an infectious disease caused by Mycobacterium tuberculosis (Mtb) infection. Cuproptosis is a novel cell death mechanism correlated with various diseases. This study sought to elucidate the role of cuproptosis-related genes (CRGs) in TB.

Methods

Based on the GSE83456 dataset, we analyzed the expression profiles of CRGs and immune cell infiltration in TB. Based on CRGs, the molecular clusters and related immune cell infiltration were explored using 92 TB samples. The Weighted Gene Co-expression Network Analysis (WGCNA) algorithm was utilized to identify the co-expression modules and cluster-specific differentially expressed genes. Subsequently, the optimal machine learning model was determined by comparing the performance of the random forest (RF), support vector machine (SVM), generalized linear model (GLM), and eXtreme Gradient Boosting (XGB). The predictive performance of the machine learning model was assessed by generating calibration curves and decision curve analysis and validated in an external dataset.

Results

11 CRGs were identified as differentially expressed cuproptosis genes. Significant differences in immune cells were observed in TB patients. Two cuproptosis-related molecular clusters expressed genes were identified. Distinct clusters were identified based on the differential expression of CRGs and immune cells. Besides, significant differences in biological functions and pathway activities were observed between the two clusters. A nomogram was generated to facilitate clinical implementation. Next, calibration curves were generated, and decision curve analysis was conducted to validate the accuracy of our model in predicting TB subtypes. XGB machine learning model yielded the best performance in distinguishing TB patients with different clusters. The top five genes from the XGB model were selected as predictor genes. The XGB model exhibited satisfactory performance during validation in an external dataset. Further analysis revealed that these five model-related genes were significantly associated with latent and active TB.

Conclusion

Our study provided hitherto undocumented evidence of the relationship between cuproptosis and TB and established an optimal machine learning model to evaluate the TB subtypes and latent and active TB patients.

Interaction molecular QTL mapping discovers cellular and environmental modifiers of genetic regulatory effects

Bulk tissue molecular quantitative trait loci (QTLs) have been the starting point for interpreting disease-associated variants, while context-specific QTLs show particular relevance for disease. Here, we present the results of mapping interaction QTLs (iQTLs) for cell type, age, and other phenotypic variables in multi-omic, longitudinal data from blood of individuals of diverse ancestries. By modeling the interaction between genotype and estimated cell type proportions, we demonstrate that cell type iQTLs could be considered as proxies for cell type-specific QTL effects. The interpretation of age iQTLs, however, warrants caution as the moderation effect of age on the genotype and molecular phenotype association may be mediated by changes in cell type composition. Finally, we show that cell type iQTLs contribute to cell type-specific enrichment of diseases that, in combination with additional functional data, may guide future functional studies. Overall, this study highlights iQTLs to gain insights into the context-specificity of regulatory effects.

Comparative Analysis of Cell Mixtures Deconvolution and Gene Signatures Generated for Blood, Immune and Cancer Cells

In the last two decades, many detailed full transcriptomic studies on complex biological samples have been published and included in large gene expression repositories. These studies primarily provide a bulk expression signal for each sample, including multiple cell-types mixed within the global signal. The cellular heterogeneity in these mixtures does not allow the activity of specific genes in specific cell types to be identified. Therefore, inferring relative cellular composition is a very powerful tool to achieve a more accurate molecular profiling of complex biological samples. In recent decades, computational techniques have been developed to solve this problem by applying deconvolution methods, designed to decompose cell mixtures into their cellular components and calculate the relative proportions of these elements. Some of them only calculate the cell proportions (supervised methods), while other deconvolution algorithms can also identify the gene signatures specific for each cell type (unsupervised methods). In these work, five deconvolution methods (CIBERSORT, FARDEEP, DECONICA, LINSEED and ABIS) were implemented and used to analyze blood and immune cells, and also cancer cells, in complex mixture samples (using three bulk expression datasets). Our study provides three analytical tools (corrplots, cell-signature plots and bar-mixture plots) that allow a thorough comparative analysis of the cell mixture data. The work indicates that CIBERSORT is a robust method optimized for the identification of immune cell-types, but not as efficient in the identification of cancer cells. We also found that LINSEED is a very powerful unsupervised method that provides precise and specific gene signatures for each of the main immune cell types tested: neutrophils and monocytes (of the myeloid lineage), B-cells, NK cells and T-cells (of the lymphoid lineage), and also for cancer cells.

Peripheral blood cellular dynamics of rheumatoid arthritis treatment informs about efficacy of response to disease modifying drugs

Rheumatoid arthritis (RA) is an autoimmune disease characterized by systemic inflammation and is mediated by multiple immune cell types. In this work, we aimed to determine the relevance of changes in cell proportions in peripheral blood mononuclear cells (PBMCs) during the development of disease and following treatment. Samples from healthy blood donors, newly diagnosed RA patients, and established RA patients that had an inadequate response to MTX and were about to start tumor necrosis factor inhibitors (TNFi) treatment were collected before and after 3 months of treatment. We used in parallel a computational deconvolution approach based on RNA expression and flow cytometry to determine the relative cell-type frequencies. Cell-type frequencies from deconvolution of gene expression indicate that monocytes (both classical and non-classical) and CD4⁺ cells (T_h1 and T_h2) were increased in RA patients compared to controls, while NK cells and B cells (naïve and mature) were significantly decreased in RA patients. Treatment with MTX caused a decrease in B cells (memory and plasma cell), and a decrease in CD4 T_h cells (T_h1 and T_h17), while treatment with TNFi resulted in a significant increase in the population of B cells. Characterization of the RNA expression patterns found that most of the differentially expressed genes in RA subjects after treatment can be explained by changes in cell frequencies (98% and 74% respectively for MTX and TNFi).

Revealing the immune cell subtype reconstitution profile in patients from the CLARITY study using deconvolution algorithms after cladribine tablets treatment

Immune Cell Deconvolution methods utilizing gene expression profiling to quantify immune cells in tissues and blood are an appealing alternative to flow cytometry. Our objective was to investigate the applicability of deconvolution approaches in clinical trial settings to better investigate the mode of action of drugs for autoimmune diseases. Popular deconvolution methods CIBERSORT and xCell were validated using gene expression from the publicly available GSE93777 dataset that has comprehensive matching flow cytometry. As shown in the online tool, ~ 50% of signatures show strong correlation (r > 0.5) with the remainder showing moderate correlation, or in a few cases, no correlation. Deconvolution methods were then applied to gene expression data from the phase III CLARITY study (NCT00213135) to evaluate the immune cell profile of relapsing multiple sclerosis patients treated with cladribine tablets. At 96 weeks after treatment, deconvolution scores showed the following changes vs placebo: naïve, mature, memory CD4⁺ and CD8⁺ T cells, non-class switched, and class switched memory B cells and plasmablasts were significantly reduced, naïve B cells and M2 macrophages were more abundant. Results confirm previously described changes in immune cell composition following cladribine tablets treatment and reveal immune homeostasis of pro- vs anti-inflammatory immune cell subtypes, potentially supporting long-term efficacy.

Widespread sex-dimorphism across single-cell transcriptomes of adult African turquoise killifish tissues

The African turquoise killifish (Nothobranchius furzeri), the shortest-lived vertebrate that can be bred in captivity, is an emerging model organism to study vertebrate aging. Here we describe the first multi-tissue, single-cell gene expression atlas of female and male turquoise killifish tissues comprising immune and metabolic cells from the blood, kidney, liver, and spleen. We were able to annotate 22 distinct cell types, define associated marker genes, and infer differentiation trajectories. Using this dataset, we found pervasive sex-dimorphic gene expression across cell types, especially in the liver. Sex-dimorphic genes tended to be involved in processes related to lipid metabolism, and indeed, we observed clear differences in lipid storage in female vs. male turquoise killifish livers. Importantly, we use machine-learning to predict sex using single-cell gene expression in our atlas and identify potential transcriptional markers for molecular sex identity in this species. As proof-of-principle, we show that our atlas can be used to deconvolute existing liver bulk RNA-seq data in this species to obtain accurate estimates of cell type proportions across biological conditions. We believe that this single-cell atlas can be a resource to the community that could notably be leveraged to identify cell type-specific genes for cell type-specific expression in transgenic animals.

Qureshi A, Li J, Mahmood T. On the Analyses of Medical Images Using Traditional Machine Learning Techniques and Convolutional Neural Networks

Convolutional neural network (CNN) has shown dissuasive accomplishment on different areas especially Object Detection, Segmentation, Reconstruction (2D and 3D), Information Retrieval, Medical Image Registration, Multi-lingual translation, Local language Processing, Anomaly Detection on video and Speech Recognition. CNN is a special type of Neural Network, which has compelling and effective learning ability to learn features at several steps during augmentation of the data. Recently, different interesting and inspiring ideas of Deep Learning (DL) such as different activation functions, hyperparameter optimization, regularization, momentum and loss functions has improved the performance, operation and execution of CNN Different internal architecture innovation of CNN and different representational style of CNN has significantly improved the performance. This survey focuses on internal taxonomy of deep learning, different models of vonvolutional neural network, especially depth and width of models and in addition CNN components, applications and current challenges of deep learning.

Involvement of RNA methylation modification patterns mediated by m7G, m6A, m5C and m1A regulators in immune microenvironment regulation of Sjögren's syndrome

Keratoconjunctivitis is the most common complication of Sjögren's syndrome (SS). It has always been a hot research topic due to its complex pathogenesis. A further understanding of keratoconjunctiva xerosis can be obtained by studying the primary diseases. 7-Methylguanine (m7G), N6-methyladenosine (m6A), 5-methylcytosine (m5C), and N1-methyladenosine (m1A) are newly discovered epigenetic mechanisms involved in the development of SS. This study aimed to investigate the effects of m7G, m6A, m5C, and m1A modifications on the immune microenvironment of SS. Three microarray datasets were downloaded from the Gene Omnibus Expression (GEO) database, including 56 SS samples and 35 normal samples. Then, genes with m7G, m6A, m5C, and m1A methylation were explored, and the RNA modification patterns mediated by 59 m7G, m6A, m5C, and m1A regulators were summarized. The effects of m7G, m6A, m5C, and m1A modifications on immune infiltrating cells were discussed. Eukaryotic translation initiation factor 3 subunit D(EIF3D) was closely related to monocytes, and the expression of EIF3D was higher in SS with less monocytes. Two distinct patterns of RNA modification mediated by the 59 m7G, m6A, m5C, and m1A regulators were also identified, which infiltrated immune cells differently. Moreover, the two distinct RNA patterns were enriched in different signaling pathways, and their biological functions were explored. The findings revealed that m7G, m6A, m5C, and m1A modifications played vital roles in the diversity and complexity of the immune microenvironment in SS.

Placental cell type deconvolution reveals that cell proportions drive preeclampsia gene expression differences

The placenta mediates adverse pregnancy outcomes, including preeclampsia, which is characterized by gestational hypertension and proteinuria. Placental cell type heterogeneity in preeclampsia is not well-understood and limits mechanistic interpretation of bulk gene expression measures. We generated single-cell RNA-sequencing samples for integration with existing data to create the largest deconvolution reference of 19 fetal and 8 maternal cell types from placental villous tissue (n = 9 biological replicates) at term (n = 40,494 cells). We deconvoluted eight published microarray case-control studies of preeclampsia (n = 173 controls, 157 cases). Preeclampsia was associated with excess extravillous trophoblasts and fewer mesenchymal and Hofbauer cells. Adjustment for cellular composition reduced preeclampsia-associated differentially expressed genes (log2 fold-change cutoff = 0.1, FDR < 0.05) from 1154 to 0, whereas downregulation of mitochondrial biogenesis, aerobic respiration, and ribosome biogenesis were robust to cell type adjustment, suggesting direct changes to these pathways. Cellular composition mediated a substantial proportion of the association between preeclampsia and FLT1 (37.8%, 95% CI [27.5%, 48.8%]), LEP (34.5%, 95% CI [26.0%, 44.9%]), and ENG (34.5%, 95% CI [25.0%, 45.3%]) overexpression. Our findings indicate substantial placental cellular heterogeneity in preeclampsia contributes to previously observed bulk gene expression differences. This deconvolution reference lays the groundwork for cellular heterogeneity-aware investigation into placental dysfunction and adverse birth outcomes.

Elevated THBS3 predicts poor overall survival for clear cell renal cell carcinoma and identifies LncRNA/RBP/THBS3 mRNA networks

This study was used to assess THBS3's overall survival (OS) prognostic values in clear cell renal cell carcinoma (ccRCC) as well as to determine the LncRNA/RNA binding protein (RBP)/THBS3 interactions. Clinical data and RNA sequencing data were gathered from the TCGA dataset. Significant pathways associated with THBS3 were identified by gene set enrichment analysis (GSEA). Cox regression analyses, both univariate and multivariate, were applied to assess factors with independent prognostic abilities. We also discussed THBS3's relationship to immunity. We discovered that THBS3 expression was increased in ccRCC samples, as well as shorter OS in the TCGA dataset (P<0.05). External verification results in GSE6344, ICGC, ArrayExpress, UALCAN datasets, and qRT-PCR remained consistent (all P<0.05). Cox regression analyses, both univariate and multivariate, identified THBS3 as a factor with independent prognostic ability (both P<0.001). THBS3 expression as well as several clinicopathological variables were included in the nomogram OS prognosis prediction method as well. GSEA identified four THBS3-related signal pathways and THBS3 was revealed to be significantly associated with MSI, TMB, neoantigen, and immunity (all P<0.05). We also identified several LncRNA/RBP/THBS3 mRNA networks as potentially THBS3-related mechanisms. For THBS3-related drug sensitivities, THBS3 was negatively associated with Actinomycin D, Cobimetinib, Eribulin mesilate, Geldanamycin analog, and Vinblastine, while it was positively related to Erlotinib drug sensitivity. In addition to being an independent prognostic factor for ccRCC, THBS3 had a close connection to immunity, with identifying LncRNA/RBPs/THBS3 mRNA networks. Verifications of our findings in vivo and in vitro should be done in the future.

Bioinformatics Analysis of Immune Cell Infiltration and Diagnostic Biomarkers between Ankylosing Spondylitis and Inflammatory Bowel Disease

Background

Ankylosing spondylitis (AS) and inflammatory bowel disease (IBD) are both autoimmune diseases, and they often occur together in clinical practice, but the pathogenesis is unclear. This study is aimed at identifying the hub genes and explore the related immune molecular mechanisms between AS and IBD by bioinformatics analysis.

Methods

From the public Gene Expression Omnibus (GEO) database, the AS and IBD datasets (GSE73754, GSE59071, GSE25101, and GSE36807) were obtained. The immune cell infiltration in the peripheral blood tissues of GSE73754 and GSE59071 was assessed using the CIBERSORT algorithm. Then, we used the Weighted Gene Coexpression Network Analysis (WGCNA) to identify the Differentially Expressed Genes (DEGs) related to AS and IBD. Then, the immune genes from the ImmPort database intersected with the DEGs to obtain hub genes. The Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyzed the functional correlation of hub genes. Then, hub genes were verified in GSE25101 and GSE36807. The clusterProfiler software and Gene Set Enrichment Analysis (GSEA) were used to conduct functional enrichment and pathway enrichment studies. Finally, the diagnostic efficacy was assessed using Receiver Operating Characteristic (ROC) curve analysis.

Results

The analysis of immune characteristics showed that both AS and IBD were related to immunity, and neutrophils were positively correlated in both diseases. Nine coexpressed genes, including FCGRT, S100A11, IFNGR1, NFKBIZ, JAK2, LYN, PLAUR, ADM, and IL1RN, were linked to immune cells. The GO and KEGG analyses results showed that enrichment analysis was mainly related to cell transport and migration. Finally, the ROC curve was verified with the validation set, and it was found that PLAUR has clinical diagnostic significance and the most excellent specificity and sensitivity, respectively.

Conclusions

PLAUR (uPAR) is a promising biomarker and will be an underlying genetic biomarker for diagnosing AS comorbid IBD. Inflammation and immunological modulation mediated by neutrophil infiltration were important in the development of AS and IBD and may be diagnostic and therapeutic targets.

A single-nucleus transcriptomics study of alcohol use disorder in the nucleus accumbens

Gene expression studies offer promising opportunities to better understand the processes underlying alcohol use disorder (AUD). As cell types differ in their function, gene expression profiles will typically vary across cell types. When studying bulk tissue, failure to account for this cellular diversity has a detrimental impact on the ability to detect disease associations. We therefore assayed the transcriptomes of 32,531 individual nuclei extracted from the nucleus accumbens (NAc) of nine donors with AUD and nine controls (72% male). Our study identified 17 clearly delineated cell types. We detected 26 transcriptome-wide significant differentially expressed genes (DEGs) that mainly involved medium spiny neurons with both D1-type and D2-type dopamine receptors, microglia (MGL) and oligodendrocytes. A higher than expected number of DEGs replicated in an existing single nucleus gene expression study of alcohol dependence in the prefrontal cortex (enrichment ratio 1.91, p value 0.019) with two genes remaining significant after a Bonferroni correction. Our most compelling result involved CD53 in MGL that replicated in the same cell type in the prefrontal cortex and was previously implicated in studies of DNA methylation, bulk gene expression and genetic variants. Several DEGs were previously reported to be associated with AUD (e.g., PER1 and MGAT5). The DEGs for MSN.3 seemed involved in neurodegeneration, disruption of circadian rhythms, alterations in glucose metabolism and changes in synaptic plasticity. For MGL, the DEGs implicated neuroinflammation and immune-related processes and for OLI, disruptions in myelination. This identification of the specific cell-types from which the association signals originate will be key for designing proper follow-up experiments and, eventually, novel clinical interventions.

Role of immune cell infiltration and small molecule drugs in adhesive capsulitis: Novel exploration based on bioinformatics analyses

Background

Adhesive capsulitis (AC) is a type of arthritis that causes shoulder joint pain, stiffness, and limited mobility. The pathogenesis of AC is still controversial. This study aims to explore the role of immune related factors in the occurrence and development of AC.

Methods

The AC dataset was downloaded from Gene Expression Omnibus (GEO) data repository. Differentially expressed immune-related genes (DEIRGs) were obtained based on R package "DESeq2" and Immport database. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) were performed to explore the functional correlation of DEIRGs. MCC method and Least Absolute Shrinkage and Selection Operator (LASSO) regression were conducted to identify the hub genes. The immune cell infiltration in shoulder joint capsule between AC and control was evaluated by CIBERSORTx, and the relationship between hub genes and infiltrating immune cells was analyzed by Spearman's rank correlation. Finally, potential small molecule drugs for AC were screened by the Connectivity Map database (CMap) and further verified by molecular docking.

Results

A total of 137 DEIRGs and eight significantly different types of infiltrating immune cells (M0 macrophages, M1 macrophages, regulatory T cells, Tfh cells, monocytes, activated NK cells, memory resting CD4+T cells and resting dendritic cells) were screened between AC and control tissues. MMP9, FOS, SOCS3, and EGF were identified as potential targets for AC. MMP9 was negatively correlated with memory resting CD4+T cells and activated NK cells, but positively correlated with M0 macrophages. SOCS3 was positively correlated with M1 macrophages. FOS was positively correlated with M1 macrophages. EGF was positively correlated with monocytes. Additionally, dactolisib (ranked first) was identified as a potential small-molecule drug for the targeted therapy of AC.

Conclusions

This is the first study on immune cell infiltration analysis in AC, and these findings may provide a new idea for the diagnosis and treatment of AC.

Integrative transcriptomic analysis in human and mouse model of anaphylaxis identifies gene signatures associated with cell movement, migration and neuroinflammatory signalling

Background

Anaphylaxis is an acute life-threatening allergic reaction and a concern at a global level; therefore, further progress in understanding the underlying mechanisms and more effective strategies for diagnosis, prevention and management are needed.

Objective

We sought to identify the global architecture of blood transcriptomic features of anaphylaxis by integrating expression data from human patients and mouse model of anaphylaxis.

Methods

Bulk RNA-sequencings of peripheral whole blood were performed in: i) 14 emergency department (ED) patients with acute anaphylaxis, predominantly to Hymenoptera venom, ii) 11 patients with peanut allergy undergoing double-blind, placebo-controlled food challenge (DBPCFC) to peanut, iii) murine model of IgE-mediated anaphylaxis. Integrative characterisation of differential gene expression, immune cell-type-specific gene expression profiles, and functional and pathway analysis was undertaken.

Results

1023 genes were commonly and significantly dysregulated during anaphylaxis in ED and DBPCFC patients; of those genes, 29 were also dysregulated in the mouse model. Cell-type-specific gene expression profiles showed a rapid downregulation of blood basophil and upregulation of neutrophil signature in ED and DBPCFC patients and the mouse model, but no consistent and/or significant differences were found for other blood cells. Functional and pathway analysis demonstrated that human and mouse blood transcriptomic signatures of anaphylaxis follow trajectories of upregulation of cell movement, migration and neuroinflammatory signalling, and downregulation of lipid activating nuclear receptors signalling.

Conclusion

Our study highlights the matched and extensive blood transcriptomic changes and suggests the involvement of discrete cellular components and upregulation of migration and neuroinflammatory pathways during anaphylaxis.

A systematic assessment of cell type deconvolution algorithms for DNA methylation data

We performed systematic assessment of computational deconvolution methods that play an important role in the estimation of cell type proportions from bulk methylation data. The proposed framework methylDeConv (available as an R package) integrates several deconvolution methods for methylation profiles (Illumina HumanMethylation450 and MethylationEPIC arrays) and offers different cell-type-specific CpG selection to construct the extended reference library which incorporates the main immune cell subsets, epithelial cells and cell-free DNAs. We compared the performance of different deconvolution algorithms via simulations and benchmark datasets and further investigated the associations of the estimated cell type proportions to cancer therapy in breast cancer and subtypes in melanoma methylation case studies. Our results indicated that the deconvolution based on the extended reference library is critical to obtain accurate estimates of cell proportions in non-blood tissues.

GEOMETRIC STRUCTURE GUIDED MODEL AND ALGORITHMS FOR COMPLETE DECONVOLUTION OF GENE EXPRESSION DATA

Complete deconvolution analysis for bulk RNA-seq data is important and helpful to distinguish whether the differences of disease-associated GEPs (gene expression profiles) in tissues of patients and normal controls are due to changes in cellular composition of tissue samples, or due to GEPs changes in specific cells. One of the major techniques to perform complete deconvolution is nonnegative matrix factorization (NMF), which also has a wide-range of applications in the machine learning community. However, the NMF is a well-known strongly ill-posed problem, so a direct application of NMF to RNA-seq data will suffer severe difficulties in the interpretability of solutions. In this paper, we develop an NMF-based mathematical model and corresponding computational algorithms to improve the solution identifiability of deconvoluting bulk RNA-seq data. In our approach, we combine the biological concept of marker genes with the solvability conditions of the NMF theories, and develop a geometric structures guided optimization model. In this strategy, the geometric structure of bulk tissue data is first explored by the spectral clustering technique. Then, the identified information of marker genes is integrated as solvability constraints, while the overall correlation graph is used as manifold regularization. Both synthetic and biological data are used to validate the proposed model and algorithms, from which solution interpretability and accuracy are significantly improved.

DNA methylation signatures of childhood trauma predict psychiatric disorders and other adverse outcomes 17 years after exposure

Childhood trauma is robustly linked to a broad range of adverse outcomes with consequences persisting far into adulthood. We conducted a prospective longitudinal study to predict psychiatric disorders and other adverse outcomes from trauma-related methylation changes 16.9 years after trauma exposure in childhood. Methylation was assayed using a sequencing-based approach that provides near-complete coverage of all 28 million sites in the blood methylome. Methylation data involved 673 assays from 489 participants aged 13.6 years (SD = 1.9) with outcomes measures collected at age 30.4 (SD = 2.26). For a subset of 303 participants we also generated methylation data in adulthood. Trauma-related methylation risk scores (MRSs) significantly predicted adult depression, externalizing problems, nicotine dependence, alcohol use disorder, serious medical problems, social problems and poverty. The predictive power of the MRSs was higher than that of reported trauma and could not be explained by the reported trauma, correlations with demographic variables, or a continuity of the predicted health problems from childhood to adulthood. Rather than measuring the occurrence of traumatic events, the MRSs seemed to capture the subject-specific impact of trauma. The majority of predictive sites did not remain associated with the outcomes suggesting the signatures of trauma do not become biologically embedded in the blood methylome. Instead, the long-term effects of trauma therefore seemed more consistent with a developmental mechanism where the initial subject-specific impacts of trauma are magnified over time. The MRSs have the potential to be a novel clinical biomarker for the assessment of trauma-related health risks.

Mature neutrophils and a NF-κB-to-IFN transition determine the unifying disease recovery dynamics in COVID-19

Disease recovery dynamics are often difficult to assess, as patients display heterogeneous recovery courses. To model recovery dynamics, exemplified by severe COVID-19, we apply a computational scheme on longitudinally sampled blood transcriptomes, generating recovery states, which we then link to cellular and molecular mechanisms, presenting a framework for studying the kinetics of recovery compared with non-recovery over time and long-term effects of the disease. Specifically, a decrease in mature neutrophils is the strongest cellular effect during recovery, with direct implications on disease outcome. Furthermore, we present strong indications for global regulatory changes in gene programs, decoupled from cell compositional changes, including an early rise in T cell activation and differentiation, resulting in immune rebalancing between interferon and NF-κB activity and restoration of cell homeostasis. Overall, we present a clinically relevant computational framework for modeling disease recovery, paving the way for future studies of the recovery dynamics in other diseases and tissues.

CPDR: An R Package of Recommending Personalized Drugs for Cancer Patients by Reversing the Individual’s Disease-Related Signature

Due to cancer heterogeneity, only some patients can benefit from drug therapy. The personalized drug usage is important for improving the treatment response rate of cancer patients. The value of the transcriptome of patients has been recently demonstrated in guiding personalized drug use, and the Connectivity Map (CMAP) is a reliable computational approach for drug recommendation. However, there is still no personalized drug recommendation tool based on transcriptomic profiles of patients and CMAP. To fill this gap, here, we proposed such a feasible workflow and a user-friendly R package—Cancer-Personalized Drug Recommendation (CPDR). CPDR has three features. 1) It identifies the individual disease signature by using the patient subgroup with transcriptomic profiles similar to those of the input patient. 2) Transcriptomic profile purification is supported for the subgroup with high infiltration of non-cancerous cells. 3) It supports in silico drug efficacy assessment using drug sensitivity data on cancer cell lines. We demonstrated the workflow of CPDR with the aid of a colorectal cancer dataset from GEO and performed the in silico validation of drug efficacy. We further assessed the performance of CPDR by a pancreatic cancer dataset with clinical response to gemcitabine. The results showed that CPDR can recommend promising therapeutic agents for the individual patient. The CPDR R package is available at https://github.com/AllenSpike/CPDR.

Transcriptome-wide association study for postpartum depression implicates altered B-cell activation and insulin resistance

Postpartum depression (PPD) affects 1 in 7 women and has negative mental health consequences for both mother and child. However, the precise biological mechanisms behind the disorder are unknown. Therefore, we performed the largest transcriptome-wide association study (TWAS) for PPD (482 cases, 859 controls) to date using RNA-sequencing in whole blood and deconvoluted cell types. No transcriptional changes were observed in whole blood. B-cells showed a majority of transcriptome-wide significant results (891 transcripts representing 789 genes) with pathway analyses implicating altered B-cell activation and insulin resistance. Integration of other data types revealed cell type-specific DNA methylation loci and disease-associated eQTLs (deQTLs), but not hormones/neuropeptides (estradiol, progesterone, oxytocin, BDNF), serve as regulators for part of the transcriptional differences between cases and controls. Further, deQTLs were enriched for several brain region-specific eQTLs, but no overlap with MDD risk loci was observed. Altogether, our results constitute a convergence of evidence for pathways most affected in PPD with data across different biological mechanisms.

A Hypoxia-Related Signature for Predicting Prognosis, Cellular Processes, Immune Microenvironment and Targeted Compounds in Lung Squamous Cell Carcinoma

Background

Lung squamous cell carcinoma (LUSC) is a malignant tumour of the lung epithelium. A hypoxic environment can promote tumour cell proliferation and invasion. Therefore, this study aims to explore hypoxia-related genes and construct reliable models to predict the prognosis, cellular processes, immune microenvironment and target compounds of lung squamous carcinoma.

Methods

The transcriptome data and matched clinical information of LUSC were retrieved from The Cancer Genome Atlas (TCGA) database. The GSVA algorithm calculated each LUSC patient’s hypoxia score, and all LUSC patients were divided into the high hypoxia score group and low hypoxia score group. Weighted gene co-expression network analysis (WGCNA) and differential expression analysis were performed to screen out differentially expressed hypoxia-related genes (DE-HRGs) in LUSC microenvironment, and the underlying regulatory mechanism of DE-HRGs in LUSC was explored through Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses. Hereafter, we established a prognosis-related genetic signature for DE-HRGs using univariate and multivariate Cox regression analyses. The relationship between gene signature and immune cells was further evaluated. Finally, the Comparative Toxicogenomics Database (CTD) was utilized to predict the targeted drugs for the prognostic genes.

Results

We obtained 376 DE-HRGs. Functional enrichment analysis indicated that the DE-HRGs were involved in the cell cycle-related regulatory processes. Next, we developed and validated 3 HRGs-based prognostic signature for LUSC, including HELLS, GPRIN1, and FAM83A. Risk score is an independent prognostic factor for LUSC. Functional enrichment analysis and immune landscape analysis suggested that the risk scoring system might be involved in altering the immune microenvironment of LUSC patients to influence patient outcomes. Ultimately, a total of 92 potential compounds were predicted for the three prognostic genes.

Conclusion

In summary, we developed and validated a hypoxia-related model for LUSC, reflecting the cellular processes and immune microenvironment characteristics and predicting the prognostic outcomes and targeted compounds.

Dual methylation and hydroxymethylation study of alcohol use disorder

Using an integrative, multi-tissue design, we sought to characterize methylation and hydroxymethylation changes in blood and brain associated with alcohol use disorder (AUD). First, we used epigenomic deconvolution to perform cell-type-specific methylome-wide association studies within subpopulations of granulocytes/T-cells/B-cells/monocytes in 1132 blood samples. Blood findings were then examined for overlap with AUD-related associations with methylation and hydroxymethylation in 50 human post-mortem brain samples. Follow-up analyses investigated if overlapping findings mediated AUD-associated transcription changes in the same brain samples. Lastly, we replicated our blood findings in an independent sample of 412 individuals and aimed to replicate published alcohol methylation findings using our results. Cell-type-specific analyses in blood identified methylome-wide significant associations in monocytes and T-cells. The monocyte findings were significantly enriched for AUD-related methylation and hydroxymethylation in brain. Hydroxymethylation in specific sites mediated AUD-associated transcription in the same brain samples. As part of the most comprehensive methylation study of AUD to date, this work involved the first cell-type-specific methylation study of AUD conducted in blood, identifying and replicating a finding in DLGAP1 that may be a blood-based biomarker of AUD. In this first study to consider the role of hydroxymethylation in AUD, we found evidence for a novel mechanism for cognitive deficits associated with AUD. Our results suggest promising new avenues for AUD research.

Contribution of classification based on ferroptosis-related genes to the heterogeneity of MAFLD

Background

Metabolic dysfunction-associated fatty liver disease (MAFLD) is a highly heterogeneous disease and its heterogeneity might be associated with ferroptosis because ferroptosis plays an important role in the development of MAFLD. We aimed to perform integrative analysis of ferroptosis related genes and MAFLD subtypes using bioinformatics.

Methods

A differential expression analysis was performed to identify key ferroptosis-related genes associated with the clinical characteristics of MAFLD. Furthermore, consensus k clustering was utilized to distinguish ferroptosis-related clinical subtypes of MAFLD and assess the association of ferroptosis-related gene expression and clinical features between patients with different subtypes of MAFLD. Moreover, the variation in the immune status and regulatory relationship of ferroptosis-related genes in individuals with MAFLD was also explored using single sample gene set enrichment analysis, weighted gene coexpression network analysis and enrichment analyses.

Results

Eight ferroptosis-related genes were identified as closely associated with both the hepatic steatosis grade and non-alcoholic fatty liver disease activity score. Two subtypes of MAFLD based on ferroptosis-related genes were identified by consensus clustering. They exhibited significantly different clinical features, immune statuses, biological processes and outcomes. The progression of the two subtypes was associated with immunity.

Conclusions

Two highly heterogeneous subtypes of MAFLD with significantly distinct clinical features, biological processes and immune statuses were identified based on ferroptosis-associated genes, which strongly supports the hypothesis that ferroptosis plays an important role in the development of MAFLD.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12876-022-02137-9.

Identification of CCL19 as a Novel Immune-Related Biomarker in Diabetic Nephropathy

Diabetic nephropathy (DN) is one of the major microvascular complications in diabetic patients and the leading cause of end-stage renal disease (ESRD). Previous studies found that immune-related genes and immune cell infiltration play important roles in the pathogenesis and development of DN. Therefore, this study aimed to explore immune-related biomarkers in DN. In this research, three microarray datasets that included 18 DN and 28 healthy tubule samples were downloaded and integrated as the training set to identify differentially expressed immune-related genes (DEIGs). A total of 63 DEIGs were identified, and most upregulated DEIGs were primarily involved in the inflammatory response and chemokine-mediated signaling pathways. The Microenvironment Cell Populations-counter (MCP-counter) algorithm was then used to estimate the abundance of infiltrated immune and stromal cell populations. According to DEIG, weighted gene coexpression network and protein–protein network analyses, CCL19 was identified as the hub immune-related biomarker. Moreover, the upregulated level of CCL19 was confirmed in other independent datasets as well as in in vitro experiments with high glucose. In summary, this study provides novel insights into the pathogenesis of diabetic nephropathy and identifies CCL19 as a potential critical gene of DN.