An association test of the spatial distribution of rare missense variants within protein structures identifies Alzheimer's disease-related patterns.

More than 90% of genetic variants are rare in most modern sequencing studies, such as the Alzheimer's Disease Sequencing Project (ADSP) whole-exome sequencing (WES) data. Furthermore, 54% of the rare variants in ADSP WES are singletons. However, both single variant and unit-based tests are limited in their statistical power to detect an association between rare variants and phenotypes. To best use missense rare variants and investigate their biological effect, we examine their association with phenotypes in the context of protein structures. We developed a protein structure-based approach, protein optimized kernel evaluation of missense nucleotides (POKEMON), which evaluates rare missense variants based on their spatial distribution within a protein rather than their allele frequency. The hypothesis behind this test is that the three-dimensional spatial distribution of variants within a protein structure provides functional context to power an association test. POKEMON identified three candidate genes (TREM2, SORL1, and EXOC3L4) and another suggestive gene from the ADSP WES data. For TREM2 and SORL1, two known Alzheimer's disease (AD) genes, the signal from the spatial cluster is stable even if we exclude known AD risk variants, indicating the presence of additional low-frequency risk variants within these genes. EXOC3L4 is a novel AD risk gene that has a cluster of variants primarily shared by case subjects around the Sec6 domain. This cluster is also validated in an independent replication data set and a validation data set with a larger sample size.

Sex-specific genetic architecture of late-life memory performance.

BACKGROUND: Women demonstrate a memory advantage when cognitively healthy yet lose this advantage to men in Alzheimer's disease. However, the genetic underpinnings of this sex difference in memory performance remain unclear. METHODS: We conducted the largest sex-aware genetic study on late-life memory to date (Nmales  = 11,942; Nfemales  = 15,641). Leveraging harmonized memory composite scores from four cohorts of cognitive aging and AD, we performed sex-stratified and sex-interaction genome-wide association studies in 24,216 non-Hispanic White and 3367 non-Hispanic Black participants. RESULTS: We identified three sex-specific loci (rs67099044-CBLN2, rs719070-SCHIP1/IQCJ-SCHIP), including an X-chromosome locus (rs5935633-EGL6/TCEANC/OFD1), that associated with memory. Additionally, we identified heparan sulfate signaling as a sex-specific pathway and found sex-specific genetic correlations between memory and cardiovascular, immune, and education traits. DISCUSSION: This study showed memory is highly and comparably heritable across sexes, as well as highlighted novel sex-specific genes, pathways, and genetic correlations that related to late-life memory. HIGHLIGHTS: Demonstrated the heritable component of late-life memory is similar across sexes. Identified two genetic loci with a sex-interaction with baseline memory. Identified an X-chromosome locus associated with memory decline in females. Highlighted sex-specific candidate genes and pathways associated with memory. Revealed sex-specific shared genetic architecture between memory and complex traits.

Longitudinal change in memory performance as a strong endophenotype for Alzheimer's disease.

INTRODUCTION: Although large-scale genome-wide association studies (GWAS) have been conducted on AD, few have been conducted on continuous measures of memory performance and memory decline. METHODS: We conducted a cross-ancestry GWAS on memory performance (in 27,633 participants) and memory decline (in 22,365 participants; 129,201 observations) by leveraging harmonized cognitive data from four aging cohorts. RESULTS: We found high heritability for two ancestry backgrounds. Further, we found a novel ancestry locus for memory decline on chromosome 4 (rs6848524) and three loci in the non-Hispanic Black ancestry group for memory performance on chromosomes 2 (rs111471504), 7 (rs4142249), and 15 (rs74381744). In our gene-level analysis, we found novel genes for memory decline on chromosomes 1 (SLC25A44), 11 (BSX), and 15 (DPP8). Memory performance and memory decline shared genetic architecture with AD-related traits, neuropsychiatric traits, and autoimmune traits. DISCUSSION: We discovered several novel loci, genes, and genetic correlations associated with late-life memory performance and decline. HIGHLIGHTS: Late-life memory has high heritability that is similar across ancestries. We discovered four novel variants associated with late-life memory. We identified four novel genes associated with late-life memory. Late-life memory shares genetic architecture with psychiatric/autoimmune traits.

