Gene Expression Quantitative Trait Loci Research Articles

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.

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.

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.

Identification of drug targets for Sjögren's syndrome: multi-omics Mendelian randomization and colocalization analyses.

Targeted therapy for Sjögren's syndrome (SS) has become an important focus for clinicians. Multi-omics-wide Mendelian randomization (MR) analyses have provided new ideas for identifying potential drug targets. We conducted summary-data-based Mendelian randomization (SMR) analysis to evaluate therapeutic targets associated with SS by integrating DNA methylation, gene expression and protein quantitative trait loci (mQTL, eQTL, and pQTL, respectively). Genetic associations with SS were derived from the FinnGen study (discovery) and the GWAS catalog (replication). Colocalization analyses were employed to determine whether two potentially relevant phenotypes share the same genetic factors in a given region. Moreover, to delve deeper into potential regulation among DNA methylation, gene expression, and protein abundance, we conducted MR analysis to explore the causal relationship between candidate gene methylation and expression, as well as between gene expression and protein abundance. Drug prediction and molecular docking were further employed to validate the pharmacological activity of the candidate drug targets. Upon integrating the multi-omics data, we identified three genes associated with SS risk: TNFAIP3, BTN3A1, and PLAU. The methylation of cg22068371 in BTN3A1 was positively associated with protein levels, consistent with the negative effect of cg22068371 methylation on the risk of SS. Additionally, positive correlations were observed between the gene methylation of PLAU (cg04939496) and expression, as well as between expression and protein levels. This consistency elucidates the promotional effects of PLAU on SS risk at the DNA methylation, gene expression, and protein levels. At the protein level, genetically predicted TNFAIP3 (OR 2.47, 95% CI 1.56-3.92) was positively associated with SS risk, while BTN3A1 (OR 2.96E-03, 95% CI 2.63E-04-3.33E-02) was negatively associated with SS risk. Molecular docking showed stable binding for candidate drugs and target proteins. Our study reveals promising therapeutic targets for the treatment of SS, providing valuable insights into targeted therapy for SS. However, further validation through future experiments is warranted.

Systematic druggable genome-wide Mendelian randomization identifies therapeutic targets for lung cancer

BackgroundDrug repurposing provides a cost-effective approach to address the need for lung cancer prevention and therapeutics. We aimed to identify actionable druggable targets using Mendelian randomization (MR).MethodsSummary-level data of gene expression quantitative trait loci (eQTLs) were sourced from the eQTLGen resource. We procured genetic associations with lung cancer and its subtypes from the TRICL, ILCCO studies (discovery) and the FinnGen study (replication). We implemented Summary-data-based Mendelian Randomization analysis to identify potential therapeutic targets for lung cancer. Colocalization analysis was further conducted to assess whether the identified signal pairs shared a causal genetic variant.FindingsIn the main analysis dataset, we identified 55 genes that demonstrate a causal relationship with lung cancer and its subtypes. However, in the replication cohort, only three genes were found to have such a causal association with lung cancer and its subtypes, and of these, HYKK (also known as AGPHD1) was consistently present in both the primary analysis dataset and the replication cohort. Following HEIDI tests and colocalization analyses, it was revealed that HYKK (AGPHD1) is associated with an increased risk of squamous cell carcinoma of the lung, with an odds ratio and confidence interval of OR = 1.28,95%CI = 1.24 to 1.33.InterpretationWe have found that the HYKK (AGPHD1) gene is associated with an increased risk of squamous cell carcinoma of the lung, suggesting that this gene may represent a potential therapeutic target for both the prevention and treatment of lung squamous cell carcinoma.

TWAS facilitates gene-scale trait genetic dissection through gene expression, structural variations, and alternative splicing in soybean

