The X-linked epigenetic regulator UTX controls NK cell-intrinsic sex differences.

Viral infection outcomes are sex biased, with males generally more susceptible than females. Paradoxically, the numbers of antiviral natural killer (NK) cells are increased in males. We demonstrate that while numbers of NK cells are increased in male mice, they display decreased effector function compared to females in mice and humans. These differences were not solely dependent on gonadal hormones, because they persisted in gonadectomized mice. Kdm6a (which encodes the protein UTX), an epigenetic regulator that escapes X inactivation, was lower in male NK cells, while NK cell-intrinsic UTX deficiency in female mice increased NK cell numbers and reduced effector responses. Furthermore, mice with NK cell-intrinsic UTX deficiency showed increased lethality to mouse cytomegalovirus. Integrative multi-omics analysis revealed a critical role for UTX in regulating chromatin accessibility and gene expression critical for NK cell homeostasis and effector function. Collectively, these data implicate UTX as a critical molecular determinant of sex differences in NK cells.

Type V Collagen in Scar Tissue Regulates the Size of Scar after Heart Injury.

Scar tissue size following myocardial infarction is an independent predictor of cardiovascular outcomes, yet little is known about factors regulating scar size. We demonstrate that collagen V, a minor constituent of heart scars, regulates the size of heart scars after ischemic injury. Depletion of collagen V led to a paradoxical increase in post-infarction scar size with worsening of heart function. A systems genetics approach across 100 in-bred strains of mice demonstrated that collagen V is a critical driver of postinjury heart function. We show that collagen V deficiency alters the mechanical properties of scar tissue, and altered reciprocal feedback between matrix and cells induces expression of mechanosensitive integrins that drive fibroblast activation and increase scar size. Cilengitide, an inhibitor of specific integrins, rescues the phenotype of increased post-injury scarring in collagen-V-deficient mice. These observations demonstrate that collagen V regulates scar size in an integrin-dependent manner.

The transcription factor Fli1 restricts the formation of memory precursor NK cells during viral infection.

Natural killer (NK) cells are innate lymphocytes that possess traits of adaptive immunity, such as memory formation. However, the molecular mechanisms by which NK cells persist to form memory cells are not well understood. Using single-cell RNA sequencing, we identified two distinct effector NK cell (NKeff) populations following mouse cytomegalovirus infection. Ly6C- memory precursor (MP) NK cells showed enhanced survival during the contraction phase in a Bcl2-dependent manner, and differentiated into Ly6C+ memory NK cells. MP NK cells exhibited distinct transcriptional and epigenetic signatures compared with Ly6C+ NKeff cells, with a core epigenetic signature shared with MP CD8+ T cells enriched in ETS1 and Fli1 DNA-binding motifs. Fli1 was induced by STAT5 signaling ex vivo, and increased levels of the pro-apoptotic factor Bim in early effector NK cells following viral infection. These results suggest that a NK cell-intrinsic checkpoint controlled by the transcription factor Fli1 limits MP NK formation by regulating early effector NK cell fitness during viral infection.

Single-cell sequencing of human white adipose tissue identifies new cell states in health and obesity.

White adipose tissue (WAT) is an essential regulator of energy storage and systemic metabolic homeostasis. Regulatory networks consisting of immune and structural cells are necessary to maintain WAT metabolism, which can become impaired during obesity in mammals. Using single-cell transcriptomics and flow cytometry, we unveil a large-scale comprehensive cellular census of the stromal vascular fraction of healthy lean and obese human WAT. We report new subsets and developmental trajectories of adipose-resident innate lymphoid cells, dendritic cells and monocyte-derived macrophage populations that accumulate in obese WAT. Analysis of cell-cell ligand-receptor interactions and obesity-enriched signaling pathways revealed a switch from immunoregulatory mechanisms in lean WAT to inflammatory networks in obese WAT. These results provide a detailed and unbiased cellular landscape of homeostatic and inflammatory circuits in healthy human WAT.

The cellular architecture of the antimicrobial response network in human leprosy granulomas.

Granulomas are complex cellular structures composed predominantly of macrophages and lymphocytes that function to contain and kill invading pathogens. Here, we investigated the single-cell phenotypes associated with antimicrobial responses in human leprosy granulomas by applying single-cell and spatial sequencing to leprosy biopsy specimens. We focused on reversal reactions (RRs), a dynamic process whereby some patients with disseminated lepromatous leprosy (L-lep) transition toward self-limiting tuberculoid leprosy (T-lep), mounting effective antimicrobial responses. We identified a set of genes encoding proteins involved in antimicrobial responses that are differentially expressed in RR versus L-lep lesions and regulated by interferon-γ and interleukin-1ß. By integrating the spatial coordinates of the key cell types and antimicrobial gene expression in RR and T-lep lesions, we constructed a map revealing the organized architecture of granulomas depicting compositional and functional layers by which macrophages, T cells, keratinocytes and fibroblasts can each contribute to the antimicrobial response.

