Decreased Brain Serotonin in rbfox1 Mutant Zebrafish and Partial Reversion of Behavioural Alterations by the SSRI Fluoxetine.

RBFOX1 functions as a master regulator of thousands of genes, exerting a pleiotropic effect on numerous neurodevelopmental and psychiatric disorders. A potential mechanism by which RBFOX1 may impact these disorders is through its modulation of serotonergic neurotransmission, a common target for pharmacological intervention in psychiatric conditions linked to RBFOX1. However, the precise effects of RBFOX1 on the serotonergic system remain largely unexplored. Here we show that homozygous rbfox1sa15940 zebrafish, which express a shorter, aberrant rbfox1 mRNA, have significantly reduced serotonin levels in telencephalon and diencephalon. We observed that the acute administration of fluoxetine partially reverses the associated behavioural alterations. The hyperactive phenotype and altered shoaling behaviour of the rbfox1sa15940/sa15940 zebrafish could be reversed with acute fluoxetine exposure in the Open Field and the Shoaling test, respectively. However, in the other paradigms, hyperactivity was not diminished, suggesting a distinct intrinsic motivation for locomotion in the different paradigms. Acute fluoxetine exposure did not reverse the alterations observed in the aggression and social novelty tests, suggesting the involvement of other neurological mechanisms in these behaviours. These findings underscore the importance of investigating the intricate working mechanisms of RBFOX1 in neurodevelopmental and psychiatric disorders to gain a better understanding of the associated disorders along with their pharmacological treatment.

The Adaptive Value of Chromosomal Inversions and Climatic Change-Studies on the Natural Populations of Drosophila subobscura from the Balkans.

The adaptive value of the Drosophila subobscura chromosomal inversion polymorphism with regard to environmental effects is well-known. However, the specific details of the inversion adaptations to the global warming scenario deserve to be analyzed. Toward this aim, polymorphism and karyotypes were studied in 574 individuals from Petnica (Serbia) in annual samples taken in June for the period 2019-2022. Comparing the results of Petnica (Cfa: humid subtropical climate) with those from Avala (Serbia: Cfb, temperate oceanic climate) and Font Groga (Barcelona, Spain; Csa: hot-summer Mediterranean climate), significant differences were observed for their chromosomal polymorphism. In Petnica, inversions from U and E chromosomes mainly reacted significantly with regard to temperature, humidity, and rainfall. Moreover, the inversion polymorphism from Petnica (2019-2022) was compared with that from 1995. In this period, a significant increase in mean and maximum temperature was observed. However, to properly explain the observed variations of inversions over time, it was necessary to carefully analyze annual seasonal changes and particular heat wave episodes. Interestingly, yearly fluctuations of U chromosome 'warm'-adapted inversions corresponded with opposite changes in 'non-thermal' inversions. Perhaps these types of inversions were not correctly defined with regard to thermal adaptation, or these fluctuations were also due to adaptations to other physical and/or biological variables. Finally, a joint study of chromosomal inversion polymorphism from many Balkan populations of D. subobscura indicated that different climatic regions presented distinct composition, including thermal-adapted inversions.

Long-term changes in the inversion chromosomal polymorphism: Drosophila subobscura population from Rasht (north of Iran).

Drosophila subobscura is characterized by a rich chromosomal polymorphism for inversions. Many inversions are adaptive to global warming and can be classified as 'warm' or 'cold' adapted. However, most studies were carried out from European populations located in the central area of the species distribution or from American colonizing populations. For this reason, we aimed to analyse the isolated and marginal Rasht population, located in the Hyrcanian forests area (Iran). The chromosomal polymorphism for inversions was compared with the previous Rasht samples (Rasht I and II) obtained 57 years ago. This polymorphism has changed based on the inversion composition and frequencies. Interestingly, the polymorphism for inversions was scarce and similar to that of Madeira, an isolated Atlantic island. Likely, this similarity is a consequence of the marginal location and isolation of the Rasht population. Also, the chromosomal thermal index (CTI) was 0.445, showing a significant increase over those from Rasht I (0.184) and II (0.210). All these observations were in agreement with the global warming expectations. Moreover, the CTI was also computed for Russian Caucasus and Turkish populations collected more than 40 years ago to better understand the adaptive potential of D. subobscura and to study the similarity between populations of different geographic areas. In summary, the inversions of D. subobscura also changed in marginal and isolated populations in agreement with the global warming expectations, and an open question is to know where is the threshold for this evolutionary change.

