Food allergy (FA) arises from a complex interplay between an individual's genetic predisposition and environmental factors, and its prevalence is increasing. Genome-wide association studies to date ha Show more
Food allergy (FA) arises from a complex interplay between an individual's genetic predisposition and environmental factors, and its prevalence is increasing. Genome-wide association studies to date have been hindered by small sample sizes and varying FA definitions. We sought to identify novel FA risk loci by conducting a genome-wide association study meta-analysis in children and adults by using a multiphenotype approach to ensure a good trade-off between sufficient sample size and valid FA definitions. Analyses were conducted separately in children and adults on the basis of the following FA phenotypes: self-report, doctor diagnosis, food-specific sensitization, and doctor diagnosis plus food-specific sensitization. A meta-analysis was performed of genome-wide association studies from up to 16 cohorts of people of European ancestry including 229,426 adults and 14,234 children. Models were adjusted for sex, age, principal components, and, if applicable, further study-specific confounders. Sensitivity models were additionally adjusted for hay fever. Replication was conducted in additional external cohorts and a validation in oral food challenge-defined FA cases. Thirty-seven single nucleotide polymorphisms met suggestive significance (P < 1 × 10 This study identified 37 single nucleotide polymorphisms suggestively associated with FA and demonstrated genetic differences across phenotypes. It highlights the need for a unified FA definition and sheds light on FA's shared genetic architecture with allergies. Show less
Identifying genetic determinants of reproductive success may highlight mechanisms underlying fertility and identify alleles under present-day selection. Using data in 785,604 individuals of European a Show more
Identifying genetic determinants of reproductive success may highlight mechanisms underlying fertility and identify alleles under present-day selection. Using data in 785,604 individuals of European ancestry, we identified 43 genomic loci associated with either number of children ever born (NEB) or childlessness. These loci span diverse aspects of reproductive biology, including puberty timing, age at first birth, sex hormone regulation, endometriosis and age at menopause. Missense variants in ARHGAP27 were associated with higher NEB but shorter reproductive lifespan, suggesting a trade-off at this locus between reproductive ageing and intensity. Other genes implicated by coding variants include PIK3IP1, ZFP82 and LRP4, and our results suggest a new role for the melanocortin 1 receptor (MC1R) in reproductive biology. As NEB is one component of evolutionary fitness, our identified associations indicate loci under present-day natural selection. Integration with data from historical selection scans highlighted an allele in the FADS1/2 gene locus that has been under selection for thousands of years and remains so today. Collectively, our findings demonstrate that a broad range of biological mechanisms contribute to reproductive success. Show less
Glycosaminoglycans are ubiquitously expressed polysaccharides that are attached to proteoglycans. Here, we showed that ablation of the heparan sulfate (HS) polymerase Ext1 in retinal progenitor cells Show more
Glycosaminoglycans are ubiquitously expressed polysaccharides that are attached to proteoglycans. Here, we showed that ablation of the heparan sulfate (HS) polymerase Ext1 in retinal progenitor cells did not affect initial progression of retinal angiogenesis, but it disrupted the pruning of blood vessels and establishment of arterioles and venules. In the absence of retinal HS, blood vessels were also vulnerable to high oxygen tension in early postnatal stages, which could be rescued by exogenous vascular endothelial growth factor (VEGF), consistent with the role of retinal HS in the fine-tuning of VEGF signaling. Furthermore, we observed that the retinal inner limiting membrane (ILM) was disrupted by deletion of Ext1 in a timing-specific manner, suggesting that retinal HS is required for the assembly but not the maintenance of the basement membrane. Lastly, we showed that further deletion of C4st1, a chondroitin sulfate (CS) sulfation enzyme, did not affect the assembly of the ILM but, when combined with Ext1 deletion, it aggravated the retinal permeability by disrupting the retinal glycocalyx. These results demonstrate an important role of CS and HS in establishing the barrier function of the extracellular matrix. Show less
The HERMES (HEart failure Molecular Epidemiology for Therapeutic targetS) consortium aims to identify the genomic and molecular basis of heart failure. The consortium currently includes 51 studies fro Show more
The HERMES (HEart failure Molecular Epidemiology for Therapeutic targetS) consortium aims to identify the genomic and molecular basis of heart failure. The consortium currently includes 51 studies from 11 countries, including 68 157 heart failure cases and 949 888 controls, with data on heart failure events and prognosis. All studies collected biological samples and performed genome-wide genotyping of common genetic variants. The enrolment of subjects into participating studies ranged from 1948 to the present day, and the median follow-up following heart failure diagnosis ranged from 2 to 116 months. Forty-nine of 51 individual studies enrolled participants of both sexes; in these studies, participants with heart failure were predominantly male (34-90%). The mean age at diagnosis or ascertainment across all studies ranged from 54 to 84 years. Based on the aggregate sample, we estimated 80% power to genetic variant associations with risk of heart failure with an odds ratio of ≥1.10 for common variants (allele frequency ≥ 0.05) and ≥1.20 for low-frequency variants (allele frequency 0.01-0.05) at P < 5 × 10 HERMES is a global collaboration aiming to (i) identify the genetic determinants of heart failure; (ii) generate insights into the causal pathways leading to heart failure and enable genetic approaches to target prioritization; and (iii) develop genomic tools for disease stratification and risk prediction. Show less
Osteoarthritis (OA) is a progressive degenerative disease of the articular cartilage caused by an unbalanced activity of proteases, cytokines and other secreted proteins. Since heparan sulfate (HS) de Show more
Osteoarthritis (OA) is a progressive degenerative disease of the articular cartilage caused by an unbalanced activity of proteases, cytokines and other secreted proteins. Since heparan sulfate (HS) determines the activity of many extracellular factors, we investigated its role in OA progression. To analyze the role of the HS level, OA was induced by anterior cruciate ligament transection (ACLT) in transgenic mice carrying a loss-of-function allele of Ext1 in clones of chondrocytes (Col2-rtTA-Cre;Ext1 All investigated mouse strains showed reduced OA scores (Col2-rtTA-Cre;Ext1 A decreased HS content or a reduced sulfation level protect against OA progression by regulating protease activity rather than expression. Show less
Heparan sulfate (HS) regulates the activity of many signaling molecules critical for the development of endochondral bones. Even so, mice with a genetically altered HS metabolism display a relatively Show more
Heparan sulfate (HS) regulates the activity of many signaling molecules critical for the development of endochondral bones. Even so, mice with a genetically altered HS metabolism display a relatively mild skeletal phenotype compared to the defects observed in other tissues and organs pointing to a reduced HS dependency of growth-factor signaling in chondrocytes. To understand this difference, we have investigated the glycosaminoglycan (GAG) composition in two mouse lines that produce either reduced levels of HS (Ext1 Show less
Cofactors are essential for driving recombinant prion protein into pathogenic conformers. Polyanions promote prion aggregation in vitro, yet the cofactors that modulate prion assembly in vivo remain l Show more
Cofactors are essential for driving recombinant prion protein into pathogenic conformers. Polyanions promote prion aggregation in vitro, yet the cofactors that modulate prion assembly in vivo remain largely unknown. Here we report that the endogenous glycosaminoglycan, heparan sulfate (HS), impacts prion propagation kinetics and deposition sites in the brain. Exostosin-1 haploinsufficient (Ext1 Show less
Binding and uptake of triglyceride-rich lipoproteins (TRLs) in mice depend on heparan sulfate and the hepatic proteoglycan, syndecan-1 (SDC1). Alteration of glucosamine N-sulfation by deletion of gluc Show more
Binding and uptake of triglyceride-rich lipoproteins (TRLs) in mice depend on heparan sulfate and the hepatic proteoglycan, syndecan-1 (SDC1). Alteration of glucosamine N-sulfation by deletion of glucosamine N-deacetylase-N-sulfotransferase 1 (Ndst1) and 2-O-sulfation of uronic acids by deletion of uronyl 2-O-sulfotransferase (Hs2st) led to diminished lipoprotein metabolism, whereas inactivation of glucosaminyl 6-O-sulfotransferase 1 (Hs6st1), which encodes one of the three 6-O-sulfotransferases, had little effect on lipoprotein binding. However, other studies have suggested that 6-O-sulfation may be important for TRL binding and uptake. In order to explain these discrepant findings, we used CRISPR/Cas9 gene editing to create a library of mutants in the human hepatoma cell line, Hep3B. Inactivation of EXT1 encoding the heparan sulfate copolymerase, NDST1 and HS2ST dramatically reduced binding of TRLs. Inactivation of HS6ST1 had no effect, but deletion of HS6ST2 reduced TRL binding. Compounding mutations in HS6ST1 and HS6ST2 did not exacerbate this effect indicating that HS6ST2 is the dominant 6-O-sulfotransferase and that binding of TRLs indeed depends on 6-O-sulfation of glucosamine residues. Uptake studies showed that TRL internalization was also affected in 6-O-sulfation deficient cells. Interestingly, genetic deletion of SDC1 only marginally impacted binding of TRLs but reduced TRL uptake to the same extent as treating the cells with heparin lyases. These findings confirm that SDC1 is the dominant endocytic proteoglycan receptor for TRLs in human Hep3B cells and that binding and uptake of TRLs depend on SDC1 and N- and 2-O-sulfation as well as 6-O-sulfation of heparan sulfate chains catalyzed by HS6ST2. Show less
Genome-wide association studies (GWAS) have identified >250 loci for body mass index (BMI), implicating pathways related to neuronal biology. Most GWAS loci represent clusters of common, noncoding var Show more
Genome-wide association studies (GWAS) have identified >250 loci for body mass index (BMI), implicating pathways related to neuronal biology. Most GWAS loci represent clusters of common, noncoding variants from which pinpointing causal genes remains challenging. Here we combined data from 718,734 individuals to discover rare and low-frequency (minor allele frequency (MAF) < 5%) coding variants associated with BMI. We identified 14 coding variants in 13 genes, of which 8 variants were in genes (ZBTB7B, ACHE, RAPGEF3, RAB21, ZFHX3, ENTPD6, ZFR2 and ZNF169) newly implicated in human obesity, 2 variants were in genes (MC4R and KSR2) previously observed to be mutated in extreme obesity and 2 variants were in GIPR. The effect sizes of rare variants are ~10 times larger than those of common variants, with the largest effect observed in carriers of an MC4R mutation introducing a stop codon (p.Tyr35Ter, MAF = 0.01%), who weighed ~7 kg more than non-carriers. Pathway analyses based on the variants associated with BMI confirm enrichment of neuronal genes and provide new evidence for adipocyte and energy expenditure biology, widening the potential of genetically supported therapeutic targets in obesity. Show less
Migraine and major depressive disorder (MDD) are common brain disorders that frequently co-occur. Despite epidemiological evidence that migraine and MDD share a genetic basis, their overlap at the mol Show more
Migraine and major depressive disorder (MDD) are common brain disorders that frequently co-occur. Despite epidemiological evidence that migraine and MDD share a genetic basis, their overlap at the molecular genetic level has not been thoroughly investigated. Using single-nucleotide polymorphism (SNP) and gene-based analysis of genome-wide association study (GWAS) genotype data, we found significant genetic overlap across the two disorders. LD Score regression revealed a significant SNP-based heritability for both migraine (h Show less
Multiple hereditary exostoses (MHE) is an autosomal dominant disorder that affects about 1 in 50,000 children worldwide. MHE, also known as hereditary multiple exostoses (HME) or multiple osteochondro Show more
Multiple hereditary exostoses (MHE) is an autosomal dominant disorder that affects about 1 in 50,000 children worldwide. MHE, also known as hereditary multiple exostoses (HME) or multiple osteochondromas (MO), is characterized by cartilage-capped outgrowths called osteochondromas that develop adjacent to the growth plates of skeletal elements in young patients. These benign tumors can affect growth plate function, leading to skeletal growth retardation, or deformations, and can encroach on nerves, tendons, muscles, and other surrounding tissues and cause motion impairment, chronic pain, and early onset osteoarthritis. In about 2-5% of patients, the osteochondromas can become malignant and life threatening. Current treatments consist of surgical removal of the most symptomatic tumors and correction of the major skeletal defects, but physical difficulties and chronic pain usually continue and patients may undergo multiple surgeries throughout life. Thus, there is an urgent need to find new treatments to prevent or reverse osteochondroma formation. The 2016 International MHE Research Conference was convened to provide a forum for the presentation of the most up-to-date and advanced clinical and basic science data and insights in MHE and related fields; to stimulate the forging of new perspectives, collaborations, and venues of research; and to publicize key scientific findings within the biomedical research community and share insights and relevant information with MHE patients and their families. This report provides a description, review, and assessment of all the exciting and promising studies presented at the Conference and delineates a general roadmap for future MHE research targets and goals. Show less
We screened variants on an exome-focused genotyping array in >300,000 participants (replication in >280,000 participants) and identified 444 independent variants in 250 loci significantly associated w Show more
We screened variants on an exome-focused genotyping array in >300,000 participants (replication in >280,000 participants) and identified 444 independent variants in 250 loci significantly associated with total cholesterol (TC), high-density-lipoprotein cholesterol (HDL-C), low-density-lipoprotein cholesterol (LDL-C), and/or triglycerides (TG). At two loci (JAK2 and A1CF), experimental analysis in mice showed lipid changes consistent with the human data. We also found that: (i) beta-thalassemia trait carriers displayed lower TC and were protected from coronary artery disease (CAD); (ii) excluding the CETP locus, there was not a predictable relationship between plasma HDL-C and risk for age-related macular degeneration; (iii) only some mechanisms of lowering LDL-C appeared to increase risk for type 2 diabetes (T2D); and (iv) TG-lowering alleles involved in hepatic production of TG-rich lipoproteins (TM6SF2 and PNPLA3) tracked with higher liver fat, higher risk for T2D, and lower risk for CAD, whereas TG-lowering alleles involved in peripheral lipolysis (LPL and ANGPTL4) had no effect on liver fat but decreased risks for both T2D and CAD. Show less
Genome-wide association studies have so far identified 56 loci associated with risk of coronary artery disease (CAD). Many CAD loci show pleiotropy; that is, they are also associated with other diseas Show more
Genome-wide association studies have so far identified 56 loci associated with risk of coronary artery disease (CAD). Many CAD loci show pleiotropy; that is, they are also associated with other diseases or traits. This study sought to systematically test if genetic variants identified for non-CAD diseases/traits also associate with CAD and to undertake a comprehensive analysis of the extent of pleiotropy of all CAD loci. In discovery analyses involving 42,335 CAD cases and 78,240 control subjects we tested the association of 29,383 common (minor allele frequency >5%) single nucleotide polymorphisms available on the exome array, which included a substantial proportion of known or suspected single nucleotide polymorphisms associated with common diseases or traits as of 2011. Suggestive association signals were replicated in an additional 30,533 cases and 42,530 control subjects. To evaluate pleiotropy, we tested CAD loci for association with cardiovascular risk factors (lipid traits, blood pressure phenotypes, body mass index, diabetes, and smoking behavior), as well as with other diseases/traits through interrogation of currently available genome-wide association study catalogs. We identified 6 new loci associated with CAD at genome-wide significance: on 2q37 (KCNJ13-GIGYF2), 6p21 (C2), 11p15 (MRVI1-CTR9), 12q13 (LRP1), 12q24 (SCARB1), and 16q13 (CETP). Risk allele frequencies ranged from 0.15 to 0.86, and odds ratio per copy of the risk allele ranged from 1.04 to 1.09. Of 62 new and known CAD loci, 24 (38.7%) showed statistical association with a traditional cardiovascular risk factor, with some showing multiple associations, and 29 (47%) showed associations at p < 1 × 10 We identified 6 loci associated with CAD at genome-wide significance. Several CAD loci show substantial pleiotropy, which may help us understand the mechanisms by which these loci affect CAD risk. Show less
Circulating blood cell counts and indices are important indicators of hematopoietic function and a number of clinical parameters, such as blood oxygen-carrying capacity, inflammation, and hemostasis. Show more
Circulating blood cell counts and indices are important indicators of hematopoietic function and a number of clinical parameters, such as blood oxygen-carrying capacity, inflammation, and hemostasis. By performing whole-exome sequence association analyses of hematologic quantitative traits in 15,459 community-dwelling individuals, followed by in silico replication in up to 52,024 independent samples, we identified two previously undescribed coding variants associated with lower platelet count: a common missense variant in CPS1 (rs1047891, MAF = 0.33, discovery + replication p = 6.38 × 10(-10)) and a rare synonymous variant in GFI1B (rs150813342, MAF = 0.009, discovery + replication p = 1.79 × 10(-27)). By performing CRISPR/Cas9 genome editing in hematopoietic cell lines and follow-up targeted knockdown experiments in primary human hematopoietic stem and progenitor cells, we demonstrate an alternative splicing mechanism by which the GFI1B rs150813342 variant suppresses formation of a GFI1B isoform that preferentially promotes megakaryocyte differentiation and platelet production. These results demonstrate how unbiased studies of natural variation in blood cell traits can provide insight into the regulation of human hematopoiesis. Show less
Adeno-associated virus 2 (AAV2) and adenovirus 5 (Ad5) are promising gene therapy vectors. Both display liver tropism and are currently thought to enter hepatocytes in vivo through cell surface hepara Show more
Adeno-associated virus 2 (AAV2) and adenovirus 5 (Ad5) are promising gene therapy vectors. Both display liver tropism and are currently thought to enter hepatocytes in vivo through cell surface heparan sulfate proteoglycans (HSPGs). To test directly this hypothesis, we created mice that lack Ext1, an enzyme required for heparan sulfate biosynthesis, in hepatocytes. Ext1(HEP) mutant mice exhibit an 8-fold reduction of heparan sulfate in primary hepatocytes and a 5-fold reduction of heparan sulfate in whole liver tissue. Conditional hepatocyte Ext1 gene deletion greatly reduced AAV2 liver transduction following intravenous injection. Ad5 transduction requires blood coagulation factor X (FX); FX binds to the Ad5 capsid hexon protein and bridges the virus to HSPGs on the cell surface. Ad5.FX transduction was abrogated in primary hepatocytes from Ext1(HEP) mice. However, in contrast to the case with AAV2, Ad5 transduction was not significantly reduced in the livers of Ext1(HEP) mice. FX remained essential for Ad5 transduction in vivo in Ext1(HEP) mice. We conclude that while AAV2 requires HSPGs for entry into mouse hepatocytes, HSPGs are dispensable for Ad5 hepatocyte transduction in vivo. This study reopens the question of how adenovirus enters cells in vivo. Our understanding of how viruses enter cells, and how they can be used as therapeutic vectors to manage disease, begins with identification of the cell surface receptors to which viruses bind and which mediate viral entry. Both adeno-associated virus 2 and adenovirus 5 are currently thought to enter hepatocytes in vivo through heparan sulfate proteoglycans (HSPGs). However, direct evidence for these conclusions is lacking. Experiments presented herein, in which hepatic heparan sulfate synthesis was genetically abolished, demonstrated that HSPGs are not likely to function as hepatocyte Ad5 receptors in vivo. The data also demonstrate that HSPGs are required for hepatocyte transduction by AAV2. These results reopen the question of the identity of the Ad5 receptor in vivo and emphasize the necessity of demonstrating the nature of the receptor by genetic means, both for understanding Ad5 entry into cells in vivo and for optimization of Ad5 vectors as therapeutic agents. Show less
White blood cells play diverse roles in innate and adaptive immunity. Genetic association analyses of phenotypic variation in circulating white blood cell (WBC) counts from large samples of otherwise Show more
White blood cells play diverse roles in innate and adaptive immunity. Genetic association analyses of phenotypic variation in circulating white blood cell (WBC) counts from large samples of otherwise healthy individuals can provide insights into genes and biologic pathways involved in production, differentiation, or clearance of particular WBC lineages (myeloid, lymphoid) and also potentially inform the genetic basis of autoimmune, allergic, and blood diseases. We performed an exome array-based meta-analysis of total WBC and subtype counts (neutrophils, monocytes, lymphocytes, basophils, and eosinophils) in a multi-ancestry discovery and replication sample of ∼157,622 individuals from 25 studies. We identified 16 common variants (8 of which were coding variants) associated with one or more WBC traits, the majority of which are pleiotropically associated with autoimmune diseases. Based on functional annotation, these loci included genes encoding surface markers of myeloid, lymphoid, or hematopoietic stem cell differentiation (CD69, CD33, CD87), transcription factors regulating lineage specification during hematopoiesis (ASXL1, IRF8, IKZF1, JMJD1C, ETS2-PSMG1), and molecules involved in neutrophil clearance/apoptosis (C10orf54, LTA), adhesion (TNXB), or centrosome and microtubule structure/function (KIF9, TUBD1). Together with recent reports of somatic ASXL1 mutations among individuals with idiopathic cytopenias or clonal hematopoiesis of undetermined significance, the identification of a common regulatory 3' UTR variant of ASXL1 suggests that both germline and somatic ASXL1 mutations contribute to lower blood counts in otherwise asymptomatic individuals. These association results shed light on genetic mechanisms that regulate circulating WBC counts and suggest a prominent shared genetic architecture with inflammatory and autoimmune diseases. Show less
Heparan sulfate (HS) is a highly modified glycosaminoglycan (GAG) bound to a core protein to form heparan sulfate proteoglycans (HSPGs) that are vital in many cellular processes ranging from developme Show more
Heparan sulfate (HS) is a highly modified glycosaminoglycan (GAG) bound to a core protein to form heparan sulfate proteoglycans (HSPGs) that are vital in many cellular processes ranging from development to adult physiology, as well as in disease, through interactions with various protein ligands. This study aimed to elucidate the role of HS in corneal epithelial homeostasis and wound healing. An inducible quadruple transgenic mouse model was generated to excise Ext1 and Ndst1, which encode the critical HS chain elongation enzyme and N-deacetylase/N-sulfotransferase, respectively, in keratin 14-positive cells upon doxycycline induction. EXT(Δ/ΔCEpi) mice (deletion of Ext1 in corneal epithelium) induced at P20 presented progressive thinning of the corneal epithelium with a significant loss in the number of epithelial layers by P55. EXT(Δ/ΔCEpi) mice presented tight junction disruption, loss of cell-basement membrane adhesion complexes, and impaired wound healing. Interestingly, EXT(Δ/ΔCEpi) and NDST(Δ/ΔCEpi) mice presented an increase in cell proliferation, which was assayed by both Ki67 staining and 5-ethynyl-2'-deoxyuridine (EdU) incorporation. Moreover, EXT(Δ/ΔCEpi) mice presented compromised epithelial stratification 7 days after a debridement wound. The conditional knockout of HS from keratocytes using the keratocan promoter led to no corneal abnormalities or any disruption in wound healing. Corneal epithelial cells require HS for maintaining corneal homeostasis, and the loss of epithelial HS leads to both impaired wound healing and impaired corneal stratification. Show less
Elevated nonfasting TG-rich lipoprotein levels are a risk factor for CVD. To further evaluate the relevance of LDL-receptor (LDLr) pathway and heparan sulfate proteoglycans (HSPGs) in TG homeostasis, Show more
Elevated nonfasting TG-rich lipoprotein levels are a risk factor for CVD. To further evaluate the relevance of LDL-receptor (LDLr) pathway and heparan sulfate proteoglycans (HSPGs) in TG homeostasis, we analyzed fasting and postprandial TG levels in mice bearing combined heterozygous mutations in both Exostosin (Ext) 1 and Ldlr, in subjects with hereditary multiple exostosis (HME) due to a heterozygous loss-of-function mutation in EXT1 or EXT2 (N = 13), and in patients with heterozygous mutations in LDLR [familial hypercholesterolemia (FH)] and SNPs in major HSPG-related genes (n = 22). Mice bearing a homozygous mutation in hepatic Ext1 exhibited elevated plasma TGs similar to mice lacking other key enzymes involved in HSPG assembly. Compound heterozygous mice lacking Ldlr and Ext1 showed synergy on plasma TG accumulation and postprandial clearance. In human subjects, a trend was observed in HME patients toward reduced postprandial TG clearance with a concomitant reduction in chylomicron clearance [area under the curve (AUC)-retinyl ester (RE) HME, 844 ± 127 vs. controls, 646 ± 119 nM/h, P = 0.09]. Moreover, in FH subjects with a high HSPG gene score, retinyl palmitate excursions were higher (AUC-RE, 2,377 ± 293 vs. 1,565 ± 181 nM/h, P < 0.05). Incremental AUC-apoB48 was similar between the groups. In conclusion, the data are supportive for a minor yet additive role of HSPG in human postprandial TG clearance, and further studies are warranted. Show less
Coffee, a major dietary source of caffeine, is among the most widely consumed beverages in the world and has received considerable attention regarding health risks and benefits. We conducted a genome- Show more
Coffee, a major dietary source of caffeine, is among the most widely consumed beverages in the world and has received considerable attention regarding health risks and benefits. We conducted a genome-wide (GW) meta-analysis of predominately regular-type coffee consumption (cups per day) among up to 91,462 coffee consumers of European ancestry with top single-nucleotide polymorphisms (SNPs) followed-up in ~30 062 and 7964 coffee consumers of European and African-American ancestry, respectively. Studies from both stages were combined in a trans-ethnic meta-analysis. Confirmed loci were examined for putative functional and biological relevance. Eight loci, including six novel loci, met GW significance (log10Bayes factor (BF)>5.64) with per-allele effect sizes of 0.03-0.14 cups per day. Six are located in or near genes potentially involved in pharmacokinetics (ABCG2, AHR, POR and CYP1A2) and pharmacodynamics (BDNF and SLC6A4) of caffeine. Two map to GCKR and MLXIPL genes related to metabolic traits but lacking known roles in coffee consumption. Enhancer and promoter histone marks populate the regions of many confirmed loci and several potential regulatory SNPs are highly correlated with the lead SNP of each. SNP alleles near GCKR, MLXIPL, BDNF and CYP1A2 that were associated with higher coffee consumption have previously been associated with smoking initiation, higher adiposity and fasting insulin and glucose but lower blood pressure and favorable lipid, inflammatory and liver enzyme profiles (P<5 × 10(-8)).Our genetic findings among European and African-American adults reinforce the role of caffeine in mediating habitual coffee consumption and may point to molecular mechanisms underlying inter-individual variability in pharmacological and health effects of coffee. Show less
Exotosin (EXT) proteins are involved in the chain elongation step of heparan sulfate (HS) biosynthesis, which is intricately involved in organ development. Loss of function mutations (LOF) in EXT1 and Show more
Exotosin (EXT) proteins are involved in the chain elongation step of heparan sulfate (HS) biosynthesis, which is intricately involved in organ development. Loss of function mutations (LOF) in EXT1 and EXT2 result in hereditary exostoses (HME). Interestingly, HS plays a role in pancreas development and beta-cell function, and genetic variations in EXT2 are associated with an increased risk for type 2 diabetes mellitus. We hypothesized that loss of function of EXT1 or EXT2 in subjects with hereditary multiple exostoses (HME) affects pancreatic insulin secretion capacity and development. We performed an oral glucose tolerance test (OGTT) followed by hyperglycemic clamps to investigate first-phase glucose-stimulated insulin secretion (GSIS) in HME patients and age and gender matched non-affected relatives. Pancreas volume was assessed with magnetic resonance imaging (MRI). OGTT did not reveal significant differences in glucose disposal, but there was a markedly lower GSIS in HME subjects during hyperglycemic clamp (iAUC HME: 0.72 [0.46-1.16] vs. controls 1.53 [0.69-3.36] nmol·l-1·min-1, p<0.05). Maximal insulin response following arginine challenge was also significantly attenuated (iAUC HME: 7.14 [4.22-10.5] vs. controls 10.2 [7.91-12.70] nmol·l-1·min-1 p<0.05), indicative of an impaired beta-cell reserve. MRI revealed a significantly smaller pancreatic volume in HME subjects (HME: 72.0±15.8 vs. controls 96.5±26.0 cm3 p = 0.04). In conclusion, loss of function of EXT proteins may affect beta-cell mass and insulin secretion capacity in humans, and render subjects at a higher risk of developing type 2 diabetes when exposed to environmental risk factors. Show less
Heparanase is the major enzyme involved in degradation of endothelial heparan sulfates, which is associated with impaired endothelial nitric oxide synthesis. However, the effect of heparan sulfate cha Show more
Heparanase is the major enzyme involved in degradation of endothelial heparan sulfates, which is associated with impaired endothelial nitric oxide synthesis. However, the effect of heparan sulfate chain length in relation to endothelial function and nitric oxide availability has never been investigated. We studied the effect of heterozygous mutations in heparan sulfate elongation genes EXT1 and EXT2 on endothelial function in vitro as well as in vivo. Flow-mediated dilation, a marker of nitric oxide bioavailability, was studied in Ext1(+/-) and Ext2(+/-) mice versus controls (n=7 per group), as well as in human subjects with heterozygous loss of function mutations in EXT1 and EXT2 (n=13 hereditary multiple exostoses and n=13 controls). Endothelial function was measured in microvascular endothelial cells under laminar flow with or without siRNA targeting EXT1 or EXT2. Endothelial glycocalyx and maximal arteriolar dilatation were significantly altered in Ext1(+/-) and Ext2(+/-) mice compared to wild-type littermates (glycocalyx: wild-type 0.67±0.1 μm, Ext1(+/-) 0.28±0.1 μm and Ext2(+/-) 0.25±0.1 μm, P<0.01, maximal arteriolar dilation during reperfusion: wild-type 11.3±1.0%), Ext1(+/-) 15.2±1.4% and Ext2(+/-) 13.8±1.6% P<0.05). In humans, brachial artery flow-mediated dilation was significantly increased in hereditary multiple exostoses patients (hereditary multiple exostoses 8.1±0.8% versus control 5.6±0.7%, P<0.05). In line, silencing of microvascular endothelial cell EXT1 and EXT2 under flow led to significant upregulation of endothelial nitric oxide synthesis and phospho-endothelial nitric oxide synthesis protein expression. Our data implicate that heparan sulfate elongation genes EXT1 and EXT2 are involved in maintaining endothelial homeostasis, presumably via increased nitric oxide bioavailability. Show less
Hair follicle (HF) morphogenesis and cycling are a result of intricate autonomous epithelial-mesenchymal interactions. Once the first HF cycle is complete it repeatedly undergoes cyclic transformation Show more
Hair follicle (HF) morphogenesis and cycling are a result of intricate autonomous epithelial-mesenchymal interactions. Once the first HF cycle is complete it repeatedly undergoes cyclic transformations. Heparan sulfate (HS) proteoglycans are found on the cell surface and in the extracellular matrix where they influence a variety of biological processes by interacting with physiologically important proteins, such as growth factors. Inhibition of heparanase (an HS endoglycosidase) in in vitro cultured HFs has been shown to induce a catagen-like process. Therefore, this study aimed to elucidate the precise role of HS in HF morphogenesis and cycling. An inducible tetratransgenic mouse model was generated to excise exostosin glycosyltransferase 1 (Ext1) in keratin 14-positive cells from P21. Interestingly, EXT1(StEpiΔ/StEpiΔ) mice presented solely anagen HFs. Moreover, waxing the fur to synchronize the HFs revealed accelerated hair regrowth in the EXT1(StEpiΔ/StEpiΔ) mice and hindered cycling into catagen. The ablation of HS in the interfollicular epidermal cells of mature skin led to the spontaneous formation of new HFs and an increase in Sonic Hedgehog expression resembling wild-type mice at P0, thereby indicating that the HS/Sonic Hedgehog signaling pathway regulates HF formation during embryogenesis and prevents HF formation in mature skin. Finally, the knock-out of HS also led to the morphogenesis and hyperplasia of sebaceous glands and sweat glands in mature mice, leading to exacerbated sebum production and accumulation on the skin surface. Therefore, our findings clearly show that an intricate control of HS levels is required for HF, sebaceous gland, and sweat gland morphogenesis and HF cycling. Show less
The heparan sulfate proteoglycan (HSPG) syndecan-1 (SDC1) acts as a major receptor for triglyceride-rich lipoprotein (TRL) clearance in the liver. We sought to identify the relevant apolipoproteins on Show more
The heparan sulfate proteoglycan (HSPG) syndecan-1 (SDC1) acts as a major receptor for triglyceride-rich lipoprotein (TRL) clearance in the liver. We sought to identify the relevant apolipoproteins on TRLs that mediate binding to SDC1 and determine their clinical relevance. Evidence supporting ApoE as a major determinant arose from its enrichment in TRLs from mice defective in hepatic heparan sulfate (Ndst1f/fAlbCre⁺ mice), decreased binding of ApoE-deficient TRLs to HSPGs on human hepatoma cells, and decreased clearance of ApoE-deficient [³H]TRLs in vivo. Evidence for a second ligand was suggested by the faster clearance of ApoE-deficient TRLs after injection into WT Ndst1f/fAlbCre⁻ versus mutant Ndst1f/fAlbCre⁺ mice and elevated fasting and postprandial plasma triglycerides in compound Apoe⁻/⁻Ndst1f/fAlbCre⁺ mice compared with either single mutant. ApoAV emerged as a candidate based on 6-fold enrichment of ApoAV in TRLs accumulating in Ndst1f/fAlbCre⁺ mice, decreased binding of TRLs to proteoglycans after depletion of ApoAV or addition of anti-ApoAV mAb, and decreased heparan sulfate-dependent binding of ApoAV-deficient particles to hepatocytes. Importantly, disruption of hepatic heparan sulfate-mediated clearance increased atherosclerosis. We conclude that clearance of TRLs by hepatic HSPGs is atheroprotective and mediated by multivalent binding to ApoE and ApoAV. Show less
Genetic loci for body mass index (BMI) in adolescence and young adulthood, a period of high risk for weight gain, are understudied, yet may yield important insight into the etiology of obesity and ear Show more
Genetic loci for body mass index (BMI) in adolescence and young adulthood, a period of high risk for weight gain, are understudied, yet may yield important insight into the etiology of obesity and early intervention. To identify novel genetic loci and examine the influence of known loci on BMI during this critical time period in late adolescence and early adulthood, we performed a two-stage meta-analysis using 14 genome-wide association studies in populations of European ancestry with data on BMI between ages 16 and 25 in up to 29 880 individuals. We identified seven independent loci (P < 5.0 × 10⁻⁸) near FTO (P = 3.72 × 10⁻²³), TMEM18 (P = 3.24 × 10⁻¹⁷), MC4R (P = 4.41 × 10⁻¹⁷), TNNI3K (P = 4.32 × 10⁻¹¹), SEC16B (P = 6.24 × 10⁻⁹), GNPDA2 (P = 1.11 × 10⁻⁸) and POMC (P = 4.94 × 10⁻⁸) as well as a potential secondary signal at the POMC locus (rs2118404, P = 2.4 × 10⁻⁵ after conditioning on the established single-nucleotide polymorphism at this locus) in adolescents and young adults. To evaluate the impact of the established genetic loci on BMI at these young ages, we examined differences between the effect sizes of 32 published BMI loci in European adult populations (aged 18-90) and those observed in our adolescent and young adult meta-analysis. Four loci (near PRKD1, TNNI3K, SEC16B and CADM2) had larger effects and one locus (near SH2B1) had a smaller effect on BMI during adolescence and young adulthood compared with older adults (P < 0.05). These results suggest that genetic loci for BMI can vary in their effects across the life course, underlying the importance of evaluating BMI at different ages. Show less
Inherited defects in the ability to catabolize glycosaminoglycans result in lysosomal storage disorders known as mucopolysaccharidoses (MPS), causing severe pathology, particularly in the brain. Enzym Show more
Inherited defects in the ability to catabolize glycosaminoglycans result in lysosomal storage disorders known as mucopolysaccharidoses (MPS), causing severe pathology, particularly in the brain. Enzyme replacement therapy has been used to treat mucopolysaccharidoses; however, neuropathology has remained refractory to this approach. To test directly whether substrate reduction might be feasible for treating MPS disease, we developed a genetic model for substrate reduction therapy by crossing MPS IIIa mice with animals partially deficient in heparan sulfate biosynthesis due to heterozygosity in Ext1 and Ext2, genes that encode the copolymerase required for heparan sulfate chain assembly. Reduction of heparan sulfate by 30-50% using this genetic strategy ameliorated the amount of disease-specific biomarker and pathology in multiple tissues, including the brain. In addition, we were able to demonstrate that substrate reduction therapy can improve the efficacy of enzyme replacement therapy in cell culture and in mice. These results provide proof of principle that targeted inhibition of heparan sulfate biosynthetic enzymes together with enzyme replacement might prove beneficial for treating mucopolysaccharidoses. Show less
Chronic kidney disease (CKD) is an important public health problem with a genetic component. We performed genome-wide association studies in up to 130,600 European ancestry participants overall, and s Show more
Chronic kidney disease (CKD) is an important public health problem with a genetic component. We performed genome-wide association studies in up to 130,600 European ancestry participants overall, and stratified for key CKD risk factors. We uncovered 6 new loci in association with estimated glomerular filtration rate (eGFR), the primary clinical measure of CKD, in or near MPPED2, DDX1, SLC47A1, CDK12, CASP9, and INO80. Morpholino knockdown of mpped2 and casp9 in zebrafish embryos revealed podocyte and tubular abnormalities with altered dextran clearance, suggesting a role for these genes in renal function. By providing new insights into genes that regulate renal function, these results could further our understanding of the pathogenesis of CKD. Show less
Specific interactions of growth factors with heparan sulfate may function as "switches" to regulate stages of branching morphogenesis in developing mammalian organs, such as breast, lung, salivary gla Show more
Specific interactions of growth factors with heparan sulfate may function as "switches" to regulate stages of branching morphogenesis in developing mammalian organs, such as breast, lung, salivary gland and kidney, but the evidence derives mostly from studies of explanted tissues or cell culture (Shah et al., 2004). We recently provided in vivo evidence that inactivation of Ndst1, the predominant N-deacetylase/N-sulfotransferase gene essential for the formation of mature heparan sulfate, results in a highly specific defect in murine lobuloalveolar development (Crawford et al., 2010). Here, we demonstrate a highly penetrant dramatic defect in primary branching by mammary epithelial-specific inactivation of Ext1, a subunit of the copolymerase complex that catalyzes the formation of the heparan sulfate chain. In contrast to Ext1 deletion, inactivation of Hs2st (which encodes an enzyme required for 2-O-sulfation of uronic acids in heparan sulfate) did not inhibit ductal formation but displayed markedly decreased secondary and ductal side-branches as well as fewer bifurcated terminal end buds. Targeted conditional deletion of c-Met, the receptor for HGF, in mammary epithelial cells showed similar defects in secondary and ductal side-branching, but did not result in any apparent defect in bifurcation of terminal end buds. Although there is published evidence indicating a role for 2-O sulfation in HGF binding, primary epithelial cells isolated from Hs2st conditional deletions were able to activate Erk in the presence of HGF and there appeared to be only a slight reduction in HGF-mediated c-Met phosphorylation in these cells compared to control. Thus, both c-Met and Hs2st play important, but partly independent, roles in secondary and ductal side-branching. When considered together with previous studies of Ndst1-deficient glands, the data presented here raise the possibility of partially-independent regulation by heparan sulfate-dependent pathways of primary ductal branching, terminal end bud bifurcation, secondary branching, ductal side-branching and lobuloalveolar formation. Show less
Multiple Hereditary Exostoses (MHE) syndrome is caused by haploinsufficiency in Golgi-associated heparan sulfate polymerases EXT1 or EXT2 and is characterized by formation of exostoses next to growing Show more
Multiple Hereditary Exostoses (MHE) syndrome is caused by haploinsufficiency in Golgi-associated heparan sulfate polymerases EXT1 or EXT2 and is characterized by formation of exostoses next to growing long bones and other skeletal elements. Recent mouse studies have indicated that formation of stereotypic exostoses requires a complete loss of Ext expression, suggesting that a similar local loss of EXT function may underlie exostosis formation in patients. To further test this possibility and gain greater insights into pathogenic mechanisms, we created heterozygous Ext1(+/-) and compound Ext1(+/-)/Ext2(+/-) mice. Like Ext2(+/-) mice described previously (Stickens et al. Development 132:5055), Ext1(+/-) mice displayed rib-associated exostosis-like outgrowths only. However, compound heterozygous mice had nearly twice as many outgrowths and, more importantly, displayed stereotypic growth plate-like exostoses along their long bones. Ext1(+/-)Ext2(+/-) exostoses contained very low levels of immuno-detectable heparan sulfate, and Ext1(+/-)Ext2(+/-) chondrocytes, endothelial cells and fibroblasts in vitro produced shortened heparan sulfate chains compared to controls and responded less vigorously to exogenous factors such as FGF-18. We also found that rib outgrowths formed in Ext1(f/+)Col2Cre and Ext1(f/+)Dermo1Cre mice, suggesting that ectopic skeletal tissue can be induced by conditional Ext ablation in local chondrogenic and/or perichondrial cells. The study indicates that formation of stereotypic exostoses requires a significant, but not complete, loss of Ext expression and that exostosis incidence and phenotype are intimately sensitive to, and inversely related to, Ext expression. The data also indicate that the nature and organization of ectopic tissue may be influenced by site-specific anatomical cues and mechanisms. Show less
Obesity is globally prevalent and highly heritable, but its underlying genetic factors remain largely elusive. To identify genetic loci for obesity susceptibility, we examined associations between bod Show more
Obesity is globally prevalent and highly heritable, but its underlying genetic factors remain largely elusive. To identify genetic loci for obesity susceptibility, we examined associations between body mass index and ∼ 2.8 million SNPs in up to 123,865 individuals with targeted follow up of 42 SNPs in up to 125,931 additional individuals. We confirmed 14 known obesity susceptibility loci and identified 18 new loci associated with body mass index (P < 5 × 10⁻⁸), one of which includes a copy number variant near GPRC5B. Some loci (at MC4R, POMC, SH2B1 and BDNF) map near key hypothalamic regulators of energy balance, and one of these loci is near GIPR, an incretin receptor. Furthermore, genes in other newly associated loci may provide new insights into human body weight regulation. Show less
To identify loci for age at menarche, we performed a meta-analysis of 32 genome-wide association studies in 87,802 women of European descent, with replication in up to 14,731 women. In addition to the Show more
To identify loci for age at menarche, we performed a meta-analysis of 32 genome-wide association studies in 87,802 women of European descent, with replication in up to 14,731 women. In addition to the known loci at LIN28B (P = 5.4 × 10⁻⁶⁰) and 9q31.2 (P = 2.2 × 10⁻³³), we identified 30 new menarche loci (all P < 5 × 10⁻⁸) and found suggestive evidence for a further 10 loci (P < 1.9 × 10⁻⁶). The new loci included four previously associated with body mass index (in or near FTO, SEC16B, TRA2B and TMEM18), three in or near other genes implicated in energy homeostasis (BSX, CRTC1 and MCHR2) and three in or near genes implicated in hormonal regulation (INHBA, PCSK2 and RXRG). Ingenuity and gene-set enrichment pathway analyses identified coenzyme A and fatty acid biosynthesis as biological processes related to menarche timing. Show less