Genomic Relationships, Novel Loci, and Pleiotropic Mechanisms across Eight Psychiatric Disorders.

Genetic influences on psychiatric disorders transcend diagnostic boundaries, suggesting substantial pleiotropy of contributing loci. However, the nature and mechanisms of these pleiotropic effects remain unclear. We performed analyses of 232,964 cases and 494,162 controls from genome-wide studies of anorexia nervosa, attention-deficit/hyperactivity disorder, autism spectrum disorder, bipolar disorder, major depression, obsessive-compulsive disorder, schizophrenia, and Tourette syndrome. Genetic correlation analyses revealed a meaningful structure within the eight disorders, identifying three groups of inter-related disorders. Meta-analysis across these eight disorders detected 109 loci associated with at least two psychiatric disorders, including 23 loci with pleiotropic effects on four or more disorders and 11 loci with antagonistic effects on multiple disorders. The pleiotropic loci are located within genes that show heightened expression in the brain throughout the lifespan, beginning prenatally in the second trimester, and play prominent roles in neurodevelopmental processes. These findings have important implications for psychiatric nosology, drug development, and risk prediction.

shaPRS: Leveraging shared genetic effects across traits or ancestries improves accuracy of polygenic scores.

We present shaPRS, a method that leverages widespread pleiotropy between traits or shared genetic effects across ancestries, to improve the accuracy of polygenic scores. The method uses genome-wide summary statistics from two diseases or ancestries to improve the genetic effect estimate and standard error at SNPs where there is homogeneity of effect between the two datasets. When there is significant evidence of heterogeneity, the genetic effect from the disease or population closest to the target population is maintained. We show via simulation and a series of real-world examples that shaPRS substantially enhances the accuracy of polygenic risk scores (PRSs) for complex diseases and greatly improves PRS performance across ancestries. shaPRS is a PRS pre-processing method that is agnostic to the actual PRS generation method, and as a result, it can be integrated into existing PRS generation pipelines and continue to be applied as more performant PRS methods are developed over time.

Novel insight into the etiology of ischemic stroke gained by integrative multiome-wide association study.

Stroke, characterized by sudden neurological deficits, is the second leading cause of death worldwide. Although genome-wide association studies (GWAS) have successfully identified many genomic regions associated with ischemic stroke (IS), the genes underlying risk and their regulatory mechanisms remain elusive. Here, we integrate a large-scale GWAS (N = 1 296 908) for IS together with molecular QTLs data, including mRNA, splicing, enhancer RNA (eRNA), and protein expression data from up to 50 tissues (total N = 11 588). We identify 136 genes/eRNA/proteins associated with IS risk across 60 independent genomic regions and find IS risk is most enriched for eQTLs in arterial and brain-related tissues. Focusing on IS-relevant tissues, we prioritize 9 genes/proteins using probabilistic fine-mapping TWAS analyses. In addition, we discover that blood cell traits, particularly reticulocyte cells, have shared genetic contributions with IS using TWAS-based pheWAS and genetic correlation analysis. Lastly, we integrate our findings with a large-scale pharmacological database and identify a secondary bile acid, deoxycholic acid, as a potential therapeutic component. Our work highlights IS risk genes/splicing-sites/enhancer activity/proteins with their phenotypic consequences using relevant tissues as well as identify potential therapeutic candidates for IS.

Blood-based genome-wide DNA methylation correlations across body-fat- and adiposity-related biochemical traits.