A genome-wide association study (GWAS) identifies trait-associated loci, but identifying the causal genes can be a bottleneck, due in part to slow decay of linkage disequilibrium (LD). A transcriptome-wide association study (TWAS) addresses this issue by identifying gene expression–phenotype associations or integrating gene expression quantitative trait loci with GWAS results. Here, we used self-pollinated soybean (Glycine max [L.] Merr.) as a model to evaluate the application of TWAS to the genetic dissection of traits in plant species with slow LD decay. We generated RNA sequencing data for a soybean diversity panel and identified the genetic expression regulation of 29 286 soybean genes. Different TWAS solutions were less affected by LD and were robust to the source of expression, identifing known genes related to traits from different tissues and developmental stages. The novel pod-color gene L2 was identified via TWAS and functionally validated by genome editing. By introducing a new exon proportion feature, we significantly improved the detection of expression variations that resulted from structural variations and alternative splicing. As a result, the genes identified through our TWAS approach exhibited a diverse range of causal variations, including SNPs, insertions or deletions, gene fusion, copy number variations, and alternative splicing. Using this approach, we identified genes associated with flowering time, including both previously known genes and novel genes that had not previously been linked to this trait, providing insights complementary to those from GWAS. In summary, this study supports the application of TWAS for candidate gene identification in species with low rates of LD decay.

#2347 Multi-omics analysis identifies traditional Chinese medicine targets for chronic kidney disease

Abstract Background and Aims As a major public health problem, chronic kidney disease (CKD) has substantial comorbidities and disease burden. It is necessary to discover new therapeutic targets considering its irreversibly. Traditional Chinese medicine (TCM) has become indispensable for the integrated treatment of CKD. The current study aims to explore potential new TCM targets in high-throughput experiment- and reference-guided database (HERB) that are related to renal phenotypes utilizing multi-omics. Method We conducted comprehensive Mendelian Randomization and Bayesian colocalization on more than 2,000 HERB target integrating genomes and transcriptomes, using summary statistics of genome-wide association study and cis gene expression quantitative trait loci. The results were validated through differential gene expression in public database, previous literature reports and clinical kidney biopsy specimens. Results In all, we discovered 34 targets associated with renal phenotypes among the 2,421 HERB TCM targets in the data mining part. Among them, most of the targets can be validated by differential gene expression analysis in GEO and Nephroseq, and have been reported by basic research, thus confirming the reliability of this method. In addition, PRKCI, KNG1, PKN3, GGT7, GBA, APEH and GSTA2 are expected to be potential new therapeutic targets for CKD. Immunohistochemical staining using clinical kidney biopsy specimens of PKN3 and KNG1, two targets with the highest evidence level, further confirm the results. Conclusion Integration of multi-omics can extensively and efficiently screen potential drug targets related to CKD. The new targets discovered in the current study, especially PKN3 and KNG1, need to be further validated in basic research.

Determinants of gastric cancer immune escape identified from non-coding immune-landscape quantitative trait loci

The landscape of non-coding mutations in cancer progression and immune evasion is largely unexplored. Here, we identify transcrptome-wide somatic and germline 3′ untranslated region (3′-UTR) variants from 375 gastric cancer patients from The Cancer Genome Atlas. By performing gene expression quantitative trait loci (eQTL) and immune landscape QTL (ilQTL) analysis, we discover 3′-UTR variants with cis effects on expression and immune landscape phenotypes, such as immune cell infiltration and T cell receptor diversity. Using a massively parallel reporter assay, we distinguish between causal and correlative effects of 3′-UTR eQTLs in immune-related genes. Our approach identifies numerous 3′-UTR eQTLs and ilQTLs, providing a unique resource for the identification of immunotherapeutic targets and biomarkers. A prioritized ilQTL variant signature predicts response to immunotherapy better than standard-of-care PD-L1 expression in independent patient cohorts, showcasing the untapped potential of non-coding mutations in cancer.

Identification and characterization of whole blood gene expression and splicing quantitative trait loci during early to mid-lactation of dairy cattle

