Retinoic acid is dispensable for meiotic initiation but required for spermiogenesis in the mammalian testis.

Retinoic acid (RA) is the proposed mammalian 'meiosis inducing substance'. However, evidence for this role comes from studies in the fetal ovary, where germ cell differentiation and meiotic initiation are temporally inseparable. In the postnatal testis, these events are separated by more than 1 week. Exploiting this difference, we discovered that, although RA is required for spermatogonial differentiation, it is dispensable for the subsequent initiation, progression and completion of meiosis. Indeed, in the absence of RA, the meiotic transcriptome program in both differentiating spermatogonia and spermatocytes entering meiosis was largely unaffected. Instead, transcripts encoding factors required during spermiogenesis were aberrant during preleptonema, and the subsequent spermatid morphogenesis program was disrupted such that no sperm were produced. Taken together, these data reveal a RA-independent model for male meiotic initiation.

Modeling mammalian spermatogonial differentiation and meiotic initiation in vitro.

In mammalian testes, premeiotic spermatogonia respond to retinoic acid by completing an essential lengthy differentiation program before initiating meiosis. The molecular and cellular changes directing these developmental processes remain largely undefined. This wide gap in knowledge is due to two unresolved technical challenges: (1) lack of robust and reliable in vitro models to study differentiation and meiotic initiation; and (2) lack of methods to isolate large and pure populations of male germ cells at each stage of differentiation and at meiotic initiation. Here, we report a facile in vitro differentiation and meiotic initiation system that can be readily manipulated, including the use of chemical agents that cannot be safely administered to live animals. In addition, we present a transgenic mouse model enabling fluorescence-activated cell sorting-based isolation of millions of spermatogonia at specific developmental stages as well as meiotic spermatocytes.

The germ cell-specific RNA binding protein RBM46 is essential for spermatogonial differentiation in mice.

Control over gene expression is exerted, in multiple stages of spermatogenesis, at the post-transcriptional level by RNA binding proteins (RBPs). We identify here an essential role in mammalian spermatogenesis and male fertility for 'RNA binding protein 46' (RBM46). A highly evolutionarily conserved gene, Rbm46 is also essential for fertility in both flies and fish. We found Rbm46 expression was restricted to the mouse germline, detectable in males in the cytoplasm of premeiotic spermatogonia and meiotic spermatocytes. To define its requirement for spermatogenesis, we generated Rbm46 knockout (KO, Rbm46-/-) mice; although male Rbm46-/- mice were viable and appeared grossly normal, they were infertile. Testes from adult Rbm46-/- mice were small, with seminiferous tubules containing only Sertoli cells and few undifferentiated spermatogonia. Using genome-wide unbiased high throughput assays RNA-seq and 'enhanced crosslinking immunoprecipitation' coupled with RNA-seq (eCLIP-seq), we discovered RBM46 could bind, via a U-rich conserved consensus sequence, to a cohort of mRNAs encoding proteins required for completion of differentiation and subsequent meiotic initiation. In summary, our studies support an essential role for RBM46 in regulating target mRNAs during spermatogonia differentiation prior to the commitment to meiosis in mice.

Lipid metabolism in dopaminergic neurons influences light entrainment.

