Associations between fasting glucose rate-of-change and the missense variant, rs373863828, in an adult Samoan cohort.

BACKGROUND: The A allele of rs373863828 in CREB3 regulatory factor is associated with high Body Mass Index, but lower odds of type 2 diabetes. These associations have been replicated elsewhere, but to date all studies have been cross-sectional. Our aims were (1) to describe the development of type 2 diabetes and change in fasting glucose between 2010 and 2018 among a longitudinal cohort of adult Samoans without type 2 diabetes or who were not using diabetes medications at baseline, and (2) to examine associations between fasting glucose rate-of-change (mmol/L per year) and the A allele of rs373863828. METHODS: We describe and test differences in fasting glucose, the development of type 2 diabetes, body mass index, age, smoking status, physical activity, urbanicity of residence, and household asset scores between 2010 and 2018 among a cohort of n = 401 adult Samoans, selected to have a ~2:2:1 ratio of GG:AG: AA rs373863828 genotypes. Multivariate linear regression was used to test whether fasting glucose rate-of-change was associated with rs373863828 genotype, and other baseline variables. RESULTS: By 2018, fasting glucose and BMI significantly increased among all genotype groups, and a substantial portion of the sample developed type 2 diabetes mellitus. The A allele was associated with a lower fasting glucose rate-of-change (ß = -0.05 mmol/L/year per allele, p = 0.058 among women; ß = -0.004 mmol/L/year per allele, p = 0.863 among men), after accounting for baseline variables. Mean fasting glucose and mean BMI increased over an eight-year period and a substantial number of individuals developed type 2 diabetes by 2018. However, fasting glucose rate-of-change, and type 2 diabetes development was lower among females with AG and AA genotypes. CONCLUSIONS: Further research is needed to understand the effect of the A allele on fasting glucose and type 2 diabetes development. Based on our observations that other risk factors increased over time, we advocate for the continued promotion for diabetes prevention and treatment programming, and the reduction of modifiable risk factors, in this setting.

Validation of human telomere length multi-ancestry meta-analysis association signals identifies POP5 and KBTBD6 as human telomere length regulation genes.

Genome-wide association studies (GWAS) have become well-powered to detect loci associated with telomere length. However, no prior work has validated genes nominated by GWAS to examine their role in telomere length regulation. We conducted a multi-ancestry meta-analysis of 211,369 individuals and identified five novel association signals. Enrichment analyses of chromatin state and cell-type heritability suggested that blood/immune cells are the most relevant cell type to examine telomere length association signals. We validated specific GWAS associations by overexpressing KBTBD6 or POP5 and demonstrated that both lengthened telomeres. CRISPR/Cas9 deletion of the predicted causal regions in K562 blood cells reduced expression of these genes, demonstrating that these loci are related to transcriptional regulation of KBTBD6 and POP5. Our results demonstrate the utility of telomere length GWAS in the identification of telomere length regulation mechanisms and validate KBTBD6 and POP5 as genes affecting telomere length regulation.

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.

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.

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 Maori 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 (Maori 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.

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.

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.

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.

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.

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.

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.

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/m2 in males (n = 220, p = 1.3 ×10-8) and from 35.9 to 37.8 kg/m2 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/m2/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/m2/year more than of those with the GG genotype.

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.

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, age2, 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.

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.

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.