Novel Pan-ERR Agonists Ameliorate Heart Failure Through Enhancing Cardiac Fatty Acid Metabolism and Mitochondrial Function.

BACKGROUND: Cardiac metabolic dysfunction is a hallmark of heart failure (HF). Estrogen-related receptors ERRα and ERRγ are essential regulators of cardiac metabolism. Therefore, activation of ERR could be a potential therapeutic intervention for HF. However, in vivo studies demonstrating the potential usefulness of ERR agonist for HF treatment are lacking, because compounds with pharmacokinetics appropriate for in vivo use have not been available. METHODS: Using a structure-based design approach, we designed and synthesized 2 structurally distinct pan-ERR agonists, SLU-PP-332 and SLU-PP-915. We investigated the effect of ERR agonist on cardiac function in a pressure overload-induced HF model in vivo. We conducted comprehensive functional, multi-omics (RNA sequencing and metabolomics studies), and genetic dependency studies both in vivo and in vitro to dissect the molecular mechanism, ERR isoform dependency, and target specificity. RESULTS: Both SLU-PP-332 and SLU-PP-915 significantly improved ejection fraction, ameliorated fibrosis, and increased survival associated with pressure overload-induced HF without affecting cardiac hypertrophy. A broad spectrum of metabolic genes was transcriptionally activated by ERR agonists, particularly genes involved in fatty acid metabolism and mitochondrial function. Metabolomics analysis showed substantial normalization of metabolic profiles in fatty acid/lipid and tricarboxylic acid/oxidative phosphorylation metabolites in the mouse heart with 6-week pressure overload. ERR agonists increase mitochondria oxidative capacity and fatty acid use in vitro and in vivo. Using both in vitro and in vivo genetic dependency experiments, we show that ERRγ is the main mediator of ERR agonism-induced transcriptional regulation and cardioprotection and definitively demonstrated target specificity. ERR agonism also led to downregulation of cell cycle and development pathways, which was partially mediated by E2F1 in cardiomyocytes. CONCLUSIONS: ERR agonists maintain oxidative metabolism, which confers cardiac protection against pressure overload-induced HF in vivo. Our results provide direct pharmacologic evidence supporting the further development of ERR agonists as novel HF therapeutics.

Adipose tissue eQTL meta-analysis reveals the contribution of allelic heterogeneity to gene expression regulation and cardiometabolic traits.

Complete characterization of the genetic effects on gene expression is needed to elucidate tissue biology and the etiology of complex traits. Here, we analyzed 2,344 subcutaneous adipose tissue samples and identified 34K conditionally distinct expression quantitative trait locus (eQTL) signals in 18K genes. Over half of eQTL genes exhibited at least two eQTL signals. Compared to primary signals, non-primary signals had lower effect sizes, lower minor allele frequencies, and less promoter enrichment; they corresponded to genes with higher heritability and higher tolerance for loss of function. Colocalization of eQTL with conditionally distinct genome-wide association study signals for 28 cardiometabolic traits identified 3,605 eQTL signals for 1,861 genes. Inclusion of non-primary eQTL signals increased colocalized signals by 46%. Among 30 genes with ≥2 pairs of colocalized signals, 21 showed a mediating gene dosage effect on the trait. Thus, expanded eQTL identification reveals more mechanisms underlying complex traits and improves understanding of the complexity of gene expression regulation.

Ventilatory efficiency measured during sub-maximal cardiopulmonary exercise testing predicts postoperative outcomes following heart transplantation.

