Skip to main content

Main menu

  • HOME
  • CONTENT
    • Early Release
    • Featured
    • Current Issue
    • Issue Archive
    • Collections
    • Podcast
  • ALERTS
  • FOR AUTHORS
    • Information for Authors
    • Fees
    • Journal Clubs
    • eLetters
    • Submit
  • EDITORIAL BOARD
  • ABOUT
    • Overview
    • Advertise
    • For the Media
    • Rights and Permissions
    • Privacy Policy
    • Feedback
  • SUBSCRIBE

User menu

  • Log in
  • My Cart

Search

  • Advanced search
Journal of Neuroscience
  • Log in
  • My Cart
Journal of Neuroscience

Advanced Search

Submit a Manuscript
  • HOME
  • CONTENT
    • Early Release
    • Featured
    • Current Issue
    • Issue Archive
    • Collections
    • Podcast
  • ALERTS
  • FOR AUTHORS
    • Information for Authors
    • Fees
    • Journal Clubs
    • eLetters
    • Submit
  • EDITORIAL BOARD
  • ABOUT
    • Overview
    • Advertise
    • For the Media
    • Rights and Permissions
    • Privacy Policy
    • Feedback
  • SUBSCRIBE
PreviousNext
Research Articles, Cellular/Molecular

Sperm Transcriptional State Associated with Paternal Transmission of Stress Phenotypes

Ashley M. Cunningham, Deena M. Walker, Aarthi Ramakrishnan, Marie A. Doyle, Rosemary C. Bagot, Hannah M. Cates, Catherine J. Peña, Orna Issler, Casey K. Lardner, Caleb Browne, Scott J. Russo, Li Shen and Eric J. Nestler
Journal of Neuroscience 21 July 2021, 41 (29) 6202-6216; DOI: https://doi.org/10.1523/JNEUROSCI.3192-20.2021
Ashley M. Cunningham
Icahn School of Medicine at Mount Sinai, Nash Family Department of Neuroscience and Friedman Brain Institute, New York, New York 10029
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Ashley M. Cunningham
Deena M. Walker
Icahn School of Medicine at Mount Sinai, Nash Family Department of Neuroscience and Friedman Brain Institute, New York, New York 10029
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Aarthi Ramakrishnan
Icahn School of Medicine at Mount Sinai, Nash Family Department of Neuroscience and Friedman Brain Institute, New York, New York 10029
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Aarthi Ramakrishnan
Marie A. Doyle
Icahn School of Medicine at Mount Sinai, Nash Family Department of Neuroscience and Friedman Brain Institute, New York, New York 10029
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Rosemary C. Bagot
Icahn School of Medicine at Mount Sinai, Nash Family Department of Neuroscience and Friedman Brain Institute, New York, New York 10029
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Hannah M. Cates
Icahn School of Medicine at Mount Sinai, Nash Family Department of Neuroscience and Friedman Brain Institute, New York, New York 10029
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Catherine J. Peña
Icahn School of Medicine at Mount Sinai, Nash Family Department of Neuroscience and Friedman Brain Institute, New York, New York 10029
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Orna Issler
Icahn School of Medicine at Mount Sinai, Nash Family Department of Neuroscience and Friedman Brain Institute, New York, New York 10029
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Casey K. Lardner
Icahn School of Medicine at Mount Sinai, Nash Family Department of Neuroscience and Friedman Brain Institute, New York, New York 10029
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Caleb Browne
Icahn School of Medicine at Mount Sinai, Nash Family Department of Neuroscience and Friedman Brain Institute, New York, New York 10029
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Caleb Browne
Scott J. Russo
Icahn School of Medicine at Mount Sinai, Nash Family Department of Neuroscience and Friedman Brain Institute, New York, New York 10029
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Scott J. Russo
Li Shen
Icahn School of Medicine at Mount Sinai, Nash Family Department of Neuroscience and Friedman Brain Institute, New York, New York 10029
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Eric J. Nestler
Icahn School of Medicine at Mount Sinai, Nash Family Department of Neuroscience and Friedman Brain Institute, New York, New York 10029
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Eric J. Nestler
  • Article
  • Figures & Data
  • Info & Metrics
  • eLetters
  • PDF
Loading

Abstract

Paternal stress can induce long-lasting changes in germ cells potentially propagating heritable changes across generations. To date, no studies have investigated differences in transmission patterns between stress-resilient and stress-susceptible mice. We tested the hypothesis that transcriptional alterations in sperm during chronic social defeat stress (CSDS) transmit increased susceptibility to stress phenotypes to the next generation. We demonstrate differences in offspring from stressed fathers that depend on paternal category (resilient vs susceptible) and offspring sex. Importantly, artificial insemination (AI) reveals that sperm mediates some of the behavioral phenotypes seen in offspring. Using RNA-sequencing (RNA-seq), we report substantial and distinct changes in the transcriptomic profiles of sperm following CSDS in susceptible versus resilient fathers, with alterations in long noncoding RNAs (lncRNAs) predominating especially in susceptibility. Correlation analysis revealed that these alterations were accompanied by a loss of regulation of protein-coding genes by lncRNAs in sperm of susceptible males. We also identify several co-expression gene modules that are enriched in differentially expressed genes (DEGs) in sperm from either resilient or susceptible fathers. Taken together, these studies advance our understanding of intergenerational epigenetic transmission of behavioral experience.

SIGNIFICANCE STATEMENT This manuscript contributes to the complex factors that influence the paternal transmission of stress phenotypes. By leveraging the segregation of males exposed to chronic social defeat stress (CSDS) into either resilient or susceptible categories we were able to identify the phenotypic differences in the paternal transmission of stress phenotypes across generations between the two lineages. Importantly, this work also alludes to the significance of both long noncoding RNAs (lncRNAs) and protein coding genes (PCGs) mediating the paternal transmission of stress. The knowledge gained from these data are of particular interest in understanding the risk for the development of psychiatric disorders such as anxiety and depression.

  • anxiety
  • intergenerational stress
  • RNA-sequencing
  • sperm

Introduction

Exposure to adverse experiences has been shown to induce long-lasting, heritable alterations in the epigenome of germ cells, potentially relaying information about the paternal environment to offspring across multiple generations (Franklin et al., 2010; Dietz et al., 2011; Bale, 2015; de Castro Barbosa et al., 2015; Cunningham et al., 2019; Manners et al., 2019; Gapp et al., 2020). Paternal exposure to a variety of stressors has been shown to influence offspring behavioral and metabolic phenotypes (Rodgers et al., 2013; Dickson et al., 2018; Manners et al., 2019). Altered offspring behavioral phenotypes are accompanied by changes in the transcriptional landscape of sperm and, importantly, experimental manipulation of small noncoding RNAs in zygotes of control offspring recapitulates altered behavioral phenotypes seen in offspring from stress lineages (Franklin et al., 2010; Gapp et al., 2014, 2020; Rodgers et al., 2015). Although long noncoding RNAs (lncRNAs) have been shown to be involved in cellular and molecular stages of development such as imprinting (Sleutels et al., 2002) and X chromosome inactivation (Brockdorff et al., 1992), only one study to date has investigated the functional role of lncRNAs in transgenerational stress transmission (Gapp et al., 2020).

Our laboratory has reported previously that paternal exposure to chronic social defeat stress (CSDS) increases baseline anxiety-like and depression-like behaviors in male and female offspring and enhanced stress sensitivity in male offspring (Dietz et al., 2011). Although there has been extensive research establishing the distinct brain circuitry (Wook Koo et al., 2016) and transcriptome (Bagot et al., 2016, 2017; Krishnan et al., 2007; Lorsch et al., 2018; Nasca et al., 2019) of resilient versus susceptible animals, how CSDS exposure differentially affects the transmission of stress phenotypes to resilient and susceptible lineages is not well understood. Here, we demonstrate that key aspects of anxiety-like behaviors can be transmitted to the offspring of socially defeated fathers, dependent in part on whether the father is resilient or susceptible. Moreover, using artificial insemination (AI), we show that sperm partly mediates the transmission patterns seen in both male and female offspring.

To further understand the molecular mechanisms underlying the intergenerational transmission of these stress phenotypes, we used RNA-sequencing (RNA-seq) to examine the transcriptomic profiles of sperm genome-wide before and after exposure to CSDS. Before exposure to CSDS, we found subtle differences in the sperm of resilient versus susceptible males. By contrast, following CSDS there are robust differences in the sperm transcriptome, with limited overlap in the biological pathways disrupted by stress exposure in the two subgroups of animals. Interestingly, after CSDS in susceptible sperm, there was a dramatic increase in the proportion of differentially expressed lncRNAs. Co-expression gene network analysis predicted lncRNAs and certain protein-coding hub regulators of these stress-specific effects on gene expression patterns in sperm.

Overall, we present novel insight into how the paternal transmission of stress phenotypes occurs in part directly via sperm. These data suggest that persistent changes in the make-up of RNAs in sperm that are transmitted to the offspring before versus after stress in stress-susceptible fathers could relay different information to offspring. Our data also suggest that lncRNAs, heretofore mostly unexplored in this context, may be an important epigenetic mechanism contributing to the intergenerational transmission of stress phenotypes.

Materials and Methods

Animals

All experimental C57BL/6J male and female mice were obtained from The Jackson Laboratory. Retired CD1 breeders used as the aggressors for the social defeat paradigm were obtained from Charles River. All mice were maintained in a temperature-controlled and humidity-controlled facility on a 12/12 h light/dark cycle (lights on at 7 A.M.) with food and water ad libitum. All procedures were conducted in accordance with the Institutional Animal Care and Use Committee guidelines of the Icahn School of Medicine at Mount Sinai.

CSDS and social interaction (SI) testing

CSDS was performed as published previously (Berton et al., 2006; Krishnan et al., 2007; Lorsch et al., 2018; Peña et al., 2019a). Briefly, F0 male C57BL/6J experimental mice were exposed to a novel aggressive CD1 male mouse for 10 min/d for 10 d. After each exposure the mice were separated for the remainder of the day by a Plexiglas barrier that allows for sensory contact without physical interaction. Non-defeated controls were run in parallel housed two animals per cage under the same conditions as their experimental counterparts but in the presence of another control mouse rather than an aggressive CD1 mouse. Twenty-four hours after the final tenth day of defeat (experimental day 11), initial SI testing was performed under red light as previously described (Bagot et al., 2016, 2017) to assess social avoidance. F0 mice were placed into an open arena (44 × 44 cm) with an empty wire cage at one side (interaction zone). Mice were given 2.5 min to explore the arena and then removed. A novel CD1 aggressor was placed in the small wire cage and the procedure was repeated for an additional 2.5 min. Time in the interaction zone was recorded automatically with video tracking software (Ethovision 10.0, Noldus). Data were analyzed as time spent in the interaction zone when the aggressor was absent (no target) compared with time spent in the interaction zone when the aggressor was present (target) to obtain a SI ratio. The segregation of defeated mice into susceptible and resilient subpopulations was performed as previously described, work which has demonstrated SI testing as a highly reliable and consistent way to differentiate susceptible versus resilient mice after CSDS (Krishnan et al., 2007; Wilkinson et al., 2009; Vialou et al., 2010). Defeated mice with SI ratio scores >1.1 were categorized as resilient. Defeated mice with SI ratio scores <0.9 were categorized as susceptible. Defeated and non-defeated males were singly housed for 28 d to allow sperm to mature and put through a second SI test (experimental day 39) to confirm their categories, as singly housing males can be stressful (Ieraci et al., 2016).

Natural mating (NM) and Artificial Insemination (AI)

On experimental day 40 (30 d after the last exposure to CSDS), F0 males were allowed to naturally mate with a naive female to produce a NM litter. Approximately 4–7 d later (experimental days 44–50), sperm from the same F0s males was used for AI of control female. We performed AI as outlined by (Stone et al., 2015). Ovulation was induced in naive females by injection with pregnant mare's serum gonadotropin (1 IU) intraperitoneal in the evening. Two days later, females were injected with human chorionic gonadotropin (1 IU) intraperitoneal and 2 d later sperm transfer was performed. Briefly, males were euthanized by cervical dislocation 1 h before the planned sperm transfer. The vas deferens and cauda epididymis were dissected and placed into 500 µl of human tubal fluid (HTF) medium (with 4 mg/ml bovine serum albumin) at 37°C in an in vitro fertilization (IVF) dish. The cauda was dissected to free the sperm. The sperm cells were allowed to swim for a few minutes, all tissue was removed, and the sperm was incubated in a CO2 incubator for ∼30 min to capacitate (the process that enables sperm to fertilize an egg). A total volume of 40 µl of sperm was delivered to the uterine horn of two females to increase the likelihood of litter production using the non-surgical embryo transfer procedure. Females were immediately paired with a vasectomized (VASEX) male overnight and copulation was confirmed by the presence of a copulation plug. After 2–3 d, the female and VASEX male were separated. VASEX males were rotated across all paternal categories to confirm that behavior of the VASEX male had no effect on maternal behavior. On embryonic day (E)12, pregnancy check was confirmed via weight gain using an 8- to 12-g increase as an indication of litter expectancy.

