The AIM2 inflammasome exacerbates atherosclerosis in clonal haematopoiesis.

Clonal haematopoiesis, which is highly prevalent in older individuals, arises from somatic mutations that endow a proliferative advantage to haematopoietic cells. Clonal haematopoiesis increases the risk of myocardial infarction and stroke independently of traditional risk factors1. Among the common genetic variants that give rise to clonal haematopoiesis, the JAK2V617F (JAK2VF) mutation, which increases JAK-STAT signalling, occurs at a younger age and imparts the strongest risk of premature coronary heart disease1,2. Here we show increased proliferation of macrophages and prominent formation of necrotic cores in atherosclerotic lesions in mice that express Jak2VF selectively in macrophages, and in chimeric mice that model clonal haematopoiesis. Deletion of the essential inflammasome components caspase 1 and 11, or of the pyroptosis executioner gasdermin D, reversed these adverse changes. Jak2VF lesions showed increased expression of AIM2, oxidative DNA damage and DNA replication stress, and Aim2 deficiency reduced atherosclerosis. Single-cell RNA sequencing analysis of Jak2VF lesions revealed a landscape that was enriched for inflammatory myeloid cells, which were suppressed by deletion of Gsdmd. Inhibition of the inflammasome product interleukin-1ß reduced macrophage proliferation and necrotic formation while increasing the thickness of fibrous caps, indicating that it stabilized plaques. Our findings suggest that increased proliferation and glycolytic metabolism in Jak2VF macrophages lead to DNA replication stress and activation of the AIM2 inflammasome, thereby aggravating atherosclerosis. Precise application of therapies that target interleukin-1ß or specific inflammasomes according to clonal haematopoiesis status could substantially reduce cardiovascular risk.

Atherosclerosis Is a Smooth Muscle Cell-Driven Tumor-Like Disease.

BACKGROUND: Atherosclerosis, a leading cause of cardiovascular disease, involves the pathological activation of various cell types, including immunocytes (eg, macrophages and T cells), smooth muscle cells (SMCs), and endothelial cells. Accumulating evidence suggests that transition of SMCs to other cell types, known as phenotypic switching, plays a central role in atherosclerosis development and complications. However, the characteristics of SMC-derived cells and the underlying mechanisms of SMC transition in disease pathogenesis remain poorly understood. Our objective is to characterize tumor cell-like behaviors of SMC-derived cells in atherosclerosis, with the ultimate goal of developing interventions targeting SMC transition for the prevention and treatment of atherosclerosis. METHODS: We used SMC lineage tracing mice and human tissues and applied a range of methods, including molecular, cellular, histological, computational, human genetics, and pharmacological approaches, to investigate the features of SMC-derived cells in atherosclerosis. RESULTS: SMC-derived cells in mouse and human atherosclerosis exhibit multiple tumor cell-like characteristics, including genomic instability, evasion of senescence, hyperproliferation, resistance to cell death, invasiveness, and activation of comprehensive cancer-associated gene regulatory networks. Specific expression of the oncogenic mutant KrasG12D in SMCs accelerates phenotypic switching and exacerbates atherosclerosis. Furthermore, we provide proof of concept that niraparib, an anticancer drug targeting DNA damage repair, attenuates atherosclerosis progression and induces regression of lesions in advanced disease in mouse models. CONCLUSIONS: Our findings demonstrate that atherosclerosis is an SMC-driven tumor-like disease, advancing our understanding of its pathogenesis and opening prospects for innovative precision molecular strategies aimed at preventing and treating atherosclerotic cardiovascular disease.

High-Dimensional Single-Cell Multimodal Landscape of Human Carotid Atherosclerosis.

