Cross-tissue and cross-species analysis of gene expression in skeletal muscle and electric organ of African weakly-electric fish (Teleostei; Mormyridae). BMC Genomics. 2015;16:668. [PMID: 26335922 DOI: 10.1186/s12864-015-1858-9]

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 chromosome-level genome of electric catfish (Malapterurus electricus) provided new insights into order Siluriformes evolution

The electric catfish (Malapterurus electricus), belonging to the family Malapteruridae, order Siluriformes (Actinopterygii: Ostariophysi), is one of the six branches that has independently evolved electrical organs. We assembled a 796.75 Mb M. electricus genome and anchored 88.72% sequences into 28 chromosomes. Gene family analysis revealed 295 expanded gene families that were enriched on functions related to glutamate receptors. Convergent evolutionary analyses of electric organs among different lineage of electric fishes further revealed that the coding gene of rho guanine nucleotide exchange factor 4-like (arhgef4), which is associated with G-protein coupled receptor (GPCR) signaling pathway, underwent adaptive parallel evolution. Gene identification suggests visual degradation in catfishes, and an important role for taste in environmental adaptation. Our findings fill in the genomic data for a branch of electric fish and provide a relevant genetic basis for the adaptive evolution of Siluriformes.

Supplementary Information

The online version contains supplementary material available at 10.1007/s42995-023-00197-8.

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.

Phylogenomics of bonytongue fishes (Osteoglossomorpha) shed light on the craniofacial evolution and biogeography of the weakly electric clade Mormyridae

Bonytongues (Osteoglossomorpha) constitute an ancient clade of teleost fishes distributed in freshwater habitats throughout the world. The group includes well-known species such as arowanas, featherbacks, pirarucus, and the weakly electric fishes in the family Mormyridae. Their disjunct distribution, extreme morphologies, and electrolocating capabilities (Gymnarchidae and Mormyridae) have attracted much scientific interest, but a comprehensive phylogenetic framework for comparative analysis is missing, especially for the species-rich family Mormyridae. Of particular interest are disparate craniofacial morphologies among mormyrids which might constitute an exceptional model system to study convergent evolution. We present a phylogenomic analysis based on 546 exons of 179 species (out of 260), 28 out of 29 genera, and all six families of extant bonytongues. Based on a recent reassessment of the fossil record of osteoglossomorphs, we inferred dates of divergence among trans-continental clades and the major groups. The estimated ages of divergence among extant taxa (e.g., Osteoglossomorpha, Osteoglossiformes, Mormyroidea) are older than previous reports, but most of the divergence dates obtained for clades on separate continents are too young to be explained by simple vicariance hypotheses. Biogeographic analysis of mormyrids indicates that their high species diversity in the Congo Basin is a consequence of range reductions of previously widespread ancestors and that the highest diversity of craniofacial morphologies among mormyrids originated in this basin. Special emphasis on a taxon-rich representation for mormyrids revealed pervasive misalignment between our phylogenomic results and mormyrid taxonomy due to repeated instances of convergence for extreme craniofacial morphologies. Estimation of ancestral phenotypes revealed contingent evolution of snout elongation and unique projections from the lower jaw to form the distinctive Schnauzenorgan. Synthesis of comparative analyses suggests that the remarkable craniofacial morphologies of mormyrids evolved convergently due to niche partitioning, likely enabled by interactions between their exclusive morphological and electrosensory adaptations.

Genomic Evidence for Convergent Molecular Adaptation in Electric Fishes

Fishes have independently evolved electric organs (EOs) at least six times, and the electric fields are used for communication, defense, and predation. However, the genetic basis of convergent evolution of EOs remains unclear. In this study, we conducted comparative genomic analyses to detect genes showing signatures of positive selection and convergent substitutions in electric fishes from three independent lineages (Mormyroidea, Siluriformes, and Gymnotiformes). Analysis of 4,657 orthologs between electric fishes and their corresponding control groups identified consistent evidence for accelerated evolution in electric fish lineages. A total of 702 positively selected genes (PSGs) were identified in electric fishes, and many of these genes corresponded to cell membrane structure, ion channels, and transmembrane transporter activity. Comparative genomic analyses revealed that widespread convergent amino acid substitutions occurred along the electric fish lineages. The overlap of convergent genes and PSGs was identified as adaptive convergence, and a subset of genes was putatively associated with electrical and muscular activities, especially scn4aa (a voltage-gated sodium channel gene). Our results provide hints to the genetic basis for the independent evolution of EOs during millions of years of 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.

