System Genetics Including Causal Inference Identify Immune Targets for Coronary Artery Disease and the Lifespan

Deep learning-derived splenic radiomics, genomics, and coronary artery disease

Background

Despite advances in managing traditional risk factors, coronary artery disease (CAD) remains the leading cause of mortality. Circulating hematopoietic cells influence risk for CAD, but the role of a key regulating organ, spleen, is unknown. The understudied spleen is a 3-dimensional structure of the hematopoietic system optimally suited for unbiased radiologic investigations toward novel mechanistic insights.

Methods

Deep learning-based image segmentation and radiomics techniques were utilized to extract splenic radiomic features from abdominal MRIs of 42,059 UK Biobank participants. Regression analysis was used to identify splenic radiomics features associated with CAD. Genome-wide association analyses were applied to identify loci associated with these radiomics features. Overlap between loci associated with CAD and the splenic radiomics features was explored to understand the underlying genetic mechanisms of the role of the spleen in CAD.

Results

We extracted 107 splenic radiomics features from abdominal MRIs, and of these, 10 features were associated with CAD. Genome-wide association analysis of CAD-associated features identified 219 loci, including 35 previously reported CAD loci, 7 of which were not associated with conventional CAD risk factors. Notably, variants at 9p21 were associated with splenic features such as run length non-uniformity.

Conclusions

Our study, combining deep learning with genomics, presents a new framework to uncover the splenic axis of CAD. Notably, our study provides evidence for the underlying genetic connection between the spleen as a candidate causal tissue-type and CAD with insight into the mechanisms of 9p21, whose mechanism is still elusive despite its initial discovery in 2007. More broadly, our study provides a unique application of deep learning radiomics to non-invasively find associations between imaging, genetics, and clinical outcomes.

Mendelian randomization reveals interactions of the blood proteome and immunome in mitral valve prolapse

BACKGROUND

Mitral valve prolapse (MVP) is a common heart disorder characterized by an excessive production of proteoglycans and extracellular matrix in mitral valve leaflets. Large-scale genome-wide association study (GWAS) underlined that MVP is heritable. The molecular underpinnings of the disease remain largely unknown.

METHODS

We interrogated cross-modality data totaling more than 500,000 subjects including GWAS, 4809 molecules of the blood proteome, and genome-wide expression of mitral valves to identify candidate drivers of MVP. Data were investigated through Mendelian randomization, network analysis, ligand-receptor inference and digital cell quantification.

RESULTS

In this study, Mendelian randomization identify that 33 blood proteins, enriched in networks for immunity, are associated with the risk of MVP. MVP- associated blood proteins are enriched in ligands for which their cognate receptors are differentially expressed in mitral valve leaflets during MVP and enriched in cardiac endothelial cells and macrophages. MVP-associated blood proteins are involved in the renewal-polarization of macrophages and regulation of adaptive immune response. Cytokine activity profiling and digital cell quantification show in MVP a shift toward cytokine signature promoting M2 macrophage polarization. Assessment of druggability identify CSF1R, CX3CR1, CCR6, IL33, MMP8, ENPEP and angiotensin receptors as actionable targets in MVP.

CONCLUSIONS

Hence, integrative analysis identifies networks of candidate molecules and cells involved in immune control and remodeling of the extracellular matrix, which drive the risk of MVP.

The landscape of the methodology in drug repurposing using human genomic data: a systematic review

The process of drug development is expensive and time-consuming. In contrast, drug repurposing can be introduced to clinical practice more quickly and at a reduced cost. Over the last decade, there has been a significant expansion of large biobanks that link genomic data to electronic health record data, public availability of various databases containing biological and clinical information and rapid development of novel methodologies and algorithms in integrating different sources of data. This review aims to provide a thorough summary of different strategies that utilize genomic data to seek drug-repositioning opportunities. We searched MEDLINE and EMBASE databases to identify eligible studies up until 1 May 2023, with a total of 102 studies finally included after two-step parallel screening. We summarized commonly used strategies for drug repurposing, including Mendelian randomization, multi-omic-based and network-based studies and illustrated each strategy with examples, as well as the data sources implemented. By leveraging existing knowledge and infrastructure to expedite the drug discovery process and reduce costs, drug repurposing potentially identifies new therapeutic uses for approved drugs in a more efficient and targeted manner. However, technical challenges when integrating different types of data and biased or incomplete understanding of drug interactions are important hindrances that cannot be disregarded in the pursuit of identifying novel therapeutic applications. This review offers an overview of drug repurposing methodologies, providing valuable insights and guiding future directions for advancing drug repurposing studies.