BackgroundCharacterization of regulatory variants (e.g., gene expression quantitative trait loci, eQTL; gene splicing QTL, sQTL) is crucial for biologically interpreting molecular mechanisms underlying loci associated with complex traits. However, regulatory variants in dairy cattle, particularly in specific biological contexts (e.g., distinct lactation stages), remain largely unknown. In this study, we explored regulatory variants in whole blood samples collected during early to mid-lactation (22–150 days after calving) of 101 Holstein cows and analyzed them to decipher the regulatory mechanisms underlying complex traits in dairy cattle.ResultsWe identified 14,303 genes and 227,705 intron clusters expressed in the white blood cells of 101 cattle. The average heritability of gene expression and intron excision ratio explained by cis-SNPs is 0.28 ± 0.13 and 0.25 ± 0.13, respectively. We identified 23,485 SNP-gene expression pairs and 18,166 SNP-intron cluster pairs in dairy cattle during early to mid-lactation. Compared with the 2,380,457 cis-eQTLs reported to be present in blood in the Cattle Genotype-Tissue Expression atlas (CattleGTEx), only 6,114 cis-eQTLs (P < 0.05) were detected in the present study. By conducting colocalization analysis between cis-e/sQTL and the results of genome-wide association studies (GWAS) from four traits, we identified a cis-e/sQTL (rs109421300) of the DGAT1 gene that might be a key marker in early to mid-lactation for milk yield, fat yield, protein yield, and somatic cell score (PP4 > 0.6). Finally, transcriptome-wide association studies (TWAS) revealed certain genes (e.g., FAM83H and TBC1D17) whose expression in white blood cells was significantly (P < 0.05) associated with complex traits.ConclusionsThis study investigated the genetic regulation of gene expression and alternative splicing in dairy cows during early to mid-lactation and provided new insights into the regulatory mechanisms underlying complex traits of economic importance.

Integrated analysis of single cell-RNA sequencing and Mendelian randomization identifies lactate dehydrogenase B as a target of melatonin in ischemic stroke.

Despite the success of single-cell RNA sequencing in identifying cellular heterogeneity in ischemic stroke, clarifying the mechanisms underlying these associations of differently expressed genes remains challenging. Several studies that integrate gene expression and gene expression quantitative trait loci (eQTLs) with genome wide-association study (GWAS) data to determine their causal role have been proposed. Here, we combined Mendelian randomization (MR) framework and single cell (sc) RNA sequencing to study how differently expressed genes (DEGs) mediating the effect of gene expression on ischemic stroke. The hub gene was further validated in the invitro model. We identified 2339 DEGs in 10 cell clusters. Among these DEGs, 58 genes were associated with the risk of ischemic stroke. After external validation with eQTL dataset, lactate dehydrogenase B (LDHB) is identified to be positively associated with ischemic stroke. The expression of LDHB has also been validated in sc RNA-seq with dominant expression in microglia and astrocytes, and melatonin is able to reduce the LDHB expression and activity invitro ischemic models. Our study identifies LDHB as a novel biomarker for ischemic stroke via combining the sc RNA-seq and MR analysis.

Epigenetic variation impacts individual differences in the transcriptional response to influenza infection.

Humans display remarkable interindividual variation in their immune response to identical challenges. Yet, our understanding of the genetic and epigenetic factors contributing to such variation remains limited. Here we performed in-depth genetic, epigenetic and transcriptional profiling on primary macrophages derived from individuals of European and African ancestry before and after infection with influenza A virus. We show that baseline epigenetic profiles are strongly predictive of the transcriptional response to influenza A virus across individuals. Quantitative trait locus (QTL) mapping revealed highly coordinated genetic effects on gene regulation, with many cis-acting genetic variants impacting concomitantly gene expression and multiple epigenetic marks. These data reveal that ancestry-associated differences in the epigenetic landscape can be genetically controlled, even more than gene expression. Lastly, among QTL variants that colocalized with immune-disease loci, only 7% were gene expression QTL, while the remaining genetic variants impact epigenetic marks, stressing the importance of considering molecular phenotypes beyond gene expression in disease-focused studies.

Circadian clock-related genome-wide mendelian randomization identifies putatively genes for ulcerative colitis and its comorbidity

