Nothing about us without us: Sharing results with communities that provide genomic data.

Sharing genetic and other study results with the communities who participate in research falls under benefit-sharing and capacity-building initiatives that underpin a more equitable biomedical research relationship. Yet, which results to return and how remain fundamental challenges that persist in the absence of practical guidance and institutional policies. Here, we discuss how the return of results can be implemented across different geographies, study designs, and project budgets.

Assessment of rosacea symptom severity by genome-wide association study and expression analysis highlights immuno-inflammatory and skin pigmentation genes.

Rosacea is a common, chronic skin disease of variable severity with limited treatment options. The cause of rosacea is unknown, but it is believed to be due to a combination of hereditary and environmental factors. Little is known about the genetics of the disease. We performed a genome-wide association study (GWAS) of rosacea symptom severity with data from 73 265 research participants of European ancestry from the 23andMe customer base. Seven loci had variants associated with rosacea at the genome-wide significance level (P < 5 × 10-8). Further analyses highlighted likely gene regions or effector genes including IRF4 (P = 1.5 × 10-17), a human leukocyte antigen (HLA) region flanked by PSMB9 and HLA-DMB (P = 2.2 × 10-15), HERC2-OCA2 (P = 4.2 × 10-12), SLC45A2 (P = 1.7 × 10-10), IL13 (P = 2.8 × 10-9), a region flanked by NRXN3 and DIO2 (P = 4.1 × 10-9), and a region flanked by OVOL1and SNX32 (P = 1.2 × 10-8). All associations with rosacea were novel except for the HLA locus. Two of these loci (HERC-OCA2 and SLC45A2) and another precedented variant (rs1805007 in melanocortin 1 receptor) with an association P value just below the significance threshold (P = 1.3 × 10-7) have been previously associated with skin phenotypes and pigmentation, two of these loci are linked to immuno-inflammation phenotypes (IL13 and PSMB9-HLA-DMA) and one has been associated with both categories (IRF4). Genes within three loci (PSMB9-HLA-DMA, HERC-OCA2 and NRX3-DIO2) were differentially expressed in a previously published clinical rosacea transcriptomics study that compared lesional to non-lesional samples. The identified loci provide specificity of inflammatory mechanisms in rosacea, and identify potential pathways for therapeutic intervention.

Parent-of-origin-specific allelic associations among 106 genomic loci for age at menarche.

Age at menarche is a marker of timing of puberty in females. It varies widely between individuals, is a heritable trait and is associated with risks for obesity, type 2 diabetes, cardiovascular disease, breast cancer and all-cause mortality. Studies of rare human disorders of puberty and animal models point to a complex hypothalamic-pituitary-hormonal regulation, but the mechanisms that determine pubertal timing and underlie its links to disease risk remain unclear. Here, using genome-wide and custom-genotyping arrays in up to 182,416 women of European descent from 57 studies, we found robust evidence (P < 5 × 10(-8)) for 123 signals at 106 genomic loci associated with age at menarche. Many loci were associated with other pubertal traits in both sexes, and there was substantial overlap with genes implicated in body mass index and various diseases, including rare disorders of puberty. Menarche signals were enriched in imprinted regions, with three loci (DLK1-WDR25, MKRN3-MAGEL2 and KCNK9) demonstrating parent-of-origin-specific associations concordant with known parental expression patterns. Pathway analyses implicated nuclear hormone receptors, particularly retinoic acid and γ-aminobutyric acid-B2 receptor signalling, among novel mechanisms that regulate pubertal timing in humans. Our findings suggest a genetic architecture involving at least hundreds of common variants in the coordinated timing of the pubertal transition.

Exome-chip meta-analysis identifies association between variation in ANKRD26 and platelet aggregation.

