A mechanistic framework for cardiometabolic and coronary artery diseases

Single-Cell Gene-Regulatory Networks of Advanced Symptomatic Atherosclerosis

BACKGROUND

While our understanding of the single-cell gene expression patterns underlying the transformation of vascular cell types during the progression of atherosclerosis is rapidly improving, the clinical and pathophysiological relevance of these changes remains poorly understood.

METHODS

Single-cell RNA sequencing data generated with SmartSeq2 (≈8000 genes/cell) in 16 588 single cells isolated during atherosclerosis progression in Ldlr-/-Apob100/100 mice with human-like plasma lipoproteins and from humans with asymptomatic and symptomatic carotid plaques was clustered into multiple subtypes. For clinical and pathophysiological context, the advanced-stage and symptomatic subtype clusters were integrated with 135 tissue-specific (atherosclerotic aortic wall, mammary artery, liver, skeletal muscle, and visceral and subcutaneous, fat) gene-regulatory networks (GRNs) inferred from 600 coronary artery disease patients in the STARNET (Stockholm-Tartu Atherosclerosis Reverse Network Engineering Task) study.

RESULTS

Advanced stages of atherosclerosis progression and symptomatic carotid plaques were largely characterized by 3 smooth muscle cells (SMCs), and 3 macrophage subtype clusters with extracellular matrix organization/osteogenic (SMC), and M1-type proinflammatory/Trem2-high lipid-associated (macrophage) phenotypes. Integrative analysis of these 6 clusters with STARNET revealed significant enrichments of 3 arterial wall GRNs: GRN33 (macrophage), GRN39 (SMC), and GRN122 (macrophage) with major contributions to coronary artery disease heritability and strong associations with clinical scores of coronary atherosclerosis severity. The presence and pathophysiological relevance of GRN39 were verified in 5 independent RNAseq data sets obtained from the human coronary and aortic artery, and primary SMCs and by targeting its top-key drivers, FRZB and ALCAM in cultured human coronary artery SMCs.

CONCLUSIONS

By identifying and integrating the most gene-rich single-cell subclusters of atherosclerosis to date with a coronary artery disease framework of GRNs, GRN39 was identified and independently validated as being critical for the transformation of contractile SMCs into an osteogenic phenotype promoting advanced, symptomatic atherosclerosis.

Omics Approaches Unveiling the Biology of Human Atherosclerotic Plaques

Atherosclerosis is a chronic inflammatory disease of the arterial wall, characterized by the buildup of plaques with the accumulation and transformation of lipids, immune cells, vascular smooth muscle cells, and necrotic cell debris. Plaques with collagen-poor thin fibrous caps infiltrated by macrophages and lymphocytes are considered unstable because they are at the greatest risk of rupture and clinical events. However, the current histologic definition of plaque types may not fully capture the complex molecular nature of atherosclerotic plaque biology and the underlying mechanisms contributing to plaque progression, rupture, and erosion. The advances in omics technologies have changed the understanding of atherosclerosis plaque biology, offering new possibilities to improve risk prediction and discover novel therapeutic targets. Genomic studies have shed light on the genetic predisposition to atherosclerosis, and integrative genomic analyses expedite the translation of genomic discoveries. Transcriptomic, proteomic, metabolomic, and lipidomic studies have refined the understanding of the molecular signature of atherosclerotic plaques, aiding in data-driven hypothesis generation for mechanistic studies and offering new prospects for biomarker discovery. Furthermore, advancements in single-cell technologies and emerging spatial analysis techniques have unveiled the heterogeneity and plasticity of plaque cells. This review discusses key omics-based discoveries that have advanced the understanding of human atherosclerotic plaque biology, focusing on insights derived from omics profiling of human atherosclerotic vascular specimens.

Mining cholesterol genes from thousands of mouse livers identifies aldolase C as a regulator of cholesterol biosynthesis

The availability of genome-wide transcriptomic and proteomic datasets is ever-increasing and often not used beyond initial publication. Here, we applied module-based coexpression network analysis to a comprehensive catalog of 35 mouse genome-wide liver expression datasets (encompassing more than 3800 mice) with the goal of identifying and validating unknown genes involved in cholesterol metabolism. From these 35 datasets, we identified a conserved module of genes enriched with cholesterol biosynthetic genes. Using a systematic approach across the 35 datasets, we identified three genes (Rdh11, Echdc1, and Aldoc) with no known role in cholesterol metabolism. We then performed functional validation studies and show that each gene is capable of regulating cholesterol metabolism. For the glycolytic gene, Aldoc, we demonstrate that it contributes to de novo cholesterol biosynthesis and regulates cholesterol and triglyceride levels in mice. As Aldoc is located within a genome-wide significant genome-wide association studies locus for human plasma cholesterol levels, our studies establish Aldoc as a causal gene within this locus. Through our work, we develop a framework for leveraging mouse genome-wide liver datasets for identifying and validating genes involved in cholesterol metabolism.

Leveraging inter-individual transcriptional correlation structure to infer discrete signaling mechanisms across metabolic tissues

Inter-organ communication is a vital process to maintain physiologic homeostasis, and its dysregulation contributes to many human diseases. Given that circulating bioactive factors are stable in serum, occur naturally, and are easily assayed from blood, they present obvious focal molecules for therapeutic intervention and biomarker development. Recently, studies have shown that secreted proteins mediating inter-tissue signaling could be identified by 'brute force' surveys of all genes within RNA-sequencing measures across tissues within a population. Expanding on this intuition, we reasoned that parallel strategies could be used to understand how individual genes mediate signaling across metabolic tissues through correlative analyses of gene variation between individuals. Thus, comparison of quantitative levels of gene expression relationships between organs in a population could aid in understanding cross-organ signaling. Here, we surveyed gene-gene correlation structure across 18 metabolic tissues in 310 human individuals and 7 tissues in 103 diverse strains of mice fed a normal chow or high-fat/high-sucrose (HFHS) diet. Variation of genes such as FGF21, ADIPOQ, GCG, and IL6 showed enrichments which recapitulate experimental observations. Further, similar analyses were applied to explore both within-tissue signaling mechanisms (liver PCSK9) and genes encoding enzymes producing metabolites (adipose PNPLA2), where inter-individual correlation structure aligned with known roles for these critical metabolic pathways. Examination of sex hormone receptor correlations in mice highlighted the difference of tissue-specific variation in relationships with metabolic traits. We refer to this resource as gene-derived correlations across tissues (GD-CAT) where all tools and data are built into a web portal enabling users to perform these analyses without a single line of code (gdcat.org). This resource enables querying of any gene in any tissue to find correlated patterns of genes, cell types, pathways, and network architectures across metabolic organs.

Multi-ancestry genetic analysis of gene regulation in coronary arteries prioritizes disease risk loci

Genome-wide association studies (GWASs) have identified hundreds of risk loci for coronary artery disease (CAD). However, non-European populations are underrepresented in GWASs, and the causal gene-regulatory mechanisms of these risk loci during atherosclerosis remain unclear. We incorporated local ancestry and haplotypes to identify quantitative trait loci for expression (eQTLs) and splicing (sQTLs) in coronary arteries from 138 ancestrally diverse Americans. Of 2,132 eQTL-associated genes (eGenes), 47% were previously unreported in coronary artery; 19% exhibited cell-type-specific expression. Colocalization revealed subgroups of eGenes unique to CAD and blood pressure GWAS. Fine-mapping highlighted additional eGenes, including TBX20 and IL5. We also identified sQTLs for 1,690 genes, among which TOR1AIP1 and ULK3 sQTLs demonstrated the importance of evaluating splicing to accurately identify disease-relevant isoform expression. Our work provides a patient-derived coronary artery eQTL resource and exemplifies the need for diverse study populations and multifaceted approaches to characterize gene regulation in disease processes.