Genetic heterogeneity: Challenges, impacts, and methods through an associative lens.

Genetic heterogeneity describes the occurrence of the same or similar phenotypes through different genetic mechanisms in different individuals. Robustly characterizing and accounting for genetic heterogeneity is crucial to pursuing the goals of precision medicine, for discovering novel disease biomarkers, and for identifying targets for treatments. Failure to account for genetic heterogeneity may lead to missed associations and incorrect inferences. Thus, it is critical to review the impact of genetic heterogeneity on the design and analysis of population level genetic studies, aspects that are often overlooked in the literature. In this review, we first contextualize our approach to genetic heterogeneity by proposing a high-level categorization of heterogeneity into "feature," "outcome," and "associative" heterogeneity, drawing on perspectives from epidemiology and machine learning to illustrate distinctions between them. We highlight the unique nature of genetic heterogeneity as a heterogeneous pattern of association that warrants specific methodological considerations. We then focus on the challenges that preclude effective detection and characterization of genetic heterogeneity across a variety of epidemiological contexts. Finally, we discuss systems heterogeneity as an integrated approach to using genetic and other high-dimensional multi-omic data in complex disease research.

A haptoglobin (HP) structural variant alters the effect of APOE alleles on Alzheimer's disease.

BACKGROUND: Haptoglobin (HP) is an antioxidant of apolipoprotein E (APOE), and previous reports have shown HP binds with APOE and amyloid beta (Aß) to aid its clearance. A common structural variant of the HP gene distinguishes it into two alleles: HP1 and HP2. METHODS: HP genotypes were imputed in 29 cohorts from the Alzheimer's Disease Genetics Consortium (N = 20,512). Associations between the HP polymorphism and Alzheimer's disease (AD) risk and age of onset through APOE interactions were investigated using regression models. RESULTS: The HP polymorphism significantly impacts AD risk in European-descent individuals (and in meta-analysis with African-descent individuals) by modifying both the protective effect of APOE ε2 and the detrimental effect of APOE ε4. The effect is particularly significant among APOE ε4 carriers. DISCUSSION: The effect modification of APOE by HP suggests adjustment and/or stratification by HP genotype is warranted when APOE risk is considered. Our findings also provided directions for further investigations on potential mechanisms behind this association.

Quality control and integration of genotypes from two calling pipelines for whole genome sequence data in the Alzheimer's disease sequencing project.

The Alzheimer's Disease Sequencing Project (ADSP) performed whole genome sequencing (WGS) of 584 subjects from 111 multiplex families at three sequencing centers. Genotype calling of single nucleotide variants (SNVs) and insertion-deletion variants (indels) was performed centrally using GATK-HaplotypeCaller and Atlas V2. The ADSP Quality Control (QC) Working Group applied QC protocols to project-level variant call format files (VCFs) from each pipeline, and developed and implemented a novel protocol, termed "consensus calling," to combine genotype calls from both pipelines into a single high-quality set. QC was applied to autosomal bi-allelic SNVs and indels, and included pipeline-recommended QC filters, variant-level QC, and sample-level QC. Low-quality variants or genotypes were excluded, and sample outliers were noted. Quality was assessed by examining Mendelian inconsistencies (MIs) among 67 parent-offspring pairs, and MIs were used to establish additional genotype-specific filters for GATK calls. After QC, 578 subjects remained. Pipeline-specific QC excluded ~12.0% of GATK and 14.5% of Atlas SNVs. Between pipelines, ~91% of SNV genotypes across all QCed variants were concordant; 4.23% and 4.56% of genotypes were exclusive to Atlas or GATK, respectively; the remaining ~0.01% of discordant genotypes were excluded. For indels, variant-level QC excluded ~36.8% of GATK and 35.3% of Atlas indels. Between pipelines, ~55.6% of indel genotypes were concordant; while 10.3% and 28.3% were exclusive to Atlas or GATK, respectively; and ~0.29% of discordant genotypes were. The final WGS consensus dataset contains 27,896,774 SNVs and 3,133,926 indels and is publicly available.