Previous genome-wide association studies (GWAS) have identified several variants associated with platelet function phenotypes; however, the proportion of variance explained by the identified variants is mostly small. Rare coding variants, particularly those with high potential for impact on protein structure/function, may have substantial impact on phenotype but are difficult to detect by GWAS. The main purpose of this study was to identify low frequency or rare variants associated with platelet function using genotype data from the Illumina HumanExome Bead Chip. Three family-based cohorts of European ancestry, including ~4,000 total subjects, comprised the discovery cohort and two independent cohorts, one of European and one of African American ancestry, were used for replication. Optical aggregometry in platelet-rich plasma was performed in all the discovery cohorts in response to adenosine diphosphate (ADP), epinephrine, and collagen. Meta-analyses were performed using both gene-based and single nucleotide variant association methods. The gene-based meta-analysis identified a significant association (P = 7.13 × 10-7) between rare genetic variants in ANKRD26 and ADP-induced platelet aggregation. One of the ANKRD26 SNVs - rs191015656, encoding a threonine to isoleucine substitution predicted to alter protein structure/function, was replicated in Europeans. Aggregation increases of ~20-50% were observed in heterozygotes in all cohorts. Novel genetic signals in ABCG1 and HCP5 were also associated with platelet aggregation to ADP in meta-analyses, although only results for HCP5 could be replicated. The SNV in HCP5 intersects epigenetic signatures in CD41+ megakaryocytes suggesting a new functional role in platelet biology for HCP5. This is the first study to use gene-based association methods from SNV array genotypes to identify rare variants related to platelet function. The molecular mechanisms and pathophysiological relevance for the identified genetic associations requires further study.

TM6SF2 rs58542926 impacts lipid processing in liver and small intestine.

The transmembrane 6 superfamily member 2 (TM6SF2) loss-of-function variant rs58542926 is a genetic risk factor for nonalcoholic fatty liver disease and progression to fibrosis but is paradoxically associated with lower levels of hepatically derived triglyceride-rich lipoproteins. TM6SF2 is expressed predominantly in liver and small intestine, sites for triglyceride-rich lipoprotein biogenesis and export. In light of this, we hypothesized that TM6SF2 may exhibit analogous effects on both liver and intestine lipid homeostasis. To test this, we genotyped rs58542926 in 983 bariatric surgery patients from the Geisinger Medical Center for Nutrition and Weight Management, Geisinger Health System, in Pennsylvania and from 3,556 study participants enrolled in the Amish Complex Disease Research Program. Although these two cohorts have different metabolic profiles, carriers in both cohorts had improved fasting lipid profiles. Importantly, following a high-fat challenge, carriers in the Amish Complex Disease Research Program cohort exhibited significantly lower postprandial serum triglycerides, suggestive of a role for TM6SF2 in the small intestine. To gain further insight into this putative role, effects of TM6SF2 deficiency were studied in a zebrafish model and in cultured human Caco-2 enterocytes. In both systems TM6SF2 deficiency resulted in defects in small intestine metabolism in response to dietary lipids, including significantly increased lipid accumulation, decreased lipid clearance, and increased endoplasmic reticulum stress. CONCLUSIONS: These data strongly support a role of TM6SF2 in the regulation of postprandial lipemia, potentially through a similar function for TM6SF2 in the lipidation and/or export of both hepatically and intestinally derived triglyceride-rich lipoproteins. (Hepatology 2017;65:1526-1542).

FTO genotype is associated with phenotypic variability of body mass index.

There is evidence across several species for genetic control of phenotypic variation of complex traits, such that the variance among phenotypes is genotype dependent. Understanding genetic control of variability is important in evolutionary biology, agricultural selection programmes and human medicine, yet for complex traits, no individual genetic variants associated with variance, as opposed to the mean, have been identified. Here we perform a meta-analysis of genome-wide association studies of phenotypic variation using ∼170,000 samples on height and body mass index (BMI) in human populations. We report evidence that the single nucleotide polymorphism (SNP) rs7202116 at the FTO gene locus, which is known to be associated with obesity (as measured by mean BMI for each rs7202116 genotype), is also associated with phenotypic variability. We show that the results are not due to scale effects or other artefacts, and find no other experiment-wise significant evidence for effects on variability, either at loci other than FTO for BMI or at any locus for height. The difference in variance for BMI among individuals with opposite homozygous genotypes at the FTO locus is approximately 7%, corresponding to a difference of ∼0.5 kilograms in the standard deviation of weight. Our results indicate that genetic variants can be discovered that are associated with variability, and that between-person variability in obesity can partly be explained by the genotype at the FTO locus. The results are consistent with reported FTO by environment interactions for BMI, possibly mediated by DNA methylation. Our BMI results for other SNPs and our height results for all SNPs suggest that most genetic variants, including those that influence mean height or mean BMI, are not associated with phenotypic variance, or that their effects on variability are too small to detect even with samples sizes greater than 100,000.

