👤 David Clarke

Annotation of regulatory elements is essential for understanding mechanisms underlying gene regulation, particularly tissue-specific regulation in human and animals. Here, we characterize 274,682 enhancers and 25,975 promoters across 24 tissues from an adult female sheep using ChIP-seq, ATAC-seq, CAGE-seq, RRBS, WGBS, and RNA-seq. We identify seven neural development-related genes with over 10 enhancers in brain tissues, highlighting the role of tissue-specific regulation. Cis-regulatory enhancer-promoter combinations provide insights into tissue-specific enhancers, such as the cerebellum-specific enhancer (chr15: 57390520-57390685) regulating BDNF, which is expressed in both the cerebellum and cerebral cortex. Comparative analysis of enhancer-promoter combinations in human, mouse, pig, cattle, and sheep reveals ruminant-specific pathways, including pentose catabolism and long-chain fatty acid import regulation. A milk fat yield quantitative trait locus (QTL) identified within an enhancer interacts with the fat metabolism-related gene COMMD1, and a birth weight-associated QTL detected within a cerebellum-specific enhancer regulates XKR4. This study provides a robust framework for exploring cis-regulatory mechanisms and tissue-specific regulation, advancing the functional annotation of the sheep reference genome. Show less

Little is known about the importance of blood lipids for risk of myocardial infarction (MI) in Chinese vs. European populations. We compared the associations with MI of apolioprotein B (ApoB) vs. low-density lipoprotein cholesterol (LDL-C) and remnant-cholesterol (remnant-C) vs. triglycerides in the China Kadoorie Biobank (CKB) and UK Biobank (UKB). Plasma levels of LDL-C, high-density lipoprotein-cholesterol (HDL-C), apolipoprotein B (ApoB), apolipoprotein A1 (ApoA1), non-HDL-C, remnant-C, LDL-C/ApoB, and HDL-C/ApoA1 ratios were measured in a nested case-control study of MI (948 cases, 6101 controls) in CKB and a prospective study (5344 cases in 279 989 participants) in UKB. Associations of lipids with MI were assessed using logistic regression in CKB and Cox regression in UKB after adjustment for confounders and correction for regression dilution. The mean levels of LDL-C were about 30% lower in CKB than in UKB [2.3 (0.6) vs. 3.7 (0.8) mmol/L], but mean levels of HDL-C were comparable [1.3 (0.3) vs. 1.5 (0.4) mmol/L], as were those for triglycerides [1.8 (1.1) vs. 1.7 (1.1) mmol/L]. While the rate ratios (RRs) of MI for 1 SD higher usual levels of LDL-C in Chinese were about half those in Europeans (1.27; 1.13-1.44 vs. 1.55; 1.49-1.61), the corresponding RRs for ApoB or non-HDL with MI were comparable between Chinese and Europeans. The findings reinforce current guidelines for primary prevention of atherosclerotic cardiovascular disease (ASCVD) in China that advocate initiation of statin treatment in individuals at high-risk of ASCVD rather than high levels of LDL-C. Show less

Integrated analyses of plasma proteomics and genetic data in prospective studies can help assess the causal relevance of proteins, improve risk prediction, and discover novel protein drug targets for type 2 diabetes (T2D). We measured plasma levels of 2,923 proteins using Olink Explore among ∼2,000 randomly selected participants from China Kadoorie Biobank (CKB) without prior diabetes at baseline. Cox regression assessed associations of individual protein with incident T2D (n = 92 cases). Proteomic-based risk models were developed with discrimination, calibration, reclassification assessed using area under the curve (AUC), calibration plots, and net reclassification index (NRI), respectively. Two-sample Mendelian randomization (MR) analyses using cis-protein quantitative trait loci identified in a genome-wide association study of CKB and UK Biobank for specific proteins were conducted to assess their causal relevance for T2D, along with colocalization analyses to examine shared causal variants between proteins and T2D. Overall, 33 proteins were significantly associated (false discovery rate <0.05) with risk of incident T2D, including IGFBP1, GHR, and amylase. The addition of these 33 proteins to a conventional risk prediction model improved AUC from 0.77 (0.73-0.82) to 0.88 (0.85-0.91) and NRI by 38%, with predicted risks well calibrated with observed risks. MR analyses provided support for the causal relevance for T2D of ENTR1, LPL, and PON3, with replication of ENTR1 and LPL in Europeans using different genetic instruments. Moreover, colocalization analyses showed strong evidence (pH4 > 0.6) of shared genetic variants of LPL and PON3 with T2D. Proteomic analyses in Chinese adults identified novel associations of multiple proteins with T2D with strong genetic evidence supporting their causal relevance and potential as novel drug targets for prevention and treatment of T2D. Show less

Recognition of the importance of conventional lipid measures and the advent of novel lipid-lowering medications have prompted the need for more comprehensive lipid panels to guide use of emerging treatments for the prevention of coronary heart disease (CHD). This report assessed the relevance of 13 apolipoproteins measured using a single mass-spectrometry assay for risk of CHD in the PROCARDIS case-control study of CHD (941 cases/975 controls). The associations of apolipoproteins with CHD were assessed after adjustment for established risk factors and correction for statin use. Apolipoproteins were grouped into 4 lipid-related classes [lipoprotein(a), low-density lipoprotein cholesterol, high-density lipoprotein cholesterol, and triglycerides] and their associations with CHD were adjusted for established CHD risk factors and conventional lipids. Analyses of these apolipoproteins in a subset of the ASCOT trial (Anglo-Scandinavian Cardiac Outcomes Trial) were used to assess their within-person variability and to estimate a correction for statin use. The findings in the PROCARDIS study were compared with those for incident cardiovascular disease in the Bruneck prospective study (n=688), including new measurements of Apo(a). Triglyceride-carrying apolipoproteins (ApoC1, ApoC3, and ApoE) were most strongly associated with the risk of CHD (2- to 3-fold higher odds ratios for top versus bottom quintile) independent of conventional lipid measures. Likewise, ApoB was independently associated with a 2-fold higher odds ratios of CHD. Lipoprotein(a) was measured using peptides from the Apo(a)-kringle repeat and Apo(a)-constant regions, but neither of these associations differed from the association with conventionally measured lipoprotein(a). Among HDL-related apolipoproteins, ApoA4 and ApoM were inversely related to CHD, independent of conventional lipid measures. The disease associations with all apolipoproteins were directionally consistent in the PROCARDIS and Bruneck studies, with the exception of ApoM. Apolipoproteins were associated with CHD independent of conventional risk factors and lipids, suggesting apolipoproteins could help to identify patients with residual lipid-related risk and guide personalized approaches to CHD risk reduction. Show less

