A genomic association study revealing subphenotypes of childhood steroid-sensitive nephrotic syndrome in a larger genomic sequencing cohort

Dissecting the genetic components that contribute to the two main subphenotypes of steroid-sensitive nephrotic syndrome (SSNS) using genome-wide association studies (GWAS) strategy is important for understanding the disease. We conducted a multicenter cohort study (360 patients and 1835 controls) combined with a GWAS strategy to identify susceptibility variants associated with the following two subphenotypes of SSNS: steroid-sensitive nephrotic syndrome without relapse (SSNSWR, 181 patients) and steroid-dependent/frequent relapse nephrotic syndrome (SDNS/FRNS, 179 patients). The distribution of two single-nucleotide polymorphisms (SNPs) in ANKRD36 and ALPG was significant between SSNSWR and healthy controls, and that of two SNPs in GAD1 and HLA-DQA1 was significant between SDNS/FRNS and healthy controls. Interestingly, rs1047989 in HLA-DQA1 was a candidate locus for SDNS/FRNS but not for SSNSWR. No significant SNPs were observed between SSNSWR and SDNS/FRNS. Meanwhile, chromosome 2:171713702 in GAD1 was associated with a greater steroid dose (>0.75 mg/kg/d) upon relapse to first remission in patients with SDNS/FRNS (odds ratio = 3.14; 95% confidence interval, 0.97-9.87; P = 0.034). rs117014418 in APOL4 was significantly associated with a decrease in eGFR of greater than 20% compared with the baseline in SDNS/FRNS patients (P = 0.0001). Protein-protein intersection network construction suggested that HLA-DQA1 and HLA-DQB1 function together through GSDMA. Thus, SSNSWR belongs to non-HLA region-dependent nephropathy, and the HLA-DQA/DQB region is likely strongly associated with disease relapse, especially in SDNS/FRNS. The study provides a novel approach for the GWAS strategy of SSNS and contributes to our understanding of the pathological mechanisms of SSNSWR and SDNS/FRNS.

Multispectral and thermal infrared data, visual scores for severity of common rust symptoms, and genotypic single nucleotide polymorphism data of three F2-derived biparental doubled-haploid maize populations

Three F2-derived biparental doubled haploid (DH) maize populations were generated for genetic mapping of resistance to common rust. Each of the three populations has the same susceptible parent, but a different resistance donor parent. Population 1 and 3 consist of 320 lines each, population 2 consists of 260 lines. The DH lines were evaluated for their susceptibility to common rust in two years and with two replications in each year. For phenotyping, a visual score (VS) for susceptibility was assigned. Additionally, unmanned aerial vehicle (UAV) derived multispectral and thermal infrared data was recorded and combined in different vegetation indices ("remote sensing", RS). The DH lines were genotyped with the DarTseq method, to obtain data on single nucleotide polymorphisms (SNPs). After quality control, 9051 markers remained. Missing values were "imputed" by the empirical mean of the marker scores of the respective locus. We used the data for comparison of genome-wide association studies and genomic prediction when based on different phenotyping methods, that is either VS or RS data. The data may be interesting for reuse for instance for benchmarking genomic prediction models, for phytopathological studies addressing common rust, or for specifications of vegetation indices.

Enrichment of the Local Synaptic Translatome for Genetic Risk Associated With Schizophrenia and Autism Spectrum Disorder

BACKGROUND

Genes that encode synaptic proteins or messenger RNA targets of the RNA-binding protein FMRP (fragile X messenger ribonucleoprotein) have been linked to schizophrenia and autism spectrum disorder (ASD) through the enrichment of genetic variants that confer risk for these disorders. FMRP binds many transcripts with synaptic functions and is thought to regulate their local translation, a process that enables rapid and compartmentalized protein synthesis required for development and plasticity.

METHODS

We used summary statistics from large-scale genome-wide association studies of schizophrenia (74,776 cases, 101,023 controls) and ASD (18,381 cases, 27,969 controls) to test the hypothesis that the subset of synaptic genes that encode localized transcripts is more strongly associated with each disorder than nonlocalized transcripts. We also postulated that this subset of synaptic genes is responsible for associations attributed to FMRP targets.

RESULTS

Schizophrenia associations were enriched in genes encoding localized synaptic transcripts compared to the remaining synaptic genes or to the remaining localized transcripts; this also applied to ASD associations, although only for transcripts observed after stimulation by fear conditioning. The genetic associations with either disorder captured by these gene sets were independent of those derived from FMRP targets. Schizophrenia association was related to FMRP interactions with messenger RNAs in somata, but not in dendrites, while ASD association was related to FMRP binding in either compartment.

CONCLUSIONS

Our data suggest that synaptic transcripts capable of local translation are particularly relevant to the pathogenesis of schizophrenia and ASD, but they do not characterize the associations attributed to current sets of FMRP targets.

Causal effect of serum 25 hydroxyvitamin D concentration on cardioembolic stroke: Evidence from two-sample Mendelian randomization

BACKGROUND AND AIMS

The putative association between serum 25-hydroxyvitamin D concentration [25(OH)D] and the risk of cardioembolic stroke (CES) has been examined in observational studies, which indicate controversial findings. We performed Mendelian randomization (MR) analysis to determine the causal relationship of serum 25(OH)D with the risk of CES.

METHODS AND RESULTS

The summary statistics dataset on the genetic variants related to 25(OH)D was used from the published GWAS of European descent participants in the UK Biobank, including 417,580 subjects, yielding 143 independent loci in 112 1-Mb regions. GWAS summary data of CES was obtained from GIGASTROKE Consortium, which included European individuals (10,804 cases, 1,234,808 controls). Our results unveiled a causal relationship between 25(OH)D and CES using IVW [OR = 0.82, 95% CI: 0.67-0.98, p = 0.037]. Horizontal pleiotropy was not seen [MR-Egger intercept = 0.001; p = 0.792], suggesting an absence of horizontal pleiotropy. Cochrane's Q [Q = 78.71, p-value = 0.924], Rucker's Q [Q = 78.64, p-value = 0.913], and I2 = 0.0% (95% CI: 0.0%, 24.6%) statistic suggested no heterogeneity. This result remained consistent using different MR methods and sensitivity analyses, including Maximum likelihood [OR = 0.82, 95%CI: 0.67-0.98, p-value = 0.036], Constrained maximum likelihood [OR = 0.76, 95%CI: 0.64-0.90, p-value = 0.002], Debiased inverse-variance weighted [OR = 0.82, 95%CI: 0.68-0.99, p-value = 0.002], MR-PRESSO [OR = 0.82, 95%CI 0.77-0.87, p-value = 0.022], RAPS [OR = 0.82, 95%CI 0.67-0.98, p-value = 0.038], MR-Lasso [OR = 0.82, 95%CI 0.68-0.99, p-value = 0.037].

CONCLUSION

Our MR analysis provides suggestive evidence that increased 25(OH)D levels may play a protective role in the development of cardioembolic stroke. Determining the role of 25(OH)D in stroke subtypes has important clinical and public health implications.

Genetic evidence strengthens the bidirectional connection between oral health status and psychiatric disorders: A two-sample Mendelian randomization study

BACKGROUND

Observational studies cannot accurately infer the causal associations between oral health status and psychiatric disorders.

METHODS

We conducted univariate and multivariate Mendelian randomization (MR) analyses using single nucleotide polymorphisms (SNPs) associated with eight oral health statuses (periodontitis, DMFS, Nteeth, toothache, loose teeth, painful gums, bleeding gums, and mouth ulcers) and four psychiatric disorders (Schizophrenia, Major Depressive Disorder (MDD), anxiety and stress-related disorder (ASRD), and Bipolar Disorder (BIP)) as instrumental variables. Genetic data were sourced from the Gene-lifestyle interactions in dental endpoints (GLIDE), UK Biobank, Psychiatric Genomics Consortium (PGC), and Lundbeck Foundation Initiative for Integrative Psychiatric Research (iPSYCH). The inverse variance-weighted (IVW) approach, supported by a comprehensive sensitivity analysis, was employed.

RESULTS

Genetically predicted mouth ulcers were significantly linked to higher MDD (OR = 2.17, 95 % CI: 1.33--3.54, P< 0.01) and BIP risks (OR = 2.25, 95 % CI: 1.22-4.15, P = 0.01). BIP heightened bleeding gums risk (OR = 1.01, 95 % CI: 1.00-1.01, P < 0.01). These associations were adjusted for smoking status and alcohol consumption. Painful gums were significantly associated with MDD risk (OR = 96.48, 95 % CI: 2.66-3495.28, P = 0.01), while MDD raised periodontitis risk (OR = 2.15, 95 % CI: 1.24-3.75, P = 0.01), both confounded by smoking and alcohol. Relatively small effects between several variables, while others could not withstand correction for multiple tests.

LIMITATIONS

The sample size and limitation to European populations limits the study generalizability.

CONCLUSIONS

This study provide evidence of possible causal relationships between several oral health conditions and mental illness. Focusing on oral health and valuing mental health are important for each other and overall health.

Causal relationship between diabetes and depression: A bidirectional Mendelian randomization study

OBJECTIVE

This study explores the causal relationship between diabetes and depression using a two-sample Mendelian Randomization (TSMR) method.

METHODS

The study selected single nucleotide polymorphisms (SNPs) closely associated with diabetes and depression in European populations from the Genome-Wide Association Study (GWAS) database, to serve as instrumental variables (IVs). The main evaluation method was inverse variance weighted analysis (IVW), supplemented by verification using Weighted median, Weighted mode, and MR Egger methods. The Odds Ratio (OR) and 95 % Confidence Interval (CI) were used as the main evaluation indicators, along with sensitivity analysis.

RESULTS

This study found a negative correlation between diabetes and depression, suggesting that diabetes may reduce the risk of depression [IVW(FE): OR: 0.901, 95 % CI: 0.823 to 0.987; P = 0.025 < 0.05]. This finding was further confirmed by the Weighted median [OR: 0.844, 95 % CI: 0.730 to 0.974; P = 0.021 < 0.05] and Weighted mode method [OR: 0.766, 95 % CI: 0.637 to 0.921; P = 0.006 < 0.05]. However, the reverse showed no causal relationship between depression and diabetes (P > 0.05). Sensitivity analysis found no pleiotropy, and there were no large influences from individual SNPs on the result's robustness; the results are stable and reliable.

CONCLUSION

For the first time, this study using TSMR analysis found a negative correlation between diabetes and the risk of depression onset in European populations, suggesting that diabetes might reduce the risk of depression. But as the mechanisms are still unclear, these findings warrant further study.

Genome-wide analysis to uncover how Pocillopora acuta survives the challenging intertidal environment

