Whole-genome expression profile in zebrafish embryos after chronic exposure to morphine: identification of new genes associated with neuronal function and mu opioid receptor expression

Aniquinazoline B, a Fungal Natural Product, Activates the μ-Opioid Receptor

The development of new μ-opioid receptor (MOR) agonists without the undesirable side effects, such as addiction or respiratory depression, has been a difficult challenge over the years. In the search for new compounds, we screened our chemical database of over 40.000 substances and further assessed the best 100 through molecular docking. We selected the top 10 compounds and evaluated them for their biological activity and potential to influence cyclic adenosine monophosphate (cAMP) levels. From the tested compounds, compound 7, called aniquinazoline B, belonging to the quinazolinone alkaloids class and isolated from the marine fungus Aspergillus nidulans, showed promising results, by inhibiting cAMP levels and in vitro binding to MOR, verified through microscale thermophoresis. Transcriptomic data investigation profiled the genes affected by compound 7 and discovered activation of different pathways compared to opioids. The western blot analysis revealed compound 7 as a balanced ligand, activating both p-ERK1/2 and β-arrestin1/2 pathways, showing this is a favorable candidate to be further tested.

Delving into the significance of the His289Tyr single-nucleotide polymorphism in the glutamate ionotropic receptor kainate-1 (Grik1) gene of a genetically audiogenic seizure model

Genetic abnormalities affecting glutamate receptors are central to excitatory overload-driven neuronal mechanisms that culminate in seizures, making them pivotal targets in epilepsy research. Increasingly used to advance this field, the genetically audiogenic seizure hamster from Salamanca (GASH/Sal) exhibits generalized seizures triggered by high-intensity acoustic stimulation and harbors significant genetic variants recently identified through whole-exome sequencing. Here, we addressed the influence of the missense single-nucleotide polymorphism (C9586732T, p.His289Tyr) in the glutamate receptor ionotropic kainate-1 (Grik1) gene and its implications for the GASH/Sal seizure susceptibility. Using a protein 3D structure prediction, we showed a potential effect of this sequence variation, located in the amino-terminal domain, on the stability and/or conformation of the kainate receptor subunit-1 protein (GluK1). We further employed a multi-technique approach, encompassing gene expression analysis (RT-qPCR), Western blotting, and immunohistochemistry in bright-field and confocal fluorescence microscopy, to investigate critical seizure-associated brain regions in GASH/Sal animals under seizure-free conditions compared to matched wild-type controls. We detected disruptions in the transcriptional profile of the Grik1 gene within the audiogenic seizure-associated neuronal network. Alterations in GluK1 protein levels were also observed in various brain structures, accompanied by an unexpected lower molecular weight band in the inferior and superior colliculi. This correlated with substantial disparities in GluK1-immunolabeling distribution across multiple brain regions, including the cerebellum, hippocampus, subdivisions of the inferior and superior colliculi, and the prefrontal cortex. Notably, the diffuse immunolabeling accumulated within perikarya, axonal fibers and terminals, exhibiting a prominent concentration in proximity to the cell nucleus. This suggests potential disturbances in the GluK1-trafficking mechanism, which could subsequently affect glutamate synaptic transmission. Overall, our study sheds light on the genetic underpinnings of seizures and underscores the importance of investigating the molecular mechanisms behind synaptic dysfunction in epileptic neural networks, laying a crucial foundation for future research and therapeutic strategies targeting GluK1-containing kainate receptors.

Etazene induces developmental toxicity in vivo Danio rerio and in silico studies of new synthetic opioid derivative

Synthetic opioids are gaining more and more popularity among recreational users as well as regular abusers. One of such novel psychoactive substance, is etazene, which is the most popular opioid drug in the darknet market nowadays. Due to limited information available concerning its activity, we aimed to characterize its developmental toxicity, including cardiotoxicity with the use of in vivo Danio rerio and in silico tools. Moreover, we aimed, for the first time, to characterize the metabolite of etazene, which could become a potential marker of its use for future forensic analysis. The results of our study proved severe dose-dependent developmental toxicity of etazene (applied concentrations 10-300 µM), including an increase in mortality, developmental malformations, and serious cardiotoxic effects, as compared with well-known and used opioid-morphine (applied concentrations 1-50 mM). In silico findings indicate the high toxic potential of etazene which may lead to drug-drug interactions and accumulation of substances. Furthermore, phase I metabolite of etazene resulting from N-dealkylation reaction was identified, and therefore it should be considered as a target for toxicological screening. Nonetheless, the exact mechanism of observed effects in response to etazene should be further examined.