BackgroundCircadian rhythm is crucial to the function of the immune system. Disorders of the circadian rhythm can contribute to inflammatory diseases such as Ulcerative colitis (UC). This Mendelian Randomization (MR) analysis applies genetic tools to represent the aggregated statistical results of exposure to circadian rhythm disorders and UC and its comorbidities, allowing for causal inferences.MethodsSummary statistics of protein, DNA methylation and gene expression quantitative trait loci in individuals of European ancestry (pQTL, mQTL, and eQTL, respectively) were used. Genetic variants located within or near 152 circadian clock-related genes and closely related to circadian rhythm disorders were selected as instrumental variables. Causal relationships with UC and its comorbidities were then estimated through employed Summary data-based Mendelian Randomization (SMR) and Inverse-Variance-Weighted MR (IVW-MR).ResultsThrough preliminary SMR analysis, we identified a potential causal relationship between circadian clock-related genes and UC along with its comorbidities, which was further confirmed by IVW-MR analysis. Our study identified strong evidence of positive correlation involving seven overlapping genes (CSNK1E, OPRL1, PIWIL2, RORC, MAX, PPP5C, and AANAT) through MWAS and TWAS in UC, four overlapping genes (OPRL1, CHRNB2, FBXL17, and SIRT1) in UC with PSC, and three overlapping genes (ARNTL, USP7, and KRAS) in UC with arthropathy.ConclusionsThis SMR study demonstrates the causal effect of circadian rhythm disorders in UC and its comorbidities. Furthermore, our investigation pinpointed candidate genes that could potentially serve as drug targets.

Liver regulatory mechanisms of noncoding variants at lipid and metabolic trait loci

Genome-wide association studies (GWASs) have identified hundreds of risk loci for liver disease and lipid-related metabolic traits, although identifying their target genes and molecular mechanisms remains challenging. We predicted target genes at GWAS signals by integrating them with molecular quantitative trait loci for liver gene expression (eQTL) and liver chromatin accessibility QTL (caQTL). We predicted specific regulatory caQTL variants at four GWAS signals located near EFHD1, LITAF, ZNF329, and GPR180. Using transcriptional reporter assays, we determined that caQTL variants rs13395911, rs11644920, rs34003091, and rs9556404 exhibit allelic differences in regulatory activity. We also performed a protein binding assay for rs13395911 and found that FOXA2 differentially interacts with the alleles of rs13395911. For variants rs13395911 and rs11644920 in putative enhancer regulatory elements, we used CRISPRi to demonstrate that repression of the enhancers altered the expression of the predicted target and/or nearby genes. Repression of the element at rs13395911 reduced the expression of EFHD1, and repression of the element at rs11644920 reduced the expression of LITAF, SNN, and TXNDC11. Finally, we showed that EFHD1 is a metabolically active gene in HepG2 cells. Together, these results provide key steps to connect genetic variants with cellular mechanisms and help elucidate the causes of liver disease.

Osteoporosis and osteoarthritis: a bi-directional Mendelian randomization study

ObjectiveTo investigate the causal relationship between low bone mineral density (BMD) and osteoarthritis (OA) using Mendelian randomization (MR) design.MethodsTwo-sample bi-directional MR analyses were performed using summary-level information on OA traits from UK Biobank and arcOGEN. Sensitivity analyses including MR-Egger, simple median, weighted median, MR pleiotropy residual sum, and outlier approaches were utilized in conjunction with inverse variance weighting (IVW). Gene ontology (GO) enrichment analyses and expression quantitative trait locus (eQTL) colocalization analyses were used to investigate the potential mechanism and shared genes between osteoporosis (OP) and OA.ResultsThe IVW method revealed that genetically predicted low femoral neck BMD was significantly linked with hip (β = 0.105, 95% CI: 0.023–0.188) and knee OA (β = 0.117, 95% CI: 0.049–0.184), but not with other site-specific OA. Genetically predicted low lumber spine BMD was significantly associated with OA at any sites (β = 0.048, 95% CI: 0.011–0.085), knee OA (β = 0.101, 95% CI: 0.045–0.156), and hip OA (β = 0.150, 95% CI: 0.077–0.224). Only hip OA was significantly linked with genetically predicted reduced total bone BMD (β = 0.092, 95% CI: 0.010–0.174). In the reverse MR analyses, no evidence for a causal effect of OA on BMD was found. GO enrichment analysis and eQTL analysis illustrated that DDN and SMAD-3 were the most prominent co-located genes.ConclusionsThese findings suggested that OP may be causally linked to an increased risk of OA, indicating that measures to raise BMD may be effective in preventing OA. More research is required to determine the underlying processes via which OP causes OA.

Identifying actionable druggable targets for breast cancer: Mendelian randomization and population-based analyses

