Deconvolution of bulk blood eQTL effects into immune cell subpopulations

Antibody signatures against viruses and microbiome reflect past and chronic exposures and associate with aging and inflammation

Encounters with pathogens and other molecules can imprint long-lasting effects on our immune system, influencing future physiological outcomes. Given the wide range of microbes to which humans are exposed, their collective impact on health is not fully understood. To explore relations between exposures and biological aging and inflammation, we profiled an antibody-binding repertoire against 2,815 microbial, viral, and environmental peptides in a population cohort of 1,443 participants. Utilizing antibody-binding as a proxy for past exposures, we investigated their impact on biological aging, cell composition, and inflammation. Immune response against cytomegalovirus (CMV), rhinovirus, and gut bacteria relates with telomere length. Single-cell expression measurements identified an effect of CMV infection on the transcriptional landscape of subpopulations of CD8 and CD4 T-cells. This examination of the relationship between microbial exposures and biological aging and inflammation highlights a role for chronic infections (CMV and Epstein-Barr virus) and common pathogens (rhinoviruses and adenovirus C).

Quantitative and Longitudinal Assessment of Systemic Innate Immunity in Health and Disease Using a 2D Gene Model

Dysregulation of innate immunity is deeply involved in infectious and autoimmune diseases. For a better understanding of pathogenesis and improved management of these diseases, it is of vital importance to implement convenient monitoring of systemic innate immunity. Built upon our previous works on the host transcriptional response to infection in peripheral blood, we proposed a 2D gene model for the simultaneous assessment of two major components of systemic innate immunity, including VirSig as the signature of the host response to viral infection and BacSig as the signature of the host response to bacterial infection. The revelation of dysregulation in innate immunity by this 2D gene model was demonstrated with a wide variety of transcriptome datasets. In acute infection, distinctive patterns of VirSig and BacSig activation were observed in viral and bacterial infection. In comparison, both signatures were restricted to a defined range in the vast majority of healthy adults, regardless of age. In addition, BacSig showed significant elevation during pregnancy and an upward trend during development. In tuberculosis (TB), elevation of BacSig and VirSig was observed in a significant portion of active TB patients, and abnormal BacSig was also associated with a longer treatment course. In cystic fibrosis (CF), abnormal BacSig was observed in a subset of patients, and no overall change in BacSig abnormality was observed after the drug treatment. In systemic sclerosis-associated interstitial lung disease (SSc-ILD), significant elevation of VirSig and BacSig was observed in some patients, and treatment with a drug led to the further deviation of BacSig from the control level. In systemic lupus erythematosus (SLE), positivity for the anti-Ro autoantibody was associated with significant elevation of VirSig in SLE patients, and the additive effect of VirSig/BacSig activation was also observed in SLE patients during pregnancy. Overall, these data demonstrated that the 2D gene model can be used to assess systemic innate immunity in health and disease, with the potential clinical applications including patient stratification, prescription of antibiotics, understanding of pathogenesis, and longitudinal monitoring of treatment response.

Adenylate cyclase activates the cAMP signalling pathway to enhance platelet-rich plasma-treated Achilles tendon disease, a theoretical bioinformatics-based study

The effectiveness of platelet-rich plasma (PRP) for the treatment of Achilles tendon disorders still needs to be evaluated through a series of prospective studies, but genomic analysis can reveal the existence of complementary PRP treatment options. Based on the 96 platelet activation-related genes in the Kyoto Encyclopedia of Genes and Genomes (KEGG) database, we performed Gene Ontology functional enrichment analysis and KEGG enrichment analysis, pathway correlation analysis, and enrichment mapping to determine the enrichment results of the gene set enrichment analysis and found that the cAMP signalling pathway may be the key to enhancing the effectiveness of PRP treatment. The cAMP signalling pathway interacts with the Rap1 signalling pathway and cGMP-PKG signalling pathway to mediate the entire pathophysiological process of Achilles tendon disease. Moreover, ADCY1-9 may be the key to the activation of the cAMP signalling network. Further based on the data in the Gene Expression Omnibus database, it was found that ADCY4 and ADCY7 may be the players that play a major role, associated with the STAT4-ADCY4-LAMA5 axis and the GRbeta-ADCY7-SEMA3C axis, which is expected to be a complementary target for enhancing the efficacy of PRP in the treatment of Achilles tendon disease.