Dietary lipids, particularly omega-3 polyunsaturated fatty acids, are speculated to impact behaviors linked to the dopaminergic system, such as movement and control of circadian rhythms. However, the ability to draw a direct link between dopaminergic omega-3 fatty acid metabolism and behavioral outcomes has been limited to the use of diet-based approaches, which are confounded by systemic effects. Here, neuronal lipid metabolism was targeted in a diet-independent manner by manipulation of long-chain acyl-CoA synthetase 6 (ACSL6) expression. ACSL6 performs the initial reaction for cellular fatty acid metabolism and prefers the omega-3 polyunsaturated fatty acid, docosahexaenoic acid (DHA). The loss of Acsl6 in mice (Acsl6-/- ) depletes neuronal membranes of DHA content and results in phenotypes linked to dopaminergic control, such as hyperlocomotion, impaired short-term spatial memory, and imbalances in dopamine neurochemistry. To investigate the role of dopaminergic ACSL6 on these outcomes, a dopaminergic neuron-specific ACSL6 knockout mouse was generated (Acsl6DA-/- ). Acsl6DA-/- mice demonstrated hyperlocomotion and imbalances in striatal dopamine neurochemistry. Circadian rhythms of both the Acsl6-/- and the Acsl6DA-/- mice were similar to control mice under basal conditions. However, upon light entrainment, a mimetic of jet lag, both the complete knockout of ACSL6 and the dopaminergic-neuron-specific loss of ACSL6 resulted in a longer recovery to entrainment compared to control mice. In conclusion, these data demonstrate that ACSL6 in dopaminergic neurons alters dopamine metabolism and regulation of light entrainment suggesting that DHA metabolism mediated by ACSL6 plays a role in dopamine neuron biology.

Granulocyte-colony stimulating factor does not prevent in vitro cisplatin-induced germ cell reduction in immature human and mouse testis.

BACKGROUND: Currently there are no established fertility preservation options for pre-pubertal boys facing cancer treatment. Granulocyte-colony stimulating factor (G-CSF) treatment has been proposed to be chemoprotective against spermatogonial cell loss in an alkylating chemotherapy model of busulfan treated adult mice. Having previously shown that exposure to the alkylating-like chemotherapy cisplatin resulted in a reduction in germ cell numbers in immature human testicular tissues, we here investigate whether G-CSF would prevent cisplatin-induced germ cell loss in immature human and mouse (fetal and pre-pubertal) testicular tissues. METHODS: Organotypic in vitro culture systems were utilised to determine the effects of clinically-relevant concentrations of G-CSF in cisplatin-exposed immature testicular tissues. Human fetal (n = 14 fetuses) and mouse pre-pubertal (n = 4 litters) testicular tissue pieces were cultured and exposed to cisplatin or vehicle control for 24 hrs and analysed at 72 and 240 hrs post-exposure. Combined G-CSF and cisplatin exposure groups explored varying concentrations and duration of G-CSF supplementation to the culture medium (including groups receiving G-CSF before, during and after cisplatin exposure). In addition, effects of G-CSF supplementation alone were investigated. Survival of total germ cell and sub-populations were identified by expression of AP2γ and MAGE-A4 for human gonocytes and (pre)spermatogonia, respectively, and MVH and PLZF, for mouse germ cells and putative spermatogonial stem cells (SSCs) respectively, were quantified. RESULTS: Exposure to cisplatin resulted in a reduced germ cell number in human fetal and mouse pre-pubertal testicular tissues at 240 hrs post-exposure. Germ cell number was not preserved by combined exposure with G-CSF using any of the exposure regimens (prior to, during or after cisplatin exposure). Continuous supplementation with G-CSF alone for 14 days did not change the germ cell composition in either human or mouse immature testicular tissues. CONCLUSIONS: This study demonstrates that exposure to G-CSF does not prevent cisplatin-induced germ cell loss in immature human and mouse testicular tissues in an in vitro system.

Conserved Transcriptome Features Define Prepubertal Primate Spermatogonial Stem Cells as Adark Spermatogonia and Identify Unique Regulators.

