Protective immunoglobulin A (IgA) responses to oral antigens are usually orchestrated by gut dendritic cells (DCs). Here, we show that lung CD103(+) and CD24(+)CD11b(+) DCs induced IgA class-switch recombination (CSR) by activating B cells through T cell-dependent or -independent pathways. Compared with lung DCs (LDC), lung CD64(+) macrophages had decreased expression of B cell activation genes and induced significantly less IgA production. Microbial stimuli, acting through Toll-like receptors, induced transforming growth factor-β (TGF-β) production by LDCs and exerted a profound influence on LDC-mediated IgA CSR. After intranasal immunization with inactive cholera toxin (CT), LDCs stimulated retinoic acid-dependent up-regulation of α4β7 and CCR9 gut-homing receptors on local IgA-expressing B cells. Migration of these B cells to the gut resulted in IgA-mediated protection against an oral challenge with active CT. However, in germ-free mice, the levels of LDC-induced, CT-specific IgA in the gut are significantly reduced. Herein, we demonstrate an unexpected role of the microbiota in modulating the protective efficacy of intranasal vaccination through their effect on the IgA class-switching function of LDCs.

Whole genome methylation profiling at a single cytosine resolution is now feasible due to the advent of high-throughput sequencing techniques together with bisulfite treatment of the DNA. To obtain the methylation value of each individual cytosine, the bisulfite-treated sequence reads are first aligned to a reference genome, and then the profiling of the methylation levels is done from the alignments. A huge effort has been made to quickly and correctly align the reads and many different algorithms and programs to do this have been created. However, the second step is just as crucial and non-trivial, but much less attention has been paid to the final inference of the methylation states. Important error sources do exist, such as sequencing errors, bisulfite failure, clonal reads, and single nucleotide variants. We developed MethylExtract, a user friendly tool to: i) generate high quality, whole genome methylation maps and ii) detect sequence variation within the same sample preparation. The program is implemented into a single script and takes into account all major error sources. MethylExtract detects variation (SNVs - Single Nucleotide Variants) in a similar way to VarScan, a very sensitive method extensively used in SNV and genotype calling based on non-bisulfite-treated reads. The usefulness of MethylExtract is shown by means of extensive benchmarking based on artificial bisulfite-treated reads and a comparison to a recently published method, called Bis-SNP. MethylExtract is able to detect SNVs within High-Throughput Sequencing experiments of bisulfite treated DNA at the same time as it generates high quality methylation maps. This simultaneous detection of DNA methylation and sequence variation is crucial for many downstream analyses, for example when deciphering the impact of SNVs on differential methylation. An exclusive feature of MethylExtract, in comparison with existing software, is the possibility to assess the bisulfite failure in a statistical way. The source code, tutorial and artificial bisulfite datasets are available at http://bioinfo2.ugr.es/MethylExtract/ and http://sourceforge.net/projects/methylextract/, and also permanently accessible from 10.5281/zenodo.7144.

Cancer cell dissemination during very early stages of breast cancer proceeds through poorly understood mechanisms. Here we show, in a mouse model of HER2+ breast cancer, that a previously described sub-population of early-evolved cancer cells requires macrophages for early dissemination. Depletion of macrophages specifically during pre-malignant stages reduces early dissemination and also results in reduced metastatic burden at end stages of cancer progression. Mechanistically, we show that, in pre-malignant lesions, CCL2 produced by cancer cells and myeloid cells attracts CD206+/Tie2+ macrophages and induces Wnt-1 upregulation that in turn downregulates E-cadherin junctions in the HER2+ early cancer cells. We also observe macrophage-containing tumor microenvironments of metastasis structures in the pre-malignant lesions that can operate as portals for intravasation. These data support a causal role for macrophages in early dissemination that affects long-term metastasis development much later in cancer progression. A pilot analysis on human specimens revealed intra-epithelial macrophages and loss of E-cadherin junctions in ductal carcinoma in situ, supporting a potential clinical relevance.

Microglia are embryonically seeded macrophages that contribute to brain development, homeostasis, and pathologies. It is thus essential to decipher how microglial properties are temporally regulated by intrinsic and extrinsic factors, such as sexual identity and the microbiome. Here, we found that microglia undergo differentiation phases, discernable by transcriptomic signatures and chromatin accessibility landscapes, which can diverge in adult males and females. Remarkably, the absence of microbiome in germ-free mice had a time and sexually dimorphic impact both prenatally and postnatally: microglia were more profoundly perturbed in male embryos and female adults. Antibiotic treatment of adult mice triggered sexually biased microglial responses revealing both acute and long-term effects of microbiota depletion. Finally, human fetal microglia exhibited significant overlap with the murine transcriptomic signature. Our study shows that microglia respond to environmental challenges in a sex- and time-dependent manner from prenatal stages, with major implications for our understanding of microglial contributions to health and disease.

