Duplication and Losses of Opsin Genes in Lophotrochozoan Evolution.

Opsins are G-coupled receptors playing a key role in metazoan visual processes. While many studies enriched our understanding of opsin diversity in several animal clades, the opsin evolution in Lophotrochozoa, one of the major metazoan groups, remains poorly understood. Using recently developed phylogenetic approaches, we investigated the opsin evolution in 74 lophotrochozoan genomes. We found that the common ancestor of Lophotrochozoa possessed at least seven opsin paralog groups that underwent divergent evolutionary history in the different phyla. Furthermore, we showed for the first time opsin-related molecules in Bilateria that we named pseudopsins, which may prove critical in uncovering opsin evolution.

Amphioxus functional genomics and the origins of vertebrate gene regulation.

Vertebrates have greatly elaborated the basic chordate body plan and evolved highly distinctive genomes that have been sculpted by two whole-genome duplications. Here we sequence the genome of the Mediterranean amphioxus (Branchiostoma lanceolatum) and characterize DNA methylation, chromatin accessibility, histone modifications and transcriptomes across multiple developmental stages and adult tissues to investigate the evolution of the regulation of the chordate genome. Comparisons with vertebrates identify an intermediate stage in the evolution of differentially methylated enhancers, and a high conservation of gene expression and its cis-regulatory logic between amphioxus and vertebrates that occurs maximally at an earlier mid-embryonic phylotypic period. We analyse regulatory evolution after whole-genome duplications, and find that-in vertebrates-over 80% of broadly expressed gene families with multiple paralogues derived from whole-genome duplications have members that restricted their ancestral expression, and underwent specialization rather than subfunctionalization. Counter-intuitively, paralogues that restricted their expression increased the complexity of their regulatory landscapes. These data pave the way for a better understanding of the regulatory principles that underlie key vertebrate innovations.

Digital Marine: An online platform for blended learning in a marine experimental biology module, the Schmid Training Course.

For over 20 years, the Schmid Training Course (STC) has offered unique opportunities for marine biology students from European universities to learn about marine model organisms. While the topics of the course have continuously changed over the years with the advent of new research techniques and discoveries, the pedagogical approach has remained largely the same - a combination of lectures, lab practicals, and field excursions. Several life science researchers, who have taught in the STC for many years, sought to bring the course's pedagogical approach into the 21st century, and with the support of Erasmus+ Programme of the European Community funding, the Digital Marine project was developed. Digital Marine began in 2018 as an international partnership between the six research centers from which the STC instructors hail, and its main objective was to introduce a flipped, blended approach to learning and teaching with respect to established and emerging marine biological model systems. The Digital Marine platform, which covers 12 marine model organisms, is now publicly available.

Hidden species diversity and mito-nuclear discordance within the Mediterranean cone snail, Lautoconus ventricosus.

The Mediterranean cone snail, Lautoconus ventricosus, is currently considered a single species inhabiting the whole Mediterranean basin and the adjacent Atlantic coasts. Yet, no population genetic study has assessed its taxonomic status. Here, we collected 245 individuals from 75 localities throughout the Mediterranean Sea and used cox1 barcodes, complete mitochondrial genomes, and genome skims to test whether L. ventricosus represents a complex of cryptic species. The maximum likelihood phylogeny based on complete mitochondrial genomes recovered six main clades (hereby named blue, brown, green, orange, red, and violet) with sufficient sequence divergence to be considered putative species. On the other hand, phylogenomic analyses based on 437 nuclear genes only recovered four out of the six clades: blue and orange clades were thoroughly mixed and the brown one was not recovered. This mito-nuclear discordance revealed instances of incomplete lineage sorting and introgression, and may have caused important differences in the dating of main cladogenetic events. Species delimitation tests proposed the existence of at least three species: green, violet, and red + blue + orange (i.e., cyan). Green plus cyan (with sympatric distributions) and violet, had West and East Mediterranean distributions, respectively, mostly separated by the Siculo-Tunisian biogeographical barrier. Morphometric analyses of the shell using species hypotheses as factor and shell length as covariate showed that the discrimination power of the studied parameters was only 70.2%, reinforcing the cryptic nature of the uncovered species, and the importance of integrative taxonomic approaches considering morphology, ecology, biogeography, and mitochondrial and nuclear population genetic variation.

Nitric Oxide Function and Nitric Oxide Synthase Evolution in Aquatic Chordates.

