Protein phosphatase 2A is essential to maintain meiotic arrest, and to prevent Ca2+ burst at spawning and eventual parthenogenesis in the larvacean Oikopleura dioica.

Unfertilized eggs of most animals are arrested at a certain point in the meiotic cell cycles. Reinitiation of meiosis and the start of embryogenesis are triggered by fertilization. This arrest is essential for preventing parthenogenetic activation and for promoting proper initiation of development by fertilization. In the larvacean Oikopleura dioica, which is a simple model organism for studies of chordate development, the unfertilized egg is arrested at metaphase of meiosis I. We show here that protein phosphatase 2A (PP2A) is essential for maintenance of meiotic arrest after spawning of oocytes. Knockdown (KD) of the maternal PP2A catalytic subunit, which was found in functional screening of maternal factors, caused unfertilized eggs to spontaneously release polar bodies after spawning, and then start pseudo-cleavages without fertilization, namely, parthenogenesis. Parthenogenetic embryos failed to undergo proper mitosis and cytokinesis because of lack of a centrosome, which is to be brought into the egg by a sperm. Activation of the KD oocytes was triggered by possible rise of ambient and intracellular pH upon their release from the gonad into seawater at spawning. Live recording of intracellular calcium level of the KD oocytes indicated that the pH rise caused an aberrant Ca2+ burst, which mimicked the Ca2+ burst that occurs at fertilization. Then, the aberrant Ca2+ burst triggered meiosis resumption through Calcium/calmodulin-dependent protein kinase (CaMK II). Therefore, PP2A is essential for maintenance of meiotic arrest and prevention of parthenogenesis by suppressing the aberrant Ca2+ burst at spawning.

Histone deacetylase activity is required for Botrylloides leachii whole-body regeneration.

The colonial tunicate Botrylloides leachii is exceptional at regenerating from a piece of vascular tunic after loss of all adults from the colony. Previous transcriptome analyses indicate a brief period of healing before regeneration of a new adult (zooid) in as little as 8-10 days. However, there is little understanding of how the resulting changes to gene expression, required to drive regeneration, are initiated and how the overall process is regulated. Rapid changes to transcription often occur in response to chromatin changes, mediated by histone modifications such as histone acetylation. Here, we investigated a group of key epigenetic modifiers, histone deacetylases (HDAC), which are known to play an important role in many biological processes such as development, healing and regeneration. Through our transcriptome data, we identified and quantified the expression levels of HDAC and histone acetyltransferase enzymes during whole-body regeneration (WBR). To determine whether HDAC activity is required for WBR, we inhibited its action using valproic acid and trichostatin A. HDAC inhibition prevented the final morphological changes normally associated with WBR and resulted in aberrant gene expression. Botrylloides leachii genes including Slit2, TGF-ß, Piwi and Fzd4 all showed altered mRNA levels upon HDAC inhibition in comparison with the control samples. Additionally, atypical expression of Bl_Piwi was found in immunocytes upon HDAC inhibition. Together, these results show that HDAC function, specifically HDAC I/IIa class enzymes, are vital for B. leachii to undergo WBR successfully.

Different Conformations of Surface Cellulose Molecules in Native Cellulose Microfibrils Revealed by Layer-by-Layer Peeling.

Layer-by-layer peeling of surface molecules of native cellulose microfibrils was performed using a repeated sequential process of 2,2,6,6-tetramethylpiperidine-1-oxyl radical-mediated oxidation followed by hot alkali extraction. Both highly crystalline algal and tunicate celluloses and low-crystalline cotton and wood celluloses were investigated. Initially, the C6-hydroxy groups of the outermost surface molecules of each algal cellulose microfibril facing the exterior had the gauche-gauche (gg) conformation, whereas those facing the interior had the gauche-trans (gt) conformation. All the other C6-hydroxy groups of the cellulose molecules inside the microfibrils contributing to crystalline cellulose I had the trans-gauche (tg) conformation. After surface peeling, the originally second-layer molecules from the microfibril surface became the outermost surface molecules, and the original tg conformation changed to gg and gt conformations. The plant cellulose microfibrils likely had disordered structures for both the outermost surface and second-layer molecules, as demonstrated using the same layer-by-layer peeling technique.

DNA barcoding of a colonial ascidian, Lissoclinum fragile (Van Name, 1902).