Comparison of spleen transcriptomes of two wild rodent species reveals differences in the immune response against Borrelia afzelii

Different host species often differ considerably in susceptibility to a given pathogen, but the causes of such differences are rarely known. The natural hosts of the tick-transmitted bacterium Borrelia afzelii, which is one of causative agents of Lyme borreliosis in humans, include a variety of small mammals like voles and mice. Previous studies have shown that B. afzelii-infected bank voles (Myodes glareolus) have about ten times higher bacterial load than infected yellow-necked mice (Apodemus flavicollis), indicating that these two species differ in resistance. In this study, we compared the immune response to B. afzelii infection in these host species by using RNA sequencing to quantify gene expression in spleen. Gene set enrichment analysis (GSEA) showed that several immune pathways were down-regulated in infected animals in both bank voles and yellow-necked mice. Moreover, IFNα response was up-regulated in B. afzelii-infected yellow-necked mice, while IL6 signaling and the complement pathway were down-regulated in infected bank voles; differences in regulation of these three pathways between bank voles and yellow-necked mice could thus contribute to the difference in resistance to B. afzelii between the species. This study provides knowledge of gene expression induced by a zoonotic pathogen in its natural host, and possible species-specific regulation of immune responses associated with resistance.

Genomic Survey of Tyrosine Kinases Repertoire in Electrophorus electricus With an Emphasis on Evolutionary Conservation and Diversification

Tyrosine kinases (TKs) play key roles in the regulation of multicellularity in organisms and involved primarily in cell growth, differentiation, and cell-to-cell communication. Genome-wide characterization of TKs has been conducted in many metazoans; however, systematic information regarding this superfamily in Electrophorus electricus (electric eel) is still lacking. In this study, we identified 114 TK genes in the E electricus genome and investigated their evolution, molecular features, and domain architecture using phylogenetic profiling to gain a better understanding of their similarities and specificity. Our results suggested that the electric eel TK (EeTK) repertoire was shaped by whole-genome duplications (WGDs) and tandem duplication events. Compared with other vertebrate TKs, gene members in Jak, Src, and EGFR subfamily duplicated specifically, but with members lost in Eph, Axl, and Ack subfamily in electric eel. We also conducted an exhaustive survey of TK genes in genomic databases, identifying 1674 TK proteins in 31 representative species covering all the main metazoan lineages. Extensive evolutionary analysis indicated that TK repertoire in vertebrates tended to be remarkably conserved, but the gene members in each subfamily were very variable. Comparative expression profile analysis showed that electric organ tissues and muscle shared a similar pattern with specific highly expressed TKs (ie, epha7, musk, jak1, and pdgfra), suggesting that regulation of TKs might play an important role in specifying an electric organ identity from its muscle precursor. We further identified TK genes exhibiting tissue-specific expression patterns, indicating that members in TKs participated in subfunctionalization representing an evolutionary divergence required for the performance of different tissues. This work generates valuable information for further gene function analysis and identifying candidate TK genes reflecting their unique tissue-function specializations in electric eel.

Intragenus F1-hybrids of African weakly electric fish (Mormyridae: Campylomormyrus tamandua ♂ × C. compressirostris ♀) are fertile

Hybridization is widespread in fish and constitutes an important mechanism in fish speciation. There is, however, little knowledge about hybridization in mormyrids. F1-interspecies hybrids between Campylomormyrus tamandua ♂ × C. compressirostris ♀ were investigated concerning: (1) fertility; (2) survival of F2-fish and (3) new gene combinations in the F2-generation concerning the structure of the electric organ and features of the electric organ discharge. These F1-hybrids achieved sexual maturity at about 12–13.5 cm total length. A breeding group comprising six males and 13 females spawned 28 times naturally proving these F1-fish to be fertile. On average 228 eggs were spawned, the average fertilization rate was 47.8%. Eggs started to hatch 70–72 h after fertilization, average hatching rate was 95.6%. Average mortality rate during embryonic development amounted to 2.3%. Average malformation rate during the free embryonic stage was 27.7%. Exogenous feeding started on day 11. In total, we raised 353 normally developed larvae all of which died consecutively, the oldest specimen reaching an age of 5 months. During survival, the activities of the larval and adult electric organs were recorded and the structure of the adult electric organ was investigated histologically.

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.

