Mapping Variation in Cellular and Transcriptional Response to 1,25-Dihydroxyvitamin D3 in Peripheral Blood Mononuclear Cells

Trans -eQTL hotspots shape complex traits by modulating cellular states

Regulatory genetic variation shapes gene expression, providing an important mechanism connecting DNA variation and complex traits. The causal relationships between gene expression and complex traits remain poorly understood. Here, we integrated transcriptomes and 46 genetically complex growth traits in a large cross between two strains of the yeast Saccharomyces cerevisiae . We discovered thousands of genetic correlations between gene expression and growth, suggesting functional connections. Local regulatory variation was a minor source of these genetic correlations. Instead, genetic correlations tended to arise from multiple independent trans -acting regulatory loci. Trans -acting hotspots that affect the expression of numerous genes accounted for particularly large fractions of genetic growth variation and of genetic correlations between gene expression and growth. Genes with genetic correlations were enriched for similar biological processes across traits, but with heterogeneous direction of effect. Our results reveal how trans -acting regulatory hotspots shape complex traits by altering cellular states.

Vitamin D enhances type I IFN signaling in COVID-19 patients

The ability of Vitamin D (VitD) to modulate antiviral responses through induction of antimicrobial peptide is well established. However, the effect of VitD on host responses to SARS-CoV-2 is not well investigated. We here report the ability of VitD to enhance host IFN-alpha/beta (a/β) signaling both in vitro and among severe COVID-19 patients treated with VitD. Blood and saliva specimens were obtained from severe COVID-19 patients treated (43 patients), or not (37 patients), with vitD, during their stay in intensive care unit. Patients were followed up to 29 days following admission, and patient survival outcomes were collected. Higher activity levels of RIG-1/MDA-5 and JAK-STAT signaling pathways were observed with significantly higher gene and protein levels of antiviral interferon stimulating genes (ISGs) such as MX-1 and ISG-15; both in vitro, following treatment of PBMCs with vitD, and in whole blood and saliva specimens of VitD treated patients. Moreover, VitD treated patients had lower risk of all-cause mortality by day 29 compared to untreated patients (adjusted hazard ratio, 0.37, 95% confidence interval of 0.14-0.94; P = 0.038). The herein uncovered regulatory role of VitD on type I IFNs suggests the importance of insuring a normal level of VitD for the prevention and probably treatment of SARS-CoV-2 infection. Additional mechanistic studies, however, are needed to fully elucidate the antiviral effects of VitD particularly in the setting of COVID-19 infection.

Vitamin D boosts immune response of macrophages through a regulatory network of microRNAs and mRNAs

Vitamin D is associated with the stimulation of innate immunity, inflammation, and host defense against pathogens. Macrophages express receptors of Vitamin D, regulating transcription of genes related to immune processes. However, the transcriptional and post-transcriptional strategies controlling gene expression in differentiated macrophages, and how they are influenced by Vitamin D are not well understood. We studied whether Vitamin D enhances immune response by regulating the expression of microRNAs and mRNAs. Analysis of the transcriptome showed differences in expression of 199 genes, of which 68% were up-regulated, revealing the cell state of monocyte-derived macrophages differentiated with Vitamin D (D3-MDMs) as compared to monocyte-derived macrophages (MDMs). The differentially expressed genes appear to be associated with pathophysiological processes, including inflammatory responses, and cellular stress. Transcriptional motifs in promoter regions of up- or down-regulated genes showed enrichment of VDR motifs, suggesting possible roles of transcriptional activator or repressor in gene expression. Further, microRNA-Seq analysis indicated that there were 17 differentially expressed miRNAs, of which, 7 were up-regulated and 10 down-regulated, suggesting that Vitamin D plays a critical role in the regulation of miRNA expression during macrophages differentiation. The miR-6501-3p, miR-1273h-5p, miR-665, miR-1972, miR-1183, miR-619-5p were down-regulated in D3-MDMs compared to MDMs. The integrative analysis of miRNA and mRNA expression profiles predict that miR-1972, miR-1273h-5p, and miR-665 regulate genes PDCD1LG2, IL-1B, and CD274, which are related to the inflammatory response. Results suggest an essential role of Vitamin D in macrophage differentiation that modulates host response against pathogens, inflammation, and cellular stress.

