A Genomic Survey of SCPP Family Genes in Fishes Provides Novel Insights into the Evolution of Fish Scales

Genomic signatures associated with recurrent scale loss in cyprinid fish

Scale morphology represents a fundamental feature of fish and a key evolutionary trait underlying fish diversification. Despite frequent and recurrent scale loss throughout fish diversification, comprehensive genome-wide analyses of the genomic signatures associated with scale loss in divergent fish lineages remain scarce. In the current study, we investigated genome-wide signatures, specifically convergent protein-coding gene loss, amino acid substitutions, and cis-regulatory sequence changes, associated with recurrent scale loss in two divergent Cypriniformes lineages based on large-scale genomic, transcriptomic, and epigenetic data. Results demonstrated convergent changes in many genes related to scale formation in divergent scaleless fish lineages, including loss of P/Q-rich scpp genes (e.g. scpp6 and scpp7), accelerated evolution of non-coding elements adjacent to the fgf and fgfr genes, and convergent amino acid changes in genes (e.g. snap29) under relaxed selection. Collectively, these findings highlight the existence of a shared genetic architecture underlying recurrent scale loss in divergent fish lineages, suggesting that evolutionary outcomes may be genetically repeatable and predictable in the convergence of scale loss in fish.

Genomics comparisons of three chromosome-level mudskipper genome assemblies reveal molecular clues for water-to-land evolution and adaptation

INTRODUCTION

Mudskippers are a large group of amphibious fishes that have developed many morphological and physiological capacities to live on land. Genomics comparisons of chromosome-level genome assemblies of three representative mudskippers, Boleophthalmus pectinirostris (BP), Periophthalmus magnuspinnatus (PM) and P. modestus (PMO), may be able to provide novel insights into the water-to-land evolution and adaptation.

METHODS

Two chromosome-level genome assemblies for BP and PM were respectively sequenced by an integration of PacBio, Nanopore and Hi-C sequencing. A series of standard assembly and annotation pipelines were subsequently performed for both mudskippers. We also re-annotated the PMO genome, downloaded from NCBI, to obtain a redundancy-reduced annotation. Three-way comparative analyses of the three mudskipper genomes in a large scale were carried out to discover detailed genomic differences, such as different gene sizes, and potential chromosomal fission and fusion events. Comparisons of several representative gene families among the three amphibious mudskippers and some other teleosts were also performed to find some molecular clues for terrestrial adaptation.

RESULTS

We obtained two high-quality haplotype genome assemblies with 23 and 25 chromosomes for BP and PM respectively. We also found two specific chromosome fission events in PM. Ancestor chromosome analysis has discovered a common fusion event in mudskipper ancestor. This fusion was then retained in all the three mudskipper species. A loss of some SCPP (secretory calcium-binding phosphoprotein) genes were identified in the three mudskipper genomes, which could lead to reduction of scales for a part-time terrestrial residence. The loss of aanat1a gene, encoding an important enzyme (arylalkylamine N-acetyltransferase 1a, AANAT1a) for dopamine metabolism and melatonin biosynthesis, was confirmed in PM but not in PMO (as previously reported existence in BP), suggesting a better air vision of PM than both PMO and BP. Such a tiny variation within the genus Periophthalmus exemplifies to prove a step-by-step evolution for the mudskippers' water-to-land adaptation.

CONCLUSION

These high-quality mudskipper genome assemblies will become valuable genetic resources for in-depth discovery of genomic evolution for the terrestrial adaptation of amphibious fishes.

Expression of secretory calcium-binding phosphoprotein (scpp) genes in medaka during the formation and replacement of pharyngeal teeth

BACKGROUND

Analyses of tooth families and tooth-forming units in medaka with regard to tooth replacement cycles and the localization of odontogenic stem cell niches in the pharyngeal dentition clearly indicate that continuous tooth replacement is maintained. The secretory calcium-binding phosphoprotein (scpp) gene cluster is involved in the formation of mineralized tissues, such as dental and bone tissues, and the genes encoding multiple SCPPs are conserved in fish, amphibians, reptiles, and mammals. In the present study, we examined the expression patterns of several scpp genes in the pharyngeal teeth of medaka to elucidate their roles during tooth formation and replacement.

METHODS

