Population-scale tissue transcriptomics maps long non-coding RNAs to complex disease.

Long non-coding RNA (lncRNA) genes have well-established and important impacts on molecular and cellular functions. However, among the thousands of lncRNA genes, it is still a major challenge to identify the subset with disease or trait relevance. To systematically characterize these lncRNA genes, we used Genotype Tissue Expression (GTEx) project v8 genetic and multi-tissue transcriptomic data to profile the expression, genetic regulation, cellular contexts, and trait associations of 14,100 lncRNA genes across 49 tissues for 101 distinct complex genetic traits. Using these approaches, we identified 1,432 lncRNA gene-trait associations, 800 of which were not explained by stronger effects of neighboring protein-coding genes. This included associations between lncRNA quantitative trait loci and inflammatory bowel disease, type 1 and type 2 diabetes, and coronary artery disease, as well as rare variant associations to body mass index.

Integration of rare expression outlier-associated variants improves polygenic risk prediction.

Polygenic risk scores (PRSs) quantify the contribution of multiple genetic loci to an individual's likelihood of a complex trait or disease. However, existing PRSs estimate this likelihood with common genetic variants, excluding the impact of rare variants. Here, we report on a method to identify rare variants associated with outlier gene expression and integrate their impact into PRS predictions for body mass index (BMI), obesity, and bariatric surgery. Between the top and bottom 10%, we observed a 20.8% increase in risk for obesity (p = 3 × 10-14), 62.3% increase in risk for severe obesity (p = 1 × 10-6), and median 5.29 years earlier onset for bariatric surgery (p = 0.008), as a function of expression outlier-associated rare variant burden when controlling for common variant PRS. We show that these predictions were more significant than integrating the effects of rare protein-truncating variants (PTVs), observing a mean 19% increase in phenotypic variance explained with expression outlier-associated rare variants when compared with PTVs (p = 2 × 10-15). We replicated these findings by using data from the Million Veteran Program and demonstrated that PRSs across multiple traits and diseases can benefit from the inclusion of expression outlier-associated rare variants identified through population-scale transcriptome sequencing.

Clonally expanding smooth muscle cells promote atherosclerosis by escaping efferocytosis and activating the complement cascade.

Atherosclerosis is the process underlying heart attack and stroke. Despite decades of research, its pathogenesis remains unclear. Dogma suggests that atherosclerotic plaques expand primarily via the accumulation of cholesterol and inflammatory cells. However, recent evidence suggests that a substantial portion of the plaque may arise from a subset of "dedifferentiated" vascular smooth muscle cells (SMCs) which proliferate in a clonal fashion. Herein we use multicolor lineage-tracing models to confirm that the mature SMC can give rise to a hyperproliferative cell which appears to promote inflammation via elaboration of complement-dependent anaphylatoxins. Despite being extensively opsonized with prophagocytic complement fragments, we find that this cell also escapes immune surveillance by neighboring macrophages, thereby exacerbating its relative survival advantage. Mechanistic studies indicate this phenomenon results from a generalized opsonin-sensing defect acquired by macrophages during polarization. This defect coincides with the noncanonical up-regulation of so-called don't eat me molecules on inflamed phagocytes, which reduces their capacity for programmed cell removal (PrCR). Knockdown or knockout of the key antiphagocytic molecule CD47 restores the ability of macrophages to sense and clear opsonized targets in vitro, allowing for potent and targeted suppression of clonal SMC expansion in the plaque in vivo. Because integrated clinical and genomic analyses indicate that similar pathways are active in humans with cardiovascular disease, these studies suggest that the clonally expanding SMC may represent a translational target for treating atherosclerosis.

Genetic Regulatory Mechanisms of Smooth Muscle Cells Map to Coronary Artery Disease Risk Loci.

Coronary artery disease (CAD) is the leading cause of death globally. Genome-wide association studies (GWASs) have identified more than 95 independent loci that influence CAD risk, most of which reside in non-coding regions of the genome. To interpret these loci, we generated transcriptome and whole-genome datasets using human coronary artery smooth muscle cells (HCASMCs) from 52 unrelated donors, as well as epigenomic datasets using ATAC-seq on a subset of 8 donors. Through systematic comparison with publicly available datasets from GTEx and ENCODE projects, we identified transcriptomic, epigenetic, and genetic regulatory mechanisms specific to HCASMCs. We assessed the relevance of HCASMCs to CAD risk using transcriptomic and epigenomic level analyses. By jointly modeling eQTL and GWAS datasets, we identified five genes (SIPA1, TCF21, SMAD3, FES, and PDGFRA) that may modulate CAD risk through HCASMCs, all of which have relevant functional roles in vascular remodeling. Comparison with GTEx data suggests that SIPA1 and PDGFRA influence CAD risk predominantly through HCASMCs, while other annotated genes may have multiple cell and tissue targets. Together, these results provide tissue-specific and mechanistic insights into the regulation of a critical vascular cell type associated with CAD in human populations.

