Next Article in Journal
Handling Dataset with Geophysical and Geological Variables on the Bolivian Andes by the GMT Scripts
Previous Article in Journal
Quality Control Impacts on Total Precipitation Gauge Records for Montane Valley and Ridge Sites in SW Alberta, Canada
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Data Descriptor

Longitudinal RNA Sequencing of Skin and DRG Neurons in Mice with Paclitaxel-Induced Peripheral Neuropathy

by
Anthony M. Cirrincione
1,†,
Cassandra A. Reimonn
2,†,
Benjamin J. Harrison
2,* and
Sandra Rieger
1,3,*
1
Department of Biology, University of Miami, Coral Gables, FL 33146, USA
2
Department of Biomedical Sciences, University of New England, Biddeford, ME 04005, USA
3
Sylvester Comprehensive Cancer Center, Miller School of Medicine, University of Miami, Miami, FL 33136, USA
*
Authors to whom correspondence should be addressed.
These authors contributed equally to this work.
Submission received: 22 April 2022 / Revised: 26 May 2022 / Accepted: 26 May 2022 / Published: 30 May 2022

Abstract

:
Paclitaxel-induced peripheral neuropathy is a condition of nerve degeneration induced by chemotherapy, which afflicts up to 70% of treated patients. Therapeutic interventions are unavailable due to an incomplete understanding of the underlying mechanisms. We previously discovered that major physiological changes in the skin underlie paclitaxel-induced peripheral neuropathy in zebrafish and rodents. The precise molecular mechanisms are only incompletely understood. For instance, paclitaxel induces the upregulation of MMP-13, which, when inhibited, prevents axon degeneration. To better understand other gene regulatory changes induced by paclitaxel, we induced peripheral neuropathy in mice following intraperitoneal injection either with vehicle or paclitaxel every other day four times total. Skin and dorsal root ganglion neurons were collected based on distinct behavioural responses categorised as “pain onset” (d4), “maximal pain” (d7), “beginning of pain resolution” (d11), and “recovery phase” (d23) for comparative longitudinal RNA sequencing. The generated datasets validate previous discoveries and reveal additional gene expression changes that warrant further validation with the goal to aid in the development of drugs that prevent or reverse paclitaxel-induced peripheral neuropathy.
Dataset: Research Square pre-print: https://doi.org/10.21203/rs.3.rs-960602/v1.
Data Deposition: The data discussed in this publication have been deposited in NCBI’s Gene Expression Omnibus (Cirrincione, et al. (2021) and are accessible through GEO series accession number GSE185084 (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE185084).
Data License: CC0

1. Summary

Paclitaxel-induced peripheral neuropathy is a condition of nerve degeneration induced by chemotherapy, which afflicts up to 70% of treated patients. Therapeutic interventions are unavailable due to an incomplete understanding of the underlying mechanisms. We previously discovered that major physiological changes in the skin underlie paclitaxel-induced peripheral neuropathy in zebrafish and rodents. The precise molecular mechanisms are only incompletely understood. For instance, paclitaxel induces the upregulation of MMP-13, which, when inhibited, prevents axon degeneration. To better understand other gene regulatory changes induced by paclitaxel, we induced peripheral neuropathy in mice following intraperitoneal injection either with vehicle or paclitaxel every other day four times total. Skin and dorsal root ganglion neurons were collected based on distinct behavioural responses categorised as “pain onset” (d4), “maximal pain” (d7), “beginning of pain resolution” (d11), and “recovery phase” (d23) for comparative longitudinal RNA sequencing. The generated datasets validate previous discoveries and reveal additional gene expression changes that warrant further validation with the goal to aid in development of drugs that prevent or reverse paclitaxel-induced peripheral neuropathy.

2. Data Description

2.1. Background and Summary

Peripheral neuropathy is a common side effect of chemotherapy characterised by paraesthesia (tingling), numbness, pain, temperature sensitivity, and motor weakness. Paclitaxel (Taxol) is one of the most widely used chemotherapeutic agents, which primarily affects the somatosensory neurons innervating the skin [1,2,3]. Pathological examinations suggest that intraepidermal unmyelinated axons are the first to degenerate upon paclitaxel treatment [4,5,6,7,8]. Thus, understanding the genetic mechanisms underlying the earliest manifestations of the disease will be essential to develop therapies that allow chemotherapy patients to complete cancer treatment without disruption, and prevent irreversible long-term symptoms.
Few studies have established expression profiles of chemotherapy-induced peripheral neuropathy. In one study, parallel gene expression profiles from dorsal root ganglion (DRG) neurons in mice were established following injection with the chemotherapeutic agents, oxaliplatin, vincristine, and cisplatin [9]. This comparative study revealed that only few genes were common among these data sets, suggesting that fundamental differences in the aetiology of chemotherapy-induced peripheral neuropathy (CIPN) must exist. This may not be surprising given the differences in the mechanisms of action for each of these chemotherapeutic agents, leading to potentially different off-target effects. In addition, the investigation of dorsal root ganglion (DRG) neurons may have obscured common upstream mechanisms. For instance, we previously showed that sensory neurons are secondarily affected by earlier epidermal damage, which promotes the degeneration of intraepidermal nerve endings in zebrafish, rats, and mice treated with paclitaxel. Epidermal keratinocytes are damaged due to increased reactive oxygen species formation and upregulation of matrix metalloproteinases, such as MMP-13, leading to extracellular matrix damage that ultimately affects the axons, leading to their degeneration [8]. Epidermal damage can be prevented when animals are treated with pharmacological MMP-13 inhibitors [10]. Therefore, the skin plays a crucial role in sensory axon homeostasis.
Existing genomic studies have focused on single time point analyses and single-cell types [11,12]. For instance, RNA sequencing was used to analyse blood samples of breast cancer survivors who suffered from long-term paclitaxel-induced peripheral neuropathy and these samples were compared to breast cancer survivors without neuropathy [13]. This study identified changes in mitochondrial genes that had been previously identified in preclinical CIPN models as differentially regulated, validating the importance of these models in studying human pathology [13]. Potentially, mitochondrial dysfunction might play a role in the deficiency of some patients to resolve their neuropathy symptoms. Despite these findings, longitudinal studies to detect changes in affected tissues over prolonged time periods have not been conducted, and currently, no data is available on gene expression changes prior to the onset of neuropathic symptoms. This information, however, will be critical to understanding the molecular gene expression networks involved in the onset, progression, and resolution of neuropathy.
To address this need, we performed a comprehensive RNAseq study using skin and DRG neuron samples of vehicle and paclitaxel-treated mice. We compared gene expression profiles according to pain profiles generated in these mice. Mice were injected four times every other day with either vehicle or paclitaxel and subsequently underwent a recovery period between Day 7 and Day 23. Tissues were collected and analysed during these time points, which were categorised as “pain onset” on Day 4, “maximal pain sensitivity” on Day 7, “beginning of pain resolution” on Day 11, and “post pain” on Day 23 (Figure 1). The generated data sets will be useful for the research community to further validate the genes implicated in paclitaxel-induced peripheral neuropathy.

2.2. Technical Validation

2.2.1. Differential Expression Analysis

Gene counts from STAR alignment were imported into the R environment using Rstudio (R Version 4.0.3; R studio version 1.3.1093; Peter Dalgaard; Frederiksberg, Denmark) and the package base version 4.0.3. The package dplyr version 1.0.5 was used to construct individual raw gene count matrices per tissue type, skin, and dorsal root ganglion neurons (DRGs), each consisting of 32 samples combined from both treatment groups. Normalised gene expression for both tissue data sets was calculated using the default parameters in the Deseq2 package version 1.3.1 (Michael Love; Chapel Hill, NC, USA).

2.2.2. Expression of Housekeeping and Validation Genes

Normalised gene counts were plotted using Prism 9 (GraphPad Prism Version 9.0.0; San Diego, CA, USA). Previously validated housekeeping genes for DRG and skin were plotted in addition to genes shown to be overexpressed during paclitaxel-induced peripheral neuropathy in either one. To validate DRG-specific housekeeping genes, we plotted TATA box binding protein (Tbp), RPTOR independent companion of MTOR, complex 2 (Rictor), and Ankyrin repeat domain-containing protein 27 (Ankrd27), as these are known to be unaffected in their expression levels during nerve injury and pain [14] (Figure 2a). Next, we plotted genes that have been established to be overexpressed in DRG neurons in the presence of paclitaxel, including Monocyte chemoattractant protein 1 (Ccl2). This gene mediates macrophage recruitment and promotes peripheral neuropathy in the presence of paclitaxel [15]. We further analysed the expression of Itgb1, which protects DRG neurons from paclitaxel-induced axon damage [16]. This revealed its upregulation, potentially a compensatory mechanism due to altered Itgb1 trafficking in DRG neurons of paclitaxel-treated mice [16]. Since our research showed that increased MMP-13 activity in the epidermis promotes the development of paclitaxel-induced peripheral neuropathy, we further determined whether Mmp13 expression also changes in DRG neurons, which revealed no significant increase in these neurons.
To validate skin-specific gene expression profiles, we choose to analyse the known keratinocyte housekeeping genes, TATA box binding protein (Tbp), Ribosomal protein large P0 (Rplp0), and Phosphoglycerate kinase 1 (Pgk1) [17,18], which as expected did not significantly vary in their expression levels regardless of treatment (Figure 2b). To validate genes that we expected to display a change in expression, we first analysed Mmp13, which displayed a paclitaxel-dependent significant increase in expression in the skin, in line with our previous findings [2,8]. Next, we analysed the expression of Tissue inhibitor of metalloproteinase-3 (Timp3), a known antagonist of MMP-13 [19], which displayed a downward regulation trend; however, this was not significant. Given the altered cell adhesion in the skin of paclitaxel-treated animals induced by increased MMP activity, we also analysed the tight junction protein, Claudin 22 (Cldn22) [2,8]. Tight junctions, for instance, have been shown to be decreased upon MMP-13 activation in intestinal epithelia upon LPS stimulation [20]. Consistent with this finding, we detected a significant decrease in Cldn22 expression in the skin [21] following paclitaxel treatment. This confirms that expected gene expression trends can be detected in our gene expression data sets, thus validating our approach. These data will be useful for the research community for further analysis. Table 1 provides a description for each of the individual raw files, their corresponding sample ID, and information about the time points, treatment, and tissue type.

3. Methods

3.1. Animals

All animal procedures were approved by the University of New England Institutional Animal Care and Use Committee. Adult male C57BL6/J mice (JAX) weighing 20–25 g were purchased from the Jackson Laboratory. Upon arrival, mice were housed 4/cage and allowed to acclimate to the facility for 7 days. All animals were kept on a 12 h light/dark cycle with ad libitum access to food and water.

3.2. Paclitaxel Treatment

The experimental design, time points, and downstream analyses are depicted in Figure 1a. Paclitaxel was administered on Days 0, 2, 4, and 6. Paclitaxel (Sigma-Aldrich, St. Louis, USA) was dissolved in (1:1 Cremophor: Ethanol) and further diluted in 0.9% NaCl to make a final concentration of 0.4 mg/mL. Six-week-old mice were injected intraperitoneally with either vehicle or paclitaxel at a volume of 10 mL/kg bodyweight to make a final concentration of 4 mg/kg (cumulative 16 mg/kg).

3.3. Analysis of Tactile Thresholds

Tactile allodynia was quantified (in the mornings prior to injections when assessed on injection days) by measuring the hind paw withdrawal threshold to von Frey filament stimulation, using the up–down method as previously reported [22]. Results are shown in Figure 1b. Animals were placed in a clear Plexiglas chamber and allowed to habituate for ~60 min. Touch-Test filaments (North Coast Medical, Morgan Hill, CA, USA) were used for all testing. The filament range was: 2.44, 2.83, 3.22, 3.61, 4.08, 4.31, 4.56, starting with 3.61. Withdrawal thresholds were determined by sequentially increasing and decreasing the stimulus intensity (“up and down” method). This up–down method was stopped four measures after the first positive response. The response threshold was subsequently analysed by using a Dixon nonparametric test and expressed as the paw withdrawal threshold in gram force values [23]. A clear paw withdrawal, shaking, or licking was considered as a positive or painful response.

3.4. Tissue Collection

Animals were exsanguinated by intracardial perfusion of 100 mL ice-cold heparinised phosphate-buffered saline (PBS). Tissues were then immediately harvested on ice before flash freezing and stored at −80 °C. DRG neurons were harvested from left lumbar segments L3, L4, and L5. Plantar skin was taken from the left hind paw. Tissues were harvested at time points consistent with the development of paclitaxel-induced neuropathy in rodents characterised as “pain onset” on Day 4, “maximal pain sensitivity” on Day 7, “beginning of pain resolution” on Day 11, and “post pain” on Day 23 [10].

3.5. RNA Extraction and Quality Control

DRG were homogenised using a glass-Teflon homogeniser on ice for 2 min with 350 µL buffer RLT plus 2-Mercaptoethanol (Qiagen, Boston, MA, USA). RNA was extracted using Qiagen RNeasy microcolumns according to the manufacturer’s protocol and eluted with 30 µL ultrapure RNase-free water. Skin was homogenised with a rotor–stator homogeniser for 2 min on ice with TRIzol reagent (Thermo, St. Louis, USA) before phase separation with phase lock gel microcentrifuge tubes (Eppendorf, Framington, MA, USA) to prevent guanidine salt contamination, according to the manufacturer’s protocol. RNA purity was confirmed using a nanodrop UV spectrometer (260/280 ratio > 2), and RNA quality was assessed with the Agilent 2100 Bioanalyzer (RIN > 7).

3.6. cDNA Library Preparation and Sequencing

Ribosomal RNA-depleted cDNA libraries were prepared using KAPA RNA HyperPrep kits with RiboErase (KAPA Biosystems) according to the manufacturer’s protocol. Libraries were sequenced with an Illumina Hi-seq 3000 (single-end, 1 × 100 bp).

3.7. Sequence Alignment, Counts, and Quality Control

Remaining adapters and poor-quality sequences were trimmed with Trimmomatic (version 0.38) [24] with the following parameters: ILLUMINACLIP:TruSeq3-SE.fa:2:30:10 SLIDINGWINDOW:4:15 MINLEN:15. Sequence quality was then verified using FastQC and summarised with MultiQC (version 1.9) (Figure 1c). Using STAR [25] (version 2.7.5a_2020-06-19), sequences were aligned to mouse genome assembly m38 with Ensembl version 100 annotations using the quantMode GeneCounts parameter. Mapping quality was verified using MultiQC (Figure 1d,e).

3.8. Code Availability

The R code used to analyse and process the raw count data from skin and DRG samples (Table 1) using DESeq2 package version on R version 4.0.3 is publicly available at https://github.com/acirrincio/PaclitaxelRNAseq/ (accessed on 25 May 2022).

Author Contributions

A.M.C. and C.A.R. contributed equally to the data analysis. B.J.H. and S.R. jointly supervised this work. S.R. wrote the manuscript. All authors edited the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This work was sponsored by NIH CTSI Miami Pilot-FY19-05 and NIH/NCI 1R01CA215973 funds awarded to S.R., and NIGMS P20GM103643 awarded to B.J.H.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All raw and normalised gene count matrices produced by DESeq2 were deposited in the NCBI Gene Expression Omnibus (GEO) database under accession number GSE185084. This GEO entry includes links to the raw count data in .txt.bz2 format. Each processed matrix data file of raw or normalised counts contains a header row corresponding to the sample followed by 52,391 rows, each corresponding to a unique transcript.

Acknowledgments

We thank the University of New England Behaviour Core staff, Denise Giuvelis and Kerry Tucker, for their assistance. We further thank Ramaz Geguchadze and Eliza Grlickova-Duzevik at the University of New England for assisting with the isolation of the DRG neurons.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Boyette-Davis, J.A.; Walters, E.T.; Dougherty, P.M. Mechanisms involved in the development of chemotherapy-induced neuropathy. Pain Manag. 2015, 5, 285–296. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Flatters, S.J.L.; Dougherty, P.M.; Colvin, L.A. Clinical and preclinical perspectives on Chemotherapy-Induced Peripheral Neuropathy (CIPN): A narrative review. Br. J. Anaesth. 2017, 119, 737–749. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Gornstein, E.; Schwarz, T.L. The paradox of paclitaxel neurotoxicity: Mechanisms and unanswered questions. Neuropharmacology 2014, 76 Pt A, 175–183. [Google Scholar] [CrossRef]
  4. Boyette-Davis, J.; Xin, W.; Zhang, H.; Dougherty, P.M. Intraepidermal nerve fiber loss corresponds to the development of taxol-induced hyperalgesia and can be prevented by treatment with minocycline. Pain 2011, 152, 308–313. [Google Scholar] [CrossRef]
  5. Bennett, G.J.; Liu, G.K.; Xiao, W.H.; Jin, H.W.; Siau, C. Terminal arbor degeneration—A novel lesion produced by the antineoplastic agent paclitaxel. Eur. J. Neurosci. 2011, 33, 1667–1676. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Jin, H.W.; Flatters, S.J.; Xiao, W.H.; Mulhern, H.L.; Bennett, G.J. Prevention of paclitaxel-evoked painful peripheral neuropathy by acetyl-L-carnitine: Effects on axonal mitochondria, sensory nerve fiber terminal arbors, and cutaneous Langerhans cells. Exp. Neurol. 2008, 210, 229–237. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Ko, M.H.; Hu, M.E.; Hsieh, Y.L.; Lan, C.T.; Tseng, T.J. Peptidergic intraepidermal nerve fibers in the skin contribute to the neuropathic pain in paclitaxel-induced peripheral neuropathy. Neuropeptides 2014, 48, 109–117. [Google Scholar] [CrossRef]
  8. Lisse, T.S.; Middleton, L.J.; Pellegrini, A.D.; Martin, P.B.; Spaulding, E.L.; Lopes, O.; Brochu, E.A.; Carter, E.V.; Waldron, A.; Rieger, S. Paclitaxel-induced epithelial damage and ectopic MMP-13 expression promotes neurotoxicity in zebrafish. Proc. Natl. Acad. Sci. USA 2016, 113, 2189–2198. [Google Scholar] [CrossRef] [Green Version]
  9. Starobova, H.; Mueller, A.; Deuis, J.R.; Carter, D.A.; Vetter, I. Inflammatory and Neuropathic Gene Expression Signatures of Chemotherapy-Induced Neuropathy Induced by Vincristine, Cisplatin, and Oxaliplatin in C57BL/6J Mice. J. Pain 2020, 21, 182–194. [Google Scholar] [CrossRef]
  10. Cirrincione, A.M.; Pellegrini, A.D.; Dominy, J.R.; Benjamin, M.E.; Utkina-Sosunova, I.; Lotti, F.; Jergova, S.; Sagen, J.; Rieger, S. Paclitaxel-induced peripheral neuropathy is caused by epidermal ROS and mitochondrial damage through conserved MMP-13 activation. Sci. Rep. 2020, 10, 3970. [Google Scholar] [CrossRef] [Green Version]
  11. Sivanesan, E.; Stephens, K.E.; Huang, Q.; Chen, Z.; Ford, N.C.; Duan, W.; He, S.Q.; Gao, X.; Linderoth, B.; Raja, S.N.; et al. Spinal cord stimulation prevents paclitaxel-induced mechanical and cold hypersensitivity and modulates spinal gene expression in rats. Pain Rep. 2019, 4, 785. [Google Scholar] [CrossRef] [PubMed]
  12. Li, Y.; Yin, C.; Liu, B.; Nie, H.; Wang, J.; Zeng, D.; Chen, R.; He, X.; Fang, J.; Du, J.; et al. Transcriptome profiling of long noncoding RNAs and mRNAs in spinal cord of a rat model of paclitaxel-induced peripheral neuropathy identifies potential mechanisms mediating neuroinflammation and pain. J. Neuroinflammation 2021, 18, 48. [Google Scholar] [CrossRef] [PubMed]
  13. Kober, K.M.; Olshen, A.; Conley, Y.P.; Schumacher, M.; Topp, K.; Smoot, B.; Mazor, M.; Chesney, M.; Hammer, M.; Paul, S.M.; et al. Expression of mitochondrial dysfunction-related genes and pathways in paclitaxel-induced peripheral neuropathy in breast cancer survivors. Mol. Pain 2018, 14, 1744806918816462. [Google Scholar] [CrossRef] [PubMed]
  14. Wang, Y.; Shan, Q.; Meng, Y.; Pan, J.; Yi, S. Mrpl10 and Tbp Are Suitable Reference Genes for Peripheral Nerve Crush Injury. Int. J. Mol. Sci. 2017, 18, 263. [Google Scholar] [CrossRef] [Green Version]
  15. Zhang, H.; Li, Y.; de Carvalho-Barbosa, M.; Kavelaars, A.; Heijnen, C.J.; Albrecht, P.J.; Dougherty, P.M. Dorsal Root Ganglion Infiltration by Macrophages Contributes to Paclitaxel Chemotherapy-Induced Peripheral Neuropathy. J. Pain 2016, 17, 775–786. [Google Scholar] [CrossRef] [Green Version]
  16. Shin, G.J.; Pero, M.E.; Hammond, L.A.; Burgos, A.; Kumar, A.; Galindo, S.E.; Lucas, T.; Bartolini, F.; Grueber, W.B. Integrins protect sensory neurons in models of paclitaxel-induced peripheral sensory neuropathy. Proc. Natl. Acad. Sci. USA 2021, 118, 2006050118. [Google Scholar] [CrossRef]
  17. Minner, F.; Poumay, Y. Candidate housekeeping genes require evaluation before their selection for studies of human epidermal keratinocytes. J. Investig. Dermatol. 2009, 129, 770–773. [Google Scholar] [CrossRef] [Green Version]
  18. Riemer, A.B.; Keskin, D.B.; Reinherz, E.L. Identification and validation of reference genes for expression studies in human keratinocyte cell lines treated with and without interferon-γ—A method for qRT-PCR reference gene determination. Exp. Derm. 2012, 21, 625–629. [Google Scholar] [CrossRef] [Green Version]
  19. Carreca, A.P.; Pravatà, V.M.; Markham, M.; Bonelli, S.; Murphy, G.; Nagase, H.; Troeberg, L.; Scilabra, S.D. TIMP-3 facilitates binding of target metalloproteinases to the endocytic receptor LRP-1 and promotes scavenging of MMP-1. Sci. Rep. 2020, 10, 12067. [Google Scholar] [CrossRef]
  20. Vandenbroucke, R.E.; Dejonckheere, E.; Van Hauwermeiren, F.; Lodens, S.; De Rycke, R.; Van Wonterghem, E.; Staes, A.; Gevaert, K.; López-Otin, C.; Libert, C. Matrix metalloproteinase 13 modulates intestinal epithelial barrier integrity in inflammatory diseases by activating TNF. EMBO Mol. Med. 2013, 5, 932–948. [Google Scholar] [CrossRef] [Green Version]
  21. Klymenko, T.; Gu, Q.; Herbert, I.; Stevenson, A.; Iliev, V.; Watkins, G.; Pollock, C.; Bhatia, R.; Cuschieri, K.; Herzyk, P.; et al. RNA-Seq Analysis of Differentiated Keratinocytes Reveals a Massive Response to Late Events during Human Papillomavirus 16 Infection, Including Loss of Epithelial Barrier Function. J. Virol. 2017, 91, 01001–01017. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Dixon, W.J. Efficient analysis of experimental observations. Annu. Rev. Pharm. Toxicol. 1980, 20, 441–462. [Google Scholar] [CrossRef] [PubMed]
  23. Chaplan, S.R.; Bach, F.W.; Pogrel, J.W.; Chung, J.M.; Yaksh, T.L. Quantitative assessment of tactile allodynia in the rat paw. J. Neurosci. Methods 1994, 53, 55–63. [Google Scholar] [CrossRef]
  24. Bolger, A.M.; Lohse, M.; Usadel, B. Trimmomatic: A flexible trimmer for Illumina sequence data. Bioinformatics 2014, 30, 2114–2120. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Dobin, A.; Davis, C.A.; Schlesinger, F.; Drenkow, J.; Zaleski, C.; Jha, S.; Batut, P.; Chaisson, M.; Gingeras, T.R. STAR: Ultrafast universal RNA-seq aligner. Bioinformatics 2013, 29, 15–21. [Google Scholar] [CrossRef]
Figure 1. Experimental overview of longitudinal RNAseq analysis comparing dorsal root ganglion neurons and skin from paclitaxel- and vehicle-treated mice. (a) Summary of experimental workflow from treatment and sample collection to data trimming and analysis. (b) Behavioural analysis using von Frey filaments to assess the touch response showing the development of paclitaxel-induced neuropathy in mice evident by reduced pain threshold at 4, 7, and 11 days with a subsequent recovery (n = 4 per treatment group per time point). (c) Read-quality scores for DRG neurons (left) and skin (right). (d,e) Total and uniquely mapped reads among DRG (left) and skin (right) samples from paclitaxel- and vehicle-treated animals. Abbreviations: VEH: vehicle; PCTX: Paclitaxel.
Figure 1. Experimental overview of longitudinal RNAseq analysis comparing dorsal root ganglion neurons and skin from paclitaxel- and vehicle-treated mice. (a) Summary of experimental workflow from treatment and sample collection to data trimming and analysis. (b) Behavioural analysis using von Frey filaments to assess the touch response showing the development of paclitaxel-induced neuropathy in mice evident by reduced pain threshold at 4, 7, and 11 days with a subsequent recovery (n = 4 per treatment group per time point). (c) Read-quality scores for DRG neurons (left) and skin (right). (d,e) Total and uniquely mapped reads among DRG (left) and skin (right) samples from paclitaxel- and vehicle-treated animals. Abbreviations: VEH: vehicle; PCTX: Paclitaxel.
Data 07 00072 g001
Figure 2. Technical validation of genes. Expression graphs of normalised gene counts showing housekeeping genes (no expected expression changes) and validation genes (expected expression changes) for (a) DRG neurons and (b) skin. (a) DRG validation genes include Ccl2 and Itgb1, known to be overexpressed in DRG neurons in the presence of paclitaxel, and Mmp13, only known to play a skin-specific role in paclitaxel-induced peripheral neuropathy. (b) Known housekeeping genes in skin include Tbp, RplpO, and Pgk1. Validation genes in skin include Mmp13, Timp3, a known inhibitor of MMP-13, and Cldn22, a tight junction protein in the epidermis. p-values from a two-sided t-test are represented above brackets for each time point. Abbreviations: VEH: vehicle; PCTX: Paclitaxel.
Figure 2. Technical validation of genes. Expression graphs of normalised gene counts showing housekeeping genes (no expected expression changes) and validation genes (expected expression changes) for (a) DRG neurons and (b) skin. (a) DRG validation genes include Ccl2 and Itgb1, known to be overexpressed in DRG neurons in the presence of paclitaxel, and Mmp13, only known to play a skin-specific role in paclitaxel-induced peripheral neuropathy. (b) Known housekeeping genes in skin include Tbp, RplpO, and Pgk1. Validation genes in skin include Mmp13, Timp3, a known inhibitor of MMP-13, and Cldn22, a tight junction protein in the epidermis. p-values from a two-sided t-test are represented above brackets for each time point. Abbreviations: VEH: vehicle; PCTX: Paclitaxel.
Data 07 00072 g002
Table 1. Raw data descriptor.
Table 1. Raw data descriptor.
File nameSample NameOrganismStrainAgeSexTissueTreatment_Timepoint
sRieger_RiegerHarrisonTaxol11202019_201926297-01_S_1_1SKIN05Mus musculusC57B6/J6 weeks maleSkinVehicle_day4
sRieger_RiegerHarrisonTaxol11202019_201926297-01_S_2_1SKIN05Mus musculusC57B6/J6 weeks maleSkinVehicle_day4
sRieger_RiegerHarrisonTaxol11202019_201926301-01_S_1_1SKIN07Mus musculusC57B6/J6 weeks maleSkinVehicle_day4
sRieger_RiegerHarrisonTaxol11202019_201926301-01_S_2_1SKIN07Mus musculusC57B6/J6 weeks maleSkinVehicle_day4
sRieger_RiegerHarrisonTaxol11202019_201926344-01_S_1_1SKIN06Mus musculusC57B6/J6 weeks maleSkinVehicle_day4
sRieger_RiegerHarrisonTaxol11202019_201926344-01_S_2_1SKIN06Mus musculusC57B6/J6 weeks maleSkinVehicle_day4
sRieger_RiegerHarrisonTaxol11202019_201926350-01_S_1_1SKIN08Mus musculusC57B6/J6 weeks maleSkinVehicle_day4
sRieger_RiegerHarrisonTaxol11202019_201926350-01_S_2_1SKIN08Mus musculusC57B6/J6 weeks maleSkinVehicle_day4
sRieger_RiegerHarrisonTaxol11202019_201926321-01_S_2_1SKIN19Mus musculusC57B6/J6 weeks maleSkinVehicle_day7
sRieger_RiegerHarrisonTaxol11202019_201926335-01_S_1_1SKIN18Mus musculusC57B6/J6 weeks maleSkinVehicle_day7
sRieger_RiegerHarrisonTaxol11202019_201926335-01_S_2_1SKIN18Mus musculusC57B6/J6 weeks maleSkinVehicle_day7
sRieger_RiegerHarrisonTaxol11202019_201926349-01_S_1_1SKIN20Mus musculusC57B6/J6 weeks maleSkinVehicle_day7
sRieger_RiegerHarrisonTaxol11202019_201926355-01_S_1_1SKIN17Mus musculusC57B6/J6 weeks maleSkinVehicle_day7
sRieger_RiegerHarrisonTaxol11202019_201926355-01_S_2_1SKIN17Mus musculusC57B6/J6 weeks maleSkinVehicle_day7
sRieger_RiegerHarrisonTaxol11202019_201926324-01_S_1_1SKIN35Mus musculusC57B6/J6 weeks maleSkinVehicle_day11
sRieger_RiegerHarrisonTaxol11202019_201926329-01_S_1_1SKIN33Mus musculusC57B6/J6 weeks maleSkinVehicle_day11
sRieger_RiegerHarrisonTaxol11202019_201926329-01_S_2_1SKIN33Mus musculusC57B6/J6 weeks maleSkinVehicle_day11
sRieger_RiegerHarrisonTaxol11202019_201926337-01_S_1_1SKIN34Mus musculusC57B6/J6 weeks maleSkinVehicle_day11
sRieger_RiegerHarrisonTaxol11202019_201926337-01_S_2_1SKIN34Mus musculusC57B6/J6 weeks maleSkinVehicle_day11
sRieger_RiegerHarrisonTaxol11202019_201926346-01_S_1_1SKIN36Mus musculusC57B6/J6 weeks maleSkinVehicle_day11
sRieger_RiegerHarrisonTaxol11202019_201926346-01_S_2_1SKIN36Mus musculusC57B6/J6 weeks maleSkinVehicle_day11
sRieger_RiegerHarrisonTaxol11202019_201926295-01_S_1_1SKIN52Mus musculusC57B6/J6 weeks maleSkinVehicle_day23
sRieger_RiegerHarrisonTaxol11202019_201926295-01_S_2_1SKIN52Mus musculusC57B6/J6 weeks maleSkinVehicle_day23
sRieger_RiegerHarrisonTaxol11202019_201926304-01_S_1_1SKIN49Mus musculusC57B6/J6 weeks maleSkinVehicle_day23
sRieger_RiegerHarrisonTaxol11202019_201926304-01_S_2_1SKIN49Mus musculusC57B6/J6 weeks maleSkinVehicle_day23
sRieger_RiegerHarrisonTaxol11202019_201926328-01_S_1_1SKIN50Mus musculusC57B6/J6 weeks maleSkinVehicle_day23
sRieger_RiegerHarrisonTaxol11202019_201926328-01_S_2_1SKIN50Mus musculusC57B6/J6 weeks maleSkinVehicle_day23
sRieger_RiegerHarrisonTaxol11202019_201926354-01_S_1_1SKIN51Mus musculusC57B6/J6 weeks maleSkinVehicle_day23
sRieger_RiegerHarrisonTaxol11202019_201926354-01_S_2_1SKIN51Mus musculusC57B6/J6 weeks maleSkinVehicle_day23
sRieger_RiegerHarrisonTaxol11202019_201926296-01_S_1_1SKIN09Mus musculusC57B6/J6 weeks maleSkinPaclitaxel_day4
sRieger_RiegerHarrisonTaxol11202019_201926296-01_S_2_1SKIN09Mus musculusC57B6/J6 weeks maleSkinPaclitaxel_day4
sRieger_RiegerHarrisonTaxol11202019_201926306-01_S_1_1SKIN10Mus musculusC57B6/J6 weeks maleSkinPaclitaxel_day4
sRieger_RiegerHarrisonTaxol11202019_201926306-01_S_2_1SKIN10Mus musculusC57B6/J6 weeks maleSkinPaclitaxel_day4
sRieger_RiegerHarrisonTaxol11202019_201926326-01_S_1_1SKIN12Mus musculusC57B6/J6 weeks maleSkinPaclitaxel_day4
sRieger_RiegerHarrisonTaxol11202019_201926339-01_S_1_1SKIN11Mus musculusC57B6/J6 weeks maleSkinPaclitaxel_day4
sRieger_RiegerHarrisonTaxol11202019_201926339-01_S_2_1SKIN11Mus musculusC57B6/J6 weeks maleSkinPaclitaxel_day4
sRieger_RiegerHarrisonTaxol11202019_201926294-01_S_1_1SKIN25Mus musculusC57B6/J6 weeks maleSkinPaclitaxel_day7
sRieger_RiegerHarrisonTaxol11202019_201926294-01_S_2_1SKIN25Mus musculusC57B6/J6 weeks maleSkinPaclitaxel_day7
sRieger_RiegerHarrisonTaxol11202019_201926334-01_S_1_1SKIN26Mus musculusC57B6/J6 weeks maleSkinPaclitaxel_day7
sRieger_RiegerHarrisonTaxol11202019_201926334-01_S_2_1SKIN26Mus musculusC57B6/J6 weeks maleSkinPaclitaxel_day7
sRieger_RiegerHarrisonTaxol11202019_201926336-01_S_1_1SKIN27Mus musculusC57B6/J6 weeks maleSkinPaclitaxel_day7
sRieger_RiegerHarrisonTaxol11202019_201926336-01_S_2_1SKIN27Mus musculusC57B6/J6 weeks maleSkinPaclitaxel_day7
sRieger_RiegerHarrisonTaxol11202019_201926338-01_S_1_1SKIN28Mus musculusC57B6/J6 weeks maleSkinPaclitaxel_day7
sRieger_RiegerHarrisonTaxol11202019_201926338-01_S_2_1SKIN28Mus musculusC57B6/J6 weeks maleSkinPaclitaxel_day7
sRieger_RiegerHarrisonTaxol11202019_201926298-01_S_1_1SKIN43Mus musculusC57B6/J6 weeks maleSkinPaclitaxel_day11
sRieger_RiegerHarrisonTaxol11202019_201926298-01_S_2_1SKIN43Mus musculusC57B6/J6 weeks maleSkinPaclitaxel_day11
sRieger_RiegerHarrisonTaxol11202019_201926310-01_S_1_1SKIN42Mus musculusC57B6/J6 weeks maleSkinPaclitaxel_day11
sRieger_RiegerHarrisonTaxol11202019_201926310-01_S_2_1SKIN42Mus musculusC57B6/J6 weeks maleSkinPaclitaxel_day11
sRieger_RiegerHarrisonTaxol11202019_201926331-01_S_1_1SKIN44Mus musculusC57B6/J6 weeks maleSkinPaclitaxel_day11
sRieger_RiegerHarrisonTaxol11202019_201926331-01_S_2_1SKIN44Mus musculusC57B6/J6 weeks maleSkinPaclitaxel_day11
sRieger_RiegerHarrisonTaxol11202019_201926332-01_S_1_1SKIN41Mus musculusC57B6/J6 weeks maleSkinPaclitaxel_day11
sRieger_RiegerHarrisonTaxol11202019_201926299-01_S_1_1SKIN60Mus musculusC57B6/J6 weeks maleSkinPaclitaxel_day23
sRieger_RiegerHarrisonTaxol11202019_201926299-01_S_2_1SKIN60Mus musculusC57B6/J6 weeks maleSkinPaclitaxel_day23
sRieger_RiegerHarrisonTaxol11202019_201926313-01_S_1_1SKIN57Mus musculusC57B6/J6 weeks maleSkinPaclitaxel_day23
sRieger_RiegerHarrisonTaxol11202019_201926313-01_S_2_1SKIN57Mus musculusC57B6/J6 weeks maleSkinPaclitaxel_day23
sRieger_RiegerHarrisonTaxol11202019_201926319-01_S_2_1SKIN58Mus musculusC57B6/J6 weeks maleSkinPaclitaxel_day23
sRieger_RiegerHarrisonTaxol11202019_201926343-01_S_1_1SKIN59Mus musculusC57B6/J6 weeks maleSkinPaclitaxel_day23
sRieger_RiegerHarrisonTaxol11202019_201926343-01_S_2_1SKIN59Mus musculusC57B6/J6 weeks maleSkinPaclitaxel_day23
sRieger_RiegerHarrisonTaxol11202019_201926362-01_S_1_1DRG08Mus musculusC57B6/J6 weeks maleDRGVehicle_day4
sRieger_RiegerHarrisonTaxol11202019_201926362-01_S_2_1DRG08Mus musculusC57B6/J6 weeks maleDRGVehicle_day4
sRieger_RiegerHarrisonTaxol11202019_201926366-01_S_1_1DRG07Mus musculusC57B6/J6 weeks maleDRGVehicle_day4
sRieger_RiegerHarrisonTaxol11202019_201926366-01_S_2_1DRG07Mus musculusC57B6/J6 weeks maleDRGVehicle_day4
sRieger_RiegerHarrisonTaxol11202019_201926387-01_S_1_1DRG06Mus musculusC57B6/J6 weeks maleDRGVehicle_day4
sRieger_RiegerHarrisonTaxol11202019_201926387-01_S_2_1DRG06Mus musculusC57B6/J6 weeks maleDRGVehicle_day4
sRieger_RiegerHarrisonTaxol11202019_201926413-01_S_1_1DRG05Mus musculusC57B6/J6 weeks maleDRGVehicle_day4
sRieger_RiegerHarrisonTaxol11202019_201926413-01_S_2_1DRG05Mus musculusC57B6/J6 weeks maleDRGVehicle_day4
sRieger_RiegerHarrisonTaxol11202019_201926361-01_S_1_1DRG17Mus musculusC57B6/J6 weeks maleDRGVehicle_day7
sRieger_RiegerHarrisonTaxol11202019_201926361-01_S_2_1DRG17Mus musculusC57B6/J6 weeks maleDRGVehicle_day7
sRieger_RiegerHarrisonTaxol11202019_201926396-01_S_1_1DRG20Mus musculusC57B6/J6 weeks maleDRGVehicle_day7
sRieger_RiegerHarrisonTaxol11202019_201926396-01_S_2_1DRG20Mus musculusC57B6/J6 weeks maleDRGVehicle_day7
sRieger_RiegerHarrisonTaxol11202019_201926400-01_S_1_1DRG19Mus musculusC57B6/J6 weeks maleDRGVehicle_day7
sRieger_RiegerHarrisonTaxol11202019_201926414-01_S_1_1DRG18Mus musculusC57B6/J6 weeks maleDRGVehicle_day7
sRieger_RiegerHarrisonTaxol11202019_201926414-01_S_2_1DRG18Mus musculusC57B6/J6 weeks maleDRGVehicle_day7
sRieger_RiegerHarrisonTaxol11202019_201926358-01_S_1_1DRG34Mus musculusC57B6/J6 weeks maleDRGVehicle_day11
sRieger_RiegerHarrisonTaxol11202019_201926358-01_S_2_1DRG34Mus musculusC57B6/J6 weeks maleDRGVehicle_day11
sRieger_RiegerHarrisonTaxol11202019_201926364-01_S_1_1DRG35Mus musculusC57B6/J6 weeks maleDRGVehicle_day11
sRieger_RiegerHarrisonTaxol11202019_201926364-01_S_2_1DRG35Mus musculusC57B6/J6 weeks maleDRGVehicle_day11
sRieger_RiegerHarrisonTaxol11202019_201926408-01_S_1_1DRG36Mus musculusC57B6/J6 weeks maleDRGVehicle_day11
sRieger_RiegerHarrisonTaxol11202019_201926408-01_S_2_1DRG36Mus musculusC57B6/J6 weeks maleDRGVehicle_day11
sRieger_RiegerHarrisonTaxol11202019_201926411-01_S_1_1DRG33Mus musculusC57B6/J6 weeks maleDRGVehicle_day11
sRieger_RiegerHarrisonTaxol11202019_201926411-01_S_2_1DRG33Mus musculusC57B6/J6 weeks maleDRGVehicle_day11
sRieger_RiegerHarrisonTaxol11202019_201926370-01_S_1_1DRG50Mus musculusC57B6/J6 weeks maleDRGVehicle_day23
sRieger_RiegerHarrisonTaxol11202019_201926370-01_S_2_1DRG50Mus musculusC57B6/J6 weeks maleDRGVehicle_day23
sRieger_RiegerHarrisonTaxol11202019_201926401-01_S_1_1DRG49Mus musculusC57B6/J6 weeks maleDRGVehicle_day23
sRieger_RiegerHarrisonTaxol11202019_201926401-01_S_2_1DRG49Mus musculusC57B6/J6 weeks maleDRGVehicle_day23
sRieger_RiegerHarrisonTaxol11202019_201926403-01_S_1_1DRG52Mus musculusC57B6/J6 weeks maleDRGVehicle_day23
sRieger_RiegerHarrisonTaxol11202019_201926403-01_S_2_1DRG52Mus musculusC57B6/J6 weeks maleDRGVehicle_day23
sRieger_RiegerHarrisonTaxol11202019_201926405-01_S_1_1DRG51Mus musculusC57B6/J6 weeks maleDRGVehicle_day23
sRieger_RiegerHarrisonTaxol11202019_201926405-01_S_2_1DRG51Mus musculusC57B6/J6 weeks maleDRGVehicle_day23
sRieger_RiegerHarrisonTaxol11202019_201926356-01_S_1_1DRG09Mus musculusC57B6/J6 weeks maleDRGPaclitaxel_day4
sRieger_RiegerHarrisonTaxol11202019_201926356-01_S_2_1DRG09Mus musculusC57B6/J6 weeks maleDRGPaclitaxel_day4
sRieger_RiegerHarrisonTaxol11202019_201926399-01_S_1_1DRG12Mus musculusC57B6/J6 weeks maleDRGPaclitaxel_day4
sRieger_RiegerHarrisonTaxol11202019_201926399-01_S_2_1DRG12Mus musculusC57B6/J6 weeks maleDRGPaclitaxel_day4
sRieger_RiegerHarrisonTaxol11202019_201926402-01_S_1_1DRG10Mus musculusC57B6/J6 weeks maleDRGPaclitaxel_day4
sRieger_RiegerHarrisonTaxol11202019_201926402-01_S_2_1DRG10Mus musculusC57B6/J6 weeks maleDRGPaclitaxel_day4
sRieger_RiegerHarrisonTaxol11202019_201926410-01_S_1_1DRG11Mus musculusC57B6/J6 weeks maleDRGPaclitaxel_day4
sRieger_RiegerHarrisonTaxol11202019_201926410-01_S_2_1DRG11Mus musculusC57B6/J6 weeks maleDRGPaclitaxel_day4
sRieger_RiegerHarrisonTaxol11202019_201926360-01_S_1_1DRG28Mus musculusC57B6/J6 weeks maleDRGPaclitaxel_day7
sRieger_RiegerHarrisonTaxol11202019_201926360-01_S_2_1DRG28Mus musculusC57B6/J6 weeks maleDRGPaclitaxel_day7
sRieger_RiegerHarrisonTaxol11202019_201926367-01_S_1_1DRG27Mus musculusC57B6/J6 weeks maleDRGPaclitaxel_day7
sRieger_RiegerHarrisonTaxol11202019_201926367-01_S_2_1DRG27Mus musculusC57B6/J6 weeks maleDRGPaclitaxel_day7
sRieger_RiegerHarrisonTaxol11202019_201926415-01_S_1_1DRG26Mus musculusC57B6/J6 weeks maleDRGPaclitaxel_day7
sRieger_RiegerHarrisonTaxol11202019_201926415-01_S_2_1DRG26Mus musculusC57B6/J6 weeks maleDRGPaclitaxel_day7
sRieger_RiegerHarrisonTaxol11202019_201926417-01_S_1_1DRG25Mus musculusC57B6/J6 weeks maleDRGPaclitaxel_day7
sRieger_RiegerHarrisonTaxol11202019_201926417-01_S_2_1DRG25Mus musculusC57B6/J6 weeks maleDRGPaclitaxel_day7
sRieger_RiegerHarrisonTaxol11202019_201926369-01_S_1_1DRG43Mus musculusC57B6/J6 weeks maleDRGPaclitaxel_day11
sRieger_RiegerHarrisonTaxol11202019_201926369-01_S_2_1DRG43Mus musculusC57B6/J6 weeks maleDRGPaclitaxel_day11
sRieger_RiegerHarrisonTaxol11202019_201926375-01_S_1_1DRG42Mus musculusC57B6/J6 weeks maleDRGPaclitaxel_day11
sRieger_RiegerHarrisonTaxol11202019_201926375-01_S_2_1DRG42Mus musculusC57B6/J6 weeks maleDRGPaclitaxel_day11
sRieger_RiegerHarrisonTaxol11202019_201926382-01_S_1_1DRG41Mus musculusC57B6/J6 weeks maleDRGPaclitaxel_day11
sRieger_RiegerHarrisonTaxol11202019_201926382-01_S_2_1DRG41Mus musculusC57B6/J6 weeks maleDRGPaclitaxel_day11
sRieger_RiegerHarrisonTaxol11202019_201926386-01_S_1_1DRG44Mus musculusC57B6/J6 weeks maleDRGPaclitaxel_day11
sRieger_RiegerHarrisonTaxol11202019_201926365-01_S_1_1DRG60Mus musculusC57B6/J6 weeks maleDRGPaclitaxel_day23
sRieger_RiegerHarrisonTaxol11202019_201926365-01_S_2_1DRG60Mus musculusC57B6/J6 weeks maleDRGPaclitaxel_day23
sRieger_RiegerHarrisonTaxol11202019_201926380-01_S_1_1DRG59Mus musculusC57B6/J6 weeks maleDRGPaclitaxel_day23
sRieger_RiegerHarrisonTaxol11202019_201926380-01_S_2_1DRG59Mus musculusC57B6/J6 weeks maleDRGPaclitaxel_day23
sRieger_RiegerHarrisonTaxol11202019_201926381-01_S_1_1DRG58Mus musculusC57B6/J6 weeks maleDRGPaclitaxel_day23
sRieger_RiegerHarrisonTaxol11202019_201926381-01_S_2_1DRG58Mus musculusC57B6/J6 weeks maleDRGPaclitaxel_day23
sRieger_RiegerHarrisonTaxol11202019_201926393-01_S_1_1DRG57Mus musculusC57B6/J6 weeks maleDRGPaclitaxel_day23
sRieger_RiegerHarrisonTaxol11202019_201926393-01_S_2_1DRG57Mus musculusC57B6/J6 weeks maleDRGPaclitaxel_day23
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Cirrincione, A.M.; Reimonn, C.A.; Harrison, B.J.; Rieger, S. Longitudinal RNA Sequencing of Skin and DRG Neurons in Mice with Paclitaxel-Induced Peripheral Neuropathy. Data 2022, 7, 72. https://doi.org/10.3390/data7060072

AMA Style

Cirrincione AM, Reimonn CA, Harrison BJ, Rieger S. Longitudinal RNA Sequencing of Skin and DRG Neurons in Mice with Paclitaxel-Induced Peripheral Neuropathy. Data. 2022; 7(6):72. https://doi.org/10.3390/data7060072

Chicago/Turabian Style

Cirrincione, Anthony M., Cassandra A. Reimonn, Benjamin J. Harrison, and Sandra Rieger. 2022. "Longitudinal RNA Sequencing of Skin and DRG Neurons in Mice with Paclitaxel-Induced Peripheral Neuropathy" Data 7, no. 6: 72. https://doi.org/10.3390/data7060072

APA Style

Cirrincione, A. M., Reimonn, C. A., Harrison, B. J., & Rieger, S. (2022). Longitudinal RNA Sequencing of Skin and DRG Neurons in Mice with Paclitaxel-Induced Peripheral Neuropathy. Data, 7(6), 72. https://doi.org/10.3390/data7060072

Article Metrics

Back to TopTop