Human Genome-wide Association Studies Research Articles

Identification of biomarkers and molecular mechanisms implicated in genetic variations underlying Alzheimer's disease pathogenesis

We analyzed data from human genome-wide association studies (GWASs) to identify genetic variants and biological pathways linked to Alzheimer's disease (AD). Ten AD biomarkers (APOE, NECTIN2, APOC1, APOC1P1, TOMM40, RNU4-67P, KRAS, Y_RNA, THORLNC, LINC01956) were found across studies, including six central genetic variants (MAPT (rs242557-A), GRIN2B (rs74442473-G), APOE (rs438811-T), ANK3 (rs438811-T), BIN1 (rs744373-G), and BDNF (rs7481773-A)). ANK3 (rs438811-T) and GRIN2B (rs74442473-G) were essential hub biomarkers for amyloid plaques, while MAPT (rs242557-A) and BIN1 (rs744373-G) were crucial for neurofibrillary tangles (NFTs). Higher-risk AD biomarkers were associated with increased protein-lipid complex formation, while lower-risk AD biomarkers were correlated with improved synaptic function. Six essential miRNAs (hsa-miR-124–3p, 15a-5p, 16–5p, 204–5p, 520g-3p, 520h) and three transcription factors (ZMAT4, ZBED6, FOXG1) emerged as possible candidates to reveal the genetic differences that lead to amyloid plaques, NFTs, and ultimately AD. These findings serve as a basis for potential AD treatments and offer new avenues for therapeutic approaches to directly target the genetic variations and processes associated with the disease.

Agonistic anti-DCIR antibody inhibits ITAM-mediated inflammatory signaling and promotes immune resolution.

Dendritic cell inhibitory receptor (DCIR) is a C-type lectin receptor selectively expressed on myeloid cells, including monocytes, macrophage, dendritic cells, and neutrophils. Its role in immune regulation has been implicated in murine models and human genome-wide association studies (GWAS), suggesting defective DCIR function associates with increased susceptibility to autoimmune diseases such as rheumatoid arthritis, lupus and Sjogren's syndrome. However, little is known about the mechanisms underlying DCIR activation to dampen inflammation. Here, we developed anti-DCIR agonistic antibodies that promote phosphorylation on DCIR's immune receptor tyrosine-based inhibitory motifs (ITIM) and recruitment of SH2 containing protein tyrosine phosphatase-2 (SHP2) for reducing inflammation. We also explored the inflammation resolution by depleting DCIR+ cells with antibodies. Utilizing a human DCIR knock-in mouse model, we validated the anti-inflammatory properties of the agonistic anti-DCIR antibody in experimental peritonitis and colitis. These findings provide critical evidence for targeting DCIR to develop transformative therapies for inflammatory diseases.

Lack of evidence for GWAS signals of exfoliation glaucoma working via monogenic loss-of-function mutation in the nearest gene.

Exfoliation syndrome (XFS) is a systemic disease of elastin-rich tissues involving a deposition of fibrillar exfoliative material (XFM) in the anterior chamber of the eye, which can promote glaucoma. The purpose of this study was to create mice with CRISPR/Cas9-induced variations in candidate genes identified from human genome-wide association studies (GWAS) and screen them for indices of XFS. Variants predicted to be deleterious were sought in the Agpat1, Cacna1a, Loxl1, Pomp, Rbms3, Sema6a, and Tlcd5 genes of C57BL/6J mice using CRISPR/Cas9-based gene editing. Strains were phenotyped by slit-lamp, SD-OCT imaging, and fundus exams at 1-5 mos of age. Smaller cohorts of 12-mos-old mice were also studied. Deleterious variants were identified in six targets; Pomp was recalcitrant to targeting. Multiple alleles of some targets were isolated, yielding 12 strains. Across all genotypes and ages, 277 mice were assessed by 902 slit-lamp exams, 928 SD-OCT exams, and 358 fundus exams. Homozygosity for Agpat1 or Cacna1a mutations led to early lethality; homozygosity for Loxl1 mutations led to pelvic organ prolapse, preventing aging. Loxl1 homozygotes exhibited a conjunctival phenotype of potential relevance to XFS. Multiple other genotype-specific phenotypes were variously identified. XFM was not observed in any mice. This study did not detect XFM in any of the strains. This may have been due to species-specific differences, background dependence, or insufficient aging. Alternatively, it is possible that the current candidates, selected based on proximity to GWAS signals, are not effectors acting via monogenic loss-of-function mechanisms.

