eQTL studies: from bulk tissues to single cells

Optimal variable identification for accurate detection of causal expression Quantitative Trait Loci with applications in heart-related diseases

Gene expression plays a pivotal role in various diseases, contributing significantly to their mechanisms. Most GWAS risk loci are in non-coding regions, potentially affecting disease risk by altering gene expression in specific tissues. This expression is notably tissue-specific, with genetic variants substantially influencing it. However, accurately detecting the expression Quantitative Trait Loci (eQTL) is challenging due to limited heritability in gene expression, extensive linkage disequilibrium (LD), and multiple causal variants. The single variant association approach in eQTL analysis is limited by its susceptibility to capture the combined effects of multiple variants, and a bias towards common variants, underscoring the need for a more robust method to accurately identify causal eQTL variants. To address this, we developed an algorithm, CausalEQTL, which integrates L 0 +L 1 penalized regression with an ensemble approach to localize eQTL, thereby enhancing prediction performance precisely. Our results demonstrate that CausalEQTL outperforms traditional models, including LASSO, Elastic Net, Ridge, in terms of power and overall performance. Furthermore, analysis of heart tissue data from the GTEx project revealed that eQTL sites identified by our algorithm provide deeper insights into heart-related tissue eQTL detection. This advancement in eQTL mapping promises to improve our understanding of the genetic basis of tissue-specific gene expression and its implications in disease. The source code and identified causal eQTLs for CausalEQTL are available on GitHub: https://github.com/zhc-moushang/CausalEQTL.

Multi-Omics Analysis Revealed the rSNPs Potentially Involved in T2DM Pathogenic Mechanism and Metformin Response

The goal of our study was to identify and assess the functionally significant SNPs with potentially important roles in the development of type 2 diabetes mellitus (T2DM) and/or their effect on individual response to antihyperglycemic medication with metformin. We applied a bioinformatics approach to identify the regulatory SNPs (rSNPs) associated with allele-asymmetric binding and expression events in our paired ChIP-seq and RNA-seq data for peripheral blood mononuclear cells (PBMCs) of nine healthy individuals. The rSNP outcomes were analyzed using public data from the GWAS (Genome-Wide Association Studies) and Genotype-Tissue Expression (GTEx). The differentially expressed genes (DEGs) between healthy and T2DM individuals (GSE221521), including metformin responders and non-responders (GSE153315), were searched for in GEO RNA-seq data. The DEGs harboring rSNPs were analyzed using the Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG). We identified 14,796 rSNPs in the promoters of 5132 genes of human PBMCs. We found 4280 rSNPs to associate with both phenotypic traits (GWAS) and expression quantitative trait loci (eQTLs) from GTEx. Between T2DM patients and controls, 3810 rSNPs were detected in the promoters of 1284 DEGs. Based on the protein-protein interaction (PPI) network, we identified 31 upregulated hub genes, including the genes involved in inflammation, obesity, and insulin resistance. The top-ranked 10 enriched KEGG pathways for these hubs included insulin, AMPK, and FoxO signaling pathways. Between metformin responders and non-responders, 367 rSNPs were found in the promoters of 131 DEGs. Genes encoding transcription factors and transcription regulators were the most widely represented group and many were shown to be involved in the T2DM pathogenesis. We have formed a list of human rSNPs that add functional interpretation to the T2DM-association signals identified in GWAS. The results suggest candidate causal regulatory variants for T2DM, with strong enrichment in the pathways related to glucose metabolism, inflammation, and the effects of metformin.

Exploring the causal association between uric acid and lung cancer in east Asian and European populations: a mendelian randomization study

BACKGROUND

Lung cancer still ranks first in the mortality rate of cancer. Uric acid is a product of purine metabolism in humans. Its presence in the serum is controversial; some say that its high levels have a protective effect against tumors, others say the opposite, that is, high levels increase the risk of cancer. Therefore, the aim of this study was to investigate the potential causal association between serum uric acid levels and lung cancer.

METHODS

Mendelian randomization was used to achieve our aim. Sensitivity analyses was performed to validate the reliability of the results, followed by reverse Mendelian analyses to determine a potential reverse causal association.

RESULTS