Offspring stress

One male and one female from each litter was exposed to stress and examined on a battery of behavioral paradigms to examine the effects of paternal exposure to CSDS on offspring stress sensitivity. The remainder of the litter was used to examine offspring behaviors at baseline (i.e., unstressed).

Subchronic variable stress (SCVS) was used to assess susceptibility to stress in male and female offspring. SCVS was performed as described previously (Bodnoff et al., 1989; Hodes et al., 2015), which consists of three different stressors with one stress given every day over 6 d. The stressors are alternated during the 6-d period to prevent stress habituation. Stressors were administered in the following order: 100 random mild foot shocks at 0.45 mA for 1 h, a tail suspension stress for 1 h and restraint stress, placed inside a 50-ml falcon tube, for 1 h within the home cage. The three stressors were given two times in the same. On day 7, F1 offspring were started on the behavioral paradigm outlined.

Modified CSDS (mCSDS) was used to assess susceptibility to a more robust stress in male offspring produced by NM and AI. Male offspring previously examined in behavioral paradigms at baseline were exposed to mCSDS. mCSDS was conducted in the same manner as standard CSDS, only performed daily for 8 d rather than the standard 10 d. Specific methods on standard CSDS can be found above.

Offspring behavioral testing

Male and female F1 offspring were examined at postnatal day (P)60 in a battery of behavioral tests in the following order: novelty-suppressed feeding (NSF), open field (OF), and SI with animals put through one behavioral test per day. Offspring were examined at baseline and following exposure to SCVS. Offspring examined at baseline were run in parallel with stressed animals. All behavioral testing was performed 1 d apart and conducted between 10 A.M. and 6 P.M.

NSF was modified as described (Peña al., 2019b). Briefly, NSF was performed after 24 h of food deprivation. On the day of testing, mice habituated to the testing room for at least 1 h. Mice were placed in a novel arena with corncob bedding and a single piece of food in the center of the arena. Time to feed was recorded manually under red light. Mice were given a maximum of 10 min to eat, after which the trial was ended, and latency of 600 s was recorded. After the mouse ate in the novel arena, the mouse was returned to the home cage where a single piece of food was located in the center, and time to eat in the home cage was recorded. Data were analyzed as latency to eat in the novel arena and latency to eat in the home cage.

Exploration of an OF arena (44 × 44 cm) was assessed during a 5-min test. A video tracking software (Ethovision 10.0, Noldus) was used to measured locomotor activity, as well as the time spent in the middle and periphery of the test arena as an index of anxiety.

SI test was performed in a similar fashion as described above on all offspring. When SI was measured in male and female animals, a same-sex C57Bl6/J mouse was used as a measure of social novelty/social avoidance. In the subset of male offspring that were exposed to mCSDS, a novel CD1 aggressor was used to measure SI after defeat exposure.

Statistics for behavior

All statistics were conducted in SPSS (version 25) or Prism GraphPad (version 8.3.1.). Outliers were removed using a Grubbs outlier test. When performing ANOVAs, if the residuals were not normally distributed, the data were transformed until a normal distribution could be achieved; where data were normally distributed (Kolmogorov–Smirnov test for normality of residuals >0.10) one-way or two-way ANOVA was used. If transformation was not successful in generating a normal distribution, a Kruskall–Wallis non-parametric test was used. If an interaction was identified, a Bonferroni post hoc analysis was conducted to determine specific differences between groups. If an effect was identified via a Kruskall–Wallis test, follow-up Dunn's multiple comparisons test were run to determine specific differences between groups. χ2 tests were used to determine differences in the distribution of RNA-seq gene biotypes. Statistical significance was set to a p ≤ 0.05 and statistically trending was set to a p value between 0.05 and 0.10.

Sperm collection and CSDS for RNA-seq

Mature sperm were collected from surgically dissected left cauda epididymis from 8- to 10-week-old test mice on experimental day 1. Under ketamine (100 mg/kg, i.p.) and xylazine (10 mg/kg, i.p.) anesthesia a 1 cm vertical incision was made just lateral to the spine, and following blunt dissection, the left testis was identified. A suture was placed just below the vas deferens, and the cauda epididymis was removed with scissors and placed in HTF medium in a cell culture plate, punctured with a needle, and the sperm cells were allowed to swim for 10–15 min. Sperm suspension and tissue were transferred to a 1.5-ml microcentrifuge tube and incubated for 30 min in a 37°C heater. The sperm suspension was pelleted by centrifugation at 10,000 rpm for 10 min at 4°C and the supernatant was pipetted off. The supernatant was examined under a microscope to confirm no sperm cells were found. All sperm samples were snap frozen on dry ice and stored at −80°C until use. The incision was closed using sterile sutures and the animals were allowed to recover from surgery for 7 d. On experimental day 8, the same surgically manipulated animals were put through a 10-d CSDS paradigm with non-surgery control animals run in parallel (detailed methods for CSDS above); 28 d after exposure to CSDS (experimental day 46) animals were put through a second SI test to confirm their SI category. One day later (experimental day 47), animals were euthanized via cervical dislocation and mature sperm was collected from the right cauda epididymis, snap frozen on dry ice, and stored at −80°C until use.

RNA extraction and RNA-seq library preparation

Total RNA was isolated from frozen mature sperm from adult male mice using a protocol combining TRIzol/chloroform extraction with RNeasy MiniElute Cleanup kit. In short, 200-µl TRIzol was used to homogenize the tissue. An additional 500-µl TRIzol was added along with 140-µl chloroform. The lysate was centrifuged at 12,000 × g for 15 min and the aqueous layer was removed and added to the RNeasy mini spin column. The rest of the protocol was completed as outlined by the manufacturer. Samples were treated with DNase to rid samples of genomic DNA. After extraction, purity and concentration were measured on an Agilent 2100 Bioanalyzer (Agilent) using an RNA 6000 Nano assay. To account for batch effects libraries were made in three batches with all batches containing equal samples from each group. Libraries were prepared using the TruSeq Stranded Total RNA Library Prep Gold protocol (Illumina). Briefly, cDNA was synthesized from ribo-depleted RNA, followed by end-repair and ligation with sequencing adaptors. The libraries were then size selected and purified using AMPure XP beads (Beckman Coulter). Bar code bases (8 bp) were introduced at one end of the adaptors during PCR amplification steps. Library size and concentration were measured by Bioanalyzer or Tape Station (Life Technologies) before sequencing. All libraries were pooled for multiplexing and sequenced on the Illumina HiSeq 2500 System using V3 chemistry with 150 base pair double-end reads at the sequenced on a HighSeq2500 System using V4 chemistry with 50 base pair single-end reads at GeneWiz LLC. QC revealed an average of 40 million reads per sample (minimum = 20 million; maximum = 50 million) with an average mapping rate of 71.3%. The number of independent tissue samples included in the final analysis was between six and eight per group.

RNA-seq differential expression analysis

Data were aligned to mm10 genome and filtered for low abundance transcripts by keeping only genes expressed at CPM ≥5 in at least 85% of samples before further analysis. After filtering, pair-wise differential expression comparisons using DEseq2 were performed to identify differentially expressed genes (DEGs) in comparisons of preresilient and postresilient versus control and susceptible versus control and a nominal significance threshold of fold change (FC) > 2 and p < 0.025 was applied. Functional annotation for gene ontology (GO) of biological processes and molecular functions was performed using EnrichR with gene identities of DEGs (Chen et al., 2013) with an adjusted p ≤ 0.05 and more than or equal to five genes reporting only non-redundant categories. A paired analysis was conducted to compare predefeat and postdefeat susceptible mice. To account for RNA-seq samples obtained from the same mice at predefeat and postdefeat time points, mouse ID was included as a covariate to the DESeq2 model.

Upstream regulator and pathway analysis

Predicted upstream regulators and molecular pathways were identified using Ingenuity Pathway Analysis Software (QIAGEN). These determinations were based on the log FC of DEGs in resilient or susceptible animals compared with controls.

Rank-rank hypergeometric overlap (RRHO) RNA-seq analysis

Threshold-free reads per kilobase million (RPKM) lists were ranked by their -log10(RPKM) values. RRHO was then applied to these gene expression values in the sperm of control animals predefeat versus postdefeat to examine the effect of time on the dynamic transcriptome of sperm.

Full threshold-free RNA expression lists were ranked by the −log10(p value) multiplied by the sign of the FC from the DESeq2 analysis and filtered to all genes with CPM ≥ 5 in at least 85% of samples. RRHO was used to evaluate the overall patterns of gene expression in resilient and susceptible sperm before and after defeat. To evaluate the overlap of gene expression between resilient and susceptible sperm before and after defeat, threshold-free genome-wide transcriptomic overlap analysis was conducted using RRHO (RRHO2; Seney et al., 2018). Briefly, the RRHO2 ranks gene lists by a signed p value, that is, the –log10(p value) multiplied by the sign of the FC, and generates matrices of overlapped genes between RNA-seq lists of interest independent of whether individual RNAs are differentially expressed or not.

lncRNA correlation analysis

Data filtered for low abundance transcripts (only genes expressed at CPM ≥ 5 in at least 85% of samples) were used to conduct correlation analysis between lncRNAs and protein coding genes (PCGs). Only samples from either susceptible or control groups and from the post defeat time point were used to calculate the RPKM values. The cor function was used to compute the r (significant: r≥±0.9) and p values (p ≤0.05) between PCGs and lncRNAs, with a correlation between genes indicated by both of the qualifications being met.

Weighted gene co-expression analysis

All animals were analyzed together through weighted gene co-expression network analysis (WGCNA) to examine the network structure of the expressed transcriptome following stress.

Using the Cytoscape software platform, we conducted a comprehensive analysis of the relationship between nodes using the maximal clique centrality (MCC) function, which has reported to be the most effective method of identifying hub nodes in a co-expression network (Chin et al., 2014).

Enrichment analysis

Fisher's exact tests (FETs) were conducted using the super exact test package in R as described previously to determine module preservation as well as enrichment of patterns and other gene lists (Wang et al., 2015)

Results

Sperm mediates a subset of behavioral phenotypes in offspring from stressed fathers

Our previous study found that offspring of CSDS-exposed fathers generated by IVF failed to show most of the phenotypes observed with NM (Dietz et al., 2011). However, it is known that dramatic epigenetic reprogramming can occur in the context of IVF, perhaps confounding the interpretation of this experiment. Therefore, to more accurately test whether sperm per se can be the active vehicle for intergenerational transmission of stress phenotypes, we used AI to investigate how sperm directly contributes to altered baseline and stress-induced phenotypes seen in litters generated from susceptible, resilient, and control fathers. Founder generation (F0) males were subjected to 10 d of CSDS followed by SI testing on day 11 to categorize them as either resilient (R) or susceptible (S; Krishnan et al., 2007; Bagot et al., 2015). Mice were singly housed for 28 d after SI testing, and SI phenotype (R or S) was again confirmed on day 39 with a matched group of controls run in parallel (Fig. 1B). These results show the consistency of the susceptible and resilient phenotypes over time after CSDS as has been demonstrated previously (Krishnan et al., 2007; Wilkinson et al., 2009). Mice were then bred with naive females the next day (at ∼13.5 weeks old; experimental day 40), and 4 d later (experimental days 44–50) sperm from the same mice was used to impregnate hormonally primed females of the same strain (Fig. 1A). Anxiety-like phenotypes were examined in offspring produced by both NM (n = 8−11 offspring/group; 7–11 litters/group) and AI (n = 8−11 offspring/group; 6–10 litters/group) using NSF (Reeb et al., 1992), SI testing, and OF (Christmas and Maxwell, 1970; Prut and Belzung, 2003; n = 7–12 offspring/group) at baseline and after subthreshold chronic variable stress (SCVS; Fig. 1C; NM: n = 7−12 offspring/group; AI: n = 7−10 offspring/group for each sex).