Ascidians (tunicates) are marine benthic organisms possessing various pharmacological activities, including anti-oxidant, anti-tumour, antimicrobial, etc. They also play a key role as model organisms to study various neurobehavioral disorders. Ascidian diversity is reportedly less in India due to lack of taxonomists as well as the limitations in morphology based taxonomy. Molecular taxonomy, comprising the sequencing of cytochrome c oxidase 1 gene (barcode region) otherwise known as DNA barcoding reduces these bottlenecks. Since several species of the family Didemnidae closely resemble in morphology, the present study was aimed to develop DNA barcodes of a colonial ascidian, Lissoclinum fragile belonging to the family Didemnidae. CO1 gene of L. fragile from Thoothukudi, Mandapam, and Vizhinjam waters were sequenced and submitted in GenBank, NCBI through Barcode submission tool. BLAST results showed maximum identity (97-100%) for L. fragile collected from different stations. The pairwise genetic distances within species and genera were calculated using Kimura two parameter (K2P) and the phylogenetic tree was constructed using Neighbour-Joining Tree.

Significant reduction in allergenicity of ovalbumin from chicken egg white following treatment with ascidian viscera N-acetylglucosaminidase.

Ovalbumin (OA) is the most abundant ingredient of chicken egg-white allergenic proteins. In the present study we investigated the possibility of reducing OA allergenicity by treatment with a natural protein exhibiting N-acetylglucosaminidase (NA) activity. Ascidian is cultivated as a food resource in northeast Asia. The ascidian viscera NA (AVNA) with almost no other exoglycosidases or proteolytic enzymes was isolated by applying size-exclusion chromatography to a protein precipitate of ascidian viscera. Intact OA was mixed with AVNA containing 0.2, 1.0, and 5.0 Units of NA. Anion-exchange chromatography was then used to isolate OA from AVNA-treated OA. The electrophoretic patterns and N-glycans of each isolated OA from AVNA-treated OA (iOA) were analyzed, and the terminal N-acetylglucosamines of iOA were selectively cleaved with no other degradation occurring. A competitive indirect enzyme-linked immunosorbent assay using rabbit anti-OA sera was performed to investigate the allergenicity of iOA, which was found to be significantly reduced depending on the increased NA activity compared to that of intact OA. These results indicate that OA allergenicity was reduced using a simple and mild treatment process with AVNA, and suggest that ascidian NA is an efficient natural protein for reducing the allergenicity of OA without requiring the use of harsh physical treatments or chemical conjugation.

Follicle cell trypsin-like protease HrOvochymase: Its cDNA cloning, localization, and involvement in the late stage of oogenesis in the ascidian Halocynthia roretzi.

We previously reported that the sperm trypsin-like protease HrAcrosin and its precursor HrProacrosin participate in fertilization of the ascidian Halocynthia roretzi. The HrProacrosin gene is annotated in the H. roretzi genome database as Harore.CG.MTP2014.S89.g15383; our previously reported sequence of HrProacrosin gene appeared to include four nucleotides inserted near the 3'-end of HrProacrosin, resulting in a frame-shift mutation and a premature termination codon. The gene architecture of HrProacrosin and Harore.CG.MTP2014.S89.g15383 resembles that of Xenopus laevis ovochymase-1/OVCH1 and ovochymase-2/OVCH2, which encode egg extracellular polyproteases. Considering these new observations, we evaluated the cDNA cloning, expression, localization, and function of Harore.CG.MTP2014.S89.g15383, herein designated as HrOvochymase/HrOVCH. We found that HrOVCH cDNA consists of a single open reading frame of 1,575 amino acids, containing a signal peptide, three trypsin-like protease domains, and six CUB domains. HrOVCH was transcribed by the testis and ovary, but the majority of protein exists in ovarian follicle cells surrounding eggs. An anti-HrOVCH antibody inhibited elevation of the vitelline coat at a late stage of oogenesis, during the period when self-sterility is acquired. As trypsin inhibitors are reported to block the acquisition of self-sterility during oogenesis, whereas trypsin induces the acquisition of self-sterility and elevation of the vitelline coat in defolliculated ovarian eggs, we propose that HrOVCH may play a role in the acquisition of self-sterility by late-stage H. roretzi oocytes.

Molecular identification of ascidians from the Palk Bay Region, Southeast coast of India.