Single-cell and spatial transcriptomics: Advances in heart development and disease applications

Current transcriptomics technologies, including bulk RNA-seq, single-cell RNA sequencing (scRNA-seq), single-nucleus RNA-sequencing (snRNA-seq), and spatial transcriptomics (ST), provide novel insights into the spatial and temporal dynamics of gene expression during cardiac development and disease processes. Cardiac development is a highly sophisticated process involving the regulation of numerous key genes and signaling pathways at specific anatomical sites and developmental stages. Exploring the cell biological mechanisms involved in cardiogenesis also contributes to congenital heart disease research. Meanwhile, the severity of distinct heart diseases, such as coronary heart disease, valvular disease, cardiomyopathy, and heart failure, is associated with cellular transcriptional heterogeneity and phenotypic alteration. Integrating transcriptomic technologies in the clinical diagnosis and treatment of heart diseases will aid in advancing precision medicine. In this review, we summarize applications of scRNA-seq and ST in the cardiac field, including organogenesis and clinical diseases, and provide insights into the promise of single-cell and spatial transcriptomics in translational research and precision medicine.

Enhancer promoter interactome and Mendelian randomization identify network of druggable vascular genes in coronary artery disease

Coronary artery disease (CAD) is a multifactorial disorder, which is partly heritable. Herein, we implemented a mapping of CAD-associated candidate genes by using genome-wide enhancer-promoter conformation (H3K27ac-HiChIP) and expression quantitative trait loci (eQTL). Enhancer-promoter anchor loops from human coronary artery smooth muscle cells (HCASMC) explained 22% of the heritability for CAD. 3D enhancer-promoter genome mapping of CAD-genes in HCASMC was enriched in vascular eQTL genes. By using colocalization and Mendelian randomization analyses, we identified 58 causal candidate vascular genes including some druggable targets (MAP3K11, CAMK1D, PDGFD, IPO9 and CETP). A network analysis of causal candidate genes was enriched in TGF beta and MAPK pathways. The pharmacologic inhibition of causal candidate gene MAP3K11 in vascular SMC reduced the expression of athero-relevant genes and lowered cell migration, a cardinal process in CAD. Genes connected to enhancers are enriched in vascular eQTL and druggable genes causally associated with CAD.

Positive Selection and Enhancer Evolution Shaped Lifespan and Body Mass in Great Apes

Within primates, the great apes are outliers both in terms of body size and lifespan, since they include the largest and longest-lived species in the order. Yet, the molecular bases underlying such features are poorly understood. Here, we leveraged an integrated approach to investigate multiple sources of molecular variation across primates, focusing on over 10,000 genes, including approximately 1,500 previously associated with lifespan, and additional approximately 9,000 for which an association with longevity has never been suggested. We analyzed dN/dS rates, positive selection, gene expression (RNA-seq), and gene regulation (ChIP-seq). By analyzing the correlation between dN/dS, maximum lifespan, and body mass, we identified 276 genes whose rate of evolution positively correlates with maximum lifespan in primates. Further, we identified five genes, important for tumor suppression, adaptive immunity, metastasis, and inflammation, under positive selection exclusively in the great ape lineage. RNA-seq data, generated from the liver of six species representing all the primate lineages, revealed that 8% of approximately 1,500 genes previously associated with longevity are differentially expressed in apes relative to other primates. Importantly, by integrating RNA-seq with ChIP-seq for H3K27ac (which marks active enhancers), we show that the differentially expressed longevity genes are significantly more likely than expected to be located near a novel "ape-specific" enhancer. Moreover, these particular ape-specific enhancers are enriched for young transposable elements, and specifically SINE-Vntr-Alus. In summary, we demonstrate that multiple evolutionary forces have contributed to the evolution of lifespan and body size in primates.

Detecting Rewiring Events in Protein-Protein Interaction Networks Based on Transcriptomic Data

Proteins rarely carry out their cellular functions in isolation. Instead, eukaryotic proteins engage in about six interactions with other proteins on average. The aggregated protein interactome of an organism forms a “hairy ball”-type protein-protein interaction (PPI) network. Yet, in a typical human cell, only about half of all proteins are expressed at a particular time. Hence, it has become common practice to prune the full PPI network to the subset of expressed proteins. If RNAseq data is available, one can further resolve the specific protein isoforms present in a cell or tissue. Here, we review various approaches, software tools and webservices that enable users to construct context-specific or tissue-specific PPI networks and how these are rewired between two cellular conditions. We illustrate their different functionalities on the example of the interactions involving the human TNR6 protein. In an outlook, we describe how PPI networks may be integrated with epigenetic data or with data on the activity of splicing factors.