Figure 1.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 1.

Paternal exposure to CSDS alters stress sensitivity of male offspring in adulthood. A, Schematic timeline of paternal stress, mating, and offspring birth. B, Within animal comparison of F0 males used in NM and AI before (day 11) and after (day 39) CSDS in SI testing. C, Schematic of male offspring stress and behavior paradigms. D, NSF ratio was affected by offspring stress in F1 sons produced by NM (F(1,54) = 16.32, p = 0.0002) or AI (F(1,48) = 4.072, p = 0.0492) but not paternal stress (p > 0.05). Bonferroni's post hoc test revealed sons from susceptible lineages generated by NM following exposure to SVCS showed an increased NSF ratio when compared with non-stressed littermates (t(1,54) = 4.004, p = 0.0006). E, Offspring stress affected the percent time spent in the middle of an OF arena (F(1,55) = 13.23, p = 0.0006) in male offspring generated by NM. Bonferroni's post hoc test revealed F1 sons from control (t(1,55) = 2.597, p = 0.0361) and resilient (t(1,55) = 2.700, p = 0.0276) lineages produced by NM spent less time in the middle of the arena when compared with unstress littermates from the same paternal lineage. Neither paternal nor offspring stress had an effect on behaviors in OF test in AI sons. F, Schematic of male offspring CSDS. G, In sons produced by NM, there was a main effect of CSDS (F(1,24) = 19.75, p = 0.0002) in SI testing. Control (t(1,24) = 3.150, p = 0.013) and resilient (t(1,24) = 2.822, p = 0.0283) sons had a higher SI ratio after CSDS exposure, whereas susceptible sons did not show this effect. In sons produced by AI, offspring stress and paternal category affected SI ratio (interaction: F(1,19) = 7.328, p = 0.014). Similar to results with NM, AI resilient sons had a higher SI ratio (t(1,19) = 3.032, p = 0.0204) after CSDS, whereas susceptible sons did not. Error bars indicate mean ± SEM. Significant (p ≤0.05) two-way ANOVA interaction noted on graph: φ, main effect of offspring stress; δ, main effect of paternal category; σ, interaction of paternal category X offspring stress. For direct comparison of two groups: *p ≤ 0.05, **p ≤ 0.01, ***p ≤ 0.001, ****p ≤ 0.0001, #p ≤ 0.1. A table summary for all statistical test results can be found in Extended Data Figure 1-1 and Figure 1-2.

Extended Data Figure 1-1

Extended Data Figure 1. Download Figure 1-1, DOCX file.

Extended Data Figure 1-2

Extended Data Figure 1-2. Download Figure 1-2, DOCX file.

We hypothesized that, if sperm mediates aspects of the paternal transmission of stress phenotypes, we will see distinct behavioral phenotypes in offspring from stress lineages compared with controls and the patterns of behavior will be similar in NM and AI litters. Unstressed offspring from stressed fathers, whether S or R and whether generated by NM or AI, were indistinguishable from offspring from control (C) lineages in all behavioral tests. We exposed male offspring to a SCVS procedure, which typically does not produce anxiety-like phenotypes in males (Hodes et al., 2015). There was no effect of male offspring exposure to SCVS in SI testing (data not shown; Extended Data Fig. 1-1 and Fig. 1-2). However, after exposure to SCVS, NM sons from S lineages showed a higher NSF ratio compared with sons from control (C) lineages (t(54) = 4.004, p = 0.0006), indicative of increased anxiety-like behavior, an effect not apparent in AI sons (t(44) = 0.07265, p >0.9999; Fig. 1D). In the OF, there was a main effect of offspring stress (F(1,55) = 13.23, p = 0.0006), with both C and R sons from NM decreasing the percent time spent in the middle of the arena when compared with unstressed littermates (t(55) = 2.597, p = 0.0361 and t(55) = 2.700, p = 0.0276, respectively). By contrast, NM sons from S lineages showed no change in the percent time spent in the middle of the arena between unstressed and stressed offspring, suggesting a blunted stress response (t(55) = 0.9344, p > 0.9999; Fig. 1E). These data provide the first evidence that male offspring from stress lineages are sensitive to stressors besides subthreshold CSDS (Dietz et al., 2011). However, these phenotypes where not seen in sons produced by AI, which suggests that SCVS may not be a robust enough stressor to produce altered stress-induced anxiety-like behaviors in offspring produced by AI.

View this table:
  • View inline
  • View popup
Table 1.

GO term summary

View this table:
  • View inline
  • View popup

To examine whether a more robust stressor would induce altered stress-responses, we next studied F1 males produced by NM or AI in an SI test before and after exposure to a mCSDS. To do this, we took the unstressed offspring from the previous experiment and exposed them to an 8-d CSDS paradigm (mCSDS; Fig. 1F). After mCSDS, NM sons displayed a main effect of stress (F(1,24) = 19.75, p = 0.0002) with sons from C (t(24) = 3.15, p = 0.013) and R (t(24) = 2.822, p = 0.0283) lineages, but not S lineage (t(24) = 1.789, p = 0.2585), showing an increase in SI ratio. Importantly, the same distinct phenotypes were seen in R and S sons produced by AI, indicating that sperm directly contributes to the intergenerational transmission of this stress phenotype (t(19) = 3.032, p = 0.0204 and t(19) = 0.1693, p = 0.9977, respectively; Fig. 1G).

There are a limited number of intergenerational stress studies that have examined how paternal exposure to stress affects female offspring. Given that there are dramatic sex differences in stress responses between males and females (LaPlant et al., 2009; Maeng and Milad, 2015; Labonté et al., 2017), we investigated potential sex differences in stress responses in offspring from stress lineages. We examined the phenotypes of female offspring from stress lineages at baseline and after exposure to SCVS (Fig. 2A). In female offspring produced by NM, at baseline daughters from S lineages showed an increased SI ratio compared with C, suggesting increased social investigation (t(49) = 3.388, p = 0.0028; Fig. 2B). This effect was lost in AI daughters after exposure to SCVS (t(49) = 0.1623, p > 0.9999). In OF, daughters produced by NM at baseline displayed a main effect of offspring stress (F(1,54) = 5.323, p = 0.0249) and post hoc analysis revealed that S daughters were more sensitive to stress with increased anxiety-like behavior, as measured by a decreased percent time spent in the middle of the arena compared with unstressed littermates (t(54) = 2.612, p = 0.0349; Fig. 2C). The same phenotype was seen in offspring produced by AI, demonstrating that sperm directly contributes to aspects of the intergeneration transmission of these behavioral phenotypes seen in daughters of S but not R lineages (t(53) = 2.067, p = 0.0582). No effect of paternal category was seen in NSF at baseline or after SCVS in NM or AI daughters (data not shown; F(2,44) = 0.5224, p = 0.5967 and F(2,54) = 0.04,218, p = 0.9587, respectively; Extended Data Fig. 2-1 and 2-2).

Figure 2.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 2.

AI with sperm from stressed fathers increases stress sensitivity of female offspring in adulthood. A, Schematic of female offspring stress and behavior paradigm. Stressed female offspring were exposed to a 6-d subthreshold chronic variable stress. B, Neither paternal category nor offspring stress affected SI ratio in NM F1 daughters. In F1 daughters produced by AI, there was a main effect of offspring stress (F(1,49) = 6.168, p = 0.0165). AI F1 daughters from susceptible lineages showed an increased SI ratio at baseline compared with unstressed controls (t(1,49) = 3.388, p = 0.0028), an effect reversed by stress exposure (t(1,49) = 2.885, p = 0.0174). C, Offspring stress affected the percent time spent in the middle of an OF arena (F(1,118) = 6.854, p = 0.01) in female offspring. F1 stressed daughters, generated by NM or by AI, from susceptible lineages spent less time in the center of the OF arena when compared with unstressed female littermates (F(1,54) = 9.376, p = 0.0034). Error bars indicate mean ± SEM. Significant (p ≤0.05) two-way ANOVA interaction noted on graph: φ, main effect of offspring stress; δ, main effect of paternal category; σ, interaction of paternal category X offspring stress. For direct comparison of two groups: *p ≤ 0.05, **p ≤ 0.01, ***p ≤ 0.001, ****p ≤ 0.0001, #p ≤ 0.1. A table summary for all statistical test results can be found in Extended Data Figure 2-1 and Figure 2-2.

Extended Data Figure 2-1

Extended Data Figure 2-1. Download Figure 2-1, DOCX file.

Extended Data Figure 2-2

Extended Data Figure 2-2. Download Figure 2-2, DOCX file.

Distinct transcriptional patterns in sperm in resilient and susceptible mice

Having established that sperm can serve as the active vehicle to transmit certain effects of paternal stress to male and female offspring, we next used RNA-seq to investigate stress-induced changes in the transcriptome of sperm that may contribute to the behavioral phenotypes seen in F1 offspring produced by AI. Given that our group has previously shown that mating male mice before their exposure to CSDS does not result in altered phenotypes in offspring (Dietz et al., 2011), we sought to examine how the sperm transcriptome changes in response to stress. First, mature sperm was surgically collected from the left caudal epididymis of male mice before CSDS (Fig. 4A). The same surgically-manipulated males were put through 10 d of CSDS and 1 d later an SI test was used to categorize them as either R or S, with C F0 animals run in parallel to account for differences in age (Fig. 3A). The mice were then singly housed for 28 d and put through another SI test (experimental day 46) and 1 d later (experimental day 47; 30 d after the last defeat exposure) mice were euthanized, and mature sperm was collected from the right caudal epididymis of C, R, and S mice. Six to seven samples per group were used for RNA-seq analysis.

Figure 3.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 3.

The effect of time on the transcriptome of sperm. A, Scheme of paradigm and sperm collection time points in control animals. Note that, while the early time point is denoted “predefeat” and the later time point is denoted “postdefeat,” all animals in this experiment were controls and not exposed to CSDS. B, Threshold-free comparison of gene expression by RRHO. Pixels represent the overlap between the transcriptome of each comparison as noted, with the significance of overlap (−log10(RPKM) of a hypergeometric test) color coded. C, Heatmaps show all predefeat and postdefeat RNAs ranked by RPKM values seeded by predefeat values. The third row of RNAs indicates DEGs that are expressed at higher (yellow) or lower (blue) levels at the later versus earlier time point.

Figure 4.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 4.

Distinct transcriptional regulation in sperm of resilient and susceptible mice. A, Schematic timeline of F0 stress and sperm collection. B, Within animal comparison of F0 males used in NM and AI before (day 18) and after (day 46) CSDS in SI testing. C, Overlap of the gene list from these sequencing datasets with published gene lists from RNA-seq of epididymis versus sperm of mice and rats shows a higher overlap of genes from our gene list with the published sperm gene list than the epididymis gene list. Moreover, only ∼8% of the genes in our sperm dataset were not seen in the published sperm dataset, whereas ∼27% of our sperm genes were not seen in the published epididymis dataset. D, Threshold-free comparison of gene by RRHO. Pixels represent the overlap between the transcriptome of each comparison as noted, with the significance of overlap (−log10(p value) of a hypergeometric test) color coded. Lower left quadrant includes co-upregulated genes, upper right quadrant includes co-downregulated genes, and upper left and lower right quadrants include oppositely regulated genes (up-down and down-up, respectively). Genes along each axis are sorted in the same order. Note the dramatic overlap of transcriptional patterns between resilient versus susceptible mice predefeat and a virtually complete loss of such overlap postdefeat. E, Heatmaps show predefeat DEGs with at least a nominal p≤0.025 and a log2 FC ≥ 2. Heatmaps are seeded by either resilient or susceptible DEGs. F, Top three biological function GO terms enriched in resilient (orange) or susceptible (pink) versus control. G, Heatmaps show DEGs with at least a nominal p ≤ 0.025 and a log2 FC≥ 2. Heatmaps are seeded by either resilient (top) or susceptible (bottom) DEGs. H, Top biological function GO terms enriched in susceptible (pink) versus control. I, Upstream regulator analysis was conducted on both resilient and susceptible sperm DEGs and comparison analysis was used to examine the overlap of regulators. Canonical pathway analysis was conducted on both resilient and susceptible sperm DEGs and comparison analysis was used to look at the overlap. Activation Z scores in heatmaps: positive (yellow) = overrepresentation of targets activated by regulator; negative (blue) = overrepresentation of targets repressed by regulator; no direction (black) = no significant enrichment of activated versus repressed targets; white = not a predicted upstream regulator; * = statistically significant enrichment (p ≤ 0.05). A list of genes DEGs for the various comparisons described can be found in Extended Data Figure 4-1.