The recent increase in obesity levels across many countries is likely to be driven by nongenetic factors. The epigenetic modification DNA methylation (DNAm) may help to explore this, as it is sensitive to both genetic and environmental exposures. While the relationship between DNAm and body-fat traits has been extensively studied, there is limited literature on the shared associations of DNAm variation across such traits. Akin to genetic correlation estimates, here, we introduce an approach to evaluate the similarities in DNAm associations between traits: DNAm correlations. As DNAm can be both a cause and consequence of complex traits, DNAm correlations have the potential to provide insights into trait relationships above that currently obtained from genetic and phenotypic correlations. Utilizing 7,519 unrelated individuals from Generation Scotland with DNAm from the EPIC array, we calculated DNAm correlations between body-fat- and adiposity-related traits by using the bivariate OREML framework in the OSCA software. For each trait, we also estimated the shared contribution of DNAm between sexes. We identified strong, positive DNAm correlations between each of the body-fat traits (BMI, body-fat percentage, and waist-to-hip ratio, ranging from 0.96 to 1.00), finding larger associations than those identified by genetic and phenotypic correlations. We identified a significant deviation from 1 in the DNAm correlations for BMI between males and females, with sex-specific DNAm changes associated with BMI identified at eight DNAm probes. Employing genome-wide DNAm correlations to evaluate the similarities in the associations of DNAm with complex traits has provided insight into obesity-related traits beyond that provided by genetic correlations.

Impact of cross-ancestry genetic architecture on GWASs in admixed populations.

Genome-wide association studies (GWASs) have identified thousands of variants for disease risk. These studies have predominantly been conducted in individuals of European ancestries, which raises questions about their transferability to individuals of other ancestries. Of particular interest are admixed populations, usually defined as populations with recent ancestry from two or more continental sources. Admixed genomes contain segments of distinct ancestries that vary in composition across individuals in the population, allowing for the same allele to induce risk for disease on different ancestral backgrounds. This mosaicism raises unique challenges for GWASs in admixed populations, such as the need to correctly adjust for population stratification. In this work we quantify the impact of differences in estimated allelic effect sizes for risk variants between ancestry backgrounds on association statistics. Specifically, while the possibility of estimated allelic effect-size heterogeneity by ancestry (HetLanc) can be modeled when performing a GWAS in admixed populations, the extent of HetLanc needed to overcome the penalty from an additional degree of freedom in the association statistic has not been thoroughly quantified. Using extensive simulations of admixed genotypes and phenotypes, we find that controlling for and conditioning effect sizes on local ancestry can reduce statistical power by up to 72%. This finding is especially pronounced in the presence of allele frequency differentiation. We replicate simulation results using 4,327 African-European admixed genomes from the UK Biobank for 12 traits to find that for most significant SNPs, HetLanc is not large enough for GWASs to benefit from modeling heterogeneity in this way.

mtPGS: Leverage multiple correlated traits for accurate polygenic score construction.

Accurate polygenic scores (PGSs) facilitate the genetic prediction of complex traits and aid in the development of personalized medicine. Here, we develop a statistical method called multi-trait assisted PGS (mtPGS), which can construct accurate PGSs for a target trait of interest by leveraging multiple traits relevant to the target trait. Specifically, mtPGS borrows SNP effect size similarity information between the target trait and its relevant traits to improve the effect size estimation on the target trait, thus achieving accurate PGSs. In the process, mtPGS flexibly models the shared genetic architecture between the target and the relevant traits to achieve robust performance, while explicitly accounting for the environmental covariance among them to accommodate different study designs with various sample overlap patterns. In addition, mtPGS uses only summary statistics as input and relies on a deterministic algorithm with several algebraic techniques for scalable computation. We evaluate the performance of mtPGS through comprehensive simulations and applications to 25 traits in the UK Biobank, where in the real data mtPGS achieves an average of 0.90%-52.91% accuracy gain compared to the state-of-the-art PGS methods. Overall, mtPGS represents an accurate, fast, and robust solution for PGS construction in biobank-scale datasets.

Improving multi-population genomic prediction accuracy using multi-trait GBLUP models which incorporate global or local genetic correlation information.