Antineoplastic treatments for cancer and other non-malignant disorders can result in long-term or permanent male infertility by ablating spermatogonial stem cells (SSCs). SSC transplantation using testicular tissue harvested before a sterilizing treatment is a promising approach for restoring male fertility in these cases, but a lack of exclusive biomarkers to unequivocally identify prepubertal SSCs limits their therapeutic potential. To address this, we performed single-cell RNA-seq on testis cells from immature baboons and macaques and compared these cells with published data from prepubertal human testis cells and functionally-defined mouse SSCs. While we found discrete groups of human spermatogonia, baboon and rhesus spermatogonia appeared less heterogenous. A cross-species analysis revealed cell types analogous to human SSCs in baboon and rhesus germ cells, but a comparison with mouse SSCs revealed significant differences with primate SSCs. Primate-specific SSC genes were enriched for components and regulators of the actin cytoskeleton and participate in cell-adhesion, which may explain why the culture conditions for rodent SSCs are not appropriate for primate SSCs. Furthermore, correlating the molecular definitions of human SSC, progenitor and differentiating spermatogonia with the histological definitions of Adark/Apale spermatogonia indicates that both SSCs and progenitor spermatogonia are Adark, while Apale spermatogonia appear biased towards differentiation. These results resolve the molecular identity of prepubertal human SSCs, define novel pathways that could be leveraged for advancing their selection and propagation in vitro, and confirm that the human SSC pool resides entirely within Adark spermatogonia.

Differential RA responsiveness directs formation of functionally distinct spermatogonial populations at the initiation of spermatogenesis in the mouse.

In the mammalian testis, sustained spermatogenesis relies on spermatogonial stem cells (SSCs); their progeny either remain as stem cells (self-renewal) or proliferate and differentiate to enter meiosis in response to retinoic acid (RA). Here, we sought to uncover elusive mechanisms regulating a key switch fundamental to spermatogonial fate: the capacity of spermatogonia to respond to RA. Using the developing mouse testis as a model, we found that spermatogonia and precursor prospermatogonia exhibit a heterogeneous capacity to respond to RA with at least two underlying causes. First, progenitor spermatogonia are prevented from responding to RA by catabolic activity of cytochrome P450 family 26 enzymes. Second, a smaller subset of undifferentiated spermatogonia enriched for SSCs exhibit catabolism-independent RA insensitivity. Moreover, for the first time, we observed that precursor prospermatogonia are heterogeneous and comprise subpopulations that exhibit the same differential RA responsiveness found in neonatal spermatogonia. We propose a novel model by which mammalian prospermatogonial and spermatogonial fates are regulated by their intrinsic capacity to respond (or not) to the differentiation signal provided by RA before, and concurrent with, the initiation of spermatogenesis.

Dynamic cytoplasmic projections connect mammalian spermatogonia in vivo.

Throughout the male reproductive lifespan, spermatogonial stem cells (SSCs) produce committed progenitors that proliferate and then remain physically connected in growing clones via short cylindrical intercellular bridges (ICBs). These ICBs, which enlarge in meiotic spermatocytes, have been demonstrated to provide a conduit for postmeiotic haploid spermatids to share sex chromosome-derived gene products. In addition to ICBs, spermatogonia exhibit multiple thin cytoplasmic projections. Here, we have explored the nature of these projections in mice and find that they are dynamic, span considerable distances from their cell body (≥25 µm), either terminate or physically connect multiple adjacent spermatogonia, and allow for sharing of macromolecules. Our results extend the current model that subsets of spermatogonia exist as isolated cells or clones, and support a model in which spermatogonia of similar developmental fates are functionally connected through a shared dynamic cytoplasm mediated by thin cytoplasmic projections.

Differential RA responsiveness among subsets of mouse late progenitor spermatogonia.

