Whole genome-wide sequence analysis of long-lived families (Long-Life Family Study) identifies MTUS2 gene associated with late-onset Alzheimer's disease.

INTRODUCTION: Late-onset Alzheimer's disease (LOAD) has a strong genetic component. Participants in Long-Life Family Study (LLFS) exhibit delayed onset of dementia, offering a unique opportunity to investigate LOAD genetics. METHODS: We conducted a whole genome sequence analysis of 3475 LLFS members. Genetic associations were examined in six independent studies (N = 14,260) with a wide range of LOAD risk. Association analysis in a sub-sample of the LLFS cohort (N = 1739) evaluated the association of LOAD variants with beta amyloid (Aß) levels. RESULTS: We identified several single nucleotide polymorphisms (SNPs) in tight linkage disequilibrium within the MTUS2 gene associated with LOAD (rs73154407, p = 7.6 × 10-9). Association of MTUS2 variants with LOAD was observed in the five independent studies and was significantly stronger within high levels of Aß42/40 ratio compared to lower amyloid. DISCUSSION: MTUS2 encodes a microtubule associated protein implicated in the development and function of the nervous system, making it a plausible candidate to investigate LOAD biology. HIGHLIGHTS: Long-Life Family Study (LLFS) families may harbor late onset Alzheimer's dementia (LOAD) variants. LLFS whole genome sequence analysis identified MTUS2 gene variants associated with LOAD. The observed LLFS variants generalized to cohorts with wide range of LOAD risk. The association of MTUS2 with LOAD was stronger within high levels of beta amyloid. Our results provide evidence for MTUS2 gene as a novel LOAD candidate locus.

Slower Decline in Processing Speed Is Associated with Familial Longevity.

INTRODUCTION: Cross-sectional analyses have associated familial longevity with better cognitive function and lower risk of cognitive impairment in comparison with individuals without familial longevity. The extent to which long-lived families also demonstrate slower rates of cognitive aging (i.e., change in cognition over time) is unknown. This study examined longitudinally collected data among 2 generations of the Long Life Family Study (LLFS) to compare rates of cognitive change across relatives and spouse controls. METHODS: We analyzed change in 6 neuropsychological test scores collected approximately 8 years apart among LLFS family members (n = 3,972) versus spouse controls (n = 1,092) using a Bayesian hierarchical model that included age, years of follow-up, sex, education, generation, and field center and all possible pairwise interactions. RESULTS: At a mean age of 88 years at enrollment in the older generation and 60 years in the younger generation, LLFS family members performed better than their spouses on the Digit Symbol Substitution Test (DSST) and the Logical Memory test. At follow-up, family members in the younger generation also showed slower decline than spouses on the DSST, whereas rates of change of Digit Span, fluency, and memory were similar between the 2 groups. DISCUSSION/CONCLUSION: Individuals in families with longevity appear to have better cognitive performance than their spouses for cognitive processes including psychomotor processing, episodic memory, and retrieval. Additionally, they demonstrate longer cognitive health spans with a slower decline on a multifactorial test of processing speed, a task requiring the integration of processes including organized visual search, working and incidental memory, and graphomotor ability. Long-lived families may be a valuable cohort for studying resilience to cognitive aging.

Association of fatigue, inflammation, and physical activity on gait speed: the Long Life Family Study.

