Cellular and transcriptional diversity over the course of human lactation.

Human breast milk (hBM) is a dynamic fluid that contains millions of cells, but their identities and phenotypic properties are poorly understood. We generated and analyzed single-cell RNA-sequencing (scRNA-seq) data to characterize the transcriptomes of cells from hBM across lactational time from 3 to 632 d postpartum in 15 donors. We found that the majority of cells in hBM are lactocytes, a specialized epithelial subset, and that cell-type frequencies shift over the course of lactation, yielding greater epithelial diversity at later points. Analysis of lactocytes reveals a continuum of cell states characterized by transcriptional changes in hormone-, growth factor-, and milk production-related pathways. Generalized additive models suggest that one subcluster, LC1 epithelial cells, increases as a function of time postpartum, daycare attendance, and the use of hormonal birth control. We identify several subclusters of macrophages in hBM that are enriched for tolerogenic functions, possibly playing a role in protecting the mammary gland during lactation. Our description of the cellular components of breast milk, their association with maternal­infant dyad metadata, and our quantification of alterations at the gene and pathway levels provide a detailed longitudinal picture of hBM cells across lactational time. This work paves the way for future investigations of how a potential division of cellular labor and differential hormone regulation might be leveraged therapeutically to support healthy lactation and potentially aid in milk production.

Maternal SARS-CoV-2 impacts fetal placental macrophage programs and placenta-derived microglial models of neurodevelopment.

BACKGROUND: The SARS-CoV-2 virus activates maternal and placental immune responses. Such activation in the setting of other infections during pregnancy is known to impact fetal brain development. The effects of maternal immune activation on neurodevelopment are mediated at least in part by fetal brain microglia. However, microglia are inaccessible for direct analysis, and there are no validated non-invasive surrogate models to evaluate in utero microglial priming and function. We have previously demonstrated shared transcriptional programs between microglia and Hofbauer cells (HBCs, or fetal placental macrophages) in mouse models. METHODS AND RESULTS: We assessed the impact of maternal SARS-CoV-2 on HBCs isolated from 24 term placentas (N = 10 SARS-CoV-2 positive cases, 14 negative controls). Using single-cell RNA-sequencing, we demonstrated that HBC subpopulations exhibit distinct cellular programs, with specific subpopulations differentially impacted by SARS-CoV-2. Assessment of differentially expressed genes implied impaired phagocytosis, a key function of both HBCs and microglia, in some subclusters. Leveraging previously validated models of microglial synaptic pruning, we showed that HBCs isolated from placentas of SARS-CoV-2 positive pregnancies can be transdifferentiated into microglia-like cells (HBC-iMGs), with impaired synaptic pruning behavior compared to HBC models from negative controls. CONCLUSION: These findings suggest that HBCs isolated at birth can be used to create personalized cellular models of offspring microglial programming.

An ex vivo ovulation system enables the discovery of novel ovulatory pathways and nonhormonal contraceptive candidates†.

Ovulation is an integral part of women's menstrual cycle and fertility. Understanding the mechanisms of ovulation has broad implications for the treatment of anovulatory diseases and the development of novel contraceptives. Now, few studies have developed effective models that both faithfully recapitulate the hallmarks of ovulation and possess scalability. We established a three-dimensional encapsulated in vitro follicle growth (eIVFG) system that recapitulates folliculogenesis and produces follicles that undergo ovulation in a controlled manner. Here, we determined whether ex vivo ovulation preserves molecular signatures of ovulation and demonstrated its use in discovering novel ovulatory pathways and nonhormonal contraceptive candidates through a high-throughput ovulation screening. Mature murine follicles from eIVFG were induced to ovulate ex vivo using human chorionic gonadotropin and collected at 0, 1, 4, and 8 hours post-induction. Phenotypic analyses confirmed key ovulatory events, including cumulus expansion, oocyte maturation, follicle rupture, and luteinization. Single-follicle RNA-sequencing analysis revealed the preservation of ovulatory genes and dynamic transcriptomic profiles and signaling. Soft clustering identified distinct gene expression patterns and new pathways that may critically regulate ovulation. We further used this ex vivo ovulation system to screen 21 compounds targeting established and newly identified ovulatory pathways. We discovered that proprotein convertases activate gelatinases to sustain follicle rupture and do not regulate luteinization and progesterone secretion. Together, our ex vivo ovulation system preserves molecular signatures of ovulation, presenting a new powerful tool for studying ovulation and anovulatory diseases as well as for establishing a high-throughput ovulation screening to identify novel nonhormonal contraceptives for women.