Himedaka (Japanese medaka, Oryzias latipes) of both sexes (body length: 28 to 33 mm) were used in this study. Real-time quantitative reverse transcription-polymerase chain reaction (PCR) (qPCR) data were evaluated using one-way analysis of variance for multi-group comparisons, and the significance of differences was determined by Tukey's comparison test. The expression of scpp genes was examined using in situ hybridization (ISH) with a digoxigenin-labeled, single-stranded antisense probe.

RESULTS

qPCR results showed that several scpp genes were strongly expressed in pharyngeal tissues. ISH analysis revealed specific expression of scpp1, scpp5, and sparc in tooth germ, and scpp5 was continually expressed in the odontoblasts of teeth attached to pedicles, but not in the osteoblasts of pedicles. In addition, many scpp genes were expressed in inner dental epithelium (ide), but not in odontoblasts, and scpp2 consistently showed epithelial-specific expression in the functional teeth. Taken together, these data indicate that specific expression of scpp2 and scpp5 may play a critical role in pharyngeal tooth formation in medaka.

CONCLUSION

We characterized changes in the expression patterns of scpp genes in medaka during the formation and replacement of pharyngeal teeth.

Whole genome assembly of the armored loricariid catfish Ancistrus triradiatus highlights herbivory signatures

The catfish Ancistrus triradiatus belongs to the species-rich family Loricariidae. Loricariids display remarkable traits such as herbivory, a benthic lifestyle, the absence of scales but the presence of dermal bony plates. They are exported as ornamental fish worldwide, with escaped fishes becoming a threat locally. Although genetic and phylogenetic studies are continuously increasing and developmental genetic investigations are underway, no genome assembly has been formally proposed for Loricariidae yet. We report a high-quality genome assembly of Ancistrus triradiatus using long and short reads, and a newly assembled transcriptome. The genome assembly is composed of 9530 scaffolds, including 85.6% of ray-finned fish BUSCOs, and 26,885 predicted protein-coding genes. The genomic GC content is higher than in other catfishes, reflecting the higher metabolism associated with herbivory. The examination of the SCPP gene family indicates that the genes presumably triggering scale loss when absent, are present in the scaleless A. triradiatus, questioning their explanatory role. The analysis of the opsin gene repertoire revealed that gene losses associated to the nocturnal lifestyle of catfishes were not entirely found in A. triradiatus, as the UV-sensitive opsin 5 is present. Finally, most gene family expansions were related to immunity except the gamma crystallin gene family which controls pupil shape and sub-aquatic vision. Thus, the genome of A. triradiatus reveals that fish herbivory may be related to the photic zone habitat, conditions metabolism, photoreception and visual functions. This genome is the first for the catfish suborder Loricarioidei and will serve as backbone for future genetic, developmental and conservation studies.

Spatiotemporal characteristics of the pharyngeal teeth in interspecific distant hybrids of cyprinid fish: Phylogeny and expression of the initiation marker genes

As an important feeding organ and taxonomical characteristic, the pharyngeal teeth of cyprinid fish have very high morphological diversity and exhibit species-specific numbers and arrangements. Many genes have been verified to regulate the pharyngeal teeth development and act as the initiation marker for teeth. Six initiation marker genes for pharyngeal teeth were used as RNA probes to investigate the expression pattern, and these genes were further used to construct a phylogenetic tree for cyprinid fish including some distant hybrids. The results from in situ hybridization showed that similarities and differences existed in the expression of dlx2b, dlx4b, dlx5a, pitx2, fth1b, and scpp5 in the pharyngeal region of the hybrids (BT) by the crosses of blunt snout bream (BSB, ♀) × topmouth culter (TC, ♂). Particularly, we found a high specificity marker gene scpp5 for the early development of pharyngeal teeth. The Scpp5 expression pattern established a clear graphic representation on the spatiotemporal characteristics of the early morphogenesis of pharyngeal teeth in BT and BSB. Our results suggested that the scpp5 expression in 4V¹, 3V¹, and 5V¹ in BT occurred earlier than that in BSB, while the replacement rate of pharyngeal teeth (4V², 3V², and 5V²) was faster in BSB. Phylogenetic analysis revealed that the six marker genes were highly conserved and could be used as the molecular marker for identifying the parents of the distant hybrids in cyprinid fish. The expression patterns of the scpp5 gene was examined in various tissues, including the brain, gill, heart, liver, muscle, skin, fins, gonad, eye, and kidney, showing that the scpp5 gene was ubiquitously expressed, indicating its important role in cyprinid fish.

