Viruses of lower vertebrates recently became a field of interest to the public due to increasing epizootics and economic losses of poikilothermic animals. These were reported worldwide from both wildlife and collections of aquatic poikilothermic animals. Several RNA and DNA viruses infecting fish, amphibians and reptiles have been studied intensively during the last 20 years. Many of these viruses induce diseases resulting in important economic losses of lower vertebrates, especially in fish aquaculture. In addition, some of the DNA viruses seem to be emerging pathogens involved in the worldwide decline in wildlife. Irido-, herpes- and polyomavirus infections may be involved in the reduction in the numbers of endangered amphibian and reptile species. In this context the knowledge of several important RNA viruses such as orthomyxo-, paramyxo-, rhabdo-, retro-, corona-, calici-, toga-, picorna-, noda-, reo- and birnaviruses, and DNA viruses such as parvo-, irido-, herpes-, adeno-, polyoma- and poxviruses, is described in this review.

The matrix metalloproteinase (MMP) family belongs to the metzincin clan of zinc-dependent metallopeptidases. Due to their enormous implications in physiology and disease, MMPs have mainly been studied in vertebrates. They are engaged in extracellular protein processing and degradation, and present extensive paralogy, with 23 forms in humans. One characteristic of MMPs is a ~165-residue catalytic domain (CD), which has been structurally studied for 14 MMPs from human, mouse, rat, pig and the oral-microbiome bacterium Tannerella forsythia. These studies revealed close overall coincidence and characteristic structural features, which distinguish MMPs from other metzincins and give rise to a sequence pattern for their identification. Here, we reviewed the literature available on MMPs outside vertebrates and performed database searches for potential MMP CDs in invertebrates, plants, fungi, viruses, protists, archaea and bacteria. These and previous results revealed that MMPs are widely present in several copies in Eumetazoa and higher plants (Tracheophyta), but have just token presence in eukaryotic algae. A few dozen sequences were found in Ascomycota (within fungi) and in double-stranded DNA viruses infecting invertebrates (within viruses). In contrast, a few hundred sequences were found in archaea and >1000 in bacteria, with several copies for some species. Most of the archaeal and bacterial phyla containing potential MMPs are present in human oral and gut microbiomes. Overall, MMP-like sequences are present across all kingdoms of life, but their asymmetric distribution contradicts the vertical descent model from a eubacterial or archaeal ancestor. This article is part of a Special Issue entitled: Matrix Metalloproteinases edited by Rafael Fridman.

The pituitary is a master endocrine gland that developed early in vertebrate evolution and therefore exists in all modern vertebrate classes. The last decade has transformed our view of this key organ. Traditionally, the pituitary has been viewed as a randomly organized collection of cells that respond to hypothalamic stimuli by secreting their content. However, recent studies have established that pituitary cells are organized in tightly wired large-scale networks that communicate with each other in both homo and heterotypic manners, allowing the gland to quickly adapt to changing physiological demands. These networks functionally decode and integrate the hypothalamic and systemic stimuli and serve to optimize the pituitary output into the generation of physiologically meaningful hormone pulses. The development of 3D imaging methods and transgenic models have allowed us to expand the research of functional pituitary networks into several vertebrate classes. Here we review the establishment of pituitary cell networks throughout vertebrate evolution and highlight the main perspectives and future directions needed to decipher the way by which pituitary networks serve to generate hormone pulses in vertebrates.

The central roles of luteinizing hormone (LH), progestin and their receptors for initiating ovulation have been well established. However, signaling pathways and downstream targets such as proteases that are essential for the rupture of follicular cells are still unclear. Recently, we found anovulation in nuclear progestin receptor (Pgr) knockout (Pgr-KO) zebrafish, which offers a new model for examining genes and pathways that are important for ovulation and fertility. In this study, we examined expression of all transcripts using RNA-Seq in preovulatory follicular cells collected following the final oocyte maturation, but prior to ovulation, from wild-type (WT) or Pgr-KO fish. Differential expression analysis revealed 3567 genes significantly differentially expressed between WT and Pgr-KO fish (fold change⩾2, p<0.05). Among those, 1543 gene transcripts were significantly more expressed, while 2024 genes were significantly less expressed, in WT than those in Pgr-KO. We then retrieved and compared transcriptional data from online databases and further identified 661 conserved genes in fish, mice, and humans that showed similar levels of high (283 genes) or low (387) expression in animals that were ovulating compared to those with no ovulation. For the first time, ovulatory genes and their involved biological processes and pathways were also visualized using Enrichment Map and Cytoscape. Intriguingly, enrichment analysis indicated that the genes with higher expression were involved in multiple ovulatory pathways and processes such as inflammatory response, angiogenesis, cytokine production, cell migration, chemotaxis, MAPK, focal adhesion, and cytoskeleton reorganization. In contrast, the genes with lower expression were mainly involved in DNA replication, DNA repair, DNA methylation, RNA processing, telomere maintenance, spindle assembling, nuclear acid transport, catabolic processes, and nuclear and cell division. Our results indicate that a large set of genes (>3000) is differentially regulated in the follicular cells in zebrafish prior to ovulation, terminating programs such as growth and proliferation, and beginning processes including the inflammatory response and apoptosis. Further studies are required to establish relationships among these genes and an ovulatory circuit in the zebrafish model.