In the application of genomic prediction, a situation often faced is that there are multiple populations in which genomic prediction (GP) need to be conducted. A common way to handle the multi-population GP is simply to combine the multiple populations into a single population. However, since these populations may be subject to different environments, there may exist genotype-environment interactions which may affect the accuracy of genomic prediction. In this study, we demonstrated that multi-trait genomic best linear unbiased prediction (MTGBLUP) can be used for multi-population genomic prediction, whereby the performances of a trait in different populations are regarded as different traits, and thus multi-population prediction is regarded as multi-trait prediction by employing the between-population genetic correlation. Using real datasets, we proved that MTGBLUP outperformed the conventional multi-population model that simply combines different populations together. We further proposed that MTGBLUP can be improved by partitioning the global between-population genetic correlation into local genetic correlations (LGC). We suggested two LGC models, LGC-model-1 and LGC-model-2, which partition the genome into regions with and without significant LGC (LGC-model-1) or regions with and without strong LGC (LGC-model-2). In analysis of real datasets, we demonstrated that the LGC models could increase universally the prediction accuracy and the relative improvement over MTGBLUP reached up to 163.86% (25.64% on average).

Contraception ends the genetic maintenance of human same-sex sexual behavior.

Because human same-sex sexual behavior (SSB) is heritable and leads to fewer offspring, it is puzzling why SSB-associated alleles have not been selectively purged. Current evidence supports the antagonistic pleiotropy hypothesis that SSB-associated alleles benefit individuals exclusively performing opposite-sex sexual behavior by increasing their number of sexual partners and consequently their number of offspring. However, by analyzing the UK Biobank, here, we show that having more sexual partners no longer predicts more offspring since the availability of oral contraceptives in the 1960s and that SSB is now genetically negatively correlated with the number of offspring, suggesting a loss of SSB's genetic maintenance in modern societies.

Fast estimation of genetic correlation for biobank-scale data.

Genetic correlation is an important parameter in efforts to understand the relationships among complex traits. Current methods that analyze individual genotype data for estimating genetic correlation are challenging to scale to large datasets. Methods that analyze summary data, while being computationally efficient, tend to yield estimates of genetic correlation with reduced precision. We propose SCORE (scalable genetic correlation estimator), a randomized method of moments estimator of genetic correlation that is both scalable and accurate. SCORE obtains more precise estimates of genetic correlations relative to summary-statistic methods that can be applied at scale; it achieves a 44% reduction in standard error relative to LD-score regression (LDSC) and a 20% reduction relative to high-definition likelihood (HDL) (averaged over all simulations). The efficiency of SCORE enables computation of genetic correlations on the UK Biobank dataset, consisting of ≈300 K individuals and ≈500 K SNPs, in a few h (orders of magnitude faster than methods that analyze individual data, such as GCTA). Across 780 pairs of traits in 291,273 unrelated white British individuals in the UK Biobank, SCORE identifies significant genetic correlation between 200 additional pairs of traits over LDSC (beyond the 245 pairs identified by both).

Genetic interactions drive heterogeneity in causal variant effect sizes for gene expression and complex traits.

Despite the growing number of genome-wide association studies (GWASs), it remains unclear to what extent gene-by-gene and gene-by-environment interactions influence complex traits in humans. The magnitude of genetic interactions in complex traits has been difficult to quantify because GWASs are generally underpowered to detect individual interactions of small effect. Here, we develop a method to test for genetic interactions that aggregates information across all trait-associated loci. Specifically, we test whether SNPs in regions of European ancestry shared between European American and admixed African American individuals have the same causal effect sizes. We hypothesize that in African Americans, the presence of genetic interactions will drive the causal effect sizes of SNPs in regions of European ancestry to be more similar to those of SNPs in regions of African ancestry. We apply our method to two traits: gene expression in 296 African Americans and 482 European Americans in the Multi-Ethnic Study of Atherosclerosis (MESA) and low-density lipoprotein cholesterol (LDL-C) in 74K African Americans and 296K European Americans in the Million Veteran Program (MVP). We find significant evidence for genetic interactions in our analysis of gene expression; for LDL-C, we observe a similar point estimate, although this is not significant, most likely due to lower statistical power. These results suggest that gene-by-gene or gene-by-environment interactions modify the effect sizes of causal variants in human complex traits.