No evidence of a causal association of type 2 diabetes and glucose metabolism with atrial fibrillation.

AIMS/HYPOTHESIS: Several epidemiological studies have shown an increased risk of atrial fibrillation in individuals with type 2 diabetes or milder forms of dysglycaemia. We aimed to assess whether this relation is causal using a Mendelian randomisation approach. METHODS: Two-sample Mendelian randomisation was used to obtain estimates of the influence of type 2 diabetes, fasting blood glucose (FBG), and HbA1c on the risk of atrial fibrillation. Instrumental variables were constructed using available summary statistics from meta-analyses of genome-wide association studies (GWAS) for type 2 diabetes and associated phenotypes. Pleiotropic SNPs were excluded from the analyses. The most recent GWAS meta-analysis summary statistics for atrial fibrillation, which included over 1 million individuals (approximately 60,000 individuals with atrial fibrillation) was used for outcome analysis. RESULTS: Neither type 2 diabetes (OR 1.01 [95% CI 0.98, 1.03]; p = 0.37), nor FBG (OR 0.95 [95% CI 0.82, 1.09] per mmol/l; p = 0.49) or HbA1c (OR 1.01 [95% CI, 0.85, 1.17] per mmol/mol [%]; p = 0.88) were associated with atrial fibrillation in Mendelian randomisation analyses. We had >80% statistical power to detect ORs of 1.08, 1.06 and 1.09 or larger for type 2 diabetes, FBG and HbA1c, respectively, for associations with atrial fibrillation. CONCLUSIONS/INTERPRETATION: This Mendelian randomisation analysis does not support a causal role of clinical significance between genetically programmed type 2 diabetes, FBG or HbA1c and development of atrial fibrillation. These data suggest that drug treatment to reduce dysglycaemia is unlikely to be an effective strategy for atrial fibrillation prevention. DATA AVAILABILITY: The datasets analysed during the current study are available from the following repository: Nielsen JB, Thorolfsdottir RB, Fritsche LG, et al (2018) GWAS summary statistics for AF (N=60,620 AF cases and 970,216 controls). Center for Statistical Genetics: http://csg.sph.umich.edu/willer/public/afib2018/nielsen-thorolfsdottir-willer-NG2018-AFib-gwas-summary-statistics.tbl.gz.

Identification of 22 novel loci associated with urinary biomarkers of albumin, sodium, and potassium excretion.

Urine biomarkers reflecting kidney function and handling of dietary sodium and potassium are strongly associated with several common diseases including chronic kidney disease, cardiovascular disease, and diabetes mellitus. Knowledge about the genetic determinants of these biomarkers may shed light on pathophysiological mechanisms underlying the development of these diseases. We performed genome-wide association studies of urinary albumin: creatinine ratio (UACR), urinary potassium: creatinine ratio (UK/UCr), urinary sodium: creatinine ratio (UNa/UCr) and urinary sodium: potassium ratio (UNa/UK) in up to 218,450 (discovery) and 109,166 (replication) unrelated individuals of European ancestry from the UK Biobank. Further, we explored genetic correlations, tissue-specific gene expression, and possible genes implicated in the regulation of these biomarkers. After replication, we identified 19 genome-wide significant independent loci associated with UACR, 6 each with UK/UCr and UNa/UCr, and 4 with UNa/UK. In addition to 22 novel associations, we confirmed several established associations, including between the CUBN locus and microalbuminuria. We detected high pairwise genetic correlation across the urinary biomarkers, and between their levels and several physiological measurements. We highlight GIPR, a potential diabetes drug target, as possibly implicated in the genetic control of urinary potassium excretion, and NRBP1, a locus associated with gout, as plausibly involved in sodium and albumin excretion. Overall, we identified 22 novel genome-wide significant associations with urinary biomarkers and confirmed several previously established associations, providing new insights into the genetic basis of these traits and their connection to chronic diseases.

Polygenic risk scores in coronary artery disease.

PURPOSE OF REVIEW: Large genome-wide association studies (GWAS) have identified variants accounting for a substantial portion of the heritable risk for coronary artery disease (CAD). These studies have catalyzed drug discovery and generated the possibility of improved risk prediction and stratification. Here, we review the current state-of-the art in polygenic risk scores (PRSs) and look to the future, as these scores move towards clinical application. RECENT FINDINGS: Over the last decade, multilocus PRSs for CAD have expanded to include millions of variants and demonstrated strong association with CAD outcomes, even when adjusted for traditional risk factors. Recently, PRSs have shown better prediction of CAD outcomes than any single traditional risk factor alone. Advances in statistical methods used to generate PRSs have improved their predictive ability and transferability between populations with varied ancestries. Initial clinical studies have also demonstrated the potential of genetic information to impact shared decision-making between patients and providers, leading to improved outcomes. SUMMARY: PRSs can improve risk stratification for CAD especially in white/European populations and have the potential to alter routine clinical care. However, unlocking this potential will require additional research in PRSs in nonwhite populations and substantial investment in clinical implementation studies.

