A statistical framework for powerful multi-trait rare variant analysis in large-scale whole-genome sequencing studies

Large-scale whole-genome sequencing (WGS) studies have improved our understanding of the contributions of coding and noncoding rare variants to complex human traits. Leveraging association effect sizes across multiple traits in WGS rare variant association analysis can improve statistical power over single-trait analysis, and also detect pleiotropic genes and regions. Existing multi-trait methods have limited ability to perform rare variant analysis of large-scale WGS data. We propose MultiSTAAR, a statistical framework and computationally-scalable analytical pipeline for functionally-informed multi-trait rare variant analysis in large-scale WGS studies. MultiSTAAR accounts for relatedness, population structure and correlation among phenotypes by jointly analyzing multiple traits, and further empowers rare variant association analysis by incorporating multiple functional annotations. We applied MultiSTAAR to jointly analyze three lipid traits (low-density lipoprotein cholesterol, high-density lipoprotein cholesterol and triglycerides) in 61,861 multi-ethnic samples from the Trans-Omics for Precision Medicine (TOPMed) Program. We discovered new associations with lipid traits missed by single-trait analysis, including rare variants within an enhancer of NIPSNAP3A and an intergenic region on chromosome 1.

Improving imputation quality in Samoans through the integration of population-specific sequences into existing reference panels

Genotype imputation is fundamental to association studies, and yet even gold standard panels like TOPMed are limited in the populations for which they yield good imputation. Specifically, Pacific Islanders are poorly represented in extant panels. To address this, we constructed an imputation reference panel using 1,285 Samoan individuals with whole-genome sequencing, combined with 1000 Genomes (1000G) samples, to create a reference panel that better represents Pacific Islander, specifically Samoan, genetic variation. We compared this panel to 1000G and TOPMed panels based on imputed variants using genotyping array data for 1,834 Samoan participants who were not part of the panels. The 1000G + 1285 Samoan panel yielded up to 2.25-2.76 times more well-imputed (r 2 ≥ 0.80) variants than TOPMed and 1000G. There was improved imputation accuracy across the minor allele frequency (MAF) spectrum, although it was more pronounced for variants with 0.01 ≤ MAF ≤ 0.05. Imputation accuracy (r 2 ) was greater for population-specific variants (high fixation index, F ST ) and those from larger haplotypes (high LD score). The gain in imputation accuracy over TOPMed was largest for small haplotypes (low LD score), reflecting the Samoan panel's ability to capture population-specific variation not well tagged by other panels. We also augmented the 1000G reference panel with varying numbers of Samoan samples and found that panels with 48 or more Samoans included outperformed TOPMed for all variants with MAF ≥ 0.001. This study identifies variants with improved imputation using population-specific reference panels and provides a framework for constructing other population-specific reference panels.

Rare variants in long non-coding RNAs are associated with blood lipid levels in the TOPMed whole-genome sequencing study

Long non-coding RNAs (lncRNAs) are known to perform important regulatory functions in lipid metabolism. Large-scale whole-genome sequencing (WGS) studies and new statistical methods for variant set tests now provide an opportunity to assess more associations between rare variants in lncRNA genes and complex traits across the genome. In this study, we used high-coverage WGS from 66,329 participants of diverse ancestries with measurement of blood lipids and lipoproteins (LDL-C, HDL-C, TC, and TG) in the National Heart, Lung, and Blood Institute (NHLBI) Trans-Omics for Precision Medicine (TOPMed) program to investigate the role of lncRNAs in lipid variability. We aggregated rare variants for 165,375 lncRNA genes based on their genomic locations and conducted rare-variant aggregate association tests using the STAAR (variant-set test for association using annotation information) framework. We performed STAAR conditional analysis adjusting for common variants in known lipid GWAS loci and rare-coding variants in nearby protein-coding genes. Our analyses revealed 83 rare lncRNA variant sets significantly associated with blood lipid levels, all of which were located in known lipid GWAS loci (in a ±500-kb window of a Global Lipids Genetics Consortium index variant). Notably, 61 out of 83 signals (73%) were conditionally independent of common regulatory variation and rare protein-coding variation at the same loci. We replicated 34 out of 61 (56%) conditionally independent associations using the independent UK Biobank WGS data. Our results expand the genetic architecture of blood lipids to rare variants in lncRNAs.

WHOLE GENOME SEQUENCING ANALYSIS OF BODY MASS INDEX IDENTIFIES NOVEL AFRICAN ANCESTRY-SPECIFIC RISK ALLELE

Obesity is a major public health crisis associated with high mortality rates. Previous genome-wide association studies (GWAS) investigating body mass index (BMI) have largely relied on imputed data from European individuals. This study leveraged whole-genome sequencing (WGS) data from 88,873 participants from the Trans-Omics for Precision Medicine (TOPMed) Program, of which 51% were of non-European population groups. We discovered 18 BMI-associated signals (P < 5 × 10-9). Notably, we identified and replicated a novel low frequency single nucleotide polymorphism (SNP) in MTMR3 that was common in individuals of African descent. Using a diverse study population, we further identified two novel secondary signals in known BMI loci and pinpointed two likely causal variants in the POC5 and DMD loci. Our work demonstrates the benefits of combining WGS and diverse cohorts in expanding current catalog of variants and genes confer risk for obesity, bringing us one step closer to personalized medicine.

A Polynesian-specific missense CETP variant alters the lipid profile

Identifying population-specific genetic variants associated with disease and disease-predisposing traits is important to provide insights into the genetic determinants of health and disease between populations, as well as furthering genomic justice. Various common pan-population polymorphisms at CETP associate with serum lipid profiles and cardiovascular disease. Here, sequencing of CETP identified a missense variant rs1597000001 (p.Pro177Leu) specific to Māori and Pacific people that associates with higher HDL-C and lower LDL-C levels. Each copy of the minor allele associated with higher HDL-C by 0.236 mmol/L and lower LDL-C by 0.133 mmol/L. The rs1597000001 effect on HDL-C is comparable with CETP Mendelian loss-of-function mutations that result in CETP deficiency, consistent with our data, which shows that rs1597000001 lowers CETP activity by 27.9%. This study highlights the potential of population-specific genetic analyses for improving equity in genomics and health outcomes for population groups underrepresented in genomic studies.

Association of rs9939609 in FTO with BMI among Polynesian peoples living in Aotearoa New Zealand and other Pacific nations

The fat mass and obesity associated (FTO) locus consistently associates with higher body mass index (BMI) across diverse ancestral groups. However, previous small studies of people of Polynesian ancestries have failed to replicate the association. In this study, we used Bayesian meta-analysis to test rs9939609, the most replicated FTO variant, for association with BMI with a large sample (n = 6095) of Aotearoa New Zealanders of Polynesian (Māori and Pacific) ancestry and of Samoan people living in the Independent State of Samoa and in American Samoa. We did not observe statistically significant association within each separate Polynesian subgroup. Bayesian meta-analysis of the Aotearoa New Zealand Polynesian and Samoan samples resulted in a posterior mean effect size estimate of +0.21 kg/m2, with a 95% credible interval [+0.03 kg/m2, +0.39 kg/m2]. While the Bayes Factor (BF) of 0.77 weakly favors the null, the BF = 1.4 Bayesian support interval is [+0.04, +0.20]. These results suggest that rs9939609 in FTO may have a similar effect on mean BMI in people of Polynesian ancestries as previously observed in other ancestral groups.

Rare variants in long non-coding RNAs are associated with blood lipid levels in the TOPMed Whole Genome Sequencing Study

Long non-coding RNAs (lncRNAs) are known to perform important regulatory functions. Large-scale whole genome sequencing (WGS) studies and new statistical methods for variant set tests now provide an opportunity to assess the associations between rare variants in lncRNA genes and complex traits across the genome. In this study, we used high-coverage WGS from 66,329 participants of diverse ancestries with blood lipid levels (LDL-C, HDL-C, TC, and TG) in the National Heart, Lung, and Blood Institute (NHLBI) Trans-Omics for Precision Medicine (TOPMed) program to investigate the role of lncRNAs in lipid variability. We aggregated rare variants for 165,375 lncRNA genes based on their genomic locations and conducted rare variant aggregate association tests using the STAAR (variant-Set Test for Association using Annotation infoRmation) framework. We performed STAAR conditional analysis adjusting for common variants in known lipid GWAS loci and rare coding variants in nearby protein coding genes. Our analyses revealed 83 rare lncRNA variant sets significantly associated with blood lipid levels, all of which were located in known lipid GWAS loci (in a ±500 kb window of a Global Lipids Genetics Consortium index variant). Notably, 61 out of 83 signals (73%) were conditionally independent of common regulatory variations and rare protein coding variations at the same loci. We replicated 34 out of 61 (56%) conditionally independent associations using the independent UK Biobank WGS data. Our results expand the genetic architecture of blood lipids to rare variants in lncRNA, implicating new therapeutic opportunities.