Interannual fluctuations in connectivity among crab populations (Liocarcinus depurator) along the Atlantic-Mediterranean transition.

An interesting evolutionary question that still remains open is the connectivity between marine populations. Marine currents can favour the dispersal of larvae or adults, but they can also produce eddies and gyres generating oceanographic fronts, thus limiting gene flow. To address this subject, we selected the Atlantic-Mediterranean transition, where several fronts are located: Gibraltar Strait (GS), Almeria-Oran Front (AOF) and Ibiza Channel (IC). Seven populations of the marine crab Liocarcinus depurator (Cadiz, West and East Alboran, Alacant, Valencia, Ebro Delta and North Catalonia) located along this transition were analysed in six consecutive years (2014-2019) using a fragment of the COI (Cytochrome Oxidase subunit I) gene. All sequences (966) belonged to two well defined haplogroups: ATL (most abundant in Atlantic waters) and MED (predominant in Mediterranean waters). Following a geographic variation, the frequency of ATL decreased significantly from Cadiz to North Catalonia. However, this variation presented steps due to the effect of oceanographic restrictions/fronts. Significant effects were recorded for GS (2015, 2017, 2018 and 2019), AOF (all years except 2018) and IC (2016). The intensity and precise location of these fronts changed over time. Multivariate analyses distinguished three main population groups: Cadiz, Alboran Sea and the remaining Mediterranean populations. These findings could be relevant to properly define Marine Protected Areas and for conservation and fisheries policies.

Deficiency of the ywhaz gene, involved in neurodevelopmental disorders, alters brain activity and behaviour in zebrafish.

Genetic variants in YWHAZ contribute to psychiatric disorders such as autism spectrum disorder and schizophrenia, and have been related to an impaired neurodevelopment in humans and mice. Here, we have used zebrafish to investigate the mechanisms by which YWHAZ contributes to neurodevelopmental disorders. We observed that ywhaz expression was pan-neuronal during developmental stages and restricted to Purkinje cells in the adult cerebellum, cells that are described to be reduced in number and size in autistic patients. We then performed whole-brain imaging in wild-type and ywhaz CRISPR/Cas9 knockout (KO) larvae and found altered neuronal activity and connectivity in the hindbrain. Adult ywhaz KO fish display decreased levels of monoamines in the hindbrain and freeze when exposed to novel stimuli, a phenotype that can be reversed with drugs that target monoamine neurotransmission. These findings suggest an important role for ywhaz in establishing neuronal connectivity during development and modulating both neurotransmission and behaviour in adults.

Adaptation of Drosophila subobscura chromosomal inversions to climatic variables: the Balkan natural population of Avala.

The adaptive value of chromosomal inversions continues raising relevant questions in evolutionary biology. In many species of the Drosophila genus, different inversions have been recognized to be related to thermal adaptation, but it is necessary to determine to which specific climatic variables the inversions are adaptive. With this aim, the behavior of thermal adapted inversions of Drosophila subobscura regarding climatic variables was studied in the natural population of Avala (Serbia) during the 2014-2017 period. The results obtained were compared with those previously reported in the Font Groga (Barcelona, Spain) population, which presents different climatic and environmental conditions. In both populations, it was observed that most thermal adapted inversions were significantly associated with the first, second or both principal components, which were related with maximum, minimum and mean temperatures. Moreover, a significant increase over years (2004-2017) for the minimum temperature was detected. In parallel, a significant variation over time in Avala was only observed for the frequencies of 'warm' and 'non-thermal' adapted inversions of the U chromosome. However, stability in the chromosomal inversion polymorphism was observed for the 2014-2017 period which might result from the temporal span of the study and/or selective process acting on the population.

DDC expression is not regulated by NFAT5 (TonEBP) in dopaminergic neural cell lines.