Genome and epigenome wide association studies identified variants in carnitine palmitoyltransferase 1a (CPT1a) that associate with lipid traits. The goal of this study was to determine the role of liver-specific CPT1a on hepatic lipid metabolism. Male and female liver-specific knockout (LKO) and littermate controls were placed on a low-fat or high-fat diet (60% kcal fat) for 15 weeks. Mice were necropsied after a 16 h fast, and tissues were collected for lipidomics, matrix-assisted laser desorption ionization mass spectrometry imaging, kinome analysis, RNA-sequencing, and protein expression by immunoblotting. Female LKO mice had increased serum alanine aminotransferase levels which were associated with greater deposition of hepatic lipids, while male mice were not affected by CPT1a deletion relative to male control mice. Mice with CPT1a deletion had reductions in DHA-containing phospholipids at the expense of monounsaturated fatty acids (MUFA)-containing phospholipids in whole liver and at the level of the lipid droplet (LD). Male and female LKO mice increased RNA levels of genes involved in LD lipolysis (Plin2, Cidec, G0S2) and in polyunsaturated fatty acid metabolism (Elovl5, Fads1, Elovl2), while only female LKO mice increased genes involved in inflammation (Ly6d, Mmp12, Cxcl2). Kinase profiling showed decreased protein kinase A activity, which coincided with increased PLIN2, PLIN5, and G0S2 protein levels and decreased triglyceride hydrolysis in LKO mice. Liver-specific deletion of CPT1a promotes sexually dimorphic steatotic liver disease (SLD) in mice, and here we have identified new mechanisms by which females are protected from HFD-induced liver injury. Show less

