Long-range LD can confound genome scans in admixed populations

Sequence variants influencing the regulation of serum IgG subclass levels

Immunoglobulin G (IgG) is the main isotype of antibody in human blood. IgG consists of four subclasses (IgG1 to IgG4), encoded by separate constant region genes within the Ig heavy chain locus (IGH). Here, we report a genome-wide association study on blood IgG subclass levels. Across 4334 adults and 4571 individuals under 18 years, we discover ten new and identify four known variants at five loci influencing IgG subclass levels. These variants also affect the risk of asthma, autoimmune diseases, and blood traits. Seven variants map to the IGH locus, three to the Fcγ receptor (FCGR) locus, and two to the human leukocyte antigen (HLA) region, affecting the levels of all IgG subclasses. The most significant associations are observed between the G1m (f), G2m(n) and G3m(b*) allotypes, and IgG1, IgG2 and IgG3, respectively. Additionally, we describe selective associations with IgG4 at 16p11.2 (ITGAX) and 17q21.1 (IKZF3, ZPBP2, GSDMB, ORMDL3). Interestingly, the latter coincides with a highly pleiotropic signal where the allele associated with lower IgG4 levels protects against childhood asthma but predisposes to inflammatory bowel disease. Our results provide insight into the regulation of antibody-mediated immunity that can potentially be useful in the development of antibody based therapeutics.

Genetic and functional analysis of Raynaud's syndrome implicates loci in vasculature and immunity

Raynaud's syndrome is a dysautonomia where exposure to cold causes vasoconstriction and hypoxia, particularly in the extremities. We performed meta-analysis in four cohorts and discovered eight loci (ADRA2A, IRX1, NOS3, ACVR2A, TMEM51, PCDH10-DT, HLA, and RAB6C) where ADRA2A, ACVR2A, NOS3, TMEM51, and IRX1 co-localized with expression quantitative trait loci (eQTLs), particularly in distal arteries. CRISPR gene editing further showed that ADRA2A and NOS3 loci modified gene expression and in situ RNAscope clarified the specificity of ADRA2A in small vessels and IRX1 around small capillaries in the skin. A functional contraction assay in the cold showed lower contraction in ADRA2A-deficient and higher contraction in ADRA2A-overexpressing smooth muscle cells. Overall, our study highlights the power of genome-wide association testing with functional follow-up as a method to understand complex diseases. The results indicate temperature-dependent adrenergic signaling through ADRA2A, effects at the microvasculature by IRX1, endothelial signaling by NOS3, and immune mechanisms by the HLA locus in Raynaud's syndrome.

Unveiling Promising Neuroimaging Biomarkers for Schizophrenia Through Clinical and Genetic Perspectives

Schizophrenia is a complex and serious brain disorder. Neuroscientists have become increasingly interested in using magnetic resonance-based brain imaging-derived phenotypes (IDPs) to investigate the etiology of psychiatric disorders. IDPs capture valuable clinical advantages and hold biological significance in identifying brain abnormalities. In this review, we aim to discuss current and prospective approaches to identify potential biomarkers for schizophrenia using clinical multimodal neuroimaging and imaging genetics. We first described IDPs through their phenotypic classification and neuroimaging genomics. Secondly, we discussed the applications of multimodal neuroimaging by clinical evidence in observational studies and randomized controlled trials. Thirdly, considering the genetic evidence of IDPs, we discussed how can utilize neuroimaging data as an intermediate phenotype to make association inferences by polygenic risk scores and Mendelian randomization. Finally, we discussed machine learning as an optimum approach for validating biomarkers. Together, future research efforts focused on neuroimaging biomarkers aim to enhance our understanding of schizophrenia.

Ancestry-aligned polygenic scores combined with conventional risk factors improve prediction of cardiometabolic outcomes in African populations

BACKGROUND

Cardiovascular diseases (CVD) are a major health concern in Africa. Improved identification and treatment of high-risk individuals can reduce adverse health outcomes. Current CVD risk calculators are largely unvalidated in African populations and overlook genetic factors. Polygenic scores (PGS) can enhance risk prediction by measuring genetic susceptibility to CVD, but their effectiveness in genetically diverse populations is limited by a European-ancestry bias. To address this, we developed models integrating genetic data and conventional risk factors to assess the risk of developing cardiometabolic outcomes in African populations.

METHODS

We used summary statistics from a genome-wide association meta-analysis (n = 14,126) in African populations to derive novel genome-wide PGS for 14 cardiometabolic traits in an independent African target sample (Africa Wits-INDEPTH Partnership for Genomic Research (AWI-Gen), n = 10,603). Regression analyses assessed relationships between each PGS and corresponding cardiometabolic trait, and seven CVD outcomes (CVD, heart attack, stroke, diabetes mellitus, dyslipidaemia, hypertension, and obesity). The predictive utility of the genetic data was evaluated using elastic net models containing multiple PGS (MultiPGS) and reference-projected principal components of ancestry (PPCs). An integrated risk prediction model incorporating genetic and conventional risk factors was developed. Nested cross-validation was used when deriving elastic net models to enhance generalisability.

RESULTS

Our African-specific PGS displayed significant but variable within- and cross- trait prediction (max.R2 = 6.8%, p = 1.86 × 10-173). Significantly associated PGS with dyslipidaemia included the PGS for total cholesterol (logOR = 0.210, SE = 0.022, p = 2.18 × 10-21) and low-density lipoprotein (logOR =  - 0.141, SE = 0.022, p = 1.30 × 10-20); with hypertension, the systolic blood pressure PGS (logOR = 0.150, SE = 0.045, p = 8.34 × 10-4); and multiple PGS associated with obesity: body mass index (max. logOR = 0.131, SE = 0.031, p = 2.22 × 10-5), hip circumference (logOR = 0.122, SE = 0.029, p = 2.28 × 10-5), waist circumference (logOR = 0.013, SE = 0.098, p = 8.13 × 10-4) and weight (logOR = 0.103, SE = 0.029, p = 4.89 × 10-5). Elastic net models incorporating MultiPGS and PPCs significantly improved prediction over MultiPGS alone. Models including genetic data and conventional risk factors were more predictive than conventional risk models alone (dyslipidaemia: R2 increase = 2.6%, p = 4.45 × 10-12; hypertension: R2 increase = 2.6%, p = 2.37 × 10-13; obesity: R2 increase = 5.5%, 1.33 × 10-34).

CONCLUSIONS

In African populations, CVD and associated cardiometabolic trait prediction models can be improved by incorporating ancestry-aligned PGS and accounting for ancestry. Combining PGS with conventional risk factors further enhances prediction over traditional models based on conventional factors. Incorporating data from target populations can improve the generalisability of international predictive models for CVD and associated traits in African populations.

Human genetic structure in Northwest France provides new insights into West European historical demography

The demographical history of France remains largely understudied despite its central role toward understanding modern population structure across Western Europe. Here, by exploring publicly available Europe-wide genotype datasets together with the genomes of 3234 present-day and six newly sequenced medieval individuals from Northern France, we found extensive fine-scale population structure across Brittany and the downstream Loire basin and increased population differentiation between the northern and southern sides of the river Loire, associated with higher proportions of steppe vs. Neolithic-related ancestry. We also found increased allele sharing between individuals from Western Brittany and those associated with the Bell Beaker complex. Our results emphasise the need for investigating local populations to better understand the distribution of rare (putatively deleterious) variants across space and the importance of common genetic legacy in understanding the sharing of disease-related alleles between Brittany and people from western Britain and Ireland.

Evidence for dynastic succession among early Celtic elites in Central Europe

The early Iron Age (800 to 450 BCE) in France, Germany and Switzerland, known as the 'West-Hallstattkreis', stands out as featuring the earliest evidence for supra-regional organization north of the Alps. Often referred to as 'early Celtic', suggesting tentative connections to later cultural phenomena, its societal and population structure remain enigmatic. Here we present genomic and isotope data from 31 individuals from this context in southern Germany, dating between 616 and 200 BCE. We identify multiple biologically related groups spanning three elite burials as far as 100 km apart, supported by trans-regional individual mobility inferred from isotope data. These include a close biological relationship between two of the richest burial mounds of the Hallstatt culture. Bayesian modelling points to an avuncular relationship between the two individuals, which may suggest a practice of matrilineal dynastic succession in early Celtic elites. We show that their ancestry is shared on a broad geographic scale from Iberia throughout Central-Eastern Europe, undergoing a decline after the late Iron Age (450 BCE to ~50 CE).

Associations Between Temporomandibular Disorders and Brain Imaging-Derived Phenotypes