The nuclear factor of activated T-cells 5 (NFAT5), also known as tonicity-responsive enhancer-binding protein (TonEBP), is a transcription factor that regulates osmoadaptive response in multiple tissues and is highly expressed in the developing central nervous system. A former study reported that NFAT5 activation through hypertonic stress increases the expression of the dopa decarboxylase enzyme (DDC), also known as aromatic-l-amino-acid decarboxylase (AADC), in human renal proximal tubule cells, leading to an increase of dopamine synthesis. In a previous study, we identified NFAT5 as a candidate gene for cocaine dependence, a complex psychiatric disorder in which dopaminergic neurotransmission plays an important role. Therefore, to test the hypothesis that NFAT5 may also affect dopamine levels in the nervous system through the regulation of DDC expression, we examined this regulation using two neural dopaminergic cell lines, SH-SY5Y and PC12. The effect of NFAT5 on the expression of the neuronal isoform of DDC was evaluated by qRT-PCR. Upon hypertonic stress, NFAT5 was activated and accumulated into the nuclei and, subsequently, the expression of NFAT5 and its known targets sodium/myo-inositol cotransporter 1 (SMIT) and sodium chloride/taurine cotransporter (TAUT) increased, as expected. However, the expression of DDC decreased. When silencing the expression of NFAT5 with a specific shRNA we observed that the downregulation of DDC is independent from NFAT5 in both cell lines and is due to hypertonic stress. In conclusion, NFAT5 does not regulate the expression of the neuronal isoform of DDC in neural dopaminergic cell lines and, consequently, it does not modulate dopamine synthesis through DDC.

MiR-9, miR-153 and miR-124 are down-regulated by acute exposure to cocaine in a dopaminergic cell model and may contribute to cocaine dependence.

Cocaine is one of the most used psychostimulant drugs worldwide. MicroRNAs are post-transcriptional regulators of gene expression that are highly expressed in brain, and several studies have shown that cocaine can alter their expression. In a previous study, we identified several protein-coding genes that are differentially expressed in a dopaminergic neuron-like model after an acute exposure to cocaine. Now, we used the prediction tool WebGestalt to identify miRNA molecules potentially involved in the regulation of these genes. Using the same cellular model, we found that seven of these miRNAs are down-regulated by cocaine: miR-124-3p, miR-124-5p, miR-137, miR-101-3p, miR-9-5p, miR-369-3p and miR-153-3p, the last three not previously related to cocaine. Furthermore, we found that three of the miRNA genes that are differentially expressed in our model (hsa-miR-9-1, hsa-miR-153-1 and hsa-miR-124-3) are nominally associated with cocaine dependence in a case-control study (2,085 cases and 4,293 controls). In summary, we highlighted novel miRNAs that may be involved in those cocaine-induced changes of gene expression that underlie addiction. Moreover, we identified genetic variants that contribute to cocaine dependence in three of these miRNA genes, supporting the idea that genes differentially expressed under cocaine may play an important role in the susceptibility to cocaine dependence.

Association of the PLCB1 gene with drug dependence.

Genetic factors involved in the susceptibility to drug addiction still remain largely unknown. MiRNAs seem to play key roles in the drug-induced plasticity of the brain that likely drives the emergence of addiction. In this work we explored the role of miRNAs in drug addiction. With this aim, we selected 62 SNPs located in the 3'UTR of target genes that are predicted to alter the binding of miRNA molecules and performed a case-control association study in a Spanish sample of 735 cases (mainly cocaine-dependent subjects with multiple drug dependencies) and 739 controls. We found an association between rs1047383 in the PLCB1 gene and drug dependence that was replicated in an independent sample (663 cases and 667 controls). Then we selected 9 miRNAs predicted to bind the rs1047383 region, but none of them showed any effect on PLCB1 expression. We also assessed two miRNAs binding a region that contains a SNP in linkage disequilibrium with rs1047383, but although one of them, hsa-miR-582, was found to downregulate PLCB1, no differences were observed between alleles. Finally, we explored the possibility that PLCB1 expression is altered by cocaine and we observed a significant upregulation of the gene in the nucleus accumbens of cocaine abusers and in human dopaminergic-like neurons after cocaine treatment. Our results, together with previous studies, suggest that PLCB1 participates in the susceptibility to drug dependence.

