Deleting interleukin-10 from myeloid cells exacerbates atherosclerosis in Apoe-/- mice.

Atherosclerosis is initiated by subendothelial retention of lipoproteins and cholesterol, which triggers a non-resolving inflammatory process that over time leads to plaque progression in the artery wall. Myeloid cells and in particular macrophages are the primary drivers of the inflammatory response and plaque formation. Several immune cells including macrophages, T cells and B cells secrete the anti-inflammatory cytokine IL-10, known to be essential for the atherosclerosis protection. The cellular source of IL-10 in natural atherosclerosis progression is unknown. This study aimed to determine the main IL10-producing cell type in atherosclerosis. To do so, we crossed VertX mice, in which IRES-green fluorescent protein (eGFP) was placed downstream of exon 5 of the Il10 gene, with atherosclerosis-prone Apoe-/- mice. We found that myeloid cells express high levels of IL-10 in VertX Apoe-/- mice in both chow and western-diet fed mice. By single cell RNA sequencing and flow cytometry analysis, we identified resident and inflammatory macrophages in atherosclerotic plaques as the main IL-10 producers. To address whether IL-10 secreted by myeloid cells is essential for the protection, we utilized LyzMCre+Il10fl/fl mice crossed into the Apoe-/- background and confirmed that macrophages were unable to secrete IL-10. Chow and western diet-fed LyzMCre+Il10fl/fl Apoe-/- mice developed significantly larger atherosclerotic plaques as measured by en face morphometry than LyzMCre-Il10 fl/flApoe-/-. Flow cytometry and cytokine measurements suggest that the depletion of IL-10 in myeloid cells increases Th17 cells with elevated CCL2, and TNFα in blood plasma. We conclude that macrophage-derived IL-10 is critical for limiting atherosclerosis in mice.

Titration of 124 antibodies using CITE-Seq on human PBMCs.

Single-cell RNA-sequencing (scRNA-Seq) is widely used to characterize immune cell populations. However, mRNA levels correlate poorly with expression of surface proteins, which are well established to define immune cell types. CITE-Seq (cellular indexing of transcriptomes and epitopes by sequencing) utilizes oligonucleotide-tagged antibodies to simultaneously analyze surface phenotypes and transcriptomes. Considering the high costs of adding surface phenotyping to scRNA-Seq, we aimed to determine which of 188 tested CITE-Seq antibodies can detect their antigens on human peripheral blood mononuclear cells (PBMCs), a commonly interrogated cell population in immunology, and find the optimal concentration for staining. The recommended concentration was optimal for 76 antibodies, whereas staining quality of 7 antibodies improved when the concentration was doubled. 33 and 8 antibodies still worked well when the concentration was reduced to 1/5 or 1/25, respectively. 64 antigens were not detected at any antibody concentration. Optimizing the antibody panel by removing antibodies not able to detect their target antigens and adjusting concentrations of the remaining antibodies will improve the analysis and may reduce costs. In conclusion, our data are a resource for building an informative and cost-effective panel of CITE-Seq antibodies and use them at their optimal concentrations in future CITE-seq experiments on human PBMCs.

Sex Differences in Coronary Artery Disease and Diabetes Revealed by scRNA-Seq and CITE-Seq of Human CD4+ T Cells.

Despite the decades-old knowledge that males and people with diabetes mellitus (DM) are at increased risk for coronary artery disease (CAD), the reasons for this association are only partially understood. Among the immune cells involved, recent evidence supports a critical role of T cells as drivers and modifiers of CAD. CD4+ T cells are commonly found in atherosclerotic plaques. We aimed to understand the relationship of CAD with sex and DM by single-cell RNA (scRNA-Seq) and antibody sequencing (CITE-Seq) of CD4+ T cells. Peripheral blood mononuclear cells (PBMCs) of 61 men and women who underwent cardiac catheterization were interrogated by scRNA-Seq combined with 49 surface markers (CITE-Seq). CAD severity was quantified using Gensini scores, with scores above 30 considered CAD+ and below 6 considered CAD-. Four pairs of groups were matched for clinical and demographic parameters. To test how sex and DM changed cell proportions and gene expression, we compared matched groups of men and women, as well as diabetic and non-diabetic subjects. We analyzed 41,782 single CD4+ T cell transcriptomes for sex differences in 16 women and 45 men with and without coronary artery disease and with and without DM. We identified 16 clusters in CD4+ T cells. The proportion of cells in CD4+ effector memory cluster 8 (CD4T8, CCR2+ Em) was significantly decreased in CAD+, especially among DM+ participants. This same cluster, CD4T8, was significantly decreased in female participants, along with two other CD4+ T cell clusters. In CD4+ T cells, 31 genes showed significant and coordinated upregulation in both CAD and DM. The DM gene signature was partially additive to the CAD gene signature. We conclude that (1) CAD and DM are clearly reflected in PBMC transcriptomes, and (2) significant differences exist between women and men and (3) between subjects with DM and non-DM.