Extended Data Figure 4-1

These tables list DEGs for the various comparisons described in this study. Download Figure 4-1, XLSX file.

Surgery did not affect the proportion of animals that were later categorized as R or S (χ2 = 4.524, df = 6; p = 0.6061). To ensure the specificity of our RNA-seq analysis for sperm rather than epithelial cells potentially collected during extraction, we compared our gene list to published RNA-seq gene lists of the epididymis (Song et al., 2019) and sperm (Soumillon et al., 2013). We found that there was more overlap of our gene list with the published sperm sequencing gene list (∼91%) than with the epididymis gene list (∼72%; Fig. 4C). Additionally, to ensure that our control animals showed typical behavior in SI tests, we examined their SI ratios and, as in Figure 1B, again found no difference between SI tests 1 and 2 (p = 0.6061 and p = 0.3323, respectively; Fig. 4B).

Previous studies have found that small non-coding RNAs in sperm are dynamic and impacted by age in both mice (Chan et al., 2020) and humans (Morgan et al., 2020). However, no studies to date have examined such age-dependent patterns in sperm PCGs and lncRNAs. We used RPKM values as a measure of gene expression in the sperm of surgery-control predefeat and postdefeat mice to perform a within-animal comparison to examine the impact of time on the sperm transcriptome. Using rank-rank hypergeometric overlay (RRHO) to compare patterns and strength of genome-wide overlap in a threshold-free manner (Cahill et al., 2018), we observed robustly overlapping patterns of gene expression in predefeat versus postdefeat control animals (Fig. 3B). However, when we directly examined DEGs in postdefeat compared with predefeat control animals, we found 929 genes that were significantly different between the two groups (Fig. 3C). This suggests that, although the sperm transcriptome is mostly the same across age, a subset of transcripts show altered expression as a function of time, thus mirroring previous results reported for small non-coding RNAs.

In predefeat sperm, we expected to see minimal transcriptomic differences between R and S sperm, since all animals examined are essentially normal C57BL/6J mice without any perturbations. We again used RRHO to compare patterns and strength of genome-wide overlap in a threshold-free manner. RRHO analyses showed that, before defeat, there is strong co-regulation of genes in the same direction in R and S mice compared with C mice. CSDS disrupts this pattern with opposite regulation observed for many of the same genes (Fig. 4D). The latter finding indicates that CSDS produces a striking divergence in transcriptional regulation in sperm between R versus S mice.

To identify DEGs in S and R sperm both predefeat and postdefeat, the significance of DEGs was set at a nominal p < 0.025 and a log2 FC≥2. Predefeat there were 873 R-specific DEGs and 400 S-specific DEGs, with PCGs making up the majority of the biotypes (R: ∼92%, S: ∼94%) followed distantly by lncRNAs (R: ∼5%, S: ∼6%) and pseudogenes (R: ∼3%, S: ∼0.25%; Fig. 5A; Extended Data Fig. 4-1). This distribution of RNA biotypes is similar to what has been reported previously for sperm (Soumillon et al., 2013). Predefeat, heatmaps seeded by either R or S DEGs revealed that both R-specific and S-specific DEGs tend to be similarly dysregulated in the other group (Fig. 4E), confirming our RRHO analysis. To gain insight into the molecular and biological functions of these DEGs, we used Enrichr GO analysis. Many of the top GO terms in R sperm (∼66%) are related to cell cycle processes, whereas all top S sperm GO terms are related to cell migration (Fig. 4F; Table 1).

View this table:
  • View inline
  • View popup
Figure 5.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 5.

CSDS alters the biotype composition of RNAs in susceptible sperm. A, Pie chart represents the biotypes making up resilient and susceptible DEGs in sperm before and after CSDS in the same animals. There is a difference in the proportion of genes in each biotype composing DEGs in susceptible (χ2 = 48.97, df = 2; p < 0.0001) but not resilient (χ2 = 2.494, df = 2; p = 0.2874) sperm following CSDS. Following CSDS the proportion of lncRNAs making up susceptible DEGs increases (χ2 = 37.61, df = 2; p < 0.0001). B, Stacked bar graphs represent the proportion of specific lncRNA biotypes making up DEGs in resilient and susceptible sperm predefeat and postdefeat. C, Heatmap showing susceptible lncRNA DEGs with at least a nominal p < 0.025 and a log2 FC≥ 2 seeded by susceptible DEGs. Distribution plot of the number of highly correlated PCGs with lncRNAs in sperm of (D) control, (E) resilient, and (F) susceptible animals. A table with the summary of the statistics from the comparison of DEGs biotypes can be found in Extended Data Figure 5-1.

Extended Data Figure 5-1

Extended Data Figure 5-1. Download Figure 5-1, DOCX file.

Before turning to DEG analysis of sperm after CSDS, we examined the effect of surgery per se on the sperm transcriptome by performing RNA-seq on sperm from surgically-manipulated control mice compared with non-surgically-manipulated control mice. There were 1527 surgery-induced DEGs in sperm with biological functions related to flagellum function, ncRNA metabolic processes, and spermatogenesis (Table 2). In subsequent analyses, surgery DEGs were filtered out to examine stress-induced transcriptional differences in the sperm of R and S animals with the cut off for DEGs set at a nominal p < 0.025 and a log2 FC ≥ 2.

Postdefeat, differential expression analysis revealed only 62 R DEGs, with PCGs making up the large majority of the biotype (∼93%) followed by lncRNAs (∼7%; Fig. 5A). This finding illustrates a very limited effect of CSDS on the sperm transcriptome of R fathers. By contrast, there were 1460 S DEGs, with PCGs making up a smaller majority of the biotype (∼79%) followed by lncRNAs (∼18%) and pseudogenes (∼3%; Fig. 5A; Extended Data Fig. 4-1). Heatmaps of the R DEGs revealed that many of these genes are similarly regulated in S animals (Fig. 4G). However, heatmaps of S DEGs revealed that the vast majority appear unchanged in R animals (Fig. 4G).

GO analysis revealed no enrichment of postdefeat R DEGs for biological functions presumably because of the very small number of DEGs. The majority of the top biological function GO terms for S sperm (∼87%) were related to flagellum function (Fig. 4H; Table 2), with one term involved in cell projection functioning. Dysregulation of flagellum processes can impact fertility (Lehti and Sironen, 2017); however, there was no difference in litter size or proportion of male or female offspring between control and stress lineages (F(2,19) = 0.2637, p = 0.7710). This finding suggests that fertility may not be affected dramatically, but rather that genes involved in flagellum function might relay paternal environmental information to offspring.

We next examined predicted upstream regulators of stress-induced transcriptomic changes in sperm of R and S mice. In R mice, there were five predicted upstream regulators with predicted inactivation [estrogen receptor 1 (ESR1), nuclear factor erythroid 2 like 2 (NFE2L2), interleukin 1-family (IL1), IL1-β, and CCAAT enhancer binding protein α (CEBPA)] and one molecule with predicted activation [Kruppel-like factor 3 (KLF3); Fig. 4I]. There were three overlapping predicted upstream regulators between predefeat and postdefeat in R mice (ESR1, NFE2L2, and KLF3), however, there was no change in the activation state of these molecules following stress.

In S sperm, there were seven molecules predicted as upstream regulators: vascular endothelial growth factor (VEGF), Rab-like protein 6 (RABL6), stratifin, aryl hydrocarbon receptor (AHR), estrogen receptor (ER), prostaglandin ER2 (PTGER2), and tumor protein p53 (Fig. 4I). There were four overlapping predicted upstream regulators in predefeat and postdefeat S sperm (VEGF, SHR, ER, and TP53). Notably, only one molecule, VEGF, went from an activated to inhibited state after CSDS. VEGF is particularly interesting as a previous intergenerational study performed by our laboratory found that plasma levels of VEGF were lower in offspring from fathers after exposure to CSDS (Dietz et al., 2011), and numerous other studies have linked VEGF to depression in both rodent (Deyama et al., 2019) and human studies (Clark-Raymond and Halaris, 2013). Of note, ESR1 is a major subtype of ER, which indicates that ER signaling was a predicted upstream regulator of stress effects in both R and S sperm.

Regulation of lncRNA expression in sperm of resilient and susceptible mice

Analysis of the distributions of RNA biotype among DEGs revealed a significant difference in predefeat versus postdefeat in S but not R sperm. Post hoc tests revealed that this was driven by an increase in the proportion of lncRNA DEGs in S animals only, with 277 lncRNA DEGs in S mice postdefeat compared with only four lncRNA DEGs in R mice postdefeat and 24 and 37 lncRNAs in S and R mice predefeat, respectively (Fig. 5A,B). Heatmaps of only the lncRNAs regulated in S mice postdefeat revealed that the vast majority are unaffected in R mice (Fig. 5C).

To further assess the impact of altered expression of lncRNAs in sperm from S fathers, we performed correlation analysis to relate the expression of lncRNAs with the expression of PCGs (Issler et al., 2020). lncRNAs and PCGs that show correlated expression are potentially involved in similar biological functions (Rinn and Chang, 2012). Histogram plots of the distribution of highly correlated PCGs (significant = r ≥±0.9 and p ≤ 0.05) in each group show distinct dysregulation in S sperm specifically, with a striking loss of correlation of lncRNAs with PCGs in S sperm but not in R or C sperm (Fig. 5D–F).

Co-expression analysis identifies gene modules in sperm of resilient and susceptible mice

We next used WGCNA (Langfelder and Horvath, 2008) to construct co-expression modules using the expression data of sperm after exposure to CSDS. WGCNA created a network that consisted of 25 discrete modules with between 33 and 5807 genes in each module (with each module assigned an arbitrary color name). We next identified gene modules differentially affected by CSDS in R and S mice by determining whether modules were enriched for DEGs under either condition. We found three modules (grey60, salmon, and medium orchid) that were enriched in R DEGs and five modules (royal blue, bisque4, dark magenta, dark orange, and dark orange2) that were enriched in S DEGs. Interestingly, no modules were enriched for both R and S DEGs, suggesting distinct regulation of genes in the sperm of resilient and susceptible animals only after stress (Fig. 6A).

Figure 6.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 6.

Identification of co-expression gene networks in sperm differentially affected by CSDS in susceptible versus resilient fathers. A, Circos plot for all 25 WGCNA modules. Module names (an arbitrary color) and sizes are indicated on the outside of the Circos plot. Enrichment for either resilient or susceptible DEGs or enrichment for lncRNA genes are indicated by the semicircle colors within each module, with increasing warm colors indicating increasing −log10(p value). Arachne plots show the co-expression networks of (B) the susceptible-enriched and lncRNA-enriched module dark magenta, and (C) the resilient-enriched module medium orchid. Node size indicates the strength of stress-specific hub genes with increased size indicating top hub rank. Pink nodes indicate differential expression for that gene in susceptible animals. Orange nodes indicate differential expression for that gene in susceptible animals. Black halos indicate that gene is a lncRNA. Edges indicate significant co-expression between two particular genes.

Using the Cytoscape software platform, we performed comprehensive analysis of the relationship between nodes using the MCC function, reported to be the most effective method of finding hub genes in a co-expression network (Chin et al., 2014). Given the loss of correlation between lncRNA and PCG DEGs in S sperm, we looked for modules enriched in lncRNAs and R or S DEGs. Of particular interest, the S-enriched module, dark magenta, was also enriched in lncRNAs and contained a lncRNA, Gm27211, as a top hub gene. Additionally, Gm27211 was differentially expressed in S animals after CSDS (Fig. 6B). Gm27211 is expressed from E10 to E17 in the central nervous system, and in male and female reproductive tracts from P4 to adulthood; however, the specific function of this lncRNA is largely unknown (Brunskill et al., 2014; Bao et al., 2016; Kistler et al., 2015; Huntley et al., 2016).

Interestingly there was no overlap in R-enriched and lncRNA-enriched modules. The R-enriched module, medium orchid, had two hub genes that were differentially expressed, sarcospan (Sspn) and Rho GTPase activating protein 21 (Arhgap21), both of which are PCGs involved in cell adhesion and actin dynamics, respectively (Parvatiyar et al., 2015, 2019; Fig. 6C).

Discussion

