Chenyi Xue

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.).

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

Significance The circadian clock allows an organism to anticipate daily changes imposed by the environment. The response to LPS is altered depending on time of day; however, the molecular mechanisms underlying this are unclear. We find that the clock in myeloid cells plays a role in LPS-induced sepsis by altering NF-κB activity and the induction of the microRNA miR-155. LPS causes the repression of BMAL1 via the targeting of miR-155 to two seed sequences in the 3′-untranslated region of Bmal1. Lack of miR-155 has profound effects on circadian function and circadian induction of cytokines by LPS. Thus, the molecular clock is using miR-155 as an important regulatory component to control inflammation variably across the circadian day in myeloid cells. 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.

Tissue-Specific RNA-Seq in Human Evoked Inflammation Identifies Blood and Adipose LincRNA Signatures of Cardiometabolic Diseases

Objective—Inappropriate transcriptional activation of innate immunity is a pathological feature of several cardiometabolic disorders, but little is known about inflammatory modulation of long intergenic noncoding RNAs (lincRNAs) in disease-relevant human tissues. Approach and Results—We applied deep RNA sequencing (>500 million filtered reads per sample) to blood and adipose during low-dose experimental endotoxemia (lipopolysaccharide) in a healthy human, with targeted replication in separate individuals undergoing endotoxemia (n=6), to identify inflammatory lincRNAs. A subset of these lincRNAs was examined for expression in adipocytes and monocytes, modulation in adipose of obese humans, and overlap with genome-wide association study signals for inflammatory and cardiometabolic traits. Of a stringent set of 4284 lincRNAs, ≈11% to 22% were expressed with 201 and 56 lincRNAs modulated by lipopolysaccharide in blood or adipose, respectively. Tissue-specific expression of a subset of 6 lipopolysaccharide-lincRNAs was replicated with lipopolysaccharide modulation confirmed for all 3 expressed in blood and 2 of 4 expressed in adipose. The broader generalizability of findings in blood of subject A was confirmed by RNA sequencing in 7 additional subjects. We confirmed adipocytes and monocytes as potential cell-sources of selective lipopolysaccharide-regulated lincRNAs, and 2 of these, linc-DMRT2 (P=0.002) and linc-TP53I13 (P=0.01), were suppressed in adipose of obese humans. Finally, we provide examples of lipopolysaccharide-modulated lincRNAs that overlap single nucleotide polymorphisms that are associated with cardiometabolic traits. Conclusions—Our findings provide novel insights into tissue-level, inflammatory transcriptome regulation in cardiometabolic diseases. These are complementary to more usual approaches limited to interrogation of DNA variations.

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.

Evaluating the Impact of Sequencing Depth on Transcriptome Profiling in Human Adipose

Recent advances in RNA sequencing (RNA-Seq) have enabled the discovery of novel transcriptomic variations that are not possible with traditional microarray-based methods. Tissue and cell specific transcriptome changes during pathophysiological stress in disease cases versus controls and in response to therapies are of particular interest to investigators studying cardiometabolic diseases. Thus, knowledge on the relationships between sequencing depth and detection of transcriptomic variation is needed for designing RNA-Seq experiments and for interpreting results of analyses. Using deeply sequenced Illumina HiSeq 2000 101 bp paired-end RNA-Seq data derived from adipose of a healthy individual before and after systemic administration of endotoxin (LPS), we investigated the sequencing depths needed for studies of gene expression and alternative splicing (AS). In order to detect expressed genes and AS events, we found that ∼100 to 150 million (M) filtered reads were needed. However, the requirement on sequencing depth for the detection of LPS modulated differential expression (DE) and differential alternative splicing (DAS) was much higher. To detect 80% of events, ∼300 M filtered reads were needed for DE analysis whereas at least 400 M filtered reads were necessary for detecting DAS. Although the majority of expressed genes and AS events can be detected with modest sequencing depths (∼100 M filtered reads), the estimated gene expression levels and exon/intron inclusion levels were less accurate. We report the first study that evaluates the relationship between RNA-Seq depth and the ability to detect DE and DAS in human adipose. Our results suggest that a much higher sequencing depth is needed to reliably identify DAS events than for DE genes.

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.

A functional synonymous coding variant in the IL1RN gene is associated with survival in septic shock.