Genome-wide analysis of clopidogrel active metabolite levels identifies novel variants that influence antiplatelet response.

Clopidogrel is one of the most commonly used therapeutics for the secondary prevention of cardiovascular events in patients with acute coronary syndromes. However, considerable interindividual variation in clopidogrel response has been documented, resulting in suboptimal therapy and an increased risk of recurrent events for some patients. In this investigation, we carried out the first genome-wide association study of circulating clopidogrel active metabolite levels in 513 healthy participants to directly measure clopidogrel pharmacokinetics. We observed that the CYP2C19 locus was the strongest genetic determinant of active metabolite formation (P=9.5×10). In addition, we identified novel genome-wide significant variants on chromosomes 3p25 (rs187941554, P=3.3×10) and 17q11 (rs80343429, P=1.3×10), as well as six additional loci that showed suggestive evidence of association (P≤1.0×10). Four of these loci showed nominal associations with on-clopidogrel ADP-stimulated platelet aggregation (P≤0.05). Evaluation of clopidogrel active metabolite concentration may help identify novel genetic determinants of clopidogrel response, which has implications for the development of novel therapeutics and improved antiplatelet treatment for at-risk patients in the future.

Null mutation in hormone-sensitive lipase gene and risk of type 2 diabetes.

BACKGROUND: Lipolysis regulates energy homeostasis through the hydrolysis of intracellular triglycerides and the release of fatty acids for use as energy substrates or lipid mediators in cellular processes. Genes encoding proteins that regulate energy homeostasis through lipolysis are thus likely to play an important role in determining susceptibility to metabolic disorders. METHODS: We sequenced 12 lipolytic-pathway genes in Old Order Amish participants whose fasting serum triglyceride levels were at the extremes of the distribution and identified a novel 19-bp frameshift deletion in exon 9 of LIPE, encoding hormone-sensitive lipase (HSL), a key enzyme for lipolysis. We genotyped the deletion in DNA from 2738 Amish participants and performed association analyses to determine the effects of the deletion on metabolic traits. We also obtained biopsy specimens of abdominal subcutaneous adipose tissue from 2 study participants who were homozygous for the deletion (DD genotype), 10 who were heterozygous (ID genotype), and 7 who were noncarriers (II genotype) for assessment of adipose histologic characteristics, lipolysis, enzyme activity, cytokine release, and messenger RNA (mRNA) and protein levels. RESULTS: Carriers of the mutation had dyslipidemia, hepatic steatosis, systemic insulin resistance, and diabetes. In adipose tissue from study participants with the DD genotype, the mutation resulted in the absence of HSL protein, small adipocytes, impaired lipolysis, insulin resistance, and inflammation. Transcription factors responsive to peroxisome-proliferator-activated receptor γ (PPAR-γ) and downstream target genes were down-regulated in adipose tissue from participants with the DD genotype, altering the regulation of pathways influencing adipogenesis, insulin sensitivity, and lipid metabolism. CONCLUSIONS: These findings indicate the physiological significance of HSL in adipocyte function and the regulation of systemic lipid and glucose homeostasis and underscore the severe metabolic consequences of impaired lipolysis. (Funded by the National Institutes of Health and others).

Multistage genome-wide association meta-analyses identified two new loci for bone mineral density.