In the present study, we provide several new lines of evidence concerning the transmission of paternal stress to the next generation. We show that susceptible and resilient fathers, after a course of CSDS, transmit different degrees of stress susceptibility to their offspring, with differences in phenotypes seen between male versus female offspring. We go on to show that a subset of these effects of paternal stress occur in offspring generated by AI, demonstrating that sperm per se is an active vehicle in the intergenerational transmission of stress phenotypes. Finally, we characterize the transcriptome of sperm of susceptible, resilient, and control fathers by RNA-seq and identify prominent regulation of lncRNAs in susceptibility selectively, thus implicating these molecules in next-generation consequences of paternal stress.

Paternal transmission studies typically capitalize on the fact that most male rodents are not involved in offspring rearing, which allows observed phenotypes in offspring to be attributed to male gamete contribution. We previously reported only minimal intergenerational transmission of stress vulnerability when offspring were produced via IVF, suggesting that behavioral or hormonal changes rather than epigenetic mechanisms may explain the dramatic alterations in offspring phenotypes (Dietz et al., 2011). However, the lack of robust behavioral phenotypes observed in this earlier study may be because of the fact that IVF alters the epigenome of sperm and egg cells (Khosla et al., 2001; Reik and Walter, 2001; Calle et al., 2012; Ventura-Juncá et al., 2015), which may wash out stress-induced epigenetic changes. In the present study, therefore, we used AI, which removes the need for many of the IVF factors that are believed to influence the epigenome. We report now for the first time that male and female offspring produced by AI using sperm from CSDS-exposed males show altered stress sensitivity, thereby establishing that epigenetic alterations in sperm do indeed mediate some of the intergenerational transmission of altered stress phenotypes observed with NM.

More specifically, we found that paternal CSDS alters how sons produced by NM respond to SCVS but does not affect such behavioral responses in sons produced by AI. This suggests the involvement of other factors such as maternal care or hormonal signals during mating in transmitting this phenotype rather than epigenetic changes in sperm. However, we hypothesized that a more robust offspring stressor would reveal stress phenotypes transmitted by sperm. We therefore exposed male offspring to a mCSDS procedure. Sons produced by NM or by AI and exposed to mCSDS display equivalent differences when compared with sons of control mice, thus demonstrating the intergenerational transmission of stress sensitivity by sperm.

Female offspring from stress lineages produced by NM and AII showed a similar heightened response to stress in OF testing. Surprisingly, daughters from susceptible lineages produced by AI but not NM showed an increase in social investigation. Although it may be surprising that female offspring from S lineages produced by AI and not by NM showed increased social investigation at baseline, previous studies have reported altered maternal care when dames are allowed to mate directly with stressed males (Mashoodh et al., 2018). These studies suggest that maternal interaction with an affected father may lead to maternal masking and underscore the importance of using gamete manipulation techniques for disentangling how male gametes and epigenetic changes to them contribute to the transmission of stress phenotypes across generations overall hinting at potential maternal masking of altered phenotypes in our NM experiments.

Studies investigating potential molecular mechanisms that contribute to the transmission of stress phenotypes across multiple generations in rodent models have identified differences in microRNAs (Rodgers et al., 2013), lncRNAs (Gapp et al., 2020), and DNA methylation (Franklin et al., 2010). However, no study to date has identified how the RNA profiles in sperm are altered by paternal stress in a manner that distinguishes between those fathers that are susceptible versus resilient to stress exposure. Here, we provide an assessment of genome-wide transcriptomic changes in sperm before and after exposure to CSDS within the same F0 males who are susceptible or resilient to that stress.

A surprising finding of this experiment was the detection of differences in gene expression in the sperm of susceptible versus resilient mice before exposure to CSDS. This is surprising since, at this time point, the susceptible, resilient, and control mice represented essentially equivalent groups of wild-type C57BL/6J mice and we previously reported that mating males before stress exposure does not result in the transmission of stress phenotypes to offspring (Dietz et al., 2011). Indeed, RRHO analysis demonstrated broad convergence of gene expression patterns in the sperm across all three groups of animals, confirming very similar transcriptomes. Nevertheless, DEG analysis identified numerous genes that show different expression levels in mice that will become susceptible versus resilient when each was compared with the control group. Together, these analyses demonstrate that, within a group of inbred C57BL/6J with equivalent environmental experiences, there is a sizable degree of variation in sperm gene expression patterns, which correlates with the inherent vulnerability of an individual mouse to be either susceptible or resilient to social stress later in life. Understanding the factors responsible for such baseline variability in sperm gene expression is an interesting topic for future studies.

There is increasing evidence that stress alters the epigenome of sperm (Franklin et al., 2010; Gapp et al., 2014, 2020), however, there is little consistency across different paternal stressors and no studies have investigated whether similar changes are seen after CSDS. In contrast to the subtle predefeat differences in sperm gene expression, exposure to CSDS dramatically restructured the transcriptome of sperm from susceptible versus resilient fathers in very different ways. RRHO and DEG analyses revealed virtually no overlap in DEGs between the two subsets of fathers following exposure to CSDS. In susceptible fathers, there was a striking transcriptional reorganization of gene expression, with 1460 genes being differentially expressed. This contrasts dramatically with the resilient fathers, which show alterations in ∼20-fold fewer genes, demonstrating that resilient fathers escape the vast majority of the stress-induced changes in sperm gene expression seen in susceptible fathers. The finding that the sperm of inbred C57BL/6J mice display such remarkable differences in transcriptomic responses to the same stress underscores the tremendous potential for diversity even in the context of constant genetic background. Future studies are needed to understand the molecular mechanisms underlying these dramatic individual differences in stress responses in sperm.

Given that lncRNAs have been established as important for sperm functioning (Zhang et al., 2010, 2019; Wen et al., 2016), have been associated with depression GWAS studies (Roy et al., 2018; Wray et al., 2018), have been seen to be dysregulated in the sperm of diabetic mice (Jiang et al., 2016), and have also been identified as a partial mediator of the paternal transmission of stress phenotypes (Gapp et al., 2020), we examined this RNA biotype in sperm DEGs. We found a ∼10-fold increase in the number of stress induced differentially expressed lncRNAs in sperm from susceptible when compared with both resilient sperm after stress and prestress sperm independent of group. lncRNAs have been shown to serve as key drivers in cancer (Ashouri et al., 2016; Lanzós et al., 2017), and manipulation of lncRNAs in adult brain have been shown to induce depression-related behavioral phenotypes (Issler et al., 2020), suggesting that lncRNA expression can have broad impacts on a variety of phenotypes. To provide insight into the function of lncRNAs in our dataset, we used correlation analysis to relate lncRNA expression with PCG expression. We found similar patterns of lncRNA-PCG relationships in sperm of control and resilient fathers, but a dramatic loss of correlation in susceptible fathers. We next used WGCNA to characterize the connectivity of genes within our sequencing dataset and identified several gene modules that differentially enrich for either resilient or susceptible DEGs, with ∼80% of the susceptible-enriched modules also enriched for lncRNAs. These data are consistent with the view that lncRNAs contribute to the paternal transmission of stress phenotypes via regulation of PCGs. We looked for lncRNA hub genes in those modules and identified one, Gm27211, which now warrants further investigation. A key challenge for future research is to now relate altered levels of Gm27211 and presumably other regulatory genes in sperm to altered gene transcription early in development and eventually in the developing and adult brain, and ultimately to alterations in stress susceptibility observed in adult offspring.

lncRNAs like Gm27211 might also drive lasting phenotypes through interactions with microRNAs. It has been established that microRNAs regulate lncRNA expression in other tissues (da Rocha et al., 2008; Royo and Cavaillé, 2008), and studies suggest that lncRNAs may regulate miRNA function, creating the possibility of multiple feedback loops between long and short noncoding RNAs (Feng et al., 2006; Jeck and Sharpless, 2014; Ballantyne et al., 2016). This is consistent with previous work that has implicated certain microRNAs in sperm of stressed fathers as also contributing to the intergenerational transmission of stress susceptibility. However, efforts to understand the role of lncRNAs in paternal transmission of stress experience to their offspring are impaired by our still very limited knowledge of the biological function of the vast majority of lncRNAs (Bassett et al., 2014).

In summary, results from the present study show that male mice exposed to CSDS categorized as either susceptible or resilient transmit partly different behavioral responses to their offspring and that sperm is one active ingredient involved in this intergenerational transmission. Moreover, the behavioral changes observed in the offspring are accompanied by stress-induced disruption of the transcriptional state in sperm. This study thereby adds valuable insight into how stress interacts across generations and how the mechanisms of this transmission might differ between resilient and susceptible animals. This work thus suggests novel potential targets through which epigenetic mechanisms contribute to inherited phenotypes. At the same time, our findings raise two crucial questions for future analysis. The complex mechanisms underlying how behavioral experience induces long-lasting changes in sperm gene expression remain largely elusive. While epididymis extracellular vesicles have been shown to be one potential mechanism (Chan et al., 2020), many others may be at play but have not yet been investigated. As well, it will be essential to delineate the precise mechanisms by which changes in sperm gene expression lead to life-long changes in brain structure and function that mediate the altered behavioral phenotypes seen in adult offspring of stressed fathers.

Footnotes

  • This work was supported by National Instituties of Health Grants P50MH096890 and R01MH051399, R00 DA042100 and T32MH087004, and the Hope for Depression Research Foundation.

  • The authors declare no competing financial interests.

  • Correspondence should be addressed to Deena M. Walker at walkerde{at}ohsu.edu or Eric J. Nestler at eric.nestler{at}mssm.edu

SfN exclusive license.