Genome and epigenome wide association studies identified variants in carnitine palmitoyltransferase 1a (CPT1a) that associate with lipid traits. The goal of this study was to determine the impact by which liver-specific CPT1a deletion impacts hepatic lipid metabolism. Six-to-eight-week old male and female liver-specific knockout (LKO) and littermate controls were placed on a low-fat or high-fat diet (HFD; 60% kcal fat) for 15 weeks. Mice were necropsied after a 16 hour fast, and tissues were collected for lipidomics, matrix-assisted laser desorption ionization mass spectrometry imaging (MALDI-MSI), kinome analysis, RNA-sequencing, and protein expression by immunoblotting. Female LKO mice had increased serum alanine aminotransferase (ALT) levels which were associated with greater deposition of hepatic lipids, while male mice were not affected by CPT1a deletion relative to male control mice. Mice with CPT1a deletion had reductions in DHA-containing phospholipids at the expense of monounsaturated fatty acids (MUFA)-containing phospholipids in both whole liver and at the level of the lipid droplet (LD). Male and female LKO mice increased RNA levels of genes involved in LD lipolysis ( Liver-specific deletion of CPT1a promotes sexually dimorphic steatotic liver disease (SLD) in mice, and here we have identified new mechanisms by which females are protected from HFD-induced liver injury. Show less

Maintaining a healthy lifestyle to reduce type 2 diabetes (T2D) risk is challenging and additional strategies for T2D prevention are needed. We evaluated several lipid control medications as potential therapeutic options for T2D prevention using tissue-specific predicted gene expression summary statistics in a two-sample Mendelian randomisation (MR) design. Large-scale European genome-wide summary statistics for lipids and T2D were leveraged in our multi-stage analysis to estimate changes in either lipid levels or T2D risk driven by tissue-specific predicted gene expression. We incorporated tissue-specific predicted gene expression summary statistics to proxy therapeutic effects of three lipid control medications [i.e., statins, icosapent ethyl (IPE), and proprotein convertase subtilisin/kexin type-9 inhibitors (PCSK-9i)] on T2D susceptibility using two-sample Mendelian randomisation (MR). IPE, as proxied via increased FADS1 expression, was predicted to lower triglycerides and was associated with a 53% reduced risk of T2D. Statins and PCSK-9i, as proxied by reduced HMGCR and PCSK9 expression, respectively, were predicted to lower LDL-C levels but were not associated with T2D susceptibility. Triglyceride lowering via IPE may reduce the risk of developing T2D in populations of European ancestry. However, experimental validation using animal models is needed to substantiate our results and to motivate randomized control trials (RCTs) for IPE as putative treatment for T2D prevention. Only summary statistics were used in this analysis. Funding information is detailed under Acknowledgments. Show less

With the widespread application of next-generation sequencing, the genetic landscape of uterine mesenchymal neoplasms has been evolving rapidly to include several recently identified fusion genes. Although chromosomal rearrangements involving the 10q22 and 17q21.31 loci have been reported in occasional uterine leiomyomas decades ago, the corresponding KAT6B::KANSL1 fusion has been only recently identified in 2 uterine tumors diagnosed as leiomyoma and leiomyosarcoma. We herein describe 13 uterine stromal neoplasms carrying a KAT6B::KANSL1 (n=11) and KAT6A::KANSL1 (n=2) fusion. Patient ages ranged from 33 to 81 years (median, 49 y). Tumor size was 2.6 to 23.5 cm (median, 8.2 cm). Nine tumors were myometrium-centered, and 3 had an intracavitary component. Original diagnoses were mostly low-grade endometrial stromal sarcoma (LG-ESS; 10 cases) with atypical features (limited CD10 expression, sex cord-like features, pericytic vasculature, and frequent myxoid changes). Treatment was hysterectomy±bilateral salpingo-oophorectomy (10), myomectomy (1), and curettage (2). Five patients were disease-free at 6 to 34 months, 3 (27%) died of disease at 2 to 47 months, and 3 were alive with disease at 2, 17, and 17 years. Histologically, most tumors showed variable overlap with LG-ESS, but they were generally well-circumscribed lacking the extensive permeative and angioinvasive growth typical of LG-ESS. They were composed of monotonous medium-sized oval and spindle cells arranged into diffuse sheets with prominent spiral-type arterioles and frequent pericytoma-like vascular pattern. Variable myxoid stromal changes were frequent. Mitotic activity ranged from 1 to >20 in 10 HPFs. Immunohistochemistry showed variable expression of CD10 (12/13), estrogen receptor (8/11), progesterone receptor (8/11), smooth muscle actin (9/11), desmin (4/12), h-caldesmon (2/10), calretinin (3/8), inhibin (1/7), WT1 (4/7), cyclin D1 (5/11; diffuse in only 1 case), and pankeratin (5/10). This series characterizes a KAT6B/A::KANSL1 fusion-positive uterine stromal neoplasm within the morphologic spectrum of LG-ESS but with atypical features. The relationship of these neoplasms to genuine LG-ESS remains unclear. This molecular subtype of uterine endometrial stromal sarcoma has the potential for an unfavorable clinical course despite the absence of widely invasive growth; nevertheless, analysis of more cases is necessary to delineate the phenotypic spectrum and biological potential of this tumor. Show less

The Melanocortin-4 Receptor (MC4R) plays a pivotal role in energy homeostasis. We used human MC4R mutations associated with an increased or decreased risk of obesity to dissect mechanisms that regulate MC4R function. Most obesity-associated mutations impair trafficking to the plasma membrane (PM), whereas obesity-protecting mutations either accelerate recycling to the PM or decrease internalization, resulting in enhanced signaling. MC4R mutations that do not affect canonical Gα Show less

Aortic stenosis (AS) has no approved medical treatment. Identifying etiological pathways for AS could identify pharmacological targets. To identify novel genetic loci and pathways associated with AS. This genome-wide association study used a case-control design to evaluate 44 703 participants (3469 cases of AS) of self-reported European ancestry from the Genetic Epidemiology Research on Adult Health and Aging (GERA) cohort (from January 1, 1996, to December 31, 2015). Replication was performed in 7 other cohorts totaling 256 926 participants (5926 cases of AS), with additional analyses performed in 6942 participants from the Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE) Consortium. Follow-up biomarker analyses with aortic valve calcium (AVC) were also performed. Data were analyzed from May 1, 2017, to December 5, 2019. Genetic variants (615 643 variants) and polyunsaturated fatty acids (ω-6 and ω-3) measured in blood samples. Aortic stenosis and aortic valve replacement defined by electronic health records, surgical records, or echocardiography and the presence of AVC measured by computed tomography. The mean (SD) age of the 44 703 GERA participants was 69.7 (8.4) years, and 22 019 (49.3%) were men. The rs174547 variant at the FADS1/2 locus was associated with AS (odds ratio [OR] per C allele, 0.88; 95% CI, 0.83-0.93; P = 3.0 × 10-6), with genome-wide significance after meta-analysis with 7 replication cohorts totaling 312 118 individuals (9395 cases of AS) (OR, 0.91; 95% CI, 0.88-0.94; P = 2.5 × 10-8). A consistent association with AVC was also observed (OR, 0.91; 95% CI, 0.83-0.99; P = .03). A higher ratio of arachidonic acid to linoleic acid was associated with AVC (OR per SD of the natural logarithm, 1.19; 95% CI, 1.09-1.30; P = 6.6 × 10-5). In mendelian randomization, increased FADS1 liver expression and arachidonic acid were associated with AS (OR per unit of normalized expression, 1.31 [95% CI, 1.17-1.48; P = 7.4 × 10-6]; OR per 5-percentage point increase in arachidonic acid for AVC, 1.23 [95% CI, 1.01-1.49; P = .04]; OR per 5-percentage point increase in arachidonic acid for AS, 1.08 [95% CI, 1.04-1.13; P = 4.1 × 10-4]). Variation at the FADS1/2 locus was associated with AS and AVC. Findings from biomarker measurements and mendelian randomization appear to link ω-6 fatty acid biosynthesis to AS, which may represent a therapeutic target. Show less

The Iberian Peninsula stands out as having variable levels of population admixture and isolation, making Spain an interesting setting for studying the genetic architecture of neurodegenerative diseases. To perform the largest PD genome-wide association study restricted to a single country. We performed a GWAS for both risk of PD and age at onset in 7,849 Spanish individuals. Further analyses included population-specific risk haplotype assessments, polygenic risk scoring through machine learning, Mendelian randomization of expression, and methylation data to gain insight into disease-associated loci, heritability estimates, genetic correlations, and burden analyses. We identified a novel population-specific genome-wide association study signal at PARK2 associated with age at onset, which was likely dependent on the c.155delA mutation. We replicated four genome-wide independent signals associated with PD risk, including SNCA, LRRK2, KANSL1/MAPT, and HLA-DQB1. A significant trend for smaller risk haplotypes at known loci was found compared to similar studies of non-Spanish origin. Seventeen PD-related genes showed functional consequence by two-sample Mendelian randomization in expression and methylation data sets. Long runs of homozygosity at 28 known genes/loci were found to be enriched in cases versus controls. Our data demonstrate the utility of the Spanish risk haplotype substructure for future fine-mapping efforts, showing how leveraging unique and diverse population histories can benefit genetic studies of complex diseases. The present study points to PARK2 as a major hallmark of PD etiology in Spain. © 2019 International Parkinson and Movement Disorder Society. Show less