Validity of anthropometric equation-based estimators of fat mass in Samoan adults

INTRODUCTION

In 1999, a set of highly accurate Polynesian-specific equations to estimate adult body fat from non-invasive field measures of age, sex, height, and weight (Equation 1), age, sex, height, weight, and bioelectrical impedance analysis (BIA) resistance (Equation 2), and age, sex, height, weight, and the sum of two skinfold thicknesses (Equation 3) were published. The purpose of this study was to evaluate the performance of the equation-based estimators in a sample of Samoan adults recruited 20 years later between 2017 and 2019.

METHODS

Age, sex, height, weight, BIA resistance, skinfold thickness, and fat mass as measured using dual energy x-ray absorptiometry (DXA) were available for 432 Samoan adults (mean age 50.9 years, 56% female) seen in 2017/2019. We compared equation-derived fat mass and DXA-derived fat mass using scatterplots and Pearson correlation coefficients. We then updated the equation coefficient estimates in a training set (2/3 of the sample) and evaluated the performance of the updated equations in a testing set (the remaining 1/3 of the sample).

RESULTS

Equation-derived fat mass was strongly correlated with DXA-derived fat mass for Equation (1) (r2  = 0.95, n = 432), Equation (2) (r2  = 0.97, n = 425), and Equation (3) (r2  = 0.95, n = 426). Updating the equation coefficient estimates resulted in mostly similar coefficients and nearly identical testing set performance for Equation (1) (r2  = 0.96, n = 153), Equation (2) (r2  = 0.98, n = 150), and Equation (3) (r2  = 0.96, n = 150).

CONCLUSIONS

The Polynesian-specific body fat estimation equations remained stable despite changing social and environmental factors and marked increase in obesity prevalence in Samoa.

Compound heterozygous splicing variants expand the genotypic spectrum of EMC1-related disorders

EMC1 encodes subunit 1 of the endoplasmic reticulum (ER) membrane protein complex (EMC), a transmembrane domain insertase involved in membrane protein biosynthesis. Variants in EMC1 are described as a cause of global developmental delay, hypotonia, cortical visual impairment, and commonly, cerebral atrophy on MRI scan. We report an individual with severe global developmental delay and progressive cerebellar atrophy in whom exome sequencing identified a heterozygous essential splice-site variant in intron-3 of EMC1 (NM_015047.3:c.287-1G>A). Whole genome sequencing (WGS) identified a deep intronic variant in intron-20 of EMC1 (NM_015047.3:c.2588-771C>G) that was poorly predicted by in silico programs to disrupt pre-mRNA splicing. Reverse Transcription-PCR (RT-PCR) revealed stochastic activation of a pseudo-exon associated with the c.2588-771C>G variant and mis-splicing arising from the c.287-1G>A variant. This case highlights the utility of WGS and RNA studies to identify and assess likely pathogenicity of deep intronic variants and expands the genotypic and phenotypic spectrum of EMC1-related disorders.

Multi-ancestry transcriptome-wide association analyses yield insights into tobacco use biology and drug repurposing

Most transcriptome-wide association studies (TWASs) so far focus on European ancestry and lack diversity. To overcome this limitation, we aggregated genome-wide association study (GWAS) summary statistics, whole-genome sequences and expression quantitative trait locus (eQTL) data from diverse ancestries. We developed a new approach, TESLA (multi-ancestry integrative study using an optimal linear combination of association statistics), to integrate an eQTL dataset with a multi-ancestry GWAS. By exploiting shared phenotypic effects between ancestries and accommodating potential effect heterogeneities, TESLA improves power over other TWAS methods. When applied to tobacco use phenotypes, TESLA identified 273 new genes, up to 55% more compared with alternative TWAS methods. These hits and subsequent fine mapping using TESLA point to target genes with biological relevance. In silico drug-repurposing analyses highlight several drugs with known efficacy, including dextromethorphan and galantamine, and new drugs such as muscle relaxants that may be repurposed for treating nicotine addiction.

Multivariate analysis of a missense variant in CREBRF reveals associations with measures of adiposity in people of Polynesian ancestries

The minor allele of rs373863828, a missense variant in CREB3 Regulatory Factor, is associated with several cardiometabolic phenotypes in Polynesian peoples. To better understand the variant, we tested the association of rs373863828 with a panel of correlated phenotypes (body mass index [BMI], weight, height, HDL cholesterol, triglycerides, and total cholesterol) using multivariate Bayesian association and network analyses in a Samoa cohort (n = 1632), Aotearoa New Zealand cohort (n = 1419), and combined cohort (n = 2976). An expanded set of phenotypes (adding estimated fat and fat-free mass, abdominal circumference, hip circumference, and abdominal-hip ratio) was tested in the Samoa cohort (n = 1496). In the Samoa cohort, we observed significant associations (log10 Bayes Factor [BF] ≥ 5.0) between rs373863828 and the overall phenotype panel (8.81), weight (8.30), and BMI (6.42). In the Aotearoa New Zealand cohort, we observed suggestive associations (1.5 < log10 BF < 5) between rs373863828 and the overall phenotype panel (4.60), weight (3.27), and BMI (1.80). In the combined cohort, we observed concordant signals with larger log10 BFs. In the Samoa-specific expanded phenotype analyses, we also observed significant associations between rs373863828 and fat mass (5.65), abdominal circumference (5.34), and hip circumference (5.09). Bayesian networks provided evidence for a direct association of rs373863828 with weight and indirect associations with height and BMI.

A stop-gain variant in BTNL9 is associated with atherogenic lipid profiles

Current understanding of lipid genetics has come mainly from studies in European-ancestry populations; limited effort has focused on Polynesian populations, whose unique population history and high prevalence of dyslipidemia may provide insight into the biological foundations of variation in lipid levels. Here, we performed an association study to fine map a suggestive association on 5q35 with high-density lipoprotein cholesterol (HDL-C) seen in Micronesian and Polynesian populations. Fine-mapping analyses in a cohort of 2,851 Samoan adults highlighted an association between a stop-gain variant (rs200884524; c.652C>T, p.R218∗; posterior probability = 0.9987) in BTNL9 and both lower HDL-C and greater triglycerides (TGs). Meta-analysis across this and several other cohorts of Polynesian ancestry from Samoa, American Samoa, and Aotearoa New Zealand confirmed the presence of this association (βHDL-C = -1.60 mg/dL, p HDL-C = 7.63 × 10-10; βTG = 12.00 mg/dL, p TG = 3.82 × 10-7). While this variant appears to be Polynesian specific, there is also evidence of association from other multiancestry analyses in this region. This work provides evidence of a previously unexplored contributor to the genetic architecture of lipid levels and underscores the importance of genetic analyses in understudied populations.

The protective effect of rs373863828 on type 2 diabetes does not operate through a body composition pathway in adult Samoans

OBJECTIVE

The aim of this study was to understand whether the paradoxical association of missense variant rs373863828 in CREB3 regulatory factor (CREBRF) with higher BMI but lower odds of diabetes is explained by either metabolically favorable body fat distribution or greater fat-free mass.

METHODS

This study explored the association of the minor allele with dual-energy x-ray absorptiometry-derived body composition in n = 421 Samoans and used path analysis to examine the mediating role of fat and fat-free mass on the relationship between rs373863828 and fasting glucose.

RESULTS

Among females, the rs373863828 minor A allele was associated with greater BMI. There was no association of genotype with percent body fat, visceral adiposity, or fat distribution in either sex. In both females and males, lean mass was greater with each A allele: 2.16 kg/copy (p = 0.0001) and 1.73 kg/copy (p = 0.02), respectively. Path analysis showed a direct negative effect of rs373863828 genotype on fasting glucose (p = 0.004) consistent with previous findings, but also an indirect positive effect on fasting glucose operating through fat-free mass (p = 0.027).

CONCLUSIONS

The protective effect of rs373863828 in CREBRF, common among Pacific Islanders, on type 2 diabetes does not operate through body composition. Rather, the variant's effects on body size/composition and fasting glucose likely operate via different, tissue-specific mechanisms.

Rare Variants in Genes Encoding Subunits of the Epithelial Na+ Channel Are Associated With Blood Pressure and Kidney Function

BACKGROUND