Revelation of candidate genes and molecular mechanism of reproductive seasonality in female rohu (Labeo rohita Ham.) by RNA sequencing

BACKGROUND

Carp fish, rohu (Labeo rohita Ham.) is important freshwater aquaculture species of South-East Asia having seasonal reproductive rhythm. There is no holistic study at transcriptome level revealing key candidate genes involved in such circannual rhythm regulated by biological clock genes (BCGs). Seasonality manifestation has two contrasting phases of reproduction, i.e., post-spawning resting and initiation of gonadal activity appropriate for revealing the associated candidate genes. It can be deciphered by RNA sequencing of tissues involved in BPGL (Brain-Pituitary-Gonad-Liver) axis controlling seasonality. How far such BCGs of this fish are evolutionarily conserved across different phyla is unknown. Such study can be of further use to enhance fish productivity as seasonality restricts seed production beyond monsoon season.

RESULT

A total of ~ 150 Gb of transcriptomic data of four tissues viz., BPGL were generated using Illumina TruSeq. De-novo assembled BPGL tissues revealed 75,554 differentially expressed transcripts, 115,534 SSRs, 65,584 SNPs, 514 pathways, 5379 transcription factors, 187 mature miRNA which regulates candidate genes represented by 1576 differentially expressed transcripts are available in the form of web-genomic resources. Findings were validated by qPCR. This is the first report in carp fish having 32 BCGs, found widely conserved in fish, amphibian, reptile, birds, prototheria, marsupials and placental mammals. This is due to universal mechanism of rhythmicity in response to environment and earth rotation having adaptive and reproductive significance.

CONCLUSION

This study elucidates evolutionary conserved mechanism of photo-periodism sensing, neuroendocrine secretion, metabolism and yolk synthesis in liver, gonadal maturation, muscular growth with sensory and auditory perception in this fish. Study reveals fish as a good model for research on biological clock besides its relevance in reproductive efficiency enhancement.

Cocaine self-administration induces sex-dependent protein expression in the nucleus accumbens

Substance use disorder (SUD) is a chronic neuropsychiatric condition characterized by long-lasting alterations in the neural circuitry regulating reward and motivation. Substantial work has focused on characterizing the molecular substrates that underlie these persistent changes in neural function and behavior. However, this work has overwhelmingly focused on male subjects, despite mounting clinical and preclinical evidence that females demonstrate dissimilar progression to SUD and responsivity to stimulant drugs of abuse, such as cocaine. Here, we show that sex is a critical biological variable that defines drug-induced plasticity in the nucleus accumbens (NAc). Using quantitative mass spectrometry, we assessed the protein expression patterns induced by cocaine self-administration and demonstrated unique molecular profiles between males and females. We show that 1. Cocaine self-administration induces non-overlapping protein expression patterns in significantly regulated proteins in males and females and 2. Critically, cocaine-induced protein regulation differentially interacts with sex to eliminate basal sexual dimorphisms in the proteome. Finally, eliminating these baseline differences in the proteome is concomitant with the elimination of sex differences in behavior for non-drug rewards. Together, these data suggest that cocaine administration is capable of rewriting basal proteomic function and reward-associated behaviors.

COPB2: a transport protein with multifaceted roles in cancer development and progression

The Coatomer protein complex subunit beta 2 (COPB2) is involved in the formation of the COPI coatomer protein complex and is responsible for the transport of vesicles between the Golgi apparatus and the endoplasmic reticulum. It plays an important role in maintaining the integrity of these cellular organelles, as well as in maintaining cell homeostasis. More importantly, COPB2 plays key roles in embryonic development and tumor progression. COPB2 is regarded as a vital oncogene in several cancer types and has been implicated in tumor cell proliferation, survival, invasion, and metastasis. Here, we summarize the current knowledge on the roles of COPB2 in cancer development and progression in the context of the hallmarks of cancer.

Zebrafish (Danio rerio) larvae show behavioral and embryonic development defects when exposed to opioids at embryo stage

