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.

Exploratory assessment of the effect of systemic administration of soluble glycoprotein 130 on cognitive performance and chemokine levels in a mouse model of experimental traumatic brain injury.

Uncontrolled neuroinflammation mediates traumatic brain injury (TBI) pathology and impairs recovery. Interleukin-6 (IL-6), a pleiotropic inflammatory regulator, is associated with poor clinical TBI outcomes. IL-6 operates via classical-signaling through membrane-bound IL-6 receptor (IL-6R) and trans-signaling through soluble IL-6 receptor (s)IL-6R. IL-6 trans-signaling specifically contributes to neuropathology, making it a potential precision therapeutic TBI target. Soluble glycoprotein 130 (sgp130) prevents IL-6 trans-signaling, sparing classical signaling, thus is a possible treatment. Mice received either controlled cortical impact (CCI) (6.0 ± 0.2 m/s; 2 mm; 50-60ms) or sham procedures. Vehicle (VEH) or sgp130-Fc was subcutaneously administered to sham (VEH or 1 µg) and CCI (VEH, 0.25 µg or 1 µg) mice on days 1, 4, 7, 10 and 13 post-surgery to assess effects on cognition [Morris Water Maze (MWM)] and ipsilateral hemisphere IL-6 related biomarkers (day 21 post-surgery). CCI + sgp130-Fc groups (0.25 µg and 1 µg) were combined for analysis given similar behavior/biomarker outcomes. CCI + VEH mice had longer latencies and path lengths to the platform and increased peripheral zone time versus Sham + VEH and Sham + sgp130-Fc mice, suggesting injury-induced impairments in learning and anxiety. CCI + sgp130-Fc mice had shorter platform latencies and path lengths and had decreased peripheral zone time, indicating a therapeutic benefit of sgp130-Fc after injury on learning and anxiety. Interestingly, Sham + sgp130-Fc mice had shorter platform latencies, path lengths and peripheral zone times than Sham + VEH mice, suggesting a beneficial effect of sgp130-Fc, independent of injury. CCI + VEH mice had increased brain IL-6 and decreased sgp130 levels versus Sham + VEH and Sham + sgp130-Fc mice. There was no treatment effect on IL-6, sIL6-R or sgp130 in Sham + VEH versus Sham + sgp130-Fc mice. There was also no treatment effect on IL-6 in CCI + VEH versus CCI + sgp130-Fc mice. However, CCI + sgp130-Fc mice had increased sIL-6R and sgp130 versus CCI + VEH mice, demonstrating sgp130-Fc treatment effects on brain biomarkers. Inflammatory chemokines (MIP-1ß, IP-10, MIG) were increased in CCI + VEH mice versus Sham + VEH and Sham + sgp130-Fc mice. However, CCI + sgp130-Fc mice had decreased chemokine levels versus CCI + VEH mice. IL-6 positively correlated, while sgp130 negatively correlated, with chemokine levels. Overall, we found that systemic sgp130-Fc treatment after CCI improved learning, decreased anxiety and reduced CCI-induced brain chemokines. Future studies will explore sex-specific dosing and treatment mechanisms for sgp130-Fc therapy.

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.

Genome-wide Interaction Study Implicates VGLL2 and Alcohol Exposure and PRL and Smoking in Orofacial Cleft Risk.

Non-syndromic cleft lip with or without cleft palate (NSCL/P) is a common birth defect, affecting approximately 1 in 700 births. NSCL/P has complex etiology including several known genes and environmental factors; however, known genetic risk variants only account for a small fraction of the heritability of NSCL/P. It is commonly suggested that gene-by-environment (G×E) interactions may help explain some of the "missing" heritability of NSCL/P. We conducted a genome-wide G×E interaction study in cases and controls of European ancestry with three common maternal exposures during pregnancy: alcohol, smoking, and vitamin use using a two-stage design. After selecting 127 loci with suggestive 2df tests for gene and G x E effects, 40 loci showed significant G x E effects after correcting for multiple tests. Notable interactions included SNPs of 6q22 near VGLL2 with alcohol and 6p22.3 near PRL with smoking. These interactions could provide new insights into the etiology of CL/P and new opportunities to modify risk through behavioral changes.

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.

FAT4 identified as a potential modifier of orofacial cleft laterality.