Cardiopulmonary exercise testing is commonly used to evaluate patients for heart transplantation. We assessed the utility of ventilatory efficiency (VE/VCO2 ) to predict perioperative outcomes following heart transplantation. We retrospectively reviewed all patients undergoing cardiopulmonary exercise testing prior to heart transplantation at our center. Spearman's coefficient showed a correlation between VE/VCO2 and ICU free days in the first 30-days post-transplant (R = -.37, p < .01). A VE /VCO2 cut-off >35 was associated with significantly lower median ICU-free days (23.0 vs. 27 days; p < .01) and a higher likelihood of postoperative morbidity (OR = 5.64, 95% CI = 1.75-18.16; p < .01). Multiple regression analysis controlling for peak oxygen consumption and right heart catheter parameters showed VE/VCO2 >35 is independently associated with lower ICU-free days (p < .01) and postoperative morbidity (p = .02). Peak oxygen consumption <15 mL/min/kg was not associated with higher ICU or hospital-free days. VE/VCO2 >35 independently predicts early postoperative morbidity in patients undergoing heart transplantation.

The Type 2 Diabetes Knowledge Portal: An open access genetic resource dedicated to type 2 diabetes and related traits.

Associations between human genetic variation and clinical phenotypes have become a foundation of biomedical research. Most repositories of these data seek to be disease-agnostic and therefore lack disease-focused views. The Type 2 Diabetes Knowledge Portal (T2DKP) is a public resource of genetic datasets and genomic annotations dedicated to type 2 diabetes (T2D) and related traits. Here, we seek to make the T2DKP more accessible to prospective users and more useful to existing users. First, we evaluate the T2DKP's comprehensiveness by comparing its datasets with those of other repositories. Second, we describe how researchers unfamiliar with human genetic data can begin using and correctly interpreting them via the T2DKP. Third, we describe how existing users can extend their current workflows to use the full suite of tools offered by the T2DKP. We finally discuss the lessons offered by the T2DKP toward the goal of democratizing access to complex disease genetic results.

Synthetic ERRα/ß/γ Agonist Induces an ERRα-Dependent Acute Aerobic Exercise Response and Enhances Exercise Capacity.

Repetitive physical exercise induces physiological adaptations in skeletal muscle that improves exercise performance and is effective for the prevention and treatment of several diseases. Genetic evidence indicates that the orphan nuclear receptors estrogen receptor-related receptors (ERRs) play an important role in skeletal muscle exercise capacity. Three ERR subtypes exist (ERRα, ß, and γ), and although ERRß/γ agonists have been designed, there have been significant difficulties in designing compounds with ERRα agonist activity. Additionally, there are limited synthetic agonists that can be used to target ERRs in vivo. Here, we report the identification of a synthetic ERR pan agonist, SLU-PP-332, that targets all three ERRs but has the highest potency for ERRα. Additionally, SLU-PP-332 has sufficient pharmacokinetic properties to be used as an in vivo chemical tool. SLU-PP-332 increases mitochondrial function and cellular respiration in a skeletal muscle cell line. When administered to mice, SLU-PP-332 increased the type IIa oxidative skeletal muscle fibers and enhanced exercise endurance. We also observed that SLU-PP-332 induced an ERRα-specific acute aerobic exercise genetic program, and the ERRα activation was critical for enhancing exercise endurance in mice. These data indicate the feasibility of targeting ERRα for the development of compounds that act as exercise mimetics that may be effective in the treatment of numerous metabolic disorders and to improve muscle function in the aging.

Improved ventilatory response during exercise over time after concussion: A case report.

This case study reports the exercise ventilatory responses of a 17-year-old female who sustained a traumatic brain injury (TBI) which resulted in loss of consciousness. Subsequently, she suffered from post-concussion syndrome accompanied by orthostatic and exercise intolerance. A cardiopulmonary exercise test (CPET) was performed 2 years post-TBI. The results demonstrated significant hypoventilation with elevated PetCO2 and ventilatory equivalents, progressive desaturation, and pre-syncope symptoms limiting exercise capacity. A repeat CPET 4 years post-TBI demonstrated a similar exercise capacity limited by pre-syncope symptoms. However, there was a marked improvement in the ventilatory response, with appropriate ventilation, PetCO2, and ventilatory equivalents, and only minor desaturation near peak exercise. Hypoventilation during exercise has been reported in subjects with post-concussion syndrome; however, with time the exercise ventilatory response can potentially normalize as observed in this subject.