Functional reprogramming of regulatory T cells in the absence of Foxp3.

Regulatory T cells (Treg cells) deficient in the transcription factor Foxp3 lack suppressor function and manifest an effector T (Teff) cell-like phenotype. We demonstrate that Foxp3 deficiency dysregulates metabolic checkpoint kinase mammalian target of rapamycin (mTOR) complex 2 (mTORC2) signaling and gives rise to augmented aerobic glycolysis and oxidative phosphorylation. Specific deletion of the mTORC2 adaptor gene Rictor in Foxp3-deficient Treg cells ameliorated disease in a Foxo1 transcription factor-dependent manner. Rictor deficiency re-established a subset of Treg cell genetic circuits and suppressed the Teff cell-like glycolytic and respiratory programs, which contributed to immune dysregulation. Treatment of Treg cells from patients with FOXP3 deficiency with mTOR inhibitors similarly antagonized their Teff cell-like program and restored suppressive function. Thus, regulatory function can be re-established in Foxp3-deficient Treg cells by targeting their metabolic pathways, providing opportunities to restore tolerance in Treg cell disorders.

Pioneer transcription factors target partial DNA motifs on nucleosomes to initiate reprogramming.

Pioneer transcription factors (TFs) access silent chromatin and initiate cell-fate changes, using diverse types of DNA binding domains (DBDs). FoxA, the paradigm pioneer TF, has a winged helix DBD that resembles linker histone and thereby binds its target sites on nucleosomes and in compacted chromatin. Herein, we compare the nucleosome and chromatin targeting activities of Oct4 (POU DBD), Sox2 (HMG box DBD), Klf4 (zinc finger DBD), and c-Myc (bHLH DBD), which together reprogram somatic cells to pluripotency. Purified Oct4, Sox2, and Klf4 proteins can bind nucleosomes in vitro, and in vivo they preferentially target silent sites enriched for nucleosomes. Pioneer activity relates simply to the ability of a given DBD to target partial motifs displayed on the nucleosome surface. Such partial motif recognition can occur by coordinate binding between factors. Our findings provide insight into how pioneer factors can target naive chromatin sites.

Mapping human haematopoietic stem cells from haemogenic endothelium to birth.

The ontogeny of human haematopoietic stem cells (HSCs) is poorly defined owing to the inability to identify HSCs as they emerge and mature at different haematopoietic sites1. Here we created a single-cell transcriptome map of human haematopoietic tissues from the first trimester to birth and found that the HSC signature RUNX1+HOXA9+MLLT3+MECOM+HLF+SPINK2+ distinguishes HSCs from progenitors throughout gestation. In addition to the aorta-gonad-mesonephros region, nascent HSCs populated the placenta and yolk sac before colonizing the liver at 6 weeks. A comparison of HSCs at different maturation stages revealed the establishment of HSC transcription factor machinery after the emergence of HSCs, whereas their surface phenotype evolved throughout development. The HSC transition to the liver marked a molecular shift evidenced by suppression of surface antigens reflecting nascent HSC identity, and acquisition of the HSC maturity markers CD133 (encoded by PROM1) and HLA-DR. HSC origin was tracked to ALDH1A1+KCNK17+ haemogenic endothelial cells, which arose from an IL33+ALDH1A1+ arterial endothelial subset termed pre-haemogenic endothelial cells. Using spatial transcriptomics and immunofluorescence, we visualized this process in ventrally located intra-aortic haematopoietic clusters. The in vivo map of human HSC ontogeny validated the generation of aorta-gonad-mesonephros-like definitive haematopoietic stem and progenitor cells from human pluripotent stem cells, and serves as a guide to improve their maturation to functional HSCs.

Challenges and perspectives in computational deconvolution of genomics data.

Deciphering cell-type heterogeneity is crucial for systematically understanding tissue homeostasis and its dysregulation in diseases. Computational deconvolution is an efficient approach for estimating cell-type abundances from a variety of omics data. Despite substantial methodological progress in computational deconvolution in recent years, challenges are still outstanding. Here we enlist four important challenges related to computational deconvolution: the quality of the reference data, generation of ground truth data, limitations of computational methodologies, and benchmarking design and implementation. Finally, we make recommendations on reference data generation, new directions of computational methodologies, and strategies to promote rigorous benchmarking.