BACKGROUND: Fatigue, inflammation, and physical activity (PA) are all independently associated with gait speed, but their directionality is not fully elucidated. AIMS: Evaluate the bidirectional associations amongst fatigue, inflammation, and PA on gait speed. METHODS: This cross sectional study included probands (n = 1280, aged 49-105) and offspring (n = 2772, aged 24-88) in the Long Life Family Study. We assessed gait speed, fatigue with the question "I could not get going", inflammation using fasting interleukin-6 (IL-6) and high sensitivity C-reactive protein (CRP), and self-reported PA as walking frequency in the past two weeks. The two generations were examined separately using linear mixed modeling. RESULTS: Lower fatigue, lower IL-6, and greater PA were all associated with faster gait speed in both generations (all p < 0.05); lower CRP was only associated with faster gait speed in the offspring. PA explained the association of fatigue and gait speed via a 16.1% (95% CI 9.7%, 26.7%) attenuation of the direct associations for the probands and 9.9% (95% CI 6.3%, 18.8%) in the offspring. In addition, IL-6 explained more of the association of fatigue and gait speed than the association between PA and gait speed, via a 14.9% (95% CI 9.2%, 23.4%) attenuation of the direct association in the offspring only. DISCUSSION: Results revealed a potential directionality from fatigue to IL-6 to PA that may lead to faster gait speed. Future work should examine these relationships longitudinally to establish temporality and causality. CONCLUSIONS: Our findings support a signal that lowering fatigue and inflammation and increasing physical activity may delay functional decline.

A Genome-Wide Association Study of 2304 Extreme Longevity Cases Identifies Novel Longevity Variants.

We performed a genome-wide association study (GWAS) of human extreme longevity (EL), defined as surviving past the 99th survival percentile, by aggregating data from four centenarian studies. The combined data included 2304 EL cases and 5879 controls. The analysis identified a locus in CDKN2B-AS1 (rs6475609, p = 7.13 × 10-8) that almost reached genome-wide significance and four additional loci that were suggestively significant. Among these, a novel rare variant (rs145265196) on chromosome 11 had much higher longevity allele frequencies in cases of Ashkenazi Jewish and Southern Italian ancestry compared to cases of other European ancestries. We also correlated EL-associated SNPs with serum proteins to link our findings to potential biological mechanisms that may be related to EL and are under genetic regulation. The findings from the proteomic analyses suggested that longevity-promoting alleles of significant genetic variants either provided EL cases with more youthful molecular profiles compared to controls or provided some form of protection from other illnesses, such as Alzheimer's disease, and disease progressions.

PopCluster: an algorithm to identify genetic variants with ethnicity-dependent effects.

MOTIVATION: Over the last decade, more diverse populations have been included in genome-wide association studies. If a genetic variant has a varying effect on a phenotype in different populations, genome-wide association studies applied to a dataset as a whole may not pinpoint such differences. It is especially important to be able to identify population-specific effects of genetic variants in studies that would eventually lead to development of diagnostic tests or drug discovery. RESULTS: In this paper, we propose PopCluster: an algorithm to automatically discover subsets of individuals in which the genetic effects of a variant are statistically different. PopCluster provides a simple framework to directly analyze genotype data without prior knowledge of subjects' ethnicities. PopCluster combines logistic regression modeling, principal component analysis, hierarchical clustering and a recursive bottom-up tree parsing procedure. The evaluation of PopCluster suggests that the algorithm has a stable low false positive rate (∼4%) and high true positive rate (>80%) in simulations with large differences in allele frequencies between cases and controls. Application of PopCluster to data from genetic studies of longevity discovers ethnicity-dependent heterogeneity in the association of rs3764814 (USP42) with the phenotype. AVAILABILITY AND IMPLEMENTATION: PopCluster was implemented using the R programming language, PLINK and Eigensoft software, and can be found at the following GitHub repository: https://github.com/gurinovich/PopCluster with instructions on its installation and usage. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Genetic Analysis Reveals a Longevity-Associated Protein Modulating Endothelial Function and Angiogenesis.