Nitric oxide (NO) is a key signaling molecule in almost all organisms and is active in a variety of physiological and pathological processes. Our understanding of the peculiarities and functions of this simple gas has increased considerably by extending studies to non-mammal vertebrates and invertebrates. In this review, we report the nitric oxide synthase (Nos) genes so far characterized in chordates and provide an extensive, detailed, and comparative analysis of the function of NO in the aquatic chordates tunicates, cephalochordates, teleost fishes, and amphibians. This comprehensive set of data adds new elements to our understanding of Nos evolution, from the single gene commonly found in invertebrates to the three genes present in vertebrates.

Individual differences and knockout in zebrafish reveal similar cognitive effects of BDNF between teleosts and mammals.

The remarkable similarities in cognitive performance between teleosts and mammals suggest that the underlying cognitive mechanisms might also be similar in these two groups. We tested this hypothesis by assessing the effects of the brain-derived neurotrophic factor (BDNF), which is critical for mammalian cognitive functioning, on fish's cognitive abilities. We found that individual differences in zebrafish's learning abilities were positively correlated with bdnf expression. Moreover, a CRISPR/Cas9 mutant zebrafish line that lacks the BDNF gene (bdnf-/-) showed remarkable learning deficits. Half of the mutants failed a colour discrimination task, whereas the remaining mutants learned the task slowly, taking three times longer than control bdnf+/+ zebrafish. The mutants also took twice as long to acquire a T-maze task compared to control zebrafish and showed difficulties exerting inhibitory control. An analysis of habituation learning revealed that cognitive impairment in mutants emerges early during development, but could be rescued with a synthetic BDNF agonist. Overall, our study indicates that BDNF has a similar activational effect on cognitive performance in zebrafish and in mammals, supporting the idea that its function is conserved in vertebrates.

Evolution of the nitric oxide synthase family in vertebrates and novel insights in gill development.

Nitric oxide (NO) is an ancestral key signalling molecule essential for life and has enormous versatility in biological systems, including cardiovascular homeostasis, neurotransmission and immunity. Although our knowledge of NO synthases (Nos), the enzymes that synthesize NO in vivo, is substantial, the origin of a large and diversified repertoire of nos gene orthologues in fishes with respect to tetrapods remains a puzzle. The recent identification of nos3 in the ray-finned fish spotted gar, which was considered lost in this lineage, changed this perspective. This finding prompted us to explore nos gene evolution, surveying vertebrate species representing key evolutionary nodes. This study provides noteworthy findings: first, nos2 experienced several lineage-specific gene duplications and losses. Second, nos3 was found to be lost independently in two different teleost lineages, Elopomorpha and Clupeocephala. Third, the expression of at least one nos paralogue in the gills of developing shark, bichir, sturgeon, and gar, but not in lamprey, suggests that nos expression in this organ may have arisen in the last common ancestor of gnathostomes. These results provide a framework for continuing research on nos genes' roles, highlighting subfunctionalization and reciprocal loss of function that occurred in different lineages during vertebrate genome duplications.

Brain Proteome and Behavioural Analysis in Wild Type, BDNF+/- and BDNF-/- Adult Zebrafish (Danio rerio) Exposed to Two Different Temperatures.

Experimental evidence suggests that environmental stress conditions can alter the expression of BDNF and that the expression of this neurotrophin influences behavioural responses in mammalian models. It has been recently demonstrated that exposure to 34 °C for 21 days alters the brain proteome and behaviour in zebrafish. The aim of this work was to investigate the role of BDNF in the nervous system of adult zebrafish under control and heat treatment conditions. For this purpose, zebrafish from three different genotypes (wild type, heterozygous BDNF+/- and knock out BDNF-/-) were kept for 21 days at 26 °C or 34 °C and then euthanized for brain molecular analyses or subjected to behavioural tests (Y-maze test, novel tank test, light and dark test, social preference test, mirror biting test) for assessing behavioural aspects such as boldness, anxiety, social preference, aggressive behaviour, interest for the novel environment and exploration. qRT-PCR analysis showed the reduction of gene expression of BDNF and its receptors after heat treatment in wild type zebrafish. Moreover, proteomic analysis and behavioural tests showed genotype- and temperature-dependent effects on brain proteome and behavioural responding. Overall, the absent expression of BDNF in KO alters (1) the brain proteome by reducing the expression of proteins involved in synapse functioning and neurotransmitter-mediated transduction; (2) the behaviour, which can be interpreted as bolder and less anxious and (3) the cellular and behavioural response to thermal treatment.

