Evidence for Non-neutral Evolution in a Sodium Channel Gene in African Weakly Electric Fish (Campylomormyrus, Mormyridae)

Gene and Allele-Specific Expression Underlying the Electric Signal Divergence in African Weakly Electric Fish

In the African weakly electric fish genus Campylomormyrus, electric organ discharge signals are strikingly different in shape and duration among closely related species, contribute to prezygotic isolation, and may have triggered an adaptive radiation. We performed mRNA sequencing on electric organs and skeletal muscles (from which the electric organs derive) from 3 species with short (0.4 ms), medium (5 ms), and long (40 ms) electric organ discharges and 2 different cross-species hybrids. We identified 1,444 upregulated genes in electric organ shared by all 5 species/hybrid cohorts, rendering them candidate genes for electric organ-specific properties in Campylomormyrus. We further identified several candidate genes, including KCNJ2 and KLF5, and their upregulation may contribute to increased electric organ discharge duration. Hybrids between a short (Campylomormyrus compressirostris) and a long (Campylomormyrus rhynchophorus) discharging species exhibit electric organ discharges of intermediate duration and showed imbalanced expression of KCNJ2 alleles, pointing toward a cis-regulatory difference at this locus, relative to electric organ discharge duration. KLF5 is a transcription factor potentially balancing potassium channel gene expression, a crucial process for the formation of an electric organ discharge. Unraveling the genetic basis of the species-specific modulation of the electric organ discharge in Campylomormyrus is crucial for understanding the adaptive radiation of this emerging model taxon of ecological (perhaps even sympatric) speciation.

A new genome assembly of an African weakly electric fish (Campylomormyrus compressirostris, Mormyridae) indicates rapid gene family evolution in Osteoglossomorpha

BACKGROUND

Teleost fishes comprise more than half of the vertebrate species. Within teleosts, most phylogenies consider the split between Osteoglossomorpha and Euteleosteomorpha/Otomorpha as basal, preceded only by the derivation of the most primitive group of teleosts, the Elopomorpha. While Osteoglossomorpha are generally species poor, the taxon contains the African weakly electric fish (Mormyroidei), which have radiated into numerous species. Within the mormyrids, the genus Campylomormyrus is mostly endemic to the Congo Basin. Campylomormyrus serves as a model to understand mechanisms of adaptive radiation and ecological speciation, especially with regard to its highly diverse species-specific electric organ discharges (EOD). Currently, there are few well-annotated genomes available for electric fish in general and mormyrids in particular. Our study aims at producing a high-quality genome assembly and to use this to examine genome evolution in relation to other teleosts. This will facilitate further understanding of the evolution of the osteoglossomorpha fish in general and of electric fish in particular.

RESULTS

A high-quality weakly electric fish (C. compressirostris) genome was produced from a single individual with a genome size of 862 Mb, consisting of 1,497 contigs with an N50 of 1,399 kb and a GC-content of 43.69%. Gene predictions identified 34,492 protein-coding genes, which is a higher number than in the two other available Osteoglossomorpha genomes of Paramormyrops kingsleyae and Scleropages formosus. A Computational Analysis of gene Family Evolution (CAFE5) comparing 33 teleost fish genomes suggests an overall faster gene family turnover rate in Osteoglossomorpha than in Otomorpha and Euteleosteomorpha. Moreover, the ratios of expanded/contracted gene family numbers in Osteoglossomorpha are significantly higher than in the other two taxa, except for species that had undergone an additional genome duplication (Cyprinus carpio and Oncorhynchus mykiss). As potassium channel proteins are hypothesized to play a key role in EOD diversity among species, we put a special focus on them, and manually curated 16 Kv1 genes. We identified a tandem duplication in the KCNA7a gene in the genome of C. compressirostris.

CONCLUSIONS