RATIONALE: Long living individuals show delay of aging, which is characterized by the progressive loss of cardiovascular homeostasis, along with reduced endothelial nitric oxide synthase activity, endothelial dysfunction, and impairment of tissue repair after ischemic injury. OBJECTIVE: Exploit genetic analysis of long living individuals to reveal master molecular regulators of physiological aging and new targets for treatment of cardiovascular disease. METHODS AND RESULTS: We show that the polymorphic variant rs2070325 (Ile229Val) in bactericidal/permeability-increasing fold-containing-family-B-member-4 (BPIFB4) associates with exceptional longevity, under a recessive genetic model, in 3 independent populations. Moreover, the expression of BPIFB4 is instrumental to maintenance of cellular and vascular homeostasis through regulation of protein synthesis. BPIFB4 phosphorylation/activation by protein-kinase-R-like endoplasmic reticulum kinase induces its complexing with 14-3-3 and heat shock protein 90, which is facilitated by the longevity-associated variant. In isolated vessels, BPIFB4 is upregulated by mechanical stress, and its knock-down inhibits endothelium-dependent vasorelaxation. In hypertensive rats and old mice, gene transfer of longevity-associated variant-BPIFB4 restores endothelial nitric oxide synthase signaling, rescues endothelial dysfunction, and reduces blood pressure levels. Furthermore, BPIFB4 is implicated in vascular repair. BPIFB4 is abundantly expressed in circulating CD34(+) cells of long living individuals, and its knock-down in endothelial progenitor cells precludes their capacity to migrate toward the chemoattractant SDF-1. In a murine model of peripheral ischemia, systemic gene therapy with longevity-associated variant-BPIFB4 promotes the recruitment of hematopoietic stem cells, reparative vascularization, and reperfusion of the ischemic muscle. CONCLUSIONS: Longevity-associated variant-BPIFB4 may represent a novel therapeutic tool to fight endothelial dysfunction and promote vascular reparative processes.

ICC-dementia (International Centenarian Consortium - dementia): an international consortium to determine the prevalence and incidence of dementia in centenarians across diverse ethnoracial and sociocultural groups.

BACKGROUND: Considerable variability exists in international prevalence and incidence estimates of dementia. The accuracy of estimates of dementia in the oldest-old and the controversial question of whether dementia incidence and prevalence decline at very old age will be crucial for better understanding the dynamics between survival to extreme old age and the occurrence and risk for various types of dementia and comorbidities. International Centenarian Consortium - Dementia (ICC-Dementia) seeks to harmonise centenarian and near-centenarian studies internationally to describe the cognitive and functional profiles of exceptionally old individuals, and ascertain the trajectories of decline and thereby the age-standardised prevalence and incidence of dementia in this population. The primary goal of the ICC-Dementia is to establish a large and thorough heterogeneous sample that has the power to answer epidemiological questions that small, separate studies cannot. A secondary aim is to examine cohort-specific effects and differential survivorship into very old age. We hope to lay the foundation for further investigation into risk and protective factors for dementia and healthy exceptional brain ageing in centenarians across diverse ethnoracial and sociocultural groups. METHODS: Studies focusing on individuals aged ≥95 years (approximately the oldest 1 percentile for men, oldest 5th percentile for women), with a minimum sample of 80 individuals, including assessment of cognition and functional status, are invited to participate. There are currently seventeen member or potential member studies from Asia, Europe, the Americas, and Oceania. Initial attempts at harmonising key variables are in progress. DISCUSSION: General challenges facing large, international consortia like ICC-Dementia include timely and effective communication among member studies, ethical and practical issues relating to human subject studies and data sharing, and the challenges related to data harmonisation. A specific challenge for ICC-Dementia relates to the concept and definition of'abnormal' in this exceptional group of individuals who are rarely free of physical, sensory and/or cognitive impairments.

A longevity-specific bank of induced pluripotent stem cells from centenarians and their offspring.

Centenarians provide a unique lens through which to study longevity, healthy aging, and resiliency. Moreover, models of human aging and resilience to disease that allow for the testing of potential interventions are virtually non-existent. We obtained and characterized over 50 centenarian and offspring peripheral blood samples including those connected to functional independence data highlighting resistance to disability and cognitive impairment. Targeted methylation arrays were used in molecular aging clocks to compare and contrast differences between biological and chronological age in these specialized subjects. Isolated peripheral blood mononuclear cells (PBMCs) were then successfully reprogrammed into high-quality induced pluripotent stem cell (iPSC) lines which were functionally characterized for pluripotency, genomic stability, and the ability to undergo directed differentiation. The result of this work is a one-of-a-kind resource for studies of human longevity and resilience that can fuel the discovery and validation of novel therapeutics for aging-related disease.

