Translational genetics identifies a phosphorylation switch in CARD9 required for innate inflammatory responses.

Population genetics continues to identify genetic variants associated with diseases of the immune system and offers a unique opportunity to discover mechanisms of immune regulation. Multiple genetic variants linked to severe fungal infections and autoimmunity are associated with caspase recruitment domain-containing protein 9 (CARD9). We leverage the CARD9 R101C missense variant to uncover a biochemical mechanism of CARD9 activation essential for antifungal responses. We demonstrate that R101C disrupts a critical signaling switch whereby phosphorylation of S104 releases CARD9 from an autoinhibited state to promote inflammatory responses in myeloid cells. Furthermore, we show that CARD9 R101C exerts dynamic effects on the skin cellular contexture during fungal infection, corrupting inflammatory signaling and cell-cell communication circuits. Card9 R101C mice fail to control dermatophyte infection in the skin, resulting in high fungal burden, yet show minimal signs of inflammation. Together, we demonstrate how translational genetics reveals molecular and cellular mechanisms of innate immune regulation.

Akkermansia muciniphila phospholipid induces homeostatic immune responses.

Multiple studies have established associations between human gut bacteria and host physiology, but determining the molecular mechanisms underlying these associations has been challenging1-3. Akkermansia muciniphila has been robustly associated with positive systemic effects on host metabolism, favourable outcomes to checkpoint blockade in cancer immunotherapy and homeostatic immunity4-7. Here we report the identification of a lipid from A. muciniphila's cell membrane that recapitulates the immunomodulatory activity of A. muciniphila in cell-based assays8. The isolated immunogen, a diacyl phosphatidylethanolamine with two branched chains (a15:0-i15:0 PE), was characterized through both spectroscopic analysis and chemical synthesis. The immunogenic activity of a15:0-i15:0 PE has a highly restricted structure-activity relationship, and its immune signalling requires an unexpected toll-like receptor TLR2-TLR1 heterodimer9,10. Certain features of the phospholipid's activity are worth noting: it is significantly less potent than known natural and synthetic TLR2 agonists; it preferentially induces some inflammatory cytokines but not others; and, at low doses (1% of EC50) it resets activation thresholds and responses for immune signalling. Identifying both the molecule and an equipotent synthetic analogue, its non-canonical TLR2-TLR1 signalling pathway, its immunomodulatory selectivity and its low-dose immunoregulatory effects provide a molecular mechanism for a model of A. muciniphila's ability to set immunological tone and its varied roles in health and disease.

Functional analyses and single cell immunoprofiling uncover sex-specific differences in SARS-CoV2 immune memory development.

SARS-CoV-2 infection leads to a broad range of outcomes and immune responses, with the development of neutralizing antibodies generally correlated with protection against reinfection. Here, we have characterized both neutralizing activity and T cell responses in a cluster of subjects with mild disease linked to a single spreading event. Surprisingly, we observed sex-specific associations between spike- and particularly nucleoprotein-specific T cell responses and neutralization, with pro-inflammatory cytokines being linked to higher titers only in males. Using single cell immunoprofiling, which provided matched transcriptome and T-cell receptor (TCR) profiles in restimulated CD4 + and CD8 + cells from these subjects, we identified differences in type I IFN signaling that may underlie this difference in antibody generation. Finally, we also identified several TCRs associated with cytokine producing T cells. Altogether, our work maps the breadth of immunological outcomes of SARS-CoV2 infections and highlight the potential role of sex-specific feedback loops during the generation of neutralizing antibodies.

Lyme Disease, Borrelia burgdorferi, and Lipid Immunogens.

The human immune system detects potentially pathogenic microbes with receptors that respond to microbial metabolites. While the overall immune signaling pathway is known in considerable detail, the initial molecular signals, the microbially produced immunogens, for important diseases like Lyme disease (LD) are often not well-defined. The immunogens for LD are produced by the spirochete Borrelia burgdorferi, and a galactoglycerolipid (1) has been identified as a key trigger for the inflammatory immune response that characterizes LD. This report corrects the original structural assignment of 1 to 3, a change of an α-galactopyranose to an α-galactofuranose headgroup. The seemingly small change has important implications for the diagnosis, prevention, and treatment of LD.