Exercise induced bronchodilation: a phenomenon more common, greater magnitude and more prolonged in older adults than in adolescents.

OBJECTIVE: There are few studies in clinically healthy subjects describing and quantifying exercise-induced bronchodilation (EIBd). This study aimed to describe and compare the magnitude and time course changes in post-exercise forced expired volume at the first second (FEV1) in healthy adolescents, younger adults, and older adults. METHODS: Adolescent (n = 73, aged 10-17 years), younger adult (n = 35, aged 18-25 years), and older adult (n = 25, aged 35-66 years) subjects with normal spirometry z-scores completed a maximal cardiopulmonary exercise test using the standardized exponential exercise test protocol on a cycle ergometer performed at stable temperature and humidity. Spirometry was performed pre-exercise and at 1-, 3-, 5-, and 10-minutes post-exercise to determine the percentage change in FEV1 compared to baseline. EIBd was defined as a ≥ 5% increase in post-exercise FEV1. RESULTS: Increases in FEV1 at one-minute post-exercise were observed in the adolescents (1.3%) and young adults (6.0%) with FEV1 returning to baseline after ten minutes. Compared to the adolescents, the older adults showed significantly greater and sustained increases in FEV1 at 1-, 3-, 5-, and 10-minutes post-exercise (6.4, 4.6, 4.7, and 3.8%, p < 0.05). At 1-minute post exercise a significantly greater proportion of younger adults (54%, p < 0.01) and older adults (64%, p < 0.01) demonstrated EIBd compared to the adolescent group (15%). CONCLUSION: Healthy older adults had a higher prevalence, greater magnitude and more prolonged EIBd compared to healthy adolescent and young adult subjects.

Integrating transcriptomics, metabolomics, and GWAS helps reveal molecular mechanisms for metabolite levels and disease risk.

Transcriptomics data have been integrated with genome-wide association studies (GWASs) to help understand disease/trait molecular mechanisms. The utility of metabolomics, integrated with transcriptomics and disease GWASs, to understand molecular mechanisms for metabolite levels or diseases has not been thoroughly evaluated. We performed probabilistic transcriptome-wide association and locus-level colocalization analyses to integrate transcriptomics results for 49 tissues in 706 individuals from the GTEx project, metabolomics results for 1,391 plasma metabolites in 6,136 Finnish men from the METSIM study, and GWAS results for 2,861 disease traits in 260,405 Finnish individuals from the FinnGen study. We found that genetic variants that regulate metabolite levels were more likely to influence gene expression and disease risk compared to the ones that do not. Integrating transcriptomics with metabolomics results prioritized 397 genes for 521 metabolites, including 496 previously identified gene-metabolite pairs with strong functional connections and suggested 33.3% of such gene-metabolite pairs shared the same causal variants with genetic associations of gene expression. Integrating transcriptomics and metabolomics individually with FinnGen GWAS results identified 1,597 genes for 790 disease traits. Integrating transcriptomics and metabolomics jointly with FinnGen GWAS results helped pinpoint metabolic pathways from genes to diseases. We identified putative causal effects of UGT1A1/UGT1A4 expression on gallbladder disorders through regulating plasma (E,E)-bilirubin levels, of SLC22A5 expression on nasal polyps and plasma carnitine levels through distinct pathways, and of LIPC expression on age-related macular degeneration through glycerophospholipid metabolic pathways. Our study highlights the power of integrating multiple sets of molecular traits and GWAS results to deepen understanding of disease pathophysiology.

Extent to which array genotyping and imputation with large reference panels approximate deep whole-genome sequencing.