Initially, children were thought to be spared from disease caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). However, a month into the epidemic, a novel multisystem inflammatory syndrome in children (MIS-C) emerged. Herein, we report on the immune profiles of nine MIS-C cases. All MIS-C patients had evidence of prior SARS-CoV-2 exposure, mounting an antibody response with intact neutralization capability. Cytokine profiling identified elevated signatures of inflammation (IL-18 and IL-6), lymphocytic and myeloid chemotaxis and activation (CCL3, CCL4, and CDCP1), and mucosal immune dysregulation (IL-17A, CCL20, and CCL28). Immunophenotyping of peripheral blood revealed reductions of non-classical monocytes, and subsets of NK and T lymphocytes, suggesting extravasation to affected tissues. Finally, profiling the autoantigen reactivity of MIS-C plasma revealed both known disease-associated autoantibodies (anti-La) and novel candidates that recognize endothelial, gastrointestinal, and immune-cell antigens. All patients were treated with anti-IL-6R antibody and/or IVIG, which led to rapid disease resolution.

Genetic variation, gene expression and DNA methylation influence each other in a complex way. To study the impact of sequence variation and DNA methylation on gene expression, we generated geno5mC, a database that contains statistically significant SNP-CpG associations that are biologically classified either through co-localization with known regulatory regions (promoters and enhancers), or through known correlations with the expression levels of nearby genes. The SNP rs727563 can be used to illustrate the usefulness of this approach. This SNP has been associated with inflammatory bowel disease through GWAS, but it is not located near any gene related to this phenotype. However, geno5mC reveals that rs727563 is associated with the methylation state of several CpGs located in promoter regions of genes reported to be involved in inflammatory processes. This case exemplifies how geno5mC can be used to infer relevant and previously unknown interactions between described disease-associated SNPs and their functional targets.

Primary Sjögren's syndrome (SS) is a systemic autoimmune disease characterized by lymphocytic infiltration and damage of exocrine salivary and lacrimal glands. The etiology of SS is complex with environmental triggers and genetic factors involved. By conducting an integrated multi-omics study, we confirmed a vast coordinated hypomethylation and overexpression effects in IFN-related genes, what is known as the IFN signature. Stratified and conditional analyses suggest a strong interaction between SS-associated HLA genetic variation and the presence of Anti-Ro/SSA autoantibodies in driving the IFN epigenetic signature and determining SS. We report a novel epigenetic signature characterized by increased DNA methylation levels in a large number of genes enriched in pathways such as collagen metabolism and extracellular matrix organization. We identified potential new genetic variants associated with SS that might mediate their risk by altering DNA methylation or gene expression patterns, as well as disease-interacting genetic variants that exhibit regulatory function only in the SS population. Our study sheds new light on the interaction between genetics, autoantibody profiles, DNA methylation and gene expression in SS, and contributes to elucidate the genetic architecture of gene regulation in an autoimmune population.

Given that gastrointestinal (GI) symptoms are a prominent extrapulmonary manifestation of COVID-19, we investigated intestinal infection with SARS-CoV-2, its effect on pathogenesis, and clinical significance.

Huntington's disease (HD) is a progressive neurodegenerative disorder caused by a CAG trinucleotide repeat expansion in the huntingtin (HTT) gene, which encodes a polyglutamine tract in the HTT protein. We found that peroxisome proliferator-activated receptor delta (PPAR-δ) interacts with HTT and that mutant HTT represses PPAR-δ-mediated transactivation. Increased PPAR-δ transactivation ameliorated mitochondrial dysfunction and improved cell survival of neurons from mouse models of HD. Expression of dominant-negative PPAR-δ in the central nervous system of mice was sufficient to induce motor dysfunction, neurodegeneration, mitochondrial abnormalities and transcriptional alterations that recapitulated HD-like phenotypes. Expression of dominant-negative PPAR-δ specifically in the striatum of medium spiny neurons in mice yielded HD-like motor phenotypes, accompanied by striatal neuron loss. In mouse models of HD, pharmacologic activation of PPAR-δ using the agonist KD3010 improved motor function, reduced neurodegeneration and increased survival. PPAR-δ activation also reduced HTT-induced neurotoxicity in vitro and in medium spiny-like neurons generated from stem cells derived from individuals with HD, indicating that PPAR-δ activation may be beneficial in HD and related disorders.