BACKGROUND: Atherosclerotic plaques are complex tissues composed of a heterogeneous mixture of cells. However, our understanding of the comprehensive transcriptional and phenotypic landscape of the cells within these lesions is limited. METHODS: To characterize the landscape of human carotid atherosclerosis in greater detail, we combined cellular indexing of transcriptomes and epitopes by sequencing and single-cell RNA sequencing to classify all cell types within lesions (n=21; 13 symptomatic) to achieve a comprehensive multimodal understanding of the cellular identities of atherosclerosis and their association with clinical pathophysiology. RESULTS: We identified 25 cell populations, each with a unique multiomic signature, including macrophages, T cells, NK (natural killer) cells, mast cells, B cells, plasma cells, neutrophils, dendritic cells, endothelial cells, fibroblasts, and smooth muscle cells (SMCs). Among the macrophages, we identified 2 proinflammatory subsets enriched in IL-1B (interleukin-1B) or C1Q expression, 2 TREM2-positive foam cells (1 expressing inflammatory genes), and subpopulations with a proliferative gene signature and SMC-specific gene signature with fibrotic pathways upregulated. Further characterization revealed various subsets of SMCs and fibroblasts, including SMC-derived foam cells. These foamy SMCs were localized in the deep intima of coronary atherosclerotic lesions. Utilizing cellular indexing of transcriptomes and epitopes by sequencing data, we developed a flow cytometry panel, using cell surface proteins CD29, CD142, and CD90, to isolate SMC-derived cells from lesions. Lastly, we observed reduced proportions of efferocytotic macrophages, classically activated endothelial cells, and contractile and modulated SMC-derived cells, while inflammatory SMCs were enriched in plaques of clinically symptomatic versus asymptomatic patients. CONCLUSIONS: Our multimodal atlas of cell populations within atherosclerosis provides novel insights into the diversity, phenotype, location, isolation, and clinical relevance of the unique cellular composition of human carotid atherosclerosis. These findings facilitate both the mapping of cardiovascular disease susceptibility loci to specific cell types and the identification of novel molecular and cellular therapeutic targets for the treatment of the disease.

Human Macrophage Long Intergenic Noncoding RNA, SIMALR, Suppresses Inflammatory Macrophage Apoptosis via NTN1 (Netrin-1).

BACKGROUND: Long noncoding RNAs (lncRNAs) have emerged as novel regulators of macrophage biology and inflammatory cardiovascular diseases. However, studies focused on lncRNAs in human macrophage subtypes, particularly human lncRNAs that are not conserved in rodents, are limited. METHODS: Through RNA-sequencing of human monocyte-derived macrophages, we identified suppressor of inflammatory macrophage apoptosis lncRNA (SIMALR). Lipopolysaccharide/IFNγ (interferon γ) stimulated human macrophages were treated with SIMALR antisense oligonucleotides and subjected to RNA-sequencing to investigate the function of SIMALR. Western blots, luciferase assay, and RNA immunoprecipitation were performed to validate function and potential mechanism of SIMALR. RNAscope was performed to identify SIMALR expression in human carotid atherosclerotic plaques. RESULTS: RNA-sequencing of human monocyte-derived macrophages identified SIMALR, a human macrophage-specific long intergenic noncoding RNA that is highly induced in lipopolysaccharide/IFNγ-stimulated macrophages. SIMALR knockdown in lipopolysaccharide/IFNγ stimulated THP1 human macrophages induced apoptosis of inflammatory macrophages, as shown by increased protein expression of cleaved PARP (poly[ADP-ribose] polymerase), caspase 9, caspase 3, and Annexin V+. RNA-sequencing of control versus SIMALR knockdown in lipopolysaccharide/IFNγ-stimulated macrophages showed Netrin-1 (NTN1) to be significantly decreased upon SIMALR knockdown. We confirmed that NTN1 knockdown in lipopolysaccharide/IFNγ-stimulated macrophages induced apoptosis. The SIMALR knockdown-induced apoptotic phenotype was rescued by adding recombinant NTN1. NTN1 promoter-luciferase reporter activity was increased in HEK293T (human embryonic kidney 293) cells treated with lentiviral overexpression of SIMALR. NTN1 promoter activity is known to require HIF1α (hypoxia-inducible factor 1 subunit alpha), and our studies suggest that SIMALR may interact with HIF1α to regulate NTN1 transcription, thereby regulating macrophages apoptosis. SIMALR was found to be expressed in macrophages in human carotid atherosclerotic plaques of symptomatic patients. CONCLUSIONS: SIMALR is a nonconserved, human macrophage lncRNA expressed in atherosclerosis that suppresses macrophage apoptosis. SIMALR partners with HIF1α (hypoxia-inducible factor 1 subunit alpha) to regulate NTN1, which is a known macrophage survival factor. This work illustrates the importance of interrogating the functions of human lncRNAs and exploring their translational and therapeutic potential in human atherosclerosis.