Investigating the shared genetic architecture of uterine leiomyoma and breast cancer: A genome-wide cross-trait analysis.

Little is known regarding the shared genetic architecture or causality underlying the phenotypic association observed for uterine leiomyoma (UL) and breast cancer (BC). Leveraging summary statistics from the hitherto largest genome-wide association study (GWAS) conducted in each trait, we investigated the genetic overlap and causal associations of UL with BC overall, as well as with its subtypes defined by the status of estrogen receptor (ER). We observed a positive genetic correlation between UL and BC overall (rg = 0.09, p = 6.00 × 10-3), which was consistent in ER+ subtype (rg = 0.06, p = 0.01) but not in ER- subtype (rg = 0.06, p = 0.08). Partitioning the whole genome into 1,703 independent regions, local genetic correlation was identified at 22q13.1 for UL with BC overall and with ER+ subtype. Significant genetic correlation was further discovered in 9 out of 14 functional categories, with the highest estimates observed in coding, H3K9ac, and repressed regions. Cross-trait meta-analysis identified 9 novel loci shared between UL and BC. Mendelian randomization demonstrated a significantly increased risk of BC overall (OR = 1.09, 95% CI = 1.01-1.18) and ER+ subtype (OR = 1.09, 95% CI = 1.01-1.17) for genetic liability to UL. No reverse causality was found. Our comprehensive genome-wide cross-trait analysis demonstrates a shared genetic basis, pleiotropic loci, as well as a putative causal relationship between UL and BC, highlighting an intrinsic link underlying these two complex female diseases.

Benchmarking of local genetic correlation estimation methods using summary statistics from genome-wide association studies.

Local genetic correlation evaluates the correlation of additive genetic effects between different traits across the same genetic variants at a genomic locus. It has been proven informative for understanding the genetic similarities of complex traits beyond that captured by global genetic correlation calculated across the whole genome. Several summary-statistics-based approaches have been developed for estimating local genetic correlation, including $\rho$-hess, SUPERGNOVA and LAVA. However, there has not been a comprehensive evaluation of these methods to offer practical guidelines on the choices of these methods. In this study, we conduct benchmark comparisons of the performance of these three methods through extensive simulation and real data analyses. We focus on two technical difficulties in estimating local genetic correlation: sample overlaps across traits and local linkage disequilibrium (LD) estimates when only the external reference panels are available. Our simulations suggest the likelihood of incorrectly identifying correlated regions and local correlation estimation accuracy are highly dependent on the estimation of the local LD matrix. These observations are corroborated by real data analyses of 31 complex traits. Overall, our findings illuminate the distinct results yielded by different methods applied in post-genome-wide association studies (post-GWAS) local correlation studies. We underscore the sensitivity of local genetic correlation estimates and inferences to the precision of local LD estimation. These observations accentuate the vital need for ongoing refinement in methodologies.

Integrating multiple traits for improving polygenic risk prediction in disease and pharmacogenomics GWAS.