Amyloid β accelerates age-related proteome-wide protein insolubility.

Loss of proteostasis is a highly conserved feature of aging across model organisms and results in the accumulation of insoluble protein aggregates. Protein insolubility is also a unifying feature of major age-related neurodegenerative diseases, including Alzheimer's Disease (AD), in which hundreds of insoluble proteins associate with aggregated amyloid beta (Aβ) in senile plaques. Despite the connection between aging and AD risk, therapeutic approaches to date have overlooked aging-driven generalized protein insolubility as a contributing factor. However, proteins that become insoluble during aging in model organisms are capable of accelerating Aβ aggregation in vitro and lifespan in vivo. Here, using an unbiased proteomics approach, we questioned the relationship between Aβ and age-related protein insolubility. Specifically, we uncovered that Aβ expression drives proteome-wide protein insolubility in C. elegans, even in young animals, and this insoluble proteome is highly similar to the insoluble proteome driven by normal aging, this vulnerable sub-proteome we term the core insoluble proteome (CIP). We show that the CIP is enriched with proteins that modify Aβ toxicity in vivo, suggesting the possibility of a vicious feedforward cycle in the context of AD. Importantly, using human genome-wide association studies (GWAS), we show that the CIP is replete with biological processes implicated not only in neurodegenerative diseases but also across a broad array of chronic, age-related diseases (CARDs). This provides suggestive evidence that age-related loss of proteostasis could play a role in general CARD risk. Finally, we show that the geroprotective, gut-derived metabolite, Urolithin A, relieves Aβ toxicity, supporting its use in clinical trials for dementia and age-related diseases.

Reduction of APOE accounts for neurobehavioral deficits in fetal alcohol spectrum disorders.

A hallmark of fetal alcohol spectrum disorders (FASD) is neurobehavioral deficits that still do not have effective treatment. Here, we present that reduction of Apolipoprotein E (APOE) is critically involved in neurobehavioral deficits in FASD. We show that prenatal alcohol exposure (PAE) changes chromatin accessibility of Apoe locus, and causes reduction of APOE levels in both the brain and peripheral blood in postnatal mice. Of note, postnatal administration of an APOE receptor agonist (APOE-RA) mitigates motor learning deficits and anxiety in those mice. Several molecular and electrophysiological properties essential for learning, which are altered by PAE, are restored by APOE-RA. Our human genome-wide association study further reveals that the interaction of PAE and a single nucleotide polymorphism in the APOE enhancer which chromatin is closed by PAE in mice is associated with lower scores in the delayed matching-to-sample task in children. APOE in the plasma is also reduced in PAE children, and the reduced level is associated with their lower cognitive performance. These findings suggest that controlling the APOE level can serve as an effective treatment for neurobehavioral deficits in FASD.

P124 Elucidating the role of risk variants in the BIM gene enhancer in T-Cell dysregulation and rheumatoid arthritis pathogenesis

