Sub-Functionalization of Duplicated Genes in the Evolution of Nine-Spined Stickleback Hatching Enzyme

Genome of the Rio Pearlfish (Nematolebias whitei), a bi-annual killifish model for Eco-Evo-Devo in extreme environments

The Rio Pearlfish, Nematolebias whitei, is a bi-annual killifish species inhabiting seasonal pools in the Rio de Janeiro region of Brazil that dry twice per year. Embryos enter dormant diapause stages in the soil, waiting for the inundation of the habitat which triggers hatching and commencement of a new life cycle. Rio Pearlfish represents a convergent, independent origin of annualism from other emerging killifish model species. While some transcriptomic datasets are available for Rio Pearlfish, thus far, a sequenced genome has been unavailable. Here, we present a high quality, 1.2 Gb chromosome-level genome assembly, genome annotations, and a comparative genomic investigation of the Rio Pearlfish as representative of a vertebrate clade that evolved environmentally cued hatching. We show conservation of 3D genome structure across teleost fish evolution, developmental stages, tissues, and cell types. Our analysis of mobile DNA shows that Rio Pearlfish, like other annual killifishes, possesses an expanded transposable element profile with implications for rapid aging and adaptation to harsh conditions. We use the Rio Pearlfish genome to identify its hatching enzyme gene repertoire and the location of the hatching gland, a key first step in understanding the developmental genetic control of hatching. The Rio Pearlfish genome expands the comparative genomic toolkit available to study convergent origins of seasonal life histories, diapause, and rapid aging phenotypes. We present the first set of genomic resources for this emerging model organism, critical for future functional genetic, and multiomic explorations of “Eco-Evo-Devo” phenotypes of resilience and adaptation to extreme environments.

Effects of hatching enzymes on egg envelope digestion in the male-brooding seahorse

In the present study, we aimed to evaluate the effects of hatching enzymes on the egg envelope digestion during the hatching period in the male brooding seahorse. The complementary DNAs encoding two hatching-enzyme genes, high choriolytic enzyme (HCE) and low choriolytic enzyme (LCE), were cloned and functionally characterized from the lined seahorse (Hippocampus erectus). The genomic-synteny analysis confirmed that teleosts shared LCE gene synteny. In contrast, the genomic location of HCE was found to be conserved with pipefish, but not other teleosts, suggesting that translocation into a novel genomic location occurred. Whole-mount in situ hybridization showed that HCE and LCE mRNAs were expressed in hatching gland cells. To determine the digestion mechanisms of HCE and LCE in hatching, recombinant HCE and LCE were generated and their enzyme activities were examined using fertilized egg envelopes and synthetic peptides. Seahorse HCE and LCE independently digested and softened the egg envelopes of the lined seahorse. Although the egg envelope was digested more following HCE and LCE co-treatment, envelope solubilization was not observed. Indeed, both HCE and LCE showed similar substrate specificities toward four different synthetic peptides designed from the cleavage sites of egg envelope proteins. HCE and LCE proteins from other euteleostean fishes showed different specificities, and the egg envelope was solubilized by the cooperative action of HCE and LCE. These results suggest that the function of LCE was degenerated in the lined seahorse. Our results imply a digestion mechanism for evolutionary adaptation in ovoviviparous fish with male pregnancy.

Genomic and phylogenetic analysis of choriolysins, and biological activity of hatching liquid in the flatfish Senegalese sole

The hatching enzymes or choriolysins are key proteases in fish life cycle controlling the release of larvae to surrounding environment that have been suggested as target for novel biotechnological uses. Due to the large amounts of eggs released by the flatfish Solea senegalensis, during the spawning season, the hatching liquid properties and choriolysin-encoding genes were investigated in this species. A genomic analysis identified four putative genes referred to as SseHCEa, SseHCEb, SseLCE and SseHE. The phylogenetic analysis classified these paralogs into two clades, the clade I containing SseHCE paralogs and the clade II containing two well-supported subclades named as HE and LCE. The two SseHCE paralogs were intron-less and both genes were tandemly arrayed very close in the genome. The synteny and gene rearrangement identified in the flatfish lineage indicated that the duplication of these two paralogs occurred recently and they are under divergent evolution. The genes SseHE and SseLCE were structured in 8 exons and 7 introns and the synteny was conserved in teleosts. Expression studies confirmed that the four genes were expressed in the hatching gland cells and they migrate co-ordinately from the head to around the yolk sac close to the hatch with specific temporal and intensity expression profiles. Although the mRNA levels of the four genes peaked in the hours previous to larval hatching, the SseHCE and SseLCE paralogs kept a longer expression than SseHE after hatching. These expression patterns were consistent even when larvae were incubated at different temperatures that modified hatching times. The analysis of hatching-liquid using SDS-PAGE and zymography analyses of hatching liquid identified a major band of expected choriolysin size. The optimal pH for protease activity was 8.5 and inhibition assays using EDTA demonstrated that most of the activity in the hatching liquid was due to metalloproteases with Ca²⁺ ions acting as the most effective metal to restore the activity. All these data provide new clues about the choriolysin evolution and function in flatfish with impact in the aquaculture and the blue cosmetic industry.

