rs953413 Regulates Polyunsaturated Fatty Acid Metabolism by Modulating ELOVL2 Expression

Dietary docosahexaenoic acid (DHA) downregulates liver DHA synthesis by inhibiting eicosapentaenoic acid elongation

DHA is abundant in the brain where it regulates cell survival, neurogenesis, and neuroinflammation. DHA can be obtained from the diet or synthesized from alpha-linolenic acid (ALA; 18:3n-3) via a series of desaturation and elongation reactions occurring in the liver. Tracer studies suggest that dietary DHA can downregulate its own synthesis, but the mechanism remains undetermined and is the primary objective of this manuscript. First, we show by tracing 13C content (δ13C) of DHA via compound-specific isotope analysis, that following low dietary DHA, the brain receives DHA synthesized from ALA. We then show that dietary DHA increases mouse liver and serum EPA, which is dependant on ALA. Furthermore, by compound-specific isotope analysis we demonstrate that the source of increased EPA is slowed EPA metabolism, not increased DHA retroconversion as previously assumed. DHA feeding alone or with ALA lowered liver elongation of very long chain (ELOVL2, EPA elongation) enzyme activity despite no change in protein content. To further evaluate the role of ELOVL2, a liver-specific Elovl2 KO was generated showing that DHA feeding in the presence or absence of a functional liver ELOVL2 yields similar results. An enzyme competition assay for EPA elongation suggests both uncompetitive and noncompetitive inhibition by DHA depending on DHA levels. To translate our findings, we show that DHA supplementation in men and women increases EPA levels in a manner dependent on a SNP (rs953413) in the ELOVL2 gene. In conclusion, we identify a novel feedback inhibition pathway where dietary DHA downregulates its liver synthesis by inhibiting EPA elongation.

FOXA1 and FOXA2: the regulatory mechanisms and therapeutic implications in cancer

FOXA1 (Forkhead Box A1) and FOXA2 (Forkhead Box A2) serve as pioneering transcription factors that build gene expression capacity and play a central role in biological processes, including organogenesis and differentiation, glycolipid metabolism, proliferation, migration and invasion, and drug resistance. Notably, FOXA1 and FOXA2 may exert antagonistic, synergistic, or complementary effects in the aforementioned biological processes. This article focuses on the molecular mechanisms and clinical relevance of FOXA1 and FOXA2 in steroid hormone-induced malignancies and highlights potential strategies for targeting FOXA1 and FOXA2 for cancer therapy. Furthermore, the article describes the prospect of targeting upstream regulators of FOXA1/FOXA2 to regulate its expression for cancer therapy because of the drug untargetability of FOXA1/FOXA2.

Human-genome single nucleotide polymorphisms affecting transcription factor binding and their role in pathogenesis

Single nucleotide polymorphisms (SNPs) are the most common type of variation in the human genome. The vast majority of SNPs identified in the human genome do not have any effect on the phenotype; however, some can lead to changes in the function of a gene or the level of its expression. Most SNPs associated with certain traits or pathologies are mapped to regulatory regions of the genome and affect gene expression by changing transcription factor binding sites. In recent decades, substantial effort has been invested in searching for such regulatory SNPs (rSNPs) and understanding the mechanisms by which they lead to phenotypic differences, primarily to individual differences in susceptibility to diseases and in sensitivity to drugs. The development of the NGS (next-generation sequencing) technology has contributed not only to the identification of a huge number of SNPs and to the search for their association (genome-wide association studies, GWASs) with certain diseases or phenotypic manifestations, but also to the development of more productive approaches to their functional annotation. It should be noted that the presence of an association does not allow one to identify a functional, truly disease-associated DNA sequence variant among multiple marker SNPs that are detected due to linkage disequilibrium. Moreover, determination of associations of genetic variants with a disease does not provide information about the functionality of these variants, which is necessary to elucidate the molecular mechanisms of the development of pathology and to design effective methods for its treatment and prevention. In this regard, the functional analysis of SNPs annotated in the GWAS catalog, both at the genome-wide level and at the level of individual SNPs, became especially relevant in recent years. A genome-wide search for potential rSNPs is possible without any prior knowledge of their association with a trait. Thus, mapping expression quantitative trait loci (eQTLs) makes it possible to identify an SNP for which - among transcriptomes of homozygotes and heterozygotes for its various alleles - there are differences in the expression level of certain genes, which can be located at various distances from the SNP. To predict rSNPs, approaches based on searches for allele-specific events in RNA-seq, ChIP-seq, DNase-seq, ATAC-seq, MPRA, and other data are also used. Nonetheless, for a more complete functional annotation of such rSNPs, it is necessary to establish their association with a trait, in particular, with a predisposition to a certain pathology or sensitivity to drugs. Thus, approaches to finding SNPs important for the development of a trait can be categorized into two groups: (1) starting from data on an association of SNPs with a certain trait, (2) starting from the determination of allele-specific changes at the molecular level (in a transcriptome or regulome). Only comprehensive use of strategically different approaches can considerably enrich our knowledge about the role of genetic determinants in the molecular mechanisms of trait formation, including predisposition to multifactorial diseases.