BACKGROUND: More than 50 ascidian species distributed in the Palk Bay, Southeast coat of India. Up to a very few molecular work has been performed to determine the evolutionary relationships of ascidians in the world. OBJECTIVES: Present study explored the value of mtDNA data in assessing phylogenetic relationships within the family Ascidiacea, Didemnidae, Styelidae, and molecular identification of ascidians from the Palk Bay, Southeast coast of India. MATERIALS AND METHODS: The phylogeny analysis was executed based on mitochondrial cytochrome oxidase subunit I (COI) sequences of ascidian species. A BLASTN search can be run to determine the similarity of an unknown DNA sequence (query) with the collection of all known DNA sequences in GenBank. RESULT: The BLASTN results showed that Didemnum candidum, Ascidia ahodori and Styela clava match the other listings for these species in GenBank. Mitochondrial COI gene sequences of collected ascidians were submitted to GenBank and obtained the accession numbers. CONCLUSION: This study is the first report of molecular identification of ascidian species of the Palk Bay, Southeast coast of India. This information about Palk Bay ascidian communities provides a baseline of general biodiversity of that ecosystem.

Release from meiotic arrest in ascidian eggs requires the activity of two phosphatases but not CaMKII.

The fertilising sperm triggers a transient Ca(2+) increase that releases eggs from cell cycle arrest in the vast majority of animal eggs. In vertebrate eggs, Erp1, an APC/C(cdc20) inhibitor, links release from metaphase II arrest with the Ca(2+) transient and its degradation is triggered by the Ca(2+)-induced activation of CaMKII. By contrast, many invertebrate groups have mature eggs that arrest at metaphase I, and these species do not possess the CaMKII target Erp1 in their genomes. As a consequence, it is unknown exactly how cell cycle arrest at metaphase I is achieved and how the fertilisation Ca(2+) transient overcomes the arrest in the vast majority of animal species. Using live-cell imaging with a novel cyclin reporter to study cell cycle arrest and its release in urochordate ascidians, the closest living invertebrate group to the vertebrates, we have identified a new signalling pathway for cell cycle resumption in which CaMKII plays no part. Instead, we find that the Ca(2+)-activated phosphatase calcineurin (CN) is required for egg activation. Moreover, we demonstrate that parthenogenetic activation of metaphase I-arrested eggs by MEK inhibition, independent of a Ca(2+) increase, requires the activity of a second egg phosphatase: PP2A. Furthermore, PP2A activity, together with CN, is required for normal egg activation during fertilisation. As ascidians are a sister group of the vertebrates, we discuss these findings in relation to cell cycle arrest and egg activation in chordates.

Senescence-associated superoxide dismutase influences mitochondrial gene expression in budding tunicates.

A recent study has shown that in the budding tunicate Polyandrocarpa misakiensis, the mitochondrial respiratory chain (MRC) dramatically attenuates the gene activity during senescence. In this study, we examined the possible involvement of superoxide dismutase (SOD) in the attenuation of gene expression of cytochrome c oxidase subunit 1 (COX1) in aged zooids. By RT-PCR and in situ hybridization, Cu/Zn-SOD (SOD1) was found to be expressed in most cells and tissues of buds and juvenile zooids but showed a conspicuous decline in senescent adult zooids, except in the gonad tissue in which the cytoplasm of juvenile oocytes was stained heavily. This expression pattern of SOD1 was similar to that of COX1. In contrast to SOD1, Mn-SOD (SOD2) was expressed constitutively in both somatic and germline tissues of buds, juvenile zooids, and senescent adult zooids. Knockdown of SOD1 by RNAi diminished the gene activity of not only SOD1 but also of COX1. The resultant zooids had transient deficiencies in growth and budding, and they recovered from these deficiencies approximately 1 month later. Our results indicate that in P. misakiensis, SOD1 is a senescence-associated nuclear gene and that the experimental decline in SOD1 gene expression accompanies the attenuation of MRC gene activity. Although it is uncertain how SOD1 is downregulated during tunicate senescence, the decreased SOD1 activity could be one of the main causes of MRC gene attenuation during normal senescence.

Metal ion selectivity of the vanadium(V)-reductase Vanabin2.