The BRCA2 interactor, centrobin, is a centrosomal protein that has been implicated in centriole duplication and microtubule stability. We used genome editing to ablate CNTROB in hTERT-RPE1 cells and observed an increased frequency of monocentriolar and acentriolar cells. Using a novel monoclonal antibody, we found that centrobin primarily localizes to daughter centrioles but also associates with mother centrioles upon serum starvation. Strikingly, centrobin loss abrogated primary ciliation upon serum starvation. Ultrastructural analysis of centrobin nulls revealed defective axonemal extension after mother centriole docking. Ciliogenesis required a C-terminal portion of centrobin that interacts with CP110 and tubulin. We also depleted centrobin in zebrafish embryos to explore its roles in an entire organism. Centrobin-depleted embryos showed microcephaly, with curved and shorter bodies, along with marked defects in laterality control, morphological features that indicate ciliary dysfunction. Our data identify new roles for centrobin as a positive regulator of vertebrate ciliogenesis.

Forkhead box, class O (FoxO) belongs to the large family of forkhead transcription factors that are characterized by a conserved forkhead box DNA-binding domain. To date, the FoxO group has four mammalian members: FoxO1, FoxO3a, FoxO4 and FoxO6, which are orthologs of DAF16, an insulin-responsive transcription factor involved in regulating longevity of worms and flies. The degree of homology between these four members is high, especially in the forkhead domain, which contains the DNA-binding interface. Yet, mouse FoxO knockouts have revealed that each FoxO gene has its unique role in the physiological process. Whether the functional divergences are primarily due to adaptive selection pressure or relaxed selective constraint remains an open question. As such, this study aims to address the evolutionary mode of FoxO, which may lead to the functional divergence.

Amelogenesis imperfecta is a group of disorders causing abnormalities in enamel formation in various phenotypes. Many mutations in the FAM83H gene have been identified to result in autosomal dominant hypocalcified amelogenesis imperfecta in different populations. However, the structure and function of FAM83H and its pathological mechanism have yet to be further explored. Evolutionary analysis is an alternative for revealing residues or motifs that are important for protein function. In the present study, we chose 50 vertebrate species in public databases representative of approximately 230 million years of evolution, including 1 amphibian, 2 fishes, 7 sauropsidas and 40 mammals, and we performed evolutionary analysis on the FAM83H protein. By sequence alignment, conserved residues and motifs were indicated, and the loss of important residues and motifs of five special species (Malayan pangolin, platypus, minke whale, nine-banded armadillo and aardvark) was discovered. A phylogenetic time tree showed the FAM83H divergent process. Positive selection sites in the C-terminus suggested that the C-terminus of FAM83H played certain adaptive roles during evolution. The results confirmed some important motifs reported in previous findings and identified some new highly conserved residues and motifs that need further investigation. The results suggest that the C-terminus of FAM83H contain key conserved regions critical to enamel formation and calcification.

The growth hormone secretagogue-receptor (GHS-R) was discovered in humans and pigs in 1996. The endogenous ligand, ghrelin, was discovered 3 years later, in 1999, and our understanding of the physiological significance of the ghrelin system in vertebrates has grown steadily since then. Although the ghrelin system in non-mammalian vertebrates is a subject of great interest, protein sequence data for the receptor in non-mammalian vertebrates has been limited until recently, and related biological information has not been well organized. In this review, we summarize current information related to the ghrelin receptor in non-mammalian vertebrates.