B cell genomics behind cross-neutralization of SARS-CoV-2 variants and SARS-CoV.

Monoclonal antibodies (mAbs) are a focus in vaccine and therapeutic design to counteract severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and its variants. Here, we combined B cell sorting with single-cell VDJ and RNA sequencing (RNA-seq) and mAb structures to characterize B cell responses against SARS-CoV-2. We show that the SARS-CoV-2-specific B cell repertoire consists of transcriptionally distinct B cell populations with cells producing potently neutralizing antibodies (nAbs) localized in two clusters that resemble memory and activated B cells. Cryo-electron microscopy structures of selected nAbs from these two clusters complexed with SARS-CoV-2 spike trimers show recognition of various receptor-binding domain (RBD) epitopes. One of these mAbs, BG10-19, locks the spike trimer in a closed conformation to potently neutralize SARS-CoV-2, the recently arising mutants B.1.1.7 and B.1.351, and SARS-CoV and cross-reacts with heterologous RBDs. Together, our results characterize transcriptional differences among SARS-CoV-2-specific B cells and uncover cross-neutralizing Ab targets that will inform immunogen and therapeutic design against coronaviruses.

Cell Type- and Stimulation-Dependent Transcriptional Programs Regulated by Atg16L1 and Its Crohn's Disease Risk Variant T300A.

Genome-wide association studies have identified common genetic variants impacting human diseases; however, there are indications that the functional consequences of genetic polymorphisms can be distinct depending on cell type-specific contexts, which produce divergent phenotypic outcomes. Thus, the functional impact of genetic variation and the underlying mechanisms of disease risk are modified by cell type-specific effects of genotype on pathological phenotypes. In this study, we extend these concepts to interrogate the interdependence of cell type- and stimulation-specific programs influenced by the core autophagy gene Atg16L1 and its T300A coding polymorphism identified by genome-wide association studies as linked with increased risk of Crohn's disease. We applied a stimulation-based perturbational profiling approach to define Atg16L1 T300A phenotypes in dendritic cells and T lymphocytes. Accordingly, we identified stimulus-specific transcriptional signatures revealing T300A-dependent functional phenotypes that mechanistically link inflammatory cytokines, IFN response genes, steroid biosynthesis, and lipid metabolism in dendritic cells and iron homeostasis and lysosomal biogenesis in T lymphocytes. Collectively, these studies highlight the combined effects of Atg16L1 genetic variation and stimulatory context on immune function.

Fundamental limits on dynamic inference from single-cell snapshots.

Single-cell expression profiling reveals the molecular states of individual cells with unprecedented detail. Because these methods destroy cells in the process of analysis, they cannot measure how gene expression changes over time. However, some information on dynamics is present in the data: the continuum of molecular states in the population can reflect the trajectory of a typical cell. Many methods for extracting single-cell dynamics from population data have been proposed. However, all such attempts face a common limitation: for any measured distribution of cell states, there are multiple dynamics that could give rise to it, and by extension, multiple possibilities for underlying mechanisms of gene regulation. Here, we describe the aspects of gene expression dynamics that cannot be inferred from a static snapshot alone and identify assumptions necessary to constrain a unique solution for cell dynamics from static snapshots. We translate these constraints into a practical algorithmic approach, population balance analysis (PBA), which makes use of a method from spectral graph theory to solve a class of high-dimensional differential equations. We use simulations to show the strengths and limitations of PBA, and then apply it to single-cell profiles of hematopoietic progenitor cells (HPCs). Cell state predictions from this analysis agree with HPC fate assays reported in several papers over the past two decades. By highlighting the fundamental limits on dynamic inference faced by any method, our framework provides a rigorous basis for dynamic interpretation of a gene expression continuum and clarifies best experimental designs for trajectory reconstruction from static snapshot measurements.

