kruX:Matrix-based non-parametric eQTL discovery

The Kruskal-Wallis test is a popular non-parametric statistical test for identifying expression quantitative trait loci (eQTLs) from genome-wide data due to its robustness against variations in the underlying genetic model and expression trait distribution, but testing billions of marker-trait combinations one-by-one can become computationally prohibitive. We developed kruX, an algorithm implemented in Matlab, Python and R that uses matrix multiplications to simultaneously calculate the Kruskal-Wallis test statistic for several millions of marker-trait combinations at once. KruX is more than ten thousand times faster than computing associations one-by-one on a typical human dataset. We used kruX and a dataset of more than 500k SNPs and 20k expression traits measured in 102 human blood samples to compare eQTLs detected by the Kruskal-Wallis test to eQTLs detected by the parametric ANOVA and linear model methods. We found that the Kruskal-Wallis test is more robust against data outliers and heterogeneous genotype group sizes and detects a higher proportion of non-linear associations, but is more conservative for calling additive linear associations. In summary, kruX enables the use of robust non-parametric methods for massive eQTL mapping without the need for a high-performance computing infrastructure.Comment: minor revision; 6 pages, 5 figures; software available at http://krux.googlecode.co

The SARS-CoV-2 receptor ACE2 is expressed in mouse pericytes but not endothelial cells : Implications for COVID-19 vascular research

Humanized mouse models and mouse-adapted SARS-CoV-2 virus are increasingly used to study COVID-19 pathogenesis, so it is impor-tant to learn where the SARS-CoV-2 receptor ACE2 is expressed. Here we mapped ACE2 expression during mouse postnatal development and in adulthood. Pericytes in the CNS, heart, and pancreas express ACE2 strongly, as do perineurial and adrenal fibroblasts, whereas endothelial cells do not at any location analyzed. In a number of other organs, pericytes do not express ACE2, including in the lung where ACE2 instead is expressed in bronchial epithelium and alveolar type II cells. The onset of ACE2 expression is organ specific: in bronchial epithelium already at birth, in brain pericytes before, andin heart pericytes after postnatal day 10.5. Establishing the vascular localization of ACE2 expression is central to correctly interpret data from modeling COVID-19 in the mouse and may shed light on the cause of vascular COVID-19 complications.Peer reviewe

Cystatin C is glucocorticoid responsive, directs recruitment of Trem2+ macrophages, and predicts failure of cancer immunotherapy

Cystatin C (CyC), a secreted cysteine protease inhibitor, has unclear biological functions. Many patients exhibit elevated plasma CyC levels, particularly during glucocorticoid (GC) treatment. This study links GCs with CyC’s systemic regulation by utilizing genome-wide association and structural equation modeling to determine CyC production genetics in the UK Biobank. Both CyC production and a polygenic score (PGS) capturing predisposition to CyC production were associated with increased all-cause and cancer-specific mortality. We found that the GC receptor directly targets CyC, leading to GC-responsive CyC secretion in macrophages and cancer cells. CyC-knockout tumors displayed significantly reduced growth and diminished recruitment of TREM2+ macrophages, which have been connected to cancer immunotherapy failure. Furthermore, the CyC-production PGS predicted checkpoint immunotherapy failure in 685 patients with metastatic cancer from combined clinical trial cohorts. In conclusion, CyC may act as a GC effector pathway via TREM2+ macrophage recruitment and may be a potential target for combination cancer immunotherapy.publishedVersio

Variation in the SERPINA6SERPINA1 locusalters morning plasma cortisol, hepatic corticosteroid binding globulin expression, gene expressionin peripheral tissues, and risk of cardiovascular disease