CRISPR-Cas9-mediated knockout of SPRY2 in human hepatocytes leads to increased glucose uptake and lipid droplet accumulation.

BACKGROUND: The prevalence of obesity and its comorbidities, including type 2 diabetes mellitus (T2DM), is dramatically increasing throughout the world; however, the underlying aetiology is incompletely understood. Genome-wide association studies (GWAS) have identified hundreds of genec susceptibility loci for obesity and T2DM, although the causal genes and mechanisms are largely unknown. SPRY2 is a candidate gene identified in GWAS of body fat percentage and T2DM, and has recently been linked to insulin production in pancreatic ß-cells. In the present study, we aimed to further understand SPRY2 via functional characterisation in HepG2 cells, an in vitro model of human hepatocytes widely used to investigate T2DM and insulin resistance. METHODS: CRISPR-Cas9 genome editing was used to target SPRY2 in HepG2 cells, and the functional consequences of SPRY2 knockout (KO) and overexpression subsequently assessed using glucose uptake and lipid droplet assays, measurement of protein kinase phosphorylation and RNA sequencing. RESULTS: The major functional consequence of SPRY2 KO was a significant increase in glucose uptake, along with elevated lipid droplet accumulation. These changes were attenuated, but not reversed, in cells overexpressing SPRY2. Phosphorylation of protein kinases across key signalling pathways (including Akt and mitogen activated protein kinases) was not altered after SPRY2 KO. Transcriptome profiling in SPRY2 KO and mock (control) cells revealed a number of differentially expressed genes related to cholesterol biosynthesis, cell cycle regulation and cellular signalling pathways. Phospholipase A2 group IIA (PLA2G2A) mRNA level was subsequently validated as significantly upregulated following SPRY2 KO, highlighting this as a potential mediator downstream of SPRY2. CONCLUSION: These findings suggest a role for SPRY2 in glucose and lipid metabolism in hepatocytes and contribute to clarifying the function of this gene in the context of metabolic diseases.

Single-cell transcriptome dataset of human and mouse in vitro adipogenesis models.

Adipogenesis is a process in which fat-specific progenitor cells (preadipocytes) differentiate into adipocytes that carry out the key metabolic functions of the adipose tissue, including glucose uptake, energy storage, and adipokine secretion. Several cell lines are routinely used to study the molecular regulation of adipogenesis, in particular the immortalized mouse 3T3-L1 line and the primary human Simpson-Golabi-Behmel syndrome (SGBS) line. However, the cell-to-cell variability of transcriptional changes prior to and during adipogenesis in these models is not well understood. Here, we present a single-cell RNA-Sequencing (scRNA-Seq) dataset collected before and during adipogenic differentiation of 3T3-L1 and SGBS cells. To minimize the effects of experimental variation, we mixed 3T3-L1 and SGBS cells and used computational analysis to demultiplex transcriptomes of mouse and human cells. In both models, adipogenesis results in the appearance of three cell clusters, corresponding to preadipocytes, early and mature adipocytes. These data provide a groundwork for comparative studies on these widely used in vitro models of human and mouse adipogenesis, and on cell-to-cell variability during this process.

Single-cell transcriptome dataset of human and mouse in vitro adipogenesis models.

Adipogenesis is a process in which fat-specific progenitor cells (preadipocytes) differentiate into adipocytes that carry out the key metabolic functions of the adipose tissue, including glucose uptake, energy storage, and adipokine secretion. Several cell lines are routinely used to study the molecular regulation of adipogenesis, in particular the immortalized mouse 3T3-L1 line and the primary human Simpson-Golabi-Behmel syndrome (SGBS) line. However, the cell-to-cell variability of transcriptional changes prior to and during adipogenesis in these models is not well understood. Here, we present a single-cell RNA-Sequencing (scRNA-Seq) dataset collected before and during adipogenic differentiation of 3T3-L1 and SGBS cells. To minimize the effects of experimental variation, we mixed 3T3-L1 and SGBS cells and used computational analysis to demultiplex transcriptomes of mouse and human cells. In both models, adipogenesis results in the appearance of three cell clusters, corresponding to preadipocytes, early and mature adipocytes. These data provide a groundwork for comparative studies on human and mouse adipogenesis, as well as on cell-to-cell variability in gene expression during this process.