Single-cell transcriptomics of staged oocytes and somatic cells reveal novel regulators of follicle activation.

In brief: Proper development of ovarian follicles, comprised of an oocyte and surrounding somatic cells, is essential to support female fertility and endocrine health. Here, we describe a method to isolate single oocytes and somatic cells from the earliest stage follicles, called primordial follicles, and we characterize signals that drive their activation. Abstract: Primordial follicles are the first class of follicles formed in the mammalian ovary and are comprised of an oocyte surrounded by a layer of squamous pre-granulosa cells. This developmental class remains in a non-growing state until individual follicles activate to initiate folliculogenesis. What regulates the timing of follicle activation and the upstream signals that govern these processes are major unanswered questions in ovarian biology. This is partly due to the paucity of data on staged follicle cells since isolating and manipulating individual oocytes and somatic cells from early follicle stages are challenging. To date, most studies on isolated primordial follicles have been conducted on cells collected from animal-age- or oocyte size-specific samples, which encompass multiple follicular stages. Here, we report a method for collecting primordial follicles and their associated oocytes and somatic cells from neonatal murine ovaries using liberase, DNase I, and Accutase. This methodology allows for the identification and collection of follicles immediately post-activation enabling unprecedented interrogation of the primordial-to-primary follicle transition. Molecular profiling by single-cell RNA sequencing revealed that processes including organelle disassembly and cadherin binding were enriched in oocytes and somatic cells as they transitioned from primordial to the primary follicle stage. Furthermore, targets including WNT4, TGFB1, FOXO3, and a network of transcription factors were identified in the transitioning oocytes and somatic cells as potential upstream regulators that collectively may drive follicle activation. Taken together, we have developed a more precise characterization and selection method for studying staged-follicle cells, revealing several novel regulators of early folliculogenesis.

Divergent Functions of Tissue-Resident and Blood-Derived Macrophages in the Hemorrhagic Brain.

Background and Purpose: Brain tissue-resident microglia and monocyte-derived macrophages (MDMs) are innate immune cells that contribute to the inflammatory response, phagocytosis of debris, and tissue repair after injury. We have previously reported that both microglia and MDMs transition from proinflammatory to reparative phenotypes over days after an intracerebral hemorrhage (ICH). However, their individual functional properties in the brain remain largely unknown. Here we characterized the differences between microglia and MDMs and further elucidate their distinct activation states and functional contributions to the pathophysiology and recovery after ICH. Methods: Autologous blood injection was used to model ICH in mice. Longitudinal transcriptomic analyses on isolated microglia and MDMs from mice at days 1, 3, 7 and 10 after ICH and naive controls identified core transcriptional programs that distinguish these cells. Imaging flow cytometry and in vivo phagocytosis assays were used to study phagocytic ability of microglia and MDMs. Antigen presentation was evaluated by ovalbumin-OTII CD4 T-cell proliferation assays with bone marrow­derived macrophages and primary microglia cultures. Results: MDMs had higher phagocytic activity and higher erythrophagocytosis in the ICH brain. Differential gene expression revealed distinct transcriptional signatures in the MDMs and microglia after ICH. MDMs had higher expression of MHCII (major histocompatibility complex class II) genes than microglia at all time points and greater ability to induce antigen-specific T-cell proliferation. Conclusions: The different ontogeny of microglia and MDMs lead to divergent responses and functions in the inflamed brain as these 2 cell populations differ in phagocytic functions and antigen-presenting capabilities in the brain after ICH.

Diverse antiviral IgG effector activities are predicted by unique biophysical antibody features.

BACKGROUND: The critical role of antibody Fc-mediated effector functions in immune defense has been widely reported in various viral infections. These effector functions confer cellular responses through engagement with innate immune cells. The precise mechanism(s) by which immunoglobulin G (IgG) Fc domain and cognate receptors may afford protection are poorly understood, however, in the context of HIV/SHIV infections. Many different in vitro assays have been developed and utilized to measure effector functions, but the extent to which these assays capture distinct antibody activities has not been fully elucidated. RESULTS: In this study, six Fc-mediated effector function assays and two biophysical antibody profiling assays were performed on a common set of samples from HIV-1 infected and vaccinated subjects. Biophysical antibody profiles supported robust prediction of diverse IgG effector functions across distinct Fc-mediated effector function assays. While a number of assays showed correlated activities, supervised machine learning models indicated unique antibody features as primary contributing factors to the associated effector functions. Additional experiments established the mechanistic relevance of relationships discovered using this unbiased approach. CONCLUSIONS: In sum, this study provides better resolution on the diversity and complexity of effector function assays, offering a clearer perspective into this family of antibody mechanisms of action to inform future HIV-1 treatment and vaccination strategies.