Combined protein and transcript single-cell RNA sequencing in human peripheral blood mononuclear cells.

BACKGROUND: Cryopreserved peripheral blood mononuclear cells (PBMCs) are frequently collected and provide disease- and treatment-relevant data in clinical studies. Here, we developed combined protein (40 antibodies) and transcript single-cell (sc)RNA sequencing (scRNA-seq) in PBMCs. RESULTS: Among 31 participants in the Women's Interagency HIV Study (WIHS), we sequenced 41,611 cells. Using Boolean gating followed by Seurat UMAPs (tool for visualizing high-dimensional data) and Louvain clustering, we identified 50 subsets among CD4+ T, CD8+ T, B, NK cells, and monocytes. This resolution was superior to flow cytometry, mass cytometry, or scRNA-seq without antibodies. Combined protein and transcript scRNA-seq allowed for the assessment of disease-related changes in transcriptomes and cell type proportions. As a proof-of-concept, we showed such differences between healthy and matched individuals living with HIV with and without cardiovascular disease. CONCLUSIONS: In conclusion, combined protein and transcript scRNA sequencing is a suitable and powerful method for clinical investigations using PBMCs.

Immunodominant MHC-II (Major Histocompatibility Complex II) Restricted Epitopes in Human Apolipoprotein B.

BACKGROUND: CD (cluster of differentiation) 4+ T-cell responses to APOB (apolipoprotein B) are well characterized in atherosclerotic mice and detectable in humans. CD4+ T cells recognize antigenic peptides displayed on highly polymorphic HLA (human leukocyte antigen)-II. Immunogenicity of individual APOB peptides is largely unknown in humans. Only 1 HLA-II-restricted epitope was validated using the DRB1*07:01-APOB3036-3050 tetramer. We hypothesized that human APOB may contain discrete immunodominant CD4+ T-cell epitopes that trigger atherosclerosis-related autoimmune responses in donors with diverse HLA alleles. METHODS: We selected 20 APOB-derived peptides (APOB20) from an in silico screen and experimentally validated binding to the most commonly occurring human HLA-II alleles. We optimized a restimulation-based workflow to evaluate antigenicity of multiple candidate peptides in HLA-typed donors. This included activation-induced marker assay, intracellular cytokine staining, IFNγ (interferon gamma) enzyme-linked immunospot and cytometric bead array. High-throughput sequencing revealed TCR (T-cell receptor) clonalities of APOB-reactive CD4+ T cells. RESULTS: Using stringent positive, negative, and crossover stimulation controls, we confirmed specificity of expansion-based protocols to detect CD4+ T cytokine responses to the APOB20 pool. Ex vivo assessment of AIM+CD4+ T cells revealed a statistically significant autoimmune response to APOB20 but not to a ubiquitously expressed negative control protein, actin. Resolution of CD4+ T responses to the level of individual peptides using IFNγ enzyme-linked immunospot led to the discovery of 6 immunodominant epitopes (APOB6) that triggered robust CD4+ T activation in most donors. APOB6-specific responding CD4+ T cells were enriched in unique expanded TCR clonotypes and preferentially expressed memory markers. Cytometric bead array analysis detected APOB6-induced secretion of both proinflammatory and regulatory cytokines. In clinical samples from patients with angiographically verified coronary artery disease, APOB6 stimulation induced higher activation and memory phenotypes and augmented secretion of proinflammatory cytokines TNF (tumor necrosis factor) and IFNγ, compared with patients with low coronary artery disease. CONCLUSIONS: Using 3 cohorts, each with ≈20 donors, we discovered and validated 6 immunodominant, HLA-II-restricted APOB epitopes. The immune response to these APOB epitopes correlated with coronary artery disease severity.

Single-Cell Antibody Sequencing in Atherosclerosis Research.