Contribution of Extracellular Vesicles in Rebuilding Injured Muscles

Skeletal myofibers are injured due to mechanical stresses experienced during physical activity, or due to myofiber fragility caused by genetic diseases. The injured myofiber needs to be repaired or regenerated to restore the loss in muscle tissue function. Myofiber repair and regeneration requires coordinated action of various intercellular signaling factors-including proteins, inflammatory cytokines, miRNAs, and membrane lipids. It is increasingly being recognized release and transmission of these signaling factors involves extracellular vesicle (EV) released by myofibers and other cells in the injured muscle. Intercellular signaling by these EVs alters the phenotype of their target cells either by directly delivering the functional proteins and lipids or by modifying longer-term gene expression. These changes in the target cells activate downstream pathways involved in tissue homeostasis and repair. The EVs are heterogeneous with regards to their size, composition, cargo, location, as well as time-course of genesis and release. These differences impact on the subsequent repair and regeneration of injured skeletal muscles. This review focuses on how intracellular vesicle production, cargo packaging, and secretion by injured muscle, modulates specific reparative, and regenerative processes. Insights into the formation of these vesicles and their signaling properties offer new understandings of the orchestrated response necessary for optimal muscle repair and regeneration.

Comparative Transcriptome Analyses Reveal the Role of Conserved Function in Electric Organ Convergence Across Electric Fishes

The independent origins of multiple electric organs (EOs) of fish are fascinating examples of convergent evolution. However, comparative transcriptomics of different electric fish lineages are scarce. In this study, we found that the gene expression of EOs and skeletal muscles from three lineages (Mormyroidea, Siluriformes, and Gymnotiformes) tended to cluster together based on the species of origin, irrespective of the organ from which they are derived. A pairwise comparison of differentially expressed genes (DEGs) revealed that no less than half of shared DEGs exhibited parallel expression differentiation, indicating conserved directionality of differential expression either in or between lineages, but only a few shared DEGs were identified across all focal species. Nevertheless, the functional enrichment analysis of DEGs indicated that there were more parallel gene expression changes at the level of pathways and biological functions. Therefore, we may conclude that there is no parallel evolution of the entire transcriptomes of EOs among different lineages. Further, our results support the hypothesis that it is not different genes but conserved biological functions that play a crucial role in the convergence of complex phenotypes. This study provides insight into the genetic basis underlying the EO convergent evolution; however, more studies in different cases will be needed to demonstrate whether this pattern can be extended to other cases to derive a general rule for convergent evolution.

Electroreception, electrogenesis and electric signal evolution

Electroreception, the capacity to detect external underwater electric fields with specialised receptors, is a phylogenetically widespread sensory modality in fishes and amphibians. In passive electroreception, a capacity possessed by c. 16% of fish species, an animal uses low-frequency-tuned ampullary electroreceptors to detect microvolt-range bioelectric fields from prey, without the need to generate its own electric field. In active electroreception (electrolocation), which occurs only in the teleost lineages Mormyroidea and Gymnotiformes, an animal senses its surroundings by generating a weak (< 1 V) electric-organ discharge (EOD) and detecting distortions in the EOD-associated field using high-frequency-tuned tuberous electroreceptors. Tuberous electroreceptors also detect the EODs of neighbouring fishes, facilitating electrocommunication. Several other groups of elasmobranchs and teleosts generate weak (< 10 V) or strong (> 50 V) EODs that facilitate communication or predation, but not electrolocation. Approximately 1.5% of fish species possess electric organs. This review has two aims. First, to synthesise our knowledge of the functional biology and phylogenetic distribution of electroreception and electrogenesis in fishes, with a focus on freshwater taxa and with emphasis on the proximate (morphological, physiological and genetic) bases of EOD and electroreceptor diversity. Second, to describe the diversity, biogeography, ecology and electric signal diversity of the mormyroids and gymnotiforms and to explore the ultimate (evolutionary) bases of signal and receptor diversity in their convergent electrogenic-electrosensory systems. Four sets of potential drivers or moderators of signal diversity are discussed. First, selective forces of an abiotic (environmental) nature for optimal electrolocation and communication performance of the EOD. Second, selective forces of a biotic nature targeting the communication function of the EOD, including sexual selection, reproductive interference from syntopic heterospecifics and selection from eavesdropping predators. Third, non-adaptive drift and, finally, phylogenetic inertia, which may arise from stabilising selection for optimal signal-receptor matching.

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.