OBJECTIVE

Temporomandibular disorders (TMD) affect the temporomandibular joint and associated structures. Despite its prevalence and impact on quality of life, the underlying mechanisms of TMD remain unclear. Magnetic resonance imaging studies suggest brain abnormalities in patients with TMD. However, these lines of evidence are essentially observational and cannot infer a causal relationship. This study employs Mendelian randomisation (MR) to probe causal relationships between TMD and brain changes.

METHODS

Genome-wide association study (GWAS) summary statistics for TMD were collected, along with brain imaging-derived phenotypes (IDPs). Instrumental variables were selected from the GWAS summary statistics and used in bidirectional 2-sample MR analyses. The inverse-variance weighted analysis was chosen as the primary method. In addition, false discovery rate (FDR) correction of P value was used.

RESULTS

Eleven IDPs related to brain imaging alterations showed significant causal associations with TMD (P-FDR < .05), validated through sensitivity analysis. In forward MR, the mean thickness of left caudal middle frontal gyrus (OR, 0.76; 95% CI, 0.67-0.87; P-FDR = 1.15  ×  10-2) and the volume of right superior frontal gyrus (OR, 1.24; 95% CI, 1.10-1.39; P-FDR = 2.26  ×  10-2) exerted significant causal effects on TMD. In the reverse MR analysis, TMD exerted a significant causal effect on 9 IDPs, including the mean thickness of the left medial orbitofrontal cortex (β = -0.10; 95% CI, -0.13 to -0.08; P-FDR = 2.06 × 10-11), the volume of the left magnocellular nucleus (β = -0.15; 95% CI, -0.22 to -0.09; P-FDR = 3.26 × 10-4), the mean intensity of the right inferior-lateral ventricle (β = -0.09; 95% CI, -0.14 to -0.04; P-FDR = 2.23 × 10-2), the volume of grey matter in the anterior division of the left superior temporal gyrus (β = 0.09; 95% CI,  0.04-0.14; P-FDR = 1.69 × 10-2), and so forth.

CONCLUSIONS

This study provides genetic evidence supporting the bidirectional causal associations between TMD and brain IDPs, shedding light on potential neurobiological mechanisms underlying TMD development and its relationship with brain structure.

Chromosome X-Wide Common Variant Association Study (XWAS) in Autism Spectrum Disorder

Autism Spectrum Disorder (ASD) displays a notable male bias in prevalence. Research into rare (<0.1) genetic variants on the X chromosome has implicated over 20 genes in ASD pathogenesis, such as MECP2, DDX3X, and DMD. The "female protective effect" in ASD suggests that females may require a higher genetic burden to manifest similar symptoms as males, yet the mechanisms remain unclear. Despite technological advances in genomics, the complexity of the biological nature of sex chromosomes leave them underrepresented in genome-wide studies. Here, we conducted an X chromosome-wide association study (XWAS) using whole-genome sequencing data from 6,873 individuals with ASD (82% males) across Autism Speaks MSSNG, Simons Simplex Cohort SSC, and Simons Foundation Powering Autism Research SPARK, alongside 8,981 population controls (43% males). We analyzed 418,652 X-chromosome variants, identifying 59 associated with ASD (p-values 7.9×10-6 to 1.51×10-5), surpassing Bonferroni-corrected thresholds. Key findings include significant regions on chrXp22.2 (lead SNP=rs12687599, p=3.57×10-7) harboring ASB9/ASB11, and another encompassing DDX53/PTCHD1-AS long non-coding RNA (lead SNP=rs5926125, p=9.47×10-6). When mapping genes within 10kb of the 59 most significantly associated SNPs, 91 genes were found, 17 of which yielded association with ASD (GRPR, AP1S2, DDX53, HDAC8, PCDH19, PTCHD1, PCDH11X, PTCHD1-AS, DMD, SYAP1, CNKSR2, GLRA2, OFD1, CDKL5, GPRASP2, NXF5, SH3KBP1). FGF13 emerged as a novel X-linked ASD candidate gene, highlighted by sex-specific differences in minor allele frequencies. These results reveal significant new insights into X chromosome biology in ASD, confirming and nominating genes and pathways for further investigation.

Evaluation of polygenic scoring methods in five biobanks shows larger variation between biobanks than methods and finds benefits of ensemble learning

Methods of estimating polygenic scores (PGSs) from genome-wide association studies are increasingly utilized. However, independent method evaluation is lacking, and method comparisons are often limited. Here, we evaluate polygenic scores derived via seven methods in five biobank studies (totaling about 1.2 million participants) across 16 diseases and quantitative traits, building on a reference-standardized framework. We conducted meta-analyses to quantify the effects of method choice, hyperparameter tuning, method ensembling, and the target biobank on PGS performance. We found that no single method consistently outperformed all others. PGS effect sizes were more variable between biobanks than between methods within biobanks when methods were well tuned. Differences between methods were largest for the two investigated autoimmune diseases, seropositive rheumatoid arthritis and type 1 diabetes. For most methods, cross-validation was more reliable for tuning hyperparameters than automatic tuning (without the use of target data). For a given target phenotype, elastic net models combining PGS across methods (ensemble PGS) tuned in the UK Biobank provided consistent, high, and cross-biobank transferable performance, increasing PGS effect sizes (β coefficients) by a median of 5.0% relative to LDpred2 and MegaPRS (the two best-performing single methods when tuned with cross-validation). Our interactively browsable online-results and open-source workflow prspipe provide a rich resource and reference for the analysis of polygenic scoring methods across biobanks.

Phenotypic and ancestry-related assortative mating in autism

BACKGROUND

Positive assortative mating (AM) in several neuropsychiatric traits, including autism, has been noted. However, it is unknown whether the pattern of AM is different in phenotypically defined autism subgroups [e.g., autism with and without intellectually disability (ID)]. It is also unclear what proportion of the phenotypic AM can be explained by the genetic similarity between parents of children with an autism diagnosis, and the consequences of AM on the genetic structure of the population.

METHODS

To address these questions, we analyzed two family-based autism collections: the Simons Foundation Powering Autism Research for Knowledge (SPARK) (1575 families) and the Simons Simplex Collection (SSC) (2283 families).

RESULTS

We found a similar degree of phenotypic and ancestry-related AM in parents of children with an autism diagnosis regardless of the presence of ID. We did not find evidence of AM for autism based on autism polygenic scores (PGS) (at a threshold of |r|> 0.1). The adjustment of ancestry-related AM or autism PGS accounted for only 0.3-4% of the fractional change in the estimate of the phenotypic AM. The ancestry-related AM introduced higher long-range linkage disequilibrium (LD) between single nucleotide polymorphisms (SNPs) on different chromosomes that are highly ancestry-informative compared to SNPs that are less ancestry-informative (D2 on the order of 1 × 10-5).

LIMITATIONS

We only analyzed participants of European ancestry, limiting the generalizability of our results to individuals of non-European ancestry. SPARK and SSC were both multicenter studies. Therefore, there could be ancestry-related AM in SPARK and SSC due to geographic stratification. The study participants from each site were unknown, so we were unable to evaluate for geographic stratification.

CONCLUSIONS

This study showed similar patterns of AM in autism with and without ID, and demonstrated that the common genetic influences of autism are likely relevant to both autism groups. The adjustment of ancestry-related AM and autism PGS accounted for < 5% of the fractional change in the estimate of the phenotypic AM. Future studies are needed to evaluate if the small increase of long-range LD induced by ancestry-related AM has impact on the downstream analysis.

Genetic evidence for the causal relationships between migraine, dementia, and longitudinal brain atrophy

BACKGROUND

Migraine is a neurological disease with a significant genetic component and is characterized by recurrent and prolonged episodes of headache. Previous epidemiological studies have reported a higher risk of dementia in migraine patients. Neuroimaging studies have also shown structural brain atrophy in regions that are common to migraine and dementia. However, these studies are observational and cannot establish causality. The present study aims to explore the genetic causal relationship between migraine and dementia, as well as the mediation roles of brain structural changes in this association using Mendelian randomization (MR).

METHODS

We collected the genome-wide association study (GWAS) summary statistics of migraine and its two subtypes, as well as four common types of dementia, including Alzheimer's disease (AD), vascular dementia, frontotemporal dementia, and Lewy body dementia. In addition, we collected the GWAS summary statistics of seven longitudinal brain measures that characterize brain structural alterations with age. Using these GWAS, we performed Two-sample MR analyses to investigate the causal effects of migraine and its two subtypes on dementia and brain structural changes. To explore the possible mediation of brain structural changes between migraine and dementia, we conducted a two-step MR mediation analysis.