Cell-type deconvolution of bulk-blood RNA-seq reveals biological insights into neuropsychiatric disorders

Genome-wide association studies (GWASs) have uncovered susceptibility loci associated with psychiatric disorders such as bipolar disorder (BP) and schizophrenia (SCZ). However, most of these loci are in non-coding regions of the genome, and the causal mechanisms of the link between genetic variation and disease risk is unknown. Expression quantitative trait locus (eQTL) analysis of bulk tissue is a common approach used for deciphering underlying mechanisms, although this can obscure cell-type-specific signals and thus mask trait-relevant mechanisms. Although single-cell sequencing can be prohibitively expensive in large cohorts, computationally inferred cell-type proportions and cell-type gene expression estimates have the potential to overcome these problems and advance mechanistic studies. Using bulk RNA-seq from 1,730 samples derived from whole blood in a cohort ascertained from individuals with BP and SCZ, this study estimated cell-type proportions and their relation with disease status and medication. For each cell type, we found between 2,875 and 4,629 eGenes (genes with an associated eQTL), including 1,211 that are not found on the basis of bulk expression alone. We performed a colocalization test between cell-type eQTLs and various traits and identified hundreds of associations that occur between cell-type eQTLs and GWASs but that are not detected in bulk eQTLs. Finally, we investigated the effects of lithium use on the regulation of cell-type expression loci and found examples of genes that are differentially regulated according to lithium use. Our study suggests that applying computational methods to large bulk RNA-seq datasets of non-brain tissue can identify disease-relevant, cell-type-specific biology of psychiatric disorders and psychiatric medication.

Interaction molecular QTL mapping discovers cellular and environmental modifiers of genetic regulatory effects

Bulk-tissue molecular quantitative trait loci (QTLs) have been the starting point for interpreting disease-associated variants, and context-specific QTLs show particular relevance for disease. Here, we present the results of mapping interaction QTLs (iQTLs) for cell type, age, and other phenotypic variables in multi-omic, longitudinal data from the blood of individuals of diverse ancestries. By modeling the interaction between genotype and estimated cell-type proportions, we demonstrate that cell-type iQTLs could be considered as proxies for cell-type-specific QTL effects, particularly for the most abundant cell type in the tissue. The interpretation of age iQTLs, however, warrants caution because the moderation effect of age on the genotype and molecular phenotype association could be mediated by changes in cell-type composition. Finally, we show that cell-type iQTLs contribute to cell-type-specific enrichment of diseases that, in combination with additional functional data, could guide future functional studies. Overall, this study highlights the use of iQTLs to gain insights into the context specificity of regulatory effects.

eQTL studies: from bulk tissues to single cells

An expression quantitative trait locus (eQTL) is a chromosomal region where genetic variants are associated with the expression levels of specific genes that can be both nearby or distant. The identifications of eQTLs for different tissues, cell types, and contexts have led to a better understanding of the dynamic regulations of gene expressions and implications of functional genes and variants for complex traits and diseases. Although most eQTL studies have been performed on data collected from bulk tissues, recent studies have demonstrated the importance of cell-type-specific and context-dependent gene regulations in biological processes and disease mechanisms. In this review, we discuss statistical methods that have been developed to enable the detection of cell-type-specific and context-dependent eQTLs from bulk tissues, purified cell types, and single cells. We also discuss the limitations of the current methods and future research opportunities.

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.

The role of DNA methylation in personalized medicine for immune-related diseases

Epigenetics functions as a bridge between host genetic & environmental factors, aiding in human health and diseases. Many immune-related diseases, including infectious and allergic diseases, have been linked to epigenetic mechanisms, particularly DNA methylation. In this review, we summarized an updated overview of DNA methylation and its importance in personalized medicine, and demonstrated that DNA methylation has excellent potential for disease prevention, diagnosis, and treatment in a personalized manner. The future implications and limitations of the DNA methylation study have also been well-discussed.

Interpreting non-coding disease-associated human variants using single-cell epigenomics

Genome-wide association studies (GWAS) have linked hundreds of thousands of sequence variants in the human genome to common traits and diseases. However, translating this knowledge into a mechanistic understanding of disease-relevant biology remains challenging, largely because such variants are predominantly in non-protein-coding sequences that still lack functional annotation at cell-type resolution. Recent advances in single-cell epigenomics assays have enabled the generation of cell type-, subtype- and state-resolved maps of the epigenome in heterogeneous human tissues. These maps have facilitated cell type-specific annotation of candidate cis-regulatory elements and their gene targets in the human genome, enhancing our ability to interpret the genetic basis of common traits and diseases.