The Genome and Adult Somatic Transcriptome of the Mormyrid Electric Fish Paramormyrops kingsleyae

Several studies have begun to elucidate the genetic and developmental processes underlying major vertebrate traits. Few of these traits have evolved repeatedly in vertebrates, preventing the analysis of molecular mechanisms underlying these traits comparatively. Electric organs have evolved multiple times among vertebrates, presenting a unique opportunity to understand the degree of constraint and repeatability of the evolutionary processes underlying novel vertebrate traits. As there is now a completed genome sequence representing South American electric eels, we were motivated to obtain genomic sequence from a linage that independently evolved electric organs to facilitate future comparative analyses of the evolution and development of electric organs. We report here the sequencing and de novo assembly of the genome of the mormyrid Paramormyrops kingsleyae using short-read sequencing. In addition, we have completed a somatic transcriptome from 11 tissues to construct a gene expression atlas of predicted genes from this assembly, enabling us to identify candidate housekeeping genes as well as genes differentially expressed in the major somatic tissues of the mormyrid electric fish. We anticipate that this resource will greatly facilitate comparative studies on the evolution and development of electric organs and electroreceptors.

Karyotype description of the African weakly electric fish Campylomormyrus compressirostris in the context of chromosome evolution in Osteoglossiformes

Karyotyping is a basic method to investigate chromosomal evolution and genomic rearrangements. Sixteen genera within the basal teleost order Osteoglossiformes are currently described cytogenetically. Our study adds information to this chromosomal dataset by determining the karyotype of Campylomormyrus compressirostris, a genus of African weakly electric fish that has not been previously examined. Our results indicate a diploid chromosome number of 2n=48 (4sm+26m+18a) with a fundamental number of FN=72. This chromosome number is identical to the number documented for the sister taxon of the genus Campylomormyrus, i.e., Gnathonemus petersii (2n=48). These results support the close relationship of Campylomormyrus and Gnathonemus. However, the karyotype formula of C. compressirostris is different from Gnathonemus petersii, thereby confirming the high variability of karyotype formulae within the Mormyridae. We infer that the differences in chromosome number and formula of Campylomormyrus relative to other mormyrids may be caused by Robertsonian fusion and pericentric inversion. In addition to the karyotype description and classification of Campylomormyrus, a ChromEvol analysis was used to determine the ancestral haploid chromosome number of osteoglossiform taxa. Our results indicate a relatively conservative haploid chromosome number of n=24 for the most recent common ancestor of Osteoglossiformes and for most of the internal nodes of osteoglossiform phylogeny. Hence, we presume that the high chromosome variability evolved recently on multiple independent occasions. Furthermore, we suggest that the most likely ancestral chromosome number of Mormyridae is either n=24 or n=25. To the best of our knowledge this is the first attempt to determine and classify the karyotype of the weakly electric fish genus Campylomormyrus and to analyze chromosomal evolution within the Osteoglossiformes based on Maximum Likelihood and Bayesian Inference analyses.

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.

Electric fish genomics: Progress, prospects, and new tools for neuroethology

Electric fish have served as a model system in biology since the 18th century, providing deep insight into the nature of bioelectrogenesis, the molecular structure of the synapse, and brain circuitry underlying complex behavior. Neuroethologists have collected extensive phenotypic data that span biological levels of analysis from molecules to ecosystems. This phenotypic data, together with genomic resources obtained over the past decades, have motivated new and exciting hypotheses that position the weakly electric fish model to address fundamental 21^st century biological questions. This review article considers the molecular data collected for weakly electric fish over the past three decades, and the insights that data of this nature has motivated. For readers relatively new to molecular genetics techniques, we also provide a table of terminology aimed at clarifying the numerous acronyms and techniques that accompany this field. Next, we pose a research agenda for expanding genomic resources for electric fish research over the next 10years. We conclude by considering some of the exciting research prospects for neuroethology that electric fish genomics may offer over the coming decades, if the electric fish community is successful in these endeavors.

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