Polygenic risk score (PRS) has been recently developed for predicting complex traits and drug responses. It remains unknown whether multi-trait PRS (mtPRS) methods, by integrating information from multiple genetically correlated traits, can improve prediction accuracy and power for PRS analysis compared with single-trait PRS (stPRS) methods. In this paper, we first review commonly used mtPRS methods and find that they do not directly model the underlying genetic correlations among traits, which has been shown to be useful in guiding multi-trait association analysis in the literature. To overcome this limitation, we propose a mtPRS-PCA method to combine PRSs from multiple traits with weights obtained from performing principal component analysis (PCA) on the genetic correlation matrix. To accommodate various genetic architectures covering different effect directions, signal sparseness and across-trait correlation structures, we further propose an omnibus mtPRS method (mtPRS-O) by combining P values from mtPRS-PCA, mtPRS-ML (mtPRS based on machine learning) and stPRSs using Cauchy Combination Test. Our extensive simulation studies show that mtPRS-PCA outperforms other mtPRS methods in both disease and pharmacogenomics (PGx) genome-wide association studies (GWAS) contexts when traits are similarly correlated, with dense signal effects and in similar effect directions, and mtPRS-O is consistently superior to most other methods due to its robustness under various genetic architectures. We further apply mtPRS-PCA, mtPRS-O and other methods to PGx GWAS data from a randomized clinical trial in the cardiovascular domain and demonstrate performance improvement of mtPRS-PCA in both prediction accuracy and patient stratification as well as the robustness of mtPRS-O in PRS association test.

A comprehensive map of genetic relationships among diagnostic categories based on 48.6 million relative pairs from the Danish genealogy.

For more than half a century, Denmark has maintained population-wide demographic, health care, and socioeconomic registers that provide detailed information on the interaction between all residents and the extensive national social services system. We leverage this resource to reconstruct the genealogy of the entire nation based on all individuals legally residing in Denmark since 1968. We cross-reference 6,691,426 individuals with nationwide health care registers to estimate heritability and genetic correlations of 10 broad diagnostic categories involving all major organs and systems. Heritability estimates for mental disorders were consistently the highest across demographic cohorts (average h2 = 0.406, 95% CI = [0.403, 0.408]), whereas estimates for cancers were the lowest (average h2 = 0.130, 95% CI = [0.125, 0.134]). The average genetic correlation of each of the 10 diagnostic categories with the other nine was highest for gastrointestinal conditions (average rg = 0.567, 95% CI = [0.566, 0.567]) and lowest for urogenital conditions (average rg = 0.386, 95% CI = [0.385, 0.388]). Mental, pulmonary, gastrointestinal, and neurological conditions had similar genetic correlation profiles.

Identification of genetic risk loci and prioritization of genes and pathways for myasthenia gravis: a genome-wide association study.

Myasthenia gravis is a chronic autoimmune disease characterized by autoantibody-mediated interference of signal transmission across the neuromuscular junction. We performed a genome-wide association study (GWAS) involving 1,873 patients diagnosed with acetylcholine receptor antibody-positive myasthenia gravis and 36,370 healthy individuals to identify disease-associated genetic risk loci. Replication of the discovered loci was attempted in an independent cohort from the UK Biobank. We also performed a transcriptome-wide association study (TWAS) using expression data from skeletal muscle, whole blood, and tibial nerve to test the effects of disease-associated polymorphisms on gene expression. We discovered two signals in the genes encoding acetylcholine receptor subunits that are the most common antigenic target of the autoantibodies: a GWAS signal within the cholinergic receptor nicotinic alpha 1 subunit (CHRNA1) gene and a TWAS association with the cholinergic receptor nicotinic beta 1 subunit (CHRNB1) gene in normal skeletal muscle. Two other loci were discovered on 10p14 and 11q21, and the previous association signals at PTPN22, HLA-DQA1/HLA-B, and TNFRSF11A were confirmed. Subgroup analyses demonstrate that early- and late-onset cases have different genetic risk factors. Genetic correlation analysis confirmed a genetic link between myasthenia gravis and other autoimmune diseases, such as hypothyroidism, rheumatoid arthritis, multiple sclerosis, and type 1 diabetes. Finally, we applied Priority Index analysis to identify potentially druggable genes/proteins and pathways. This study provides insight into the genetic architecture underlying myasthenia gravis and demonstrates that genetic factors within the loci encoding acetylcholine receptor subunits contribute to its pathogenesis.

Effect of case and control definitions on genome-wide association study (GWAS) findings.