Abstract Background/Aims The balance governing T-cell responses is pivotal for immune homeostasis. Disruptions to this balance may lead to aberrant behaviours, potentially resulting in autoimmunity. Rheumatoid Arthritis (RA), an autoimmune disorder primarily manifesting in the joints, is characterised by a defective immune response, resulting in chronic inflammation. Enhanced T-cell apoptosis resistance in this context can contribute to sustained activation and abnormal proliferation in a pro-inflammatory environment. Understanding these disease mechanisms could lead to improvements in clinical care. Genome-wide association studies (GWAS) have identified hundreds of single nucleotide polymorphisms (SNPs) associated with RA, preferentially enriched in non-coding, T-cell-specific putative enhancers. Enhancers are DNA elements that can regulate the transcription of specific genes, often from many kilobases away. Our previous work identified BCL2L11, or BIM, as the probable causal gene within the RA locus rs13396472, despite the GWAS variants' intronic positioning relative to ACOXL. The RA risk variant disrupts a distant enhancer regulating the BIM gene, pivotal for apoptosis. Given the role of apoptosis in immune tolerance, disruption of this pathway driven by RA risk variants could pave the way for autoimmune responses. Our aim was to investigate the effect of this enhancer on T-cell proliferation and function in an in vivo model. Methods We mapped the human RA-locus to the mouse genome and utilised CRISPR-Cas9 to delete this putative RA enhancer in mice. Two sgRNAs, combined with Cas9, were microinjected into C57BL6/J zygote pronuclei, followed by embryo transfer to pseudo-pregnant mice. PCR and Sanger sequencing confirmed enhancer deletion in pups, establishing an RA-enhancer deletion colony. To gauge cell proportion differences, murine splenocytes or thymocytes were stained with specific fluorochrome markers and analysed via a flow cytometer. For the antigen-induced arthritis (AIA) model, wildtype and RA-enhancer deletion mice were subjected to arthritis induction via intra-articular injection of methylated bovine serum albumin, following prior immunisation, and monitored for arthritis signs. Results We observed an increased CD4+CD8+ double positive (DP) T-cell population in aged mice homozygous and heterozygous for the RA-enhancer deletion, but not in younger and wild-type controls. Knockout mice exhibited visibly enlarged spleens in response to the immunisation step prior to AIA induction. Further work will establish the mechanisms underlying this and disease progression after arthritis induction. These observations underscore a potential impact this enhancer might have on immune responses and survival upon immune system alteration. Conclusion Our research indicates that perturbing mouse non-coding sequences syntenic to human GWAS loci can provide a physiological model to elucidate novel mechanisms underlying RA, potentially suggesting therapeutic targets. Future research will explore the mechanism by which the enhancer deletion leads to an increase in DP T-cells and the altered phenotype, and the implications of such a change for T- and B-cell biology and RA. Disclosure S. Rossi: None. E. Richards: None. J. Ding: None. C. Shi: None. M. Menon: None. J. Gibbs: None. R. Grencis: None. A. Adamson: None. G. Orozco: None.

Genetic loci regulate Sarbecovirus pathogenesis: A comparison across mice and humans

Coronavirus (CoV) cause considerable morbidity and mortality in humans and other mammals, as evidenced by the emergence of Severe Acute Respiratory CoV (SARS-CoV) in 2003, Middle East Respiratory CoV (MERS-CoV) in 2012, and SARS-CoV-2 in 2019. Although poorly characterized, natural genetic variation in human and other mammals modulate virus pathogenesis, as reflected by the spectrum of clinical outcomes ranging from asymptomatic infections to lethal disease. Using multiple human epidemic and zoonotic Sarbecoviruses, coupled with murine Collaborative Cross genetic reference populations, we identify several dozen quantitative trait loci that regulate SARS-like group-2B CoV pathogenesis and replication. Under a Chr4 QTL, we deleted a candidate interferon stimulated gene, Trim14 which resulted in enhanced SARS-CoV titers and clinical disease, suggesting an antiviral role during infection. Importantly, about 60 % of the murine QTL encode susceptibility genes identified as priority candidates from human genome-wide association studies (GWAS) studies after SARS-CoV-2 infection, suggesting that similar selective forces have targeted analogous genes and pathways to regulate Sarbecovirus disease across diverse mammalian hosts. These studies provide an experimental platform in rodents to investigate the molecular-genetic mechanisms by which potential cross mammalian susceptibility loci and genes regulate type-specific and cross-SARS-like group 2B CoV replication, immunity, and pathogenesis in rodent models. Our study also provides a paradigm for identifying susceptibility loci for other highly heterogeneous and virulent viruses that sporadically emerge from zoonotic reservoirs to plague human and animal populations.

Monkey multi-organ cell atlas exposed to estrogen