References

  1. ↵
    1. Ashouri A,
    2. Sayin VI,
    3. Van den Eynden J,
    4. Singh SX,
    5. Papagiannakopoulos T,
    6. Larsson E
    (2016) Pan-cancer transcriptomic analysis associates long non-coding RNAs with key mutational driver events. Nat Commun 7:13197. doi:10.1038/ncomms13197 pmid:28959951
    OpenUrlCrossRefPubMed
  2. ↵
    1. Bagot RC,
    2. Parise EM,
    3. Peña CJ,
    4. Zhang HX,
    5. Maze I,
    6. Chaudhury D,
    7. Persaud B,
    8. Cachope R,
    9. Bolaños-Guzmán CA,
    10. Cheer JF,
    11. Cheer J,
    12. Deisseroth K,
    13. Han MH,
    14. Nestler EJ
    (2015) Ventral hippocampal afferents to the nucleus accumbens regulate susceptibility to depression. Nat Commun 6:7062. doi:10.1038/ncomms8062 pmid:25952660
    OpenUrlCrossRefPubMed
  3. ↵
    1. Bagot RC,
    2. Cates HM,
    3. Purushothaman I,
    4. Lorsch ZS,
    5. Walker DM,
    6. Wang J,
    7. Huang X,
    8. Schlüter OM,
    9. Maze I,
    10. Peña CJ,
    11. Heller EA,
    12. Issler O,
    13. Wang M,
    14. Song WM,
    15. Stein JL,
    16. Liu X,
    17. Doyle MA,
    18. Scobie KN,
    19. Sun HS,
    20. Neve RL, et al
    . (2016) Circuit-wide transcriptional profiling reveals brain region-specific gene networks regulating depression susceptibility. Neuron 90:969–983. doi:10.1016/j.neuron.2016.04.015 pmid:27181059
    OpenUrlCrossRefPubMed
  4. ↵
    1. Bagot RC,
    2. Cates HM,
    3. Purushothaman I,
    4. Vialou V,
    5. Heller EA,
    6. Yieh L,
    7. LaBonté B,
    8. Peña CJ,
    9. Shen L,
    10. Wittenberg GM,
    11. Nestler EJ
    (2017) Ketamine and imipramine reverse transcriptional signatures of susceptibility and induce resilience-specific gene expression profiles. Biol Psychiatry 81:285–295. doi:10.1016/j.biopsych.2016.06.012 pmid:27569543
    OpenUrlCrossRefPubMed
  5. ↵
    1. Bale TL
    (2015) Epigenetic and transgenerational reprogramming of brain development. Nat Rev Neurosci 16:332–344. doi:10.1038/nrn3818 pmid:25921815
    OpenUrlCrossRefPubMed
  6. ↵
    1. Ballantyne MD,
    2. McDonald RA,
    3. Baker AH
    (2016) lncRNA/MicroRNA interactions in the vasculature. Clin Pharmacol Ther 99:494–501. doi:10.1002/cpt.355 pmid:26910520
    OpenUrlCrossRefPubMed
  7. ↵
    1. Bao J,
    2. Vitting-Seerup K,
    3. Waage J,
    4. Tang C,
    5. Ge Y,
    6. Porse BT,
    7. Yan W
    (2016) UPF2-Dependent nonsense-mediated mRNA decay pathway is essential for spermatogenesis by selectively eliminating longer 3'UTR transcripts. PLoS Genet 12:e1005863. doi:10.1371/journal.pgen.1005863 pmid:27149259
    OpenUrlCrossRefPubMed
  8. ↵
    1. Bassett AR,
    2. Akhtar A,
    3. Barlow DP,
    4. Bird AP,
    5. Brockdorff N,
    6. Duboule D,
    7. Ephrussi A,
    8. Ferguson-Smith AC,
    9. Gingeras TR,
    10. Haerty W,
    11. Higgs DR,
    12. Miska EA,
    13. Ponting CP
    (2014) Considerations when investigating lncRNA function in vivo. Elife 3:e03058. doi:10.7554/eLife.03058 pmid:25124674
    OpenUrlCrossRefPubMed
  9. ↵
    1. Berton O,
    2. McClung CA,
    3. Dileone RJ,
    4. Krishnan V,
    5. Renthal W,
    6. Russo SJ,
    7. Graham D,
    8. Tsankova NM,
    9. Bolanos CA,
    10. Rios M,
    11. Monteggia LM,
    12. Self DW,
    13. Nestler EJ
    (2006) Essential role of BDNF in the mesolimbic dopamine pathway in social defeat stress. Science 311:864–868. doi:10.1126/science.1120972 pmid:16469931
    OpenUrlAbstract/FREE Full Text
  10. ↵
    1. Bodnoff SR,
    2. Suranyi-Cadotte B,
    3. Quirion R,
    4. Meaney MJ
    (1989) A comparison of the effects of diazepam versus several typical and atypical anti-depressant drugs in an animal model of anxiety. Psychopharmacology (Berl) 97:277–279. doi:10.1007/BF00442264 pmid:2567028
    OpenUrlCrossRefPubMed
  11. ↵
    1. Brockdorff N,
    2. Ashworth A,
    3. Kay GF,
    4. McCabe VM,
    5. Norris DP,
    6. Cooper PJ,
    7. Swift S,
    8. Rastan S
    (1992) The product of the mouse Xist gene is a 15 kb inactive X-specific transcript containing no conserved ORF and located in the nucleus. Cell 71:515–526. doi:10.1016/0092-8674(92)90519-I pmid:1423610
    OpenUrlCrossRefPubMed
  12. ↵
    1. Brunskill EW,
    2. Potter AS,
    3. Distasio A,
    4. Dexheimer P,
    5. Plassard A,
    6. Aronow BJ,
    7. Potter SS
    (2014) A gene expression atlas of early craniofacial development. Dev Biol 391:133–146. doi:10.1016/j.ydbio.2014.04.016 pmid:24780627
    OpenUrlCrossRefPubMed
  13. ↵
    1. Cahill KM,
    2. Huo Z,
    3. Tseng GC,
    4. Logan RW,
    5. Seney ML
    (2018) Improved identification of concordant and discordant gene expression signatures using an updated rank-rank hypergeometric overlap approach. Sci Rep 8:9588. doi:10.1038/s41598-018-27903-2 pmid:29942049
    OpenUrlCrossRefPubMed
  14. ↵
    1. Calle A,
    2. Fernandez-Gonzalez R,
    3. Ramos-Ibeas P,
    4. Laguna-Barraza R,
    5. Perez-Cerezales S,
    6. Bermejo-Alvarez P,
    7. Ramirez MA,
    8. Gutierrez-Adan A
    (2012) Long-term and transgenerational effects of in vitro culture on mouse embryos. Theriogenology 77:785–793. doi:10.1016/j.theriogenology.2011.07.016 pmid:21855990
    OpenUrlCrossRefPubMed
  15. ↵
    1. Chan JC,
    2. Morgan CP,
    3. Adrian Leu N,
    4. Shetty A,
    5. Cisse YM,
    6. Nugent BM,
    7. Morrison KE,
    8. Jašarević E,
    9. Huang W,
    10. Kanyuch N,
    11. Rodgers AB,
    12. Bhanu NV,
    13. Berger DS,
    14. Garcia BA,
    15. Ament S,
    16. Kane M,
    17. Neill Epperson C,
    18. Bale TL
    (2020) Reproductive tract extracellular vesicles are sufficient to transmit intergenerational stress and program neurodevelopment. Nat Commun 11:1499. doi:10.1038/s41467-020-15305-w pmid:32198406
    OpenUrlCrossRefPubMed
  16. ↵
    1. Chen EY,
    2. Tan CM,
    3. Kou Y,
    4. Duan Q,
    5. Wang Z,
    6. Meirelles GV,
    7. Clark NR,
    8. Ma'ayan A
    (2013) Enrichr: interactive and collaborative HTML5 gene list enrichment analysis tool. BMC Bioinformatics 14:128. doi:10.1186/1471-2105-14-128 pmid:23586463
    OpenUrlCrossRefPubMed
  17. ↵
    1. Chin CH,
    2. Chen SH,
    3. Wu HH,
    4. Ho CW,
    5. Ko MT,
    6. Lin CY
    (2014) cytoHubba: identifying hub objects and sub-networks from complex interactome. BMC Syst Biol 8:S11. doi:10.1186/1752-0509-8-S4-S11
    OpenUrlCrossRefPubMed
  18. ↵
    1. Christmas AJ,
    2. Maxwell DR
    (1970) A comparison of the effects of some benzodiazepines and other drugs on aggressive and exploratory behaviour in mice and rats. Neuropharmacology 9:17–29. doi:10.1016/0028-3908(70)90044-4 pmid:5464000
    OpenUrlCrossRefPubMed
  19. ↵
    1. Clark-Raymond A,
    2. Halaris A
    (2013) VEGF and depression: a comprehensive assessment of clinical data. J Psychiatr Res 47:1080–1087. doi:10.1016/j.jpsychires.2013.04.008 pmid:23684549
    OpenUrlCrossRefPubMed
  20. ↵
    1. Cunningham AM,
    2. Santos TL,
    3. Gutzeit VA,
    4. Hamilton H,
    5. Hen R,
    6. Donaldson ZR
    (2019) Functional interrogation of a depression-related serotonergic single nucleotide polymorphism, rs6295, using a humanized mouse model. ACS Chem Neurosci 10:3197–3206. doi:10.1021/acschemneuro.8b00638 pmid:30694044
    OpenUrlCrossRefPubMed
  21. ↵
    1. da Rocha ST,
    2. Edwards CA,
    3. Ito M,
    4. Ogata T,
    5. Ferguson-Smith AC
    (2008) Genomic imprinting at the mammalian Dlk1-Dio3 domain. Trends Genet 24:306–316. doi:10.1016/j.tig.2008.03.011 pmid:18471925
    OpenUrlCrossRefPubMed
  22. ↵
    1. de Castro Barbosa T,
    2. Ingerslev LR,
    3. Alm PS,
    4. Versteyhe S,
    5. Massart J,
    6. Rasmussen M,
    7. Donkin I,
    8. Sjögren R,
    9. Mudry JM,
    10. Vetterli L,
    11. Gupta S,
    12. Krook A,
    13. Zierath JR,
    14. Barrès R
    (2015) High-fat diet reprograms the epigenome of rat spermatozoa and transgenerationally affects metabolism of the offspring. Mol Metab 5:184–197. doi:10.1016/j.molmet.2015.12.002 pmid:26977389
    OpenUrlCrossRefPubMed
  23. ↵
    1. Deyama S,
    2. Bang E,
    3. Wohleb ES,
    4. Li XY,
    5. Kato T,
    6. Gerhard DM,
    7. Dutheil S,
    8. Dwyer JM,
    9. Taylor SR,
    10. Picciotto MR,
    11. Duman RS
    (2019) Role of neuronal VEGF signaling in the prefrontal cortex in the rapid antidepressant effects of ketamine. Am J Psychiatry 176:388–400. doi:10.1176/appi.ajp.2018.17121368 pmid:30606046
    OpenUrlCrossRefPubMed
  24. ↵
    1. Dickson DA,
    2. Paulus JK,
    3. Mensah V,
    4. Lem J,
    5. Saavedra-Rodriguez L,
    6. Gentry A,
    7. Pagidas K,
    8. Feig LA
    (2018) Reduced levels of miRNAs 449 and 34 in sperm of mice and men exposed to early life stress. Transl Psychiatry 8:101. doi:10.1038/s41398-018-0146-2 pmid:29795112
    OpenUrlCrossRefPubMed
  25. ↵
    1. Dietz DM,
    2. Laplant Q,
    3. Watts EL,
    4. Hodes GE,
    5. Russo SJ,
    6. Feng J,
    7. Oosting RS,
    8. Vialou V,
    9. Nestler EJ
    (2011) Paternal transmission of stress-induced pathologies. Biol Psychiatry 70:408–414. doi:10.1016/j.biopsych.2011.05.005 pmid:21679926
    OpenUrlCrossRefPubMed
  26. ↵
    1. Feng J,
    2. Bi C,
    3. Clark BS,
    4. Mady R,
    5. Shah P,
    6. Kohtz JD
    (2006) The Evf-2 noncoding RNA is transcribed from the Dlx-5/6 ultraconserved region and functions as a Dlx-2 transcriptional coactivator. Genes Dev 20:1470–1484. doi:10.1101/gad.1416106 pmid:16705037
    OpenUrlAbstract/FREE Full Text
  27. ↵
    1. Franklin TB,
    2. Russig H,
    3. Weiss IC,
    4. Gräff J,
    5. Linder N,
    6. Michalon A,
    7. Vizi S,
    8. Mansuy IM
    (2010) Epigenetic transmission of the impact of early stress across generations. Biol Psychiatry 68:408–415. doi:10.1016/j.biopsych.2010.05.036 pmid:20673872
    OpenUrlCrossRefPubMed
  28. ↵
    1. Gapp K,
    2. Jawaid A,
    3. Sarkies P,
    4. Bohacek J,
    5. Pelczar P,
    6. Prados J,
    7. Farinelli L,
    8. Miska E,
    9. Mansuy IM
    (2014) Implication of sperm RNAs in transgenerational inheritance of the effects of early trauma in mice. Nat Neurosci 17:667–669. doi:10.1038/nn.3695 pmid:24728267
    OpenUrlCrossRefPubMed
  29. ↵
    1. Gapp K,
    2. van Steenwyk G,
    3. Germain PL,
    4. Matsushima W,
    5. Rudolph KLM,
    6. Manuella F,
    7. Roszkowski M,
    8. Vernaz G,
    9. Ghosh T,
    10. Pelczar P,
    11. Mansuy IM,
    12. Miska EA
    (2020) Alterations in sperm long RNA contribute to the epigenetic inheritance of the effects of postnatal trauma. Mol Psychiatry 25:2162–2174.
    OpenUrl
  30. ↵
    1. Hodes GE,
    2. Pfau ML,
    3. Purushothaman I,
    4. Ahn HF,
    5. Golden SA,
    6. Christoffel DJ,
    7. Magida J,
    8. Brancato A,
    9. Takahashi A,
    10. Flanigan ME,
    11. Ménard C,
    12. Aleyasin H,
    13. Koo JW,
    14. Lorsch ZS,
    15. Feng J,
    16. Heshmati M,
    17. Wang M,
    18. Turecki G,
    19. Neve R,
    20. Zhang B, et al
    . (2015) Sex differences in nucleus accumbens transcriptome profiles associated with susceptibility versus resilience to subchronic variable stress. J Neurosci 35:16362–16376. doi:10.1523/JNEUROSCI.1392-15.2015 pmid:26674863
    OpenUrlAbstract/FREE Full Text
  31. ↵
    1. Huntley MA,
    2. Lou M,
    3. Goldstein LD,
    4. Lawrence M,
    5. Dijkgraaf GJ,
    6. Kaminker JS,
    7. Gentleman R
    (2016) Complex regulation of ADAR-mediated RNA-editing across tissues. BMC Genomics 17:61. doi:10.1186/s12864-015-2291-9 pmid:26768488
    OpenUrlCrossRefPubMed
  32. ↵
    1. Ieraci A,
    2. Mallei A,
    3. Popoli M
    (2016) Social isolation stress induces anxious-depressive-like behavior and alterations of neuroplasticity-related genes in adult male mice. Neural Plast 2016:6212983. doi:10.1155/2016/6212983 pmid:26881124
    OpenUrlCrossRefPubMed
  33. ↵
    1. Issler O,
    2. van der Zee YY,
    3. Ramakrishnan A,
    4. Wang J,
    5. Tan C,
    6. Loh YHE,
    7. Purushothaman I,
    8. Walker DM,
    9. Lorsch ZS,
    10. Hamilton PJ,
    11. Peña CJ,
    12. Flaherty E,
    13. Hartley BJ,
    14. Torres-Berrío A,
    15. Parise EM,
    16. Kronman H,
    17. Duffy JE,
    18. Estill MS,
    19. Calipari ES,
    20. Labonté B, et al
    . (2020) Sex-specific role for the long non-coding RNA LINC00473 in depression. Neuron 106:912–926. e915. doi:10.1016/j.neuron.2020.03.023
    OpenUrlCrossRef
  34. ↵
    1. Jeck WR,
    2. Sharpless NE
    (2014) Detecting and characterizing circular RNAs. Nat Biotechnol 32:453–461. doi:10.1038/nbt.2890 pmid:24811520
    OpenUrlCrossRefPubMed
  35. ↵
    1. Jiang GJ,
    2. Zhang T,
    3. An T,
    4. Zhao DD,
    5. Yang XY,
    6. Zhang DW,
    7. Zhang Y,
    8. Mu QQ,
    9. Yu N,
    10. Ma XS,
    11. Gao SH
    (2016) Differential expression of long noncoding RNAs between sperm samples from diabetic and non-diabetic mice. PLoS One 11:e0154028. doi:10.1371/journal.pone.0154028 pmid:27119337
    OpenUrlCrossRefPubMed
  36. ↵
    1. Khosla S,
    2. Dean W,
    3. Brown D,
    4. Reik W,
    5. Feil R
    (2001) Culture of preimplantation mouse embryos affects fetal development and the expression of imprinted genes. Biol Reprod 64:918–926. doi:10.1095/biolreprod64.3.918 pmid:11207209
    OpenUrlCrossRefPubMed
  37. ↵
    1. Kistler WS,
    2. Baas D,
    3. Lemeille S,
    4. Paschaki M,
    5. Seguin-Estevez Q,
    6. Barras E,
    7. Ma W,
    8. Duteyrat J-L,
    9. Morlé L,
    10. Durand B,
    11. Reith W
    (2015) RFX2 is a major transcriptional regulator of spermiogenesis. PLoS Genet 11:e1005368. doi:10.1371/journal.pgen.1005368 pmid:26162102
    OpenUrlCrossRefPubMed
  38. ↵
    1. Krishnan V,
    2. Han MH,
    3. Graham DL,
    4. Berton O,
    5. Renthal W,
    6. Russo SJ,
    7. Laplant Q,
    8. Graham A,
    9. Lutter M,
    10. Lagace DC,
    11. Ghose S,
    12. Reister R,
    13. Tannous P,
    14. Green TA,
    15. Neve RL,
    16. Chakravarty S,
    17. Kumar A,
    18. Eisch AJ,
    19. Self DW,
    20. Lee FS, et al
    . (2007) Molecular adaptations underlying susceptibility and resistance to social defeat in brain reward regions. Cell 131:391–404. doi:10.1016/j.cell.2007.09.018 pmid:17956738
    OpenUrlCrossRefPubMed
  39. ↵
    1. Labonté B,
    2. Engmann O,
    3. Purushothaman I,
    4. Menard C,
    5. Wang J,
    6. Tan C,
    7. Scarpa JR,
    8. Moy G,
    9. Loh YHE,
    10. Cahill M,
    11. Lorsch ZS,
    12. Hamilton PJ,
    13. Calipari ES,
    14. Hodes GE,
    15. Issler O,
    16. Kronman H,
    17. Pfau M,
    18. Obradovic ALJ,
    19. Dong Y,
    20. Neve RL, et al
    . (2017) Sex-specific transcriptional signatures in human depression. Nat Med 23:1102–1111. doi:10.1038/nm.4386 pmid:28825715
    OpenUrlCrossRefPubMed
  40. ↵
    1. Langfelder P,
    2. Horvath S
    (2008) WGCNA: an R package for weighted correlation network analysis. BMC Bioinformatics 9:559. doi:10.1186/1471-2105-9-559
    OpenUrlCrossRefPubMed
  41. ↵
    1. Lanzós A,
    2. Carlevaro-Fita J,
    3. Mularoni L,
    4. Reverter F,
    5. Palumbo E,
    6. Guigó R,
    7. Johnson R
    (2017) Discovery of cancer driver long noncoding RNAs across 1112 tumour genomes: new candidates and distinguishing features. Sci Rep 7:41544. doi:10.1038/srep41544 pmid:28128360
    OpenUrlCrossRefPubMed
  42. ↵
    1. LaPlant Q,
    2. Chakravarty S,
    3. Vialou V,
    4. Mukherjee S,
    5. Koo JW,
    6. Kalahasti G,
    7. Bradbury KR,
    8. Taylor SV,
    9. Maze I,
    10. Kumar A,
    11. Graham A,
    12. Birnbaum SG,
    13. Krishnan V,
    14. Truong H-T,
    15. Neve RL,
    16. Nestler EJ,
    17. Russo SJ
    (2009) Role of nuclear factor kappaB in ovarian hormone-mediated stress hypersensitivity in female mice. Biol Psychiatry 65:874–880. doi:10.1016/j.biopsych.2009.01.024 pmid:19251249
    OpenUrlCrossRefPubMed
  43. ↵
    1. Lehti MS,
    2. Sironen A
    (2017) Formation and function of sperm tail structures in association with sperm motility defects. Biol Reprod 97:522–536. doi:10.1093/biolre/iox096 pmid:29024992
    OpenUrlCrossRefPubMed
  44. ↵
    1. Lorsch ZS,
    2. Loh YHE,
    3. Purushothaman I,
    4. Walker DM,
    5. Parise EM,
    6. Salery M,
    7. Cahill ME,
    8. Hodes GE,
    9. Pfau ML,
    10. Kronman H,
    11. Hamilton PJ,
    12. Issler O,
    13. Labonté B,
    14. Symonds AE,
    15. Zucker M,
    16. Zhang TY,
    17. Meaney MJ,
    18. Russo SJ,
    19. Shen L,
    20. Bagot RC, et al
    . (2018) Estrogen receptor alpha drives pro-resilient transcription in mouse models of depression. Nat Commun 9:1116. doi:10.1038/s41467-018-03567-4 pmid:29549264
    OpenUrlCrossRefPubMed
  45. ↵
    1. Maeng LY,
    2. Milad MR
    (2015) Sex differences in anxiety disorders: interactions between fear, stress, and gonadal hormones. Horm Behav 76:106–117. doi:10.1016/j.yhbeh.2015.04.002 pmid:25888456
    OpenUrlCrossRefPubMed
  46. ↵
    1. Manners MT,
    2. Yohn NL,
    3. Lahens NF,
    4. Grant GR,
    5. Bartolomei MS,
    6. Blendy JA
    (2019) Transgenerational inheritance of chronic adolescent stress: effects of stress response and the amygdala transcriptome. Genes Brain Behav 18:e12493. doi:10.1111/gbb.12493 pmid:29896789
    OpenUrlCrossRefPubMed
  47. ↵
    1. Mashoodh R,
    2. Habrylo IB,
    3. Gudsnuk KM,
    4. Pelle G,
    5. Champagne FA
    (2018) Maternal modulation of paternal effects on offspring development. Proc Biol Sci 285:20180118.
    OpenUrlCrossRefPubMed
  48. ↵
    1. Morgan CP,
    2. Shetty AC,
    3. Chan JC,
    4. Berger DS,
    5. Ament SA,
    6. Epperson CN,
    7. Bale TL
    (2020) Repeated sampling facilitates within- and between-subject modeling of the human sperm transcriptome to identify dynamic and stress-responsive sncRNAs. Sci Rep 10:17498. doi:10.1038/s41598-020-73867-7 pmid:33060642
    OpenUrlCrossRefPubMed
  49. ↵
    1. Nasca C,
    2. Menard C,
    3. Hodes G,
    4. Bigio B,
    5. Pena C,
    6. Lorsch Z,
    7. Zelli D,
    8. Ferris A,
    9. Kana V,
    10. Purushothaman I,
    11. Dobbin J,
    12. Nassim M,
    13. DeAngelis P,
    14. Merad M,
    15. Rasgon N,
    16. Meaney M,
    17. Nestler EJ,
    18. McEwen BS,
    19. Russo SJ
    (2019) Multidimensional predictors of susceptibility and resilience to social defeat stress. Biol Psychiatry 86:483–491. doi:10.1016/j.biopsych.2019.06.030 pmid:31466563
    OpenUrlCrossRefPubMed
  50. ↵
    1. Parvatiyar MS,
    2. Marshall JL,
    3. Nguyen RT,
    4. Jordan MC,
    5. Richardson VA,
    6. Roos KP,
    7. Crosbie-Watson RH
    (2015) Sarcospan regulates cardiac isoproterenol response and prevents Duchenne muscular dystrophy-associated cardiomyopathy. J Am Heart Assoc 4:e002481.
    OpenUrlAbstract/FREE Full Text
  51. ↵
    1. Parvatiyar MS,
    2. Brownstein AJ,
    3. Kanashiro-Takeuchi RM,
    4. Collado JR,
    5. Dieseldorff Jones KM,
    6. Gopal J,
    7. Hammond KG,
    8. Marshall JL,
    9. Ferrel A,
    10. Beedle AM,
    11. Chamberlain JS,
    12. Renato Pinto J,
    13. Crosbie RH
    (2019) Stabilization of the cardiac sarcolemma by sarcospan rescues DMD-associated cardiomyopathy. JCI Insight 5:e123855. doi:10.1172/jci.insight.123855
    OpenUrlCrossRef
  52. ↵
    1. Peña CJ,
    2. Nestler EJ,
    3. Bagot RC
    (2019a) Environmental programming of susceptibility and resilience to stress in adulthood in male mice. Front Behav Neurosci 13:40.
    OpenUrl
  53. ↵
    1. Peña CJ,
    2. Smith M,
    3. Ramakrishnan A,
    4. Cates HM,
    5. Bagot RC,
    6. Kronman HG,
    7. Patel B,
    8. Chang AB,
    9. Purushothaman I,
    10. Dudley J,
    11. Morishita H,
    12. Shen L,
    13. Nestler EJ
    (2019b) Early life stress alters transcriptomic patterning across reward circuitry in male and female mice. Nat Commun 10:5098. doi:10.1038/s41467-019-13085-6 pmid:31704941
    OpenUrlCrossRefPubMed
  54. ↵
    1. Prut L,
    2. Belzung C
    (2003) The open field as a paradigm to measure the effects of drugs on anxiety-like behaviors: a review. Eur J Pharmacol 463:3–33. doi:10.1016/S0014-2999(03)01272-X
    OpenUrlCrossRefPubMed
  55. ↵
    1. Reeb RM,
    2. Rogers BP,
    3. Richardson B
    (1992) Research: the bridge to collaboration. Nursingconnections 5:31–36.
    OpenUrlPubMed
  56. ↵
    1. Reik W,
    2. Walter J
    (2001) Genomic imprinting: parental influence on the genome. Nat Rev Genet 2:21–32. doi:10.1038/35047554 pmid:11253064
    OpenUrlCrossRefPubMed
  57. ↵
    1. Rinn JL,
    2. Chang HY
    (2012) Genome regulation by long noncoding RNAs. Annu Rev Biochem 81:145–166. doi:10.1146/annurev-biochem-051410-092902 pmid:22663078
    OpenUrlCrossRefPubMed
  58. ↵
    1. Rodgers AB,
    2. Morgan CP,
    3. Bronson SL,
    4. Revello S,
    5. Bale TL
    (2013) Paternal stress exposure alters sperm microRNA content and reprograms offspring HPA stress axis regulation. J Neurosci 33:9003–9012. doi:10.1523/JNEUROSCI.0914-13.2013 pmid:23699511
    OpenUrlAbstract/FREE Full Text
  59. ↵
    1. Rodgers AB,
    2. Morgan CP,
    3. Leu NA,
    4. Bale TL
    (2015) Transgenerational epigenetic programming via sperm microRNA recapitulates effects of paternal stress. Proc Natl Acad Sci USA 112:13699–13704. doi:10.1073/pnas.1508347112 pmid:26483456
    OpenUrlAbstract/FREE Full Text
  60. ↵
    1. Roy B,
    2. Wang Q,
    3. Dwivedi Y
    (2018) Long noncoding RNA-associated transcriptomic changes in resiliency or susceptibility to depression and response to antidepressant treatment. Int J Neuropsychopharmacol 21:461–472. doi:10.1093/ijnp/pyy010 pmid:29390069
    OpenUrlCrossRefPubMed
  61. ↵
    1. Royo H,
    2. Cavaillé J
    (2008) Non-coding RNAs in imprinted gene clusters. Biol Cell 100:149–166. doi:10.1042/BC20070126 pmid:18271756
    OpenUrlCrossRefPubMed
  62. ↵
    1. Seney ML,
    2. Huo Z,
    3. Cahill K,
    4. French L,
    5. Puralewski R,
    6. Zhang J,
    7. Logan RW,
    8. Tseng G,
    9. Lewis DA,
    10. Sibille E
    (2018) Opposite molecular signatures of depression in men and women. Biol Psychiatry 84:18–27. doi:10.1016/j.biopsych.2018.01.017 pmid:29548746
    OpenUrlCrossRefPubMed
  63. ↵
    1. Sleutels F,
    2. Zwart R,
    3. Barlow DP
    (2002) The non-coding Air RNA is required for silencing autosomal imprinted genes. Nature 415:810–813. doi:10.1038/415810a pmid:11845212
    OpenUrlCrossRefPubMed
  64. ↵
    1. Song X,
    2. Lin NH,
    3. Wang YL,
    4. Chen B,
    5. Wang HX,
    6. Hu K
    (2019) Comprehensive transcriptome analysis based on RNA sequencing identifies critical genes for lipopolysaccharide-induced epididymitis in a rat model. Asian J Androl 21:605–611. doi:10.4103/aja.aja_21_19 pmid:31044753
    OpenUrlCrossRefPubMed
  65. ↵
    1. Soumillon M,
    2. Necsulea A,
    3. Weier M,
    4. Brawand D,
    5. Zhang X,
    6. Gu H,
    7. Barthès P,
    8. Kokkinaki M,
    9. Nef S,
    10. Gnirke A,
    11. Dym M,
    12. de Massy B,
    13. Mikkelsen TS,
    14. Kaessmann H
    (2013) Cellular source and mechanisms of high transcriptome complexity in the mammalian testis. Cell Rep 3:2179–2190. doi:10.1016/j.celrep.2013.05.031 pmid:23791531
    OpenUrlCrossRefPubMed
  66. ↵
    1. Stone BJ,
    2. Steele KH,
    3. Fath-Goodin A
    (2015) A rapid and effective nonsurgical artificial insemination protocol using the NSET device for sperm transfer in mice without anesthesia. Transgenic Res 24:775–781. doi:10.1007/s11248-015-9887-3 pmid:26065409
    OpenUrlCrossRefPubMed
  67. ↵
    1. Ventura-Juncá P,
    2. Irarrázaval I,
    3. Rolle AJ,
    4. Gutiérrez JI,
    5. Moreno RD,
    6. Santos MJ
    (2015) In vitro fertilization (IVF) in mammals: epigenetic and developmental alterations. Scientific and bioethical implications for IVF in humans. Biol Res 48:68. doi:10.1186/s40659-015-0059-y pmid:26683055
    OpenUrlCrossRefPubMed
  68. ↵
    1. Vialou V,
    2. Robison AJ,
    3. Laplant QC,
    4. Covington HE 3rd.,
    5. Dietz DM,
    6. Ohnishi YN,
    7. Mouzon E,
    8. Rush AJ 3rd.,
    9. Watts EL,
    10. Wallace DL,
    11. Iñiguez SD,
    12. Ohnishi YH,
    13. Steiner MA,
    14. Warren BL,
    15. Krishnan V,
    16. Bolaños CA,
    17. Neve RL,
    18. Ghose S,
    19. Berton O,
    20. Tamminga CA, et al
    . (2010) DeltaFosB in brain reward circuits mediates resilience to stress and antidepressant responses. Nat Neurosci 13:745–752. doi:10.1038/nn.2551 pmid:20473292
    OpenUrlCrossRefPubMed
  69. ↵
    1. Wang M,
    2. Zhao Y,
    3. Zhang B
    (2015) Efficient test and visualization of multi-set intersections. Sci Rep 5:16923. doi:10.1038/srep16923 pmid:26603754
    OpenUrlCrossRefPubMed
  70. ↵
    1. Wen K,
    2. Yang L,
    3. Xiong T,
    4. Di C,
    5. Ma D,
    6. Wu M,
    7. Xue Z,
    8. Zhang X,
    9. Long L,
    10. Zhang W,
    11. Zhang J,
    12. Bi X,
    13. Dai J,
    14. Zhang Q,
    15. Lu ZJ,
    16. Gao G
    (2016) Critical roles of long noncoding RNAs in Drosophila spermatogenesis. Genome Res 26:1233–1244. doi:10.1101/gr.199547.115 pmid:27516619
    OpenUrlAbstract/FREE Full Text
  71. ↵
    1. Wilkinson MB,
    2. Xiao G,
    3. Kumar A,
    4. LaPlant Q,
    5. Renthal W,
    6. Sikder D,
    7. Kodadek TJ,
    8. Nestler EJ
    (2009) Imipramine treatment and resiliency exhibit similar chromatin regulation in the mouse nucleus accumbens in depression models. J Neurosci 29:7820–7832. doi:10.1523/JNEUROSCI.0932-09.2009 pmid:19535594
    OpenUrlAbstract/FREE Full Text
  72. ↵
    1. Wook Koo J,
    2. Labonté B,
    3. Engmann O,
    4. Calipari ES,
    5. Juarez B,
    6. Lorsch Z,
    7. Walsh JJ,
    8. Friedman AK,
    9. Yorgason JT,
    10. Han M-H,
    11. Nestler EJ
    (2016) Essential role of mesolimbic brain-derived neurotrophic factor in chronic social stress-induced depressive behaviors. Biol Psychiatry 80:469–478. doi:10.1016/j.biopsych.2015.12.009 pmid:26858215
    OpenUrlCrossRefPubMed
  73. ↵
    1. Wray NR,
    2. Ripke S,
    3. Mattheisen M,
    4. Trzaskowski M,
    5. Byrne EM,
    6. Abdellaoui A,
    7. Adams MJ,
    8. Agerbo E,
    9. Air TM,
    10. Andlauer TMF,
    11. Bacanu SA,
    12. Bækvad-Hansen M,
    13. Beekman AFT,
    14. Bigdeli TB,
    15. Binder EB,
    16. Blackwood DRH,
    17. Bryois J,
    18. Buttenschøn HN,
    19. Bybjerg-Grauholm J,
    20. Cai N, et al
    . (2018) Genome-wide association analyses identify 44 risk variants and refine the genetic architecture of major depression. Nat Genet 50:668–681. doi:10.1038/s41588-018-0090-3 pmid:29700475
    OpenUrlCrossRefPubMed
  74. ↵
    1. Zhang L,
    2. Lu H,
    3. Xin D,
    4. Cheng H,
    5. Zhou R
    (2010) A novel ncRNA gene from mouse chromosome 5 trans-splices with Dmrt1 on chromosome 19. Biochem Biophys Res Commun 400:696–700. doi:10.1016/j.bbrc.2010.08.130 pmid:20816665
    OpenUrlCrossRefPubMed
  75. ↵
    1. Zhang X,
    2. Zhang P,
    3. Song D,
    4. Xiong S,
    5. Zhang H,
    6. Fu J,
    7. Gao F,
    8. Chen H,
    9. Zeng X
    (2019) Expression profiles and characteristics of human lncRNA in normal and asthenozoospermia sperm. Biol Reprod 100:982–993. doi:10.1093/biolre/ioy253 pmid:30517597
    OpenUrlCrossRefPubMed