Atherosclerotic plaque development in mice is enhanced by myeloid ZEB1 downregulation

Accumulation of lipid-laden macrophages within the arterial neointima is a critical step in atherosclerotic plaque formation. Here, we show that reduced levels of the cellular plasticity factor ZEB1 in macrophages increase atherosclerotic plaque formation and the chance of cardiovascular events. Compared to control counterparts (Zeb1WT/ApoeKO), male mice with Zeb1 ablation in their myeloid cells (Zeb1∆M/ApoeKO) have larger atherosclerotic plaques and higher lipid accumulation in their macrophages due to delayed lipid traffic and deficient cholesterol efflux. Zeb1∆M/ApoeKO mice display more pronounced systemic metabolic alterations than Zeb1WT/ApoeKO mice, with higher serum levels of low-density lipoproteins and inflammatory cytokines and larger ectopic fat deposits. Higher lipid accumulation in Zeb1∆M macrophages is reverted by the exogenous expression of Zeb1 through macrophage-targeted nanoparticles. In vivo administration of these nanoparticles reduces atherosclerotic plaque formation in Zeb1∆M/ApoeKO mice. Finally, low ZEB1 expression in human endarterectomies is associated with plaque rupture and cardiovascular events. These results set ZEB1 in macrophages as a potential target in the treatment of atherosclerosis.

Gene Regulatory Networks in Coronary Artery Disease

PURPOSE OF REVIEW

Coronary artery disease is a complex disorder and the leading cause of mortality worldwide. As technologies for the generation of high-throughput multiomics data have advanced, gene regulatory network modeling has become an increasingly powerful tool in understanding coronary artery disease. This review summarizes recent and novel gene regulatory network tools for bulk tissue and single cell data, existing databases for network construction, and applications of gene regulatory networks in coronary artery disease.

RECENT FINDINGS

New gene regulatory network tools can integrate multiomics data to elucidate complex disease mechanisms at unprecedented cellular and spatial resolutions. At the same time, updates to coronary artery disease expression data in existing databases have enabled researchers to build gene regulatory networks to study novel disease mechanisms. Gene regulatory networks have proven extremely useful in understanding CAD heritability beyond what is explained by GWAS loci and in identifying mechanisms and key driver genes underlying disease onset and progression. Gene regulatory networks can holistically and comprehensively address the complex nature of coronary artery disease. In this review, we discuss key algorithmic approaches to construct gene regulatory networks and highlight state-of-the-art methods that model specific modes of gene regulation. We also explore recent applications of these tools in coronary artery disease patient data repositories to understand disease heritability and shared and distinct disease mechanisms and key driver genes across tissues, between sexes, and between species.

Improving liver transplant outcomes with transplant-omics and network biology

PURPOSE OF REVIEW

Molecular omics data is increasingly ubiquitous throughout medicine. In organ transplantation, recent large-scale research efforts are generating the 'transplant-ome' - the entire set of molecular omics data, including the genome, transcriptome, proteome, and metabolome. Importantly, early studies in anesthesiology have demonstrated how perioperative interventions alter molecular profiles in various patient populations. The next step for anesthesiologists and intensivists will be to tailor perioperative care to the transplant-ome of individual liver transplant patients.

RECENT FINDINGS

In liver transplant patients, elements of the transplant-ome predict complications and point to novel interventions. Importantly, molecular profiles of both the donor organ and recipient contribute to this risk, and interventions like normothermic machine perfusion influence these profiles. As we can now measure various omics molecules simultaneously, we can begin to understand how these molecules interact to form molecular networks and emerging technologies offer noninvasive and continuous ways to measure these networks throughout the perioperative period. Molecules that regulate these networks are likely mediators of complications and actionable clinical targets throughout the perioperative period.

SUMMARY

The transplant-ome can be used to tailor perioperative care to the individual liver transplant patient. Monitoring molecular networks continuously and noninvasively would provide new opportunities to optimize perioperative management.

Profiling the genome and proteome of metabolic dysfunction-associated steatotic liver disease identifies potential therapeutic targets

BACKGROUND & AIMS

Metabolic dysfunction-associated steatotic liver disease (MASLD) affects over 25% of the population and currently has no effective treatments. Plasma proteins with causal evidence may represent promising drug targets. We aimed to identify plasma proteins in the causal pathway of MASLD and explore their interaction with obesity.

METHODS

We analysed 2,941 plasma proteins in 43,978 European participants from UK Biobank. We performed genome-wide association study (GWAS) for all MASLD-associated proteins and created the largest MASLD GWAS (109,885 cases/1,014,923 controls). We performed Mendelian Randomization (MR) and integrated proteins and their encoding genes in MASLD ranges to identify candidate causal proteins. We then validated them through independent replication, exome sequencing, liver imaging, bulk and single-cell gene expression, liver biopsies, pathway, and phenome-wide data. We explored the role of obesity by MR and multivariable MR across proteins, body mass index, and MASLD.

RESULTS

We found 929 proteins associated with MASLD, reported five novel genetic loci associated with MASLD, and identified 17 candidate MASLD protein targets. We identified four novel targets for MASLD (CD33, GRHPR, HMOX2, and SCG3), provided protein evidence supporting roles of AHCY, FCGR2B, ORM1, and RBKS in MASLD, and validated nine previously known targets. We found that CD33, FCGR2B, ORM1, RBKS, and SCG3 mediated the association of obesity and MASLD, and HMOX2, ORM1, and RBKS had effect on MASLD independent of obesity.

CONCLUSIONS

This study identified new protein targets in the causal pathway of MASLD, providing new insights into the multi-omics architecture and pathophysiology of MASLD. These findings advise further therapeutic interventions for MASLD.

Integration of transcriptomes of senescent cell models with multi-tissue patient samples reveals reduced COL6A3 as an inducer of senescence

Senescent cells are a major contributor to age-dependent cardiovascular tissue dysfunction, but knowledge of their in vivo cell markers and tissue context is lacking. To reveal tissue-relevant senescence biology, we integrate the transcriptomes of 10 experimental senescence cell models with a 224 multi-tissue gene co-expression network based on RNA-seq data of seven tissues biopsies from ∼600 coronary artery disease (CAD) patients. We identify 56 senescence-associated modules, many enriched in CAD GWAS genes and correlated with cardiometabolic traits-which supports universality of senescence gene programs across tissues and in CAD. Cross-tissue network analyses reveal 86 candidate senescence-associated secretory phenotype (SASP) factors, including COL6A3. Experimental knockdown of COL6A3 induces transcriptional changes that overlap the majority of the experimental senescence models, with cell-cycle arrest linked to modulation of DREAM complex-targeted genes. We provide a transcriptomic resource for cellular senescence and identify candidate biomarkers, SASP factors, and potential drivers of senescence in human tissues.

Integrative single-cell meta-analysis reveals disease-relevant vascular cell states and markers in human atherosclerosis