Short-term exposure to the simple polyphenolic compound gallic acid induces neuronal hyperactivity in zebrafish larvae.

A growing body of evidence suggests that the biological effects of polyphenols are not restricted to antioxidant activity, but they exert a wide range of modulatory effects on metabolic pathways, cellular signaling and gene expression. In this study, we tested the minimum safe concentration of gallic acid (GA) in 72 hpf zebrafish larvae in order to evaluate the effects on the central nervous system and the behavioral response. We showed that a short exposure (30 min) induces the depletion of the two main excitatory and inhibitory neurotransmitters, Glu and GABA, respectively, in the larval nervous system. The acute impairment of GABAergic-glutamatergic balance was paralleled by an increase of the fosab neuronal activity marker in specific brain areas, such as the forebrain, olfactory bulbs, pallial area, ventral midbrain, tegmentum, and the medulla oblongata ventral area. The neuronal excitation was mirrored by the increased cumulative motor response. The inhibition of the olfactory epithelium with brief cadmium exposition suggests a direct involvement of olfaction in the larvae response to GA. Our results demonstrate that a brief exposure to GA induces motoneuronal hyperexcitability in zebrafish. The behavioral response was probably elicited through the activation of an odorous, or chemical, stimulus. The specificity of the activated neuronal territories suggests the involvement of additional signaling pathways. Although the underlying molecular mechanisms remain to be elucidated, our data support the hypothesis that GA acts as an excitatory molecule, capable of inducing a specific nerve response. These results offer a new vision on potential effects of GA.

Functional conserved non-coding elements among tunicates and chordates.

Non-coding regions with dozens to several hundred base pairs of extreme conservation have been found in all metazoan genomes. The distribution of these conserved non-coding elements (CNE) within and across genomes has suggested that many of them may have roles as transcriptional regulatory elements. A combination of bioinformatics and experimental approaches can be used to identify CNEs with regulatory activity in phylogenetically distant species. Nevertheless, the high divergent rate of genomic sequences of several organisms, such as tunicates, complicates the characterization of these conserved elements and very few examples really may prove their functional activity. We used a comparative approach to facilitate the identification of CNEs among distantly related or highly divergent species and experimentally demonstrated the functional significance of these novel CNEs. We first experimentally tested, in C. robusta and D. rerio transgenic embryos, the regulatory activity of conserved elements associated to genes involved in developmental control among different chordates (Homo sapiens and Danio rerio for vertebrates, Ciona robusta and Ciona savignyi for tunicates and Branchiostoma floridae for cephalochordates). Once demonstrated the cross-species functional conservation of these CNEs, the same gene loci were used as references to locate homologous regions and possible CNEs in available tunicate genomes. Comparison of tunicate-specific and chordate-specific CNEs revealed absence of conservation of the regulatory elements in spite of conservation of regulatory patterns, likely due to evolutionary specification of the respective developmental networks. This result highlights the importance of an integrative in-silico/in-vivo approach to CNEs investigation, encompassing both bioinformatics, essential for putative CNEs identification, and laboratory experiments, pivotal for the understanding of CNEs functionality.

Transphyletic conservation of nitric oxide synthase regulation in cephalochordates and tunicates.

Nitric oxide synthase is ubiquitously present in metazoans and is involved in a wide range of biological processes. Three distinct Nos genes have been so far identified in vertebrates exhibiting a complex expression pattern and transcriptional regulation. Nevertheless, although independent events of Nos duplication have been observed in several taxa, only few studies described the regulatory mechanisms responsible for their activation in non-vertebrate animals. To shed light on the mechanisms underlying neuronal-type Nos expression, we focused on two non-vertebrate chordates: the cephalochordate Branchiostoma lanceolatum and the tunicate Ciona robusta. Here, throughout transphyletic and transgenic approaches, we identified genomic regions in both species acting as Nos functional enhancers during development. In vivo analyses of Nos genomic fragments revealed their ability to recapitulate the endogenous expression territories. Therefore, our results suggest the existence of evolutionary conserved mechanisms responsible for neuronal-type Nos regulation in non-vertebrate chordates. In conclusion, this study paves the way for future characterization of conserved transcriptional logic underlying the expression of neuronal-type Nos genes in chordates.

