Search

The epithelial NaWe 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.We found that variants inOur results suggest that variants in extrarenal ENaCs, in addition to ENaCs expressed in kidneys, influence blood pressure and kidney function.

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

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

Background: The availability of whole-genome sequencing data in large studies has enabled the assessment of coding and noncoding variants across the allele frequency spectrum for their associations with blood pressure. Methods: We conducted a multiancestry whole-genome sequencing analysis of blood pressure among 51 456 Trans-Omics for Precision Medicine and Centers for Common Disease Genomics program participants (stage-1). Stage-2 analyses leveraged array data from UK Biobank (N=383 145), Million Veteran Program (N=318 891), and Reasons for Geographic and Racial Differences in Stroke (N=10 643) participants, along with whole-exome sequencing data from UK Biobank (N=199 631) participants. Results: Two blood pressure signals achieved genome-wide significance in meta-analyses of stage-1 and stage-2 single variant findings ( P -8 ). Among them, a rare intergenic variant at novel locus, LOC100506274 , was associated with lower systolic blood pressure in stage-1 (beta [SE]=−32.6 [6.0]; P =4.99×10 -8 ) but not stage-2 analysis ( P =0.11). Furthermore, a novel common variant at the known INSR locus was suggestively associated with diastolic blood pressure in stage-1 (beta [SE]=−0.36 [0.07]; P =4.18×10 -7 ) and attained genome-wide significance in stage-2 (beta [SE]=−0.29 [0.03]; P =7.28×10 -23 ). Nineteen additional signals suggestively associated with blood pressure in meta-analysis of single and aggregate rare variant findings ( P -6 and P -4 , respectively). Discussion: We report one promising but unconfirmed rare variant for blood pressure and, more importantly, contribute insights for future blood pressure sequencing studies. Our findings suggest promise of aggregate analyses to complement single variant analysis strategies and the need for larger, diverse samples, and family studies to enable robust rare variant identification.

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.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.Mean BMI increased from 32.1 to 33.5 kg/mIn Samoans, the minor A allele of rs373863828 is associated with an increased rate-of-change in BMI in females. On average, BMI of females with the AA genotype increased 0.30 kg/m

