- 1Qi Huang Chinese Medicine Academy, Jiangxi University of Chinese Medicine, Nanchang, Jiangxi, China
- 2School of Clinical Medicine, Jiangxi University of Chinese Medicine, Nanchang, Jiangxi, China
- 3Pulmonary Disease Department, Affiliated Hospital of Jiangxi University of Traditional Chinese Medicine, Nanchang, Jiangxi, China
Background: Abdominal aortic aneurysm (AAA) and Sjögren's syndrome (SS) frequently coexist, suggesting shared pathogenesis, but their genetic-immunological links are unclear.
Objective: Investigate the shared genetic architecture and immune pathways between SS and AAA.
Methods: Using European GWAS summary statistics (SS: 585 cases/1,546 controls; AAA: 4,083 cases/420,324 controls), we applied complementary genomics approaches: LDSC (genetic correlation), S-LDSC (tissue heritability), PLACO/FUMA (pleiotropic loci), MAGMA/Metascape (pathways), SMR (druggable targets), and HyPrColoc (immune cells).
Key results: A significant positive genetic correlation exists (rg = 0.32, PFDR = 0.021). We identified 8 shared risk SNPs and 6 pleiotropic genes (e.g., HLA-B, HLA-DQB2, LSM2) within loci 6p22.2-21.32. Pathway analyses revealed significant enrichment for MHC class II antigen presentation (P = 3.1 × 10−12) and U6 snRNA binding/spliceosome (P = 2.8 × 10−7). Tissue-specific heritability enrichment occurred in artery/aorta, kidney, and secretory tissues (all pS-LDSC < 9.3 × 10−4). Immune co-localization implicated myeloid dendritic cells expressing HLA-DR (rs9272318) in convergent dysregulation. HLA-B emerged as a prioritised druggable target (pSMR = 1.65 × 10−8).
Conclusion: This study establishes a shared genetic and immunological basis for SS and AAA, driven primarily by dysregulated HLA-mediated antigen presentation (HLA-B/HLA-DR), spliceosome dysfunction, and NK cell impairment. These findings provide mechanistic insights for early AAA detection in SS patients and support developing immunotherapies targeting HLA pathways.
1 Introduction
Abdominal aortic aneurysm (AAA) is a focal dilatation of the arterial wall due to weakness caused by atherosclerosis (AS) (1). It is typified by immune cell infiltration, enhanced proteolytic activity, and continuing extracellular matrix component degradation, leading to the aorta wall to dilate (2). AAA is usually asymptomatic and challenging to identify quickly, but once it ruptures, the surgical repair mortality rate can reach 31%–70% (1). Sjogren's Syndrome (SS) is an autoimmune disease (AID) marked by lymphocytic infiltration of salivary and lacrimal glands, leading to glandular dysfunction, where CD4-positive helper T cells (TH) and their cytokines are crucial in the pathogenesis (3).
There are several biological effect mechanisms and possible correlations between AAA and autoimmune-related substances based on evidence from modern medical research.There are significant variations in the quantity of five immune cell types and approximately 40 AID-related bioactive peptides in circulation between AAA and non-AAA patients (4). In AAA patients, IL-1 levels surpass 90 pg/ml (5). Other bioactive peptides, such as myeloperoxidase (MPO) and platelet factor 4 (PF4), additionally display an increasing tendency, surpassing 250 ng/ml and 300 ng/ml, respectively (6). These findings imply that the occurrence of AAA is pathologically impacted by autoimmune diseases (7). Disrupting the connection between Th17 cells and IL-17 can effectively delay the course of AAA, according to studies conducted using mouse models generated by Ang II (8). Additionally, the occurrence of AAA is significantly influenced by neutrophil extracellular traps (NETs), which are secreted by neutrophils. Myeloperoxidase and citrullinated histone H3 (H3Cit), which are indicators of NETs, are markedly increased in AAA patients and have a positive correlation with both MPO and elastase (P < 0.001) (9). According to epidemiological research, there is an association between certain AID and AAA.For example, 309 cases (2.0%) of 14,816 patients who underwent elective AAA repair were diagnosed with Rheumatoid Arthritis (RA) during a 13-year cohort study (10). In addition, a 20-year retrospective study discovered that 38.9% of Systemic Lupus Erythematosus (SLE) patients on long-term prednisone treatment may have accelerated AS progression, which might lead significantly to the appearance and death of AAA (11). However, the co-morbidity between SS and AAA has not been studied; hence the purpose of this study was to look at the genetic relationship between SS and AAA.
The concurrent occurrence of AAA and AID consistently encourages research into their shared genetic foundation. Focusing on the proposed, the aim of this study was to investigate the genetic correlation between SS in AID and AAA. However, the impact of AID on immune-related bioactive peptides is the main topic of current study. Based on Lu et al., the induction of autoimmune antigens may be linked to the incidence of AAA (12, 13). The pathological progress in AAA is impacted by this process, which causes the intima to vanish and collagen and elastic fibers to rupture, degrade, and be damaged (14, 15). The effect of inherited telomere length on the correlation between AAA and AID, however, was only investigated by Philip C. Haycock et al. (16). Zhang et al. did not explore the genetic interaction pathways between AAA and SLE, but they did examine the immunological dysregulation based on by genetic parameters, which may be a possible explanation for the co-occurrence of both diseases (17). Relevant routes in AID that are meaningful for the presence, grading, progression, or rupture of AAA (5). It is clear that abnormalities in these processes cause AID and AAA to develop. The urgent need to precisely identify the shared genetic risk loci between SS and AAA. Further clarifying their possible immune-mediated connections is highlighted by considering that there are still substantial knowledge gaps regarding the pleiotropic mechanisms and bidirectional causal relationships between these two diseases. Thus, the purpose of this study was to look into the risk genetic structural loci that SS and AAA share. It is necessary to be aware that the statistical validity of clinical or epidemiological research may face challenges.
The following six series of studies were conducted in order to systematically evaluate the genetic overlap, shared susceptibility genes, and potential effector connections between SS and AAA: (1) assessing if there is a genetic association between the two diseases using linkage disequilibrium (LD) score regression (LDSC) based on GWAS summary data (18). (2) The genetic enrichment of SNPs in particular tissues and organs was evaluated using a stratified LDSC (S-LDSC) regression approach. (3) Pleiotropic genetic loci were identified at the SNP level using “PLACO,” risk loci were identified using FUMA, and shared pleiotropic loci were identified using eQTL mapping (19, 20). (4) FUMA, PLACO, MAGMA, and summary data-based MR (SMR) are used to clarify the underlying mechanisms (21). (5) Using immune co-localization to probe immune-mediated processes (22). The various forms of pleiotropy (horizontal or vertical) were then partially explained, and pairwise causal connections were evaluated using the randomization analysis. The flowchart of our study is shown in Figure 1.