Opioid abuse continues to plague society, and in recent years, there has been an epidemic, leading to increased addiction and death. It is poorly understood how prenatal opioid use affects the lives of children. The aim of this work was to evaluate the effect of early embryonic codeine or morphine exposure in zebrafish (Danio rerio), examining gastrulation progression (epiboly), teratogenic effects, mortality and locomotor behavior response to light/dark cycles. Zebrafish embryos were exposed to codeine or morphine (designated C or M) at 1, 5 or 10 mg/L (designated 01, 05 or 10, respectively) from 3 to 24 h postfertilization (hpf) or from 3 to 48 hpf (designated -24 or - 48 for 1 or 2 days of exposure, respectively). The C10-24, C01-48, C05-48 and C10-48 groups showed significantly smaller eyes than control larvae at 7 days postfertilization (dpf). Locomotor behavior of control larvae in light/dark cycles showed greater swimming time and distance in dark cycles. Two-day codeine exposure produced strong effects, showing no significant response due to light/dark cycles in distance moved. Morphine exposed groups showed similar effects as observed in 2-day codeine exposed groups, showing less large movement activity and also no significant difference between inactive duration in response to light/dark cycles. In conclusion, we observed low teratogenic effects and mortality effects. Animals exposed to high levels and higher exposure times of opioids were hypoactive, relative to controls, in the dark period. Future studies will be needed to understand the neural defects producing behavior changes.

Dry Needling of a Healthy Rat Achilles Tendon Increases Its Gene Expressions: A Pilot Study

BACKGROUND

Tendon dry needling is a potential treatment for tendinopathies. Several hypotheses have been proposed to explain its underlying mechanisms. No studies (to the best of our knowledge) have investigated changes in gene expression.

OBJECTIVE

To investigate histological and gene expression changes after the application of dry needling to the healthy Achilles tendons of rats.

METHODS

Six Sprague-Dawley male rats were randomly divided into two groups: no intervention or dry needling. Dry needling consisted of three sessions (once per week) to the Achilles tendon. Molecular expression of several genes involved in tendon repair and remodeling (e.g., Cox2, Mmp2, Mmp9, Col1a1, Col3a1, Vefg, and Scx) was assessed 7 days after the last needling session (day 28) or 28 days after for the no-intervention group. Histological tissue changes were determined with hematoxylin-eosin analyses.

RESULTS

The hematoxylin-eosin-stained images revealed no substantial differences in collagen structure or the presence of inflammatory cells between the dry needling and no-intervention groups. A significant increase in the molecular expression of Cox2, Mmp2, Col3a1, and Scx genes was observed in Achilles tendons treated with dry needling when compared with the no-intervention group.

CONCLUSION

This animal pilot study found that the application of dry needling to the healthy Achilles tendons of rats is able to increase the expression of genes associated with collagen regeneration and tissue remodeling of the extracellular matrix with no further histological damage to the tendon.

Changes in Gene Expression Associated with Collagen Regeneration and Remodeling of Extracellular Matrix after Percutaneous Electrolysis on Collagenase-Induced Achilles Tendinopathy in an Experimental Animal Model: A Pilot Study

Percutaneous electrolysis is an emerging intervention proposed for the management of tendinopathies. Tendon pathology is characterized by a significant cell response to injury and gene expression. No study investigating changes in expression of those genes associated with collagen regeneration and remodeling of extracellular matrix has been conducted. The aim of this pilot study was to investigate gene expression changes after the application of percutaneous electrolysis on experimentally induced Achilles tendinopathy with collagenase injection in an animal model. Fifteen Sprague Dawley male rats were randomly divided into three different groups (no treatment vs. percutaneous electrolysis vs. needling). Achilles tendinopathy was experimentally induced with a single bolus of collagenase injection. Interventions consisted of 3 sessions (one per week) of percutaneous electrolysis or just needling. The rats were euthanized, and molecular expression of genes involved in tendon repair and remodeling, e.g., Cox2, Mmp2, Mmp9, Col1a1, Col3a1, Vegf and Scx, was examined at 28 days after injury. Histological tissue changes were determined with hematoxylin–eosin and safranin O analyses. The images of hematoxylin–eosin and Safranin O tissue images revealed that collagenase injection induced histological changes compatible with a tendinopathy. No further histological changes were observed after the application of percutaneous electrolysis or needling. A significant increase in molecular expression of Cox2, Mmp9 and Vegf genes was observed in Achilles tendons treated with percutaneous electrolysis to a greater extent than after just needling. The expression of Mmp2, Col1a1, Col3a1, or Scx genes also increased, but did not reach statistical significance. This animal study demonstrated that percutaneous electrolysis applied on an experimentally induced Achilles tendinopathy model could increase the expression of some genes associated with collagen regeneration and remodeling of extracellular matrix. The observed gene overexpression was higher with percutaneous electrolysis than with just needling.

