Single Cell Multi-Omics of an iPSC Model of Human Sinoatrial Node Development Reveals Genetic Determinants of Heart Rate and Arrhythmia Susceptibility.

Cellular heterogeneity within the sinoatrial node (SAN) is functionally important but has been difficult to model in vitro , presenting a major obstacle to studies of heart rate regulation and arrhythmias. Here we describe a scalable method to derive sinoatrial node pacemaker cardiomyocytes (PCs) from human induced pluripotent stem cells that recapitulates differentiation into distinct PC subtypes, including SAN Head, SAN Tail, transitional zone cells, and sinus venosus myocardium. Single cell (sc) RNA-sequencing, sc-ATAC-sequencing, and trajectory analyses were used to define epigenetic and transcriptomic signatures of each cell type, and to identify novel transcriptional pathways important for PC subtype differentiation. Integration of our multi-omics datasets with genome wide association studies uncovered cell type-specific regulatory elements that associated with heart rate regulation and susceptibility to atrial fibrillation. Taken together, these datasets validate a novel, robust, and realistic in vitro platform that will enable deeper mechanistic exploration of human cardiac automaticity and arrhythmia.

Phagocytosis increases an oxidative metabolic and immune suppressive signature in tumor macrophages.

Phagocytosis is a key macrophage function, but how phagocytosis shapes tumor-associated macrophage (TAM) phenotypes and heterogeneity in solid tumors remains unclear. Here, we utilized both syngeneic and novel autochthonous lung tumor models in which neoplastic cells express the fluorophore tdTomato (tdTom) to identify TAMs that have phagocytosed neoplastic cells in vivo. Phagocytic tdTompos TAMs upregulated antigen presentation and anti-inflammatory proteins, but downregulated classic proinflammatory effectors compared to tdTomneg TAMs. Single-cell transcriptomic profiling identified TAM subset-specific and common gene expression changes associated with phagocytosis. We uncover a phagocytic signature that is predominated by oxidative phosphorylation (OXPHOS), ribosomal, and metabolic genes, and this signature correlates with worse clinical outcome in human lung cancer. Expression of OXPHOS proteins, mitochondrial content, and functional utilization of OXPHOS were increased in tdTompos TAMs. tdTompos tumor dendritic cells also display similar metabolic changes. Our identification of phagocytic TAMs as a distinct myeloid cell state links phagocytosis of neoplastic cells in vivo with OXPHOS and tumor-promoting phenotypes.

Blinatumomab-induced T cell activation at single cell transcriptome resolution.

BACKGROUND: Bi-specific T-cell engager (BiTE) antibody is a class of bispecific antibodies designed for cancer immunotherapy. Blinatumomab is the first approved BiTE to treat acute B cell lymphoblastic leukemia (B-ALL). It brings killer T and target B cells into close proximity, activating patient's autologous T cells to kill malignant B cells via mechanisms such as cytolytic immune synapse formation and inflammatory cytokine production. However, the activated T-cell subtypes and the target cell-dependent T cell responses induced by blinatumomab, as well as the mechanisms of resistance to blinatumomab therapy are largely unknown. RESULTS: In this study, we performed single-cell sequencing analysis to identify transcriptional changes in T cells following blinatumomab-induced T cell activation using single cells from both, a human cell line model and a patient-derived model of blinatumomab-mediated cytotoxicity. In total, the transcriptome of 17,920 single T cells from the cell line model and 2271 single T cells from patient samples were analyzed. We found that CD8+ effector memory T cells, CD4+ central memory T cells, naïve T cells, and regulatory T cells were activated after blinatumomab treatment. Here, blinatumomab-induced transcriptional changes reflected the functional immune activity of the blinatumomab-activated T cells, including the upregulation of pathways such as the immune system, glycolysis, IFNA signaling, gap junctions, and IFNG signaling. Co-stimulatory (TNFRSF4 and TNFRSF18) and co-inhibitory (LAG3) receptors were similarly upregulated in blinatumomab-activated T cells, indicating ligand-dependent T cell functions. Particularly, B-ALL cell expression of TNFSF4, which encodes the ligand of T cell co-stimulatory receptor TNFRSF4, was found positively correlated with the response to blinatumomab treatment. Furthermore, recombinant human TNFSF4 protein enhanced the cytotoxic activity of blinatumomab against B-ALL cells. CONCLUSION: These results reveal a target cell-dependent mechanism of T-cell activation by blinatumomab and suggest that TNFSF4 may be responsible for the resistant mechanism and a potential target for combination therapy with blinatumomab, to treat B-ALL or other B-cell malignancies.