Peanut allergy (PA) is a complex disease with both environmental and genetic risk factors. Previously, PA loci were identified in filaggrin (FLG) and HLA in candidate gene studies, and loci in HLA were identified in a genome-wide association study and meta-analysis. We sought to investigate genetic susceptibility to PA. Eight hundred fifty cases and 926 hyper-control subjects and more than 7.8 million genotyped and imputed single nucleotide polymorphisms (SNPs) were analyzed in a genome-wide association study to identify susceptibility variants for PA in the Canadian population. A meta-analysis of 2 phenotypes (PA and food allergy) was conducted by using 7 studies from the Canadian, American (n = 2), Australian, German, and Dutch (n = 2) populations. An SNP near integrin α6 (ITGA6) reached genome-wide significance with PA (P = 1.80 × 10 This study identifies multiple novel loci as risk factors for PA and food allergy and establishes C11orf30 as a risk locus for both PA and food allergy. Multiple genes (C11orf30/EMSY, SKAP1, and CTNNA3) identified by this study are involved in epigenetic regulation of gene expression. Show less

Increasing levels of high-density lipoprotein (HDL) cholesterol through pharmacologic inhibition of cholesteryl ester transfer protein (CETP) is a potentially important strategy for prevention and treatment of cardiovascular disease (CVD). To use genetic variants in the CETP gene to assess potential risks and benefits of lifelong lower CETP activity on CVD and other outcomes. This prospective biobank study included 151 217 individuals aged 30 to 79 years who were enrolled from 5 urban and 5 rural areas of China from June 25, 2004, through July 15, 2008. All participants had baseline genotype data, 17 854 of whom had lipid measurements and 4657 of whom had lipoprotein particle measurements. Median follow-up of 9.2 years (interquartile range, 8.2-10.1 years) was completed January 1, 2016, through linkage to health insurance records and death and disease registries. Five CETP variants, including an East Asian loss-of-function variant (rs2303790), combined in a genetic score weighted to associations with HDL cholesterol levels. Baseline levels of lipids and lipoprotein particles, cardiovascular risk factors, incidence of carotid plaque and predefined major vascular and nonvascular diseases, and a phenome-wide range of diseases. Among the 151 217 individuals included in this study (58.4% women and 41.6% men), the mean (SD) age was 52.3 (10.9) years. Overall, the mean (SD) low-density lipoprotein (LDL) cholesterol level was 91 (27) mg/dL; HDL cholesterol level, 48 (12) mg/dL. CETP variants were strongly associated with higher concentrations of HDL cholesterol (eg, 6.1 [SE, 0.4] mg/dL per rs2303790-G allele; P = 9.4 × 10-47) but were not associated with lower LDL cholesterol levels. Within HDL particles, cholesterol esters were increased and triglycerides reduced, whereas within very low-density lipoprotein particles, cholesterol esters were reduced and triglycerides increased. When scaled to 10-mg/dL higher levels of HDL cholesterol, the CETP genetic score was not associated with occlusive CVD (18 550 events; odds ratio [OR], 0.98; 95% CI, 0.91-1.06), major coronary events (5767 events; OR, 1.08; 95% CI, 0.95-1.22), myocardial infarction (3118 events; OR, 1.14; 95% CI, 0.97-1.35), ischemic stroke (13 759 events; OR, 0.94; 95% CI, 0.86-1.02), intracerebral hemorrhage (6532 events; OR, 0.94; 95% CI, 0.83-1.06), or other vascular diseases or carotid plaque. Similarly, rs2303790 was not associated with any vascular diseases or plaque. No associations with nonvascular diseases were found other than an increased risk for eye diseases with rs2303790 (4090 events; OR, 1.43; 95% CI, 1.13-1.80; P = .003). CETP variants were associated with altered HDL metabolism but did not lower LDL cholesterol levels and had no significant association with risk for CVD. These results suggest that in the absence of reduced LDL cholesterol levels, increasing HDL cholesterol levels by inhibition of CETP may not confer significant benefits for CVD. Show less

Deregulation of the cyclin-dependent kinases (CDKs) has been implicated in the pathogenesis of multiple cancer types. Consequently, CDKs have garnered intense interest as therapeutic targets for the treatment of cancer. We describe herein the molecular and cellular effects of CCT068127, a novel inhibitor of CDK2 and CDK9. Optimized from the purine template of seliciclib, CCT068127 exhibits greater potency and selectivity against purified CDK2 and CDK9 and superior antiproliferative activity against human colon cancer and melanoma cell lines. X-ray crystallography studies reveal that hydrogen bonding with the DFG motif of CDK2 is the likely mechanism of greater enzymatic potency. Commensurate with inhibition of CDK activity, CCT068127 treatment results in decreased retinoblastoma protein (RB) phosphorylation, reduced phosphorylation of RNA polymerase II, and induction of cell cycle arrest and apoptosis. The transcriptional signature of CCT068127 shows greatest similarity to other small-molecule CDK and also HDAC inhibitors. CCT068127 caused a dramatic loss in expression of DUSP6 phosphatase, alongside elevated ERK phosphorylation and activation of MAPK pathway target genes. MCL1 protein levels are rapidly decreased by CCT068127 treatment and this associates with synergistic antiproliferative activity after combined treatment with CCT068127 and ABT263, a BCL2 family inhibitor. These findings support the rational combination of this series of CDK2/9 inhibitors and BCL2 family inhibitors for the treatment of human cancer. Show less

Spinal cord injury (SCI) results in loss of movement, sensibility, and autonomic control at the level of the lesion and at lower parts of the body. Several experimental strategies have been used in attempts to increase endogenous mechanisms of neuroprotection, neuroplasticity, and repair, but with limited success. It is known that glucose-dependent insulinotropic peptide (GIP) and its receptor (GIPR) can enhance synaptic plasticity, neurogenesis, and axonal outgrowth. However, their role in the injury has never been studied. The aim of this study was to evaluate the changes in expression levels of both GIP and GIPR in acute and chronic phases of SCI in rats. Following SCI (2 to 24 h after damage), the rat spinal cord showed a lesion in which the epicenter had a cavity with hemorrhage and necrosis. Furthermore, the lesion cavity also showed ballooned cells 14 and 28 days after injury. We found that SCI induced increases in GIPR expression in areas neighboring the site of injury at 6 h and 28 days after the injury. Moreover, higher GIP expression was observed in these regions on day 28. Neuronal projections from the injury epicenter showed an increase in GIP immunoreactivity 24 h and 14 and 28 days after SCI. Interestingly, GIP was also found in progenitor cells at the spinal cord canal 24 h after injury, whereas both GIP and GIPR were present in progenitor cells at the injury epicenter 14 days after in SCI animals. These results suggest that GIP and its receptor might be implicated with neurogenesis and the repair process after SCI. Show less

Molecular evolution has focused on the divergence of molecular functions, yet we know little about how structurally distinct protein folds emerge de novo. We characterized the evolutionary trajectories and selection forces underlying emergence of β-propeller proteins, a globular and symmetric fold group with diverse functions. The identification of short propeller-like motifs (<50 amino acids) in natural genomes indicated that they expanded via tandem duplications to form extant propellers. We phylogenetically reconstructed 47-residue ancestral motifs that form five-bladed lectin propellers via oligomeric assembly. We demonstrate a functional trajectory of tandem duplications of these motifs leading to monomeric lectins. Foldability, i.e., higher efficiency of folding, was the main parameter leading to improved functionality along the entire evolutionary trajectory. However, folding constraints changed along the trajectory: initially, conflicts between monomer folding and oligomer assembly dominated, whereas subsequently, upon tandem duplication, tradeoffs between monomer stability and foldability took precedence. Show less

APOE ɛ4, the most significant genetic risk factor for Alzheimer disease (AD), may mask effects of other loci. We re-analyzed genome-wide association study (GWAS) data from the International Genomics of Alzheimer's Project (IGAP) Consortium in APOE ɛ4+ (10 352 cases and 9207 controls) and APOE ɛ4- (7184 cases and 26 968 controls) subgroups as well as in the total sample testing for interaction between a single-nucleotide polymorphism (SNP) and APOE ɛ4 status. Suggestive associations (P<1 × 10(-4)) in stage 1 were evaluated in an independent sample (stage 2) containing 4203 subjects (APOE ɛ4+: 1250 cases and 536 controls; APOE ɛ4-: 718 cases and 1699 controls). Among APOE ɛ4- subjects, novel genome-wide significant (GWS) association was observed with 17 SNPs (all between KANSL1 and LRRC37A on chromosome 17 near MAPT) in a meta-analysis of the stage 1 and stage 2 data sets (best SNP, rs2732703, P=5·8 × 10(-9)). Conditional analysis revealed that rs2732703 accounted for association signals in the entire 100-kilobase region that includes MAPT. Except for previously identified AD loci showing stronger association in APOE ɛ4+ subjects (CR1 and CLU) or APOE ɛ4- subjects (MS4A6A/MS4A4A/MS4A6E), no other SNPs were significantly associated with AD in a specific APOE genotype subgroup. In addition, the finding in the stage 1 sample that AD risk is significantly influenced by the interaction of APOE with rs1595014 in TMEM106B (P=1·6 × 10(-7)) is noteworthy, because TMEM106B variants have previously been associated with risk of frontotemporal dementia. Expression quantitative trait locus analysis revealed that rs113986870, one of the GWS SNPs near rs2732703, is significantly associated with four KANSL1 probes that target transcription of the first translated exon and an untranslated exon in hippocampus (P ⩽ 1.3 × 10(-8)), frontal cortex (P ⩽ 1.3 × 10(-9)) and temporal cortex (P⩽1.2 × 10(-11)). Rs113986870 is also strongly associated with a MAPT probe that targets transcription of alternatively spliced exon 3 in frontal cortex (P=9.2 × 10(-6)) and temporal cortex (P=2.6 × 10(-6)). Our APOE-stratified GWAS is the first to show GWS association for AD with SNPs in the chromosome 17q21.31 region. Replication of this finding in independent samples is needed to verify that SNPs in this region have significantly stronger effects on AD risk in persons lacking APOE ɛ4 compared with persons carrying this allele, and if this is found to hold, further examination of this region and studies aimed at deciphering the mechanism(s) are warranted. Show less

Myocardial infarction (MI), a leading cause of death around the world, displays a complex pattern of inheritance. When MI occurs early in life, genetic inheritance is a major component to risk. Previously, rare mutations in low-density lipoprotein (LDL) genes have been shown to contribute to MI risk in individual families, whereas common variants at more than 45 loci have been associated with MI risk in the population. Here we evaluate how rare mutations contribute to early-onset MI risk in the population. We sequenced the protein-coding regions of 9,793 genomes from patients with MI at an early age (≤50 years in males and ≤60 years in females) along with MI-free controls. We identified two genes in which rare coding-sequence mutations were more frequent in MI cases versus controls at exome-wide significance. At low-density lipoprotein receptor (LDLR), carriers of rare non-synonymous mutations were at 4.2-fold increased risk for MI; carriers of null alleles at LDLR were at even higher risk (13-fold difference). Approximately 2% of early MI cases harbour a rare, damaging mutation in LDLR; this estimate is similar to one made more than 40 years ago using an analysis of total cholesterol. Among controls, about 1 in 217 carried an LDLR coding-sequence mutation and had plasma LDL cholesterol > 190 mg dl(-1). At apolipoprotein A-V (APOA5), carriers of rare non-synonymous mutations were at 2.2-fold increased risk for MI. When compared with non-carriers, LDLR mutation carriers had higher plasma LDL cholesterol, whereas APOA5 mutation carriers had higher plasma triglycerides. Recent evidence has connected MI risk with coding-sequence mutations at two genes functionally related to APOA5, namely lipoprotein lipase and apolipoprotein C-III (refs 18, 19). Combined, these observations suggest that, as well as LDL cholesterol, disordered metabolism of triglyceride-rich lipoproteins contributes to MI risk. Show less

There is great interest in the relationship between polyunsaturated fatty acids and health. Yet, the combinatory effect of factors such as sex, ethnicity, genetic polymorphisms and hormonal contraceptives (HC) on the concentrations of these fatty acids is unknown. Therefore, we sought to determine the effects of FADS polymorphisms, and HC use in females, on aggregate desaturase indices (ADI), and plasma docosahexaenoic acid (DHA) concentrations in Caucasian and East Asian males and females. Fasting plasma samples were collected from subjects (Caucasian males: 113 and females: 298; East Asian males: 98 and females: 277) from the Toronto Nutrigenomics and Health Study. Fatty acid concentrations were measured by gas chromatography. ADI were estimated by dividing concentrations of arachidonic acid by linoleic acid (n-6 ADI) and eicosapentaenoic acid (EPA) by α-linolenic acid (n-3 ADI). [DHA/EPA] desaturase index was used to determine effects of FADS2 polymorphisms and HC use on EPA conversion to DHA. In Caucasians, associations between n-6 ADI and multiple SNP (FADS1 rs174547, FADS2 rs174576, and rs174611 in males; FADS1 rs174547, FADS2 rs174570, rs174576, rs174679, rs174611, rs174593, rs174626, rs2072114, rs2845573, and rs2851682 in females) withstood multiple testing. In East Asian females, 5 SNP-n-6 ADI associations (FADS2 rs174602, rs174626, rs2072114, rs2845573, and rs2851682) withstood multiple testing. One FADS2 SNP was associated with altered [DHA/EPA] desaturase index in Caucasian females only (rs174576, p < 0.0001). HC use had a significant effect on DHA concentrations in Caucasian females only (P < 0.0001). We demonstrate ethnic- and sex-specific effects of FADS polymorphisms on desaturase indices, and ethnic-specific effect of HC use on plasma DHA concentrations. Show less

DIRAS3 is an imprinted tumor suppressor gene that is downregulated in 60% of human ovarian cancers. Re-expression of DIRAS3 at physiological levels inhibits proliferation, decreases motility, induces autophagy, and regulates tumor dormancy. Functional inhibition of autophagy with choroquine in dormant xenografts that express DIRAS3 significantly delays tumor regrowth after DIRAS3 levels are reduced, suggesting that autophagy sustains dormant ovarian cancer cells. This study documents a newly discovered role for DIRAS3 in forming the autophagosome initiation complex (AIC) that contains BECN1, PIK3C3, PIK3R4, ATG14, and DIRAS3. Participation of BECN1 in the AIC is inhibited by binding of BECN1 homodimers to BCL2. DIRAS3 binds BECN1, disrupting BECN1 homodimers and displacing BCL2. Binding of DIRAS3 to BECN1 increases the association of BECN1 with PIK3C3 and ATG14, facilitating AIC activation. Amino acid starvation of cells induces DIRAS3 expression, reduces BECN1-BCL2 interaction and promotes autophagy, whereas DIRAS3 depletion blocks amino acid starvation-induced autophagy. In primary ovarian cancers, punctate expression of DIRAS3, BECN1, and the autophagic biomarker MAP1LC3 are highly correlated (P<0.0001), underlining the clinical relevance of these mechanistic studies. Punctate expression of DIRAS3 and MAP1LC3 was detected in only 21-23% of primary ovarian cancers but in 81-84% of tumor nodules found on the peritoneal surface at second-look operations following primary chemotherapy. This reflects a 4-fold increase (P<0.0001) in autophagy between primary disease and post-treatment recurrence. We suggest that DIRAS3 not only regulates the AIC, but induces autophagy in dormant, nutrient-deprived ovarian cancer cells that remain after conventional chemotherapy, facilitating their survival. Show less

The strong observational association between total homocysteine (tHcy) concentrations and risk of coronary artery disease (CAD) and the null associations in the homocysteine-lowering trials have prompted the need to identify genetic variants associated with homocysteine concentrations and risk of CAD. We tested whether common genetic polymorphisms associated with variation in tHcy are also associated with CAD. We conducted a meta-analysis of genome-wide association studies (GWAS) on tHcy concentrations in 44,147 individuals of European descent. Polymorphisms associated with tHcy (P < 10(⁻⁸) were tested for association with CAD in 31,400 cases and 92,927 controls. Common variants at 13 loci, explaining 5.9% of the variation in tHcy, were associated with tHcy concentrations, including 6 novel loci in or near MMACHC (2.1 × 10⁻⁹), SLC17A3 (1.0 × 10⁻⁸), GTPB10 (1.7 × 10⁻⁸), CUBN (7.5 × 10⁻¹⁰), HNF1A (1.2 × 10⁻¹²)), and FUT2 (6.6 × 10⁻⁹), and variants previously reported at or near the MTHFR, MTR, CPS1, MUT, NOX4, DPEP1, and CBS genes. Individuals within the highest 10% of the genotype risk score (GRS) had 3-μmol/L higher mean tHcy concentrations than did those within the lowest 10% of the GRS (P = 1 × 10⁻³⁶). The GRS was not associated with risk of CAD (OR: 1.01; 95% CI: 0.98, 1.04; P = 0.49). We identified several novel loci that influence plasma tHcy concentrations. Overall, common genetic variants that influence plasma tHcy concentrations are not associated with risk of CAD in white populations, which further refutes the causal relevance of moderately elevated tHcy concentrations and tHcy-related pathways for CAD. Show less

Congenital hypogonadotropic hypogonadism (CHH) and its anosmia-associated form (Kallmann syndrome [KS]) are genetically heterogeneous. Among the >15 genes implicated in these conditions, mutations in FGF8 and FGFR1 account for ~12% of cases; notably, KAL1 and HS6ST1 are also involved in FGFR1 signaling and can be mutated in CHH. We therefore hypothesized that mutations in genes encoding a broader range of modulators of the FGFR1 pathway might contribute to the genetics of CHH as causal or modifier mutations. Thus, we aimed to (1) investigate whether CHH individuals harbor mutations in members of the so-called "FGF8 synexpression" group and (2) validate the ability of a bioinformatics algorithm on the basis of protein-protein interactome data (interactome-based affiliation scoring [IBAS]) to identify high-quality candidate genes. On the basis of sequence homology, expression, and structural and functional data, seven genes were selected and sequenced in 386 unrelated CHH individuals and 155 controls. Except for FGF18 and SPRY2, all other genes were found to be mutated in CHH individuals: FGF17 (n = 3 individuals), IL17RD (n = 8), DUSP6 (n = 5), SPRY4 (n = 14), and FLRT3 (n = 3). Independently, IBAS predicted FGF17 and IL17RD as the two top candidates in the entire proteome on the basis of a statistical test of their protein-protein interaction patterns to proteins known to be altered in CHH. Most of the FGF17 and IL17RD mutations altered protein function in vitro. IL17RD mutations were found only in KS individuals and were strongly linked to hearing loss (6/8 individuals). Mutations in genes encoding components of the FGF pathway are associated with complex modes of CHH inheritance and act primarily as contributors to an oligogenic genetic architecture underlying CHH. Show less

Mutations in components of the Wnt signaling pathway, including β-catenin and AXIN1, are found in more than 50% of human hepatocellular carcinomas (HCCs). Disruption of Axin1 causes embryonic lethality in mice. We generated mice with conditional disruption of Axin1 to study its function specifically in adult liver. Mice with a LoxP-flanked allele of Axin1 were generated by homologous recombination. Mice homozygous for the Axin1fl/fl allele were crossed with AhCre mice; in offspring, Axin1 was disrupted in liver following injection of β-naphthoflavone (Axin1fl/fl/Cre mice). Liver tissues were collected and analyzed by quantitative real-time polymerase chain reaction and immunoprecipitation, histology, and immunoblot assays. Deletion of Axin1 from livers of adult mice resulted in an acute and persistent increase in hepatocyte cell volume, proliferation, and transcription of genes that induce the G(2)/M transition in the cell cycle and cytokinesis. A subset of Wnt target genes was activated, including Axin2, c-Myc, and cyclin D1. However, loss of Axin1 did not increase nuclear levels of β-catenin or cause changes in liver zonation that have been associated with loss of the adenomatous polyposis coli (APC) or constitutive activation of β-catenin. After 1 year, 5 of 9 Axin1fl/fl/Cre mice developed liver tumors with histologic features of HCC. Hepatocytes from adult mice with conditional disruption of Axin1 in liver have a transcriptional profile that differs from that associated with loss of APC or constitutive activation of β-catenin. It might be similar to a proliferation profile observed in a subset of human HCCs with mutations in AXIN1. Axin1fl/fl mice could be a useful model of AXIN1-associated tumorigenesis and HCC. Show less

Recent evidence indicates that genetic variation in fatty acid desaturases 1 and 2 (FADS1 and FADS2) is associated with changes in plasma fatty acid profiles; however, the association with altered desaturase activity has not been examined in different ethnic populations. The present study examined whether genetic variation in the FADS gene cluster regulates desaturase activity in two populations of young Canadian adults (Caucasian and Asian) and whether altered desaturase activity was reflected in both n-3 and n-6 fatty acid profiles. FADS1 and FADS2 were genotyped in a random subset of participants (Caucasian, n=78; Asian, n=69) from the Toronto Nutrigenomics and Health study using MALDI-TOF mass spectrometry, and plasma fatty acids were measured by gas chromatography. Desaturase activities were estimated using the following fatty acid ratios: γ-linoleic acid to linoleic acid (GLA:LA), arachidonic acid to linoleic acid (AA:LA), arachidonic acid to dihomo-γ-linoleic acid (AA:DGLA), and eicosapentaneoic acid to α-linolenic acid (EPA:ALA). Nineteen single nucleotide polymorphisms (SNPs) were examined, and several SNPs (9 in Caucasians and 8 in Asians) were associated with various desaturase activities. The most significant association detected was between the FADS1 rs174547 SNP and AA:LA in both Caucasians (p=4.0 × 10(-8)) and Asians (p=5.0 × 10(-5)). Although the minor allele for this SNP differed between Caucasians (T) and Asians (C), carriers of the C allele had a lower desaturase activity than carriers of the T allele in both groups. To determine whether rs174547 was a dominant SNP in the FADS gene cluster, we constructed an additional model which included this SNP as a covariate. Only one SNP (rs498793 in FADS2) remained associated with the EPA:ALA ratio (p=1.1 × 10(-5)) in Asians. This study shows that genetic variation in the FADS gene cluster (in particular rs174547) can alter desaturase activity in subjects of Caucasians and Asian descent. Show less

Proinsulin is a precursor of mature insulin and C-peptide. Higher circulating proinsulin levels are associated with impaired β-cell function, raised glucose levels, insulin resistance, and type 2 diabetes (T2D). Studies of the insulin processing pathway could provide new insights about T2D pathophysiology. We have conducted a meta-analysis of genome-wide association tests of ∼2.5 million genotyped or imputed single nucleotide polymorphisms (SNPs) and fasting proinsulin levels in 10,701 nondiabetic adults of European ancestry, with follow-up of 23 loci in up to 16,378 individuals, using additive genetic models adjusted for age, sex, fasting insulin, and study-specific covariates. Nine SNPs at eight loci were associated with proinsulin levels (P < 5 × 10(-8)). Two loci (LARP6 and SGSM2) have not been previously related to metabolic traits, one (MADD) has been associated with fasting glucose, one (PCSK1) has been implicated in obesity, and four (TCF7L2, SLC30A8, VPS13C/C2CD4A/B, and ARAP1, formerly CENTD2) increase T2D risk. The proinsulin-raising allele of ARAP1 was associated with a lower fasting glucose (P = 1.7 × 10(-4)), improved β-cell function (P = 1.1 × 10(-5)), and lower risk of T2D (odds ratio 0.88; P = 7.8 × 10(-6)). Notably, PCSK1 encodes the protein prohormone convertase 1/3, the first enzyme in the insulin processing pathway. A genotype score composed of the nine proinsulin-raising alleles was not associated with coronary disease in two large case-control datasets. We have identified nine genetic variants associated with fasting proinsulin. Our findings illuminate the biology underlying glucose homeostasis and T2D development in humans and argue against a direct role of proinsulin in coronary artery disease pathogenesis. 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 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

Levels of circulating glucose are tightly regulated. To identify new loci influencing glycemic traits, we performed meta-analyses of 21 genome-wide association studies informative for fasting glucose, fasting insulin and indices of beta-cell function (HOMA-B) and insulin resistance (HOMA-IR) in up to 46,186 nondiabetic participants. Follow-up of 25 loci in up to 76,558 additional subjects identified 16 loci associated with fasting glucose and HOMA-B and two loci associated with fasting insulin and HOMA-IR. These include nine loci newly associated with fasting glucose (in or near ADCY5, MADD, ADRA2A, CRY2, FADS1, GLIS3, SLC2A2, PROX1 and C2CD4B) and one influencing fasting insulin and HOMA-IR (near IGF1). We also demonstrated association of ADCY5, PROX1, GCK, GCKR and DGKB-TMEM195 with type 2 diabetes. Within these loci, likely biological candidate genes influence signal transduction, cell proliferation, development, glucose-sensing and circadian regulation. Our results demonstrate that genetic studies of glycemic traits can identify type 2 diabetes risk loci, as well as loci containing gene variants that are associated with a modest elevation in glucose levels but are not associated with overt diabetes. Show less

The concept that optic nerve fiber loss might be reduced by neuroprotection arose in the mid 1990s. The subsequent research effort, focused mainly on rodent models, has not yet transformed into a successful clinical trial, but provides mechanistic understanding of retinal ganglion cell death and points to potential therapeutic strategies. This review highlights advances made over the last year. In excitotoxicity and axotomy models retinal ganglion cell death has been shown to result from a complex interaction between retinal neurons and Müller glia, which release toxic molecules including tumor necrosis factor alpha. This counteracts neuroprotection by neurotrophins such as nerve growth factor, which bind to p75NTR receptors on Müller glia stimulating the toxic release. Another negative effect against neurotrophin-mediated protection involves the action of LINGO-1 at trkB brain-derived neurotrophic factor (BDNF) receptors, and BDNF neuroprotection is enhanced by an antagonist to LINGO-1. As an alternative to pharmacotherapy, retinal defences can be stimulated by exposure to infrared radiation. The mechanisms involved in glaucoma and other optic nerve disorders are being clarified in rodent models, focusing on retrograde degeneration following axonal damage, excitotoxicity and inflammatory/autoimmune mechanisms. Neuroprotective strategies are being refined in the light of the mechanistic understanding. Show less

Blood low-density lipoprotein (LDL) cholesterol, high-density lipoprotein (HDL) cholesterol and triglyceride levels are risk factors for cardiovascular disease. To dissect the polygenic basis of these traits, we conducted genome-wide association screens in 19,840 individuals and replication in up to 20,623 individuals. We identified 30 distinct loci associated with lipoprotein concentrations (each with P < 5 x 10(-8)), including 11 loci that reached genome-wide significance for the first time. The 11 newly defined loci include common variants associated with LDL cholesterol near ABCG8, MAFB, HNF1A and TIMD4; with HDL cholesterol near ANGPTL4, FADS1-FADS2-FADS3, HNF4A, LCAT, PLTP and TTC39B; and with triglycerides near AMAC1L2, FADS1-FADS2-FADS3 and PLTP. The proportion of individuals exceeding clinical cut points for high LDL cholesterol, low HDL cholesterol and high triglycerides varied according to an allelic dosage score (P < 10(-15) for each trend). These results suggest that the cumulative effect of multiple common variants contributes to polygenic dyslipidemia. Show less

While conventional LDL-C, HDL-C, and triglyceride measurements reflect aggregate properties of plasma lipoprotein fractions, NMR-based measurements more accurately reflect lipoprotein particle concentrations according to class (LDL, HDL, and VLDL) and particle size (small, medium, and large). The concentrations of these lipoprotein sub-fractions may be related to risk of cardiovascular disease and related metabolic disorders. We performed a genome-wide association study of 17 lipoprotein measures determined by NMR together with LDL-C, HDL-C, triglycerides, ApoA1, and ApoB in 17,296 women from the Women's Genome Health Study (WGHS). Among 36 loci with genome-wide significance (P<5x10(-8)) in primary and secondary analysis, ten (PCCB/STAG1 (3q22.3), GMPR/MYLIP (6p22.3), BTNL2 (6p21.32), KLF14 (7q32.2), 8p23.1, JMJD1C (10q21.3), SBF2 (11p15.4), 12q23.2, CCDC92/DNAH10/ZNF664 (12q24.31.B), and WIPI1 (17q24.2)) have not been reported in prior genome-wide association studies for plasma lipid concentration. Associations with mean lipoprotein particle size but not cholesterol content were found for LDL at four loci (7q11.23, LPL (8p21.3), 12q24.31.B, and LIPG (18q21.1)) and for HDL at one locus (GCKR (2p23.3)). In addition, genetic determinants of total IDL and total VLDL concentration were found at many loci, most strongly at LIPC (15q22.1) and APOC-APOE complex (19q13.32), respectively. Associations at seven more loci previously known for effects on conventional plasma lipid measures reveal additional genetic influences on lipoprotein profiles and bring the total number of loci to 43. Thus, genome-wide associations identified novel loci involved with lipoprotein metabolism-including loci that affect the NMR-based measures of concentration or size of LDL, HDL, and VLDL particles-all characteristics of lipoprotein profiles that may impact disease risk but are not available by conventional assay. Show less