Marijuana and Opioid Use during Pregnancy: Using Zebrafish to Gain Understanding of Congenital Anomalies Caused by Drug Exposure during Development

Marijuana and opioid addictions have increased alarmingly in recent decades, especially in the United States, posing threats to society. When the drug user is a pregnant mother, there is a serious risk to the developing baby. Congenital anomalies are associated with prenatal exposure to marijuana and opioids. Here, we summarize the current data on the prevalence of marijuana and opioid use among the people of the United States, particularly pregnant mothers. We also summarize the current zebrafish studies used to model and understand the effects of these drug exposures during development and to understand the behavioral changes after exposure. Zebrafish experiments recapitulate the drug effects seen in human addicts and the birth defects seen in human babies prenatally exposed to marijuana and opioids. Zebrafish show great potential as an easy and inexpensive model for screening compounds for their ability to mitigate the drug effects, which could lead to new therapeutics.

Understanding the neurobiological effects of drug abuse: Lessons from zebrafish models

Drug abuse and brain disorders related to drug comsumption are public health problems with harmful individual and social consequences. The identification of therapeutic targets and precise pharmacological treatments to these neuropsychiatric conditions associated with drug abuse are urgently needed. Understanding the link between neurobiological mechanisms and behavior is a key aspect of elucidating drug abuse-related targets. Due to various molecular, biochemical, pharmacological, and physiological features, the zebrafish (Danio rerio) has been considered a suitable vertebrate for modeling complex processes involved in drug abuse responses. In this review, we discuss how the zebrafish has been successfully used for modeling neurobehavioral phenotypes related to drug abuse and review the effects of opioids, cannabinoids, alcohol, nicotine, and psychedelic drugs on the central nervous system (CNS). Moreover, we summarize recent advances in zebrafish-based studies and outline potential advantages and limitations of the existing zebrafish models to explore the neurochemical bases of drug abuse and addiction. Finally, we discuss how the use of zebrafish models may present fruitful approaches to provide valuable clinically translatable data.

The Yin and Yang of Autosomal Recessive Primary Microcephaly Genes: Insights from Neurogenesis and Carcinogenesis

The stem cells of neurogenesis and carcinogenesis share many properties, including proliferative rate, an extensive replicative potential, the potential to generate different cell types of a given tissue, and an ability to independently migrate to a damaged area. This is also evidenced by the common molecular principles regulating key processes associated with cell division and apoptosis. Autosomal recessive primary microcephaly (MCPH) is a neurogenic mitotic disorder that is characterized by decreased brain size and mental retardation. Until now, a total of 25 genes have been identified that are known to be associated with MCPH. The inactivation (yin) of most MCPH genes leads to neurogenesis defects, while the upregulation (yang) of some MCPH genes is associated with different kinds of carcinogenesis. Here, we try to summarize the roles of MCPH genes in these two diseases and explore the underlying mechanisms, which will help us to explore new, attractive approaches to targeting tumor cells that are resistant to the current therapies.

Review on the use of zebrafish embryos to study the effects of anesthetics during early development

Over the years, the potential toxicity of anesthetics has raised serious concerns about its safe use during pregnancy. As evidence emerged from research in animal models, showing that some anesthetic drugs are potential teratogenic, the determination of the risk of exposures to anesthetic drugs at early life stages became mandatory. However, due to inaccessibility and ethical constrains related to experimental conditions, the use of early life stages in mammalian models is limited. In this regard, some animal and nonanimal models have been suggested to surpass mammalian use in experimentation. Among them, the zebrafish embryo test has been recognized as a promising alternative in toxicology research, as well as an inexpensive and practical test. Substantial information collected from developmental research following compounds exposure, has contributed to the application of zebrafish assays in research, although only a few studies have focused on the use of early life stages of zebrafish to evaluate the developmental effects of anesthetics. Based on the recent advances of science and technology, there is a clear potential for zebrafish early life stages to provide new insights into anesthetics teratogenicity. This review provides an overview of recent anesthesia research using zebrafish embryos, demonstrating its usefulness to the anesthesia field, discussing the recent findings on various aspects related to the effects of anesthetics during early life development and the strengths and limitations of this model system.