A significant causal association was found between serum uric acid levels and lung cancer in East Asian and European populations. Further sublayer analysis revealed a significant causal association between uric acid and small cell lung cancer, while no potential association was observed between uric acid and non-small cell lung cancer, squamous lung cancer, and lung adenocarcinoma. The sensitivity analyses confirmed the reliability of the results. Reverse Mendelian analysis showed no reverse causal association between uric acid and lung cancer.

CONCLUSIONS

The results of this study suggested that serum uric acid levels were negatively associated with lung cancer, with uric acid being a potential protective factor for lung cancer. In addition, uric acid level monitoring was simple and inexpensive. Therefore, it might be used as a biomarker for lung cancer, promoting its wide use clinical practice.

Functional Genomics in Psoriasis

Psoriasis is an autoimmune cutaneous condition that significantly impacts quality of life and represents a burden on society due to its prevalence. Genome-wide association studies (GWASs) have pinpointed several psoriasis-related risk loci, underlining the disease's complexity. Functional genomics is paramount to unveiling the role of such loci in psoriasis and disentangling its complex nature. In this review, we aim to elucidate the main findings in this field and integrate our discussion with gold-standard techniques in molecular biology-i.e., Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR)-and high-throughput technologies. These tools are vital to understanding how disease risk loci affect gene expression in psoriasis, which is crucial in identifying new targets for personalized treatments in advanced precision medicine.

Integrated analysis of single-cell RNA-seq, bulk RNA-seq, Mendelian randomization, and eQTL reveals T cell-related nomogram model and subtype classification in rheumatoid arthritis

Objective

Rheumatoid arthritis (RA) is a systemic disease that attacks the joints and causes a heavy economic burden on humans worldwide. T cells regulate RA progression and are considered crucial targets for therapy. Therefore, we aimed to integrate multiple datasets to explore the mechanisms of RA. Moreover, we established a T cell-related diagnostic model to provide a new method for RA immunotherapy.

Methods

scRNA-seq and bulk-seq datasets for RA were obtained from the Gene Expression Omnibus (GEO) database. Various methods were used to analyze and characterize the T cell heterogeneity of RA. Using Mendelian randomization (MR) and expression quantitative trait loci (eQTL), we screened for potential pathogenic T cell marker genes in RA. Subsequently, we selected an optimal machine learning approach by comparing the nine types of machine learning in predicting RA to identify T cell-related diagnostic features to construct a nomogram model. Patients with RA were divided into different T cell-related clusters using the consensus clustering method. Finally, we performed immune cell infiltration and clinical correlation analyses of T cell-related diagnostic features.

Results

By analyzing the scRNA-seq dataset, we obtained 10,211 cells that were annotated into 7 different subtypes based on specific marker genes. By integrating the eQTL from blood and RA GWAS, combined with XGB machine learning, we identified a total of 8 T cell-related diagnostic features (MIER1, PPP1CB, ICOS, GADD45A, CD3D, SLFN5, PIP4K2A, and IL6ST). Consensus clustering analysis showed that RA could be classified into two different T-cell patterns (Cluster 1 and Cluster 2), with Cluster 2 having a higher T-cell score than Cluster 1. The two clusters involved different pathways and had different immune cell infiltration states. There was no difference in age or sex between the two different T cell patterns. In addition, ICOS and IL6ST were negatively correlated with age in RA patients.

Conclusion

Our findings elucidate the heterogeneity of T cells in RA and the communication role of these cells in an RA immune microenvironment. The construction of T cell-related diagnostic models provides a resource for guiding RA immunotherapeutic strategies.

The effect of metabolism-related lifestyle and clinical risk factors on digestive system cancers in East Asian populations: a two-sample Mendelian randomization analysis