Drug repurposing provides a cost-effective approach to address the need for breast cancer prevention and therapeutics. We aimed to identify actionable druggable targets using Mendelian randomization (MR) and then validate the candidate drugs using population-based analyses. We identified genetic instruments for 1406 actionable targets of approved non-oncological drugs based on gene expression, DNA methylation, and protein expression quantitative trait loci (eQTL, mQTL, and pQTL, respectively). Genome-wide association study (GWAS) summary statistics were obtained from the Breast Cancer Association Consortium (122,977 cases, 105,974 controls). We further conducted a nested case-control study using data retrieved from Swedish registers to validate the candidate drugs that were identified from MR analyses. We identified six significant MR associations with gene expression levels (TUBB, MDM2, CSK, ULK3, MC1R and KCNN4) and two significant associations with gene methylation levels across 21 CpG islands (RPS23 and MAPT). Results from the nested case-control study showed that the use of raloxifene (targeting MAPT) was associated with 35% reduced breast cancer risk (odds ratio, OR, 0.65; 95% confidence interval, CI, 0.51-0.83). However, usage of estradiol, tolterodine, and nitrofurantoin (also targeting MAPT) was associated with increased breast cancer risk, with adjusted ORs and 95% CI of 1.10 (1.07-1.13), 1.16 (1.09-1.24), and 1.09 (1.05-1.13), respectively. The effect of raloxifene and nitrofurantoin lost significance in further validation analyses using active-comparator and new-user design. This large-scale MR analysis, combined with population-based validation, identified eight druggable target genes for breast cancer and suggested that raloxifene is an effective chemoprevention against breast cancer. Swedish Research Council, Cancerfonden, Crafoordska Stiftelsen, Allmänna Sjukhusets i Malmö Stiftelsen för bekämpande av cancer, 111 Project and MAS cancer.

Causality between depression and ankylosing spondylitis in a European population: Results from a Mendelian randomization analysis.

The aim of this study was to explore the application of Mendelian randomization (MR) Egger and inverse variance weighted (IVW) in a causal effect on depression and ankylosing spondylitis (AS). Instrumental variables (IVs) were determined using genome-wide association studies. The 2-sample MR analysis was conducted by MR Egger to test the causal effect between depression and AS. The pleiotropy of potential instrumental variables was evaluated. The results of MR Egger and IVW were further compared. A total of 3 single nucleotide polymorphisms as the construct IVs were included. IVW results showed a significant causal effect between depression and AS (P < .001). Depression could promote the risk of AS (odds ratio = 1.060, 95% confidence interval: 1.026-1.094). However, the MR Egger showed no causal effect (P = .311). Heterogeneity statistics suggested that no heterogeneity was existed (P > .05). It was also suggested that there was no horizontal pleiotropy in IVs (MR Egger intercept: -0.0004, P = .471). Reverse MR analysis suggested that there was no causal effect between AS and depression (P > .05). Gene expression quantitative trait locus (QTLs) suggested that rs2517601 and RNF39 were positively correlated (beta = 1.066, P < .001). Depression may be one of the causes of AS by MR analysis in a European population. We can estimate the causal effect based on IVW when horizontal pleiotropy is very tiny.

EQTL Catalogue 2023: New datasets, X chromosome QTLs, and improved detection and visualisation of transcript-level QTLs.

The eQTL Catalogue is an open database of uniformly processed human molecular quantitative trait loci (QTLs). We are continuously updating the resource to further increase its utility for interpreting genetic associations with complex traits. Over the past two years, we have increased the number of uniformly processed studies from 21 to 31 and added X chromosome QTLs for 19 compatible studies. We have also implemented Leafcutter to directly identify splice-junction usage QTLs in all RNA sequencing datasets. Finally, to improve the interpretability of transcript-level QTLs, we have developed static QTL coverage plots that visualise the association between the genotype and average RNA sequencing read coverage in the region for all 1.7 million fine mapped associations. To illustrate the utility of these updates to the eQTL Catalogue, we performed colocalisation analysis between vitamin D levels in the UK Biobank and all molecular QTLs in the eQTL Catalogue. Although most GWAS loci colocalised both with eQTLs and transcript-level QTLs, we found that visual inspection could sometimes be used to distinguish primary splicing QTLs from those that appear to be secondary consequences of large-effect gene expression QTLs. While these visually confirmed primary splicing QTLs explain just 6/53 of the colocalising signals, they are significantly less pleiotropic than eQTLs and identify a prioritised causal gene in 4/6 cases.