Figure 1. Study workflow: integrated pipeline for shared genetics and immune mechanisms between SS and AAA. (Key steps: LDSC was used to assess genetic correlation; S-LDSC and MAGMA identified cell/tissue and gene-level enrichment; PLACO detected pleiotropic loci, with FUMA annotation; SMR integrated with GTEx prioritized drug targets; HyPrColoc was applied for multi-trait immune colocalization across 731 immune cell traits, identifying 2,193 region–cell colocalization events and highlighting regulatory pathways linking SNPs to immune cell function. AAA, abdominal aortic aneurysm; SS, Sjögren's syndrome; LDSC, linkage disequilibrium score regression; S-LDSC, stratified LDSC; MHC, major histocompatibility complex; MAF, minor allele frequency; PLACO, pleiotropy analysis under composite null hypothesis; FUMA, functional mapping and annotation; eQTL, expression quantitative trait locus; GTEx, Genotype-Tissue Expression; SMR, summary-data-based Mendelian randomization; GO, Gene Ontology; GWAS, genome-wide association study; HyPrColoc, Hypothesis Prioritisation in multi-trait Colocalization).
2 Method
2.1 Dataset
GWAS summary statistics for Sjögren's syndrome (SS) were obtained from open-access source: The FinnGen release R11 dataset (4,083 cases and 420,324 controls, all of European ancestry). Abdominal aortic aneurysm (AAA) summary statistics were derived from the GWAS Catalog (GCST012796; 585 cases, 1,546 controls) (23). Genotyping was performed using high-throughput platforms (whole-genome sequencing or SNP microarrays), and standard quality control procedures were applied to remove low-quality samples and variants.
Beagle 4.1 software and the Finnish population-specific SISu v3 reference panel were used for AAA genotype imputation, and the TOPMed reference panel for SS. Eagle 2.3.5 was used for phasing. For eQTL analyses, we used (1) eQTLGen Consortium whole blood cis-eQTL summary statistics (n ≈ 31,684) for immune-related expression profiles, and (2) GTEx v8 multi-tissue cis-eQTL data (54 tissues) for tissue-specific regulation. All SNP positions were aligned to GRCh37(hg 19). To avoid potential sample overlap between the SS and AAA cohorts, we verified cohort independence by (i) confirming distinct recruitment sources (FinnGen vs. independent GWAS Catalog datasets) and (ii) cross-referencing with eQTL sample lists. No individual-level overlap was identified between the GWAS and eQTL datasets (Table 1).
Stringent quality control (QC) procedures were applied to ensure data reliability. Firstly, variants with minor allele frequency (MAF) < 0.01, sample call rate <95%, or SNP call rate <99% were excluded. Post-imputation QC required an INFO score >0.9.
For SMR analyses, cis-eQTLs within ±1 Mb of transcription start sites were tested, with significance defined as p_SMR < 0.05 and p_HEIDI > 0.05 to exclude linkage. MAGMA mapped SNPs to genes ±50 kb using the 1000 Genomes Phase 3 EUR reference panel (24), and gene-based tests were Bonferroni-corrected. PLACO was implemented under the composite null hypothesis using HapMap3 SNPs after LD pruning (r2 < 0.1), with significance threshold p < 5 × 10−8. No individual-level sample overlap existed between the AAA and SS datasets, ensuring unbiased cross-trait analyses.
2.2 Quality control standards
We employed the LD Score regression (LDSC) to estimate the genetic correlation between AAA and SS. LD scores were precomputed from HapMap3 SNPs using the 1000 Genomes Project Phase 3 European reference panel (24), matched to the ancestry of the GWAS summary statistics (25). By using the block jackknife method in LDSC, the standard error (SE) was calculated and then adjusted for attenuation bias. Summary statistics were taken from independent sources (AAA: FinnGen; SS: GWAS Catalog), and we confirmed no overlapping individuals across datasets. It is crucial to mention in our study that SS and AAA studies did not exactly overlap in population, which improves the validity of the conclusions. When multiple correlations were tested, Benjamini–Hochberg FDR control was applied with fdr < 0.05.
2.3 Organ-level association analysis
We assessed AAA–SS heritability enrichment in specific tissues using Stratified LD Score regression (S-LDSC) (26) with baseline LD model v2.2 and 54 GTEx v8 tissue annotations (27). Enrichment Z-scores were computed using default settings with 1000 Genomes EUR LD reference. A Bonferroni-corrected significance threshold of pbon < 0.05/54 was used.
2.4 Gene-level exploration and analysis
We mapped overlapping genes according to the lead SNP in each locus in order to investigate the shared mechanisms across AAA and SS. Multi-marker effects were then assessed with Multi-marker Analysis of GenoMic Annotation (MAGMA) (28), mapping SNPs to genes within ±50 kb and using LD from the 1000 Genomes Phase 3 European reference panel. Gene-level significance was controlled by Bonferroni correction across 17,008 genes (pbon < 0.05/17,008 = 2.94 × 10−6). For gene-set enrichment, we used MSigDB v7.5.1 curated sets (c2.all) (28, 29, 30). Multiple testing was controlled by Bonferroni across 10,678 tested sets. Pathway enrichment of mapped genes was performed with Metascape (metascape.org) to aid functional interpretation (31). Genome-wide pleiotropic loci identified by PLACO underwent tissue-specific enrichment using GTEx v8 annotations across 54 tissues. We also conducted tissue-specific expression profiling for pleiotropic genes by computing their average transformed expression in each GTEx tissue and classifying genes as up- or down regulated by the sign of the corresponding t-statistics.
2.5 Pleiotropy analysis and locus validation
We applied Pleiotropy under the Composite null hypothesis (PLACO), a pleiotropy analysis framework under the composite null hypothesis (32), to detect SNPs jointly associated with AAA and SS using a genome-wide significance threshold of p < 5 × 10−8. Variants were LD-pruned at r2 < 0.1 based on the 1000 Genomes Project European reference panel. Genome-wide significant SNPs were mapped to genomic risk loci using Functional Mapping and Annotation (FUMA) (33), and shared loci were further evaluated by Bayesian colocalization analysis (34) to assess whether the same causal variant was likely driving both traits. Compared with traditional univariate genome-wide association studies, PLACO increases power to identify shared genetic variants by jointly modeling the null hypotheses for each phenotype, enabling detection of pleiotropic effects even when single-trait associations do not reach genome-wide significance. This approach effectively controls false positives while maintaining sensitivity to loci with modest, concordant effects across traits (32).
2.6 Potential exploration of drug targets in European populations
We applied summary-data–based Mendelian randomization (SMR) to integrate summary statistics from genome-wide association studies (GWAS) with expression quantitative trait locus (eQTL) data, in order to identify genes whose expression levels may causally or pleiotropically influence AAA and SS (35). GWAS data were combined with eQTLGen whole-blood and GTEx v8 cis-eQTL datasets using the parameters p_eqtl_smr 5 × 10−8 and p_HEIDI 0.05. Only the top cis-eQTL per gene was retained, and associations with HEIDI p < 0.05 were excluded. The HEIDI (Heterogeneity in Dependent Instruments) test distinguishes linkage from pleiotropy by assessing heterogeneity among SNP–gene expression associations within a locus (36). This approach increases confidence that observed overlap between GWAS and eQTL signals is due to a shared causal variant, rather than distinct but correlated variants (27). Candidate genes passing these filters were considered supported by multiple lines of evidence, providing insights into regulatory mechanisms underlying the traits.
2.7 Immune co-localization analysis
We enhanced the co-localization analysis to find comparable immune-regulatory pathways and processes between SS and AAA by concurrently integrating a broad range of immunological GWAS data (AAA, SS, and 731 immune cell phenotypes) using “HyPrColoc v1.0” (Hypothesis Prioritisation in multi-trait Colocalization) with default priors (prior_prob_coloc = 0.02, prior_prob_single = 1 × 10−4) (37). Immune cell GWAS data were sourced from the “GWAS Catalog”. Co-localization significance: “regional probability >0.7” (36, 37). This approach enables hypothesis-free exploration of shared immune mechanisms, although statistical power for rare immune cell subsets (e.g., myeloid dendritic cells) was limited by modest cohort sizes.
3 Result
3.1 Shared genetic structure between SS and AAA
First, we evaluated genetic correlations between AAA and selected AIDs. Table 2 presents LDSC results. Genetic correlation coefficients (rg) were FDR-corrected. Significant associations (rg > 0, pFDR < 0.05) were reported, accounting for genetic covariance between traits. Results revealed significant genetic associations (p < 0.05) between AAA and three AIDs: SS, RA, and primary biliary cholangitis (PBC).
3.2 Results of organ-level association analysis
We employed S-LDSC to assess tissue-specific enrichment of SNP heritability for SS and AAA. First, we obtained GTEx datasets containing gene sets and baseline models for 54 human tissues. Using regression coefficient Z-scores and P-values, we applied S-LDSC to evaluate tissue-specific genetic enrichment significance. Tissue-specific analysis revealed significant enrichment in multiple tissues including artery, aorta, bladder, kidney, prostate, testis, uterus, and vagina (Figure 2).
3.3 Results of gene-level enrichment analysis
Using FUMA for MAGMA gene enrichment analysis, we identified 388 significantly enriched genes associated with SS and AAA. Six genes showed significant enrichment (pfdr><0.05) (Figure 3). These genes participate in key pathways including: allograft rejection; NK cell-mediated cytotoxicity protection; MHC class I/II peptide assembly; antigen processing/presentation; and U6 snRNA binding (Figures 4, 5). Gene set analysis further implicated these genes in pro-inflammatory cytokine secretion, macrophage activation, lysosomal degradation regulation, and snRNP assembly via protein binding. Significant enrichment occurred in secretory tissues including thyroid, stomach, pituitary, pancreas, kidney medulla, and kidney cortex, suggesting SS-mediated disruption of secretory functions. However, tissue specificity did not survive Bonferroni correction for 54 tests (pb > 0.05/54). These results demonstrate shared SS-AAA genetic mechanisms.

Figure 3. Results of gene enrichment [The x-axis represents the chromosomal position, and the y axis represents the uncorrected −log10 (P) from two-sided z-tests for SNP associations with SS. The red line represents a significance threshold of P = 5 × 10−6].

Figure 4. Results of MAGMA gene-set enrichment analysis (the red dotted line indicates the multiple testing-corrected significance threshold of 0.05. Significantly enriched pathways (highlighted in red) include U6 snRNA binding and MHC class II complex, suggesting shared RNA processing and immune activation mechanisms).

Figure 5. Results of MAGMA tissue-specific enrichment analysis (the red dotted line indicates the multiple testing-corrected significance threshold of 0.05).
3.4 Results of SNP-level pleiotropy evaluation and risk loci analysis
Based on LDSC-identified shared genetics, we applied pleiotropy analysis (PLACO) to detect pleiotropic loci. Using PLACO results, we employed FUMA to identify AAA-associated SS risk loci (Table 3). Notably, shared pleiotropic regions (e.g., 6p22.2, 6p21.33, 6p21.32) were identified for both SS and AAA (Supplementary Table S1).
3.5 Drug targets in European populations
First, we identified three potential drug targets using SMR (p_SMR < 0.05, p_HEIDI > 0.05, p_b < 0.05/15,638 = 3.2 × 10−6) (Figure 6). While HLA-B mediates organ transplantation compatibility (38), HLA-DQB2 is implicated in metastatic lung adenocarcinoma (39), and LSM2 in hepatocellular carcinoma (HCC) (40), these loci lack established roles in AAA or SS pathogenesis. We integrated PLACO, FUMA, MAGMA, and SMR results to identify pleiotropic genes associated with multiple traits (Table 4). These genes showed significant tissue-specific genetic signals. eQTL and SMR analyses confirmed their pleiotropic effects across traits and provided chromosomal localization.

Figure 6. European population drug target (signals represent gene hits across trait pairs. eQTL, expression quantitative trait locus; SMR, summary-data-based Mendelian randomization).
3.6 Results of immuno-colocalization
Immune-mediated AAA pathogenesis is suggested by shared dysfunction in secretory tissues (thyroid, stomach, pituitary, pancreas, renal medulla/cortex). We applied HyPrColoc to identify immune cells showing multi-trait co-localization. This analysis targeted critical immune cells (Supplementary Table S2). One immune cell type showed significant contribution to SS, supported by a multi-effector locus with shared causal variants (p > 0.6). We confirmed conventional dendritic cells (cDCs) significantly influence specific immune subsets. Notably, myeloid dendritic cells showed AAA-associated HLA-DR expression driven by rs9272318, suggesting mechanistic links between secretory dysfunction and immune processes (Figure 7).

Figure 7. Immune colocalization analysis of AAA and SS risk loci with immune cell traits using hyPrColoc. (HyPrColoc (Hypothesis Prioritisation in multi-trait Colocalization) was applied to three pleiotropic loci identified by PLACO for AAA and SS, testing colocalization with 731 immune cell phenotypes. The top panel shows the SNP-level posterior probability of being the shared causal variant for each candidate SNP. The bottom panel compares regional posterior probability (x-axis) and SNP-level posterior probability (y-axis) across all locus–cell pairs. The vertical reference line at a regional posterior probability of 0.7 indicates the threshold for significant colocalization. Points to the right of this line represent AAA–SS–immune trait combinations with strong evidence of a shared causal variant).
4 Discussion
SS is a systemic autoimmune disorder affecting multiple organ systems. AAA involvement in SS may lead to severe complications. We applied complementary analytical approaches with distinct modeling frameworks to comprehensively investigate pleiotropic relationships.
LDSC identified significant genetic associations between AIDs and AAA, specifically SS, RA, and PBC. Although RA and PBC comorbidity with AAA is established (10, 41, 42), SS-AAA evidence was previously confined to a single case report (43). We provide the first population-level genetic evidence of shared SS-AAA genetic architecture. This shared risk likely arises from immune dysregulation, demonstrated by T cell, B cell, and antigen-presenting cell infiltration in AAA lesions (12)—key features of SS pathogenesis.
4.1 Multi-omics approach revealed key shared biological mechanisms
Antigen Presentation & MHC Pathways: Gene enrichment analysis revealed MHC class II pathway involvement (Figure 4). Immunolocalization identified HLA-DR+ myeloid dendritic cells as a key shared signature (regional probability >0.6; Figure 7). This convergence implies aberrant antigen presentation activates autoreactive T cells in both diseases (44, 45), evidenced by HLA-DR+ dendritic cell enrichment in AAA lesions (12) and SS salivary glands (46).
Spliceosome Dysregulation: U6 snRNA binding pathway enrichment (Figure 4) indicates spliceosome dysfunction as a novel link. This connects autoimmunity to ECM degradation—a feature of SS (glandular disruption) (47) and AAA (aortic rupture via collagen/elastin defects) (2, 48).
NK Cell Dysfunction: Enrichment of NK cell-mediated cytotoxicity (Figure 4) and lead SNP rs2076030 in the KIR locus (6p22.2; Table 3) implicate NK dysfunction as a shared mechanism. In SS, impaired NK surveillance permits autoantibody production (46), whereas in AAA, defective cytotoxicity fails to control macrophage-driven inflammation (12). This KIR/HLA imbalance constitutes a regional pathogenic risk.
Shared genomic loci (6p21.32-33, 6p22.2; Table 3, Supplementary Table S1) overlap regions associated with autoimmune diseases, including rheumatoid arthritis (RA; rs2076030 at 6p22.2) (49). These loci contain genes regulating immune activation, antigen presentation, and inflammation (50, 51). Tissue-specific heritability enrichment (S-LDSC; Figure 2) occurred in immune-relevant (artery, aorta) and secretory tissues (bladder, kidney, prostate, testis, uterus, vagina). This parallels SS-related secretory gland and renal/bladder dysfunction (52), suggesting systemic inflammation (53, 54, 55) or immune-metabolic injury (56, 57) links SS secretory defects to AAA vascular vulnerability. Artery/aorta enrichment indicates direct vascular immune involvement, potentially through hemodynamic stress- or IFN-γ-induced vascular smooth muscle cell (VSMC) inflammation (58).
4.2 Clinical translation and novel therapeutic insights
Drug target analysis identified HLA-B, HLA-DQB2, and LSM2 as druggable targets (Table 4, Figure 6). HLA-B emerges as a particularly promising candidate. Its role in presenting antigenic peptides to activate self-reactive T cells (59, 60) is central to SS pathogenesis and is now implicated in AAA susceptibility. However, the direct therapeutic targeting of specific HLA molecules remains challenging due to their extreme polymorphism and fundamental role in adaptive immunity. Consequently, no clinical trials are currently running that directly target HLA-B or HLA-DR alleles for SS or AAA. However, broader immunomodulatory strategies that indirectly impact HLA-mediated antigen presentation, such as T-cell co-stimulation blockade (e.g., with Abatacept) or Janus kinase (JAK) inhibition, are being explored in related autoimmune contexts and represent a promising indirect avenue for translation (61, 62). The colocalization of rs9272318 with HLA-DR on myeloid DCs further implicates antigen presentation as a convergent and therapeutically targetable node. These findings provide a strong genetic rationale for exploring immunomodulatory approaches (63). This could involve repurposing existing therapies or developing novel biologics targeting these shared pathways, enabling early AAA intervention and concurrent management of SS-AAA comorbidity (3).
4.3 Molecular mechanistic conclusion
In summary, our integrated genetic analyses delineate a convergent pathogenic model whereby shared risk loci (e.g., HLA-B, HLA-DQB2) primarily drive the comorbidity of SS and AAA through a core mechanism of dysregulated antigen presentation. This process activates autoreactive T cells and initiates a chronic inflammatory cascade (64). While fundamentally immune-mediated, this inflammation secondarily triggers the activation of vascular endothelial and smooth muscle cells within the aortic wall, culminating in the extracellular matrix degradation and remodeling that defines AAA progression (65). Thus, the interplay between innate immune dysregulation and subsequent vascular cell activation forms the critical bridge connecting systemic autoimmunity to local aortic wall instability (66, 67).
5 Limitation
While our study identified shared genetic architecture between SS and AAA in European populations, the generalizability of these findings to diverse populations requires further validation. Future studies should prioritize multi-ethnic cohorts to evaluate the relevance of identified risk loci (including HLA-B and LSM2) and immune pathways across populations. Population stratification analyses were limited by our use of summary-level data. GWAS data for specific immune cell subsets (e.g., HLA-DR + myeloid dendritic cells) had limited sample sizes, potentially reducing the robustness of colocalization signals such as rs9272318-HLA-DR. Leveraging larger immunogenomic resources (e.g., UK Biobank) and performing single-cell eQTL mapping in diseased tissues will be essential. Concurrently, single-cell multi-omics expansion strategies and AI-enhanced phenotypic analysis (developing deep phenotypic inference models: predicting the proportions of 78 immune cell subpopulations using conventional blood test parameters; establishing a transfer learning framework: utilizing mouse model data to enhance human GWAS signal detection) can improve statistical power and result reliability. Three-dimensional genome-guided localization (integrating Hi-C/ChIA-PET data to construct immune cell-specific chromatin interaction networks) and cross-omics Bayesian integration models can further achieve closed-loop validation from genetic signals to functional mechanisms. Although we propose mechanistic links (HLA-DR+ dendritic cells, U6 snRNA processing, KIR-HLA interactions) supported by genetic colocalization and pathway enrichment, experimental validation is required: 1) validate spatially resolved HLA-DR expression and U6 snRNA-associated splicing factors in tissues from patients with SS-AAA comorbidity; 2) assess NK cell function in carriers of identified KIR variants with both SS and AAA.
6 Conclusion
Our study reveals significant genetic and immunological crosstalk between SS and AAA. We identified shared risk loci in chromosomal regions 6p22.2, 6p21.33, and 6p21.32, implicating three convergent mechanisms: 1) dysregulated antigen presentation by HLA-DR+ myeloid dendritic cells; 2) spliceosome dysfunction involving U6 snRNA-binding factors; 3) impaired natural killer (NK) cell function associated with KIR-HLA interactions. Druggable targets, particularly HLA-B, were identified, supporting development of novel immunotherapies concurrently targeting both disorders. These findings provide critical insights for: 1) early AAA detection in SS cohorts; 2) concurrent management of SS-AAA comorbidity.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.
Ethics statement
All data analyzed in this study are publicly available, de-identified GWAS summary statistics. The collection of the original data was conducted in compliance with ethical guidelines and with informed consent obtained in the respective studies (FinnGen and GWAS Catalog). Therefore, this study itself does not require new ethical approval.
Author contributions
YZ: Methodology, Writing – original draft. HX: Data curation, Software, Writing – original draft. QL: Data curation, Writing – review & editing. ZD: Methodology, Writing – review & editing. JH: Supervision, Writing – review & editing. LW: Funding acquisition, Supervision, Writing – review & editing.
Funding
The author(s) declare that financial support was received for the research and/or publication of this article. This work was supported by the Key Laboratory of TCM Pulmonary Science of Jiangxi Province (No. 2024SSY06321), the Key Research and Development Project of Jiangxi Province, the field of Social Development, the development and application of the Yin-Yang Attribute Breath Recognition Instrument for chronic obstructive pulmonary Disease (20232BBG70021), the Fifth Batch of National Traditional Chinese Medicine Excellent Clinical Talents Training Project. (Announcement from the Personnel and Education Department of the National Administration of Traditional Chinese Medicine. No. 2022-1).
Acknowledgments
We thank all the studies for making the summary association statistics data publicly.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Generative AI statement
The author(s) declare that no Generative AI was used in the creation of this manuscript.
Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.
Publisher's note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcvm.2025.1554991/full#supplementary-material
Abbreviations
AAA, abdominal aortic aneurysm; AS, atherosclerosis; AID, sutoimmune disease; CI, confidence interval; DEGs, differentially expressed genes; eQTL, expression quantitative trait loci; FUMA, functional mapping and annotation; GWAS, genome-wide association study; HEIDI, heterogeneity in dependent instrument; HLA, human leukocyte antigen; HyPrColoc, hypothesis prioritization for multi-trait colocalization; H3Cit, citrullinated histone H3; KIR, killer cell immunoglobulin-like receptor; LD, linkage disequilibrium; LDSC; linkage disequilibrium score regression; MAF, minor allele frequency; MHC, major histocompatibility complex; MPO, myeloperoxidase; MAGMA, multi-marker analysis of genomic annotation; MSigDB, molecular signatures database; NETs, neutrophil extracellular traps; NK, natural killer cells; OR, odds ratio; PBC, primary biliary cholangitis; PF4, platelet factor 4; PLACO, pleiotropic analysis under composite null hypothesis; RA, rheumatoid arthritis; SCG2, secretory granule protein II; SE, standard error; SLE, systemic lupus erythematosus; SMR, summary-based mendelian randomization; SNP, single nucleotide polymorphism; SS, Sjogren's syndrome; S-LDSC, stratified-linkage disequilibrium score regression; TH, helper T cells; VSMC, vascular smooth muscle cell.
References
1. Sakalihasan N, Limet R, Defawe OD. Abdominal aortic aneurysm. Lancet. (2005) 365:1577–89. doi: 10.1016/S0140-6736(05)66459-8
2. Golledge J, Thanigaimani S, Powell JT, Tsao PS. Pathogenesis and management of abdominal aortic aneurysm. Eur Heart J. (2023) 44:2682–97. doi: 10.1093/eurheartj/ehad386
3. Chi X, Huang M, Tu H, Zhang B, Lin X, Xu H, et al. Innate and adaptive immune abnormalities underlying autoimmune diseases: the genetic connections. Sci China Life Sci. (2023) 66:1482–517. doi: 10.1007/s11427-021-2187-3
4. Huanggu H, Yang D, Zheng Y. Blood immunological profile of abdominal aortic aneurysm based on autoimmune injury. Autoimmun Rev. (2023) 22:103258. doi: 10.1016/j.autrev.2022.103258
5. Liu P, Sun Z, Zhang Y, Guo W. Myeloid related protein 8/14 is a new candidate biomarker and therapeutic target for abdominal aortic aneurysm. Biomed Pharmacother. (2019) 118:109229. doi: 10.1016/j.biopha.2019.109229
6. Wang SK, Green LA, Gutwein AR, Drucker NA, Motaganahalli RL, Gupta AK, et al. Description of human AAA by cytokine and immune cell aberrations compared to risk-factor matched controls. Surgery. (2018) 164:354–8. doi: 10.1016/j.surg.2018.03.002
7. Juvonen J, Surcel HM, Satta J, Teppo AM, Bloigu A, Syrjälä H, et al. Elevated circulating levels of inflammatory cytokines in patients with abdominal aortic aneurysm. Arterioscler Thromb Vasc Biol. (1997) 17:2843–7. doi: 10.1161/01.atv.17.11.2843
8. Wei Z, Wang Y, Zhang K, Liao Y, Ye P, Wu J, et al. Inhibiting the Th17/IL-17A-related inflammatory responses with digoxin confers protection against experimental abdominal aortic aneurysm. Arterioscler Thromb Vasc Biol. (2014) 34:2429–38. doi: 10.1161/ATVBAHA.114.304435
9. Brandau A, Ibrahim N, Klopf J, Hayden H, Ozsvar-Kozma M, Afonyushkin T, et al. Association of lipoproteins with neutrophil extracellular traps in patients with abdominal aortic aneurysm. Biomedicines. (2022) 10:217. doi: 10.3390/biomedicines10020217
10. Salata K, Almaghlouth I, Hussain MA, de Mestral C, Greco E, Aljabri BA, et al. Outcomes of abdominal aortic aneurysm repair among patients with rheumatoid arthritis. J Vasc Surg. (2021) 73:1261–8.e5. doi: 10.1016/j.jvs.2020.08.134
11. Nishimoto A, Haraguchi S, Kishi T, Fujisaki T, Ikeda Y, Sakemi T, et al. A case of systemic lupus erythematosus patient who died of rupture of abdominal aortic aneurysm. Nihon Jinzo Gakkai Shi. (2003) 45:701–5.14631766
12. Lu S, White JV, Nwaneshiudu I, Nwaneshiudu A, Monos DS, Solomides CC, et al. Human abdominal aortic aneurysm (AAA): evidence for an autoimmune antigen-driven disease. Autoimmun Rev. (2022) 21:103164. doi: 10.1016/j.autrev.2022.103164
13. Ibrahim N, Eilenberg W, Neumayer C, Brostjan C. Neutrophil extracellular traps in cardiovascular and aortic disease: a narrative review on molecular mechanisms and therapeutic targeting. Int J Mol Sci. (2024) 25:3983. doi: 10.3390/ijms25073983
14. Dubis J, Niepiekło-Miniewska W, Jędruchniewicz N, Sobczyński M, Witkiewicz W, Zapotoczny N, et al. Associations of genes for killer cell immunoglobulin-like receptors and their human leukocyte antigen-A/B/C ligands with abdominal aortic aneurysm. Cells. (2021) 10:3357. doi: 10.3390/cells10123357
15. Lei C, Yang D, Chen S, Chen W, Sun X, Wu X, et al. Patterns of immune infiltration in stable and raptured abdominal aortic aneurysms: a gene-expression-based retrospective study. Gene. (2020) 762:145056. doi: 10.1016/j.gene.2020.145056
16. Telomeres Mendelian Randomization Collaboration, Haycock PC, Burgess S, Nounu A, Zheng J, Okoli GN, et al. Association between telomere length and risk of cancer and non-neoplastic diseases: a Mendelian randomization study. JAMA Oncol. (2017) 3:636–51. doi: 10.1001/jamaoncol.2016.5945
17. Zhang L, Li Q, Zhou C, Zhang Z, Zhang J, Qin X. Immune-dysregulated neutrophils characterized by upregulation of CXCL1 may be a potential factor in the pathogenesis of abdominal aortic aneurysm and systemic lupus erythematosus. Heliyon. (2023) 9:e18037. doi: 10.1016/j.heliyon.2023.e18037
18. Kappelmann N, Arloth J, Georgakis MK, Czamara D, Rost N, Ligthart S, et al. Dissecting the association between inflammation, metabolic dysregulation, and specific depressive symptoms: a genetic correlation and 2-sample Mendelian randomization study. JAMA Psychiatry. (2021) 78:161–70. doi: 10.1001/jamapsychiatry.2020.3436
19. Lu H, Qiao J, Shao Z, Wang T, Huang S, Zeng P. A comprehensive gene-centric pleiotropic association analysis for 14 psychiatric disorders with GWAS summary statistics. BMC Med. (2021) 19:314. doi: 10.1186/s12916-021-02186-z
20. Sun S, Liu Y, Li L, Xiong L, Jiao M, Yang J, et al. Unveiling the shared genetic architecture between testosterone and polycystic ovary syndrome. Sci Rep. (2024) 14:23931. doi: 10.1038/s41598-024-75816-0
21. Zhang H, Zheng R, Yu B, Yu Y, Luo X, Yin S, et al. Dissecting shared genetic architecture between depression and body mass index. BMC Med. (2024) 22:455. doi: 10.1186/s12916-024-03681-9
22. Vanlerberghe C, Jourdain AS, Frenois F, Ait-Yahya E, Bamshad M, Dieux A, et al. Functional characterization vs in silico prediction for TBX5 missense and splice variants in holt-oram syndrome. Genet Med. (2024) 26:101267. doi: 10.1016/j.gim.2024.101267
23. Taylor KE, Wong Q, Levine DM, McHugh C, Laurie C, Doheny K, et al. Genome-Wide association analysis reveals genetic heterogeneity of sjögren’s syndrome according to ancestry. Arthritis Rheumatol. (2017) 69:1294–305. doi: 10.1002/art.40040
24. PubMed. A global reference for human genetic variation—PubMed. Available online at: https://pubmed.ncbi.nlm.nih.gov/26432245/ (Accessed September 9, 2025).
25. Bulik-Sullivan BK, Loh P-R, Finucane HK, Ripke S, Yang J, Schizophrenia Working Group of the Psychiatric Genomics Consortium, et al. LD Score regression distinguishes confounding from polygenicity in genome-wide association studies. Nat Genet. (2015) 47:291–5. doi: 10.1038/ng.3211
26. Park S, Kim S, Kim B, Kim DS, Kim J, Ahn Y, et al. Multivariate genomic analysis of 5 million people elucidates the genetic architecture of shared components of the metabolic syndrome. Nat Genet. (2024) 56:2380–91. doi: 10.1038/s41588-024-01933-1
27. GTEx Consortium. The GTEx consortium atlas of genetic regulatory effects across human tissues. Science. (2020) 369:1318–30. doi: 10.1126/science.aaz1776
28. de Leeuw CA, Mooij JM, Heskes T, Posthuma D. MAGMA: generalized gene-set analysis of GWAS data. PLoS Comput Biol. (2015) 11:e1004219. doi: 10.1371/journal.pcbi.1004219
29. Liberzon A, Birger C, Thorvaldsdóttir H, Ghandi M, Mesirov JP, Tamayo P. The molecular signatures database (MSigDB) hallmark gene set collection. Cell Syst. (2015) 1:417–25. doi: 10.1016/j.cels.2015.12.004
30. Huang DW, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. (2009) 4:44–57. doi: 10.1038/nprot.2008.211
31. Zhou Y, Zhou B, Pache L, Chang M, Khodabakhshi AH, Tanaseichuk O, et al. Metascape provides a biologist-oriented resource for the analysis of systems-level datasets. Nat Commun. (2019) 10:1523. doi: 10.1038/s41467-019-09234-6
32. Ray D, Chatterjee N. A powerful method for pleiotropic analysis under composite null hypothesis identifies novel shared loci between type 2 diabetes and prostate cancer. PLoS Genet. (2020) 16:e1009218. doi: 10.1371/journal.pgen.1009218
33. Watanabe K, Taskesen E, van Bochoven A, Posthuma D. Functional mapping and annotation of genetic associations with FUMA. Nat Commun. (2017) 8:1826. doi: 10.1038/s41467-017-01261-5
34. Giambartolomei C, Vukcevic D, Schadt EE, Franke L, Hingorani AD, Wallace C, et al. Bayesian Test for colocalisation between pairs of genetic association studies using summary statistics. PLoS Genet. (2014) 10:e1004383. doi: 10.1371/journal.pgen.1004383
35. Zhu Z, Zhang F, Hu H, Bakshi A, Robinson MR, Powell JE, et al. Integration of summary data from GWAS and eQTL studies predicts complex trait gene targets. Nat Genet. (2016) 48:481–7. doi: 10.1038/ng.3538
36. Foley CN, Staley JR, Breen PG, Sun BB, Kirk PDW, Burgess S, et al. A fast and efficient colocalization algorithm for identifying shared genetic risk factors across multiple traits. Nat Commun. (2021) 12:764. doi: 10.1038/s41467-020-20885-8
37. Orrù V, Steri M, Sidore C, Marongiu M, Serra V, Olla S, et al. Complex genetic signatures in immune cells underlie autoimmunity and inform therapy. Nat Genet. (2020) 52:1036–45. doi: 10.1038/s41588-020-0684-4
38. Li L, Tian W. The full-length genomic sequence of the HLA-B*46:01:25 allele. HLA. (2024) 104:e15770. doi: 10.1111/tan.15770
39. Dias E Silva D, Mambetsariev I, Fricke J, Babikian R, Dingal ST, Mazdisnian F, et al. A novel HLA-DQB2::MET gene fusion variant in lung adenocarcinoma with prolonged response to tepotinib: a case report. Transl Lung Cancer Res. (2024) 13:1163–8. doi: 10.21037/tlcr-24-34
40. Qin P, Huang H, Wang J, Jiang T, Zeng N, Wang Q, et al. The mechanism of LSM2 in the progression of live hepatocellular carcinoma was analyzed based on bioinformatics. Med Oncol. (2023) 40:276. doi: 10.1007/s12032-023-02152-0
41. Kaneko K, Nonomura Y, Watanabe K, Koike R, Kubota T, Harigai M, et al. Infected abdominal aortic aneurysm caused by nontyphoid Salmonella in an immunocompromised patient with rheumatoid arthritis. J Infect Chemother. (2009) 15:312–5. doi: 10.1007/s10156-009-0699-3
42. Tin K, Sobani ZA, Horovitz J, Rahmani R. Aortic and splanchnic artery aneurysms: unusual causes of biliary obstruction—a retrospective cohort from literature. Int J Surg. (2017) 39:163–8. doi: 10.1016/j.ijsu.2017.01.108
43. Taglietti L, Pontoglio S, Di Flumeri GR, Re P, Vettoretto N, Ghilardi G, et al. Acute thrombosis of abdominal aorta and hypercoagulable disorders. Int Angiol. (2008) 27:157–65.18427402
44. Perricone C, Bruno L, Cafaro G, Latini A, Ceccarelli F, Borgiani P, et al. Sjogren’s syndrome: everything you always wanted to know about genetic and epigenetic factors. Autoimmun Rev. (2024) 23:103673. doi: 10.1016/j.autrev.2024.103673
45. Moñux G, Serrano FJ, Vigil P, De la Concha EG. Role of HLA-DR in the pathogenesis of abdominal aortic aneurysm. Eur J Vasc Endovasc Surg. (2003) 26:211–4. doi: 10.1053/ejvs.2002.1897
46. Zong Y, Yang Y, Zhao J, Li L, Luo D, Hu J, et al. Characterisation of macrophage infiltration and polarisation based on integrated transcriptomic and histological analyses in primary sjögren’s syndrome. Front Immunol. (2023) 14:1292146. doi: 10.3389/fimmu.2023.1292146
47. Meissner J, Eysmont K, Matylla-Kulińska K, Konarska MM. Characterization of Cwc2, U6 snRNA, and Prp8 interactions destabilized by Prp16 ATPase at the transition between the first and second steps of splicing. RNA. (2024) 30:1199–212. doi: 10.1261/rna.079886.123
48. Atkinson R, Georgiou M, Yang C, Szymanska K, Lahat A, Vasconcelos EJR, et al. PRPF8-mediated Dysregulation of hBrr2 helicase disrupts human spliceosome kinetics and 5´-splice-site selection causing tissue-specific defects. Nat Commun. (2024) 15:3138. doi: 10.1038/s41467-024-47253-0
49. Lee H-S, Lee AT, Criswell LA, Seldin MF, Amos CI, Carulli JP, et al. Several regions in the major histocompatibility complex confer risk for anti-CCP-antibody positive rheumatoid arthritis, independent of the DRB1 locus. Mol Med. (2008) 14:293–300. doi: 10.2119/2007-00123.Lee
50. Wang Y, Broderick P, Webb E, Wu X, Vijayakrishnan J, Matakidou A, et al. Common 5p15.33 and 6p21.33 variants influence lung cancer risk. Nat Genet. (2008) 40:1407–9. doi: 10.1038/ng.273
51. Zhou X-J, Su T, Xie J, Xie Q-H, Wang L-Z, Hu Y, et al. Genome-Wide association study in acute tubulointerstitial nephritis. J Am Soc Nephrol. (2023) 34:895–908. doi: 10.1681/ASN.0000000000000091
52. Erdik N, Erdik A, Kizil D, Yavuzbilge G, Turk SM, Guneysu CA, et al. The relationship between neuropathic pain and lower urinary tract symptom scores in patients with primary sjögren’s syndrome. BMC Urol. (2024) 24:257. doi: 10.1186/s12894-024-01650-1
53. Moreno-Santos I, Macías-González M, Porras-Martín C, Castellano-Castillo D, Sánchez-Espín G, Gómez-Doblas JJ, et al. Role of epicardial adipose tissue NPR-C in acute coronary syndrome. Atherosclerosis. (2019) 286:79–87. doi: 10.1016/j.atherosclerosis.2019.05.010
54. Pang R, Ma X, Guo H, Qi X. Meta-analysis of mortality-associated factors in primary sjögren’s syndrome patients with interstitial lung disease. Clin Rheumatol. (2025) 44:23–31. doi: 10.1007/s10067-024-07191-0
55. Yatim KM, Lakkis FG. A brief journey through the immune system. Clin J Am Soc Nephrol. (2015) 10:1274–81. doi: 10.2215/CJN.10031014
56. Libby P. Inflammation in atherosclerosis. Arterioscler Thromb Vasc Biol. (2012) 32:2045–51. doi: 10.1161/ATVBAHA.108.179705
57. Zhou J, Dong C, Tan J, Wang G, Li Z, Li S, et al. Promoting effect and immunologic role of secretogranin II on bladder cancer progression via regulating MAPK and NF-κB pathways. Apoptosis. (2024) 29:121–41. doi: 10.1007/s10495-023-01898-2
58. Tellides G, Pober JS. Inflammatory and immune responses in the arterial media. Circ Res. (2015) 116:312–22. doi: 10.1161/CIRCRESAHA.116.301312
59. Saleh R, Taha RZ, Sasidharan Nair V, Toor SM, Alajez NM, Elkord E. Transcriptomic profiling of circulating HLA-DR- myeloid cells, compared with HLA-DR+ myeloid antigen-presenting cells. Immunol Invest. (2021) 50:952–63. doi: 10.1080/08820139.2020.1795875
60. Yang X, Garner LI, Zvyagin IV, Paley MA, Komech EA, Jude KM, et al. Autoimmunity-associated T cell receptors recognize HLA-B*27-bound peptides. Nature. (2022) 612:771–7. doi: 10.1038/s41586-022-05501-7
61. Zaimoku Y, Patel BA, Adams SD, Shalhoub R, Groarke EM, Lee AAC, et al. HLA Associations, somatic loss of HLA expression, and clinical outcomes in immune aplastic anemia. Blood. (2021) 138:2799–809. doi: 10.1182/blood.2021012895
62. Samra S, Bergerson JRE, Freeman AF, Turvey SE. JAK-STAT signaling pathway, immunodeficiency, inflammation, immune dysregulation, and inborn errors of immunity. J Allergy Clin Immunol. (2025) 155:357–67. doi: 10.1016/j.jaci.2024.09.020
63. Larson HJ, Eagan RL, Bowers AP, Smilingyte I, de Figueiredo A, Coulter A. Building public trust in preparation for future health shocks: a research agenda. Br Med J. (2024) 387:e078464. doi: 10.1136/bmj-2023-078464
64. Sarrand J, Soyfoo MS. Involvement of epithelial-mesenchymal transition (EMT) in autoimmune diseases. Int J Mol Sci. (2023) 24:14481. doi: 10.3390/ijms241914481
65. Peng Z, Lv S-J, Chen H, Rao H, Guo Z, Wan Q, et al. Disruption of PCSK9 suppresses inflammation and attenuates abdominal aortic aneurysm formation. Arterioscler Thromb Vasc Biol. (2025) 45:e1–14. doi: 10.1161/ATVBAHA.123.320391
66. Loppnow H, Buerke M, Werdan K, Rose-John S. Contribution of vascular cell-derived cytokines to innate and inflammatory pathways in atherogenesis. J Cell Mol Med. (2011) 15:484–500. doi: 10.1111/j.1582-4934.2010.01245.x
Keywords: abdominal aortic aneurysm, Sjögren's syndrome, genetic correlations, pleiotropic loci, immune colocalization
Citation: Zhou Y, Xiong H, Luo Q, Ding Z, He J and Wang L (2025) Genetic relationship between Sjögren's syndrome and abdominal aortic aneurysm: insights from a European population's genome-wide association analysis. Front. Cardiovasc. Med. 12:1554991. doi: 10.3389/fcvm.2025.1554991
Received: 3 January 2025; Accepted: 12 September 2025;
Published: 24 September 2025.
Edited by:
Jirko Kühnisch, Charité University Medicine Berlin, GermanyReviewed by:
Lingfeng Luo, Stanford University, United StatesYanhong Li, First Affiliated Hospital of Xinjiang Medical University, China
Copyright: © 2025 Zhou, Xiong, Luo, Ding, He and Wang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Jun He, NjQ5OTc3NDM3QHFxLmNvbQ==; Lihua Wang, MjcyMTMyMjE4MkBxcS5jb20=
†These authors have contributed equally to this work and share first authorship