Natterin-like depletion by CRISPR/Cas9 impairs zebrafish (Danio rerio) embryonic development

Background The Natterin protein family was first discovered in the venom of the medically significant fish Thalassophryne nattereri, and over the last decade natterin-like genes have been identified in various organisms, notably performing immune-related functions. Previous findings support natterin-like genes as effector defense molecules able to activate multiprotein complexes driving the host innate immune response, notably due to the pore-forming function of the aerolysin superfamily members. Herein, employing a combination of the CRISPR/Cas9 depletion system, phenotype-based screening, and morphometric methods, we evaluated the role of one family member, LOC795232, in the embryonic development of zebrafish since it might be implicated in multiple roles and characterization of the null mutant is central for analysis of gene activity. Results Multiple sequence alignment revealed that the candidate natterin-like has the highest similarity to zebrafish aep1, a putative and better characterized fish-specific defense molecule from the same family. Compared to other species, zebrafish have many natterin-like copies. Whole-mount in situ hybridization confirmed the knockout and mutant embryos exhibited epiboly delay, growth retardation, yolk sac and heart edema, absent or diminished swim bladder, spinal defects, small eyes and head, heart dysfunction, and behavioral impairment. As previously demonstrated, ribonucleoproteins composed of Cas9 and duplex guide RNAs are effective at inducing mutations in the F0 zebrafish. Conclusions The considerably high natterin-like copies in zebrafish compared to other species might be due to the teleost-specific whole genome duplication and followed by subfunctionalization or neofunctionalization. In the present work, we described some of the natterin-like features in the zebrafish development and infer that natterin-like proteins potentially contribute to the embryonary development and immune response. Supplementary Information The online version contains supplementary material available at 10.1186/s12864-022-08369-z.


Background
Over the past 15 years, significant progress has been made in expanding the knowledge on the Natterin family since the first discovery of the four founder members in the venom of the medically significant toadfish Thalassophryne nattereri (VTn) [1]. One of the main symptoms resulting from T. nattereri envenomation is the immediate and intense pain that persists over 24 h. Erythema and edema are also shortly noticed with the efflorescence of bubbles with serous content. These lesions progress to long-remaining necrosis with a delayed healing process [2].
The presence of Natterin family members in evolutionarily divergent non-venomous species suggests an adaptive value consistent with the functions plurality, including immunity, signaling, and development rather than its feature only as a toxin. It is believed that Natterins have a substantial role as effector defense molecules. Their ability to activate cells and cytoplasmic multiprotein complexes that drive the host innate immune response might be due to the fact that all Natterin-like proteins have pore-forming domains and are generally classified as part of the aerolysin superfamily [3][4][5][6][7].
The aep1, the best-known natterin-like protein from zebrafish, was suggested as a novel defense molecule and the first one to have its structures unveiled by crystal and electron microscopy [13], along with additional functions reported in other models; then it has been used as a reference to compare with other zebrafish natterin-like. Considering that, to expand the knowledge over this protein family, we decided to focus on another zebrafish protein that lacks functional description in the literature and has the highest sequence similarity to aep1, i.e., LOC795232.
During the last three decades, zebrafish has emerged as an alternative vertebrate model for molecular and cellular development studies. Compared with other traditional vertebrate models, it offers advantages like a high reproduction rate, easy use and low cost, and high homology of genes known to be associated with human diseases. The embryo transparency and the easy prospect for genetic edition, including the combined system CRISPR (clustered regularly interspaced short palindromic repeats) with Cas9 nuclease [16], makes it possible for a wide exploration of the molecular mechanisms involved in different physiological and developmental processes.
Together with the discovery of the wide expression of natterin-like in different organisms, these advances provide a broader perspective of mechanisms by which pathways featuring natterin-like may protect against infections and provide specificity in physiological processes. Therefore, we hypothesize that in addition to Natterin playing a role as toxins and as effector molecule in immunity, it might have a role in fish development as other examples of critical genes mentioned in the literature that also have multiple roles, notably in immune response and development, such as β-defensins and HOX, respectively [17,18].
To test the hypothesis that members of the Natterin family act in key functions, we evaluate in zebrafish a combination of RNA duplexes of chemically synthesized crRNA and trans-activating crRNA (tracrRNA) molecules complexed with Cas9 protein to form duplex guide ribonucleoproteins (dgRNPs) that act as a depletion system. Phenotype-based screening and morphometric methods were applied to evaluate the specific function of the natterin-like LOC795232, present on zebrafish chromosome 7, in the early embryonic stages of development.
• Natterin-like might be a novel fish-specific defense molecule. • Natterin-like proteins are thought to be pore-forming molecules. • Reverse genetic study and phenotypic characterization suggests natterin-like genes may have roles in zebrafish development.

