A chromosome-level genome of electric catfish (Malapterurus electricus) provided new insights into order Siluriformes evolution

Abstract

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.

Keywords

Electric fish / Karyotype analysis / Convergent evolution / Environmental adaptation / Historical effective population size

References

[] Aguilar C, Miller MJ, Loaiza JR, González R, Krahe R, De León LF. Tempo and mode of allopatric divergence in the weakly electric fish Sternopygus dariensis in the Isthmus of Panama. Sci Rep, 2019, 9: 18828.

[] Arias CF, Dikow RB, McMillan WO, De León LF. De novo genome assembly of the electric fish Brachyhypopomus occidentalis (Hypopomidae, Gymnotiformes). Genome Biol Evol, 2021, 13: 1-7.

[] Avdeyev P, Jiang S, Aganezov S, Hu F, Alekseyev MA. Reconstruction of ancestral genomes in presence of gene gain and loss. J Comput Biol, 2016, 23: 150-164.

[] Bairoch A, Apweiler R. The SWISS-PROT protein sequence database and its supplement TrEMBL in 2000. Nucleic Acids Res, 2000, 28: 45-48.

[] Bao W, Kojima KK, Kohany O. Repbase Update, a database of repetitive elements in eukaryotic genomes. Mob DNA, 2015, 6: 11.

[] Bauer R. Untersuchungen zur Entladungstätigkeit und zum Beutefangverhalten des Zitterwelses Malapterurus electricus Gmelin 1789 (Siluroidea, Malapteruridae, Lacep. 1803). Z Vgl Physiol, 1968, 59: 371-402.

[] Benson G. Tandem repeats finder: a program to analyze DNA sequences. Nucleic Acids Res, 1999, 27: 573-580.

[] Bigiani A, Ghiaroni V, Fieni F. Channels as taste receptors in vertebrates. Prog Biophys Mol Biol, 2003, 83: 193-225.

[] Birney E, Clamp M, Durbin R. Genewise and genomewise. Genome Res, 2004, 14: 988-995.

[] Chen Y, Chen Y, Shi C, Huang Z, Zhang Y, Li S, Li Y, Ye J, Yu C, Li Z, Zhang X, Wang J, Yang H, Fang L, Chen Q. SOAPnuke: a MapReduce acceleration-supported software for integrated quality control and preprocessing of high-throughput sequencing data. Gigascience, 2018, 7: 1-6.

[] Clark PU, Dyke AS, Shakun JD, Carlson AE, Clark J, Wohlfarth B, Mitrovica JX, Hostetler SW, McCabe AM. The Last Glacial Maximum. Science, 2009, 325: 710-714.

[] Danecek P, Bonfield JK, Liddle J, Marshall J, Ohan V, Pollard MO, Whitwham A, Keane T, McCarthy SA, Davies RM, Li H. Twelve years of SAMtools and BCFtools. Gigascience, 2021, 10: 1-4.

[] Ding W, Zhang X, Zhao X, Jing W, Cao Z, Li J, Huang Y, You X, Wang M, Shi Q, Bing X. A chromosome-level genome assembly of the mandarin fish (Siniperca chuatsi). Front Genet, 2021, 12: 1-15.

[] Dudchenko O, Batra SS, Omer AD, Nyquist SK, Hoeger M, Durand NC, Shamim MS, Machol I, Lander ES, Aiden AP, Aiden EL. De novo assembly of the Aedes aegypti genome using Hi-C yields chromosome-length scaffolds. Science, 2017, 356: 92-95.

[] Durand NC, Shamim MS, Machol I, Rao SSP, Huntley MH, Lander ES, Aiden EL. Juicer provides a one-click system for analyzing loop-resolution Hi-C experiments. Cell Syst, 2016, 3: 95-98.

[] Ebbesson SOE, Donnel DO. Retinal projections in the electric catfish (Malapterurus electricus). Cell Tissue Res, 1980, 213: 497-503.