The relationship between 11 different polygenic longevity scores, parental lifespan, and disease diagnosis in the UK Biobank.

Large-scale genome-wide association studies (GWAS) strongly suggest that most traits and diseases have a polygenic component. This observation has motivated the development of disease-specific "polygenic scores (PGS)" that are weighted sums of the effects of disease-associated variants identified from GWAS that correlate with an individual's likelihood of expressing a specific phenotype. Although most GWAS have been pursued on disease traits, leading to the creation of refined "Polygenic Risk Scores" (PRS) that quantify risk to diseases, many GWAS have also been pursued on extreme human longevity, general fitness, health span, and other health-positive traits. These GWAS have discovered many genetic variants seemingly protective from disease and are often different from disease-associated variants (i.e., they are not just alternative alleles at disease-associated loci) and suggest that many health-positive traits also have a polygenic basis. This observation has led to an interest in "polygenic longevity scores (PLS)" that quantify the "risk" or genetic predisposition of an individual towards health. We derived 11 different PLS from 4 different available GWAS on lifespan and then investigated the properties of these PLS using data from the UK Biobank (UKB). Tests of association between the PLS and population structure, parental lifespan, and several cancerous and non-cancerous diseases, including death from COVID-19, were performed. Based on the results of our analyses, we argue that PLS are made up of variants not only robustly associated with parental lifespan, but that also contribute to the genetic architecture of disease susceptibility, morbidity, and mortality.

A Novel Gene ARHGAP44 for Longitudinal Changes in Glycated Hemoglobin (HbA1c) in Subjects without Type 2 Diabetes: Evidence from the Long Life Family Study (LLFS) and the Framingham Offspring Study (FOS).

Glycated hemoglobin (HbA1c) indicates average glucose levels over three months and is associated with insulin resistance and type 2 diabetes (T2D). Longitudinal changes in HbA1c (ΔHbA1c) are also associated with aging processes, cognitive performance, and mortality. We analyzed ΔHbA1c in 1,886 non-diabetic Europeans from the Long Life Family Study to uncover gene variants influencing ΔHbA1c. Using growth curve modeling adjusted for multiple covariates, we derived ΔHbA1c and conducted linkage-guided sequence analysis. Our genome-wide linkage scan identified a significant locus on 17p12. In-depth analysis of this locus revealed a variant rs56340929 (explaining 27% of the linkage peak) in the ARHGAP44 gene that was significantly associated with ΔHbA1c. RNA transcription of ARHGAP44 was associated with ΔHbA1c. The Framingham Offspring Study data further supported these findings on the gene level. Together, we found a novel gene ARHGAP44 for ΔHbA1c in family members without T2D. Follow-up studies using longitudinal omics data in large independent cohorts are warranted.

Cross-platform proteomics signatures of extreme old age.

In previous work we used a Somalogic platform targeting approximately 5000 proteins to generate a serum protein signature of centenarians that we validated in independent studies that used the same technology. We set here to validate and possibly expand the results by profiling the serum proteome of a subset of individuals included in the original study using liquid chromatography tandem mass spectrometry (LC-MS/MS). Following pre-processing, the LC-MS/MS data provided quantification of 398 proteins, with only 266 proteins shared by both platforms. At 1% FDR statistical significance threshold, the analysis of LC-MS/MS data detected 44 proteins associated with extreme old age, including 23 of the original analysis. To identify proteins for which associations between expression and extreme-old age were conserved across platforms, we performed inter-study conservation testing of the 266 proteins quantified by both platforms using a method that accounts for the correlation between the results. From these tests, a total of 80 proteins reached 5% FDR statistical significance, and 26 of these proteins had concordant pattern of gene expression in whole blood. This signature of 80 proteins points to blood coagulation, IGF signaling, extracellular matrix (ECM) organization, and complement cascade as important pathways whose protein level changes provide evidence for age-related adjustments that distinguish centenarians from younger individuals.