Population snapshots predict early haematopoietic and erythroid hierarchies.

The formation of red blood cells begins with the differentiation of multipotent haematopoietic progenitors. Reconstructing the steps of this differentiation represents a general challenge in stem-cell biology. Here we used single-cell transcriptomics, fate assays and a theory that allows the prediction of cell fates from population snapshots to demonstrate that mouse haematopoietic progenitors differentiate through a continuous, hierarchical structure into seven blood lineages. We uncovered coupling between the erythroid and the basophil or mast cell fates, a global haematopoietic response to erythroid stress and novel growth factor receptors that regulate erythropoiesis. We defined a flow cytometry sorting strategy to purify early stages of erythroid differentiation, completely isolating classically defined burst-forming and colony-forming progenitors. We also found that the cell cycle is progressively remodelled during erythroid development and during a sharp transcriptional switch that ends the colony-forming progenitor stage and activates terminal differentiation. Our work showcases the utility of linking transcriptomic data to predictive fate models, and provides insights into lineage development in vivo.

High-throughput single-cell fate potential assay of murine hematopoietic progenitors in vitro.

The advent of single-cell transcriptomics has led to the proposal of a number of novel high-resolution models for the hematopoietic system. Testing the predictions generated by such models requires cell fate potential assays of matching, single-cell resolution. Here we detail the development of an in vitro high-throughput single-cell culture assay using flow cytometrically sorted single murine bone marrow progenitors, which measures their differentiation into any of five myeloid lineages. We identify critical parameters for single-cell culture outcome, including the choice of sorter nozzle size and pressure, culture media, and the coating of culture dishes with extracellular matrix proteins. Further, we find that accurate assay readout requires the titration of antibodies specifically for their use under low-cell-number conditions. Our approach may be used as a template for the development of single-cell fate potential assays for a variety of blood cell progenitors.

Setd1a and NURF mediate chromatin dynamics and gene regulation during erythroid lineage commitment and differentiation.

The modulation of chromatin structure is a key step in transcription regulation in mammalian cells and eventually determines lineage commitment and differentiation. USF1/2, Setd1a and NURF complexes interact to regulate chromatin architecture in erythropoiesis, but the mechanistic basis for this regulation is hitherto unknown. Here we showed that Setd1a and NURF complexes bind to promoters to control chromatin structural alterations and gene activation in a cell context dependent manner. In human primary erythroid cells USF1/2, H3K4me3 and the NURF complex were significantly co-enriched at transcription start sites of erythroid genes, and their binding was associated with promoter/enhancer accessibility that resulted from nucleosome repositioning. Mice deficient for Setd1a, an H3K4 trimethylase, in the erythroid compartment exhibited reduced Ter119/CD71 positive erythroblasts, peripheral blood RBCs and hemoglobin levels. Loss of Setd1a led to a reduction of promoter-associated H3K4 methylation, inhibition of gene transcription and blockade of erythroid differentiation. This was associated with alterations in NURF complex occupancy at erythroid gene promoters and reduced chromatin accessibility. Setd1a deficiency caused decreased associations between enhancer and promoter looped interactions as well as reduced expression of erythroid genes such as the adult ß-globin gene. These data indicate that Setd1a and NURF complexes are specifically targeted to and coordinately regulate erythroid promoter chromatin dynamics during erythroid lineage differentiation.

Setd1a regulates progenitor B-cell-to-precursor B-cell development through histone H3 lysine 4 trimethylation and Ig heavy-chain rearrangement.