The epithelial Na+ channel (ENaC) is intrinsically linked to fluid volume homeostasis and blood pressure. Specific rare mutations in SCNN1A, SCNN1B, and SCNN1G, genes encoding the α, β, and γ subunits of ENaC, respectively, are associated with extreme blood pressure phenotypes. No associations between blood pressure and SCNN1D, which encodes the δ subunit of ENaC, have been reported. A small number of sequence variants in ENaC subunits have been reported to affect functional transport in vitro or blood pressure. The effects of the vast majority of rare and low-frequency ENaC variants on blood pressure are not known.

METHODS

We explored the association of low frequency and rare variants in the genes encoding ENaC subunits, with systolic blood pressure, diastolic blood pressure, mean arterial pressure, and pulse pressure. Using whole-genome sequencing data from 14 studies participating in the Trans-Omics in Precision Medicine Whole-Genome Sequencing Program, and sequence kernel association tests.

RESULTS

We found that variants in SCNN1A and SCNN1B were associated with diastolic blood pressure and mean arterial pressure (P<0.00625). Although SCNN1D is poorly expressed in human kidney tissue, SCNN1D variants were associated with systolic blood pressure, diastolic blood pressure, mean arterial pressure, and pulse pressure (P<0.00625). ENaC variants in 2 of the 4 subunits (SCNN1B and SCNN1D) were also associated with estimated glomerular filtration rate (P<0.00625), but not with stroke.

CONCLUSIONS

Our results suggest that variants in extrarenal ENaCs, in addition to ENaCs expressed in kidneys, influence blood pressure and kidney function.

Insights From a Large-Scale Whole-Genome Sequencing Study of Systolic Blood Pressure, Diastolic Blood Pressure, and Hypertension

BACKGROUND

The availability of whole-genome sequencing data in large studies has enabled the assessment of coding and noncoding variants across the allele frequency spectrum for their associations with blood pressure.

METHODS

We conducted a multiancestry whole-genome sequencing analysis of blood pressure among 51 456 Trans-Omics for Precision Medicine and Centers for Common Disease Genomics program participants (stage-1). Stage-2 analyses leveraged array data from UK Biobank (N=383 145), Million Veteran Program (N=318 891), and Reasons for Geographic and Racial Differences in Stroke (N=10 643) participants, along with whole-exome sequencing data from UK Biobank (N=199 631) participants.

RESULTS

Two blood pressure signals achieved genome-wide significance in meta-analyses of stage-1 and stage-2 single variant findings (P<5×10^-8). Among them, a rare intergenic variant at novel locus, LOC100506274, was associated with lower systolic blood pressure in stage-1 (beta [SE]=-32.6 [6.0]; P=4.99×10^-8) but not stage-2 analysis (P=0.11). Furthermore, a novel common variant at the known INSR locus was suggestively associated with diastolic blood pressure in stage-1 (beta [SE]=-0.36 [0.07]; P=4.18×10^-7) and attained genome-wide significance in stage-2 (beta [SE]=-0.29 [0.03]; P=7.28×10^-23). Nineteen additional signals suggestively associated with blood pressure in meta-analysis of single and aggregate rare variant findings (P<1×10^-6 and P<1×10^-4, respectively).

DISCUSSION

We report one promising but unconfirmed rare variant for blood pressure and, more importantly, contribute insights for future blood pressure sequencing studies. Our findings suggest promise of aggregate analyses to complement single variant analysis strategies and the need for larger, diverse samples, and family studies to enable robust rare variant identification.

Whole genome sequence association analysis of fasting glucose and fasting insulin levels in diverse cohorts from the NHLBI TOPMed program

The genetic determinants of fasting glucose (FG) and fasting insulin (FI) have been studied mostly through genome arrays, resulting in over 100 associated variants. We extended this work with high-coverage whole genome sequencing analyses from fifteen cohorts in NHLBI's Trans-Omics for Precision Medicine (TOPMed) program. Over 23,000 non-diabetic individuals from five race-ethnicities/populations (African, Asian, European, Hispanic and Samoan) were included. Eight variants were significantly associated with FG or FI across previously identified regions MTNR1B, G6PC2, GCK, GCKR and FOXA2. We additionally characterize suggestive associations with FG or FI near previously identified SLC30A8, TCF7L2, and ADCY5 regions as well as APOB, PTPRT, and ROBO1. Functional annotation resources including the Diabetes Epigenome Atlas were compiled for each signal (chromatin states, annotation principal components, and others) to elucidate variant-to-function hypotheses. We provide a catalog of nucleotide-resolution genomic variation spanning intergenic and intronic regions creating a foundation for future sequencing-based investigations of glycemic traits.

A loss-of-function IFNAR1 allele in Polynesia underlies severe viral diseases in homozygotes

Globally, autosomal recessive IFNAR1 deficiency is a rare inborn error of immunity underlying susceptibility to live attenuated vaccine and wild-type viruses. We report seven children from five unrelated kindreds of western Polynesian ancestry who suffered from severe viral diseases. All the patients are homozygous for the same nonsense IFNAR1 variant (p.Glu386*). This allele encodes a truncated protein that is absent from the cell surface and is loss-of-function. The fibroblasts of the patients do not respond to type I IFNs (IFN-α2, IFN-ω, or IFN-β). Remarkably, this IFNAR1 variant has a minor allele frequency >1% in Samoa and is also observed in the Cook, Society, Marquesas, and Austral islands, as well as Fiji, whereas it is extremely rare or absent in the other populations tested, including those of the Pacific region. Inherited IFNAR1 deficiency should be considered in individuals of Polynesian ancestry with severe viral illnesses.

The missense variant, rs373863828, in CREBRF plays a role in longitudinal changes in body mass index in Samoans

OBJECTIVE

A missense variant, rs373863828, in CREBRF is associated with obesity in Polynesians. We investigate whether rs373863828 and other factors are associated with body mass index (BMI) rate-of-change between 2010 and 2017-19 in Samoans.

METHODS

We used sex-stratified models to test whether BMI rate-of-change was associated with rs373863828, baseline BMI, age, residence, physical activity, and household asset score in a cohort study of 480 Samoan adults measured in both 2010 (mean age 43.8 years) and 2017-19.

RESULTS

Mean BMI increased from 32.1 to 33.5 kg/m² in males (n = 220, p = 1.3 ×10^-8) and from 35.9 to 37.8 kg/m² in females (n = 260, p = 1.2 ×10^-13). In females, the A allele was associated with a higher rate-of-change (0.150 kg/m²/year/allele, p = 1.7 ×10^-4). Across 10-year age groups, mean BMI rate-of-change was lower in older participants. The BMI rate of change differed by genotype: it was, in females with AA genotype, approximately half that seen in GG and AG participants. In females lower baseline household asset scores were associated with a higher rate-of-change (p = 0.002).

CONCLUSIONS

In Samoans, the minor A allele of rs373863828 is associated with an increased rate-of-change in BMI in females. On average, BMI of females with the AA genotype increased 0.30 kg/m²/year more than of those with the GG genotype.

Mendelian randomization supports bidirectional causality between telomere length and clonal hematopoiesis of indeterminate potential

Human genetic studies support an inverse causal relationship between leukocyte telomere length (LTL) and coronary artery disease (CAD), but directionally mixed effects for LTL and diverse malignancies. Clonal hematopoiesis of indeterminate potential (CHIP), characterized by expansion of hematopoietic cells bearing leukemogenic mutations, predisposes both hematologic malignancy and CAD. TERT (which encodes telomerase reverse transcriptase) is the most significantly associated germline locus for CHIP in genome-wide association studies. Here, we investigated the relationship between CHIP, LTL, and CAD in the Trans-Omics for Precision Medicine (TOPMed) program (n = 63,302) and UK Biobank (n = 47,080). Bidirectional Mendelian randomization studies were consistent with longer genetically imputed LTL increasing propensity to develop CHIP, but CHIP then, in turn, hastens to shorten measured LTL (mLTL). We also demonstrated evidence of modest mediation between CHIP and CAD by mLTL. Our data promote an understanding of potential causal relationships across CHIP and LTL toward prevention of CAD.

CREBRF missense variant rs373863828 has both direct and indirect effects on type 2 diabetes and fasting glucose in Polynesian peoples living in Samoa and Aotearoa New Zealand

INTRODUCTION

The minor allele of a missense variant, rs373863828, in CREBRF is associated with higher body mass index (BMI), lower fasting glucose, and lower odds of type 2 diabetes. rs373863828 is common in Pacific Island populations (minor allele frequency (MAF) 0.096-0.259) but rare in non-Pacific Island populations (MAF <0.001). We examined the cross-sectional associations between BMI and rs373863828 in type 2 diabetes and fasting glucose with a large sample of adults of Polynesian ancestries from Samoa, American Samoa, and Aotearoa New Zealand, and estimated the direct and indirect (via BMI) effects of rs373863828 on type 2 diabetes and fasting glucose.