RESULTS

The MR analysis demonstrated a significant association between genetically predicted migraine and an increased risk of AD (OR = 1.097, 95% CI = [1.040, 1.158], p = 7.03 × 10- 4). Moreover, migraine significantly accelerated annual atrophy of the total cortical surface area (-65.588 cm2 per year, 95% CI = [-103.112, -28.064], p = 6.13 × 10- 4) and thalamic volume (-9.507 cm3 per year, 95% CI = [-15.512, -3.502], p = 1.91 × 10- 3). The migraine without aura (MO) subtype increased the risk of AD (OR = 1.091, 95% CI = [1.059, 1.123], p = 6.95 × 10- 9) and accelerated annual atrophy of the total cortical surface area (-31.401 cm2 per year, 95% CI = [-43.990, -18.811], p = 1.02 × 10- 6). The two-step MR mediation analysis revealed that thalamic atrophy partly mediated the causal effect of migraine on AD, accounting for 28.2% of the total effect.

DISCUSSION

This comprehensive MR study provided genetic evidence for the causal effect of migraine on AD and identified longitudinal thalamic atrophy as a potential mediator in this association. These findings may inform brain intervention targets to prevent AD risk in migraine patients.

MultiSuSiE improves multi-ancestry fine-mapping in All of Us whole-genome sequencing data

Leveraging data from multiple ancestries can greatly improve fine-mapping power due to differences in linkage disequilibrium and allele frequencies. We propose MultiSuSiE, an extension of the sum of single effects model (SuSiE) to multiple ancestries that allows causal effect sizes to vary across ancestries based on a multivariate normal prior informed by empirical data. We evaluated MultiSuSiE via simulations and analyses of 14 quantitative traits leveraging whole-genome sequencing data in 47k African-ancestry and 94k European-ancestry individuals from All of Us. In simulations, MultiSuSiE applied to Afr47k+Eur47k was well-calibrated and attained higher power than SuSiE applied to Eur94k; interestingly, higher causal variant PIPs in Afr47k compared to Eur47k were entirely explained by differences in the extent of LD quantified by LD 4th moments. Compared to very recently proposed multi-ancestry fine-mapping methods, MultiSuSiE attained higher power and/or much lower computational costs, making the analysis of large-scale All of Us data feasible. In real trait analyses, MultiSuSiE applied to Afr47k+Eur94k identified 579 fine-mapped variants with PIP > 0.5, and MultiSuSiE applied to Afr47k+Eur47k identified 44% more fine-mapped variants with PIP > 0.5 than SuSiE applied to Eur94k. We validated MultiSuSiE results for real traits via functional enrichment of fine-mapped variants. We highlight several examples where MultiSuSiE implicates well-studied or biologically plausible fine-mapped variants that were not implicated by other methods.

Breast Cancer Polygenic Risk Score Validation and Effects of Variable Imputation

Breast cancer (BC) is a complex disease affecting one in eight women in the USA. Advances in population genomics have led to the development of polygenic risk scores (PRSs) with the potential to augment current risk models, but replication is often limited. We evaluated 2 robust PRSs with 313 and 3820 SNPs and the effects of multiple genotype imputation replications in BC cases and control populations. Biological samples from BC cases and cancer-free controls were drawn from three European ancestry cohorts. Genotyping on the Illumina Global Screening Array was followed by stringent quality control measures and 20 genotype imputation replications. A total of 468 unrelated cases and 4337 controls were scored, revealing significant differences in mean PRS percentiles between cases and controls (p < 0.001) for both SNP sets (313-SNP PRS: 52.81 and 48.07; 3820-SNP PRS: 55.45 and 49.81), with receiver operating characteristic curve analysis showing area under the curve values of 0.596 and 0.603 for the 313-SNP and 3820-SNP PRS, respectively. PRS fluctuations (from ~2-3% up to 9%) emerged across imputation iterations. Our study robustly reaffirms the predictive capacity of PRSs for BC by replicating their performance in an independent BC population and showcases the need to average imputed scores for reliable outcomes.

Improved multi-ancestry fine-mapping identifies cis-regulatory variants underlying molecular traits and disease risk

Multi-ancestry statistical fine-mapping of cis-molecular quantitative trait loci (cis-molQTL) aims to improve the precision of distinguishing causal cis-molQTLs from tagging variants. However, existing approaches fail to reflect shared genetic architectures. To solve this limitation, we present the Sum of Shared Single Effects (SuShiE) model, which leverages LD heterogeneity to improve fine-mapping precision, infer cross-ancestry effect size correlations, and estimate ancestry-specific expression prediction weights. We apply SuShiE to mRNA expression measured in PBMCs (n=956) and LCLs (n=814) together with plasma protein levels (n=854) from individuals of diverse ancestries in the TOPMed MESA and GENOA studies. We find SuShiE fine-maps cis-molQTLs for 16% more genes compared with baselines while prioritizing fewer variants with greater functional enrichment. SuShiE infers highly consistent cis-molQTL architectures across ancestries on average; however, we also find evidence of heterogeneity at genes with predicted loss-of-function intolerance, suggesting that environmental interactions may partially explain differences in cis-molQTL effect sizes across ancestries. Lastly, we leverage estimated cis-molQTL effect-sizes to perform individual-level TWAS and PWAS on six white blood cell-related traits in AOU Biobank individuals (n=86k), and identify 44 more genes compared with baselines, further highlighting its benefits in identifying genes relevant for complex disease risk. Overall, SuShiE provides new insights into the cis-genetic architecture of molecular traits.

Adjusting for principal components can induce spurious associations in genome-wide association studies in admixed populations

Principal component analysis (PCA) is widely used to control for population structure in genome-wide association studies (GWAS). Top principal components (PCs) typically reflect population structure, but challenges arise in deciding how many PCs are needed and ensuring that PCs do not capture other artifacts such as regions with atypical linkage disequilibrium (LD). In response to the latter, many groups suggest performing LD pruning or excluding known high LD regions prior to PCA. However, these suggestions are not universally implemented and the implications for GWAS are not fully understood, especially in the context of admixed populations. In this paper, we investigate the impact of pre-processing and the number of PCs included in GWAS models in African American samples from the Women's Women's Health Initiative SNP Health Association Resource and two Trans-Omics for Precision Medicine Whole Genome Sequencing Project contributing studies (Jackson Heart Study and Genetic Epidemiology of Chronic Obstructive Pulmonary Disease Study). In all three samples, we find the first PC is highly correlated with genome-wide ancestry whereas later PCs often capture local genomic features. The pattern of which, and how many, genetic variants are highly correlated with individual PCs differs from what has been observed in prior studies focused on European populations and leads to distinct downstream consequences: adjusting for such PCs yields biased effect size estimates and elevated rates of spurious associations due to the phenomenon of collider bias. Excluding high LD regions identified in previous studies does not resolve these issues. LD pruning proves more effective, but the optimal choice of thresholds varies across datasets. Altogether, our work highlights unique issues that arise when using PCA to control for ancestral heterogeneity in admixed populations and demonstrates the importance of careful pre-processing and diagnostics to ensure that PCs capturing multiple local genomic features are not included in GWAS models.

Association between rheumatoid arthritis and autoimmune thyroid disease: evidence from complementary genetic methods

OBJECTIVES

To assess the causal association of Rheumatoid Arthritis (RA) with Autoimmune thyroid disease (AITD).

METHOD

Complementary genetic approaches, including genetic correlation, Mendelian randomization (MR) and colocalization analysis, were conducted to assess the potential causal association between RA and AITD using summary statistics from large-scale genome-wide association studies (GWASs). Various sensitivity analyses had been conducted to assess the robustness and the consistency of the findings.

RESULTS

The linkage disequilibrium score regression revealed a shared genetic structure between RA and AITD, with the significant genetic correlation between RA and autoimmune hyperthyroidism and autoimmune hypothyroidism estimated to be 0.3945 (P = 2.83 × 10-6) and 0.2771 (P = 1.04 × 10-6) respectively. The results of MR analysis showed that RA had a positive causal relationship with autoimmune hypothyroidism and autoimmune hyperthyroidism. The odds ratio (OR) were 1.29 (95% CI, 1.17-1.42; P = 1.08 × 10-7) and 1.47 (95% CI, 1.25-1.72; P = 1.85 × 10-6), respectively. In reverse MR analysis, autoimmune hypothyroidism had a positive causal relationship with RA, OR was 1.51 (95% CI, 1.37-1.66; P = 1.10 × 10-16); autoimmune hyperthyroidism had no causal relationship with RA relationship (P = 0.22). Similar results were found using different MR methods. In addition, colocalization analysis suggested that shared causal variants existed between RA and AITD.