Interaction molecular QTL mapping discovers cellular and environmental modifiers of genetic regulatory effects

Bulk tissue molecular quantitative trait loci (QTLs) have been the starting point for interpreting disease-associated variants, while context-specific QTLs show particular relevance for disease. Here, we present the results of mapping interaction QTLs (iQTLs) for cell type, age, and other phenotypic variables in multi-omic, longitudinal data from blood of individuals of diverse ancestries. By modeling the interaction between genotype and estimated cell type proportions, we demonstrate that cell type iQTLs could be considered as proxies for cell type-specific QTL effects. The interpretation of age iQTLs, however, warrants caution as the moderation effect of age on the genotype and molecular phenotype association may be mediated by changes in cell type composition. Finally, we show that cell type iQTLs contribute to cell type-specific enrichment of diseases that, in combination with additional functional data, may guide future functional studies. Overall, this study highlights iQTLs to gain insights into the context-specificity of regulatory effects.

Sex-biased gene expression at single-cell resolution: cause and consequence of sexual dimorphism

Gene expression differences between males and females are thought to be key for the evolution of sexual dimorphism, and sex-biased genes are often used to study the molecular footprint of sex-specific selection. However, gene expression is often measured from complex aggregations of diverse cell types, making it difficult to distinguish between sex differences in expression that are due to regulatory rewiring within similar cell types and those that are simply a consequence of developmental differences in cell-type abundance. To determine the role of regulatory versus developmental differences underlying sex-biased gene expression, we use single-cell transcriptomic data from multiple somatic and reproductive tissues of male and female guppies, a species that exhibits extensive phenotypic sexual dimorphism. Our analysis of gene expression at single-cell resolution demonstrates that nonisometric scaling between the cell populations within each tissue and heterogeneity in cell-type abundance between the sexes can influence inferred patterns of sex-biased gene expression by increasing both the false-positive and false-negative rates. Moreover, we show that, at the bulk level, the subset of sex-biased genes that are the product of sex differences in cell-type abundance can significantly confound patterns of coding-sequence evolution. Taken together, our results offer a unique insight into the effects of allometry and cellular heterogeneity on perceived patterns of sex-biased gene expression and highlight the power of single-cell RNA-sequencing in distinguishing between sex-biased genes that are the result of regulatory change and those that stem from sex differences in cell-type abundance, and hence are a consequence rather than a cause of sexual dimorphism.

A computational method for cell type-specific expression quantitative trait loci mapping using bulk RNA-seq data

Mapping cell type-specific gene expression quantitative trait loci (ct-eQTLs) is a powerful way to investigate the genetic basis of complex traits. A popular method for ct-eQTL mapping is to assess the interaction between the genotype of a genetic locus and the abundance of a specific cell type using a linear model. However, this approach requires transforming RNA-seq count data, which distorts the relation between gene expression and cell type proportions and results in reduced power and/or inflated type I error. To address this issue, we have developed a statistical method called CSeQTL that allows for ct-eQTL mapping using bulk RNA-seq count data while taking advantage of allele-specific expression. We validated the results of CSeQTL through simulations and real data analysis, comparing CSeQTL results to those obtained from purified bulk RNA-seq data or single cell RNA-seq data. Using our ct-eQTL findings, we were able to identify cell types relevant to 21 categories of human traits.

Identification of genetic variants that impact gene co-expression relationships using large-scale single-cell data

BACKGROUND

Expression quantitative trait loci (eQTL) studies show how genetic variants affect downstream gene expression. Single-cell data allows reconstruction of personalized co-expression networks and therefore the identification of SNPs altering co-expression patterns (co-expression QTLs, co-eQTLs) and the affected upstream regulatory processes using a limited number of individuals.

RESULTS