Genome-wide association studies (GWAS) have significantly advanced our understanding of the genetic underpinnings of diseases, but case and control cohort definitions for a given disease can vary between different published studies. For example, two GWAS for the same disease using the UK Biobank data set might use different data sources (i.e., self-reported questionnaires, hospital records, etc.) or different levels of granularity (i.e., specificity of inclusion criteria) to define cases and controls. The extent to which this variability in cohort definitions impacts the end-results of a GWAS study is unclear. In this study, we systematically evaluated the effect of the data sources used for case and control definitions on GWAS findings. Using the UK Biobank, we selected three diseases-glaucoma, migraine, and iron-deficiency anemia. For each disease, we designed 13 GWAS, each using different combinations of data sources to define cases and controls, and then calculated the pairwise genetic correlations between all GWAS for each disease. We found that the data sources used to define cases for a given disease can have a significant impact on GWAS end-results, but the extent of this depends heavily on the disease in question. This suggests the need for greater scrutiny on how case cohorts are defined for GWAS.

Genetic correlation between smoking behavior and gastroesophageal reflux disease: insights from integrative multi-omics data.

BACKGROUND: Observational studies have preliminarily revealed an association between smoking and gastroesophageal reflux disease (GERD). However, little is known about the causal relationship and shared genetic architecture between the two. This study aims to explore their common genetic correlations by leveraging genome-wide association studies (GWAS) of smoking behavior-specifically, smoking initiation (SI), never smoking (NS), ever smoking (ES), cigarettes smoked per day (CPD), age of smoking initiation(ASI) and GERD. METHODS: Firstly, we conducted global cross-trait genetic correlation analysis and heritability estimation from summary statistics (HESS) to explore the genetic correlation between smoking behavior and GERD. Then, a joint cross-trait meta-analysis was performed to identify shared "pleiotropic SNPs" between smoking behavior and GERD, followed by co-localization analysis. Additionally, multi-marker analyses using annotation (MAGMA) were employed to explore the degree of enrichment of single nucleotide polymorphism (SNP) heritability in specific tissues, and summary data-based Mendelian randomization (SMR) was further utilized to investigate potential functional genes. Finally, Mendelian randomization (MR) analysis was conducted to explore the causal relationship between the smoking behavior and GERD. RESULTS: Consistent genetic correlations were observed through global and local genetic correlation analyses, wherein SI, ES, and CPD showed significantly positive genetic correlations with GERD, while NS and ASI showed significantly negative correlations. HESS analysis also identified multiple significantly associated loci between them. Furthermore, three novel "pleiotropic SNPs" (rs4382592, rs200968, rs1510719) were identified through cross-trait meta-analysis and co-localization analysis to exist between SI, NS, ES, ASI, and GERD, mapping the genes MED27, HIST1H2BO, MAML3 as new pleiotropic genes between SI, NS, ES, ASI, and GERD. Moreover, both smoking behavior and GERD were found to be co-enriched in multiple brain tissues, with GMPPB, RNF123, and RBM6 identified as potential functional genes co-enriched in Cerebellar Hemisphere, Cerebellum, Cortex/Nucleus accumbens in SI and GERD, and SUOX identified in Caudate nucleus, Cerebellum, Cortex in NS and GERD. Lastly, consistent causal relationships were found through MR analysis, indicating that SI, ES, and CPD increase the risk of GERD, while NS and higher ASI decrease the risk. CONCLUSION: We identified genetic loci associated with smoking behavior and GERD, as well as brain tissue sites of shared enrichment, prioritizing three new pleiotropic genes and four new functional genes. Finally, the causal relationship between smoking behavior and GERD was demonstrated, providing insights for early prevention strategies for GERD.

Methodological artefacts cause counter-intuitive evolutionary conclusions in a simulation study.

In their simulation study, Garcia-Costoya et al. (2023) conclude that evolutionary constraints might aid populations facing climate change. However, we are concerned that this conclusion is largely a consequence of the simulated temperature variation being too small, and, most importantly, that uneven limitations to standing variation disadvantage unconstrained populations.