Functional Characterization of Genetic Variant Effects on Expression

Thousands of common genetic variants in the human population have been associated with disease risk and phenotypic variation by genome-wide association studies (GWAS). However, the majority of GWAS variants fall into noncoding regions of the genome, complicating our understanding of their regulatory functions, and few molecular mechanisms of GWAS variant effects have been clearly elucidated. Here, we set out to review genetic variant effects, focusing on expression quantitative trait loci (eQTLs), including their utility in interpreting GWAS variant mechanisms. We discuss the interrelated challenges and opportunities for eQTL analysis, covering determining causal variants, elucidating molecular mechanisms of action, and understanding context variability. Addressing these questions can enable better functional characterization of disease-associated loci and provide insights into fundamental biological questions of the noncoding genetic regulatory code and its control of gene expression. Expected final online publication date for the Annual Review of Biomedical Data Science, Volume 5 is August 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Single-cell sequencing reveals lineage-specific dynamic genetic regulation of gene expression during human cardiomyocyte differentiation

Dynamic and temporally specific gene regulatory changes may underlie unexplained genetic associations with complex disease. During a dynamic process such as cellular differentiation, the overall cell type composition of a tissue (or an in vitro culture) and the gene regulatory profile of each cell can both experience significant changes over time. To identify these dynamic effects in high resolution, we collected single-cell RNA-sequencing data over a differentiation time course from induced pluripotent stem cells to cardiomyocytes, sampled at 7 unique time points in 19 human cell lines. We employed a flexible approach to map dynamic eQTLs whose effects vary significantly over the course of bifurcating differentiation trajectories, including many whose effects are specific to one of these two lineages. Our study design allowed us to distinguish true dynamic eQTLs affecting a specific cell lineage from expression changes driven by potentially non-genetic differences between cell lines such as cell composition. Additionally, we used the cell type profiles learned from single-cell data to deconvolve and re-analyze data from matched bulk RNA-seq samples. Using this approach, we were able to identify a large number of novel dynamic eQTLs in single cell data while also attributing dynamic effects in bulk to a particular lineage. Overall, we found that using single cell data to uncover dynamic eQTLs can provide new insight into the gene regulatory changes that occur among heterogeneous cell types during cardiomyocyte differentiation.

Whole transcriptome analysis of high and low IFN-α producers reveals differential response patterns following rhinovirus stimulation

Objectives

Viral respiratory infections cause considerable morbidity and economic loss. While rhinoviruses (RV) typically cause little more than the common cold, they can produce severe infections and disease exacerbations in susceptible individuals, such as those with asthma. Variations in the regulation of key antiviral cytokines, particularly type I interferon (IFN‐α and IFN‐β), may contribute to RV susceptibility. To understand this variability, we compared the transcriptomes of high and low type I IFN producers.

Methods

Blood mononuclear cells from 238 individuals with or without asthma were cultured in the presence or absence of RV. Those samples demonstrating high or low RV‐stimulated IFN‐α production (N = 75) underwent RNA‐sequencing.

Results

Gene expression patterns were similar in samples from healthy participants and those with asthma. At baseline, the high IFN‐α producer group showed higher expression of genes associated with plasmacytoid dendritic cells, the innate immune response and vitamin D activation, but lower expression of oxidative stress pathways than the low IFN‐α producer group. After RV stimulation, the high IFN‐α producer group showed higher expression of genes found in immune response biological pathways and lower expression of genes linked to developmental and catabolic processes when compared to the low IFN‐α producer group.

Conclusions

These differences suggest that the high IFN‐α group has a higher level of immune system readiness, resulting in a more intense and perhaps more focussed pathogen‐specific immune response. These results contribute to a better understanding of the variability in type I IFN production between individuals.