We conduct a co-eQTL meta-analysis across four scRNA-seq peripheral blood mononuclear cell datasets using a novel filtering strategy followed by a permutation-based multiple testing approach. Before the analysis, we evaluate the co-expression patterns required for co-eQTL identification using different external resources. We identify a robust set of cell-type-specific co-eQTLs for 72 independent SNPs affecting 946 gene pairs. These co-eQTLs are replicated in a large bulk cohort and provide novel insights into how disease-associated variants alter regulatory networks. One co-eQTL SNP, rs1131017, that is associated with several autoimmune diseases, affects the co-expression of RPS26 with other ribosomal genes. Interestingly, specifically in T cells, the SNP additionally affects co-expression of RPS26 and a group of genes associated with T cell activation and autoimmune disease. Among these genes, we identify enrichment for targets of five T-cell-activation-related transcription factors whose binding sites harbor rs1131017. This reveals a previously overlooked process and pinpoints potential regulators that could explain the association of rs1131017 with autoimmune diseases.

CONCLUSION

Our co-eQTL results highlight the importance of studying context-specific gene regulation to understand the biological implications of genetic variation. With the expected growth of sc-eQTL datasets, our strategy and technical guidelines will facilitate future co-eQTL identification, further elucidating unknown disease mechanisms.

Brain expression quantitative trait locus and network analyses reveal downstream effects and putative drivers for brain-related diseases

Identification of therapeutic targets from genome-wide association studies (GWAS) requires insights into downstream functional consequences. We harmonized 8,613 RNA-sequencing samples from 14 brain datasets to create the MetaBrain resource and performed cis- and trans-expression quantitative trait locus (eQTL) meta-analyses in multiple brain region- and ancestry-specific datasets (n ≤ 2,759). Many of the 16,169 cortex cis-eQTLs were tissue-dependent when compared with blood cis-eQTLs. We inferred brain cell types for 3,549 cis-eQTLs by interaction analysis. We prioritized 186 cis-eQTLs for 31 brain-related traits using Mendelian randomization and co-localization including 40 cis-eQTLs with an inferred cell type, such as a neuron-specific cis-eQTL (CYP24A1) for multiple sclerosis. We further describe 737 trans-eQTLs for 526 unique variants and 108 unique genes. We used brain-specific gene-co-regulation networks to link GWAS loci and prioritize additional genes for five central nervous system diseases. This study represents a valuable resource for post-GWAS research on central nervous system diseases.

eQTL Studies: from Bulk Tissues to Single Cells

An expression quantitative trait locus (eQTL) is a chromosomal region where genetic variants are associated with the expression levels of certain genes that can be both nearby or distant. The identifications of eQTLs for different tissues, cell types, and contexts have led to better understanding of the dynamic regulations of gene expressions and implications of functional genes and variants for complex traits and diseases. Although most eQTL studies to date have been performed on data collected from bulk tissues, recent studies have demonstrated the importance of cell-type-specific and context-dependent gene regulations in biological processes and disease mechanisms. In this review, we discuss statistical methods that have been developed to enable the detections of cell-type-specific and context-dependent eQTLs from bulk tissues, purified cell types, and single cells. We also discuss the limitations of the current methods and future research opportunities.

An immune dysfunction score for stratification of patients with acute infection based on whole-blood gene expression

Dysregulated host responses to infection can lead to organ dysfunction and sepsis, causing millions of global deaths each year. To alleviate this burden, improved prognostication and biomarkers of response are urgently needed. We investigated the use of whole-blood transcriptomics for stratification of patients with severe infection by integrating data from 3149 samples from patients with sepsis due to community-acquired pneumonia or fecal peritonitis admitted to intensive care and healthy individuals into a gene expression reference map. We used this map to derive a quantitative sepsis response signature (SRSq) score reflective of immune dysfunction and predictive of clinical outcomes, which can be estimated using a 7- or 12-gene signature. Last, we built a machine learning framework, SepstratifieR, to deploy SRSq in adult and pediatric bacterial and viral sepsis, H1N1 influenza, and COVID-19, demonstrating clinically relevant stratification across diseases and revealing some of the physiological alterations linking immune dysregulation to mortality. Our method enables early identification of individuals with dysfunctional immune profiles, bringing us closer to precision medicine in infection.

Limited evidence for blood eQTLs in human sexual dimorphism

BACKGROUND

The genetic underpinning of sexual dimorphism is very poorly understood. The prevalence of many diseases differs between men and women, which could be in part caused by sex-specific genetic effects. Nevertheless, only a few published genome-wide association studies (GWAS) were performed separately in each sex. The reported enrichment of expression quantitative trait loci (eQTLs) among GWAS-associated SNPs suggests a potential role of sex-specific eQTLs in the sex-specific genetic mechanism underlying complex traits.