This study investigated long-term substitution rate differences using three calibration points, divergences between lobe-finned vertebrates and ray-finned fish, between mammals and sauropsids, and between holosteans (gar and bowfin) and teleost fish with amino acid sequence data of 625 genes for 25 bony vertebrates. The result showed that the substitution rate was two to three times higher in the stem branches of lobe-finned vertebrates before the mammal-sauropsid divergence than in amniotes. The rate in the stem branch of ray-finned fish before the holostean-teleost fish divergence was also a few times higher than the holostean rate, whereas it was similar to or somewhat slower than the teleost fish rate. The phylogenetic relationship of coelacanth and lungfish with tetrapod was difficult to determine because of the short interval of the divergences. Considering the high rate in the stem branches, the divergences of coelacanth and lungfish from the stem branch were estimated as 408-427 Ma and 399-414 Ma, respectively, with the interval of 9-13 Myr. With the external calibration of the mammal-sauropsid split, the estimated times for ordinal divergences within eutherian mammals tend to be smaller than those in previous studies that used the calibration points within the lineage, with deeper divergences before the Cretaceous-Paleogene boundary and shallower ones after the boundary. In contrast the estimated times within birds were larger than those of previous studies, with the divergence between Galliformes and Anseriformes ∼80 Ma and that between Galloanserae and Neoaves 110 Ma.

Generation of mature red blood cells, consisting mainly of hemoglobin, is a remarkable example of coordinated action of various signaling networks. Chromatin condensation is an essential step for terminal erythroid differentiation and subsequent nuclear expulsion in mammals. Here, we profiled 3D genome organization in the blood cells from ten species belonging to different vertebrate classes. Our analysis of contact maps revealed a striking absence of such 3D interaction patterns as loops or TADs in blood cells of all analyzed representatives. We also detect large-scale chromatin rearrangements in blood cells from mammals, birds, reptiles and amphibians: their contact maps display strong second diagonal pattern, representing an increased frequency of long-range contacts, unrelated to TADs or compartments. This pattern is completely atypical for interphase chromosome structure. We confirm that these principles of genome organization are conservative in vertebrate erythroid cells.

All biological functions in vertebrates are synchronized with daily and seasonal changes in the environment by the time keeping hormone melatonin. Its nocturnal surge is primarily due to the rhythmic activity of the arylalkylamine N-acetyl transferase AANAT, which thus became the focus of many investigations regarding its evolution and function. Various vertebrate isoforms have been reported from cartilaginous fish to mammals but their origin has not been clearly established. Using phylogeny and synteny, we took advantage of the increasing number of available genomes in order to test whether the various rounds of vertebrate whole genome duplications were responsible for the diversification of AANAT. We highlight a gene secondary loss of the AANAT2 in the Sarcopterygii, revealing for the first time that the AAANAT1/2 duplication occurred before the divergence between Actinopterygii (bony fish) and Sarcopterygii (tetrapods, lobe-finned fish, and lungfish). We hypothesize the teleost-specific whole genome duplication (WDG) generated the appearance of the AANAT1a/1b and the AANAT2/2'paralogs, the 2' isoform being rapidly lost in the teleost common ancestor (ray-finned fish). We also demonstrate the secondary loss of the AANAT1a in a Paracantopterygii (Atlantic cod) and of the 1b in some Ostariophysi (zebrafish and cave fish). Salmonids present an even more diverse set of AANATs that may be due to their specific WGD followed by secondary losses. We propose that vertebrate AANAT diversity resulted from 3 rounds of WGD followed by previously uncharacterized secondary losses. Extant isoforms show subfunctionalized localizations, enzyme activities and affinities that have increased with time since their emergence.

Hemostasis is a defense mechanism which protects the organism in the event of injury to stop bleeding. Recently, we established that all the known major mammalian hemostatic factors are conserved in early vertebrates. However, since their highly vascularized gills experience high blood pressure and are exposed to the environment, even very small injuries could be fatal to fish. Since trypsins are forerunners for coagulation proteases and are expressed by many extrapancreatic cells such as endothelial cells and epithelial cells, we hypothesized that trypsin or trypsin-like proteases from gill epithelial cells may protect these animals from gill bleeding following injuries. In this paper we identified the release of three different trypsins from fish gills into water under stress or injury, which have tenfold greater serine protease activity compared to bovine trypsin. We found that these trypsins activate the thrombocytes and protect the fish from gill bleeding. We found 27 protease-activated receptors (PARs) by analyzing zebrafish genome and classified them into five groups, based on tethering peptides, and two families, PAR1 and PAR2, based on homologies. We also found a canonical member of PAR2 family, PAR2-21A which is activated more readily by trypsin, and PAR2-21A tethering peptide stops gill bleeding just as trypsin. This finding provides evidence that trypsin cleaves a PAR2 member on thrombocyte surface. In conclusion, we believe that the gills are evolutionarily selected to produce trypsin to activate PAR2 on thrombocyte surface and protect the gills from bleeding. We also speculate that trypsin may also protect the fish from bleeding from other body injuries due to quick contact with the thrombocytes. Thus, this finding provides evidence for the role of trypsins in primary hemostasis in early vertebrates.

