SkateBase, an elasmobranch genome project and collection of molecular resources for chondrichthyan fishes

Gene expression, evolution, and the genetics of electrosensing in the smalltooth sawfish, Pristis pectinata

Sawfishes (Pristidae) are large, highly threatened rays named for their tooth-studded rostrum, which is used for prey sensing and capture. Of all five species, the smalltooth sawfish, Pristis pectinata, has experienced the greatest decline in range, currently found in only ~20% of its historic range. To better understand the genetic underpinnings of these taxonomically and morphologically unique animals, we collected transcriptomic data from several tissue types, mapped them to the recently completed reference genome, and contrasted the patterns observed with comparable data from other elasmobranchs. Evidence of positive selection was detected in 79 genes in P. pectinata, several of which are involved in growth factor/receptor tyrosine kinase signaling and body symmetry and may be related to the unique morphology of sawfishes. Changes in these genes may impact cellular responses to environmental conditions such as temperature, dissolved oxygen, and salinity. Data acquired also allow for examination of the molecular components of P. pectinata electrosensory systems, which are highly developed in sawfishes and have likely been influential in their evolutionary success.

Teeth outside the mouth: The evolution and development of shark denticles

Vertebrate skin appendages are incredibly diverse. This diversity, which includes structures such as scales, feathers, and hair, likely evolved from a shared anatomical placode, suggesting broad conservation of the early development of these organs. Some of the earliest known skin appendages are dentine and enamel-rich tooth-like structures, collectively known as odontodes. These appendages evolved over 450 million years ago. Elasmobranchs (sharks, skates, and rays) have retained these ancient skin appendages in the form of both dermal denticles (scales) and oral teeth. Despite our knowledge of denticle function in adult sharks, our understanding of their development and morphogenesis is less advanced. Even though denticles in sharks appear structurally similar to oral teeth, there has been limited data directly comparing the molecular development of these distinct elements. Here, we chart the development of denticles in the embryonic small-spotted catshark (Scyliorhinus canicula) and characterize the expression of conserved genes known to mediate dental development. We find that shark denticle development shares a vast gene expression signature with developing teeth. However, denticles have restricted regenerative potential, as they lack a sox2+ stem cell niche associated with the maintenance of a dental lamina, an essential requirement for continuous tooth replacement. We compare developing denticles to other skin appendages, including both sensory skin appendages and avian feathers. This reveals that denticles are not only tooth-like in structure, but that they also share an ancient developmental gene set that is likely common to all epidermal appendages.

Little skate genome provides insights into genetic programs essential for limb-based locomotion

The little skate Leucoraja erinacea, a cartilaginous fish, displays pelvic fin driven walking-like behavior using genetic programs and neuronal subtypes similar to those of land vertebrates. However, mechanistic studies on little skate motor circuit development have been limited, due to a lack of high-quality reference genome. Here, we generated an assembly of the little skate genome, with precise gene annotation and structures, which allowed post-genome analysis of spinal motor neurons (MNs) essential for locomotion. Through interspecies comparison of mouse, skate and chicken MN transcriptomes, shared and divergent gene expression profiles were identified. Comparison of accessible chromatin regions between mouse and skate MNs predicted shared transcription factor (TF) motifs with divergent ones, which could be used for achieving differential regulation of MN-expressed genes. A greater number of TF motif predictions were observed in MN-expressed genes in mouse than in little skate. These findings suggest conserved and divergent molecular mechanisms controlling MN development of vertebrates during evolution, which might contribute to intricate gene regulatory networks in the emergence of a more sophisticated motor system in tetrapods.

Combined phylogeny and new classification of catsharks (Chondrichthyes: Elasmobranchii: Carcharhiniformes)

This is the first study to combine morphological and molecular characters to infer the phylogenetic relationships among catsharks. All currently valid genera classified in the family Scyliorhinidae s.l. and representatives of other carcharhinoid families plus one lamnoid and two orectoloboids were included as terminal taxa. A total of 143 morphological characters and 44 NADH2 sequences were analysed. Parsimony analyses under different weighting schemes and strengths were used to generate hypotheses of phylogenetic relationships. The phylogenetic analysis of 78 terminal taxa, using the combined dataset and weighting each column separately (SEP; k = 3) resulted in one most-parsimonious cladogram of 4441 steps with the greatest internal resolution of clades and strongest support. The main changes in nomenclature and classification are the revised definition and scope of Scyliorhinidae, Apristurus and Pentanchus and the revalidation of Atelomycteridae. The monophyly of Pentanchidae is supported, as is that of most catshark genera. Two new subfamilies of the family Pentanchidae are defined: Halaelurinae subfam. nov. and Galeinae subfam. nov. Our analysis emphasizes the relevance of morphological characters in the inference of evolutionary history of carcharhinoids and sheds light on the taxonomic status of some genera in need of further exploration.

Origin and evolutionary landscape of Nr2f transcription factors across Metazoa

BACKGROUND

Nuclear Receptor Subfamily 2 Group F (Nr2f) orphan nuclear hormone transcription factors (TFs) are fundamental regulators of many developmental processes in invertebrates and vertebrates. Despite the importance of these TFs throughout metazoan development, previous work has not clearly outlined their evolutionary history.

RESULTS

We integrated molecular phylogeny with comparisons of intron/exon structure, domain architecture, and syntenic conservation to define critical evolutionary events that distinguish the Nr2f gene family in Metazoa. Our data indicate that a single ancestral eumetazoan Nr2f gene predated six main Bilateria subfamilies, which include single Nr2f homologs, here referred to as Nr2f1/2/5/6, that are present in invertebrate protostomes and deuterostomes, Nr2f1/2 homologs in agnathans, and Nr2f1, Nr2f2, Nr2f5, and Nr2f6 orthologs that are found in gnathostomes. Four cnidarian Nr2f1/2/5/6 and three agnathan Nr2f1/2 members are each due to independent expansions, while the vertebrate Nr2f1/Nr2f2 and Nr2f5/Nr2f6 members each form paralogous groups that arose from the established series of whole-genome duplications (WGDs). Nr2f6 members are the most divergent Nr2f subfamily in gnathostomes. Interestingly, in contrast to the other gnathostome Nr2f subfamilies, Nr2f5 has been independently lost in numerous vertebrate lineages. Furthermore, our analysis shows there are differential expansions and losses of Nr2f genes in teleosts following their additional rounds of WGDs.

CONCLUSION

Overall, our analysis of Nr2f gene evolution helps to reveal the origins and previously unrecognized relationships of this ancient TF family, which may allow for greater insights into the conservation of Nr2f functions that shape Metazoan body plans.

Bacterial community dynamics during embryonic development of the little skate (Leucoraja erinacea)

Background

Microbial transmission from parent to offspring is hypothesized to be widespread in vertebrates. However, evidence for this is limited as many evolutionarily important clades remain unexamined. There is currently no data on the microbiota associated with any Chondrichthyan species during embryonic development, despite the global distribution, ecological importance, and phylogenetic position of this clade. In this study, we take the first steps towards filling this gap by investigating the microbiota associated with embryonic development in the little skate, Leucoraja erinacea, a common North Atlantic species and popular system for chondrichthyan biology.

Methods

To assess the potential for bacterial transmission in an oviparous chondrichthyan, we used 16S rRNA amplicon sequencing to characterize the microbial communities associated with the skin, gill, and egg capsule of the little skate, at six points during ontogeny. Community composition was analyzed using the QIIME2 pipeline and microbial continuity between stages was tracked using FEAST.

Results

We identify site-specific and stage-specific microbiota dominated by the bacterial phyla Proteobacteria and Bacteroidetes. This composition is similar to, but distinct from, that of previously published data on the adult microbiota of other chondrichthyan species. Our data reveal that the skate egg capsule harbors a highly diverse bacterial community–particularly on the internal surface of the capsule–and facilitates intergenerational microbial transfer to the offspring. Embryonic skin and external gill tissues host similar bacterial communities; the skin and gill communities later diverge as the internal gills and skin denticles develop.

Conclusions

Our study is the first exploration of the chondrichthyan microbiota throughout ontogeny and provides the first evidence of vertical transmission in this group.

Supplementary Information

The online version contains supplementary material available at 10.1186/s42523-021-00136-x.

Contingency and chance erase necessity in the experimental evolution of ancestral proteins