Genome-wide association meta-analysis of neuropathologic features of Alzheimer's disease and related dementias.

Alzheimer's disease (AD) and related dementias are a major public health challenge and present a therapeutic imperative for which we need additional insight into molecular pathogenesis. We performed a genome-wide association study and analysis of known genetic risk loci for AD dementia using neuropathologic data from 4,914 brain autopsies. Neuropathologic data were used to define clinico-pathologic AD dementia or controls, assess core neuropathologic features of AD (neuritic plaques, NPs; neurofibrillary tangles, NFTs), and evaluate commonly co-morbid neuropathologic changes: cerebral amyloid angiopathy (CAA), Lewy body disease (LBD), hippocampal sclerosis of the elderly (HS), and vascular brain injury (VBI). Genome-wide significance was observed for clinico-pathologic AD dementia, NPs, NFTs, CAA, and LBD with a number of variants in and around the apolipoprotein E gene (APOE). GalNAc transferase 7 (GALNT7), ATP-Binding Cassette, Sub-Family G (WHITE), Member 1 (ABCG1), and an intergenic region on chromosome 9 were associated with NP score; and Potassium Large Conductance Calcium-Activated Channel, Subfamily M, Beta Member 2 (KCNMB2) was strongly associated with HS. Twelve of the 21 non-APOE genetic risk loci for clinically-defined AD dementia were confirmed in our clinico-pathologic sample: CR1, BIN1, CLU, MS4A6A, PICALM, ABCA7, CD33, PTK2B, SORL1, MEF2C, ZCWPW1, and CASS4 with 9 of these 12 loci showing larger odds ratio in the clinico-pathologic sample. Correlation of effect sizes for risk of AD dementia with effect size for NFTs or NPs showed positive correlation, while those for risk of VBI showed a moderate negative correlation. The other co-morbid neuropathologic features showed only nominal association with the known AD loci. Our results discovered new genetic associations with specific neuropathologic features and aligned known genetic risk for AD dementia with specific neuropathologic changes in the largest brain autopsy study of AD and related dementias.

Systems biology approach to late-onset Alzheimer's disease genome-wide association study identifies novel candidate genes validated using brain expression data and Caenorhabditis elegans experiments.

INTRODUCTION: We sought to determine whether a systems biology approach may identify novel late-onset Alzheimer's disease (LOAD) loci. METHODS: We performed gene-wide association analyses and integrated results with human protein-protein interaction data using network analyses. We performed functional validation on novel genes using a transgenic Caenorhabditis elegans Aß proteotoxicity model and evaluated novel genes using brain expression data from people with LOAD and other neurodegenerative conditions. RESULTS: We identified 13 novel candidate LOAD genes outside chromosome 19. Of those, RNA interference knockdowns of the C. elegans orthologs of UBC, NDUFS3, EGR1, and ATP5H were associated with Aß toxicity, and NDUFS3, SLC25A11, ATP5H, and APP were differentially expressed in the temporal cortex. DISCUSSION: Network analyses identified novel LOAD candidate genes. We demonstrated a functional role for four of these in a C. elegans model and found enrichment of differentially expressed genes in the temporal cortex.

Genomic variants, genes, and pathways of Alzheimer's disease: An overview.

Alzheimer's disease (AD) (MIM: 104300) is a highly heritable disease with great complexity in its genetic contributors, and represents the most common form of dementia. With the gradual aging of the world's population, leading to increased prevalence of AD, and the substantial cost of care for those afflicted, identifying the genetic causes of disease represents a critical effort in identifying therapeutic targets. Here we provide a comprehensive review of genomic studies of AD, from the earliest linkage studies identifying monogenic contributors to early-onset forms of AD to the genome-wide and rare variant association studies of recent years that are being used to characterize the mosaic of genetic contributors to late-onset AD (LOAD), and which have identified approximately ∼20 genes with common variants contributing to LOAD risk. In addition, we explore studies employing alternative approaches to identify genetic contributors to AD, including studies of AD-related phenotypes and multi-variant association studies such as pathway analyses. Finally, we introduce studies of next-generation sequencing, which have recently helped identify multiple low-frequency and rare variant contributors to AD, and discuss on-going efforts with next-generation sequencing studies to develop statistically well- powered and comprehensive genomic studies of AD. Through this review, we help uncover the many insights the genetics of AD have provided into the pathways and pathophysiology of AD. © 2016 Wiley Periodicals, Inc.