We present the fourth genome of an electric fish and the third well-annotated genome for Osteoglossomorpha, enabling us to compare gene family evolution among major teleost lineages. Osteoglossomorpha appear to exhibit rapid gene family evolution, with more gene family expansions than contractions. The curated Kv1 gene family showed seven gene clusters, which is more than in other analyzed fish genomes outside Osteoglossomorpha. The KCNA7a, encoding for a potassium channel central for EOD production and modulation, is tandemly duplicated which may related to the diverse EOD observed among Campylomormyrus species.

Molecular evolution of the ependymin-related gene epdl2 in African weakly electric fish

Gene duplication and subsequent molecular evolution can give rise to taxon-specific gene specializations. In previous work, we found evidence that African weakly electric fish (Mormyridae) may have as many as three copies of the epdl2 gene, and the expression of two epdl2 genes is correlated with electric signal divergence. Epdl2 belongs to the ependymin-related family (EPDR), a functionally diverse family of secretory glycoproteins. In this study, we first describe vertebrate EPDR evolution and then present a detailed evolutionary history of epdl2 in Mormyridae with emphasis on the speciose genus Paramormyrops. Using Sanger sequencing, we confirm three apparently functional epdl2 genes in Paramormyrops kingsleyae. Next, we developed a nanopore-based amplicon sequencing strategy and bioinformatics pipeline to obtain and classify full-length epdl2 gene sequences (N = 34) across Mormyridae. Our phylogenetic analysis proposes three or four epdl2 paralogs dating from early Paramormyrops evolution. Finally, we conducted selection tests which detected positive selection around the duplication events and identified ten sites likely targeted by selection in the resulting paralogs. These sites' locations in our modeled 3D protein structure involve four sites in ligand binding and six sites in homodimer formation. Together, these findings strongly imply an evolutionary mechanism whereby epdl2 genes underwent selection-driven functional specialization after tandem duplications in the rapidly speciating Paramormyrops. Considering previous evidence, we propose that epdl2 may contribute to electric signal diversification in mormyrids, an important aspect of species recognition during mating.

Divergent cis-regulatory evolution underlies the convergent loss of sodium channel expression in electric fish

South American and African weakly electric fish independently evolved electric organs from muscle. In both groups, a voltage-gated sodium channel gene independently lost expression from muscle and gained it in the electric organ, allowing the channel to become specialized for generating electric signals. It is unknown how this voltage-gated sodium channel gene is targeted to muscle in any vertebrate. We describe an enhancer that selectively targets sodium channel expression to muscle. Next, we demonstrate how the loss of this enhancer, but not trans-activating factors, drove the loss of sodium channel gene expression from muscle in South American electric fish. While this enhancer is also altered in African electric fish, key transcription factor binding sites and enhancer activity are retained, suggesting that the convergent loss of sodium channel expression from muscle in these two electric fish lineages occurred via different processes.

Vocal and Electric Fish: Revisiting a Comparison of Two Teleost Models in the Neuroethology of Social Behavior

The communication behaviors of vocal fish and electric fish are among the vertebrate social behaviors best understood at the level of neural circuits. Both forms of signaling rely on midbrain inputs to hindbrain pattern generators that activate peripheral effectors (sonic muscles and electrocytes) to produce pulsatile signals that are modulated by frequency/repetition rate, amplitude and call duration. To generate signals that vary by sex, male phenotype, and social context, these circuits are responsive to a wide range of hormones and neuromodulators acting on different timescales at multiple loci. Bass and Zakon (2005) reviewed the behavioral neuroendocrinology of these two teleost groups, comparing how the regulation of their communication systems have both converged and diverged during their parallel evolution. Here, we revisit this comparison and review the complementary developments over the past 16 years. We (a) summarize recent work that expands our knowledge of the neural circuits underlying these two communication systems, (b) review parallel studies on the action of neuromodulators (e.g., serotonin, AVT, melatonin), brain steroidogenesis (via aromatase), and social stimuli on the output of these circuits, (c) highlight recent transcriptomic studies that illustrate how contemporary molecular methods have elucidated the genetic regulation of social behavior in these fish, and (d) describe recent studies of mochokid catfish, which use both vocal and electric communication, and that use both vocal and electric communication and consider how these two systems are spliced together in the same species. Finally, we offer avenues for future research to further probe how similarities and differences between these two communication systems emerge over ontogeny and evolution.