Dynamic effects of genetic variation on gene expression revealed following hypoxic stress in cardiomyocytes

One life-threatening outcome of cardiovascular disease is myocardial infarction, where cardiomyocytes are deprived of oxygen. To study inter-individual differences in response to hypoxia, we established an in vitro model of induced pluripotent stem cell-derived cardiomyocytes from 15 individuals. We measured gene expression levels, chromatin accessibility, and methylation levels in four culturing conditions that correspond to normoxia, hypoxia, and short- or long-term re-oxygenation. We characterized thousands of gene regulatory changes as the cells transition between conditions. Using available genotypes, we identified 1,573 genes with a cis expression quantitative locus (eQTL) in at least one condition, as well as 367 dynamic eQTLs, which are classified as eQTLs in at least one, but not in all conditions. A subset of genes with dynamic eQTLs is associated with complex traits and disease. Our data demonstrate how dynamic genetic effects on gene expression, which are likely relevant for disease, can be uncovered under stress.

Common and personal target genes of the micronutrient vitamin D in primary immune cells from human peripheral blood

Vitamin D is essential for the function of the immune system. In this study, we treated peripheral blood mononuclear cells (PBMCs) of healthy adults with the biologically active form of vitamin D₃, 1α,25-dihydroxyvitamin D₃ (1,25(OH)₂D₃) using two different approaches: single repeats with PBMCs obtained from a cohort of 12 individuals and personalized analysis based on triplicates of five study participants. This identified 877 (cohort approach) and 3951 (personalized approach) genes that significantly (p < 0.05) changed their expression 24 h after 1,25(OH)₂D₃ stimulation. From these, 333 and 1232 were classified as supertargets, a third of which were identified as novel. Individuals differed largely in their vitamin D response not only by the magnitude of expression change but also by their personal selection of (super)target genes. Functional analysis of the target genes suggested the overarching role of vitamin D in the regulation of metabolism, proliferation and differentiation, but in particular in the control of functions mediated by the innate and adaptive immune system, such as responses to infectious diseases and chronic inflammatory disorders. In conclusion, immune cells are an important target of vitamin D and common genes may serve as biomarkers for personal responses to the micronutrient.

In vivo transcriptome changes of human white blood cells in response to vitamin D

In the vitamin D intervention study VitDbol (NCT02063334) blood samples were drawn directly before an oral bolus (2000 μg vitamin D₃) and 24 h later. The focus of phase II of VitDbol was the transcriptome-wide analysis of the effects of vitamin D gene expression in human peripheral blood mononuclear cells (PBMCs). All five participants responded in an individual fashion to the bolus by increases in serum levels of the vitamin D metabolites 25-hydroxyvitamin D₃ (25(OH)D₃) and 1α,25-dihydroxyvitamin D₃ (1,25(OH)₂D₃). RNA sequencing identified 15.040 commonly expressed genes in PBMCs, 702 (4,7%) of which were significantly (p < 0,05) affected by the vitamin D₃ bolus. KEGG pathway analysis suggested that these genes are involved in general protein translation, monocyte differentiation and cellular growth control. Previously published transcriptome-wide studies in comparable cell systems confirmed 234 of the 702 vitamin D target genes, leaving many genes, such as HLA-A and HLA-C, as novel discoveries. Interestingly, in vivo stimulated PBMCs of this study showed a larger number of common vitamin D target genes with the monocytic cell line THP-1 than with in vitro stimulated PBMCs. The expression pattern of vitamin D target genes differed significantly between individuals and the average expression change can serve as a marker for vitamin D responsiveness. In conclusion, this study demonstrates that under in vivo conditions changes in 25(OH)D₃ and 1,25(OH)₂D₃ serum concentrations alter the expression of more than 700 vitamin D target genes in human leukocytes.

Association between preoperative peripheral blood mononuclear cell gene expression profiles, early postoperative organ function recovery potential and long-term survival in advanced heart failure patients undergoing mechanical circulatory support

BACKGROUND