SUCLG2 identified as both a determinator of CSF Aß1-42 levels and an attenuator of cognitive decline in Alzheimer's disease.

Cerebrospinal fluid amyloid-beta 1-42 (Aß1-42) and phosphorylated Tau at position 181 (pTau181) are biomarkers of Alzheimer's disease (AD). We performed an analysis and meta-analysis of genome-wide association study data on Aß1-42 and pTau181 in AD dementia patients followed by independent replication. An association was found between Aß1-42 level and a single-nucleotide polymorphism in SUCLG2 (rs62256378) (P = 2.5×10(-12)). An interaction between APOE genotype and rs62256378 was detected (P = 9.5 × 10(-5)), with the strongest effect being observed in APOE-ε4 noncarriers. Clinically, rs62256378 was associated with rate of cognitive decline in AD dementia patients (P = 3.1 × 10(-3)). Functional microglia experiments showed that SUCLG2 was involved in clearance of Aß1-42.

Parkinsonism and distinct dementia patterns in a family with the MAPT R406W mutation.

BACKGROUND: The Arg406Trp (R406W) missense mutation in the microtubule-associated protein-tau gene (MAPT) is a known cause of early-onset dementia. Various dementia phenotypes have been described, including frontotemporal dementia (FTD), FTD with parkinsonism, and early-onset Alzheimer disease (EOAD)-like presentations. METHODS: Using whole-exome capture with subsequent sequencing, we identified the R406W mutation in a family with multiple individuals with clinically diagnosed EOAD, in a pattern suggesting autosomal dominant inheritance. We reevaluated all available family members clinically. RESULTS: Each of the affected individuals had a course meeting clinical criteria for EOAD. Two distinct disease trajectories were apparent: one rapidly progressive, and the other long and gradual. Four of five affected individuals also manifested parkinsonian symptoms. FTD features were not prominent and, when present, appeared only late in the course of dementia. CONCLUSIONS: The MAPT R406W mutation is associated with EOAD-like symptoms and parkinsonism without FTD, as well as distinct cognitive courses.

Brain and Blood Transcriptome-Wide Association Studies Identify Five Novel Genes Associated with Alzheimer's Disease.

INTRODUCTION: Transcriptome-wide Association Studies (TWAS) extend genome-wide association studies (GWAS) by integrating genetically-regulated gene expression models. We performed the most powerful AD-TWAS to date, using summary statistics from cis -eQTL meta-analyses and the largest clinically-adjudicated Alzheimer's Disease (AD) GWAS. METHODS: We implemented the OTTERS TWAS pipeline, leveraging cis -eQTL data from cortical brain tissue (MetaBrain; N=2,683) and blood (eQTLGen; N=31,684) to predict gene expression, then applied these models to AD-GWAS data (Cases=21,982; Controls=44,944). RESULTS: We identified and validated five novel gene associations in cortical brain tissue ( PRKAG1 , C3orf62 , LYSMD4 , ZNF439 , SLC11A2 ) and six genes proximal to known AD-related GWAS loci (Blood: MYBPC3 ; Brain: MTCH2 , CYB561 , MADD , PSMA5 , ANXA11 ). Further, using causal eQTL fine-mapping, we generated sparse models that retained the strength of the AD-TWAS association for MTCH2 , MADD , ZNF439 , CYB561 , and MYBPC3 . DISCUSSION: Our comprehensive AD-TWAS discovered new gene associations and provided insights into the functional relevance of previously associated variants.

An X Chromosome Transcriptome Wide Association Study Implicates ARMCX6 in Alzheimer's Disease.