[] Edgar RC. MUSCLE: Multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res, 2004, 32: 1792-1797.

[] Enright JM, Toomey MB, Sato SY, Temple SE, Allen JR, Fujiwara R, Kramlinger VM, Nagy LD, Johnson KM, Xiao Y, How MJ, Johnson SL, Roberts NW, Kefalov VJ, Peter Guengerich F, Corbo JC. Cyp27c1 red-shifts the spectral sensitivity of photoreceptors by converting vitamin A1 into A2. Curr Biol, 2015, 25: 3048-3057.

[] Fagbenro OA, Adedire CO, Aiyegbeni ML. Food composition and digestive enzymes in the gut of the African electric catfish, Malapterurus electricus (Gmelin 1789) (Malapteruridae). Trop Zool, 2001, 14: 1-6.

[] Feng P, Zheng J, Rossiter SJ, Wang D, Zhao H. Massive losses of taste receptor genes in toothed and baleen whales. Genome Biol Evol, 2014, 6: 1254-1265.

[] Fricke R, Eschmeyer WN, Fong JD (2021) Species by family/subfamily. http://researcharchive.calacademy.org/research/ichthyology/catalog/SpeciesByFamily.asp. Accessed 29 June 2022

[] Gallant JR, O’Connell LA. Studying convergent evolution to relate genotype to behavioral phenotype. J Exp Biol, 2020, 223: jeb213447.

[] Gallant JR, Traeger LL, Volkening JD, Moffett H, Chen PH, Novina CD, Phillips GN Jr Anand R, Wells GB, Pinch M, Güth R, Unguez GA, Albert JS, Zakon HH, Samanta MP, Sussman MR. Nonhuman genetics. Genomic basis for the convergent evolution of electric organs. Science, 2014, 344: 1522-1525.

[] Goli Bi BEP, Kamelan TM, Berté S, Kien KB, Ndiaye S, Kouamelan EP. The first data on the population parameters of the main fish species in manmade Lake Buyo (River Sassandra, Côte d’Ivoire). Egypt J Aquat Biol Fish, 2019, 23: 585-597.

[] Götz S, García-Gómez JM, Terol J, Williams TD, Nagaraj SH, Nueda MJ, Robles M, Talón M, Dopazo J, Conesa A. High-throughput functional annotation and data mining with the Blast2GO suite. Nucleic Acids Res, 2008, 36: 3420-3435.

[] Han MV, Thomas GWC, Lugo-Martinez J, Hahn MW. Estimating gene gain and loss rates in the presence of error in genome assembly and annotation using CAFE 3. Mol Biol Evol, 2013, 30: 1987-1997.

[] Harris RS (2007) Improved pairwise alignment of genomic DNA. The Pennsylvania State University

[] Hasegawa M, Kishino H, Yano T-a. Dating of the human-ape splitting by a molecular clock of mitochondrial DNA. J Mol Evol, 1985, 22: 160-174.

[] Hoff KJ, Stanke M. Predicting genes in single genomes with AUGUSTUS. Curr Protoc Bioinforma, 2019, 65.

[] Hu Y, Wu Q, Ma S, Ma T, Shan L, Wang X, Nie Y, Ning Z, Yan L, Xiu Y, Wei F. Comparative genomics reveals convergent evolution between the bamboo-eating giant and red pandas. Proc Natl Acad Sci USA, 2017, 114: 1081-1086.

[] Janetzko A, Zimmermann H, Volknandt W. The electromotor system of the electric catfish (Malapterurus electricus): a fine-structural analysis. Cell Tissue Res, 1987, 247: 613-624.

[] Jones DT, Taylor WR, Thornton JM. The rapid generation of mutation data matrices from protein sequences. Comput Appl Biosci, 1992, 8: 275-282.

[] Jones P, Binns D, Chang H-Y, Fraser M, Li W, McAnulla C, McWilliam H, Maslen J, Mitchell A, Nuka G, Pesseat S, Quinn AF, Sangrador-Vegas A, Scheremetjew M, Yong S-Y, Lopez R, Hunter S. InterProScan 5: genome-scale protein function classification. Bioinformatics, 2014, 30: 1236-1240.