Two Forms of Sexual Dimorphism in Gene Expression in Drosophila melanogaster: Their Coincidence and Evolutionary Genetics.

Phenotypic sexual dimorphism can be mediated by sex differences in gene expression. We examine two forms of sexual dimorphism in gene expression in Drosophila melanogaster: 1) sex-biased gene expression (SBGE) in which the sexes differ in the amount a gene is expressed and 2) sexual dimorphism in isoform usage, that is, sex-specific splicing (SSS). In whole body (but not the head) expression, we find a negative association between SBGE and SSS, possibly suggesting that these are alternate routes to resolving sexual antagonistic selection. Next, we evaluate whether expression dimorphism contributes to the heterogeneity among genes in rmf, the intersexual genetic correlation in body expression that constrains the extent to which a gene's expression can evolve independently between the sexes. We find lower rmf values for genes with than without SSS. We find higher rmf values for male- than female-biased genes (except genes with extreme male bias), even though male-biased genes are known to have greater evolutionary divergence in expression. Finally, we examine population genetic patterns in relation to SBGE and SSS because genes with expression dimorphism have likely experienced a history of sex differences in selection. SSS is associated with reduced values of Tajima's D and elevated direction of selection (DoS) values, suggestive of higher rates of adaptive evolution. Though DoS is highly elevated for genes with extreme male bias, DoS otherwise tends to decline from female-biased to unbiased to male-biased genes. Collectively, the results indicate that SBGE and SSS are differentially distributed across the genome and are associated with different forms of selection.

Phenotypic and genetic effect of carotid intima-media thickness on the risk of stroke.

While carotid intima-media thickness (cIMT) as a noninvasive surrogate measure of atherosclerosis is widely considered a risk factor for stroke, the intrinsic link underlying cIMT and stroke has not been fully understood. We aimed to evaluate the clinical value of cIMT in stroke through the investigation of phenotypic and genetic relationships between cIMT and stroke. We evaluated phenotypic associations using observational data from UK Biobank (N = 21,526). We then investigated genetic relationships leveraging genomic data conducted in predominantly European ancestry for cIMT (N = 45,185) and any stroke (AS, Ncase/Ncontrol=40,585/406,111). Observational analyses suggested an increased hazard of stroke per one standard deviation increase in cIMT (cIMTmax-AS: hazard ratio (HR) = 1.39, 95%CI = 1.09-1.79; cIMTmean-AS: HR = 1.39, 95%CI = 1.09-1.78; cIMTmin-AS: HR = 1.32, 95%CI = 1.04-1.68). A positive global genetic correlation was observed (cIMTmax-AS: [Formula: see text]=0.23, P=9.44 × 10-5; cIMTmean-AS: [Formula: see text]=0.21, P=3.00 × 10-4; cIMTmin-AS: [Formula: see text]=0.16, P=6.30 × 10-3). This was further substantiated by five shared independent loci and 15 shared expression-trait associations. Mendelian randomization analyses suggested no causal effect of cIMT on stroke (cIMTmax-AS: odds ratio (OR)=1.12, 95%CI=0.97-1.28; cIMTmean-AS: OR=1.09, 95%CI=0.93-1.26; cIMTmin-AS: OR=1.03, 95%CI = 0.90-1.17). A putative association was observed for genetically predicted stroke on cIMT (AS-cIMTmax: beta=0.07, 95%CI = 0.01-0.13; AS-cIMTmean: beta=0.08, 95%CI = 0.01-0.15; AS-cIMTmin: beta = 0.08, 95%CI = 0.01-0.16) in the reverse direction MR, which attenuated to non-significant in sensitivity analysis. Our work does not find evidence supporting causal associations between cIMT and stroke. The pronounced cIMT-stroke association is intrinsic, and mostly attributed to shared genetic components. The clinical value of cIMT as a surrogate marker for stroke risk in the general population is likely limited.