Coronary artery disease (CAD) is characterized by atherosclerotic plaque formation in the arterial wall. CAD progression involves complex interactions and phenotypic plasticity among vascular and immune cell lineages. Single-cell RNA-seq (scRNA-seq) studies have highlighted lineage-specific transcriptomic signatures, but human cell phenotypes remain controversial. Here, we perform an integrated meta-analysis of 22 scRNA-seq libraries to generate a comprehensive map of human atherosclerosis with 118,578 cells. Besides characterizing granular cell-type diversity and communication, we leverage this atlas to provide insights into smooth muscle cell (SMC) modulation. We integrate genome-wide association study data and uncover a critical role for modulated SMC phenotypes in CAD, myocardial infarction, and coronary calcification. Finally, we identify fibromyocyte/fibrochondrogenic SMC markers (LTBP1 and CRTAC1) as proxies of atherosclerosis progression and validate these through omics and spatial imaging analyses. Altogether, we create a unified atlas of human atherosclerosis informing cell state-specific mechanistic and translational studies of cardiovascular diseases.

Systems genetics approaches for understanding complex traits with relevance for human disease

Quantitative traits are often complex because of the contribution of many loci, with further complexity added by environmental factors. In medical research, systems genetics is a powerful approach for the study of complex traits, as it integrates intermediate phenotypes, such as RNA, protein, and metabolite levels, to understand molecular and physiological phenotypes linking discrete DNA sequence variation to complex clinical and physiological traits. The primary purpose of this review is to describe some of the resources and tools of systems genetics in humans and rodent models, so that researchers in many areas of biology and medicine can make use of the data.

Sex differences in the genetic and molecular mechanisms of coronary artery disease

Sex differences in coronary artery disease (CAD) presentation, risk factors and prognosis have been widely studied. Similarly, studies on atherosclerosis have shown prominent sex differences in plaque biology. Our understanding of the underlying genetic and molecular mechanisms that drive these differences remains fragmented and largely understudied. Through reviewing genetic and epigenetic studies, we identified more than 40 sex-differential candidate genes (13 within known CAD loci) that may explain, at least in part, sex differences in vascular remodeling, lipid metabolism and endothelial dysfunction. Studies with transcriptomic and single-cell RNA sequencing data from atherosclerotic plaques highlight potential sex differences in smooth muscle cell and endothelial cell biology. Especially, phenotypic switching of smooth muscle cells seems to play a crucial role in female atherosclerosis. This matches the known sex differences in atherosclerotic phenotypes, with men being more prone to lipid-rich plaques, while women are more likely to develop fibrous plaques with endothelial dysfunction. To unravel the complex mechanisms that drive sex differences in CAD, increased statistical power and adjustments to study designs and analysis strategies are required. This entails increasing inclusion rates of women, performing well-defined sex-stratified analyses and the integration of multi-omics data.

FSTL3 partially mediates the association of increased nonalcoholic fatty liver disease fibrosis risk with acute myocardial infarction in patients with type 2 diabetes mellitus

BACKGROUND

The study aimed to investigate an association of increased liver fibrosis with acute myocardial infarction (AMI), and to investigate the mediating effect of serum follistatin-like protein 3 (FSTL3) on the association in patients with type 2 diabetes mellitus (T2DM).

METHOD

A total of 1424 participants were included in this study, and were firstly divided into two groups: 429 T2DM patients and 995 T2DM patients with NAFLD to assess the association of NAFLD and AMI. Then 995 T2DM co-existent NAFLD patients were categorized by NAFLD fibrosis risk to explore the association between NAFLD fibrosis risk and AMI. Immunohistochemistry staining and semi-quantitative analysis of liver FSTL3 were performed in 60 patients with NAFLD. There were 323 individuals (191 without AMI and 132 with AMI) in T2DM co-existent NAFLD patients who had serum samples, and serum FSTL3 was tested and mediation effect of FSTL3 in association of NAFLD fibrosis and AMI was performed.

RESULTS

First, increased NAFLD fibrosis risk was an independent risk factor for AMI in patients with T2DM and co-existent NAFLD. In addition, analysis of Gene Expression Omnibus (GEO) database and immunohistochemical staining confirmed the increased expression of FSTL3 in the liver of NAFLD patients with fibrosis. Serum FSTL3 significantly increased in patients with high NAFLD fibrosis risk and AMI, and closely associated with NAFLD fibrosis and AMI severity in T2DM patients with co-existent NAFLD. Most importantly, analysis of the level of mediation revealed that increased serum FSTL3 partially mediated the association of increased NAFLD fibrosis risk with AMI in T2DM patients with co-existent NAFLD.

CONCLUSIONS

NAFLD fibrosis was closely associated with AMI in T2DM patients. FSTL3 expression was enriched in the liver of NAFLD patients with significant and advanced fibrosis, and serum FSTL3 partially mediated the association of increased liver fibrosis risk with AMI in T2DM patients.

Evolution of the search for a common mechanism of congenital risk of coronary heart disease and type 2 diabetes mellitus in the chromosomal locus 9p21.3

9.21.3 chromosomal locus predisposes to coronary heart disease (CHD) and type 2 diabetes mellitus (DM2), but their overall pathological mechanism and clinical applicability remain unclear. The review uses publications of the study results of 9.21.3 chromosomal locus in association with CHD and DM2, which are important for changing the focus of clinical practice. The eligibility criteria are full-text articles published in the PubMed database (MEDLINE) up to December 31, 2022. A total of 56 publications were found that met the inclusion criteria. Using the examples of the progressive stages in understanding the role of the chromosomal locus 9p.21.3, scientific ideas were grouped, from a fragmentary study of independent pathological processes to a systematic study of the overall development of CHD and DM2. The presented review can become a source of new scientific hypotheses for further studies, the results of which can determine the general mechanism of the congenital risk of CHD and DM2 and change the focus of clinical practice.

Leveraging inter-individual transcriptional correlation structure to infer discrete signaling mechanisms across metabolic tissues

Inter-organ communication is a vital process to maintain physiologic homeostasis, and its dysregulation contributes to many human diseases. Beginning with the discovery of insulin over a century ago, characterization of molecules responsible for signal between tissues has required careful and elegant experimentation where these observations have been integral to deciphering physiology and disease. Given that circulating bioactive factors are stable in serum, occur naturally, and are easily assayed from blood, they present obvious focal molecules for therapeutic intervention and biomarker development. For example, physiologic dissection of the actions of soluble proteins such as proprotein convertase subtilisin/kexin type 9 (PCSK9) and glucagon-like peptide 1 (GLP1) have yielded among the most promising therapeutics to treat cardiovascular disease and obesity, respectively1–4. A major obstacle in the characterization of such soluble factors is that defining their tissues and pathways of action requires extensive experimental testing in cells and animal models. Recently, studies have shown that secreted proteins mediating inter-tissue signaling could be identified by “brute-force” surveys of all genes within RNA-sequencing measures across tissues within a population5–9. Expanding on this intuition, we reasoned that parallel strategies could be used to understand how individual genes mediate signaling across metabolic tissues through correlative analyses of gene variation between individuals. Thus, comparison of quantitative levels of gene expression relationships between organs in a population could aid in understanding cross-organ signaling. Here, we surveyed gene-gene correlation structure across 18 metabolic tissues in 310 human individuals and 7 tissues in 103 diverse strains of mice fed a normal chow or HFHS diet. Variation of genes such asFGF21, ADIPOQ, GCGandIL6showed enrichments which recapitulate experimental observations. Further, similar analyses were applied to explore both within-tissue signaling mechanisms (liverPCSK9) as well as genes encoding enzymes producing metabolites (adiposePNPLA2), where inter-individual correlation structure aligned with known roles for these critical metabolic pathways. Examination of sex hormone receptor correlations in mice highlighted the difference of tissue-specific variation in relationships with metabolic traits. We refer to this resource asGene-DerivedCorrelationsAcrossTissues (GD-CAT) where all tools and data are built into a web portal enabling users to perform these analyses without a single line of code (gdcat.org). This resource enables querying of any gene in any tissue to find correlated patterns of genes, cell types, pathways and network architectures across metabolic organs.