In a previous study, Vanabin2, a member of a family of V(IV)-binding proteins, or Vanabins, was shown to act as a V(V)-reductase. The current study assesses the ability of Vanabin2 to reduce various transition metal ions in vitro. An NADPH-coupled oxidation assay yielded no evidence of reduction activity with the hexavalent transition metal anions, Mo(VI)O4(2-) and W(VI)O4(2-), or with three divalent cations, Mn(II), Ni(II), and Co(II). Although Cu(II) is readily reduced by glutathione and is gradually oxidized in air, this process was not affected by the presence of Vanabin2. In the experiments conducted thus far, Vanabin2 acts only as a V(V)-reductase. This high selectivity may account for the metal ion selectivity of vanadium accumulation in ascidians.

The Pictet-Spengler mechanism involved in the biosynthesis of tetrahydroisoquinoline antitumor antibiotics: a novel function for a nonribosomal peptide synthetase.

The core scaffold of microbial tetrahydroisoquinoline antitumor antibiotics is biosynthesized by a nonribosomal peptide synthetase (NRPS) with novel functions, which catalyzes a highly unusual seven-step transformation involving multiple reductions of thioester intermediates and two rounds of the Pictet-Spengler reaction. The reaction mechanism of saframycin NRPS SfmC has been firmly established by a series of in vitro experiments using various substrate analogs, SfmC domain-deletion mutants and (2)H-labeled NADH and NADPH. The Pictet-Spengler reaction found in the biosynthesis of saframycin heavily relies on the chain length of the cryptic long acyl chain in the peptide substrates. This chapter describes protocols for biochemical characterization of the saframycin NRPS SfmC. They include (1) bioinformatic analysis of related gene clusters, (2) synthesis of intermediate analogs, and (3) enzymatic reactions for both analytical and preparative scale.

Urochordate ascidians possess a single isoform of Aurora kinase that localizes to the midbody via TPX2 in eggs and cleavage stage embryos.

Aurora kinases are key proteins found throughout the eukaryotes that control mitotic progression. Vertebrate Aurora-A and B kinases are thought to have evolved from a single Aurora-kinase isoform closest to that found in present day urochordates. In urochordate ascidians Aurora binds both TPX2 (a vertebrate AURKA partner) and INCENP (a vertebrate AURKB partner) and localizes to centrosomes and spindle microtubules as well as chromosomes and midbody during both meiosis and mitosis. Ascidian Aurora also displays this localization pattern during mitosis in echinoderms, strengthening the idea that non-vertebrate deuterostomes such as the urochordates and echinoderms possess a single form of Aurora kinase that has properties of vertebrate Aurora-kinase A and B. In the ascidian, TPX2 localizes to the centrosome and the spindle poles also as in vertebrates. However, we were surprised to find that TPX2 also localized strongly to the midbody in ascidian eggs and embryos. We thus examined more closely Aurora localization to the midbody by creating two separate point mutations of ascidian Aurora predicted to perturb binding to TPX2. Both forms of mutated Aurora behaved as predicted: neither localized to spindle poles where TPX2 is enriched. Interestingly, neither form of mutated Aurora localized to the midbody where TPX2 is also enriched, suggesting that ascidian Aurora midbody localization required TPX2 binding in ascidians. Functional analysis revealed that inhibition of Aurora kinase with a pharmacological inhibitor or with a dominant negative kinase dead form of Aurora caused cytokinesis failure and perturbed midbody formation during polar body extrusion. Our data support the view that vertebrate Aurora-A and B kinases evolved from a single non-vertebrate deuterostome ancestor. Moreover, since TPX2 localizes to the midbody in ascidian eggs and cleavage stage embryos it may be worthwhile re-assessing whether Aurora A kinase or TPX2 localize to the midbody in eggs and cleavage stage embryos.

Diversity of nonribosomal peptide synthetase genes in the microbial metagenomes of marine sponges.

Genomic mining revealed one major nonribosomal peptide synthetase (NRPS) phylogenetic cluster in 12 marine sponge species, one ascidian, an actinobacterial isolate and seawater. Phylogenetic analysis predicts its taxonomic affiliation to the actinomycetes and hydroxy-phenyl-glycine as a likely substrate. Additionally, a phylogenetically distinct NRPS gene cluster was discovered in the microbial metagenome of the sponge Aplysina aerophoba, which shows highest similarities to NRPS genes that were previously assigned, by ways of single cell genomics, to a Chloroflexi sponge symbiont. Genomic mining studies such as the one presented here for NRPS genes, contribute to on-going efforts to characterize the genomic potential of sponge-associated microbiota for secondary metabolite biosynthesis.