Back to top

In this issue

The Journal of Neuroscience: 41 (29)
Journal of Neuroscience
Vol. 41, Issue 29
21 Jul 2021
  • Table of Contents
  • Table of Contents (PDF)
  • About the Cover
  • Index by author
  • Ed Board (PDF)
Email

Thank you for sharing this Journal of Neuroscience article.

NOTE: We request your email address only to inform the recipient that it was you who recommended this article, and that it is not junk mail. We do not retain these email addresses.

Enter multiple addresses on separate lines or separate them with commas.
Sperm Transcriptional State Associated with Paternal Transmission of Stress Phenotypes
(Your Name) has forwarded a page to you from Journal of Neuroscience
(Your Name) thought you would be interested in this article in Journal of Neuroscience.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Print
View Full Page PDF
Citation Tools
Sperm Transcriptional State Associated with Paternal Transmission of Stress Phenotypes
Ashley M. Cunningham, Deena M. Walker, Aarthi Ramakrishnan, Marie A. Doyle, Rosemary C. Bagot, Hannah M. Cates, Catherine J. Peña, Orna Issler, Casey K. Lardner, Caleb Browne, Scott J. Russo, Li Shen, Eric J. Nestler
Journal of Neuroscience 21 July 2021, 41 (29) 6202-6216; DOI: 10.1523/JNEUROSCI.3192-20.2021

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Respond to this article
Request Permissions
Share
Sperm Transcriptional State Associated with Paternal Transmission of Stress Phenotypes
Ashley M. Cunningham, Deena M. Walker, Aarthi Ramakrishnan, Marie A. Doyle, Rosemary C. Bagot, Hannah M. Cates, Catherine J. Peña, Orna Issler, Casey K. Lardner, Caleb Browne, Scott J. Russo, Li Shen, Eric J. Nestler
Journal of Neuroscience 21 July 2021, 41 (29) 6202-6216; DOI: 10.1523/JNEUROSCI.3192-20.2021
del.icio.us logo Digg logo Reddit logo Twitter logo Facebook logo Google logo Mendeley logo
  • Tweet Widget
  • Facebook Like
  • Google Plus One