The transcriptomic information obtained by single cell RNA sequencing (scRNA-seq) can be supplemented by information on the cell surface phenotype by using oligonucleotide-tagged monoclonal antibodies (scAb-Seq). This is of particular importance in immune cells, where the correlation between mRNA and cell surface expression is very weak. scAb-Seq is facilitated by the availability of commercial antibodies and antibody mixes. Now panels of up to 200 antibodies are available for human and mouse cells. Proteins are detected by antibodies conjugated to a tripartite DNA sequence that contains a primer for amplification and sequencing, a unique oligonucleotide that acts as an antibody barcode and a poly(dA) sequence, simultaneously detecting extension of antibody-specific DNA sequences and cDNAs in the same poly(dT)-primed reaction. For each cell, surface protein expression is captured and sequenced along with the cell's transcriptome. Here, we list the steps needed to produce antibody sequencing data from tissue or blood cells.

Thymus-Derived CD4+CD8+ Cells Reside in Mediastinal Adipose Tissue and the Aortic Arch.

Double-positive CD4+CD8αß+ (DP) cells are thought to reside as T cell progenitors exclusively within the thymus. We recently discovered an unexpected CD4+ and CD8αß+ immune cell population in healthy and atherosclerotic mice by single-cell RNA sequencing. Transcriptomically, these cells resembled thymic DPs. Flow cytometry and three-dimensional whole-mount imaging confirmed DPs in thymus, mediastinal adipose tissue, and aortic adventitia, but nowhere else. Deep transcriptional profiling revealed differences between DP cells isolated from the three locations. All DPs were dependent on RAG2 expression and the presence of the thymus. Mediastinal adipose tissue DPs resided in close vicinity to invariant NKT cells, which they could activate in vitro. Thymus transplantation failed to reconstitute extrathymic DPs, and frequencies of extrathymic DPs were unaltered by pharmacologic inhibition of S1P1, suggesting that their migration may be locally confined. Our results define two new, transcriptionally distinct subsets of extrathymic DPs that may play a role in aortic vascular homeostasis.

Longitudinal clonal tracking in humanized mice reveals sustained polyclonal repopulation of gene-modified human-HSPC despite vector integration bias.

BACKGROUND: Current understanding of hematopoiesis is largely derived from mouse models that are physiologically distant from humans. Humanized mice provide the most physiologically relevant small animal model to study human diseases, most notably preclinical gene therapy studies. However, the clonal repopulation dynamics of human hematopoietic stem and progenitor cells (HSPC) in these animal models is only partially understood. Using a new clonal tracking methodology designed for small sample volumes, we aim to reveal the underlying clonal dynamics of human cell repopulation in a mouse environment. METHODS: Humanized bone marrow-liver-thymus (hu-BLT) mice were generated by transplanting lentiviral vector-transduced human fetal liver HSPC (FL-HSPC) in NOD.Cg-PrkdcscidIl2rgtm1Wjl/SzJ (NSG) mice implanted with a piece of human fetal thymus. We developed a methodology to track vector integration sites (VIS) in a mere 25 µl of mouse blood for longitudinal and quantitative clonal analysis of human HSPC repopulation in mouse environment. We explored transcriptional and epigenetic features of human HSPC for possible VIS bias. RESULTS: A total of 897 HSPC clones were longitudinally tracked in hu-BLT mice-providing a first-ever demonstration of clonal dynamics and coordinated expansion of therapeutic and control vector-modified human cell populations simultaneously repopulating in the same humanized mice. The polyclonal repopulation stabilized at 19 weeks post-transplant and the contribution of the largest clone doubled within 4 weeks. Moreover, 550 (~ 60%) clones persisted over 6 weeks and were highly shared between different organs. The normal clonal profiles confirmed the safety of our gene therapy vectors. Multi-omics analysis of human FL-HSPC revealed that 54% of vector integrations in repopulating clones occurred within ± 1 kb of H3K36me3-enriched regions. CONCLUSIONS: Human repopulation in mice is polyclonal and stabilizes more rapidly than that previously observed in humans. VIS preference for H3K36me3 has no apparent negative effects on HSPC repopulation. Our study provides a methodology to longitudinally track clonal repopulation in small animal models extensively used for stem cell and gene therapy research and with lentiviral vectors designed for clinical applications. Results of this study provide a framework for understanding the clonal behavior of human HPSC repopulating in a mouse environment, critical for translating results from humanized mice models to the human settings.