METHODS

To explore this scenario, we combined sex-specific whole blood RNA-seq eQTL data from 3447 European individuals included in BIOS Consortium and GWAS data from UK Biobank. Next, to test the presence of sex-biased causal effect of gene expression on complex traits, we performed sex-specific transcriptome-wide Mendelian randomization (TWMR) analyses on the two most sexually dimorphic traits, waist-to-hip ratio (WHR) and testosterone levels. Finally, we performed power analysis to calculate the GWAS sample size needed to observe sex-specific trait associations driven by sex-biased eQTLs.

RESULTS

Among 9 million SNP-gene pairs showing sex-combined associations, we found 18 genes with significant sex-biased cis-eQTLs (FDR 5%). Our phenome-wide association study of the 18 top sex-biased eQTLs on >700 traits unraveled that these eQTLs do not systematically translate into detectable sex-biased trait-associations. In addition, we observed that sex-specific causal effects of gene expression on complex traits are not driven by sex-specific eQTLs. Power analyses using real eQTL- and causal-effect sizes showed that millions of samples would be necessary to observe sex-biased trait associations that are fully driven by sex-biased cis-eQTLs. Compensatory effects may further hamper their detection.

CONCLUSIONS

Our results suggest that sex-specific eQTLs in whole blood do not translate to detectable sex-specific trait associations of complex diseases, and vice versa that the observed sex-specific trait associations cannot be explained by sex-specific eQTLs.

Detecting signatures of selection on gene expression

A substantial amount of phenotypic diversity results from changes in gene expression levels and patterns. Understanding how the transcriptome evolves is therefore a key priority in identifying mechanisms of adaptive change. However, in contrast to powerful models of sequence evolution, we lack a consensus model of gene expression evolution. Furthermore, recent work has shown that many of the comparative approaches used to study gene expression are subject to biases that can lead to false signatures of selection. Here we first outline the main approaches for describing expression evolution and their inherent biases. Next, we bridge the gap between the fields of phylogenetic comparative methods and transcriptomics to reinforce the main pitfalls of inferring selection on expression patterns and use simulation studies to show that shifts in tissue composition can heavily bias inferences of selection. We close by highlighting the multi-dimensional nature of transcriptional variation and identifying major unanswered questions in disentangling how selection acts on the transcriptome.

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.

Immune disease variants modulate gene expression in regulatory CD4+ T cells

Identifying cellular functions dysregulated by disease-associated variants could implicate novel pathways for drug targeting or modulation in cell therapies. However, follow-up studies can be challenging if disease-relevant cell types are difficult to sample. Variants associated with immune diseases point toward the role of CD4+ regulatory T cells (Treg cells). We mapped genetic regulation (quantitative trait loci [QTL]) of gene expression and chromatin activity in Treg cells, and we identified 133 colocalizing loci with immune disease variants. Colocalizations of immune disease genome-wide association study (GWAS) variants with expression QTLs (eQTLs) controlling the expression of CD28 and STAT5A, involved in Treg cell activation and interleukin-2 (IL-2) signaling, support the contribution of Treg cells to the pathobiology of immune diseases. Finally, we identified seven known drug targets suitable for drug repurposing and suggested 63 targets with drug tractability evidence among the GWAS signals that colocalized with Treg cell QTLs. Our study is the first in-depth characterization of immune disease variant effects on Treg cell gene expression modulation and dysregulation of Treg cell function.

Dissecting the transcriptome in cardiovascular disease

The human transcriptome comprises a complex network of coding and non-coding RNAs implicated in a myriad of biological functions. Non-coding RNAs exhibit highly organized spatial and temporal expression patterns and are emerging as critical regulators of differentiation, homeostasis, and pathological states, including in the cardiovascular system. This review defines the current knowledge gaps, unmet methodological needs, and describes the challenges in dissecting and understanding the role and regulation of the non-coding transcriptome in cardiovascular disease. These challenges include poor annotation of the non-coding genome, determination of the cellular distribution of transcripts, assessment of the role of RNA processing and identification of cell-type specific changes in cardiovascular physiology and disease. We highlight similarities and differences in the hurdles associated with the analysis of the non-coding and protein-coding transcriptomes. In addition, we discuss how the lack of consensus and absence of standardized methods affect reproducibility of data. These shortcomings should be defeated in order to make significant scientific progress and foster the development of clinically applicable non-coding RNA-based therapeutic strategies to lessen the burden of cardiovascular disease.