SETD1A is a member of trithorax-related histone methyltransferases that methylate lysine 4 at histone H3 (H3K4). We showed previously that Setd1a is required for mesoderm specification and hematopoietic lineage differentiation in vitro. However, it remains unknown whether or not Setd1a controls specific hematopoietic lineage commitment and differentiation during animal development. Here, we reported that homozygous Setd1a knockout (KO) mice are embryonic lethal. Loss of the Setd1a gene in the hematopoietic compartment resulted in a blockage of the progenitor B-cell-to-precursor B-cell development in bone marrow (BM) and B-cell maturation in spleen. The Setd1a-cKO (conditional knockout) mice exhibited an enlarged spleen with disrupted spleen architecture and leukocytopenia. Mechanistically, Setd1a deficiency in BM reduced the levels of H3K4me3 at critical B-cell gene loci, including Pax5 and Rag1/2, which are critical for the IgH (Ig heavy-chain) locus contractions and rearrangement. Subsequently, the differential long-range looped interactions of the enhancer Eµ with proximal 5' DH region and 3' regulatory regions as well as with Pax5-activated intergenic repeat elements and 5' distal VH genes were compromised by the Setd1a-cKO. Together, our findings revealed a critical role of Setd1a and its mediated epigenetic modifications in regulating the IgH rearrangement and B-cell development.

The H3K27 demethylase JMJD3 is required for maintenance of the embryonic respiratory neuronal network, neonatal breathing, and survival.

JMJD3 (KDM6B) antagonizes Polycomb silencing by demethylating lysine 27 on histone H3. The interplay of methyltransferases and demethylases at this residue is thought to underlie critical cell fate transitions, and the dynamics of H3K27me3 during neurogenesis posited for JMJD3 a critical role in the acquisition of neural fate. Despite evidence of its involvement in early neural commitment, however, its role in the emergence and maturation of the mammalian CNS remains unknown. Here, we inactivated Jmjd3 in the mouse and found that its loss causes perinatal lethality with the complete and selective disruption of the pre-Bötzinger complex (PBC), the pacemaker of the respiratory rhythm generator. Through genetic and electrophysiological approaches, we show that the enzymatic activity of JMJD3 is selectively required for the maintenance of the PBC and controls critical regulators of PBC activity, uncovering an unanticipated role of this enzyme in the late structuring and function of neuronal networks.

The lymphoma-associated NPM-ALK oncogene elicits a p16INK4a/pRb-dependent tumor-suppressive pathway.

Oncogene-induced senescence (OIS) is a barrier for tumor development. Oncogene-dependent DNA damage and activation of the ARF/p53 pathway play a central role in OIS and, accordingly, ARF and p53 are frequently mutated in human cancer. A number of leukemia/lymphoma-initiating oncogenes, however, inhibit ARF/p53 and only infrequently select for ARF or p53 mutations, suggesting the involvement of other tumor-suppressive pathways. We report that NPM-ALK, the initiating oncogene of anaplastic large cell lymphomas (ALCLs), induces DNA damage and irreversibly arrests the cell cycle of primary fibroblasts and hematopoietic progenitors. This effect is associated with inhibition of p53 and is caused by activation of the p16INK4a/pRb tumor-suppressive pathway. Analysis of NPM-ALK lymphomagenesis in transgenic mice showed p16INK4a-dependent accumulation of senescent cells in premalignant lesions and decreased tumor latency in the absence of p16INK4a. Accordingly, human ALCLs showed no expression of either p16INK4a or pRb. Up-regulation of the histone-demethylase Jmjd3 and de-methylation at the p16INK4a promoter contributed to the effect of NPM-ALK on p16INK4a, which was transcriptionally regulated. These data demonstrate that p16INK4a/pRb may function as an alternative pathway of oncogene-induced senescence, and suggest that the reactivation of p16INK4a expression might be a novel strategy to restore the senescence program in some tumors.

A large fraction of extragenic RNA pol II transcription sites overlap enhancers.