RESEARCH DESIGN AND METHODS

We regressed type 2 diabetes and fasting glucose on BMI and rs373863828 stratified by obesity, regressed type 2 diabetes and fasting glucose on BMI stratified by rs373863828 genotype, and assessed the effects of rs373863828 on type 2 diabetes and fasting glucose with path analysis. The regression analyses were completed separately in four samples that were recruited during different time periods between 1990 and 2010 and then the results were meta-analyzed. All samples were pooled for the path analysis.

RESULTS

Association of BMI with type 2 diabetes and fasting glucose may be greater in those without obesity (OR=7.77, p=0.015 and β=0.213, p=9.53×10^-5, respectively) than in those with obesity (OR=5.01, p=1.12×10^-9 and β=0.162, p=5.63×10^-6, respectively). We did not observe evidence of differences in the association of BMI with type 2 diabetes or fasting glucose by genotype. In the path analysis, the minor allele has direct negative (lower odds of type 2 diabetes and fasting glucose) and indirect positive (higher odds of type 2 diabetes and fasting glucose) effects on type 2 diabetes risk and fasting glucose, with the indirect effects mediated through a direct positive effect of rs373863828 on BMI.

CONCLUSIONS

There may be a stronger effect of BMI on fasting glucose in Polynesian individuals without obesity than in those with obesity. Carrying the rs373863828 minor allele does not decouple higher BMI from higher odds of type 2 diabetes.

Genetic determinants of telomere length from 109,122 ancestrally diverse whole-genome sequences in TOPMed

Genetic studies on telomere length are important for understanding age-related diseases. Prior GWAS for leukocyte TL have been limited to European and Asian populations. Here, we report the first sequencing-based association study for TL across ancestrally-diverse individuals (European, African, Asian and Hispanic/Latino) from the NHLBI Trans-Omics for Precision Medicine (TOPMed) program. We used whole genome sequencing (WGS) of whole blood for variant genotype calling and the bioinformatic estimation of telomere length in n=109,122 individuals. We identified 59 sentinel variants (p-value <5×10-9) in 36 loci associated with telomere length, including 20 newly associated loci (13 were replicated in external datasets). There was little evidence of effect size heterogeneity across populations. Fine-mapping at OBFC1 indicated the independent signals colocalized with cell-type specific eQTLs for OBFC1 (STN1). Using a multi-variant gene-based approach, we identified two genes newly implicated in telomere length, DCLRE1B (SNM1B) and PARN. In PheWAS, we demonstrated our TL polygenic trait scores (PTS) were associated with increased risk of cancer-related phenotypes.

Iron homeostasis pathway DNA methylation trajectories reveal a role for STEAP3 metalloreductase in patient outcomes after aneurysmal subarachnoid hemorrhage

BACKGROUND

Following aneurysmal subarachnoid hemorrhage (aSAH), the brain is susceptible to ferroptosis, a type of iron-dependent cell death. Therapeutic intervention targeting the iron homeostasis pathway shows promise for mitigating ferroptosis and improving recovery in animal models, but little work has been conducted in humans. DNA methylation (DNAm) plays a key role in gene expression and brain function, plasticity, and injury recovery, making it a potentially useful biomarker of outcomes or therapeutic target for intervention. Therefore, in this longitudinal, observational study, we examined the relationships between trajectories of DNAm in candidate genes related to iron homeostasis and acute (cerebral vasospasm and delayed cerebral ischemia) and long-term (Glasgow Outcome Scale [GOS, unfavorable = 1-3] and death) patient outcomes after aSAH.

RESULTS

Longitudinal, genome-wide DNAm data were generated from DNA extracted from post-aSAH cerebrospinal fluid (n = 260 participants). DNAm trajectories of 637 CpG sites in 36 candidate genes related to iron homeostasis were characterized over 13 days post-aSAH using group-based trajectory analysis, an unsupervised clustering method. Significant associations were identified between inferred DNAm trajectory groups at several CpG sites and acute and long-term outcomes. Among our results, cg25713625 in the STEAP3 metalloreductase gene (STEAP3) stood out. Specifically, in comparing the highest cg25713625 DNAm trajectory group with the lowest, we observed significant associations (i.e., based on p-values less than an empirical significance threshold) with unfavorable GOS at 3 and 12 months (OR = 11.7, p = 0.0006 and OR = 15.6, p = 0.0018, respectively) and death at 3 and 12 months (OR = 19.1, p = 0.0093 and OR = 12.8, p = 0.0041, respectively). These results were replicated in an independent sample (n = 100 participants) observing significant associations with GOS at 3 and 12 months (OR = 8.2, p = 0.001 and OR = 6.3, p = 0.0.0047, respectively) and death at 3 months (OR = 2.3, p = 0.008) and a suggestive association (i.e., p-value < 0.05 not meeting an empirical significance threshold) with death at 12 months (OR = 2.0, p = 0.0272). In both samples, an additive effect of the DNAm trajectory group was observed as the percentage of participants with unfavorable long-term outcomes increased substantially with higher DNAm trajectory groups.

CONCLUSION

Our results support a role for DNAm of cg25713625/STEAP3 in recovery following aSAH. Additional research is needed to further explore the role of DNAm of cg25713625/STEAP3 as a biomarker of unfavorable outcomes, or therapeutic target to improve outcomes, to translate these findings clinically.

A murine model of the human CREBRFR457Q obesity-risk variant does not influence energy or glucose homeostasis in response to nutritional stress

Obesity and diabetes have strong heritable components, yet the genetic contributions to these diseases remain largely unexplained. In humans, a missense variant in Creb3 regulatory factor (CREBRF) [rs373863828 (p.Arg457Gln); CREBRFR457Q] is strongly associated with increased odds of obesity but decreased odds of diabetes. Although virtually nothing is known about CREBRF's mechanism of action, emerging evidence implicates it in the adaptive transcriptional response to nutritional stress downstream of TORC1. The objectives of this study were to generate a murine model with knockin of the orthologous variant in mice (CREBRFR458Q) and to test the hypothesis that this CREBRF variant promotes obesity and protects against diabetes by regulating energy and glucose homeostasis downstream of TORC1. To test this hypothesis, we performed extensive phenotypic analysis of CREBRFR458Q knockin mice at baseline and in response to acute (fasting/refeeding), chronic (low- and high-fat diet feeding), and extreme (prolonged fasting) nutritional stress as well as with pharmacological TORC1 inhibition, and aging to 52 weeks. The results demonstrate that the murine CREBRFR458Q model of the human CREBRFR457Q variant does not influence energy/glucose homeostasis in response to these interventions, with the exception of possible greater loss of fat relative to lean mass with age. Alternative preclinical models and/or studies in humans will be required to decipher the mechanisms linking this variant to human health and disease.

Pleiotropic effects between cardiovascular disease risk factors and measures of cognitive and physical function in long-lived adults

Cardiovacular disease (CVD) is the leading cause of death among older adults and is often accompanied by functional decline. It is unclear what is driving this co-occurrence, but it may be behavioral, environmental and/or genetic. We used a family-based study to estimate the phenotypic and shared genetic correlation between CVD risk factors and physical and cognitive functional measures. Participants (n = 1,881) were from the Long Life Family Study, which enrolled families based on their exceptional longevity (sample mean age = 69.4 years, 44% female). Cardiovascular disease risk factors included carotid vessel measures [intima-media thickness and inter-adventitial diameter], obesity [body mass index (BMI) and waist circumference], and hypertension [systolic and diastolic blood pressures]. Function was measured in the physical [gait speed, grip strength, chair stand] and cognitive [digital symbol substitution test, retained and working memory, semantic fluency, and trail making tests] domains. We used SOLAR to estimate the genetic, environmental, and phenotypic correlation between each pair adjusting for age, age², sex, field center, smoking, height, and weight. There were significant phenotypic correlations (range |0.05–0.22|) between CVD risk factors and physical and cognitive function (all P < 0.05). Most significant genetic correlations (range |0.21–0.62|) were between CVD risk factorsand cognitive function, although BMI and waist circumference had significant genetic correlation with gait speed and chair stand time (range |0.29–0.53|; all P < 0.05). These results suggest that CVD risk factors may share a common genetic-and thus, biologic-basis with both cognitive and physical function. This is particularly informative for research into the genetic determinants of chronic disease.

The impact of global and local Polynesian genetic ancestry on complex traits in Native Hawaiians