ASEP: Gene-based detection of allele-specific expression across individuals in a population by RNA sequencing.

Allele-specific expression (ASE) analysis, which quantifies the relative expression of two alleles in a diploid individual, is a powerful tool for identifying cis-regulated gene expression variations that underlie phenotypic differences among individuals. Existing methods for gene-level ASE detection analyze one individual at a time, therefore failing to account for shared information across individuals. Failure to accommodate such shared information not only reduces power, but also makes it difficult to interpret results across individuals. However, when only RNA sequencing (RNA-seq) data are available, ASE detection across individuals is challenging because the data often include individuals that are either heterozygous or homozygous for the unobserved cis-regulatory SNP, leading to sample heterogeneity as only those heterozygous individuals are informative for ASE, whereas those homozygous individuals have balanced expression. To simultaneously model multi-individual information and account for such heterogeneity, we developed ASEP, a mixture model with subject-specific random effect to account for multi-SNP correlations within the same gene. ASEP only requires RNA-seq data, and is able to detect gene-level ASE under one condition and differential ASE between two conditions (e.g., pre- versus post-treatment). Extensive simulations demonstrated the convincing performance of ASEP under a wide range of scenarios. We applied ASEP to a human kidney RNA-seq dataset, identified ASE genes and validated our results with two published eQTL studies. We further applied ASEP to a human macrophage RNA-seq dataset, identified genes showing evidence of differential ASE between M0 and M1 macrophages, and confirmed our findings by results from cardiometabolic trait-relevant genome-wide association studies. To the best of our knowledge, ASEP is the first method for gene-level ASE detection at the population level that only requires the use of RNA-seq data. With the growing adoption of RNA-seq, we believe ASEP will be well-suited for various ASE studies for human diseases.

Cis-regulated expression of non-conserved lincRNAs associates with cardiometabolic related traits.

Many complex disease risk loci map to intergenic regions containing long intergenic noncoding RNAs (lincRNAs). The majority of these is not conserved outside humans, raising the question whether genetically regulated expression of non-conserved and conserved lincRNAs has similar rates of association with complex traits. Here we leveraged data from the Genotype-Tissue Expression (GTEx) project and multiple public genome-wide association study (GWAS) resources. Using an established transcriptome-wide association study (TWAS) tool, FUSION, we interrogated the associations between cis-regulated expression of lincRNAs and multiple cardiometabolic traits. We found that cis-regulated expression of non-conserved lincRNAs had a strikingly similar trend of association with complex cardiometabolic traits as conserved lincRNAs. This finding challenges the conventional notion of conservation that has led to prioritization of conserved loci for functional studies and calls attention to the need to develop comprehensive strategies to study the large number of non-conserved human lincRNAs that may contribute to human disease.

Nonconserved Long Intergenic Noncoding RNAs Associate With Complex Cardiometabolic Disease Traits.

OBJECTIVE: Transcriptome profiling of human tissues has revealed thousands of long intergenic noncoding RNAs (lincRNAs) at loci identified through large-scale genome-wide studies for complex cardiometabolic traits. This raises the question of whether genetic variation at nonconserved lincRNAs has any systematic association with complex disease, and if so, how different this pattern is from conserved lincRNAs. We evaluated whether the associations between nonconserved lincRNAs and 8 complex cardiometabolic traits resemble or differ from the pattern of association for conserved lincRNAs. Approach and Results: Our investigation of over 7000 lincRNA annotations from GENCODE Release 33-GRCh38.p13 for complex trait genetic associations leveraged several large, established meta-analyses genome-wide association study summary data resources, including GIANT (Genetic Investigation of Anthropometric Traits), UK Biobank, GLGC (Global Lipids Genetics Consortium), Cardiogram (Coronary Artery Disease Genome Wide Replication and Meta-Analysis), and DIAGRAM (Diabetes Genetics Replication and Meta-Analysis)/DIAMANTE (Diabetes Meta-Analysis of Trans-Ethnic Association Studies). These analyses revealed that (1) nonconserved lincRNAs associate with a range of cardiometabolic traits at a rate that is generally consistent with conserved lincRNAs; (2) these findings persist across different definitions of conservation; and (3) overall across all cardiometabolic traits, approximately one-third of genome-wide association study-associated lincRNAs are nonconserved, and this increases to about two-thirds using a more stringent definition of conservation. CONCLUSIONS: These findings suggest that the traditional notion of conservation driving prioritization for functional and translational follow-up of complex cardiometabolic genomic discoveries may need to be revised in the context of the abundance of nonconserved long noncoding RNAs in the human genome and their apparent predilection to associate with complex cardiometabolic traits.