Understanding the genetic basis of human diseases and traits is dependent on the identification and accurate genotyping of genetic variants. Deep whole-genome sequencing (WGS), the gold standard technology for SNP and indel identification and genotyping, remains very expensive for most large studies. Here, we quantify the extent to which array genotyping followed by genotype imputation can approximate WGS in studies of individuals of African, Hispanic/Latino, and European ancestry in the US and of Finnish ancestry in Finland (a population isolate). For each study, we performed genotype imputation by using the genetic variants present on the Illumina Core, OmniExpress, MEGA, and Omni 2.5M arrays with the 1000G, HRC, and TOPMed imputation reference panels. Using the Omni 2.5M array and the TOPMed panel, ≥90% of bi-allelic single-nucleotide variants (SNVs) are well imputed (r2 > 0.8) down to minor-allele frequencies (MAFs) of 0.14% in African, 0.11% in Hispanic/Latino, 0.35% in European, and 0.85% in Finnish ancestries. There was little difference in TOPMed-based imputation quality among the arrays with >700k variants. Individual-level imputation quality varied widely between and within the three US studies. Imputation quality also varied across genomic regions, producing regions where even common (MAF > 5%) variants were consistently not well imputed across ancestries. The extent to which array genotyping and imputation can approximate WGS therefore depends on reference panel, genotype array, sample ancestry, and genomic location. Imputation quality by variant or genomic region can be queried with our new tool, RsqBrowser, now deployed on the Michigan Imputation Server.

ACE2 expression in adipose tissue is associated with cardio-metabolic risk factors and cell type composition-implications for COVID-19.

BACKGROUND: COVID-19 severity varies widely. Although some demographic and cardio-metabolic factors, including age and obesity, are associated with increasing risk of severe illness, the underlying mechanism(s) are uncertain. SUBJECTS/METHODS: In a meta-analysis of three independent studies of 1471 participants in total, we investigated phenotypic and genetic factors associated with subcutaneous adipose tissue expression of Angiotensin I Converting Enzyme 2 (ACE2), measured by RNA-Seq, which acts as a receptor for SARS-CoV-2 cellular entry. RESULTS: Lower adipose tissue ACE2 expression was associated with multiple adverse cardio-metabolic health indices, including type 2 diabetes (T2D) (P = 9.14 × 10-6), obesity status (P = 4.81 × 10-5), higher serum fasting insulin (P = 5.32 × 10-4), BMI (P = 3.94 × 10-4), and lower serum HDL levels (P = 1.92 × 10-7). ACE2 expression was also associated with estimated proportions of cell types in adipose tissue: lower expression was associated with a lower proportion of microvascular endothelial cells (P = 4.25 × 10-4) and higher proportion of macrophages (P = 2.74 × 10-5). Despite an estimated heritability of 32%, we did not identify any proximal or distal expression quantitative trait loci (eQTLs) associated with adipose tissue ACE2 expression. CONCLUSIONS: Our results demonstrate that individuals with cardio-metabolic features known to increase risk of severe COVID-19 have lower background ACE2 levels in this highly relevant tissue. Reduced adipose tissue ACE2 expression may contribute to the pathophysiology of cardio-metabolic diseases, as well as the associated increased risk of severe COVID-19.

Genome-wide association studies of metabolites in Finnish men identify disease-relevant loci.

Few studies have explored the impact of rare variants (minor allele frequency < 1%) on highly heritable plasma metabolites identified in metabolomic screens. The Finnish population provides an ideal opportunity for such explorations, given the multiple bottlenecks and expansions that have shaped its history, and the enrichment for many otherwise rare alleles that has resulted. Here, we report genetic associations for 1391 plasma metabolites in 6136 men from the late-settlement region of Finland. We identify 303 novel association signals, more than one third at variants rare or enriched in Finns. Many of these signals identify genes not previously implicated in metabolite genome-wide association studies and suggest mechanisms for diseases and disease-related traits.

Emerging Role of Nuclear Receptors for the Treatment of NAFLD and NASH.