Epidemiological studies of obesity, Type-2 diabetes (T2D), cardiovascular diseases and several common cancers have revealed an increased risk in Native Hawaiians compared to European- or Asian-Americans living in the Hawaiian islands. However, there remains a gap in our understanding of the genetic factors that affect the health of Native Hawaiians. To fill this gap, we studied the genetic risk factors at both the chromosomal and sub-chromosomal scales using genome-wide SNP array data on ~4,000 Native Hawaiians from the Multiethnic Cohort. We estimated the genomic proportion of Native Hawaiian ancestry ("global ancestry," which we presumed to be Polynesian in origin), as well as this ancestral component along each chromosome ("local ancestry") and tested their respective association with binary and quantitative cardiometabolic traits. After attempting to adjust for non-genetic covariates evaluated through questionnaires, we found that per 10% increase in global Polynesian genetic ancestry, there is a respective 8.6%, and 11.0% increase in the odds of being diabetic (P = 1.65×10-4) and having heart failure (P = 2.18×10-4), as well as a 0.059 s.d. increase in BMI (P = 1.04×10-10). When testing the association of local Polynesian ancestry with risk of disease or biomarkers, we identified a chr6 region associated with T2D. This association was driven by an uniquely prevalent variant in Polynesian ancestry individuals. However, we could not replicate this finding in an independent Polynesian cohort from Samoa due to the small sample size of the replication cohort. In conclusion, we showed that Polynesian ancestry, which likely capture both genetic and lifestyle risk factors, is associated with an increased risk of obesity, Type-2 diabetes, and heart failure, and that larger cohorts of Polynesian ancestry individuals will be needed to replicate the putative association on chr6 with T2D.

ECHS1 disease in two unrelated families of Samoan descent: Common variant - rare disorder

Mutations in the short-chain enoyl-CoA hydratase (SCEH) gene, ECHS1, cause a rare autosomal recessive disorder of valine catabolism. Patients usually present with developmental delay, regression, dystonia, feeding difficulties, and abnormal MRI with bilateral basal ganglia involvement. We present clinical, biochemical, molecular, and functional data for four affected patients from two unrelated families of Samoan descent with identical novel compound heterozygous mutations. Family 1 has three affected boys while Family 2 has an affected daughter, all with clinical and MRI findings of Leigh syndrome and intermittent episodes of acidosis and ketosis. WES identified a single heterozygous variant in ECHS1 at position c.832G > A (p.Ala278Thr). However, western blot revealed significantly reduced ECHS1 protein for all affected family members. Decreased SCEH activity in fibroblasts and a mild increase in marker metabolites in urine further supported ECHS1 as the underlying gene defect. Additional investigations at the DNA (aCGH, WGS) and RNA (qPCR, RT-PCR, RNA-Seq, RNA-Array) level identified a silent, common variant at position c.489G > A (p.Pro163=) as the second mutation. This substitution, present at high frequency in the Samoan population, is associated with decreased levels of normally spliced mRNA. To our understanding, this is the first report of a novel, hypomorphic allele c.489G > A (p.Pro163=), associated with SCEH deficiency.

Genetic Variability in the Iron Homeostasis Pathway and Patient Outcomes After Aneurysmal Subarachnoid Hemorrhage

BACKGROUND/OBJECTIVE

Iron can be detrimental to most tissues both in excess and in deficiency. The brain in particular is highly susceptible to the consequences of excessive iron, especially during blood brain barrier disruption after injury. Preliminary evidence suggests that iron homeostasis is important during recovery after neurologic injury; therefore, the exploration of genetic variability in genes involved in iron homeostasis is an important area of patient outcomes research. The purpose of this study was to examine the relationship between tagging single nucleotide polymorphisms (SNPs) in candidate genes related to iron homeostasis and acute and long-term patient outcomes after aneurysmal subarachnoid hemorrhage (aSAH).

METHODS

This study was a longitudinal, observational, candidate gene association study of participants with aSAH that used a two-tier design including tier 1 (discovery, n = 197) and tier 2 (replication, n = 277). Participants were followed during the acute outcome phase for development of cerebral vasospasm and delayed cerebral ischemia (DCI) and during the long-term outcome phase for death and gross functional outcome using the Glasgow Outcome Scale (GOS; poor = 1-3). Genetic association analyses were performed using a logistic regression model adjusted for age, sex, and Fisher grade. Approximate Bayes factors (ABF) and Bayesian false discovery probabilities (BFDP) were used to prioritize and interpret results.

RESULTS

In tier 1, 235 tagging SNPs in 28 candidate genes were available for analysis and 26 associations (20 unique SNPs in 12 genes) were nominated for replication in tier 2. In tier 2, we observed an increase in evidence of association for three associations in the ceruloplasmin (CP) and cubilin (CUBN) genes. We observed an association of rs17838831 (CP) with GOS at 3 months (tier 2 results, odds ratio [OR] = 2.10, 95% confidence interval [CI] = 1.14-3.86, p = 0.018, ABF = 0.52, and BFDP = 70.8%) and GOS at 12 months (tier 2 results, OR = 1.86, 95% CI 0.98-3.52, p = 0.058, ABF = 0.72, and BFDP = 77.3%) as well as rs10904850 (CUBN) with DCI (tier 2 results, OR = 0.70, 95% CI 0.48-1.02, p = 0.064, ABF = 0.59, and BFDP = 71.8%).

CONCLUSIONS

Among the genes examined, our findings support a role for CP and CUBN in patient outcomes after aSAH. In an effort to translate these findings into clinical utility and improve outcomes after aSAH, additional research is needed to examine the functional roles of these genes after aSAH.

Relationship Between Serum IGF-1 and BMI Differs by Age

BACKGROUND

Serum levels of insulin-like growth factor 1 (IGF-1) and body mass index (BMI) are both associated with susceptibility to age-related diseases. Reports on the correlation between them have been conflicting, with both positive to negative correlations reported. However, the age ranges of the participants varied widely among these studies.

METHODS

Using data on 4241 participants (aged 24-110) from the Long Life Family Study, we investigated the relationship between IGF-1 and BMI by age groups using regression analysis.

RESULTS

When stratified by age quartile, the relationship between IGF-1 and BMI varied: in the first quartile (Q1, 20-58 years) the relationship was negative (β = -0.2, p = .002); in Q2 (58-66 years) and Q3 (67-86 years) the relationship was negative (β = -0.07, β = -0.01, respectively) but nonsignificant; and in Q4 (87-110 years) the relationship was positive (β = 0.31, p = .0002). This pattern did not differ by sex. We observed a similar age-related pattern between IGF-1 and BMI among participants in the third National Health and Nutritional Examination Survey.

CONCLUSIONS

Our results that the relationship between IGF-1 and BMI differs by age may explain some of the inconsistency in reports about their relationship and encourage additional studies to understand the mechanisms underlying it.

Genome-wide association studies in Samoans give insight into the genetic architecture of fasting serum lipid levels

The current understanding of the genetic architecture of lipids has largely come from genome-wide association studies (GWAS). To date, few GWAS have examined the genetic architecture of lipids in Polynesians, and none have in Samoans, whose unique population history, including many population bottlenecks, may provide insight into the biological foundations of variation in lipid levels. Here we performed a GWAS of four fasting serum lipid levels: total cholesterol (TC), high-density lipoprotein (HDL), low-density lipoprotein (LDL), and triglycerides (TG) in a sample of 2849 Samoans, with validation genotyping for associations in a replication cohort comprising 1798 Samoans and American Samoans. We identified multiple genome-wide significant associations (P < 5 × 10^-8) previously seen in other populations-APOA1 with TG, CETP with HDL, and APOE with TC and LDL-and several suggestive associations (P < 1 × 10^-5), including an association of variants downstream of MGAT1 and RAB21 with HDL. However, we observed different association signals for variants near APOE than what has been previously reported in non-Polynesian populations. The association with several known lipid loci combined with the newly identified associations with variants near MGAT1 and RAB21 suggest that while some of the genetic architecture of lipids is shared between Samoans and other populations, part of the genetic architecture may be Polynesian-specific.

Exploring the Paradoxical Relationship of a Creb 3 Regulatory Factor Missense Variant With Body Mass Index and Diabetes Among Samoans: Protocol for the Soifua Manuia (Good Health) Observational Cohort Study

BACKGROUND

The prevalence of obesity and diabetes in Samoa, like many other Pacific Island nations, has reached epidemic proportions. Although the etiology of these conditions can be largely attributed to the rapidly changing economic and nutritional environment, a recently identified genetic variant, rs373863828 (CREB 3 regulatory factor, CREBRF: c.1370G>A p.[R457Q]) is associated with increased odds of obesity, but paradoxically, decreased odds of diabetes.