Initiation of spermatogonial differentiation in the mouse testis begins with the response to retinoic acid (RA) characterized by activation of KIT and STRA8 expression. In the adult, spermatogonial differentiation is spatiotemporally coordinated by a pulse of RA every 8.6 days that is localized to stages VII-VIII of the seminiferous epithelial cycle. Dogmatically, progenitor spermatogonia that express retinoic acid receptor gamma (RARG) at these stages will differentiate in response to RA, but this has yet to be tested functionally. Previous single-cell RNA-seq data identified phenotypically and functionally distinct subsets of spermatogonial stem cells (SSCs) and progenitor spermatogonia, where late progenitor spermatogonia were defined by expression of RARG and Dppa3. Here, we found late progenitor spermatogonia (RARGhigh KIT-) were further divisible into two subpopulations based on Dppa3 reporter expression (Dppa3-ECFP or Dppa3-EGFP) and were observed across all stages of the seminiferous epithelial cycle. However, nearly all Dppa3+ spermatogonia were differentiating (KIT+) late in the seminiferous epithelial cycle (stages X-XII), while Dppa3- late progenitors remained abundant, suggesting that Dppa3+ and Dppa3- late progenitors differentially responded to RA. Following acute RA treatment (2-4 h), significantly more Dppa3+ late progenitors induced KIT, including at the midpoint of the cycle (stages VI-IX), than Dppa3- late progenitors. Subsequently, single-cell analyses indicated a subset of Dppa3+ late progenitors expressed higher levels of Rxra, which we confirmed by RXRA whole-mount immunostaining. Together, these results indicate RARG alone is insufficient to initiate a spermatogonial response to RA in the adult mouse testis and suggest differential RXRA expression may discriminate responding cells.

Acyl-CoA synthetase 6 enriches seminiferous tubules with the ω-3 fatty acid docosahexaenoic acid and is required for male fertility in the mouse.

Docosahexaenoic acid (DHA) is an ω-3 dietary-derived polyunsaturated fatty acid of marine origin enriched in testes and necessary for normal fertility, yet the mechanisms regulating the enrichment of DHA in the testes remain unclear. Long-chain ACSL6 (acyl-CoA synthetase isoform 6) activates fatty acids for cellular anabolic and catabolic metabolism by ligating a CoA to a fatty acid, is highly expressed in testes, and has high preference for DHA. Here, we investigated the role of ACSL6 for DHA enrichment in the testes and its requirement for male fertility. Acsl6-/- males were severely subfertile with smaller testes, reduced cauda epididymal sperm counts, germ cell loss, and disorganization of the seminiferous epithelium. Total fatty acid profiling of Acsl6-/- testes revealed reduced DHA and increased ω-6 arachidonic acid, a fatty acid profile also reflected in phospholipid composition. Strikingly, lipid imaging demonstrated spatial redistribution of phospholipids in Acsl6-/- testes. Arachidonic acid-containing phospholipids were predominantly interstitial in control testes but diffusely localized across Acsl6-/- testes. In control testes, DHA-containing phospholipids were predominantly within seminiferous tubules, which contain Sertoli cells and spermatogenic cells but relocalized to the interstitium in Acsl6-/- testes. Taken together, these data demonstrate that ACSL6 is an initial driving force for germ cell DHA enrichment and is required for normal spermatogenesis and male fertility.

The regulatory repertoire of PLZF and SALL4 in undifferentiated spermatogonia.

Spermatogonial stem cells (SSCs) maintain spermatogenesis throughout adulthood through balanced self-renewal and differentiation, yet the regulatory logic of these fate decisions is poorly understood. The transcription factors Sal-like 4 (SALL4) and promyelocytic leukemia zinc finger (PLZF; also known as ZBTB16) are known to be required for normal SSC function, but their targets are largely unknown. ChIP-seq in mouse THY1(+) spermatogonia identified 4176 PLZF-bound and 2696 SALL4-bound genes, including 1149 and 515 that were unique to each factor, respectively, and 1295 that were bound by both factors. PLZF and SALL4 preferentially bound gene promoters and introns, respectively. Motif analyses identified putative PLZF and SALL4 binding sequences, but rarely both at shared sites, indicating significant non-autonomous binding in any given cell. Indeed, the majority of PLZF/SALL4 shared sites contained only PLZF motifs. SALL4 also bound gene introns at sites containing motifs for the differentiation factor DMRT1. Moreover, mRNA levels for both unique and shared target genes involved in both SSC self-renewal and differentiation were suppressed following SALL4 or PLZF knockdown. Together, these data reveal the full profile of PLZF and SALL4 regulatory targets in undifferentiated spermatogonia, including SSCs, which will help elucidate mechanisms controlling the earliest cell fate decisions in spermatogenesis.