Non-alcoholic fatty liver (NAFLD) over the past years has become a metabolic pandemic linked to a collection of metabolic diseases. The nuclear receptors ERRs, REV-ERBs, RORs, FXR, PPARs, and LXR are master regulators of metabolism and liver physiology. The characterization of these nuclear receptors and their biology has promoted the development of synthetic ligands. The possibility of targeting these receptors to treat NAFLD is promising, as several compounds including Cilofexor, thiazolidinediones, and Saroglitazar are currently undergoing clinical trials. This review focuses on the latest development of the pharmacology of these metabolic nuclear receptors and how they may be utilized to treat NAFLD and subsequent comorbidities.

Quantifying tidal expiratory flow limitation using a vector-based analysis technique.

Objective. Tidal expiratory flow limitation (EFLT) is commonly identified by tidal breaths exceeding the forced vital capacity (FVC) loop. This technique, known as the Hyatt method, is limited by the difficulties in defining the FVC and tidal flow-volume (TV) loops. The vector-based analysis (VBA) technique described and piloted in this manuscript identifies and quantifies EFLTas tidal breaths that conform to the contour of the FVC loop.Approach. The FVC and TV loops are interpolated to generate uniformly spaced plots. VBA is performed to determine the smallest vector difference between each point on the FVC and TV curves, termed the flow reserve vector (FRV). From the FVC point yielding the lowest FRV, the tangential angles of the FVC and TV segments are recorded. If the TV and FVC loops become parallel, the difference between the tangential angles tends towards zero. We infer EFLTas parallel TV and FVC segments where the FRV is < 0.1 and the tangential angle is within ±18 degrees for ≥5% of TV. EFLTis quantified by the percent of TV loop fulfilling these criteria. We compared the presence and degree of EFLTat rest and during peak exercise using the Hyatt method and our VBA technique in 25 healthy subjects and 20 subjects with moderate-severe airflow obstruction.Main results. Compared to the Hyatt method, our VBA technique reported a significantly lower degree of EFLTin healthy subjects during peak exercise, and in obstructed subjects at rest and during peak exercise. In contrast to the Hyatt method, our VBA technique re-classified five subjects (one in the healthy group and four in the obstructed group) as demonstrating EFLT.Significance.Our VBA technique provides an alternative approach to determine and quantify EFLTwhich may reduce the overestimation of the degree EFLTand more accurately identify subjects experiencing EFLT.

Novel application of Poincaré analysis to detect and quantify exercise oscillatory ventilation.

Objective.Exercise oscillatory ventilation (EOV) is frequently observed in individuals with cardiac disease. Assessment of EOV relies on pattern recognition and this subjectivity and lack of quantification limits the widespread clinical use of EOV as a prognostic marker. Poincaré analysis quantifies the short (SD1) and long-term (SD2) variability of a signal and may provide an alternative means to identify and quantify unstable exercise breathing patterns. This study aimed to determine if Poincaré analysis can distinguish between the breathing patterns of healthy control subjects and individuals being assessed for heart transplantation with and without EOV.Approach.Thirty-nine subjects performed a cardiopulmonary exercise test as part of heart transplant assessment and were subjectively classified into two groups according to the presence of EOV: non-EOV (n = 19) and EOV (n = 20). The control group (n = 24) consisted of healthy adults. Poincaré analysis (SD1 and SD2) was performed for minute ventilation (V̇E) and tidal volume (VT) normalized to forced vital capacity (V̇EnandV̇Tn), and breathing frequency (BF) for breath-by-breath data over the 10-15 ml · min-1 · kg-1V̇O2range.Main results.Poincaré analysis showed similar exercise ventilatory responses between the non-EOV and control group. BF was found to discriminate between subjects with stable and unstable ventilation. BF SD1 was significantly higher in the EOV group compared to the non-EOV (7.9 versus 4.6,p < 0.01) and control (7.9 versus 4.2,p < 0.01) groups. The EOV group had significantly greater BF SD2 compared to the non-EOV (5.7 versus 3.5,p < 0.01) and control (5.7 versus 3.5,p < 0.01) groups.Significance.We demonstrated that this novel application of Poincaré analysis can objectively distinguish and quantify unstable from stable breathing patterns during exercise. In subjects being assessed for heart transplantation the presence of EOV is associated with greater BF variability. Poincaré analysis provides an objective measure to identify and quantify EOV.Summary at a glance.As EOV may indicate abnormal ventilatory control, there is a need for an objective measure to identify and quantify unstable from stable ventilation during exercise. We developed a method of quantifying BF variation by the application of Poincaré analysis and demonstrated higher than normal variability of BF in subjects being assessed for heart transplantation who demonstrated EOV.