Hypomethylated, CpG-rich DNA segments (CpG islands, CGIs) are epigenome markers involved in key biological processes. Aberrant methylation is implicated in the appearance of several disorders as cancer, immunodeficiency, or centromere instability. Furthermore, methylation differences at promoter regions between human and chimpanzee strongly associate with genes involved in neurological/psychological disorders and cancers. Therefore, the evolutionary comparative analyses of CGIs can provide insights on the functional role of these epigenome markers in both health and disease. Given the lack of specific tools, we developed CpGislandEVO. Briefly, we first compile a database of statistically significant CGIs for the best assembled mammalian genome sequences available to date. Second, by means of a coupled browser front-end, we focus on the CGIs overlapping orthologous genes extracted from OrthoDB, thus ensuring the comparison between CGIs located on truly homologous genome segments. This allows comparing the main compositional features between homologous CGIs. Finally, to facilitate nucleotide comparisons, we lifted genome coordinates between assemblies from different species, which enables the analysis of sequence divergence by direct count of nucleotide substitutions and indels occurring between homologous CGIs. The resulting CpGislandEVO database, linking together CGIs and single-cytosine DNA methylation data from several mammalian species, is freely available at our website.

The zinc finger transcription factor EGR-2 has been shown to play an important role in the induction of T cell anergy and the regulation of peripheral T cell tolerance. In vitro, a prior study has show that T cells deficient in EGR-2 are hyperproliferative to IL-2 and produce elevated levels of the effector cytokine IFN-γ. EGR-2 deficient mice have increased levels of CD44(high) T cells in peripheral lymphoid organs, and with age, develop autoimmune-like features.

Caloric restriction is known to improve inflammatory and autoimmune diseases. However, the mechanisms by which reduced caloric intake modulates inflammation are poorly understood. Here we show that short-term fasting reduced monocyte metabolic and inflammatory activity and drastically reduced the number of circulating monocytes. Regulation of peripheral monocyte numbers was dependent on dietary glucose and protein levels. Specifically, we found that activation of the low-energy sensor 5'-AMP-activated protein kinase (AMPK) in hepatocytes and suppression of systemic CCL2 production by peroxisome proliferator-activator receptor alpha (PPARα) reduced monocyte mobilization from the bone marrow. Importantly, we show that fasting improves chronic inflammatory diseases without compromising monocyte emergency mobilization during acute infectious inflammation and tissue repair. These results reveal that caloric intake and liver energy sensors dictate the blood and tissue immune tone and link dietary habits to inflammatory disease outcome.

Recent technical advances highlight that to understand mammalian development and human disease we need to consider transcriptional and epigenetic cell-to-cell differences within cell populations. This is particularly important in key areas of biomedicine like stem cell differentiation and intratumor heterogeneity. The recently developed nucleosome occupancy and methylome (NOMe) assay facilitates the simultaneous study of DNA methylation and nucleosome positioning on the same DNA strand. NOMe-treated DNA can be sequenced by sanger (NOMe-PCR) or high throughput approaches (NOMe-seq). NOMe-PCR provides information for a single locus at the single molecule while NOMe-seq delivers genome-wide data that is usually interrogated to obtain population-averaged measures. Here, we have developed a bioinformatic tool that allow us to easily obtain locus-specific information at the single molecule using genome-wide NOMe-seq datasets obtained from bulk populations. We have used NOMePlot to study mouse embryonic stem cells and found that polycomb-repressed bivalent gene promoters coexist in two different epigenetic states, as defined by the nucleosome binding pattern detected around their transcriptional start site.

While conventional pathogenic protists have been extensively studied, there is an underappreciated constitutive protist microbiota that is an integral part of the vertebrate microbiome. The impact of these species on the host and their potential contributions to mucosal immune homeostasis remain poorly studied. Here, we show that the protozoan Tritrichomonas musculis activates the host epithelial inflammasome to induce IL-18 release. Epithelial-derived IL-18 promotes dendritic cell-driven Th1 and Th17 immunity and confers dramatic protection from mucosal bacterial infections. Along with its role as a "protistic" antibiotic, colonization with T. musculis exacerbates the development of T-cell-driven colitis and sporadic colorectal tumors. Our findings demonstrate a novel mutualistic host-protozoan interaction that increases mucosal host defenses at the cost of an increased risk of inflammatory disease.

Central nervous system (CNS) macrophages comprise microglia and border-associated macrophages (BAMs) residing in the meninges, the choroid plexus, and the perivascular spaces. Most CNS macrophages emerge during development, with the exception of choroid plexus and dural macrophages, which are replaced by monocytes in adulthood. Whether microglia and BAMs share a developmental program or arise from separate lineages remains unknown. Here, we identified two phenotypically, transcriptionally, and locally distinct brain macrophages throughout development, giving rise to either microglia or BAMs. Two macrophage populations were already present in the yolk sac suggesting an early segregation. Fate-mapping models revealed that BAMs mostly derived from early erythro-myeloid progenitors in the yolk sac. The development of microglia was dependent on TGF-β, whereas the genesis of BAMs occurred independently of this cytokine. Collectively, our data show that developing parenchymal and non-parenchymal brain macrophages are separate entities in terms of ontogeny, gene signature, and requirement for TGF-β.