Signatures of Neuropsychological Test Results in the Long Life Family Study: A Cluster Analysis.

BACKGROUND: Discovering patterns of cognitive domains and characterizing how these patterns associate with other risk factors and biomarkers can improve our understanding of the determinants of cognitive aging. OBJECTIVE: To discover patterns of cognitive domains using neuropsychological test results in Long Life Family Study (LLFS) and characterize how these patterns associate with aging markers. METHODS: 5,086 LLFS participants were administered neuropsychological tests at enrollment. We performed a cluster analysis of six baseline neuropsychological test scores and tested the association between the identified clusters and various clinical variables, biomarkers, and polygenic risk scores using generalized estimating equations and the Chi-square test. We used Cox regression to correlate the clusters with the hazard of various medical events. We investigated whether the cluster information could enhance the prediction of cognitive decline using Bayesian beta regression. RESULTS: We identified 12 clusters with different cognitive signatures that represent profiles of performance across multiple neuropsychological tests. These signatures significantly correlated with 26 variables including polygenic risk scores, physical and pulmonary functions, and blood biomarkers and were associated with the hazard of mortality (p < 0.01), cardiovascular disease (p = 0.03), dementia (p = 0.01), and skin cancer (p = 0.03). CONCLUSION: The identified cognitive signatures capture multiple domains simultaneously and provide a holistic vision of cognitive function, showing that different patterns of cognitive function can coexist in aging individuals. Such patterns can be used for clinical intervention and primary care.

Rare genetic variants correlate with better processing speed.

We conducted a genome-wide association study of Digit Symbol Substitution Test scores administered in 4207 family members of the Long Life Family Study (LLFS). Genotype data were imputed to the HRC panel of 64,940 haplotypes resulting in ∼15M genetic variants with a quality score > 0.7. The results were replicated using genetic data imputed to the 1000 Genomes phase 3 reference panel from 2 Danish twin cohorts: the study of Middle Aged Danish Twins and the Longitudinal Study of Aging Danish Twins. The genome-wide association study in LLFS discovered 18 rare genetic variants (minor allele frequency (MAF) < 1.0%) that reached genome-wide significance (p-value < 5 × 10-8). Among these, 17 rare variants in chromosome 3 had large protective effects on the processing speed, including rs7623455, rs9821776, rs9821587, rs78704059, which were replicated in the combined Danish twin cohort. These SNPs are located in/near 2 genes, THRB and RARB, that belonged to the thyroid hormone receptors family that may influence the speed of metabolism and cognitive aging. The gene-level tests in LLFS confirmed that these 2 genes are associated with processing speed.

Mosaic Chromosomal Alterations and Human Longevity.

Mosaic chromosomal alterations (mCAs) are structural alterations associated with aging, cancer, cardiovascular disease, infectious diseases, and mortality. The distribution of mCAs in centenarians and individuals with familial longevity is poorly understood. We used MOsaic CHromosomal Alteration (MoChA) to discover mCAs in 2050 centenarians, offspring, and 248 controls from the New England Centenarian Study (NECS) and in 3 642 subjects with familial longevity and 920 spousal controls from the Long-Life Family Study (LLFS). We analyzed study-specific associations of somatic mCAs with age, familial longevity, the incidence of age-related diseases, and mortality and aggregated the results by meta-analysis. We show that the accumulation of mCAs > 100 KB increased to 102 years and plateaued at older ages. Centenarians and offspring accumulated fewer autosomal mCAs compared with controls (relative risk 0.637, p = .0147). Subjects with the APOE E4 allele had a 35.3% higher risk of accumulating autosomal mCAs (p = .002). Males were at higher risk for mCAs compared to females (male relative risk 1.36, p = 5.15e-05). mCAs were associated with increased hazard for cancer (hazard ratio 1.2) and dementia (hazard ratio 1.259) at a 10% false discovery rate. We observed a borderline significant association between mCAs and risk for mortality (hazard ratio 1.07, p = .0605). Our results show that the prevalence of individuals with mCAs does not continue to increase at ages >102 years and factors promoting familial longevity appear to confer protections from mCAs. These results suggest that limited mCA accumulation could be an important mechanism for extreme human longevity that needs to be investigated.