Abstract Awareness of estrogen’s effects on health is broadening rapidly. The effects of long-term high levels of estrogen on the body involve multiple organs. Here, we used both single-cell chromatin accessibility and RNA sequencing data to analyze the potential effect of estrogen on major organs. The integrated cell map enabled in-depth dissection and comparison of molecular dynamics, cell-type compositions, and cellular heterogeneity across multiple tissues and organs under estrogen stimulation. We also inferred pseudotime cell trajectories and cell–cell communications to uncover key molecular signatures underlying their cellular processes in major organs in response to estrogen. For example, estrogen could induce the differentiation of IFIT3+ neutrophils into S100A9+ neutrophils involved in the function of endosome-to-lysosome transport and the multivesicular body sorting pathway in liver tissues. Furthermore, through integration with human genome-wide association study data, we further identified a subset of risk genes during disease development that were induced by estrogen, such as AKT1 (related to endometrial cancer), CCND1 (related to breast cancer), HSPH1 (related to colorectal cancer), and COVID-19 and asthma-related risk genes. Our work uncovers the impact of estrogen on the major organs, constitutes a useful resource, and reveals the contribution and mechanism of estrogen to related diseases.

BridGE: a pathway-based analysis tool for detecting genetic interactions from GWAS.

Genetic interactions have the potential to modulate phenotypes, including human disease. In principle, genome-wide association studies (GWAS) provide a platform for detecting genetic interactions; however, traditional methods for identifying them, which tend to focus on testing individual variant pairs, lack statistical power. In this protocol, we describe a novel computational approach, called Bridging Gene sets with Epistasis (BridGE), for discovering genetic interactions between biological pathways from GWAS data. We present a Python-based implementation of BridGE along with instructions for its application to a typical human GWAS cohort. The major stages include initial data processing and quality control, construction of a variant-level genetic interaction network, measurement of pathway-level genetic interactions, evaluation of statistical significance using sample permutations and generation of results in a standardized output format. The BridGE software pipeline includes options for running the analysis on multiple cores and multiple nodes for users who have access to computing clusters or a cloud computing environment. In a cluster computing environment with 10 nodes and 100 GB of memory per node, the method can be run in less than 24 h for typical human GWAS cohorts. Using BridGE requires knowledge of running Python programs and basic shell script programming experience.

Genome wide association study in Swedish Labrador retrievers identifies genetic loci associated with hip dysplasia and body weight

Genome wide association studies (GWAS) have been utilized to identify genetic risk loci associated with both simple and complex inherited disorders. Here, we performed a GWAS in Labrador retrievers to identify genetic loci associated with hip dysplasia and body weight. Hip dysplasia scores were available for 209 genotyped dogs. We identified a significantly associated locus for hip dysplasia on chromosome 24, with three equally associated SNPs (p = 4.3 × 10–7) in complete linkage disequilibrium located within NDRG3, a gene which in humans has been shown to be differentially expressed in osteoarthritic joint cartilage. Body weight, available for 85 female dogs, was used as phenotype for a second analysis. We identified two significantly associated loci on chromosome 10 (p = 4.5 × 10–7) and chromosome 31 (p = 2.5 × 10–6). The most associated SNPs within these loci were located within the introns of the PRKCE and CADM2 genes, respectively. PRKCE has been shown to play a role in regulation of adipogenesis whilst CADM2 has been associated with body weight in multiple human GWAS. In summary, we identified credible candidate loci explaining part of the genetic inheritance for hip dysplasia and body weight in Labrador retrievers with strong candidate genes in each locus previously implicated in the phenotypes investigated.

Divergent role of Mitochondrial Amidoxime Reducing Component 1 (MARC1) in human and mouse.

Recent human genome-wide association studies have identified common missense variants in MARC1, p.Ala165Thr and p.Met187Lys, associated with lower hepatic fat, reduction in liver enzymes and protection from most causes of cirrhosis. Using an exome-wide association study we recapitulated earlier MARC1 p.Ala165Thr and p.Met187Lys findings in 540,000 individuals from five ancestry groups. We also discovered novel rare putative loss of function variants in MARC1 with a phenotype similar to MARC1 p.Ala165Thr/p.Met187Lys variants. In vitro studies of recombinant human MARC1 protein revealed Ala165Thr substitution causes protein instability and aberrant localization in hepatic cells, suggesting MARC1 inhibition or deletion may lead to hepatoprotection. Following this hypothesis, we generated Marc1 knockout mice and evaluated the effect of Marc1 deletion on liver phenotype. Unexpectedly, our study found that whole-body Marc1 deficiency in mouse is not protective against hepatic triglyceride accumulation, liver inflammation or fibrosis. In attempts to explain the lack of the observed phenotype, we discovered that Marc1 plays only a minor role in mouse liver while its paralogue Marc2 is the main Marc family enzyme in mice. Our findings highlight the major difference in MARC1 physiological function between human and mouse.