Causal inference in drug discovery and development

To discover new drugs is to seek and to prove causality. As an emerging approach leveraging human knowledge and creativity, data, and machine intelligence, causal inference holds the promise of reducing cognitive bias and improving decision-making in drug discovery. Although it has been applied across the value chain, the concepts and practice of causal inference remain obscure to many practitioners. This article offers a nontechnical introduction to causal inference, reviews its recent applications, and discusses opportunities and challenges of adopting the causal language in drug discovery and development.

Lowering of Circulating Sclerostin May Increase Risk of Atherosclerosis and Its Risk Factors: Evidence From a Genome-Wide Association Meta-Analysis Followed by Mendelian Randomization

OBJECTIVE

In this study, we aimed to establish the causal effects of lowering sclerostin, target of the antiosteoporosis drug romosozumab, on atherosclerosis and its risk factors.

METHODS

A genome-wide association study meta-analysis was performed of circulating sclerostin levels in 33,961 European individuals. Mendelian randomization (MR) was used to predict the causal effects of sclerostin lowering on 15 atherosclerosis-related diseases and risk factors.

RESULTS

We found that 18 conditionally independent variants were associated with circulating sclerostin. Of these, 1 cis signal in SOST and 3 trans signals in B4GALNT3, RIN3, and SERPINA1 regions showed directionally opposite signals for sclerostin levels and estimated bone mineral density. Variants with these 4 regions were selected as genetic instruments. MR using 5 correlated cis-SNPs suggested that lower sclerostin increased the risk of type 2 diabetes mellitus (DM) (odds ratio [OR] 1.32 [95% confidence interval (95% CI) 1.03-1.69]) and myocardial infarction (MI) (OR 1.35 [95% CI 1.01-1.79]); sclerostin lowering was also suggested to increase the extent of coronary artery calcification (CAC) (β = 0.24 [95% CI 0.02-0.45]). MR using both cis and trans instruments suggested that lower sclerostin increased hypertension risk (OR 1.09 [95% CI 1.04-1.15]), but otherwise had attenuated effects.

CONCLUSION

This study provides genetic evidence to suggest that lower levels of sclerostin may increase the risk of hypertension, type 2 DM, MI, and the extent of CAC. Taken together, these findings underscore the requirement for strategies to mitigate potential adverse effects of romosozumab treatment on atherosclerosis and its related risk factors.

Identifying shared transcriptional risk patterns between atherosclerosis and cancer

Cancer and cardiovascular disease (CVD) are the leading causes of death worldwide. Numerous overlapping pathophysiologic mechanisms have been hypothesized to drive the development of both diseases. Further investigation of these common pathways could allow for the identification of mutually detrimental processes and therapeutic targeting to derive mutual benefit. In this study, we intersect transcriptomic datasets correlated with disease severity or patient outcomes for both cancer and atherosclerotic CVD. These analyses confirmed numerous pathways known to underlie both diseases, such as inflammation and hypoxia, but also identified several novel shared pathways. We used these to explore common translational targets by applying the drug prediction software, OCTAD, to identify compounds that simultaneously reverse the gene expression signature for both diseases. These analyses suggest that certain tumor-specific therapeutic approaches may be implemented so that they avoid cardiovascular consequences, and in some cases may even be used to simultaneously target co-prevalent cancer and atherosclerosis.

The Genetics of Coronary Artery Disease: A Vascular Perspective

Genome-wide association studies (GWAS) have identified a large number of genetic loci for coronary artery disease (CAD), with many located close to genes associated with traditional CAD risk pathways, such as lipid metabolism and inflammation. It is becoming evident with recent CAD GWAS meta-analyses that vascular pathways are also highly enriched and present an opportunity for novel therapeutics. This review examines GWAS-enriched vascular gene loci, the pathways involved and their potential role in CAD pathogenesis. The functionality of variants is explored from expression quantitative trait loci, massively parallel reporter assays and CRISPR-based gene-editing tools. We discuss how this research may lead to novel therapeutic tools to treat cardiovascular disorders.

Cell-type-specific co-expression inference from single cell RNA-sequencing data

The advancement of single cell RNA-sequencing (scRNA-seq) technology has enabled the direct inference of co-expressions in specific cell types, facilitating our understanding of cell-type-specific biological functions. For this task, the high sequencing depth variations and measurement errors in scRNA-seq data present two significant challenges, and they have not been adequately addressed by existing methods. We propose a statistical approach, CS-CORE, for estimating and testing cell-type-specific co-expressions, that explicitly models sequencing depth variations and measurement errors in scRNA-seq data. Systematic evaluations show that most existing methods suffered from inflated false positives as well as biased co-expression estimates and clustering analysis, whereas CS-CORE gave accurate estimates in these experiments. When applied to scRNA-seq data from postmortem brain samples from Alzheimer's disease patients/controls and blood samples from COVID-19 patients/controls, CS-CORE identified cell-type-specific co-expressions and differential co-expressions that were more reproducible and/or more enriched for relevant biological pathways than those inferred from existing methods.

Causal relationship between osteoarthritis with atrial fibrillation and coronary atherosclerosis: a bidirectional Mendelian randomization study of European ancestry

Background

Osteoarthritis (OA) is a degenerative disease with high prevalence. Some observational studies have shown that patients with osteoarthritis often have co-existing cardiovascular diseases (CVD) such as atrial fibrillation (AF) and coronary atherosclerosis (CA). However, there is still a lack of stronger evidence confirming the association between osteoarthritis and cardiovascular disease. In this study, we used a bidirectional two-sample Mendelian randomization study to investigate the relationship between OA with AF and CA.

Methods

OA data from the UK Biobank and arcOGEN (Arthritis Research UK Osteoarthritis Genetics, a study that aimed to find genetic determinants of osteoarthritis and elucidate the genetic architecture of the disease) integration were selected for the study (n = 417,596), AF data were obtained from six studies (n = 1,030,836), and coronary atherosclerosis data were derived from the FinnGen (n = 218,792). MR analysis was performed primarily using the Inverse variance weighted (IVW) method, with MR Egger, weighted median, simple mode, weighted mode as supplements, sensitivity analysis was performed using Cochran Q statistic, and leave-one-out analysis.

Results

We found that OA and AF were positively associated [IVW: OR (95% CI): 1.11 (1.04, 1.19), P = 0.002], while OA and CA were negatively associated [IVW: OR (95% CI): 0.88 (0.79, 0.98), P = 0.02]. In the reverse MR analysis, no effect of AF on OA was found [IVW: OR (95% CI): 1.00 (0.97, 1.03), P = 0.84], meanwhile, CA and OA were found to be associated negatively [IVW: OR (95% CI): 0.95 (0.92, 0.99), P = 0.01]. No violations of MR assumptions were found in the sensitivity analysis.

Conclusion

This research confirms that OA is a risk factor for AF, and there is a mutual protective factor between OA and CA. However, further studies are still necessary to elucidate the underlying mechanisms.

Leukocytes carrying Clonal Hematopoiesis of Indeterminate Potential (CHIP) Mutations invade Human Atherosclerotic Plaques

Background

Leukocyte progenitors derived from clonal hematopoiesis of undetermined potential (CHIP) are associated with increased cardiovascular events. However, the prevalence and functional relevance of CHIP in coronary artery disease (CAD) are unclear, and cells affected by CHIP have not been detected in human atherosclerotic plaques.