The roles of chance, contingency, and necessity in evolution are unresolved because they have never been assessed in a single system or on timescales relevant to historical evolution. We combined ancestral protein reconstruction and a new continuous evolution technology to mutate and select proteins in the B-cell lymphoma-2 (BCL-2) family to acquire protein–protein interaction specificities that occurred during animal evolution. By replicating evolutionary trajectories from multiple ancestral proteins, we found that contingency generated over long historical timescales steadily erased necessity and overwhelmed chance as the primary cause of acquired sequence variation; trajectories launched from phylogenetically distant proteins yielded virtually no common mutations, even under strong and identical selection pressures. Chance arose because many sets of mutations could alter specificity at any timepoint; contingency arose because historical substitutions changed these sets. Our results suggest that patterns of variation in BCL-2 sequences – and likely other proteins, too – are idiosyncratic products of a particular and unpredictable course of historical events.

One of the most fundamental and unresolved questions in evolutionary biology is whether the outcomes of evolution are predictable. Is the diversity of life we see today the expected result of organisms adapting to their environment throughout history (also known as natural selection) or the product of random chance? Or did chance events early in history shape the paths that evolution could take next, determining the biological forms that emerged under natural selection much later?

These questions are hard to study because evolution happened only once, long ago. To overcome this barrier, Xie, Pu, Metzger et al. developed an experimental approach that can evolve reconstructed ancestral proteins that existed deep in the past. Using this method, it is possible to replay evolution multiple times, from various historical starting points, under conditions similar to those that existed long ago. The end products of the evolutionary trajectories can then be compared to determine how predictable evolution actually is.

Xie, Pu, Metzger et al. studied proteins belonging to the BCL-2 family, which originated some 800 million years ago. These proteins have diversified greatly over time in both their genetic sequences and their ability to bind to specific partner proteins called co-regulators. Xie, Pu, Metzger et al. synthesized BCL-2 proteins that existed at various times in the past. Each ancestral protein was then allowed to evolve repeatedly under natural selection to acquire the same co-regulator binding functions that evolved during history.

At the end of each evolutionary trajectory, the genetic sequence of the resulting BCL-2 proteins was recorded. This revealed that the outcomes of evolution were almost completely unpredictable: trajectories initiated from the same ancestral protein produced proteins with very different sequences, and proteins launched from different ancestral starting points were even more dissimilar.

Further experiments identified the mutations in each trajectory that caused changes in coregulator binding. When these mutations were introduced into other ancestral proteins, they did not yield the same change in function. This suggests that early chance events influenced each protein’s evolution in an unpredictable way by opening and closing the paths available to it in the future.

This research expands our understanding of evolution on a molecular level whilst providing a new experimental approach for studying evolutionary drivers in more detail. The results suggest that BCL-2 proteins, in all their various forms, are unique products of a particular, unpredictable course of history set in motion by ancient chance events.

The dystrophia myotonica WD repeat-containing protein DMWD and WDR20 differentially regulate USP12 deubiquitinase

Despite its potential clinical relevance, the product of the DMWD (dystrophia myotonica, WD repeat containing) gene is a largely uncharacterized protein. The DMWD amino acid sequence is similar to that of WDR20, a known regulator of the USP12 and USP46 deubiquitinases (DUBs). Here, we apply a combination of in silico and experimental methods to investigate several aspects of DMWD biology. Molecular evolution and phylogenetic analyses reveal that WDR20 and DMWD, similar to USP12 and USP46, arose by duplication of a common ancestor during the whole genome duplication event in the vertebrate ancestor lineage. The analysis of public human gene expression datasets suggests that DMWD expression is positively correlated with USP12 expression in normal tissues and negatively correlated with WDR20 expression in tumors. Strikingly, a survey of the annotated interactome for DMWD and WDR20 reveals a largely nonoverlapping set of interactors for these proteins. Experimentally, we first confirmed that DMWD binds both USP12 and USP46 through direct coimmunoprecipitation of epitope-tagged proteins. We found that DMWD and WDR20 share the same binding interface in USP12, suggesting that their interaction with the DUB may be mutually exclusive. Finally, we show that both DMWD and WDR20 promote USP12 enzymatic activity, but they differentially modulate the subcellular localization of the DUB. Altogether, our findings suggest a model whereby mutually exclusive binding of DMWD and WDR20 to USP12 may lead to formation of deubiquitinase complexes with distinct subcellular localization, potentially targeting different substrate repertoires.

New Species Can Broaden Myelin Research: Suitability of Little Skate, Leucoraja erinacea

Although myelinated nervous systems are shared among 60,000 jawed vertebrates, studies aimed at understanding myelination have focused more and more on mice and zebrafish. To obtain a broader understanding of the myelination process, we examined the little skate, Leucoraja erinacea. The reasons behind initiating studies at this time include: the desire to study a species belonging to an out group of other jawed vertebrates; using a species with embryos accessible throughout development; the availability of genome sequences; and the likelihood that mammalian antibodies recognize homologs in the chosen species. We report that the morphological features of myelination in a skate hatchling, a stage that supports complex behavioral repertoires needed for survival, are highly similar in terms of: appearances of myelinating oligodendrocytes (CNS) and Schwann cells (PNS); the way their levels of myelination conform to axon caliber; and their identity in terms of nodal and paranodal specializations. These features provide a core for further studies to determine: axon-myelinating cell communication; the structures of the proteins and lipids upon which myelinated fibers are formed; the pathways used to transport these molecules to sites of myelin assembly and maintenance; and the gene regulatory networks that control their expressions.

Connexins during 500 Million Years-From Cyclostomes to Mammals

It was previously shown that the connexin gene family had relatively similar subfamily structures in several vertebrate groups. Still, many details were left unclear. There are essentially no data between tunicates, which have connexins that cannot be divided into the classic subfamilies, and teleosts, where the subfamilies are easily recognized. There are also relatively few data for the groups that diverged between the teleosts and mammals. As many of the previously analyzed genomes have been improved, and many more genomes are available, we reanalyzed the connexin gene family and included species from all major vertebrate groups. The major results can be summarized as follows: (i) The same connexin subfamily structures are found in all Gnathostomata (jawed vertebrates), with some variations due to genome duplications, gene duplications and gene losses. (ii) In contrast to previous findings, birds do not have a lower number of connexins than other tetrapods. (iii) The cyclostomes (lampreys and hagfishes) possess genes in the alpha, beta, gamma and delta subfamilies, but only some of the genes show a phylogenetic affinity to specific genes in jawed vertebrates. Thus, two major evolutionary transformations have occurred in this gene family, from tunicates to cyclostomes and from cyclostomes to jawed vertebrates.

Development and regeneration of the crushing dentition in skates (Rajidae)

Sharks and rays (elasmobranchs) have the remarkable capacity to continuously regenerate their teeth. The polyphyodont system is considered the ancestral condition of the gnathostome dentition. Despite this shared regenerative ability, sharks and rays exhibit dramatic interspecific variation in their tooth morphology. Ray (batoidea) teeth typically constitute crushing pads of flattened teeth, whereas shark teeth are pointed, multi-cuspid units. Although recent research has addressed the molecular development of the shark dentition, little is known about that of the ray. Furthermore, how dental diversity within the elasmobranch lineage is achieved remains unknown. Here, we examine dental development and regeneration in two Batoid species: the thornback skate (Raja clavata) and the little skate (Leucoraja erinacea). Using in situ hybridization and immunohistochemistry, we examine the expression of a core gnathostome dental gene set during early development of the skate dentition and compare it to development in the shark. Elasmobranch tooth development is highly conserved, with sox2 likely playing an important role in the initiation and regeneration of teeth. Alterations to conserved genes expressed in an enamel knot-like signalling centre may explain the morphological diversity of elasmobranch teeth, thereby enabling sharks and rays to occupy diverse dietary and ecological niches.

The importance of a single amino acid substitution in reduced red blood cell carbonic anhydrase function of early-diverging fish