Background: The X chromosome is often omitted in disease association studies despite containing thousands of genes that may provide insight into well-known sex differences in the risk of Alzheimer's disease (AD). Objective: To model the expression of X chromosome genes and evaluate their impact on AD risk in a sex-stratified manner. Methods: Using elastic net, we evaluated multiple modeling strategies in a set of 175 whole blood samples and 126 brain cortex samples, with whole genome sequencing and RNA-seq data. SNPs (MAF > 0.05) within the cis-regulatory window were used to train tissue-specific models of each gene. We apply the best models in both tissues to sex-stratified summary statistics from a meta-analysis of Alzheimer's Disease Genetics Consortium (ADGC) studies to identify AD-related genes on the X chromosome. Results: Across different model parameters, sample sex, and tissue types, we modeled the expression of 217 genes (95 genes in blood and 135 genes in brain cortex). The average model R2 was 0.12 (range from 0.03 to 0.34). We also compared sex-stratified and sex-combined models on the X chromosome. We further investigated genes that escaped X chromosome inactivation (XCI) to determine if their genetic regulation patterns were distinct. We found ten genes associated with AD at p < 0.05, with only ARMCX6 in female brain cortex (p = 0.008) nearing the significance threshold after adjusting for multiple testing (α = 0.002). Conclusions: We optimized the expression prediction of X chromosome genes, applied these models to sex-stratified AD GWAS summary statistics, and identified one putative AD risk gene, ARMCX6.

Human whole-exome genotype data for Alzheimer's disease.

The heterogeneity of the whole-exome sequencing (WES) data generation methods present a challenge to a joint analysis. Here we present a bioinformatics strategy for joint-calling 20,504 WES samples collected across nine studies and sequenced using ten capture kits in fourteen sequencing centers in the Alzheimer's Disease Sequencing Project. The joint-genotype called variant-called format (VCF) file contains only positions within the union of capture kits. The VCF was then processed specifically to account for the batch effects arising from the use of different capture kits from different studies. We identified 8.2 million autosomal variants. 96.82% of the variants are high-quality, and are located in 28,579 Ensembl transcripts. 41% of the variants are intronic and 1.8% of the variants are with CADD > 30, indicating they are of high predicted pathogenicity. Here we show our new strategy can generate high-quality data from processing these diversely generated WES samples. The improved ability to combine data sequenced in different batches benefits the whole genomics research community.

Genetic factors in nonsmokers with age-related macular degeneration revealed through genome-wide gene-environment interaction analysis.

Relatively little is known about the interaction between genes and environment in the complex etiology of age-related macular degeneration (AMD). This study aimed to identify novel factors associated with AMD by analyzing gene-smoking interactions in a genome-wide association study of 1207 AMD cases and 686 controls of Caucasian background with genotype data on 668,238 single nucleotide polymorphisms (SNPs) after quality control. Participants' history of smoking at least 100 cigarettes lifetime was determined by a self-administered questionnaire. SNP associations modeled the effect of the minor allele additively on AMD using logistic regression, with adjustment for age, sex, and ever/never smoking. Joint effects of SNPs and smoking were examined comparing a null model containing only age, sex, and smoking against an extended model including genotypic and interaction terms. Genome-wide significant main effects were detected at three known AMD loci: CFH (P = 7.51×10(-30) ), ARMS2 (P = 1.94×10(-23) ), and RDBP/CFB/C2 (P = 4.37×10(-10) ), while joint effects analysis revealed three genomic regions with P < 10(-5) . Analyses stratified by smoking found genetic associations largely restricted to nonsmokers, with one notable exception: the chromosome 18q22.1 intergenic SNP rs17073641 (between SERPINB8 and CDH7), more strongly associated in nonsmokers (OR = 0.57, P = 2.73 × 10(-5) ), with an inverse association among smokers (OR = 1.42, P = 0.00228), suggesting that smoking modifies the effect of some genetic polymorphisms on AMD risk.

Gene-environment interactions and obesity traits among postmenopausal African-American and Hispanic women in the Women's Health Initiative SHARe Study.

