The current study used two inbred mouse strains, C57BL/6 J and A/J, to investigate the genetics of behavioral responses to fentanyl. Mice were tested for conditioned place preference and fentanyl-induced locomotor activity. C57BL/6J mice formed a conditioned place preference to fentanyl injections and fentanyl increased their activity. Neither effect was noted in A/J mice. We conducted RNA-sequencing on the nucleus accumbens of mice used for fentanyl-induced locomotor activity. Surprisingly, we noted few differentially expressed genes using treatment as the main factor. However many genes differed between strains.
Authors:
Samuel J Harp, Mariangela Martini, Will Rosenow, Larry D Mesner, Hugh Johnson, Charles R Farber, Emilie F Rissman
The current study used two inbred mouse strains, C57BL/6 J and A/J, to investigate the genetics of behavioral responses to fentanyl. Mice were tested for conditioned place preference and fentanyl-induced locomotor activity. C57BL/6J mice formed a conditioned place preference to fentanyl injections and fentanyl increased their activity. Neither effect was noted in A/J mice. We conducted RNA-sequencing on the nucleus accumbens of mice used for fentanyl-induced locomotor activity. Surprisingly, we noted few differentially expressed genes using treatment as the main factor. However many genes differed between strains.
Authors:
Samuel J Harp, Mariangela Martini, Will Rosenow, Larry D Mesner, Hugh Johnson, Charles R Farber, Emilie F Rissman
The current study used two inbred mouse strains, C57BL/6 J and A/J, to investigate the genetics of behavioral responses to fentanyl. Mice were tested for conditioned place preference and fentanyl-induced locomotor activity. C57BL/6J mice formed a conditioned place preference to fentanyl injections and fentanyl increased their activity. Neither effect was noted in A/J mice. We conducted RNA-sequencing on the nucleus accumbens of mice used for fentanyl-induced locomotor activity. Surprisingly, we noted few differentially expressed genes using treatment as the main factor. However many genes differed between strains.
Authors:
Samuel J Harp, Mariangela Martini, Will Rosenow, Larry D Mesner, Hugh Johnson, Charles R Farber, Emilie F Rissman
Differential gene expression in nucleus accumbens somatostatin interneurons_cocaine_mice_logFC
Description:
To characterize transcriptional alterations that cocaine induces in these cells, we perform cell type-specific RNA-sequencing on FACS-isolated nuclei of somatostatin interneurons and identified 1100 DETs enriched for processes related to neural plasticity. To profile the entire (non poly-A selected) transcriptome of NAc somatostatin interneurons, we generated a transgenic reporter line (SST-TLG498 mice) to label the nuclei of these cells with a modified form of EGFP that is retained in the nuclear membrane (EGFP-F)22, enabling their isolation from NAc dissections using FACS. We succeeded in FACS-isolating nuclei suitable for RNA-sequencing from individual SST-TLG498 mice. We proceeded with differential expression analysis of the RNA-sequencing data to identify differentially expressed transcripts (DETs) in NAc somatostatin interneurons in response to repeated cocaine exposure: 778 transcripts were upregulated by cocaine and 322 were downregulated.
Authors:
Efrain A Ribeiro, Marine Salery, Joseph R Scarpa, Erin S Calipari, Peter J Hamilton, Stacy M Ku, Hope Kronman, Immanuel Purushothaman, Barbara Juarez, Mitra Heshmati, Marie Doyle, Casey Lardner, Dominicka Burek, Ana Strat, Stephen Pirpinias, Ezekiell Mouzon, Ming-Hu Han, Rachael L Neve, Rosemary C Bagot, Andrew Kasarskis, Ja Wook Koo, Eric J Nestler
DEG male mouse forebrain PND14 morphine vs saline_pvalue
Description:
To examine forebrain transcriptomic changes that might elucidate mechanisms of withdrawal, delayed development, and any long-term behavior changes, we generated transcriptomic signatures following our “3-trimester” exposure model (3-Tri). In addition, we also examined transcriptomes from animals that received opioids only during the gestational period (PND1) or only during the last trimester from PND 1–14 (PND 14). We sought to determine whether transcriptomic signatures vary based on the window of exposure, perhaps contributing to the discrepancies in the literature regarding acute and long-term outcomes. Brains were dissected from PND 1 pups 6 h after discovery. Brains were dissected from post-natal exposure only (PND 14) or 3-trimester exposure (3-tri) 6 h after the last morphine or saline injection. The number of animals per group was similar (N = 5–7 animals, male and female C57Bl/6NTac mice), and the quality controls, library construction and sequence parameters were also identical across all groups. Libraries were sequenced on a NovaSeq 6000 at a depth of 30 million total reads/sample using paired-end sequencing of 150 base pairs (PE150), to a depth of 30 million total reads/sample. Reads were then mapped to the mouse reference genome (Mus Musculus, GRCm38/mm10) using HISAT2 (version 2.2.1), and duplicated fragments were removed using Picard MarkDuplicates. Differential expression analysis between two conditions (e.g., Morphine and Saline) was performed in R (version 4.1.1) with DESeq2 (v1.32.0) package. Genes were assigned by the authors as differentially expressed if the (adjusted) (nominal) p-value < 0.05. All genes/scores are presented here.
Authors:
Amelia D Dunn, Shivon A Robinson, Chiso Nwokafor, Molly Estill, Julia Ferrante, Li Shen, Crystal O Lemchi, Jordi Creus-Muncunill, Angie Ramirez, Juliet Mengaziol, Julia K Brynildsen, Mark Leggas, Jamie Horn, Michelle E Ehrlich, Julie A Blendy
DEG female mouse forebrain PND1 morphine vs saline_pvalue
Description:
To examine forebrain transcriptomic changes that might elucidate mechanisms of withdrawal, delayed development, and any long-term behavior changes, we generated transcriptomic signatures following our “3-trimester” exposure model (3-Tri). In addition, we also examined transcriptomes from animals that received opioids only during the gestational period (PND1) or only during the last trimester from PND 1–14 (PND 14). We sought to determine whether transcriptomic signatures vary based on the window of exposure, perhaps contributing to the discrepancies in the literature regarding acute and long-term outcomes. Brains were dissected from PND 1 pups 6 h after discovery. Brains were dissected from post-natal exposure only (PND 14) or 3-trimester exposure (3-tri) 6 h after the last morphine or saline injection. The number of animals per group was similar (N = 5–7 animals, male and female C57Bl/6NTac mice), and the quality controls, library construction and sequence parameters were also identical across all groups. Libraries were sequenced on a NovaSeq 6000 at a depth of 30 million total reads/sample using paired-end sequencing of 150 base pairs (PE150), to a depth of 30 million total reads/sample. Reads were then mapped to the mouse reference genome (Mus Musculus, GRCm38/mm10) using HISAT2 (version 2.2.1), and duplicated fragments were removed using Picard MarkDuplicates. Differential expression analysis between two conditions (e.g., Morphine and Saline) was performed in R (version 4.1.1) with DESeq2 (v1.32.0) package. Genes were assigned by the authors as differentially expressed if the (adjusted) (nominal) p-value < 0.05. All genes/scores are presented here.
Authors:
Amelia D Dunn, Shivon A Robinson, Chiso Nwokafor, Molly Estill, Julia Ferrante, Li Shen, Crystal O Lemchi, Jordi Creus-Muncunill, Angie Ramirez, Juliet Mengaziol, Julia K Brynildsen, Mark Leggas, Jamie Horn, Michelle E Ehrlich, Julie A Blendy
Perinatal fentanyl exposed DEG in male mouse VTA_pvalue
Description:
We identified individual genes with significant gene expression changes following perinatal fentanyl exposure. Pregnant C57BL/6 dams received 10μg/ml fentanyl in their drinking water from embryonic day 0 (E0) through birth until weaning at postnatal day 21 (P21). Liquid consumption by dams all through gestation to weaning was comparable between water and fentanyl in water group. Post-weaning, at P35, we performed an unbiased transcriptomic analysis in perinatal fentanyl exposed mice using high-throughput RNAseq of bulk-tissue punches from VTA, NAc, PrL, S1 and VBT from four juvenile (P35) mice per sex/treatment. For RNA sequencing, only samples with RNA integrity numbers >8 were used. Samples were sequenced on an Illumina HiSeq 4000 with a 75bp paired-end read. 64–100 million reads were obtained for each sample. Reads were aligned to the mouse genome (Mus musculus. GRCm38) using TopHat2 (version 2.0.8; maximum number of mismatches=2; segment length=30; maximum multi-hits per read=25; maximum intron length=50,000). A False Discovery Rate (FDR)<0.05 was used to define differentially expressed genes.
Authors:
Jimmy Olusakin, Gautam Kumar, Mahashweta Basu, Cali A Calarco, Megan E Fox, Jason B Alipio, Catherine Haga, Makeda D Turner, Asaf Keller, Seth A Ament, Mary Kay Lobo
Perinatal fentanyl exposed DEG in male mouse VTA_qvalue
Description:
We identified individual genes with significant gene expression changes following perinatal fentanyl exposure. Pregnant C57BL/6 dams received 10μg/ml fentanyl in their drinking water from embryonic day 0 (E0) through birth until weaning at postnatal day 21 (P21). Liquid consumption by dams all through gestation to weaning was comparable between water and fentanyl in water group. Post-weaning, at P35, we performed an unbiased transcriptomic analysis in perinatal fentanyl exposed mice using high-throughput RNAseq of bulk-tissue punches from VTA, NAc, PrL, S1 and VBT from four juvenile (P35) mice per sex/treatment. For RNA sequencing, only samples with RNA integrity numbers >8 were used. Samples were sequenced on an Illumina HiSeq 4000 with a 75bp paired-end read. 64–100 million reads were obtained for each sample. Reads were aligned to the mouse genome (Mus musculus. GRCm38) using TopHat2 (version 2.0.8; maximum number of mismatches=2; segment length=30; maximum multi-hits per read=25; maximum intron length=50,000). A False Discovery Rate (FDR)<0.05 was used to define differentially expressed genes.
Authors:
Jimmy Olusakin, Gautam Kumar, Mahashweta Basu, Cali A Calarco, Megan E Fox, Jason B Alipio, Catherine Haga, Makeda D Turner, Asaf Keller, Seth A Ament, Mary Kay Lobo
Average rotarod training latency Chr# 5 rs13478110 (9741228) with right flanking marker rs13478092(3595407) and left marker rs3718776 (150393227). This was mapped in 300 + (b6x129)F2 mice.
Genes identified as expressed higher (up) in the AJ strain than in the AJ strain. Differentially expressed genes had a Q-value < 0.05 following the Benjamini-Hochberg methodology for false discovery rates in the limma+voom pipeline within edgeR. Q-value is reported from the topTable function.
Genes identified as expressed lower (down) in the AJ strain than in the CAST strain. Differentially expressed genes had a Q-value < 0.05 following the Benjamini-Hochberg methodology for false discovery rates in the limma+voom pipeline within edgeR. Q-value is reported from the topTable function.
Genes identified as expressed higher (up) in the AJ strain than in the NZO strain. Differentially expressed genes had a Q-value < 0.05 following the Benjamini-Hochberg methodology for false discovery rates in the limma+voom pipeline within edgeR. Q-value is reported from the topTable function.
Genes identified as expressed lower (down) in the AJ strain than in the PWK strain. Differentially expressed genes had a Q-value < 0.05 following the Benjamini-Hochberg methodology for false discovery rates in the limma+voom pipeline within edgeR. Q-value is reported from the topTable function.
Genes identified as expressed higher (up) in the AJ strain than in the S129 strain. Differentially expressed genes had a Q-value < 0.05 following the Benjamini-Hochberg methodology for false discovery rates in the limma+voom pipeline within edgeR. Q-value is reported from the topTable function.
Genes identified as expressed higher (up) in the AJ strain than in the WSB strain. Differentially expressed genes had a Q-value < 0.05 following the Benjamini-Hochberg methodology for false discovery rates in the limma+voom pipeline within edgeR. Q-value is reported from the topTable function.
Genes identified as expressed lower (down) in the AJ strain than in the CAST strain. Differentially expressed genes had a Q-value < 0.05 following the Benjamini-Hochberg methodology for false discovery rates in the limma+voom pipeline within edgeR. Q-value is reported from the topTable function.
Genes identified as expressed higher (up) in the AJ strain than in the AJ strain. Differentially expressed genes had a Q-value < 0.05 following the Benjamini-Hochberg methodology for false discovery rates in the limma+voom pipeline within edgeR. Q-value is reported from the topTable function.
Genes identified as expressed higher (up) in the AJ strain than in the NOD strain. Differentially expressed genes had a Q-value < 0.05 following the Benjamini-Hochberg methodology for false discovery rates in the limma+voom pipeline within edgeR. Q-value is reported from the topTable function.
Genes identified as expressed higher (up) in the AJ strain than in the NOD strain. Differentially expressed genes had a Q-value < 0.05 following the Benjamini-Hochberg methodology for false discovery rates in the limma+voom pipeline within edgeR. Q-value is reported from the topTable function.
Genes identified as expressed higher (up) in the AJ strain than in the NZO strain. Differentially expressed genes had a Q-value < 0.05 following the Benjamini-Hochberg methodology for false discovery rates in the limma+voom pipeline within edgeR. Q-value is reported from the topTable function.
Genes identified as expressed higher (up) in the AJ strain than in the NZO strain. Differentially expressed genes had a Q-value < 0.05 following the Benjamini-Hochberg methodology for false discovery rates in the limma+voom pipeline within edgeR. Q-value is reported from the topTable function.
Genes identified as expressed lower (down) in the AJ strain than in the PWK strain. Differentially expressed genes had a Q-value < 0.05 following the Benjamini-Hochberg methodology for false discovery rates in the limma+voom pipeline within edgeR. Q-value is reported from the topTable function.
Genes identified as expressed higher (up) in the AJ strain than in the S129 strain. Differentially expressed genes had a Q-value < 0.05 following the Benjamini-Hochberg methodology for false discovery rates in the limma+voom pipeline within edgeR. Q-value is reported from the topTable function.
Genes identified as expressed higher (up) in the AJ strain than in the S129 strain. Differentially expressed genes had a Q-value < 0.05 following the Benjamini-Hochberg methodology for false discovery rates in the limma+voom pipeline within edgeR. Q-value is reported from the topTable function.
GeneWeaver