Mammalian genomes are pervasively transcribed outside mapped protein-coding genes. One class of extragenic transcription products is represented by long non-coding RNAs (lncRNAs), some of which result from Pol_II transcription of bona-fide RNA genes. Whether all lncRNAs described insofar are products of RNA genes, however, is still unclear. Here we have characterized transcription sites located outside protein-coding genes in a highly regulated response, macrophage activation by endotoxin. Using chromatin signatures, we could unambiguously classify extragenic Pol_II binding sites as belonging to either canonical RNA genes or transcribed enhancers. Unexpectedly, 70% of extragenic Pol_II peaks were associated with genomic regions with a canonical chromatin signature of enhancers. Enhancer-associated extragenic transcription was frequently adjacent to inducible inflammatory genes, was regulated in response to endotoxin stimulation, and generated very low abundance transcripts. Moreover, transcribed enhancers were under purifying selection and contained binding sites for inflammatory transcription factors, thus suggesting their functionality. These data demonstrate that a large fraction of extragenic Pol_II transcription sites can be ascribed to cis-regulatory genomic regions. Discrimination between lncRNAs generated by canonical RNA genes and products of transcribed enhancers will provide a framework for experimental approaches to lncRNAs and help complete the annotation of mammalian genomes.

Genotyping results of Iranian PCG families suggests one or more PCG locus other than GCL3A, GCL3B, and GCL3C exist.

PURPOSE: To assess whether loci other than GLC3A, GLC3B, and GLC3C are linked to primary congenital glaucoma (PCG). METHODS: The gene CYP1B1 at GLC3A was screened in 19 Iranian PCG probands who had been recruited mostly from among individuals of Turkish ethnicity and individuals from central and eastern Iran. The gene MYOC was screened in patients from this cohort who lacked CYP1B1 mutations and in ten patients previously shown not to carry CYP1B1 mutations. Family members of 19 probands without mutations in either of these genes were recruited for assessment of linkage to GLC3B and GLC3C by genotyping microsatellite markers. In total, 127 individuals, including 35 affected with PCG, were genotyped. RESULTS: Eleven (57.9%) of the newly recruited PCG patients did not carry disease-associated mutations in CYP1B1. Disease-associated MYOC mutations were not observed in any of the patients screened. Inheritance of PCG in all the families was consistent with an autosomal recessive pattern. Linkage to GLC3B and GLC3C was ruled out in nine of the families on the basis of autozygosity mapping and haplotype analysis. CONCLUSIONS: Observation of the absence of linkage to GLC3B and GLC3C in at least nine families without CYP1B1 mutations suggests that at least one PCG-causing locus other than GLC3A, GLC3B, and GLC3C may exist. The disease-causing gene or genes in the novel locus or loci may account for PCG in a notable fraction of Iranian patients.

Jmjd3 contributes to the control of gene expression in LPS-activated macrophages.

Jmjd3, a JmjC family histone demethylase, is induced by the transcription factor NF-kB in response to microbial stimuli. Jmjd3 erases H3K27me3, a histone mark associated with transcriptional repression and involved in lineage determination. However, the specific contribution of Jmjd3 induction and H3K27me3 demethylation to inflammatory gene expression remains unknown. Using chromatin immunoprecipitation-sequencing we found that Jmjd3 is preferentially recruited to transcription start sites characterized by high levels of H3K4me3, a marker of gene activity, and RNA polymerase II (Pol_II). Moreover, 70% of lipopolysaccharide (LPS)-inducible genes were found to be Jmjd3 targets. Although most Jmjd3 target genes were unaffected by its deletion, a few hundred genes, including inducible inflammatory genes, showed moderately impaired Pol_II recruitment and transcription. Importantly, most Jmjd3 target genes were not associated with detectable levels of H3K27me3, and transcriptional effects of Jmjd3 absence in the window of time analysed were uncoupled from measurable effects on this histone mark. These data show that Jmjd3 fine-tunes the transcriptional output of LPS-activated macrophages in an H3K27 demethylation-independent manner.

Sex Bias in Primary Congenital Glaucoma Patients with and without CYP1B1 Mutations.