OBJECTIVE

The overarching goal of the Soifua Manuia (Good Health) study was to precisely characterize the association of the CREBRF variant with metabolic (body composition and glucose homeostasis) and behavioral traits (dietary intake, physical activity, sleep, and weight control behaviors) that influence energy homeostasis in 500 adults.

METHODS

A cohort of adult Samoans who participated in a genome-wide association study of adiposity in Samoa in 2010 was followed up, based on the presence or absence of the CREBRF variant, between August 2017 and March 2019. Over a period of 7-10 days, each participant completed the main study protocol, which consisted of anthropometric measurements (weight, height, circumferences, and skinfolds), body composition assessment (bioelectrical impedance and dual-energy x-ray absorptiometry), point-of-care glycated hemoglobin measurement, a fasting blood draw and oral glucose tolerance test, urine collection, blood pressure measurement, hand grip strength measurement, objective physical activity and sleep apnea monitoring, and questionnaire measures (eg, health interview, cigarette and alcohol use, food frequency questionnaire, socioeconomic position, stress, social support, food and water insecurity, sleep, body image, and dietary preferences). In January 2019, a subsample of the study participants (n=118) completed a buttock fat biopsy procedure to collect subcutaneous adipose tissue samples.

RESULTS

Enrollment of 519 participants was completed in March 2019. Data analyses are ongoing, with results expected in 2020 and 2021.

CONCLUSIONS

While the genetic variant rs373863828, in CREBRF, has the largest known effect size of any identified common obesity gene, very little is currently understood about the mechanisms by which it confers increased odds of obesity but paradoxically lowered odds of type 2 diabetes. The results of this study will provide insights into how the gene functions on a whole-body level, which could provide novel targets to prevent or treat obesity, diabetes, and associated metabolic disorders. This study represents the human arm of a comprehensive and integrated approach involving humans as well as preclinical models that will provide novel insights into metabolic disease.

INTERNATIONAL REGISTERED REPORT IDENTIFIER (IRRID)

RR1-10.2196/17329.

Evolutionary history of modern Samoans

Archaeological studies estimate the initial settlement of Samoa at 2,750 to 2,880 y ago and identify only limited settlement and human modification to the landscape until about 1,000 to 1,500 y ago. At this point, a complex history of migration is thought to have begun with the arrival of people sharing ancestry with Near Oceanic groups (i.e., Austronesian-speaking and Papuan-speaking groups), and was then followed by the arrival of non-Oceanic groups during European colonialism. However, the specifics of this peopling are not entirely clear from the archaeological and anthropological records, and is therefore a focus of continued debate. To shed additional light on the Samoan population history that this peopling reflects, we employ a population genetic approach to analyze 1,197 Samoan high-coverage whole genomes. We identify population splits between the major Samoan islands and detect asymmetrical gene flow to the capital city. We also find an extreme bottleneck until about 1,000 y ago, which is followed by distinct expansions across the islands and subsequent bottlenecks consistent with European colonization. These results provide for an increased understanding of Samoan population history and the dynamics that inform it, and also demonstrate how rapid demographic processes can shape modern genomes.

A missense variant in CREBRF is associated with taller stature in Samoans

OBJECTIVES

Studies have demonstrated that rs373863828, a missense variant in CREBRF, is associated with a number of anthropometric traits including body mass index (BMI), obesity, percent body fat, hip circumference, and abdominal circumference. Given the biological relationship between height and adiposity, we hypothesized that the effect of this variant on BMI might be due in part to an association of this variant with height.

METHODS

We tested the hypothesis that minor allele of rs373863828 is associated with height in a Samoan population in two adult cohorts and in a separate cohort of children (age 5-18 years old) using linear mixed modeling.

RESULTS

We found evidence of a strong relationship between rs373863828 and greater mean height in Samoan adults (0.77 cm greater average height for each copy of the minor allele) with the same direction of effect in Samoan children.

CONCLUSIONS

These results suggest that the missense variant rs373863828 in CREBRF, first identified through an association with larger BMI, may be related to an underlying biological mechanism affecting overall body size including stature.

Genome-wide linkage analysis of carotid artery traits in exceptionally long-lived families

BACKGROUND AND AIMS

Atherosclerosis develops with age and is partially controlled by genetics. Research to date has identified common variants with small effects on atherosclerosis related traits. We aimed to use family-based genome-wide linkage analysis to identify chromosomal regions potentially harboring rare variants with larger effects for atherosclerosis related traits.

METHODS

Participants included 2205 individuals from the Long Life Family Study (LLFS), which recruited families with exceptional longevity from Boston, New York, Pittsburgh, and Denmark. Participants underwent B-mode ultrasonography of the carotid arteries to measure intima-media thickness (IMT), inter-adventitial diameter (IAD), and plaque presence and severity. We conducted residual heritability and genome-wide linkage analyses adjusted for age, age², sex, and field center using pedigree-based maximum-likelihood methods in SOLAR.

RESULTS

All carotid traits were significantly heritable with a range of 0.68 for IAD to 0.38 for IMT. We identified three chromosomal regions with linkage to IAD (3q13; max LOD 5.3), plaque severity (17q22-q23, max LOD 3.2), and plaque presence (17q24, max LOD 3.1). No common allelic variants within these linkage peaks were associated with the carotid artery traits.

CONCLUSIONS

We identified three chromosomal regions with evidence of linkage to carotid artery diameter and atherosclerotic plaque in exceptionally long-lived families. Since common allelic variants within our linkage peaks did not account for our findings, future follow-up resequencing of these regions in LLFS families should help advance our understanding of atherosclerosis, CVD, and healthy vascular aging.

Discordant association of the CREBRF rs373863828 A allele with increased BMI and protection from type 2 diabetes in Māori and Pacific (Polynesian) people living in Aotearoa/New Zealand

AIMS/HYPOTHESIS

The A (minor) allele of CREBRF rs373863828 has been associated with increased BMI and reduced risk of type 2 diabetes in the Samoan populations of Samoa and American Samoa. Our aim was to test rs373863828 for associations with BMI and the odds of type 2 diabetes, gout and chronic kidney disease (CKD) in Māori and Pacific (Polynesian) people living in Aotearoa/New Zealand.

METHODS

Linear and logistic regression models were used to analyse the association of the A allele of CREBRF rs373863828 with BMI, log-transformed BMI, waist circumference, type 2 diabetes, gout and CKD in 2286 adults. The primary analyses were adjusted for age, sex, the first four genome-wide principal components and (where appropriate) BMI, waist circumference and type 2 diabetes. The primary analysis was conducted in ancestrally defined groups and association effects were combined using meta-analysis.

RESULTS

For the A allele of rs373863828, the effect size was 0.038 (95% CI 0.022, 0.055, p = 4.8 × 10-6) for log-transformed BMI, with OR 0.59 (95% CI 0.47, 0.73, p = 1.9 × 10-6) for type 2 diabetes. There was no evidence for an association of genotype with variance in BMI (p = 0.13), and nor was there evidence for associations with serum urate (β = 0.012 mmol/l, pcorrected = 0.10), gout (OR 1.00, p = 0.98) or CKD (OR 0.91, p = 0.59).

CONCLUSIONS/INTERPRETATION

Our results in New Zealand Polynesian adults replicate, with very similar effect sizes, the association of the A allele of rs373863828 with higher BMI but lower odds of type 2 diabetes among Samoan adults living in Samoa and American Samoa.

A novel healthy blood pressure phenotype in the Long Life Family Study

BACKGROUND

Hypertension tends to run in families and has both genetic and environmental determinants. We assessed the hypothesis that a novel healthy blood pressure (BP) phenotype is also familial and sought to identify its associated factors.

METHODS

We developed a healthy BP phenotype in the Long Life Family Study, a cohort of two-generation families selected for longevity. Participants from the offspring generation (n = 2211, ages 32-88) were classified as having healthy BP if their age-adjusted and sex-adjusted SBP z-score was between -1.5 and -0.5. Offspring on antihypertensive medications were classified as not having healthy BP. Families with at least two offspring (n = 419 families) were defined as meeting the healthy BP phenotype if at least two and at least 50% of their offspring had healthy BP.

RESULTS

Among 2211 offspring, 476 (21.5%) met the healthy BP phenotype. When examining the 419 families, only 44 (10.5%) families met the criteria for the healthy BP phenotype. Both offspring and probands from families with healthy BP performed better on neuropsychological tests that place demands on complex attention and executive function when compared with offspring and probands from remaining families. Among families with the healthy BP phenotype compared with families without, a higher proportion of offspring met the American Heart Association definition of ideal cardiovascular health (10.8 versus 3.8%, respectively; driven by BP, smoking status, and BMI components).