Single-Cell Genomics Reveals a Novel Cell State During Smooth Muscle Cell Phenotypic Switching and Potential Therapeutic Targets for Atherosclerosis in Mouse and Human.

BACKGROUND: Smooth muscle cells (SMCs) play significant roles in atherosclerosis via phenotypic switching, a pathological process in which SMC dedifferentiation, migration, and transdifferentiation into other cell types. Yet how SMCs contribute to the pathophysiology of atherosclerosis remains elusive. METHODS: To reveal the trajectories of SMC transdifferentiation during atherosclerosis and to identify molecular targets for disease therapy, we combined SMC fate mapping and single-cell RNA sequencing of both mouse and human atherosclerotic plaques. We also performed cell biology experiments on isolated SMC-derived cells, conducted integrative human genomics, and used pharmacological studies targeting SMC-derived cells both in vivo and in vitro. RESULTS: We found that SMCs transitioned to an intermediate cell state during atherosclerosis, which was also found in human atherosclerotic plaques of carotid and coronary arteries. SMC-derived intermediate cells, termed "SEM" cells (stem cell, endothelial cell, monocyte), were multipotent and could differentiate into macrophage-like and fibrochondrocyte-like cells, as well as return toward the SMC phenotype. Retinoic acid (RA) signaling was identified as a regulator of SMC to SEM cell transition, and RA signaling was dysregulated in symptomatic human atherosclerosis. Human genomics revealed enrichment of genome-wide association study signals for coronary artery disease in RA signaling target gene loci and correlation between coronary artery disease risk alleles and repressed expression of these genes. Activation of RA signaling by all-trans RA, an anticancer drug for acute promyelocytic leukemia, blocked SMC transition to SEM cells, reduced atherosclerotic burden, and promoted fibrous cap stability. CONCLUSIONS: Integration of cell-specific fate mapping, single-cell genomics, and human genetics adds novel insights into the complexity of SMC biology and reveals regulatory pathways for therapeutic targeting of SMC transitions in atherosclerotic cardiovascular disease.

Genome-wide interrogation reveals hundreds of long intergenic noncoding RNAs that associate with cardiometabolic traits.

Long intergenic noncoding RNAs (lincRNAs) play important roles in disease, but the vast majority of these transcripts remain uncharacterized. We defined a set of 54 944 human lincRNAs by drawing on four publicly available lincRNA datasets, and annotated ∼2.5 million single nucleotide polymorphisms (SNPs) from each of 15 cardiometabolic genome-wide association study datasets into these lincRNAs. We identified hundreds of lincRNAs with at least one trait-associated SNP: 898 SNPs in 343 unique lincRNAs at 5% false discovery rate, and 469 SNPs in 146 unique lincRNAs meeting Bonferroni-corrected P < 0.05. An additional 64 trait-associated lincRNAs were identified using a class-level testing strategy at Bonferroni-corrected P < 0.05. To better understand the genomic context and prioritize trait-associated lincRNAs, we examined the pattern of linkage disequilibrium between SNPs in the lincRNAs and SNPs that met genome-wide-significance in the region (±500 kb of lincRNAs). A subset of the lincRNA-trait association findings was replicated in independent Genome-wide association studies data from the Pakistan Risk of Myocardial Infarction Study study. For trait-associated lincRNAs, we also investigated synteny and conservation relative to mouse, expression patterns in five cardiometabolic-relevant tissues, and allele-specific expression in RNA sequencing data for adipose tissue and leukocytes. Finally, we revealed a functional role in human adipocytes for linc-NFE2L3-1, which is expressed in adipose and is associated with waist-hip ratio adjusted for BMI. This comprehensive profile of trait-associated lincRNAs provides novel insights into disease mechanism and serves as a launching point for interrogation of the biology of specific lincRNAs in cardiometabolic disease.