What has single-cell RNA-seq taught us about mammalian spermatogenesis?

Mammalian spermatogenesis is a complex developmental program that transforms mitotic testicular germ cells (spermatogonia) into mature male gametes (sperm) for production of offspring. For decades, it has been known that this several-weeks-long process involves a series of highly ordered and morphologically recognizable cellular changes as spermatogonia proliferate, spermatocytes undertake meiosis, and spermatids develop condensed nuclei, acrosomes, and flagella. Yet, much of the underlying molecular logic driving these processes has remained opaque because conventional characterization strategies often aggregated groups of cells to meet technical requirements or due to limited capability for cell selection. Recently, a cornucopia of single-cell transcriptome studies has begun to lift the veil on the full compendium of gene expression phenotypes and changes underlying spermatogenic development. These datasets have revealed the previously obscured molecular heterogeneity among and between varied spermatogenic cell types and are reinvigorating investigation of testicular biology. This review describes the extent of available single-cell RNA-seq profiles of spermatogenic and testicular somatic cells, how those data were produced and evaluated, their present value for advancing knowledge of spermatogenesis, and their potential future utility at both the benchtop and bedside.

Granulocyte colony-stimulating factor (G-CSF) promotes spermatogenic regeneration from surviving spermatogonia after high-dose alkylating chemotherapy.

BACKGROUND: The lifesaving chemotherapy and radiation treatments that allow patients to survive cancer can also result in a lifetime of side-effects, including male infertility. Infertility in male cancer survivors is thought to primarily result from killing of the spermatogonial stem cells (SSCs) responsible for producing spermatozoa since SSCs turn over slowly and are thereby sensitive to antineoplastic therapies. We previously demonstrated that the cytokine granulocyte colony-stimulating factor (G-CSF) can preserve spermatogenesis after alkylating chemotherapy (busulfan). METHODS: Male mice were treated with G-CSF or controls before and/or after sterilizing busulfan treatment and evaluated immediately or 10-19 weeks later for effects on spermatogenesis. RESULTS: We demonstrated that the protective effect of G-CSF on spermatogenesis was stable for at least 19 weeks after chemotherapy, nearly twice as long as previously shown. Further, G-CSF treatment enhanced spermatogenic measures 10 weeks after treatment in the absence of a cytotoxic insult, suggesting G-CSF acts as a mitogen in steady-state spermatogenesis. In agreement with this conclusion, G-CSF treatment for 3 days before busulfan treatment exacerbated the loss of spermatogenesis observed with G-CSF alone. Reciprocally, spermatogenic recovery was modestly enhanced in mice treated with G-CSF for 4 days after busulfan. These results suggested that G-CSF promoted spermatogonial proliferation, leading to enhanced spermatogenic regeneration from surviving SSCs. Similarly, there was a significant increase in proportion of PLZF+ undifferentiated spermatogonia that were Ki67+ (proliferating) 1 day after G-CSF treatment. CONCLUSIONS: Together, these results clarify that G-CSF protects spermatogenesis after alkylating chemotherapy by stimulating proliferation of surviving spermatogonia, and indicate it may be useful as a retrospective fertility-restoring treatment.

TSPAN8 Expression Distinguishes Spermatogonial Stem Cells in the Prepubertal Mouse Testis.

Precise separation of spermatogonial stem cells (SSCs) from progenitor spermatogonia that lack stem cell activity and are committed to differentiation remains a challenge. To distinguish between these spermatogonial subtypes, we identified genes that exhibited bimodal mRNA levels at the single-cell level among undifferentiated spermatogonia from Postnatal Day 6 mouse testes, including Tspan8, Epha2, and Pvr, each of which encode cell surface proteins useful for cell selection. Transplantation studies provided definitive evidence that a TSPAN8-high subpopulation is enriched for SSCs. RNA-seq analyses identified genes differentially expressed between TSPAN8-high and -low subpopulations that clustered into multiple biological pathways potentially involved in SSC renewal or differentiation, respectively. Methyl-seq analysis identified hypomethylated domains in the promoters of these genes in both subpopulations that colocalized with peaks of histone modifications defined by ChIP-seq analysis. Taken together, these results demonstrate functional heterogeneity among mouse undifferentiated spermatogonia and point to key biological characteristics that distinguish SSCs from progenitor spermatogonia.