Background: Plasma proteins are critical mediators of cardiovascular processes and are the targets of many drugs. Previous efforts to characterize the genetic architecture of the plasma proteome have been limited by a focus on individuals of European descent and leveraged genotyping arrays and imputation. Here we describe whole genome sequence analysis of the plasma proteome in individuals with greater African ancestry, increasing our power to identify novel genetic determinants. Methods: Proteomic profiling of 1301 proteins was performed in 1852 Black adults from the Jackson Heart Study using aptamer-based proteomics (SomaScan). Whole genome sequencing association analysis was ascertained for all variants with minor allele count ≥5. Results were validated using an alternative, antibody-based, proteomic platform (Olink) as well as replicated in the Multi-Ethnic Study of Atherosclerosis and the HERITAGE Family Study (Health, Risk Factors, Exercise Training and Genetics). Results: We identify 569 genetic associations between 479 proteins and 438 unique genetic regions at a Bonferroni-adjusted significance level of 3.8×10 -11 . These associations include 114 novel locus-protein relationships and an additional 217 novel sentinel variant-protein relationships. Novel cardiovascular findings include new protein associations at the APOE gene locus including ZAP70 (sentinel single nucleotide polymorphism [SNP] rs7412-T, β=0.61±0.05, P =3.27×10 -30 ) and MMP-3 (β=-0.60±0.05, P =1.67×10 -32 ), as well as a completely novel pleiotropic locus at the HPX gene, associated with 9 proteins. Further, the associations suggest new mechanisms of genetically mediated cardiovascular disease linked to African ancestry; we identify a novel association between variants linked to APOL1-associated chronic kidney and heart disease and the protein CKAP2 (rs73885319-G, β=0.34±0.04, P =1.34×10 -17 ) as well as an association between ATTR amyloidosis and RBP4 levels in community-dwelling individuals without heart failure. Conclusions: Taken together, these results provide evidence for the functional importance of variants in non-European populations, and suggest new biological mechanisms for ancestry-specific determinants of lipids, coagulation, and myocardial function.

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.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).Equation-derived fat mass was strongly correlated with DXA-derived fat mass for Equation (1) (rThe Polynesian-specific body fat estimation equations remained stable despite changing social and environmental factors and marked increase in obesity prevalence in Samoa.

Motivation Allele-specific differences in molecular traits can be obtained from next-generation sequencing data and could potentially improve testing power, but such information is usually overlooked in association studies. Furthermore, the variation of molecular quantitative traits (e.g. gene expression) could result from the interaction effect of genotypes and phenotypes, but it is challenging to identify such interaction signals in complex disease studies in humans due to small genetic effect sizes and/or small sample sizes. Results We develop a novel statistical method, the combined haplotype interaction test (CHIT), which tests for association between molecular quantitative traits and phenotype–genotype interactions by modeling the total read counts and allele-specific reads in a target region. CHIT can be used as a supplementary analysis to the regular linear interaction regression. In our simulations, CHIT obtains non-inflated type I error rates, and it has higher power than a standard interaction quantitative trait locus approach based on linear regression models. Finally, we illustrate CHIT by testing associations between gene expression obtained by RNA-seq and the interaction of SNPs and atopy status from a study of childhood asthma in Puerto Ricans, and results demonstrate that CHIT could be more powerful than a standard linear interaction expression quantitative trait loci approach. Availability and implementation The CHIT algorithm has been implemented in Python. The source code and documentation are available and can be downloaded from https://github.com/QiYanPitt/CHIT. Supplementary information Supplementary data are available at Bioinformatics online.

Epigenomic modifications of the brain-derived neurotrophic factor (BDNF) gene have been postulated to underlie the pathogenesis of neurodevelopmental, psychiatric, and neurological conditions. This systematic review summarizes current evidence investigating the association of BDNF epigenomic modifications (DNA methylation, non-coding RNA, histone modifications) with brain-related phenotypes in humans. A novel contribution is our creation of an open access web-based application, the BDNF DNA Methylation Map, to interactively visualize specific positions of CpG sites investigated across all studies for which relevant data were available. Our literature search of four databases through September 27, 2021 returned 1,701 articles, of which 153 met inclusion criteria. Our review revealed exceptional heterogeneity in methodological approaches, hindering the identification of clear patterns of robust and/or replicated results. We summarize key findings and provide recommendations for future epigenomic research. The existing literature appears to remain in its infancy and requires additional rigorous research to fulfill its potential to explain BDNF-linked risk for brain-related conditions and improve our understanding of the molecular mechanisms underlying their pathogenesis.Highlights153 articles examined BDNF epigenomic modifications and brain-related phenotypesNovel BDNF DNA Methylation Map allows users to interactively visualize CpGsBDNF epigenomics lack robust/replicated results due to methodological heterogeneity

BackgroundThe Pacific-specific minor allele of rs373863828, a missense variant in CREB3 Regulatory Factor (CREBRF), is associated with several cardiometabolic phenotypes in Polynesian peoples, but the variant’s function remains poorly understood. To broaden our understanding of this variant, we used joint multivariate and network analyses to examine the relationships between rs373863828 and a panel of correlated anthropometric and lipid phenotypes.MethodsWe tested the association of rs373863828 with a panel of phenotypes (body mass index [BMI], weight, height, HDL cholesterol, triglycerides, and total cholesterol) under a multivariate Bayesian association model in a cohort from Samoa (N = 1 632), a Māori and Pacific Island (Polynesian) cohort from Aotearoa New Zealand (N = 1 419), and the combined cohort (N = 2 976). An expanded set of phenotypes (adding estimated fat and fat-free mass, abdominal circumference, hip circumference, and abdominal-hip ratio) was also tested in the Samoa cohort (N = 1 496). Bayesian networks were learned to further understand the structure of the relationships.ResultsIn the Samoa cohort, significant associations (log10 Bayes factor ≥5.0) were found between rs373863828 and the overall phenotype panel (7.97), weight (8.35) and BMI (6.39). In the Aotearoa New Zealand cohort, suggestive associations (log10 Bayes factor between 1.5 and 5) were found between rs373863828 and the overall phenotype panel (3.64), weight (3.30), and BMI (1.79). In the combined cohort, concordant signals with stronger magnitudes were observed. In the expanded phenotype analyses among the Samoa cohort, significant associations were also observed between rs373863828 and fat mass (5.68), abdominal circumference (5.37), and hip circumference (5.15).Bayesian networks provided evidence for a direct association of rs373863828 with weight and indirect associations with height and BMI.ConclusionsWhen correlation structures were considered, multivariate Bayesian analyses provided additional evidence of rs373863828’s pleiotropic effects and highlighted a strong direct effect only on weight.

Scientific data visualization is a critical aspect of fully understanding data patterns and trends. To date, the majority of data visualizations in nursing research – as with other biomedical fields – have been static. The availability of electronic scientific journal articles (which are quickly becoming the norm) has created new opportunities for dynamic and interactive data visualization which carry added cognitive benefits and support the ability to understand data more fully. Therefore, here we highlight the benefits of R, an open-source programming language, for scientific data visualization, with a specific focus on creating dynamic, interactive figures using the R shiny package. For R users, we have included a tutorial with example code to create three increasingly complex shiny applications. For individuals more interested in understanding the potential of R shiny as an innovative tool to interact with research data, we have included links to online versions of the examples that do not require any programming or R experience. We believe that widespread adoption of dynamic and interactive scientific data visualization will further support nurse scientists’ higher-level mission of advancing our understanding of health and wellness of individuals and communities.

Background. Epigenetic biomarkers have the potential to explain outcome heterogeneity following traumatic brain injury (TBI) but are largely unexplored. Objective. This exploratory pilot study characterized brain-derived neurotrophic factor ( BDNF) DNA methylation trajectories following severe TBI. Methods. Brain-derived neurotrophic factor DNA methylation trajectories in cerebrospinal fluid (CSF) over the first 5 days following severe TBI in 112 adults were examined in association with 3- and 12-month outcomes. Results. Group-based trajectory analysis revealed low and high DNA methylation groups at two BDNF cytosine-phosphate-guanine (CpG) targets that showed suggestive associations ( P < .05) with outcomes. Membership in the high DNA methylation groups was associated with better outcomes after controlling for age, sex, and injury severity. Associations of age × trajectory group interactions with outcomes at a third CpG site revealed a pattern of the same or better outcomes with higher ages in the high DNA methylation group and worse outcomes with higher ages in the low DNA methylation group. Conclusions. Although no observed associations met the empirical significance threshold after correcting for multiple comparisons, suggestive associations of the main effect models were consistent in their direction of effect and were observed across two CpG sites and two outcome time points. Results suggest that higher acute CSF BDNF DNA methylation may promote recovery following severe TBI in adults, and this effect may be more robust with higher age. While the results require replication in larger and racially diverse independent samples, BDNF DNA methylation may serve as an early postinjury biomarker helping to explain outcome heterogeneity following TBI.

Preeclampsia is a leading cause of maternal and neonatal morbidity and mortality. Chronological age and race are associated with preeclampsia, but the role of these factors is not entirely understood. We hypothesized that DNA methylation age, a measure of biological age, would be higher in individuals with preeclampsia than in individuals with normotensive pregnancy and that DNA methylation age would differ by race across pregnancy. This was a longitudinal, exploratory study of 56 pregnant individuals (n = 28 preeclampsia cases and n = 28 normotensive controls). Genome-wide DNA methylation data were generated from trimester-specific peripheral blood samples. DNA methylation age was estimated using the "Improved Precision" clock, and ∆age, the difference between DNA methylation age and chronological age, was computed. DNA methylation age was compared with chronological age using Pearson correlations. The relationships between ∆age and preeclampsia status, self-reported race, and covariates were tested using multiple linear regression and performed both with and without consideration of cell-type heterogeneity. We observed strong correlation between chronological age and DNA methylation age across pregnancy, with significantly stronger correlation observed in White participants than in Black participants. We observed no association between ∆age and preeclampsia status. However, ∆age was higher in participants with higher pre-pregnancy body mass index in trimester 1 and lower in Black participants than in White participants in trimesters 2 and 3. Observations were largely consistent when controlling for cell-type heterogeneity. Our findings in a small sample support the need for additional studies to investigate the relationship between race and biological age, which could provide further insight into racial disparities across pregnancy. However, this study does not support an association between ∆age and preeclampsia status.

The Trans-Omics for Precision Medicine (TOPMed) programme seeks to elucidate the genetic architecture and biology of heart, lung, blood and sleep disorders, with the ultimate goal of improving diagnosis, treatment and prevention of these diseases. The initial phases of the programme focused on whole-genome sequencing of individuals with rich phenotypic data and diverse backgrounds. Here we describe the TOPMed goals and design as well as the available resources and early insights obtained from the sequence data. The resources include a variant browser, a genotype imputation server, and genomic and phenotypic data that are available through dbGaP (Database of Genotypes and Phenotypes)1. In the first 53,831 TOPMed samples, we detected more than 400 million single-nucleotide and insertion or deletion variants after alignment with the reference genome. Additional previously undescribed variants were detected through assembly of unmapped reads and customized analysis in highly variable loci. Among the more than 400 million detected variants, 97% have frequencies of less than 1% and 46% are singletons that are present in only one individual (53% among unrelated individuals). These rare variants provide insights into mutational processes and recent human evolutionary history. The extensive catalogue of genetic variation in TOPMed studies provides unique opportunities for exploring the contributions of rare and noncoding sequence variants to phenotypic variation. Furthermore, combining TOPMed haplotypes with modern imputation methods improves the power and reach of genome-wide association studies to include variants down to a frequency of approximately 0.01%., The goals, resources and design of the NHLBI Trans-Omics for Precision Medicine (TOPMed) programme are described, and analyses of rare variants detected in the first 53,831 samples provide insights into mutational processes and recent human evolutionary history.

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.

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 (TG). 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, pHDL-C = 7.63 × 10−10; βTG = 12.00 mg/dL, pTG = 3.82 × 10−7). While this variant appears to be Polynesian-specific, there is also evidence of association from other multi-ancestry 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.

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.652CT, p.R218∗; posterior probability = 0.9987) in