Genomic diversity and genome-wide association analysis related to yield and fatty acid composition of wild American oil palm.

Existing Elaeis guineensis cultivars lack sufficient genetic diversity due to extensive breeding. Harnessing variation in wild crop relatives is necessary to expand the breadth of agronomically valuable traits. Using RAD sequencing, we examine the natural diversity of wild American oil palm populations (Elaeis oleifera), a sister species of the cultivated Elaeis guineensis oil palm. We genotyped 192 wild E. oleifera palms collected from seven Latin American countries along with four cultivated E. guineensis palms. Honduras, Costa Rica, Panama and Colombia palms are panmictic and genetically similar. Genomic patterns of diversity suggest that these populations likely originated from the Amazon Basin. Despite evidence of a genetic bottleneck and high inbreeding observed in these populations, there is considerable genetic and phenotypic variation for agronomically valuable traits. Genome-wide association revealed several candidate genes associated with fatty acid composition along with vegetative and yield-related traits. These observations provide valuable insight into the geographic distribution of diversity, phenotypic variation and its genetic architecture that will guide choices of wild genotypes for crop improvement.

Global-level population genomics reveals differential effects of geography and phylogeny on horizontal gene transfer in soil bacteria.

Although microorganisms are known to dominate Earth's biospheres and drive biogeochemical cycling, little is known about the geographic distributions of microbial populations or the environmental factors that pattern those distributions. We used a global-level hierarchical sampling scheme to comprehensively characterize the evolutionary relationships and distributional limitations of the nitrogen-fixing bacterial symbionts of the crop chickpea, generating 1,027 draft whole-genome sequences at the level of bacterial populations, including 14 high-quality PacBio genomes from a phylogenetically representative subset. We find that diverse Mesorhizobium taxa perform symbiosis with chickpea and have largely overlapping global distributions. However, sampled locations cluster based on the phylogenetic diversity of Mesorhizobium populations, and diversity clusters correspond to edaphic and environmental factors, primarily soil type and latitude. Despite long-standing evolutionary divergence and geographic isolation, the diverse taxa observed to nodulate chickpea share a set of integrative conjugative elements (ICEs) that encode the major functions of the symbiosis. This symbiosis ICE takes 2 forms in the bacterial chromosome-tripartite and monopartite-with tripartite ICEs confined to a broadly distributed superspecies clade. The pairwise evolutionary relatedness of these elements is controlled as much by geographic distance as by the evolutionary relatedness of the background genome. In contrast, diversity in the broader gene content of Mesorhizobium genomes follows a tight linear relationship with core genome phylogenetic distance, with little detectable effect of geography. These results illustrate how geography and demography can operate differentially on the evolution of bacterial genomes and offer useful insights for the development of improved technologies for sustainable agriculture.

Between-species differences in gene copy number are enriched among functions critical for adaptive evolution in Arabidopsis halleri.

BACKGROUND: Gene copy number divergence between species is a form of genetic polymorphism that contributes significantly to both genome size and phenotypic variation. In plants, copy number expansions of single genes were implicated in cultivar- or species-specific tolerance of high levels of soil boron, aluminium or calamine-type heavy metals, respectively. Arabidopsis halleri is a zinc- and cadmium-hyperaccumulating extremophile species capable of growing on heavy-metal contaminated, toxic soils. In contrast, its non-accumulating sister species A. lyrata and the closely related reference model species A. thaliana exhibit merely basal metal tolerance. RESULTS: For a genome-wide assessment of the role of copy number divergence (CND) in lineage-specific environmental adaptation, we conducted cross-species array comparative genome hybridizations of three plant species and developed a global signal scaling procedure to adjust for sequence divergence. In A. halleri, transition metal homeostasis functions are enriched twofold among the genes detected as copy number expanded. Moreover, biotic stress functions including mostly disease Resistance (R) gene-related genes are enriched twofold among genes detected as copy number reduced, when compared to the abundance of these functions among all genes. CONCLUSIONS: Our results provide genome-wide support for a link between evolutionary adaptation and CND in A. halleri as shown previously for Heavy metal ATPase4. Moreover our results support the hypothesis that elemental defences, which result from the hyperaccumulation of toxic metals, allow the reduction of classical defences against biotic stress as a trade-off.