Unraveling inflammatory responses using systems genetics and gene-environment interactions in macrophages.

Many common diseases have an important inflammatory component mediated in part by macrophages. Here we used a systems genetics strategy to examine the role of common genetic variation in macrophage responses to inflammatory stimuli. We examined genome-wide transcript levels in macrophages from 92 strains of the Hybrid Mouse Diversity Panel. We exposed macrophages to control media, bacterial lipopolysaccharide (LPS), or oxidized phospholipids. We performed association mapping under each condition and identified several thousand expression quantitative trait loci (eQTL), gene-by-environment interactions, and eQTL "hot spots" that specifically control LPS responses. We used siRNA knockdown of candidate genes to validate an eQTL hot spot in chromosome 8 and identified the gene 2310061C15Rik as a regulator of inflammatory responses in macrophages. We have created a public database where the data presented here can be used as a resource for understanding many common inflammatory traits that are modeled in the mouse and for the dissection of regulatory relationships between genes.

Scl represses cardiomyogenesis in prospective hemogenic endothelium and endocardium.

Endothelium in embryonic hematopoietic tissues generates hematopoietic stem/progenitor cells; however, it is unknown how its unique potential is specified. We show that transcription factor Scl/Tal1 is essential for both establishing the hematopoietic transcriptional program in hemogenic endothelium and preventing its misspecification to a cardiomyogenic fate. Scl(-/-) embryos activated a cardiac transcriptional program in yolk sac endothelium, leading to the emergence of CD31+Pdgfrα+ cardiogenic precursors that generated spontaneously beating cardiomyocytes. Ectopic cardiogenesis was also observed in Scl(-/-) hearts, where the disorganized endocardium precociously differentiated into cardiomyocytes. Induction of mosaic deletion of Scl in Scl(fl/fl)Rosa26Cre-ER(T2) embryos revealed a cell-intrinsic, temporal requirement for Scl to prevent cardiomyogenesis from endothelium. Scl(-/-) endothelium also upregulated the expression of Wnt antagonists, which promoted rapid cardiomyocyte differentiation of ectopic cardiogenic cells. These results reveal unexpected plasticity in embryonic endothelium such that loss of a single master regulator can induce ectopic cardiomyogenesis from endothelial cells.

Adenovirus small E1A directs activation of Alu transcription at YAP/TEAD- and AP-1-bound enhancers through interactions with the EP400 chromatin remodeler.

Alu retrotransposons, which form the largest family of mobile DNA elements in the human genome, have recently come to attention as a potential source of regulatory novelties, most notably by participating in enhancer function. Even though Alu transcription by RNA polymerase III is subjected to tight epigenetic silencing, their expression has long been known to increase in response to various types of stress, including viral infection. Here we show that, in primary human fibroblasts, adenovirus small e1a triggered derepression of hundreds of individual Alus by promoting TFIIIB recruitment by Alu-bound TFIIIC. Epigenome profiling revealed an e1a-induced decrease of H3K27 acetylation and increase of H3K4 monomethylation at derepressed Alus, making them resemble poised enhancers. The enhancer nature of e1a-targeted Alus was confirmed by the enrichment, in their upstream regions, of the EP300/CBP acetyltransferase, EP400 chromatin remodeler and YAP1 and FOS transcription factors. The physical interaction of e1a with EP400 was critical for Alu derepression, which was abrogated upon EP400 ablation. Our data suggest that e1a targets a subset of enhancer Alus whose transcriptional activation, which requires EP400 and is mediated by the e1a-EP400 interaction, may participate in the manipulation of enhancer activity by adenoviruses.

Single-stranded pre-methylated 5mC adapters uncover the methylation profile of plasma ultrashort Single-stranded cell-free DNA.

