Genomic atlas of the human plasma proteome

The effect of histo-blood group ABO system transferase (BGAT) on pregnancy related outcomes：A Mendelian randomization study

Protein level of Histo-Blood Group ABO System Transferase (BGAT) has been reported to be associated with cardiometabolic diseases. But its effect on pregnancy related outcomes still remains unclear. Here we conducted a two-sample Mendelian randomization (MR) study to ascertain the putative causal roles of protein levels of BGAT in pregnancy related outcomes. Cis-acting protein quantitative trait loci (pQTLs) robustly associated with protein level of BGAT (P < 5 ×10-8) were used as instruments to proxy the BGAT protein level (N = 35,559, data from deCODE), with two additional pQTL datasets from Fenland (N = 10,708) and INTERVAL (N = 3301) used as validation exposures. Ten pregnancy related diseases and complications were selected as outcomes. We observed that a higher protein level of BGAT showed a putative causal effect on venous complications and haemorrhoids in pregnancy (VH) (odds ratio [OR]=1.19, 95% confidence interval [95% CI]=1.12-1.27, colocalization probability=91%), which was validated by using pQTLs from Fenland and INTERVAL. The Mendelian randomization results further showed effects of the BGAT protein on gestational hypertension (GH) (OR=0.97, 95% CI=0.96-0.99), despite little colocalization evidence to support it. Sensitivity analyses, including proteome-wide Mendelian randomization of the cis-acting BGAT pQTLs, showed little evidence of horizontal pleiotropy. Correctively, our study prioritised BGAT as a putative causal protein for venous complications and haemorrhoids in pregnancy. Future epidemiology and clinical studies are needed to investigate whether BGAT can be considered as a drug target to prevent adverse pregnancy outcomes.

Associations between genetically predicted concentrations of circulating inflammatory cytokines and the risk of ten pregnancy-related adverse outcomes: A two-sample Mendelian randomization study

BACKGROUND

Evidence from increasing observational studies indicates that systemic inflammation plays a role in pregnancy-related adverse events. However, the causal associations between them are largely unclear. To investigate the potential causal effects of genetically regulated concentrations of inflammatory cytokines on the risk of adverse pregnancy outcomes, we performed a Mendelian randomization (MR) analysis.

METHODS

The cis-protein quantitative trait loci for the 47 inflammatory cytokines derived from the latest genome-wide association studies (GWASs) consisting of 31,112 European individuals were used as the instrumental variables. The latest GWAS summary data for the ten adverse pregnancy events were obtained from the FinnGen project (samples ranging from 141,014 to 190,879). The inverse-variance weighted regression or Ward ratio was used as the primary MR analysis method. Sensitivity analyses based on the other five methods were performed to verify MR results. A replication MR analysis was conducted to further clarify the significant associations using data from the UK Biobank.

RESULTS

Twenty-three of the 220 associations were nominally significant (P < 0.05). Among them, seven robust associations survived the Bonferroni correction and passed sensitivity analyses, including positive associations of soluble intercellular adhesion molecule (sICAM-1) with the risk of excessive vomiting in pregnancy, preeclampsia (PE), and pregnancy hypertension (PH), vascular endothelial growth factor with the risk of medical abortion, macrophage colony-stimulating factor (MCSF) with the risk of spontaneous abortion (SA), and an inverse association of macrophage inflammatory protein-1α with the risk of medical abortion. The associations of MCSF with SA, and sICAM-1 with both PE and PH were further confirmed in the replication analysis.

CONCLUSIONS

This study provides further evidence of the role of systemic inflammation, especially endothelial dysfunction in the pathology of adverse pregnancy events, and the identified cytokines warrant in-depth research to explore their underlying mechanisms of action and to evaluate their potential as targets for disease screening, prevention, and treatment in the future.

Matrix Metalloproteinase and Aortic Aneurysm: A Two-sample Mendelian Randomization Study

BACKGROUND

Studies have linked matrix metalloproteinases (MMPs) to both thoracic aortic aneurysm and abdominal aortic aneurysm (TAA and AAA). The precise MMPs entailed in this procedure, however, were still unknown. This study used a two-sample Mendelian randomization (MR) analysis to look into the causal relationship between MMPs and the risk of TAA and AAA.

METHODS

Eight MMPs, including MMP-1, MMP-2, MMP-3, MMP-8, MMP-9, MMP-10, MMP-12, and MMP-13, were found among people of European ancestry with accessible Genome-Wide Association Studies (GWAS). We employed the findings from Genome-Wide Association Studies (GWAS) for 8 MMPs, and TAA and AAA from the FinnGen consortiums (3,201 cases and 317,899 controls, respectively) were used in a two-sample MR analysis. The primary method of analysis for MR was the inverse variance weighted (IVW) method, along with analyses of heterogeneity and horizontal pleiotropy. 31 single-nucleotide polymorphisms connected to MMP were retrieved.

RESULTS

IVW demonstrated a negative causal association between TAA and AAA and serum MMP-12 levels. The incidence of TAA decreased by 1.031% for every 1 ng/mL increase in serum MMP-12 [odds ratio (OR) = 0.897, 95% confidence interval (CI): 0.831-0.968, P = 0.005]. The incidence of AAA fell by 1.653% (OR = 0.835, 95% CI: 0.752-0.926, P = 0.001) for every 1 ng/mL increase in serum MMP-12. There was no horizontal pleiotropy or heterogeneity in the MR data (P > 0.05).

CONCLUSIONS

The levels of TAA and AAA and serum MMP-12 are causally related. MMP-12 is a factor that reduces the risk of AAA and TTA. Our study suggested that MMP-12 level is causally associated with a decreased risk of TAA and AAA.

Mitochondrial biological function and risk of atrial fibrillation and atrial flutter: A 2-sample Mendelian randomization study

Current research suggests that mitochondrial dysfunction can be a contributing factor in the development of cardiac arrhythmias. In pursuit of elucidating the causal link between the biological functions of mitochondria and the occurrence of atrial fibrillation/flutter, we conducted a 2-sample Mendelian randomization (MR) study. Mitochondrial proteins were selected for exposure in this study. To enhance the accuracy of our study, we selected data on AF/AFL from the FinnGen study and the UK Biobank for MR analysis, respectively. The inverse variance-weighted method was utilized as the primary analysis technique for MR. In addition, we performed a series of sensitivity analyses to detect heterogeneity and horizontal pleiotropy. MR results indicated a significant positive association between NAD-dependent protein deacylase sirtuin-5 and AF/AFL (odds ratio = 1.084, 95% confidence interval: 1.037-1.133, P = 3.679 × 10-4, Adjusted P = .024), with consistent outcomes observed in replication analysis (odds ratio = 1.002, 95% confidence interval: 1.001-1.003, P = 4.808 × 10-4, Adjusted P = .032). NAD-dependent protein deacylase sirtuin-5 can significantly promote the occurrence of AF/AFL, and its specific mechanisms warrant further investigation.

The EAT-Lancet Diet Index, Plasma Proteins, and Risk of Heart Failure in a Population-Based Cohort

BACKGROUND

The landmark EAT-Lancet Commission proposed that a planetary health diet is comprised mainly of plant-based foods. However, studies examining whether this diet is associated with heart failure (HF) are currently lacking. In addition, the potential proteomics mechanism on the association between diet and HF warrants further elucidation.

OBJECTIVES

This study aims to both examine the association between the EAT-Lancet diet index and risk of HF and identify plasma proteins underlying such an association.

METHODS

This prospective cohort study included 23,260 participants. HF cases during the follow-up were identified through the Swedish national register. An EAT-Lancet diet index (score range: 0-42) was created to assess adherence to the EAT-Lancet reference diet. In a subcohort (n = 4,742), fasting plasma proteins were quantified.

RESULTS

During a median follow-up of 25.0 years, 1,768 incident HF cases were documented. After adjusting for sociodemographic, lifestyle, diabetes, hypertension, use of lipid-lowering drugs, and body mass index, the HR per 3-point increase of the EAT-Lancet diet index was 0.93 (95% CI: 0.88-0.97). This association was robust in several sensitivity analyses. Among the included 136 plasma proteins, a total of 8 proteins (AM, GDF15, IL-6, TIM, CTSD, CCL20, FS, and FUR) were both inversely associated with the EAT-Lancet diet index and positively associated with risk of HF; the overall proteomic score mediated 9.4% (95% CI: 2.2%-32.1%) of the association.

CONCLUSIONS

Higher adherence to the EAT-Lancet diet was associated with a lower risk of HF. The identified eight plasma proteins provide information on potential pathways mediating such an association.

Mendelian randomization of circulating proteome identifies IFN-γ as a druggable target in aplastic anemia

BACKGROUND

Aplastic anemia (AA) is a kind of bone marrow failure (BMF) characterized by pancytopenia with hypoplasia/aplasia of bone marrow. Immunosuppressive therapy and bone marrow transplantation are effective methods to treat severe aplastic anemia. However, the efficacy is limited by complications and the availability of suitable donors. This study aimed to determine whether any circulating druggable protein levels may have causal effects on AA and provide potential novel drug targets for AA.

METHODS