Opioid Neurobiology, Neurogenetics and Neuropharmacology in Zebrafish

Despite the high prevalence of medicinal use and abuse of opioids, their neurobiology and mechanisms of action are not fully understood. Experimental (animal) models are critical for improving our understanding of opioid effects in vivo. As zebrafish (Danio rerio) are increasingly utilized as a powerful model organism in neuroscience research, mounting evidence suggests these fish as a useful tool to study opioid neurobiology. Here, we discuss the zebrafish opioid system with specific focus on opioid gene expression, existing genetic models, as well as its pharmacological and developmental regulation. As many human brain diseases involve pain and aberrant reward, we also summarize zebrafish models relevant to opioid regulation of pain and addiction, including evidence of functional interplay between the opioid system and central dopaminergic and other neurotransmitter mechanisms. Additionally, we critically evaluate the limitations of zebrafish models for translational opioid research and emphasize their developing utility for improving our understanding of evolutionarily conserved mechanisms of pain-related, addictive, affective and other behaviors, as well as for fostering opioid-related drug discovery.

The Transcriptome of the Zebrafish Embryo After Chemical Exposure: A Meta-Analysis

Numerous studies have been published in the past years investigating the transcriptome of the zebrafish embryo (ZFE) upon being subjected to chemical stress. Aiming at a more mechanistic understanding of the results of such studies, knowledge about commonalities of transcript regulation in response to chemical stress is needed. Thus, our goal in this study was to identify and interpret genes and gene sets constituting a general response to chemical exposure. Therefore, we aggregated and reanalyzed published toxicogenomics data obtained with the ZFE. We found that overlap of differentially transcribed genes in response to chemical stress across independent studies is generally low and the most commonly differentially transcribed genes appear in less than 50% of all treatments across studies. However, effect size analysis revealed several genes showing a common trend of differential expression, among which genes related to calcium homeostasis emerged as key, especially in exposure settings up to 24 h post-fertilization. Additionally, we found that these and other downregulated genes are often linked to anatomical regions developing during the respective exposure period. Genes showing a trend of increased expression were, among others, linked to signaling pathways (e.g., Wnt, Fgf) as well as lysosomal structures and apoptosis. The findings of this study help to increase the understanding of chemical stress responses in the developing zebrafish embryo and provide a starting point to improve experimental designs for this model system. In future, improved time- and concentration-resolved experiments should offer better understanding of stress response patterns and access to mechanistic information.

Identifying novel members of the Wntless interactome through genetic and candidate gene approaches

Wnt signaling is an important pathway that regulates several aspects of embryogenesis, stem cell maintenance, and neural connectivity. We have recently determined that opioids decrease Wnt secretion, presumably by inhibiting the recycling of the Wnt trafficking protein Wntless (Wls). This effect appears to be mediated by protein-protein interaction between Wls and the mu-opioid receptor (MOR), the primary cellular target of opioid drugs. The goal of this study was to identify novel protein interactors of Wls that are expressed in the brain and may also play a role in reward or addiction. Using genetic and candidate gene approaches, we show that among a variety of protein, Wls interacts with the dopamine transporter (target of cocaine), cannabinoid receptors (target of THC), Adenosine A2A receptor (target of caffeine), and SGIP1 (endocytic regulator of cannabinoid receptors). Our study shows that aside from opioid receptors, Wntless interacts with additional proteins involved in reward and/or addiction. Future studies will determine whether Wntless and WNT signaling play a more universal role in these processes.

Reduction in activity by noxious chemical stimulation is ameliorated by immersion in analgesic drugs in zebrafish