Blood handling and leukocyte isolation methods impact the global transcriptome of immune cells.

BACKGROUND: Transcriptional profiling with ultra-low input methods can yield valuable insights into disease, particularly when applied to the study of immune cells using RNA-sequencing. The advent of these methods has allowed for their use in profiling cells collected in clinical trials and other studies that involve the coordination of human-derived material. To date, few studies have sought to quantify what effects that collection and handling of this material can have on resulting data. RESULTS: We characterized the global effects of blood handling, methods for leukocyte isolation, and preservation media on low numbers of immune cells isolated from blood. We found overall that storage/shipping temperature of blood prior to leukocyte isolation and sorting led to global changes in both CD8+ T cells and monocytes, including alterations in immune-related gene sets. We found that the use of a leukocyte filtration system minimized these alterations and we applied this method to generate high-quality transcriptional data from sorted immune cells isolated from the blood of intracerebral hemorrhage patients and matched healthy controls. CONCLUSIONS: Our data underscore the necessity of processing samples with comparably defined protocols prior to transcriptional profiling and demonstrate that a filtration method can be applied to quickly isolate immune cells of interest while minimizing transcriptional bias.

Correction to: blood handling and leukocyte isolation methods impact the global transcriptome of immune cells.

It has been highlighted that the original article [1] contained a typesetting mistake in the middle name of Arthur F. Steinschneider.

Diversity of Antiviral IgG Effector Activities Observed in HIV-Infected and Vaccinated Subjects.

Diverse Ab effector functions mediated by the Fc domain have been commonly associated with reduced risk of infection in a growing number of nonhuman primate and human clinical studies. This study evaluated the anti-HIV Ab effector activities in polyclonal serum samples from HIV-infected donors, VAX004 vaccine recipients, and healthy HIV-negative subjects using a variety of primary and cell line-based assays, including Ab-dependent cellular cytotoxicity (ADCC), Ab-dependent cell-mediated viral inhibition, and Ab-dependent cellular phagocytosis. Additional assay characterization was performed with a panel of Fc-engineered variants of mAb b12. The goal of this study was to characterize different effector functions in the study samples and identify assays that might most comprehensively and dependably capture Fc-mediated Ab functions mediated by different effector cell types and against different viral targets. Deployment of such assays may facilitate assessment of functionally unique humoral responses and contribute to identification of correlates of protection with potential mechanistic significance in future HIV vaccine studies. Multivariate and correlative comparisons identified a set of Ab-dependent cell-mediated viral inhibition and phagocytosis assays that captured different Ab activities and were distinct from a group of ADCC assays that showed a more similar response profile across polyclonal serum samples. The activities of a panel of b12 monoclonal Fc variants further identified distinctions among the ADCC assays. These results reveal the natural diversity of Fc-mediated Ab effector responses among vaccine recipients in the VAX004 trial and in HIV-infected subjects, and they point to the potential importance of polyfunctional Ab responses.

Development of a High-Throughput Functional Screen Using Nanowell-Assisted Cell Patterning.

Living-cell-based screens can facilitate lead discovery of functional therapeutics of interest. A versatile and scalable method is reported that uses dense arrays of nanowells for imparting defined patterns on monolayers of cells. It is shown that this approach can coordinate a multi-component biological assay by designing and implementing a high-throughput, functional nanoliter-scale neutralization assay to identify neutralizing antibodies against HIV.

Automated pipeline for rapid production and screening of HIV-specific monoclonal antibodies using pichia pastoris.

Monoclonal antibodies (mAbs) that bind and neutralize human pathogens have great therapeutic potential. Advances in automated screening and liquid handling have resulted in the ability to discover antigen-specific antibodies either directly from human blood or from various combinatorial libraries (phage, bacteria, or yeast). There remain, however, bottlenecks in the cloning, expression and evaluation of such lead antibodies identified in primary screens that hinder high-throughput screening. As such, "hit-to-lead identification" remains both expensive and time-consuming. By combining the advantages of overlap extension PCR (OE-PCR) and a genetically stable yet easily manipulatable microbial expression host Pichia pastoris, we have developed an automated pipeline for the rapid production and screening of full-length antigen-specific mAbs. Here, we demonstrate the speed, feasibility and cost-effectiveness of our approach by generating several broadly neutralizing antibodies against human immunodeficiency virus (HIV).