Methods

CHIP mutations in blood and tissues were identified by targeted deep-DNA-sequencing (DNAseq: coverage >3,000) and whole-genome-sequencing (WGS: coverage >35). CHIP-mutated leukocytes were visualized in human atherosclerotic plaques by mutaFISH™. Functional relevance of CHIP mutations was studied by RNAseq.

Results

DNAseq of whole blood from 540 deceased CAD patients of the Munich cardIovaScular StudIes biObaNk (MISSION) identified 253 (46.9%) CHIP mutation carriers (mean age 78.3 years). DNAseq on myocardium, atherosclerotic coronary and carotid arteries detected identical CHIP mutations in 18 out of 25 mutation carriers in tissue DNA. MutaFISH™ visualized individual macrophages carrying DNMT3A CHIP mutations in human atherosclerotic plaques. Studying monocyte-derived macrophages from Stockholm-Tartu Atherosclerosis Reverse Networks Engineering Task (STARNET; n=941) by WGS revealed CHIP mutations in 14.2% (mean age 67.1 years). RNAseq of these macrophages revealed that expression patterns in CHIP mutation carriers differed substantially from those of non-carriers. Moreover, patterns were different depending on the underlying mutations, e.g. those carrying TET2 mutations predominantly displayed upregulated inflammatory signaling whereas ASXL1 mutations showed stronger effects on metabolic pathways.

Conclusions

Deep-DNA-sequencing reveals a high prevalence of CHIP mutations in whole blood of CAD patients. CHIP-affected leukocytes invade plaques in human coronary arteries. RNAseq data obtained from macrophages of CHIP-affected patients suggest that pro-atherosclerotic signaling differs depending on the underlying mutations. Further studies are necessary to understand whether specific pathways affected by CHIP mutations may be targeted for personalized treatment.

Sirtuin 2 deficiency aggravates ageing-induced vascular remodelling in humans and mice

AIMS

The mechanisms underlying ageing-induced vascular remodelling remain unclear. This study investigates the role and underlying mechanisms of the cytoplasmic deacetylase sirtuin 2 (SIRT2) in ageing-induced vascular remodelling.

METHODS AND RESULTS

Transcriptome and quantitative real-time PCR data were used to analyse sirtuin expression. Young and old wild-type and Sirt2 knockout mice were used to explore vascular function and pathological remodelling. RNA-seq, histochemical staining, and biochemical assays were used to evaluate the effects of Sirt2 knockout on the vascular transcriptome and pathological remodelling and explore the underlying biochemical mechanisms. Among the sirtuins, SIRT2 had the highest levels in human and mouse aortas. Sirtuin 2 activity was reduced in aged aortas, and loss of SIRT2 accelerated vascular ageing. In old mice, SIRT2 deficiency aggravated ageing-induced arterial stiffness and constriction-relaxation dysfunction, accompanied by aortic remodelling (thickened vascular medial layers, breakage of elastin fibres, collagen deposition, and inflammation). Transcriptome and biochemical analyses revealed that the ageing-controlling protein p66Shc and metabolism of mitochondrial reactive oxygen species (mROS) contributed to SIRT2 function in vascular ageing. Sirtuin 2 repressed p66Shc activation and mROS production by deacetylating p66Shc at lysine 81. Elimination of reactive oxygen species by MnTBAP repressed the SIRT2 deficiency-mediated aggravation of vascular remodelling and dysfunction in angiotensin II-challenged and aged mice. The SIRT2 coexpression module in aortas was reduced with ageing across species and was a significant predictor of age-related aortic diseases in humans.

CONCLUSION

The deacetylase SIRT2 is a response to ageing that delays vascular ageing, and the cytoplasm-mitochondria axis (SIRT2-p66Shc-mROS) is important for vascular ageing. Therefore, SIRT2 may serve as a potential therapeutic target for vascular rejuvenation.

Exploring the genetic basis of coronary artery disease using functional genomics

Genome-wide Association Studies (GWAS) have identified more than 300 loci associated with coronary artery disease (CAD), defining the genetic risk map of the disease. However, the translation of the association signals into biological-pathophysiological mechanisms constitute a major challenge. Through a group of examples of studies focused on CAD, we discuss the rationale, basic principles and outcomes of the main methodologies implemented to prioritize and characterize causal variants and their target genes. Additionally, we highlight the strategies as well as the current methods that integrate association and functional genomics data to dissect the cellular specificity underlying the complexity of disease mechanisms. Despite the limitations of existing approaches, the increasing knowledge generated through functional studies helps interpret GWAS maps and opens novel avenues for the clinical usability of association data.

Multiomic Analysis and CRISPR Perturbation Screens Identify Endothelial Cell Programs and Novel Therapeutic Targets for Coronary Artery Disease

Endothelial cells (EC) are an important mediator of atherosclerosis and vascular disease. Their exposure to atherogenic risk factors such as hypertension and serum cholesterol leads to endothelial dysfunction and many disease-associated processes. Identifying which of these multiple EC functions is causally related to disease risk has been challenging. There is evidence from in vivo models and human sequencing studies that dysregulation of nitric oxide production directly affects risk of coronary artery disease. Human genetics can help prioritize the other EC functions with causal relationships because germline mutations are acquired at birth and serve as a randomized test of which pathways affect disease risk. Though several coronary artery disease risk variants have been linked to EC function, this process has been slow and laborious. Unbiased analyses of EC dysfunction using multiomic approaches promise to identify the causal genetic mechanisms responsible for vascular disease. Here, we review the data from genomic, epigenomic, and transcriptomic studies that prioritize EC-specific causal pathways. New methods that CRISPR (clustered regularly interspaced short palindromic repeats) perturbation technology with genomic, epigenomic, and transcriptomic analysis promise to speed up the characterization of disease-associated genetic variation. We summarize several recent studies in ECs which use high-throughput genetic perturbation to identify disease-relevant pathways and novel mechanisms of disease. These genetically validated pathways can accelerate the identification of drug targets for the prevention and treatment of atherosclerosis.

FHL5 Controls Vascular Disease-Associated Gene Programs in Smooth Muscle Cells

BACKGROUND

Genome-wide association studies have identified hundreds of loci associated with common vascular diseases, such as coronary artery disease, myocardial infarction, and hypertension. However, the lack of mechanistic insights for many GWAS loci limits their translation into the clinic. Among these loci with unknown functions is UFL1-four-and-a-half LIM (LIN-11, Isl-1, MEC-3) domain 5 (FHL5; chr6q16.1), which reached genome-wide significance in a recent coronary artery disease/ myocardial infarction GWAS meta-analysis. UFL1-FHL5 is also associated with several vascular diseases, consistent with the widespread pleiotropy observed for GWAS loci.

METHODS

We apply a multimodal approach leveraging statistical fine-mapping, epigenomic profiling, and ex vivo analysis of human coronary artery tissues to implicate FHL5 as the top candidate causal gene. We unravel the molecular mechanisms of the cross-phenotype genetic associations through in vitro functional analyses and epigenomic profiling experiments in coronary artery smooth muscle cells.

RESULTS

We prioritized FHL5 as the top candidate causal gene at the UFL1-FHL5 locus through expression quantitative trait locus colocalization methods. FHL5 gene expression was enriched in the smooth muscle cells and pericyte population in human artery tissues with coexpression network analyses supporting a functional role in regulating smooth muscle cell contraction. Unexpectedly, under procalcifying conditions, FHL5 overexpression promoted vascular calcification and dysregulated processes related to extracellular matrix organization and calcium handling. Lastly, by mapping FHL5 binding sites and inferring FHL5 target gene function using artery tissue gene regulatory network analyses, we highlight regulatory interactions between FHL5 and downstream coronary artery disease/myocardial infarction loci, such as FOXL1 and FN1 that have roles in vascular remodeling.