Looking for putative phenoloxidases of compound ascidians: haemocyanin-like proteins in Polyandrocarpa misakiensis and Botryllus schlosseri.

Phenoloxidases (POs) and haemocyanins constitute a family of copper-containing proteins widely distributed among invertebrates. Both of them are able, under appropriate conditions, to convert polyphenols to quinones and induce cytotoxicity through the production of reactive oxygen species, a fundamental event in many immune responses. In ascidians, PO activity has been described and studied in both solitary and colonial species and the enzyme is involved in inflammatory and cytotoxic reactions against foreign cells or molecules, and in the formation of the cytotoxic foci which characterise the nonfusion reaction of botryllids. Expressed genes for two putative POs (CiPO1 and CiPO2) have been recently identified in C. intestinalis. In the present study, we determined the cDNA sequences of two haemocyanin-like proteins from two colonial ascidians: Botryllus schlosseri from the Mediterranean Sea and Polyandrocarpa misakiensis from Japan. Multiple sequence alignments evidenced the similarity between the above sequences and crustacean proPOs whereas the analysis of the three-dimensional structure reveals high similarity with arthropod haemocyanins which share common precursors with arthropod proPOs. Botryllus HLP grouped in the same cluster with Ciona POs, whereas Polyandrocarpa HLP clustered with arthropod haemocyanins; all of them share the full conservation of the six histidines at the two copper-binding sites as well as of other motifs, also found in arthropod haemocyanin subunits, involved in the regulation of enzyme activity. In situ hybridisation indicated that the genes are transcribed inside morula cells, a characteristic haemocyte type in ascidians where PO activity is located, at the beginning of their differentiation. These results represent a first attempt to identify candidate molecules responsible of the PO activity in compound ascidians.

Increased inter-colony fusion rates are associated with reduced COI haplotype diversity in an invasive colonial ascidian Didemnum vexillum.

Considerable progress in our understanding of the population genetic changes associated with biological invasions has been made over the past decade. Using selectively neutral loci, it has been established that reductions in genetic diversity, reflecting founder effects, have occurred during the establishment of some invasive populations. However, some colonial organisms may actually gain an ecological advantage from reduced genetic diversity because of the associated reduction in inter-colony conflict. Here we report population genetic analyses, along with colony fusion experiments, for a highly invasive colonial ascidian, Didemnum vexillum. Analyses based on mitochondrial cytochrome oxidase I (COI) partial coding sequences revealed two distinct D. vexillum clades. One COI clade appears to be restricted to the probable native region (i.e., north-west Pacific Ocean), while the other clade is present in widely dispersed temperate coastal waters around the world. This clade structure was supported by 18S ribosomal DNA (rDNA) sequence data, which revealed a one base-pair difference between the two clades. Recently established populations of D. vexillum in New Zealand displayed greatly reduced COI genetic diversity when compared with D. vexillum in Japan. In association with this reduction in genetic diversity was a significantly higher inter-colony fusion rate between randomly paired New Zealand D. vexillum colonies (80%, standard deviation ±18%) when compared with colonies found in Japan (27%, standard deviation ±15%). The results of this study add to growing evidence that for colonial organisms reductions in population level genetic diversity may alter colony interaction dynamics and enhance the invasive potential of newly colonizing species.

Pattern of cell proliferation during budding in the colonial ascidian Diplosoma listerianum.

Many invertebrates reproduce asexually by budding, but morphogenesis and the role of cell proliferation in this diverse and nonconserved regeneration-like process are generally poorly understood and particularly little investigated in didemnid ascidians. We here analyzed cell proliferation patterns and telomerase activity during budding in the colonial didemnid ascidian Diplosoma listerianum, with special focus on the thoracic bud where a new brain develops de novo. To help define developmental stages of the thoracic bud, the distribution of acetylated tubulin was also examined. We found extensive cell proliferation in both the thoracic and abdominal buds of D. listerianum as well as higher telomerase activity in bud tissue compared to adult tissues. In the parent adult, proliferation was found in various tissues, but was especially intense in the adult esophagus and epicardial structures that protrude into the proliferating and developing buds, confirming these tissues as the primary source of the cells that form the buds. The neural complex in the thoracic bud forms from a hollow tube that appears to separate into the neural gland and the cerebral ganglion. Whereas most of the bud undergoes proliferation, including the hollow tube and the neural gland, the cerebral ganglion shows little or no proliferation. Pulse-chase labeling experiments indicate that the ganglion, as well as the myocardium, in adult zooids are instead composed of postmitotic cells.