RATIONALE Death from infection is a highly heritable trait, yet there are few genetic variants with known mechanism influencing survival during septic shock. OBJECTIVES We hypothesized that a synonymous coding variant in the IL-1 receptor antagonist gene (IL1RN), rs315952, previously associated with reduced risk for acute respiratory distress syndrome, would be functional and associate with improved survival in septic shock. METHODS We used a human endotoxin (LPS) model of evoked inflammatory stress to measure plasma IL-1 receptor antagonist (IL1RA) following low-dose Food and Drug Administration-grade LPS injection (1 ng/kg) in 294 human volunteers. RNA sequencing of adipose tissue pre- and post-LPS was used to test for allelic imbalance at rs315952. In the Vasopressin and Septic Shock Trial cohort, we performed a genetic association study for survival, mortality, and organ failure-free days. MEASUREMENTS AND MAIN RESULTS Adipose tissue displayed significant allelic imbalance favoring the rs315952C allele in subjects of European ancestry. Consistent with this, carriers of rs315952C had slightly higher plasma IL1RA at baseline (0.039) and higher evoked IL1RA post-LPS (0.011). In the Vasopressin and Septic Shock Trial cohort, rs315952C associated with improved survival (P = 0.028), decreased adjusted 90-day mortality (P = 0.044), and faster resolution of shock (P = 0.029). CONCLUSIONS In European ancestry subjects, the IL1RN variant rs315952C is preferentially transcribed and associated with increased evoked plasma IL1RA and with improved survival from septic shock. It may be that genetically determined IL1RA levels influence survival from septic shock.

The Long Noncoding RNA Landscape in Hypoxic and Inflammatory Renal Epithelial Injury

Long noncoding RNAs (lncRNAs) are emerging as key species-specific regulators of cellular and disease processes. To identify potential lncRNAs relevant to acute and chronic renal epithelial injury, we performed unbiased whole transcriptome profiling of human proximal tubular epithelial cells (PTECs) in hypoxic and inflammatory conditions. RNA sequencing revealed that the protein-coding and noncoding transcriptomic landscape differed between hypoxia-stimulated and cytokine-stimulated human PTECs. Hypoxia- and inflammation-modulated lncRNAs were prioritized for focused followup according to their degree of induction by these stress stimuli, their expression in human kidney tissue, and whether exposure of human PTECs to plasma of critically ill sepsis patients with acute kidney injury modulated their expression. For three lncRNAs (MIR210HG, linc-ATP13A4-8, and linc-KIAA1737-2) that fulfilled our criteria, we validated their expression patterns, examined their loci for conservation and synteny, and defined their associated epigenetic marks. The lncRNA landscape characterized here provides insights into novel transcriptomic variations in the renal epithelial cell response to hypoxic and inflammatory stress.

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.

Adipose tissue RNASeq reveals novel gene–nutrient interactions following n-3 PUFA supplementation and evoked inflammation in humans

Dietary consumption of long-chain omega-3 polyunsaturated fatty acids (n-3 PUFA) may protect against cardiometabolic disease through modulation of systemic and adipose inflammation. However, it is often difficult to detect the subtle effects of n-3 PUFA on inflammatory biomarkers in traditional intervention studies. We aimed to identify novel n-3 PUFA modulated gene expression using unbiased adipose transcriptomics during evoked endotoxemia in a clinical trial of n-3 PUFA supplementation. We analyzed adipose gene expression using RNA sequencing in the fenofibrate and omega-3 fatty acid modulation of endotoxemia (FFAME) trial of healthy individuals at three timepoints: before and after n-3 PUFA supplementation (n=8; 3600mg/day EPA/DHA) for 6weeks compared with placebo (n=6), as well as during a subsequent evoked inflammatory challenge (lipopolysaccharide 0.6ng/kg i.v.). As expected, supplementation with n-3 PUFA vs. placebo alone had only modest effects on adipose tissue gene expression, e.g., increased expression of immediate early response IER2. In contrast, the transcriptomic response to evoked endotoxemia was significantly modified by n-3 PUFA supplementation, with several genes demonstrating significant n-3 PUFA gene-nutrient interactions, e.g., enhanced transcriptional responses in specific immune genes IER5L, HES1, IL1RN, CCL18, IL1RN, IL7R, IL8, CCL3 and others. These data highlight potential mechanisms whereby n-3 PUFA consumption may enhance the immune response to an inflammatory challenge. In conclusion, unbiased transcriptomics during evoked inflammation reveals novel immune modulating functions of n-3 PUFA nutritional intervention in a dynamic pathophysiological setting.