Whole-genome bisulfite sequencing (BS-Seq) measures cytosine methylation changes at single-base resolution and can be used to profile cell-free DNA (cfDNA). In plasma, ultrashort single-stranded cfDNA (uscfDNA, ∼50 nt) has been identified together with 167 bp double-stranded mononucleosomal cell-free DNA (mncfDNA). However, the methylation profile of uscfDNA has not been described. Conventional BS-Seq workflows may not be helpful because bisulfite conversion degrades larger DNA into smaller fragments, leading to erroneous categorization as uscfDNA. We describe the '5mCAdpBS-Seq' workflow in which pre-methylated 5mC (5-methylcytosine) single-stranded adapters are ligated to heat-denatured cfDNA before bisulfite conversion. This method retains only DNA fragments that are unaltered by bisulfite treatment, resulting in less biased uscfDNA methylation analysis. Using 5mCAdpBS-Seq, uscfDNA had lower levels of DNA methylation (∼15%) compared to mncfDNA and was enriched in promoters and CpG islands. Hypomethylated uscfDNA fragments were enriched in upstream transcription start sites (TSSs), and the intensity of enrichment was correlated with expressed genes of hemopoietic cells. Using tissue-of-origin deconvolution, we inferred that uscfDNA is derived primarily from eosinophils, neutrophils, and monocytes. As proof-of-principle, we show that characteristics of the methylation profile of uscfDNA can distinguish non-small cell lung carcinoma from non-cancer samples. The 5mCAdpBS-Seq workflow is recommended for any cfDNA methylation-based investigations.

Heterogeneous pdgfrb+ cells regulate coronary vessel development and revascularization during heart regeneration.

Endothelial cells emerge from the atrioventricular canal to form coronary blood vessels in juvenile zebrafish hearts. We find that pdgfrb is first expressed in the epicardium around the atrioventricular canal and later becomes localized mainly in the mural cells. pdgfrb mutant fish show severe defects in mural cell recruitment and coronary vessel development. Single-cell RNA sequencing analyses identified pdgfrb+ cells as epicardium-derived cells (EPDCs) and mural cells. Mural cells associated with coronary arteries also express cxcl12b and smooth muscle cell markers. Interestingly, these mural cells remain associated with coronary arteries even in the absence of Pdgfrß, although smooth muscle gene expression is downregulated. We find that pdgfrb expression dynamically changes in EPDCs of regenerating hearts. Differential gene expression analyses of pdgfrb+ EPDCs and mural cells suggest that they express genes that are important for regeneration after heart injuries. mdka was identified as a highly upregulated gene in pdgfrb+ cells during heart regeneration. However, pdgfrb but not mdka mutants show defects in heart regeneration after amputation. Our results demonstrate that heterogeneous pdgfrb+ cells are essential for coronary development and heart regeneration.

Non-coding RNA yREX3 from human extracellular vesicles exerts macrophage-mediated cardioprotection via a novel gene-methylating mechanism.

BACKGROUND AND AIMS: Extracellular vesicles (EVs) secreted by cardiosphere-derived cells exert immunomodulatory effects through the transmission of small non-coding RNAs. METHODS: The mechanism and role of yREX3, a small Y RNA abundant in EVs in myocardial injury, was investigated. RESULTS: yREX3 attenuates cardiac ischaemic injury by selective DNA methylation. Synthetic yREX3 encapsulated in lipid nanoparticles triggers broad transcriptomic changes in macrophages, localizes to the nucleus, and mediates epigenetic silencing of protein interacting with C kinase-1 (Pick1) through methylation of upstream CpG sites. Moreover, yREX3 interacts with polypyrimidine tract binding protein 3 (PTBP3) to methylate the Pick1 gene locus in a DNA methyltransferase-dependent manner. Suppression of Pick1 in macrophages potentiates Smad3 signalling and enhances efferocytosis, minimizing heart necrosis in rats with myocardial infarction. Adoptive transfer of Pick1-deficient macrophages recapitulates the cardioprotective effects of yREX3 in vivo. CONCLUSIONS: These findings highlight the role of a small Y RNA mined from EVs with a novel gene-methylating mechanism.

Single-cell RNA sequencing of batch Chlamydomonas cultures reveals heterogeneity in their diurnal cycle phase.

The photosynthetic unicellular alga Chlamydomonas (Chlamydomonas reinhardtii) is a versatile reference for algal biology because of its ease of culture in the laboratory. Genomic and systems biology approaches have previously described transcriptome responses to environmental changes using bulk data, thus representing the average behavior from pools of cells. Here, we apply single-cell RNA sequencing (scRNA-seq) to probe the heterogeneity of Chlamydomonas cell populations under three environments and in two genotypes differing by the presence of a cell wall. First, we determined that RNA can be extracted from single algal cells with or without a cell wall, offering the possibility to sample natural algal communities. Second, scRNA-seq successfully separated single cells into nonoverlapping cell clusters according to their growth conditions. Cells exposed to iron or nitrogen deficiency were easily distinguished despite a shared tendency to arrest photosynthesis and cell division to economize resources. Notably, these groups of cells not only recapitulated known patterns observed with bulk RNA-seq but also revealed their inherent heterogeneity. A substantial source of variation between cells originated from their endogenous diurnal phase, although cultures were grown in constant light. We exploited this result to show that circadian iron responses may be conserved from algae to land plants. We document experimentally that bulk RNA-seq data represent an average of typically hidden heterogeneity in the population.