Phylogeny of NF-YA trans-activation splicing isoforms in vertebrate evolution

NF-Y is a trimeric pioneer Transcription Factor (TF) whose target sequence -the CCAAT box- is present in ~25% of mammalian promoters. We reconstruct the phylogenetic history of the regulatory NF-YA subunit in vertebrates. We find that in addition to the remarkable conservation of the subunits-interaction and DNA-binding parts, the Transcriptional Activation Domain (TAD) is also conserved (>90% identity among bony vertebrates). We infer the phylogeny of the alternatively spliced exon-3 and partial splicing events of exon-7 -7N and 7C- revealing independent clade-specific losses of these regions. These isoforms shape the TAD. Absence of exon-3 in basal deuterostomes, cartilaginous fishes and hagfish, but not in lampreys, suggests that the "short" isoform is primordial, with emergence of exon-3 in chordates. Exon 7N was present in the vertebrate common ancestor, while 7C is a molecular innovation of teleost fishes. RNA-seq analysis in several species confirms expression of all these isoforms. We identify 3 blocks of amino acids in the TAD shared across deuterostomes, yet structural predictions and sequence analyses suggest an evolutionary drive for maintenance of an Intrinsically Disordered Region -IDR- within the TAD. Overall, these data help reconstruct the logic for alternative splicing of this essential eukaryotic TF.

Parallel evolution of ameloblastic scpp genes in bony and cartilaginous vertebrates

In bony vertebrates, skeletal mineralization relies on the secretory calcium-binding phosphoproteins (Scpp) family whose members are acidic extracellular proteins posttranslationally regulated by the Fam20°C kinase. As scpp genes are absent from the elephant shark genome, they are currently thought to be specific to bony fishes (osteichthyans). Here, we report a scpp gene present in elasmobranchs (sharks and rays) that evolved from local tandem duplication of sparc-L 5′ exons and show that both genes experienced recent gene conversion in sharks. The elasmobranch scpp is remarkably similar to the osteichthyan scpp members as they share syntenic and gene structure features, code for a conserved signal peptide, tyrosine-rich and aspartate/glutamate-rich regions, and harbor putative Fam20°C phosphorylation sites. In addition, the catshark scpp is coexpressed with sparc-L and fam20°C in tooth and scale ameloblasts, similarly to some osteichthyan scpp genes. Despite these strong similarities, molecular clock and phylogenetic data demonstrate that the elasmobranch scpp gene originated independently from the osteichthyan scpp gene family. Our study reveals convergent events at the sparc-L locus in the two sister clades of jawed vertebrates, leading to parallel diversification of the skeletal biomineralization toolkit. The molecular evolution of sparc-L and its coexpression with fam20°C in catshark ameloblasts provides a unifying genetic basis that suggests that all convergent scpp duplicates inherited similar features from their sparc-L precursor. This conclusion supports a single origin for the hypermineralized outer odontode layer as produced by an ancestral developmental process performed by Sparc-L, implying the homology of the enamel and enameloid tissues in all vertebrates.

Convergent losses of SCPP genes and ganoid scales among non-teleost actinopterygians

