SCEBE: an efficient and scalable algorithm for genome-wide association studies on longitudinal outcomes with mixed-effects modeling

A framework for conducting GWAS using repeated measures data with an application to childhood BMI

Genetic effects on changes in human traits over time are understudied and may have important pathophysiological impact. We propose a framework that enables data quality control, implements mixed models to evaluate trajectories of change in traits, and estimates phenotypes to identify age-varying genetic effects in GWAS. Using childhood BMI as an example trait, we included 71,336 participants from six cohorts and estimated the slope and area under the BMI curve within four time periods (infancy, early childhood, late childhood and adolescence) for each participant, in addition to the age and BMI at the adiposity peak and the adiposity rebound. GWAS of the 12 estimated phenotypes identified 28 genome-wide significant variants at 13 loci, one of which (in DAOA) has not been previously associated with childhood or adult BMI. Genetic studies of changes in human traits over time could uncover unique biological mechanisms influencing quantitative traits.

Detecting time-varying genetic effects in Alzheimer's disease using a longitudinal genome-wide association studies model

INTRODUCTION

The development and progression of Alzheimer's disease (AD) is a complex process, during which genetic influences on phenotypes may also change. Incorporating longitudinal phenotypes in genome-wide association studies (GWAS) could unmask these genetic loci.

METHODS

We conducted a longitudinal GWAS using a varying coefficient test to identify age-dependent single nucleotide polymorphisms (SNPs) in AD. Data from 1877 Alzheimer's Neuroimaging Data Initiative participants, including impairment status and amyloid positron emission tomography (PET) scan standardized uptake value ratio (SUVR) scores, were analyzed using a retrospective varying coefficient mixed model association test (RVMMAT).

RESULTS

RVMMAT identified 244 SNPs with significant time-varying effects on AD impairment status, with 12 SNPs on chromosome 19 validated using National Alzheimer's Coordinating Center data. Age-stratified analyses showed these SNPs' effects peaked between 70 and 80 years. Additionally, 73 SNPs were linked to longitudinal amyloid accumulation changes. Pathway analyses implicated immune and neuroinflammation-related disruptions.

DISCUSSION

Our findings demonstrate that longitudinal GWAS models can uncover time-varying genetic signals in AD.

Highlights

Identify time-varying genetic effects using a longitudinal GWAS model in AD.Illustrate age-dependent genetic effects on both diagnoses and amyloid accumulation.Replicate time-varying effect of APOE in a second dataset.

A framework for conducting time-varying genome-wide association studies: An application to body mass index across childhood in six multiethnic cohorts

Genetic effects on changes in human traits over time are understudied and may have important pathophysiological impact. We propose a framework that enables data quality control, implements mixed models to evaluate trajectories of change in traits, and estimates phenotypes to identify age-varying genetic effects in genome-wide association studies (GWASs). Using childhood body mass index (BMI) as an example, we included 71,336 participants from six cohorts and estimated the slope and area under the BMI curve within four time periods (infancy, early childhood, late childhood and adolescence) for each participant, in addition to the age and BMI at the adiposity peak and the adiposity rebound. GWAS on each of the estimated phenotypes identified 28 genome-wide significant variants at 13 loci across the 12 estimated phenotypes, one of which was novel (in DAOA) and had not been previously associated with childhood or adult BMI. Genetic studies of changes in human traits over time could uncover novel biological mechanisms influencing quantitative traits.

Detecting time-varying genetic effects in Alzheimer's disease using a longitudinal GWAS model

Background

The development and progression of Alzheimer's disease (AD) is a complex process that can change over time, during which genetic influences on phenotypes may also fluctuate. Incorporating longitudinal phenotypes in genome wide association studies (GWAS) could help unmask genetic loci with time-varying effects. In this study, we incorporated a varying coefficient test in a longitudinal GWAS model to identify single nucleotide polymorphisms (SNPs) that may have time- or age-dependent effects in AD.

Methods

Genotype data from 1,877 participants in the Alzheimer's Neuroimaging Data Initiative (ADNI) were imputed using the Haplotype Reference Consortium (HRC) panel, resulting in 9,573,130 SNPs. Subjects' longitudinal impairment status at each visit was considered as a binary and clinical phenotype. Participants' composite standardized uptake value ratio (SUVR) derived from each longitudinal amyloid PET scan was considered as a continuous and biological phenotype. The retrospective varying coefficient mixed model association test (RVMMAT) was used in longitudinal GWAS to detect time-varying genetic effects on the impairment status and SUVR measures. Post-hoc analyses were performed on genome-wide significant SNPs, including 1) pathway analyses; 2) age-stratified genotypic comparisons and regression analyses; and 3) replication analyses using data from the National Alzheimer's Coordinating Center (NACC).

Results

Our model identified 244 genome-wide significant SNPs that revealed time-varying genetic effects on the clinical impairment status in AD; among which, 12 SNPs on chromosome 19 were successfully replicated using data from NACC. Post-hoc age-stratified analyses indicated that for most of these 244 SNPs, the maximum genotypic effect on impairment status occurred between 70 to 80 years old, and then declined with age. Our model further identified 73 genome-wide significant SNPs associated with the temporal variation of amyloid accumulation. For these SNPs, an increasing genotypic effect on PET-SUVR was observed as participants' age increased. Functional pathway analyses on significant SNPs for both phenotypes highlighted the involvement and disruption of immune responses- and neuroinflammation-related pathways in AD.

Conclusion

We demonstrate that longitudinal GWAS models with time-varying coefficients can boost the statistical power in AD-GWAS. In addition, our analyses uncovered potential time-varying genetic variants on repeated measurements of clinical and biological phenotypes in AD.

Genome-wide Analysis of Motor Progression in Parkinson Disease

Background and Objectives

The genetic basis of Parkinson disease (PD) motor progression is largely unknown. Previous studies of the genetics of PD progression have included small cohorts and shown a limited overlap with genetic PD risk factors from case-control studies. Here, we have studied genomic variation associated with PD motor severity and early-stage progression in large longitudinal cohorts to help to define the biology of PD progression and potential new drug targets.

Methods

We performed a GWAS meta-analysis of early PD motor severity and progression up to 3 years from study entry. We used linear mixed-effect models with additive effects, corrected for age at diagnosis, sex, and the first 5 genetic principal components to assess variability in axial, limb, and total Movement Disorder Society-Unified Parkinson's Disease Rating Scale (MDS-UPDRS) III scores.

Results

We included 3,572 unrelated European ancestry patients with PD from 5 observational cohorts and 1 drug trial. The average AAO was 62.6 years (SD = 9.83), and 63% of participants were male. We found an average increase in the total MDS-UPDRS III score of 2.3 points/year. We identified an association between PD axial motor progression and variation at the GJA5 locus at 1q12 (β = -0.25, SE = 0.04, p = 3.4e-10). Exploration of the regulation of gene expression in the region (cis-expression quantitative trait loci [eQTL] analysis) showed that the lead variant was associated with expression of ACP6, a lysophosphatidic acid phosphatase that regulates mitochondrial lipid biosynthesis (cis-eQTL p-values in blood and brain RNA expression data sets: <10-14 in eQTLGen and 10-7 in PsychEncode).

Discussion

Our study highlights the potential role of mitochondrial lipid homeostasis in the progression of PD, which may be important in establishing new drug targets that might modify disease progression.