The induction of preterm labor in rhesus macaques is determined by the strength of immune response to intrauterine infection.

Intrauterine infection/inflammation (IUI) is a major contributor to preterm labor (PTL). However, IUI does not invariably cause PTL. We hypothesized that quantitative and qualitative differences in immune response exist in subjects with or without PTL. To define the triggers for PTL, we developed rhesus macaque models of IUI driven by lipopolysaccharide (LPS) or live Escherichia coli. PTL did not occur in LPS challenged rhesus macaques, while E. coli-infected animals frequently delivered preterm. Although LPS and live E. coli both caused immune cell infiltration, E. coli-infected animals showed higher levels of inflammatory mediators, particularly interleukin 6 (IL-6) and prostaglandins, in the chorioamnion-decidua and amniotic fluid (AF). Neutrophil infiltration in the chorio-decidua was a common feature to both LPS and E. coli. However, neutrophilic infiltration and IL6 and PTGS2 expression in the amnion was specifically induced by live E. coli. RNA sequencing (RNA-seq) analysis of fetal membranes revealed that specific pathways involved in augmentation of inflammation including type I interferon (IFN) response, chemotaxis, sumoylation, and iron homeostasis were up-regulated in the E. coli group compared to the LPS group. Our data suggest that the intensity of the host immune response to IUI may determine susceptibility to PTL.

Exponential dynamics of DNA methylation with age.

The association of DNA methylation with age has been extensively studied. Previous work has investigated the trajectories of methylation with age, and developed predictive biomarkers of age. However, we still have a limited understanding of the functional form of methylation-age dynamics. To address this we present a theoretical framework to model the dynamics of DNA methylation at single sites. We show that this model leads to convergence to a steady-state methylation level at an exponential rate. By fitting the model to a dataset that measures changes in DNA methylation in the brain from birth to old age, we show that the timescales of this exponential convergence are heterogeneous across sites. To model this heterogeneity we generated a simulation of CpG Methylation changes with time and investigated the functional form of the dynamics of methylation with age under the empirical distribution of timescales estimated from the dataset. The resulting dynamics of the average methylation of the system were characterized and were found to closely follow an exponential trajectory. We conclude that DNA methylation can be modeled as a system that starts out of equilibrium at birth and approaches equilibrium with age in an exponential fashion. These insights illustrate the importance of accounting for nonlinear dynamics when utilizing age associated DNA methylation changes for constructing biomarkers of aging. Thus DNA methylation, along with the exponentially increasing risk of mortality with age, further establishes the exponential nature of aging.

A complex network of interactions governs DNA methylation at telomeric regions.

DNA methylation modulates telomere function. In Arabidopsis thaliana, telomeric regions have a bimodal chromatin organization with unmethylated telomeres and methylated subtelomeres. To gain insight into this organization we have generated TAIR10-Tel, a modified version of the Arabidopsis reference genome with additional sequences at most chromosome ends. TAIR10-Tel has allowed us to analyse DNA methylation at nucleotide resolution level in telomeric regions. We have analysed the wild-type strain and mutants that encode inactive versions of all currently known relevant methyltransferases involved in cytosine methylation. These analyses have revealed that subtelomeric DNA methylation extends 1 to 2 kbp from Interstitial Telomeric Sequences (ITSs) that abut or are very near to telomeres. However, DNA methylation drops at the telomeric side of the telomere-subtelomere boundaries and disappears at the inner part of telomeres. We present a comprehensive and integrative model for subtelomeric DNA methylation that should help to decipher the mechanisms that govern the epigenetic regulation of telomeres. This model involves a complex network of interactions between methyltransferases and subtelomeric DNA sequences.

The IntAct database: efficient access to fine-grained molecular interaction data.

The IntAct molecular interaction database (https://www.ebi.ac.uk/intact) is a curated resource of molecular interactions, derived from the scientific literature and from direct data depositions. As of August 2021, IntAct provides more than one million binary interactions, curated by twelve global partners of the International Molecular Exchange consortium, for which the IntAct database provides a shared curation and dissemination platform. The IMEx curation policy has always emphasised a fine-grained data and curation model, aiming to capture the relevant experimental detail essential for the interpretation of the provided molecular interaction data. Here, we present recent curation focus and progress, as well as a completely redeveloped website which presents IntAct data in a much more user-friendly and detailed way.