Multiple alignment analysis of Natterin-like in zebrafish
Many genes were duplicated during teleost evolution after divergence from the tetrapod lineage [19]. In T. nattereri four genes are responsible for encoding all members of the Natterin family (natterin 1 to 4). A bioinformatic survey aiming to know how many natterin genes are in zebrafish confirmed that its genome contains 10 natterin-like genes encoding 11 proteins as described in Table 1. When compared in terms of structure, all of them have an aerolysin-like motif in the C-terminal region (~ 116-360 aa) and a jacalin-like module in the N-terminal region (~ 13-186 aa), except XP_021325134.1, the largest natterin-like protein of zebrafish with 368 amino acids (aa), encoded by LOC101882550, which has the aerolysin-like domain in C-terminal and unknown module in N-terminal portion. The comparison among zebrafish natterin-like proteins supported the selection of the XP_017212453.1 encoded by LOC795232 as the candidate of our study since it presents the highest identity with well-known similar proteins. All zebrafish natterin-like [except XP_021325134.1, percentage of identity (PID) 10.1%] had a PID ranging from 50.6 to 60.6% with the candidate gene (Table 1). When the sequences were compared in the conserved β portion of the aerolysin domain, the similarity to our subject increased and ranged from 54.2 to 66.7%. In contrast, the NP_001013322.1 protein, encoded by another putative zebrafish natterin-like (aep1), showed the higher identity for either the complete sequence or conserved aerolysin domain comparison (Fig. 1, Table 1).
Then we found the conserved residues of the Natterin family, i.e., AGIP (ala-gly-ile-pro), in most zebrafish proteins, except XP_021325134.1. Unlike all four natterins found in T. nattereri, the zebrafish proteins also showed in the N-terminal region the YPT (tyr-pro-thr) conserved residues, a galactose-binding site [20], with a 90% consensus, and other conserved residues with unknown functions, like VLTNVN, FDLPYTG, AND LQYETKG_KGV, within the aerolysin domains, except for XP_021325134.1 (Fig. 1). These differences in the founder proteins from a venomous fish to those in zebrafish may highlight the variability in the function they exert in the latter, which seems to be more related to immune and developmental processes.

Natterin-like gene LOC795232 is expressed during embryogenesis, and CRISPR/Cas9 successfully knocked it out
We followed the embryos from 24 to 144 h post-fertilization (hpf ) to confirm the expression of LOC795232 transcripts by qualitative whole-mount in situ hybridization (WISH) analysis, which is widely used for spatial and temporal detection of the expression throughout development [21]. The WISH suggests broadest expression at 24 hpf, and the transcripts were visualized in both the tailbud region and in the head, notably between the gill and heart regions of WT embryos, probably comprising the heart during the early stage of embryo development (24 to 144 hpf, Fig. 2). After this period, the expression decreased. However, from the hatching longpec stage (48 hpf ) to the early larval period (72 to 144 hpf ), the natterin-like transcripts persisted at low levels restricted to the head region (data not shown). Further, WISH analysis of KO displayed a mild expression in the embryos compared to WT.
To investigate the natterin-like role in zebrafish development, we used CRISPR/Cas9 approach to generate a deletion using the synthesized duplex guide RNA (dgRNA) composed by crRNA and tracrRNA, identified according to Bae et al. [22], in exon 2 of chromosome 7 (genomic identity: ID 795232, XM_017356964, 1298 bp). The dgRNA-Cas RNP complex was microinjected in 0 hpf embryos as 50 ng.μL − 1 of dgRNA with Cas9 enzyme The zebrafish natterin-like protein sequences were obtained from National Center for Biotechnology Information (NCBI) and aligned through Clustal Omega (European Molecular Biology Laboratory -EMBL/European Bioinformatics Institute -EBI). The alignment is displayed in MView evidencing the coverage (cov), percentage of identity (pid), and consensus among proteins. The highlight in yellow→indicates the reference natterin protein studied herein coded by LOC795232 and the Natterin family conserved residues (AGIP) are highlighted in the red square at 250 ng.μL − 1 . In this experimental model, CRISPR/ Cas9-mediated mutagenesis is considerably improved by using a two-RNA component (crRNA:tracrRNA) version of the CRISPR system and it had been shown previously to be effective at inducing mutations in zebrafish F0 embryos [23][24][25]. Considering the robustness of the system and validity in F0 embryos, no stable adult lines were generated in the study.