Characterisation of genomic variation among corals can help uncover variants underlying trait differences and contribute towards genotype prioritisation in coastal restoration projects. For example, there is growing interest in identifying resilient genotypes for transplantation, and to better understand the genetic processes that allow some individuals to survive in specific conditions better than others. The coral species Pocillopora acuta is known to survive in a wide range of habitats, from reefs artificial coastal defences, suggesting its potential use as a starter species for ecological engineering efforts involving coral transplantation onto intertidal seawalls. However, the intertidal section of coastal armour is a challenging environment for corals, with conditions during periods of emersion being particularly stressful. Here, we scanned the entire genome of P. acuta corals to identify the regions harbouring single nucleotide polymorphisms (SNPs) and copy number variations (CNVs) that separate intertidal colonies (n = 18) from those found in subtidal areas (n = 21). Findings revealed 74,391 high quality SNPs distributed across 386 regions of the P. acuta genome. While the majority of the detected SNPs were in non-coding regions, 12% were identified in exons (i.e. coding regions). Functional SNPs that were significantly associated with intertidal colonies were found in overrepresented genomic regions linked to cellular homeostasis, metabolism, and signalling processes, which may represent local environmental adaptation in the intertidal. Interestingly, regions that exhibited CNVs were also associated with metabolic and signalling processes, suggesting P. acuta corals living in the intertidal have a high capacity to perform biological functions critical for survival in extreme environments.

Genome-wide association study identifies DRAM1 associated with Tourette syndrome in Taiwan

BACKGROUND

Tourette syndrome (TS) is a neurodevelopmental disorder characterized by motor and vocal tics. Several susceptibility loci associated with TS have been identified previously in populations of European descent using genome-wide association studies (GWAS). However, the exact pathogenic mechanism underlying TS is unknown; additionally, the results of previous GWAS for TS were based on Western populations, which may not translate to other populations. Therefore, we conducted a GWAS in Taiwanese patients with TS and chronic tic disorders (CTDs), with an aim to elucidate the genetic basis and potential risk factors for TS in this population.

METHODS

GWAS was performed on a Taiwanese TS/CTDs cohort with a sample size of 1,007 patients with TS and 25,522 ancestry-matched controls. Additionally, polygenic risk score was calculated and assessed.

RESULTS

Genome-wide significant locus, rs12313062 (p=1.43 × 10-8) and other 9 single nucleotide polymorphisms, were identified in chromosomes 12q23.2, associated with DRAM1 and was a novel susceptibility locus identified in TS/CTDs group. DRAM1, a lysosomal transmembrane protein regulated by p53, modulates autophagy and apoptosis, with potential implications for neuropsychiatric conditions associated with autophagy disruption.

CONCLUSIONS

This study conducted the first GWAS for TS in a Taiwanese population, identifying a significant locus on chromosome 12q23.2 associated with DRAM1. These findings provide novel insights into the neurobiology of TS and potential directions for future research in this area.

From genomic insights to clinical hope: Targeting NEU1 in IgA nephropathy

BACKGROUND

IgA Nephropathy (IgAN), the primary form of glomerulonephritis, presents significant clinical challenges due to its obscure pathogenesis and lack of targeted treatments. We conducted a proteome-wide Mendelian randomization (MR) study to identify therapeutic targets for IgAN.

METHODS

Utilizing a plasma proteome dataset comprising 4907 blood plasma proteins as the exposure variable, and renal biopsy-confirmed IgAN cases as the outcome, this study employed MR to pinpoint proteins potentially pathogenic to IgAN. The robustness of our findings was affirmed through external dataset validation, reverse causation testing, and Bayesian colocalization analysis. Additionally, we conducted phenotypic scanning and analyzed downstream metabolites to investigate candidate proteins's biological function.

RESULTS

In our study, a significant association was identified between an increase in neuraminidase 1 (NEU1) expression and the risk of IgAN. Specifically, a one standard deviation increase in NEU1 expression was associated with an odds ratio of 11.80 for the development of IgAN (95% confidence interval: 4.03-34.54). This association was substantiated across various statistical models and external validations. Colocalization analysis indicated a shared causal variant between NEU1 expression and IgAN. Furthermore, an increased influenza risk associated with NEU1 was observed, supporting the therapeutic potential of NEU1 inhibitors for IgAN. However, our study found no significant role for neuraminic acid-related metabolites in IgAN's development, suggesting an independent pathway for NEU1's influence.

CONCLUSION

This study identifies NEU1 as a promising therapeutic target for IgAN, backed by robust genetic evidence. Future research should explore NEU1's therapeutic potential in diverse populations and clinical scenarios, further establishing its role in IgAN.

Genetic influences on testosterone and PTSD

Females are twice as likely to experience PTSD as compared to males. Although sex differences in prevalence are well-established, little is known about why such sex differences occur. Biological factors that vary with sex, including sex hormone production, may contribute to these differences. Considerable evidence links sex hormones, such as testosterone, to PTSD risk though less is known about the shared genetic underpinnings. The objective of the present study was to test for genetic relationships between testosterone and PTSD. To do so, we used summary statistics from large, publicly available genetic consortia to conduct linkage disequilibrium score regression to estimate the genetic correlations between PTSD and testosterone in males and females, and two-sample, bi-directional Mendelian randomization to examine potential causal relationships of testosterone on PTSD and the reverse. Heritability estimates of testosterone were significantly higher in males (0.17, SE = 0.02) than females (0.11, SE = 0.01; z = 2.46, p = 00.01). The correlation between testosterone and PTSD was negative in males (rg = -0.11, SE = 0.02, p = 6.7 x 10-6), but not significant in females (rg = 0.002, SE = 0.03, p = 0.95). MR analyses found no evidence of a causal effect of testosterone on PTSD or the reverse. Findings are consistent with phenotypic literature suggesting a relationship between testosterone and PTSD that may be sex-specific. This work provides early evidence of a relationship between testosterone and PTSD genotypically and suggests an avenue for future research that will enable a better understanding of disparities in PTSD.

Causal effect of gut microbiota on the risk of prostatitis: a two-sample Mendelian randomization study

BACKGROUND

Recent studies demonstrated that chronic prostatitis (CP) is closely related to the gut microbiota (GM). Nevertheless, the causal relationship between GM and CP has not been fully elucidated. Therefore, the two-sample Mendelian randomization (MR) analysis was employed to investigate this association.

METHODS

The summary data of gut microbiota derived from a genome-wide association study (GWAS) involving 18,340 individuals in the MiBioGen study served as the exposure, and the corresponding summary statistics for CP risk, representing the outcome, were obtained from the FinnGen databases (R9). The causal effects between GM and CP were estimated using the inverse-variance weighted (IVW) method supplemented with MR-Egger, weighted median, weighted mode, and simple mode methods. Additionally, the false discovery rate (FDR) correction was performed to adjust results. The detection and quantification of heterogeneity and pleiotropy were accomplished through the MR pleiotropy residual sum and outlier method, Cochran's Q statistics, and MR-Egger regression.

RESULTS

The IVW estimates indicated that a total of 11 GM taxa were related to the risk of CP. Seven of them was correlated with an increased risk of CP, while the remained linked with a decreased risk of CP. However, only Methanobacteria (OR 0.86; 95% CI 0.74-0.99), Methanobacteriales (OR 0.86; 95% CI 0.74-0.99), NB1n (OR 1.16; 95% CI 1.16-1.34), Methanobacteriaceae (OR 0.86; 95% CI 0.74-0.99), Odoribactergenus Odoribacter (OR 1.43; 95% CI 1.05-1.94), and Sutterellagenus Sutterella (OR 1.33; 95% CI 1.01-1.76) still maintain significant association with CP after FDR correction. Consistent directional effects for all analyses were observed in the supplementary methods. Subsequently, sensitivity analyses indicated the absence of heterogeneity, directional pleiotropy, or outliers concerning the causal effect of specific gut microbiota on CP (p > 0.05).

CONCLUSION

Our study demonstrated a gut microbiota-prostate axis, offering crucial data supporting the promising use of the GM as a candidate target for CP prevention, diagnosis, and treatment. There is a necessity for randomized controlled trials to validate the protective effect of the linked GM against the risk of CP, and to further investigate the underlying mechanisms involved.

The casual relationship between autoimmune diseases and multiple myeloma: a Mendelian randomization study

Observational studies showed possible associations between systemic lupus erythematosus and multiple myeloma. However, whether there is a casual relationship between different types of autoimmune diseases (type 1 diabetes mellitus, rheumatoid arthritis, systemic lupus erythematosus, psoriasis, multiple sclerosis, primary sclerosing cholangitis, primary biliary cirrhosis, and juvenile idiopathic arthritis) and multiple myeloma (MM) is not well known. We performed a two-sample Mendelian randomization (MR) study to estimate the casual relationship. Summary-level data of autoimmune diseases were gained from published genome-wide association studies while data of MM was obtained from UKBiobank. The Inverse-Variance Weighted (IVW) method was used as the primary analysis method to interpret the study results, with MR-Egger and weighted median as complementary methods of analysis. There is causal relationship between primary sclerosing cholangitis [OR = 1.00015, 95% CI 1.000048-1.000254, P = 0.004] and MM. Nevertheless, no similar causal relationship was found between the remaining seven autoimmune diseases and MM. Considering the important role of age at recruitment and body mass index (BMI) in MM, we excluded these relevant instrument variables, and similar results were obtained. The accuracy and robustness of these findings were confirmed by sensitivity tests. Overall, MR analysis suggests that genetic liability to primary sclerosing cholangitis could be causally related to the increasing risk of MM. This finding may serve as a guide for clinical attention to patients with autoimmune diseases and their early screening for MM.

Genome-wide Association Study of Susceptibility Loci for Self-Reported Atopic Dermatitis and Allergic Rhinitis in the Korean Population

BACKGROUND

Allergic diseases include atopic dermatitis (AD) and allergic rhinitis (AR), which are chronic, relapsing inflammatory disorders of the skin or mucosa that usually accompany immunoglobulin E-mediated immune responses. They are complex, multifactorial diseases with an etiology involving interactions between genetic and environmental factors.

OBJECTIVE

We performed a genome-wide association study (GWAS) to identify single nucleotide polymorphisms (SNPs) associated with allergic diseases in the Korean population.

METHODS

A total of 8,840 samples were obtained from the Korean Association Resource Consortium dataset of the Korean Genome and Epidemiology Study Ansan-Anseong cohort. The allergic disease phenotype was determined based on self-reported physician diagnoses. After quality control, 8,823 subjects with 877,242 variants remained for the final analysis. The GWAS was performed using logistic regression analysis in an additive model adjusted for age and sex.

RESULTS

A total of 636 patients with allergic disease and 8,176 controls were analyzed. Three SNPs were associated with allergic disease at a level of genome-wide suggestive significance (p<1.0×10-5) in the Korean population: rs7275360, located in neural cell adhesion molecule 2; rs698195; and rs3750552, located in family with sequence similarity 189, member A2. These polymorphisms were on chromosomes 21q21.1, 7q31.1, and 9q21.12, respectively.

CONCLUSION

We identified 3 novel SNPs significantly associated with allergic diseases in the Korean population. Further research is required to confirm the association between these novel SNPs and allergic disease in the Korean population and in other ethnicities.