CONCLUSION

In this cohort of familial longevity, few families had a novel healthy BP phenotype in multiple members. Families with this healthy BP phenotype may represent a specific pathway to familial longevity.

Genomewide Association Scan of a Mortality Associated Endophenotype for a Long and Healthy Life in the Long Life Family Study

Background

Identification of genes or fundamental biological pathways that regulate aging phenotypes and longevity could lead to possible interventions to increase healthy longevity.

Methods

Using data from the Long Life Family Study, we performed genomewide association analyses on an endophenotype construct, LF1, comprising a linear combination of traits across health domains. LF1 primarily reflected traits from the pulmonary and physical activity domains.

Results

We detected a significant association between LF1 and a locus on chromosome 10p15 (p-value = 4.65 × 10-8) and suggestive evidence (p-value < 5 × 10-6) for association on chromosomes 1, 2, 8, 12, 15, 18, and 22. Using data from the Health, Aging and Body Composition Study, we subsequently replicated the association for the 1p13 region near the NBPF6 locus (p-value = 3.65 × 10-4).

Conclusions

Our analyses indicate that loci influencing a healthy aging endophenotype construct predominantly comprised of pulmonary and physical function domains may be located on chromosome 1p13 near the NBPF6 locus. Further investigation of this possible locus and other suggestive loci may reveal novel biological pathways that influence healthy aging.

Association Between Mortality and Heritability of the Scale of Aging Vigor in Epidemiology

OBJECTIVES

To investigate the association between mortality and heritability of a rescaled Fried frailty index, the Scale of Aging Vigor in Epidemiology (SAVE), to determine its value for genetic analyses.

DESIGN

Longitudinal, community-based cohort study.

SETTING

The Long Life Family Study (LLFS) in the United States and Denmark.

PARTICIPANTS

Long-lived individuals (N = 4,875, including 4,075 genetically related individuals) and their families (N = 551).

MEASUREMENTS

The SAVE was administered to 3,599 participants and included weight change, weakness (grip strength), fatigue (questionnaire), physical activity (days walked in prior 2 weeks), and slowness (gait speed); each component was scored 0, 1, or 2 using approximate tertiles, and summed (range 0 (vigorous) to 10 (frail)). Heritability was determined using a variance component-based family analysis using a polygenic model. Association with mortality in the proband generation (N = 1,421) was calculated using Cox proportional hazards mixed-effect models.

RESULTS

Heritability of the SAVE was 0.23 (P < .001) overall (n = 3,599), 0.31 (P < .001) in probands (n = 1,479), and 0.26 (P < .001) in offspring (n = 2,120). In adjusted models, higher SAVE scores were associated with higher mortality (score 5-6: hazard ratio (HR) = 2.83, 95% confidence interval (CI) = 1.46-5.51; score 7-10: HR = 3.40, 95% CI = 1.72-6.71) than lower scores (0-2).

CONCLUSION

The SAVE was associated with mortality and was moderately heritable in the LLFS, suggesting a genetic component to age-related vigor and frailty and supporting its use for further genetic analyses.

Genome-Wide Association Study and Linkage Analysis of the Healthy Aging Index

BACKGROUND

The Healthy Aging Index (HAI) is a tool for measuring the extent of health and disease across multiple systems.

METHODS

We conducted a genome-wide association study and a genome-wide linkage analysis to map quantitative trait loci associated with the HAI and a modified HAI weighted for mortality risk in 3,140 individuals selected for familial longevity from the Long Life Family Study. The genome-wide association study used the Long Life Family Study as the discovery cohort and individuals from the Cardiovascular Health Study and the Framingham Heart Study as replication cohorts.

RESULTS

There were no genome-wide significant findings from the genome-wide association study; however, several single-nucleotide polymorphisms near ZNF704 on chromosome 8q21.13 were suggestively associated with the HAI in the Long Life Family Study (p < 10(-) (6)) and nominally replicated in the Cardiovascular Health Study and Framingham Heart Study. Linkage results revealed significant evidence (log-odds score = 3.36) for a quantitative trait locus for mortality-optimized HAI in women on chromosome 9p24-p23. However, results of fine-mapping studies did not implicate any specific candidate genes within this region of interest.

CONCLUSIONS

ZNF704 may be a potential candidate gene for studies of the genetic underpinnings of longevity.

Genetic variants associated with lung function: the long life family study

Background

Reduced forced expiratory volume in 1 second (FEV₁) and the ratio of FEV₁ to forced vital capacity (FVC) are strong predictors of mortality and lung function is higher among individuals with exceptional longevity. However, genetic factors associated with lung function in individuals with exceptional longevity have not been identified.

Method

We conducted a genome wide association study (GWAS) to identify novel genetic variants associated with lung function in the Long Life Family Study (LLFS) (n = 3,899). Replication was performed using data from the CHARGE/SpiroMeta consortia. The association between SNPs and FEV₁ and FEV₁/FVC was analyzed using a linear mixed effects model adjusted for age, age², sex, height, field center, ancestry principal components and kinship structure to adjust for family relationships separately for ever smokers and never smokers. In the linkage analysis, we used the residuals of the FEV₁ and FEV₁/FVC, adjusted for age, sex, height, ancestry principal components (PCs), smoking status, pack-years, and field center.

Results

We identified nine SNPs in strong linkage disequilibrium in the CYP2U1 gene to be associated with FEV₁ and a novel SNP (rs889574) associated with FEV₁/FVC, none of which were replicated in the CHARGE/SpiroMeta consortia. Using linkage analysis, we identified a novel linkage peak in chromosome 2 at 219 cM for FEV₁/FVC (LOD: 3.29) and confirmed a previously reported linkage peak in chromosome 6 at 28 cM (LOD: 3.33) for FEV₁.

Conclusion

Future studies need to identify the rare genetic variants underlying the linkage peak in chromosome 6 for FEV₁.

Electronic supplementary material

The online version of this article (doi:10.1186/s12931-014-0134-x) contains supplementary material, which is available to authorized users.

Stratified randomization controls better for batch effects in 450K methylation analysis: a cautionary tale

BACKGROUND

Batch effects in DNA methylation microarray experiments can lead to spurious results if not properly handled during the plating of samples.

METHODS

Two pilot studies examining the association of DNA methylation patterns across the genome with obesity in Samoan men were investigated for chip- and row-specific batch effects. For each study, the DNA of 46 obese men and 46 lean men were assayed using Illumina's Infinium HumanMethylation450 BeadChip. In the first study (Sample One), samples from obese and lean subjects were examined on separate chips. In the second study (Sample Two), the samples were balanced on the chips by lean/obese status, age group, and census region. We used methylumi, watermelon, and limma R packages, as well as ComBat, to analyze the data. Principal component analysis and linear regression were, respectively, employed to identify the top principal components and to test for their association with the batches and lean/obese status. To identify differentially methylated positions (DMPs) between obese and lean males at each locus, we used a moderated t-test.

RESULTS

Chip effects were effectively removed from Sample Two but not Sample One. In addition, dramatic differences were observed between the two sets of DMP results. After "removing" batch effects with ComBat, Sample One had 94,191 probes differentially methylated at a q-value threshold of 0.05 while Sample Two had zero differentially methylated probes. The disparate results from Sample One and Sample Two likely arise due to the confounding of lean/obese status with chip and row batch effects.

CONCLUSION

Even the best possible statistical adjustments for batch effects may not completely remove them. Proper study design is vital for guarding against spurious findings due to such effects.

Prevalence of adiposity and associated cardiometabolic risk factors in the Samoan genome-wide association study

OBJECTIVE

To describe the prevalence of obesity-related noncommunicable diseases (NCDs) and associated risk factors in a sample of Samoan adults studied in 2010 as part of a genome-wide assocation study (GWAS) for obesity related traits.

METHODS

Anthropometric and biochemical data collected from n = 3,475 participants (n = 1,437 male; n = 2,038 female) aged 24.5 to <65 years were used to describe the prevalence of obesity, diabetes, hypertension, and dyslipidemia within the study sample. One way analysis of variance, χ(2) tests, and binary logistic regression were used to identify differences in disease and risk factor prevalence by 10-year age group, gender, or by census region of residence.

RESULTS

