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

Powerful, scalable and resource-efficient meta-analysis of rare variant associations in large whole genome sequencing studies.

Meta-analysis of whole genome sequencing/whole exome sequencing (WGS/WES) studies provides an attractive solution to the problem of collecting large sample sizes for discovering rare variants associated with complex phenotypes. Existing rare variant meta-analysis approaches are not scalable to biobank-scale WGS data. Here we present MetaSTAAR, a powerful and resource-efficient rare variant meta-analysis framework for large-scale WGS/WES studies. MetaSTAAR accounts for relatedness and population structure, can analyze both quantitative and dichotomous traits and boosts the power of rare variant tests by incorporating multiple variant functional annotations. Through meta-analysis of four lipid traits in 30,138 ancestrally diverse samples from 14 studies of the Trans Omics for Precision Medicine (TOPMed) Program, we show that MetaSTAAR performs rare variant meta-analysis at scale and produces results comparable to using pooled data. Additionally, we identified several conditionally significant rare variant associations with lipid traits. We further demonstrate that MetaSTAAR is scalable to biobank-scale cohorts through meta-analysis of TOPMed WGS data and UK Biobank WES data of ~200,000 samples.

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.

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.

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.

Rare coding variants in RCN3 are associated with blood pressure.

BACKGROUND: While large genome-wide association studies have identified nearly one thousand loci associated with variation in blood pressure, rare variant identification is still a challenge. In family-based cohorts, genome-wide linkage scans have been successful in identifying rare genetic variants for blood pressure. This study aims to identify low frequency and rare genetic variants within previously reported linkage regions on chromosomes 1 and 19 in African American families from the Trans-Omics for Precision Medicine (TOPMed) program. Genetic association analyses weighted by linkage evidence were completed with whole genome sequencing data within and across TOPMed ancestral groups consisting of 60,388 individuals of European, African, East Asian, Hispanic, and Samoan ancestries. RESULTS: Associations of low frequency and rare variants in RCN3 and multiple other genes were observed for blood pressure traits in TOPMed samples. The association of low frequency and rare coding variants in RCN3 was further replicated in UK Biobank samples (N = 403,522), and reached genome-wide significance for diastolic blood pressure (p = 2.01 × 10- 7). CONCLUSIONS: Low frequency and rare variants in RCN3 contributes blood pressure variation. This study demonstrates that focusing association analyses in linkage regions greatly reduces multiple-testing burden and improves power to identify novel rare variants associated with blood pressure traits.

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.

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.

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.

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.

Associations between socioeconomic resources and adiposity traits in adults: Evidence from Samoa.

In low- and middle-income countries, earlier in economic development, obesity tends to be more prevalent in high socioeconomic resource groups compared to low. Later in development, the distribution of obesity tends to show the opposite pattern, becoming more prevalent in those with low socioeconomic resources. This shift in obesity prevalence tends to occur between a gross national income per capita (GNI) of US$1,000 to $4,000 dollars. Whether a similar pattern occurs in Pacific Island countries has not been well documented. In Samoa, the GNI rose to US$3,200 dollars in 2010 at which time over 80% of adults were overweight or obese. We aimed to understand the association of socioeconomic resources, assessed by household assets, with adult body mass index (BMI) and abdominal circumference (AC) in Samoa. Data were from a genome-wide association study for obesity among 3,370 Samoans aged 24.5-<65 years in 2010. Household asset scores were calculated based on ownership of consumer durables, housing construction, and access to basic services. Sex-stratified multivariate linear regressions were used to assess adiposity trait differences by household asset ownership, after controlling for age, education, and household urbanicity. Higher asset ownership was associated with higher BMI and AC and the positive relationship remained robust after controlling for potential confounders. Despite significant economic growth preceding the year 2010 in Samoa, the obesity burden had not shifted to low socioeconomic groups in a similar way that has been observed in countries further along in the economic transition. The mechanism by which socioeconomic resources influence adiposity is complex and may be particularly complicated in Samoa by migrant remittances received both as cash and household assets. Social and physical environments may constrain the positive health behavior change necessary to reduce obesity even in the context of high socioeconomic position, a situation that requires further investigation.

Tobacco smoking patterns in Samoa in 2010: Implications for interventions.

INTRODUCTION: Tobacco use in Samoa has been described over time by age, sex and education, but little work exists on other sociodemographic factors associated with tobacco use. We describe current smoking and daily tobacco use in adults from Samoa, with a focus on sex and age stratified analyses of the influence of occupation, education, census region, household asset ownership and alcohol use in order to help develop potential targeted interventions. METHODS: In 2010, a nationwide survey of 3745 adults aged 25-65 years from 33 villages was completed in Samoa. Current smoking status, daily tobacco use, as well as current alcohol use, and a variety of sociodemographic factors were assessed by interview. Bivariate and multivariable models, and sex and age group stratified analyses, were performed to determine different patterns of correlates. RESULTS: More than half of all men (51.3%) and 21.8% of women were current tobacco smokers. Men and women smoked on average 10.9 and 8.7 cigarettes/day, respectively. Twenty per cent of men consumed ≥20 cigarettes/day. In men, being married, a subsistence-farmer/laborer, an alcohol user, and having low household assets, were independently associated with being a tobacco smoker (all p<0.01). Among women, not completing secondary education, being 25-34 years, residing in urban Apia, and being an alcohol user, were independently associated with being a tobacco smoker (all p<0.01). CONCLUSIONS: Tobacco use in Samoa remains high and correlates of smoking suggest that interventions for cessation need to be developed within the contexts of sex, age, education, and household socioeconomic status.

A thrifty variant in CREBRF strongly influences body mass index in Samoans.

Samoans are a unique founder population with a high prevalence of obesity, making them well suited for identifying new genetic contributors to obesity. We conducted a genome-wide association study (GWAS) in 3,072 Samoans, discovered a variant, rs12513649, strongly associated with body mass index (BMI) (P = 5.3 × 10(-14)), and replicated the association in 2,102 additional Samoans (P = 1.2 × 10(-9)). Targeted sequencing identified a strongly associated missense variant, rs373863828 (p.Arg457Gln), in CREBRF (meta P = 1.4 × 10(-20)). Although this variant is extremely rare in other populations, it is common in Samoans (frequency of 0.259), with an effect size much larger than that of any other known common BMI risk variant (1.36-1.45 kg/m(2) per copy of the risk-associated allele). In comparison to wild-type CREBRF, the Arg457Gln variant when overexpressed selectively decreased energy use and increased fat storage in an adipocyte cell model. These data, in combination with evidence of positive selection of the allele encoding p.Arg457Gln, support a 'thrifty' variant hypothesis as a factor in human obesity.