Genome-wide association studies (GWAS) of obesity measures have identified associations with single nucleotide polymorphisms (SNPs). However, no large-scale evaluation of gene-environment interactions has been performed. We conducted a search of gene-environment (G × E) interactions in post-menopausal African-American and Hispanic women from the Women's Health Initiative SNP Health Association Resource GWAS study. Single SNP linear regression on body mass index (BMI) and waist-to-hip circumference ratio (WHR) adjusted for multidimensional-scaling-derived axes of ancestry and age was run in race-stratified data with 871,512 SNPs available from African-Americans (N = 8,203) and 786,776 SNPs from Hispanics (N = 3,484). Tests of G × E interaction at all SNPs for recreational physical activity (m h/week), dietary energy intake (kcal/day), alcohol intake (categorical), cigarette smoking years, and cigarette smoking (ever vs. never) were run in African-Americans and Hispanics adjusted for ancestry and age at interview, followed by meta-analysis of G × E interaction terms. The strongest evidence for concordant G × E interactions in African-Americans and Hispanics was for smoking and marker rs10133840 (Q statistic P = 0.70, beta = -0.01, P = 3.81 × 10(-7)) with BMI as the outcome. The strongest evidence for G × E interaction within a cohort was in African-Americans with WHR as outcome for dietary energy intake and rs9557704 (SNP × kcal = -0.04, P = 2.17 × 10(-7)). No results exceeded the Bonferroni-corrected statistical significance threshold.

Association of SCARB1 variants with subclinical atherosclerosis and incident cardiovascular disease: the multi-ethnic study of atherosclerosis.

OBJECTIVE: We previously reported a statistically significant association of SCARB1 intronic single nucleotide polymorphism (SNP) rs10846744 with common carotid intimal-medial artery thickness in each of the 4 Multi-Ethnic Study of Atherosclerosis racial/ethnic groups (white, Chinese, black, and Hispanic). METHODS AND RESULTS: Using an expanded sample of 7936 Multi-Ethnic Study of Atherosclerosis participants, phenotyped for measures of subclinical atherosclerosis, incident myocardial infarction, and cardiovascular disease, and genotyped through the SNP Health Association Resource project, we have now examined the genetic association of these phenotypes with 126 genotyped and imputed SCARB1 SNPs. We also performed stratified analyses to examine whether SCARB1 SNP effects differed by sex. Our analysis of the full Multi-Ethnic Study of Atherosclerosis cohort provides strong evidence for the association of rs10846744 with common carotid intimal-medial thickness (P=1.04E-4 in combined analysis of all 4 Multi-Ethnic Study of Atherosclerosis racial/ethnic groups). In sex-stratified analysis, we observed statistically significant association of rs10846744 with incident cardiovascular disease events in males (P=0.01). Examining analytical results from the Myocardial Infarction Genetics Consortium for replication, we observed further support for the association of rs10846744 with myocardial infarction. CONCLUSIONS: The SCARB1 SNP, rs10846744, exerts a major effect on subclinical atherosclerosis and incident cardiovascular disease in humans.

Dementia revealed: novel chromosome 6 locus for late-onset Alzheimer disease provides genetic evidence for folate-pathway abnormalities.

Genome-wide association studies (GWAS) of late-onset Alzheimer disease (LOAD) have consistently observed strong evidence of association with polymorphisms in APOE. However, until recently, variants at few other loci with statistically significant associations have replicated across studies. The present study combines data on 483,399 single nucleotide polymorphisms (SNPs) from a previously reported GWAS of 492 LOAD cases and 496 controls and from an independent set of 439 LOAD cases and 608 controls to strengthen power to identify novel genetic association signals. Associations exceeding the experiment-wide significance threshold (alpha=1.03x10(-7)) were replicated in an additional 1,338 cases and 2,003 controls. As expected, these analyses unequivocally confirmed APOE's risk effect (rs2075650, P=1.9x10(-36)). Additionally, the SNP rs11754661 at 151.2 Mb of chromosome 6q25.1 in the gene MTHFD1L (which encodes the methylenetetrahydrofolate dehydrogenase (NADP+ dependent) 1-like protein) was significantly associated with LOAD (P=4.70x10(-8); Bonferroni-corrected P=0.022). Subsequent genotyping of SNPs in high linkage disequilibrium (r2>0.8) with rs11754661 identified statistically significant associations in multiple SNPs (rs803424, P=0.016; rs2073067, P=0.03; rs2072064, P=0.035), reducing the likelihood of association due to genotyping error. In the replication case-control set, we observed an association of rs11754661 in the same direction as the previous association at P=0.002 (P=1.90x10(-10) in combined analysis of discovery and replication sets), with associations of similar statistical significance at several adjacent SNPs (rs17349743, P=0.005; rs803422, P=0.004). In summary, we observed and replicated a novel statistically significant association in MTHFD1L, a gene involved in the tetrahydrofolate synthesis pathway. This finding is noteworthy, as MTHFD1L may play a role in the generation of methionine from homocysteine and influence homocysteine-related pathways and as levels of homocysteine are a significant risk factor for LOAD development.