Various secretory calcium-binding phosphoprotein (SCPP) genes are expressed in the skin and jaw during the formation of bone, teeth, and scales in osteichthyans (bony vertebrates). Among these mineralized skeletal units is the ganoid scale, found in many fossil actinopterygians (ray-finned fish) but confirmed only in Polypteriformes (bichirs, reedfish) and Lepisosteiformes (gars) among extant clades. Here, we examined SCPP genes in the genome of seven non-teleost actinopterygian species that possess or do not possess ganoid scales. As a result, 39-43 SCPP genes were identified in Polypteriformes and Lepisosteiformes, whereas 22-24 SCPP genes were found in Acipenseriformes (sturgeons, paddlefish) and Amiiformes (bowfin). Most of these genes form two clusters in the genome of Polypteriformes, Lepisosteiformes, and Amiiformes, and these two clusters are duplicated in Acipenseriformes. Despite their distant phylogenetic relationship, Polypteriformes and Lepisosteiformes retain many orthologous SCPP genes. These results imply that common ancestors of extant actinopterygians possessed a large repertoire of SCPP genes, and that many SCPP genes were lost independently in Acipenseriformes and Amiiformes. Notably, most SCPP genes originally located in one of the two SCPP gene clusters are retained in Polypteriformes and Lepisosteiformes but were secondarily lost in Acipenseriformes and Amiiformes. In Lepisosteiformes, orthologs of these lost genes show high or detectable expression levels in the skin but not in the jaw. We thus hypothesize that many SCPP genes located in this cluster are involved in the formation of ganoid scales in Polypteriformes and Lepisosteiformes, and that their orthologs and ganoid scales were convergently lost in Acipenseriformes and Amiiformes.

From scales to armor: Scale losses and trunk bony plate gains in ray-finned fishes

Actinopterygians (ray‐finned fishes) are the most diversified group of vertebrates and are characterized by a variety of protective structures covering their integument, the evolution of which has intrigued biologists for decades. Paleontological records showed that the first mineralized vertebrate skeleton was composed of dermal bony plates covering the body, including odontogenic and skeletogenic components. Later in evolution, the exoskeleton of actinopterygian's trunk was composed of scale structures. Although scales are nowadays a widespread integument cover, some contemporary lineages do not have scales but bony plates covering their trunk, whereas other lineages are devoid of any such structures. To understand the evolution of the integument coverage and particularly the transition between different structures, we investigated the pattern of scale loss events along with actinopterygian evolution and addressed the functional relationship between the scaleless phenotype and the ecology of fishes. Furthermore, we examined whether the emergence of trunk bony plates was dependent over the presence or absence of scales. To this aim, we used two recently published actinopterygian phylogenies, one including >11,600 species, and by using stochastic mapping and Bayesian methods, we inferred scale loss events and trunk bony plate acquisitions. Our results reveal that a scaled integument is the most frequent state in actinopterygians, but multiple independent scale loss events occurred along their phylogeny with essentially no scale re‐acquisition. Based on linear mixed models, we found evidence supporting that after a scale loss event, fishes tend to change their ecology and adopt a benthic lifestyle. Furthermore, we show that trunk bony plates appeared independently multiple times along the phylogeny. By using fitted likelihood models for character evolution, we show that trunk bony plate acquisitions were dependent on a previous scale loss event. Overall, our findings support the hypothesis that integument cover is a key evolutionary trait underlying actinopterygian radiation.

Coevolution of enamel, ganoin, enameloid, and their matrix SCPP genes in osteichthyans

We resolve debate over the evolution of vertebrate hypermineralized tissues through analyses of matrix protein-encoding secretory calcium-binding phosphoprotein (SCPP) genes and phylogenetic inference of hypermineralized tissues. Among these genes, AMBN and ENAM are found in both sarcopterygians and actinopterygians, whereas AMEL and SCPP5 are found only in sarcopterygians and actinopterygians, respectively. Actinopterygian AMBN, ENAM, and SCPP5 are expressed during the formation of hypermineralized tissues on scales and teeth: ganoin, acrodin, and collar enamel in gar, and acrodin and collar enameloid in zebrafish. Our phylogenetic analyses indicate the emergence of an ancestral enamel in stem-osteichthyans, whereas ganoin emerged in stem-actinopterygians and true enamel in stem-sarcopterygians. Thus, AMBN and ENAM originated in concert with ancestral enamel, SCPP5 evolved in association with ganoin, and AMEL evolved with true enamel. Shifts in gene expression domain and timing explain the evolution of different hypermineralized tissues. We propose that hypermineralized tissues in osteichthyans coevolved with matrix SCPP genes.

•

Ganoin emerged in actinopterygians; true enamel arose in sarcopterygians

•

Dental enamel, acrodin, and enameloid in actinopterygians are related to ganoin

•

SCPP5 evolved in association with ganoin, whereas AMEL evolved with true enamel

•

Shifts in SCPP gene expression explain the evolution of hypermineralized tissues

Comparative transcriptome analysis of scaled and scaleless skins in Gymnocypris eckloni provides insights into the molecular mechanism of scale degeneration