The EJC component Magoh in non-vertebrate chordates.

Earliest craniates possess a newly enlarged, elaborated forebrain with new cell types and neuronal networks. A key question in vertebrate evolution is when and how this cerebral expansion took place. The exon-junction complex (EJC) plays an essential role in mRNA processing of all Eukarya. Recently, it has been proposed that the EJC represses recursive RNA splicing in Deuterostomes, with implication in human brain diseases like microcephaly and depression. However, the EJC or EJC subunit contribution to brain development in non-vertebrate Deuterostomes remained unknown. Being interested in the evolution of chordate characters, we focused on the model species, Branchiostoma lanceolatum (Cephalochordata) and Ciona robusta (Tunicata), with the aim to investigate the ancestral and the derived expression state of Magoh orthologous genes. This study identifies that Magoh is part of a conserved syntenic group exclusively in vertebrates and suggests that Magoh has experienced duplication and loss events in mammals. During early development in amphioxus and ascidian, maternal contribution and zygotic expression of Magoh genes in various types of progenitor cells and tissues are consistent with the condition observed in other Bilateria. Later in development, we also show expression of Magoh in the brain of cephalochordate and ascidian larvae. Collectively, these results provide a basis to further define what functional role(s) Magoh exerted during nervous system development and evolution.

Into the bloom: Molecular response of pelagic tunicates to fluctuating food availability.

The planktonic tunicates appendicularians and thaliaceans are highly efficient filter feeders on a wide range of prey size including bacteria and have shorter generation times than any other marine grazers. These traits allow some tunicate species to reach high population densities and ensure their success in a favourable environment. However, there are still few studies focusing on which genes and gene pathways are associated with responses of pelagic tunicates to environmental variability. Herein, we present the effect of food availability increase on tunicate community and gene expression at the Marquesas Islands (South-East Pacific Ocean). By using data from the Tara Oceans expedition, we show that changes in phytoplankton density and composition trigger the success of a dominant larvacean species (an undescribed appendicularian). Transcriptional signature to the autotroph bloom suggests key functions in specific physiological processes, i.e., energy metabolism, muscle contraction, membrane trafficking, and proteostasis. The relative abundance of reverse transcription-related Pfams was lower at bloom conditions, suggesting a link with adaptive genetic diversity in tunicates in natural ecosystems. Downstream of the bloom, pelagic tunicates were outcompeted by copepods. Our work represents the first metaomics study of the biological effects of phytoplankton bloom on a key zooplankton taxon.

The evolutionary landscape of the Rab family in chordates.

Intracellular traffic amongst organelles represents a key feature for eukaryotes and is orchestrated principally by members of Rab family, the largest within Ras superfamily. Given that variations in Rab repertoire have been fundamental in animal diversification, we provided the most exhaustive survey regarding the Rab toolkit of chordates. Our findings reveal the existence of 42 metazoan conserved subfamilies exhibiting a univocal intron/exon structure preserved from cnidarians to vertebrates. Since the current view does not capture the Rab complexity, we propose a new Rab family classification in three distinct monophyletic clades. The Rab complement of chordates shows a dramatic diversification due to genome duplications and independent gene duplications and losses with sharp differences amongst cephalochordates, tunicates and gnathostome vertebrates. Strikingly, the analysis of the domain architecture of this family highlighted the existence of chimeric calcium-binding Rabs, which are animal novelties characterized by a complex evolutionary history in gnathostomes and whose role in cellular metabolism is obscure. This work provides novel insights in the knowledge of Rab family: our hypothesis is that chordates represent a hotspot of Rab variability, with many events of gene gains and losses impacting intracellular traffic capabilities. Our results help to elucidate the role of Rab members in the transport amongst endomembranes and shed light on intracellular traffic routes in vertebrates. Then, since the predominant role of Rabs in the molecular communication between different cellular districts, this study paves to way to comprehend inherited or acquired human disorders provoked by dysfunctions in Rab genes.

Diatom-Derived Polyunsaturated Aldehydes Activate Similar Cell Death Genes in Two Different Systems: Sea Urchin Embryos and Human Cells.