Systematic comparison of high-throughput single-cell RNA-seq methods for immune cell profiling.

BACKGROUND: Elucidation of immune populations with single-cell RNA-seq has greatly benefited the field of immunology by deepening the characterization of immune heterogeneity and leading to the discovery of new subtypes. However, single-cell methods inherently suffer from limitations in the recovery of complete transcriptomes due to the prevalence of cellular and transcriptional dropout events. This issue is often compounded by limited sample availability and limited prior knowledge of heterogeneity, which can confound data interpretation. RESULTS: Here, we systematically benchmarked seven high-throughput single-cell RNA-seq methods. We prepared 21 libraries under identical conditions of a defined mixture of two human and two murine lymphocyte cell lines, simulating heterogeneity across immune-cell types and cell sizes. We evaluated methods by their cell recovery rate, library efficiency, sensitivity, and ability to recover expression signatures for each cell type. We observed higher mRNA detection sensitivity with the 10x Genomics 5' v1 and 3' v3 methods. We demonstrate that these methods have fewer dropout events, which facilitates the identification of differentially-expressed genes and improves the concordance of single-cell profiles to immune bulk RNA-seq signatures. CONCLUSION: Overall, our characterization of immune cell mixtures provides useful metrics, which can guide selection of a high-throughput single-cell RNA-seq method for profiling more complex immune-cell heterogeneity usually found in vivo.

Dynamic changes in the regulatory T-cell heterogeneity and function by murine IL-2 mutein.

The therapeutic expansion of Foxp3+ regulatory T cells (Tregs) shows promise for treating autoimmune and inflammatory disorders. Yet, how this treatment affects the heterogeneity and function of Tregs is not clear. Using single-cell RNA-seq analysis, we characterized 31,908 Tregs from the mice treated with a half-life extended mutant form of murine IL-2 (IL-2 mutein, IL-2M) that preferentially expanded Tregs, or mouse IgG Fc as a control. Cell clustering analysis revealed that IL-2M specifically expands multiple sub-states of Tregs with distinct expression profiles. TCR profiling with single-cell analysis uncovered Treg migration across tissues and transcriptional changes between clonally related Tregs after IL-2M treatment. Finally, we identified IL-2M-expanded Tnfrsf9+Il1rl1+ Tregs with superior suppressive function, highlighting the potential of IL-2M to expand highly suppressive Foxp3+ Tregs.

Affinity Maturation of the Anti-Citrullinated Protein Antibody Paratope Drives Epitope Spreading and Polyreactivity in Rheumatoid Arthritis.

OBJECTIVE: Anti-citrullinated protein antibodies (ACPAs) are a hallmark of rheumatoid arthritis (RA). While epitope spreading of the serum ACPA response is believed to contribute to RA pathogenesis, little is understood regarding how this phenomenon occurs. This study was undertaken to analyze the antibody repertoires of individuals with RA to gain insight into the mechanisms leading to epitope spreading of the serum ACPA response in RA. METHODS: Plasmablasts from the blood of 6 RA patients were stained with citrullinated peptide tetramers to identify ACPA-producing B cells by flow cytometry. Plasmablasts were single-cell sorted and sequenced to obtain antibody repertoires. Sixty-nine antibodies were recombinantly expressed, and their anticitrulline reactivities were characterized using a cyclic citrullinated peptide enzyme-linked immuosorbent assay and synovial antigen arrays. Thirty-six mutated antibodies designed either to represent ancestral antibodies or to test paratope residues critical for binding, as determined from molecular modeling studies, were also tested for anticitrulline reactivities. RESULTS: Clonally related monoclonal ACPAs and their shared ancestral antibodies each exhibited differential reactivity against citrullinated antigens. Molecular modeling identified residues within the complementarity-determining region loops and framework regions predicted to be important for citrullinated antigen binding. Affinity maturation resulted in mutations of these key residues, which conferred binding to different citrullinated epitopes and/or increased polyreactivity to citrullinated epitopes. CONCLUSION: These results demonstrate that the different somatic hypermutations accumulated by clonally related B cells during affinity maturation alter the antibody paratope to mediate epitope spreading and polyreactivity of the ACPA response in RA, suggesting that these may be key properties that likely contribute to the pathogenicity of ACPAs.