CONCLUSIONS

Taken together, these studies provide mechanistic insights into the pleiotropic genetic associations of UFL1-FHL5. We show that FHL5 mediates vascular disease risk through transcriptional regulation of downstream vascular remodeling gene programs. These transacting mechanisms may explain a portion of the heritable risk for complex vascular diseases.

Multi-omic molecular profiling and network biology for precision anaesthesiology: a narrative review

Technological advancement, data democratisation, and decreasing costs have led to a revolution in molecular biology in which the entire set of DNA, RNA, proteins, and various other molecules - the 'multi-omic' profile - can be measured in humans. Sequencing 1 million bases of human DNA now costs US$0.01, and emerging technologies soon promise to reduce the cost of sequencing the whole genome to US$100. These trends have made it feasible to sample the multi-omic profile of millions of people, much of which is publicly available for medical research. Can anaesthesiologists use these data to improve patient care? This narrative review brings together a rapidly growing literature in multi-omic profiling across numerous fields that points to the future of precision anaesthesiology. Here, we discuss how DNA, RNA, proteins, and other molecules interact in molecular networks that can be used for preoperative risk stratification, intraoperative optimisation, and postoperative monitoring. This literature provides evidence for four fundamental insights: (1) Clinically similar patients have different molecular profiles and, as a consequence, different outcomes. (2) Vast, publicly available, and rapidly growing molecular datasets have been generated in chronic disease patients and can be repurposed to estimate perioperative risk. (3) Multi-omic networks are altered in the perioperative period and influence postoperative outcomes. (4) Multi-omic networks can serve as empirical, molecular measurements of a successful postoperative course. With this burgeoning universe of molecular data, the anaesthesiologist-of-the-future will tailor their clinical management to an individual's multi-omic profile to optimise postoperative outcomes and long-term health.

The potential of integrating human and mouse discovery platforms to advance our understanding of cardiometabolic diseases

Cardiometabolic diseases encompass a range of interrelated conditions that arise from underlying metabolic perturbations precipitated by genetic, environmental, and lifestyle factors. While obesity, dyslipidaemia, smoking, and insulin resistance are major risk factors for cardiometabolic diseases, individuals still present in the absence of such traditional risk factors, making it difficult to determine those at greatest risk of disease. Thus, it is crucial to elucidate the genetic, environmental, and molecular underpinnings to better understand, diagnose, and treat cardiometabolic diseases. Much of this information can be garnered using systems genetics, which takes population-based approaches to investigate how genetic variance contributes to complex traits. Despite the important advances made by human genome-wide association studies (GWAS) in this space, corroboration of these findings has been hampered by limitations including the inability to control environmental influence, limited access to pertinent metabolic tissues, and often, poor classification of diseases or phenotypes. A complementary approach to human GWAS is the utilisation of model systems such as genetically diverse mouse panels to study natural genetic and phenotypic variation in a controlled environment. Here, we review mouse genetic reference panels and the opportunities they provide for the study of cardiometabolic diseases and related traits. We discuss how the post-GWAS era has prompted a shift in focus from discovery of novel genetic variants to understanding gene function. Finally, we highlight key advantages and challenges of integrating complementary genetic and multi-omics data from human and mouse populations to advance biological discovery.

RNF130 Regulates LDLR Availability and Plasma LDL Cholesterol Levels

BACKGROUND

Removal of circulating plasma low-density lipoprotein cholesterol (LDL-C) by the liver relies on efficient endocytosis and intracellular vesicle trafficking. Increasing the availability of hepatic LDL receptors (LDLRs) remains a major clinical target for reducing LDL-C levels. Here, we describe a novel role for RNF130 (ring finger containing protein 130) in regulating plasma membrane availability of LDLR.

METHODS

We performed a combination of gain-of-function and loss-of-function experiments to determine the effect of RNF130 on LDL-C and LDLR recycling. We overexpressed RNF130 and a nonfunctional mutant RNF130 in vivo and measured plasma LDL-C and hepatic LDLR protein levels. We performed in vitro ubiquitination assays and immunohistochemical staining to measure levels and cellular distribution of LDLR. We supplement these experiments with 3 separate in vivo models of RNF130 loss-of-function where we disrupted Rnf130 using either ASO (antisense oligonucleotides), germline deletion, or AAV CRISPR (adeno-associated virus clustered regularly interspaced short palindromic repeats) and measured hepatic LDLR and plasma LDL-C.

RESULTS

We demonstrate that RNF130 is an E3 ubiquitin ligase that ubiquitinates LDLR resulting in redistribution of the receptor away from the plasma membrane. Overexpression of RNF130 decreases hepatic LDLR and increases plasma LDL-C levels. Further, in vitro ubiquitination assays demonstrate RNF130-dependent regulation of LDLR abundance at the plasma membrane. Finally, in vivo disruption of Rnf130 using ASO, germline deletion, or AAV CRISPR results in increased hepatic LDLR abundance and availability and decreased plasma LDL-C levels.

CONCLUSIONS

Our studies identify RNF130 as a novel posttranslational regulator of LDL-C levels via modulation of LDLR availability, thus providing important insight into the complex regulation of hepatic LDLR protein levels.

Integrative multi-ancestry genetic analysis of gene regulation in coronary arteries prioritizes disease risk loci

Genome-wide association studies (GWAS) have identified hundreds of genetic risk loci for coronary artery disease (CAD). However, non-European populations are underrepresented in GWAS and the causal gene-regulatory mechanisms of these risk loci during atherosclerosis remain unclear. We incorporated local ancestry and haplotype information to identify quantitative trait loci (QTL) for gene expression and splicing in coronary arteries obtained from 138 ancestrally diverse Americans. Of 2,132 eQTL-associated genes (eGenes), 47% were previously unreported in coronary arteries and 19% exhibited cell-type-specific expression. Colocalization analysis with GWAS identified subgroups of eGenes unique to CAD and blood pressure. Fine-mapping highlighted additional eGenes of interest, including TBX20 and IL5 . Splicing (s)QTLs for 1,690 genes were also identified, among which TOR1AIP1 and ULK3 sQTLs demonstrated the importance of evaluating splicing events to accurately identify disease-relevant gene expression. Our work provides the first human coronary artery eQTL resource from a patient sample and exemplifies the necessity of diverse study populations and multi-omic approaches to characterize gene regulation in critical disease processes.

Study Design Overview

Roles of Skeletal Muscle in Development: A Bioinformatics and Systems Biology Overview

The ability to assess various cellular events consequent to perturbations, such as genetic mutations, disease states and therapies, has been recently revolutionized by technological advances in multiple "omics" fields. The resulting deluge of information has enabled and necessitated the development of tools required to both process and interpret the data. While of tremendous value to basic researchers, the amount and complexity of the data has made it extremely difficult to manually draw inference and identify factors key to the study objectives. The challenges of data reduction and interpretation are being met by the development of increasingly complex tools that integrate disparate knowledge bases and synthesize coherent models based on current biological understanding. This chapter presents an example of how genomics data can be integrated with biological network analyses to gain further insight into the developmental consequences of genetic perturbations. State of the art methods for conducting similar studies are discussed along with modern methods used to analyze and interpret the data.

Epicardial fat volume evaluated with multidetector computed tomography and other risk factors for prevalence of three-vessel coronary lesions