Multiorgan dysfunction syndrome contributes to adverse outcomes in advanced heart failure (AdHF) patients after mechanical circulatory support (MCS) implantation and is associated with aberrant leukocyte activity. We tested the hypothesis that preoperative peripheral blood mononuclear cell (PBMC) gene expression profiles (GEP) can predict early postoperative improvement or non-improvement in patients undergoing MCS implantation. We believe this information may be useful in developing prognostic biomarkers.

METHODS & DESIGN

We conducted a study with 29 patients undergoing MCS-surgery in a tertiary academic medical center from 2012 to 2014. PBMC samples were collected one day before surgery (day -1). Clinical data was collected on day -1 and day 8 postoperatively. Patients were classified by Sequential Organ Failure Assessment score and Model of End-stage Liver Disease Except INR score (measured eight days after surgery): Group I = improving (both scores improved from day -1 to day 8, n = 17) and Group II = not improving (either one or both scores did not improve from day -1 to day 8, n = 12). RNA-sequencing was performed on purified mRNA and analyzed using Next Generation Sequencing Strand. Differentially expressed genes (DEGs) were identified by Mann-Whitney test with Benjamini-Hochberg correction. Preoperative DEGs were used to construct a support vector machine algorithm to predict Group I vs. Group II membership.

RESULTS

Out of 28 MCS-surgery patients alive 8 days postoperatively, one-year survival was 88% in Group I and 27% in Group II. We identified 28 preoperative DEGs between Group I and II, with an average 93% prediction accuracy. Out of 105 DEGs identified preoperatively between year 1 survivors and non-survivors, 12 genes overlapped with the 28 predictive genes.

CONCLUSIONS

In AdHF patients following MCS implantation, preoperative PBMC-GEP predicts early changes in organ function scores and correlates with long-term outcomes. Therefore, gene expression lends itself to outcome prediction and warrants further studies in larger longitudinal cohorts.

Colonic transcriptional response to 1α,25(OH)2 vitamin D3 in African- and European-Americans

Colorectal cancer (CRC) is a significant health burden especially among African Americans (AA). Epidemiological studies have correlated low serum vitamin D with CRC risk, and, while hypovitaminosis D is more common and more severe in AA, the mechanisms by which vitamin D modulates CRC risk and how these differ by race are not well understood. Active vitamin D (1α,25(OH)₂D₃) has chemoprotective effects primarily through transcriptional regulation of target genes in the colon. We hypothesized that transcriptional response to 1α,25(OH)₂D₃ differs between AA and European Americans (EA) irrespective of serum vitamin D and that regulatory variants could impact transcriptional response. We treated ex vivo colon cultures from 34 healthy subjects (16 AA and 18 EA) with 0.1μM 1α,25(OH)₂D₃ or vehicle control for 6h and performed genome-wide transcriptional profiling. We found 8 genes with significant differences in transcriptional response to 1α,25(OH)₂D₃ between AA and EA with definitive replication of inter-ethnic differences for uridine phosphorylase 1 (UPP1) and zinc finger-SWIM containing 4 (ZSWIM4). We performed expression quantitative trait loci (eQTL) mapping and identified response cis-eQTLs for ZSWIM4 as well as for histone deacetylase 3 (HDAC3), the latter of which showed a trend toward significant inter-ethnic differences in transcriptional response. Allele frequency differences of eQTLs for ZSWIM4 and HDAC3 accounted for observed transcriptional differences between populations. Taken together, our results demonstrate that transcriptional response to 1α,25(OH)₂D₃ differs between AA and EA independent of serum 25(OH)D levels. We provide evidence in support of a genetic regulatory mechanism underlying transcriptional differences between populations for ZSWIM4 and HDAC3. Further work is needed to elucidate how response eQTLs modify vitamin D response and whether genotype and/or transcriptional response correlate with chemopreventive effects. Relevant biomarkers, such as tissue-specific 1α,25(OH)₂D₃ transcriptional response, could identify individuals likely to benefit from vitamin D for CRC prevention as well as elucidate basic mechanisms underlying CRC disparities.