Targeting postsynaptic glutamate receptor scaffolding proteins PSD-95 and PICK1 for obesity treatment.

Human genome-wide association studies (GWAS) suggest a functional role for central glutamate receptor signaling and plasticity in body weight regulation. Here, we use UK Biobank GWAS summary statistics of body mass index (BMI) and body fat percentage (BF%) to identify genes encoding proteins known to interact with postsynaptic α-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid (AMPA) and N-methyl-d-aspartate (NMDA) receptors. Loci in/near discs large homolog 4 (DLG4) and protein interacting with C kinase 1 (PICK1) reached genome-wide significance (P < 5 × 10-8) for BF% and/or BMI. To further evaluate the functional role of postsynaptic density protein-95 (PSD-95; gene name: DLG4) and PICK1 in energy homeostasis, we used dimeric PSD-95/disc large/ZO-1 (PDZ) domain-targeting peptides of PSD-95 and PICK1 to demonstrate that pharmacological inhibition of PSD-95 and PICK1 induces prolonged weight-lowering effects in obese mice. Collectively, these data demonstrate that the glutamate receptor scaffolding proteins, PICK1 and PSD-95, are genetically linked to obesity and that pharmacological targeting of their PDZ domains represents a promising therapeutic avenue for sustained weight loss.

Inactivation of phosphodiesterase-4B gene in rat nucleus accumbens shell by CRISPR/Cas9 or positive allosteric modulation of the protein affects the motivation to chronically self-administer nicotine

Large scale human genome wide association studies (GWAS) have identified a growing pool of genes associated with cigarette smoking. One of the most prominent, phosphodiesterase-4B (PDE4B), has been associated with multiple smoking phenotypes. Although PDE4B modulates the half-life of neuronal cAMP, its precise role in smoking behaviors is unknown. To address this knowledge gap, we used a reverse translational approach. We inactivated PDE4B in bilateral medial nucleus accumbens shell (NAcs) neurons by injecting AAV containing a specific gRNA in female transgenic Cas9+ Long Evans rats. These rats then were given 23-h chronic access to nicotine intravenous self-administration (IVSA) under a schedule of increasing fixed ratios (FR). With the increased effort required at FR7, nicotine SA (i.e. active presses and drug infusions) declined significantly in controls, whereas it was maintained in the mutagenized group. A progressive ratio (PR) study also showed significantly greater cumulative nicotine infusions in the PDE4B-edited group. Hence, we hypothesized that enhanced PDE4B protein activity would reduce nicotine IVSA. A positive allosteric modulator, 2-(3-(4-chloro-3-fluorophenyl)-5-ethyl-1H-1,2,4-triazol-1-yl)-N-(3,5-dichlorobenzyl)acetamide (MR-L2), was microinfused into NAcs bilaterally at FR3 or FR5; in both cohorts, MR-L2 acutely reduced nicotine IVSA. In summary, these studies show that the activity of PDE4B regulates the capacity of NAcs to maintain nicotine IVSA in face of the cost of increasing work. This finding and the results of the PR study indicate that PDE4B affects the motivation to obtain nicotine. These reverse translational studies in rats provide insight into the motivational effects of NAcs PDE4B that advance our understanding of the smoking behaviors mapped in human GWAS.

Genetic effects of sequence-conserved enhancer-like elements on human complex traits.