Voltage-gated sodium channels, Nav1, play a crucial role in the generation and propagation of action potentials and substantially contribute to the shape of their rising phase. The electric organ discharge (EOD) of African weakly electric fish (Mormyroidea) is the sum of action potentials fired from all electrocytes of the electric organ at the same time and hence voltage-gated sodium channels are one factor-together with the electrocyte's morphology and innervation pattern-that determines the properties of these EODs. Due to the fish-specific genome duplication, teleost fish possess eight copies of sodium channel genes (SCN), which encode for Nav1 channels. In mormyroids, SCN4aa is solely expressed in the electrocytes of the adult electric organ. In this study, we compared entire SCN4aa sequences of six species of the genus Campylomormyrus and identified nonsynonymous substitutions among them. SCN4aa in Campylomormyrus exhibits a much higher evolutionary rate compared to its paralog SCN4ab, whose expression is not restricted to the electric organ. We also found evidence for strong positive selection on the SCN4aa gene within Mormyridae and along the lineage ancestral to the Mormyridae. We have identified sites at which all nonelectric teleosts are monomorphic in their amino acid, but mormyrids have different amino acids. Our findings confirm the crucial role of SCN4aa in EOD evolution among mormyrid weakly electric fish. The inferred positive selection within Mormyridae makes this gene a prime candidate for further investigation of the divergent evolution of pulse-type EODs among closely related species.

Screening differentially expressed genes in hepatocellular carcinoma by cross-species and cross-carcinogenic factors strategy based on RNA-Seq

AIM: To identify key molecules involved in progression of hepatocellular carcinoma (HCC) by cross-species (human and tupaia) and cross-carcinogenic factors (HBV and AFB1) strategy based on RNA sequencing (RNA-Seq).

METHODS: The transcripts in human and tupaia HCC, tumor adjacent liver tissue (para-HCC) and normal liver tissue were thoroughly analyzed by RNA sequencing. Tupaia HCC was induced by HBV infection or aflatoxin B1 (AFB1). Differentially expressed genes were collected between these tissues, and common differential genes which cross human and Tupaia as well as HBV and AFB1 carcinogenic factors were identified.

RESULTS: Compared to human para-HCC and normal liver tissues, 68 differential genes were screened in human HCC, among which 14 were up-regulated and 54 down-regulated. Compared to Tupaia para-HCC and normal liver tissues, 314 differential genes were screened in HBV induced Tupaia HCC, and 20 were screened in AFB1 induced Tupaia HCC. There were 11 common differential genes between HBV and AFB1 induced Tupaia HCC, all of which were down-regulated. There were 2 common differential genes between human HCC and Tupaia HCC, and they were apolipoprotein F (APOF) and insulin-like growth factor binding protein, acid labile subunit (IGFALS), both of which were down-regulated in HCC.

CONCLUSION: The cross-species, cross-carcinogenic factors screening strategy based on RNA-Seq may promote the process of identifying key molecules for human HCC. APOF and IGFALS may be important factors for HCC.

The myogenic electric organ of Sternopygus macrurus: a non-contractile tissue with a skeletal muscle transcriptome

In most electric fish species, the electric organ (EO) derives from striated muscle cells that suppress many muscle properties. In the gymnotiform Sternopygus macrurus, mature electrocytes, the current-producing cells of the EO, do not contain sarcomeres, yet they continue to make some cytoskeletal and sarcomeric proteins and the muscle transcription factors (MTFs) that induce their expression. In order to more comprehensively examine the transcriptional regulation of genes associated with the formation and maintenance of the contractile sarcomere complex, results from expression analysis using qRT-PCR were informed by deep RNA sequencing of transcriptomes and miRNA compositions of muscle and EO tissues from adult S. macrurus. Our data show that: (1) components associated with the homeostasis of the sarcomere and sarcomere-sarcolemma linkage were transcribed in EO at levels similar to those in muscle; (2) MTF families associated with activation of the skeletal muscle program were not differentially expressed between these tissues; and (3) a set of microRNAs that are implicated in regulation of the muscle phenotype are enriched in EO. These data support the development of a unique and highly specialized non-contractile electrogenic cell that emerges from a striated phenotype and further differentiates with little modification in its transcript composition. This comprehensive analysis of parallel mRNA and miRNA profiles is not only a foundation for functional studies aimed at identifying mechanisms underlying the transcription-independent myogenic program in S. macrurus EO, but also has important implications to many vertebrate cell types that independently activate or suppress specific features of the skeletal muscle program.