CRISPR/Cas9-Mediated Gene Editing in Human iPSC-Derived Macrophage Reveals Lysosomal Acid Lipase Function in Human Macrophages-Brief Report.

OBJECTIVE: To gain mechanistic insights into the role of LIPA (lipase A), the gene encoding LAL (lysosomal acid lipase) protein, in human macrophages. APPROACH AND RESULTS: We used CRISPR (clustered regularly interspaced short palindromic repeats)/Cas9 (CRISPR-associated protein 9) technology to knock out LIPA in human induced pluripotent stem cells and then differentiate to macrophage (human-induced pluripotent stem cells-derived macrophage [IPSDM]) to explore the human macrophage LIPA loss-of-function phenotypes. LIPA was abundantly expressed in monocyte-derived macrophages and was markedly induced on IPSDM differentiation to comparable levels as in human monocyte-derived macrophage. IPSDM with knockout of LIPA (LIPA-/-) had barely detectable LAL enzymatic activity. Control and LIPA-/- IPSDM were loaded with [3H]-cholesteryl oleate-labeled AcLDL (acetylated low-density lipoprotein) followed by efflux to apolipoprotein A-I. Efflux of liberated [3H]-cholesterol to apolipoprotein A-I was abolished in LIPA-/- IPSDM, indicating deficiency in LAL-mediated lysosomal cholesteryl ester hydrolysis. In cells loaded with [3H]-cholesterol-labeled AcLDL, [3H]-cholesterol efflux was, however, not different between control and LIPA-/- IPSDM. ABCA1 (ATP-binding cassette, subfamily A, member 1) expression was upregulated by AcLDL loading but to a similar extent between control and LIPA-/- IPSDM. In nonlipid loaded state, LIPA-/- IPSDM had high levels of cholesteryl ester mass compared with minute amounts in control IPSDM. Yet, with AcLDL loading, overall cholesteryl ester mass was increased to similar levels in both control and LIPA-/- IPSDM. LIPA-/- did not impact lysosomal apolipoprotein-B degradation or expression of IL1B, IL6, and CCL5. CONCLUSIONS: LIPA-/- IPSDM reveals macrophage-specific hallmarks of LIPA deficiency. CRISPR/Cas9 and IPSDM provide important tools to study human macrophage biology and more broadly for future studies of disease-associated LIPA genetic variation in human macrophages.

Circadian control of innate immunity in macrophages by miR-155 targeting Bmal1.

The response to an innate immune challenge is conditioned by the time of day, but the molecular basis for this remains unclear. In myeloid cells, there is a temporal regulation to induction by lipopolysaccharide (LPS) of the proinflammatory microRNA miR-155 that correlates inversely with levels of BMAL1. BMAL1 in the myeloid lineage inhibits activation of NF-κB and miR-155 induction and protects mice from LPS-induced sepsis. Bmal1 has two miR-155-binding sites in its 3'-UTR, and, in response to LPS, miR-155 binds to these two target sites, leading to suppression of Bmal1 mRNA and protein in mice and humans. miR-155 deletion perturbs circadian function, gives rise to a shorter circadian day, and ablates the circadian effect on cytokine responses to LPS. Thus, the molecular clock controls miR-155 induction that can repress BMAL1 directly. This leads to an innate immune response that is variably responsive to challenges across the circadian day.

De novo RNA sequence assembly during in vivo inflammatory stress reveals hundreds of unannotated lincRNAs in human blood CD14+ monocytes and in adipose tissue.