Multisystem inflammatory syndrome in children (MIS-C) presents with fever, inflammation and multiple organ involvement in individuals under 21 years following severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection. To identify genes, pathways and cell types driving MIS-C, we sequenced the blood transcriptomes of MIS-C cases, pediatric cases of coronavirus disease 2019, and healthy controls. We define a MIS-C transcriptional signature partially shared with the transcriptional response to SARS-CoV-2 infection and with the signature of Kawasaki disease, a clinically similar condition. By projecting the MIS-C signature onto a co-expression network, we identified disease gene modules and found genes downregulated in MIS-C clustered in a module enriched for the transcriptional signatures of exhausted CD8 + T-cells and CD56 dim CD57 + NK cells. Bayesian network analyses revealed nine key regulators of this module, including TBX21 , a central coordinator of exhausted CD8 + T-cell differentiation. Together, these findings suggest dysregulated cytotoxic lymphocyte response to SARS-Cov-2 infection in MIS-C.

There is currently no approved treatment for primary Sjögren's syndrome, a disease that primarily affects adult women. The difficulty in developing effective therapies is -in part- because of the heterogeneity in the clinical manifestation and pathophysiology of the disease. Finding common molecular signatures among patient subgroups could improve our understanding of disease etiology, and facilitate the development of targeted therapeutics. Here, we report, in a cross-sectional cohort, a molecular classification scheme for Sjögren's syndrome patients based on the multi-omic profiling of whole blood samples from a European cohort of over 300 patients, and a similar number of age and gender-matched healthy volunteers. Using transcriptomic, genomic, epigenetic, cytokine expression and flow cytometry data, combined with clinical parameters, we identify four groups of patients with distinct patterns of immune dysregulation. The biomarkers we identify can be used by machine learning classifiers to sort future patients into subgroups, allowing the re-evaluation of response to treatments in clinical trials.

The Cancer Immune Monitoring and Analysis Centers - Cancer Immunologic Data Commons (CIMAC-CIDC) Network is supported by the NCI to identify biomarkers of response to cancer immunotherapies across clinical trials using state-of-the-art assays. A primary platform for CIMAC-CIDC studies is cytometry by time of flight (CyTOF), performed at all CIMAC laboratories. To ensure the ability to generate comparable CyTOF data across labs, a multistep cross-site harmonization effort was undertaken.

Multisystem inflammatory syndrome in children (MIS-C) presents with fever, inflammation and pathology of multiple organs in individuals under 21 years of age in the weeks following severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection. Although an autoimmune pathogenesis has been proposed, the genes, pathways and cell types causal to this new disease remain unknown. Here we perform RNA sequencing of blood from patients with MIS-C and controls to find disease-associated genes clustered in a co-expression module annotated to CD56dimCD57+ natural killer (NK) cells and exhausted CD8+ T cells. A similar transcriptome signature is replicated in an independent cohort of Kawasaki disease (KD), the related condition after which MIS-C was initially named. Probing a probabilistic causal network previously constructed from over 1,000 blood transcriptomes both validates the structure of this module and reveals nine key regulators, including TBX21, a central coordinator of exhausted CD8+ T cell differentiation. Together, this unbiased, transcriptome-wide survey implicates downregulation of NK cells and cytotoxic T cell exhaustion in the pathogenesis of MIS-C.

Plasmacytoid dendritic cells (pDCs) can rapidly produce interferons and other soluble factors in response to extracellular viruses or virus mimics such as CpG-containing DNA. pDCs can also recognize live cells infected with certain RNA viruses, but the relevance and functional consequences of such recognition remain unclear. We studied the response of primary DCs to the prototypical persistent DNA virus, human cytomegalovirus (CMV). Human pDCs produced high amounts of type I interferon (IFN-I) when incubated with live CMV-infected fibroblasts but not with free CMV; the response involved integrin-mediated adhesion, transfer of DNA-containing virions to pDCs, and the recognition of DNA through TLR9. Compared with transient polyfunctional responses to CpG or free influenza virus, pDC response to CMV-infected cells was long-lasting, dominated by the production of IFN-I and IFN-III, and lacked diversification into functionally distinct populations. Similarly, pDC activation by influenza-infected lung epithelial cells was highly efficient, prolonged, and dominated by interferon production. Prolonged pDC activation by CMV-infected cells facilitated the activation of natural killer cells critical for CMV control. Last, patients with CMV viremia harbored phenotypically activated pDCs and increased circulating IFN-I and IFN-III. Thus, recognition of live infected cells is a mechanism of virus detection by pDCs that elicits a unique antiviral immune response.