Transcriptome-wide single nucleotide polymorphisms related to electric organ discharge differentiation among African weakly electric fish species

African weakly electric fish of the mormyrid genus Campylomormyrus generate pulse-type electric organ discharges (EODs) for orientation and communication. Their pulse durations are species-specific and elongated EODs are a derived trait. So far, differential gene expression among tissue-specific transcriptomes across species with different pulses and point mutations in single ion channel genes indicate a relation of pulse duration and electrocyte geometry/excitability. However, a comprehensive assessment of expressed Single Nucleotide Polymorphisms (SNPs) throughout the entire transcriptome of African weakly electric fish, with the potential to identify further genes influencing EOD duration, is still lacking. This is of particular value, as discharge duration is likely based on multiple cellular mechanisms and various genes. Here we provide the first transcriptome-wide SNP analysis of African weakly electric fish species (genus Campylomormyrus) differing by EOD duration to identify candidate genes and cellular mechanisms potentially involved in the determination of an elongated discharge of C. tshokwe. Non-synonymous substitutions specific to C. tshokwe were found in 27 candidate genes with inferred positive selection among Campylomormyrus species. These candidate genes had mainly functions linked to transcriptional regulation, cell proliferation and cell differentiation. Further, by comparing gene annotations between C. compressirostris (ancestral short EOD) and C. tshokwe (derived elongated EOD), we identified 27 GO terms and 2 KEGG pathway categories for which C. tshokwe significantly more frequently exhibited a species-specific expressed substitution than C. compressirostris. The results indicate that transcriptional regulation as well cell proliferation and differentiation take part in the determination of elongated pulse durations in C. tshokwe. Those cellular processes are pivotal for tissue morphogenesis and might determine the shape of electric organs supporting the observed correlation between electrocyte geometry/tissue structure and discharge duration. The inferred expressed SNPs and their functional implications are a valuable resource for future investigations on EOD durations.

The transcriptional correlates of divergent electric organ discharges in Paramormyrops electric fish

Background

Understanding the genomic basis of phenotypic diversity can be greatly facilitated by examining adaptive radiations with hypervariable traits. In this study, we focus on a rapidly diverged species group of mormyrid electric fish in the genus Paramormyrops, which are characterized by extensive phenotypic variation in electric organ discharges (EODs). The main components of EOD diversity are waveform duration, complexity and polarity. Using an RNA-sequencing based approach, we sought to identify gene expression correlates for each of these EOD waveform features by comparing 11 specimens of Paramormyrops that exhibit variation in these features.

Results

Patterns of gene expression among Paramormyrops are highly correlated, and 3274 genes (16%) were differentially expressed. Using our most restrictive criteria, we detected 145–183 differentially expressed genes correlated with each EOD feature, with little overlap between them. The predicted functions of several of these genes are related to extracellular matrix, cation homeostasis, lipid metabolism, and cytoskeletal and sarcomeric proteins. These genes are of significant interest given the known morphological differences between electric organs that underlie differences in the EOD waveform features studied.

Conclusions

In this study, we identified plausible candidate genes that may contribute to phenotypic differences in EOD waveforms among a rapidly diverged group of mormyrid electric fish. These genes may be important targets of selection in the evolution of species-specific differences in mate-recognition signals.

Electrostatic Tuning of a Potassium Channel in Electric Fish

Molecular variation contributes to the evolution of adaptive phenotypes, though it is often difficult to understand precisely how. The adaptively significant electric organ discharge behavior of weakly electric fish is the direct result of biophysical membrane properties set by ion channels. Here, we describe a voltage-gated potassium-channel gene in African electric fishes that is under positive selection and highly expressed in the electric organ. The channel produced by this gene shortens electric organ action potentials by activating quickly and at hyperpolarized membrane potentials. The source of these properties is a derived patch of negatively charged amino acids in an extracellular loop near the voltage sensor. We demonstrate that this negative patch acts by contributing to the global surface charge rather than by local interactions with specific amino acids in the channel's extracellular face. We suggest a more widespread role for this loop in the evolutionary tuning of voltage-dependent channels.