International audience; 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.

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

BACKGROUND/OBJECTIVES: In Samoa, where 80% of the adult population is living with obesity, understanding the determinants of adiposity and growth during infancy may inform prevention efforts. We examined the association of a missense variant, rs373863828, in the CREBRF gene with body composition in Samoan infants. Adults with one or more copies of the rs373863828 minor allele (A) have higher odds of obesity, based on body-mass index (BMI), but paradoxically decreased odds of diabetes compared to those without the allele. Our study may offer novel insight into the natural history and pathogenesis of this unexpected relationship. SUBJECTS/METHODS: In a prospective study, we measured body composition in early infancy, and at 2- and 4-months of age using anthropometry and Dual-Energy X-Ray Absorptiometry (DXA). We genotyped subjects at the CREBRF rs373863828 locus and compared infants with (AA/AG) and without (GG) the variant. In longitudinal analyses, we calculated the absolute change in each outcome from the early infant to the 4-month assessment, adjusting for baseline and other covariates. RESULTS: In cross-sectional analyses, there was no significant difference in infant BMI or fat mass by genotype. After adjusting for covariates, infants with the variant had 4.0 ± 1.8g more bone mass (p=0.026) and 210.9 ±79.6g more lean mass (p=0.009) at 4-months and accumulated 176.9 ± 73.0g more lean mass between the early infant and 4-month assessment (p=0.017). CONCLUSIONS: The CREBRF rs373863828 minor allele (A) was not associated with increased BMI or adiposity in Samoan infants, but instead with increased lean and bone mass. Our findings suggest that lean (i.e. muscle) and bone mass accretion should be explored as pathways to explain the ‘protective’ effect of the CREBRF variant against diabetes.