The stress hormone cortisol modulates fuel metabolism, cardiovascular homoeostasis, mood, inflammation and cognition. The CORtisol NETwork (CORNET) consortium previously identified a single locus associated with morning plasma cortisol. Identifying additional genetic variants that explain more of the variance in cortisol could provide new insights into cortisol biology and provide statistical power to test the causative role of cortisol in common diseases. The CORNET consortium extended its genome-wide association meta-analysis for morning plasma cortisol from 12,597 to 25,314 subjects and from ~2.2 M to ~7 M SNPs, in 17 population-based cohorts of European ancestries. We confirmed the genetic association with SERPINA6/SERPINA1. This locus contains genes encoding corticosteroid binding globulin (CBG) and α1-antitrypsin. Expression quantitative trait loci (eQTL) analyses undertaken in the STARNET cohort of 600 individuals showed that specific genetic variants within the SERPINA6/SERPINA1 locus influence expression of SERPINA6 rather than SERPINA1 in the liver. Moreover, trans-eQTL analysis demonstrated effects on adipose tissue gene expression, suggesting that variation

Variation in the SERPINA6/SERPINA1 locus alters morning plasma cortisol, hepatic corticosteroid binding globulin expression, gene expression in peripheral tissues, and risk of cardiovascular disease

The stress hormone cortisol modulates fuel metabolism, cardiovascular homoeostasis, mood, inflammation and cognition. The CORtisol NETwork (CORNET) consortium previously identified a single locus associated with morning plasma cortisol. Identifying additional genetic variants that explain more of the variance in cortisol could provide new insights into cortisol biology and provide statistical power to test the causative role of cortisol in common diseases. The CORNET consortium extended its genome-wide association meta-analysis for morning plasma cortisol from 12,597 to 25,314 subjects and from similar to 2.2 M to similar to 7 M SNPs, in 17 population-based cohorts of European ancestries. We confirmed the genetic association with SERPINA6/SERPINA1. This locus contains genes encoding corticosteroid binding globulin (CBG) and alpha 1-antitrypsin. Expression quantitative trait loci (eQTL) analyses undertaken in the STARNET cohort of 600 individuals showed that specific genetic variants within the SERPINA6/SERPINA1 locus influence expression of SERPINA6 rather than SERPINA1 in the liver. Moreover, trans-eQTL analysis demonstrated effects on adipose tissue gene expression, suggesting that variations in CBG levels have an effect on delivery of cortisol to peripheral tissues. Two-sample Mendelian randomisation analyses provided evidence that each genetically-determined standard deviation (SD) increase in morning plasma cortisol was associated with increased odds of chronic ischaemic heart disease (0.32, 95% CI 0.06-0.59) and myocardial infarction (0.21, 95% CI 0.00-0.43) in UK Biobank and similarly in CARDIoGRAMplusC4D. These findings reveal a causative pathway for CBG in determining cortisol action in peripheral tissues and thereby contributing to the aetiology of cardiovascular disease.</p

Reverse Engineering of Gene Networks with LASSO and Nonlinear Basis Functions

The quest to determine cause from effect is often referred to as reverse engineering in the context of cellular networks. Here we propose and evaluate an algorithm for reverse engineering a gene regulatory network from time-series kind steady-state data. Our algorithmic pipeline, which is rather standard in its parts but not in its integrative composition, combines ordinary differential equations, parameter estimations by least angle regression, and cross-validation procedures for determining the in-degrees and selection of nonlinear transfer functions. The result of the algorithm is a complete directed net-work, in which each edge has been assigned a score front it bootstrap procedure. To evaluate the performance, we submitted the outcome of the algorithm to the reverse engineering assessment competition DREAM2, where we used the data corresponding to the InSillico1 and InSilico2 networks as input. Our algorithm outperformed all other algorithms when inferring one of the directed gene-to-gene networks.This is the authors’ version of the following article:Mika Gustafsson, Michael Hörnquist, Jesper Lundstrom, Johan Bjorkegren and Jesper Tegnér, Reverse Engineering of Gene Networks with LASSO and Nonlinear Basis Functions, 2009, Annals of the New York Academy of Sciences, Volume 1158 Issue, The Challenges of Systems Biology Community Efforts to Harness Biological Complexity, 265-275.which has been published in final form at: http://dx.doi.org/10.1111/j.1749-6632.2008.03764.xCopyright: Blackwell Publishing Ltdhttp://www.blackwellpublishing.com

Genetic regulation of the placental transcriptome underlies birth weight and risk of childhood obesity.