Long intergenic noncoding RNAs (lincRNAs) have emerged as key regulators of cellular functions and physiology. Yet functional lincRNAs often have low, context-specific and tissue-specific expression. We hypothesized that many human monocyte and adipose lincRNAs would be absent in current public annotations due to lincRNA tissue specificity, modest sequencing depth in public data, limitations of transcriptome assembly algorithms, and lack of dynamic physiological contexts. Deep RNA sequencing (RNA-Seq) was performed in peripheral blood CD14+ monocytes (monocytes; average ~247 million reads per sample) and adipose tissue (average ~378 million reads per sample) collected before and after human experimental endotoxemia, an in vivo inflammatory stress, to identify tissue-specific and clinically relevant lincRNAs. Using a stringent filtering pipeline, we identified 109 unannotated lincRNAs in monocytes and 270 unannotated lincRNAs in adipose. Most unannotated lincRNAs are not conserved in rodents and are tissue specific, while many have features of regulated expression and are enriched in transposable elements. Specific subsets have enhancer RNA characteristics or are expressed only during inflammatory stress. A subset of unannotated lincRNAs was validated and replicated for their presence and inflammatory induction in independent human samples and for their monocyte and adipocyte origins. Through interrogation of public genome-wide association data, we also found evidence of specific disease association for selective unannotated lincRNAs. Our findings highlight the critical need to perform deep RNA-Seq in a cell-, tissue-, and context-specific manner to annotate the full repertoire of human lincRNAs for a complete understanding of lincRNA roles in dynamic cell functions and in human disease.

Integrative genomics identifies 7p11.2 as a novel locus for fever and clinical stress response in humans.

Fever predicts clinical outcomes in sepsis, trauma and during cardiovascular stress, yet the genetic determinants are poorly understood. We used an integrative genomics approach to identify novel genomic determinants of the febrile response to experimental endotoxemia. We highlight multiple integrated lines of evidence establishing the clinical relevance of this novel fever locus. Through genome-wide association study (GWAS) of evoked endotoxemia (lipopolysaccharide (LPS) 1 ng/kg IV) in healthy subjects of European ancestry we discovered a locus on chr7p11.2 significantly associated with the peak febrile response to LPS (top single nucleotide polymorphism (SNP) rs7805622, P = 2.4 × 10(-12)), as well as with temperature fluctuation over time. We replicated this association in a smaller independent LPS study (rs7805622, P = 0.03). In clinical translation, this locus was also associated with temperature and mortality in critically ill patients with trauma or severe sepsis. The top GWAS SNPs are not located within protein-coding genes, but have significant cis-expression quantitative trait loci (eQTL) associations with expression of a cluster of genes ∼400 kb upstream, several of which (SUMF2, CCT6A, GBAS) are regulated by LPS in vivo in blood cells. LPS- and cold-treatment of adipose stromal cells in vitro suggest genotype-specific modulation of eQTL candidate genes (PSPH). Several eQTL genes were up-regulated in brown and white adipose following cold exposure in mice, highlighting a potential role in thermogenesis. Thus, through genomic interrogation of experimental endotoxemia, we identified and replicated a novel fever locus on chr7p11.2 that modulates clinical responses in trauma and sepsis, and highlight integrated in vivo and in vitro evidence for possible novel cis candidate genes conserved across human and mouse.

Whole-exome sequencing identifies rare and low-frequency coding variants associated with LDL cholesterol.

Elevated low-density lipoprotein cholesterol (LDL-C) is a treatable, heritable risk factor for cardiovascular disease. Genome-wide association studies (GWASs) have identified 157 variants associated with lipid levels but are not well suited to assess the impact of rare and low-frequency variants. To determine whether rare or low-frequency coding variants are associated with LDL-C, we exome sequenced 2,005 individuals, including 554 individuals selected for extreme LDL-C (>98(th) or <2(nd) percentile). Follow-up analyses included sequencing of 1,302 additional individuals and genotype-based analysis of 52,221 individuals. We observed significant evidence of association between LDL-C and the burden of rare or low-frequency variants in PNPLA5, encoding a phospholipase-domain-containing protein, and both known and previously unidentified variants in PCSK9, LDLR and APOB, three known lipid-related genes. The effect sizes for the burden of rare variants for each associated gene were substantially higher than those observed for individual SNPs identified from GWASs. We replicated the PNPLA5 signal in an independent large-scale sequencing study of 2,084 individuals. In conclusion, this large whole-exome-sequencing study for LDL-C identified a gene not known to be implicated in LDL-C and provides unique insight into the design and analysis of similar experiments.