[] Kanehisa M, Goto S. KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res, 2000, 28: 27-30.

[] Kasahara M, Naruse K, Sasaki S, Nakatani Y, Qu W, Ahsan B, Yamada T, Nagayasu Y, Doi K, Kasai Y, Jindo T, Kobayashi D, Shimada A, Toyoda A, Kuroki Y, Fujiyama A, Sasaki T, Shimizu A, Asakawa S, Shimizu N, et al. The medaka draft genome and insights into vertebrate genome evolution. Nature, 2007, 447: 714-719.

[] Katoh K, Standley DM. MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Mol Biol Evol, 2013, 30: 772-780.

[] Keilwagen J, Hartung F, Grau J. GeMoMa: homology-based gene prediction utilizing intron position conservation and RNA-seq data. Methods Mol Biol, 2019, 1962: 161-177.

[] Kent WJ. BLAT —The BLAST -Like Alignment Tool. Genome Res, 2002, 12: 656-664.

[] Kent WJ, Baertsch R, Hinrichs A, Miller W, Haussler D. Evolution’s cauldron: duplication, deletion, and rearrangement in the mouse and human genomes. Proc Natl Acad Sci USA, 2003, 100: 11484-11489.

[] Kim D, Pertea G, Trapnell C, Pimentel H, Kelley R, Salzberg SL. TopHat2: accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions. Genome Biol, 2013, 14: R36.

[] Kim J, Farré M, Auvil L, Capitanu B, Larkin DM, Ma J, Lewin HA. Reconstruction and evolutionary history of eutherian chromosomes. Proc Natl Acad Sci USA, 2017, 114: E5379-E5388.

[] Kim OTP, Nguyen PT, Shoguchi E, Hisata K, Vo TTB, Inoue J, Shinzato C, Le BTN, Nishitsuji K, Kanda M, Nguyen VH, Nong HV, Satoh N. A draft genome of the striped catfish, Pangasianodon hypophthalmus, for comparative analysis of genes relevant to development and a resource for aquaculture improvement. BMC Genomics, 2018, 19: 733.

[] Kitts PA, Church DM, Thibaud-Nissen F, Choi J, Hem V, Sapojnikov V, Smith RG, Tatusova T, Xiang C, Zherikov A, DiCuccio M, Murphy TD, Pruitt KD, Kimchi A. Assembly: a resource for assembled genomes at NCBI. Nucleic Acids Res, 2016, 44: D73-80.

[] Krijgsman W, Tesakov A, Yanina T, Lazarev S, Danukalova G, Van Baak CGC, Agustí J, Alçiçek MC, Aliyeva E, Bista D, Bruch A, Büyükmeriç Y, Bukhsianidze M, Flecker R, Frolov P, Hoyle TM, Jorissen EL, Kirscher U, Koriche SA, Kroonenberg SB, et al. Quaternary time scales for the Pontocaspian domain: interbasinal connectivity and faunal evolution. Earth-Sci Rev, 2019, 188: 1-40.

[] Kumazawa T, Brand JG, Teeter JH. Amino acid-activated channels in the catfish taste system. Biophys J, 1998, 75: 2757-2766.

[] Lamanna F, Kirschbaum F, Waurick I, Dieterich C, Tiedemann R. 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: 1-17.

[] Lee AM, Myhre AM, Markussen SS, Engen S, Solberg EJ, Haanes H, Røed K, Herfindal I, Heim M, Saether B-E. Decomposing demographic contributions to the effective population size with moose as a case study. Mol Ecol, 2020, 29: 56-70.

[] LeGrande WH. Chromosomal evolution in North American catfishes (Siluriformes: Ictaluridae) with particular emphasis on the madtoms, Noturus. Copeia, 1981, 1981: 33-52.