Programmed cell death, such as apoptosis and autophagy, are key processes that are activated early on during development, leading to remodelling in embryos and homeostasis in adult organisms. Genomic conservation of death factors has been largely investigated in the animal and plant kingdoms. In this study, we analysed, for the first time, the expression profile of 11 genes involved in apoptosis (extrinsic and intrinsic pathways) and autophagy in sea urchin Paracentrotus lividus embryos exposed to antiproliferative polyunsaturated aldehydes (PUAs), and we compared these results with those obtained on the human cell line A549 treated with the same molecules. We found that sea urchins and human cells activated, at the gene level, a similar cell death response to these compounds. Despite the evolutionary distance between sea urchins and humans, we observed that the activation of apoptotic and autophagic genes in response to cytotoxic compounds is a conserved process. These results give first insight on death mechanisms of P. lividus death mechanisms, also providing additional information for the use of this marine organism as a useful in vitro model for the study of cell death signalling pathways activated in response to chemical compounds.

Detection of long non-coding RNA homology, a comparative study on alignment and alignment-free metrics.

BACKGROUND: Long non-coding RNAs (lncRNAs) represent a novel class of non-coding RNAs having a crucial role in many biological processes. The identification of long non-coding homologs among different species is essential to investigate such roles in model organisms as homologous genes tend to retain similar molecular and biological functions. Alignment-based metrics are able to effectively capture the conservation of transcribed coding sequences and then the homology of protein coding genes. However, unlike protein coding genes the poor sequence conservation of long non-coding genes makes the identification of their homologs a challenging task. RESULTS: In this study we compare alignment-based and alignment-free string similarity metrics and look at promoter regions as a possible source of conserved information. We show that promoter regions encode relevant information for the conservation of long non-coding genes across species and that such information is better captured by alignment-free metrics. We perform a genome wide test of this hypothesis in human, mouse, and zebrafish. CONCLUSIONS: The obtained results persuaded us to postulate the new hypothesis that, unlike protein coding genes, long non-coding genes tend to preserve their regulatory machinery rather than their transcribed sequence. All datasets, scripts, and the prediction tools adopted in this study are available at https://github.com/bioinformatics-sannio/lncrna-homologs .

Evolution of selenophosphate synthetases: emergence and relocation of function through independent duplications and recurrent subfunctionalization.

Selenoproteins are proteins that incorporate selenocysteine (Sec), a nonstandard amino acid encoded by UGA, normally a stop codon. Sec synthesis requires the enzyme Selenophosphate synthetase (SPS or SelD), conserved in all prokaryotic and eukaryotic genomes encoding selenoproteins. Here, we study the evolutionary history of SPS genes, providing a map of selenoprotein function spanning the whole tree of life. SPS is itself a selenoprotein in many species, although functionally equivalent homologs that replace the Sec site with cysteine (Cys) are common. Many metazoans, however, possess SPS genes with substitutions other than Sec or Cys (collectively referred to as SPS1). Using complementation assays in fly mutants, we show that these genes share a common function, which appears to be distinct from the synthesis of selenophosphate carried out by the Sec- and Cys- SPS genes (termed SPS2), and unrelated to Sec synthesis. We show here that SPS1 genes originated through a number of independent gene duplications from an ancestral metazoan selenoprotein SPS2 gene that most likely already carried the SPS1 function. Thus, in SPS genes, parallel duplications and subsequent convergent subfunctionalization have resulted in the segregation to different loci of functions initially carried by a single gene. This evolutionary history constitutes a remarkable example of emergence and evolution of gene function, which we have been able to trace thanks to the singular features of SPS genes, wherein the amino acid at a single site determines unequivocally protein function and is intertwined to the evolutionary fate of the entire selenoproteome.

Zebrafish, a Novel Model System to Study Uremic Toxins: The Case for the Sulfur Amino Acid Lanthionine.

The non-proteinogenic amino acid lanthionine is a byproduct of hydrogen sulfide biosynthesis: the third endogenous vasodilator gas, after nitric oxide and carbon monoxide. While hydrogen sulfide is decreased in uremic patients on hemodialysis, lanthionine is increased and has been proposed as a new uremic toxin, since it is able to impair hydrogen sulfide production in hepatoma cells. To characterize lanthionine as a uremic toxin, we explored its effects during the early development of the zebrafish (Danio rerio), a widely used model to study the organ and tissue alterations induced by xenobiotics. Lanthionine was employed at concentrations reproducing those previously detected in uremia. Light-induced visual motor response was also studied by means of the DanioVision system. Treatment of zebrafish embryos with lanthionine determined acute phenotypical alterations, on heart organogenesis (disproportion in cardiac chambers), increased heart beating, and arrhythmia. Lanthionine also induced locomotor alterations in zebrafish embryos. Some of these effects could be counteracted by glutathione. Lanthionine exerted acute effects on transsulfuration enzymes and the expression of genes involved in inflammation and metabolic regulation, and modified microRNA expression in a way comparable with some alterations detected in uremia. Lanthionine meets the criteria for classification as a uremic toxin. Zebrafish can be successfully used to explore uremic toxin effects.