CONCLUSIONS

Our study suggested a potentially causal effect of genetically predicted RA on autoimmune hyperthyroidism and a bidirectional causal relationship between RA and autoimmune hypothyroidism was also observed with complementary genetic approaches, which supports the importance and necessity of thyroid function screening and monitoring in RA patient management in clinical practice.

Thalamocortical organoids enable in vitro modeling of 22q11.2 microdeletion associated with neuropsychiatric disorders

Thalamic dysfunction has been implicated in multiple psychiatric disorders. We sought to study the mechanisms by which abnormalities emerge in the context of the 22q11.2 microdeletion, which confers significant genetic risk for psychiatric disorders. We investigated early stages of human thalamus development using human pluripotent stem cell-derived organoids and show that the 22q11.2 microdeletion underlies widespread transcriptional dysregulation associated with psychiatric disorders in thalamic neurons and glia, including elevated expression of FOXP2. Using an organoid co-culture model, we demonstrate that the 22q11.2 microdeletion mediates an overgrowth of thalamic axons in a FOXP2-dependent manner. Finally, we identify ROBO2 as a candidate molecular mediator of the effects of FOXP2 overexpression on thalamic axon overgrowth. Together, our study suggests that early steps in thalamic development are dysregulated in a model of genetic risk for schizophrenia and contribute to neural phenotypes in 22q11.2 deletion syndrome.

Common infections and neuroimaging markers of dementia in three UK cohort studies

INTRODUCTION

We aimed to investigate associations between common infections and neuroimaging markers of dementia risk (brain volume, hippocampal volume, white matter lesions) across three population-based studies.

METHODS

We tested associations between serology measures (pathogen serostatus, cumulative burden, continuous antibody responses) and outcomes using linear regression, including adjustments for total intracranial volume and scanner/clinic information (basic model), age, sex, ethnicity, education, socioeconomic position, alcohol, body mass index, and smoking (fully adjusted model). Interactions between serology measures and apolipoprotein E (APOE) genotype were tested. Findings were meta-analyzed across cohorts (Nmain  = 2632; NAPOE-interaction  = 1810).

RESULTS

Seropositivity to John Cunningham virus associated with smaller brain volumes in basic models (β = -3.89 mL [-5.81, -1.97], Padjusted  < 0.05); these were largely attenuated in fully adjusted models (β = -1.59 mL [-3.55, 0.36], P = 0.11). No other relationships were robust to multiple testing corrections and sensitivity analyses, but several suggestive associations were observed.

DISCUSSION

We did not find clear evidence for relationships between common infections and markers of dementia risk. Some suggestive findings warrant testing for replication.

Genetic influences on circulating retinol and its relationship to human health

Retinol is a fat-soluble vitamin that plays an essential role in many biological processes throughout the human lifespan. Here, we perform the largest genome-wide association study (GWAS) of retinol to date in up to 22,274 participants. We identify eight common variant loci associated with retinol, as well as a rare-variant signal. An integrative gene prioritisation pipeline supports novel retinol-associated genes outside of the main retinol transport complex (RBP4:TTR) related to lipid biology, energy homoeostasis, and endocrine signalling. Genetic proxies of circulating retinol were then used to estimate causal relationships with almost 20,000 clinical phenotypes via a phenome-wide Mendelian randomisation study (MR-pheWAS). The MR-pheWAS suggests that retinol may exert causal effects on inflammation, adiposity, ocular measures, the microbiome, and MRI-derived brain phenotypes, amongst several others. Conversely, circulating retinol may be causally influenced by factors including lipids and serum creatinine. Finally, we demonstrate how a retinol polygenic score could identify individuals more likely to fall outside of the normative range of circulating retinol for a given age. In summary, this study provides a comprehensive evaluation of the genetics of circulating retinol, as well as revealing traits which should be prioritised for further investigation with respect to retinol related therapies or nutritional intervention.

A genome-wide association study for survival from a multi-centre European study identified variants associated with COVID-19 risk of death

The clinical manifestations of SARS-CoV-2 infection vary widely among patients, from asymptomatic to life-threatening. Host genetics is one of the factors that contributes to this variability as previously reported by the COVID-19 Host Genetics Initiative (HGI), which identified sixteen loci associated with COVID-19 severity. Herein, we investigated the genetic determinants of COVID-19 mortality, by performing a case-only genome-wide survival analysis, 60 days after infection, of 3904 COVID-19 patients from the GEN-COVID and other European series (EGAS00001005304 study of the COVID-19 HGI). Using imputed genotype data, we carried out a survival analysis using the Cox model adjusted for age, age2, sex, series, time of infection, and the first ten principal components. We observed a genome-wide significant (P-value < 5.0 × 10-8) association of the rs117011822 variant, on chromosome 11, of rs7208524 on chromosome 17, approaching the genome-wide threshold (P-value = 5.19 × 10-8). A total of 113 variants were associated with survival at P-value < 1.0 × 10-5 and most of them regulated the expression of genes involved in immune response (e.g., CD300 and KLR genes), or in lung repair and function (e.g., FGF19 and CDH13). Overall, our results suggest that germline variants may modulate COVID-19 risk of death, possibly through the regulation of gene expression in immune response and lung function pathways.

Connecting genomic results for psychiatric disorders to human brain cell types and regions reveals convergence with functional connectivity

Understanding the temporal and spatial brain locations etiological for psychiatric disorders is essential for targeted neurobiological research. Integration of genomic insights from genome-wide association studies with single-cell transcriptomics is a powerful approach although past efforts have necessarily relied on mouse atlases. Leveraging a comprehensive atlas of the adult human brain, we prioritized cell types via the enrichment of SNP-heritabilities for brain diseases, disorders, and traits, progressing from individual cell types to brain regions. Our findings highlight specific neuronal clusters significantly enriched for the SNP-heritabilities for schizophrenia, bipolar disorder, and major depressive disorder along with intelligence, education, and neuroticism. Extrapolation of cell-type results to brain regions reveals important patterns for schizophrenia with distinct subregions in the hippocampus and amygdala exhibiting the highest significance. Cerebral cortical regions display similar enrichments despite the known prefrontal dysfunction in those with schizophrenia highlighting the importance of subcortical connectivity. Using functional MRI connectivity from cases with schizophrenia and neurotypical controls, we identified brain networks that distinguished cases from controls that also confirmed involvement of the central and lateral amygdala, hippocampal body, and prefrontal cortex. Our findings underscore the value of single-cell transcriptomics in decoding the polygenicity of psychiatric disorders and offer a promising convergence of genomic, transcriptomic, and brain imaging modalities toward common biological targets.

Longitudinal machine learning uncouples healthy aging factors from chronic disease risks

To understand human longevity, inherent aging processes must be distinguished from known etiologies leading to age-related chronic diseases. Such deconvolution is difficult to achieve because it requires tracking patients throughout their entire lives. Here, we used machine learning to infer health trajectories over the entire adulthood age range using extrapolation from electronic medical records with partial longitudinal coverage. Using this approach, our model tracked the state of patients who were healthy and free from known chronic disease risk and distinguished individuals with higher or lower longevity potential using a multivariate score. We showed that the model and the markers it uses performed consistently on data from Israeli, British and US populations. For example, mildly low neutrophil counts and alkaline phosphatase levels serve as early indicators of healthy aging that are independent of risk for major chronic diseases. We characterize the heritability and genetic associations of our longevity score and demonstrate at least 1 year of extended lifespan for parents of high-scoring patients compared to matched controls. Longitudinal modeling of healthy individuals is thereby established as a tool for understanding healthy aging and longevity.

Assortative mating and parental genetic relatedness contribute to the pathogenicity of variably expressive variants