Translocation of promoter-conserved hatching enzyme genes with intron-loss provides a new insight in the role of retrocopy during teleostean evolution

The hatcing enzyme gene (HE) encodes a protease that is indispensable for the hatching process and is conserved during vertebrate evolution. During teleostean evolution, it is known that HE experienced a drastic transfiguration of gene structure, namely, losing all of its introns. However, these facts are contradiction with each other, since intron-less genes typically lose their original promoter because of duplication via mature mRNA, called retrocopy. Here, using a comparative genomic assay, we showed that HEs have changed their genomic location several times, with the evolutionary timings of these translocations being identical to those of intron-loss. We further showed that HEs maintain the promoter sequence upstream of them after translocation. Therefore, teleostean HEs are unique genes which have changed intra- (exon-intron) and extra-genomic structure (genomic loci) several times, although their indispensability for the reproductive process of hatching implies that HE genes are translocated by retrocopy with their promoter sequence.

Neofunctionalization of a duplicate hatching enzyme gene during the evolution of teleost fishes

BACKGROUND

Duplication and subsequent neofunctionalization of the teleostean hatching enzyme gene occurred in the common ancestor of Euteleostei and Otocephala, producing two genes belonging to different phylogenetic clades (clade I and II). In euteleosts, the clade I enzyme inherited the activity of the ancestral enzyme of swelling the egg envelope by cleavage of the N-terminal region of egg envelope proteins. The clade II enzyme gained two specific cleavage sites, N-ZPd and mid-ZPd but lost the ancestral activity. Thus, euteleostean clade II enzymes assumed a new function; solubilization of the egg envelope by the cooperative action with clade I enzyme. However, in Otocephala, the clade II gene was lost during evolution. Consequently, in a late group of Otocephala, only the clade I enzyme is present to swell the egg envelope. We evaluated the egg envelope digestion properties of clade I and II enzymes in Gonorynchiformes, an early diverging group of Otocephala, using milkfish, and compared their digestion with those of other fishes. Finally, we propose a hypothesis of the neofunctionalization process.

RESULTS

The milkfish clade II enzyme cleaved N-ZPd but not mid-ZPd, and did not cause solubilization of the egg envelope. We conclude that neofunctionalization is incomplete in the otocephalan clade II enzymes. Comparison of clade I and clade II enzyme characteristics implies that the specificity of the clade II enzymes gradually changed during evolution after the duplication event, and that a change in substrate was required for the addition of the mid-ZPd site and loss of activity at the N-terminal region.

CONCLUSIONS

We infer the process of neofunctionalization of the clade II enzyme after duplication of the gene. The ancestral clade II gene gained N-ZPd cleavage activity in the common ancestral lineage of the Euteleostei and Otocephala. Subsequently, acquisition of cleavage activity at the mid-ZPd site and loss of cleavage activity in the N-terminal region occurred during the evolution of Euteleostei, but not of Otocephala. The clade II enzyme provides an example of the development of a neofunctional gene for which the substrate, the egg envelope protein, has adapted to a gradual change in the specificity of the corresponding enzyme.

Molecular co-evolution of a protease and its substrate elucidated by analysis of the activity of predicted ancestral hatching enzyme

Background

Hatching enzyme is a protease that digests the egg envelope, enabling hatching of the embryo. We have comprehensively studied the molecular mechanisms of the enzyme action to its substrate egg envelope, and determined the gene/protein structure and phylogenetic relationships. Because the hatching enzyme must have evolved while maintaining its ability to digest the egg envelope, the hatching enzyme-egg envelope protein pair is a good model for studying molecular co-evolution of a protease and its substrate.

Results

Hatching enzymes from medaka (Oryzias latipes) and killifish (Fundulus heteroclitus) showed species-specific egg envelope digestion. We found that by introducing four medaka-type residue amino acid substitutions into recombinant killifish hatching enzyme, the mutant killifish hatching enzyme could digest medaka egg envelope. Further, we studied the participation of the cleavage site of the substrate in the species-specificity of hatching enzyme. A P2-site single amino acid substitution was responsible for the species-specificity. Estimation of the activity of the predicted ancestral enzymes towards various types of cleavage sites along with prediction of the evolutionary timing of substitutions allowed prediction of a possible evolutionary pathway, as follows: ancestral hatching enzyme, which had relatively strict substrate specificity, developed broader specificity as a result of four amino acid substitutions in the active site cleft of the enzyme. Subsequently, a single substitution occurred within the cleavage site of the substrate, and the recent feature of species-specificity was established in the hatching enzyme-egg envelope system.

Conclusions

The present study clearly provides an ideal model for protease-substrate co-evolution. The evolutionary process giving rise to species-specific egg envelope digestion of hatching enzyme was initiated by amino acid substitutions in the enzyme, resulting in altered substrate specificity, which later allowed an amino acid substitution in the substrate.