Molecular phylogeny and divergence times of deuterostome animals

Comparative Analysis of Bivalve and Sea Urchin Genetics and Development: Investigating the Dichotomy in Bilateria

This comprehensive review presents a comparative analysis of early embryogenesis in Protostomia and Deuterostomia, the first of which exhibit a mosaic pattern of development, where cells are fated deterministically, while Deuterostomia display a regulatory pattern of development, where the fate of cells is indeterminate. Despite these fundamental differences, there are common transcriptional mechanisms that underline their evolutionary linkages, particularly in the field of functional genomics. By elucidating both conserved and unique regulatory strategies, this review provides essential insights into the comparative embryology and developmental dynamics of these groups. The objective of this review is to clarify the shared and distinctive characteristics of transcriptional regulatory mechanisms. This will contribute to the extensive areas of functional genomics, evolutionary biology and developmental biology, and possibly lay the foundation for future research and discussion on this seminal topic.

Morphology of the wings and attachment apparatus in the evolution of the family Hippoboscidae (Diptera)

Using a complex analysis of the molecular genetics, morphological, and ecological characteristics of Hippoboscidae flies, the phylogenetic structure and trends in the evolution of morphological characters that contribute to the ectoparasitic lifestyle of hippoboscid flies of the north of Eurasia were studied for the first time. The research was carried out on 26 Palearctic species from 10 genera. The analysis of molecular phylogeny revealed the levels of clustering of the family with the species predominantly parasitizing mammals or birds, the time of cluster formation, and the divergence of species in the Palearctic conditions. An independent adaptation to birds occurred in the genera Icosta, Pseudolynchia, Ornithoica, and others. Bird parasites are characterized by bifid tarsal claws, long hooks on pulvilli, and long empodium setae (except genus Ornithoica). Mammalian parasites are characterized by simple tarsal claws, short lobes of hooks on pulvilli, and zones on empodium with short setae. Specialization in empodium and pulvillus morphotypes and wing reduction are higher diverged in mammalian parasites than in bird parasites. The decrease of flight ability and wing reduction independently arose in different subfamilies of Hippoboscidae flies. Our results assume that the tribe Ornithomyini is a paraphyletic group, since, according to the complex of morphological features, the genus Ornithoica can be considered a separate lineage of evolution.

Exploring the evolution and function of Canoe's intrinsically disordered region in linking cell-cell junctions to the cytoskeleton during embryonic morphogenesis

One central question for cell and developmental biologists is defining how epithelial cells can change shape and move during embryonic development without tearing tissues apart. This requires robust yet dynamic connections of cells to one another, via the cell-cell adherens junction, and of junctions to the actin and myosin cytoskeleton, which generates force. The last decade revealed that these connections involve a multivalent network of proteins, rather than a simple linear pathway. We focus on Drosophila Canoe, homolog of mammalian Afadin, as a model for defining the underlying mechanisms. Canoe and Afadin are complex, multidomain proteins that share multiple domains with defined and undefined binding partners. Both also share a long carboxy-terminal intrinsically disordered region (IDR), whose function is less well defined. IDRs are found in many proteins assembled into large multiprotein complexes. We have combined bioinformatic analysis and the use of a series of canoe mutants with early stop codons to explore the evolution and function of the IDR. Our bioinformatic analysis reveals that the IDRs of Canoe and Afadin differ dramatically in sequence and sequence properties. When we looked over shorter evolutionary time scales, we identified multiple conserved motifs. Some of these are predicted by AlphaFold to be alpha-helical, and two correspond to known protein interaction sites for alpha-catenin and F-actin. We next identified the lesions in a series of eighteen canoe mutants, which have early stop codons across the entire protein coding sequence. Analysis of their phenotypes are consistent with the idea that the IDR, including the conserved motifs in the IDR, are critical for protein function. These data provide the foundation for further analysis of IDR function.

Expansion and collapse of VEGF diversity in major clades of the animal kingdom

Together with the platelet-derived growth factors (PDGFs), the vascular endothelial growth factors (VEGFs) form the PDGF/VEGF subgroup among cystine knot growth factors. The evolutionary relationships within this subgroup have not been examined thoroughly to date. Here, we comprehensively analyze the PDGF/VEGF growth factors throughout all animal phyla and propose a phylogenetic tree. Vertebrate whole-genome duplications play a role in expanding PDGF/VEGF diversity, but several limited duplications are necessary to account for the temporal pattern of emergence. The phylogenetically oldest PDGF/VEGF-like growth factor likely featured a C-terminus with a BR3P signature, a hallmark of the modern-day lymphangiogenic growth factors VEGF-C and VEGF-D. Some younger VEGF genes, such as VEGFB and PGF, appeared completely absent in important vertebrate clades such as birds and amphibia, respectively. In contrast, individual PDGF/VEGF gene duplications frequently occurred in fish on top of the known fish-specific whole-genome duplications. The lack of precise counterparts for human genes poses limitations but also offers opportunities for research using organisms that diverge considerably from humans. Sources for the graphical abstract: 326 MYA and older [1]; 72-240 MYA [2]; 235-65 MYA [3].

Chromosome-scale assembly and quantitative trait locus mapping for major economic traits of the Culter alburnus genome using Illumina and PacBio sequencing with Hi-C mapping information

Topmouth culter (Culter alburnus) is an economically important freshwater fish with high nutritional value. However, its potential genetic advantages have not been fully exploited. Therefore, we aimed to determine the genome sequence of C. alburnus and examine quantitative trait loci (QTLs) related to major economic traits. The results showed that 24 pseudochromosomes were anchored by 914.74 Mb of the C. alburnus genome sequence. De novo sequencing identified 31,279 protein-coding genes with an average length of 8507 bp and average coding sequ ence of 1115 bp. In addition, a high-density genetic linkage map consisting of 24 linkage groups was constructed based on 353,532 high-quality single nucleotide polymorphisms and 4,710 bin markers. A total of 28 QTLs corresponding to 11 genes, 26 QTLs corresponding to 11 genes, and 12 QTLs corresponding to 5 genes were identified for sex, intermuscular spine number and body weight traits, respectively. In this study, we assembled an accurate and nearly complete genome of C. alburnus by combining Illumina, PacBio, and high-throughput Chromosome conformation capture (Hi-C) technologies. In addition, we identified QTLs that explained variances in intermuscular spine number, body weight, and sex differences in C. alburnus. These genetic markers or candidate genes associated with growth traits provide a basis for marker-assisted selection in C. alburnus.

Genomic Analysis of Amphioxus Reveals a Wide Range of Fragments Homologous to Viral Sequences

Amphioxus species are considered living fossils and are important in the evolutionary study of chordates and vertebrates. To explore viral homologous sequences, a high-quality annotated genome of the Beihai amphioxus (Branchiostoma belcheri beihai) was examined using virus sequence queries. In this study, 347 homologous fragments (HFs) of viruses were identified in the genome of B. belcheri beihai, of which most were observed on 21 genome assembly scaffolds. HFs were preferentially located within protein-coding genes, particularly in their CDS regions and promoters. A range of amphioxus genes with a high frequency of HFs is proposed, including histone-related genes that are homologous to the Histone H4 or Histone H2B domains of viruses. Together, this comprehensive analysis of viral HFs provides insights into the neglected role of viral integration in the evolution of amphioxus.

Exploring the evolution and function of Canoe’s intrinsically disordered region in linking cell-cell junctions to the cytoskeleton during embryonic morphogenesis

One central question for cell and developmental biologists is defining how epithelial cells can change shape and move during embryonic development without tearing tissues apart. This requires robust yet dynamic connections of cells to one another, via the cell-cell adherens junction, and of junctions to the actin and myosin cytoskeleton, which generates force. The last decade revealed that these connections involve a multivalent network of proteins, rather than a simple linear pathway. We focus on Drosophila Canoe, homolog of mammalian Afadin, as a model for defining the underlying mechanisms. Canoe and Afadin are complex, multidomain proteins that share multiple domains with defined and undefined binding partners. Both also share a long carboxy-terminal intrinsically disordered region (IDR), whose function is less well defined. IDRs are found in many proteins assembled into large multiprotein complexes. We have combined bioinformatic analysis and the use of a series of canoe mutants with early stop codons to explore the evolution and function of the IDR. Our bioinformatic analysis reveals that the IDRs of Canoe and Afadin differ dramatically in sequence and sequence properties. When we looked over shorter evolutionary time scales, we identified multiple conserved motifs. Some of these are predicted by AlphaFold to be alpha-helical, and two correspond to known protein interaction sites for alpha-catenin and F-actin. We next identified the lesions in a series of eighteen canoe mutants, which have early stop codons across the entire protein coding sequence. Analysis of their phenotypes are consistent with the idea that the IDR, including its C-terminal conserved motifs, are important for protein function. These data provide the foundation for further analysis of IDR function.