Research has recently demonstrated that larval zebrafish show similar molecular responses to nociception to those of adults. Our study explored whether unprotected larval zebrafish exhibited altered behaviour after exposure to noxious chemicals and screened a range of analgesic drugs to determine their efficacy to reduce these responses. This approach aimed to validate larval zebrafish as a reliable replacement for adults as well as providing a high-throughput means of analysing behavioural responses. Zebrafish at 5 days post-fertilization were exposed to known noxious stimuli: acetic acid (0.01%, 0.1% and 0.25%) and citric acid (0.1%, 1% and 5%). The behavioural response of each was recorded and analysed using novel tracking software that measures time spent active in 25 larvae at one time. Subsequently, the efficacy of aspirin, lidocaine, morphine and flunixin as analgesics after exposure to 0.1% acetic acid was tested. Larvae exposed to 0.1% and 0.25% acetic acid spent less time active, whereas those exposed to 0.01% acetic acid and 0.1–5% citric acid showed an increase in swimming activity. Administration of 2.5 mg l−1 aspirin, 5 mg l−1 lidocaine and 48 mg l−1 morphine prevented the behavioural changes induced by acetic acid. These results suggest that larvae respond to a noxious challenge in a similar way to adult zebrafish and other vertebrates and that the effect of nociception on activity can be ameliorated by using analgesics. Therefore, adopting larval zebrafish could represent a direct replacement of a protected adult fish with a non-protected form in pain- and nociception-related research.

Endogenous Opiates and Behavior: 2015

This paper is the thirty-eighth consecutive installment of the annual review of research concerning the endogenous opioid system. It summarizes papers published during 2015 that studied the behavioral effects of molecular, pharmacological and genetic manipulation of opioid peptides, opioid receptors, opioid agonists and opioid antagonists. The particular topics that continue to be covered include the molecular-biochemical effects and neurochemical localization studies of endogenous opioids and their receptors related to behavior, and the roles of these opioid peptides and receptors in pain and analgesia, stress and social status, tolerance and dependence, learning and memory, eating and drinking, drug abuse and alcohol, sexual activity and hormones, pregnancy, development and endocrinology, mental illness and mood, seizures and neurologic disorders, electrical-related activity and neurophysiology, general activity and locomotion, gastrointestinal, renal and hepatic functions, cardiovascular responses, respiration and thermoregulation, and immunological responses.

Early avoidance of a heroin-paired taste-cue and subsequent addiction-like behavior in rats

The ability to predict individual vulnerability to substance abuse would allow for a better understanding of the progression of the disease and development of better methods for prevention and/or early intervention. Here we use drug-induced devaluation of a saccharin cue in an effort to predict later addiction-like behavior in a model akin to that used by Deroche-Gamonet et al. (2004) and seek to link such vulnerability to changes in expression of various mu opioid receptor and D2 receptor-interacting proteins in brain. The results show that the greatest heroin-induced suppression of intake of a saccharin cue is associated with the greatest vulnerability to later addiction-like behavior and to differences in the expression of WLS, β-catenin, and NCS-1 in brain compared to rats that exhibited the least suppression of intake of the heroin-paired cue and/or saline controls. Finally, because the self-administration model employed produced no significant differences in drug intake between groups, overall, the resultant changes in protein expression can be more closely linked to individual differences in motivation for drug.

Effects of 17 α-methyltestosterone on transcriptome, gonadal histology and sex steroid hormones in rare minnow Gobiocypris rarus

The 17α-methyltestosterone (MT), a synthetic androgen, is known for its interference effects on the endocrine system. Aiming to investigate the transcriptome profiling of gonads induced by MT and to understand the molecular mechanism by which MT causes adverse effects in fish, transcriptome profiling of gonads, gonadal histology and the sex steroid hormones in response to MT were analyzed in Gobiocypris rarus. Eight libraries, 4 from the ovary and 4 from the testis, were constructed and sequenced and then a total number of clean reads per sample ranging from 7.03 to 9.99 million were obtained. In females, a total of 191 transcripts were differentially regulated by MT, consisting of 102 up-regulated transcripts and 89 down-regulated transcripts. In males, 268 differentially expressed genes with 108 up-regulated and 160 down-regulated were detected upon MT exposure. Testosterone serves as the major sex steroid hormone content in G. rarus of both sexes. The concentrations of 17β-estradiol, testosterone and 11-ketotestosterone were significantly increased in females and decreased in males after MT exposure. Interestingly, MT caused a decreased number of vitellogenic oocytes in the ovary and spermatozoa in the testis. After MT exposure, four differentially expressed genes (ndufa4, slc1a3a, caskin-2 and rpt3) were found in G. rarus of both sexes. Overall, we suggest that MT seemed to affect genes involved in pathways related to physiological processes in the gonads of G. rarus. These processes include the electron transfer of Complex IV, endothelial cell activation, axon growth and guidance, and proteasome assembly and glutamate transport metabolic.