GWAS identified variants associated with birth weight (BW), childhood obesity (CO) and childhood BMI (CBMI), and placenta is a critical organ for fetal development and postnatal health. We examined the role of placental transcriptome and eQTLs in mediating the genetic causes for BW, CO and CBMI, and applied integrative analysis (Colocalization and MetaXcan). GWAS loci associated with BW, CO, and CBMI were substantially enriched for placenta eQTLs (6.76, 4.83 and 2.26 folds, respectively). Importantly, compared to eQTLs of adult tissues, only placental eQTLs contribute significantly to both anthropometry outcomes at birth (BW) and childhood phenotypes (CO/CBMI). Eight, six and one transcripts colocalized with BW, CO and CBMI risk loci, respectively. Our study reveals that placental transcription in utero likely plays a key role in determining postnatal body size, and as such may hold new possibilities for therapeutic interventions to prevent childhood obesity

A single-cell transcriptomic inventory of murine smooth muscle cells

Smooth muscle cells (SMCs) execute important physiological functions in numerous vital organ systems, including the vascular, gastrointestinal, respiratory, and urogenital tracts. SMC differ morphologically and functionally at these different anatomical locations, but the molecular underpinnings of the differences remain poorly understood. Here, using deep single-cell RNA sequencing combined with in situ gene and pro-tein expression analysis in four murine organs-heart, aorta, lung, and colon-we identify a molecular basis for high-level differences among vascular, visceral, and airway SMC, as well as more subtle differences between, for example, SMC in elastic and muscular arteries and zonation of elastic artery SMC along the direction of blood flow. Arterial SMC exhibit extensive organotypic heterogeneity, whereas venous SMC are similar across organs. We further identify a specific SMC subtype within the pulmonary vasculature. This comparative SMC cross-organ resource offers insight into SMC subtypes and their specific functions

Single-cell analysis uncovers fibroblast heterogeneity and criteria for fibroblast and mural cell identification and discrimination

Many important cell types in adult vertebrates have a mesenchymal origin, including fibroblasts and vascular mural cells. Although their biological importance is undisputed, the level of mesenchymal cell heterogeneity within and between organs, while appreciated, has not been analyzed in detail. Here, we compare single-cell transcriptional profiles of fibroblasts and vascular mural cells across four murine muscular organs: heart, skeletal muscle, intestine and bladder. We reveal gene expression signatures that demarcate fibroblasts from mural cells and provide molecular signatures for cell subtype identification. We observe striking inter- and intra-organ heterogeneity amongst the fibroblasts, primarily reflecting differences in the expression of extracellular matrix components. Fibroblast subtypes localize to discrete anatomical positions offering novel predictions about physiological function(s) and regulatory signaling circuits. Our data shed new light on the diversity of poorly defined classes of cells and provide a foundation for improved understanding of their roles in physiological and pathological processes. To define and distinguish fibroblasts from vascular mural cells have remained challenging. Here, using single-cell RNA sequencing and tissue imaging, the authors provide a molecular basis for cell type classification and reveal inter- and intra-organ diversity of these cell types.Correction in: NATURE COMMUNICATIONS, Volume: 11, Issue: 1, Article Number: 4493, DOI: 10.1038/s41467-020-18511-8</p

Integrative functional genomics identifies regulatory mechanisms at coronary artery disease loci

Coronary artery disease (CAD) is the leading cause of mortality and morbidity, driven by both genetic and environmental risk factors. Meta-analyses of genome-wide association studies have identified 4150 loci associated with CAD and myocardial infarction susceptibility in humans. A majority of these variants reside in non-coding regions and are co-inherited with hundreds of candidate regulatory variants, presenting a challenge to elucidate their functions. Herein, we use integrative genomic, epigenomic and transcriptomic profiling of perturbed human coronary artery smooth muscle cells and tissues to begin to identify causal regulatory variation and mechanisms responsible for CAD associations. Using these genome-wide maps, we prioritize 64 candidate variants and perform allele-specific binding and expression analyses at seven top candidate loci: 9p21.3, SMAD3, PDGFD, IL6R, BMP1, CCDC97/TGFB1 and LMOD1. We validate our findings in expression quantitative trait loci cohorts, which together reveal new links between CAD associations and regulatory function in the appropriate disease context