Genetic variants strongly associated with circulating druggable protein levels to perform Mendelian randomization (MR) analyses were used. The effect of these druggable protein levels on AA risk was measured using the summary statistics from a large-scale proteomic genome-wide association study (GWAS) and FinnGen database ( https://www.finngen.fi/en/access_results ). Multivariable MR analyses were performed to statistically adjust for potential confounders, including platelet counts, reticulocyte counts, neutrophil counts, and proportions of hematopoietic stem cells.

RESULTS

The data showed that higher level of circulating IFN-γ levels was causally associated with AA susceptibility. The causal effects of circulating IFN-γ levels on the AA were broadly consistent, when adjusted for platelet counts, reticulocyte counts, neutrophil counts and proportions of hematopoietic stem cells.

CONCLUSIONS

High levels of circulating IFN-γ levels might increase the risk of AA and might provide a potential novel target for AA.

Sequential immunotherapy: towards cures for autoimmunity

Despite major progress in the treatment of autoimmune diseases in the past two decades, most therapies do not cure disease and can be associated with increased risk of infection through broad suppression of the immune system. However, advances in understanding the causes of autoimmune disease and clinical data from novel therapeutic modalities such as chimeric antigen receptor T cell therapies provide evidence that it may be possible to re-establish immune homeostasis and, potentially, prolong remission or even cure autoimmune diseases. Here, we propose a 'sequential immunotherapy' framework for immune system modulation to help achieve this ambitious goal. This framework encompasses three steps: controlling inflammation; resetting the immune system through elimination of pathogenic immune memory cells; and promoting and maintaining immune homeostasis via immune regulatory agents and tissue repair. We discuss existing drugs and those in development for each of the three steps. We also highlight the importance of causal human biology in identifying and prioritizing novel immunotherapeutic strategies as well as informing their application in specific patient subsets, enabling precision medicine approaches that have the potential to transform clinical care.

Impaired GK-GKRP interaction rather than direct GK activation worsens lipid profiles and contributes to long-term complications: a Mendelian randomization study

BACKGROUND

Glucokinase (GK) plays a key role in glucose metabolism. In the liver, GK is regulated by GK regulatory protein (GKRP) with nuclear sequestration at low plasma glucose level. Some GK activators (GKAs) disrupt GK-GKRP interaction which increases hepatic cytoplasmic GK level. Excess hepatic GK activity may exceed the capacity of glycogen synthesis with excess triglyceride formation. It remains uncertain whether hypertriglyceridemia associated with some GKAs in previous clinical trials was due to direct GK activation or impaired GK-GKRP interaction.

METHODS

Using publicly available genome-wide association study summary statistics, we selected independent genetic variants of GCKR and GCK associated with fasting plasma glucose (FPG) as instrumental variables, to mimic the effects of impaired GK-GKRP interaction and direct GK activation, respectively. We applied two-sample Mendelian Randomization (MR) framework to assess their causal associations with lipid-related traits, risks of metabolic dysfunction-associated steatotic liver disease (MASLD) and cardiovascular diseases. We verified these findings in one-sample MR analysis using individual-level statistics from the Hong Kong Diabetes Register (HKDR).

RESULTS

Genetically-proxied impaired GK-GKRP interaction increased plasma triglycerides, low-density lipoprotein cholesterol and apolipoprotein B levels with increased odds ratio (OR) of 14.6 (95% CI 4.57-46.4) per 1 mmol/L lower FPG for MASLD and OR of 2.92 (95% CI 1.78-4.81) for coronary artery disease (CAD). Genetically-proxied GK activation was associated with decreased risk of CAD (OR 0.69, 95% CI 0.54-0.88) and not with dyslipidemia. One-sample MR validation in HKDR showed consistent results.

CONCLUSIONS

Impaired GK-GKRP interaction, rather than direct GK activation, may worsen lipid profiles and increase risks of MASLD and CAD. Development of future GKAs should avoid interfering with GK-GKRP interaction.

Proteome-wide association study using cis and trans variants and applied to blood cell and lipid-related traits in the Women's Health Initiative study

In most Proteome-Wide Association Studies (PWAS), variants near the protein-coding gene (±1 Mb), also known as cis single nucleotide polymorphisms (SNPs), are used to predict protein levels, which are then tested for association with phenotypes. However, proteins can be regulated through variants outside of the cis region. An intermediate GWAS step to identify protein quantitative trait loci (pQTL) allows for the inclusion of trans SNPs outside the cis region in protein-level prediction models. Here, we assess the prediction of 540 proteins in 1002 individuals from the Women's Health Initiative (WHI), split equally into a GWAS set, an elastic net training set, and a testing set. We compared the testing r2 between measured and predicted protein levels using this proposed approach, to the testing r2 using only cis SNPs. The two methods usually resulted in similar testing r2, but some proteins showed a significant increase in testing r2 with our method. For example, for cartilage acidic protein 1, the testing r2 increased from 0.101 to 0.351. We also demonstrate reproducible findings for predicted protein association with lipid and blood cell traits in WHI participants without proteomics data and in UK Biobank utilizing our PWAS weights.

Visual environment in schools and student depressive symptoms: Insights from a prospective study across multiple cities in eastern China

OBJECTIVE

To investigate the effects of the school visual environment on depressive symptoms in children and adolescents based on cohort study in eastern China. The school visual environment-related indicators included in this study comprise personal factors (visual impairment) and school-related factors (classroom lighting, school green spaces and school air quality).

METHOD

The follow-up cohort comprises 15,348 students from 283 primary and secondary schools in eastern China. This represents the one-year outcomes of a school-based myopia-mental health cohort study. Data collection includes basic demographics (age, gender, region, etc.), physical examination indicators, behavioral indicators, and school visual environment-related indicators.

RESULT

After a one-year follow-up, we found that compared to the more severe vision impairment group (≤4.0), healthy vision group (≥5.0) had a positive effect against the occurrence of depressive symptoms during consecutive follow-ups, with an RR value of 0.61 (95% CI: 0.57-0.66). Higher values of blackboard illumination appear to be associated with greater positive effects, with an RR (Q75%∼Q100% range) value of 0.87(95% CI: 0.81-0.93). School green spaces seem to exhibit relatively good positive effects when in the Q25%∼Q75% range. The combination of physical activity (Weekly high-intensity exercise) with school air quality(PM2.5≤50%)showed a better positive effect, with an RR value of 0.51(95%CI:0.48-0.55).

CONCLUSION

When addressing students' depressive symptoms, it is crucial to improve the visual environment both at the school level and in students' personal level. Paying appropriate attention to modifiable behaviors, like regular participation in high-intensity exercise sessions, can help alleviate students' depressive symptoms.

Plasma proteome mediate the impact of PM2.5 on stroke: A 2-step Mendelian randomization study

The objectives of this study were to measure the mediation effect of plasma proteins and to clarify their mediating role in the relationship between stroke risk and particulate matter 2.5 (PM2.5) exposure. The possible mediating role of plasma proteins on the causative link between PM2.5 exposure and stroke incidence were examined using a two-step Mendelian randomization (MR) approach based on two-sample Mendelian randomization (TSMR). The findings revealed a significant positive causal relationship between PM2.5 exposure and stroke, with an inverse variance weighted odds ratio of 1.219 (95 % CI: 1.002 - 1.482, P < 0.05). Additionally, a positive causal association was identified between PM2.5 exposure and several plasma proteins, including FAM134B, SAP, ITGB7, Elafin, and DCLK3. Among these, FAM134B, ITGB7, Elafin, and DCLK3 also demonstrated a positive causal association with stroke, whereas only SAP was found to be negatively causally associated with stroke. Remarkably, four plasma proteins, namely DCLK3, FAM134B, Elafin, and ITGB7, were identified as mediators, accounting for substantial proportions (14.5 %, 13.6 %, 11.1 %, and 9.9 %) of the causal association between PM2.5 and stroke. These results remained robust across various sensitivity analyses. Consequently, the study highlights the significant and independent impact of PM2.5 on stroke risk and identifies specific plasma proteins as potential targets for preventive interventions against PM2.5-induced stroke.

Assessment of circulating proteins in thyroid cancer: Proteome-wide Mendelian randomization and colocalization analysis

The causality between circulating proteins and thyroid cancer (TC) remains unclear. We employed five large-scale circulating proteomic genome-wide association studies (GWASs) with up to 100,000 participants and a TC meta-GWAS (nCase = 3,418, nControl = 292,703) to conduct proteome-wide Mendelian randomization (MR) and Bayesian colocalization analysis. Protein and gene expressions were validated in thyroid tissue. Through MR analysis, we identified 26 circulating proteins with a putative causal relationship with TCs, among which NANS protein passed multiple corrections (P BH = 3.28e-5, 0.05/1,525). These proteins were involved in amino acids and organic acid synthesis pathways. Colocalization analysis further identified six proteins associated with TCs (VCAM1, LGMN, NPTX1, PLEKHA7, TNFAIP3, and BMP1). Tissue validation confirmed BMP1, LGMN, and PLEKHA7's differential expression between normal and TC tissues. We found limited evidence for linking circulating proteins and the risk of TCs. Our study highlighted the contribution of proteins, particularly those involved in amino acid metabolism, to TCs.

Towards Healthy Longevity: Comprehensive Insights from Molecular Targets and Biomarkers to Biological Clocks

Aging is a complex and time-dependent decline in physiological function that affects most organisms, leading to increased risk of age-related diseases. Investigating the molecular underpinnings of aging is crucial to identify geroprotectors, precisely quantify biological age, and propose healthy longevity approaches. This review explores pathways that are currently being investigated as intervention targets and aging biomarkers spanning molecular, cellular, and systemic dimensions. Interventions that target these hallmarks may ameliorate the aging process, with some progressing to clinical trials. Biomarkers of these hallmarks are used to estimate biological aging and risk of aging-associated disease. Utilizing aging biomarkers, biological aging clocks can be constructed that predict a state of abnormal aging, age-related diseases, and increased mortality. Biological age estimation can therefore provide the basis for a fine-grained risk stratification by predicting all-cause mortality well ahead of the onset of specific diseases, thus offering a window for intervention. Yet, despite technological advancements, challenges persist due to individual variability and the dynamic nature of these biomarkers. Addressing this requires longitudinal studies for robust biomarker identification. Overall, utilizing the hallmarks of aging to discover new drug targets and develop new biomarkers opens new frontiers in medicine. Prospects involve multi-omics integration, machine learning, and personalized approaches for targeted interventions, promising a healthier aging population.

The role of serum interleukins in Cancer: A Multi-center Mendelian Randomization study

The occurrence of cancer is often accompanied by immune evasion and tumor-promoting inflammation, with interleukins (IL) playing a pivotal role in the immune-inflammatory mechanism. However, the precise contribution of serum interleukins in cancer remains elusive. We obtained GWAS summary data for 35 interleukins from eight independent large-scale serum proteome studies of European ancestry populations and for 23 common cancers from the FinnGen Consortium. We then conducted a multicenter Mendelian Randomization (MR) study to explore the relationship between systemic inflammatory status and cancers. 24 causal associations between interleukins and cancers were supported by multicenter data, 18 of which were reported for the first time. Our results indicated that IL-1α (Hodgkin lymphoma), IL-5 (bladder cancer), IL-7 (prostate cancer), IL-11 (bone malignant tumor), IL-16 (lung cancer), IL-17A (pancreatic cancer), IL-20 (bladder cancer), IL-22 (lymphocytic leukemia), IL-34 (breast cancer), IL-36β (prostate cancer), and IL-36γ (liver cancer) were risk factors for related cancers. Conversely, IL-9 (malignant neoplasms of the corpus uteri), IL-17C (liver cancer), and IL-31 (colorectal cancer, bladder cancer, pancreatic cancer, and cutaneous melanoma) exhibited protective effects against related cancers. Notably, the dual effects of serum interleukins were also observed. IL-18 acted as a risk factor for prostate cancer, however, was a protective factor against laryngeal cancer. Similarly, IL-19 promoted the development of lung cancer and myeloid leukemia, while conferring protection against Breast, cervical, and thyroid cancers. Our study confirmed the genetic association between multiple serum interleukins and cancers. Immune and anti-inflammatory strategies targeting these associations provide opportunities for prevention and treatment.

Plasma proteome analysis implicates novel proteins as potential therapeutic targets for chronic kidney disease: A proteome-wide association study

Chronic kidney disease (CKD) is prevalent globally with limited therapeutic drugs available. To systemically identify novel proteins involved in the pathogenesis of CKD and possible therapeutic targets, we integrated human plasma proteomes with the genome-wide association studies (GWASs) of CKD, estimated glomerular filtration rate (eGFR) and blood urea nitrogen (BUN) to perform proteome-wide association study (PWAS), Mendelian Randomization and Bayesian colocalization analyses. The single-cell RNA sequencing data of healthy human and mouse kidneys were analyzed to explore the cell-type specificity of identified genes. Functional enrichment analysis was conducted to investigate the involved signaling pathways. The PWAS identified 22 plasma proteins significantly associated with CKD. Of them, the significant associations of three proteins (INHBC, LMAN2, and SNUPN) were replicated in the GWASs of eGFR, and BUN. Mendelian Randomization analyses showed that INHBC and SNUPN were causally associated with CKD, eGFR, and BUN. The Bayesian colocalization analysis identified shared causal variants for INHBC in CKD, eGFR, and BUN (all PP4 > 0.75). The single-cell RNA sequencing revealed that the INHBC gene was sparsely scattered within the kidney cells. This proteomic study revealed that INHBC, LMAN2, and SNUPN may be involved in the pathogenesis of CKD, which represent novel therapeutic targets and warrant further exploration in future research.

Plasma Proteomics of Exercise Blood Pressure and Incident Hypertension

Importance

Blood pressure response during acute exercise (exercise blood pressure [EBP]) is associated with the future risk of hypertension and cardiovascular disease (CVD). Biochemical characterization of EBP could inform disease biology and identify novel biomarkers of future hypertension.

Objective

To identify protein markers associated with EBP and test their association with incident hypertension.

Design, Setting, and Participants

This study assayed 4977 plasma proteins in 681 healthy participants (from 763 assessed) of the Health, Risk Factors, Exercise Training and Genetics (HERITAGE; data collection from January 1993 to December 1997 and plasma proteomics from January 2019 to January 2020) Family Study at rest who underwent 2 cardiopulmonary exercise tests. Individuals were free of CVD at the time of recruitment. Individuals with resting SBP ≥160 mm Hg or DBP ≥100 mm Hg or taking antihypertensive drug therapy were excluded from the study. The association between resting plasma protein levels to both resting BP and EBP was evaluated. Proteins associated with EBP were analyzed for their association with incident hypertension in the Framingham Heart Study (FHS; n = 1177) and validated in the Jackson Heart Study (JHS; n = 772) and Multi-Ethnic Study of Atherosclerosis (MESA; n = 1367). Proteins associated with incident hypertension were tested for putative causal links in approximately 700 000 individuals using cis-protein quantitative loci mendelian randomization (cis-MR). Data were analyzed from January 2023 to January 2024.

Exposures

Plasma proteins.

Main Outcomes and Measures

EBP was defined as the BP response during a fixed workload (50 W) on a cycle ergometer. Hypertension was defined as BP ≥140/90 mm Hg or taking antihypertensive medication.

Results

Among the 681 participants in the HERITAGE Family Study, the mean (SD) age was 34 (13) years; 366 participants (54%) were female; 238 (35%) were self-reported Black and 443 (65%) were self-reported White. Proteomic profiling of EBP revealed 34 proteins that would not have otherwise been identified through profiling of resting BP alone. Transforming growth factor β receptor 3 (TGFBR3) and prostaglandin D2 synthase (PTGDS) had the strongest association with exercise systolic BP (SBP) and diastolic BP (DBP), respectively (TGFBR3: exercise SBP, β estimate, -3.39; 95% CI, -4.79 to -2.00; P = 2.33 × 10-6; PTGDS: exercise DBP β estimate, -2.50; 95% CI, -3.29 to -1.70; P = 1.18 × 10-9). In fully adjusted models, TGFBR3 was inversely associated with incident hypertension in FHS, JHS, and MESA (hazard ratio [HR]: FHS, 0.86; 95% CI, 0.75-0.97; P = .01; JHS, 0.87; 95% CI, 0.77-0.97; P = .02; MESA, 0.84; 95% CI, 0.71-0.98; P = .03; pooled cohort, 0.86; 95% CI, 0.79-0.92; P = 6 × 10-5). Using cis-MR, genetically predicted levels of TGFBR3 were associated with SBP, hypertension, and CVD events (SBP: β, -0.38; 95% CI, -0.64 to -0.11; P = .006; hypertension: odds ratio [OR], 0.99; 95% CI, 0.98-0.99; P < .001; heart failure with hypertension: OR, 0.86; 95% CI, 0.77-0.97; P = .01; CVD: OR, 0.84; 95% CI, 0.77-0.92; P = 8 × 10-5; cerebrovascular events: OR, 0.77; 95% CI, 0.70-0.85; P = 5 × 10-7).

Conclusions and Relevance

Plasma proteomic profiling of EBP identified a novel protein, TGFBR3, which may protect against elevated BP and long-term CVD outcomes.

A Mendelian randomization study investigating the causal relationships between inflammation and immunoglobulin A nephropathy

BACKGROUND

Immunoglobulin A nephropathy (IgAN) is an autoimmune disease characterized by the production of galactose‑deficient IgA1 (Gd‑IgA1) and the deposition of immune complexes in the kidney. Exploring the landscape of immune dysregulation in IgAN is valuable for pathogenesis and disease treatment. We conducted Mendelian randomization (MR) to assess the causal correlations between inflammation and IgAN.

METHODS

Based on available genetic datasets, we investigated potential causal links between inflammation and the risk of IgAN using two-sample MR. We used genome-wide association study (GWAS) summary statistics of 5 typical inflammation markers, 41 inflammatory cytokines, and 731 immune cell signatures, accessed from the public GWAS Catalog. The primary method employed for MR analysis was Inverse Variance Weighted (IVW). To confirm consistency across results, four supplementary MR methods were also conducted: MR-Egger, Weighted Median, Weighted Mode, and Simple Mode. To assess pleiotropy, we used the MR-Egger regression intercept test and Mendelian Randomization Pleiotropy RESidual Sum and Outlier (MR-PRESSO) test. Cochrane's Q statistic was applied to evaluate heterogeneity. Additionally, the stability of the MR findings was verified through the leave-one-out sensitivity analysis.

RESULTS

This study revealed that interleukin-7 (IL-7) and stem cell growth factor beta (SCGF-β) were possibly associated with the risk of IgAN according to the IVW approach, with estimated odds ratios (OR) of 1.059 (95 % confidence interval [CI] 1.015 to 1.104, P = 0.008) and 1.043 (95 % CI 1.002 to 1.085, P = 0.037). Five immune traits were identified that might be linked to IgAN risk, each with P-values below 0.01, including natural killer T %T cell (OR = 1.058, 95 % CI: 1.020 to 1.097, P = 0.002), natural killer T %lymphocyte (OR = 1.055, 95 % CI: 1.016 to 1.096, P = 0.006), CD25++ CD8+ T cell %T cell (OR = 1.057, 95 % CI: 1.016 to 1.099, P = 0.006), CD3 on effector memory CD4+ T cell (OR = 1.045, 95 % CI: 1.019 to 1.071, P = 0.001), and CD3 on CD28+ CD45RA+ CD8+ T cell (OR = 1.042, 95 % CI: 1.016 to 1.068, P = 0.001). CD4 on central memory CD4+ T cell might be a protective factor for IgAN (OR = 0.922, 95 % CI: 0.875 to 0.971, P = 0.002). Moreover, IgAN may be implicated in a high risk of elevated granulocyte colony-stimulating factor (G-CSF) (OR = 1.114, 95 % CI 1.002 to 1.239, P = 0.046).

CONCLUSION

Our study revealed exposures among typical inflammation markers, inflammatory cytokines, and immune cell signatures that may potentially linked to IgAN risk by MR analysis. This insight may advance our understanding of the etiology of IgAN and support the development of targeted therapeutic strategies.

Integrating plasma protein-centric multi-omics to identify potential therapeutic targets for pancreatic cancer

BACKGROUND

Deciphering the role of plasma proteins in pancreatic cancer (PC) susceptibility can aid in identifying novel targets for diagnosis and treatment.

METHODS

We examined the relationship between genetically determined levels of plasma proteins and PC through a systemic proteome-wide Mendelian randomization (MR) analysis utilizing cis-pQTLs from multiple centers. Rigorous sensitivity analyses, colocalization, reverse MR, replications with varying instrumental variable selections and additional datasets, as well as subsequent meta-analysis, were utilized to confirm the robustness of significant findings. The causative effect of corresponding protein-coding genes' expression and their expression pattern in single-cell types were then investigated. Enrichment analysis, between-protein interaction and causation, knock-out mice models, and mediation analysis with established PC risk factors were applied to indicate the pathogenetic pathways. These candidate targets were ultimately prioritized upon druggability and potential side effects predicted by a phenome-wide MR.

RESULTS

Twenty-one PC-related circulating proteins were identified in the exploratory phase with no evidence for horizontal pleiotropy or reverse causation. Of these, 11 were confirmed in a meta-analysis integrating external validations. The causality at a transcription level was repeated for neutrophil elastase, hydroxyacylglutathione hydrolase, lipase member N, protein disulfide-isomerase A5, xyloside xylosyltransferase 1. The carbohydrate sulfotransferase 11 and histo-blood group ABO system transferase exhibited high-support genetic colocalization evidence and were found to affect PC carcinogenesis partially through modulating body mass index and type 2 diabetes, respectively. Approved drugs have been established for eight candidate targets, which could potentially be repurposed for PC therapies. The phenome-wide investigation revealed 12 proteins associated with 51 non-PC traits, and interference on protein disulfide-isomerase A5 and cystatin-D would increase the risk of other malignancies.

CONCLUSIONS

By employing comprehensive methodologies, this study demonstrated a genetic predisposition linking 21 circulating proteins to PC risk. Our findings shed new light on the PC etiology and highlighted potential targets as priorities for future efforts in early diagnosis and therapeutic strategies of PC.

Exploring the role of mitochondrial proteins SIRT5 and MRPL33 through Mendelian randomization in primary biliary cholangitis

BACKGROUND

Primary biliary cholangitis (PBC) is an autoimmune liver disease characterized by elevated serum antimitochondrial antibody levels in 90-95 % of cases. However, the exact causal relationship between mitochondrial proteins and PBC remains unclear. This study aims to investigate and clarify this relationship.

METHODS

Genome-wide association data for mitochondrial proteins and PBC were obtained from public databases. The assessment of causal relationships between exposures and outcomes employed the Inverse Variance Weighted (IVW) method, MR Egger regression, and Weighted Median. Sensitivity analyses were systematically carried out to appraise the robustness of the Mendelian Randomization (MR) findings.

RESULTS

The analysis revealed two mitochondrial proteins exhibiting a causal relationship with PBC. Elevated SIRT5 levels demonstrated a positive correlation with an augmented susceptibility to PBC in the IVW approach (odds ratio, OR: 1.2907, 95 % CI: 1.062-1.568, p = 0.0102). Conversely, increased MRPL33 levels were associated with a decreased risk of PBC (OR: 0.8957, 95 % CI: 0.807-0.993, p = 0.0376). Sensitivity analysis corroborated these findings consistently.

CONCLUSION

This investigation advances the notion of a potential causal association between elevated SIRT5 levels and an increased risk of PBC, alongside a decreased risk of PBC linked to elevated MRPL33 levels. The identified mitochondrial proteins may serve as viable biomarkers, offering pertinent insights for the understanding and addressing of PBC.

Genetic analyses of inflammatory polyneuropathy and chronic inflammatory demyelinating polyradiculoneuropathy identified candidate genes

Chronic inflammatory demyelinating polyneuropathy (CIDP) is a rare, immune-mediated disorder in which an aberrant immune response causes demyelination and axonal damage of the peripheral nerves. Genetic contribution to CIDP is unclear and no genome-wide association study (GWAS) has been reported so far. In this study, we aimed to identify CIDP-related risk loci, genes, and pathways. We first focused on CIDP, and 516 CIDP cases and 403,545 controls were included in the GWAS analysis. We also investigated genetic risk for inflammatory polyneuropathy (IP), in which we performed a GWAS study using FinnGen data and combined the results with GWAS from the UK Biobank using a fixed-effect meta-analysis. A total of 1,261 IP cases and 823,730 controls were included in the analysis. Stratified analyses by gender were performed. Mendelian randomization (MR), colocalization, and transcriptome-wide association study (TWAS) analyses were performed to identify associated genes. Gene-set analyses were conducted to identify associated pathways. We identified one genome-wide significant locus at 20q13.33 for CIDP risk among women, the top variant located at the intron region of gene CDH4. Sex-combined MR, colocalization, and TWAS analyses identified three candidate pathogenic genes for CIDP and five genes for IP. MAGMA gene-set analyses identified a total of 18 pathways related to IP or CIDP. Sex-stratified analyses identified three genes for IP among males and two genes for IP among females. Our study identified suggestive risk genes and pathways for CIDP and IP. Functional analyses should be conducted to further confirm these associations.

Potential Biomarkers in Cerebrospinal Fluid and Plasma for Dementia

Background

The identification of biomarkers for different dementias in plasma and cerebrospinal fluid (CSF) has made substantial progress. However, they are observational studies, and there remains a lack of research on dementias with low incidence rates.

Objective

We performed a comprehensive Mendelian randomization to identify potential biomarkers for different dementia type.

Methods

The summary-level datasets encompassed 734 plasma and 154 cerebrospinal fluid proteins sourced from recently published genome-wide association studies (GWAS). Summary statistics for different dementias, including any dementia (refering to any type of dementia symptoms, 218,792 samples), Alzheimer's disease (AD, 63,926 samples), vascular dementia (212,389 samples), frontotemporal dementia (3,024 samples), dementia with Lewy bodies (DLB, 6,618 samples), and dementia in Parkinson's disease (216,895 samples), were collected from large GWAS. The primary method is inverse variance weighting, with additional sensitivity analyses conducted to ensure the robustness of the findings.

Results

The molecules released into CSF, namely APOE2 for any dementia, APOE2 and Siglec-3 for AD, APOE2 for vascular dementia, and APOE2 for DLB, might be potential biomarkers. CD33 for AD and SNCA for DLB in plasma could be promising biomarkers.

Conclusions

This is the first study to integrate plasma and CSF proteins to identify potential biomarkers for different dementias.

Bidirectional Two-Sample Mendelian Randomization Study of Immunoglobulin G N-Glycosylation and Senescence-Associated Secretory Phenotype

Observational studies revealed changes in Immunoglobulin G (IgG) N-glycosylation during the aging process. However, it lacks causal insights and remains unclear in which direction causal relationships exist. The two-sample bidirectional Mendelian randomization (MR) design was adopted to explore causal associations between IgG N-glycans and the senescence-associated secretory phenotype (SASP). Inverse variance weighted (IVW) and Wald ratio methods were used as the main analyses, supplemented by sensitivity analyses. Forward MR analyses revealed causal associations between the glycan peak (GP) and SASP, including GP6 (odds ratio [OR] = 0.428, 95% confidence interval [CI] = 0.189-0.969) and GP17 (OR = 0.709, 95%CI = 0.504-0.995) with growth/differentiation factor 15 (GDF15), GP19 with an advanced glycosylation end-product-specific receptor (RAGE) (OR = 2.142, 95%  CI  = 1.384-3.316), and GP15 with matrix metalloproteinase 2 (MMP2) (OR = 1.136, 95%  CI =1.008-1.282). The reverse MR indicated that genetic liability to RAGE was associated with increased levels of GP17 (OR = 1.125, 95%  CI  = 1.003-1.261) and GP24 (OR = 1.222, 95%  CI  = 1.046-1.428), while pulmonary and activation-regulated chemokines (PARC) exhibited causal associations with GP10 (OR = 1.269, 95%  CI  = 1.048-1.537) and GP15 (OR = 1.297, 95%  CI = 1.072-1.570). The findings provided suggested evidence on the bidirectional causality between IgG N-glycans and SASP, which might reveal potential regulatory mechanisms.

Association between cathepsins and benign prostate diseases: a bidirectional two-sample Mendelian randomization study

Objectives

The relationship between cathepsins and prostate cancer (PCa) has been reported. However, there is a lack of research on cathepsins and benign prostate diseases (BPDs). This study investigated the potential genetic link between cathepsins and BPDs through the utilization of Mendelian randomization (MR) analysis to determine if a causal relationship exists.

Methods

Publicly accessible summary statistics on BPDs were obtained from FinnGen Biobank. The data comprised 149,363 individuals, with 30,066 cases and 119,297 controls for BPH, and 123,057 individuals, with 3,760 cases and 119,297 controls for prostatitis. The IEU OpenGWAS provided the Genome-wide association data on ten cathepsins. To evaluate the causal relationship between BPDs and cathepsins, five distinct MR analyses were employed, with the primary method being the inverse variance weighted (IVW) approach. Additionally, sensitivity analyses were conducted to examine the horizontal pleiotropy and heterogeneity of the findings.

Results

The examination of IVW MR findings showed that cathepsin O had a beneficial effect on BPH (IVW OR=0.94, 95% CI 0.89-0.98, P=0.0055), while cathepsin X posed a threat to prostatitis (IVW OR=1.08, 95% CI 1.00-1.16, P=0.047). Through reverse MR analysis, it was revealed that prostatitis had an adverse impact on cathepsin V (IVW OR=0.89, 95% CI 0.80-0.99, P=0.035), while no favorable association was observed between BPH and cathepsins. The results obtained from MR-Egger, weighted median, simple mode, and weighted mode methods were consistent with the findings of the IVW approach. Based on sensitivity analyses, heterogeneity, and horizontal pleiotropy are unlikely to distort the results.

Conclusion

This study offers the initial evidence of a genetic causal link between cathepsins and BPDs. Our findings revealed that cathepsin O was beneficial in preventing BPH, whereas cathepsin X posed a potential threat to prostatitis. Additionally, prostatitis negatively affected cathepsin V level. These three cathepsins could be targets of diagnosis and treatment for BPDs, which need further research.

Independent causal effect of migraines on Alzheimer's disease risk: a multivariate Mendelian randomization study

Background

The observational studies investigated the impact of migraine on Alzheimer's Disease (AD). However, these findings were limited by confounding factors and reverse causation, leading to contradictory results.

Methods

We utilized Univariable Mendelian Randomization (UVMR) to explore the link between migraine (13,971 cases/470,627 controls) and AD risk (Bellenguez et al., 39,106 cases/46,828 controls; FinnGen, 111,471 cases/111,471 controls). Meta-analysis was performed for comprehensive synthesis. Employing Multivariable Mendelian Randomization (MVMR), we created models incorporating migraine and 35 potential AD risk factors, examining migraine's independent impact on AD onset risk under considering these factors.

Results

The meta-analysis of inverse variance weighted MR results, combining data from Bellenguez et al. (odds ratio (OR) [95% confidence interval (CI)]: 1.5717 [1.1868-2.0814], p = 0.0016) and FinnGen (OR [95% CI]: 1.2904 [0.5419-3.0730], p = 0.5646), provided evidence for a causal relationship between genetically predicted migraine and the heightened risk of AD occurrence (OR [95% CI]: 1.54 [1.18, 2.00], p < 0.01). After adjusting for Diastolic blood pressure (OR [95% CI]: 1.4120 [0.8487-2.3493], p = 0.1840) and Tumor necrosis factor alpha (OR [95% CI]: 1.2411 [0.8352-1.8443], p = 0.2852), no discernible association was detected between migraine and the risk of AD.

Conclusion

This study offers compelling evidence indicating a significant correlation between genetically predicted migraine and an elevated risk of AD.

Cathepsins and Parkinson's disease: insights from Mendelian randomization analyses

Background

Parkinson's disease (PD), the second most prevalent neurodegenerative condition, has a multifaceted etiology. Cathepsin-cysteine proteases situated within lysosomes participate in a range of physiological and pathological processes, including the degradation of harmful proteins. Prior research has pointed towards a potential link between cathepsins and PD; however, the precise causal relationship between the cathepsin family and PD remains unclear.

Methods

This study employed univariate and multivariate Mendelian randomization (MR) analyses to explore the causal relationship between the nine cathepsins and Parkinson's disease (PD) risk. For the primary analysis, genome-wide association study (GWAS) summary statistics for the plasma levels of the nine cathepsins and PD was obtained from the INTERVAL study and the International Parkinson's Disease Genomics Consortium. GWAS for PD replication analysis were obtained from the FinnGen consortium, and a meta-analysis was performed for the primary and replication analyses to evaluate the association between genetically predicted cathepsin plasma levels and PD risk. After identifying significant MR estimates, genetic co-localization analyses were conducted to determine whether shared or distinct causal variants influenced both cathepsins and PD.

Results

Elevated cathepsin B levels were associated with a decreased risk of PD in univariate MR analysis (odds ratio [OR] = 0.890, 95% confidence interval [CI]: 0.831-0.954, pFDR = 0.009). However, there was no indication that PD affected cathepsin B levels (OR = 0.965, 95% CI: 0.858-1.087, p = 0.852). In addition, after adjusting for the remaining cathepsins, cathepsin B levels independently and significantly contributed to the reduced risk of PD in multivariate MR analysis (OR = 0.887, 95% CI: 0.823-0.957, p = 0.002). The results of the replication MR analysis with the FinnGen GWAS for PD (OR = 0.921, 95% CI: 0.860-0.987, p = 0.020) and meta-analysis (OR = 0.905, 95% CI: 0.862-0.951, p < 0.001) were consistent with those of the primary analysis. Colocalization analysis did not provide any evidence of a shared causal variant between cathepsins and PD (PP.H4.abf = 0.005).

Conclusion

This genetic investigation supports the hypothesis that cathepsin B exerts a protective effect against PD. The quantification of cathepsin B levels could potentially serve as a predictive biomarker for susceptibility to PD, providing new insights into the pathomechanisms of the disease and possible interventions.

Evaluating the Causal Effects of Circulating Proteome on the Risk of Sepsis and Related Outcomes

The current investigation deployed Mendelian randomization (MR) to elucidate the causal relationship between circulating proteins and sepsis. A rigorous two-sample MR analysis evaluated the effect of plasma proteins on the sepsis susceptibility. To affirm the integrity of MR findings, a suite of supplementary analyses, including Bayesian colocalization, Steiger filtering, the assessment of protein-altering polymorphisms, and the correlation between expression quantitative trait loci and protein quantitative trait loci (pQTLs), was employed. The study further integrated the examination of protein-protein interactions and pathway enrichment, along with the identification of pharmacologically actionable targets, to advance our comprehension and outline potential sepsis therapies. Subsequent analyses leveraging cis-pQTLs within MR studies unveiled noteworthy relationships: 94 specific proteins exhibited significant links with sepsis-related 28 day mortality, while 96 distinct proteins correlated with survival outcomes in sepsis. Furthermore, incorporating both cis- and trans-pQTLs in MR investigations revealed more comprehensive findings, associating 201 unique proteins with sepsis-related 28 day mortality and 199 distinct proteins with survival outcomes in sepsis. Markedly, colocalization analyses confirmed that eight of these proteins exhibited prominent evidence for colocalization, emphasizing their potential criticality in sepsis pathophysiology. Further in silico analyses were conducted to delineate putative regulatory networks and to highlight prospective drug targets among these proteins. Employing the MR methodology has shed light on plasma proteins implicated in the etiopathogenesis of sepsis. This novel approach unveiled numerous biomarkers and targets, providing a scientific rationale for the development of new therapeutic strategies and prophylactic measures against sepsis.

Unraveling the role of plasma proteins in dementia: insights from two cohort studies in the UK, with causal evidence from Mendelian randomization

Population-based proteomics offer a groundbreaking avenue to predict dementia onset. This study employed a proteome-wide, data-driven approach to investigate protein-dementia associations in 229 incident all-cause dementia (ACD) among 3,249 participants from the English Longitudinal Study of Ageing (ELSA) over a median 9.8-year follow-up, then validated in 1,506 incident ACD among 52,745 individuals from the UK Biobank (UKB) over median 13.7 years. NEFL and RPS6KB1 were robustly associated with incident ACD; MMP12 was associated with vascular dementia in ELSA. Additional markers EDA2R and KIM1 (HAVCR1) were identified from sensitivity analyses. Combining NEFL and RPS6KB1 with other factors yielded high predictive accuracy (area under the curve (AUC)=0.871) for incident ACD. Replication in the UKB confirmed associations between identified proteins with various dementia subtypes. Results from reverse Mendelian Randomization also supported the role of several proteins as early dementia biomarkers. These findings underscore proteomics' potential in identifying novel risk screening targets for dementia.

The impact of cannabis use on erectile dysfunction and sex hormones: a Mendelian randomization analysis

Previous study has highlighted an association between cannabis use (CU) and an increased risk of erectile dysfunction (ED), potentially due to indirect effects on sex hormonal balance. However, the evidence remains controversial, and the causal relationship is unclear. This study utilized genome-wide association study (GWAS) data to investigate the causal relationships between cannabis use disorder (CUD), lifetime cannabis use (LCU), and ED, as well as levels of sex hormones including estradiol (E2), bioavailable testosterone (BT), follicle-stimulating hormone (FSH), and luteinizing hormone (LH) through Mendelian randomization (MR) analysis. The primary method of analysis was the inverse variance weighted (IVW) method. Data from the FinnGen and UK Biobank were used for replication and meta-analysis. The results indicated no causal relationship between genetically predicted CUD (OR = 0.97, 95% CI 0.87-1.10, P = 0.66) and LCU (OR = 1.13, 95% CI 0.84-1.50, P = 0.42) with the risk of ED. The meta-analysis provided consistent evidence (P > 0.05). No causal relationships were found between CUD and LCU with E2(CUD: β = 0.00, 95% CI 0.00-0.01, P = 0.37; LCU: β = 0.00, 95% CI -0.02-0.01, P = 0.62), BT (CUD: β = 0.00, 95% CI -0.03-0.02, P = 0.90; LCU: β = 0.02, 95% CI -0.04-0.09, P = 0.46), FSH (CUD: β = 0.01, 95% CI -0.18-0.20, P = 0.92; LCU: β = 0.01, 95% CI -0.44-0.47, P = 0.95), and LH (CUD: β = 0.01, 95% CI -0.18-0.21, P = 0.90; LCU: β = 0.13, 95% CI -0.22-0.49, P = 0.46). Sensitivity analyses detected no evidence of horizontal pleiotropy or heterogeneity, ensuring the robustness of the results. In conclusion, this MR analysis did not provide evidence supporting a causal relationship between CU and ED or sex hormone levels.

Causal Association of Cytokines and Growth Factors with Stroke and Its Subtypes: a Mendelian Randomization Study

Cytokines and growth factors contribute to nerve growth and angiogenesis and are associated with the development of vascular disease. This Mendelian randomization (MR) study was designed to examine the causal relationship between factors associated with stem cell paracrine mechanisms and with stroke and its subtypes. We used pooled statistics on cytokine levels from three studies (INTERIAL, Olink Proseek CVD array, and KORA) encompassing 7795 participants in Europe. Data for stroke and its subtypes were pooled from these European populations (40,585 cases and 406,111 controls) in a multiprogenitor genome-wide association study (GWAS). MR was performed using established analytical methods, including inverse variance weighting (IVW), weighted median (WM), and MR-Egger. Genetically determined high IGF-1 levels were found to associate negatively with risk of stroke, ischemic stroke (large-artery atherosclerosis), and ischemic stroke (cardiogenic embolism). Meanwhile, high IL-13 levels had a positive causal relationship with ischemic stroke (large-artery atherosclerosis). An additional 27 cytokines were found to have a causal association with stroke or its subtypes. However, these results should be interpreted with caution given that the power efficacy was <80%. This MR study supports the concept of a causal relationship of 29 cytokines with stroke or its subtypes. Our genetic analysis provides new insights into stroke prevention and treatment by demonstrating an association of stem cell paracrine-related cytokines with stroke risk.

Novel insights into the complex interplay of immune dysregulation and inflammatory biomarkers in preeclampsia and fetal growth restriction: A two-step Mendelian randomization analysis

Background

The relationship between genetic immune dysregulation and the occurrence of preeclampsia (PE) or PE with fetal growth restriction (PE with FGR) has yielded inconsistent findings, and the underlying mediators of this association remain elusive. We aimed to explore the causal impact of genetic immune dysregulation on the risk of PE or PE with FGR and to elucidate the role of specific transcriptomes in mediating this relationship.

Methods

A two-step Mendelian randomization (MR) analysis was performed to explore the link between immune dysregulation and PE or PE with FGR, as well as to identify potential inflammatory biomarkers that act as mediators. GWAS summary data for outcomes were obtained from the FinnGen dataset. The analyses encompassed five systemic immune-associated diseases, four chronic genital inflammatory diseases, and thirty-one inflammatory biomarkers. Summary-data-based MR (SMR) and HEIDI analysis were conducted to test whether the effect size of single nucleotide polymorphisms (SNPs) on outcomes was mediated by the expression of immune-associated genes.

Results

The primary univariable analysis revealed a significant positive correlation between systemic lupus erythematosus (SLE), type 1 diabetes (T1D), type 2 diabetes (T2D), and rheumatoid arthritis (RA) with the risk of PE or PE with FGR. Surprisingly, a counterintuitive finding showed a significant negative association between endometriosis of pelvic peritoneum (EMoP) and the risk of PE with FGR. None of the inflammatory factors had a causal relationship with PE or PE with FGR. However, there was a significant association between lymphocyte count and the risk of PE with FGR. Within the lymphocyte subset, both the proportion of Natural Killer (NK) cells and absolute counts of naïve CD4+ T cells demonstrated significant effects on the risk of PE with FGR. Two-step MR analysis underscored the genetically predicted lymphocyte count as a significant mediator between T1D and PE with FGR. Additionally, SMR analysis indicated the potential involvement of SH2B3 in the occurrence of PE with FGR.

Conclusions

Our findings provided substantial evidence of the underlying causal relationship between immune dysregulation and PE or PE with FGR and some of these diseases proved to accelerate immune cells disorders and then contribute to the risk of incident PE or PE with FGR.

Comprehensive single-cell transcriptomic profiling reveals molecular subtypes and prognostic biomarkers with implications for targeted therapy in esophageal squamous cell carcinoma

BACKGROUND

Esophageal squamous cell carcinoma (ESCC) is a genetically heterogeneous disease with poor clinical outcomes. Identification of biomarkers linked to DNA replication stress may enable improved prognostic risk stratification and guide therapeutic decision making. We performed integrated single-cell RNA sequencing and computational analyses to define the molecular determinants and subtypes underlying ESCC heterogeneity.

METHODS

Single-cell RNA sequencing was performed on ESCC samples and analyzed using Seurat. Differential gene expression analysis was used to identify esophageal cell phenotypes. DNA replication stress-related genes were intersected with single-cell differential expression data to identify potential prognostic genes, which were used to generate a DNA replication stress (DRS) score. This score and associated genes were evaluated in survival analysis. Putative prognostic biomarkers were evaluated by Cox regression and consensus clustering. Mendelian randomization analyses assessed the causal role of PRKCB.

RESULTS

High DRS score associated with poor survival. Four genes (CDKN2A, NUP155, PPP2R2A, PRKCB) displayed prognostic utility. Three molecular subtypes were identified with discrete survival and immune properties. A 12-gene signature displayed robust prognostic performance. PRKCB was overexpressed in ESCC, while PRKCB knockdown reduced ESCC cell migration.

CONCLUSIONS

This integrated single-cell sequencing analysis provides new insights into the molecular heterogeneity and prognostic determinants underlying ESCC. The findings identify potential prognostic biomarkers and a gene expression signature that may enable improved patient risk stratification in ESCC. Experimental validation of the role of PRKCB substantiates the potential clinical utility of our results.

Metalloproteins as risk factors for osteoarthritis: improving and understanding causal estimates using Mendelian randomization

Osteoarthritis (OA) is one of the most prevalent musculoskeletal disorders and a primary cause of pain and disability among the elderly population. Research on the relationship between metalloproteins (MPs) and OA is limited, and causality remains unclear. Our objective is to utilize Mendelian randomization (MR) to explore the possible causal relationship between MPs and OA. The data on MPs were derived from a Genome-Wide Association Study (GWAS) analysis involving 3301 samples. The GWAS data for OA were obtained from an analysis involving 462,933 European individuals. In this study, a variety of two-sample Mendelian randomization methods (two-sample MR) to evaluate the causal effect of MPs on OA, including inverse variance weighted method (IVW), MR-Egger method, weighted median method (WM), simple mode, weight mode, and Wald ratio. The primary MR analysis using the IVW method reveals a significant negative correlation between Metallothionein-1F (MT-1F), zinc finger protein 134 (ZNF134), calcium/calmodulin-dependent protein kinase type 1D (CAMK1D), and EF-hand calcium-binding domain-containing protein 14 (EFCAB14) with the occurrence of osteoarthritis (OA) (p value < 0.05). However, no causal relationship was observed in the opposite direction between these MPs and OA. Notably, even in combined models accounting for confounding factors, the negative association between these four MPs and OA remained significant. Sensitivity analysis demonstrated no evidence of horizontal pleiotropy or heterogeneity, and leave-one-out analysis confirmed the robustness of the results. In this study, we have established a conspicuous association between four distinct MPs and OA. This discovery augments our understanding of potential avenues for the diagnosis and treatment of this condition. Key Points • The MR method was employed to assess the relationship between MPs and OA. • A total of four types of MPs have demonstrated inhibitory effects on the occurrence of OA.

Genetically predicted systemic inflammation and the risk of atrial fibrillation: A bidirectional two-sample Mendelian randomization study

Background

Systemic inflammation has been proposed to be associated with the incidence of atrial fibrillation (AF), but whether it is a cause or a consequence of AF remains uncertain. We sought to explore the causal associations between systemic inflammation and AF using bidirectional Mendelian randomization (MR) analysis.

Methods

Independent genetic variants strongly associated with AF were selected as instrumental variables from the largest genome-wide association study (GWAS) with up to 1,030,836 individuals. Regarding inflammation traits, genetic associations with 41 inflammatory cytokines and 5 inflammatory biomarkers were obtained from their corresponding GWASs databases. Effect estimates were primarily evaluated using the inverse-variance weighted (IVW) method, supplemented by sensitivity analyses using MR-Egger, weighted median, and MR-PRESSO methods.

Results

In our initial MR analyses, we observed suggestive associations of genetically predicted interleukin-17 (IL-17), interleukin-2 receptor subunit alpha (IL-2rα), and procalcitonin (PCT) with AF. One standard deviation (SD) increase in IL-17, IL-2rα, and PCT caused an increase in AF risk by 6.3 % (OR 1.063, 95 %CI 1.011---1.118, p = 0.018), 4.9 % (OR 1.049, 95 %CI 1.007---1.094, p = 0.023) and 3.4 % (OR 1.034, 95 %CI 1.005---1.064, p = 0.022), respectively. Furthermore, our reverse MR analyses indicated that genetically predicted AF contributed to a suggestive increase in the levels of macrophage inflammatory protein-1β (MIP1β) (β 0.055, 95 %CI 0.006 to 0.103, p = 0.028), while a decrease in the levels of fibrinogen (Fbg) (β -0.091, 95 %CI -0.140 to -0.041, p < 0.001), which remained significant after multiple test correction.

Conclusions

Our MR study identified several inflammatory biomarkers with suggestive causal associations regarding the upstream and downstream regulation of AF occurrence, offering new insights for therapeutic exploitation of AF. Further research is required to validate the underlying link between systemic inflammation and AF in larger cohorts.

Predicting Preterm Birth Using Proteomics

The complexity of preterm birth (PTB), both spontaneous and medically indicated, and its various etiologies and associated risk factors pose a significant challenge for developing tools to accurately predict risk. This review focuses on the discovery of proteomics signatures that might be useful for predicting spontaneous PTB or preeclampsia, which often results in PTB. We describe methods for proteomics analyses, proteomics biomarker candidates that have so far been identified, obstacles for discovering biomarkers that are sufficiently accurate for clinical use, and the derivation of composite signatures including clinical parameters to increase predictive power.

Association between oxidative stress and chronic orofacial pain and potential druggable targets: Evidence from a Mendelian randomization study

BACKGROUND

Oxidative stress indicators affect chronic orofacial pain (COFP), but how to reduce these effects is uncertain.

OBJECTIVES

11 oxidative stress biomarkers were collected as exposures, while four forms of COFP were chosen as outcomes for Mendelian randomization (MR) study.

METHODS

The effect estimates between oxidative stress and COFP were calculated using inverse variance-weighted MR (IVW-MR). Then, functional mapping and annotation (FUMA) was utilized in order to carry out SNP-based functional enrichment analyses. In addition, the IVW-MR method was applied to combine effect estimates when using genetic variants associated with oxidative stress biomarkers as an instrument for exploring potential druggable targets.

RESULTS

The results indicated that oxidative stress biomarkers (causal OR of uric acid (UA), 0.998 for myofascial pain, 95% CI 0.996-1.000, p < .05; and OR of glutathione transferase (GST), 1.002 for dentoalveolar pain, 95% CI 1.000-1.003, p < .05) were significantly linked with the probability of COFP. Functional analysis also demonstrated that UA and myofascial pain genes were prominent in nitrogen and uracil metabolism, while GST and dentoalveolar pain genes were enriched in glutathione metabolism. Also, the study provided evidence that solute carrier family 2 member 9 (SLC2A9) and glutathione S-transferase alpha 2 (GSTA2) cause discomfort in the myofascial pain (OR = 1.003, 95% CI 1.000-1.006; p < .05) and dentoalveolar region (OR = 1.001, 95% CI 1.000-1.002; p < .05), respectively.

CONCLUSIONS

In conclusion, this MR study indicates that genetically predicted myofascial pain was significantly associated with decreased UA and dentoalveolar pain was significantly associated with increased GST level. SLC2A9 inhibitor and GSTA2 inhibitor were novel chronic orofacial pain therapies and biomarkers, but clinical trials are called to examine if these oxidative biomarkers have the protective effect against orofacial pain, and further research are needed to explore the underlying mechanisms.

A genome-wide association study of mass spectrometry proteomics using the Seer Proteograph platform

Genome-wide association studies (GWAS) with proteomics are essential tools for drug discovery. To date, most studies have used affinity proteomics platforms, which have limited discovery to protein panels covered by the available affinity binders. Furthermore, it is not clear to which extent protein epitope changing variants interfere with the detection of protein quantitative trait loci (pQTLs). Mass spectrometry-based (MS) proteomics can overcome some of these limitations. Here we report a GWAS using the MS-based Seer Proteograph™ platform with blood samples from a discovery cohort of 1,260 American participants and a replication in 325 individuals from Asia, with diverse ethnic backgrounds. We analysed 1,980 proteins quantified in at least 80% of the samples, out of 5,753 proteins quantified across the discovery cohort. We identified 252 and replicated 90 pQTLs, where 30 of the replicated pQTLs have not been reported before. We further investigated 200 of the strongest associated cis-pQTLs previously identified using the SOMAscan and the Olink platforms and found that up to one third of the affinity proteomics pQTLs may be affected by epitope effects, while another third were confirmed by MS proteomics to be consistent with the hypothesis that genetic variants induce changes in protein expression. The present study demonstrates the complementarity of the different proteomics approaches and reports pQTLs not accessible to affinity proteomics, suggesting that many more pQTLs remain to be discovered using MS-based platforms.

Cathepsins and cancer risk: a Mendelian randomization study

Background

Previous observational epidemiological studies reported an association between cathepsins and cancer, however, a causal relationship is uncertain. This study evaluated the causal relationship between cathepsins and cancer using Mendelian randomization (MR) analysis.

Methods

We used publicly available genome-wide association study (GWAS) data for bidirectional MR analysis. Inverse variance weighting (IVW) was used as the primary MR method of MR analysis.

Results

After correction for the False Discovery Rate (FDR), two cathepsins were found to be significantly associated with cancer risk: cathepsin H (CTSH) levels increased the risk of lung cancer (OR = 1.070, 95% CI = 1.027-1.114, P = 0.001, PFDR = 0.009), and CTSH levels decreased the risk of basal cell carcinoma (OR = 0.947, 95% CI = 0.919-0.975, P = 0.0002, P FDR = 0.002). In addition, there was no statistically significant effect of the 20 cancers on the nine cathepsins. Some unadjusted low P-value phenotypes are worth mentioning, including a positive correlation between cathepsin O (CTSO) and breast cancer (OR = 1.012, 95% CI = 1.001-1.025, P = 0.041), cathepsin S (CTSS) and pharyngeal cancer (OR = 1.017, 95% CI = 1.001-1.034, P = 0.043), and CTSS and endometrial cancer (OR = 1.055, 95% CI = 1.012-1.101, P = 0.012); and there was a negative correlation between cathepsin Z and ovarian cancer (CTSZ) (OR = 0.970, 95% CI = 0.949-0.991, P = 0.006), CTSS and prostate cancer (OR = 0.947, 95% CI = 0.902-0.944, P = 0.028), and cathepsin E (CTSE) and pancreatic cancer (OR = 0.963, 95% CI = 0.938-0.990, P = 0.006).

Conclusion

Our MR analyses showed a causal relationship between cathepsins and cancers and may help provide new insights for further mechanistic and clinical studies of cathepsin-mediated cancer.

Exogenous and endogenous antioxidants in osteoporosis risk: causal associations unveiled by Mendelian Randomization analysis

Background

Recent epidemiological studies and animal experiments have highlighted the significant role of oxidative stress in the development of osteoporosis (OP). The provision of antioxidants is widely considered a fundamental strategy to combat free radical-induced stress, inhibit oxidative damage, and potentially reverse the adverse effects of oxidative stress on bone health. However, there is no consensus in the scientific literature regarding the practical effectiveness of antioxidants in OP prevention and treatment. Some studies have not shown a clear connection between antioxidant supplementation and decreased OP risk. Therefore, it is essential to clarify the potential causal relationship between antioxidants and the development of OP.

Methods

The study utilized the inverse variance weighted (IVW) approach as the primary analytical method in the Mendelian Randomization (MR) framework to investigate the causal effects of five exogenous and six endogenous antioxidants on the risk of OP. To thoroughly assess potential pleiotropic effects and heterogeneity among the data analyzed, the MR-Egger intercept test was employed, and Cochran's Q statistic was calculated.

Results

In the evaluation of exogenous antioxidants, single-directional two-sample MR analyses did not reveal any statistically significant relationship between these agents and the risk of OP. Regarding endogenous antioxidants, bidirectional two-sample MR analyses were conducted, which generally indicated that most genetically regulated endogenous antioxidants had no significant association with the onset risk of OP. A significant causal relationship was found between OP and serum albumin levels (β: -0.0552, 95%CI: -0.0879 to -0.0225, p < 0.0011 after Bonferroni adjustment, power = 100%).

Conclusion

The research uncovers OP as a possible determinant contributing to a decrement in serum albumin levels, and further suggests a potentially intimate relationship between the downward trajectory of serum albumin concentrations and the advancement of the OP disease process.

Genetic Associations with Placental Proteins in Maternal Serum Identify Biomarkers for Hypertension in Pregnancy

Background

Preeclampsia is a complex syndrome that accounts for considerable maternal and perinatal morbidity and mortality. Despite its prevalence, no effective disease-modifying therapies are available. Maternal serum placenta-derived proteins have been in longstanding use as markers of risk for aneuploidy and placental dysfunction, but whether they have a causal contribution to preeclampsia is unknown.

Objective

We aimed to investigate the genetic regulation of serum placental proteins in early pregnancy and their potential causal links with preeclampsia and gestational hypertension.

Study design

This study used a nested case-control design with nulliparous women enrolled in the nuMoM2b study from eight clinical sites across the United States between 2010 and 2013. The first- and second-trimester serum samples were collected, and nine proteins were measured, including vascular endothelial growth factor (VEGF), placental growth factor, endoglin, soluble fms-like tyrosine kinase-1 (sFlt-1), a disintegrin and metalloproteinase domain-containing protein 12 (ADAM-12), pregnancy-associated plasma protein A, free beta-human chorionic gonadotropin, inhibin A, and alpha-fetoprotein. This study used genome-wide association studies to discern genetic influences on these protein levels, treating proteins as outcomes. Furthermore, Mendelian randomization was used to evaluate the causal effects of these proteins on preeclampsia and gestational hypertension, and their further causal relationship with long-term hypertension, treating proteins as exposures.

Results

A total of 2,352 participants were analyzed. We discovered significant associations between the pregnancy zone protein locus and concentrations of ADAM-12 (rs6487735, P= 3.03×10 -22 ), as well as between the vascular endothelial growth factor A locus and concentrations of both VEGF (rs6921438, P= 7.94×10 -30 ) and sFlt-1 (rs4349809, P= 2.89×10 -12 ). Our Mendelian randomization analyses suggested a potential causal association between first-trimester ADAM-12 levels and gestational hypertension (odds ratio=0.78, P= 8.6×10 -4 ). We also found evidence for a potential causal effect of preeclampsia (odds ratio=1.75, P =8.3×10 -3 ) and gestational hypertension (odds ratio=1.84, P =4.7×10 -3 ) during the index pregnancy on the onset of hypertension 2-7 years later. The additional mediation analysis indicated that the impact of ADAM-12 on postpartum hypertension could be explained in part by its indirect effect through gestational hypertension (mediated effect=-0.15, P= 0.03).

Conclusions

Our study discovered significant genetic associations with placental proteins ADAM-12, VEGF, and sFlt-1, offering insights into their regulation during pregnancy. Mendelian randomization analyses demonstrated evidence of potential causal relationships between the serum levels of placental proteins, particularly ADAM-12, and gestational hypertension, potentially informing future prevention and treatment investigations.

Cathepsin-mediated regulation of alpha-synuclein in Parkinson's disease: a Mendelian randomization study

Objective

The observational association between cathepsin and Parkinson's disease (PD) has been partially explored in previous research. However, the causal relationship remains unclear. In this study, our objective is to investigate the causal link between cathepsin and PD using Mendelian randomization (MR) analysis and elucidate the underlying mechanisms governing their interaction.

Methods

Utilizing bidirectional two-sample MR and multivariable MR, we systematically investigates the causal relationship between nine cathepsins and PD. The data pertaining to cathepsins were obtained from the Integrative Epidemiology Unit (IEU) Open GWAS Project, while data related to PD were sourced from versions R9 and R10 of the FinnGen database. The primary analytical method utilized was the inverse variance weighted (IVW), with MR analysis initially conducted using PD data from R9, complemented by a series of sensitivity analyses. Subsequently, replication analysis was performed on the R10 dataset, and meta-analysis were employed to merge the findings from both datasets. To explore potential mechanisms by which Cathepsins may impact PD, MR analyses were performed on significant Cathepsins with alpha-synuclein. MR analysis and colocalization analysis were conducted on expression quantitative trait loci (eQTL) data of gene related to alpha-synuclein with PD data.

Result

Forward MR analyses revealed more cathepsin B (CTSB) associated with less PD risk (OR = 0.898, 95%CI: 0.834-0.966, p = 0.004), while more cathepsin H (CTSH) (OR = 1.076, 95%CI: 1.007-1.149, p = 0.029) and more cathepsin S (CTSS) (OR = 1.076, 95%CI: 1.007-1.150, p = 0.030) associated with increasing PD risk. Meta-analyses validated these associations. Multivariate MR Results were consistent with those before adjustment. No significant results were observed in bidirectional MR analysis. In the investigation of the underlying mechanism, our findings demonstrate that CTSB significantly reduces the levels of alpha-synuclein (OR = 0.909, 95%CI: 0.841-0.983, p = 0.017). Concurrently, a genetically determined positive correlation between alpha-synuclein and PD is illuminated by both eQTL MR and colocalization analysis.

Conclusion

In conclusion, this MR study yields robust evidence suggesting an association between elevated levels of CTSB and reduced PD risk, mediated by the downregulation of alpha-synuclein levels. Conversely, higher levels of CTSH and CTSS are associated with an increased risk of PD. These findings offer novel insights into the pathophysiological mechanisms of PD and identify potential drug targets for disease prevention and treatment warranting further clinical investigations.

Roles of blood metabolites in mediating the relationship between vitiligo and autoimmune diseases: Evidence from a Mendelian randomization study

OBJECTIVE

This study employed Mendelian Randomization (MR) to investigate the causal relationship between genetic susceptibility to vitiligo and the risk of various autoimmune diseases, along with the mediating role of blood metabolites.

METHODS

We performed two-sample MR analyses using aggregated genome-wide association studies (GWAS) data on 486 blood metabolites, vitiligo, and nine autoimmune diseases to investigate blood metabolites' causal effects on the susceptibility of vitiligo and the associations of vitiligo with nine autoimmune comorbidities. We also applied multivariable MR to unravel metabolites by which vitiligo influences the pathogenesis of autoimmune diseases.

RESULTS

Our findings indicate that vitiligo amplified the risk of several autoimmune diseases, including rheumatoid arthritis (OR 1.17; 95 % CI 1.08-1.27), psoriasis (OR 1.10; 95 % CI 1.04-1.17), type 1 diabetes (OR 1.41; 95 % CI 1.23-1.63), pernicious anemia (OR 1.23; 95 % CI 1.12-1.36), autoimmune hypothyroidism (OR 1.19; 95 % CI 1.11-1.26), alopecia areata (OR 1.22; 95 % CI 1.10-1.35), and autoimmune Addison's disease (OR 1.22; 95 % CI 1.12-1.33). Additionally, our analysis identified correlations with vitiligo for 14 known (nine risk, five protective) and seven uncharacterized serum metabolites. After adjusting for genetically predicted levels of histidine and pyruvate, the associations between vitiligo and these diseases were attenuated.

CONCLUSIONS

We substantiated vitiligo's influence on susceptibility to seven autoimmune diseases and conducted a thorough investigation of serum metabolites correlated with vitiligo. Histidine and pyruvate are potential mediators of vitiligo associated with autoimmune diseases.By combining metabolomics with genomics, we provide new perspectives on the etiology of vitiligo and its immune comorbidities.

Predicted mechanistic impacts of human protein missense variants

Genome sequencing efforts have led to the discovery of tens of millions of protein missense variants found in the human population with the majority of these having no annotated role and some likely contributing to trait variation and disease. Sequence-based artificial intelligence approaches have become highly accurate at predicting variants that are detrimental to the function of proteins but they do not inform on mechanisms of disruption. Here we combined sequence and structure-based methods to perform proteome-wide prediction of deleterious variants with information on their impact on protein stability, protein-protein interactions and small-molecule binding pockets. AlphaFold2 structures were used to predict approximately 100,000 small-molecule binding pockets and stability changes for over 200 million variants. To inform on protein-protein interfaces we used AlphaFold2 to predict structures for nearly 500,000 protein complexes. We illustrate the value of mechanism-aware variant effect predictions to study the relation between protein stability and abundance and the structural properties of interfaces underlying trans protein quantitative trait loci (pQTLs). We characterised the distribution of mechanistic impacts of protein variants found in patients and experimentally studied example disease linked variants in FGFR1.

Role of Ubiquitin-conjugating enzyme E2 (UBE2) in two immune-mediated inflammatory skin diseases: a mendelian randomization analysis

Psoriasis vulgaris (PV) and Atopic dermatitis (AD) are the two major types of immune-mediated inflammatory skin disease (IMISD). Limited studies reported the association between Ubiquitin-conjugating enzyme E2 (UBE2) and IMISD. We employed a two-sample Mendelian randomization (MR) study to assess the causality between UBE2 and PV & AD. UBE2 association genome-wide association study (GWAS) data were collected. The inverse variance weighted (IVW) method was utilized as the principal method in our Mendelian randomization (MR) study, with additional using the MR-Egger, weighted median, simple mode, and weighted mode methods. The MR-Egger intercept test, Cochran's Q test, MR-Pleiotropy RESidual Sum and Outlier (MR-PRESSO) and leave-one-out analysis were conducted to identify heterogeneity and pleiotropy, colocalization analysis was also performed. The results showed that Ubiquitin-conjugating enzyme E2 variant 1 (UBE2V1) was causally associated with PV (OR = 0.909, 95% CI: 0.830-0.996, P = 0.040), Ubiquitin-conjugating enzyme E2 L3 (UBE2L3) was causally associated with AD (OR = 0.799, 95% CI: 0.709-0.990, P < 0.001). Both UBE2V1 and UBE2L3 may play protective roles in patients with PV or AD, respectively. No other significant result has been investigated. No heterogeneity or pleiotropy was observed. This study provided new evidence of the relationship between UBE2V1 and PV, UBE2L3 and AD. Our MR suggested that UBE2V1 plays an inhibitory role in PV progression, UBE2L3 plays an inhibitory role in AD. These could be novel and effective ways to treat PV and AD.

Investigation of lactotransferrin messenger RNA expression levels as an anti-type 2 asthma biomarker

BACKGROUND

Lactotransferrin (LTF) has an immunomodulatory function, and its expression levels are associated with asthma susceptibility.

OBJECTIVES

We sought to investigate LTF messenger RNA (mRNA) expression levels in human bronchial epithelial cells (BECs) as an anti-type 2 (T2) asthma biomarker.

METHODS

Association analyses between LTF mRNA expression levels in BECs and asthma-related phenotypes were performed in the Severe Asthma Research Program (SARP) cross-sectional (n = 155) and longitudinal (n = 156) cohorts using a generalized linear model. Correlation analyses of mRNA expression levels between LTF and all other genes were performed by Spearman correlation.

RESULTS

Low LTF mRNA expression levels were associated with asthma susceptibility and severity (P < .025), retrospective and prospective asthma exacerbations, and low lung function (P < 8.3 × 10-3). Low LTF mRNA expression levels were associated with high airway T2 inflammation biomarkers (sputum eosinophils and fractional exhaled nitric oxide; P < 8.3 × 10-3) but were not associated with blood eosinophils or total serum IgE. LTF mRNA expression levels were negatively correlated with expression levels of TH2 or asthma-associated genes (POSTN, NOS2, and MUC5AC) and eosinophil-related genes (IL1RL1, CCL26, and IKZF2) and positively correlated with expression levels of TH1 and inflammation genes (IL12A, MUC5B, and CC16) and TH17-driven cytokines or chemokines for neutrophils (CXCL1, CXCL6, and CSF3) (P < 3.5 × 10-6).

CONCLUSIONS

Low LTF mRNA expression levels in BECs are associated with asthma susceptibility, severity, and exacerbations through upregulation of airway T2 inflammation. LTF is a potential anti-T2 biomarker, and its expression levels may help determine the balance of eosinophilic and neutrophilic asthma.

Identifying potential drug targets for idiopathic pulmonary fibrosis: a mendelian randomization study based on the druggable genes

BACKGROUND

Idiopathic pulmonary fibrosis (IPF) is a chronic fibrotic interstitial lung disease characterized by progressive dyspnea and decreased lung function, yet its exact etiology remains unclear. It is of great significance to discover new drug targets for IPF.

METHODS

We obtained the cis-expression quantitative trait locus (cis-eQTL) of druggable genes from eQTLGen Consortium as exposure and the genome wide association study (GWAS) of IPF from the International IPF Genetics Consortium as outcomes to simulate the effects of drugs on IPF by employing mendelian randomization analysis. Then colocalization analysis was performed to calculate the probability of both cis-eQTL of druggable genes and IPF sharing a causal variant. For further validation, we conducted protein quantitative trait locus (pQTL) analysis to reaffirm our findings.

RESULTS

The expression of 45 druggable genes was significantly associated with IPF susceptibility at FDR < 0.05. The expression of 23 and 15 druggable genes was significantly associated with decreased forced vital capacity (FVC) and diffusing capacity of the lungs for carbon monoxide (DLco) in IPF patients, respectively. IPF susceptibility and two significant genes (IL-7 and ABCB2) were likely to share a causal variant. The results of the pQTL analysis demonstrated that high levels of IL-7 in plasma are associated with a reduced risk of IPF (OR = 0.67, 95%CI: 0.47-0.97).

CONCLUSION

IL-7 stands out as the most promising potential drug target to mitigate the risk of IPF. Our study not only sheds light on potential drug targets but also provides a direction for future drug development in IPF.

Proteome-wide Mendelian randomization identifies potential therapeutic targets for nonalcoholic fatty liver diseases

Nonalcoholic fatty liver disease (NAFLD) is the predominant cause of liver pathology. Current evidence highlights plasma proteins as potential therapeutic targets. However, their mechanistic roles in NAFLD remain unclear. This study investigated the involvement of specific plasma proteins and intermediate risk factors in NAFLD progression. Two-sample Mendelian randomization (MR) analysis was conducted to examine the association between plasma proteins and NAFLD. Colocalization analysis determined the shared causal variants between the identified proteins and NAFLD. The MR analysis was applied separately to proteins, risk factors, and NAFLD. Mediator shares were computed by detecting the correlations among these elements. Phenome-wide association studies (phewas) were utilized to assess the safety implications of targeting these proteins. Among 1,834 cis-protein quantitative trait loci (cis-pQTLs), after-FDR correction revealed correlations between the plasma levels of four gene-predicted proteins (CSPG3, CILP2, Apo-E, and GCKR) and NAFLD. Colocalization analysis indicated shared causal variants for CSPG3 and GCKR in NAFLD (posterior probability > 0.8). Out of the 22 risk factors screened for MR analysis, only 8 showed associations with NAFLD (p ≤ 0.05), while 4 linked to CSPG3 and GCKR. The mediator shares for these associations were calculated separately. Additionally, reverse MR analysis was performed on the pQTLs, risk factors, and NAFLD, which exhibited a causal relationship with forward MR analysis. Finally, phewas summarized the potential side effects of associated-targeting proteins, including CSPG3 and GCKR. Our research emphasized the potential therapeutic targets for NAFLD and provided modifiable risk factors for preventing NAFLD.

Exploration of potential novel drug targets for diabetic retinopathy by plasma proteome screening

The aim of this study is to identify novel potential drug targets for diabetic retinopathy (DR). A bidirectional two-sample Mendelian randomization (MR) analysis was performed using protein quantitative trait loci (pQTL) of 734 plasma proteins as the exposures and clinically diagnosed DR as the outcome. Genetic instruments for 734 plasma proteins were obtained from recently published genome-wide association studies (GWAS), and external plasma proteome data was retrieved from the Icelandic Decoding Genetics Study and UK Biobank Pharma Proteomics Project. Summary-level data of GWAS for DR were obtained from the Finngen Consortium, comprising 14,584 cases and 202,082 population controls. Steiger filtering, Bayesian co-localization, and phenotype scanning were used to further verify the causal relationships calculated by MR. Three significant (p < 6.81 × 10-5) plasma protein-DR pairs were identified during the primary MR analysis, including CFH (OR = 0.8; 95% CI 0.75-0.86; p = 1.29 × 10-9), B3GNT8 (OR = 1.09; 95% CI 1.05-1.12; p = 5.9 × 10-6) and CFHR4 (OR = 1.11; 95% CI 1.06-1.16; p = 1.95 × 10-6). None of the three proteins showed reverse causation. According to Bayesian colocalization analysis, CFH (coloc.abf-PPH4 = 0.534) and B3GNT8 (coloc.abf-PPH4 = 0.638) in plasma shared the same variant with DR. All three identified proteins were validated in external replication cohorts. Our research shows a cause-and-effect connection between genetically determined levels of CFH, B3GNT8 and CFHR4 plasma proteins and DR. The discovery implies that these proteins hold potential as drug target in the process of developing drugs to treat DR.

Causal roles and clinical utility of cardiovascular proteins in colorectal cancer risk: a multi-modal study integrating mendelian randomization, expression profiling, and survival analysis

PURPOSE

This comprehensive investigation delved into the intricate causal interplay existing between cardiovascular-related plasma proteins and the susceptibility to colorectal cancer, leveraging the robust framework of Mendelian randomization, and employed expression profiling and survival analysis to unravel the latent clinical worth embedded within pertinent gene expressions.

METHODS

Protein quantitative trait loci (pQTLs) of 85 cardiovascular proteins were employed as instrumental variables to investigate the causal relationship between proteins and CRC risk using a Mendelian randomization approach. Causal inferences were graded as strong, intermediate or weak based on statistical checks. Drug-target MR examined VEGF receptors for their potential as therapeutic targets for colorectal cancer. Differential expression analysis, diagnostic ROC curves, and survival analyses were performed for identified proteins using RNA-seq data from The Cancer Genome Atlas (TCGA) colorectal cancer cohort.

RESULTS

Using cis-pQTLs, LOX-1, VEGF-A and OPG were associated with increased CRC risk (strong evidence), while PTX3, TNF-R2 and MMP-7 were protective (strong evidence). Pan-pQTL analysis found MMP-10 increased risk (intermediate evidence) and ADM increased risk (weak evidence). Drug-target MR found VEGF R1 may be promising therapeutic targets. Differential expression analysis revealed seven genes encoding the identified proteins were dysregulated in tumors. ROC analysis showed five gene expression had high diagnostic accuracy. KM analysis showed four genes had prognostic value.

CONCLUSIONS

This large-scale MR study implicates several cardiovascular proteins in CRC susceptibility and progression. Findings highlight roles for VEGF signaling and extracellular matrix regulation. Results nominate specific proteins as potential diagnostic biomarkers or therapeutic targets warranting further investigation.

Identification of plasma protein markers of allergic disease risk: a mendelian randomization approach to proteomic analysis

BACKGROUND

While numerous allergy-related biomarkers and targeted treatment strategies have been developed and employed, there are still signifcant limitations and challenges in the early diagnosis and targeted treatment for allegic diseases. Our study aims to identify circulating proteins causally associated with allergic disease-related traits through Mendelian randomization (MR)-based analytical framework.

METHODS

Large-scale cis-MR was employed to estimate the effects of thousands of plasma proteins on five main allergic diseases. Additional analyses including MR Steiger analyzing and Bayesian colocalisation, were performed to test the robustness of the associations; These findings were further validated utilizing meta-analytical methods in the replication analysis. Both proteome- and transcriptome-wide association studies approach was applied, and then, a protein-protein interaction was conducted to examine the interplay between the identified proteins and the targets of existing medications.

RESULTS

Eleven plasma proteins were identified with links to atopic asthma (AA), atopic dermatitis (AD), and allergic rhinitis (AR). Subsequently, these proteins were classified into four distinct target groups, with a focus on tier 1 and 2 targets due to their higher potential to become drug targets. MR analysis and extra validation revealed STAT6 and TNFRSF6B to be Tier 1 and IL1RL2 and IL6R to be Tier 2 proteins with the potential for AA treatment. Two Tier 1 proteins, CRAT and TNFRSF6B, and five Tier 2 proteins, ERBB3, IL6R, MMP12, ICAM1, and IL1RL2, were linked to AD, and three Tier 2 proteins, MANF, STAT6, and TNFSF8, to AR.

CONCLUSION

Eleven Tier 1 and 2 protein targets that are promising drug target candidates were identified for AA, AD, and AR, which influence the development of allergic diseases and expose new diagnostic and therapeutic targets.