Meta-analysis of ACE inhibitor-induced angioedema identifies novel risk locus

BACKGROUND

Angioedema is a rare but potentially life-threatening adverse drug reaction in patients receiving angiotensin-converting enzyme inhibitors (ACEis). Research suggests that susceptibility to ACEi-induced angioedema (ACEi-AE) involves both genetic and nongenetic risk factors. Genome- and exome-wide studies of ACEi-AE have identified the first genetic risk loci. However, understanding of the underlying pathophysiology remains limited.

OBJECTIVE

We sought to identify further genetic factors of ACEi-AE to eventually gain a deeper understanding of its pathophysiology.

METHODS

By combining data from 8 cohorts, a genome-wide association study meta-analysis was performed in more than 1000 European patients with ACEi-AE. Secondary bioinformatic analyses were conducted to fine-map associated loci, identify relevant genes and pathways, and assess the genetic overlap between ACEi-AE and other traits. Finally, an exploratory cross-ancestry analysis was performed to assess shared genetic factors in European and African-American patients with ACEi-AE.

RESULTS

Three genome-wide significant risk loci were identified. One of these, located on chromosome 20q11.22, has not been implicated previously in ACEi-AE. Integrative secondary analyses highlighted previously reported genes (BDKRB2 [bradykinin receptor B2] and F5 [coagulation factor 5]) as well as biologically plausible novel candidate genes (PROCR [protein C receptor] and EDEM2 [endoplasmic reticulum degradation enhancing alpha-mannosidase like protein 2]). Lead variants at the risk loci were found with similar effect sizes and directions in an African-American cohort.

CONCLUSIONS

The present results contributed to a deeper understanding of the pathophysiology of ACEi-AE by (1) providing further evidence for the involvement of bradykinin signaling and coagulation pathways and (2) suggesting, for the first time, the involvement of the fibrinolysis pathway in this adverse drug reaction. An exploratory cross-ancestry comparison implicated the relevance of the associated risk loci across diverse ancestries.

Genome-wide association and Mendelian randomisation analysis among 30,699 Chinese pregnant women identifies novel genetic and molecular risk factors for gestational diabetes and glycaemic traits

AIMS/HYPOTHESIS

Gestational diabetes mellitus (GDM) is the most common disorder in pregnancy; however, its underlying causes remain obscure. This study aimed to investigate the genetic and molecular risk factors contributing to GDM and glycaemic traits.

METHODS

We collected non-invasive prenatal test (NIPT) sequencing data along with four glycaemic and 55 biochemical measurements from 30,699 pregnant women during a 2 year period at Shenzhen Baoan Women's and Children's Hospital in China. Genome-wide association studies (GWAS) were conducted between genotypes derived from NIPTs and GDM diagnosis, baseline glycaemic levels and glycaemic levels after glucose challenges. In total, 3317 women were diagnosed with GDM, while 19,565 served as control participants. The results were replicated using two independent cohorts. Additionally, we performed one-sample Mendelian randomisation to explore potential causal associations between the 55 biochemical measurements and risk of GDM and glycaemic levels.

RESULTS

We identified four genetic loci significantly associated with GDM susceptibility. Among these, MTNR1B exhibited the highest significance (rs10830963-G, OR [95% CI] 1.57 [1.45, 1.70], p=4.42×10-29), although its effect on type 2 diabetes was modest. Furthermore, we found 31 genetic loci, including 14 novel loci, that were significantly associated with the four glycaemic traits. The replication rates of these associations with GDM, fasting plasma glucose levels and 0 h, 1 h and 2 h OGTT glucose levels were 4 out of 4, 6 out of 9, 10 out of 11, 5 out of 7 and 4 out of 4, respectively. Mendelian randomisation analysis suggested that a genetically regulated higher lymphocytes percentage and lower white blood cell count, neutrophil percentage and absolute neutrophil count were associated with elevated glucose levels and an increased risk of GDM.

CONCLUSIONS/INTERPRETATION

Our findings provide new insights into the genetic basis of GDM and glycaemic traits during pregnancy in an East Asian population and highlight the potential role of inflammatory pathways in the aetiology of GDM and variations in glycaemic levels.

DATA AVAILABILITY