Intragenus (Campylomormyrus) and intergenus hybrids in mormyrid fish: Physiological and histological investigations of the electric organ ontogeny

African weakly electric mormyrid fish show a high diversity of their electric organ discharge (EOD) both across and within genera. Thanks to a recently developed technique of artificial reproduction in mormyrid fish, we were able to perform hybridizations between different genera and within one genus (Campylomormyrus). The hybrids of intergenus hybridizations exhibited different degrees of reduced survival related to the phylogenetic distance of the parent species: hybrids of the crosses between C. rhynchophorus and its sister genus Gnathonemus survived and developed normally. Hybrids between C. rhynchophorus and a Mormyrus species (a more basal clade compared to Campylomormyrus s) survived up to 42days and developed many malformations, e.g., at the level of the unpaired fins. Hybrids between C. numenius and Hippopotamyrus pictus (a derived clade, only distantly related to Campylomormyrus) only survived for two days during embryological development. Eight different hybrid combinations among five Campylomormyrus species (C. tamandua, C. compressirostris, C. tshokwe, C. rhynchophorus, C. numenius) were performed. The aim of the hybridizations was to combine species with (1) either caudal or rostral position of the main stalk innervating the electrocytes in the electric organ and (2) short, median or long duration of their EOD. The hybrids, though they are still juveniles, show very interesting features concerning electrocyte geometry as well as EOD form and duration: the caudal position of the stalk is prevailing over the rostral position, and the penetration of the stalk is dominant over the non-penetrating feature (in the Campylomormyrus hybrids); in the hybrid between C. rhynchophorus and Gnathonemus petersii it is the opposite. When crossing species with long and short EODs, it is always the long duration EOD that is expressed in the hybrids. The F1-Hybrids of the cross C. tamandua×C. compressirostris are fertile: viable F2-fish could be obtained with artificial reproduction.

Voltage-Gated Na+ Channel Isoforms and Their mRNA Expression Levels and Protein Abundance in Three Electric Organs and the Skeletal Muscle of the Electric Eel Electrophorus electricus

This study aimed to obtain the coding cDNA sequences of voltage-gated Na+ channel (scn) α-subunit (scna) and β-subunit (scnb) isoforms from, and to quantify their transcript levels in, the main electric organ (EO), Hunter's EO, Sach's EO and the skeletal muscle (SM) of the electric eel, Electrophorus electricus, which can generate both high and low voltage electric organ discharges (EODs). The full coding sequences of two scna (scn4aa and scn4ab) and three scnb (scn1b, scn2b and scn4b) were identified for the first time (except scn4aa) in E. electricus. In adult fish, the scn4aa transcript level was the highest in the main EO and the lowest in the Sach's EO, indicating that it might play an important role in generating high voltage EODs. For scn4ab/Scn4ab, the transcript and protein levels were unexpectedly high in the EOs, with expression levels in the main EO and the Hunter's EO comparable to those of scn4aa. As the key domains affecting the properties of the channel were mostly conserved between Scn4aa and Scn4ab, Scn4ab might play a role in electrogenesis. Concerning scnb, the transcript level of scn4b was much higher than those of scn1b and scn2b in the EOs and the SM. While the transcript level of scn4b was the highest in the main EO, protein abundance of Scn4b was the highest in the SM. Taken together, it is unlikely that Scna could function independently to generate EODs in the EOs as previously suggested. It is probable that different combinations of Scn4aa/Scn4ab and various Scnb isoforms in the three EOs account for the differences in EODs produced in E. electricus. In general, the transcript levels of various scn isoforms in the EOs and the SM were much higher in adult than in juvenile, and the three EOs of the juvenile fish could be functionally indistinct.