A Highly Polymorphic Copy Number Variant in the NSF Gene is Associated with Cocaine Dependence.

Cocaine dependence is a complex psychiatric disorder involving both genetic and environmental factors. Several neurotransmitter systems mediate cocaine's effects, dependence and relapse, being the components of the neurotransmitter release machinery good candidates for the disorder. Previously, we identified a risk haplotype for cocaine dependence in the NSF gene, encoding the protein N-Ethylmaleimide-Sensitive Factor essential for synaptic vesicle turnover. Here we examined the possible contribution to cocaine dependence of a large copy number variant (CNV) that encompasses part of the NSF gene. We performed a case-control association study in a discovery sample (359 cases and 356 controls) and identified an association between cocaine dependence and the CNV (P = 0.013), that was confirmed in the replication sample (508 cases and 569 controls, P = 7.1e-03) and in a pooled analysis (P = 1.8e-04), with an over-representation of low number of copies in cases. Subsequently, we studied the functional impact of the CNV on gene expression and found that the levels of two NSF transcripts were significantly increased in peripheral blood mononuclear cells (PBMC) along with the number of copies of the CNV. These results, together with a previous study from our group, support the role of NSF in the susceptibility to cocaine dependence.

Replication study of previous migraine genome-wide association study findings in a Spanish sample of migraine with aura.

BACKGROUND: Migraine is a common disabling condition that affects approximately 15% of the population. Several genome-wide association studies have attempted to identify susceptibility variants involved in migraine, reporting several candidate loci for the disorder. METHODS: In order to replicate findings from previous genome-wide association studies, a case-control association study was performed. Twelve single nucleotide polymorphisms were genotyped in a Spanish sample of 512 migraine with aura patients and 535 migraine-free controls. RESULTS: Nominal associations were found for single nucleotide polymorphisms rs2651899 (within the PRDM16 gene), rs10166942 (near TRPM8), rs12134493 (close to TSPAN2) and rs10504861 (near MMP16) in our migraine with aura sample. CONCLUSIONS: Our study provides suggestive replication, in a Spanish migraine with aura sample, of four genome-wide association study findings previously reported in common migraine. However, larger sample sets should be explored to confirm our results.

Evolution of Chilean colonizing populations of Drosophila subobscura: lethal genes and chromosomal arrangements.

Knowledge of the frequency, distribution, and fate of lethal genes in chromosomal inversions helps to illuminate the evolution of recently founded populations. We analyze the relationship between lethal genes and inversions in two colonizing populations of D. subobscura in Chile. In the ancestral Palearctic populations of this species, lethal genes seem distributed at random on chromosomes. But in colonizing American populations, some lethal genes are associated with specific chromosomal arrangements. Some of these associated lethals were detected only during the first stages of the colonization (O( 3+4+2 )), and never thereafter, whereas others have persisted (O( 3+4+7 ) and O(5)). However, most lethal genes in American populations have been observed only once: they have arisen by novel mutation and soon disappear. Finally, recombination between different inversions has been observed in America. However, the persistence of lethal genes associated with the heterotic inversions O( 3+4+7 ) and O(5) could indicate that recombination inside these inversions is rare.

The genetic structure of Balkan populations of Drosophila subobscura.

Although Drosophila subobscura has been a model organism for European and American population geneticists, little information is available on the genetic structure of its natural populations. In this paper we report the estimates of some population parameters. We have used data from lethal allelism tests in four Balkan populations (Kamariste, Djerdap and Petnica in Serbia and Zanjic in Montenegro). In all populations, lethal genes were found to have a deleterious effect on heterozygotes. The N(e) values varied greatly from 370 (Petnica) to 19413 (Kamariste), depending on the habitat conditions and some environmental factors. Finally, gene flow between the Balkan populations was detected by the estimates of N(m) (from 4.68 in Petnica to 106.2 in Kamariste) and m (from 0.0041 in Djerdap to 0.0126 in Petnica). These results agree with those obtained in a previous study where the frequencies of allelism between populations were greater than predicted by independently arising lethal genes.