We examined more than 97,000 families from four neurodevelopmental disease cohorts and the UK Biobank to identify phenotypic and genetic patterns in parents contributing to neurodevelopmental disease risk in children. We identified within- and cross-disorder correlations between six phenotypes in parents and children, such as obsessive-compulsive disorder (R = 0.32-0.38, p < 10-126). We also found that measures of sub-clinical autism features in parents are associated with several autism severity measures in children, including biparental mean Social Responsiveness Scale scores and proband Repetitive Behaviors Scale scores (regression coefficient = 0.14, p = 3.38 × 10-4). We further describe patterns of phenotypic similarity between spouses, where spouses show correlations for six neurological and psychiatric phenotypes, including a within-disorder correlation for depression (R = 0.24-0.68, p < 0.001) and a cross-disorder correlation between anxiety and bipolar disorder (R = 0.09-0.22, p < 10-92). Using a simulated population, we also found that assortative mating can lead to increases in disease liability over generations and the appearance of "genetic anticipation" in families carrying rare variants. We identified several families in a neurodevelopmental disease cohort where the proband inherited multiple rare variants in disease-associated genes from each of their affected parents. We further identified parental relatedness as a risk factor for neurodevelopmental disorders through its inverse relationship with variant pathogenicity and propose that parental relatedness modulates disease risk by increasing genome-wide homozygosity in children (R = 0.05-0.26, p < 0.05). Our results highlight the utility of assessing parent phenotypes and genotypes toward predicting features in children who carry rare variably expressive variants and implicate assortative mating as a risk factor for increased disease severity in these families.

Whole genomes from Angola and Mozambique inform about the origins and dispersals of major African migrations

As the continent of origin for our species, Africa harbours the highest levels of diversity anywhere on Earth. However, many regions of Africa remain under-sampled genetically. Here we present 350 whole genomes from Angola and Mozambique belonging to ten Bantu ethnolinguistic groups, enabling the construction of a reference variation catalogue including 2.9 million novel SNPs. We investigate the emergence of Bantu speaker population structure, admixture involving migrations across sub-Saharan Africa and model the demographic histories of Angolan and Mozambican Bantu speakers. Our results bring together concordant views from genomics, archaeology, and linguistics to paint an updated view of the complexity of the Bantu Expansion. Moreover, we generate reference panels that better represents the diversity of African populations involved in the trans-Atlantic slave trade, improving imputation accuracy in African Americans and Brazilians. We anticipate that our collection of genomes will form the foundation for future African genomic healthcare initiatives.

Haplotype-based inference of recent effective population size in modern and ancient DNA samples

Individuals sharing recent ancestors are likely to co-inherit large identical-by-descent (IBD) genomic regions. The distribution of these IBD segments in a population may be used to reconstruct past demographic events such as effective population size variation, but accurate IBD detection is difficult in ancient DNA data and in underrepresented populations with limited reference data. In this work, we introduce an accurate method for inferring effective population size variation during the past ~2000 years in both modern and ancient DNA data, called HapNe. HapNe infers recent population size fluctuations using either IBD sharing (HapNe-IBD) or linkage disequilibrium (HapNe-LD), which does not require phasing and can be computed in low coverage data, including data sets with heterogeneous sampling times. HapNe shows improved accuracy in a range of simulated demographic scenarios compared to currently available methods for IBD-based and LD-based inference of recent effective population size, while requiring fewer computational resources. We apply HapNe to several modern populations from the 1,000 Genomes Project, the UK Biobank, the Allen Ancient DNA Resource, and recently published samples from Iron Age Britain, detecting multiple instances of recent effective population size variation across these groups.

African American/Black race, apolipoprotein L1 , and serum creatinine among persons with HIV

OBJECTIVE

Accurate estimation of kidney function is critical among persons with HIV (PWH) to avoid under-dosing of antiretroviral therapies and ensure timely referral for kidney transplantation. Existing estimation equations for kidney function include race, the appropriateness of which has been debated. Given advancements in understanding of race and the necessity of accuracy in kidney function estimation, this study aimed to examine whether race, or genetic factors, improved prediction of serum creatinine among PWH.

DESIGN

This cross-sectional study utilized data from the Center for AIDS Research Network of Integrated Clinical Systems cohort (2008-2018). The outcome was baseline serum creatinine.

METHODS

Ordinary least squares regression was used to examine whether inclusion of race or genetic factors [ apolipoprotein-L1 ( APOL1 ) variants and genetic African ancestry] improved serum creatinine prediction. A reduction in root mean squared error (RMSE) greater than 2% was a clinically relevant improvement in predictive ability.

RESULTS

There were 4183 PWH included. Among PWH whose serum creatinine was less than 1.7 mg/dl, race was significantly associated with serum creatinine ( β  = 0.06, SE = 0.01, P  < 0.001) but did not improve predictive ability. African ancestry and APOL1 variants similarly failed to improve predictive ability. Whereas, when serum creatinine was at least 1.7 mg/dl, inclusion of race reduced the RMSE by 2.1%, indicating improvement in predictive ability. APOL1 variants further improved predictive ability by reducing the RMSE by 2.9%.

CONCLUSION

These data suggest that, among PWH, inclusion of race or genetic factors may only be warranted at higher serum creatinine levels. Work eliminating existing healthcare disparities while preserving the utility of estimating equations is needed.

Joint multi-ancestry and admixed GWAS reveals the complex genetics behind human cranial vault shape

The cranial vault in humans is highly variable, clinically relevant, and heritable, yet its genetic architecture remains poorly understood. Here, we conduct a joint multi-ancestry and admixed multivariate genome-wide association study on 3D cranial vault shape extracted from magnetic resonance images of 6772 children from the ABCD study cohort yielding 30 genome-wide significant loci. Follow-up analyses indicate that these loci overlap with genomic risk loci for sagittal craniosynostosis, show elevated activity cranial neural crest cells, are enriched for processes related to skeletal development, and are shared with the face and brain. We present supporting evidence of regional localization for several of the identified genes based on expression patterns in the cranial vault bones of E15.5 mice. Overall, our study provides a comprehensive overview of the genetics underlying normal-range cranial vault shape and its relevance for understanding modern human craniofacial diversity and the etiology of congenital malformations.

Genome-wide association analysis of plasma lipidome identifies 495 genetic associations

The human plasma lipidome captures risk for cardiometabolic diseases. To discover new lipid-associated variants and understand the link between lipid species and cardiometabolic disorders, we perform univariate and multivariate genome-wide analyses of 179 lipid species in 7174 Finnish individuals. We fine-map the associated loci, prioritize genes, and examine their disease links in 377,277 FinnGen participants. We identify 495 genome-trait associations in 56 genetic loci including 8 novel loci, with a considerable boost provided by the multivariate analysis. For 26 loci, fine-mapping identifies variants with a high causal probability, including 14 coding variants indicating likely causal genes. A phenome-wide analysis across 953 disease endpoints reveals disease associations for 40 lipid loci. For 11 coronary artery disease risk variants, we detect strong associations with lipid species. Our study demonstrates the power of multivariate genetic analysis in correlated lipidomics data and reveals genetic links between diseases and lipid species beyond the standard lipids.

Air Pollution and Allergic Rhinitis: Findings from a Prospective Cohort Study

To investigate the association of long-term exposure to ambient air pollution with the risk of allergic rhinitis (AR), we performed a longitudinal analysis of 379,488 participants (47.4% women) free of AR at baseline in the UK Biobank. The annual average concentrations of PM2.5, PMcoarse, PM10, NO2, and NOx were estimated by land use regression models. Cox proportional hazard models were used to calculate hazard ratios (HRs) and 95% confidence intervals (CIs). A weighted polygenic risk score was constructed. During a median follow-up period of 12.5 years, 3095 AR cases were identified. We observed significant associations between the risk of AR and PM2.5 (HR: 1.51, 95% CI: 1.27-1.79, per 5 μg/m3), PMcoarse (HR: 1.28, 95% CI: 1.06-1.55, per 5 μg/m3), PM10 (HR: 1.45, 95% CI: 1.20-1.74, per 10 μg/m3), NO2 (HR: 1.14, 95% CI: 1.09-1.19, per 10 μg/m3), and NOx (HR: 1.10, 95% CI: 1.05-1.15, per 20 μg/m3). Moreover, participants with high air pollution combined with high genetic risk showed the highest risk of AR, although no multiplicative or additive interaction was observed. In conclusion, long-term exposure to air pollutants was associated with an elevated risk of AR, particularly in high-genetic-risk populations, emphasizing the urgent need to improve air quality.

Assessing the associations between known genetic variants and substance use in people with HIV in the United States

BACKGROUND

The prevalence of substance use in people with HIV (PWH) in the United States is higher than in the general population and is an important driver of HIV-related outcomes. We sought to assess if previously identified genetic associations that contribute to substance use are also observed in a population of PWH.

METHODS