Metabolic factors play a critical role in the development of digestive system cancers (DSCs), and East Asia has the highest incidence of malignant tumors in the digestive system. We performed a two-sample Mendelian randomization analysis to explore the associations between 19 metabolism-related lifestyle and clinical risk factors and DSCs, including esophageal, gastric, colorectal, hepatocellular, biliary tract, and pancreatic cancer. The causal association was explored for all combinations of each risk factor and each DSC. We gathered information on the instrumental variables (IVs) from various sources and retrieved outcome information from Biobank Japan (BBJ). The data were all from studies of east Asian populations. Finally, 17,572 DSCs cases and 195,745 controls were included. Our analysis found that genetically predicted alcohol drinking was a strong indicator of gastric cancer (odds ratio (OR) = 0.95; 95% confidence interval (CI): 0.93-0.98) and hepatocellular carcinoma (OR = 1.11; 95% CI: 1.05-1.18), whereas coffee consumption had a potential protective effect on hepatocellular carcinoma (OR = 0.69; 95% CI: 0.53-0.90). Triglyceride was potentially associated with a decreased risk of biliary tract cancer (OR = 0.53; 95% CI: 0.34-0.81), and uric acid was associated with pancreatic cancer risk (OR = 0.59; 95% CI: 0.37-0.96). Metabolic syndrome (MetS) was associated with esophageal and gastric cancer. Additionally, there was no evidence for a causal association between other risk factors, including body mass index, waist circumference, waist-to-hip ratio, educational levels, lipoprotein cholesterol, total cholesterol, glycine, creatinine, gout, and Graves' disease, and DSCs. The leave-one-out analysis revealed that the single nucleotide polymorphism (SNP) rs671 from the ALDH2 gene has a disproportionately high contribution to the causal association between alcohol drinking and gastric cancer and hepatocellular carcinoma, as well as the association between coffee consumption and hepatocellular carcinoma. The present study revealed multiple metabolism-related lifestyle and clinical risk factors and a valuable SNP rs671 for DSCs, highlighting the significance of metabolic factors in both the prevention and treatment of DSCs.

Cell-type deconvolution of bulk-blood RNA-seq reveals biological insights into neuropsychiatric disorders

Genome-wide association studies (GWASs) have uncovered susceptibility loci associated with psychiatric disorders such as bipolar disorder (BP) and schizophrenia (SCZ). However, most of these loci are in non-coding regions of the genome, and the causal mechanisms of the link between genetic variation and disease risk is unknown. Expression quantitative trait locus (eQTL) analysis of bulk tissue is a common approach used for deciphering underlying mechanisms, although this can obscure cell-type-specific signals and thus mask trait-relevant mechanisms. Although single-cell sequencing can be prohibitively expensive in large cohorts, computationally inferred cell-type proportions and cell-type gene expression estimates have the potential to overcome these problems and advance mechanistic studies. Using bulk RNA-seq from 1,730 samples derived from whole blood in a cohort ascertained from individuals with BP and SCZ, this study estimated cell-type proportions and their relation with disease status and medication. For each cell type, we found between 2,875 and 4,629 eGenes (genes with an associated eQTL), including 1,211 that are not found on the basis of bulk expression alone. We performed a colocalization test between cell-type eQTLs and various traits and identified hundreds of associations that occur between cell-type eQTLs and GWASs but that are not detected in bulk eQTLs. Finally, we investigated the effects of lithium use on the regulation of cell-type expression loci and found examples of genes that are differentially regulated according to lithium use. Our study suggests that applying computational methods to large bulk RNA-seq datasets of non-brain tissue can identify disease-relevant, cell-type-specific biology of psychiatric disorders and psychiatric medication.

Cell type deconvolution of bulk blood RNA-Seq to reveal biological insights of neuropsychiatric disorders

Genome-wide association studies (GWAS) have uncovered susceptibility loci associated with psychiatric disorders like bipolar disorder (BP) and schizophrenia (SCZ). However, most of these loci are in non-coding regions of the genome with unknown causal mechanisms of the link between genetic variation and disease risk. Expression quantitative trait loci (eQTL) analysis of bulk tissue is a common approach to decipher underlying mechanisms, though this can obscure cell-type specific signals thus masking trait-relevant mechanisms. While single-cell sequencing can be prohibitively expensive in large cohorts, computationally inferred cell type proportions and cell type gene expression estimates have the potential to overcome these problems and advance mechanistic studies. Using bulk RNA-Seq from 1,730 samples derived from whole blood in a cohort ascertained for individuals with BP and SCZ this study estimated cell type proportions and their relation with disease status and medication. We found between 2,875 and 4,629 eGenes for each cell type, including 1,211 eGenes that are not found using bulk expression alone. We performed a colocalization test between cell type eQTLs and various traits and identified hundreds of associations between cell type eQTLs and GWAS loci that are not detected in bulk eQTLs. Finally, we investigated the effects of lithium use on cell type expression regulation and found examples of genes that are differentially regulated dependent on lithium use. Our study suggests that computational methods can be applied to large bulk RNA-Seq datasets of non-brain tissue to identify disease-relevant, cell type specific biology of psychiatric disorders and psychiatric medication.