Single-cell and spatiotemporal profile of ovulation in the mouse ovary.

Ovulation is a spatiotemporally coordinated process that involves several tightly controlled events, including oocyte meiotic maturation, cumulus expansion, follicle wall rupture and repair, and ovarian stroma remodeling. To date, no studies have detailed the precise window of ovulation at single-cell resolution. Here, we performed parallel single-cell RNA-seq and spatial transcriptomics on paired mouse ovaries across an ovulation time course to map the spatiotemporal profile of ovarian cell types. We show that major ovarian cell types exhibit time-dependent transcriptional states enriched for distinct functions and have specific localization profiles within the ovary. We also identified gene markers for ovulation-dependent cell states and validated these using orthogonal methods. Finally, we performed cell-cell interaction analyses to identify ligand-receptor pairs that may drive ovulation, revealing previously unappreciated interactions. Taken together, our data provides a rich and comprehensive resource of murine ovulation that can be mined for discovery by the scientific community.

Organoid co-culture model of the human endometrium in a fully synthetic extracellular matrix enables the study of epithelial-stromal crosstalk.

BACKGROUND: The human endometrium undergoes recurring cycles of growth, differentiation, and breakdown in response to sex hormones. Dysregulation of epithelial-stromal communication during hormone-mediated signaling may be linked to myriad gynecological disorders for which treatments remain inadequate. Here, we describe a completely defined, synthetic extracellular matrix that enables co-culture of human endometrial epithelial and stromal cells in a manner that captures healthy and disease states across a simulated menstrual cycle. METHODS: We parsed cycle-dependent endometrial integrin expression and matrix composition to define candidate cell-matrix interaction cues for inclusion in a polyethylene glycol (PEG)-based hydrogel crosslinked with matrix metalloproteinase-labile peptides. We semi-empirically screened a parameter space of biophysical and molecular features representative of the endometrium to define compositions suitable for hormone-driven expansion and differentiation of epithelial organoids, stromal cells, and co-cultures of the two cell types. FINDINGS: Each cell type exhibited characteristic morphological and molecular responses to hormone changes when co-encapsulated in hydrogels tuned to a stiffness regime similar to the native tissue and functionalized with a collagen-derived adhesion peptide (GFOGER) and a fibronectin-derived peptide (PHSRN-K-RGD). Analysis of cell-cell crosstalk during interleukin 1B (IL1B)-induced inflammation revealed dysregulation of epithelial proliferation mediated by stromal cells. CONCLUSIONS: Altogether, we demonstrate the development of a fully synthetic matrix to sustain the dynamic changes of the endometrial microenvironment and support its applications to understand menstrual health and endometriotic diseases. FUNDING: This work was supported by The John and Karine Begg Foundation, the Manton Foundation, and NIH U01 (EB029132).

Systematic benchmarking of imaging spatial transcriptomics platforms in FFPE tissues.

Emerging imaging spatial transcriptomics (iST) platforms and coupled analytical methods can recover cell-to-cell interactions, groups of spatially covarying genes, and gene signatures associated with pathological features, and are thus particularly well-suited for applications in formalin fixed paraffin embedded (FFPE) tissues. Here, we benchmarked the performance of three commercial iST platforms on serial sections from tissue microarrays (TMAs) containing 23 tumor and normal tissue types for both relative technical and biological performance. On matched genes, we found that 10x Xenium shows higher transcript counts per gene without sacrificing specificity, but that all three platforms concord to orthogonal RNA-seq datasets and can perform spatially resolved cell typing, albeit with different false discovery rates, cell segmentation error frequencies, and with varying degrees of sub-clustering for downstream biological analyses. Taken together, our analyses provide a comprehensive benchmark to guide the choice of iST method as researchers design studies with precious samples in this rapidly evolving field.

Maternal SARS-CoV-2 impacts fetal placental macrophage programs and placenta-derived microglial models of neurodevelopment.