The hypothalamus coordinates neuroendocrine functions in vertebrates. To explore its evolutionary origin, we describe integrated transcriptome/connectome brain maps for swimming tadpoles of Ciona, which serves as an approximation of the ancestral proto-vertebrate. This map features several cell types related to different regions of the vertebrate hypothalamus, including the mammillary nucleus, the arcuate nucleus, and magnocellular neurons. Coronet cells express melanopsin and share additional properties with the saccus vasculosus, a specialized region of the hypothalamus that mediates photoperiodism in nontropical fishes. Comparative transcriptome analyses identified orthologous cell types for mechanosensory switch neurons, and VP+ and VPR+ relay neurons in different regions of the mouse hypothalamus. These observations provide evidence that the hypothalamus predates the evolution of the vertebrate brain. We discuss the possibility that switch neurons, coronet cells, and FoxP+ /VPR+ relay neurons comprise a behavioral circuit that helps trigger metamorphosis of Ciona larvae in response to twilight.

Evolutionary innovation relies partially on changes in gene regulation. While a growing body of evidence demonstrates that such innovation is generated by functional changes or translocation of regulatory elements via mobile genetic elements, the de novo generation of enhancers from non-regulatory/non-mobile sequences has, to our knowledge, not previously been demonstrated. Here we show evidence for the de novo genesis of enhancers in vertebrates. For this, we took advantage of the massive gene loss following the last whole genome duplication in teleosts to systematically identify regions that have lost their coding capacity but retain sequence conservation with mammals. We found that these regions show enhancer activity while the orthologous coding regions have no regulatory activity. These results demonstrate that these enhancers have been de novo generated in fish. By revealing that minor changes in non-regulatory sequences are sufficient to generate new enhancers, our study highlights an important playground for creating new regulatory variability and evolutionary innovation.