Aiming to identify novel genetic variants and to confirm previously identified genetic variants associated with bone mineral density (BMD), we conducted a three-stage genome-wide association (GWA) meta-analysis in 27 061 study subjects. Stage 1 meta-analyzed seven GWA samples and 11 140 subjects for BMDs at the lumbar spine, hip and femoral neck, followed by a Stage 2 in silico replication of 33 SNPs in 9258 subjects, and by a Stage 3 de novo validation of three SNPs in 6663 subjects. Combining evidence from all the stages, we have identified two novel loci that have not been reported previously at the genome-wide significance (GWS; 5.0 × 10(-8)) level: 14q24.2 (rs227425, P-value 3.98 × 10(-13), SMOC1) in the combined sample of males and females and 21q22.13 (rs170183, P-value 4.15 × 10(-9), CLDN14) in the female-specific sample. The two newly identified SNPs were also significant in the GEnetic Factors for OSteoporosis consortium (GEFOS, n = 32 960) summary results. We have also independently confirmed 13 previously reported loci at the GWS level: 1p36.12 (ZBTB40), 1p31.3 (GPR177), 4p16.3 (FGFRL1), 4q22.1 (MEPE), 5q14.3 (MEF2C), 6q25.1 (C6orf97, ESR1), 7q21.3 (FLJ42280, SHFM1), 7q31.31 (FAM3C, WNT16), 8q24.12 (TNFRSF11B), 11p15.3 (SOX6), 11q13.4 (LRP5), 13q14.11 (AKAP11) and 16q24 (FOXL1). Gene expression analysis in osteogenic cells implied potential functional association of the two candidate genes (SMOC1 and CLDN14) in bone metabolism. Our findings independently confirm previously identified biological pathways underlying bone metabolism and contribute to the discovery of novel pathways, thus providing valuable insights into the intervention and treatment of osteoporosis.

AP2 suppresses osteoblast differentiation and mineralization through down-regulation of Frizzled-1.

Transcription factor activating protein 2 (AP2) plays an important role in cellular differentiation. Although profound craniofacial and long bone developmental abnormalities have been observed in AP2-knockout mice, the molecular effects of AP2 on osteoblasts are poorly defined. We demonstrated that AP2 regulates the expression of human Frizzled 1 (FZD1), a co-receptor for the Wnt signalling pathway, in human osteoblast cell lines and primary bone marrow stromal cells (BMSCs). We also identified a putative AP2-binding site in the FZD1 proximal promoter in silico and characterized this binding element further in Saos2 in vitro by ChIP, electrophoretic mobility shift and promoter reporter assays. The transcriptional repression of the FZD1 promoter by AP2 was confirmed in normal human fetal osteoblasts (hFOB). Furthermore, overexpression of AP2 resulted in a significant reduction in both differentiation and mineralization of Saos2 cells. Knockdown of FZD1 expression before AP2 up-regulation diminished the AP2-dependent inhibition of Saos2 cell differentiation and mineralization. Similarly, overexpressing FZD1 before AP2 treatment in both Saos2 and BMSCs diminished the inhibitory effect of AP2 on osteoblast differentiation and mineralization. Taken together, these results demonstrate that AP2 is a negative regulator of osteoblast differentiation and mineralization, and its inhibitory effect may be mediated in part through down-regulation of FZD1 expression.

WNT16 influences bone mineral density, cortical bone thickness, bone strength, and osteoporotic fracture risk.