Transcriptional and translational heterogeneity among neonatal mouse spermatogonia.

Spermatogonial stem cells (SSCs) are a subset of undifferentiated spermatogonia responsible for ongoing spermatogenesis in mammalian testes. Spermatogonial stem cells arise from morphologically homogeneous prospermatogonia, but growing evidence suggests that only a subset of prospermatogonia develops into the foundational SSC pool. This predicts that subtypes of undifferentiated spermatogonia with discrete mRNA and protein signatures should be distinguishable in neonatal testes. We used single-cell quantitative RT-PCR to examine mRNA levels of 172 genes in individual spermatogonia from 6-day postnatal (P6) mouse testes. Cells enriched from P6 testes using the StaPut or THY1(+) magnetic cell sorting methods exhibited considerable heterogeneity in the abundance of specific germ cell and stem cell mRNAs, segregating into one somatic and three distinct spermatogonial clusters. However, P6 Id4-eGFP(+) transgenic spermatogonia, which are known to be enriched for SSCs, were more homogeneous in their mRNA levels, exhibiting uniform levels for the majority of genes examined (122 of 172). Interestingly, these cells displayed nonuniform (50 of 172) expression of a smaller cohort of these genes, suggesting there is substantial heterogeneity even within the Id4-eGFP(+) population. Further, although immunofluorescence staining largely demonstrated conformity between mRNA and protein levels, some proteins were observed in patterns that were disparate from those detected for the corresponding mRNAs in Id4-eGFP(+) spermatogonia (e.g., Kit, Sohlh2, Stra8), suggesting additional heterogeneity is introduced at the posttranscriptional level. Taken together, these data demonstrate the existence of multiple spermatogonial subtypes in P6 mouse testes and raise the intriguing possibility that these subpopulations may correlate with the development of functionally distinct spermatogenic cell types.

Introduction: The Why's and How's for Studying Spermatogenesis and Spermatogonial Stem Cells.

Spermatogenesis is maintained throughout adulthood by a pool of adult stem cells termed spermatogonial stem cells (SSCs). Research investigations into spermatogenesis can provide insight into the etiology of certain types of male infertility (e.g., Sertoli cell only syndrome), elucidate means of improving food animal production, reveal new therapeutic avenues to address naturally occurring defects in sperm production, mitigate iatrogenic male infertility (e.g., arising from cancer therapy), and potentially intervene for male contraception. This chapter will serve as a commentary about why studying spermatogenesis is important, including a high-level overview of spermatogonia and SSCs, and make the case for a critical need for use of stringent definitions for SSCs and experimental platforms that allow for clear distinction of the multiple types of spermatogonia that exist in testes of mammals.

Single-Molecule Fluorescence In Situ Hybridization for Spatial Detection of mRNAs in Sections of Mammalian Testes.

Single-molecule fluorescence in situ hybridization (smFISH) enables the detection and localization of individual mRNAs in tissue sections with single-molecule resolution while preserving spatial context, and thus, is a useful tool for examining gene expression in biological systems. In particular, the growing reliance on single-cell RNA sequencing (scRNA-seq) to explore cellular heterogeneity has reinvigorated this approach as a validation tool to spatially re-map mRNA expression patterns described in isolated cells to their parent tissue. While use of antibody-based methods, such as indirect immunofluorescence (IIF), remain popular as validation strategies, smFISH often affords superior specificity and maintains congruency with scRNA-seq. Here, we present a detailed protocol that combines multiplexed smFISH using the RNAscope approach with IIF to co-visualize mRNAs and proteins within sections of mouse testes. We provide step-by-step guidelines from testis preparation through visualization that enables mapping of combinations of up to four mRNA/protein targets in frozen sections on the RNAscope platform.