In a transcription-wide association study of nasal epithelium, we identify novel and previously reported susceptibility genes for atopic asthma in children and show that gene co-expression networks differ markedly between children with and without atopic asthma.

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.

Significance There are multiple archaeological debates regarding early Samoan population history. Here, we add genetic data to this discussion, which supports an initial small population size at the founding of the Samoan islands. We also indicate a major demographic change approximately 1,000 y ago that mirrors the archaeological record. In addition, we demonstrate the utility of rare genetic variants in identifying sparse population structure. These genetic results help establish a detailed demographic model for the Samoan population, which will aid in future studies of Oceanic populations for both history and disease., 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.

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 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 (CV) 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 between 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 between 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.

Identifying high-risk individuals with diseases through reliable prediction models guides screening and preventive treatment. Most complex diseases have a genetic basis influenced by multiple genes and so disease risk can be estimated using polygenic risk score (PRS) algorithms. Many PRS algorithms have been developed so far. Among them, BayesR shows good characteristics of unbiasedness, accuracy, sparseness, and robustness. It detects the associated SNPs, estimates the SNP effects, and makes prediction of disease risks based on all SNPs simultaneously. However, this method assumes that the phenotypes follow a Gaussian distribution, which cannot be met in case-control studies. Here, we made an extension of the BayesR method, called BayesRB, by adding auxiliary variables to the BayesR model. We explored the characteristics, efficacy, and accuracy of BayesRB when estimating SNP effects and predicting disease risks compared with three traditional algorithms under different conditions using both simulated data and real data from the Welcome Trust Case Control Consortium (WTCCC). For SNP effect estimation, BayesRB shows unbiasedness and sparseness for big and small effect SNPs, respectively. For disease risk prediction, BayesRB had the best performance among the methods. This study provides a theoretical basis for complex disease risk prediction and disease prevention