[] Letunic I, Bork P. Interactive tree of life (iTOL) v5: an online tool for phylogenetic tree display and annotation. Nucleic Acids Res, 2021, 49: W293-W296.

[] Li H, Durbin R. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics, 2009, 25: 1754-1760.

[] Li H, Durbin R. Inference of human population history from individual whole-genome sequences. Nature, 2011, 475: 493-496.

[] Li L, Stoeckert CJJ, Roos DS. OrthoMCL: identification of ortholog groups for eukaryotic genomes. Genome Res, 2003, 13: 2178-2189.

[] Lin JJ, Wang FY, Li WH, Wang TY. The rises and falls of opsin genes in 59 ray-finned fish genomes and their implications for environmental adaptation. Sci Rep, 2017, 7: 1-13.

[] Liu B, Shi Y, Yuan J, Hu X, Zhang H, Li N, Li Z, Chen Y, Mu D, Fan W (2013) Estimation of genomic characteristics by analyzing k-mer frequency in de novo genome projects. arXiv: https://doi.org/10.48550/arXiv.1308.2012

[] Liu D-W, Lu Y, Yan HY, Zakon HH. South American weakly electric fish (Gymnotiformes) are long-wavelength-sensitive cone monochromats. Brain Behav Evol, 2016, 88: 204-212.

[] Liu Z, Liu S, Yao J, Bao L, Zhang J, Li Y, Jiang C, Sun L, Wang R, Zhang Y, Zhou T, Zeng Q, Fu Q, Gao S, Li N, Koren S, Jiang Y, Zimin A, Xu P, Phillippy AM, et al. The channel catfish genome sequence provides insights into the evolution of scale formation in teleosts. Nat Commun, 2016, 7: 11757.

[] Lundberg JG, Friel JP (2003) Siluriformes. Catfishes. http://tolweb.org/Siluriformes/15065/2003.01.20

[] Luo R, Liu B, Xie Y, Li Z, Huang W, Yuan J, He G, Chen Y, Pan Q, Liu Y, Tang J, Wu G, Zhang H, Shi Y, Liu Y, Yu C, Wang B, Lu Y, Han C, Cheung DW, et al. SOAPdenovo2: an empirically improved memory-efficient short-read de novo assembler. Gigascience, 2012, 1: 18.

[] Manni M, Berkeley MR, Seppey M, Simão FA, Zdobnov EM. BUSCO update: novel and streamlined workflows along with broader and deeper phylogenetic coverage for scoring of eukaryotic, prokaryotic, and viral genomes. Mol Biol Evol, 2021, 38: 4647-4654.

[] Mather N, Traves SM, Ho SYW. A practical introduction to sequentially Markovian coalescent methods for estimating demographic history from genomic data. Ecol Evol, 2020, 10: 579.

[] Miller KG, Kominz MA, Browning JV, Wright JD, Mountain GS, Katz ME, Sugarman PJ, Cramer BS, Christie-Blick N, Pekar SF. The Phanerozoic record of global sea-level change. Science, 2005, 310: 1293-1298.

[] Moreno-Hagelsieb G, Latimer K. Choosing BLAST options for better detection of orthologs as reciprocal best hits. Bioinformatics, 2007, 24: 319-324.

[] Nakatani Y, McLysaght A. Genomes as documents of evolutionary history: a probabilistic macrosynteny model for the reconstruction of ancestral genomes. Bioinformatics, 2017, 33: i369-i378.

[] Naruse K, Tanaka M, Mita K, Shima A, Postlethwait J, Mitani H. A medaka gene map: the trace of ancestral vertebrate proto-chromosomes revealed by comparative gene mapping. Genome Res, 2004, 14: 820-828.

[] Nascimento RA, Shimizu MH, Venancio IM, Chiessi CM, Kuhnert H, Johnstone H, Govin A, Lessa D, Ballalai JM, Piacsek P, Mulitza S, Albuquerque ALS. Warmer western tropical South Atlantic during the Last Interglacial relative to the current interglacial period. Glob Planet Change, 2022, 215.