Bulk and Single-Cell RNA-Seq Analyses for Studies of Spermatogonia.

Robust methods have been developed that leverage next-generation sequencing (NGS) to measure abundance of all mRNAs (RNA-seq) in samples as small as individual cells in order to study the testicular transcriptome in mammals. In this chapter, we present robust options for implementing bioinformatics workflows for the analysis of bulk RNA-seq from aggregate samples of hundreds to millions of cells and single-cell RNA-seq from individual cells. We also provide detailed protocols for using the R packages DESeq2 and Seurat, important parameters for successful implementation, and considerations for drawing conclusions from the results.

Isolation of Undifferentiated Spermatogonia from Adult and Developing Mouse Testes.

In the mammalian testis, the mitotic complements of spermatogenic cells are spermatogonia, including spermatogonial stem cells (SSCs) which form the basis of life-long spermatogenesis and male fertility. Thus, investigating spermatogonia and subdivisions thereof is essential to increase our understanding of male germline development and infertility. This protocol describes the isolation of spermatogonia from both adult and developing [postnatal day 6 (P6)] mouse testes. Cell suspensions of the adult mouse testis from the Id4-Egfp transgenic mouse line are obtained through a two-step enzymatic digestion and are subjected to Percoll pre-enrichment before spermatogonia are isolated by selecting testis cells that are CD9bright and ID4-EGFP+ through FACS. For P6 mice, the testis is digested using trypsin-DNase, and spermatogonia are isolated by FACS selection of ID4-EGFP+ testis cells. In both cases, nearly pure populations of undifferentiated spermatogonia are obtained that can be further subdivided using additional parameters (e.g., EGFP intensity, cell surface protein immunostaining), and recovered for use in various downstream applications, such as biochemical analyses (e.g., transcriptome/epigenome), functional analyses by SSC transplantation or propagation in vitro.

Cellular and transcriptome signatures unveiled by single-cell RNA-Seq following ex vivo infection of murine splenocytes with Borrelia burgdorferi.

Introduction: Lyme disease, the most common tick-borne infectious disease in the US, is caused by a spirochetal pathogen Borrelia burgdorferi (Bb). Distinct host responses are observed in susceptible and resistant strains of inbred of mice following infection with Bb reflecting a subset of inflammatory responses observed in human Lyme disease. The advent of post-genomic methodologies and genomic data sets enables dissecting the host responses to advance therapeutic options for limiting the pathogen transmission and/or treatment of Lyme disease. Methods: In this study, we used single-cell RNA-Seq analysis in conjunction with mouse genomics exploiting GFP-expressing Bb to sort GFP+ splenocytes and GFP- bystander cells to uncover novel molecular and cellular signatures that contribute to early stages of immune responses against Bb. Results: These data decoded the heterogeneity of splenic neutrophils, macrophages, NK cells, B cells, and T cells in C3H/HeN mice in response to Bb infection. Increased mRNA abundance of apoptosis-related genes was observed in neutrophils and macrophages clustered from GFP+ splenocytes. Moreover, complement-mediated phagocytosis-related genes such as C1q and Ficolin were elevated in an inflammatory macrophage subset, suggesting upregulation of these genes during the interaction of macrophages with Bb-infected neutrophils. In addition, the role of DUSP1 in regulating the expression of Casp3 and pro-inflammatory cytokines Cxcl1, Cxcl2, Il1b, and Ccl5 in Bb-infected neutrophils were identified. Discussion: These findings serve as a growing catalog of cell phenotypes/biomarkers among murine splenocytes that can be exploited for limiting spirochetal burden to limit the transmission of the agent of Lyme disease to humans via reservoir hosts.