Overexpression of REDUCED WALL ACETYLATION C increases xylan acetylation and biomass recalcitrance in Populus.

Plant lignocellulosic biomass, i.e. secondary cell walls of plants, is a vital alternative source for bioenergy. However, the acetylation of xylan in secondary cell walls impedes the conversion of biomass to biofuels. Previous studies have shown that REDUCED WALL ACETYLATION (RWA) proteins are directly involved in the acetylation of xylan but the regulatory mechanism of RWAs is not fully understood. In this study, we demonstrate that overexpression of a Populus trichocarpa PtRWA-C gene increases the level of xylan acetylation and increases the lignin content and S/G ratio, ultimately yielding poplar woody biomass with reduced saccharification efficiency. Furthermore, through gene coexpression network and expression quantitative trait loci (eQTL) analysis, we found that PtRWA-C was regulated not only by the secondary cell wall hierarchical regulatory network but also by an AP2 family transcription factor HARDY (HRD). Specifically, HRD activates PtRWA-C expression by directly binding to the PtRWA-C promoter, which is also the cis-eQTL for PtRWA-C. Taken together, our findings provide insights into the functional roles of PtRWA-C in xylan acetylation and consequently saccharification and shed light on synthetic biology approaches to manipulate this gene and alter cell wall properties. These findings have substantial implications for genetic engineering of woody species, which could be used as a sustainable source of biofuels, valuable biochemicals, and biomaterials.

A computational method for cell type-specific expression quantitative trait loci mapping using bulk RNA-seq data

Mapping cell type-specific gene expression quantitative trait loci (ct-eQTLs) is a powerful way to investigate the genetic basis of complex traits. A popular method for ct-eQTL mapping is to assess the interaction between the genotype of a genetic locus and the abundance of a specific cell type using a linear model. However, this approach requires transforming RNA-seq count data, which distorts the relation between gene expression and cell type proportions and results in reduced power and/or inflated type I error. To address this issue, we have developed a statistical method called CSeQTL that allows for ct-eQTL mapping using bulk RNA-seq count data while taking advantage of allele-specific expression. We validated the results of CSeQTL through simulations and real data analysis, comparing CSeQTL results to those obtained from purified bulk RNA-seq data or single cell RNA-seq data. Using our ct-eQTL findings, we were able to identify cell types relevant to 21 categories of human traits.

Mycobacterium tuberculosis Gene Expression Associated With Fluoroquinolone Resistance and Efflux Pump Inhibition.

We evaluated the relationship between response to efflux pump inhibition in fluoroquinolone-resistant Mycobacterium tuberculosis (Mtb) isolates and differences in gene expression and expression quantitative trait loci (eQTL). We determined ofloxacin minimum inhibitory concentration (MIC) for ofloxacin-resistant and -susceptible Mtb isolates without and with the efflux pump inhibitor verapamil. We performed RNA sequencing (RNA-seq), whole genome sequencing (WGS), and eQTL analysis, focusing on efflux pump, transport, and secretion-associated genes. Of 42 ofloxacin-resistant Mtb isolates, 27 had adequate WGS coverage and acceptable RNA-seq quality. Of these 27, 7 had >2-fold reduction in ofloxacin MIC with verapamil; 6 had 2-fold reduction, and 14 had <2-fold reduction. Five genes (including Rv0191) had significantly increased expression in the MIC fold change >2 compared to <2 groups. Among regulated genes, 31 eQTLs (without ofloxacin) and 35 eQTLs (with ofloxacin) had significant allele frequency differences between MIC fold change >2 and <2 groups. Of these, Rv1410c, Rv2459, and Rv3756c (without ofloxacin) and Rv0191 and Rv3756c (with ofloxacin) have previously been associated with antituberculosis drug resistance. In this first reported eQTL analysis in Mtb, Rv0191 had increased gene expression and significance in eQTL analysis, making it a candidate for functional evaluation of efflux-mediated fluoroquinolone resistance in Mtb.