T Cell-Dependent Affinity Maturation and Innate Immune Pathways Differentially Drive Autoreactive B Cell Responses in Rheumatoid Arthritis.

OBJECTIVE: Rheumatoid arthritis (RA) is characterized by the activation of B cells that produce anti-citrullinated protein antibodies (ACPAs) and rheumatoid factors (RFs), but the mechanisms by which tolerance is broken in these B cells remain incompletely understood. We undertook this study to investigate whether ACPA+ and RF+ B cells break tolerance through distinct molecular mechanisms. METHODS: We developed antigen-tetramers to isolate ACPA+ and RF+ B cells and performed single-cell RNA sequencing on 2,349 B cells from 6 RA patients and 1 healthy donor to analyze their immunoglobulin repertoires and transcriptional programs. Prominent immunoglobulins were expressed as monoclonal antibodies and tested for autoantigen reactivity. RESULTS: ACPA+ and RF+ B cells were enriched in the peripheral blood of RA patients relative to healthy controls. Characterization of patient-derived monoclonal antibodies confirmed ACPA and RF targeting of tetramer-specific B cells at both antigen-inexperienced and affinity-matured B cell stages. ACPA+ B cells used more class-switched isotypes and exhibited more somatic hypermutations relative to RF+ B cells, and these differences were accompanied by down-regulation of CD72 and up-regulation of genes that promote class-switching and T cell-dependent responses. In contrast, RF+ B cells expressed transcriptional programs that stimulate rapid memory reactivation through multiple innate immune pathways. Coexpression analysis revealed that ACPA+ and RF+ B cell-enriched genes belong to distinct transcriptional regulatory networks. CONCLUSION: Our findings suggest that ACPA+ and RF+ B cells are imprinted with distinct transcriptional programs, which suggests that these autoantibodies associated with increased inflammation in RA arise from 2 different molecular mechanisms.

Street-experienced peripheral B cells traffic to the brain.

Deep sequencing of immunoglobulin repertoires in multiple sclerosis patients reveals dysregulation of peripheral B cell immunity (Palanichamy et al. and Stern et al., this issue).

Barcode-enabled sequencing of plasmablast antibody repertoires in rheumatoid arthritis.

OBJECTIVE: A hallmark of rheumatoid arthritis (RA) is the production of autoantibodies, including anti-citrullinated protein antibodies (ACPAs). Nevertheless, the specific targets of these autoantibodies remain incompletely defined. During an immune response, B cells specific for the inciting antigen(s) are activated and differentiate into plasmablasts, which are released into the blood. We undertook this study to sequence the plasmablast antibody repertoire to define the targets of the active immune response in RA. METHODS: We developed a novel DNA barcoding method to sequence the cognate heavy- and light-chain pairs of antibodies expressed by individual blood plasmablasts in RA. The method uses a universal 5' adapter that enables full-length sequencing of the antibodies' variable regions and recombinant expression of the paired antibody chains. The sequence data sets were bioinformatically analyzed to generate phylogenetic trees that identify clonal families of antibodies sharing heavy- and light-chain VJ sequences. Representative antibodies were expressed, and their binding properties were characterized using anti-cyclic citrullinated peptide 2 (anti-CCP-2) enzyme-linked immunosorbent assay (ELISA) and antigen microarrays. RESULTS: We used our sequencing method to generate phylogenetic trees representing the antibody repertoires of peripheral blood plasmablasts from 4 individuals with anti-CCP+ RA, and recombinantly expressed 14 antibodies that were either "singleton" antibodies or representative of clonal antibody families. Anti-CCP-2 ELISA identified 4 ACPAs, and antigen microarray analysis identified ACPAs that differentially targeted epitopes on α-enolase, citrullinated fibrinogen, and citrullinated histone H2B. CONCLUSION: Our data provide evidence that autoantibodies targeting α-enolase, citrullinated fibrinogen, and citrullinated histone H2B are produced by the ongoing activated B cell response in, and thus may contribute to the pathogenesis of, RA.