We performed genome-wide association studies (GWAS) of alcohol, smoking, and cannabis use phenotypes in a multi-ancestry population of 7,542 PWH from the Center for AIDS Research Network of Integrated Clinical Systems (CNICS). We conducted multi-ancestry GWAS for individuals of African (n = 3,748), Admixed American (n = 1,334), and European (n = 2,460) ancestry. Phenotype data were self-reported and collected using patient reported outcomes (PROs) and three questions from AUDIT-C, an alcohol screening tool. We analyzed nine phenotypes: 1) frequency of alcohol consumption, 2) typical number of drinks on a day when drinking alcohol, 3) frequency of five or more alcoholic drinks in a 30-day period, 4) smoking initiation, 5) smoking cessation, 6) cigarettes per day, 7) cannabis use initiation, 8) cannabis use cessation, 9) frequency of cannabis use during the previous 30 days. For each phenotype we considered a) variants previously identified as associated with a substance use trait and b) novel associations.

RESULTS

We observed evidence for effects of previously reported single nucleotide polymorphisms (SNPs) related to alcohol (rs1229984, p = 0.001), tobacco (rs11783093, p = 2.22E-4), and cannabis use (rs2875907, p = 0.005). We also report two novel loci (19p13.2, p = 1.3E-8; and 20p11.21, p = 2.1E-8) associated with cannabis use cessation.

CONCLUSIONS

Our analyses contribute to understanding the genetic bases of substance use in a population with relatively higher rates of use compared to the general population.

Polygenic risk for alcohol consumption and multisite chronic pain: Associations with ad lib drinking behavior

Interrelations between alcohol use disorder and chronic pain have received increasing empirical attention, and several lines of evidence support the possibility of shared genetic liability. However, research on the genetic contributions to the component processes of these complex and potentially overlapping phenotypes remains scarce. The goal of the present study was to test polygenic risk scores (PRSs) for alcohol consumption and multisite chronic pain as predictors of ad lib drinking behavior during an experimental taste test. PRSs were calculated for 209 pain-free, moderate-to-heavy drinkers (57.9% male; 63.6% White). Among White participants, the alcohol and chronic pain PRSs showed nominally significant (ps < .05) positive associations with the volume of alcohol consumed and peak blood alcohol concentration (BAC), respectively. However, associations did not survive correction for multiple comparisons. When stratifying results by experimental condition (between-subjects design: no-pain vs. pain), the alcohol PRS was significantly and negatively associated with the volume of alcohol poured, consumed, and peak BAC among Black participants randomized to the no-pain condition (all false discovery rate [FDR]p < .05). Conversely, the chronic pain PRS was significantly and positively associated with study outcomes among White participants in both the no-pain (alcohol consumed; FDRp = .037) and pain conditions (peak BAC; FDRp = .017). These findings lend partial support to the assertion that alcohol consumption in the laboratory is reflective of drinking behavior in naturalistic settings. This was also the first study to use a pain-related PRS to predict alcohol outcomes, which may be indicative of shared etiology between base and target traits. (PsycInfo Database Record (c) 2023 APA, all rights reserved).

Mexican Biobank advances population and medical genomics of diverse ancestries

Latin America continues to be severely underrepresented in genomics research, and fine-scale genetic histories and complex trait architectures remain hidden owing to insufficient data1. To fill this gap, the Mexican Biobank project genotyped 6,057 individuals from 898 rural and urban localities across all 32 states in Mexico at a resolution of 1.8 million genome-wide markers with linked complex trait and disease information creating a valuable nationwide genotype-phenotype database. Here, using ancestry deconvolution and inference of identity-by-descent segments, we inferred ancestral population sizes across Mesoamerican regions over time, unravelling Indigenous, colonial and postcolonial demographic dynamics2-6. We observed variation in runs of homozygosity among genomic regions with different ancestries reflecting distinct demographic histories and, in turn, different distributions of rare deleterious variants. We conducted genome-wide association studies (GWAS) for 22 complex traits and found that several traits are better predicted using the Mexican Biobank GWAS compared to the UK Biobank GWAS7,8. We identified genetic and environmental factors associating with trait variation, such as the length of the genome in runs of homozygosity as a predictor for body mass index, triglycerides, glucose and height. This study provides insights into the genetic histories of individuals in Mexico and dissects their complex trait architectures, both crucial for making precision and preventive medicine initiatives accessible worldwide.

CCLHunter: An efficient toolkit for cancer cell line authentication

Cancer cell lines are essential in cancer research, yet accurate authentication of these cell lines can be challenging, particularly for consanguineous cell lines with close genetic similarities. We introduce a new Cancer Cell Line Hunter (CCLHunter) method to tackle this challenge. This approach utilizes the information of single nucleotide polymorphisms, expression profiles, and kindred topology to authenticate 1389 human cancer cell lines accurately. CCLHunter can precisely and efficiently authenticate cell lines from consanguineous lineages and those derived from other tissues of the same individual. Our evaluation results indicate that CCLHunter has a complete accuracy rate of 93.27%, with an accuracy of 89.28% even for consanguineous cell lines, outperforming existing methods. Additionally, we provide convenient access to CCLHunter through standalone software and a web server at https://ngdc.cncb.ac.cn/cclhunter.

Genetic Distinctness and Diversity of American Aberdeen Cattle Compared to Common Beef Breeds in the United States