MicroRNAs (miRNAs) are a class of non-coding RNA that plays an important role in posttranscriptional regulation of mRNA. Evidence has shown that miRNA gene variability might interfere with its function resulting in phenotypic variation and disease susceptibility. A major role in miRNA target recognition is ascribed to complementarity with the miRNA seed region that can be affected by polymorphisms. In the present study, we developed an online tool for the detection of miRNA polymorphisms (miRNA SNiPer) in vertebrates (http://www.integratomics-time.com/miRNA-SNiPer) and generated a catalog of miRNA seed region polymorphisms (miR-seed-SNPs) consisting of 149 SNPs in six species. Although a majority of detected polymorphisms were due to point mutations, two consecutive nucleotide substitutions (double nucleotide polymorphisms, DNPs) were also identified in nine miRNAs. We determined that miR-SNPs are frequently located within the quantitative trait loci (QTL), chromosome fragile sites, and cancer susceptibility loci, indicating their potential role in the genetic control of various complex traits. To test this further, we performed an association analysis between the mmu-miR-717 seed SNP rs30372501, which is polymorphic in a large number of standard inbred strains, and all phenotypic traits in these strains deposited in the Mouse Phenome Database. Analysis showed a significant association between the mmu-miR-717 seed SNP and a diverse array of traits including behavior, blood-clinical chemistry, body weight size and growth, and immune system suggesting that seed SNPs can indeed have major pleiotropic effects. The bioinformatics analyses, data and tools developed in the present study can serve researchers as a starting point in testing more targeted hypotheses and designing experiments using optimal species or strains for further mechanistic studies.

Gamma-aminobutyric acid (GABA) is a widely conserved signaling molecule that in animals has been adapted as a neurotransmitter. GABA is synthesized from the amino acid glutamate by the action of glutamate decarboxylases (GADs). Two vertebrate genes, GAD1 and GAD2, encode distinct GAD proteins: GAD67 and GAD65, respectively. We have identified a third vertebrate GAD gene, GAD3. This gene is conserved in fishes as well as tetrapods. We analyzed protein sequence, gene structure, synteny, and phylogenetics to identify GAD3 as a homolog of GAD1 and GAD2. Interestingly, we found that GAD3 was lost in the hominid lineage. Because of the importance of GABA as a neurotransmitter, GAD3 may play important roles in vertebrate nervous systems.

Acoustic communication is crucial to humans and many other tetrapods, including birds, frogs, crocodilians, and mammals. However, large-scale patterns in its evolution are largely unstudied. Here, we address several fundamental questions about the origins of acoustic communication in terrestrial vertebrates (tetrapods), using phylogenetic methods. We show that origins of acoustic communication are significantly associated with nocturnal activity. We find that acoustic communication does not increase diversification rates, a surprising result given the many speciation-focused studies of frog calls and bird songs. We also demonstrate that the presence of acoustic communication is strongly conserved over time. Finally, we find that acoustic communication evolved independently in most major tetrapod groups, often with remarkably ancient origins (~100-200 million years ago). Overall, we show that the role of ecology in shaping signal evolution applies to surprisingly deep timescales, whereas the role of signal evolution in diversification may not.

Gene expression profiles in homologous tissues have been observed to be different between species, which may be due to differences between species in the gene expression program in each cell type, but may also reflect differences in cell type composition of each tissue in different species. Here, we compare expression profiles in matching primary cells in human, mouse, rat, dog, and chicken using Cap Analysis Gene Expression (CAGE) and short RNA (sRNA) sequencing data from FANTOM5. While we find that expression profiles of orthologous genes in different species are highly correlated across cell types, in each cell type many genes were differentially expressed between species. Expression of genes with products involved in transcription, RNA processing, and transcriptional regulation was more likely to be conserved, while expression of genes encoding proteins involved in intercellular communication was more likely to have diverged during evolution. Conservation of expression correlated positively with the evolutionary age of genes, suggesting that divergence in expression levels of genes critical for cell function was restricted during evolution. Motif activity analysis showed that both promoters and enhancers are activated by the same transcription factors in different species. An analysis of expression levels of mature miRNAs and of primary miRNAs identified by CAGE revealed that evolutionary old miRNAs are more likely to have conserved expression patterns than young miRNAs. We conclude that key aspects of the regulatory network are conserved, while differential expression of genes involved in cell-to-cell communication may contribute greatly to phenotypic differences between species.

Vascular mural cells (vMCs) are essential components of the vertebrate vascular system, controlling blood vessel maturation and homeostasis. Discrete molecular mechanisms have been associated with vMC development and differentiation. The function of hemodynamic forces in controlling vMC recruitment is unclear. Using transgenic lines marking developing vMCs in zebrafish embryos, we find that vMCs are recruited by arterial-fated vessels and that the process is flow dependent. We take advantage of tissue-specific CRISPR gene targeting to demonstrate that hemodynamic-dependent Notch activation and the ensuing arterial genetic program is driven by endothelial primary cilia. We also identify zebrafish foxc1b as a cilia-dependent Notch-specific target that is required within endothelial cells to drive vMC recruitment. In summary, we have identified a hemodynamic-dependent mechanism in the developing vasculature that controls vMC recruitment.

The germline mutation rate determines the pace of genome evolution and is an evolving parameter itself1. However, little is known about what determines its evolution, as most studies of mutation rates have focused on single species with different methodologies2. Here we quantify germline mutation rates across vertebrates by sequencing and comparing the high-coverage genomes of 151 parent-offspring trios from 68 species of mammals, fishes, birds and reptiles. We show that the per-generation mutation rate varies among species by a factor of 40, with mutation rates being higher for males than for females in mammals and birds, but not in reptiles and fishes. The generation time, age at maturity and species-level fecundity are the key life-history traits affecting this variation among species. Furthermore, species with higher long-term effective population sizes tend to have lower mutation rates per generation, providing support for the drift barrier hypothesis3. The exceptionally high yearly mutation rates of domesticated animals, which have been continually selected on fecundity traits including shorter generation times, further support the importance of generation time in the evolution of mutation rates. Overall, our comparative analysis of pedigree-based mutation rates provides ecological insights on the mutation rate evolution in vertebrates.