PURPOSE: To investigate variations in sex ratio among Iranian primary congenital glaucoma (PCG) patients with and without mutations in the CYP1B1 gene and to evaluate possible clinical variations associated with sex in these two groups. METHODS: Phenotypical data on 104 unrelated Iranian PCG patients who had previously been screened for CYP1B1 mutations were analyzed. Emphasis was placed on analysis of sex ratios among patients with and without CYP1B1 mutations. In addition to sex, familial and sporadic incidence and clinical features including age at onset, bilateral/unilateral involvement, corneal diameter, intraocular pressure, and cup-disc ratios were compared between these two groups. Information on phenotypical parameters was available for most but not all patients. RESULTS: Among the 93 PCG patients whose sex was recorded, 57 were male (61.3%) and 36 were female (38.7%) (P=0.03). Patients with CYP1B1 mutations included 37 male (66.1%) and 29 female (43.9%) subjects (P=0.30), while patients without the mutation included 20 (74.1%) male and 7 (25.9%) female individuals (P=0.013). Our data did not provide conclusive evidence on difference in severity of the disease between those with and without CYP1B1 mutations, nor between the two sexes. CONCLUSION: Consistent with data on PCG patients from other populations, the overall incidence of PCG in Iran seems to be higher among male subjects. The difference in incidence between the two sexes was not significant among patients whose disease was due to mutations in CYP1B1. The overall higher incidence of PCG among male subjects seems to be attributable to a higher incidence in male patients not harboring CYP1B1 mutations, suggesting that other genes or factors may be involved in manifestation of PCG phenotypes in a sex dependent manner.

Contributions of MYOC and CYP1B1 mutations to JOAG.

PURPOSE: To investigate the role of MYOC and CYP1B1 in Iranian juvenile open angle glaucoma (JOAG) patients. METHODS: Twenty-three JOAG probands, their available affected and unaffected family members, and 100 ethnically matched control individuals without history of ocular disease were recruited. Clinical examinations of the probands included slit lamp biomicroscopy, intraocular pressure (IOP) measurement, gonioscopic evaluation, fundus examination, and perimetry measurement. Familial cases were classified according to the mode of inheritance. Exons of MYOC and CYP1B1 were sequenced, and novel variations assessed in the control individuals. Potential disease-associated variations were tested for segregation with disease status in available family members. RESULTS: The mode of inheritance of the disease in the families of four probands (17.4%) appeared to be autosomal dominant and in at least eight (34.8%) to be autosomal recessive. Four patients carried MYOC mutations, and an equal number carried CYP1B1 mutations. The MYOC mutations were heterozygous; two of them (p.C8X and p.L334P) are novel, and one codes for the shortest truncated protein so far reported. Autosomal recessive inheritance was consistent with inheritance observed in families of patients carrying CYP1B1 mutations. All these patients carried homozygous mutations. CONCLUSIONS: MYOC and CYP1B1 contributed equally to the disease status of the Iranian JOAG patients studied. The contribution of the two genes appeared to be independent in that no patient carried mutations in both genes. The fraction of Iranian patients carrying MYOC mutations was comparable to previously reported populations.

CYP1B1 mutation profile of Iranian primary congenital glaucoma patients and associated haplotypes.

The mutation spectrum of CYP1B1 among 104 primary congenital glaucoma patients of the genetically heterogeneous Iranian population was investigated by sequencing. We also determined intragenic single nucleotide polymorphism (SNP) haplotypes associated with the mutations and compared these with haplotypes of other populations. Finally, the frequency distribution of the haplotypes was compared among primary congenital glaucoma patients with and without CYP1B1 mutations and normal controls. Genotype classification of six high-frequency SNPs was performed using the PHASE 2.0 software. CYP1B1 mutations in the Iranian patients were very heterogeneous. Nineteen nonconservative mutations associated with disease, and 10 variations not associated with disease were identified. Ten mutations and three variations not associated with disease were novel. The 13 novel variations make a notable contribution to the approximately 70 known variations in the gene. CYP1B1 mutations were identified in 70% of the patients. The four most common mutations were G61E, R368H, R390H, and R469W, which together constituted 76.2% of the CYP1B1 mutated alleles found. Six unique core SNP haplotypes were identified, four of which were common to the patients with and without CYP1B1 mutations and controls studied. Three SNP blocks determined the haplotypes. Comparison of haplotypes with those of other populations suggests a common origin for many of the mutations.