The vast majority of findings from human genome-wide association studies (GWAS) map to non-coding sequences, complicating their mechanistic interpretations and clinical translations. Non-coding sequences that are evolutionarily conserved and biochemically active could offer clues to the mechanisms underpinning GWAS discoveries. However, genetic effects of such sequences have not been systematically examined across a wide range of human tissues and traits, hampering progress to fully understand regulatory causes of human complex traits. Here we develop a simple yet effective strategy to identify functional elements exhibiting high levels of human-mouse sequence conservation and enhancer-like biochemical activity, which scales well to 313 epigenomic datasets across 106 human tissues and cell types. Combined with 468 GWAS of European (EUR) and East Asian (EAS) ancestries, these elements show tissue-specific enrichments of heritability and causal variants for many traits, which are significantly stronger than enrichments based on enhancers without sequence conservation. These elements also help prioritize candidate genes that are functionally relevant to body mass index (BMI) and schizophrenia but were not reported in previous GWAS with large sample sizes. Our findings provide a comprehensive assessment of how sequence-conserved enhancer-like elements affect complex traits in diverse tissues and demonstrate a generalizable strategy of integrating evolutionary and biochemical data to elucidate human disease genetics.

Central Channelopathies in Obesity.

As obesity has raised heightening awareness, researchers have attempted to identify potential targets that can be treated for therapeutic intervention. Focusing on the central nervous system (CNS), the key organ in maintaining energy balance, a plethora of ion channels that are expressed in the CNS have been inspected and determined through manipulation in different hypothalamic neural subpopulations for their roles in fine-tuning neuronal activity on energy state alterations, possibly acting as metabolic sensors. However, a remaining gap persists between human clinical investigations and mouse studies. Despite having delineated the pathways and mechanisms of how the mouse study-identified ion channels modulate energy homeostasis, only a few targets overlap with the obesity-related risk genes extracted from human genome-wide association studies. Here, we present the most recently discovered CNS-specific metabolism-correlated ion channels using reverse and forward genetics approaches in mice and humans, respectively, in the hope of illuminating the prospects for future therapeutic development.

An inverse causal relationship between serum 25-hydroxyvitamin D levels and pulmonary hypertension: A two-sample Mendelianrandomization study.

Observational studies have confirmed that 25-hydroxyvitamin D (25(OH)D) is associated with pulmonary hypertension (PH), but the causal association between each other is unclear. Therefore, Mendelian randomization (MR) method was performed to validate the causal association between PH and serum 25(OH)D levels. The summary data for 25(OH)D and PH were from the National Human Genome Research Institute-European Bioinformatics Institute. Catalog of human genome-wide association studies and FinnGen biobank consortium. MR analysis was utilized to explore the potential causal association between PH and 25(OH)D. To evaluate this association, inverse variance weighting was considered as the primary method. Cochran's Q test, MR-Egger intercept test, and "leave-one-out" sensitivity analyses were utilized to control the pleiotropy and heterogeneity in the study. Two-sample MR analysis revealed an inverse causal relationship between 25(OH)D and PH (odds ratio: 0.376, 95% confidence interval: 0.162-0.876, p = 2.334 × 10-2). There was no significant heterogeneity and pleiotropy. The present study confirmed the inverse causal relationship between 25(OH)D and PH. This pathway may provide another treatment pathway in PH. Further studies to elucidate this pathway is indicated.

Nicotinic acetylcholine receptor signaling maintains epithelial barrier integrity.

Disruption of epithelial barriers is a common disease manifestation in chronic degenerative diseases of the airways, lung, and intestine. Extensive human genetic studies have identified risk loci in such diseases, including in chronic obstructive pulmonary disease (COPD) and inflammatory bowel diseases. The genes associated with these loci have not fully been determined, and functional characterization of such genes requires extensive studies in model organisms. Here, we report the results of a screen in Drosophila melanogaster that allowed for rapid identification, validation, and prioritization of COPD risk genes that were selected based on risk loci identified in human genome-wide association studies (GWAS). Using intestinal barrier dysfunction in flies as a readout, our results validate the impact of candidate gene perturbations on epithelial barrier function in 56% of the cases, resulting in a prioritized target gene list. We further report the functional characterization in flies of one family of these genes, encoding for nicotinic acetylcholine receptor (nAchR) subunits. We find that nAchR signaling in enterocytes of the fly gut promotes epithelial barrier function and epithelial homeostasis by regulating the production of the peritrophic matrix. Our findings identify COPD-associated genes critical for epithelial barrier maintenance, and provide insight into the role of epithelial nAchR signaling for homeostasis.