American Aberdeen (AD) cattle in the USA descend from an Aberdeen Angus herd originally brought to the Trangie Agricultural Research Centre, New South Wales, AUS. Although put under specific selection pressure for yearling growth rate, AD remain genomically uncharacterized. The objective was to characterize the genetic diversity and structure of purebred and crossbred AD cattle relative to seven common USA beef breeds using available whole-genome SNP data. A total of 1140 animals consisting of 404 purebred (n = 8 types) and 736 admixed individuals (n = 10 types) was used. Genetic diversity metrics, an analysis of molecular variance, and a discriminant analysis of principal components were employed. When linkage disequilibrium was not accounted for, markers influenced basic diversity parameter estimates, especially for AD cattle. Even so, intrapopulation and interpopulation estimates separate AD cattle from other purebred types (e.g., Latter's pairwise FST ranged from 0.1129 to 0.2209), where AD cattle were less heterozygous and had lower allelic richness than other purebred types. The admixed AD-influenced cattle were intermediate to other admixed types for similar parameters. The diversity metrics separation and differences support strong artificial selection pressures during and after AD breed development, shaping the evolution of the breed and making them genomically distinct from similar breeds.

Recent natural selection conferred protection against schizophrenia by non-antagonistic pleiotropy

Schizophrenia is a debilitating psychiatric disorder associated with a reduced fertility and decreased life expectancy, yet common predisposing variation substantially contributes to the onset of the disorder, which poses an evolutionary paradox. Previous research has suggested balanced selection, a mechanism by which schizophrenia risk alleles could also provide advantages under certain environments, as a reliable explanation. However, recent studies have shown strong evidence against a positive selection of predisposing loci. Furthermore, evolutionary pressures on schizophrenia risk alleles could have changed throughout human history as new environments emerged. Here in this study, we used 1000 Genomes Project data to explore the relationship between schizophrenia predisposing loci and recent natural selection (RNS) signatures after the human diaspora out of Africa around 100,000 years ago on a genome-wide scale. We found evidence for significant enrichment of RNS markers in derived alleles arisen during human evolution conferring protection to schizophrenia. Moreover, both partitioned heritability and gene set enrichment analyses of mapped genes from schizophrenia predisposing loci subject to RNS revealed a lower involvement in brain and neuronal related functions compared to those not subject to RNS. Taken together, our results suggest non-antagonistic pleiotropy as a likely mechanism behind RNS that could explain the persistence of schizophrenia common predisposing variation in human populations due to its association to other non-psychiatric phenotypes.

GCPBayes pipeline: a tool for exploring pleiotropy at the gene level

Cross-phenotype association using gene-set analysis can help to detect pleiotropic genes and inform about common mechanisms between diseases. Although there are an increasing number of statistical methods for exploring pleiotropy, there is a lack of proper pipelines to apply gene-set analysis in this context and using genome-scale data in a reasonable running time. We designed a user-friendly pipeline to perform cross-phenotype gene-set analysis between two traits using GCPBayes, a method developed by our team. All analyses could be performed automatically by calling for different scripts in a simple way (using a Shiny app, Bash or R script). A Shiny application was also developed to create different plots to visualize outputs from GCPBayes. Finally, a comprehensive and step-by-step tutorial on how to use the pipeline is provided in our group's GitHub page. We illustrated the application on publicly available GWAS (genome-wide association studies) summary statistics data to identify breast cancer and ovarian cancer susceptibility genes. We have shown that the GCPBayes pipeline could extract pleiotropic genes previously mentioned in the literature, while it also provided new pleiotropic genes and regions that are worthwhile for further investigation. We have also provided some recommendations about parameter selection for decreasing computational time of GCPBayes on genome-scale data.

Genotyping and population characteristics of the China Kadoorie Biobank

The China Kadoorie Biobank (CKB) is a population-based prospective cohort of >512,000 adults recruited from 2004 to 2008 from 10 geographically diverse regions across China. Detailed data from questionnaires and physical measurements were collected at baseline, with additional measurements at three resurveys involving ∼5% of surviving participants. Analyses of genome-wide genotyping, for >100,000 participants using custom-designed Axiom arrays, reveal extensive relatedness, recent consanguinity, and signatures reflecting large-scale population movements from recent Chinese history. Systematic genome-wide association studies of incident disease, captured through electronic linkage to death and disease registries and to the national health insurance system, replicate established disease loci and identify 14 novel disease associations. Together with studies of candidate drug targets and disease risk factors and contributions to international genetics consortia, these demonstrate the breadth, depth, and quality of the CKB data. Ongoing high-throughput omics assays of collected biosamples and planned whole-genome sequencing will further enhance the scientific value of this biobank.

Studying the genetics of participation using footprints left on the ascertained genotypes

The trait of participating in a genetic study probably has a genetic component. Identifying this component is difficult as we cannot compare genetic information of participants with nonparticipants directly, the latter being unavailable. Here, we show that alleles that are more common in participants than nonparticipants would be further enriched in genetic segments shared by two related participants. Genome-wide analysis was performed by comparing allele frequencies in shared and not-shared genetic segments of first-degree relative pairs of the UK Biobank. In nonoverlapping samples, a polygenic score constructed from that analysis is significantly associated with educational attainment, body mass index and being invited to a dietary study. The estimated correlation between the genetic components underlying participation in UK Biobank and educational attainment is estimated to be 36.6%-substantial but far from total. Taking participation behaviour into account would improve the analyses of the study data, including those of health traits.

Assessment of bidirectional relationships between brain imaging-derived phenotypes and stroke: a Mendelian randomization study

BACKGROUND

Stroke is a major cause of mortality and long-term disability worldwide. Whether the associations between brain imaging-derived phenotypes (IDPs) and stroke are causal is uncertain.

METHODS

We performed two-sample bidirectional Mendelian randomization (MR) analyses to explore the causal associations between IDPs and stroke. Summary data of 587 brain IDPs (up to 33,224 individuals) from the UK Biobank and five stroke types (sample size range from 301,663 to 446,696, case number range from 5,386 to 40,585) from the MEGASTROKE consortium were used.

RESULTS

Forward MR indicated 14 IDPs belong to projection fibers or association fibers were associated with stroke. For example, higher genetically determined mean diffusivity (MD) in the right external capsule was causally associated with an increased risk of small vessel stroke (IVW OR = 2.76, 95% CI 2.07 to 3.68, P = 5.87 × 10^-12). Reverse MR indicated that genetically determined higher risk of any ischemic stroke was associated with increased isotropic or free water volume fraction (ISOVF) in body of corpus callosum (IVW β = 0.23, 95% CI 0.14 to 0.33, P = 3.22 × 10^-7). This IDP is a commissural fiber and it is not included in the IDPs identified by forward MR.

CONCLUSIONS

We identified 14 IDPs with statistically significant evidence of causal effects on stroke or stroke subtypes. We also identified potential causal effects of stroke on one IDP of commissural fiber. These findings might guide further work toward identifying preventative strategies at the brain imaging levels.

Genetic Ancestry Estimates within Dutch Family Units and Across Genotyping Arrays: Insights from Empirical Analysis Using Two Estimation Methods

Accurate inference of genetic ancestry is crucial for population-based association studies, accounting for population heterogeneity and structure. This study analyzes genome-wide SNP data from the Netherlands Twin Register to compare genetic ancestry estimates. The focus is on the comparison of ancestry estimates between family members and individuals genotyped on multiple arrays (Affymetrix 6.0, Affymetrix Axiom, and Illumina GSA). Two conventional methods, principal component analysis and ADMIXTURE, were implemented to estimate ancestry, each serving its specific purpose, rather than for direct comparison. The results reveal that as the degree of genetic relatedness decreases, the Euclidean distances of genetic ancestry estimates between family members significantly increase (empirical p < 0.001), regardless of the estimation method and genotyping array. Ancestry estimates among individuals genotyped on multiple arrays also show statistically significant differences (empirical p < 0.001). Additionally, this study investigates the relationship between the ancestry estimates of non-identical twin offspring with ancestrally diverse parents and those with ancestrally similar parents. The results indicate a statistically significant weak correlation between the variation in ancestry estimates among offspring and differences in ancestry estimates among parents (Spearman's rho: 0.07, p = 0.005). This study highlights the utility of current methods in inferring genetic ancestry, emphasizing the importance of reference population composition in determining ancestry estimates.

Depression pathophysiology, risk prediction of recurrence and comorbid psychiatric disorders using genome-wide analyses

Depression is a common psychiatric disorder and a leading cause of disability worldwide. Here we conducted a genome-wide association study meta-analysis of six datasets, including >1.3 million individuals (371,184 with depression) and identified 243 risk loci. Overall, 64 loci were new, including genes encoding glutamate and GABA receptors, which are targets for antidepressant drugs. Intersection with functional genomics data prioritized likely causal genes and revealed new enrichment of prenatal GABAergic neurons, astrocytes and oligodendrocyte lineages. We found depression to be highly polygenic, with ~11,700 variants explaining 90% of the single-nucleotide polymorphism heritability, estimating that >95% of risk variants for other psychiatric disorders (anxiety, schizophrenia, bipolar disorder and attention deficit hyperactivity disorder) were influencing depression risk when both concordant and discordant variants were considered, and nearly all depression risk variants influenced educational attainment. Additionally, depression genetic risk was associated with impaired complex cognition domains. We dissected the genetic and clinical heterogeneity, revealing distinct polygenic architectures across subgroups of depression and demonstrating significantly increased absolute risks for recurrence and psychiatric comorbidity among cases of depression with the highest polygenic burden, with considerable sex differences. The risks were up to 5- and 32-fold higher than cases with the lowest polygenic burden and the background population, respectively. These results deepen the understanding of the biology underlying depression, its disease progression and inform precision medicine approaches to treatment.

Cross-ancestry genome-wide association meta-analyses of hippocampal and subfield volumes

The hippocampus is critical for memory and cognition and neuropsychiatric disorders, and its subfields differ in architecture and function. Genome-wide association studies on hippocampal and subfield volumes are mainly conducted in European populations; however, other ancestral populations are under-represented. Here we conduct cross-ancestry genome-wide association meta-analyses in 65,791 individuals for hippocampal volume and 38,977 for subfield volumes, including 7,009 individuals of East Asian ancestry. We identify 339 variant-trait associations at P < 1.13 × 10^-9 for 44 hippocampal traits, including 23 new associations. Common genetic variants have similar effects on hippocampal traits across ancestries, although ancestry-specific associations exist. Cross-ancestry analysis improves the fine-mapping precision and the prediction performance of polygenic scores in under-represented populations. These genetic variants are enriched for Wnt signaling and neuron differentiation and affect cognition, emotion and neuropsychiatric disorders. These findings may provide insight into the genetic architectures of hippocampal and subfield volumes.

Assortative mating and parental genetic relatedness drive the pathogenicity of variably expressive variants

We examined more than 38,000 spouse pairs from four neurodevelopmental disease cohorts and the UK Biobank to identify phenotypic and genetic patterns in parents associated with neurodevelopmental disease risk in children. We identified correlations between six phenotypes in parents and children, including correlations of clinical diagnoses such as obsessive-compulsive disorder (R=0.31-0.49, p<0.001), and two measures of sub-clinical autism features in parents affecting several autism severity measures in children, such as bi-parental mean Social Responsiveness Scale (SRS) scores affecting proband SRS scores (regression coefficient=0.11, p=0.003). We further describe patterns of phenotypic and genetic similarity between spouses, where spouses show both within- and cross-disorder correlations for seven neurological and psychiatric phenotypes, including a within-disorder correlation for depression (R=0.25-0.72, p<0.001) and a cross-disorder correlation between schizophrenia and personality disorder (R=0.20-0.57, p<0.001). Further, these spouses with similar phenotypes were significantly correlated for rare variant burden (R=0.07-0.57, p<0.0001). We propose that assortative mating on these features may drive the increases in genetic risk over generations and the appearance of "genetic anticipation" associated with many variably expressive variants. We further identified parental relatedness as a risk factor for neurodevelopmental disorders through its inverse correlations with burden and pathogenicity of rare variants and propose that parental relatedness drives disease risk by increasing genome-wide homozygosity in children (R=0.09-0.30, p<0.001). Our results highlight the utility of assessing parent phenotypes and genotypes in predicting features in children carrying variably expressive variants and counseling families carrying these variants.

Exploring the relationship between socioeconomic deprivation index and Alzheimer's disease using summary-level data: From genetic correlation to causality

Patients with Alzheimer's disease (AD) are markedly increasing as population aging and no disease-modifying therapies are currently available for AD. Previous studies suggested a broad link between socioeconomic status and a variety of disorders, including mental illness and cognitive abilities. However, the association between socioeconomic deprivation and AD has been unknown. We here employed Townsend deprivation index (TDI) to explore such relation and found a positive genetic correlation (r̂_g=0.211, P = 8.00 × 10^-4) between the two traits with summary statistics data (N = 455,258 for TDI and N = 455,815 for AD). Then, we performed pleiotropy analysis at both variant and gene levels using a powerful method called PLACO and detected 87 distinct pleiotropic genes. Functional analysis demonstrated these genes were significantly enriched in pancreas, liver, heart, blood, brain, and muscle tissues. Using Mendelian randomization methods, we further found that one genetically predicted standard deviation elevation in TDI could lead to approximately 18.5% (95% confidence intervals 1.6- 38.2%, P = 0.031) increase of AD risk, and that the identified causal association was robust against used MR approaches, horizontal pleiotropy, and instrumental selection. Overall, this study provides deep insight into common genetic components underlying TDI and AD, and further reveals causal connection between them. It is also helpful to develop a more suitable plan for ameliorating inequities, hardship, and disadvantage, with the hope of improving health outcomes among economically disadvantaged people.

Characterization of genome-wide STR variation in 6487 human genomes

Short tandem repeats (STRs) are abundant and highly mutagenic in the human genome. Many STR loci have been associated with a range of human genetic disorders. However, most population-scale studies on STR variation in humans have focused on European ancestry cohorts or are limited by sequencing depth. Here, we depicted a comprehensive map of 366,013 polymorphic STRs (pSTRs) constructed from 6487 deeply sequenced genomes, comprising 3983 Chinese samples (~31.5x, NyuWa) and 2504 samples from the 1000 Genomes Project (~33.3x, 1KGP). We found that STR mutations were affected by motif length, chromosome context and epigenetic features. We identified 3273 and 1117 pSTRs whose repeat numbers were associated with gene expression and 3'UTR alternative polyadenylation, respectively. We also implemented population analysis, investigated population differentiated signatures, and genotyped 60 known disease-causing STRs. Overall, this study further extends the scale of STR variation in humans and propels our understanding of the semantics of STRs.

Localizing Post-Admixture Adaptive Variants with Object Detection on Ancestry-Painted Chromosomes

Gene flow between previously differentiated populations during the founding of an admixed or hybrid population has the potential to introduce adaptive alleles into the new population. If the adaptive allele is common in one source population, but not the other, then as the adaptive allele rises in frequency in the admixed population, genetic ancestry from the source containing the adaptive allele will increase nearby as well. Patterns of genetic ancestry have therefore been used to identify post-admixture positive selection in humans and other animals, including examples in immunity, metabolism, and animal coloration. A common method identifies regions of the genome that have local ancestry "outliers" compared with the distribution across the rest of the genome, considering each locus independently. However, we lack theoretical models for expected distributions of ancestry under various demographic scenarios, resulting in potential false positives and false negatives. Further, ancestry patterns between distant sites are often not independent. As a result, current methods tend to infer wide genomic regions containing many genes as under selection, limiting biological interpretation. Instead, we develop a deep learning object detection method applied to images generated from local ancestry-painted genomes. This approach preserves information from the surrounding genomic context and avoids potential pitfalls of user-defined summary statistics. We find the method is robust to a variety of demographic misspecifications using simulated data. Applied to human genotype data from Cabo Verde, we localize a known adaptive locus to a single narrow region compared with multiple or long windows obtained using two other ancestry-based methods.

Associations between genetic liabilities to smoking behavior and schizophrenia symptoms in patients with a psychotic disorder, their siblings and healthy controls

It is unknown how smoking behavior polygenic scores (PRS) relate to psychosis and psychotic symptoms. To elucidate this, genotype and phenotype data were collected from patients with schizophrenia, their unaffected siblings, and healthy controls in a six-year follow-up prospective cohort study. Associations between smoking behaviors, PRS and schizophrenia symptoms were explored using linear mixed-effect models. The mean number of cigarettes smoked per day were 18 for patients, 13 for siblings and 12 for controls. In the overall sample, PRSs-smoking initiation (i.e., ever smoking as a binary phenotype, PRS-SI) were positively associated with positive symptoms, negative symptoms, and depressive symptoms, whereas PRSs-AI (age at regular smoking initiation) were negatively associated with all symptom dimensions, with similar effect sizes. When considering groups separately, PRS were only associated with psychotic symptoms in siblings and controls. In conclusion, unaffected siblings show smoking behaviors at an intermediate level between patients and healthy controls. Additionally, PRS-SI and PRS-AI are associated with all symptom dimensions only in unaffected siblings and healthy controls, possibly owing to the dominant role of other (genetic) risk factors in patients. Future studies may examine mechanisms via which genetic risk for smoking affects mental health symptoms.

Estimating human mobility in Holocene Western Eurasia with large-scale ancient genomic data

The recent increase in openly available ancient human DNA samples allows for large-scale meta-analysis applications. Trans-generational past human mobility is one of the key aspects that ancient genomics can contribute to since changes in genetic ancestry-unlike cultural changes seen in the archaeological record-necessarily reflect movements of people. Here, we present an algorithm for spatiotemporal mapping of genetic profiles, which allow for direct estimates of past human mobility from large ancient genomic datasets. The key idea of the method is to derive a spatial probability surface of genetic similarity for each individual in its respective past. This is achieved by first creating an interpolated ancestry field through space and time based on multivariate statistics and Gaussian process regression and then using this field to map the ancient individuals into space according to their genetic profile. We apply this algorithm to a dataset of 3138 aDNA samples with genome-wide data from Western Eurasia in the last 10,000 y. Finally, we condense this sample-wise record with a simple summary statistic into a diachronic measure of mobility for subregions in Western, Central, and Southern Europe. For regions and periods with sufficient data coverage, our similarity surfaces and mobility estimates show general concordance with previous results and provide a meta-perspective of genetic changes and human mobility.

Causal relationships between migraine and microstructural white matter: a Mendelian randomization study

BACKGROUND

Migraine is a disabling neurological disorder with the pathophysiology yet to be understood. The microstructural alteration in brain white matter (WM) has been suggested to be related to migraine in recent studies, but these evidence are observational essentially and cannot infer a causal relationship. The present study aims to reveal the causal relationship between migraine and microstructural WM using genetic data and Mendelian randomization (MR).

METHODS

We collected the Genome-wide association study (GWAS) summary statistics of migraine (48,975 cases / 550,381 controls) and 360 WM imaging-derived phenotypes (IDPs) (31,356 samples) that were used to measure microstructural WM. Based on instrumental variables (IVs) selected from the GWAS summary statistics, we conducted bidirectional two-sample MR analyses to infer bidirectional causal associations between migraine and microstructural WM. In forward MR analysis, we inferred the causal effect of microstructural WM on migraine by reporting the odds ratio (OR) that quantified the risk change of migraine for per 1 standard deviation (SD) increase of IDPs. In reverse MR analysis, we inferred the causal effect of migraine on microstructural WM by reporting the β value that represented SDs of changes in IDPs were caused by migraine.

RESULTS

Three WM IDPs showed significant causal associations (p < 3.29 × 10- 4, Bonferroni correction) with migraine and were proved to be reliable via sensitivity analysis. The mode of anisotropy (MO) of left inferior fronto-occipital fasciculus (OR = 1.76, p = 6.46 × 10- 5) and orientation dispersion index (OD) of right posterior thalamic radiation (OR = 0.78, p = 1.86 × 10- 4) exerted significant causal effects on migraine. Migraine exerted a significant causal effect on the OD of left superior cerebellar peduncle (β = - 0.09, p = 2.78 × 10- 4).

CONCLUSIONS

Our findings provided genetic evidence for the causal relationships between migraine and microstructural WM, bringing new insights into brain structure for the development and experience of migraine.