PURPOSE

To retrospectively investigate the epicardial fat volume with multidetector computed tomography (MDCT) and other risk factors for the prevalence of three-vessel coronary lesion.

MATERIALS AND METHODS

MDCT was performed on 424 subjects with or without three-vessel coronary lesion. Blood was tested for triglyceride, high-density lipoprotein (HDL), low-density lipoprotein (LDL), apolipoprotein A (ApoA), apolipoprotein B (ApoB), alanine aminotransferase (ALT), aspartate aminotransferase (AST), lipoprotein a, and fasting blood glucose.

RESULTS

Among all the subjects, a significant (P < 0.05) negative linear correlation existed between age and ALT or ALT/AST. The epicardial fat had a significant (P < 0.05) negative linear correlation with HDL and Apo A but a positive correlation with age and ApoB/ApoA. The epicardial fat volume and the fasting blood glucose were significantly (P = 0.001) greater in the patients than in the control group, whereas HDL and Apo A were both significantly (P < 0.0001) smaller in the patients than in the control groups. A significant prediction value (P < 0.05) existed in age increase, male gender, epicardial fat increase, low HDL, high LDL, and elevated fasting blood glucose.

CONCLUSION

Three-vessel coronary lesions are more prevalent in subjects with greater volume of epicardial fat and in male gender.

Cell-type-specific co-expression inference from single cell RNA-sequencing data

The inference of gene co-expressions from microarray and RNA-sequencing data has led to rich insights on biological processes and disease mechanisms. However, the bulk samples analyzed in most studies are a mixture of different cell types. As a result, the inferred co-expressions are confounded by varying cell type compositions across samples and only offer an aggregated view of gene regulations that may be distinct across different cell types. The advancement of single cell RNA-sequencing (scRNA-seq) technology has enabled the direct inference of co-expressions in specific cell types, facilitating our understanding of cell-type-specific biological functions. However, the high sequencing depth variations and measurement errors in scRNA-seq data present significant challenges in inferring cell-type-specific gene co-expressions, and these issues have not been adequately addressed in the existing methods. We propose a statistical approach, CS-CORE, for estimating and testing cell-type-specific co-expressions, built on a general expression-measurement model that explicitly accounts for sequencing depth variations and measurement errors in the observed single cell data. Systematic evaluations show that most existing methods suffer from inflated false positives and biased co-expression estimates and clustering analysis, whereas CS-CORE has appropriate false positive control, unbiased co-expression estimates, good statistical power and satisfactory performance in downstream co-expression analysis. When applied to analyze scRNA-seq data from postmortem brain samples from Alzheimer’s disease patients and controls and blood samples from COVID-19 patients and controls, CS-CORE identified cell-type-specific co-expressions and differential co-expressions that were more reproducible and/or more enriched for relevant biological pathways than those inferred from other methods.

Cell-type-specific co-expression inference from single cell RNA-sequencing data

The inference of gene co-expressions from microarray and RNA-sequencing data has led to rich insights on biological processes and disease mechanisms. However, the bulk samples analyzed in most studies are a mixture of different cell types. As a result, the inferred co-expressions are confounded by varying cell type compositions across samples and only offer an aggregated view of gene regulations that may be distinct across different cell types. The advancement of single cell RNA-sequencing (scRNA-seq) technology has enabled the direct inference of co-expressions in specific cell types, facilitating our understanding of cell-type-specific biological functions. However, the high sequencing depth variations and measurement errors in scRNA-seq data present significant challenges in inferring cell-type-specific gene co-expressions, and these issues have not been adequately addressed in the existing methods. We propose a statistical approach, CS-CORE, for estimating and testing cell-type-specific co-expressions, built on a general expression-measurement model that explicitly accounts for sequencing depth variations and measurement errors in the observed single cell data. Systematic evaluations show that most existing methods suffer from inflated false positives and biased co-expression estimates and clustering analysis, whereas CS-CORE has appropriate false positive control, unbiased co-expression estimates, good statistical power and satisfactory performance in downstream co-expression analysis. When applied to analyze scRNA-seq data from postmortem brain samples from Alzheimer’s disease patients and controls and blood samples from COVID-19 patients and controls, CS-CORE identified cell-type-specific co-expressions and differential co-expressions that were more reproducible and/or more enriched for relevant biological pathways than those inferred from other methods.

Loss of soluble guanylyl cyclase in platelets contributes to atherosclerotic plaque formation and vascular inflammation

Variants in genes encoding the soluble guanylyl cyclase (sGC) in platelets are associated with coronary artery disease (CAD) risk. Here, by using histology, flow cytometry and intravital microscopy, we show that functional loss of sGC in platelets of atherosclerosis-prone Ldlr^-/- mice contributes to atherosclerotic plaque formation, particularly via increasing in vivo leukocyte adhesion to atherosclerotic lesions. In vitro experiments revealed that supernatant from activated platelets lacking sGC promotes leukocyte adhesion to endothelial cells (ECs) by activating ECs. Profiling of platelet-released cytokines indicated that reduced platelet angiopoietin-1 release by sGC-depleted platelets, which was validated in isolated human platelets from carriers of GUCY1A1 risk alleles, enhances leukocyte adhesion to ECs. I mp or ta ntly, p ha rm ac ol ogical sGC stimulation increased platelet angiopoietin-1 release in vitro and reduced leukocyte recruitment and atherosclerotic plaque formation in atherosclerosis-prone Ldlr^-/- mice. Therefore, pharmacological sGC stimulation might represent a potential therapeutic strategy to prevent and treat CAD.

Within- and cross-tissue gene regulations were disrupted by PM2.5 nitrate exposure and associated with respiratory functions

BACKGROUND

Pathogenesis of complex diseases often involves multiple organs/tissue-types. To date, the PM_2.5 exposure's toxic effects and induced disease risks were not studied at multi-tissue level.

METHODS

C57BL/6 mice (n = 40) were exposed to PM_2.5 NO₃^- and clean air, respectively, and afterwards assessed respiratory functions and transcriptome in relevant tissues: blood and lung. We constructed within- and cross-tissue gene regulation networks and identified network modules associated with exposure and respiratory functions.

RESULTS

PM_2.5 NO₃^- exposure elevated naïve B cells proportion in blood (p = 0.0028). Among the 6000 highest expressed genes in blood, 18.8 % (1133 genes) were altered by exposure at p ≤ 0.05 level, among which 763 genes were also associated with respiratory function (enrichment folds = 7.63, p = 2.7E-189). The exposure disrupted blood genes were primarily in the immunoregulation pathways. Both within- and cross-tissue gene network modules were perturbed by exposure and associated with respiratory function. An immunodeficiency related cross-tissue module of 555 genes was affected by exposure (p = 0.0023) and strongly correlated with FEV_0.05/FVC (r = 0.61 and p = 3E-5).

CONCLUSIONS

This study aims to fill in a major knowledge gap and investigated the effect of PM_2.5 exposure simultaneously in multiple tissues. We provided novel evidence that PM_2.5 NO₃^- exposure profoundly perturbed within- and cross-tissue gene regulations, and highlighted their roles in the etiology of respiratory decline.

How dysregulation of the immune system promotes diabetes mellitus and cardiovascular risk complications

