Expression Quantitative Trait Loci Research Articles

A comprehensive assessment of the association between common drugs and psychiatric disorders using Mendelian randomization and real-world pharmacovigilance database

Medications prescribed for chronic diseases can lead to short-term neuropsychiatric symptoms, but their long-term effects on brain structures and psychiatric conditions remain unclear. We comprehensively analyzed the FDA Adverse Event Reporting System database and conducted drug target Mendelian Randomization (MR) studies on six categories of common drugs, 477 brain imaging-derived phenotypes (IDPs) and eight psychiatric disorders. Genetic instruments were extracted from expression quantitative trait loci (eQTLs) in blood, brain, and other target tissues, protein quantitative trait loci (pQTLs) in blood, and genome-wide association studies (GWAS) of hemoglobin and cholesterol. Summary statistics for brain IDPs, psychiatric disorders, and gut microbiome were obtained from the BIG40, Psychiatric Genomics Consortium, and MiBioGen. A two-step MR and mediation analysis were employed to screen possible mediators of drug-IDP effects from 119 gut microbiota genera and identify their mediation proportions. Among 19 drug classes, six drugs were found to be associated with higher risks of psychiatric adverse events, while 11 drugs were associated with higher risks of gastrointestinal adverse events in the FAERS analysis. We identified ten drug-psychiatric disorder associations, 202 drug-IDP associations, 16 drug-microbiota associations, and four drug-microbiota-IDP causal links. For example, PPARG activation mediated HbA1c reduction caused a higher risk of bipolar disorder (BD) II. Genetically proxied GLP-1R agonists were significantly associated with an increase in the volume of the CA3-head of the right hippocampus and the area of the left precuneus cortex, both of which have been shown to correlate with cognition in previous studies. Common drugs may affect brain structure and risk of psychiatric disorder. Oral medications in particular may exert some of these effects by influencing gut microbiota. This study calls for greater attention to be paid to the neuropsychiatric adverse effects of drugs and encourages drug repurposing. National Natural Science Foundation of China (grant No. 82330035, 82130043, 82172685, and 82001223), National Natural Science Foundation of Hunan Province (grant No. 2021SK1010), and the Science Foundation for Distinguished Young Scholars of Changsha (grant No. kq2209006).

Mendelian randomization reveals potential causal relationships between cellular senescence-related genes and multiple cancer risks.

Cellular senescence is widely acknowledged as having strong associations with cancer. However, the intricate relationships between cellular senescence-related (CSR) genes and cancer risk remain poorly explored, with insights on causality remaining elusive. In this study, Mendelian Randomization (MR) analyses were used to draw causal inferences from 866 CSR genes as exposures and summary statistics for 18 common cancers as outcomes. We focused on genetic variants affecting gene expression, DNA methylation, and protein expression quantitative trait loci (cis-eQTL, cis-mQTL, and cis-pQTL, respectively), which were strongly linked to CSR genes alterations. Variants were selected as instrumental variables (IVs) and analyzed for causality with cancer using both summary-data-based MR (SMR) and two-sample MR (TSMR) approaches. Bayesian colocalization was used to unravel potential regulatory mechanisms underpinning risk variants in cancer, and further validate the robustness of MR results. We identified five CSR genes (CNOT6, DNMT3B, MAP2K1, TBPL1, and SREBF1), 18 DNA methylation genes, and LAYN protein expression which were all causally associated with different cancer types. Beyond causality, a comprehensive analysis of gene function, pathways, and druggability values was also conducted. These findings provide a robust foundation for unravelling CSR genes molecular mechanisms and promoting clinical drug development for cancer.

Identification of genetic variants associated with clinical features of sickle cell disease.

Sickle cell disease (SCD) is an inherited blood disorder marked by homozygosity of hemoglobin S, which is a defective hemoglobin caused by a missense mutation in the β-globin gene. However, clinical phenotypes of SCD vary among patients. To investigate genetic variants associated with various clinical phenotypes of SCD, we genotyped DNA samples from 520 SCD subjects and used a genome-wide association study (GWAS) approach to identify genetic variants associated with phenotypic features of SCD. For HbF levels, the previously reported 2p16.1 locus (BCL11A) reached genome significance (rs1427407, P = 8.58 × 10-10) in our GWAS as expected. In addition, we found a new genome-wide significance locus at 15q14 (rs8182015, P = 2.07 × 10-8) near gene EMC7. GWAS of acute chest syndrome (ACS) detected a locus (rs79915189, P = 3.70 × 10-8) near gene IDH2 at 15q26.1. The SNP, rs79915189, is also an expression quantitative trait locus (eQTL) of IDH2 in multiple tissues. For vasoocclusive episode (VOE), GWAS detected multiple significant signals at 2p25.1 (rs62118798, P = 4.27 × 10-8), 15q26.1 (rs62020555, P = 2.04 × 10-9) and 15q26.3 (rs117797325, P = 4.63 × 10-8). Our findings provide novel insights into the genetic mechanisms of SCD suggesting that common genetic variants play an important role in the presentation of the clinical phenotypes of patients with SCD.