Multi-modal profiling of peripheral blood cells across the human lifespan reveals distinct immune cell signatures of aging and longevity.

BACKGROUND: Age-related changes in immune cell composition and functionality are associated with multimorbidity and mortality. However, many centenarians delay the onset of aging-related disease suggesting the presence of elite immunity that remains highly functional at extreme old age. METHODS: To identify immune-specific patterns of aging and extreme human longevity, we analyzed novel single cell profiles from the peripheral blood mononuclear cells (PBMCs) of a random sample of 7 centenarians (mean age 106) and publicly available single cell RNA-sequencing (scRNA-seq) datasets that included an additional 7 centenarians as well as 52 people at younger ages (20-89 years). FINDINGS: The analysis confirmed known shifts in the ratio of lymphocytes to myeloid cells, and noncytotoxic to cytotoxic cell distributions with aging, but also identified significant shifts from CD4+ T cell to B cell populations in centenarians suggesting a history of exposure to natural and environmental immunogens. We validated several of these findings using flow cytometry analysis of the same samples. Our transcriptional analysis identified cell type signatures specific to exceptional longevity that included genes with age-related changes (e.g., increased expression of STK17A, a gene known to be involved in DNA damage response) as well as genes expressed uniquely in centenarians' PBMCs (e.g., S100A4, part of the S100 protein family studied in age-related disease and connected to longevity and metabolic regulation). INTERPRETATION: Collectively, these data suggest that centenarians harbor unique, highly functional immune systems that have successfully adapted to a history of insults allowing for the achievement of exceptional longevity. FUNDING: TK, SM, PS, GM, SA, TP are supported by NIH-NIAUH2AG064704 and U19AG023122. MM and PS are supported by NIHNIA Pepper center: P30 AG031679-10. This project is supported by the Flow Cytometry Core Facility at BUSM. FCCF is funded by the NIH Instrumentation grant: S10 OD021587.

A metabolomic signature of the APOE2 allele.

With the goal of identifying metabolites that significantly correlate with the protective e2 allele of the apolipoprotein E (APOE) gene, we established a consortium of five studies of healthy aging and extreme human longevity with 3545 participants. This consortium includes the New England Centenarian Study, the Baltimore Longitudinal Study of Aging, the Arivale study, the Longevity Genes Project/LonGenity studies, and the Long Life Family Study. We analyzed the association between APOE genotype groups E2 (e2e2 and e2e3 genotypes, N = 544), E3 (e3e3 genotypes, N = 2299), and E4 (e3e4 and e4e4 genotypes, N = 702) with metabolite profiles in the five studies and used fixed effect meta-analysis to aggregate the results. Our meta-analysis identified a signature of 19 metabolites that are significantly associated with the E2 genotype group at FDR < 10%. The group includes 10 glycerolipids and 4 glycerophospholipids that were all higher in E2 carriers compared to E3, with fold change ranging from 1.08 to 1.25. The organic acid 6-hydroxyindole sulfate, previously linked to changes in gut microbiome that were reflective of healthy aging and longevity, was also higher in E2 carriers compared to E3 carriers. Three sterol lipids and one sphingolipid species were significantly lower in carriers of the E2 genotype group. For some of these metabolites, the effect of the E2 genotype opposed the age effect. No metabolites reached a statistically significant association with the E4 group. This work confirms and expands previous results connecting the APOE gene to lipid regulation and suggests new links between the e2 allele, lipid metabolism, aging, and the gut-brain axis.

Digitally generated Trail Making Test data: Analysis using hidden Markov modeling.