Integration of functional genomics data to uncover cell type-specific pathways affected in Parkinson's disease

Parkinson's disease (PD) is the second most prevalent late-onset neurodegenerative disorder worldwide after Alzheimer's disease for which available drugs only deliver temporary symptomatic relief. Loss of dopaminergic neurons (DaNs) in the substantia nigra and intracellular alpha-synuclein inclusions are the main hallmarks of the disease but the events that cause this degeneration remain uncertain. Despite cell types other than DaNs such as astrocytes, microglia and oligodendrocytes have been recently associated with the pathogenesis of PD, we still lack an in-depth characterisation of PD-affected brain regions at cell-type resolution that could help our understanding of the disease mechanisms. Nevertheless, publicly available large-scale brain-specific genomic, transcriptomic and epigenomic datasets can be further exploited to extract different layers of cell type-specific biological information for the reconstruction of cell type-specific transcriptional regulatory networks. By intersecting disease risk variants within the networks, it may be possible to study the functional role of these risk variants and their combined effects at cell type- and pathway levels, that, in turn, can facilitate the identification of key regulators involved in disease progression, which are often potential therapeutic targets.

Identification of tissue-specific and common methylation quantitative trait loci in healthy individuals using MAGAR

Background

Understanding the influence of genetic variants on DNA methylation is fundamental for the interpretation of epigenomic data in the context of disease. There is a need for systematic approaches not only for determining methylation quantitative trait loci (methQTL), but also for discriminating general from cell type-specific effects.

Results