Novel genetic structures associated with adverse response to chemotherapy in breast cancer.

Introduction: The role of genetic variants in response to chemotherapy has been investigated in several studies. This study aimed to investigate genetic variants associated with response to chemotherapy in breast cancer (BC) patients. Methods: Significant variants (p < 5 × 10-8) associated with response to chemotherapy were obtained from GWA studies. Candidate variants were identified by haplotype analysis (r2 ≥ 0.9, D'≥0.9) using 1000Genome LD data. To determine the effects of the variants on gene expression, expression quantitative trait loci (eQTL) were evaluated. To compare the expression of the identified genes in tumor samples, expression levels were compared between TCGA tumor types and adjacent normal tissues. Results: Six rs3820706, rs147451859, rs4784750, rs17587029, rs16830728, and rs16972207 variants were significantly associated with response to chemotherapy in BC patients (p < 5 × 10-8). Seven novel haplotypic structures were identified to be associated with adverse response to chemotherapy in BC patients. These haplotypes formed two genetic structures associated with neutropenia, leukopenia, chemotherapy-induced cytotoxicity (GAG-TTAT), and chemotherapy-induced alopecia (CC-CAACTCCCGTTGCGG). These variants are located on PPCDC, NLRC5, STAM2, and TNFSF13B genes, and the expression of these genes significantly changed in BC tissues than normal tissues (P ≤ 0.05), also showing gene-gene correlation (P ≤ 0.05). Conclusions: These genetic variants and their associated novel haplotypic structures can predict adverse response to chemotherapy in BC patients and could potentially form BC-associated genetic panel for adverse response to chemotherapy.

Extensive co-regulation of neighbouring genes complicates the use of eQTLs in target gene prioritisation

Identifying causal genes underlying genome-wide association studies (GWAS) is a fundamental problem in human genetics. Although colocalisation with gene expression quantitative trait loci (eQTLs) is often used to prioritise GWAS target genes, systematic benchmarking has been limited due to unavailability of large ground truth datasets. Here, we re-analysed plasma protein QTL data from 3,301 individuals of the INTERVAL cohort together with 131 eQTL Catalogue datasets. Focusing on variants located within or close to the affected protein identified 793 proteins with at least one cis-pQTL where we could assume that the most likely causal gene was the gene coding for the protein. We then benchmarked the ability of cis-eQTLs to recover these causal genes by comparing three Bayesian colocalisation methods (coloc.susie, coloc.abf and CLPP) and five Mendelian randomisation (MR) approaches (three varieties of inverse-variance weighted MR, MR-RAPS, and MRLocus). We found that assigning fine-mapped pQTLs to their closest protein coding genes outperformed all colocalisation methods regarding both precision (71.9%) and recall (76.9%). Furthermore, the colocalisation method with the highest recall (coloc.susie - 46.3%) also had the lowest precision (45.1%). Combining evidence from multiple conditionally distinct colocalising QTLs with MR increased precision to 81%, but this was accompanied by a large reduction in recall to 7.1%. Furthermore, the choice of the MR method greatly affected performance, with the standard inverse-variance weighted MR often producing many false positives. Our results highlight that linking GWAS variants to target genes remains challenging with eQTL evidence alone, and prioritising novel targets requires triangulation of evidence from multiple sources.

Genetic Associations With an Amyotrophic Lateral Sclerosis Reversal Phenotype.

The term "ALS Reversal" describes patients who initially meet diagnostic criteria for amyotrophic lateral sclerosis (ALS) or had clinical features most consistent with progressive muscular atrophy (PMA) but subsequently demonstrated substantial and sustained clinical improvement. The objective of this genome-wide association study (GWAS) was to identify correlates of this unusual clinical phenotype. Participants were recruited from a previously created database of individuals with the ALS Reversal phenotype. Whole-genome sequencing (WGS) data were compared with ethnicity-matched patients with typically progressive ALS enrolled through the CReATe Consortium's Phenotype-Genotype-Biomarker (PGB) study. These results were replicated using an independent ethnically matched WGS data set from Target ALS. Significant results were further explored with available databases of genetic regulatory markers and expression quantitative trait loci (eQTL) analysis. WGS from 22 participants with documented ALS Reversals was compared with the PGB primary cohort (n = 103) and the Target ALS validation cohort (n = 140). Two genetic loci met predefined criteria for statistical significance (two-sided permutation p ≤ 0.01) and remained plausible after fine-mapping. The lead single nucleotide variant (SNV) from the first locus was rs4242007 (primary cohort GWAS OR = 12.0, 95% CI 4.1 to 34.6), which is in an IGFBP7 intron and is in near-perfect linkage disequilibrium with a SNV in the IGFBP7 promoter region. Both SNVs are associated with decreased frontal cortex IGFBP7 expression in eQTL data sets. Notably, 3 Reversals, but none of the typically progressive individuals (n = 243), were homozygous for rs4242007. The importance of the second locus, located near GRIP1, is uncertain given the absence of an associated effect on nearby gene transcription. We found a significant association between the Reversal phenotype and an IGFBP7 noncoding SNV that is associated with IGFBP7 expression. This is biologically relevant as IGFBP7 is a reported inhibitor of the insulin growth factor-1 (IGF-1) receptor that activates the possibly neuroprotective IGF-1 signaling pathway. This finding is limited by small sample size but suggests that there may be merit in further exploration of IGF-1 pathway signaling as a therapeutic mechanism for ALS. This study was registered with ClinicalTrials.gov (NCT03464903) on March 14, 2018. The first participant was enrolled on June 22, 2018.

Unraveling the molecular mechanisms of Ace2-mediated post-COVID-19 cognitive dysfunction through systems genetics approach

The dysregulation of Angiotensin-converting enzyme 2 (ACE2) in central nervous system is believed associates with COVID-19 induced cognitive dysfunction. However, the detailed mechanism remains largely unknown. In this study, we performed a comprehensive system genetics analysis on hippocampal ACE2 based on BXD mice panel. Expression quantitative trait loci (eQTLs) mapping showed that Ace2 was strongly trans-regulated, and the elevation of Ace2 expression level was significantly correlated with impaired cognitive functions. Further Gene co-expression analysis showed that Ace2 may be correlated with the membrane proteins in Calcium signaling pathway. Further, qRT-PCR confirmed that SARS-CoV-2 spike S1 protein upregulated ACE2 expression together with eight membrane proteins in Calcium Signaling pathway. Moreover, such elevation can be attenuated by recombinant ACE2. Collectively, our findings revealed a potential mechanism of Ace2 in cognitive dysfunction, which could be beneficial for COVID–19–induced cognitive dysfunction prevention and potential treatment.

Insights into preeclampsia: a bioinformatics approach to deciphering genetic and immune contributions

BackgroundPreeclampsia (PE) is a global pregnancy concern, characterized by hypertension with an unclear etiology. This study employs Mendelian randomization (MR) and single-cell RNA sequencing (scRNA-seq) to clarify its genetic and molecular roots, offering insights into diagnosis and treatment avenues.MethodsWe integrated PE-specific genome-wide association study (GWAS) data, expression and protein quantitative trait loci (eQTL and pQTL) data, and single-cell data from peripheral blood mononuclear cells (PBMCs). We identified highly variable genes using single-cell information and employed MR to determine potential causality. We also combined pQTL and GWAS data, discerned genes positively associated with PE through scRNA-seq, and leveraged the Enrichr platform to unearth drug-gene interactions.ResultsOur scRNA-seq pinpointed notable cell type distribution variances, especially in T helper cells (Th cells), between PE and control groups. We unveiled 591 highly variable genes and 6 directly PE-associated genes. Although MR revealed correlations with PE risk, pQTL analysis was inconclusive due to data constraints. Using DSigDB, 93 potential therapeutic agents, like Retinoic acid targeting core genes (IFITM3, NINJ1, COTL1, CD69, and YWHAZ), emerged as prospective multi-target treatments.ConclusionUtilizing MR and scRNA-seq, this study underscores significant cellular disparities, particularly in Th cells, and identifies crucial genes related to PE. Despite some limitations, these genes have been revealed in PE’s underlying mechanism. Potential therapeutic agents, such as Retinoic acid, suggest promising treatment pathways.

Rs12411980 single-nucleotide polymorphism related to PRTFDC1 expression is significantly associated with phantom tooth pain.

Phantom tooth pain (PTP) is one type of non-odontogenic neuropathic toothache, which rarely occurs after appropriate pulpectomy or tooth extraction. The cause of PTP is unknown. We investigated pain-related genetic factors that are associated with PTP. Four pain-associated genes, including G protein-coupled receptor 158 (GPR158) and phosphoribosyl transferase domain containing 1 (PRTFDC1), are adjacent to each other on the human genome. Some of these four genes or their genomic region may be related to PTP. We statistically analyzed associations between single-nucleotide polymorphisms (SNPs) in the genomic region and PTP in patients with PTP (PTP group), other orofacial pain (OFP group), and healthy control subjects. We then performed a database search of expression quantitative trait loci (eQTLs). For the seven SNPs that were significantly associated with PTP even after Bonferroni correction, we focused on the rs12411980 tag SNP (P = 9.42 × 10-4). Statistical analyses of the PTP group and healthy subject groups (group labels: NOC and TD) revealed that the rate of the GG genotype of the rs12411980 SNP was significantly higher in the PTP group than in the healthy subject groups (PTP group vs. NOC group: P = 2.92 × 10-4, PTP group vs. TD group: P = 5.46 × 10-4; percentage of GG: 30% in PTP group, 12% in NOC group, 11% in TD group). These results suggest that the GG genotype of the rs12411980 SNP is more susceptible to PTP. The rs2765697 SNP that is in strong linkage disequilibrium with the rs12411980 SNP is an eQTL that is associated with higher PRTFDC1 expression in the minor allele homozygotes in the healthy subject groups of the rs2765697 SNP. Thus, PRTFDC1 expression similarly increases in the minor allele homozygotes (GG genotype) in the healthy subject groups of the rs12411980 SNP, which would lead to greater susceptibility to PTP.

Integrative analysis of transcriptome-wide association study and mRNA expression profile identified risk genes for bipolar disorder

ObjectiveBipolar disorder (BD) is a debilitating neuropsychiatric disorder, which is associated with genetic variation through “vast but mixed” Genome-Wide Association Studies (GWAS). Transcriptome-Wide Association Study (TWAS) is more effective in explaining genetic factors that influence complex diseases and can help identifying risk genes more reliably. So, this study aims to identify potential BD risk genes in pedigrees with TWAS. MethodsWe conducted a TWAS analysis with expression quantitative trait loci (eQTL) analysis on extended BD pedigrees, and the BD genome-wide association study (GWAS) summary data acquired from the Psychiatric Genomics Consortium (PGC). Furthermore, the BD-associated genes identified by TWAS were validated by mRNA expression profiles from the Gene Expression Omnibus (GEO) Datasets (GSE23848 and GSE46416). Functional enrichment and annotation analysis were implemented by RStudio (version 4.2.0). ResultsTWAS identified 362 genes with P value < 0.05, and 18 genes remain significant after Bonferroni correction, such as SEMA3G (PTWAS=1.07 × 10-11), ALOX5AP (PTWAS=3.12 × 10-8), and PLEC (PTWAS=1.27 × 10-7). Further 6 overlapped genes were detected in integrative analysis, such as UQCRB (PTWAS=0.0020, PmRNA=0.0000), TMPRSS9 (PTWAS=0.0405, PmRNA=0.0032), and SNX10 (PTWAS=0.0104, PmRNA=0.0015). Using genes identified by TWAS, Gene Ontology (GO) enrichment analysis identified 40 significant GO terms, such as mitochondrial ATP synthesis coupled electron transport, mitochondrial respiratory, aerobic electron transport chain, oxidative phosphorylation, mitochondrial membrane proteins, and ubiquinone activity. The Kyoto Encyclopedia of Genes and Genomes (KEGG) Pathway enrichment analysis identified significant 15 pathways for BD, such as Oxidative phosphorylation, endocannabinoids signaling, neurodegeneration, and reactive oxide species. ConclusionsWe found a set of BD-associated genes and pathways, validating the important role of neurodevelopmental abnormalities, inflammatory responses, and mitochondrial dysfunction in the pathology of BD, offering novel information for comprehending the genetic basis of BD.

Genetic regulation of B cell receptor signaling pathway: Insights from expression quantitative trait locus analysis using a mixed model

The B cell receptor (BCR) signaling pathway regulates non-immune cellular response through various pathways like MAPK, NF-kB, and PI3K-Akt. This study aimed to identify expression quantitative trait loci (eQTL) and their regulatory functions on BCR signaling pathway genes. A mixed model was employed to analyze eQTL using RNA expression levels in lymphoblastoid from 376 Europeans in the GEUVADIS dataset. In total, 266 SNPs, including 115 cis-acting SNPs, were found for association with transcription of 13 genes (P < 5 × 10-8), revealing 19 independent signals for five genes through linkage disequilibrium analysis. Functional analysis, aligning them with DNase sensitive sites, transcription factor binding sites, histone modification, promoters/enhancers, CpG islands, and ChIA-PET, identified regulatory variants targeting SYK, VAV2, and PLCG2. Notably, rs2562397 was validated as a SYK promoter variant, and rs694505, rs636667, and rs4889409 were confirmed as enhancer variants for VAV2 and PLCG2. Their allelic differences in gene expression were also confirmed using ENCODE ChIP-seq and Sei neural network prediction. Persistent differential expression of these genes by alleles might impact the adaptive immune system, vascular development, and/or relevant diseases that have been previously associated with other variants of the genes. Comprehensive genetic architecture studies of the BCR signaling pathway, along with experiments demonstrating related mechanisms, will greatly contribute to understanding the underlying mechanisms of relevant disease development and implementing precision medicine.

Novel genetic variants in the NLRP3 inflammasome-related PANX1 and APP genes predict survival of patients with hepatitis B virus-related hepatocellular carcinoma.

The nod-like receptor protein 3 (NLRP3) is one of the most characterized inflammasomes involved in the pathogenesis of several cancers, including hepatocellular carcinoma (HCC). However, the effects of genetic variants in the NLRP3 inflammasome-related genes on survival of hepatitis B virus (HBV)-related HCC patients are unclear. We performed multivariable Cox proportional hazards regression analysis to evaluate associations between 299 single-nucleotide polymorphisms (SNPs) in 16 NLRP3 inflammasome-related genes and overall survival (OS) of 866 patients with HBV-related HCC. We further performed expression quantitative trait loci (eQTL) analysis using the data from the GTEx project and 1000 Genomes projects, and performed differential expression analysis using the TCGA dataset to explore possible molecular mechanisms underlying the observed associations. We found that two functional SNPs (PANX1 rs3020013 A > G and APP rs9976425 C > T) were significantly associated with HBV-related HCC OS with the adjusted hazard ratio (HR) of 0.83 [95% confidence interval (CI) = 0.73-0.95, P = 0.008], and 1.26 (95% CI = 1.02-1.55, P = 0.033), respectively. Moreover, the eQTL analysis revealed that the rs3020013 G allele was correlated with decreased mRNA expression levels of PANX1 in both normal liver tissues (P = 0.044) and whole blood (P < 0.001) in the GTEx dataset, and PANX1 mRNA expression levels were significantly higher in HCC samples and associated with a poorer survival of HCC patients. However, we did not observe such correlations for APP rs9976425. These results indicated that SNPs in the NLRP3 inflammasome-related genes may serve as potential biomarkers for HBV-related HCC survival, once replicated by additional larger studies.

New Insights on the Therapeutic Potential of Runt-Related Transcription Factor2 for Osteoarthritis: Evidence from Mendelian Randomization.

Research has highlighted the role of runt-related transcription factor2 (Runx2) in the development of osteoarthritis (OA); however, its causal association remains unclear. This study aimed to explore whether Runx2 expression is causally associated with OA and assess its therapeutic potential for OA. Genetic proxy instruments for Runx2 expression were obtained from gene expression quantitative trait locus (eQTLs) study of eQTLGen Consortium (n = 31,684). Aggregated genome-wide association study (GWAS) data for OA (including all OA [177,517 cases and 649,173 controls], knee OA (KOA) [62,497 cases and 333,557 controls], and hip OA (HOA) [36,445 cases and 316,943 controls]) were extracted from the Genetics of Osteoarthritis Consortium. We integrated eQTLs data with OA GWAS data to estimate their causal association and to estimate the potential of Runx2 as a drug target in the treatment of OA using summary data-based Mendelian randomization (SMR) analysis. Furthermore, different OA GWAS data (including all OA [77,052 cases and 378,169 controls], KOA [24,955 cases and 378,169 controls], and HOA [15,704 cases and 378,169 controls]) derived from the GWAS Catalog database were used for replication study. SMR analysis showed that high expression levels of Runx2 were associated with an increased risk of all OA [odds ratio (OR) 1.044, 95% confidence interval (CI) 1.023-1.067; P = 5.03 × 10-5], KOA (OR1.040, 95%CI 1.006-1.075; P = 0.021), and HOA (OR1.067, 95%CI 1.022-1.113; P = 0.003). This suggests that Runx2 inhibitors may have promising potential for the treatment of OA. Notably, the causal effects of Runx2 with all OA (OR1.053, 95%CI 1.027-1.079; P = 3.95 × 10-5) and KOA (OR1.043, 95%CI 1.001-1.087; P = 0.045) were repeated in the replication study, but limited evidence supported the association of Runx2 expression levels with HOA (OR1.045, 95%CI 0.993-1.101; P = 0.094). Our analyses indicate a positive correlation between Runx2 expression and OA risk across all three phenotypes, suggesting the potential of Runx2 inhibitors in the treatment of OA and providing evidence from a genetic perspective.

SPATS2L is a positive feedback regulator of the type I interferon signaling pathway and plays a vital role in lupus.

Through genome-wide association studies (GWAS) and integrated expression quantitative trait locus (eQTL) analyses, numerous susceptibility genes ("eGenes", whose expressions are significantly associated with common variants) associated with systemic lupus erythematosus (SLE) have been identified. Notably, a subset of these eGenes is correlated with disease activity. However, the precise mechanisms through which these genes contribute to the initiation and progression of the disease remain to be fully elucidated. In this investigation, we initially identify SPATS2L as an SLE eGene correlated with disease activity. eSignaling and transcriptomic analyses suggest its involvement in the type I interferon (IFN) pathway. We observe a significant increase in SPATS2L expression following type I IFN stimulation, and the expression levels are dependent on both the concentration and duration of stimulation. Furthermore, through dual-luciferase reporter assays, western blot analysis, and imaging flow cytometry, we confirm that SPATS2L positively modulates the type I IFN pathway, acting as a positive feedback regulator. Notably, siRNA-mediated intervention targeting SPATS2L, an interferon-inducible gene, in peripheral blood mononuclear cells (PBMCs) from patients with SLE reverses the activation of the interferon pathway. In conclusion, our research highlights the pivotal role of SPATS2L as a positive-feedback regulatory molecule within the type I IFN pathway. Our findings suggest that SPATS2L plays a critical role in the onset and progression of SLE and may serve as a promising target for disease activity assessment and intervention strategies.

Potentially functional variants of ERRFI1 in hypoxia-related genes predict survival of non-small cell lung cancer patients.

Hypoxia is often involved in tumor microenvironment, and the hypoxia-induced signaling pathways play a key role in aggressive cancer phenotypes, including angiogenesis, immune evasion, and therapy resistance. However, it is unknown what role genetic variants in the hypoxia-related genes play in survival of patients with non-small cell lung cancer (NSCLC). We evaluated the associations between 16,092 single-nucleotide polymorphisms (SNPs) in 182 hypoxia-related genes and survival outcomes of NSCLC patients. Data from the Prostate, Lung, Colorectal, and Ovarian (PLCO) Cancer Screening Trial were used as the discovery dataset, and the Harvard Lung Cancer Susceptibility (HLCS) Study served as the replication dataset. We also performed additional linkage disequilibrium analysis and a stepwise multivariable Cox proportional hazards regression analysis in the PLCO dataset. An independent SNP, ERRFI1 rs28624 A > C, was identified with an adjusted hazards ratio (HR) of 1.31 (95% CI = 1.14-1.51, p = 0.0001) for overall survival (OS). In further analyses, unfavorable genotypes AC and CC, compared with the AA genotype, were associated a worse OS (HR = 1.20, 95% CI = 1.03-1.39, p = 0.014) and disease-specific survival (HR = 1.21, 95% CI = 1.04-1.42, p = 0.016). Further expression quantitative trait loci analysis indicated that ERRFI1 rs28624C genotypes were significantly associated with higher ERRFI1 mRNA expression levels in the whole blood. Additional analysis showed that high ERRFI1 mRNA expression levels were associated with a worse OS in patients with lung adenocarcinoma. Our findings suggest that genetic variants in the hypoxia-related gene ERRFI1 may modulate NSCLC survival, potentially through their effect on the gene expression.

Prognostic Significance of VAV3 Gene Variants and Expression in Renal Cell Carcinoma.

Renal cell carcinoma (RCC) is characterized by high mortality and morbidity rates. Vav guanine nucleotide exchange factors (VAVs), crucial for signal transduction between cell membrane receptors and intracellular mediators, have been implicated in carcinogenesis. However, their potential prognostic value in RCC remains unclear. The impact of 150 common VAV polymorphisms on RCC risk and survival was investigated in a cohort of 630 individuals. Publicly available gene expression datasets were utilized to analyze VAV gene expression in relation to patient outcomes. The VAV3 rs17019888 polymorphism was significantly associated with RCC risk and overall survival after adjusting for false discovery rates. Expression quantitative trait loci analysis revealed that the risk allele of rs17019888 is linked to reduced VAV3 expression. Analysis of 19 kidney cancer gene expression datasets revealed lower VAV3 expression in RCC tissues compared to normal tissues, with higher expression correlating with better prognosis. Gene set enrichment analysis demonstrated that VAV3 negatively regulates the ubiquitin-proteasome system, extracellular matrix and membrane receptors, inflammatory responses, matrix metalloproteinases, and cell cycle pathways. Furthermore, elevated VAV3 expression was associated with increased infiltration of B cells, macrophages, and neutrophils into the RCC tumor microenvironment. Our findings suggest that VAV3 gene variants influence RCC risk and survival, contributing to a favorable prognosis in RCC.

Systems genetic analysis of lignin biosynthesis in Populus tremula.

The genetic control underlying natural variation in lignin content and composition in trees is not fully understood. We performed a systems genetic analysis to uncover the genetic regulation of lignin biosynthesis in a natural 'SwAsp' population of aspen (Populus tremula) trees. We analyzed gene expression by RNA sequencing (RNA-seq) in differentiating xylem tissues, and lignin content and composition using Pyrolysis-GC-MS in mature wood of 268 trees from 99 genotypes. Abundant variation was observed for lignin content and composition, and genome-wide association study identified proteins in the pentose phosphate pathway and arabinogalactan protein glycosylation among the top-ranked genes that are associated with these traits. Variation in gene expression and the associated genetic polymorphism was revealed through the identification of 312 705 local and 292 003 distant expression quantitative trait loci (eQTL). A co-expression network analysis suggested modularization of lignin biosynthesis and novel functions for the lignin-biosynthetic CINNAMYL ALCOHOL DEHYDROGENASE 2 and CAFFEOYL-CoA O-METHYLTRANSFERASE 3. PHENYLALANINE AMMONIA LYASE 3 was co-expressed with HOMEOBOX PROTEIN 5 (HB5), and the role of HB5 in stimulating lignification was demonstrated in transgenic trees. The systems genetic approach allowed linking natural variation in lignin biosynthesis to trees´ responses to external cues such as mechanical stimulus and nutrient availability.

Gene expression and splicing QTL analysis of blood cells in African American participants from the Jackson Heart Study.

Most gene expression and alternative splicing quantitative trait loci (eQTL/sQTL) studies have been biased toward European ancestry individuals. Here, we performed eQTL and sQTL analyses using TOPMed whole-genome sequencing-derived genotype data and RNA-sequencing data from stored peripheral blood mononuclear cells in 1,012 African American participants from the Jackson Heart Study (JHS). At a false discovery rate of 5%, we identified 17,630 unique eQTL credible sets covering 16,538 unique genes; and 24,525 unique sQTL credible sets covering 9,605 unique genes, with lead QTL at P < 5e-8. About 24% of independent eQTLs and independent sQTLs with a minor allele frequency > 1% in JHS were rare (minor allele frequency < 0.1%), and therefore unlikely to be detected, in European ancestry individuals. Finally, we created an open database, which is freely available online, allowing fast query and bulk download of our QTL results.

Causal Effects of Oxidative Stress on Diabetes Mellitus and Microvascular Complications: Insights Integrating Genome-Wide Mendelian Randomization, DNA Methylation, and Proteome.

Oxidative stress (OS) is involved in the development of diabetes, but the genetic mechanisms are not completely understood. We integrated multi-omics data in order to explore the genetic relations between OS-related genes, diabetes mellitus, and microvascular complications using Mendelian randomization and colocalization analysis. Summary-level data related to OS were acquired from respective studies of methylation, expression, and protein abundance quantitative trait loci. Genetic associations concerning diabetes, diabetic nephropathy (DN), and diabetic retinopathy (DR) were derived from the FinnGen study. Summary-data-based Mendelian randomization (SMR) analysis was conducted to evaluate the correlations between molecular features concerned with OS-related genes and diabetes mellitus, along with its microvascular complications. Additionally, we performed colocalization analysis to determine if the detected signal pairs shared a causal genetic variant. At the genetic level, we identified ten potential causal associations of oxidative stress genes with diabetes, along with microvascular complications, through SMR and colocalization analysis. After integrating the DNA methylation quantitative trait loci (mQTL) and expression QTL (eQTL) data, our analyses revealed a correlation between the methylation site cg26343298 and reduced expression of TP53INP1, supporting the protective role of cg26343298 methylation on type 2 diabetes (T2D) and diabetic nephropathy. Similarly, an inverse association was observed between gene methylation and expression in CHEK1 (cg07110182), confirming the beneficial effect of modification of CHEK1 by cg07110182 in diabetic retinopathy. In addition, upregulation of SUOX expression by cg22580629 was linked to a reduced risk of diabetic retinopathy. At circulating protein levels, genetically predicted a higher level of ICAM1 (OR 1.05, 95%CI 1.03-1.08) was positively connected with the risk of diabetic retinopathy. This SMR study elucidated that the TP53INP1 gene was putatively associated with T2D and DN risk, while the SUOX and CHEK1 genes were associated with DR risk through oxidative stress mechanisms. Additionally, our study showed a positive correlation between the ICAM-1 protein and DR. These findings may enhance our understanding of their pathogenesis and suggest new therapeutic targets for clinical practice.

Identification of immune-inflammation targets for intracranial aneurysms: A multiomics and epigenome-wide study integrating summary-data-based mendelian randomization, single-cell-type expression analysis, and DNA methylation regulation.

Dysfunction of the immune system and inflammation plays a vital role in developing intracranial aneurysms (IAs). However, the progress of genetic pathophysiology is complicated and not entirely elaborated. This study aimed to explore the genetic associations of immune- and inflammation-related genes (IIRGs) with IAs and their subtypes using Mendelian randomization, colocalization test, and integrated multiomics functional analysis. We conducted a summary-data-based Mendelian randomization (SMR) analysis using data from several genome-wide association studies of gene expression (31,684 European individuals) and protein quantitative trait loci (35,559 Icelanders), as well as information on IAs and their subtypes from The International Stroke Genetics Consortium (IGSC) for discovery phase and the FinnGen study for replication. This analysis aimed to determine the causal relationship between IIRGs and the risk of IAs and their subtypes. Further functional analyses, including DNA methylation regulation (1980, European individuals), single-cell-type expression analysis, and protein-protein interaction, were conducted to detect the specific cell type with enriched expression and discover potential drug targets. After integrating multi-omics evidence from expression quantitative trait loci (eQTL)and protein quantitative trait loci(pQTL), we found that tier 1: RELT [odds ratio (OR): 0.14, 95% confidence interval (CI): 0.04-0.50], TNFSF12 (OR: 1.24, 95% CI: 1.24-1.43), tier 3:ICAM5 (OR: 0.89, 95% CI: 0.82-0.96), and ERAP2 (OR: 1.07, 95% CI: 1.02-1.12) were associated with the risk of IAs; tier 3: RELT (OR: 0.11, 95% CI: 0.02-0.54), ERAP2 (OR: 1.08, 95% CI: 1.02-1.13), and TNFSF12 (OR: 1.24, 95% CI: 1.05-1.47) were associated with the risk of aneurysmal subarachnoid hemorrhage (aSAH); and tier 1:RELT (OR: 0.04, 95% CI: 0.01-0.30) was associated with the risk of unruptured intracranial aneurysms (uIAs). Further functional analyses showed that RELT was regulated by cg06382664 and cg18850434 and ICAM5 was regulated by cg04295144 in IAs; RELT was regulated by cg06382664, cg08770935, cg16533363, and cg18850434 in aSAH; and RELT was regulated by cg06382664 and cg21810604 in uIAs. In addition, we found that H6PD (OR: 1.13, 95% CI: 1.01-1.28), NT5M (OR: 1.91, 95% CI: 1.21-3.01), and NPTXR (OR: 1.13, 95% CI: 1.01-1.26) were associated with IAs; NT5M (OR: 2.13, 95% CI: 1.23-3.66) was associated aSAH; and AP4M1 (OR: 0.06, 95% CI: 0.01-0.42) and STX7 (OR: 3.97, 95% CI: 1.41-11.18) were related to uIAs. STX7 and TNFSF12 were mainly enriched in microglial cells, whereas H6PD, STX7, and TNFSF12 were mainly enriched in astrocytes. After integrating multi-omics evidence, we eventually identified IIRGs: RELT, TNFSF12, ICAM5 and ERAP2 were the novel therapy targets for IAs. These new results confirmed a vital role of immune and inflammation in the etiology of IAs, contributing to enhance our understanding of the immune and inflammatory mechanisms in the pathogenesis of IAs and revealing the complex genetic causality of IAs.