Jump to section

  • Article
    • Abstract
    • Introduction
    • Materials and Methods
    • Results
    • Discussion
    • Footnotes
    • References
  • Figures & Data
  • Info & Metrics
  • eLetters
  • PDF

Keywords

  • anxiety
  • intergenerational stress
  • RNA-sequencing
  • sperm

Responses to this article

Respond to this article

Jump to comment:

No eLetters have been published for this article.

Related Articles

Cited By...

More in this TOC Section

Research Articles

  • Evidence that ultrafast non-quantal transmission underlies synchronized vestibular action potential generation
  • Nfia is Critical for AII Amacrine Cell Production: Selective Bipolar Cell Dependencies and Diminished ERG
  • Multimodal Imaging for Validation and Optimization of Ion Channel-Based Chemogenetics in Nonhuman Primates
Show more Research Articles

Cellular/Molecular

  • Cleavage of VAMP2/3 Affects Oligodendrocyte Lineage Development in the Developing Mouse Spinal Cord
  • P2Y6 receptor-dependent microglial phagocytosis of synapses during development regulates synapse density and memory
  • Acute ethanol modulates synaptic inhibition in the basolateral amygdala via rapid NLRP3 inflammasome activation and regulates anxiety-like behavior in rats
Show more Cellular/Molecular
  • Home
  • Alerts
  • Visit Society for Neuroscience on Facebook
  • Follow Society for Neuroscience on Twitter
  • Follow Society for Neuroscience on LinkedIn
  • Visit Society for Neuroscience on Youtube
  • Follow our RSS feeds

Content

  • Early Release
  • Current Issue
  • Issue Archive
  • Collections

Information

  • For Authors
  • For Advertisers
  • For the Media
  • For Subscribers

About

  • About the Journal
  • Editorial Board
  • Privacy Policy
  • Contact
(JNeurosci logo)
(SfN logo)

Copyright © 2023 by the Society for Neuroscience.
JNeurosci Online ISSN: 1529-2401

The ideas and opinions expressed in JNeurosci do not necessarily reflect those of SfN or the JNeurosci Editorial Board. Publication of an advertisement or other product mention in JNeurosci should not be construed as an endorsement of the manufacturer’s claims. SfN does not assume any responsibility for any injury and/or damage to persons or property arising from or related to any use of any material contained in JNeurosci.