[] Ozerov MY, Flajshans M, Noreikiene K, Vasemägi A, Gross R, Flajšhans M, Noreikiene K, Vasemägi A, Gross R, Flajshans M, Noreikiene K, Vasemägi A, Gross R. Draft genome assembly of the freshwater apex predator wels catfish (Silurus glanis) using linked-read sequencing. G3- Genes Genom Genet, 2020, 10: 3897-3906.

[] Parey E, Louis A, Montfort J, Guiguen Y, Roest Crollius H, Berthelot C. An atlas of fish genome evolution reveals delayed rediploidization following the teleost whole-genome duplication. Genome Res, 2022, 32: 1685-1697.

[] Parker J, Tsagkogeorga G, Cotton JA, Liu Y, Provero P, Stupka E, Rossiter SJ. Genome-wide signatures of convergent evolution in echolocating mammals. Nature, 2013, 502: 228-231.

[] Ravi V, Venkatesh B. The divergent genomes of teleosts. Annu Rev Anim Biosci, 2018, 6: 47-68.

[] Sanderson MJ. R8s: inferring absolute rates of molecular evolution and divergence times in the absence of a molecular clock. Bioinformatics, 2003, 19: 301-302.

[] Schikorski T, Braun N, Zimmermann H. Cytoarchitectural organization of the electromotor system in the electric catfish (Malapterurus electricus). Cell Tissue Res, 1992, 269: 481-493.

[] Servant N, Varoquaux N, Lajoie BR, Viara E, Chen CJ, Vert JP, Heard E, Dekker J, Barillot E. HiC-Pro: an optimized and flexible pipeline for Hi-C data processing. Genome Biol, 2015, 16: 259.

[] Shi P, Zhang J. Contrasting modes of evolution between vertebrate sweet/umami receptor genes and bitter receptor genes. Mol Biol Evol, 2006, 23: 292-300.

[] Stamatakis A. RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics, 2014, 30: 1312-1313.

[] Stanke M, Keller O, Gunduz I, Hayes A, Waack S, Morgenstern B. AUGUSTUS: ab initio prediction of alternative transcripts. Nucleic Acids Res, 2006, 34: W435-W439.

[] Talavera G, Castresana J. Improvement of phylogenies after removing divergent and ambiguously aligned blocks from protein sequence alignments. Syst Biol, 2007, 56: 564-577.

[] Tarailo-Graovac M, Chen N. Using repeat masker to identify repetitive elements in genomic sequences. Curr Protoc Bioinformatics, 2009, 25: 4.10.1-4.10.14.

[] Terhorst J, Kamm JA, Song YS. Robust and scalable inference of population history from hundreds of unphased whole genomes. Nat Genet, 2017, 49: 303-309.

[] Thomas GWC, Hahn MW. Determining the null model for detecting adaptive convergence from genomic data: a case study using echolocating mammals. Mol Biol Evol, 2015, 32: 1232-1236.

[] Traeger LL, Volkening JD, Moffett H, Gallant JR, Chen PH, Novina CD, Phillips GN, Anand R, Wells GB, Pinch M, Güth R, Unguez GA, Albert JS, Zakon H, Sussman MR, Samanta MP. Unique patterns of transcript and miRNA expression in the South American strong voltage electric eel (Electrophorus electricus). BMC Genomics, 2015, 16: 243.

[] Wang Y, Tang H, Debarry JD, Tan X, Li J, Wang X, Lee TH, Jin H, Marler B, Guo H, Kissinger JC, Paterson AH. MCScanX: a toolkit for detection and evolutionary analysis of gene synteny and collinearity. Nucleic Acids Res, 2012, 40: 1-14.

[] Wang O, Chin R, Cheng X, Yan Wu MK, Mao Q, Tang J, Sun Y, Anderson E, Lam HK, Chen D, Zhou Y, Wang L, Fan F, Zou Y, Xie Y, Zhang RY, Drmanac S, Nguyen D, Xu C, Villarosa C, et al. Efficient and unique cobarcoding of second-generation sequencing reads from long DNA molecules enabling cost-effective and accurate sequencing, haplotyping, and de novo assembly. Genome Res, 2019, 29: 798-808.