In most vertebrates, red blood cell carbonic anhydrase (RBC CA) plays a critical role in carbon dioxide (CO₂) transport and excretion across epithelial tissues. Many early-diverging fishes (e.g., hagfish and chondrichthyans) are unique in possessing plasma-accessible membrane-bound CA-IV in the gills, allowing some CO₂ excretion to occur without involvement from the RBCs. However, implications of this on RBC CA function are unclear. Through homology cloning techniques, we identified the putative protein sequences for RBC CA from nine early-diverging species. In all cases, these sequences contained a modification of the proton shuttle residue His-64, and activity measurements from three early-diverging fish demonstrated significantly reduced CA activity. Site-directed mutagenesis was used to restore the His-64 proton shuttle, which significantly increased RBC CA activity, clearly illustrating the functional significance of His-64 in fish red blood cell CA activity. Bayesian analyses of 55 vertebrate cytoplasmic CA isozymes suggested that independent evolutionary events led to the modification of His-64 and thus reduced CA activity in hagfish and chondrichthyans. Additionally, in early-diverging fish that possess branchial CA-IV, there is an absence of His-64 in RBC CAs and the absence of the Root effect [where a reduction in pH reduces hemoglobin's capacity to bind with oxygen (O₂)]. Taken together, these data indicate that low-activity RBC CA may be present in all fish with branchial CA-IV, and that the high-activity RBC CA seen in most teleosts may have evolved in conjunction with enhanced hemoglobin pH sensitivity.

Estimating effective population size using RADseq: Effects of SNP selection and sample size

Effective population size (Ne ) is a key parameter of population genetics. However, N e remains challenging to estimate for natural populations as several factors are likely to bias estimates. These factors include sampling design, sequencing method, and data filtering. One issue inherent to the restriction site-associated DNA sequencing (RADseq) protocol is missing data and SNP selection criteria (e.g., minimum minor allele frequency, number of SNPs). To evaluate the potential impact of SNP selection criteria on Ne estimates (Linkage Disequilibrium method) we used RADseq data for a nonmodel species, the thornback ray. In this data set, the inbreeding coefficient F IS was positively correlated with the amount of missing data, implying data were missing nonrandomly. The precision of Ne estimates decreased with the number of SNPs. Mean Ne estimates (averaged across 50 random data sets with2000 SNPs) ranged between 237 and 1784. Increasing the percentage of missing data from 25% to 50% increased Ne estimates between 82% and 120%, while increasing the minor allele frequency (MAF) threshold from 0.01 to 0.1 decreased estimates between 71% and 75%. Considering these effects is important when interpreting RADseq data-derived estimates of effective population size in empirical studies.

Skeletal Mineralization in Association with Type X Collagen Expression Is an Ancestral Feature for Jawed Vertebrates

In order to characterize the molecular bases of mineralizing cell evolution, we targeted type X collagen, a nonfibrillar network forming collagen encoded by the Col10a1 gene. It is involved in the process of endochondral ossification in ray-finned fishes and tetrapods (Osteichthyes), but until now unknown in cartilaginous fishes (Chondrichthyes). We show that holocephalans and elasmobranchs have respectively five and six tandemly duplicated Col10a1 gene copies that display conserved genomic synteny with osteichthyan Col10a1 genes. All Col10a1 genes in the catshark Scyliorhinus canicula are expressed in ameloblasts and/or odontoblasts of teeth and scales, during the stages of extracellular matrix protein secretion and mineralization. Only one duplicate is expressed in the endoskeletal (vertebral) mineralizing tissues. We also show that the expression of type X collagen is present in teeth of two osteichthyans, the zebrafish Danio rerio and the western clawed frog Xenopus tropicalis, indicating an ancestral jawed vertebrate involvement of type X collagen in odontode formation. Our findings push the origin of Col10a1 gene prior to the divergence of osteichthyans and chondrichthyans, and demonstrate its ancestral association with mineralization of both the odontode skeleton and the endoskeleton.

Further studies of ion channels in the electroreceptor of the skate through deep sequencing, cloning and cross species comparisons

Our comparative studies seek to understand the structure and function of ion channels in cartilaginous fish that can detect very low voltage gradients in seawater. The principal channels of the electroreceptor include a calcium activated K channel whose α subunit is Kcnma1, and a voltage-dependent calcium channel, Cacna1d. It has also been suggested based on physiological and pharmacological evidence that a voltage-gated K channel is present in the basal membranes of the receptor cells which modulates synaptic transmitter release. Large conductance calcium-activated K channels (BK) are comprised of four α subunits, encoded by Kcnma1 and modulatory β subunits of the Kcnmb class. We recently cloned and published the skate Kcnma1 gene and most of Kcnmb4 using purified mRNA of homogenized electroreceptors. Bellono et al. have recently performed RNA sequencing (RNA-seq) on purified mRNA from skate electroreceptors and found several ion channels including Kcnma1. We searched the Bellono et al. RNA-seq repository for additional channels and subunits. Our most significant findings are the presence of two Shaker type voltage dependent K channel sequences which are grouped together as isoforms in the data repository. The larger of these is a skate ortholog of the voltage dependent fast potassium channel Kv1.1, which is expressed at appreciable levels. The second ortholog is similar to Kv1.5 but has fewer N-terminal amino acids than other species. The sequence for Kv1.5 in the skate is very strongly aligned with the recently reported sequence for potassium channels in the electroreceptors of the cat shark, S. retifer, which also modulate synaptic transmission. The latter channel was designated as Kv1.3 in the initial report, but we suggest that these channels are actually orthologs of each other, and that Kv1.5 is the prevailing designation. We also found a beta subunit sequence (Kcnab2) which may co-assemble with one or both of the voltage gated channels. The new channels and subunits were verified by RT-PCR and the Kv1.1 sequence was confirmed by cloning. We also searched the RNA-seq repository for accessory subunits of Kcnma1, and found a computer-generated assembly that contained a complete sequence of its β subunit, Kcnmb2. Skate Kcnmb2 has a total of 279 amino acids, with 51 novel amino acids at the N-terminus which may play a specific physiological role. This sequence was confirmed by PCR and cloning. However, skate Kcnmb2 is expressed at low levels in the electroreceptor compared to Kcnma1 and skate Kcnmb1 is absent. The evolutionary origin of the newly described K channels and their subunits was studied by alignments with mammalian sequences, including human, and also those in related fish: the whale shark (R. typus), the ghost shark (C.milii), and (S. retifer). There are also orthologous K channels of the lamprey, which has electroreceptors. Tree building and bootstrap programs were used to confirm phylogenetic inferences. Further research should focus on the subcellular locations of these channels, their gating behavior, and the effects of accessory subunits on gating.

Evolutionary Plasticity in Detoxification Gene Modules: The Preservation and Loss of the Pregnane X Receptor in Chondrichthyes Lineages

To appraise how evolutionary processes, such as gene duplication and loss, influence an organism's xenobiotic sensitivity is a critical question in toxicology. Of particular importance are gene families involved in the mediation of detoxification responses, such as members of the nuclear receptor subfamily 1 group I (NR1I), the pregnane X receptor (PXR), and the constitutive androstane receptor (CAR). While documented in multiple vertebrate genomes, PXR and CAR display an intriguing gene distribution. PXR is absent in birds and reptiles, while CAR shows a tetrapod-specific occurrence. More elusive is the presence of PXR and CAR gene orthologs in early branching and ecologically-important Chondrichthyes (chimaeras, sharks and rays). Therefore, we investigated various genome projects and use them to provide the first identification and functional characterization of a Chondrichthyan PXR from the chimaera elephant shark (Callorhinchus milii, Holocephali). Additionally, we substantiate the targeted PXR gene loss in Elasmobranchii (sharks and rays). Compared to other vertebrate groups, the chimaera PXR ortholog displays a diverse expression pattern (skin and gills) and a unique activation profile by classical xenobiotic ligands. Our findings provide insights into the molecular landscape of detoxification mechanisms and suggest lineage-specific adaptations in response to xenobiotics in gnathostome evolution.

Ezh1 arises from Ezh2 gene duplication but its function is not required for zebrafish development

Trimethylation on H3K27 mediated by Polycomb Repressive Complex 2 (PRC2) is required to control gene repression programs involved in development, regulation of tissue homeostasis or maintenance and lineage specification of stem cells. In Drosophila, the PRC2 catalytic subunit is the single protein E(z), while in mammals this function is fulfilled by two proteins, Ezh1 and Ezh2. Based on database searches, we propose that Ezh1 arose from an Ezh2 gene duplication that has occurred in the common ancestor to elasmobranchs and bony vertebrates. Expression studies in zebrafish using in situ hybridization and RT-PCR followed by the sequencing of the amplicon revealed that ezh1 mRNAs are maternally deposited. Then, ezh1 transcripts are ubiquitously distributed in the entire embryo at 24 hpf and become more restricted to anterior part of the embryo at later developmental stages. To unveil the function of ezh1 in zebrafish, a mutant line was generated using the TALEN technology. Ezh1-deficient mutant fish are viable and fertile, but the loss of ezh1 function is responsible for the earlier death of ezh2 mutant larvae indicating that ezh1 contributes to zebrafish development in absence of zygotic ezh2 gene function. Furthermore, we show that presence of ezh1 transcripts from the maternal origin accounts for the delayed lethality of ezh2-deficient larvae.

Chondrichthyan research in South America: Endocrinology overview and research trends over 50 years (1967-2016) compared to the rest of the world

The endocrine system plays a crucial role in regulating the activity of cells and organs among vertebrates, including the class Chondrichthyes. Accordingly, Chondrichthyan endocrinology publications have been steadily increasing in the global literature. However, while interest in South American Chondrichthyan research has been growing over the last 50 years, the field of endocrinology related to Chondrichthyans has been limited. Understanding the trajectory of a scientific discipline assists researchers and stakeholders in making decisions regarding which research areas require further attention. Further, visualisation techniques based on bibliometric analysis of scientific publications assist in understanding fluctuations in the trends of specific research fields over time. In this study, Chondrichthyan research publications over time were analysed by creating visualisation maps using VOSviewer bibliometric software. Trends in South America Chondrichthyan research with an emphasis on endocrinology were explored over a 50-year period (1967-2016). These trends were compared with Chondrichthyans research worldwide for the more recent 15-year period (2002-2016). The number of South America Chondrichthyan scientific publications increased from six during the 1967-1981 period to 112 in 2016. However, only eight papers were found published in the area of Chondrichthyan endocrinology research. Fisheries, reproduction and taxonomy were the dominate research areas in South America over the 50 years. For the more recent 15 years, South American publications comprised 11% of the total literature published globally. While South America research outputs fluctuated closely with global research trends, differences appeared when comparing areas of growth. This study describes the trends in Chondrichthyan research literature globally and more specifically in South America. Although South American countries may never contribute to the same scale as the wider international scientific community, the future of Chondrichthyans would strongly benefit from the contributions of the many diverse research groups around the world.

Inositol 1,4,5-triphosphate receptor is selectively expressed in cerebellum but not cerebellum-like structures of the elasmobranch fish, Leucoraja erinacea

The Inositol 1,4,5-trisphosphate receptor type 1 protein (Ip3r1) performs an essential role for the induction of cerebellar long-term depression. Here, I describe the use of RT-PCR, qPCR, western blotting and immunohistochemistry to assay Ip3r1 gene expression and localize Ip3r1 protein in the hindbrain of the elasmobranch fish, Leucoraja erinacea. Elasmobranchs are representatives of the most basal, yet extant lineage of gnathostomes, or jawed vertebrates. The cerebellum is a synapomorphy for gnathostomes and thus elasmobranch cerebellar physiology may serve as a proxy for the ancestral state of other jawed vertebrates. LeIp3r1 is selectively expressed in the cerebellum of the little skate and the resultant protein is localized to Purkinje cells. If Ip3r1 performs the same functions in the skate cerebellum as in the mammalian cerebellum, then parallel fiber-Purkinje cell long-term depression through Ip3r1 mediated intracellular calcium regulation may be a conserved feature of cerebellar physiology. Cerebellum and surrounding hindbrain regions termed cerebellum-like structures share a common developmental genetic toolkit. LeIp3r1 expression was lowly detected in cerebellum-like structures indicating that although generatively homologous, the cerebellum and cerebellum-like structures do not share a complete overlap of common expression. Because of the little skate's important phylogenetic placement, performing molecular methodologies to assay targeted gene expression and determine protein localization in the hindbrain can be valuable for our understanding of cerebellar evolution and comparative neural development.

The Chemosensory Receptor Repertoire of a True Shark Is Dominated by a Single Olfactory Receptor Family

Throughout the animal kingdom chemical senses are one of the primary means by which organisms make sense of their environment. To achieve perception of complex chemosensory stimuli large repertoires of olfactory and gustatory receptors are employed in bony vertebrates, which are characterized by high evolutionary dynamics in receptor repertoire size and composition. However, little is known about their evolution in earlier diverging vertebrates such as cartilaginous fish, which include sharks, skates, rays, and chimeras. Recently, the olfactory repertoire of a chimera, elephant shark, was found to be curiously reduced in odorant receptor number. Elephant sharks rely heavily on electroreception to localize prey; thus, it is unclear how representative their chemosensory receptor repertoire sizes would be for cartilaginous fishes in general. Here, we have mined the genome of a true shark, Scyliorhinus canicula (catshark) for olfactory and gustatory receptors, and have performed a thorough phylogenetic study to shed light on the evolution of chemosensory receptors in cartilaginous fish. We report the presence of several gustatory receptors of the TAS1R family in catshark and elephant shark, whereas TAS2R receptors are absent. The catshark olfactory repertoire is dominated by V2R receptors, with 5–8 receptors in the other three families (OR, ORA, TAAR). Species-specific expansions are mostly limited to the V2R family. Overall, the catshark chemosensory receptor repertoires are generally similar in size to those of elephant shark, if somewhat larger, showing similar evolutionary tendencies across over 400 Myr of separate evolution between catshark and elephant shark.

PIRSitePredict for protein functional site prediction using position-specific rules

Methods focused on predicting 'global' annotations for proteins (such as molecular function, biological process and presence of domains or membership in a family) have reached a relatively mature stage. Methods to provide fine-grained 'local' annotation of functional sites (at the level of individual amino acid) are now coming to the forefront, especially in light of the rapid accumulation of genetic variant data. We have developed a computational method and workflow that predicts functional sites within proteins using position-specific conditional template annotation rules (namely PIR Site Rules or PIRSRs for short). Such rules are curated through review of known protein structural and other experimental data by structural biologists and are used to generate high-quality annotations for the UniProt Knowledgebase (UniProtKB) unreviewed section. To share the PIRSR functional site prediction method with the broader scientific community, we have streamlined our workflow and developed a stand-alone Java software package named PIRSitePredict. We demonstrate the use of PIRSitePredict for functional annotation of de novo assembled genome/transcriptome by annotating uncharacterized proteins from Trinity RNA-seq assembly of embryonic transcriptomes of the following three cartilaginous fishes: Leucoraja erinacea (Little Skate), Scyliorhinus canicula (Small-spotted Catshark) and Callorhinchus milii (Elephant Shark). On average about 1200 lines of annotations were predicted for each species.

Potential Human Health Applications from Marine Biomedical Research with Elasmobranch Fishes

Members of the subclass of fishes collectively known as elasmobranchs (Class Chondrichthyes, Subclass Elasmobranchii) include sharks, skates, rays, guitarfish, and sawfish. Having diverged from the main line of vertebrate evolution some 400 million years ago, these fishes have continued to be successful in our ever-changing oceans. Much of their success must be attributed to their uncanny ability to remain healthy. Based on decades of basic research, some of their secrets may be very close to benefitting man. In this short review, some of the molecular and cellular biological areas that show promise for potential human applications are presented. With a brief background and current status of relevant research, these topics include development of new antibiotics and novel treatments for cancer, macular degeneration, viral pathogens, and Parkinson’s disease; potentially useful genomic information from shark transcriptomes; shark antibody-derived drug delivery systems; and immune cell-derived compounds as potential cancer therapeutic agents.

Phylotranscriptomics suggests the jawed vertebrate ancestor could generate diverse helper and regulatory T cell subsets

BACKGROUND

The cartilaginous fishes diverged from other jawed vertebrates ~ 450 million years ago (mya). Despite this key evolutionary position, the only high-quality cartilaginous fish genome available is for the elephant shark (Callorhinchus milii), a chimaera whose ancestors split from the elasmobranch lineage ~ 420 mya. Initial analysis of this resource led to proposals that key components of the cartilaginous fish adaptive immune system, most notably their array of T cell subsets, was primitive compared to mammals. This proposal is at odds with the robust, antigen-specific antibody responses reported in elasmobranchs following immunization. To explore this discrepancy, we generated a multi-tissue transcriptome for small-spotted catshark (Scyliorhinus canicula), a tractable elasmobranch model for functional studies. We searched this, and other newly available sequence datasets, for CD4+ T cell subset-defining genes, aiming to confirm the presence or absence of each subset in cartilaginous fishes.

RESULTS

We generated a new transcriptome based on a normalised, multi-tissue RNA pool, aiming to maximise representation of tissue-specific and lowly expressed genes. We utilized multiple transcriptomic datasets and assembly variants in phylogenetic reconstructions to unambiguously identify several T cell subset-specific molecules in cartilaginous fishes for the first time, including interleukins, interleukin receptors, and key transcription factors. Our results reveal the inability of standard phylogenetic reconstruction approaches to capture the site-specific evolutionary processes of fast-evolving immune genes but show that site-heterogeneous mixture models can adequately do so.

CONCLUSIONS

Our analyses reveal that cartilaginous fishes are capable of producing a range of CD4+ T cell subsets comparable to that of mammals. Further, that the key molecules required for the differentiation and functioning of these subsets existed in the jawed vertebrate ancestor. Additionally, we highlight the importance of considering phylogenetic diversity and, where possible, utilizing multiple datasets for individual species, to accurately infer gene presence or absence at higher taxonomic levels.

An ancient Turing-like patterning mechanism regulates skin denticle development in sharks

Vertebrates have a vast array of epithelial appendages, including scales, feathers, and hair. The developmental patterning of these diverse structures can be theoretically explained by Alan Turing's reaction-diffusion system. However, the role of this system in epithelial appendage patterning of early diverging lineages (compared to tetrapods), such as the cartilaginous fishes, is poorly understood. We investigate patterning of the unique tooth-like skin denticles of sharks, which closely relates to their hydrodynamic and protective functions. We demonstrate through simulation models that a Turing-like mechanism can explain shark denticle patterning and verify this system using gene expression analysis and gene pathway inhibition experiments. This mechanism bears remarkable similarity to avian feather patterning, suggesting deep homology of the system. We propose that a diverse range of vertebrate appendages, from shark denticles to avian feathers and mammalian hair, use this ancient and conserved system, with slight genetic modulation accounting for broad variations in patterning.

Haptoglobin Is a Divergent MASP Family Member That Neofunctionalized To Recycle Hemoglobin via CD163 in Mammals

In mammals, haptoglobin (Hp) is an acute-phase plasma protein that binds with high affinity to hemoglobin (Hb) released by intravascular hemolysis. The resultant Hp-Hb complexes are bound and cleared by the scavenger receptor CD163, limiting Hb-induced oxidative damage. In this study, we show that Hp is a divergent member of the complement-initiating MASP family of proteins, which emerged in the ancestor of jawed vertebrates. We demonstrate that Hp has been independently lost from multiple vertebrate lineages, that characterized Hb-interacting residues of mammals are poorly conserved in nonmammalian species maintaining Hp, and that the extended loop 3 region of Hp, which mediates CD163 binding, is present only in mammals. We show that the Hb-binding ability of cartilaginous fish (nurse shark, Ginglymostoma cirratum; small-spotted catshark, Scyliorhinus canicula; and thornback ray, Raja clavata) and teleost fish (rainbow trout, Oncorhynchus mykiss) Hp is species specific, and where binding does occur it is likely mediated through a different structural mechanism to mammalian Hp. The continued, high-level expression of Hp in cartilaginous fishes in which Hb binding is not evident signals that Hp has (an)other, yet unstudied, role(s) in these species. Previous work indicates that mammalian Hp also has secondary, immunomodulatory functions that are independent of Hb binding; our work suggests these may be remnants of evolutionary more ancient functions, retained after Hb removal became the primary role of Hp in mammals.

Embryonic transcriptome sequencing of the ocellate spot skate Okamejei kenojei

Chondrichthyans (cartilaginous fishes) exhibit highly variable reproductive styles, categorized as viviparity and oviparity. Among these, species with oviparity provide an enormous potential of molecular experimentation with stable sample supply which does not demand the sacrifices of live mothers. Cartilaginous fishes are divided into two subclasses, chimaeras (Holocephali) and elasmobranchs (Elasmobranchii), and the latter consists of two monophyletic groups, Batoidea (rays, skates and torpedoes) and Selachimorpha (sharks). Here we report transcriptome assemblies of the ocellate spot skate Okamejei kenojei, produced by strand-specific RNA-seq of its embryonic tissues. We obtained a total of 325 million illumina short reads from libraries prepared using four different tissue domains and assembled them all together. Our assembly result confirmed the species authenticity and high continuity of contig sequences. Also, assessment of its coverage of pre-selected one-to-one orthologs supported high diversity of transcripts in the assemblies. Our products are expected to provide a basis of comparative molecular studies encompassing other chondrichthyan species with emerging genomic and transcriptomic sequence information.

Shark genomes provide insights into elasmobranch evolution and the origin of vertebrates

Modern cartilaginous fishes are divided into elasmobranchs (sharks, rays and skates) and chimaeras, and the lack of established whole-genome sequences for the former has prevented our understanding of early vertebrate evolution and the unique phenotypes of elasmobranchs. Here we present de novo whole-genome assemblies of brownbanded bamboo shark and cloudy catshark and an improved assembly of the whale shark genome. These relatively large genomes (3.8-6.7 Gbp) contain sparse distributions of coding genes and regulatory elements and exhibit reduced molecular evolutionary rates. Our thorough genome annotation revealed Hox C genes previously hypothesized to have been lost, as well as distinct gene repertories of opsins and olfactory receptors that would be associated with adaptation to unique underwater niches. We also show the early establishment of the genetic machinery governing mammalian homoeostasis and reproduction at the jawed vertebrate ancestor. This study, supported by genomic, transcriptomic and epigenomic resources, provides a foundation for the comprehensive, molecular exploration of phenotypes unique to sharks and insights into the evolutionary origins of vertebrates.

A role for the CXCR4-CXCL12 axis in the little skate, Leucoraja erinacea

The interaction between C-X-C chemokine receptor type 4 (CXCR4) and its cognate ligand C-X-C motif chemokine ligand 12 (CXCL12) plays a critical role in regulating hematopoietic stem cell activation and subsequent cellular mobilization. Extensive studies of these genes have been conducted in mammals, but much less is known about the expression and function of CXCR4 and CXCL12 in non-mammalian vertebrates. In the present study, we identify simultaneous expression of CXCR4 and CXCL12 orthologs in the epigonal organ (the primary hematopoietic tissue) of the little skate, Leucoraja erinacea. Genetic and phylogenetic analyses were functionally supported by significant mobilization of leukocytes following administration of Plerixafor, a CXCR4 antagonist and clinically important drug. Our results provide evidence that, as in humans, Plerixafor disrupts CXCR4/CXCL12 binding in the little skate, facilitating release of leukocytes into the bloodstream. Our study illustrates the value of the little skate as a model organism, particularly in studies of hematopoiesis and potentially for preclinical research on hematological and vascular disorders.

The nanos1 gene was duplicated in early Vertebrates and the two paralogs show different gonadal expression profiles in a shark

Nanos are RNA-binding proteins playing crucial roles in germ cell development and maintenance. Based on phylogenetic and synteny analyses, this study reveals that nanos1 gene has undergone multiple duplications and gene copies losses in Vertebrates. Chondrichthyan species display two nanos1 genes (named nanos1A/1B), which were both retrieved in some Osteichthyes at basal positions in Sarcopterygii and Actinopterygii lineages. In contrast, Teleosts have lost nanos1A but duplicated nanos1B leading to the emergence of two ohnologs (nanos1Ba/1Bb), whereas Tetrapods have lost nanos1B gene. The two successive nanos gene duplications may result from the second and third whole genome duplication events at the basis of Vertebrates and Teleosts respectively. The expression profiles of nanos1A and nanos1B paralogs were characterized in the dogfish, Scyliorhinus canicula. Nanos1A was strongly expressed in brain and also localized in all germ cell types in the polarized testis. In contrast, nanos1B was detected in testis with the highest expression in the germinative zone. In addition, Nanos1B protein was predominantly located in the nuclei of male germinal cells. In the ovary, both paralogs were detected in germinal and somatic cells. Our study opens new perspectives concerning the complex evolution of nanos1 paralogs and their potential distinct roles in Vertebrates gonads.

A staging table for the embryonic development of the brownbanded bamboo shark (Chiloscyllium punctatum)

BACKGROUND

Studying cartilaginous fishes (chondrichthyans) has helped us understand vertebrate evolution and diversity. However, resources such as genome sequences, embryos, and detailed staging tables are limited for species within this clade. To overcome these limitations, we have focused on a species, the brownbanded bamboo shark (Chiloscyllium punctatum), which is a relatively common aquarium species that lays eggs continuously throughout the year. In addition, because of its relatively small genome size, this species is promising for molecular studies.

RESULTS

To enhance biological studies of cartilaginous fishes, we establish a normal staging table for the embryonic development of the brownbanded bamboo shark. Bamboo shark embryos take around 118 days to reach the hatching period at 25°C, which is approximately 1.5 times as fast as the small-spotted catshark (Scyliorhinus canicula) takes. Our staging table divides the embryonic period into 38 stages. Furthermore, we found culture conditions that allow early embryos to grow in partially opened egg cases.

CONCLUSIONS

In addition to the embryonic staging table, we show that bamboo shark embryos exhibit relatively fast embryonic growth and are amenable to culture, key characteristics that enhance their experimental utility. Therefore, the present study is a foundation for cartilaginous fish research. Developmental Dynamics 247:712-723, 2018. © 2017 Wiley Periodicals, Inc.

The Ancient Origins of Neural Substrates for Land Walking

Walking is the predominant locomotor behavior expressed by land-dwelling vertebrates, but it is unknown when the neural circuits that are essential for limb control first appeared. Certain fish species display walking-like behaviors, raising the possibility that the underlying circuitry originated in primitive marine vertebrates. We show that the neural substrates of bipedalism are present in the little skate Leucoraja erinacea, whose common ancestor with tetrapods existed ∼420 million years ago. Leucoraja exhibits core features of tetrapod locomotor gaits, including left-right alternation and reciprocal extension-flexion of the pelvic fins. Leucoraja also deploys a remarkably conserved Hox transcription factor-dependent program that is essential for selective innervation of fin/limb muscle. This network encodes peripheral connectivity modules that are distinct from those used in axial muscle-based swimming and has apparently been diminished in most modern fish. These findings indicate that the circuits that are essential for walking evolved through adaptation of a genetic regulatory network shared by all vertebrates with paired appendages. VIDEO ABSTRACT.

Evolution and Expression of Tissue Globins in Ray-Finned Fishes

The globin gene family encodes oxygen-binding hemeproteins conserved across the major branches of multicellular life. The origins and evolutionary histories of complete globin repertoires have been established for many vertebrates, but there remain major knowledge gaps for ray-finned fish. Therefore, we used phylogenetic, comparative genomic and gene expression analyses to discover and characterize canonical “non-blood” globin family members (i.e., myoglobin, cytoglobin, neuroglobin, globin-X, and globin-Y) across multiple ray-finned fish lineages, revealing novel gene duplicates (paralogs) conserved from whole genome duplication (WGD) and small-scale duplication events. Our key findings were that: (1) globin-X paralogs in teleosts have been retained from the teleost-specific WGD, (2) functional paralogs of cytoglobin, neuroglobin, and globin-X, but not myoglobin, have been conserved from the salmonid-specific WGD, (3) triplicate lineage-specific myoglobin paralogs are conserved in arowanas (Osteoglossiformes), which arose by tandem duplication and diverged under positive selection, (4) globin-Y is retained in multiple early branching fish lineages that diverged before teleosts, and (5) marked variation in tissue-specific expression of globin gene repertoires exists across ray-finned fish evolution, including several previously uncharacterized sites of expression. In this respect, our data provide an interesting link between myoglobin expression and the evolution of air breathing in teleosts. Together, our findings demonstrate great-unrecognized diversity in the repertoire and expression of nonblood globins that has arisen during ray-finned fish evolution.

Evolutionary history of the T cell receptor complex as revealed by small-spotted catshark (Scyliorhinus canicula)

In every jawed vertebrate species studied so far, the T cell receptor (TCR) complex is composed of two different TCR chains (α/β or γ/δ) and a number of CD3 subunits responsible for transmitting signals into the T cell. In this study, we characterised all of the TCR and CD3 genes of small-spotted catshark (Scyliorhinus canicula) and analysed their expression in a broad range of tissues. While the TCR complex is highly conserved across jawed vertebrates, we identified a number of differences in catshark, most notably the presence of two copies of both TCRβ and CD3γδ, and the absence of a functionally-important proline rich region from CD3ε. We also demonstrate that TCRβ has duplicated independently multiple times in jawed vertebrate evolution, bringing additional diversity to the TCR complex. This study reveals new insights about the evolutionary history of the TCR complex and raises new avenues for future exploration.

Evolution of specificity in cartilaginous fish glycoprotein hormones and receptors

Glycoprotein hormones (GpH) interact very specifically with their receptors to mediate hypothalamic-pituitary-peripheral gland endocrine signaling. Vertebrates typically have three functionally distinct GpH endocrine signaling complexes: follicle-stimulating hormone, luteinizing hormone, and thyroid-stimulating hormone, and their receptors. Each hormone consists of a common α subunit bound to one of three different β subunits. Individual hormone subunits and receptors are present in genomes of early metazoans, and a subset of hormone subunits and receptors has been recently characterized in sea lamprey. However, it remains unclear when the full complement of hormone and receptor protein families first appeared, and when specificity of interactions between GpH hormones and receptors first evolved. Here we present phylogenetic analyses showing that the elephant shark (Callorhinchus milii) genome contains sequences representing the current diversity of all hormone subunits and receptors in these co-evolving protein families. We examined specificity of hormone and receptor interactions using functional assays testing reporter gene activation by elephant shark follicle-stimulating hormone, luteinizing hormone, and thyroid-stimulating hormone receptors. We show highly specific, dose-responsive hormone interactions for all three complexes. Our results suggest that co-evolution of specificity between proteins in these endocrine signaling complexes occurred prior to the divergence of Chondrichthyes from the chordate lineage.

Developing an ancient epithelial appendage: FGF signalling regulates early tail denticle formation in sharks

Background

Vertebrate epithelial appendages constitute a diverse group of organs that includes integumentary structures such as reptilian scales, avian feathers and mammalian hair. Recent studies have provided new evidence for the homology of integumentary organ development throughout amniotes, despite their disparate final morphologies. These structures develop from conserved molecular signalling centres, known as epithelial placodes. It is not yet certain whether this homology extends beyond the integumentary organs of amniotes, as there is a lack of knowledge regarding their development in basal vertebrates. As the ancient sister lineage of bony vertebrates, extant chondrichthyans are well suited to testing the phylogenetic depth of this homology. Elasmobranchs (sharks, skates and rays) possess hard, mineralised epithelial appendages called odontodes, which include teeth and dermal denticles (placoid scales). Odontodes constitute some of the oldest known vertebrate integumentary appendages, predating the origin of gnathostomes. Here, we used an emerging model shark (Scyliorhinus canicula) to test the hypothesis that denticles are homologous to other placode-derived amniote integumentary organs. To examine the conservation of putative gene regulatory network (GRN) member function, we undertook small molecule inhibition of fibroblast growth factor (FGF) signalling during caudal denticle formation.

Results

We show that during early caudal denticle morphogenesis, the shark expresses homologues of conserved developmental gene families, known to comprise a core GRN for early placode morphogenesis in amniotes. This includes conserved expression of FGFs, sonic hedgehog (shh) and bone morphogenetic protein 4 (bmp4). Additionally, we reveal that denticle placodes possess columnar epithelial cells with a reduced rate of proliferation, a conserved characteristic of amniote skin appendage development. Small molecule inhibition of FGF signalling revealed placode development is FGF dependent, and inhibiting FGF activity resulted in downregulation of shh and bmp4 expression, consistent with the expectation from comparison to the amniote integumentary appendage GRN.

Conclusion

Overall, these findings suggest the core GRN for building vertebrate integumentary epithelial appendages has been highly conserved over 450 million years. This provides evidence for the continuous, historical homology of epithelial appendage placodes throughout jawed vertebrates, from sharks to mammals. Epithelial placodes constitute the shared foundation upon which diverse vertebrate integumentary organs have evolved.

Electronic supplementary material

The online version of this article (doi:10.1186/s13227-017-0071-0) contains supplementary material, which is available to authorized users.

The evolution and functional diversification of the deubiquitinating enzyme superfamily

Ubiquitin and ubiquitin-like molecules are attached to and removed from cellular proteins in a dynamic and highly regulated manner. Deubiquitinating enzymes are critical to this process, and the genetic catalogue of deubiquitinating enzymes expanded greatly over the course of evolution. Extensive functional redundancy has been noted among the 93 members of the human deubiquitinating enzyme (DUB) superfamily. This is especially true of genes that were generated by duplication (termed paralogs) as they often retain considerable sequence similarity. Because complete redundancy in systems should be eliminated by selective pressure, we theorized that many overlapping DUBs must have significant and unique spatiotemporal roles that can be evaluated in an evolutionary context. We have determined the evolutionary history of the entire class of deubiquitinating enzymes, including the sequence and means of duplication for all paralogous pairs. To establish their uniqueness, we have investigated cell-type specificity in developmental and adult contexts, and have investigated the coemergence of substrates from the same duplication events. Our analysis has revealed examples of DUB gene subfunctionalization, neofunctionalization, and nonfunctionalization.

Comparative transcriptomics of elasmobranchs and teleosts highlight important processes in adaptive immunity and regional endothermy

BACKGROUND

Comparative genomic and/or transcriptomic analyses involving elasmobranchs remain limited, with genome level comparisons of the elasmobranch immune system to that of higher vertebrates, non-existent. This paper reports a comparative RNA-seq analysis of heart tissue from seven species, including four elasmobranchs and three teleosts, focusing on immunity, but concomitantly seeking to identify genetic similarities shared by the two lamnid sharks and the single billfish in our study, which could be linked to convergent evolution of regional endothermy.

RESULTS

Across seven species, we identified an average of 10,877 Swiss-Prot annotated genes from an average of 32,474 open reading frames within each species' heart transcriptome. About half of these genes were shared between all species while the remainder included functional differences between our groups of interest (elasmobranch vs. teleost and endotherms vs. ectotherms) as revealed by Gene Ontology (GO) and selection analyses. A repeatedly represented functional category, in both the uniquely expressed elasmobranch genes (total of 259) and the elasmobranch GO enrichment results, involved antibody-mediated immunity, either in the recruitment of immune cells (Fc receptors) or in antigen presentation, including such terms as "antigen processing and presentation of exogenous peptide antigen via MHC class II", and such genes as MHC class II, HLA-DPB1. Molecular adaptation analyses identified three genes in elasmobranchs with a history of positive selection, including legumain (LGMN), a gene with roles in both innate and adaptive immunity including producing antigens for presentation by MHC class II. Comparisons between the endothermic and ectothermic species revealed an enrichment of GO terms associated with cardiac muscle contraction in endotherms, with 19 genes expressed solely in endotherms, several of which have significant roles in lipid and fat metabolism.

CONCLUSIONS

This collective comparative evidence provides the first multi-taxa transcriptomic-based perspective on differences between elasmobranchs and teleosts, and suggests various unique features associated with the adaptive immune system of elasmobranchs, pointing in particular to the potential importance of MHC Class II. This in turn suggests that expanded comparative work involving additional tissues, as well as genome sequencing of multiple elasmobranch species would be productive in elucidating the regulatory and genome architectural hallmarks of elasmobranchs.

Adaptive phenotypic response to climate enabled by epigenetics in a K-strategy species, the fish Leucoraja ocellata (Rajidae)

The relative importance of genetic versus epigenetic changes in adaptive evolution is a hotly debated topic, with studies showing that some species appear to be able to adapt rapidly without significant genetic change. Epigenetic mechanisms may be particularly important for the evolutionary potential of species with long maturation times and low reproductive potential ('K-strategists'), particularly when faced with rapidly changing environmental conditions. Here we study the transcriptome of two populations of the winter skate (Leucoraja ocellata), a typical 'K-strategist', in Atlantic Canada; an endemic population in the southern Gulf of St Lawrence and a large population on the Scotian Shelf. The endemic population has been able to adapt to a 10°C higher water temperature over short evolutionary time (7000 years), dramatically reducing its body size (by 45%) significantly below the minimum maturation size of Scotian Shelf and other populations of winter skate, as well as exhibiting other adaptations in life history and physiology. We demonstrate that the adaptive response to selection has an epigenetic basis, cataloguing 3653 changes in gene expression that may have enabled this species to rapidly respond to the novel environment. We argue that the epigenetic augmentation of species evolutionary potential (its regulation though gene expression) can enable K-strategists to survive and adapt to different environments, and this mechanism may be particularly important for the persistence of sharks, skates and rays in the light of future climate change.

The fin-to-limb transition as the re-organization of a Turing pattern

A Turing mechanism implemented by BMP, SOX9 and WNT has been proposed to control mouse digit patterning. However, its generality and contribution to the morphological diversity of fins and limbs has not been explored. Here we provide evidence that the skeletal patterning of the catshark Scyliorhinus canicula pectoral fin is likely driven by a deeply conserved Bmp–Sox9–Wnt Turing network. In catshark fins, the distal nodular elements arise from a periodic spot pattern of Sox9 expression, in contrast to the stripe pattern in mouse digit patterning. However, our computer model shows that the Bmp–Sox9–Wnt network with altered spatial modulation can explain the Sox9 expression in catshark fins. Finally, experimental perturbation of Bmp or Wnt signalling in catshark embryos produces skeletal alterations which match in silico predictions. Together, our results suggest that the broad morphological diversity of the distal fin and limb elements arose from the spatial re-organization of a deeply conserved Turing mechanism.

Mouse digit patterning is controlled by a Turing network of Bmp, Sox9, and Wnt. Here, Onimaru et al. show that fin patterning in the catshark, Scyliorhinus canicula, is controlled by the same network with a different spatial organization; thus, the Turing network is deeply conserved in limb development.

The Lungfish Transcriptome: A Glimpse into Molecular Evolution Events at the Transition from Water to Land

Lungfish and coelacanths are the only living sarcopterygian fish. The phylogenetic relationship of lungfish to the last common ancestor of tetrapods and their close morphological similarity to their fossil ancestors make this species uniquely interesting. However their genome size, the largest among vertebrates, is hampering the generation of a whole genome sequence. To provide a partial solution to the problem, a high-coverage lungfish reference transcriptome was generated and assembled. The present findings indicate that lungfish, not coelacanths, are the closest relatives to land-adapted vertebrates. Whereas protein-coding genes evolve at a very slow rate, possibly reflecting a “living fossil” status, transposable elements appear to be active and show high diversity, suggesting a role for them in the remarkable expansion of the lungfish genome. Analyses of single genes and gene families documented changes connected to the water to land transition and demonstrated the value of the lungfish reference transcriptome for comparative studies of vertebrate evolution.

Phylogenetic analysis and tissue distribution of elasmobranch glucose transporters and their response to feeding

Elasmobranch diets consist of high quantities of protein and lipids, but very low levels of carbohydrates including glucose. Reflecting this diet, most tissues use lipids and ketone bodies as their main metabolic fuel. However, the rectal gland has been shown to be dependent on glucose as a fuel, so we hypothesized that glucose transporters (GLUTs) would be present and upregulated in the gland during times of activation (e.g. following a meal). In this study, we searched for and identified putative class I GLUTs in three elasmobranchs and a holocephalan using transcriptomes, and used these to reconstruct a Bayesian phylogeny. We determined that each of the four species possessed three of the four class I GLUT sequences, but the identities of the isoforms present in each species differed between the elasmobranchs (GLUT1, 3 and 4) and the holocephalan (GLUT1, 2 and 3). We then used qPCR to measure mRNA levels of these GLUTs in the rectal gland, liver, intestine, and muscle of fed and starved spiny dogfish (Squalus suckleyi). The rectal gland data showed higher mRNA levels of GLUT4 in the starved relative to the fed fish. In the muscle, both GLUT1 and 4 were significantly elevated at 24 h post-feeding, as was the case for GLUT4 in the liver. In the intestine on the other hand, GLUT4 was significantly elevated by 6 h post-feeding, remaining elevated through 48 h. We suggest that GLUT4 has taken on the role of GLUT2 in elasmobranchs as the expression patterns observed in the liver and intestine are representative of GLUT2 in other vertebrates.

An ancient dental gene set governs development and continuous regeneration of teeth in sharks

The evolution of oral teeth is considered a major contributor to the overall success of jawed vertebrates. This is especially apparent in cartilaginous fishes including sharks and rays, which develop elaborate arrays of highly specialized teeth, organized in rows and retain the capacity for life-long regeneration. Perpetual regeneration of oral teeth has been either lost or highly reduced in many other lineages including important developmental model species, so cartilaginous fishes are uniquely suited for deep comparative analyses of tooth development and regeneration. Additionally, sharks and rays can offer crucial insights into the characters of the dentition in the ancestor of all jawed vertebrates. Despite this, tooth development and regeneration in chondrichthyans is poorly understood and remains virtually uncharacterized from a developmental genetic standpoint. Using the emerging chondrichthyan model, the catshark (Scyliorhinus spp.), we characterized the expression of genes homologous to those known to be expressed during stages of early dental competence, tooth initiation, morphogenesis, and regeneration in bony vertebrates. We have found that expression patterns of several genes from Hh, Wnt/β-catenin, Bmp and Fgf signalling pathways indicate deep conservation over ~450 million years of tooth development and regeneration. We describe how these genes participate in the initial emergence of the shark dentition and how they are redeployed during regeneration of successive tooth generations. We suggest that at the dawn of the vertebrate lineage, teeth (i) were most likely continuously regenerative structures, and (ii) utilised a core set of genes from members of key developmental signalling pathways that were instrumental in creating a dental legacy redeployed throughout vertebrate evolution. These data lay the foundation for further experimental investigations utilizing the unique regenerative capacity of chondrichthyan models to answer evolutionary, developmental, and regenerative biological questions that are impossible to explore in classical models.

Calcium activated K⁺ channels in the electroreceptor of the skate confirmed by cloning. Details of subunits and splicing

Elasmobranchs detect small potentials using excitable cells of the ampulla of Lorenzini which have calcium-activated K(+) channels, first described in 1974. A distinctive feature of the outward current in voltage clamped ampullae is its apparent insensitivity to voltage. The sequence of a BK channel α isoform expressed in the ampulla of the skate was characterized. A signal peptide is present at the beginning of the gene. When compared to human isoform 1 (the canonical sequence), the largest difference was absence of a 59 amino acid region from the S8-S9 intra-cellular linker that contains the strex regulatory domain. The ampulla isoform was also compared with the isoform predicted in late skate embryos where strex was also absent. The BK voltage sensors were conserved in both skate isoforms. Differences between the skate and human BK channel included alternative splicing. Alternative splicing occurs at seven previously defined sites that are characteristic for BK channels in general and hair cells in particular. Skate BK sequences were highly similar to the Australian ghost shark and several other vertebrate species. Based on alignment of known BK sequences with the skate genome and transcriptome, there are at least two isoforms of Kcnma1α expressed in the skate. One of the β subunits (β4), which is known to decrease voltage sensitivity, was also identified in the skate genome and transcriptome and in the ampulla. These studies advance our knowledge of BK channels and suggest further studies in the ampulla and other excitable tissues.

Evolution of the highly networked deubiquitinating enzymes USP4, USP15, and USP11

Background

USP4, USP15 and USP11 are paralogous deubiquitinating enzymes as evidenced by structural organization and sequence similarity. Based on known interactions and substrates it would appear that they have partially redundant roles in pathways vital to cell proliferation, development and innate immunity, and elevated expression of all three has been reported in various human malignancies. The nature and order of duplication events that gave rise to these extant genes has not been determined, nor has their functional redundancy been established experimentally at the organismal level.

Methods

We have employed phylogenetic and syntenic reconstruction methods to determine the chronology of the duplication events that generated the three paralogs and have performed genetic crosses to evaluate redundancy in mice.

Results

Our analyses indicate that USP4 and USP15 arose from whole genome duplication prior to the emergence of jawed vertebrates. Despite having lower sequence identity USP11 was generated later in vertebrate evolution by small-scale duplication of the USP4-encoding region. While USP11 was subsequently lost in many vertebrate species, all available genomes retain a functional copy of either USP4 or USP15, and through genetic crosses of mice with inactivating mutations we have confirmed that viability is contingent on a functional copy of USP4 or USP15. Loss of ubiquitin-exchange regulation, constitutive skipping of the seventh exon and neural-specific expression patterns are derived states of USP11. Post-translational modification sites differ between USP4, USP15 and USP11 throughout evolution.

Conclusions

In isolation sequence alignments can generate erroneous USP gene phylogenies. Through a combination of methodologies the gene duplication events that gave rise to USP4, USP15, and USP11 have been established. Although it operates in the same molecular pathways as the other USPs, the rapid divergence of the more recently generated USP11 enzyme precludes its functional interchangeability with USP4 and USP15. Given their multiplicity of substrates the emergence (and in some cases subsequent loss) of these USP paralogs would be expected to alter the dynamics of the networks in which they are embedded.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-015-0511-1) contains supplementary material, which is available to authorized users.