We aimed to identify genetic variants associated with cortical bone thickness (CBT) and bone mineral density (BMD) by performing two separate genome-wide association study (GWAS) meta-analyses for CBT in 3 cohorts comprising 5,878 European subjects and for BMD in 5 cohorts comprising 5,672 individuals. We then assessed selected single-nucleotide polymorphisms (SNPs) for osteoporotic fracture in 2,023 cases and 3,740 controls. Association with CBT and forearm BMD was tested for ∼2.5 million SNPs in each cohort separately, and results were meta-analyzed using fixed effect meta-analysis. We identified a missense SNP (Thr>Ile; rs2707466) located in the WNT16 gene (7q31), associated with CBT (effect size of -0.11 standard deviations [SD] per C allele, P = 6.2 × 10(-9)). This SNP, as well as another nonsynonymous SNP rs2908004 (Gly>Arg), also had genome-wide significant association with forearm BMD (-0.14 SD per C allele, P = 2.3 × 10(-12), and -0.16 SD per G allele, P = 1.2 × 10(-15), respectively). Four genome-wide significant SNPs arising from BMD meta-analysis were tested for association with forearm fracture. SNP rs7776725 in FAM3C, a gene adjacent to WNT16, was associated with a genome-wide significant increased risk of forearm fracture (OR = 1.33, P = 7.3 × 10(-9)), with genome-wide suggestive signals from the two missense variants in WNT16 (rs2908004: OR = 1.22, P = 4.9 × 10(-6) and rs2707466: OR = 1.22, P = 7.2 × 10(-6)). We next generated a homozygous mouse with targeted disruption of Wnt16. Female Wnt16(-/-) mice had 27% (P<0.001) thinner cortical bones at the femur midshaft, and bone strength measures were reduced between 43%-61% (6.5 × 10(-13)

Characterization of European ancestry nonalcoholic fatty liver disease-associated variants in individuals of African and Hispanic descent.

UNLABELLED: Nonalcoholic fatty liver disease (NAFLD) is an obesity-related condition affecting over 50% of individuals in some populations and is expected to become the number one cause of liver disease worldwide by 2020. Common, robustly associated genetic variants in/near five genes were identified for hepatic steatosis, a quantifiable component of NAFLD, in European ancestry individuals. Here we tested whether these variants were associated with hepatic steatosis in African- and/or Hispanic-Americans and fine-mapped the observed association signals. We measured hepatic steatosis using computed tomography in five African American (n = 3,124) and one Hispanic American (n = 849) cohorts. All analyses controlled for variation in age, age(2) , gender, alcoholic drinks, and population substructure. Heritability of hepatic steatosis was estimated in three cohorts. Variants in/near PNPLA3, NCAN, LYPLAL1, GCKR, and PPP1R3B were tested for association with hepatic steatosis using a regression framework in each cohort and meta-analyzed. Fine-mapping across African American cohorts was conducted using meta-analysis. African- and Hispanic-American cohorts were 33.9/37.5% male, with average age of 58.6/42.6 years and body mass index of 31.8/28.9 kg/m(2) , respectively. Hepatic steatosis was 0.20-0.34 heritable in African- and Hispanic-American families (P < 0.02 in each cohort). Variants in or near PNPLA3, NCAN, GCKR, PPP1R3B in African Americans and PNPLA3 and PPP1R3B in Hispanic Americans were significantly associated with hepatic steatosis; however, allele frequency and effect size varied across ancestries. Fine-mapping in African Americans highlighted missense variants at PNPLA3 and GCKR and redefined the association region at LYPLAL1. CONCLUSION: Multiple genetic variants are associated with hepatic steatosis across ancestries. This explains a substantial proportion of the genetic predisposition in African- and Hispanic-Americans. Missense variants in PNPLA3 and GCKR are likely functional across multiple ancestries.

Meta-analysis of genome-wide studies identifies MEF2C SNPs associated with bone mineral density at forearm.

BACKGROUND: Forearm fractures affect 1.7 million individuals worldwide each year and most occur earlier in life than hip fractures. While the heritability of forearm bone mineral density (BMD) and fracture is high, their genetic determinants are largely unknown. AIM: To identify genetic variants associated with forearm BMD and forearm fractures. METHODS: BMD at distal radius, measured by dual-energy x-ray absorptiometry, was tested for association with common genetic variants. We conducted a meta-analysis of genome-wide association studies for BMD in 5866 subjects of European descent and then selected the variants for replication in 715 Mexican American samples. Gene-based association was carried out to supplement the single-nucleotide polymorphism (SNP) association test. We then tested the BMD-associated SNPs for association with forearm fracture in 2023 cases and 3740 controls. RESULTS: We found that five SNPs in the introns of MEF2C were associated with forearm BMD at a genome-wide significance level (p<5×10(-8)) in meta-analysis (lead SNP, rs11951031[T] -0.20 SDs per allele, p=9.01×10(-9)). The gene-based association test suggested an association between MEF2C and forearm BMD (p=0.003). The association between MEF2C variants and risk of fracture did not achieve statistical significance (SNP rs12521522[A]: OR=1.14 (95% CI 0.92 to 1.35), p=0.14). Meta-analysis also revealed two genome-wide suggestive loci at CTNNA2 and 6q23.2. CONCLUSIONS: These findings demonstrate that variants at MEF2C were associated with forearm BMD, implicating this gene in the determination of BMD at forearm.

Genome-wide association analysis identifies variants associated with nonalcoholic fatty liver disease that have distinct effects on metabolic traits.

Nonalcoholic fatty liver disease (NAFLD) clusters in families, but the only known common genetic variants influencing risk are near PNPLA3. We sought to identify additional genetic variants influencing NAFLD using genome-wide association (GWA) analysis of computed tomography (CT) measured hepatic steatosis, a non-invasive measure of NAFLD, in large population based samples. Using variance components methods, we show that CT hepatic steatosis is heritable (∼26%-27%) in family-based Amish, Family Heart, and Framingham Heart Studies (n = 880 to 3,070). By carrying out a fixed-effects meta-analysis of genome-wide association (GWA) results between CT hepatic steatosis and ∼2.4 million imputed or genotyped SNPs in 7,176 individuals from the Old Order Amish, Age, Gene/Environment Susceptibility-Reykjavik study (AGES), Family Heart, and Framingham Heart Studies, we identify variants associated at genome-wide significant levels (p<5×10(-8)) in or near PNPLA3, NCAN, and PPP1R3B. We genotype these and 42 other top CT hepatic steatosis-associated SNPs in 592 subjects with biopsy-proven NAFLD from the NASH Clinical Research Network (NASH CRN). In comparisons with 1,405 healthy controls from the Myocardial Genetics Consortium (MIGen), we observe significant associations with histologic NAFLD at variants in or near NCAN, GCKR, LYPLAL1, and PNPLA3, but not PPP1R3B. Variants at these five loci exhibit distinct patterns of association with serum lipids, as well as glycemic and anthropometric traits. We identify common genetic variants influencing CT-assessed steatosis and risk of NAFLD. Hepatic steatosis associated variants are not uniformly associated with NASH/fibrosis or result in abnormalities in serum lipids or glycemic and anthropometric traits, suggesting genetic heterogeneity in the pathways influencing these traits.

Genetic variation at NCAN locus is associated with inflammation and fibrosis in non-alcoholic fatty liver disease in morbid obesity.

OBJECTIVE: Obesity-associated non-alcoholic fatty liver disease (NAFLD) may cause liver dysfunction and failure. In a previously reported genome-wide association meta-analysis, single nucleotide polymorphisms (SNPs) near PNPLA3, NCAN, GCKR, LYPLAL1 and PPP1R3B were associated with NAFLD and with distinctive serum lipid profiles. The present study examined the relevance of these variants to NAFLD in extreme obesity. METHODS: In 1,092 bariatric surgery patients, the candidate SNPs were genotyped and association analyses with liver histology and serum lipids were performed. RESULTS: We replicated the association of hepatosteatosis with PNPLA3 rs738409[G] and with NCAN rs2228603[T]. We also replicated the association of rs2228603[T] with hepatic inflammation and fibrosis. rs2228603[T] was associated with lower serum low-density lipoprotein, total cholesterol and triglycerides. After stratification by the presence or absence of NAFLD, these associations were present predominantly in the subgroup with NAFLD. CONCLUSION: NCAN rs2228603[T] is a risk factor for liver inflammation and fibrosis, suggesting that this locus is responsible for progression from steatosis to steatohepatitis. In this bariatric cohort, rs2228603[T] was associated with low serum lipids only in patients with NAFLD. This supports a NAFLD model in which the liver may sequester triglycerides as a result of either increased triglyceride uptake and/or decreased lipolysis.

Causal relationship between obesity and vitamin D status: bi-directional Mendelian randomization analysis of multiple cohorts.

BACKGROUND: Obesity is associated with vitamin D deficiency, and both are areas of active public health concern. We explored the causality and direction of the relationship between body mass index (BMI) and 25-hydroxyvitamin D [25(OH)D] using genetic markers as instrumental variables (IVs) in bi-directional Mendelian randomization (MR) analysis. METHODS AND FINDINGS: We used information from 21 adult cohorts (up to 42,024 participants) with 12 BMI-related SNPs (combined in an allelic score) to produce an instrument for BMI and four SNPs associated with 25(OH)D (combined in two allelic scores, separately for genes encoding its synthesis or metabolism) as an instrument for vitamin D. Regression estimates for the IVs (allele scores) were generated within-study and pooled by meta-analysis to generate summary effects. Associations between vitamin D scores and BMI were confirmed in the Genetic Investigation of Anthropometric Traits (GIANT) consortium (n = 123,864). Each 1 kg/m(2) higher BMI was associated with 1.15% lower 25(OH)D (p = 6.52×10⁻²7). The BMI allele score was associated both with BMI (p = 6.30×10⁻6²) and 25(OH)D (-0.06% [95% CI -0.10 to -0.02], p = 0.004) in the cohorts that underwent meta-analysis. The two vitamin D allele scores were strongly associated with 25(OH)D (p≤8.07×10⁻57 for both scores) but not with BMI (synthesis score, p = 0.88; metabolism score, p = 0.08) in the meta-analysis. A 10% higher genetically instrumented BMI was associated with 4.2% lower 25(OH)D concentrations (IV ratio: -4.2 [95% CI -7.1 to -1.3], p = 0.005). No association was seen for genetically instrumented 25(OH)D with BMI, a finding that was confirmed using data from the GIANT consortium (p≥0.57 for both vitamin D scores). CONCLUSIONS: On the basis of a bi-directional genetic approach that limits confounding, our study suggests that a higher BMI leads to lower 25(OH)D, while any effects of lower 25(OH)D increasing BMI are likely to be small. Population level interventions to reduce BMI are expected to decrease the prevalence of vitamin D deficiency.

Association of adiposity genetic variants with menarche timing in 92,105 women of European descent.

Obesity is of global health concern. There are well-described inverse relationships between female pubertal timing and obesity. Recent genome-wide association studies of age at menarche identified several obesity-related variants. Using data from the ReproGen Consortium, we employed meta-analytical techniques to estimate the associations of 95 a priori and recently identified obesity-related (body mass index (weight (kg)/height (m)(2)), waist circumference, and waist:hip ratio) single-nucleotide polymorphisms (SNPs) with age at menarche in 92,116 women of European descent from 38 studies (1970-2010), in order to estimate associations between genetic variants associated with central or overall adiposity and pubertal timing in girls. Investigators in each study performed a separate analysis of associations between the selected SNPs and age at menarche (ages 9-17 years) using linear regression models and adjusting for birth year, site (as appropriate), and population stratification. Heterogeneity of effect-measure estimates was investigated using meta-regression. Six novel associations of body mass index loci with age at menarche were identified, and 11 adiposity loci previously reported to be associated with age at menarche were confirmed, but none of the central adiposity variants individually showed significant associations. These findings suggest complex genetic relationships between menarche and overall obesity, and to a lesser extent central obesity, in normal processes of growth and development.

Association of the vitamin D metabolism gene CYP24A1 with coronary artery calcification.

OBJECTIVE: The vitamin D endocrine system is essential for calcium homeostasis, and low levels of vitamin D metabolites have been associated with cardiovascular disease risk. We hypothesized that DNA sequence variation in genes regulating vitamin D metabolism and signaling pathways might influence variation in coronary artery calcification (CAC). METHODS AND RESULTS: We genotyped single-nucleotide polymorphisms (SNPs) in GC, CYP27B1, CYP24A1, and VDR and tested their association with CAC quantity, as measured by electron beam computed tomography. Initial association studies were carried out in a discovery sample comprising 697 Amish subjects, and SNPs nominally associated with CAC quantity (4 SNPs in CYP24A1, P=0.008 to 0.00003) were then tested for association with CAC quantity in 2 independent cohorts of subjects of white European ancestry (Genetic Epidemiology Network of Arteriopathy study [n=916] and the Penn Coronary Artery Calcification sample [n=2061]). One of the 4 SNPs, rs2762939, was associated with CAC quantity in both the Genetic Epidemiology Network of Arteriopathy (P=0.007) and Penn Coronary Artery Calcification (P=0.01) studies. In all 3 populations, the rs2762939 C allele was associated with lower CAC quantity. Metaanalysis for the association of this SNP with CAC quantity across all 3 studies yielded a P value of 2.9×10(-6). CONCLUSIONS: A common SNP in the CYP24A1 gene was associated with CAC quantity in 3 independent populations. This result suggests a role for vitamin D metabolism in the development of CAC quantity.

Pleiotropic associations of heterozygosity for the SERPINA1 Z allele in the UK Biobank.

Homozygosity for the SERPINA1 Z allele causes α1-antitrypsin deficiency, a rare condition that can cause lung and liver disease. However, the effects of Z allele heterozygosity on nonrespiratory phenotypes, and on lung function in the general population, remain unclear. We conducted a large, population-based study to determine Z allele effects on >2400 phenotypes in the UK Biobank (N=303 353). Z allele heterozygosity was strongly associated with increased height (ß=1.02 cm, p=3.91×10-68), and with other nonrespiratory phenotypes including increased risk of gall bladder disease, reduced risk of heart disease and lower blood pressure, reduced risk of osteoarthritis and reduced bone mineral density, increased risk of headache and enlarged prostate, as well as with blood biomarkers of liver function. Heterozygosity was associated with higher height-adjusted forced expiratory volume in 1 s (FEV1) (ß=19.36 mL, p=9.21×10-4) and FEV1/forced vital capacity (ß=0.0031, p=1.22×10-5) in nonsmokers, whereas in smokers, this protective effect was abolished. Furthermore, we show for the first time that sex modifies the association of the Z allele on lung function. We conclude that Z allele heterozygosity and homozygosity exhibit opposing effects on lung function in the UK population, and that these associations are modified by smoking and sex. In exploratory analyses, heterozygosity for the Z allele also showed pleiotropic associations with nonrespiratory health-related traits and disease risk.

A Noncoding Variant Near PPP1R3B Promotes Liver Glycogen Storage and MetS, but Protects Against Myocardial Infarction.

CONTEXT: Glycogen storage diseases are rare. Increased glycogen in the liver results in increased attenuation. OBJECTIVE: Investigate the association and function of a noncoding region associated with liver attenuation but not histologic nonalcoholic fatty liver disease. DESIGN: Genetics of Obesity-associated Liver Disease Consortium. SETTING: Population-based. MAIN OUTCOME: Computed tomography measured liver attenuation. RESULTS: Carriers of rs4841132-A (frequency 2%-19%) do not show increased hepatic steatosis; they have increased liver attenuation indicative of increased glycogen deposition. rs4841132 falls in a noncoding RNA LOC157273 ~190 kb upstream of PPP1R3B. We demonstrate that rs4841132-A increases PPP1R3B through a cis genetic effect. Using CRISPR/Cas9 we engineered a 105-bp deletion including rs4841132-A in human hepatocarcinoma cells that increases PPP1R3B, decreases LOC157273, and increases glycogen perfectly mirroring the human disease. Overexpression of PPP1R3B or knockdown of LOC157273 increased glycogen but did not result in decreased LOC157273 or increased PPP1R3B, respectively, suggesting that the effects may not all occur via affecting RNA levels. Based on electronic health record (EHR) data, rs4841132-A associates with all components of the metabolic syndrome (MetS). However, rs4841132-A associated with decreased low-density lipoprotein (LDL) cholesterol and risk for myocardial infarction (MI). A metabolic signature for rs4841132-A includes increased glycine, lactate, triglycerides, and decreased acetoacetate and beta-hydroxybutyrate. CONCLUSIONS: These results show that rs4841132-A promotes a hepatic glycogen storage disease by increasing PPP1R3B and decreasing LOC157273. rs4841132-A promotes glycogen accumulation and development of MetS but lowers LDL cholesterol and risk for MI. These results suggest that elevated hepatic glycogen is one cause of MetS that does not invariably promote MI.