Loss-of-function mutations in APOC3, triglycerides, and coronary disease.

BACKGROUND: Plasma triglyceride levels are heritable and are correlated with the risk of coronary heart disease. Sequencing of the protein-coding regions of the human genome (the exome) has the potential to identify rare mutations that have a large effect on phenotype. METHODS: We sequenced the protein-coding regions of 18,666 genes in each of 3734 participants of European or African ancestry in the Exome Sequencing Project. We conducted tests to determine whether rare mutations in coding sequence, individually or in aggregate within a gene, were associated with plasma triglyceride levels. For mutations associated with triglyceride levels, we subsequently evaluated their association with the risk of coronary heart disease in 110,970 persons. RESULTS: An aggregate of rare mutations in the gene encoding apolipoprotein C3 (APOC3) was associated with lower plasma triglyceride levels. Among the four mutations that drove this result, three were loss-of-function mutations: a nonsense mutation (R19X) and two splice-site mutations (IVS2+1G→A and IVS3+1G→T). The fourth was a missense mutation (A43T). Approximately 1 in 150 persons in the study was a heterozygous carrier of at least one of these four mutations. Triglyceride levels in the carriers were 39% lower than levels in noncarriers (P<1×10(-20)), and circulating levels of APOC3 in carriers were 46% lower than levels in noncarriers (P=8×10(-10)). The risk of coronary heart disease among 498 carriers of any rare APOC3 mutation was 40% lower than the risk among 110,472 noncarriers (odds ratio, 0.60; 95% confidence interval, 0.47 to 0.75; P=4×10(-6)). CONCLUSIONS: Rare mutations that disrupt APOC3 function were associated with lower levels of plasma triglycerides and APOC3. Carriers of these mutations were found to have a reduced risk of coronary heart disease. (Funded by the National Heart, Lung, and Blood Institute and others.).

Functional analysis and transcriptomic profiling of iPSC-derived macrophages and their application in modeling Mendelian disease.

RATIONALE: An efficient and reproducible source of genotype-specific human macrophages is essential for study of human macrophage biology and related diseases. OBJECTIVE: To perform integrated functional and transcriptome analyses of human induced pluripotent stem cell-derived macrophages (IPSDMs) and their isogenic human peripheral blood mononuclear cell-derived macrophage (HMDM) counterparts and assess the application of IPSDM in modeling macrophage polarization and Mendelian disease. METHODS AND RESULTS: We developed an efficient protocol for differentiation of IPSDM, which expressed macrophage-specific markers and took up modified lipoproteins in a similar manner to HMDM. Like HMDM, IPSDM revealed reduction in phagocytosis, increase in cholesterol efflux capacity and characteristic secretion of inflammatory cytokines in response to M1 (lipopolysaccharide+interferon-γ) activation. RNA-Seq revealed that nonpolarized (M0) as well as M1 or M2 (interleukin-4) polarized IPSDM shared transcriptomic profiles with their isogenic HMDM counterparts while also revealing novel markers of macrophage polarization. Relative to IPSDM and HMDM of control individuals, patterns of defective cholesterol efflux to apolipoprotein A-I and high-density lipoprotein-3 were qualitatively and quantitatively similar in IPSDM and HMDM of patients with Tangier disease, an autosomal recessive disorder because of mutations in ATP-binding cassette transporter AI. Tangier disease-IPSDM also revealed novel defects of enhanced proinflammatory response to lipopolysaccharide stimulus. CONCLUSIONS: Our protocol-derived IPSDM are comparable with HMDM at phenotypic, functional, and transcriptomic levels. Tangier disease-IPSDM recapitulated hallmark features observed in HMDM and revealed novel inflammatory phenotypes. IPSDMs provide a powerful tool for study of macrophage-specific function in human genetic disorders as well as molecular studies of human macrophage activation and polarization.