Whole Genome Sequencing of the Asian Arowana (Scleropages formosus) Provides Insights into the Evolution of Ray-Finned Fishes

The Asian arowana (Scleropages formosus) is of commercial importance, conservation concern, and is a representative of one of the oldest lineages of ray-finned fish, the Osteoglossomorpha. To add to genomic knowledge of this species and the evolution of teleosts, the genome of a Malaysian specimen of arowana was sequenced. A draft genome is presented consisting of 42,110 scaffolds with a total size of 708 Mb (2.85% gaps) representing 93.95% of core eukaryotic genes. Using a k-mer-based method, a genome size of 900 Mb was also estimated. We present an update on the phylogenomics of fishes based on a total of 27 species (23 fish species and 4 tetrapods) using 177 orthologous proteins (71,360 amino acid sites), which supports established relationships except that arowana is placed as the sister lineage to all teleost clades (Bayesian posterior probability 1.00, bootstrap replicate 93%), that evolved after the teleost genome duplication event rather than the eels (Elopomorpha). Evolutionary rates are highly heterogeneous across the tree with fishes represented by both slowly and rapidly evolving lineages. A total of 94 putative pigment genes were identified, providing the impetus for development of molecular markers associated with the spectacular colored phenotypes found within this species.

Development of the Synarcual in the Elephant Sharks (Holocephali; Chondrichthyes): Implications for Vertebral Formation and Fusion