Orofacial clefts (OFCs) are common (1 in 700 births) congenital malformations that include a cleft lip (CL) and cleft lip and palate (CLP). These OFC subtypes are also heterogeneous themselves, with the CL occurring on the left, right, or both sides of the upper lip. Unilateral CL and CLP have a 2:1 bias towards left-sided clefts, suggesting a nonrandom process. Here, we performed a study of left- and right-sided clefts within the CL and CLP subtypes to better understand the genetic factors controlling cleft laterality. We conducted genome-wide modifier analyses by comparing cases that had right unilateral CL (RCL; N = 130), left unilateral CL (LCL; N = 216), right unilateral CLP (RCLP; N = 416), or left unilateral CLP (LCLP; N = 638), and identified a candidate region on 4q28, 400 kb downstream from FAT4, that approached genome-wide significance for LCL versus RCL (p = 8.4 × 10-8 ). Consistent with its potential involvement as a genetic modifier of CL, we found that Fat4 exhibits a specific domain of expression in the mesenchyme of the medial nasal processes that form the median upper lip. Overall, these results suggest that the epidemiological similarities in left- to right-sided clefts in CL and CLP are not reflected in the genetic association results.

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

BACKGROUND: The prevalence of obesity and diabetes in Samoa, like many other Pacific Island nations, has reached epidemic proportions. Although the etiology of these conditions can be largely attributed to the rapidly changing economic and nutritional environment, a recently identified genetic variant, rs373863828 (CREB 3 regulatory factor, CREBRF: c.1370G>A p.[R457Q]) is associated with increased odds of obesity, but paradoxically, decreased odds of diabetes. OBJECTIVE: The overarching goal of the Soifua Manuia (Good Health) study was to precisely characterize the association of the CREBRF variant with metabolic (body composition and glucose homeostasis) and behavioral traits (dietary intake, physical activity, sleep, and weight control behaviors) that influence energy homeostasis in 500 adults. METHODS: A cohort of adult Samoans who participated in a genome-wide association study of adiposity in Samoa in 2010 was followed up, based on the presence or absence of the CREBRF variant, between August 2017 and March 2019. Over a period of 7-10 days, each participant completed the main study protocol, which consisted of anthropometric measurements (weight, height, circumferences, and skinfolds), body composition assessment (bioelectrical impedance and dual-energy x-ray absorptiometry), point-of-care glycated hemoglobin measurement, a fasting blood draw and oral glucose tolerance test, urine collection, blood pressure measurement, hand grip strength measurement, objective physical activity and sleep apnea monitoring, and questionnaire measures (eg, health interview, cigarette and alcohol use, food frequency questionnaire, socioeconomic position, stress, social support, food and water insecurity, sleep, body image, and dietary preferences). In January 2019, a subsample of the study participants (n=118) completed a buttock fat biopsy procedure to collect subcutaneous adipose tissue samples. RESULTS: Enrollment of 519 participants was completed in March 2019. Data analyses are ongoing, with results expected in 2020 and 2021. CONCLUSIONS: While the genetic variant rs373863828, in CREBRF, has the largest known effect size of any identified common obesity gene, very little is currently understood about the mechanisms by which it confers increased odds of obesity but paradoxically lowered odds of type 2 diabetes. The results of this study will provide insights into how the gene functions on a whole-body level, which could provide novel targets to prevent or treat obesity, diabetes, and associated metabolic disorders. This study represents the human arm of a comprehensive and integrated approach involving humans as well as preclinical models that will provide novel insights into metabolic disease. INTERNATIONAL REGISTERED REPORT IDENTIFIER (IRRID): RR1-10.2196/17329.

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.

Multiethnic GWAS Reveals Polygenic Architecture of Earlobe Attachment.

The genetic basis of earlobe attachment has been a matter of debate since the early 20th century, such that geneticists argue both for and against polygenic inheritance. Recent genetic studies have identified a few loci associated with the trait, but large-scale analyses are still lacking. Here, we performed a genome-wide association study of lobe attachment in a multiethnic sample of 74,660 individuals from four cohorts (three with the trait scored by an expert rater and one with the trait self-reported). Meta-analysis of the three expert-rater-scored cohorts revealed six associated loci harboring numerous candidate genes, including EDAR, SP5, MRPS22, ADGRG6 (GPR126), KIAA1217, and PAX9. The large self-reported 23andMe cohort recapitulated each of these six loci. Moreover, meta-analysis across all four cohorts revealed a total of 49 significant (p < 5 × 10-8) loci. Annotation and enrichment analyses of these 49 loci showed strong evidence of genes involved in ear development and syndromes with auricular phenotypes. RNA sequencing data from both human fetal ear and mouse second branchial arch tissue confirmed that genes located among associated loci showed evidence of expression. These results provide strong evidence for the polygenic nature of earlobe attachment and offer insights into the biological basis of normal and abnormal ear development.

