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.

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

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.

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.

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.

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.

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.

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.

Discordant association of the CREBRF rs373863828 A allele with increased BMI and protection from type 2 diabetes in Maori 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 Maori 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.

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.

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

BACKGROUND: Reduced forced expiratory volume in 1 second (FEV1) and the ratio of FEV1 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 FEV1 and FEV1/FVC was analyzed using a linear mixed effects model adjusted for age, age2, 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 FEV1 and FEV1/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 FEV1 and a novel SNP (rs889574) associated with FEV1/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 FEV1/FVC (LOD: 3.29) and confirmed a previously reported linkage peak in chromosome 6 at 28 cM (LOD: 3.33) for FEV1. CONCLUSION: Future studies need to identify the rare genetic variants underlying the linkage peak in chromosome 6 for FEV1.

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.

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.

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.

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.