The synarcual is a structure incorporating multiple elements of two or more anterior vertebrae of the axial skeleton, forming immediately posterior to the cranium. It has been convergently acquired in the fossil group ‘Placodermi’, in Chondrichthyes (Holocephali, Batoidea), within the teleost group Syngnathiformes, and to varying degrees in a range of mammalian taxa. In addition, cervical vertebral fusion presents as an abnormal pathology in a variety of human disorders. Vertebrae develop from axially arranged somites, so that fusion could result from a failure of somite segmentation early in development, or from later heterotopic development of intervertebral bone or cartilage. Examination of early developmental stages indicates that in the Batoidea and the ‘Placodermi’, individual vertebrae developed normally and only later become incorporated into the synarcual, implying regular somite segmentation and vertebral development. Here we show that in the holocephalan Callorhinchus milii, uniform and regular vertebral segmentation also occurs, with anterior individual vertebra developing separately with subsequent fusion into a synarcual. Vertebral elements forming directly behind the synarcual continue to be incorporated into the synarcual through growth. This appears to be a common pattern through the Vertebrata. Research into human disorders, presenting as cervical fusion at birth, focuses on gene misexpression studies in humans and other mammals such as the mouse. However, in chondrichthyans, vertebral fusion represents the normal morphology, moreover, taxa such Leucoraja (Batoidea) and Callorhinchus (Holocephali) are increasingly used as laboratory animals, and the Callorhinchus genome has been sequenced and is available for study. Our observations on synarcual development in three major groups of early jawed vertebrates indicate that fusion involves heterotopic cartilage and perichondral bone/mineralised cartilage developing outside the regular skeleton. We suggest that chondrichthyans have potential as ideal extant models for identifying the genes involved in these processes, for application to human skeletal heterotopic disorders.