Identification of long non-coding transcripts with feature selection: a comparative study.

BACKGROUND: The unveiling of long non-coding RNAs as important gene regulators in many biological contexts has increased the demand for efficient and robust computational methods to identify novel long non-coding RNAs from transcripts assembled with high throughput RNA-seq data. Several classes of sequence-based features have been proposed to distinguish between coding and non-coding transcripts. Among them, open reading frame, conservation scores, nucleotide arrangements, and RNA secondary structure have been used with success in literature to recognize intergenic long non-coding RNAs, a particular subclass of non-coding RNAs. RESULTS: In this paper we perform a systematic assessment of a wide collection of features extracted from sequence data. We use most of the features proposed in the literature, and we include, as a novel set of features, the occurrence of repeats contained in transposable elements. The aim is to detect signatures (groups of features) able to distinguish long non-coding transcripts from other classes, both protein-coding and non-coding. We evaluate different feature selection algorithms, test for signature stability, and evaluate the prediction ability of a signature with a machine learning algorithm. The study reveals different signatures in human, mouse, and zebrafish, highlighting that some features are shared among species, while others tend to be species-specific. Compared to coding potential tools and similar supervised approaches, including novel signatures, such as those identified here, in a machine learning algorithm improves the prediction performance, in terms of area under precision and recall curve, by 1 to 24%, depending on the species and on the signature. CONCLUSIONS: Understanding which features are best suited for the prediction of long non-coding RNAs allows for the development of more effective automatic annotation pipelines especially relevant for poorly annotated genomes, such as zebrafish. We provide a web tool that recognizes novel long non-coding RNAs with the obtained signatures from fasta and gtf formats. The tool is available at the following url: http://www.bioinformatics-sannio.org/software/ .

Rab32 and Rab38 genes in chordate pigmentation: an evolutionary perspective.

BACKGROUND: The regulation of cellular membrane trafficking in all eukaryotes is a very complex mechanism, mostly regulated by the Rab family proteins. Among all membrane-enclosed organelles, melanosomes are the cellular site for synthesis, storage and transport of melanin granules, making them an excellent model for studies on organelle biogenesis and motility. Specific Rab proteins, as Rab32 and Rab38, have been shown to play a key role in melanosome biogenesis. We analysed the Rab32 and Rab38 genes in the teleost zebrafish and in the cephalochordate amphioxus, gaining insight on their evolutionary history following gene and genome duplications. RESULTS: We studied the molecular evolution of Rab supergroup III in deuterostomes by phylogenetic reconstruction, intron and synteny conservation. We discovered a novel amino acid stretch, named FALK, shared by three related classes belonging to Rab supergroup III: Rab7L1, Rab32LO and Rab32/Rab38. Among these, we demonstrated that the Rab32LO class, already present in the last common eukaryotic ancestor, was lost in urochordates and vertebrates. Synteny shows that one zebrafish gene, Rab38a, which is expressed in pigmented cells, retained the linkage with tyrosinase, a protein essential for pigmentation. Moreover, the chromosomal linkage of Rab32 or Rab38 with a member of the glutamate receptor metabotropic (Grm) family has been retained in all analysed gnathostomes, suggesting a conserved microsynteny in the vertebrate ancestor. Expression patterns of Rab32 and Rab38 genes in zebrafish, and Rab32/38 in amphioxus, indicate their involvement in development of pigmented cells and notochord. CONCLUSIONS: Phylogenetic, intron conservation and synteny analyses point towards an evolutionary scenario based on a duplication of a single invertebrate Rab32/38 gene giving rise to vertebrate Rab32 and Rab38. The expression patterns of Rab38 paralogues highlight sub-functionalization event. Finally, the discovery of a chromosomal linkage between the Rab32 or Rab38 gene with a Grm opens new perspectives on possible conserved bystander gene regulation across the vertebrate evolution.