miR-92b-5p negatively regulates IKK through targeting its ORF region in the innate immune responses of amphioxus (Branchiostoma belcheri)

Cephalochordate (Amphioxus), situated at a key phylogenetic position in the phylum Chordata, serves as a model organism for studying the origin and evolution of the vertebrate innate immune. In this study, five members of precursor miR-92 family (miR-92a-1, miR-92a-2, miR-92b, miR-92c and miR-92d) are identified in Branchiostoma belcheri, and their evolutionary conservation and potential molecular functions in innate immunity are analyzed. Among them, miR-92b-5p was validated in HEK293 cells to target the coding region but not classic 3'UTR of IKK (inhibitor of nuclear factor kappa-B kinase) mRNA, one integral component of MAPK and TLR4 immune signaling. Furthermore, the spatiotemporal expression patterns of miR-92b-5p and IKK were examined in different tissues or different time points (2 h, 4 h, 8 h, 12 h, 24 h and 48 h) post LPS stimulation at RNA and protein level in vivo. The seemingly inverse expression pattern between miR-92b-5p and IKK supports the involvement of miR-92b-5p in Branchiostoma belcheri innate immune response. In conclusion, our work not only illustrates the evolutionary pattern of Branchiostoma belcheri miR-92 family across chordates, but also reveals that miR-92b-5p could target IKK expression to regulate innate immune response.

EchinoDB: an update to the web-based application for genomic and transcriptomic data on echinoderms

BACKGROUND