Telomerase deficiency in a colonial ascidian after prolonged asexual propagation.

In organisms that propagate by agametic cloning, the parental body is the reproductive unit and fitness increases with clonal size, so that colonial metazoans, despite lack of experimental data, have been considered potentially immortal. Using asexual propagation rate as a measure of somatic performance, and telomerase activity and relative telomere length as molecular markers of senescence, old (7-12 years) asexual strains of a colonial ascidian, Diplosoma listerianum, were compared with their recent sexually produced progeny. We report for the first time evidence for long-term molecular senescence in asexual lineages of a metazoan, and that only passage between sexual generations provides total rejuvenation permitting indefinite propagation and growth. Thus, this colonial ascidian has not fully escaped ageing. The possibility of somatic replicative senescence also potentially helps to explain why metazoans, with the capacity for asexual propagation through agametic cloning, commonly undergo cycles of sexual reproduction in the wild.

Gene expression profile of Ci-VSP in juveniles and adult blood cells of ascidian.

VSP is a transmembrane protein whose cytoplasmic region shows significant similarity to phosphatase and tensin homolog deleted on chromosome 10 (PTEN). Notably, VSP exhibits a unique ability to transduce electrical signals into phosphoinositide turnover by coupling a transmembrane voltage sensor domain to the PTEN-like phosphoinositide phosphatase domain. Moreover, VSP gene is known to be widely conserved among deuterostome genomes, though the function of VSP in vivo remains largely unknown. In the present study, the expression pattern of ascidian VSP(Ci-VSP) was examined in embryos and juveniles of a marine invertebrate chordate, Ciona intestinalis. RT-PCR showed that Ci-VSP is expressed at the larval stage and that expression persists in juveniles. Whole mount in situ hybridization showed that Ci-VSP is expressed in cells of the stomach, intestine and blood cells of 2- to 3-week-old juveniles. Moreover, double staining blood cells from 2-month-old adults with Ci-VSP and Ci-PTEN probes showed that Ci-VSP-positive cells are a distinct population, separate from cells expressing Ci-PTEN. These findings suggest that in addition to its previously suggested roles in testis or sperm, Ci-VSP plays a key role in voltage-induced signal transduction in cells of the digestive system and blood.

Oikopleura dioica alcohol dehydrogenase class 3 provides new insights into the evolution of retinoic acid synthesis in chordates.

Enzymes that synthesize retinoic acid (RA) constitute the first level of regulation of RA action. In vertebrates, enzymes of the medium-chain alcohol dehydrogenase (MDR-Adh) family catalyze the first step of the RA synthetic pathway by oxidizing retinol. Among MDR-Adh enzymes, Adh3 is the only member present in non-vertebrates, and whether Adh3 is actually involved in RA biosynthesis remains uncertain. Here, we investigate the MDR-Adh family in Oikopleura dioica, a urochordate representing the sister group to vertebrates. Oikopleura is of special interest because it has lost the classical RA role in development, which relaxed evolutionary constraints to preserve the RA-genetic machinery, leading to the loss of RA-system components. The hypothesis that Adh3 plays a role in RA synthesis predicts that the relaxation of selection in Oikopleura should have led to the loss of Adh3, or changes in residues related to retinol oxidation. The hypothesis also predicts changes in the expression pattern of Oikopleura Adh3 compared to other chordates that preserved RA-signaling. Our results, however, revealed the presence of a highly conserved Adh3 gene in Oikopleura, with no significant changes in functional residues. Our results also revealed that the Oikopleura Adh3 expression remains unchanged in comparison to other non-vertebrate chordates, restricted to specific compartments of the digestive system. Because Adh3 has been highly conserved in an animal that has dismantled the RA system, we conclude that Adh3 preservation is not due to a conserved role in RA synthesis. Thereby, if Adh3 plays a role in RA synthesis in vertebrates, it might be a lineage-specific neofunctionalization.