Genome-wide meta-analyses of nonsyndromic orofacial clefts identify novel associations between FOXE1 and all orofacial clefts, and TP63 and cleft lip with or without cleft palate.

Nonsyndromic orofacial clefts (OFCs) are a heterogeneous group of common craniofacial birth defects with complex etiologies that include genetic and environmental risk factors. OFCs are commonly categorized as cleft lip with or without cleft palate (CL/P) and cleft palate alone (CP), which have historically been analyzed as distinct entities. Genes for both CL/P and CP have been identified via multiple genome-wide linkage and association studies (GWAS); however, altogether, known variants account for a minority of the estimated heritability in risk to these craniofacial birth defects. We performed genome-wide meta-analyses of CL/P, CP, and all OFCs across two large, multiethnic studies. We then performed population-specific meta-analyses in sub-samples of Asian and European ancestry. In addition to observing associations with known variants, we identified a novel genome-wide significant association between SNPs located in an intronic TP63 enhancer and CL/P (p = 1.16 × 10-8). Several novel loci with compelling candidate genes approached genome-wide significance on 4q21.1 (SHROOM3), 12q13.13 (KRT18), and 8p21 (NRG1). In the analysis of all OFCs combined, SNPs near FOXE1 reached genome-wide significance (p = 1.33 × 10-9). Our results support the highly heterogeneous nature of OFCs and illustrate the utility of meta-analysis for discovering new genetic risk factors.

A multi-ethnic genome-wide association study identifies novel loci for non-syndromic cleft lip with or without cleft palate on 2p24.2, 17q23 and 19q13.

Orofacial clefts (OFCs), which include non-syndromic cleft lip with or without cleft palate (CL/P), are among the most common birth defects in humans, affecting approximately 1 in 700 newborns. CL/P is phenotypically heterogeneous and has a complex etiology caused by genetic and environmental factors. Previous genome-wide association studies (GWASs) have identified at least 15 risk loci for CL/P. As these loci do not account for all of the genetic variance of CL/P, we hypothesized the existence of additional risk loci. We conducted a multiethnic GWAS in 6480 participants (823 unrelated cases, 1700 unrelated controls and 1319 case-parent trios) with European, Asian, African and Central and South American ancestry. Our GWAS revealed novel associations on 2p24 near FAM49A, a gene of unknown function (P = 4.22 × 10-8), and 19q13 near RHPN2, a gene involved in organizing the actin cytoskeleton (P = 4.17 × 10-8). Other regions reaching genome-wide significance were 1p36 (PAX7), 1p22 (ARHGAP29), 1q32 (IRF6), 8q24 and 17p13 (NTN1), all reported in previous GWASs. Stratification by ancestry group revealed a novel association with a region on 17q23 (P = 2.92 × 10-8) among individuals with European ancestry. This region included several promising candidates including TANC2, an oncogene required for development, and DCAF7, a scaffolding protein required for craniofacial development. In the Central and South American ancestry group, significant associations with loci previously identified in Asian or European ancestry groups reflected their admixed ancestry. In summary, we have identified novel CL/P risk loci and suggest new genes involved in craniofacial development, confirming the highly heterogeneous etiology of OFCs.

Identification of functional variants for cleft lip with or without cleft palate in or near PAX7, FGFR2, and NOG by targeted sequencing of GWAS loci.

Although genome-wide association studies (GWASs) for nonsyndromic orofacial clefts have identified multiple strongly associated regions, the causal variants are unknown. To address this, we selected 13 regions from GWASs and other studies, performed targeted sequencing in 1,409 Asian and European trios, and carried out a series of statistical and functional analyses. Within a cluster of strongly associated common variants near NOG, we found that one, rs227727, disrupts enhancer activity. We furthermore identified significant clusters of non-coding rare variants near NTN1 and NOG and found several rare coding variants likely to affect protein function, including four nonsense variants in ARHGAP29. We confirmed 48 de novo mutations and, based on best biological evidence available, chose two of these for functional assays. One mutation in PAX7 disrupted the DNA binding of the encoded transcription factor in an in vitro assay. The second, a non-coding mutation, disrupted the activity of a neural crest enhancer downstream of FGFR2 both in vitro and in vivo. This targeted sequencing study provides strong functional evidence implicating several specific variants as primary contributory risk alleles for nonsyndromic clefting in humans.