Transcriptome-Wide Analysis Reveals Modulation of Human Macrophage Inflammatory Phenotype Through Alternative Splicing.

OBJECTIVE: Human macrophages can shift phenotype across the inflammatory M1 and reparative M2 spectrum in response to environmental challenges, but the mechanisms promoting inflammatory and cardiometabolic disease-associated M1 phenotypes remain incompletely understood. Alternative splicing (AS) is emerging as an important regulator of cellular function, yet its role in macrophage activation is largely unknown. We investigated the extent to which AS occurs in M1 activation within the cardiometabolic disease context and validated a functional genomic cell model for studying human macrophage-related AS events. APPROACH AND RESULTS: From deep RNA-sequencing of resting, M1, and M2 primary human monocyte-derived macrophages, we found 3860 differentially expressed genes in M1 activation and detected 233 M1-induced AS events; the majority of AS events were cell- and M1-specific with enrichment for pathways relevant to macrophage inflammation. Using genetic variant data for 10 cardiometabolic traits, we identified 28 trait-associated variants within the genomic loci of 21 alternatively spliced genes and 15 variants within 7 differentially expressed regulatory splicing factors in M1 activation. Knockdown of 1 such splicing factor, CELF1, in primary human macrophages led to increased inflammatory response to M1 stimulation, demonstrating CELF1's potential modulation of the M1 phenotype. Finally, we demonstrated that an induced pluripotent stem cell-derived macrophage system recapitulates M1-associated AS events and provides a high-fidelity macrophage AS model. CONCLUSIONS: AS plays a role in defining macrophage phenotype in a cell- and stimulus-specific fashion. Alternatively spliced genes and splicing factors with trait-associated variants may reveal novel pathways and targets in cardiometabolic diseases.

Distinct functions of PPARγ isoforms in regulating adipocyte plasticity.

A better understanding of the mechanisms underlying obesity and its comorbidities is key to designing new therapies and treatments. PPARγ is a master regulator of adipocyte biology but the functions of its isoforms are poorly distinguished. Here we demonstrated that PPARγ1 is preferentially expressed in catabolic fat depots while PPARγ2 presents itself at a higher level in browning-resistant depots. PPARγ2, but not PPARγ1, responds to endogenous ligands to induce adipogenesis, and the isoforms regulate distinct sets of white and brown adipocyte genes. Moreover, PPARγ1 negatively correlates while PPARγ2 positively correlates with adiposity in human subcutaneous and visceral fat. These results together indicate that PPARγ1 and PPARγ2 have distinct functions in regulating adipocyte plasticity, and future research should take into account the binary roles of both isoforms in order to identify druggable gene targets and pathways relevant for treatment of metabolic disorders.

PennSeq: accurate isoform-specific gene expression quantification in RNA-Seq by modeling non-uniform read distribution.

Correctly estimating isoform-specific gene expression is important for understanding complicated biological mechanisms and for mapping disease susceptibility genes. However, estimating isoform-specific gene expression is challenging because various biases present in RNA-Seq (RNA sequencing) data complicate the analysis, and if not appropriately corrected, can affect isoform expression estimation and downstream analysis. In this article, we present PennSeq, a statistical method that allows each isoform to have its own non-uniform read distribution. Instead of making parametric assumptions, we give adequate weight to the underlying data by the use of a non-parametric approach. Our rationale is that regardless what factors lead to non-uniformity, whether it is due to hexamer priming bias, local sequence bias, positional bias, RNA degradation, mapping bias or other unknown reasons, the probability that a fragment is sampled from a particular region will be reflected in the aligned data. This empirical approach thus maximally reflects the true underlying non-uniform read distribution. We evaluate the performance of PennSeq using both simulated data with known ground truth, and using two real Illumina RNA-Seq data sets including one with quantitative real time polymerase chain reaction measurements. Our results indicate superior performance of PennSeq over existing methods, particularly for isoforms demonstrating severe non-uniformity. PennSeq is freely available for download at http://sourceforge.net/projects/pennseq.