The Trail Making Test (TMT) is a neuropsychological test used to assess cognitive dysfunction. The TMT consists of two parts: TMT-A requires connecting numbers 1 to 25 sequentially; TMT-B requires connecting numbers 1 to 12 and letters A to L sequentially, alternating between numbers and letters. We propose using a digitally recorded version of TMT to capture cognitive or physical functions underlying test performance. We analyzed digital versions of TMT-A and -B to derive time metrics and used Bayesian hidden Markov models to extract additional metrics. We correlated these derived metrics with cognitive and physical function scores using regression. On both TMT-A and -B, digital metrics associated with graphomotor processing test scores and gait speed. Digital metrics on TMT-B were additionally associated with episodic memory test scores and grip strength. These metrics provide additional information of cognitive state and can differentiate cognitive and physical factors affecting test performance.

Association between late maternal age and age-related endophenotypes in the Long Life Family Study.

Extended maternal age has been suggested as marker of delayed age-associated disabilities. We use the Long Life Family Study (LLFS) offspring generation to investigate the association between extended maternal age at last childbirth and healthy-aging endophenotypes. We hypothesize that women with extended maternal age at last childbirth will exhibit healthier endophenotype profiles compared to younger mothers. The association between maternal age and age-related endophenotypes previously derived in LLFS was assessed using Generalized Estimating Equations to adjust for relatedness. The quartiles of the maternal age at last childbirth were modeled as the independent variables. Univariate analyses tested the association between maternal age at last childbirth and age at clinical assessment, education, field center, Apolipoprotein E (APOE) genotype, depression, stress, smoking and successful pregnancies. Only the variables significantly associated in the univariate analyses were considered in secondary multivariate analyses. Univariate analyses showed that compared to older mothers (age at last birth ≥35), mothers 30 years old or younger at last childbirth are less educated (12 ± 3 years versus 13 ± 3 years) and have a higher frequency of smoking (9% versus 3% for maternal age ≥35). Results showed that older mothers (age at last birth ≥31-34 or ≥ 35) demonstrated significantly better cognitive profiles (p = 0.017 and p = 0.021 respectively) compared with mothers with last childbirth age ≤30. Later maternal age among women from long-life families is associated with a better cognitive profile, supporting the hypothesis that later age at childbirth may be a marker for healthy aging.

Evaluation of GENESIS, SAIGE, REGENIE and fastGWA-GLMM for genome-wide association studies of binary traits in correlated data.

Performing a genome-wide association study (GWAS) with a binary phenotype using family data is a challenging task. Using linear mixed effects models is typically unsuitable for binary traits, and numerical approximations of the likelihood function may not work well with rare genetic variants with small counts. Additionally, imbalance in the case-control ratios poses challenges as traditional statistical methods such as the Score test or Wald test perform poorly in this setting. In the last couple of years, several methods have been proposed to better approximate the likelihood function of a mixed effects logistic regression model that uses Saddle Point Approximation (SPA). SPA adjustment has recently been implemented in multiple software, including GENESIS, SAIGE, REGENIE and fastGWA-GLMM: four increasingly popular tools to perform GWAS of binary traits. We compare Score and SPA tests using real family data to evaluate computational efficiency and the agreement of the results. Additionally, we compare various ways to adjust for family relatedness, such as sparse and full genetic relationship matrices (GRM) and polygenic effect estimates. We use the New England Centenarian Study imputed genotype data and the Long Life Family Study whole-genome sequencing data and the binary phenotype of human extreme longevity to compare the agreement of the results and tools' computational performance. The evaluation suggests that REGENIE might not be a good choice when analyzing correlated data of a small size. fastGWA-GLMM is the most computationally efficient compared to the other three tools, but it appears to be overly conservative when applied to family-based data. GENESIS, SAIGE and fastGWA-GLMM produced similar, although not identical, results, with SPA adjustment performing better than Score tests. Our evaluation also demonstrates the importance of adjusting by full GRM in highly correlated datasets when using GENESIS or SAIGE.