Here, we present a two-step computational framework MAGAR (https://bioconductor.org/packages/MAGAR), which fully supports the identification of methQTLs from matched genotyping and DNA methylation data, and additionally allows for illuminating cell type-specific methQTL effects. In a pilot analysis, we apply MAGAR on data in four tissues (ileum, rectum, T cells, B cells) from healthy individuals and demonstrate the discrimination of common from cell type-specific methQTLs. We experimentally validate both types of methQTLs in an independent data set comprising additional cell types and tissues. Finally, we validate selected methQTLs located in the PON1, ZNF155, and NRG2 genes by ultra-deep local sequencing. In line with previous reports, we find cell type-specific methQTLs to be preferentially located in enhancer elements.

Conclusions

Our analysis demonstrates that a systematic analysis of methQTLs provides important new insights on the influences of genetic variants to cell type-specific epigenomic variation.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13072-021-00415-6.

Multi-Omics Approaches in Immunological Research

The immune system plays a vital role in health and disease, and is regulated through a complex interactive network of many different immune cells and mediators. To understand the complexity of the immune system, we propose to apply a multi-omics approach in immunological research. This review provides a complete overview of available methodological approaches for the different omics data layers relevant for immunological research, including genetics, epigenetics, transcriptomics, proteomics, metabolomics, and cellomics. Thereafter, we describe the various methods for data analysis as well as how to integrate different layers of omics data. Finally, we discuss the possible applications of multi-omics studies and opportunities they provide for understanding the complex regulatory networks as well as immune variation in various immune-related diseases.

Sex-specific gene expression in the blood of four primates

Blood is an important non-reproductive tissue, but little is known about the sex-specific gene expressions in the blood. Therefore, we investigated sex-specific gene expression differences in the blood tissues of four primates, rhesus macaques (Macaca mulatta), Tibetan macaques (M. thibetana), yellow baboons (Papio cynocephalus), and humans. We identified seven sex-specific differentially expressed genes (SDEGs) in each non-human primate and 31 SDEGs in humans. The four primates had only one common SDEG, MAP7D2. In humans, immune-related SDEGs were identified as up-regulated, but also down-regulated in females. We also found that most of the X-Y gene pairs had similar expression levels between species, except pair EIF1AY/EIF1AX. The expression level of X-Y gene pairs of rhesus and Tibetan macaques showed no significant differential expression levels, while humans had six significant XY-biased and three XX-biased X-Y gene pairs. Our observed sex differences in blood should increase understanding of sex differences in primate blood tissue.

Targeting of the CD80/86 proinflammatory axis as a therapeutic strategy to prevent severe COVID-19

An excessive immune response known as cytokine storm is the hallmark of severe COVID-19. The cause of this cytokine rampage is yet not known. Based on recent epidemiological evidence, we hypothesized that CD80/86 signaling is essential for this hyperinflammation, and that blocking this proinflammatory axis could be an effective therapeutic approach to protect against severe COVID-19. Here we provide exploratory evidence that abatacept, a drug that blocks CD80/86 co-stimulation, produces changes at the systemic level that are highly antagonistic of the proinflammatory processes elicited by COVID-19. Using RNA-seq from blood samples from a longitudinal cohort of n = 38 rheumatic patients treated with abatacept, we determined the immunological processes that are significantly regulated by this treatment. We then analyzed available blood RNA-seq from two COVID19 patient cohorts, a very early cohort from the epicenter of the pandemic in China (n = 3 COVID-19 cases and n = 3 controls), and a recent and larger cohort from the USA (n = 49 severe and n = 51 mild COVD-19 patients). We found a highly significant antagonism between SARS-CoV-2 infection and COVID-19 severity with the systemic response to abatacept. Analysis of previous single-cell RNA-seq data from bronchoalveolar lavage fluid from mild and severe COVID-19 patients and controls, reinforce the implication of the CD80/86 proinflammatory axis. Our functional results further support abatacept as a candidate therapeutic approach to prevent severe COVID-19.

Feasibility of predicting allele specific expression from DNA sequencing using machine learning

Allele specific expression (ASE) concerns divergent expression quantity of alternative alleles and is measured by RNA sequencing. Multiple studies show that ASE plays a role in hereditary diseases by modulating penetrance or phenotype severity. However, genome diagnostics is based on DNA sequencing and therefore neglects gene expression regulation such as ASE. To take advantage of ASE in absence of RNA sequencing, it must be predicted using only DNA variation. We have constructed ASE models from BIOS (n = 3432) and GTEx (n = 369) that predict ASE using DNA features. These models are highly reproducible and comprise many different feature types, highlighting the complex regulation that underlies ASE. We applied the BIOS-trained model to population variants in three genes in which ASE plays a clinically relevant role: BRCA2, RET and NF1. This resulted in predicted ASE effects for 27 variants, of which 10 were known pathogenic variants. We demonstrated that ASE can be predicted from DNA features using machine learning. Future efforts may improve sensitivity and translate these models into a new type of genome diagnostic tool that prioritizes candidate pathogenic variants or regulators thereof for follow-up validation by RNA sequencing. All used code and machine learning models are available at GitHub and Zenodo.

DeCompress: tissue compartment deconvolution of targeted mRNA expression panels using compressed sensing

Targeted mRNA expression panels, measuring up to 800 genes, are used in academic and clinical settings due to low cost and high sensitivity for archived samples. Most samples assayed on targeted panels originate from bulk tissue comprised of many cell types, and cell-type heterogeneity confounds biological signals. Reference-free methods are used when cell-type-specific expression references are unavailable, but limited feature spaces render implementation challenging in targeted panels. Here, we present DeCompress, a semi-reference-free deconvolution method for targeted panels. DeCompress leverages a reference RNA-seq or microarray dataset from similar tissue to expand the feature space of targeted panels using compressed sensing. Ensemble reference-free deconvolution is performed on this artificially expanded dataset to estimate cell-type proportions and gene signatures. In simulated mixtures, four public cell line mixtures, and a targeted panel (1199 samples; 406 genes) from the Carolina Breast Cancer Study, DeCompress recapitulates cell-type proportions with less error than reference-free methods and finds biologically relevant compartments. We integrate compartment estimates into cis-eQTL mapping in breast cancer, identifying a tumor-specific cis-eQTL for CCR3 (C-C Motif Chemokine Receptor 3) at a risk locus. DeCompress improves upon reference-free methods without requiring expression profiles from pure cell populations, with applications in genomic analyses and clinical settings.

Nonlinear ridge regression improves cell-type-specific differential expression analysis

BACKGROUND

Epigenome-wide association studies (EWAS) and differential gene expression analyses are generally performed on tissue samples, which consist of multiple cell types. Cell-type-specific effects of a trait, such as disease, on the omics expression are of interest but difficult or costly to measure experimentally. By measuring omics data for the bulk tissue, cell type composition of a sample can be inferred statistically. Subsequently, cell-type-specific effects are estimated by linear regression that includes terms representing the interaction between the cell type proportions and the trait. This approach involves two issues, scaling and multicollinearity.

RESULTS

First, although cell composition is analyzed in linear scale, differential methylation/expression is analyzed suitably in the logit/log scale. To simultaneously analyze two scales, we applied nonlinear regression. Second, we show that the interaction terms are highly collinear, which is obstructive to ordinary regression. To cope with the multicollinearity, we applied ridge regularization. In simulated data, nonlinear ridge regression attained well-balanced sensitivity, specificity and precision. Marginal model attained the lowest precision and highest sensitivity and was the only algorithm to detect weak signal in real data.

CONCLUSION

Nonlinear ridge regression performed cell-type-specific association test on bulk omics data with well-balanced performance. The omicwas package for R implements nonlinear ridge regression for cell-type-specific EWAS, differential gene expression and QTL analyses. The software is freely available from https://github.com/fumi-github/omicwas.

Leveraging the cell lineage to predict cell-type specificity of regulatory variation from bulk genomics

Recent computational methods have enabled the inference of the cell-type-specificity of eQTLs based on bulk transcriptomes from highly heterogeneous tissues. However, these methods are limited in their scalability to highly heterogeneous tissues and limited in their broad applicability to any cell-type specificity of eQTLs. Here we present and demonstrate Cell Lineage Genetics (CeL-Gen), a novel computational approach that allows inference of eQTLs together with the subsets of cell types in which they have an effect, from bulk transcriptome data. To obtain improved scalability and broader applicability, CeL-Gen takes as input the known cell lineage tree and relies on the observation that dynamic changes in genetic effects occur relatively infrequently during cell differentiation. CeL-Gen can therefore be used not only to tease apart genetic effects derived from different cell types but also to infer the particular differentiation steps in which genetic effects are altered.

Correction for both common and rare cell types in blood is important to identify genes that correlate with age

BACKGROUND

Aging is a multifactorial process that affects multiple tissues and is characterized by changes in homeostasis over time, leading to increased morbidity. Whole blood gene expression signatures have been associated with aging and have been used to gain information on its biological mechanisms, which are still not fully understood. However, blood is composed of many cell types whose proportions in blood vary with age. As a result, previously observed associations between gene expression levels and aging might be driven by cell type composition rather than intracellular aging mechanisms. To overcome this, previous aging studies already accounted for major cell types, but the possibility that the reported associations are false positives driven by less prevalent cell subtypes remains.

RESULTS

Here, we compared the regression model from our previous work to an extended model that corrects for 33 additional white blood cell subtypes. Both models were applied to whole blood gene expression data from 3165 individuals belonging to the general population (age range of 18-81 years). We evaluated that the new model is a better fit for the data and it identified fewer genes associated with aging (625, compared to the 2808 of the initial model; P ≤ 2.5⨯10^-6). Moreover, 511 genes (~ 18% of the 2808 genes identified by the initial model) were found using both models, indicating that the other previously reported genes could be proxies for less abundant cell types. In particular, functional enrichment of the genes identified by the new model highlighted pathways and GO terms specifically associated with platelet activity.

CONCLUSIONS

We conclude that gene expression analyses in blood strongly benefit from correction for both common and rare blood cell types, and recommend using blood-cell count estimates as standard covariates when studying whole blood gene expression.

Resolving trained immunity with systems biology

Trained immunity is characterized by long-term functional reprogramming of innate immune cells following challenge with pathogens or microbial ligands during infection or vaccination. This cellular reprogramming leads to increased responsiveness upon restimulation, and is mediated through epigenetic and metabolic modifications. In this review, we describe how molecular mechanisms underlying trained immunity, for example, induced by β-glucan or Bacille Calmette-Guérin (BCG) vaccination, can be investigated by using and integrating different layers of information including genome, epigenome, transcriptome, proteome, metabolome, microbiome, immune cell phenotyping, and function. We also describe the most commonly used experimental and computational techniques. Finally, we provide a number of examples of how a systems biology approach was applied to study trained immunity to understand interindividual variation or the complex interplay between molecular layers. In conclusion, trained immunity represents an opportunity for regulating innate immune function, and understanding the complex interplay of mechanisms that mediate trained immunity might enable us to employ it as a clinical tool in the future.