The SARS-CoV-2 virus activates maternal and placental immune responses, which in the setting of other infections occurring during pregnancy are known to impact fetal brain development. The effects of maternal immune activation on neurodevelopment are mediated at least in part by fetal brain microglia. However, microglia are inaccessible for direct analysis, and there are no validated non-invasive surrogate models to evaluate in utero microglial priming and function. We have previously demonstrated shared transcriptional programs between microglia and Hofbauer cells (HBCs, or fetal placental macrophages) in mouse models. Here, we assessed the impact of maternal SARS-CoV-2 on HBCs isolated from term placentas using single-cell RNA-sequencing. We demonstrated that HBC subpopulations exhibit distinct cellular programs, with specific subpopulations differentially impacted by SARS-CoV-2. Assessment of differentially expressed genes implied impaired phagocytosis, a key function of both HBCs and microglia, in some subclusters. Leveraging previously validated models of microglial synaptic pruning, we showed that HBCs isolated from placentas of SARS-CoV-2 positive pregnancies can be transdifferentiated into microglia-like cells, with altered morphology and impaired synaptic pruning behavior compared to HBC models from negative controls. These findings suggest that HBCs isolated at birth can be used to create personalized cellular models of offspring microglial programming.

Leukocyte dynamics after intracerebral hemorrhage in a living patient reveal rapid adaptations to tissue milieu.

Intracerebral hemorrhage (ICH) is a devastating form of stroke with a high mortality rate and few treatment options. Discovery of therapeutic interventions has been slow given the challenges associated with studying acute injury in the human brain. Inflammation induced by exposure of brain tissue to blood appears to be a major part of brain tissue injury. Here, we longitudinally profiled blood and cerebral hematoma effluent from a patient enrolled in the Minimally Invasive Surgery with Thrombolysis in Intracerebral Hemorrhage Evacuation trial, offering a rare window into the local and systemic immune responses to acute brain injury. Using single-cell RNA-Seq (scRNA-Seq), this is the first report to our knowledge that characterized the local cellular response during ICH in the brain of a living patient at single-cell resolution. Our analysis revealed shifts in the activation states of myeloid and T cells in the brain over time, suggesting that leukocyte responses are dynamically reshaped by the hematoma microenvironment. Interestingly, the patient had an asymptomatic rebleed that our transcriptional data indicated occurred prior to detection by CT scan. This case highlights the rapid immune dynamics in the brain after ICH and suggests that sensitive methods such as scRNA-Seq would enable greater understanding of complex intracerebral events.

Longitudinal transcriptomics define the stages of myeloid activation in the living human brain after intracerebral hemorrhage.

Opportunities to interrogate the immune responses in the injured tissue of living patients suffering from acute sterile injuries such as stroke and heart attack are limited. We leveraged a clinical trial of minimally invasive neurosurgery for patients with intracerebral hemorrhage (ICH), a severely disabling subtype of stroke, to investigate the dynamics of inflammation at the site of brain injury over time. Longitudinal transcriptional profiling of CD14+ monocytes/macrophages and neutrophils from hematomas of patients with ICH revealed that the myeloid response to ICH within the hematoma is distinct from that in the blood and occurs in stages conserved across the patient cohort. Initially, hematoma myeloid cells expressed a robust anabolic proinflammatory profile characterized by activation of hypoxia-inducible factors (HIFs) and expression of genes encoding immune factors and glycolysis. Subsequently, inflammatory gene expression decreased over time, whereas anti-inflammatory circuits were maintained and phagocytic and antioxidative pathways up-regulated. During this transition to immune resolution, glycolysis gene expression and levels of the potent proresolution lipid mediator prostaglandin E2 remained elevated in the hematoma, and unexpectedly, these elevations correlated with positive patient outcomes. Ex vivo activation of human macrophages by ICH-associated stimuli highlighted an important role for HIFs in production of both inflammatory and anti-inflammatory factors, including PGE2, which, in turn, augmented VEGF production. Our findings define the time course of myeloid activation in the human brain after ICH, revealing a conserved progression of immune responses from proinflammatory to proresolution states in humans after brain injury and identifying transcriptional programs associated with neurological recovery.

Single-cell transcriptomic atlas of the human retina identifies cell types associated with age-related macular degeneration.

Genome-wide association studies (GWAS) have identified genetic variants associated with age-related macular degeneration (AMD), one of the leading causes of blindness in the elderly. However, it has been challenging to identify the cell types associated with AMD given the genetic complexity of the disease. Here we perform massively parallel single-cell RNA sequencing (scRNA-seq) of human retinas using two independent platforms, and report the first single-cell transcriptomic atlas of the human retina. Using a multi-resolution network-based analysis, we identify all major retinal cell types, and their corresponding gene expression signatures. Heterogeneity is observed within macroglia, suggesting that human retinal glia are more diverse than previously thought. Finally, GWAS-based enrichment analysis identifies glia, vascular cells, and cone photoreceptors to be associated with the risk of AMD. These data provide a detailed analysis of the human retina, and show how scRNA-seq can provide insight into cell types involved in complex, inflammatory genetic diseases.