Background Following aneurysmal subarachnoid hemorrhage (aSAH), the brain is susceptible to ferroptosis, a type of iron-dependent cell death. Therapeutic intervention targeting the iron homeostasis pathway shows promise for mitigating ferroptosis and improving recovery in animal models, but little work has been conducted in humans. DNA methylation (DNAm) plays a key role in gene expression and brain function, plasticity, and injury recovery, making it a potentially useful biomarker of outcomes or therapeutic target for intervention. Therefore, in this longitudinal, observational study, we examined the relationships between trajectories of DNAm in candidate genes related to iron homeostasis and acute (cerebral vasospasm and delayed cerebral ischemia) and long-term (Glasgow Outcome Scale [GOS, unfavorable = 1–3] and death) patient outcomes after aSAH. Results Longitudinal, genome-wide DNAm data were generated from DNA extracted from post-aSAH cerebrospinal fluid (n = 260 participants). DNAm trajectories of 637 CpG sites in 36 candidate genes related to iron homeostasis were characterized over 13 days post-aSAH using group-based trajectory analysis, an unsupervised clustering method. Significant associations were identified between inferred DNAm trajectory groups at several CpG sites and acute and long-term outcomes. Among our results, cg25713625 in the STEAP3 metalloreductase gene (STEAP3) stood out. Specifically, in comparing the highest cg25713625 DNAm trajectory group with the lowest, we observed significant associations (i.e., based on p-values less than an empirical significance threshold) with unfavorable GOS at 3 and 12 months (OR = 11.7, p = 0.0006 and OR = 15.6, p = 0.0018, respectively) and death at 3 and 12 months (OR = 19.1, p = 0.0093 and OR = 12.8, p = 0.0041, respectively). These results were replicated in an independent sample (n = 100 participants) observing significant associations with GOS at 3 and 12 months (OR = 8.2, p = 0.001 and OR = 6.3, p = 0.0.0047, respectively) and death at 3 months (OR = 2.3, p = 0.008) and a suggestive association (i.e., p-value OR = 2.0, p = 0.0272). In both samples, an additive effect of the DNAm trajectory group was observed as the percentage of participants with unfavorable long-term outcomes increased substantially with higher DNAm trajectory groups. Conclusion Our results support a role for DNAm of cg25713625/STEAP3 in recovery following aSAH. Additional research is needed to further explore the role of DNAm of cg25713625/STEAP3 as a biomarker of unfavorable outcomes, or therapeutic target to improve outcomes, to translate these findings clinically.