Natterin-like gene LOC795232 loss promotes epiboly arrest and increased mortality and morphological abnormalities
We tested for the influence of natterin-like expression on development by analyzing epiboly, segmentation, survival, and occurrence of abnormalities. First, we assessed potential epiboly delay when expression was abolished in KO compared to WT, according to Bruce and Heisenberg [26]. A substantial delay in epiboly progression was observed in KO embryos from the oblong-sphere stage and throughout the doming and progression phase (Additional File 1).
Epiboly retard was first evident at 5 hpf. Then, at 11 hpf, we observed that KO embryos (Fig. 3B) did not present complete closure of the blastopore (red asterisks) at the end of the gastrulation period compared to WT embryos ( Fig. 3A). At the segmentation period (19 hpf ), when it was possible to observe anteroposteriorly elongated body with head, somites in the trunk, and tail and rudimentary organs in WT embryos (Fig. 3C), KO embryos ( Fig. 3D) were indistinguishable in appearance, but the number of somites was lower than in WT.
Moreover, we observed that KO embryos had a survival rate of 71.43% in the first 24 hpf compared to 90.12% in WT. The survival in mutants decreased over time, reaching 51.43% in 144 hpf, while the level of survival was 81.24% in WT (Fig. 4A). Next, we observed that KO showed yolk sac edema from 48 hpf (6%), reaching the highest percentage (20%) at 72 hpf (Fig. 4B). The KO abnormality incidence decreased by 96 hpf to 15% but remained high up to 120-144 hpf (12.5 and 11.4%, respectively). Still, delay in yolk absorption, which affects the efficiency of lipid metabolism, leading to a delay in the development of embryos is an example of sub-lethal abnormalities [27]. Regarding pericardial edema, the percentage of individuals with this abnormality progressively Fig. 2 Qualitative evaluation of natterin-like gene expression by whole-mount in situ Hybridization (WISH). Larvae after 24 h post-fertilization (hpf ) were fixed, dehydrated, permeabilized and hybridized with the DIG-labeled probe, complementing the natterin-like RNA sequence (Qiagen, #339500 LCD0168623-BKG). Probe hybridization was determined by binding the AP-labeled anti-DIG antibody (#11093274910, Roche Diagnostics) and revealed through chromogenic reaction by NBT/BCIP solution. The images were obtained in the stereomicroscope STEREO LUMAR-V-12 Carl Zeiss. The picture snippets indicate the area of the larva's body where the blue-purple precipitate (indicating gene product) can be observed in the non-depleted wild-type (WT) larvae (left), indicative of the presence of the natterin-like transcripts; and in the depleted (KO) larvae by CRISPR/Cas9 system (crRNA natterin-like WD07479944, Sigma) (right). The embryos body changed the shape due to long-lasting high temperature incubation in conical-bottom microtubes, not representing developmental malformations increased from 24 up to 120 hpf and stabilized at 144 hpf (7.1, 11.8, 15, 18.8, and 18.6%). As soon as the swim bladder became evident at 96 hpf [28], 100% of KOs did not present this organ. The frequency dropped to 41.3% at 120 hpf and 35.7% at 144 hpf, demonstrating a retard in its inflation. Still, Fig. 4B shows that few individuals (1.2% at 48 hpf, 1.3% at 96 hpf, and 1.3% at 120 hpf ) had pigmentation deficiency (See additional files 2 and 3 for detailed phenotype-based screening data and supplemental images of embryos).  When teratogenic abnormalities were assessed in KOs, we observed increased frequency of individuals with spinal defects from 48 to 144 hpf (5.9, 8.5, 12.5, 16.2, and 17.1%). Moreover, from 24 to 120 hpf, we observed an increase in the developmental delay (2.9, 3.5, 6.8, 10, and 11.3%), which dropped to 4.3% at 144 hpf (Fig. 4C).
General morphological screening in early stages of development comprehended phenotypic measurements of the CRISPR/Cas9 system depletion efficacy. Alongside genotypic analysis, the phenotype evaluation makes it possible to screen the role of one or multiple genes [29], and the assessment can be reliably performed due to the extensive catalog of phenotypic changes in zebrafish described in the literature [30].
While non-injected WT embryos did not show any exacerbated morphological defects (embryo controls noninjected or injected with CRISPR dilution buffer did not differ significantly, as assessed in pilot tests), the surviving mutants presented various unusual phenotypes that persisted until later developmental stages. They include abnormal head and brain development without hemorrhage but with apparent necrosis, small eyes with hypopigmentation, body axis with a spinal defect, yolk sac and pericardial edema, swim bladder loss, and tail defects. Most KOs showed several abnormalities (Fig. 5).

Natterin-like gene LOC795232 loss assessed through morphometric analysis.
To better understand the involvement of the natterin-like in embryonic development, we performed a morphometric analysis in zebrafish mutants and compared to noninjected WT embryos (20 per group; Additional file 3). First, body length was measured as an indication of growth impairment. We observed that in comparison to WT, KOs embryos showed a growth arrest from 96 hpf, resulting in a 13% reduction in the body length at 144 hpf (Fig. 6A). Although we observed that the eyes of the mutants grew over time, they ended up 15% smaller than the control ( Fig. 6B; Additional file 2).
In contrast to the continuous growth of the eyes, WT had a progressive decrease in the yolk sac area directly related to its expected absorption. However, the mutants kept larger yolk sacs over time, with an increase of 106% compared to WT, resulting not just due to the low consumption but also because of the edemas recorded (Fig. 6C). As shown in Fig. 4B, all mutants had not shown a swim bladder at 96 hpf, an important element in the larval development, directly involved in posture, fluctuation, and swimming [31]. However, those who developed it showed a 66.3% decrease in the swim bladder area at 120 hpf and 52.5% at 144 hpf compared to WT ( Fig. 6D; Additional file 2).
The head-trunk angle, a measurement of the overall head and axial skeleton development, increased primarily between the periods of segmentation and hatching (20 to 70 hpf ) as a consequence of the embryo's straightening [32]. In Fig. 6E, we observed that the natterin-like depletion modulated the detachment of the head-trunk, decreasing the angle after 72 hpf compared to WT, with a smaller angle (13%) at the end of the analysis.
In zebrafish, the smaller size of the brain may be the result of reduced neural crest cells and muscles that develop from the paraxial mesoderm [33]. Herein, while the mutants presented similar head sizes to the WT at 24 hpf, at 48 hpf KO heads increased by 14%. However, between 96 and 144 hpf, KO heads were smaller than WT (14, 10, and 8%) ( Fig. 6F; Additional files 2 and 3).

Natterin-like gene LOC795232 regulates cardiac and behavioral functions
The zebrafish linear heart tube begins synchronized contractions around 24 hpf and supports a rudimentary circulatory system with red blood cells (erythrocytes) flowing through the dorsal aorta and cardinal vein, the latter of which only fully connects to the heart later [34]. We analyzed the influence of natterin-like on cardiac function by measuring the pericardium area and heart rate. We observed an increase in the pericardial area by 69% at 144 hpf in KO embryos (Fig. 7A). Again, it might be explained by the high incidence of pericardial edema recorded.
Pronounced cardiac malformation coupled with pericardial edema (Fig. 4B) led us to evaluate whether the depletion of natterin-like modulated the heartbeat rate. The heartbeat count revealed a continuous increase from 24 to 72 hpf in WT embryos, maintained until 144 hpf (Additional file 2). However, gene depletion induced bradycardia at 48 and 72 hpf in contrast to tachycardia from 96 to 144 hpf (Fig. 7B).
Since our analysis showed that natterin-like depletion promoted changes in several structures directly related to the movement, such as swim bladder, spine, head, and eyes, along with decreased size, our last step was to evaluate the locomotor and behavioral activity. WT larvae had a normal swimming pattern after the darklight period exposition, as shown in Fig. 8. However, we observed that KO had slightly less activity recorded during the light period, followed by hyperactivity during the dark period (p = 0.0140; Fig. 8). Even during the acclimatization period, it is noticeable that the KO larvae were rather agitated and distinguishable from the movement pattern demonstrated by WT animals.

Discussion
The first members of the Natterin family were described as toxins present in a venomous fish [35]. Afterwards, natterin-like were demonstrated in non-venomous fish and other divergent species, where they were assumed to participate in the immune response.
By in silico screening on available genome databases we previously identified 331 species presenting 859 natterin or natterin-like genes, distributed throughout all kingdoms of life, including plants, fungi, and sessile marine animals with primitive anatomical organization, and teleost fish [36]. Interestingly, although fish represent the majority of species that contain natterin-like genes (109 species with 598 sequences), only four species are venomous presenting venom apparatus, i.e., Plotosus canius, Plotosus lineatus, Thalassophryne amazonica, and T. nattereri [36,37]. Members of the family may have roles in development since many of them show distinct importance in embryonic progress as well as in multifunctional fish tissues such as skin, gills, and intestines of the arctic charr (Salvelinus alpinus) [38], lumpsucker (Cyclopterus lumpus) [39,40], and ovate pompano (Trachinotus ovatus) [41]. Additionally, Cokus et al. [42] identified among numerous genes that aep1 was highly expressed in the embryonic skin of zebrafish, which was not previously reported to be expressed in the skin during early development.
Employing a combination of the CRISPR/Cas9 depletion system, phenotype-based screening, and morphometric methods, we evaluated the specific function of one family member, the natterin-like protein (XP_017212453.1) encoded by the gene LOC795232 in comparison to WT non-injected controls. This representative protein is highly similar to aep1, which is classified as a sugar-binding natterin-like protein and described for receptor recognition and pore formation. In addition, it is activated early in vertebrate's development, before the maturation of the adaptive immune system [15].
Even the gene product annotations regarding aep1 point to molecular functions strongly associated with binding mechanisms, and the cellular components expressing this gene are related to the cell periphery, membrane component, and pore complex, reinforcing its designated functions as pore-forming protein, just like all the natterin-and aerolysin-like proteins. Membrane binding may trigger drastic conformational changes of the aerolysin domain in an environmental-dependent manner, resulting in the membrane-bound octameric pore formation. Besides, deeper analysis suggests a process with a distinct activation mechanism from the previously characterized in prokaryotic members [13].
In the CRISPR/Cas9 system, the dgRNA design was achieved by bioinformatics algorithms with a 97.56% predicted specificity and 63.6% efficiency, confirming the recent findings of Naert et al. [43] that indicates a probability above 80% of a projected dgRNA induce frameshift mutations and generate at least 64% efficiency of mutant biallelic cells. Further works also demonstrate that mutagenesis by dgRNA-Cas9 RNPs is highly effective at stimulating double strand breaks (DSB)-repair-induced mutations in the zebrafish, even at target sites that appear resistant to the activity of canonical sgRNA-Cas9 RNPs [23], and greatly heritable, reaching up to 86% [44], > 80% [45], or 75-99% [46], depending on the study.
Following delivery of dgRNPs into zebrafish eggs, roughly all cells present in developing embryos harbor bi-allelic indel mutations [23,47]. It appears that dgRNPs rarely have off-target effects with phenotypic consequences in the zebrafish [23]. Preeminently, this complex can induce targeted genetic modifications in zebrafish embryos comparable to those obtained using zinc-finger nucleases (ZFNs) and transcription activator-like effector nucleases (TALENs), although presenting the lowest rate of side effects with negligible off-targets. Mutation rates at potential off-target sites are only 1.1-2.5% [44,47].
The detection of expressed mRNAs through WISH is a robust tool for assessing the spatial distribution of gene transcripts [48]. The zebrafish is especially suitable for this analysis and there are many described gene expression patterns available online at the Zebrafish Information Network (ZFIN, https:// zfin. org/ action/ expre ssion/ search/) [47,49,50].
The highest expression of natterin-like at 24 hpf coincides with the beginning of the differentiation of primitive myeloerythroid lineage within the intermediate cell mass (ICM) [51] and with the emergence of the cardiac contractions and the bloodstream [52]. In zebrafish embryos, primitive hematopoietic cells derived from lateral mesoderm arising in two regions, i.e., the posterior caudal region of the animal between the notochord and endoderm of the trunk, the ICM; and the second in the anterior lateral plate mesoderm (ALPM), a more anterior location under the head. These locations coincide with the higher site-specific expression of the natterinlike in the embryos at 24 hpf. The ICM is further divided into the anterior-trunk domain and the posterior blood islands, which form slightly later in the ventral region of the tail. Erythrocytes develop exclusively in the ICM, while the ALPM is the primary site for early macrophage precursors and other myeloid cell types origin [53,54].
It is noteworthy to point out the persistent expression of natterin-like up to at least 144 hpf in the head region, where the zebrafish's self-renewable hematopoietic stem cells (HSCs) seed the pronephros/primitive kidney, indicating induction of natterin-like during hematopoiesis and validating its contribution to this process. HSPCs embedded in the caudal hematopoietic tissue might serve as a source of embryonic macrophages, neutrophils, and monocytes. The kidney, which corresponds to the mammalian bone marrow, will produce myeloid, erythroid, thromboid, and lymphoid lineages, leaving the thymus to produce mature lymphoid T cells throughout adulthood [51].
Zebrafish epiboly demands the synchronized vegetal pole expansion of the blastoderm and yolk syncytial layer. That is used widely during animal development; thus, knowledge acquired by studying epiboly in zebrafish is likely to be relevant for understanding the process in other systems and contexts [55]. Epiboly delay was first evident in KO embryos at 5 hpf after the maternal-tozygotic transition phase [56]. Then, at the 11 hpf stage, we observed that KO embryos did not present complete blastopore closure at the end of the gastrulation period.
Since somites are generated every 30 min in zebrafish [57], the depletion of natterin-like led to a decrease in the number of somites in KO embryos during the segmentation period. Defects in the spine and shortening of the body length observed in natterin-like deficient larvae may be associated with initial defects in the formation of somites that in turn lead to skeletal and muscle deformities in zebrafish [58].
These results reveal natterin-like gene products of zygotic origin playing roles in zebrafish embryo morphogenesis and suggest that their absence may disrupt developmental processes. Most of the major signaling pathways, e.g., TGF-β, FGF, Wnt, Delta-Notch, and retinoic acid, have been identified to drive vertebrate mesoderm development, and many of their interactions have been elucidated. For example, Nodal, BMP, Wnt, and FGF pathways communicate in complex ways to specify both cell fate and cell movement during gastrulation [59]. Wnt, FGF, and Delta-Notch pathways interact with associated transcription factors to direct segmentation [60]. Finally, BMP, Notch, and Wnt pathways interact with associated transcription factors to regulate blood and vessel formation [61].
In addition, the decreased KO eye size may be attributed to the increased activity of Wnt/β-catenin promoted by the absence of natterin-like [62]. Masterblind (mbl)/axin1 mutant zebrafish that have constitutively hyperactive Wnt signaling lack or have highly reduced forebrain and retina [63]. As described by Ueno et al. [64], Wnt/β-catenin signaling in zebrafish embryos, at different times of development, promotes cardiac differentiation before the gastrulation stage and inhibits it later, which led us to associate the control of this pathway with natterin-like. Furthermore, we might link the absence of the bladders in natterin-like depleted larvae with modulation of the Hedgehog signaling pathway, which is involved in its development [28].
Furthermore, zebrafish have become increasingly used in behavioral neuroscience. Although the fish brains are smaller and simpler, the genetic, neuronal, and physiological mechanisms that drive behavioral responses to a variety of stimuli are similar to those observed in mammals [65,66]. The locomotor activity of the 144 hpf zebrafish depends on the integrity of brain function, nervous system development, and visual pathways [67]. Our data showing locomotor hyperactivity in the KO larvae indicate that the absence of the gene may also have affected this complex behavior circuit and, because of that, influence other neurological relevant functions. As a matter of fact, several behavioral responses are the result of cognitive processes, which depend upon structural, physiological, and biochemical characteristics of the central nervous system [65]. Behavioral changes and visual-motor deficiency, for example, after genome editing, have recently been reported by Zhu et al. [68] and Safarian et al. [69], among others.
In this work, it seems that even though the zebrafish have a considerable number of protein copies from the Natterin family, the knockout of one of them caused phenotypic problems during early development. Such abundance in the natterin-like representatives recorded in zebrafish might be due to the teleost-specific (Ts3R) whole-genome duplication (WGD) event. The teleost fish underwent a third genome duplication round in addition to those that occurred in the ancient vertebrate lineage [70,71]. As a consequence of the genomic enlargement and rearrangements during the Ts3R-WGD, the taxon experienced an increase in morphological complexity and innovation, and wide speciation, becoming the most diverse vertebrate group on Earth [72]. However, there are predominantly two ways duplicates may follow, i.e., subfunctionalization, where following genome duplication, the two gene copies degenerate to perform complementary functions that jointly match that of the single ancestral gene; or one copy maintains the ancestral demand while the other diverge, known as neofunctionalization [73]. The last phenomenon explains why each protein may possess a unique role in the development, and its loss results in failures. Besides, regarding the Natterin family, it might explain the broader diversity of functionalities designated to natterin-like proteins, which could be investigated in a molecular evolution study of these genes in more fishes and vertebrates in the future.
Compared to zebrafish, the other species where the four founder natterin members were described, T. nattereri, has fewer genes encoding natterin proteins. This discrepancy is probably because about 75% of the genes from the Ts3R-WGD event may revert to singletons [74]. On the other hand, the duplication-degenerationcomplementation (DDC) hypothesis [75,76] has been proposed as an explanation for the high retention of duplicate genes, which also suggests that genes with simple tissue-and time-specific regulatory elements would be more likely to revert to singletons than those with complex regulation. Thus, we infer that LOC795232 can be both involved in multifunctional processes and share roles with other natterin-like proteins, such as aep1.
Throughout evolution, mainly immune molecular families have expanded in some species, providing critical functional effects against pathogens. Proteins of the Natterin family may be evidence of this diversification. Even the Agnathans, a superclass of jawless fish and the most primitive members of the vertebrates, present a natterin family member, which implements variable lymphocyte receptors (VLRs)-activated complement cytotoxicity for antigen recognition [10]. Moreover, a natterin-like protein found in zebrafish (aep1) with high homology with the one from the lamprey was identified to be a putative immune defense molecule by forming pore-like structures like aerolysin [13]. Ultimately, Chen et al. [15] suggested that aep1 may be a pro-inflammatory protein and an innate immune molecule that triggers the antimicrobial immune responses preventing zebrafish from bacterial infection. Overall, considering the concepts discussed before, to different extents, the natterin-like genes are potentially primitive and important agents in development and immune response.

Conclusions
The present work provided the first demonstration of the natterin-like role in embryonic development using zebrafish. Depletion generated mutants with abnormal phenotypes that worsened over time and died prematurely. The severe KO phenotypic abnormalities included curved body axis with small bodies, head and eyes deformities, with lack or reduced swim bladders, frequently accompanied by pericardial and yolk sac edema. These abnormalities affected the zebrafish's physiologically relevant functions, leading to severe heart dysfunction and a locomotor hyperactivity pattern suggestive of high levels of stress/anxiety.
Taken together, these results demonstrated that programmable CRISPR/Cas9 systems provide a precise, rapid, and reliable method for altering genes in vivo. It opens the venue to using RNA-guided nucleases for genome editing remarkably in the zebrafish, which is highly validated as a powerful vertebrate model for investigating gene functions and translational research on developmental disorders, thereby demonstrating its potential as a simple, customizable and ready-to-use genome editing tool. The emergence of zebrafish genome editing offers promising new opportunities to understand the genetic basis of natterin-like roles employing targeted mutations. The study of zebrafish through morphogenesis and embryogenic development can provide compelling and broadly applicable insights into the genetic, molecular, and cellular control of the Natterin family members. The knowledge over this family is essential for a better comprehension of these proteins since they are not just overspread around fish but also in organisms from a diversity of taxons and there is increasingly growing evidence of their fundamental contributions to the development and immune response.

Multiple sequence alignment
For the multiple alignments and selection of the gene, we used all known zebrafish natterin-like protein sequences obtained from the National Center for Biotechnology Information (NCBI) aligned through the multiple sequence alignment tool from the software Clustal Omega [European Molecular Biology Laboratory -The European Bioinformatics Institute (EMBL-EBI); https:// www. ebi. ac. uk/ Tools/ msa/ clust alo/]. The alignment was visualized in the viewer MView (https:// www. ebi. ac. uk/ Tools/ msa/ mview/). Clustal aligns sequences using a heuristic method that progressively builds a multiple sequence alignment from a series of pairwise alignments. Essentially, Clustal creates multiple sequence alignments through three main steps: performing a pairwise alignment using the progressive alignment method, creating a guide tree, and using the guide tree to carry out multiple alignments.

Zebrafish husbandry
Adult zebrafish (< 18 months old) from AB strain (International Zebrafish Resource Center, Eugene, OR, US) were kept separated by sex and bred under standard conditions of temperature at 28 °C, pH 7, and light-dark cycle (14/10 h) in individual aquariums in an ALESCO (Campinas, Brazil) rack using system water (60 μg [32].

Zebrafish anesthesia, dechorionation, and euthanasia
Anesthesia was performed by immersing larvae in 2 mL of 0.5x E2 medium containing 0.4% tricaine (ethyl-3-aminobenzoate, #MS-222, Sigma Chemical Co., St. Louis, MO, US) for 2 min at room temperature before analysis. At the end of the experiments, euthanasia was obtained by immersion in 4% tricaine diluted in 0.5x E2 medium. After exposure, larvae were checked in an M205C stereomicroscope (Leica Microsystems) to ensure that their hearts were not beating before being placed in a 10% bleach solution. When suitable, 24 hpf larvae were anesthetized and dechorionated by immersion in pronase (#P5147, Sigma) at 0.02 mg.mL − 1 for 5 min.

Evaluation of embryonic survival and development
Mortality (egg coagulation) was checked daily as well as the occurrence of deformities in the surviving embryos, such as the absence of somites and non-detachment of the tail (lethality parameters); head/eye malformation, abnormal yolk sac absorption, pericardium and yolk sac edema, uninflated or absent swimming bladder, and abnormal pigmentation (sublethality parameters); curved tail or shortened tail, spinal deformity, and delayed growth (teratogenicity parameters) were analyzed according to Shah et al. [77]. The images obtained using an M205C stereomicroscope (Leica Microsystems) at 27x magnification were used for measurements using the ImageJ v.1.8.0_172.

Evaluation of zebrafish epiboly
Whole embryo images were collected from non-dechorionated animals aligned in a 1.2% agarose support and covered by 0.5x E2 medium. Images were acquired (27x magnification) every 1 h until 24 h of development in stereomicroscope Leica M205C (LAS V4.11 software) for analysis of epiboly initiation at 5 hpf when the yolk cell domes and deep cells move radially outwards, forming a cap of cells over the yolk; during the progression phase (11 hpf ) when the blastoderm continues to thin, expanding its surface area to envelop the yolk cell; and at 19 hpf when the blastoderm has covered approximately 50% of the yolk when deep cell epiboly temporarily pauses as cells begin to converge dorsally until gastrulation begins [55].

Phenotype-based screening
The WT and KO embryos were anesthetized, aligned in a glass dish in the lateral position, and photographed under an M205C stereomicroscope (Leica Microsystems). The images obtained were used according to the methods previously reported by Bilder et al. [29]. For measurement of the body length, they were evaluated horizontally from the top of the head to the tip of the tail (μm); the head size evaluated by the antero-caudal measurement of the forebrain, midbrain, and hindbrain (μm); the yolk circumference area (μm 2 ); the eye diameter (μm); the swim bladder circumference area (μm 2 ); the angle of the head (in degree, °) evaluated by the opening of the head in relation to the yolk sac; and the pericardial area (μm 2 ) using ImageJ v.1.8.0_172.

Cardiac function analysis
The anesthetized zebrafish embryos of the WT or KO groups, aligned side by side in groups of 5, were recorded for 15 s per day for 6 days in M205C stereomicroscope (Leica Microsystems; LAS V4.11 software) at 80x magnification for the heart rate evaluation by heartbeats count.

Zebrafish locomotor behavior assessment
The evaluation of locomotor activity was carried out by analyzing the swimming behavior of 144 hpf zebrafish larvae upon the dark-light transition according to the modified method of Macaulay et al. [78]. WT and KO larvae (n = 20) were transferred to 96-well plates, one larva per well in 100 μL of 0.5x E2 medium, and analyzed in a Zebrabox System (Viewpoint Life sciences, Lyon, FR). The larvae were analyzed for a total of 25 min; consisting of 20 min of acclimatization in the dark followed by 5 min of 10 cycles of 25 s in the light to induce visual and neurological stimulation followed by 5 s in the dark. Locomotor activity was quantified and analyzed by ZebraLab ™ software by Viewpoint.

Statistical analysis
All values were expressed as mean ± SEM. Experiments were performed independently two times. Parametric data were evaluated using analysis of variance, followed by the Bonferroni correction for multiple comparisons. Non-parametric data were assessed using the Mann-Whitney test. Differences were considered statistically significant at p < 0.05 using Graph-Pad Prism (Graph Pad Software, v6.02, 2013, La Jolla, CA, US).