Diabetes mellitus (DM) is a chronic metabolic disorder characterized by persistent hyperglycemia due to insulin resistance or failure to produce insulin. Patients with DM develop microvascular complications that include chronic kidney disease and retinopathy, and macrovascular complications that mainly consist in an accelerated and more severe atherosclerosis compared to the general population, increasing the risk of cardiovascular (CV) events, such as stroke or myocardial infarction by 2- to 4-fold. DM is commonly associated with a low-grade chronic inflammation that is a known causal factor in its development and its complications. Moreover, it is now well-established that inflammation and immune cells play a major role in both atherosclerosis genesis and progression, as well as in CV event occurrence. In this review, after a brief presentation of DM physiopathology and its macrovascular complications, we will describe the immune system dysregulation present in patients with type 1 or type 2 diabetes and discuss its role in DM cardiovascular complications development. More specifically, we will review the metabolic changes and aberrant activation that occur in the immune cells driving the chronic inflammation through cytokine and chemokine secretion, thus promoting atherosclerosis onset and progression in a DM context. Finally, we will discuss how genetics and recent systemic approaches bring new insights into the mechanisms behind these inflammatory dysregulations and pave the way toward precision medicine.

Identification of the Transcription Factor ATF3 as a Direct and Indirect Regulator of the LDLR

Coronary artery disease (CAD) is a complex, multifactorial disease caused, in particular, by inflammation and cholesterol metabolism. At the molecular level, the role of tissue-specific signaling pathways leading to CAD is still largely unexplored. This study relied on two main resources: (1) genes with impact on atherosclerosis/CAD, and (2) liver-specific transcriptome analyses from human and mouse studies. The transcription factor activating transcription factor 3 (ATF3) was identified as a key regulator of a liver network relevant to atherosclerosis and linked to inflammation and cholesterol metabolism. ATF3 was predicted to be a direct and indirect (via MAF BZIP Transcription Factor F (MAFF)) regulator of low-density lipoprotein receptor (LDLR). Chromatin immunoprecipitation DNA sequencing (ChIP-seq) data from human liver cells revealed an ATF3 binding motif in the promoter regions of MAFF and LDLR. siRNA knockdown of ATF3 in human Hep3B liver cells significantly upregulated LDLR expression (p < 0.01). Inflammation induced by lipopolysaccharide (LPS) stimulation resulted in significant upregulation of ATF3 (p < 0.01) and subsequent downregulation of LDLR (p < 0.001). Liver-specific expression data from human CAD patients undergoing coronary artery bypass grafting (CABG) surgery (STARNET) and mouse models (HMDP) confirmed the regulatory role of ATF3 in the homeostasis of cholesterol metabolism. This study suggests that ATF3 might be a promising treatment candidate for lowering LDL cholesterol and reducing cardiovascular risk.

The HDAC9-associated risk locus promotes coronary artery disease by governing TWIST1

Genome wide association studies (GWAS) have identified thousands of single nucleotide polymorphisms (SNPs) associated with the risk of common disorders. However, since the large majority of these risk SNPs reside outside gene-coding regions, GWAS generally provide no information about causal mechanisms regarding the specific gene(s) that are affected or the tissue(s) in which these candidate gene(s) exert their effect. The ‘gold standard’ method for understanding causal genes and their mechanisms of action are laborious basic science studies often involving sophisticated knockin or knockout mouse lines, however, these types of studies are impractical as a high-throughput means to understand the many risk variants that cause complex diseases like coronary artery disease (CAD). As a solution, we developed a streamlined, data-driven informatics pipeline to gain mechanistic insights on complex genetic loci. The pipeline begins by understanding the SNPs in a given locus in terms of their relative location and linkage disequilibrium relationships, and then identifies nearby expression quantitative trait loci (eQTLs) to determine their relative independence and the likely tissues that mediate their disease-causal effects. The pipeline then seeks to understand associations with other disease-relevant genes, disease sub-phenotypes, potential causality (Mendelian randomization), and the regulatory and functional involvement of these genes in gene regulatory co-expression networks (GRNs). Here, we applied this pipeline to understand a cluster of SNPs associated with CAD within and immediately adjacent to the gene encoding HDAC9. Our pipeline demonstrated, and validated, that this locus is causal for CAD by modulation of TWIST1 expression levels in the arterial wall, and by also governing a GRN related to metabolic function in skeletal muscle. Our results reconciled numerous prior studies, and also provided clear evidence that this locus does not govern HDAC9 expression, structure or function. This pipeline should be considered as a powerful and efficient way to understand GWAS risk loci in a manner that better reflects the highly complex nature of genetic risk associated with common disorders.

Connections for Matters of the Heart: Network Medicine in Cardiovascular Diseases

Cardiovascular diseases (CVD) are diverse disorders affecting the heart and vasculature in millions of people worldwide. Like other fields, CVD research has benefitted from the deluge of multiomics biomedical data. Current CVD research focuses on disease etiologies and mechanisms, identifying disease biomarkers, developing appropriate therapies and drugs, and stratifying patients into correct disease endotypes. Systems biology offers an alternative to traditional reductionist approaches and provides impetus for a comprehensive outlook toward diseases. As a focus area, network medicine specifically aids the translational aspect of in silico research. This review discusses the approach of network medicine and its application to CVD research.

Atherosclerosis: Recent developments

Atherosclerosis is an inflammatory disease of the large arteries that is the major cause of cardiovascular disease (CVD) and stroke. Here, we review the current understanding of the molecular, cellular, genetic, and environmental contributions to atherosclerosis, from both individual pathway and systems perspectives. We place an emphasis on recent developments, some of which have yielded unexpected biology, including previously unknown heterogeneity of inflammatory and smooth muscle cells in atherosclerotic lesions, roles for senescence and clonal hematopoiesis, and links to the gut microbiome.

The resolvin D2 - GPR18 axis is expressed in human coronary atherosclerosis and transduces atheroprotection in apolipoprotein E deficient mice

Chronic inflammation in atherosclerosis reflects a failure in the resolution of inflammation. Pro-resolving lipid mediators derived from omega-3 fatty acids reduce the development of atherosclerosis in murine models. The aim of the present study was to decipher the role of the specialized proresolving mediator (SPM) resolvin D2 (RvD2) in atherosclerosis and its signaling through the G-protein coupled receptor (GPR) 18. The ligand and receptor were detected in human coronary arteries in relation to the presence of atherosclerotic lesions and its cellular components. Importantly, RvD2 levels were significantly higher in atherosclerotic compared with healthy human coronary arteries. Furthermore, apolipoprotein E (ApoE) deficient hyperlipidemic mice were treated with either RvD2 or vehicle in the absence and presence of the GPR18 antagonist O-1918. RvD2 significantly reduced atherosclerosis, necrotic core, and pro-inflammatory macrophage marker expression. RvD2 in addition enhanced macrophage phagocytosis. The beneficial effects of RvD2 were not observed in the presence of O-1918. Taken together, these results provide evidence of atheroprotective pro-resolving signalling through the RvD2-GPR18 axis.

eQTLs as causal instruments for the reconstruction of hormone linked gene networks

Hormones act within in highly dynamic systems and much of the phenotypic response to variation in hormone levels is mediated by changes in gene expression. The increase in the number and power of large genetic association studies has led to the identification of hormone linked genetic variants. However, the biological mechanisms underpinning the majority of these loci are poorly understood. The advent of affordable, high throughput next generation sequencing and readily available transcriptomic databases has shown that many of these genetic variants also associate with variation in gene expression levels as expression Quantitative Trait Loci (eQTLs). In addition to further dissecting complex genetic variation, eQTLs have been applied as tools for causal inference. Many hormone networks are driven by transcription factors, and many of these genes can be linked to eQTLs. In this mini-review, we demonstrate how causal inference and gene networks can be used to describe the impact of hormone linked genetic variation upon the transcriptome within an endocrinology context.