LocusZoom.js: interactive and embeddable visualization of genetic association study results.

SUMMARY: LocusZoom.js is a JavaScript library for creating interactive web-based visualizations of genetic association study results. It can display one or more traits in the context of relevant biological data (such as gene models and other genomic annotation), and allows interactive refinement of analysis models (by selecting linkage disequilibrium reference panels, identifying sets of likely causal variants, or comparisons to the GWAS catalog). It can be embedded in web pages to enable data sharing and exploration. Views can be customized and extended to display other data types such as phenome-wide association study (PheWAS) results, chromatin co-accessibility, or eQTL measurements. A new web upload service harmonizes datasets, adds annotations, and makes it easy to explore user-provided result sets. AVAILABILITY AND IMPLEMENTATION: LocusZoom.js is open-source software under a permissive MIT license. Code and documentation are available at: https://github.com/statgen/locuszoom/. Installable packages for all versions are also distributed via NPM. Additional features are provided as standalone libraries to promote reuse. Use with your own GWAS results at https://my.locuszoom.org/. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Intestinal α1-2-Fucosylation Contributes to Obesity and Steatohepatitis in Mice.

BACKGROUND & AIMS: Fucosyltransferase 2 (Fut2)-mediated intestinal α1- 2-fucosylation is important for host-microbe interactions and has been associated with several diseases, but its role in obesity and hepatic steatohepatitis is not known. The aim of this study was to investigate the role of Fut2 in a Western-style diet-induced mouse model of obesity and steatohepatitis. METHODS: Wild-type (WT) and Fut2-deficient littermate mice were used and features of the metabolic syndrome and steatohepatitis were assessed after 20 weeks of Western diet feeding. RESULTS: Intestinal α1-2-fucosylation was suppressed in WT mice after Western diet feeding, and supplementation of α1-2-fucosylated glycans exacerbated obesity and steatohepatitis in these mice. Fut2-deficient mice were protected from Western diet-induced features of obesity and steatohepatitis despite an increased caloric intake. These mice have increased energy expenditure and thermogenesis, as evidenced by a higher core body temperature. Protection from obesity and steatohepatitis associated with Fut2 deficiency is transmissible to WT mice via microbiota exchange; phenotypic differences between Western diet-fed WT and Fut2-deficient mice were reduced with antibiotic treatment. Fut2 deficiency attenuated diet-induced bile acid accumulation by altered relative abundance of bacterial enzyme 7-α-hydroxysteroid dehydrogenases metabolizing bile acids and by increased fecal excretion of secondary bile acids. This also was associated with increased intestinal farnesoid X receptor/fibroblast growth factor 15 signaling, which inhibits hepatic synthesis of bile acids. Dietary supplementation of α1-2-fucosylated glycans abrogates the protective effects of Fut2 deficiency. CONCLUSIONS: α1-2-fucosylation is an important host-derived regulator of intestinal microbiota and plays an important role for the pathogenesis of obesity and steatohepatitis in mice.

Genetic variant effects on gene expression in human pancreatic islets and their implications for T2D.