[] Waterhouse RM, Seppey M, Simao FA, Manni M, Ioannidis P, Klioutchnikov G, Kriventseva EV, Zdobnov EM. BUSCO applications from quality assessments to gene prediction and phylogenomics. Mol Biol Evol, 2018, 35: 543-548.

[] Wong KHY, Levy-Sakin M, Kwok PY. De novo human genome assemblies reveal spectrum of alternative haplotypes in diverse populations. Nat Commun, 2018, 9: 3040.

[] Woods IG, Wilson C, Friedlander B, Chang P, Reyes DK, Nix R, Kelly PD, Chu F, Postlethwait JH, Talbot WS. The zebrafish gene map defines ancestral vertebrate chromosomes. Genome Res, 2005, 15: 1307-1314.

[] Xu S, He Z, Guo Z, Zhang Z, Wyckoff GJ, Greenberg A, Wu C-II, Shi S, Zhang J. Genome-wide convergence during evolution of mangroves from woody plants. Mol Biol Evol, 2017, 34: 1008-1015.

[] Xu S, Wang J, Guo Z, He Z, Shi S. Genomic convergence in the adaptation to extreme environments. Plant Commun, 2020, 1: 872-879.

[] Yang Z. PAML 4: Phylogenetic analysis by maximum likelihood. Mol Biol Evol, 2007, 24: 1586-1591.

[] Yoo K-S, Lee K, Lee Y-S, Oh W-J, Kim HK. Rho guanine nucleotide exchange factor 4 (Arhgef4) deficiency enhances spatial and object recognition memory. Exp Neurobiol, 2020, 29: 334-343.

[] Zachos J, Pagani H, Sloan L, Thomas E, Billups K. Trends, rhythms, and aberrations in global climate 65 Ma to present. Science, 2001, 292: 686-693.

[] Zakon HH, Lu Y, Zwickl DJ, Hillis DM. Sodium channel genes and the evolution of diversity in communication signals of electric fishes: convergent molecular evolution. Proc Natl Acad Sci USA, 2006, 103: 3675-3680.

[] Zhang X, Sun L, Chen W, Wu S, Li Y, Li X, Zhang B, Yao J, Wang H, Xu A. ARHGEF4-mediates the actin cytoskeleton reorganization of hepatic stellate cells in 3-dimensional collagen matrices. Cell Adh Migr, 2019, 13: 169-181.

[] Zhao H, Xu D, Zhang S, Zhang J. Genomic and genetic evidence for the loss of umami taste in bats. Genome Biol Evol, 2012, 4: 73-79.

[] Zhao H, Yang JR, Xu H, Zhang J. Pseudogenization of the umami taste receptor gene Tas1r1 in the giant panda coincided with its dietary switch to bamboo. Mol Biol Evol, 2010, 27: 2669.

[] Zheng S, Shao F, Tao W, Liu Z, Long J, Wang X, Zhang S, Zhao Q, Carleton KL, Kocher TD, Jin L, Wang Z, Peng Z, Wang D, Zhang Y. Chromosome-level assembly of southern catfish (Silurus meridionalis) provides insights into visual adaptation to nocturnal and benthic lifestyles. Mol Ecol Resour, 2021, 21: 1575-1592.

[] Zhu C, Liu H, Pan Z, Cheng L, Sun Y, Wang H, Chang G, Wu N, Ding H, Zhao H, Zhang L, Yu X. Insights into chromosomal evolution and sex determination of Pseudobagrus ussuriensis (Bagridae, Siluriformes) based on a chromosome-level genome. DNA Res, 2022, 29: 1-15.

[] Zou Z, Zhang J. No genome-wide protein sequence convergence for echolocation. Mol Biol Evol, 2015, 32: 1237-1241.