Here we release a new version of EchinoDB, EchinoDB v2.0 ( https://echinodb.uncc.edu ). EchinoDB is a database of genomic and transcriptomic data on echinoderms. The initial database consisted of groups of 749,397 orthologous and paralogous transcripts arranged in orthoclusters by sequence similarity.

RESULTS

The updated version of EchinoDB includes two new major datasets: the RNA-Seq data of the brittle star Ophioderma brevispinum and the high-quality genomic assembly data of the green sea urchin Lytechinus variegatus. In addition, we enabled keyword searches for annotated data and installed an updated version of Sequenceserver to allow Basic Local Alignment Search Tool (BLAST) searches. The data are downloadable in FASTA format. The first version of EchinoDB appeared in 2016 and was implemented in GO on a local server. The new version has been updated using R Shiny to include new features and improvements in the application. Furthermore, EchinoDB now runs entirely in the cloud for increased reliability and scaling.

CONCLUSION

EchinoDB serves a user base drawn from the fields of phylogenetics, developmental biology, genomics, physiology, neurobiology, and regeneration. As use cases, we illustrate the function of EchinoDB in retrieving components of signaling pathways involved in the tissue regeneration process of different echinoderms, including the emerging model species Ophioderma brevispinum. Moreover, we use EchinoDB to shed light on the conservation of the molecular components involved in two echinoderm-specific phenomena: spicule matrix proteins involved in the formation of stereom endoskeleton and the tensilin protein that contributes to the capacity of the connective tissues to quickly change its mechanical properties. The genes involved in the former had been previously studied in echinoids, while gene sequences involved in the latter had been previously described in holothuroids. Specifically, we ask (a) if the biomineralization-related proteins previously reported only in sea urchins are also present in other, non-echinoid, echinoderms and (b) if tensilin, the protein responsible for the control of stiffness of the mutable collagenous tissue, previously described in sea cucumbers, is conserved across the phylum.

Characterization of myogenic regulatory factors, myod and myf5 from Megalobrama amblycephala and the effect of lipopolysaccharide on satellite cells in skeletal muscle

Myod and Myf5 are muscle-specific basic helix-loop-helix (bHLH) transcription factors that play essential roles in regulating skeletal muscle development and growth. In order to investigate potential function of myod and myf5 of Megalobrama amblycephala, an economically important freshwater fish species, in the present study, we characterized the sequences and expression profiles of M. amblycephala myod and myf5. The open reading frame (ORF) sequences of myod and myf5 encoded 275 and 240 amino acids, respectively, possessing analogous structure with the highly conserved domains, bHLH and C-terminal helix III domains. Spatio-temporal expression patterns revealed that myod and myf5 were predominant in skeletal muscle with the highest expression in white muscle, and the highest at 10 days post-hatching (dph) and the segmentation period, respectively. Furthermore, we evaluated the effects of lipopolysaccharide (LPS) on the expression of muscle-related genes in white and red muscle, and proliferation and differentiation of satellite cells. The myod, myf5 and pax-7 expression generally increased and then decreased with increase of LPS concentration and treatment time in red muscle, while these genes showed inconsistent expression patterns in white muscle. In addition, LPS administration caused the frequency increase of satellite cells in red and white muscle especially at 3 and 7 days after LPS-injection.

Deep conservation and co-option of programmed cell death facilitates evolution of alternative phenotypes at multiple biological levels

Alternative phenotypes, such as polyphenisms and sexual dimorphisms, are widespread in nature and appear at all levels of biological organization, from genes and cells to morphology and behavior. Yet, our understanding of the mechanisms through which alternative phenotypes develop and how they evolve remains understudied. In this review, we explore the association between alternative phenotypes and programmed cell death, a mechanism responsible for the elimination of superfluous cells during development. We discuss the ancient origins and deep conservation of programmed cell death (its function, forms and underlying core regulatory gene networks), and propose that it was co-opted repeatedly to generate alternative phenotypes at the level of cells, tissues, organs, external morphology, and even individuals. We review several examples from across the tree of life to explore the conditions under which programmed cell death is likely to facilitate the evolution of alternative phenotypes.

Exploring the Candidate Terminal Glycan Profile in Neural Regeneration of the Sea Urchin Paracentrotus lividus, Using Lectin Blotting and Mass Spectrometry

Glycans are expressed as conjugates of glycoproteins, glycolipids, and proteoglycans. The huge diversity of glycans on glycoconjugates contributes to many biological processes, from glycan-based molecular recognition to developmental events, such as regeneration in the nervous system. Echinoderms, which have a close phylogenetic relationship with chordates, are an important group of marine invertebrates for body regeneration. Although many major roles of glycans on glycoconjugates are known, their role in the glycosylation profile of the nervous system in sea urchins is poorly understood. In this study, we aimed to determine the terminal glycan profile by lectin blotting and to quantify sialic acids by the capillary liquid chromatography electrospray ionization tandem mass spectrometry system in the nervous tissue of the sea urchin Paracentrotus lividus. We determined the N-acetyl-D-glucosamine, mannose, and sialic acids (mainly α2,3 linked) by lectin blotting and five types of sialic acids (N-glycolylneuraminic acid, N-acetylneuraminic acid, 9-O-acetyl-N-alycolylneuraminic acid, 5-N-acetyl-9-O-acetyl-N-acetylneuraminic acid, and di-O-acetylated-N-alycolylneuraminic acid) by capillary liquid chromatography electrospray ionization tandem mass spectrometry. This potential first description of the terminal glycan profile in the nervous system of the sea urchin is expected to help us understand its role in nervous system development and regeneration.

Hagfish to Illuminate the Developmental and Evolutionary Origins of the Vertebrate Retina

The vertebrate eye is a vital sensory organ that has long fascinated scientists, but the details of how this organ evolved are still unclear. The vertebrate eye is distinct from the simple photoreceptive organs of other non-vertebrate chordates and there are no clear transitional forms of the eye in the fossil record. To investigate the evolution of the eye we can examine the eyes of the most ancient extant vertebrates, the hagfish and lamprey. These jawless vertebrates are in an ideal phylogenetic position to study the origin of the vertebrate eye but data on eye/retina development in these organisms is limited. New genomic and gene expression data from hagfish and lamprey suggest they have many of the same genes for eye development and retinal neurogenesis as jawed vertebrates, but functional work to determine if these genes operate in retinogenesis similarly to other vertebrates is missing. In addition, hagfish express a marker of proliferative retinal cells (Pax6) near the margin of the retina, and adult retinal growth is apparent in some species. This finding of eye growth late into hagfish ontogeny is unexpected given the degenerate eye phenotype. Further studies dissecting retinal neurogenesis in jawless vertebrates would allow for comparison of the mechanisms of retinal development between cyclostome and gnathostome eyes and provide insight into the evolutionary origins of the vertebrate eye.

Study of the glial cytoarchitecture of the developing olfactory bulb of a shark using immunochemical markers of radial glia

During development of the olfactory bulb (OB), glial cells play key roles in axonal guiding/targeting, glomerular formation and synaptic plasticity. Studies in mammals have shown that radial glial cells and peripheral olfactory glia (olfactory ensheathing cells, OECs) are involved in the development of the OB. Most studies about the OB glia were carried out in mammals, but data are lacking in most non-mammalian vertebrates. In the present work, we studied the development of the OB glial system in the cartilaginous fish Scyliorhinus canicula (catshark) using antibodies against glial markers, such as glial fibrillary acidic protein (GFAP), brain lipid-binding protein (BLBP), and glutamine synthase (GS). These glial markers were expressed in cells with radial morphology lining the OB ventricle of embryos and this expression continues in ependymal cells (tanycytes) in early juveniles. Astrocyte-like cells were also observed in the granular layer and surrounding glomeruli. Numerous GS-positive cells were present in the primary olfactory pathway of embryos. In the developmental stages analysed, the olfactory nerve layer and the glomerular layer were the regions with higher GFAP, BLBP and GS immuno-reactivity. In addition, numerous BLBP-expressing cells (a marker of mammalian OECs) showing proliferative activity were present in the olfactory nerve layer. Our findings suggest that glial cells of peripheral and central origin coexist in the OB of catshark embryos and early juveniles. These results open the path for future studies about the differential roles of glial cells in the catshark OB during embryonic development and in adulthood.

Selective Gene Loss of Visual and Olfactory Guanylyl Cyclase Genes Following the Two Rounds of Vertebrate-Specific Whole-Genome Duplications

Photoreceptors convey visual information and come in two flavors; dim-light and bright-light dedicated rod and cones. Both cell types feature highly specialized phototransduction cascades that convert photonic energy into intracellular signals. Although a substantial amount of phototransduction gene ohnologs are expressed either in rods or cones, visual guanylyl cyclases (GCs) involved in the calcium (Ca²⁺) dependent feedback regulation of phototransduction are neither rod nor cone specific. The co-existence of visual GCs in both photoreceptor types suggests that specialization of these ohnologs occurred despite their overlapping expression.

Here, we analyze gene retention and inactivation patterns of vertebrate visual and closely related olfactory GCs following two rounds (2R) of vertebrate-specific whole-genome duplication events (2R WGD). Although eutherians generally use two visual and one olfactory GC, independent inactivation occurred in some lineages. Sauropsids (birds, lizards, snakes, turtles, and crocodiles) generally have only one visual GC (GC-E). Additionally, turtles (testodes) also lost the olfactory GC (GC-D). Pseudogenization in mammals occurred in specific species/families likely according to functional needs (i.e., many species with reduced vision only have GC-E). Likewise, some species not relying on scent marks lack GC-D, the olfactory GC enzyme. Interestingly, in the case of fish, no species can be found with fewer than three (two visual and one olfactory) genes and the teleost-specific 3R WGD can increase this number to up to five. This suggests that vision in fish now requires at least two visual GCs.

Estimating the age of the subfamily Orthocoronavirinae using host divergence times as calibration ages at two internal nodes

Viruses of the subfamily Orthocoronavirinae can cause mild to severe disease in people, including COVID-19, MERS and SARS. Their most common natural hosts are bat and bird species, which are mostly split across four virus genera. Molecular clock analyses of orthocoronaviruses suggested the most recent common ancestor of these viruses might have emerged either around 10,000 years ago or, using models accounting for selection, many millions of years. Here, we reassess the evolutionary history of these viruses. We present time-aware phylogenetic analyses of a RNA-dependent RNA polymerase locus from 123 orthocoronaviruses isolated from birds and bats, including those in New Zealand, which were geographically isolated from other bats around 35 million years ago. We used this age, as well as the age of the avian-mammals split, to calibrate the molecular clocks, under the assumption that these ages are applicable to the analyzed viruses. We found that the time to the most recent ancestor common for all orthocoronaviruses is likely 150 or more million years, supporting clock analyses that account for selection.

The Tentacular Spectacular: Evolution of Regeneration in Sea Anemones

Sea anemones vary immensely in life history strategies, environmental niches and their ability to regenerate. While the sea anemone Nematostella vectensis is the starlet of many key regeneration studies, recent work is emerging on the diverse regeneration strategies employed by other sea anemones. This manuscript will explore current molecular mechanisms of regeneration employed by non-model sea anemones Exaiptasia diaphana (an emerging model species for coral symbiosis studies) and Calliactis polypus (a less well-studied species) and examine how these species compare to the model sea anemone N. vectensis. We summarize the field of regeneration within sea anemones, within the greater context of phylum Cnidaria and in other invertebrate models of regeneration. We also address the current knowledge on two key systems that may be implemented in regeneration: the innate immune system and developmental pathways, including future aspects of work and current limitations.

Nkx2-1 and FoxE regionalize glandular (mucus-producing) and thyroid-equivalent traits in the endostyle of the chordate Oikopleura dioica

The endostyle is a ventral pharyngeal organ used for internal filter feeding of basal chordates and is considered homologous to the follicular thyroid of vertebrates. It contains mucus-producing (glandular) and thyroid-equivalent regions organized along the dorsoventral (DV) axis. Although thyroid-related genes (Nkx2-1, FoxE, and thyroid peroxidase (TPO)) are known to be expressed in the endostyle, their roles in establishing regionalization within the organ have not been demonstrated. We report that Nkx2-1 and FoxE are essential for establishing DV axial identity in the endostyle of Oikopleura dioica. Genome and expression analyses showed von Willebrand factor-like (vWFL) and TPO/dual oxidase (Duox)/Nkx2-1/FoxE as orthologs of glandular and thyroid-related genes, respectively. Knockdown experiments showed that Nkx2-1 is necessary for the expression of glandular and thyroid-related genes, whereas FoxE is necessary only for thyroid-related genes. Moreover, Nkx2-1 expression is necessary for FoxE expression in larvae during organogenesis. The results demonstrate the essential roles of Nkx2-1 and FoxE in establishing regionalization in the endostyle, including (1) the Nkx2-1-dependent glandular region, and (2) the Nkx2-1/FoxE-dependent thyroid-equivalent region. DV axial regionalization may be responsible for organizing glandular and thyroid-equivalent traits of the pharynx along the DV axis.

Evolution of the Avian Eggshell Biomineralization Protein Toolkit - New Insights From Multi-Omics

The avian eggshell is a remarkable biomineral, which is essential for avian reproduction; its properties permit embryonic development in the desiccating terrestrial environment, and moreover, are critically important to preserve unfertilized egg quality for human consumption. This calcium carbonate (CaCO₃) bioceramic is made of 95% calcite and 3.5% organic matrix; it protects the egg contents against microbial penetration and mechanical damage, allows gaseous exchange, and provides calcium for development of the embryonic skeleton. In vertebrates, eggshell occurs in the Sauropsida and in a lesser extent in Mammalia taxa; avian eggshell calcification is one of the fastest known CaCO₃ biomineralization processes, and results in a material with excellent mechanical properties. Thus, its study has triggered a strong interest from the researcher community. The investigation of eggshell biomineralization in birds over the past decades has led to detailed characterization of its protein and mineral constituents. Recently, our understanding of this process has been significantly improved using high-throughput technologies (i.e., proteomics, transcriptomics, genomics, and bioinformatics). Presently, more or less complete eggshell proteomes are available for nine birds, and therefore, key proteins that comprise the eggshell biomineralization toolkit are beginning to be identified. In this article, we review current knowledge on organic matrix components from calcified eggshell. We use these data to analyze the evolution of selected matrix proteins and underline their role in the biological toolkit required for eggshell calcification in avian species. Amongst the panel of eggshell-associated proteins, key functional domains are present such as calcium-binding, vesicle-binding and protein-binding. These technical advances, combined with progress in mineral ultrastructure analyses, have opened the way for new hypotheses of mineral nucleation and crystal growth in formation of the avian eggshell, including transfer of amorphous CaCO₃ in vesicles from uterine cells to the eggshell mineralization site. The enrichment of multi-omics datasets for bird species is critical to understand the evolutionary context for development of CaCO₃ biomineralization in metazoans, leading to the acquisition of the robust eggshell in birds (and formerly dinosaurs).

Accelerated loss of hypoxia response in zebrafish with familial Alzheimer's disease-like mutation of presenilin 1

Ageing is the major risk factor for Alzheimer's disease (AD), a condition involving brain hypoxia. The majority of early-onset familial AD (EOfAD) cases involve dominant mutations in the gene PSEN1. PSEN1 null mutations do not cause EOfAD. We exploited putative hypomorphic and EOfAD-like mutations in the zebrafish psen1 gene to explore the effects of age and genotype on brain responses to acute hypoxia. Both mutations accelerate age-dependent changes in hypoxia-sensitive gene expression supporting that ageing is necessary, but insufficient, for AD occurrence. Curiously, the responses to acute hypoxia become inverted in extremely aged fish. This is associated with an apparent inability to upregulate glycolysis. Wild-type PSEN1 allele expression is reduced in post-mortem brains of human EOfAD mutation carriers (and extremely aged fish), possibly contributing to EOfAD pathogenesis. We also observed that age-dependent loss of HIF1 stabilization under hypoxia is a phenomenon conserved across vertebrate classes.

Whole-genome resequencing reveals the pleistocene temporal dynamics of Branchiostoma belcheri and Branchiostoma floridae populations

Global climatic fluctuations governed the ancestral demographic histories of species and contributed to place the current population status into a more extensive ecological and evolutionary context. Genetic variations will leave unambiguous signatures in the patterns of intraspecific genetic variation in extant species since the genome of each individual is an imperfect mosaic of the ancestral genomes. Here, we report the genome sequences of 20 Branchiostoma individuals by whole-genome resequencing strategy. We detected over 140 million genomic variations for each Branchiostoma individual. In particular, we applied the pairwise sequentially Markovian coalescent (PSMC) method to estimate the trajectories of changes in the effective population size (N e) of Branchiostoma population during the Pleistocene. We evaluated the threshold of sequencing depth for proper inference of demographic histories using PSMC was ≥25×. The PSMC results highlight the role of historical global climatic fluctuations in the long-term population dynamics of Branchiostoma. The inferred ancestral N e of the Branchiostoma belcheri populations from Zhanjiang and Xiamen (China) seawaters was different in amplitude before the first (mutation rate = 3 × 10-9) or third glaciation (mutation rate = 9 × 10-9) of the Pleistocene, indicating that the two populations most probably started to evolve in isolation in their respective seas after the first or third glaciation of the Pleistocene. A pronounced population bottleneck coinciding with the last glacial maximum was observed in all Branchiostoma individuals, followed by a population expansion occurred during the late Pleistocene. Species that have experienced long-term declines may be especially vulnerable to recent anthropogenic activities. Recently, the industrial pollution and the exploitation of sea sand have destroyed the harmonious living environment of amphioxus species. In the future, we need to protect the habitat of Branchiostoma and make full use of these detected genetic variations to facilitate the functional study of Branchiostoma for adaptation to local environments.

An evolutionarily ancient mechanism for regulation of hemoglobin expression in vertebrate red cells

The oxygen transport function of hemoglobin (HB) is thought to have arisen ∼500 million years ago, roughly coinciding with the divergence between jawless (Agnatha) and jawed (Gnathostomata) vertebrates. Intriguingly, extant HBs of jawless and jawed vertebrates were shown to have evolved twice, and independently, from different ancestral globin proteins. This raises the question of whether erythroid-specific expression of HB also evolved twice independently. In all jawed vertebrates studied to date, one of the HB gene clusters is linked to the widely expressed NPRL3 gene. Here we show that the nprl3-linked hb locus of a jawless vertebrate, the river lamprey (Lampetra fluviatilis), shares a range of structural and functional properties with the equivalent jawed vertebrate HB locus. Functional analysis demonstrates that an erythroid-specific enhancer is located in intron 7 of lamprey nprl3, which corresponds to the NPRL3 intron 7 MCS-R1 enhancer of jawed vertebrates. Collectively, our findings signify the presence of an nprl3-linked multiglobin gene locus, which contains a remote enhancer that drives globin expression in erythroid cells, before the divergence of jawless and jawed vertebrates. Different globin genes from this ancestral cluster evolved in the current NPRL3-linked HB genes in jawless and jawed vertebrates. This provides an explanation of the enigma of how, in different species, globin genes linked to the same adjacent gene could undergo convergent evolution.

Transcriptomic analysis of sea cucumber (Holothuria leucospilota) coelomocytes revealed the echinoderm cytokine response during immune challenge

Background

The sea cucumber Holothuria leucospilota belongs to echinoderm, which is evolutionally the most primitive group of deuterostomes. Sea cucumber has a cavity between its digestive tract and the body wall that is filled with fluid and suspended coelomic cells similar to blood cells. The humoral immune response of the sea cucumber is based on the secretion of various immune factors from coelomocytes into the coelomic cavity. The aim of this study is to lay out a foundation for the immune mechanisms in echinoderms and their origins in chordates by using RNA-seq.

Results

Sea cucumber primary coelomocytes were isolated from healthy H. leucospilota and incubated with lipopolysaccharide (LPS, 10 μg/ml), polyinosinic-polycytidylic acid [Poly (I:C), 10 μg/ml] and heat-inactived Vibrio harveyi (10⁷ cell/ml) for 24 h, respectively. After high-throughput mRNA sequencing on an Illumina HiSeq2500, a de novo transcriptome was assembled and the Unigenes were annotated. Thirteen differentially expressed genes (DEGs) were selected randomly from our data and subsequently verified by using RT-qPCR. The results of RT-qPCR were consistent with those of the RNA-seq (R² = 0.61). The top 10 significantly enriched signaling pathways and immune-related pathways of the common and unique DEGs were screened from the transcriptome data. Twenty-one cytokine candidate DEGs were identified, which belong to 4 cytokine families, namely, BCL/CLL, EPRF1, IL-17 and TSP/TPO. Gene expression in response to LPS dose-increased treatment (0, 10, 20 and 50 μg/ml) showed that IL-17 family cytokines were significantly upregulated after 10 μg/ml LPS challenge for 24 h.

Conclusion

A de novo transcriptome was sequenced and assembled to generate the gene expression profiling across the sea cucumber coelomocytes treated with LPS, Poly (I:C) and V. harveyi. The cytokine genes identified in DEGs could be classified into 4 cytokine families, in which the expression of IL-17 family cytokines was most significantly induced after 10 μg/ml LPS challenge for 24 h. Our findings have laid the foundation not only for the research of molecular mechanisms related to the immune response in echinoderms but also for their origins in chordates, particularly in higher vertebrates.

Variation in base composition underlies functional and evolutionary divergence in non-LTR retrotransposons

Background

Non-LTR retrotransposons often exhibit base composition that is markedly different from the nucleotide content of their host’s gene. For instance, the mammalian L1 element is AT-rich with a strong A bias on the positive strand, which results in a reduced transcription. It is plausible that the A-richness of mammalian L1 is a self-regulatory mechanism reflecting a trade-off between transposition efficiency and the deleterious effect of L1 on its host. We examined if the A-richness of L1 is a general feature of non-LTR retrotransposons or if different clades of elements have evolved different nucleotide content. We also investigated if elements belonging to the same clade evolved towards different base composition in different genomes or if elements from different clades evolved towards similar base composition in the same genome.

Results

We found that non-LTR retrotransposons differ in base composition among clades within the same host but also that elements belonging to the same clade differ in base composition among hosts. We showed that nucleotide content remains constant within the same host over extended period of evolutionary time, despite mutational patterns that should drive nucleotide content away from the observed base composition.

Conclusions

Our results suggest that base composition is evolving under selection and may be reflective of the long-term co-evolution between non-LTR retrotransposons and their host. Finally, the coexistence of elements with drastically different base composition suggests that these elements may be using different strategies to persist and multiply in the genome of their host.

Preferential intracellular pH regulation is a common trait amongst fishes exposed to high environmental CO2

Acute (<96 h) exposure to elevated environmental CO2 (hypercarbia) induces a pH disturbance in fishes that is often compensated by concurrent recovery of intracellular and extracellular pH (pHi and pHe, respectively; coupled pH regulation). However, coupled pH regulation may be limited at CO2 partial pressure (PCO2 ) tensions far below levels that some fishes naturally encounter. Previously, four hypercarbia-tolerant fishes had been shown to completely and rapidly regulate heart, brain, liver and white muscle pHi during acute exposure to >4 kPa PCO2  (preferential pHi regulation) before pHe compensation was observed. Here, we test the hypothesis that preferential pHi regulation is a widespread strategy of acid-base regulation among fish by measuring pHi regulation in 10 different fish species that are broadly phylogenetically separated, spanning six orders, eight families and 10 genera. Contrary to previous views, we show that preferential pHi regulation is the most common strategy for acid-base regulation within these fishes during exposure to severe acute hypercarbia and that this strategy is associated with increased hypercarbia tolerance. This suggests that preferential pHi regulation may confer tolerance to the respiratory acidosis associated with hypercarbia, and we propose that it is an exaptation that facilitated key evolutionary transitions in vertebrate evolution, such as the evolution of air breathing.

Molecular evolution of GIP and Exendin and their receptors

Glucose-dependent insulinotropic polypeptide (GIP) is a product of the Gip gene and acts as an incretin hormone in mammals. Gip is most closely related to the proglucagon (Gcg) and Exendin genes and diverged from these very early in vertebrate evolution. In mammals, GIP acts through its specific receptor, encoded by the Gipr gene, which belongs to a subfamily of 7-transmembrane G-protein coupled receptor (GPCR) genes that also includes those for the proglucagon-derived peptides (Gcgr, Glp1r, and Glp2r), and the receptor for Exendin (Grlr). Gip, Gipr, Exendin, and Grlr genes are found in species from most vertebrate classes. While most species that have a Gip gene also have a Gipr gene, two classes of vertebrates, cartilaginous fish and birds, retain conserved Gip genes but lack Gipr genes. This raises the possibility the GIP signals through other receptors in some vertebrates. Exendin genes and the gene for its receptor, Grlr, are also found in diverse vertebrates, with the notable exception of mammals. Both GIP and Exendin likely have important roles in vertebrate physiology, but their roles are either dispensable or can be replaced by other hormones.

Neurons and Glia Cells in Marine Invertebrates: An Update

The nervous system (NS) of invertebrates and vertebrates is composed of two main types of cells: neurons and glia. In both types of organisms, nerve cells have similarities in biochemistry and functionality. The neurons are in charge of the synapse, and the glial cells are in charge of important functions of neuronal and homeostatic modulation. Knowing the mechanisms by which NS cells work is important in the biomedical area for the diagnosis and treatment of neurological disorders. For this reason, cellular and animal models to study the properties and characteristics of the NS are always sought. Marine invertebrates are strategic study models for the biological sciences. The sea slug Aplysia californica and the squid Loligo pealei are two examples of marine key organisms in the neurosciences field. The principal characteristic of marine invertebrates is that they have a simpler NS that consists of few and larger cells, which are well organized and have accessible structures. As well, the close phylogenetic relationship between Chordata and Echinodermata constitutes an additional advantage to use these organisms as a model for the functionality of neuronal cells and their cellular plasticity. Currently, there is great interest in analyzing the signaling processes between neurons and glial cells, both in vertebrates and in invertebrates. However, only few types of glial cells of invertebrates, mostly insects, have been studied, and it is important to consider marine organisms' research. For this reason, the objective of the review is to present an update of the most relevant information that exists around the physiology of marine invertebrate neuronal and glial cells.

Primary orthologs from local sequence context

BACKGROUND

The evolutionary history of genes serves as a cornerstone of contemporary biology. Most conserved sequences in mammalian genomes don't code for proteins, yielding a need to infer evolutionary history of sequences irrespective of what kind of functional element they may encode. Thus, sequence-, as opposed to gene-, centric modes of inferring paths of sequence evolution are increasingly relevant. Customarily, homologous sequences derived from the same direct ancestor, whose ancestral position in two genomes is usually conserved, are termed "primary" (or "positional") orthologs. Methods based solely on similarity don't reliably distinguish primary orthologs from other homologs; for this, genomic context is often essential. Context-dependent identification of orthologs traditionally relies on genomic context over length scales characteristic of conserved gene order or whole-genome sequence alignment, and can be computationally intensive.

RESULTS

We demonstrate that short-range sequence context-as short as a single "maximal" match- distinguishes primary orthologs from other homologs across whole genomes. On mammalian whole genomes not preprocessed by repeat-masker, potential orthologs are extracted by genome intersection as "non-nested maximal matches:" maximal matches that are not nested into other maximal matches. It emerges that on both nucleotide and gene scales, non-nested maximal matches recapitulate primary or positional orthologs with high precision and high recall, while the corresponding computation consumes less than one thirtieth of the computation time required by commonly applied whole-genome alignment methods. In regions of genomes that would be masked by repeat-masker, non-nested maximal matches recover orthologs that are inaccessible to Lastz net alignment, for which repeat-masking is a prerequisite. mmRBHs, reciprocal best hits of genes containing non-nested maximal matches, yield novel putative orthologs, e.g. around 1000 pairs of genes for human-chimpanzee.

CONCLUSIONS

We describe an intersection-based method that requires neither repeat-masking nor alignment to infer evolutionary history of sequences based on short-range genomic sequence context. Ortholog identification based on non-nested maximal matches is parameter-free, and less computationally intensive than many alignment-based methods. It is especially suitable for genome-wide identification of orthologs, and may be applicable to unassembled genomes. We are agnostic as to the reasons for its effectiveness, which may reflect local variation of mean mutation rate.

A new approach to molecular biosurveillance of invasive species using DNA metabarcoding

Non-indigenous species (NIS) reach every corner of the world, at times wreaking havoc on ecosystems and costing the global economy billions of dollars. A rapid and accurate biosurveillance tool tailored to a particular biogeographic region is needed to detect NIS when they are first introduced into an area as traditional detection methods are expensive and require specialized expertise. Metabarcoding of environmental and community DNA meets those biosurveillance requirements; a novel tool tailored to the Northwest Pacific Ocean is presented here using an approach that could revolutionize early detection of NIS. Eight newly designed genetic markers for multiple gene regions were implemented to meet the stringent taxonomic requirements for the detection of NIS across four major marine phyla. The tool was considered highly successful because it identified 12 known NIS in the study area and a further seven species representing potential new records. Overall community composition detected here was statistically different between substrate types; zooplankton sampling accounted for significantly higher species richness than filtered sea water in most cases, but this was dominated by mollusk and arthropod species. Both substrate types sampled were required to identify the wide taxonomic breadth of known NIS in the study area. Intensive sampling is known to be paramount for the detection of rare species, including new incursions of NIS, thus it is recommended to include diverse DNA sampling protocols based on species' life-history characteristics for broad detection capacity. Application of a metabarcoding-based molecular biosurveillance tool optimized for biogeographic regions enables rapid and accurate early detection across a wide taxonomic range to allow quick implementation of eradication or control efforts and potentially mitigate some of the devastating effects of NIS worldwide.

The stalk-eyed fly as a model for aggression - is there a conserved role for 5-HT between vertebrates and invertebrates?

Serotonin (5-HT) has largely been accepted to be inhibitory to vertebrate aggression, whereas an opposing stimulatory role has been proposed for invertebrates. Herein, we argue that critical gaps in our understanding of the nuanced role of 5-HT in invertebrate systems drove this conclusion prematurely, and that emerging data suggest a previously unrecognized level of phylogenetic conservation with respect to neurochemical mechanisms regulating the expression of aggressive behaviors. This is especially apparent when considering the interplay among factors governing 5-HT activity, many of which share functional homology across taxa. We discuss recent findings using insect models, with an emphasis on the stalk-eyed fly, to demonstrate how particular 5-HT receptor subtypes mediate the intensity of aggression with respect to discrete stages of the interaction (initiation, escalation and termination), which mirrors the complex behavioral regulation currently recognized in vertebrates. Further similarities emerge when considering the contribution of neuropeptides, which interact with 5-HT to ultimately determine contest progression and outcome. Relative to knowledge in vertebrates, much less is known about the function of 5-HT receptors and neuropeptides in invertebrate aggression, particularly with respect to sex, species and context, prompting the need for further studies. Our Commentary highlights the need to consider multiple factors when determining potential taxonomic differences, and raises the possibility of more similarities than differences between vertebrates and invertebrates with regard to the modulatory effect of 5-HT on aggression.

Evolution of Neuropeptide Precursors in Polyneoptera (Insecta)

Neuropeptides are among the structurally most diverse signaling molecules and participate in intercellular information transfer from neurotransmission to intrinsic or extrinsic neuromodulation. Many of the peptidergic systems have a very ancient origin that can be traced back to the early evolution of the Metazoa. In recent years, new insights into the evolution of these peptidergic systems resulted from the increasing availability of genome and transcriptome data which facilitated the investigation of the complete neuropeptide precursor sequences. Here we used a comprehensive transcriptome dataset of about 200 species from the 1KITE initiative to study the evolution of single-copy neuropeptide precursors in Polyneoptera. This group comprises well-known orders such as cockroaches, termites, locusts, and stick insects. Due to their phylogenetic position within the insects and the large number of old lineages, these insects are ideal candidates for studying the evolution of insect neuropeptides and their precursors. Our analyses include the orthologs of 21 single-copy neuropeptide precursors, namely ACP, allatotropin, AST-CC, AST-CCC, CCAP, CCHamide-1 and 2, CNMamide, corazonin, CRF-DH, CT-DH, elevenin, HanSolin, NPF-1 and 2, MS, proctolin, RFLamide, SIFamide, sNPF, and trissin. Based on the sequences obtained, the degree of sequence conservation between and within the different polyneopteran lineages is discussed. Furthermore, the data are used to postulate the individual neuropeptide sequences that were present at the time of the insect emergence more than 400 million years ago. The data confirm that the extent of sequence conservation across Polyneoptera is remarkably different between the different neuropeptides. Furthermore, the average evolutionary distance for the single-copy neuropeptides differs significantly between the polyneopteran orders. Nonetheless, the single-copy neuropeptide precursors of the Polyneoptera show a relatively high degree of sequence conservation. Basic features of these precursors in this very heterogeneous insect group are explained here in detail for the first time.

Characterization of testis-specific serine/threonine kinase 1-like (TSSK1-like) gene and expression patterns in diploid and triploid Pacific abalone (Haliotis discus hannai; Gastropoda; Mollusca) males

Testis-specific serine/threonine kinase 1-like (TSSK1-like), which plays important roles in late-phase spermatogenesis and male fertility, was characterized in Pacific abalone Haliotis discus hannai, an important commercial marine gastropod. Further, its expression patterns were assessed in diploid and induced triploid males showing differential degrees of testis maturation. Abalone TSSK1-like shared conserved structural features with mammalian TSSK1s and other potential metazoan orthologs, especially regarding the catalytic STKc domain. Phylogenetically, abalone TSSK1-like displayed a genetic affiliation with its molluscan TSSK1-like orthologs and human TSSK1. Additionally, abalone TSSK1-like gene showed a tetrapartite exon-intron organization, unlike the intronless structure of most amniotic tetrapodian TSSK1s. Molecular phylogenetic analysis in the metazoan lineage suggested a possible revision in the origin of the earliest ancestral TSSK1. Further, abalone TSSK1-like showed testis-predominant expression, which was significantly influenced by both age and seasonal reproductive cycles. Comparative expression analyses between diploid and triploid abalone males suggested that robust TSSK1-like expression occurred primarily at the post-meiotic stage. Additionally, RT-PCR assay indicates that mature abalone sperms retain TSSK1-like transcripts after release. Taken together, this study provides useful insights for further studies to assess male reproduction and sterility and/or partial fertility of induced male triploidy in abalone species.

The Melanin-Concentrating Hormone (MCH) System: A Tale of Two Peptides

The melanin-concentrating hormone (MCH) system is a robust integrator of exogenous and endogenous information, modulating arousal and energy balance in mammals. Its predominant function in teleosts, however, is to concentrate melanin in the scales, contributing to the adaptive color change observed in several teleost species. These contrasting functions resulted from a gene duplication that occurred after the teleost divergence, which resulted in the generation of two MCH-coding genes in this clade, which acquired distinctive sequences, distribution, and functions, examined in detail here. We also describe the distribution of MCH immunoreactivity and gene expression in a large number of species, in an attempt to identify its core elements. While initially originated as a periventricular peptide, with an intimate relationship with the third ventricle, multiple events of lateral migration occurred during evolution, making the ventrolateral and dorsolateral hypothalamus the predominant sites of MCH in teleosts and mammals, respectively. Substantial differences between species can be identified, likely reflecting differences in habitat and behavior. This observation aligns well with the idea that MCH is a major integrator of internal and external information, ensuring an appropriate response to ensure the organism’s homeostasis. New studies on the MCH system in species that have not yet been investigated will help us understand more precisely how these habitat changes are connected to the hypothalamic neurochemical circuits, paving the way to new intervention strategies that may be used with pharmacological purposes.

Steroid receptors and vertebrate evolution

Considering that life on earth evolved about 3.7 billion years ago, vertebrates are young, appearing in the fossil record during the Cambrian explosion about 542 to 515 million years ago. Results from sequence analyses of genomes from bacteria, yeast, plants, invertebrates and vertebrates indicate that receptors for adrenal steroids (aldosterone, cortisol), and sex steroids (estrogen, progesterone, testosterone) also are young, with an estrogen receptor and a 3-ketosteroid receptor first appearing in basal chordates (cephalochordates: amphioxus), which are close ancestors of vertebrates. Duplication and divergence of the 3-ketosteroid receptor yielded an ancestral progesterone receptor and an ancestral corticoid receptor, the common ancestor of the glucocorticoid and mineralocorticoid receptors, in jawless vertebrates (cyclostomes: lampreys, hagfish). This was followed by evolution of an androgen receptor, distinct glucocorticoid and mineralocorticoid receptors and estrogen receptor-α and -β in cartilaginous fishes (Chondrichthyes: sharks). Further evolution of mineralocorticoid signaling occurred with the evolution of aldosterone synthase in lungfish, a forerunner of terrestrial vertebrates. Adrenal and sex steroid receptors are not found in echinoderms and hemichordates, which are ancestors in the lineage of cephalochordates and vertebrates. The evolution of steroid receptors at key nodes in the evolution of vertebrates, in which steroid receptors act as master switches to regulate differentiation, development, reproduction, immune responses, electrolyte homeostasis and stress responses, suggests an important role for steroid receptors in the evolutionary success of vertebrates, considering that the human genome contains about 22,000 genes, which is not much larger than genomes of invertebrates, such as Caenorhabditis elegans (~18,000 genes) and Drosophila (~14,000 genes).

Effect of seasonal high temperature on the immune response in Apostichopus japonicus by transcriptome analysis

The sea cucumber Apostichopus japonicus is a flourishing aquaculture species in China. However, there are challenges for sea cucumber aquaculture, one of which is the high temperature in summer. In this study, we explored the transcriptome expression profiles with seasons (APR, JUN and JUL) in the muscle tissue of A. japonicus. The temperature of the natural coast was 13 °C, 21 °C and 25 °C respectively when sampling. Compared with APR group, changes of expression profiles were more significant in JUL group than that in JUN group. A total of 46 differential expressed genes (DEGs) involved in both innate and adaptive immunity were highlighted, including 27 up-regulated and 19 down-regulated genes. They were further grouped into 10 sub-classes: heat shock, coagulation cascades, antigen processing and presentation, inflammatory response, transporter activity, immunoglobulin, lectin C, cell adhesion, reactive oxygen species (ROS) scavenging, apoptosis and autophagy. The study will offer deep insights of the molecular mechanisms underlying the physiological responses to seasonal high temperature in A. japonicus. Particularly, knowledge about the immunological effects of seasonal temperature on the species is critical for the optimal management practices for both wild and aquaculture populations.

Genome Sequencing of the Japanese Eel (Anguilla japonica) for Comparative Genomic Studies on tbx4 and a tbx4 Gene Cluster in Teleost Fishes

Limbs originated from paired fish fins are an important innovation in Gnathostomata. Many studies have focused on limb development-related genes, of which the T-box transcription factor 4 gene (tbx4) has been considered as one of the most essential factors in the regulation of the hindlimb development. We previously confirmed pelvic fin loss in tbx4-knockout zebrafish. Here, we report a high-quality genome assembly of the Japanese eel (Anguilla japonica), which is an economically important fish without pelvic fins. The assembled genome is 1.13 Gb in size, with a scaffold N50 of 1.03 Mb. In addition, we collected 24 tbx4 sequences from 22 teleost fishes to explore the correlation between tbx4 and pelvic fin evolution. However, we observed complete exon structures of tbx4 in several pelvic-fin-loss species such as Ocean sunfish (Mola mola) and ricefield eel (Monopterus albus). More interestingly, an inversion of a special tbx4 gene cluster (brip1-tbx4-tbx2b- bcas3) occurred twice independently, which coincides with the presence of fin spines. A nonsynonymous mutation (M82L) was identified in the nuclear localization sequence (NLS) of the Japanese eel tbx4. We also examined variation and loss of hindlimb enhancer B (HLEB), which may account for pelvic fin loss in Tetraodontidae and Diodontidae. In summary, we generated a genome assembly of the Japanese eel, which provides a valuable genomic resource to study the evolution of fish tbx4 and helps elucidate the mechanism of pelvic fin loss in teleost fishes. Our comparative genomic studies, revealed for the first time a potential correlation between the tbx4 gene cluster and the evolutionary development of toxic fin spines. Because fin spines in teleosts are usually venoms, this tbx4 gene cluster may facilitate the genetic engineering of toxin-related marine drugs.

The Axial Organ and the Pharynx Are Sites of Hematopoiesis in the Sea Urchin

Background: The location of coelomocyte proliferation in adult sea urchins is unknown and speculations since the early 1800s have been based on microanatomy and tracer uptake studies. In adult sea urchins (Strongylocentrotus purpuratus) with down-regulated immune systems, coelomocyte numbers increase in response to immune challenge, and whether some or all of these cells are newly proliferated is not known. The gene regulatory network that encodes transcription factors that control hematopoiesis in embryonic and larval sea urchins has not been investigated in adults. Hence, to identify the hematopoietic tissue in adult sea urchins, cell proliferation, expression of phagocyte specific genes, and expression of genes encoding transcription factors that function in the conserved regulatory network that controls hematopoiesis in embryonic and larval sea urchins were investigated for several tissues. Results: Cell proliferation was induced in adult sea urchins either by immune challenge through injection of heat-killed Vibrio diazotrophicus or by cell depletion through aspiration of coelomic fluid. In response to either of these stimuli, newly proliferated coelomocytes constitute only about 10% of the cells in the coelomic fluid. In tissues, newly proliferated cells and cells that express SpTransformer proteins (formerly Sp185/333) that are markers for phagocytes are present in the axial organ, gonad, pharynx, esophagus, and gut with no differences among tissues. The expression level of genes encoding transcription factors that regulate hematopoiesis show that both the axial organ and the pharynx have elevated expression compared to coelomocytes, esophagus, gut, and gonad. Similarly, an RNAseq dataset shows similar results for the axial organ and pharynx, but also suggests that the axial organ may be a site for removal and recycling of cells in the coelomic cavity. Conclusions: Results presented here are consistent with previous speculations that the axial organ may be a site of coelomocyte proliferation and that it may also be a center for cellular removal and recycling. A second site, the pharynx, may also have hematopoietic activity, a tissue that has been assumed to function only as part of the intestinal tract.

The Natural History of Teneurins: A Billion Years of Evolution in Three Key Steps

The entire evolutionary history of the animal gene family, Teneurin, can be summed up in three key steps, plus three salient footnotes. In a shared ancestor of all bilaterians, the first step began with gene fusions that created a protein with an amino-terminal intracellular domain bridged via a single transmembrane helix to extracellular EGF-like domains. This first step was completed with a further gene fusion: an additional carboxy-terminal stretch of about 2000 amino acids (aa) was adopted, as-a-whole, from bacteria. The 2000 aa structure in Teneurin was recently solved in three dimensions. The 2000 aa region appears in a number of bacteria, yet was co-opted solely into Teneurin, and into no other eukaryotic proteins. Outside of bilaterian animals, no Teneurins exist, with a “Monosiga brevicollis caveat” brought below, as ‘the third footnote.” Subsequent to the “urTeneurin’s” genesis-by-fusions, all bilaterians bore a single Teneurin gene, always encoding an extraordinarily conserved Type II transmembrane protein with invariant domain content and order. The second key step was a duplication that led to an exception to singleton Teneurin genomes. A pair of Teneurin paralogs, Ten-a and Ten-m, are found in representatives of all four Arthropod sub-phyla, in: insects, crustaceans, myriapods, and chelicerates. In contrast, in every other protostome species’ genome, including those of all non-Arthropod ecdysozoan phyla, only a single Teneurin gene occurs. The closest, sister, phylum of arthropods, the Onychophorans (velvet worms), bear a singleton Teneurin. Ten-a and Ten-m therefore arose from a duplication in an urArthropod only after Arthropods split from Onychophorans, but before the splits that led to the four Arthropod sub-phyla. The third key step was a quadruplication of Teneurins at the root of vertebrate radiation. Four Teneurin paralogs (Teneurins 1 through 4) arose first by a duplication of a single chordate gene likely leading to one 1/4–type gene, and one 2/3-type gene: the two copies found in extant jawless vertebrates. Relatively soon thereafter, a second duplication round yielded the -1, -2, -3, and -4 paralog types now found in all jawed vertebrates, from sharks to humans. It is possible to assert that these duplication events correlate well to the Ohno hypothesized 2R (two round) vertebrate whole genome duplication (WGD), as refined in more recent treatments. The quadruplication can therefore be placed at approximately 400 Myr ago. Echinoderms, hemichordates, cephalochordates, and urochordates have only a single copy of Teneurin in their genomes. These deuterostomes and non-vertebrate chordates provide the anchor showing that the quadruplication happened at the root of vertebrates. A first footnote must be brought concerning some of the ‘invertebrate’ relatives of vertebrates, among Deuterostomes. A family of genes that encode 7000 aa proteins was derived from, but is distinct from, the Teneurin family. This distinct family arose early in deuterostomes, yet persists today only in hemichordate and cephalochordate genomes. They are named here TRIPs (Teneurin-related immense proteins). As a second of three ‘footnotes’: a limited number of species exist with additional Teneurin gene copies. However, these further duplications of Teneurins occur for paralog types (a, m, or 1–4) only in specific lineages within Arthropods or Vertebrates. All examples are paralog duplications that evidently arose in association with lineage specific WGDs. The increased Teneurin paralog numbers correlate with WGDs known and published in bony fish, Xenopus, plus select Chelicerates lineages and Crustaceans. The third footnote, alluded to above, is that a Teneurin occurs in one unicellular species: Monosiga brevicollis. Teneurins are solely a metazoan, bilaterian-specific family, to the exclusion of the Kingdoms of prokaryotes, plants, fungi, and protists. The single exception occurs among the unicellular, opisthokont, closest relatives of metazoans, the choanoflagellates. There is a Teneurin in Monosiga brevicollis, one species of the two fully sequenced choanoflagellate species. In contrast, outside of triploblast-bilaterians, there are no Teneurins in any diploblast genomes, including even sponges – those metazoans closest to choanoflagellates. Perhaps the ‘birth’ of the original Teneurin occurred in a shared ancestor of M. brevicollis and metazoans, then was lost in M. brevicollis’ sister species, and was serially and repeatedly lost in all diploblast metazoans. Alternatively, and as favored above, it first arose in the ‘urBilaterian,’ then was subsequently acquired from some bilaterian via horizontal transfer by a single choanoflagellate clade. The functional partnership of Teneurins and Latrophilins was discovered in rodents through the LPH1-TENM2 interaction. Recent work extends this to further members of each family. Surveying when the interacting domains of Teneurins and Latrophilins co-exist within different organisms can give an indication of how widespread their functional cooperation might be, across bilaterians. Paralog number for the two families is relatively correlated among bilaterians, and paralog numbers underwent co-increase in the WGDs mentioned above. With co-increasing paralog numbers, the possible combinatorial pairs grow factorially. This should have a significant impact for increasing nervous system complexity. The 3 key events in the ‘natural history’ of the Teneurins and their Latrophilin partners coincide with the ascendance of particularly successful metazoan clades: bilaterians; arthropods; and vertebrates. Perhaps we can attribute some of this success to the unique Teneurin family, and to its partnership with Latrophilins.

Origin and evolution of the specialized forms of proteasomes involved in antigen presentation

Proteasomes are a multi-subunit protease complex that produces peptides bound by major histocompatibility complex (MHC) class I molecules. Phylogenetic studies indicate that two specialized forms of proteasomes, immunoproteasomes and thymoproteasomes, and the proteasome activator PA28αβ emerged in a common ancestor of jawed vertebrates which acquired adaptive immunity based on the MHC, T cell receptors, and B cell receptors ~ 500 million years ago. Comparative genomics studies now provide strong evidence that the genes coding for the immunoproteasome subunits emerged by genome-wide duplication. On the other hand, the gene encoding the thymoproteasome subunit β5t emerged by tandem duplication from the gene coding for the β5 subunit. Strikingly, birds lack immunoproteasomes, thymoproteasomes, and the proteasome activator PA28αβ, raising an interesting question of whether they have evolved any compensatory mechanisms.

Transcriptomic analysis of immunoglobulin novel antigen receptor (IgNAR) heavy chain constant domains of brownbanded bamboo shark (Chiloscyllium punctatum)

Cartilaginous fish are the evolutionarily oldest group of animals which possess antibodies, T cell receptors and major histocompatibility complex (MHC). The immunoglobulin novel antigen receptor (IgNAR) found in cartilaginous fish is a heavy chain homodimer which lacks light chain. The presence of non-canonical cysteine molecules and lack of CDR2 region make it more significant. To synthesize active binding domains based on variable region of IgNAR (VNAR), knowledge on the constant region dynamics play a significant role. The IgNAR exhibit species variations in its primary sequence features; hence, this study was conducted to determine the IgNAR heavy chain constant domain of the brownbanded bamboo shark (Chiloscyllium punctatum). Peripheral blood leukocytes (PBL) isolated from adult bamboo sharks were used to synthesize a cDNA library. A total of four billion residues of two million sequences (average length 218.41 bp) were obtained. Assembled sequences were aligned with published cartilaginous fish IgNAR constant region sequences. Transcriptome analysis revealed two distinct types of IgNAR in the brownbanded bamboo shark. Also, constant-1 domain sequences displayed 13 unique sequences which may reflect the least number of IgNAR gene clusters. The phylogenetic analysis revealed the closest relationship with the nurse shark (Ginglymostoma cirratum) followed by the wobbegong shark (Orectolobus maculatus) which belong to the same order Orectolobiformes. Analysis of the constant domains of the brownbanded bamboo shark IgNAR revealed an evolutionarily conserved nature and this knowledge can be used to design primers for VNAR cloning. Furthermore, knowledge on the structural features in IgNAR constant domains that increase the stability could be useful in the process of stabilizing human immunoglobulins.

Ambulacrarians and the Ancestry of Deuterostome Nervous Systems

The evolutionary origin and history of metazoan nervous systems has been at the heart of numerous scientific debates for well over a century. This has been a particularly difficult issue to resolve within the deuterostomes, chiefly due to the distinct neural architectures observed within this group of animals. Indeed, deuterosomes feature central nervous systems, apical organs, nerve cords, and basiepidermal nerve nets. Comparative analyses investigating the anatomy and molecular composition of deuterostome nervous systems have nonetheless succeeded in identifying a number of shared and derived features. These analyses have led to the elaboration of diverse theories about the origin and evolutionary history of deuterostome nervous systems. Here, we provide an overview of these distinct theories. Further, we argue that deciphering the adult nervous systems of representatives of all deuterostome phyla, including echinoderms, which have long been neglected in this type of surveys, will ultimately provide answers to the questions concerning the ancestry and evolution of deuterostome nervous systems.

Identification of a novel antimicrobial peptide from the sea star Patiria pectinifera

Antimicrobial peptides (AMPs) are components of innate immunity found in many forms of life. However, there have been no reports of AMPs in sea star (Phylum Echinodermata). Here we report the isolation and characterization of a novel antimicrobial peptide from the coelomic epithelium extract of the sea star Patiria pectinifera. The isolated peptide comprises 38 amino acid residues, is cationic (pI 9.2), has four cysteine residues that form two disulfide bonds (C1-C3 and C2-C4), is amidated at the C-terminus, and is designated P. pectinifera cysteine-rich antimicrobial peptide (PpCrAMP). Synthetic PpCrAMP identical to the native peptide exhibited the most potent antimicrobial activity compared to analogs with different disulfide bond configurations. Expression analysis of PpCrAMP precursor transcripts revealed constitutive expression in the coelomic epithelium and tube feet of P. pectinifera. Analysis of genomic DNA and cDNA encoding the PpCrAMP precursor protein revealed that an intron splits the coding region of the mature peptide into a positively charged N-terminal domain and a C-terminal domain harboring four cysteine residues and a glycine for C-terminal amidation. No significant homology with other known AMPs was observed, while orthologs of PpCrAMP were found in other echinoderm species. These findings indicate that PpCrAMP is the prototype of a family a novel cysteine-rich AMPs that participate in mechanisms of innate immunity in echinoderms. Furthermore, the discovery of PpCrAMP may lead to the identification of related AMPs in vertebrates and protostome invertebrates.

Whole Genome Sequencing of the Pirarucu (Arapaima gigas) Supports Independent Emergence of Major Teleost Clades

The Pirarucu (Arapaima gigas) is one of the world's largest freshwater fishes and member of the superorder Osteoglossomorpha (bonytongues), one of the oldest lineages of ray-finned fishes. This species is an obligate air-breather found in the basin of the Amazon River with an attractive potential for aquaculture. Its phylogenetic position among bony fishes makes the Pirarucu a relevant subject for evolutionary studies of early teleost diversification. Here, we present, for the first time, a draft genome version of the A. gigas genome, providing useful information for further functional and evolutionary studies. The A. gigas genome was assembled with 103-Gb raw reads sequenced in an Illumina platform. The final draft genome assembly was ∼661 Mb, with a contig N50 equal to 51.23 kb and scaffold N50 of 668 kb. Repeat sequences accounted for 21.69% of the whole genome, and a total of 24,655 protein-coding genes were predicted from the genome assembly, with an average of nine exons per gene. Phylogenomic analysis based on 24 fish species supported the postulation that Osteoglossomorpha and Elopomorpha (eels, tarpons, and bonefishes) are sister groups, both forming a sister lineage with respect to Clupeocephala (remaining teleosts). Divergence time estimations suggested that Osteoglossomorpha and Elopomorpha lineages emerged independently in a period of ∼30 Myr in the Jurassic. The draft genome of A. gigas provides a valuable genetic resource for further investigations of evolutionary studies and may also offer a valuable data for economic applications.

The immunoproteasome and thymoproteasome: functions, evolution and human disease

The basic principle of adaptive immunity is to strictly discriminate between self and non-self, and a central challenge to overcome is the enormous variety of pathogens that might be encountered. In cell-mediated immunity, immunological discernment takes place at a molecular or cellular level. Central to both mechanisms of discernment is the generation of antigenic peptides associated with MHC class I molecules, which is achieved by a proteolytic complex called the proteasome. To adequately accomplish the discrimination between self and non-self that is essential for adaptive immunity and self-tolerance, two proteasome subtypes have evolved via gene duplication: the immunoproteasome and the thymoproteasome. In this Review, we describe various aspects of these immunity-dedicated proteasomes, from their discovery to recent findings.

Morphological Stasis and Proteome Innovation in Cephalochordates

Lancelets, extant representatives of basal chordates, are prototypic examples of evolutionary stasis; they preserved a morphology and body-plan most similar to the fossil chordates from the early Cambrian. Such a low level of morphological evolution is in harmony with a low rate of amino acid substitution; cephalochordate proteins were shown to evolve slower than those of the slowest evolving vertebrate, the elephant shark. Surprisingly, a study comparing the predicted proteomes of Chinese amphioxus, Branchiostoma belcheri and the Florida amphioxus, Branchiostoma floridae has led to the conclusion that the rate of creation of novel domain combinations is orders of magnitude greater in lancelets than in any other Metazoa, a finding that contradicts the notion that high rates of protein innovation are usually associated with major evolutionary innovations. Our earlier studies on a representative sample of proteins have provided evidence suggesting that the differences in the domain architectures of predicted proteins of these two lancelet species reflect annotation errors, rather than true innovations. In the present work, we have extended these studies to include a larger sample of genes and two additional lancelet species, Asymmetron lucayanum and Branchiostoma lanceolatum. These analyses have confirmed that the domain architecture differences of orthologous proteins of the four lancelet species are because of errors of gene prediction, the error rate in the given species being inversely related to the quality of the transcriptome dataset that was used to aid gene prediction.

Genome-wide organization, evolutionary diversification of the COMMD family genes of amphioxus (Branchiostoma belcheri) with the possible role in innate immunity

The COMMD (COpper Metabolism gene MURR1 Domain) gene family with ten members participates in various biological processes, such as the regulation of copper and sodium transport, NF-κB activity and cell cycle progression. However, studies on the COMMD gene family in amphioxus (Branchiostoma belcheri) are yet largely unknown. In this study, we have identified and characterized the ten COMMD family members from amphioxus (designated as AmphiCOMMDs). Firstly, we clone the full length of AmphiCOMMDs, and all AmphiCOMMD proteins contain the conserved COMM domain with two NES (Nuclear Export Signal) motifs. Secondly, the genomic structure analysis demonstrates that genes of the COMMD family have undergone intron loss and gain during the process of divergence from amphioxus to vertebrates. Thirdly, phylogenetic analysis indicates that AmphiCOMMDs are more closely related to vertebrates, implying the AmphiCOMMDs may be the ancestor of the vertebrate COMMDs. Fourthly, AmphiCOMMDs are ubiquitously and differentially expressed in five investigated tissues (muscles, gills, intestine, heaptic cecum and notochord). Finally, our results show that expression levels of AmphiCOMMD genes are fluctuating after LPS stimulation to some different extent. Taken together, our studies have elaborated the evolutionary dynamic and the innate immune role of the COMMD family genes in amphioxus.