Most signals detected by genome-wide association studies map to non-coding sequence and their tissue-specific effects influence transcriptional regulation. However, key tissues and cell-types required for functional inference are absent from large-scale resources. Here we explore the relationship between genetic variants influencing predisposition to type 2 diabetes (T2D) and related glycemic traits, and human pancreatic islet transcription using data from 420 donors. We find: (a) 7741 cis-eQTLs in islets with a replication rate across 44 GTEx tissues between 40% and 73%; (b) marked overlap between islet cis-eQTL signals and active regulatory sequences in islets, with reduced eQTL effect size observed in the stretch enhancers most strongly implicated in GWAS signal location; (c) enrichment of islet cis-eQTL signals with T2D risk variants identified in genome-wide association studies; and (d) colocalization between 47 islet cis-eQTLs and variants influencing T2D or glycemic traits, including DGKB and TCF7L2. Our findings illustrate the advantages of performing functional and regulatory studies in disease relevant tissues.

Adiponectin GWAS loci harboring extensive allelic heterogeneity exhibit distinct molecular consequences.

Loci identified in genome-wide association studies (GWAS) can include multiple distinct association signals. We sought to identify the molecular basis of multiple association signals for adiponectin, a hormone involved in glucose regulation secreted almost exclusively from adipose tissue, identified in the Metabolic Syndrome in Men (METSIM) study. With GWAS data for 9,262 men, four loci were significantly associated with adiponectin: ADIPOQ, CDH13, IRS1, and PBRM1. We performed stepwise conditional analyses to identify distinct association signals, a subset of which are also nearly independent (lead variant pairwise r2<0.01). Two loci exhibited allelic heterogeneity, ADIPOQ and CDH13. Of seven association signals at the ADIPOQ locus, two signals colocalized with adipose tissue expression quantitative trait loci (eQTLs) for three transcripts: trait-increasing alleles at one signal were associated with increased ADIPOQ and LINC02043, while trait-increasing alleles at the other signal were associated with decreased ADIPOQ-AS1. In reporter assays, adiponectin-increasing alleles at two signals showed corresponding directions of effect on transcriptional activity. Putative mechanisms for the seven ADIPOQ signals include a missense variant (ADIPOQ G90S), a splice variant, a promoter variant, and four enhancer variants. Of two association signals at the CDH13 locus, the first signal consisted of promoter variants, including the lead adipose tissue eQTL variant for CDH13, while a second signal included a distal intron 1 enhancer variant that showed ~2-fold allelic differences in transcriptional reporter activity. Fine-mapping and experimental validation demonstrated that multiple, distinct association signals at these loci can influence multiple transcripts through multiple molecular mechanisms.

ACE2 expression in adipose tissue is associated with COVID-19 cardio-metabolic risk factors and cell type composition.

COVID-19 severity has varied widely, with demographic and cardio-metabolic factors increasing risk of severe reactions to SARS-CoV-2 infection, but the underlying mechanisms for this remain uncertain. We investigated phenotypic and genetic factors associated with subcutaneous adipose tissue expression of Angiotensin I Converting Enzyme 2 ( ACE2 ), which has been shown to act as a receptor for SARS-CoV-2 cellular entry. In a meta-analysis of three independent studies including up to 1,471 participants, lower adipose tissue ACE2 expression was associated with adverse cardio-metabolic health indices including type 2 diabetes (T2D) and obesity status, higher serum fasting insulin and BMI, and lower serum HDL levels (P<5.32x10 -4 ). ACE2 expression levels were also associated with estimated proportions of cell types in adipose tissue; lower ACE2 expression was associated with a lower proportion of microvascular endothelial cells (P=4.25x10 -4 ) and higher macrophage proportion (P=2.74x10 -5 ), suggesting a link to inflammation. Despite an estimated heritability of 32%, we did not identify any proximal or distal genetic variants (eQTLs) associated with adipose tissue ACE2 expression. Our results demonstrate that at-risk individuals have lower background ACE2 levels in this highly relevant tissue. Further studies will be required to establish how this may contribute to increased COVID-19 severity.