Summary statistics for GDM; fasting plasma glucose; 0 h, 1 h and 2h OGTT; and the 55 biomarkers are available in the GWAS Atlas (study accession no.: GVP000001, https://ngdc.cncb.ac.cn/gwas/browse/GVP000001) .

Exploring the impact of gut microbiota on abdominal aortic aneurysm risk through a bidirectional Mendelian randomization analysis

OBJECTIVE

The abdominal aortic aneurysm (AAA) is associated with alterations in the composition of the gut microbiota; however, the precise causal relationship remains unclear. Elucidating this complex interplay could provide new insights into the pathogenesis of AAA.

METHODS

A bidirectional two-sample Mendelian randomization analysis was conducted using genome-wide association study summary data on the gut microbiota (n = 18,340) and AAA (n = 353,087). A total of 196 gut microbial taxa across taxonomic levels were examined for their potential causal effects on AAA risk. Conversely, the effect of AAA on these microbial taxa was also analyzed.

RESULTS

Nine microbial taxa were identified as having a causal influence on AAA risk. Specifically, increased risk were associated with genus Bilophila (odds ratio [OR], 1.359; P = .0119), genus Catenibacterium (OR, 1.348; P = .0058), genus family XIII AD3011 group (OR, 1.507; P = .004), genus Oxalobacter (OR, 1.157; P = .0449), and genus Prevotella 7 (OR, 1.194; P = .0306), whereas decreased risks were linked to class Lentisphaeria (OR, 0.829; P = .0361), order Victivallales (OR, 0.829; P = .0361), family Victivallaceae (OR, 0.814; P = .0057), and genus Anaerotruncus (OR, 0.773; P = .0497). Furthermore, AAA was found to influence the abundance of 14 microbial taxa across various taxonomic levels. Notably, bidirectional associations were observed with the class Lentisphaeria and the order Victivallales.

CONCLUSIONS

This study provides novel evidence for a reciprocal causal relationship between gut microbiota and AAA risk, thereby offering new insights into the pathogenesis of AAA. These findings also suggest promising avenues for microbiome-based therapeutic interventions.

Shared genetic architecture between hypothyroidism and rheumatoid arthritis: A large-scale cross-trait analysis

BACKGROUND

In recent years, mounting evidence has indicated a co-morbid relationship between hypothyroidism and rheumatoid arthritis (RA), however, the shared genetic factors underlying this association remain unclear. This study aims to investigate the common genetic architecture between hypothyroidism and RA.

METHODS

Genome-wide association study (GWAS) summary statistics from recently published studies were utilized to examine the genetic correlation, shared genetic loci, and potential causal relationship between hypothyroidism and RA. Statistical methods included linkage disequilibrium score regression (LDSC), high-definition likelihood (HDL), cross-trait meta-analyses, colocalization analysis, multi-marker analysis of genomic annotation (MAGMA), tissue-specific enrichment analysis (TSEA), functional enrichment analysis, and latent causal variable method (LCV).

RESULTS

Our study demonstrated a significant genetic correlation between hypothyroidism and RA(LDSC:rg=0.3803,p=7.23e-11;HDL:rg=0.3849,p=1.02e-21). Through cross-trait meta-analysis, we identified 1035 loci, including 43 novel genetic loci. By integrating colocalization analysis and the MAGMA algorithm, we found a substantial number of genes, such as PTPN22, TYK2, and CTLA-4, shared between the two diseases, which showed significant enrichment across 14 tissues. These genes were primarily associated with the regulation of alpha-beta T cell proliferation, positive regulation of T cell activation, positive regulation of leukocyte cell-cell adhesion, T cell receptor signaling pathway, and JAK-STAT signaling pathway. However, our study did not reveal a significant causal association between the two diseases using the LCV approach.

CONCLUSION

Based on these findings, there is a significant genetic correlation between hypothyroidism and RA, suggesting a shared genetic basis for these conditions.

Decoding and reconstructing disease relations between dry eye and depression: a multimodal investigation comprising meta-analysis, genetic pathways and Mendelian randomization

INTRODUCTION

The clinical presentations of dry eye disease (DED) and depression (DEP) often comanifest. However, the robustness and the mechanisms underlying this association were undetermined.

OBJECTIVES

To this end, we set up a three-segment study that employed multimodality results (meta-analysis, genome-wide association study [GWAS] and Mendelian randomization [MR]) to elucidate the association, common pathways and causality between DED and DEP.

METHODS

A meta-analysis comprising 26 case-control studies was first conducted to confirm the DED-DEP association. Next, we performed a linkage disequilibrium (LD)-adjusted GWAS and targeted phenotype association study (PheWAS) in East Asian TW Biobank (TWB) and European UK Biobank (UKB) populations. Single-nucleotide polymorphisms (SNPs) were further screened for molecular interactions and common pathways at the functional gene level. To further elucidate the activated pathways in DED and DEP, a systemic transcriptome review was conducted on RNA sequencing samples from the Gene Expression Omnibus. Finally, 48 MR experiments were implemented to examine the bidirectional causation between DED and DEP.

RESULTS

Our meta-analysis showed that DED patients are associated with an increased DEP prevalence (OR = 1.83), while DEP patients have a concurrent higher risk of DED (OR = 2.34). Notably, cross-disease GWAS analysis revealed that similar genetic architecture (rG = 0.19) and pleiotropic functional genes contributed to phenotypes in both diseases. Through protein-protein interaction and ontology convergence, we summarized the pleiotropic functional genes under the ontology of immune activation, which was further validated by a transcriptome systemic review. Importantly, the inverse variance-weighted (IVW)-MR experiments in both TWB and UKB populations (p value <0.001) supported the bidirectional exposure-outcome causation for DED-to-DEP and DEP-to-DED. Despite stringent LD-corrected instrumental variable re-selection, the bidirectional causation between DED and DEP remained.

CONCLUSION

With the multi-modal evidence combined, we consolidated the association and causation between DED and DEP.

Phenome-wide causal associations between osteoarthritis and other complex traits through the latent causal variable analysis

BACKGROUND

Individuals with osteoarthritis present with comorbidities, and the potential causal associations remain incompletely elucidated. The present study undertook a large-scale investigation about the causality between osteoarthritis and variable traits, using the summary-level data of genome-wide association studies (GWAS).

METHODS

The present study included the summary-level GWS data of knee osteoarthritis, hip osteoarthritis, hip or knee osteoarthritis, hand osteoarthritis, and other 1355 traits. Genetic correlation analysis was conducted between osteoarthritis and other traits through cross-trait bivariate linkage disequilibrium score regression. Subsequently, latent causal variable analysis was performed to explore the causal association when there was a significant genetic correlation. Genetic correlation and latent causal variable analysis were conducted on the Complex Traits Genomics Virtual Lab platform ( https://vl.genoma.io/ ).

RESULTS

We found 133 unique phenotypes showing causal relationships with osteoarthritis. Our results confirmed several well-established risk factors of osteoarthritis, such as obesity, weight, BMI, and meniscus derangement. Additionally, our findings suggested putative causal links between osteoarthritis and multiple factors. Socioeconomic determinants such as occupational exposure to dust and diesel exhaust, extended work hours exceeding 40 per week, and unemployment status were implicated. Furthermore, our analysis revealed causal associations with cardiovascular and metabolic disorders, including heart failure, deep venous thrombosis, type 2 diabetes mellitus, and elevated cholesterol levels. Soft tissue and musculoskeletal disorders, such as hallux valgus, internal derangement of the knee, and spondylitis, were also identified to be causally related to osteoarthritis. The study also identified the putative causal associations of osteoarthritis with digestive and respiratory diseases, such as Barrett's esophagus, esophagitis, and asthma, as well as psychiatric conditions including panic attacks and manic or hyperactive episodes. Additionally, we observed osteoarthritis causally related to pharmacological treatments, such as the use of antihypertensive medications, anti-asthmatic drugs, and antidepressants.

CONCLUSION

Our study uncovered a wide range of traits causally associated with osteoarthritis. Further studies are needed to validate and illustrate the detailed mechanism of those causal associations.

Characterization of cell cycle, inflammation, and oxidative stress signaling role in non-communicable diseases: Insights into genetic variants, microRNAs and pathways

Non-Communicable Diseases (NCDs) significantly impact global health, contributing to over 70% of premature deaths, as reported by the World Health Organization (WHO). These diseases have complex and multifactorial origins, involving genetic, epigenetic, environmental and lifestyle factors. While Genome-Wide Association Study (GWAS) is widely recognized as a valuable tool for identifying variants associated with complex phenotypes; the multifactorial nature of NCDs necessitates a more comprehensive exploration, encompassing not only the genetic but also the epigenetic aspect. For this purpose, we employed a bioinformatics-multiomics approach to examine the genetic and epigenetic characteristics of NCDs (i.e. colorectal cancer, coronary atherosclerosis, squamous cell lung cancer, psoriasis, type 2 diabetes, and multiple sclerosis), aiming to identify novel biomarkers for diagnosis and prognosis. Leveraging GWAS summary statistics, we pinpointed Single Nucleotide Polymorphisms (SNPs) independently associated with each NCD. Subsequently, we identified genes linked to cell cycle, inflammation and oxidative stress mechanisms, revealing shared genes across multiple diseases, suggesting common functional pathways. From an epigenetic perspective, we identified microRNAs (miRNAs) with regulatory functions targeting these genes of interest. Our findings underscore critical genetic pathways implicated in these diseases. In colorectal cancer, the dysregulation of the "Cytokine Signaling in Immune System" pathway, involving LAMA5 and SMAD7, regulated by Hsa-miR-21-5p, Hsa-miR-103a-3p, and Hsa-miR-195-5p, emerged as pivotal. In coronary atherosclerosis, the pathway associated with "binding of TCF/LEF:CTNNB1 to target gene promoters" displayed noteworthy implications, with the MYC factor controlled by Hsa-miR-16-5p as a potential regulatory factor. Squamous cell lung carcinoma analysis revealed significant pathways such as "PTK6 promotes HIF1A stabilization," regulated by Hsa-let-7b-5p. In psoriasis, the "Endosomal/Vacuolar pathway," involving HLA-C and Hsa-miR-148a-3p and Hsa-miR-148b-3p, was identified as crucial. Type 2 Diabetes implicated the "Regulation of TP53 Expression" pathway, controlled by Hsa-miR-106a-5p and Hsa-miR-106b-5p. In conclusion, our study elucidates the genetic framework and molecular mechanisms underlying NCDs, offering crucial insights into potential genetic/epigenetic biomarkers for diagnosis and prognosis. The specificity of pathways and related miRNAs in different pathologies highlights promising candidates for further clinical validation, with the potential to advance personalized treatments and alleviate the global burden of NCDs.

Causal association between cardiovascular proteins and membranous nephropathy: a bidirectional Mendelian randomization

PURPOSE

Multiple circulating proteins have been reported to participate in human diseases. However, the association between cardiovascular proteins and membranous nephropathy (MN) remained profoundly elusive.

METHODS

A bidirectional Mendelian randomization (MR) analysis was conducted to explore the causal correlation between ninety cardiovascular proteins and MN. Genome-wide association study (GWAS) data of cardiovascular proteins and MN were all from European research. Inverse variance weighted (IVW) was used as the main approach. Moreover, MR-Egger, weighted median, weighted mode, and simple mode were also performed. Cochrane's Q test, MR-Egger, and MR-PRESSO were conducted for sensitivity analysis.

RESULTS

According to IVW method, fatty acid-binding protein and thrombomodulin (TM) were identified as risk factors for MN, while a protective role was detected in tissue-type plasminogen activator. Additionally, MN was associated with an elevated level of macrophage colony-stimulating factor 1, stem cell factor, TM, and tissue factor. Reversely, MN was also correlated with a downregulated level of beta-nerve growth factor, Cathepsin D, hepatocyte growth factor, interleukin-6 receptor subunit alpha, macrophage colony-stimulating factor 1, and myeloperoxidase. In the sensitivity analysis, no significant pleiotropy and heterogeneity was detected.

CONCLUSION

This was the first study to reveal the causal association between cardiovascular proteins and MN. These specific cardiovascular proteins could be novel biomarkers for MN, and is helpful for timely identify the risk of other diseases that might result from MN. However, further clinical studies are needed to confirm our results.

Causal relationship between the immune phenotype of monocytes and myasthenia gravis: A Mendelian randomization study

Background

Monocytes play an essential role in developing autoimmune diseases; however, their association with myasthenia gravis (MG) development is unclear.

Methods

We performed a two-sample Mendelian randomization analysis to assess the causal relationship between monocyte-associated traits and MG, reviewing summary statistics of genome-wide association studies (GWAS).

Results

Using the inverse variance weighted method, the following were found to be causally associated with MG: HLA-DR on monocytes (OR, 1.363; 95% CI, 1.158-1.605; P = 2E-04), HLA-DR on CD14+ monocytes (OR, 1.324; 95% CI, 1.183-1.482; P = 1.08E-06), HLA-DR on CD14+CD16- monocytes (OR, 1.313; 95% CI, 1.177-1.465; P = 1.07E-06), CD40 on monocytes (OR, 1.135; 95% CI, 1.012-1.272; P < 0.05), CD40 on CD14+CD16- monocytes (OR, 1.142; 95% CI, 1.015-1.285; P < 0.05), CD40 on CD14+CD16+ monocytes (OR, 1.142; 95% CI, 1.021-1.278; P < 0.05), CD64 on CD14+CD16+ monocytes (OR, 1.286; 95% CI, 1.019-1.623; P < 0.05).

Conclusions

The present study suggests a causal relationship between the upregulation of CD40, HLA-DR, and CD64 on monocytes and the development of MG. Altered monocyte function may potentially be a risk factor for MG and a therapeutic target.

Causal association of obesity with epigenetic aging and telomere length: a bidirectional mendelian randomization study

BACKGROUND

In observational studies, there exists an association between obesity and epigenetic age as well as telomere length. However, varying and partially conflicting outcomes have notably arisen from distinct studies on this topic. In the present study, two-way Mendelian randomization was used to identify potential causal associations between obesity and epigenetic age and telomeres.

METHODS

A genome-wide association study was conducted using data from individuals of European ancestry to investigate bidirectional Mendelian randomization (MR) regarding the causal relationships between obesity, as indicated by three obesity indicators (body mass index or BMI, waist circumference adjusted for BMI or WCadjBMI, and waist-to-hip ratio adjusted for BMI or WHRadjBMI), and four epigenetic age measures (HannumAge, HorvathAge, GrimAge, PhenoAge), as well as telomere length. To assess these causal associations, various statistical methods were employed, including Inverse Variance Weighted (IVW), Weighted Median, MR Egger, Weighted Mode, and Simple Mode. To address the issue of multiple testing, we applied the Bonferroni correction. These methods were used to determine whether there is a causal link between obesity and epigenetic age, as well as telomere length, and to explore potential bidirectional relationships. Forest plots and scatter plots were generated to show causal associations between exposures and outcomes. For a comprehensive visualization of the results, leave-one-out sensitivity analysis plots, individual SNP-based forest plots for MR analysis, and funnel plots were included in the presentation of the results.

RESULTS

A strong causal association was identified between obesity and accelerated HannumAge, GrimAge, PhenoAge and telomere length shrinkage. The causal relationship between WCadjBMI and PhenoAge acceleration (OR: 2.099, 95%CI: 1.248-3.531, p = 0.005) was the strongest among them. However, only the p-values for the causal associations of obesity with GrimAge, PhenoAge, and telomere length met the criteria after correction using the Bonferroni multiple test. In the reverse MR analysis, there were statistically significant causal associations between HorvathAge, PhenoAge and GrimAge and BMI, but these associations exhibited lower effect sizes, as indicated by their Odds Ratios (ORs). Notably, sensitivity analysis revealed the robustness of the study results.

CONCLUSIONS

The present findings reveal a causal relationship between obesity and the acceleration of epigenetic aging as well as the reduction of telomere length, offering valuable insights for further scientific investigations aimed at developing strategies to mitigate the aging process in humans.

Epigenetic age acceleration and risk of aortic valve stenosis: a bidirectional Mendelian randomization study

BACKGROUND

Aortic valve stenosis (AVS) is the most prevalent cardiac valve lesion in developed countries, and pathogenesis is closely related to aging. DNA methylation-based epigenetic clock is now recognized as highly accurate predictor of the aging process and associated health outcomes. This study aimed to explore the causal relationship between epigenetic clock and AVS by conducting a bidirectional Mendelian randomization (MR) analysis.

METHODS

Summary genome-wide association study statistics of epigenetic clocks (HannumAge, HorvathAge, PhenoAge, and GrimAge) and AVS were obtained and assessed for significant instrumental variables from Edinburgh DataShare (n = 34,710) and FinnGen biobank (cases = 9870 and controls = 402,311). The causal association between epigenetic clock and AVS was evaluated using inverse variance weighted (IVW), weighted median (WM), and MR-Egger methods. Multiple analyses (heterogeneity analysis, pleiotropy analysis, and sensitivity analysis) were performed for quality control assessment.

RESULTS

The MR analysis showed that the epigenetic age acceleration of HorvathAge and PhenoAge was associated with an increased risk of AVS (HorvathAge: OR = 1.043, P = 0.016 by IVW, OR = 1.058, P = 0.018 by WM; PhenoAge: OR = 1.058, P = 0.005 by IVW, OR = 1.053, P = 0.039 by WM). Quality control assessment proved our findings were reliable and robust. However, there was a lack of evidence supporting a causal link from AVS to epigenetic aging.

CONCLUSION

The present MR analysis unveiled a causal association between epigenetic clocks, especially HorvathAge and PhenoAge, with AVS. Further research is required to elucidate the underlying mechanisms and develop strategies for potential interventions.

Genome-wide association studies highlight novel risk loci for septal defects and left-sided congenital heart defects

BACKGROUND

Congenital heart defects (CHD) are structural defects of the heart affecting approximately 1% of newborns. They exhibit low penetrance and non-Mendelian patterns of inheritance as varied and complex traits. While genetic factors are known to play an important role in the development of CHD, the specific genetics remain unknown for the majority of patients. To elucidate the underlying genetic risk, we performed a genome wide association study (GWAS) of CHDs in general and specific CHD subgroups using the FinnGen Release 10 (R10) (N > 393,000), followed by functional fine-mapping through eQTL and co-localization analyses using the GTEx database.

RESULTS

We discovered three genome-wide significant loci associated with general CHD. Two of them were located in chromosome 17: 17q21.32 (rs2316327, intronic: LRRC37A2, Odds ratio (OR) [95% Confidence Interval (CI)] = 1.17[1.12-1.23], p = 1.5 × 10-9) and 17q25.3 (rs1293973611, nearest: BAHCC1, OR[95%CI] = 4.48[2.80-7.17], p = 7.0 × 10-10), respectively, and in addition to general CHD, the rs1293973611 locus was associated with the septal defect subtype. The third locus was in band 1p21.2 (rs35046143, nearest: PALMD, OR[95%CI] = 1.15[1.09-1.21], p = 7.1 × 10-9), and it was associated with general CHD and left-sided lesions. In the subgroup analysis, two additional loci were associated with septal defects (rs75230966 and rs6824295), and one with left-sided lesions (rs1305393195). In the eQTL analysis the variants rs2316327 (general CHD), and rs75230966 (septal defects) both located in 17q21.32 (with a LD r2 of 0.41) were both predicted to significantly associate with the expression of WNT9B in the atrial appendage tissue category. This effect was further confirmed by co-localization analysis, which also implicated WNT3 expression in the atrial appendage. A meta-analysis of general CHD together with the UK Biobank (combined N = 881,678) provided a different genome-wide significant locus in LRRC37A2; rs16941382 (OR[95%CI] = 1.15[1.11-1.20], p = 1.5 × 10-9) which is in significant LD with rs2316327.

CONCLUSIONS

Our results of general CHD and different CHD subcategories identified a complex risk locus on chromosome 17 near BAHCC1 and LRRC37A2, interacting with the genes WNT9B, WNT3 and MYL4, may constitute potential novel CHD risk associated loci, warranting future experimental tests to determine their role.

Illuminating Shared Genetic Associations Between Oesophageal Carcinoma and Pulmonary Carcinoma Risk

Background: Lung cancer and oesophageal cancer are prevalent malignancies with rising incidence and mortality worldwide. While some environmental and behavioural risk factors for these cancers are established, the contribution of genetic factors to their pathogenesis remains incompletely defined. This study aimed to interrogate the intricate genetic relationship between lung cancer and oesophageal cancer and their potential comorbidity. Methods: We utilised linkage disequilibrium score regression (LDSC) to analyse the genetic correlation between oesophageal carcinoma and lung carcinoma. We then employed several approaches, including pleiotropic analysis under the composite null hypothesis (PLACO), multi-marker analysis of genomic annotation (MAGMA), cis-expression quantitative trait loci (eQTL) analysis, and a pan-cancer assessment to identify pleiotropic loci and genes. Finally, we performed bidirectional Mendelian randomisation (MR) to evaluate the causal relationship between these malignancies. Results: LDSC revealed a significant genetic correlation between oesophageal carcinoma and lung carcinoma. Further analysis identified shared gene loci including PGBD1, ZNF323, and WNK1 using PLACO. MAGMA identified enriched pathways and 9 pleiotropic genes including HIST1H1B, HIST1H4L, and HIST1H2BL. eQTL analysis integrating oesophageal, lung, and blood tissues revealed 26 shared genes including TERT, NKAPL, RAD52, BTN3A2, GABBR1, CLPTM1L, and TRIM27. A pan-cancer exploration of the identified genes was undertaken. MR analysis showed no evidence for a bidirectional causal relationship between oesophageal carcinoma and lung carcinoma. Conclusions: This study provides salient insights into the intricate genetic links between lung carcinoma and oesophageal carcinoma. Utilising multiple approaches for genetic correlation, locus and gene analysis, and causal assessment, we identify shared genetic susceptibilities and regulatory mechanisms. These findings reveal new leads and targets to further elucidate the genetic basis of lung and oesophageal carcinoma, aiding development of preventive and therapeutic strategies.

Genome-wide association analyses using machine learning-based phenotyping reveal genetic architecture of occupational creativity and overlap with psychiatric disorders

Creativity is known to be heritable and exhibits familial aggregation with psychiatric disorders; however, the complex nature of their relationship has not been well-established. In the present study, we demonstrate that using an expanded and validated machine learning (ML)-based phenotyping of occupational creativity (OC) can allow us to further understand the trait of creativity, which was previously difficult to define and study. We conducted the largest genome-wide association study (GWAS) on OC with 241,736 participants from the UK Biobank and identified 25 lead variants that have not yet been reported and three candidate causal genes that were previously associated with educational attainment and psychiatric disorders. We found extensive genetic overlap between OC and psychiatric disorders with mixed effect direction through various post-GWAS analyses, including the bivariate causal mixture model. In addition, we discovered a strongly genetic correlation between our original GWAS and the GWAS adjusted for education years (rg = 0.95). Our GWAS analysis via ML-based phenotyping contributes to the understanding of the genetic architecture of creativity, which may inform genetic discovery and genetic prediction in human cognition and psychiatric disorders.

Causal relationships between coffee intake, apolipoprotein B and gastric, colorectal, and esophageal cancers: univariable and multivariable Mendelian randomization

PURPOSE

Coffee intake and apolipoprotein B levels have been linked to gastric, colorectal, and esophageal cancers in numerous recent studies. However, whether these associations are all causal remains unestablished. This study aimed to assess the potential causal associations of apolipoprotein B and coffee intake with the risk of gastric, colorectal, and esophageal cancers using Mendelian randomization analysis.

METHODS

In this study, we utilized a two-sample Mendelian randomization analysis to access the causal effects of coffee intake and apolipoprotein B on gastric, colorectal, and esophageal cancers. The summary statistics of coffee intake (n = 428,860) and apolipoprotein B (n = 439,214) were obtained from the UK Biobank. In addition, the summary statistics of gastric cancer, colorectal cancer, and esophageal cancer were obtained from the FinnGen biobank (n = 218,792). Inverse variance weighted, MR-Egger, weighted median, and weighted mode were applied to examine the causal relationship between coffee intake, apolipoprotein B and gastric, colorectal, and esophageal cancers. MR-Egger intercept test, Cochran's Q test, and leave-one-out analysis were performed to evaluate possible heterogeneity and pleiotropy. Steiger filtering and bidirectional mendelian randomization analysis were performed to evaluate the possible reverse causality.

RESULTS

The result of the inverse variance weighted method indicated that apolipoprotein B levels were significantly associated with a higher risk of gastric cancer (OR = 1.392, 95% CI 1.027-1.889, P = 0.0333) and colorectal cancer (OR = 1.188, 95% CI 1.001-1.411, P = 0.0491). Furthermore, multivariable Mendelian randomization analysis also revealed a positive association between apolipoprotein B levels and colorectal cancer risk, but the effect of apolipoprotein B on gastric cancer risk disappeared after adjustment of coffee intake, body mass index or lipid-related traits. However, we did not discover any conclusive evidence linking coffee intake to gastric, colorectal, or esophageal cancers.

CONCLUSIONS

This study suggested a causal association between genetically increased apolipoprotein B levels and higher risk of colorectal cancer. No causal relationship was observed between coffee intake and gastric, colorectal, or esophageal cancers.

Variant to gene mapping for carpal tunnel syndrome risk loci implicates skeletal muscle regulatory elements

BACKGROUND

Carpal tunnel syndrome (CTS) is a common disorder caused by compression of the median nerve in the wrist, resulting in pain and numbness throughout the hand and forearm. While multiple behavioural and physiological factors influence CTS risk, a growing body of evidence supports a strong genetic contribution. Recent genome-wide association study (GWAS) efforts have reported 53 independent signals associated with CTS. While GWAS can identify genetic loci conferring risk, it does not determine which cell types drive the genetic aetiology of the trait, which variants are "causal" at a given signal, and which effector genes correspond to these non-coding variants. These obstacles limit interpretation of potential disease mechanisms.

METHODS

We analysed CTS GWAS findings in the context of chromatin conformation between gene promoters and accessible chromatin regions across cellular models of bone, skeletal muscle, adipocytes and neurons. We identified proxy variants in high LD with the lead CTS sentinel SNPs residing in promoter connected open chromatin in the skeletal muscle and bone contexts.

FINDINGS

We detected significant enrichment for heritability in skeletal muscle myotubes, as well as a weaker correlation in human mesenchymal stem cell-derived osteoblasts. In myotubes, our approach implicated 117 genes contacting 60 proxy variants corresponding to 20 of the 53 GWAS signals. In the osteoblast context we implicated 30 genes contacting 24 proxy variants coinciding with 12 signals, of which 19 genes shared. We subsequently prioritized BZW2 as a candidate effector gene in CTS and implicated it as novel gene that perturbs myocyte differentiation in vitro.

INTERPRETATION

Taken together our results suggest that the CTS genetic component influences the size, integrity, and organization of multiple tissues surrounding the carpal tunnel, in particular muscle and bone, to predispose the nerve to being compressed in this disease setting.

FUNDING

This work was supported by NIH Grant UM1 DK126194 (SFAG and WY), R01AG072705 (SFAG & KDH) and the Center for Spatial and Functional Genomics at CHOP (SFAG & ADW). SFAG is supported by the Daniel B. Burke Endowed Chair for Diabetes Research. WY is supported by the Perelman School of Medicine of the University of Pennsylvania.

Screen time in the development of cardiovascular diseases: A two-sample Mendelian randomization study

BACKGROUND AND AIMS

Coronary artery disease (CAD), heart failure (HF), and ischemic heart disease (IHD) are three common cardiovascular diseases that are closely associated with metabolic activity. The global incidence and prevalence of these conditions are on the rise, primarily due to unhealthy lifestyles, aging populations, and the increasing prevalence of obesity and diabetes. Excessive screen time has emerged as a potential risk factor for various adverse health outcomes, although limited research has explored its relationship with cardiovascular disease outcomes.

METHODS AND RESULTS

A Mendelian randomization (MR) study was conducted, employing exposure-associated genetic variants as instrumental variables to explore the causal relationship between screen time use and cardiovascular disease outcomes. Single nucleotide polymorphisms (SNPs) were utilized as pooled data for the genetic variable instrument, investigating the association between screen use duration and three types of cardiovascular diseases: coronary artery disease (CAD), heart failure (HF), and ischemic heart disease (IHD). Through the MR analysis, it was revealed that the use of mobile phones and TV screens exhibited a significant causal association with the occurrence of CAD, heart failure, and IHD. However, no significant association was observed between the use of computers and these three types of cardiovascular diseases.

CONCLUSION

Our study suggests that excessive screen time use is associated with the development of cardiovascular disease. However, it should be noted that the consequences of screen time can vary depending on the reasons and purposes for its use. Implementing reasonable control over screen time, particularly for entertainment purposes, holds promise as a potential approach to mitigating cardiovascular disease.

Intervertebral Disk Degeneration and Bone Mineral Density: A Bidirectional Mendelian Randomization Study

This study aimed to investigate the causal relationship between bone mineral density (BMD) and intervertebral disk degeneration (IVDD) using a two-sample bidirectional Mendelian randomization analysis. Summary-level data from the Genome-Wide Association Study (GWAS) were used. Instrumental variables (IVs) for IVDD were selected from the large-scale Genome-Wide Association Study (GWAS) (20,001 cases and 164,682 controls). Bone mineral density (BMD) at five different sites (heel (n = 426,824), total body (TB) (n = 56,284), forearm (FA) (n = 8143), femoral neck (FN) (n = 32,735), and lumbar spine (LS) (n = 28,498)) was used as a phenotype for OP. Bidirectional causality between IVDD and BMD was assessed using inverse variance weighting (IVW) and other methods. Related sensitivity analyses were performed. Myopia was also analyzed as a negative control result to ensure the validity of IVs. Heel bone mineral density (heel BMD), total body bone mineral density (TB-BMD), femoral neck bone mineral density (FN-BMD), and lumbar spine bone mineral density (LS-BMD) have a direct causal relationship on intervertebral disk degeneration (IVDD) [heel BMD-related analysis: beta = 0.06, p = 0.03; TB-BMD-related analysis: beta = 0.18, p = 8.72E-08; FN-BMD-related analysis: beta = 0.15, p = 4.89E-03; LS-BMD-related analysis: beta = 0.16, p = 1.43E-04]. There was no evidence of a significant causal effect of IVDD on BMD. In conclusion, our study found a significant positive causal effect of lower BMD on IVDD, and we identified significant causal effects of heel, TB-, FN-, and LS-BMD on IVDD, but there was no evidence of a significant causal effect of IVDD on BMD.

Genome-wide association study identifies genetic regulation of oestrone concentrations and association with endometrial cancer risk in postmenopausal women

BACKGROUND

Oestrone, predominantly made in fat, is the main circulating oestrogen and important for target tissue oestradiol production in women after menopause. The present study was undertaken to determine the genetic regulation of blood oestrone, measured with precision, in postmenopausal women and to explore associations between the identified genetic loci and endometrial cancer in a large, independent cohort.

METHODS

A genome-wide association study (GWAS) was undertaken in women aged at least 70 years to identify genetic associations with blood oestrone concentrations measured by liquid chromatography and tandem mass spectrometry. The GWAS included participants from the Sex Hormones in Older Women (SHOW) study, a sub-study of the longitudinal ASPREE (ASPirin in Reducing Events in the Elderly) randomised trial. Of the 6358 women providing a biobank sample at enrolment, 4951 unrelated women of European ancestry, not taking sex hormones, anti-oestrogens, anti-androgens or systemic glucocorticoids were included in the GWAS. Single nucleotide polymorphisms (SNPs) from loci identified below the genome-wide significance threshold were then tested in an independent cohort (the UK Biobank) for association with endometrial cancer risk, using logistic regression and adjusting for age, body mass index (BMI) and the top 10 genetic principal components.

FINDINGS

The median age of the 4951 women included in the GWAS was 75.9 years (range 70-94.8 years). The GWAS identified four independent SNPs associated with oestrone concentrations (p < 5 × 10-8). Among them, the effect (minor) alleles rs34670419-T, rs2846729-T and rs2414098-T were associated with lower oestrone concentrations. Carrying these effect alleles was associated with lower oestrone concentrations in a dose-dependent manner. The effect allele rs56400819-A was associated with higher oestrone concentrations. When applied to UK Biobank, carrier status for rs2414098-T associated with the CYP19A1 gene which encodes the aromatase enzyme required for oestrogen synthesis was significantly associated with lower endometrial cancer risk (adjusted odd ratio [aOR] 0.87 [95% CI 0.82-0.93]; p = 6.69 × 10-5 for women across all ages and aOR 0.89 [95% CI 0.83-0.96]; p = 0.003 for postmenopausal women). None of the models that included age, body mass index (BMI), the top 10 genetic principal components, parity and diabetes mellitus explained more than 7.6% of the variation in risk.

INTERPRETATION

We have shown genetic regulation of oestrone concentrations in postmenopausal women, and that SNPs associated with oestrone were also associated with endometrial cancer risk, independent of BMI, parity and diabetes mellitus. Although the apparent contribution was modest, the biological influence of oestrone concentrations may be greater through conversion to oestradiol in endometrial tissue.

FUNDING

The ASPREE trial was supported by the National Institute on Aging and the National Cancer Institute at the National Institutes of Health (Grant U01AG029824); the National Health and Medical Research Council (NHMRC) of Australia (Grant 34047, 1127060); Monash University (Australia); and the Victorian Cancer Agency (Australia). The ASPREE Healthy Ageing Biobank was funded by the CSIRO (Flagship Grant), the National Cancer Institute (Grant U01 AG029824) and Monash University. This analysis of sex hormones was funded by an NHMRC of Australia Project Grant (No. 1105305). SRD holds an NHMRC Investigator Grant (2016627). PL is supported by a National Heart Foundation Future Leader Fellowship (102604).

Phenotypic and genome-wide studies on dicarbonyls: major associations to glomerular filtration rate and gamma-glutamyltransferase activity

BACKGROUND

The dicarbonyl compounds methylglyoxal (MG), glyoxal (GO) and 3-deoxyglucosone (3-DG) have been linked to various diseases. However, disease-independent phenotypic and genotypic association studies with phenome-wide and genome-wide reach, respectively, have not been provided.

METHODS

MG, GO and 3-DG were measured by LC-MS in 1304 serum samples of two populations (KORA, n = 482; BiDirect, n = 822) and assessed for associations with genome-wide SNPs (GWAS) and with phenome-wide traits. Redundancy analysis (RDA) was used to identify major independent trait associations.

FINDINGS

Mutual correlations of dicarbonyls were highly significant, being stronger between MG and GO (ρ = 0.6) than between 3-DG and MG or GO (ρ = 0.4). Significant phenotypic results included associations of all dicarbonyls with sex, waist-to-hip ratio, glomerular filtration rate (GFR), gamma-glutamyltransferase (GGT), and hypertension, of MG and GO with age and C-reactive protein, of GO and 3-DG with glucose and antidiabetics, of MG with contraceptives, of GO with ferritin, and of 3-DG with smoking. RDA revealed GFR, GGT and, in case of 3-DG, glucose as major contributors to dicarbonyl variance. GWAS did not identify genome-wide significant loci. SNPs previously associated with glyoxalase activity did not reach nominal significance. When multiple testing was restricted to the lead SNPs of GWASs on the traits selected by RDA, 3-DG was found to be associated (p = 2.3 × 10-5) with rs1741177, an eQTL of NF-κB inhibitor NFKBIA.

INTERPRETATION

This large-scale, population-based study has identified numerous associations, with GFR and GGT being of pivotal importance, providing unbiased perspectives on dicarbonyls beyond the current state.

FUNDING

Deutsche Forschungsgemeinschaft, Helmholtz Munich, German Centre for Cardiovascular Research (DZHK), German Federal Ministry of Research and Education (BMBF).

RETROSPECTIVE VARYING COEFFICIENT ASSOCIATION ANALYSIS OF LONGITUDINAL BINARY TRAITS: APPLICATION TO THE IDENTIFICATION OF GENETIC LOCI ASSOCIATED WITH HYPERTENSION

Many genetic studies contain rich information on longitudinal phenotypes that require powerful analytical tools for optimal analysis. Genetic analysis of longitudinal data that incorporates temporal variation is important for understanding the genetic architecture and biological variation of complex diseases. Most of the existing methods assume that the contribution of genetic variants is constant over time and fail to capture the dynamic pattern of disease progression. However, the relative influence of genetic variants on complex traits fluctuates over time. In this study, we propose a retrospective varying coefficient mixed model association test, RVMMAT, to detect time-varying genetic effect on longitudinal binary traits. We model dynamic genetic effect using smoothing splines, estimate model parameters by maximizing a double penalized quasi-likelihood function, design a joint test using a Cauchy combination method, and evaluate statistical significance via a retrospective approach to achieve robustness to model misspecification. Through simulations we illustrated that the retrospective varying-coefficient test was robust to model misspecification under different ascertainment schemes and gained power over the association methods assuming constant genetic effect. We applied RVMMAT to a genome-wide association analysis of longitudinal measure of hypertension in the Multi-Ethnic Study of Atherosclerosis. Pathway analysis identified two important pathways related to G-protein signaling and DNA damage. Our results demonstrated that RVMMAT could detect biologically relevant loci and pathways in a genome scan and provided insight into the genetic architecture of hypertension.

Genome Wide Association Study of Neuropathic Ocular Pain

Purpose

To conduct a genome-wide association study (GWAS) of individuals with neuropathic ocular pain (NOP) symptoms to identify genomic variants that may predispose to NOP development.

Design

Prospective study of individuals with NOP.

Participants

Three hundred twenty-nine patients recruited from the Miami Veterans Affairs eye clinic.

Methods

The Neuropathic Pain Symptom Inventory modified for the eye (NPSI-Eye) was completed to calculate a NPSI-Eye-Sub-Score (summed ratings of burning and wind sensitivity) as an indicator of NOP severity. A GWAS was performed for the NPSI-Eye-Sub-Score with a significance threshold of P < 5 × 10-8. A gene-based analysis was performed using the multimarker analysis of genomic annotation software (in the functional mapping and annotation of GWAS online platform). The 13 865 778 single nucleotide polymorphisms (SNPs) from our GWAS analysis were mapped to 10 834 protein coding genes, and significant genes were run through gene set enrichment analysis.

Main Outcome Measures

Identification of SNPs and protein products that may be associated with the development of NOP.

Results

One hundred seventy-one SNPs reached a threshold of P < 10-5, of which 10 SNPs reached the suggestive level of significance of P < 5 × 10-7 and 1 SNP met our genome-wide significance threshold of P < 5 × 10-8. This lead SNP, rs140293404 (P = 1.23 × 10-8), is an intronic variant found within gene ENSG00000287251 coding for transcript ENST00000662732.1. Rs140293404 is in linkage disequilibrium with exon variant rs7926353 (r2 > 0.8) within ENSG00000279046 coding for transcript ENST00000624288.1. The most significant genes from gene-based tests were matrix metalloproteinase-19 (MMP19) (P = 1.12 × 10-5), zinc finger RNA-binding motif and serine/arginine rich-1 (ZRSR1) (P = 1.48 × 10-4), CTC-487M23.8 (P = 1.79 × 10-4), receptor expression-enhancing protein-5 (REEP5) (P = 2.36 × 10-4), and signal recognition particle-19 (SRP19) (P = 2.56 × 10-4). From gene set enrichment analysis, the sensory perception (false discovery rate = 6.57 × 10-3) and olfactory signaling (false discovery rate = 1.63 × 10-2) pathways were enriched with the most significant genes.

Conclusions

Our GWAS revealed genes with protein products that may impact sensory perception, lending biological plausibility to a role for SNPs identified by our GWAS in the development of NOP. A better understanding of the biological relevance of these genes and pathways in the pathophysiology associated with NOP may facilitate future novel mechanism-based treatments.

Financial Disclosures

Proprietary or commercial disclosure may be found in the Footnotes and Disclosures at the end of this article.

Genome-wide association studies reveal stable loci for wheat grain size under different sowing dates

Background

Wheat (Tritium aestivum L.) production is critical for global food security. In recent years, due to climate change and the prolonged growing period of rice varieties, the delayed sowing of wheat has resulted in a loss of grain yield in the area of the middle and lower reaches of the Yangtze River. It is of great significance to screen for natural germplasm resources of wheat that are resistant to late sowing and to explore genetic loci that stably control grain size and yield.

Methods

A collection of 327 wheat accessions from diverse sources were subjected to genome-wide association studies using genotyping-by-sequencing. Field trials were conducted under normal, delayed, and seriously delayed sowing conditions for grain length, width, and thousand-grain weight at two sites. Additionally, the additive main effects and multiplicative interaction (AMMI) model was applied to evaluate the stability of thousand-grain weight of 327 accessions across multiple sowing dates.

Results

Four wheat germplasm resources have been screened, demonstrating higher stability of thousand-grain weight. A total of 43, 35, and 39 significant MTAs were determined across all chromosomes except for 4D under the three sowing dates, respectively. A total of 10.31% of MTAs that stably affect wheat grain size could be repeatedly identified in at least two sowing dates, with PVE ranging from 0.03% to 38.06%. Among these, six were for GL, three for GW, and one for TGW. There were three novel and stable loci (4A_598189950, 4B_307707920, 2D_622241054) located in conserved regions of the genome, which provide excellent genetic resources for pyramid breeding strategies of superior loci. Our findings offer a theoretical basis for cultivar improvement and marker-assisted selection in wheat breeding practices.

Genetic control of abiotic stress-related specialized metabolites in sunflower

BACKGROUND

Abiotic stresses in plants include all the environmental conditions that significantly reduce yields, like drought and heat. One of the most significant effects they exert at the cellular level is the accumulation of reactive oxygen species, which cause extensive damage. Plants possess two mechanisms to counter these molecules, i.e. detoxifying enzymes and non-enzymatic antioxidants, which include many classes of specialized metabolites. Sunflower, the fourth global oilseed, is considered moderately drought resistant. Abiotic stress tolerance in this crop has been studied using many approaches, but the control of specialized metabolites in this context remains poorly understood. Here, we performed the first genome-wide association study using abiotic stress-related specialized metabolites as molecular phenotypes in sunflower. After analyzing leaf specialized metabolites of 450 hybrids using liquid chromatography-mass spectrometry, we selected a subset of these compounds based on their association with previously known abiotic stress-related quantitative trait loci. Eventually, we characterized these molecules and their associated genes.

RESULTS

We putatively annotated 30 compounds which co-localized with abiotic stress-related quantitative trait loci and which were associated to seven most likely candidate genes. A large proportion of these compounds were potential antioxidants, which was in agreement with the role of specialized metabolites in abiotic stresses. The seven associated most likely candidate genes, instead, mainly belonged to cytochromes P450 and glycosyltransferases, two large superfamilies which catalyze greatly diverse reactions and create a wide variety of chemical modifications. This was consistent with the high plasticity of specialized metabolism in plants.

CONCLUSIONS

This is the first characterization of the genetic control of abiotic stress-related specialized metabolites in sunflower. By providing hints concerning the importance of antioxidant molecules in this biological context, and by highlighting some of the potential molecular mechanisms underlying their biosynthesis, it could pave the way for novel applications in breeding. Although further analyses will be required to better understand this topic, studying how antioxidants contribute to the tolerance to abiotic stresses in sunflower appears as a promising area of research.

Psoriasis and risk of central neurological disorders in European populations: A mendelian randomization analysis

Background

People with psoriasis are at a higher risk for having central neurological problems, according to previous studies; however, it is unclear if there is a genetic link between the risk of developing psoriasis and developing central neurological disorders. In this study, the possible link between genetically predisposed psoriasis and the risk of common central nervous system disorders was comprehensively investigated.

Methods

There was no overlap in the participant populations between the psoriasis and central neurological disorders genome-wide association studies, which provide the genetic resources. Inverse variance weighting, often used as Mendelian randomization (MR) analysis, is the main method. To guarantee the accuracy of our findings, a number of sensitivity studies were carried out.

Results

MR analysis revealed that although psoriasis was reported to increase the risk of Parkinson's disease (OR = 4.42, 95%CI[-3.81～6.79], P = 0.58) and epilepsy (OR = 4.71, 95%CI[-2.20～5.30], P = 0.42) in this study, they did not reach statistical significance. At the same time, this study did not observe that psoriasis would increase the risk of multiple sclerosis (OR = 0.01, 95%CI [-12.61～3.83], P = 0.30) and migraine (OR = 0.99, 95%CI [0.94～1.05], P = 0.78), they also did not reach statistical significance. Under all sensitivity assessments, the results remained stable.

Conclusions

Psoriasis does not appear to raise the risk of migraine, Parkinson's disease, multiple sclerosis, or epilepsy, according to our study.

Plasma campesterol and ABCG5/ABCG8 gene loci on the risk of cholelithiasis and cholecystitis: evidence from Mendelian randomization and colocalization analyses

The causal relationships between plasma metabolites and cholelithiasis/cholecystitis risks remain elusive. Using two-sample Mendelian randomization, we found that genetic proxied plasma campesterol level showed negative correlation with the risk of both cholelithiasis and cholecystitis. Furthermore, the increased risk of cholelithiasis is correlating with the increased level of plasma campesterol. Lastly, genetic colocalization study showed that the leading SNP, rs4299376, which residing at the ABCG5/ABCG8 gene loci, was shared by plasma campesterol level and cholelithiasis, indicating that the aberrant transportation of plant sterol/cholesterol from the blood stream to the bile duct/gut lumen might be the key in preventing cholesterol gallstone formation.

Evaluation of penalized and machine learning methods for asthma disease prediction in the Korean Genome and Epidemiology Study (KoGES)

BACKGROUND

Genome-wide association studies have successfully identified genetic variants associated with human disease. Various statistical approaches based on penalized and machine learning methods have recently been proposed for disease prediction. In this study, we evaluated the performance of several such methods for predicting asthma using the Korean Chip (KORV1.1) from the Korean Genome and Epidemiology Study (KoGES).

RESULTS

First, single-nucleotide polymorphisms were selected via single-variant tests using logistic regression with the adjustment of several epidemiological factors. Next, we evaluated the following methods for disease prediction: ridge, least absolute shrinkage and selection operator, elastic net, smoothly clipped absolute deviation, support vector machine, random forest, boosting, bagging, naïve Bayes, and k-nearest neighbor. Finally, we compared their predictive performance based on the area under the curve of the receiver operating characteristic curves, precision, recall, F1-score, Cohen's Kappa, balanced accuracy, error rate, Matthews correlation coefficient, and area under the precision-recall curve. Additionally, three oversampling algorithms are used to deal with imbalance problems.

CONCLUSIONS

Our results show that penalized methods exhibit better predictive performance for asthma than that achieved via machine learning methods. On the other hand, in the oversampling study, randomforest and boosting methods overall showed better prediction performance than penalized methods.

The causality between systemic inflammatory regulators and chronic respiratory diseases: A bidirectional Mendelian-randomization study

INTRODUCTION

Accumulative evidence suggests the associations between systemic inflammatory regulators and chronic respiratory diseases (CRDs). However, the intrinsic causation remains implicit. Therefore, this study aimed to examine causative associations by mendelian randomization (MR) and to identify valuable active factors.

METHODS

Based on data from the GWAS database, we performed MR analyses of 41 serum cytokines from 8,293 Finnish and European descent cohorts from GBMI and UKBB for five major CRDs. We mainly applied inverse variance weighted regression, supplemented by MR-Egger regression, weighted median, maximum likelihood, weighted mode, and simple mode algorithms. Moreover, sensitivity analyses were conducted using Cochrane's Q test, MR-Egger intercept, MR-PRESSO Global test and MR-Steiger filtering. Eventually, the consistency of MR results was assessed by leave-one-out.

RESULTS

Our results suggest that 12 genetically predicted systemic inflammatory regulators probably participate in the progression of CRDs, including four risk factors (IL-1RA, IL-4, MIP-1A, PDGF-BB) and one protective factor (IL-6) in IPF, two protective factors (SCF, SDF-1A) in COPD, and two protective factors (SCF, SDF-1A) in asthma, two protective factors (GROA, IL-2RA) were also included in asthma, whereas only one factor (HGF) was protective against bronchiectasis. Additionally, two protective factors (FGF-BASIC, G-CSF) were identified in sarcoidosis. Sensitivity analyses showed no horizontal pleiotropy and significant heterogeneity. Finally, based on the findings of inverse MR analysis, no inverse causal association was uncovered, confirming the robustness of results.

CONCLUSION

Our study unearths potential associations between systemic inflammatory modulators and common CRDs, providing new insights for inflammation-mediated CRD prevention and therapeutic approaches.

Lung function and the risk of frailty in the European population: a mendelian randomization study

BACKGROUND

Epidemiological evidence has suggested a relationship between lung function and frailty, but the precise nature of the causality remains unclear. In this study, we applied a two-sample Mendelian randomization (MR) analysis to determine the causal effects of lung function on frailty.

METHODS

Single nucleotide polymorphisms (SNPs) independently related (P ≤ 5E-08) to lung function, as identified by genome-wide association study (GWAS), were applied as instrumental variables (IV). The association with frailty index (FI) was investigated using summary-level data from the latest GWAS on FI (n = 175,226). Different statistical methods were employed to evaluate the causal estimates between lung function and FI. The pleiotropy, heterogeneity, and leave-one-out analysis were applied to confirm the stability of the MR estimates.

RESULTS

Using the random-effect inverse-variance weighted approach, genetically proxied forced expiratory volume in the first second (FEV1), ratio of FEV1 on forced vital capacity (FVC) [FEV1/FVC], and peak expiratory flow (PEF) were significantly and inversely associated with FI (FEV1, β = -0.08, P = 2.03E-05; FEV1/FVC, β = -0.06, P = 9.51E-06; PEF, β = -0.07, P = 4.09E-04) with good statistical power (99.7-100%). However, no significant association was observed between FVC and FI (β = -0.01, P = 0.681). Leave-one-out analysis showed that there was no single SNP driving the bias of the estimates. There was potential heterogeneity, but no obvious pleiotropy was founded in this MR study.

CONCLUSIONS

Our findings indicate that impaired pulmonary function is closely related to the risk of frailty. Enhancing lung function in the elderly population may contribute to the prevention of frailty to a certain extent.

Identifying risk loci for FTD and shared genetic component with ALS: A large-scale multitrait association analysis

Current genome-wide association studies of frontotemporal dementia (FTD) are underpowered due to limited samples. Further, common genetic etiologies between FTD and amyotrophic lateral sclerosis (ALS) remain unknown. Using the largest summary statistics of FTD (3526 cases and 9402 controls) and ALS (27,205 cases and 110,881 controls), we found a significant genetic correlation between them (rˆg = 0.637, P = 0.032) and identified 190 FTD-related variants within 5 loci (3p22.1, 5q35.1, 9p21.2, 19p13.11, and 20q13.13). Among these, ALS and FTD had causal variants in 9p21.2 and 19p13.11. Moreover, MOBP (3p22.1), C9orf72 (9p21.2), MOB3B (9p21.2), UNC13A (19p13.11), SLC9A8 (20q13.13), SNAI1 (20q13.13), and SPATA2 (20q13.13) were discovered by both SNP- and gene-level analyses, which together discovered 15 FTD-associated genes, with 10 not detected before (IFNK, RNF114, SLC9A8, SPATA2, SNAI1, SCFD1, POLDIP2, TMEM97, G2E3, and PIGW). Functional analyses showed these genes were enriched in heart left ventricle, kidney cortex, and some brain regions. Overall, this study provides insights into genetic determinants of FTD and shared genetic etiology underlying FTD and ALS.

Genetics of Parkinson's disease heterogeneity: A genome-wide association study of clinical subtypes

INTRODUCTION

Substantial heterogeneity between individual patients in the clinical presentation of Parkinson's disease (PD) has led to the classification of distinct PD subtypes. However, genetic susceptibility factors for specific PD subtypes are not well understood. Therefore, the present study aimed to investigate the genetics of PD heterogeneity by performing a genome-wide association study (GWAS) of PD subtypes.

METHODS

A total of 799 PD patients were included and classified into tremor-dominant (TD) (N = 345), akinetic-rigid (AR) (N = 227), gait-difficulty (GD) (N = 82), and mixed (MX) (N = 145) phenotypic subtypes. After array genotyping and subsequent imputation, a total of 7,918,344 variants were assessed for association with each PD subtype using logistic regression models that were adjusted for age, sex, and the top five principal components of GWAS data.

RESULTS

We identified one genome-wide significant association (P < 5 × 10-8), which was between the MIR3976HG rs7504760 variant and the AR subtype (Odds ratio [OR] = 6.12, P = 2.57 × 10-8). Suggestive associations (P < 1 × 10-6) were observed regarding TD for RP11-497G19.3/RP11-497G19.1 rs7304254 (OR = 3.33, P = 3.89 × 10-7), regarding GD for HES2 rs111473931 (OR = 3.18, P = 6.85 × 10-7), RP11-400D2.3/CTD-2012I17.1 rs149082205 (OR = 8.96, P = 9.08 × 10-7), and RN7SL408P/SGK1 rs56161738 (OR = 2.97, P = 6.19 × 10-7), and regarding MX for MMRN2 rs112991171 (OR = 4.98, P = 1.02 × 10-7).

CONCLUSION

Our findings indicate that genetic variation may account for part of the clinical heterogeneity of PD. In particular, we found a novel genome-wide significant association between MIR3976HG variation and the AR PD subtype. Replication of these findings will be important in order to better define the genetic architecture of clinical variability in PD disease presentation.

Atopic dermatitis-associated genetic variants regulate LOC100294145 expression implicating interleukin-27 production and type 1 interferon signaling

Background

Atopic dermatitis (AD) is a complex inflammatory disease with a strong genetic component. A singular approach of genome wide association studies (GWAS) can identify AD-associated genetic variants, but is unable to explain their functional relevance in AD. This study aims to characterize AD-associated genetic variants and elucidate the mechanisms leading to AD through a multi-omics approach.

Methods

GWAS identified an association between genetic variants at 6p21.32 locus and AD. Genotypes of 6p21.32 locus variants were evaluated against LOC100294145 expression in peripheral blood mononuclear cells (PBMCs). Their influence on LOC100294145 promoter activity was measured in vitro via a dual-luciferase assay. The function of LOC100294145 was then elucidated through a combination of co-expression analyses and gene enrichment with g:Profiler. Mendelian randomization was further used to assess the causal regulatory effect of LOC100294145 on its co-expressed genes.

Results

Minor alleles of rs116160149 and rs115388857 at 6p21.32 locus were associated with increased AD risk (p = 2.175 × 10-8, OR = 1.552; p = 2.805 × 10-9, OR = 1.55) and higher LOC100294145 expression in PBMCs (adjusted p = 0.182; 8.267 × 10-12). LOC100294145 expression was also found to be increased in those with AD (adjusted p = 3.653 × 10-2). The genotype effect of 6p21.32 locus on LOC100294145 promoter activity was further validated in vitro. Co-expression analyses predicted LOC100294145 protein's involvement in interleukin-27 and type 1 interferon signaling, which was further substantiated through mendelian randomization.

Conclusion

Genetic variants at 6p21.32 locus increase AD susceptibility through raising LOC100294145 expression. A multi-omics approach enabled the deduction of its pathogenesis model comprising dysregulation of hub genes involved in type 1 interferon and interleukin 27 signaling.

Increased expression of TNFRSF14 and LIGHT in biliary epithelial cells of patients with primary sclerosing cholangitis

BACKGROUND AND AIMS

There is a lack of biliary epithelial molecular markers for primary sclerosing cholangitis (PSC). We analyzed candidates from disease susceptibility genes identified in recent genome-wide association studies (GWAS).

METHODS

Expression levels of GWAS genes were analyzed in archival liver tissues of patients with PSC and controls. Immunohistochemical analysis was performed to evaluate expression levels in the biliary epithelia of PSC (N = 45) and controls (N = 12). Samples from patients with primary biliary cholangitis (PBC) were used as disease controls (N = 20).

RESULTS

Hepatic expression levels of ATXN2, HHEX, PRDX5, MST1, and TNFRSF14 were significantly altered in the PSC group. We focused on the immune-related receptor, TNFRSF14. Immunohistochemistry revealed that high expression of TNFRSF14 in biliary epithelial cells was observed only in the PSC group. In addition, the expression of LIGHT, which encodes a TNFRSF14-activating ligand, was increased in PSC liver. Immunohistochemistry showed that high expression of LIGHT was more common in PSC biliary epithelia (53%) than in the PBC (15%) or control (0%) groups; moreover, it was positively associated with fibrotic progression, although it was not an independent prognostic factor.

CONCLUSIONS

TNFRSF14 and LIGHT are promising candidate markers for PSC.

Predictive Performance of an Updated Polygenic Risk Score for Age-Related Macular Degeneration

PURPOSE

A recent genome-wide association study of age-related macular degeneration (AMD) identified new AMD-associated risk variants. These variants now can be incorporated into an updated polygenic risk score (PRS). This study aimed to assess the performance of an updated PRS, PRS2023, in an independent cohort of older individuals with retinal imaging data and to compare performance with an older PRS, PRS2016.

DESIGN

Cross-sectional study.

PARTICIPANTS

A total of 4175 participants of European ancestry, 70 years of age or older, with genotype and retinal imaging data.

METHODS

We used logistic regression models and area under the receiver operating characteristic curve (AUC) to assess the performance of PRS2023 compared with PRS2016. AMD status and severity were graded using color fundus photography.

MAIN OUTCOME MEASURES

Association of PRS2023 and PRS2016 with AMD risk at baseline.

RESULTS

At enrollment among 4175 participants, 2605 participants (62.4%) had no AMD and 853 participants (20.4%), 671 participants (16.1%), and 46 participants (1.1%) had early, intermediate, and late-stage AMD, respectively. More than 27% of the participants with a high PRS2023 (top quartile) had intermediate or late-stage AMD, compared with < 15% for those in the middle 2 quartiles and less than 13% for those in the lowest quartile. Both PRS2023 and PRS2016 were associated significantly with AMD after adjustment for age, sex, smoking status, and lipid levels, with increasing odds ratios (ORs) for worsening AMD grades. PRS2023 outperformed PRS2016 (P = 0.03 for all AMD and P = 0.03 for late AMD, DeLong test comparing AUC). PRS2023 was associated with late-stage AMD with an adjusted OR of 5.05 (95% confidence interval [CI], 3.41-7.47) per standard deviation. The AUC of a model containing conventional or nongenetic risk factors and PRS2023 was 91% (95% CI, 87%-95%) for predicting late-stage AMD, which improved 12% over the model without the PRS (AUC, 79%; P < 0.001 for difference).

CONCLUSIONS

A new PRS, PRS2023, for AMD outperforms a previous PRS and predicts increasing risk for late-stage AMD (with stronger association for more severe imaging-confirmed AMD grades). Our findings have clinical implications for the improved prediction and risk stratification of AMD.

FINANCIAL DISCLOSURE(S)

Proprietary or commercial disclosure may be found in the Footnotes and Disclosures at the end of this article.

Genome-wide association study suggests genetic candidate loci of insulin dysregulation in Finnhorses

Equine metabolic syndrome (EMS) is a common welfare problem in horses worldwide. It is characterized by insulin dysregulation (ID), predisposition to laminitis and often obesity. EMS is multifactorial by nature, with both the environment and genetics contributing to the phenotype. Environmental factors, such as feeding and exercise, can be controlled, thus forming the basis for treatment and prevention. Genetic factors, by contrast, are less well-known and not easily controllable. The aim of this study was to identify potential genetic loci influencing ID/EMS in Finnhorses. A single-breed (Finnhorse) case-control genome-wide association study (GWAS) of ID was conducted with controls that included age-appropriate non-ID horses. ID status was determined with an oral sugar test (OST) for fasted horses. Seventy-one Finnhorses participated (n = 34 ID, n = 37 control). DNA samples (hair roots) were genotyped for 65 157 single-nucleotide polymorphisms (SNPs) with the Illumina Equine SNP70 BeadChip, and these data were analysed for association and FST outliers with genomic tools. P-values that exceeded the suggestive threshold (P = 1.00 ×10-5) were found in SNP BIEC2_383954 (P = 3.45 ×10-6) in chromosome 17 and SNP BIEC2_312374 (P = 1.89 ×10-5) in chromosome 15. Hierarchical and Bayesian FST outlier tests also detected these SNPs. Potential candidate genes associated with the ID close to SNP BIEC2_383954, with functions in carbohydrate metabolism, were Arginine and Glutamate Rich 1 (ARGLU1) and Ephrin-B2 (EFNB2).