High-throughput functional dissection of noncoding SNPs with biased allelic enhancer activity for insulin resistance-relevant phenotypes

Most of the single-nucleotide polymorphisms (SNPs) associated with insulin resistance (IR)-relevant phenotypes by genome-wide association studies (GWASs) are located in noncoding regions, complicating their functional interpretation. Here, we utilized an adapted STARR-seq to evaluate the regulatory activities of 5,987 noncoding SNPs associated with IR-relevant phenotypes. We identified 876 SNPs with biased allelic enhancer activity effects (baaSNPs) across 133 loci in three IR-relevant cell lines (HepG2, preadipocyte, and A673), which showed pervasive cell specificity and significant enrichment for cell-specific open chromatin regions or enhancer-indicative markers (H3K4me1, H3K27ac). Further functional characterization suggested several transcription factors (TFs) with preferential allelic binding to baaSNPs. We also incorporated multi-omics data to prioritize 102 candidate regulatory target genes for baaSNPs and revealed prevalent long-range regulatory effects and cell-specific IR-relevant biological functional enrichment on them. Specifically, we experimentally verified the distal regulatory mechanism at IRS1 locus, in which rs952227-A reinforces IRS1 expression by long-range chromatin interaction and preferential binding to the transcription factor HOXC6 to augment the enhancer activity. Finally, based on our STARR-seq screening data, we predicted the enhancer activity of 227,343 noncoding SNPs associated with IR-relevant phenotypes (fasting insulin adjusted for BMI, HDL cholesterol, and triglycerides) from the largest available GWAS summary statistics. We further provided an open resource (http://www.bigc.online/fnSNP-IR) for better understanding genetic regulatory mechanisms of IR-relevant phenotypes.

Insight of a lipid metabolism prognostic model to identify immune landscape and potential target for retroperitoneal liposarcoma

Introduction

The exploration of lipid metabolism dysregulation may provide novel perspectives for retroperitoneal liposarcoma (RPLS). In our study, we aimed to investigate potential targets and facilitate further understanding of immune landscape in RPLS, through lipid metabolism-associated genes (LMAGs) based prognostic model.

Methods

Gene expression profiles and corresponding clinical information of 234 cases were enrolled from two public databases and the largest retroperitoneal tumor research center of East China, including cohort-TCGA (n=58), cohort-GSE30929 (n=92), cohort-FD (n=50), cohort-scRNA-seq (n=4) and cohort-validation (n=30). Consensus clustering analysis was performed to identify lipid metabolism-associated molecular subtypes (LMSs). A prognostic risk model containing 13 LMAGs was established using LASSO algorithm and multivariate Cox analysis in cohort-TCGA. ESTIMATE, CIBERSORT, XCELL and MCP analyses were performed to visualize the immune landscape. WGCNA was used to identify three hub genes among the 13 model LMAGs, and preliminarily validated in both cohort-GSE30929 and cohort-FD. Moreover, TIMER was used to visualize the correlation between antigen-presenting cells and potential targets. Finally, single-cell RNA-sequencing (scRNA-seq) analysis of four RPLS and multiplexed immunohistochemistry (mIHC) were performed in cohort-validation to validate the discoveries of bioinformatics analysis.

Results

LMS1 and LMS2 were characterized as immune-infiltrated and -excluded tumors, with significant differences in molecular features and clinical prognosis, respectively. Elongation of very long chain fatty acids protein 2 (ELOVL2), the enzyme that catalyzed the elongation of long chain fatty acids, involved in the maintenance of lipid metabolism and cellular homeostasis in normal cells, was identified and negatively correlated with antigen-presenting cells and identified as a potential target in RPLS. Furthermore, ELOVL2 was enriched in LMS2 with significantly lower immunoscore and unfavorable prognosis. Finally, a high-resolution dissection through scRNA-seq was performed in four RPLS, revealing the entire tumor ecosystem and validated previous findings.

Discussion

The LMS subgroups and risk model based on LMAGs proposed in our study were both promising prognostic classifications for RPLS. ELOVL2 is a potential target linking lipid metabolism to immune regulations against RPLS, specifically for patients with LMS2 tumors.

Age-associated changes in circulatory fatty acids: new insights on adults and long-lived individuals

Long-lived individuals (LLIs) are considered an ideal model to study healthy human aging. Blood fatty acid (FA) profile of a cohort of LLIs (90-111 years old, n = 49) from Sicily was compared to adults (18-64 years old, n = 69) and older adults (65-89 years old, n = 54) from the same area. Genetic variants in key enzymes related to FA biosynthesis and metabolism were also genotyped to investigate a potential genetic predisposition in determining the FA profile. Gas chromatography was employed to determine the FA profile, and genotyping was performed using high-resolution melt (HRM) analysis. Blood levels of total polyunsaturated FA (PUFA) and total trans-FA decreased with age, while the levels of saturated FA (SFA) remained unchanged. Interestingly, distinctively higher circulatory levels of monounsaturated FA (MUFA) in LLIs compared to adults and older adults were observed. In addition, among LLIs, rs174537 in the FA desaturase 1/2 (FADS1/2) gene was associated with linoleic acid (LA, 18:2n-6) and docosatetraenoic acid (DTA, 22:4n-6) levels, and the rs953413 in the elongase of very long FA 2 (ELOVL2) was associated with DTA levels. We further observed that rs174579 and rs174626 genotypes in FADS1/2 significantly affect delta-6 desaturase (D6D) activity. In conclusion, our results suggest that the LLIs have a different FA profile characterized by high MUFA content, which indicates reduced peroxidation while maintaining membrane fluidity.

Increases in plasma n-3 tetracosapentaenoic acid and tetracosahexaenoic acid following 12 weeks of EPA, but not DHA, supplementation in women and men

Dietary feeding and stable isotope studies in rodents support that the 24-carbon omega-3 polyunsaturated fatty acids, tetracosapentaenoic acid (24:5n-3, TPAn-3) and tetracosahexaenoic acid (24:6n-3, THA), are immediate precursors to docosahexaenoic acid (DHA, 22:6n-3). In this study, we assessed for the first time, changes in TPAn-3 or THA levels following omega-3 PUFA supplementation in humans, providing insight into human omega-3 PUFA metabolism. In this secondary analysis of a double-blind randomized control trial, women and men (19 - 30 years, n = 10 - 14 per sex, per diet) were supplemented with 3 g/day EPA, DHA, or olive oil control for 12 weeks. Plasma TPAn-3 and THA concentrations were determined by gas chromatography-mass spectrometry to determine changes following supplementation in a sex-specific manner (sex x time). EPA supplementation significantly increased (p < 0.0001) plasma TPAn-3 by 215% (1.3 ± 0.1 - 4.1 ± 0.7, nmol/mL ± SEM) and THA by 112% (1.7 ± 0.2 - 3.6 ± 0.5, nmol/mL ± SEM). Furthermore, women had 111% and 99% higher plasma TPAn-3 and THA in the EPA supplemented group compared to men (p < 0.0001). There were no significant effects of time on plasma TPAn-3 or THA concentrations in the DHA supplemented or olive oil supplemented groups. In conclusion, EPA, but not DHA, supplementation in humans increased plasma TPAn-3 and THA levels, suggesting that THA accumulates prior to conversion to DHA in the n-3 PUFA synthesis pathway. Furthermore, women generally exhibit higher plasma TPAn-3 and THA concentrations compared with men, suggesting that women have a greater ability to accumulate 24-carbon n-3 PUFA in plasma via EPA and DPAn-3 elongation, which may explain the known higher DHA levels in women. Summary: In this secondary analysis of a double-blind randomized control trial, we assessed changes in omega-3 (n-3) tetracosapentaenoic acid (24:5n-3, TPAn-3) and tetracosahexaenoic acid (24:6n-3, THA) plasma levels in women and men (19 - 30 years, n = 10 - 14 per sex, per diet) following 12-weeks of n-3 PUFA supplementation (3 g/day EPA, DHA or olive oil). Women had higher plasma TPAn-3 in all supplementation groups and higher THA levels in the EPA and olive oil groups (p < 0.0001) compared to men. EPA supplementation increased (p < 0.0001) plasma TPAn-3 by 215% (1.3 ± 0.1 - 4.1 ± 0.7, nmol/mL ± SEM) and THA by 112% (1.7 ± 0.2 - 3.6 ± 0.5, nmol/mL ± SEM), but DHA supplementation had no effect. For the first time in humans, we show that plasma TPAn-3 and THA levels are higher in women and increased with EPA, but not DHA supplementation, suggesting an accumulation of THA prior to conversion to DHA in the n-3 PUFA synthesis pathway.

The risk variant rs11836367 contributes to breast cancer onset and metastasis by attenuating Wnt signaling via regulating NTN4 expression

Most genome-wide association study (GWAS)-identified breast cancer-associated causal variants remain uncharacterized. To provide a framework of understanding GWAS-identified variants to function, we performed a comprehensive study of noncoding regulatory variants at the NTN4 locus (12q22) and NTN4 gene in breast cancer etiology. We find that rs11836367 is the more likely causal variant, disrupting enhancer activity in both enhancer reporter assays and endogenous genome editing experiments. The protective T allele of rs11837367 increases the binding of GATA3 to the distal enhancer and up-regulates NTN4 expression. In addition, we demonstrate that loss of NTN4 gene in mice leads to tumor earlier onset, progression, and metastasis. We discover that NTN4, as a tumor suppressor, can attenuate the Wnt signaling pathway by directly binding to Wnt ligands. Our findings bridge the gaps among breast cancer-associated single-nucleotide polymorphisms, transcriptional regulation of NTN4, and breast cancer biology, which provides previously unidentified insights into breast cancer prediction and prevention.

CRISPR/Cas9-Mediated Transgenesis of the Masu Salmon (Oncorhynchus masou) elovl2 Gene Improves n-3 Fatty Acid Content in Channel Catfish (Ictalurus punctatus)

Omega-3 polyunsaturated fatty acids (n-3 PUFAs), particularly eicosapentaenoic acid (EPA, 20:5n-3) and docosahexaenoic acid (DHA, 22:6n-3), play a very important role in human health. Channel catfish (Ictalurus punctatus) is one of the leading freshwater aquaculture species in the USA, but has low levels of EPA and DHA compared to some fish such as salmon. To improve EPA and DHA content, a modification of the n-3 PUFA biosynthetic pathway was achieved through the insertion of an elovl2 transgene isolated from masu salmon (Oncorhynchus masou) driven by a carp β-actin promoter using a two-hit by gRNA and two oligos with a targeting plasmid (2H2OP) CRISPR/Cas9 approach. Integration rate of the transgene was high (37.5%) and detected in twelve different tissues of P₁ transgenic fish with tissue-specific gene expression. Liver and muscle had relative high gene expression (13.4- and 9.2-fold change, respectively). Fatty acid analysis showed DHA content in the muscle from transgenic fish was 1.62-fold higher than in non-transgenic fish (P < 0.05). Additionally, total n-3 PUFAs and omega-6 polyunsaturated fatty acids (n-6 PUFAs) increased to 1.41-fold and 1.50-fold, respectively, suggesting the β-actin-elovl2 transgene improved biosynthesis of PUFAs in channel catfish as a whole. The n-9 fatty acid level decreased in the transgenic fish compared to the control. Morphometric analysis showed that there were significant differences between injected fish with sgRNAs (including positive and negative fish) and sham-injected controls (P < 0.001). Potential off-target effects are likely the major factor responsible for morphological deformities. Optimization of sgRNA design to maximize activity and reduce off-target effects of CRISPR/Cas9 should be examined in future transgenic research, but this research shows a promising first step in the improvement of n-3 PUFAs in channel catfish.

The landscape of GWAS validation; systematic review identifying 309 validated non-coding variants across 130 human diseases

Background

The remarkable growth of genome-wide association studies (GWAS) has created a critical need to experimentally validate the disease-associated variants, 90% of which involve non-coding variants.

Methods

To determine how the field is addressing this urgent need, we performed a comprehensive literature review identifying 36,676 articles. These were reduced to 1454 articles through a set of filters using natural language processing and ontology-based text-mining. This was followed by manual curation and cross-referencing against the GWAS catalog, yielding a final set of 286 articles.

Results

We identified 309 experimentally validated non-coding GWAS variants, regulating 252 genes across 130 human disease traits. These variants covered a variety of regulatory mechanisms. Interestingly, 70% (215/309) acted through cis-regulatory elements, with the remaining through promoters (22%, 70/309) or non-coding RNAs (8%, 24/309). Several validation approaches were utilized in these studies, including gene expression (n = 272), transcription factor binding (n = 175), reporter assays (n = 171), in vivo models (n = 104), genome editing (n = 96) and chromatin interaction (n = 33).

Conclusions

This review of the literature is the first to systematically evaluate the status and the landscape of experimentation being used to validate non-coding GWAS-identified variants. Our results clearly underscore the multifaceted approach needed for experimental validation, have practical implications on variant prioritization and considerations of target gene nomination. While the field has a long way to go to validate the thousands of GWAS associations, we show that progress is being made and provide exemplars of validation studies covering a wide variety of mechanisms, target genes, and disease areas.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12920-022-01216-w.

A common variant in 11q23.3 associated with hyperlipidemia is mediated by the binding and regulation of GATA4

Large-scale genome-wide associations comprising multiple studies have identified hundreds of genetic loci commonly associated with hyperlipidemia-related phenotypes. However, single large cohort remains necessary in aiming to investigate ethnicity-specific genetic risks and mechanical insights. A community-based cohort comprising 23,988 samples that included both genotype and biochemical information was assembled for the genome-wide association analysis (GWAS) of hyperlipidemia. The analysis identified fifty genetic variants (P < 5 × 10⁻⁸) on five different chromosomes, and a subsequent validation analysis confirmed the significance of the lead variants. Integrated analysis combined with cell-based experiments of the most statistically significant locus in 11q23.3 revealed rs651821 (P = 4.52 × 10⁻⁷⁶) as the functional variant. We showed transcription factor GATA4 preferentially binds the T allele of rs651821, the protective allele for hyperlipidemia, which promoted APOA5 expression in liver cells and individuals with the TT genotype of rs651821. As GATA4-APOA5 axis maintains triglyceride homeostasis, GATA4 activation by phenylephrine implies synergism for lowering triglyceride levels in hyperlipidemia patients. Our study demonstrates that rs651821 mediates APOA5 activation via allele-specific regulation by GATA4. We suggest elevating GATA4 activity could provide a therapeutic potential for treating the development of hyperlipidemia.

Nutrigenetics and nutrigenomics-A personalized approach to nutrition

The prevalence of non-communicable diseases has been on an upward trajectory for some time and this puts an enormous burden on the healthcare expenditure. Lifestyle modifications including dietary interventions hold an immense promise to manage and prevent these diseases. Recent advances in genomic research provide evidence that focussing these efforts on individual variations in abilities to metabolize nutrients (nutrigenetics) and exploring the role of dietary compounds on gene expression (nutrigenomics and nutri-epigenomics) can lead to more meaningful personalized dietary strategies to promote optimal health. This chapter aims to provide examples on these gene-diet interactions at multiple levels to support the need of embedding targeted dietary interventions as a way forward to prevent, avoid and manage diseases.

Polymorphisms rs55710213 and rs56334587 regulate SCD1 expression by modulating HNF4A binding

The stearoyl-CoA desaturase 1 (SCD1) gene at 10q24.31 encodes the rate limiting enzyme SCD1 that catalyzes the biosynthesis of monounsaturated fatty acids (MUFAs) from saturated fatty acids (SFAs). Dysregulated SCD1 activity has been observed in many human diseases including non-alcoholic fatty liver disease (NAFLD), obesity, hypertension, hyperlipidemia, metabolic syndrome and several types of cancer. HNF4A is a central regulator of glucose and lipid metabolism and previous studies suggested that it is deeply involved in regulating the SCD1 activity in the liver. However, the underlying mechanisms on whether and how SCD1 is regulated by HNF4A have not been explored in detail. In this study, we found that HNF4A regulates SCD1 expression by directly binding to the key regulatory regions in the SCD1 locus. Knocking down of HNF4A significantly downregulated the expression of SCD1. Variants rs55710213 and rs56334587 in intron 5 of SCD1 directly reside in a canonical HNF4A binding site. The GG haplotype of rs55710213 and rs56334587 is associated with decreased SCD1 activity by disrupting the binding of HNF4A, which further decreased the enhancer activity and SCD1 expression. In conclusion, our study demonstrated that SCD1 is directly regulated by HNF4A, which may be helpful in the understanding of the altered metabolic pathways in many diseases associated with dysregulated SCD1 or HNF4A or both.

Regulatory SNPs: Altered Transcription Factor Binding Sites Implicated in Complex Traits and Diseases

The vast majority of the genetic variants (mainly SNPs) associated with various human traits and diseases map to a noncoding part of the genome and are enriched in its regulatory compartment, suggesting that many causal variants may affect gene expression. The leading mechanism of action of these SNPs consists in the alterations in the transcription factor binding via creation or disruption of transcription factor binding sites (TFBSs) or some change in the affinity of these regulatory proteins to their cognate sites. In this review, we first focus on the history of the discovery of regulatory SNPs (rSNPs) and systematized description of the existing methodical approaches to their study. Then, we brief the recent comprehensive examples of rSNPs studied from the discovery of the changes in the TFBS sequence as a result of a nucleotide substitution to identification of its effect on the target gene expression and, eventually, to phenotype. We also describe state-of-the-art genome-wide approaches to identification of regulatory variants, including both making molecular sense of genome-wide association studies (GWAS) and the alternative approaches the primary goal of which is to determine the functionality of genetic variants. Among these approaches, special attention is paid to expression quantitative trait loci (eQTLs) analysis and the search for allele-specific events in RNA-seq (ASE events) as well as in ChIP-seq, DNase-seq, and ATAC-seq (ASB events) data.

Multifaceted regulation of hepatic lipid metabolism by YY1

This study shows that YY1 regulates hepatic lipid metabolism by directly or indirectly regulating the expression of several key upstream transcription factors and their coactivators.

Recent studies suggested that dysregulated YY1 plays a pivotal role in many liver diseases. To obtain a detailed view of genes and pathways regulated by YY1 in the liver, we carried out RNA sequencing in HepG2 cells after YY1 knockdown. A rigid set of 2,081 differentially expressed genes was identified by comparing the YY1-knockdown samples (n = 8) with the control samples (n = 14). YY1 knockdown significantly decreased the expression of several key transcription factors and their coactivators in lipid metabolism. This is illustrated by YY1 regulating PPARA expression through binding to its promoter and enhancer regions. Our study further suggest that down-regulation of the key transcription factors together with YY1 knockdown significantly decreased the cooperation between YY1 and these transcription factors at various regulatory regions, which are important in regulating the expression of genes in hepatic lipid metabolism. This was supported by the finding that the expression of SCD and ELOVL6, encoding key enzymes in lipogenesis, were regulated by the cooperation between YY1 and PPARA/RXRA complex over their promoters.

Harnessing the potential of CRISPR-based platforms to advance the field of hospital medicine

INTRODUCTION

Clustered regularly interspaced short palindromic repeats (CRISPR) are segments of nucleic acid that play a role in prokaryotic defense and form the basis of a genome editing technology that allows permanent alteration of genetic material. This methodology, known as CRISPR-Cas9, is poised to revolutionize molecular biology, but no literature yet exists on how these advances will affect hospitalists.

AREAS COVERED

These specialists in inpatient medicine care for a wide variety of hospitalized patients, including those with infectious disease, cancer, cardiovascular disease, autoimmune disease, hematologic disease, and a variety of other conditions that may soon be impacted by advances in gene-modifying technology provided by CRISPR-Cas9. A Literature search was performed using PubMed [1 December 2019-17 April 2020].

EXPERT OPINION

This paper reviews the remarkable diagnostic and therapeutic potential of the CRISPR-Cas9 platform and concludes with a look at ethical issues and technical hurdles pertaining to the implementation of permanent gene modification in the practice of Hospital Medicine.

New avenues for systematically inferring cell-cell communication: through single-cell transcriptomics data

For multicellular organisms, cell-cell communication is essential to numerous biological processes. Drawing upon the latest development of single-cell RNA-sequencing (scRNA-seq), high-resolution transcriptomic data have deepened our understanding of cellular phenotype heterogeneity and composition of complex tissues, which enables systematic cell-cell communication studies at a single-cell level. We first summarize a common workflow of cell-cell communication study using scRNA-seq data, which often includes data preparation, construction of communication networks, and result validation. Two common strategies taken to uncover cell-cell communications are reviewed, e.g., physically vicinal structure-based and ligand-receptor interaction-based one. To conclude, challenges and current applications of cell-cell communication studies at a single-cell resolution are discussed in details and future perspectives are proposed.