Identifying functional anti-Staphylococcus aureus antibodies by sequencing antibody repertoires of patient plasmablasts.

Infection by Staphylococcus aureus is on the rise, and there is a need for a better understanding of host immune responses that combat S. aureus. Here we use DNA barcoding to enable deep sequencing of the paired heavy- and light-chain immunoglobulin genes expressed by individual plasmablasts derived from S. aureus-infected humans. Bioinformatic analysis of the antibody repertoires revealed clonal families of heavy-chain sequences and enabled rational selection of antibodies for recombinant expression. Of the ten recombinant antibodies produced, seven bound to S. aureus, of which four promoted opsonophagocytosis of S. aureus. Five of the antibodies bound to known S. aureus cell-surface antigens, including fibronectin-binding protein A. Fibronectin-binding protein A-specific antibodies were isolated from two independent S. aureus-infected patients and mediated neutrophil killing of S. aureus in in vitro assays. Thus, our DNA barcoding approach enabled efficient identification of antibodies involved in protective host antibody responses against S. aureus.

Interspecific germline transmission of cultured primordial germ cells.

In birds, the primordial germ cell (PGC) lineage separates from the soma within 24 h following fertilization. Here we show that the endogenous population of about 200 PGCs from a single chicken embryo can be expanded one million fold in culture. When cultured PGCs are injected into a xenogeneic embryo at an equivalent stage of development, they colonize the testis. At sexual maturity, these donor PGCs undergo spermatogenesis in the xenogeneic host and become functional sperm. Insemination of semen from the xenogeneic host into females from the donor species produces normal offspring from the donor species. In our model system, the donor species is chicken (Gallus domesticus) and the recipient species is guinea fowl (Numida meleagris), a member of a different avian family, suggesting that the mechanisms controlling proliferation of the germline are highly conserved within birds. From a pragmatic perspective, these data are the basis of a novel strategy to produce endangered species of birds using domesticated hosts that are both tractable and fecund.

Early evolution of the LIM homeobox gene family.

BACKGROUND: LIM homeobox (Lhx) transcription factors are unique to the animal lineage and have patterning roles during embryonic development in flies, nematodes and vertebrates, with a conserved role in specifying neuronal identity. Though genes of this family have been reported in a sponge and a cnidarian, the expression patterns and functions of the Lhx family during development in non-bilaterian phyla are not known. RESULTS: We identified Lhx genes in two cnidarians and a placozoan and report the expression of Lhx genes during embryonic development in Nematostella and the demosponge Amphimedon. Members of the six major LIM homeobox subfamilies are represented in the genomes of the starlet sea anemone, Nematostella vectensis, and the placozoan Trichoplax adhaerens. The hydrozoan cnidarian, Hydra magnipapillata, has retained four of the six Lhx subfamilies, but apparently lost two others. Only three subfamilies are represented in the haplosclerid demosponge Amphimedon queenslandica. A tandem cluster of three Lhx genes of different subfamilies and a gene containing two LIM domains in the genome of T. adhaerens (an animal without any neurons) indicates that Lhx subfamilies were generated by tandem duplication. This tandem cluster in Trichoplax is likely a remnant of the original chromosomal context in which Lhx subfamilies first appeared. Three of the six Trichoplax Lhx genes are expressed in animals in laboratory culture, as are all Lhx genes in Hydra. Expression patterns of Nematostella Lhx genes correlate with neural territories in larval and juvenile polyp stages. In the aneural demosponge, A. queenslandica, the three Lhx genes are expressed widely during development, including in cells that are associated with the larval photosensory ring. CONCLUSIONS: The Lhx family expanded and diversified early in animal evolution, with all six subfamilies already diverged prior to the cnidarian-placozoan-bilaterian last common ancestor. In Nematostella, Lhx gene expression is correlated with neural territories in larval and juvenile polyp stages. This pattern is consistent with a possible role in patterning the Nematostella nervous system. We propose a scenario in which Lhx genes play a homologous role in neural patterning across eumetazoans.