Background

The scale degeneration is thought to be related to the adaptation to the extreme environment with cold climate and high-altitude in schizothoracine fishes. Gymnocypris eckloni, a schizothoracine fish living in plateau waters with the elevation above 2500 m, is nearly esquamate and only covered with shoulder scales and anal scales, making it a good model species to study the molecular mechanism of scale degeneration.

Results

The transcriptomes of shoulder scaled skins (SSS), anal scaled skins (ASS) and scaleless skins (NSS) were sequenced and analyzed in G. eckloni at the age of 1 year. Histological examination showed that shoulder scale had completed its differentiation and anal scale just initiated the differentiation. A total of 578,046 unigenes were obtained from the transcriptomes, with 407,799 unigenes annotated in public databases. A total of 428 and 142 differentially expressed unigenes (DEUs) were identified between SSS and NSS, and between ASS and NSS, respectively, with 45 DEUs that were overlapped. Annotation analysis indicated that these DEUs were mainly enriched in Gene Ontology (GO) terms and KEGG pathways associated with bone and muscle formation, such as myofibril, contractile fiber, cytoskeletal protein binding, muscle structure development, cardiac muscle contraction, hypertrophic cardiomyopathy (HCM) and calcium signaling pathway.

Conclusions

Our results would provide insights into the molecular mechanisms of scale degeneration in G. eckloni and other congeneric fishes. In addition, the transcriptome data provides candidate genes and markers for future studies.

Using zebrafish to study skeletal genomics

While genome-wide association studies (GWAS) have revolutionized our understanding of the genetic architecture of skeletal diseases, animal models are required to identify causal mechanisms and to translate underlying biology into new therapies. Despite large-scale knockout mouse phenotyping efforts, the skeletal functions of most genes residing at GWAS-identified loci remain unknown, highlighting a need for complementary model systems to accelerate gene discovery. Over the past several decades, zebrafish (Danio rerio) has emerged as a powerful system for modeling the genetics of human diseases. In this review, our goal is to outline evidence supporting the utility of zebrafish for accelerating our understanding of human skeletal genomics, as well as gaps in knowledge that need to be filled for this purpose. We do this by providing a basic foundation of the zebrafish skeletal morphophysiology and phenotypes, and surveying evidence of skeletal gene homology and the use of zebrafish for post-GWAS analysis in other tissues and organs. We also outline challenges in translating zebrafish mutant phenotypes. Finally, we conclude with recommendations of future directions and how to leverage the large body of tools and knowledge of skeletal genetics in zebrafish for the needs of human skeletal genomic exploration. Due to their amenability to rapid genetic approaches, as well as the large number of conserved genetic and phenotypic features, there is a strong rationale supporting the use of zebrafish for human skeletal genomic studies.

Identification and characterization of microRNAs involved in scale biomineralization in the naked carp Gymnocypris przewalskii

The mineralized scale derived from skin plays a protective role for the fish body and also possesses important application values in the biomaterial field. However, little is known about fish scale biomineralization and related molecular regulatory mechanisms. Here, we used a comparative microRNA sequencing approach to identify and characterize differentially expressed microRNAs (DEMs) involved in scale biomineralization in the naked carp Gymnocypris przewalskii. A total of 18, 43, and 66 DEMs were obtained from skin tissues covered with initial, developing, and mature scales (IS, DS, and MS) compared with scale-uncovered skin. The target genes of these DEMs were significantly enriched in a sole biomineralization-related sphingolipid signaling pathway. Seven DEMs (dre-miR-124-3p, dre-miR-133a-2-5p, dre-miR-184, dre-miR-206-3p, novel_33, novel_56 and novel_75) were common in IS, DS, and MS. Dre-miR-124-3p, dre-miR-206-3p, and novel_33 were predicted to be able to target biomineralization-related genes. Stem-loop real-time quantitative PCR further confirmed that the common DEMs had higher expression levels in scale-covered skin tissues than that in the gill, intestine, and brain, except for dre-miR-133a-2-5p. Our results suggest that these identified microRNAs may play a role in scale biomineralization in G. przewalskii, and the obtained microRNAs are expected to be candidates in understanding the molecular mechanism of scale biomineralization in fish species.