Associations of Sex, Race, and Apolipoprotein E Alleles With Multiple Domains of Cognition Among Older Adults.

Importance: Sex differences are established in associations between apolipoprotein E (APOE) ε4 and cognitive impairment in Alzheimer disease (AD). However, it is unclear whether sex-specific cognitive consequences of APOE are consistent across races and extend to the APOE ε2 allele. Objective: To investigate whether sex and race modify APOE ε4 and ε2 associations with cognition. Design, Setting, and Participants: This genetic association study included longitudinal cognitive data from 4 AD and cognitive aging cohorts. Participants were older than 60 years and self-identified as non-Hispanic White or non-Hispanic Black (hereafter, White and Black). Data were previously collected across multiple US locations from 1994 to 2018. Secondary analyses began December 2021 and ended September 2022. Main Outcomes and Measures: Harmonized composite scores for memory, executive function, and language were generated using psychometric approaches. Linear regression assessed interactions between APOE ε4 or APOE ε2 and sex on baseline cognitive scores, while linear mixed-effect models assessed interactions on cognitive trajectories. The intersectional effect of race was modeled using an APOE × sex × race interaction term, assessing whether APOE × sex interactions differed by race. Models were adjusted for age at baseline and corrected for multiple comparisons. Results: Of 32 427 participants who met inclusion criteria, there were 19 007 females (59%), 4453 Black individuals (14%), and 27 974 White individuals (86%); the mean (SD) age at baseline was 74 years (7.9). At baseline, 6048 individuals (19%) had AD, 4398 (14%) were APOE ε2 carriers, and 12 538 (38%) were APOE ε4 carriers. Participants missing APOE status were excluded (n = 9266). For APOE ε4, a robust sex interaction was observed on baseline memory (ß = -0.071, SE = 0.014; P = 9.6 × 10-7), whereby the APOE ε4 negative effect was stronger in females compared with males and did not significantly differ among races. Contrastingly, despite the large sample size, no APOE ε2 × sex interactions on cognition were observed among all participants. When testing for intersectional effects of sex, APOE ε2, and race, an interaction was revealed on baseline executive function among individuals who were cognitively unimpaired (ß = -0.165, SE = 0.066; P = .01), whereby the APOE ε2 protective effect was female-specific among White individuals but male-specific among Black individuals. Conclusions and Relevance: In this study, while race did not modify sex differences in APOE ε4, the APOE ε2 protective effect could vary by race and sex. Although female sex enhanced ε4-associated risk, there was no comparable sex difference in ε2, suggesting biological pathways underlying ε4-associated risk are distinct from ε2 and likely intersect with age-related changes in sex biology.

ODACH: a one-shot distributed algorithm for Cox model with heterogeneous multi-center data.

We developed a One-shot Distributed Algorithm for Cox proportional-hazards model to analyze Heterogeneous multi-center time-to-event data (ODACH) circumventing the need for sharing patient-level information across sites. This algorithm implements a surrogate likelihood function to approximate the Cox log-partial likelihood function that is stratified by site using patient-level data from a lead site and aggregated information from other sites, allowing the baseline hazard functions and the distribution of covariates to vary across sites. Simulation studies and application to a real-world opioid use disorder study showed that ODACH provides estimates close to the pooled estimator, which analyzes patient-level data directly from all sites via a stratified Cox model. Compared to the estimator from meta-analysis, the inverse variance-weighted average of the site-specific estimates, ODACH estimator demonstrates less susceptibility to bias, especially when the event is rare. ODACH is thus a valuable privacy-preserving and communication-efficient method for analyzing multi-center time-to-event data.