Obesity was highly prevalent among the study sample; 64.6% of females and 41.2% of males were obese according to Polynesian cutoffs (BMI ≥ 32 kg/m(2) ). Females were less likely than males to have hypertension (31.7% vs. 36.7%) but equally likely to have diabetes (17.8% vs. 16.4%). With the exception of obesity and low HDL-cholesterol in females only, there were significant differences in the prevalence of all NCDs and associated risk factors by age group, with the oldest age group (55 to <65 years) most affected. In both sexes, residents of the Apia Urban Area were at significantly greater risk of obesity, diabetes, low HDL-cholesterol, and high triglycerides than residents of the more rural Savaii region.

CONCLUSIONS

The phenotypic characteristics of this sample provide evidence of a continuation of previously reported temporal trends toward obesity and its associated disorders. Attention must be paid to the critical NCD situation in Samoa.

Genetic determinants of heel bone properties: genome-wide association meta-analysis and replication in the GEFOS/GENOMOS consortium

Quantitative ultrasound of the heel captures heel bone properties that independently predict fracture risk and, with bone mineral density (BMD) assessed by X-ray (DXA), may be convenient alternatives for evaluating osteoporosis and fracture risk. We performed a meta-analysis of genome-wide association (GWA) studies to assess the genetic determinants of heel broadband ultrasound attenuation (BUA; n = 14 260), velocity of sound (VOS; n = 15 514) and BMD (n = 4566) in 13 discovery cohorts. Independent replication involved seven cohorts with GWA data (in silico n = 11 452) and new genotyping in 15 cohorts (de novo n = 24 902). In combined random effects, meta-analysis of the discovery and replication cohorts, nine single nucleotide polymorphisms (SNPs) had genome-wide significant (P < 5 × 10(-8)) associations with heel bone properties. Alongside SNPs within or near previously identified osteoporosis susceptibility genes including ESR1 (6q25.1: rs4869739, rs3020331, rs2982552), SPTBN1 (2p16.2: rs11898505), RSPO3 (6q22.33: rs7741021), WNT16 (7q31.31: rs2908007), DKK1 (10q21.1: rs7902708) and GPATCH1 (19q13.11: rs10416265), we identified a new locus on chromosome 11q14.2 (rs597319 close to TMEM135, a gene recently linked to osteoblastogenesis and longevity) significantly associated with both BUA and VOS (P < 8.23 × 10(-14)). In meta-analyses involving 25 cohorts with up to 14 985 fracture cases, six of 10 SNPs associated with heel bone properties at P < 5 × 10(-6) also had the expected direction of association with any fracture (P < 0.05), including three SNPs with P < 0.005: 6q22.33 (rs7741021), 7q31.31 (rs2908007) and 10q21.1 (rs7902708). In conclusion, this GWA study reveals the effect of several genes common to central DXA-derived BMD and heel ultrasound/DXA measures and points to a new genetic locus with potential implications for better understanding of osteoporosis pathophysiology.

Heritability of and mortality prediction with a longevity phenotype: the healthy aging index

BACKGROUND

Longevity-associated genes may modulate risk for age-related diseases and survival. The Healthy Aging Index (HAI) may be a subphenotype of longevity, which can be constructed in many studies for genetic analysis. We investigated the HAI's association with survival in the Cardiovascular Health Study and heritability in the Long Life Family Study.

METHODS

The HAI includes systolic blood pressure, pulmonary vital capacity, creatinine, fasting glucose, and Modified Mini-Mental Status Examination score, each scored 0, 1, or 2 using approximate tertiles and summed from 0 (healthy) to 10 (unhealthy). In Cardiovascular Health Study, the association with mortality and accuracy predicting death were determined with Cox proportional hazards analysis and c-statistics, respectively. In Long Life Family Study, heritability was determined with a variance component-based family analysis using a polygenic model.

RESULTS

Cardiovascular Health Study participants with unhealthier index scores (7-10) had 2.62-fold (95% confidence interval: 2.22, 3.10) greater mortality than participants with healthier scores (0-2). The HAI alone predicted death moderately well (c-statistic = 0.643, 95% confidence interval: 0.626, 0.661, p < .0001) and slightly worse than age alone (c-statistic = 0.700, 95% confidence interval: 0.684, 0.717, p < .0001; p < .0001 for comparison of c-statistics). Prediction increased significantly with adjustment for demographics, health behaviors, and clinical comorbidities (c-statistic = 0.780, 95% confidence interval: 0.765, 0.794, p < .0001). In Long Life Family Study, the heritability of the HAI was 0.295 (p < .0001) overall, 0.387 (p < .0001) in probands, and 0.238 (p = .0004) in offspring.

CONCLUSION

The HAI should be investigated further as a candidate phenotype for uncovering longevity-associated genes in humans.

Genome-wide association analysis of age-at-onset in Alzheimer's disease

The risk of Alzheimer's disease (AD) is strongly determined by genetic factors and recent genome-wide association studies (GWAS) have identified several genes for the disease risk. In addition to the disease risk, age-at-onset (AAO) of AD has also strong genetic component with an estimated heritability of 42%. Identification of AAO genes may help to understand the biological mechanisms that regulate the onset of the disease. Here we report the first GWAS focused on identifying genes for the AAO of AD. We performed a genome-wide meta-analysis on three samples comprising a total of 2222 AD cases. A total of ~2.5 million directly genotyped or imputed single-nucleotide polymorphisms (SNPs) were analyzed in relation to AAO of AD. As expected, the most significant associations were observed in the apolipoprotein E (APOE) region on chromosome 19 where several SNPs surpassed the conservative genome-wide significant threshold (P<5E-08). The most significant SNP outside the APOE region was located in the DCHS2 gene on chromosome 4q31.3 (rs1466662; P=4.95E-07). There were 19 additional significant SNPs in this region at P<1E-04 and the DCHS2 gene is expressed in the cerebral cortex and thus is a potential candidate for affecting AAO in AD. These findings need to be confirmed in additional well-powered samples.

Assessment of gene-by-sex interaction effect on bone mineral density

Sexual dimorphism in various bone phenotypes, including bone mineral density (BMD), is widely observed; however, the extent to which genes explain these sex differences is unclear. To identify variants with different effects by sex, we examined gene-by-sex autosomal interactions genome-wide, and performed expression quantitative trait loci (eQTL) analysis and bioinformatics network analysis. We conducted an autosomal genome-wide meta-analysis of gene-by-sex interaction on lumbar spine (LS) and femoral neck (FN) BMD in 25,353 individuals from 8 cohorts. In a second stage, we followed up the 12 top single-nucleotide polymorphisms (SNPs; p < 1 × 10(-5) ) in an additional set of 24,763 individuals. Gene-by-sex interaction and sex-specific effects were examined in these 12 SNPs. We detected one novel genome-wide significant interaction associated with LS-BMD at the Chr3p26.1-p25.1 locus, near the GRM7 gene (male effect = 0.02 and p = 3.0 × 10(-5) ; female effect = -0.007 and p = 3.3 × 10(-2) ), and 11 suggestive loci associated with either FN- or LS-BMD in discovery cohorts. However, there was no evidence for genome-wide significant (p < 5 × 10(-8) ) gene-by-sex interaction in the joint analysis of discovery and replication cohorts. Despite the large collaborative effort, no genome-wide significant evidence for gene-by-sex interaction was found to influence BMD variation in this screen of autosomal markers. If they exist, gene-by-sex interactions for BMD probably have weak effects, accounting for less than 0.08% of the variation in these traits per implicated SNP. © 2012 American Society for Bone and Mineral Research.

Association of CLU and PICALM variants with Alzheimer's disease

Two recent large genome-wide association studies have reported significant associations in the CLU (APOJ), CR1, and PICALM genes with the risk of Alzheimer's disease (AD). In order to replicate these findings, we examined 7 single nucleotide polymorphisms (SNPs) most significantly implicated by these studies in a large case-control sample comprising 2707 individuals. Principle components analysis revealed no population substructure in our sample. While no association was observed with CR1 SNPs (p = 0.30-0.457), a trend of association was seen with the PICALM (p = 0.071-0.086) and CLU (p = 0.148-0.258) SNPs. A meta-analysis of 3 studies revealed significant associations with all 3 genes. Our data from an independent and large case-control sample suggest that these gene regions should be followed up by comprehensive resequencing to find functional variants.

No association between CALHM1 variation and risk of Alzheimer disease

A polymorphism in the calcium homeostasis modulator 1 gene (CALHM1) has recently been associated with risk of late-onset Alzheimer disease. We examined this variant (rs2986017) in 945 Caucasian Americans with late-onset Alzheimer disease and 875 age-matched Caucasian American controls. No association with risk of late-onset Alzheimer disease (p=0.368 for genotypes; p=0.796 for alleles) was observed in our sample. However, a potential modest association of minor allele homozygosity (TT) with an earlier age-at-onset was seen (p=0.034).