PTPN2 Regulates Metabolic Flux to Affect β-Cell Susceptibility to Inflammatory Stress.

Protein tyrosine phosphatase N2 (Ptpn2) is a type 1 diabetes (T1D) candidate gene identified from human genome-wide association studies. PTPN2 is highly expressed in human and murine islets and becomes elevated upon inflammation and models of T1D, suggesting that PTPN2 may be important for beta cell survival in the context of T1D. To test whether PTPN2 contributed to beta cell dysfunction in an inflammatory environment, we generated a beta cell-specific deletion of Ptpn2 in mice (Ptpn2 βKO). While unstressed animals exhibit normal metabolic profiles, low and high dose streptozotocin (STZ) treated Ptpn2 βKO mice displayed hyperglycemia and accelerated death, respectively. Furthermore, cytokine treated Ptpn2 KO islets resulted in impaired glucose-stimulated insulin secretion, mitochondrial defects and reduced glucose-induced metabolic flux, suggesting beta cells lacking Ptpn2 are more susceptible to inflammatory stress associated with T1D due to maladaptive metabolic fitness. Consistent with the phenotype, proteomic analysis identified an important metabolic enzyme ATP-citrate lyase (ACLY) as a novel PTPN2 substrate.

Identification of genetic drivers of plasma lipoprotein size in the Diversity Outbred mouse population

Despite great progress in understanding lipoprotein physiology, there is still much to be learned about the genetic drivers of lipoprotein abundance, composition, and function. We used ion mobility spectrometry to survey 16 plasma lipoprotein subfractions in 500 Diversity Outbred mice maintained on a Western-style diet. We identified 21 quantitative trait loci (QTL) affecting lipoprotein abundance. To refine the QTL and link them to disease risk in humans, we asked if the human homologs of genes located at each QTL were associated with lipid traits in human genome-wide association studies. Integration of mouse QTL with human genome-wide association studies yielded candidate gene drivers for 18 of the 21 QTL. This approach enabled us to nominate the gene encoding the neutral ceramidase, Asah2, as a novel candidate driver at a QTL on chromosome 19 for large HDL particles (HDL-2b). To experimentally validate Asah2, we surveyed lipoproteins in Asah2−/− mice. Compared to wild-type mice, female Asah2−/− mice showed an increase in several lipoproteins, including HDL. Our results provide insights into the genetic regulation of circulating lipoproteins, as well as mechanisms by which lipoprotein subfractions may affect cardiovascular disease risk in humans.

Missense variants in SORT1 are associated with LDL-C in an Amish population

Common noncoding variants at the human 1p13.3 locus associated with SORT1 expression are among those most strongly associated with low-density lipoprotein cholesterol (LDL-C) in human genome-wide association studies. However, validation studies in mice and cell lines have produced variable results regarding the directionality of the effect of SORT1 on LDL-C. This, together with the fact that the 1p13.3 variants are associated with expression of several genes, has raised the question of whether SORT1 is the causal gene at this locus. Using whole exome sequencing in members of an Amish population, we identified coding variants in SORT1 that are associated with increased (rs141749679, K302E) and decreased (rs149456022, Q225H) LDL-C. Further, analysis of plasma lipoprotein particle subclasses by ion mobility in a subset of rs141749679 (K302E) carriers revealed higher levels of large LDL particles compared to noncarriers. In contrast to the effect of these variants in the Amish, the sortilin K302E mutation introduced into a C57BL6/J mouse via CRISPR/Cas9 had decreased non-high-density lipoprotein cholesterol, and the sortilin Q225H mutation did not alter cholesterol levels in mice. This is indicative of different effects of these mutations on cholesterol metabolism in the two species. To our knowledge, this is the first evidence that naturally occurring coding variants in SORT1 are associated with LDL-C, thus supporting SORT1 as the gene responsible for the association of the 1p13.3 locus with LDL-C.