Sequencing of the sea lamprey (Petromyzon marinus) genome provides insights into vertebrate evolution

Evolution of intrinsic disorder in eukaryotic proteins

Conformational flexibility conferred though regions of intrinsic structural disorder allows proteins to behave as dynamic molecules. While it is well-known that intrinsically disordered regions can undergo disorder-to-order transitions in real-time as part of their function, we also are beginning to learn more about the dynamics of disorder-to-order transitions along evolutionary time-scales. Intrinsically disordered regions endow proteins with functional promiscuity, which is further enhanced by the ability of some of these regions to undergo real-time disorder-to-order transitions. Disorder content affects gene retention after whole genome duplication, but it is not necessarily conserved. Altered patterns of disorder resulting from evolutionary disorder-to-order transitions indicate that disorder evolves to modify function through refining stability, regulation, and interactions. Here, we review the evolution of intrinsically disordered regions in eukaryotic proteins. We discuss the interplay between secondary structure and disorder on evolutionary time-scales, the importance of disorder for eukaryotic proteome expansion and functional divergence, and the evolutionary dynamics of disorder.

Evolution of the α2-adrenoreceptors in vertebrates: ADRA2D is absent in mammals and crocodiles

Evolutionary studies of genes that have been functionally characterized and whose variation has been associated with pathological conditions represent an opportunity to understand the genetic basis of pathologies. α₂-Adrenoreceptors (ADRA2) are a class of G protein-coupled receptors that regulate several physiological processes including blood pressure, platelet aggregation, insulin secretion, lipolysis, and neurotransmitter release. This gene family has been extensively studied from a molecular/physiological perspective, yet much less is known about its evolutionary history. Accordingly, the goal of this study was to investigate the evolutionary history of α₂-adrenoreceptors (ADRA2) in vertebrates. Our results show that in addition to the three well-recognized α₂-adrenoreceptor genes (ADRA2A, ADRA2B and ADRA2C), we recovered a clade that corresponds to the fourth member of the α₂-adrenoreceptor gene family (ADRA2D). We also recovered a clade that possesses two ADRA2 sequences found in two lamprey species. Furthermore, our results show that mammals and crocodiles are characterized by possessing three α₂-adrenoreceptor genes, whereas all other vertebrate groups possess the full repertoire of α₂-adrenoreceptor genes. Among vertebrates ADRA2D seems to be a dispensable gene, as it was lost two independent times during the evolutionary history of the group. Additionally, we found that most examined species possess the most common alleles described for humans; however, there are cases in which non-human mammals possess the alternative variant. Finally, transcript abundance profiles revealed that during the early evolutionary history of gnathostomes, the expression of ADRA2D in different taxonomic groups became specialized to different tissues, but in the ancestor of sarcopterygians this specialization would have been lost.

Deep ancestry of programmed genome rearrangement in lampreys

In most multicellular organisms, the structure and content of the genome is rigorously maintained over the course of development. However some species have evolved genome biologies that permit, or require, developmentally regulated changes in the physical structure and content of the genome (programmed genome rearrangement: PGR). Relatively few vertebrates are known to undergo PGR, although all agnathans surveyed to date (several hagfish and one lamprey: Petromyzon marinus) show evidence of large scale PGR. To further resolve the ancestry of PGR within vertebrates, we developed probes that allow simultaneous tracking of nearly all sequences eliminated by PGR in P. marinus and a second lamprey species (Entosphenus tridentatus). These comparative analyses reveal conserved subcellular structures (lagging chromatin and micronuclei) associated with PGR and provide the first comparative embryological evidence in support of the idea that PGR represents an ancient and evolutionarily stable strategy for regulating inherent developmental/genetic conflicts between germline and soma.

pBACode: a random-barcode-based high-throughput approach for BAC paired-end sequencing and physical clone mapping

Applications that use Bacterial Artificial Chromosome (BAC) libraries often require paired-end sequences and knowledge of the physical location of each clone in plates. To facilitate obtaining this information in high-throughput, we generated pBACode vectors: a pool of BAC cloning vectors, each with a pair of random barcodes flanking its cloning site. In a pBACode BAC library, the BAC ends and their linked barcodes can be sequenced in bulk. Barcode pairs are determined by sequencing the empty pBACode vectors, which allows BAC ends to be paired according to their barcodes. For physical clone mapping, the barcodes are used as unique markers for their linked genomic sequence. After multi-dimensional pooling of BAC clones, the barcodes are sequenced and deconvoluted to locate each clone. We generated a pBACode library of 94,464 clones for the flounder Paralichthys olivaceus and obtained paired-end sequence from 95.4% of the clones. Incorporating BAC paired-ends into the genome preassembly improved its continuity by over 10-fold. Furthermore, we were able to use the barcodes to map the physical locations of each clone in just 50 pools, with up to 11 808 clones per pool. Our physical clone mapping located 90.2% of BAC clones, enabling targeted characterization of chromosomal rearrangements.

Segmental arithmetic: summing up the Hox gene regulatory network for hindbrain development in chordates

Organization and development of the early vertebrate hindbrain are controlled by a cascade of regulatory interactions that govern the process of segmentation and patterning along the anterior-posterior axis via Hox genes. These interactions can be assembled into a gene regulatory network that provides a framework to interpret experimental data, generate hypotheses, and identify gaps in our understanding of the progressive process of hindbrain segmentation. The network can be broadly separated into a series of interconnected programs that govern early signaling, segmental subdivision, secondary signaling, segmentation, and ultimately specification of segmental identity. Hox genes play crucial roles in multiple programs within this network. Furthermore, the network reveals properties and principles that are likely to be general to other complex developmental systems. Data from vertebrate and invertebrate chordate models are shedding light on the origin and diversification of the network. Comprehensive cis-regulatory analyses of vertebrate Hox gene regulation have enabled powerful cross-species gene regulatory comparisons. Such an approach in the sea lamprey has revealed that the network mediating segmental Hox expression was present in ancestral vertebrates and has been maintained across diverse vertebrate lineages. Invertebrate chordates lack hindbrain segmentation but exhibit conservation of some aspects of the network, such as a role for retinoic acid in establishing nested Hox expression domains. These comparisons lead to a model in which early vertebrates underwent an elaboration of the network between anterior-posterior patterning and Hox gene expression, leading to the gene-regulatory programs for segmental subdivision and rhombomeric segmentation. WIREs Dev Biol 2017, 6:e286. doi: 10.1002/wdev.286 For further resources related to this article, please visit the WIREs website.

An evolutionary perspective on the systems of adaptive immunity

We propose an evolutionary perspective to classify and characterize the diverse systems of adaptive immunity that have been discovered across all major domains of life. We put forward a new function-based classification according to the way information is acquired by the immune systems: Darwinian immunity (currently known from, but not necessarily limited to, vertebrates) relies on the Darwinian process of clonal selection to 'learn' by cumulative trial-and-error feedback; Lamarckian immunity uses templated targeting (guided adaptation) to internalize heritable information on potential threats; finally, shotgun immunity operates through somatic mechanisms of variable targeting without feedback. We argue that the origin of Darwinian (but not Lamarckian or shotgun) immunity represents a radical innovation in the evolution of individuality and complexity, and propose to add it to the list of major evolutionary transitions. While transitions to higher-level units entail the suppression of selection at lower levels, Darwinian immunity re-opens cell-level selection within the multicellular organism, under the control of mechanisms that direct, rather than suppress, cell-level evolution for the benefit of the individual. From a conceptual point of view, the origin of Darwinian immunity can be regarded as the most radical transition in the history of life, in which evolution by natural selection has literally re-invented itself. Furthermore, the combination of clonal selection and somatic receptor diversity enabled a transition from limited to practically unlimited capacity to store information about the antigenic environment. The origin of Darwinian immunity therefore comprises both a transition in individuality and the emergence of a new information system - the two hallmarks of major evolutionary transitions. Finally, we present an evolutionary scenario for the origin of Darwinian immunity in vertebrates. We propose a revival of the concept of the 'Big Bang' of vertebrate immunity, arguing that its origin involved a 'difficult' (i.e. low-probability) evolutionary transition that might have occurred only once, in a common ancestor of all vertebrates. In contrast to the original concept, we argue that the limiting innovation was not the generation of somatic diversity, but the regulatory circuitry needed for the safe operation of amplifiable immune responses with somatically acquired targeting. Regulatory complexity increased abruptly by genomic duplications at the root of the vertebrate lineage, creating a rare opportunity to establish such circuitry. We discuss the selection forces that might have acted at the origin of the transition, and in the subsequent stepwise evolution leading to the modern immune systems of extant vertebrates.

Evolution of the cAMP-dependent protein kinase (PKA) catalytic subunit isoforms

The 3’,5’-cyclic adenosine monophosphate (cAMP)-dependent protein kinase, or protein kinase A (PKA), pathway is one of the most versatile and best studied signaling pathways in eukaryotic cells. The two paralogous PKA catalytic subunits Cα and Cβ, encoded by the genes PRKACA and PRKACB, respectively, are among the best understood model kinases in signal transduction research. In this work, we explore and elucidate the evolution of the alternative 5’ exons and the splicing pattern giving rise to the numerous PKA catalytic subunit isoforms. In addition to the universally conserved Cα1/Cβ1 isoforms, we find kinase variants with short N-termini in all main vertebrate classes, including the sperm-specific Cα2 isoform found to be conserved in all mammals. We also describe, for the first time, a PKA Cα isoform with a long N-terminus, paralogous to the PKA Cβ2 N-terminus. An analysis of isoform-specific variation highlights residues and motifs that are likely to be of functional importance.

Draft sequencing and assembly of the genome of the world's largest fish, the whale shark: Rhincodon typus Smith 1828

BACKGROUND

The whale shark (Rhincodon typus) has by far the largest body size of any elasmobranch (shark or ray) species. Therefore, it is also the largest extant species of the paraphyletic assemblage commonly referred to as fishes. As both a phenotypic extreme and a member of the group Chondrichthyes - the sister group to the remaining gnathostomes, which includes all tetrapods and therefore also humans - its genome is of substantial comparative interest. Whale sharks are also listed as an endangered species on the International Union for Conservation of Nature's Red List of threatened species and are of growing popularity as both a target of ecotourism and as a charismatic conservation ambassador for the pelagic ecosystem. A genome map for this species would aid in defining effective conservation units and understanding global population structure.

RESULTS

We characterised the nuclear genome of the whale shark using next generation sequencing (454, Illumina) and de novo assembly and annotation methods, based on material collected from the Georgia Aquarium. The data set consisted of 878,654,233 reads, which yielded a draft assembly of 1,213,200 contigs and 997,976 scaffolds. The estimated genome size was 3.44Gb. As expected, the proteome of the whale shark was most closely related to the only other complete genome of a cartilaginous fish, the holocephalan elephant shark. The whale shark contained a novel Toll-like-receptor (TLR) protein with sequence similarity to both the TLR4 and TLR13 proteins of mammals and TLR21 of teleosts. The data are publicly available on GenBank, FigShare, and from the NCBI Short Read Archive under accession number SRP044374.

CONCLUSIONS

This represents the first shotgun elasmobranch genome and will aid studies of molecular systematics, biogeography, genetic differentiation, and conservation genetics in this and other shark species, as well as providing comparative data for studies of evolutionary biology and immunology across the jawed vertebrate lineages.

CTCF binding landscape in jawless fish with reference to Hox cluster evolution

The nuclear protein CCCTC-binding factor (CTCF) contributes as an insulator to chromatin organization in animal genomes. Currently, our knowledge of its binding property is confined mainly to mammals. In this study, we identified CTCF homologs in extant jawless fishes and performed ChIP-seq for the CTCF protein in the Arctic lamprey. Our phylogenetic analysis suggests that the lamprey lineage experienced gene duplication that gave rise to its unique paralog, designated CTCF2, which is independent from the previously recognized duplication between CTCF and CTCFL. The ChIP-seq analysis detected comparable numbers of CTCF binding sites between lamprey, chicken, and human, and revealed that the lamprey CTCF protein binds to the two-part motif, consisting of core and upstream motifs previously reported for mammals. These findings suggest that this mode of CTCF binding was established in the last common ancestor of extant vertebrates (more than 500 million years ago). We analyzed CTCF binding inside Hox clusters, which revealed a reinforcement of CTCF binding in the region spanning Hox1-4 genes that is unique to lamprey. Our study provides not only biological insights into the antiquity of CTCF-based epigenomic regulation known in mammals but also a technical basis for comparative epigenomic studies encompassing the whole taxon Vertebrata.

Uncovering missing pieces: duplication and deletion history of arrestins in deuterostomes

BACKGROUND

The cytosolic arrestin proteins mediate desensitization of activated G protein-coupled receptors (GPCRs) via competition with G proteins for the active phosphorylated receptors. Arrestins in active, including receptor-bound, conformation are also transducers of signaling. Therefore, this protein family is an attractive therapeutic target. The signaling outcome is believed to be a result of structural and sequence-dependent interactions of arrestins with GPCRs and other protein partners. Here we elucidated the detailed evolution of arrestins in deuterostomes.

RESULTS

Identity and number of arrestin paralogs were determined searching deuterostome genomes and gene expression data. In contrast to standard gene prediction methods, our strategy first detects exons situated on different scaffolds and then solves the problem of assigning them to the correct gene. This increases both the completeness and the accuracy of the annotation in comparison to conventional database search strategies applied by the community. The employed strategy enabled us to map in detail the duplication- and deletion history of arrestin paralogs including tandem duplications, pseudogenizations and the formation of retrogenes. The two rounds of whole genome duplications in the vertebrate stem lineage gave rise to four arrestin paralogs. Surprisingly, visual arrestin ARR3 was lost in the mammalian clades Afrotheria and Xenarthra. Duplications in specific clades, on the other hand, must have given rise to new paralogs that show signatures of diversification in functional elements important for receptor binding and phosphate sensing.

CONCLUSION

The current study traces the functional evolution of deuterostome arrestins in unprecedented detail. Based on a precise re-annotation of the exon-intron structure at nucleotide resolution, we infer the gain and loss of paralogs and patterns of conservation, co-variation and selection.

BACH transcription factors in innate and adaptive immunity

BTB and CNC homology (BACH) proteins are transcriptional repressors of the basic region leucine zipper (bZIP) transcription factor family. Recent studies indicate widespread roles of BACH proteins in controlling the development and function of the innate and adaptive immune systems, including the differentiation of effector and memory cells of the B and T cell lineages, CD4+ regulatory T cells and macrophages. Here, we emphasize similarities at a molecular level in the cell-type-specific activities of BACH factors, proposing that competitive interactions of BACH proteins with transcriptional activators of the bZIP family form a common mechanistic theme underlying their diverse actions. The findings contribute to a general understanding of how transcriptional repressors shape lineage commitment and cell-type-specific functions through repression of alternative lineage programmes.

Evolution of the vertebrate insulin receptor substrate (Irs) gene family

BACKGROUND

Insulin receptor substrate (Irs) proteins are essential for insulin signaling as they allow downstream effectors to dock with, and be activated by, the insulin receptor. A family of four Irs proteins have been identified in mice, however the gene for one of these, IRS3, has been pseudogenized in humans. While it is known that the Irs gene family originated in vertebrates, it is not known when it originated and which members are most closely related to each other. A better understanding of the evolution of Irs genes and proteins should provide insight into the regulation of metabolism by insulin.

RESULTS

Multiple genes for Irs proteins were identified in a wide variety of vertebrate species. Phylogenetic and genomic neighborhood analyses indicate that this gene family originated very early in vertebrae evolution. Most Irs genes were duplicated and retained in fish after the fish-specific genome duplication. Irs genes have been lost of various lineages, including Irs3 in primates and birds and Irs1 in most fish. Irs3 and Irs4 experienced an episode of more rapid protein sequence evolution on the ancestral mammalian lineage. Comparisons of the conservation of the proteins sequences among Irs paralogs show that domains involved in binding to the plasma membrane and insulin receptors are most strongly conserved, while divergence has occurred in sequences involved in interacting with downstream effector proteins.

CONCLUSIONS

The Irs gene family originated very early in vertebrate evolution, likely through genome duplications, and in parallel with duplications of other components of the insulin signaling pathway, including insulin and the insulin receptor. While the N-terminal sequences of these proteins are conserved among the paralogs, changes in the C-terminal sequences likely allowed changes in biological function.

SCPP Genes and Their Relatives in Gar: Rapid Expansion of Mineralization Genes in Osteichthyans

Gar is an actinopterygian that has bone, dentin, enameloid, and ganoin (enamel) in teeth and/or scales. Mineralization of these tissues involves genes encoding various secretory calcium-binding phosphoproteins (SCPPs) in osteichthyans, but no SCPP genes have been identified in chondrichthyans to date. In the gar genome, we identified 38 SCPP genes, seven of which encode "acidic-residue-rich" proteins and 31 encode "Pro/Gln (P/Q) rich" proteins. These gar SCPP genes constitute the largest known repertoire, including many newly identified P/Q-rich genes expressed in teeth and/or scales. Among gar SCPP genes, six acidic and three P/Q-rich genes were identified as orthologs of sarcopterygian genes. The sarcopterygian orthologs of most of these acidic genes are involved in bone and/or dentin formation, and sarcopterygian orthologs of all three P/Q-rich genes participate in enamel formation. The finding of these genes in gar suggests that an elaborate SCPP gene-based genetic system for tissue mineralization was already present in stem osteichthyans. While SCPP genes have been thought to originate from ancient SPARCL1, SPARCL1L1 appears to be more closely related to these genes, because it established a structure similar to acidic SCPP genes probably in stem gnathostomes, perhaps at about the same time with the origin of tissue mineralization. Assuming enamel evolved in stem osteichthyans, all P/Q-rich SCPP genes likely arose within the osteichthyan lineage. Furthermore, the absence of acidic SCPP genes in chondrichthyans might be explained by the secondary loss of earliest acidic genes. It appears that many SCPP genes expanded rapidly in stem osteichthyans and in basal actinopterygians.

Lamprey neural crest migration is Snail-dependent and occurs without a differential shift in cadherin expression

The acquisition of neural crest cells was a key step in the origin of the vertebrate body plan. An outstanding question is how neural crest cells acquired their ability to undergo an epithelial-mesenchymal transition (EMT) and migrate extensively throughout the vertebrate embryo. We tested if differential regulation of classical cadherins-a highly conserved feature of neural crest EMT and migration in jawed vertebrates-mediates these cellular behaviors in lamprey, a basal jawless vertebrate. Lamprey has single copies of the type I and type II classical cadherins (CadIA and CadIIA). CadIIA is expressed in premigratory neural crest, and requires the transcription factor Snail for proper expression, yet CadIA is never expressed in the neural tube during neural crest development, suggesting that differential regulation of classical cadherin expression is not required to initiate neural crest migration in basal vertebrates. We hypothesize that neural crest cells evolved by retention of regulatory programs linking distinct mesenchymal and multipotency properties, and emigrated from the neural tube without differentially regulating type I/type II cadherins. Our results point to the coupling of mesenchymal state and multipotency as a key event facilitating the origin of migratory neural crest cells.

Lineage-specific rediploidization is a mechanism to explain time-lags between genome duplication and evolutionary diversification

BACKGROUND

The functional divergence of duplicate genes (ohnologues) retained from whole genome duplication (WGD) is thought to promote evolutionary diversification. However, species radiation and phenotypic diversification are often temporally separated from WGD. Salmonid fish, whose ancestor underwent WGD by autotetraploidization ~95 million years ago, fit such a 'time-lag' model of post-WGD radiation, which occurred alongside a major delay in the rediploidization process. Here we propose a model, 'lineage-specific ohnologue resolution' (LORe), to address the consequences of delayed rediploidization. Under LORe, speciation precedes rediploidization, allowing independent ohnologue divergence in sister lineages sharing an ancestral WGD event.

RESULTS

Using cross-species sequence capture, phylogenomics and genome-wide analyses of ohnologue expression divergence, we demonstrate the major impact of LORe on salmonid evolution. One-quarter of each salmonid genome, harbouring at least 4550 ohnologues, has evolved under LORe, with rediploidization and functional divergence occurring on multiple independent occasions >50 million years post-WGD. We demonstrate the existence and regulatory divergence of many LORe ohnologues with functions in lineage-specific physiological adaptations that potentially facilitated salmonid species radiation. We show that LORe ohnologues are enriched for different functions than 'older' ohnologues that began diverging in the salmonid ancestor.

CONCLUSIONS

LORe has unappreciated significance as a nested component of post-WGD divergence that impacts the functional properties of genes, whilst providing ohnologues available solely for lineage-specific adaptation. Under LORe, which is predicted following many WGD events, the functional outcomes of WGD need not appear 'explosively', but can arise gradually over tens of millions of years, promoting lineage-specific diversification regimes under prevailing ecological pressures.

The origin and evolution of human glutaminases and their atypical C-terminal ankyrin repeats

On the basis of tissue-specific enzyme activity and inhibition by catalytic products, Hans Krebs first demonstrated the existence of multiple glutaminases in mammals. Currently, two human genes are known to encode at least four glutaminase isoforms. However, the phylogeny of these medically relevant enzymes remains unclear, prompting us to investigate their origin and evolution. Using prokaryotic and eukaryotic glutaminase sequences, we built a phylogenetic tree whose topology suggested that the multidomain architecture was inherited from bacterial ancestors, probably simultaneously with the hosting of the proto-mitochondrion endosymbiont. We propose an evolutionary model wherein the appearance of the most active enzyme isoform, glutaminase C (GAC), which is expressed in many cancers, was a late retrotransposition event that occurred in fishes from the Chondrichthyes class. The ankyrin (ANK) repeats in the glutaminases were acquired early in their evolution. To obtain information on ANK folding, we solved two high-resolution structures of the ANK repeat-containing C termini of both kidney-type glutaminase (KGA) and GLS2 isoforms (glutaminase B and liver-type glutaminase). We found that the glutaminase ANK repeats form unique intramolecular contacts through two highly conserved motifs; curiously, this arrangement occludes a region usually involved in ANK-mediated protein-protein interactions. We also solved the crystal structure of full-length KGA and present a small-angle X-ray scattering model for full-length GLS2. These structures explain these proteins' compromised ability to assemble into catalytically active supra-tetrameric filaments, as previously shown for GAC. Collectively, these results provide information about glutaminases that may aid in the design of isoform-specific glutaminase inhibitors.

Molecular cytogenetic characterization of chromosome site-specific repetitive sequences in the Arctic lamprey (Lethenteron camtschaticum, Petromyzontidae)

All extant lamprey karyotypes are characterized by almost all dot-shaped microchromosomes. To understand the molecular basis of chromosome structure in lampreys, we performed chromosome C-banding and silver staining and chromosome mapping of the 18S–28S and 5S ribosomal RNA (rRNA) genes and telomeric TTAGGG repeats in the Arctic lamprey (Lethenteron camtschaticum). In addition, we cloned chromosome site-specific repetitive DNA sequences and characterized them by nucleotide sequencing, chromosome in situ hybridization, and filter hybridization. Three types of repetitive sequences were detected; a 200-bp AT-rich repetitive sequence, LCA-EcoRIa that co-localized with the 18S–28S rRNA gene clusters of 3 chromosomal pairs; a 364-bp AT-rich LCA-EcoRIb sequence that showed homology to the EcoRI sequence family from the sea lamprey (Petromyzon marinus), which contains short repeats as centromeric motifs; and a GC-rich 702-bp LCA-ApaI sequence that was distributed on nearly all chromosomes and showed significant homology with the integrase-coding region of a Ty3/Gypsy family long terminal repeat (LTR) retrotransposon. All three repetitive sequences are highly conserved within the Petromyzontidae or within Petromyzontidae and Mordaciidae. Molecular cytogenetic characterization of these site-specific repeats showed that they may be correlated with programed genome rearrangement (LCA-EcoRIa), centromere structure and function (LCA-EcoRIb), and site-specific amplification of LTR retroelements through homogenization between non-homologous chromosomes (LCA-ApaI).

IL17 factors are early regulators in the gut epithelium during inflammatory response to Vibrio in the sea urchin larva

IL17 cytokines are central mediators of mammalian immunity. In vertebrates, these factors derive from diverse cellular sources. Sea urchins share a molecular heritage with chordates that includes the IL17 system. Here, we characterize the role of epithelial expression of IL17 in the larval gut-associated immune response. The purple sea urchin genome encodes 10 IL17 subfamilies (35 genes) and 2 IL17 receptors. Most of these subfamilies are conserved throughout echinoderms. Two IL17 subfamilies are sequentially strongly upregulated and attenuated in the gut epithelium in response to bacterial disturbance. IL17R1 signal perturbation results in reduced expression of several response genes including an IL17 subtype, indicating a potential feedback. A third IL17 subfamily is activated in adult immune cells indicating that expression in immune cells and epithelia is divided among families. The larva provides a tractable model to investigate the regulation and consequences of gut epithelial IL17 expression across the organism.

DOI:http://dx.doi.org/10.7554/eLife.23481.001

To protect themselves from the constant invasion of harmful microbes, animals have evolved complex immune systems. The gut is one of the most active sites of the immune system and plays a key role in regulating immune responses. In mammals, cells lining the gut wall can sense the presence of harmful bacteria and communicate this information to tissues across the body by producing specialized proteins called Interleukin-17 (IL-17). IL-17 proteins are important for regulating inflammation and are thought to activate specific immune cells in an infected area.

Some aspects of immune systems are similar between different animal species, which can provide clues of how immunity evolved and how it is regulated. For example, sea urchins, which evolved 400-600 million years ago, begin life as simple larvae consisting of a few thousand cells. As oceans harbor a multitude of bacteria and viruses, sea urchin larvae need an efficient immune system to defend themselves. These larvae can respond to specific types of bacteria within a few hours after the microbes have entered their gut by modifying gene expression in distant cells. As these changes occur in cells that are removed from the bacteria, it is thought that the gut cells that initially sense the bacteria, somehow communicate this information.

Now, Buckley et al. exposed sea urchin larvae to a marine bacterium and measured the responses of the cells and their gene expression. The infection affected several types of cells, and in the first 24 hours, a subset of immune cells changed shape and started migrating to the gut wall. In addition, IL-17 gene expression changed significantly in gut cells in the early phases of the larval immune response. Buckley et al. identified three types of IL-17 proteins involved in sea urchin immunity: two that are important for the immune response in the gut during the larval stage, and a third that is only present in adults. These findings suggest that IL-17 signaling is an ancient and central element of gut-associated immune response, which even exists in animals that evolved long before humans.

These findings demonstrate that the sea urchin larva represents a unique and ideal experimental model to study immune responses in a living organism that is more closely related to mammals than some other models, like fruit flies or worms. By understanding the fundamental mechanisms that mediate gut health, this work may highlight new drug targets to treat conditions like Crohn’s disease and colon cancer.

DOI:http://dx.doi.org/10.7554/eLife.23481.002

Evolution and Development of the Inner Ear Efferent System: Transforming a Motor Neuron Population to Connect to the Most Unusual Motor Protein via Ancient Nicotinic Receptors

All craniate chordates have inner ears with hair cells that receive input from the brain by cholinergic centrifugal fibers, the so-called inner ear efferents (IEEs). Comparative data suggest that IEEs derive from facial branchial motor (FBM) neurons that project to the inner ear instead of facial muscles. Developmental data showed that IEEs develop adjacent to FBMs and segregation from IEEs might depend on few transcription factors uniquely associated with IEEs. Like other cholinergic terminals in the peripheral nervous system (PNS), efferent terminals signal on hair cells through nicotinic acetylcholine channels, likely composed out of alpha 9 and alpha 10 units (Chrna9, Chrna10). Consistent with the evolutionary ancestry of IEEs is the even more conserved ancestry of Chrna9 and 10. The evolutionary appearance of IEEs may reflect access of FBMs to a novel target, possibly related to displacement or loss of mesoderm-derived muscle fibers by the ectoderm-derived ear vesicle. Experimental transplantations mimicking this possible aspect of ear evolution showed that different motor neurons of the spinal cord or brainstem form cholinergic synapses on hair cells when ears replace somites or eyes. Transplantation provides experimental evidence in support of the evolutionary switch of FBM neurons to become IEEs. Mammals uniquely evolved a prestin related motor system to cause shape changes in outer hair cells regulated by the IEEs. In summary, an ancient motor neuron population drives in craniates via signaling through highly conserved Chrna receptors a uniquely derived cellular contractility system that is essential for hearing in mammals.

The elephant shark methylome reveals conservation of epigenetic regulation across jawed vertebrates

BACKGROUND

Methylation of CG dinucleotides constitutes a critical system of epigenetic memory in bony vertebrates, where it modulates gene expression and suppresses transposon activity. The genomes of studied vertebrates are pervasively hypermethylated, with the exception of regulatory elements such as transcription start sites (TSSs), where the presence of methylation is associated with gene silencing. This system is not found in the sparsely methylated genomes of invertebrates, and establishing how it arose during early vertebrate evolution is impeded by a paucity of epigenetic data from basal vertebrates.

METHODS

 We perform whole-genome bisulfite sequencing to generate the first genome-wide methylation profiles of a cartilaginous fish, the elephant shark Callorhinchus milii. Employing these to determine the elephant shark methylome structure and its relationship with expression, we compare this with higher vertebrates and an invertebrate chordate using published methylation and transcriptome data.  Results: Like higher vertebrates, the majority of elephant shark CG sites are highly methylated, and methylation is abundant across the genome rather than patterned in the mosaic configuration of invertebrates. This global hypermethylation includes transposable elements and the bodies of genes at all expression levels. Significantly, we document an inverse relationship between TSS methylation and expression in the elephant shark, supporting the presence of the repressive regulatory architecture shared by higher vertebrates.

CONCLUSIONS

 Our demonstration that methylation patterns in a cartilaginous fish are characteristic of higher vertebrates imply the conservation of this epigenetic modification system across jawed vertebrates separated by 465 million years of evolution. In addition, these findings position the elephant shark as a valuable model to explore the evolutionary history and function of vertebrate methylation.

Evolution of the vertebrate phototransduction cascade activation steps

We examine the molecular phylogeny of the proteins underlying the activation steps of vertebrate phototransduction, for both agnathan and jawed vertebrate taxa. We expand the number of taxa analysed and we update the alignment and tree building methodology from a previous analysis. For each of the four primary components (the G-protein transducin alpha subunit, Gα_T, the cyclic GMP phosphodiesterase, PDE6, and the alpha and beta subunits of the cGMP-gated ion channel, CNGC), the phylogenies appear consistent with expansion from an ancestral proto-vertebrate cascade during two rounds of whole-genome duplication followed by divergence of the agnathan and jawed vertebrate lineages. In each case, we consider possible scenarios for the underlying gene duplications and losses, and we apply relevant constraints to the tree construction. From tests of the topology of the resulting trees, we obtain a scenario for the expansion of each component during 2R that accurately fits the observations. Similar analysis of the visual opsins indicates that the only expansion to have occurred during 2R was the formation of Rh1 and Rh2. Finally, we propose a hypothetical scenario for the conversion of an ancestral chordate cascade into the proto-vertebrate phototransduction cascade, prior to whole-genome duplication. Together, our models provide a plausible account for the origin and expansion of the vertebrate phototransduction cascade.

Premammalian origin of the sperm-specific Slo3 channel

Slo3 is a sperm-specific potassium (K+) channel essential for male fertility. Slo3 channels have so far been considered to be specific to mammals. Through exploratory genomics, we identified the Slo3 gene in the genome of terrestrial (birds and reptiles) and aquatic (fish) vertebrates. In the case of fish, Slo3 has undergone several episodes of gene loss. Transcriptomic analysis showed that vertebrate Slo3 transcript orthologues are predominantly expressed in testis, in concordance with the mammalian Slo3. We conclude that the Slo3 gene arose during the radiation of early vertebrates, much earlier than previously thought. Our findings add to the growing evidence indicating that the phylogenetic profiles of sperm-specific channels are intermittent throughout metazoan evolution, which probably reflects the adaptation of sperm to different ionic milieus and fertilization environments.

Autopolyploidy genome duplication preserves other ancient genome duplications in Atlantic salmon (Salmo salar)

Salmonids (e.g. Atlantic salmon, Pacific salmon, and trouts) have a long legacy of genome duplication. In addition to three ancient genome duplications that all teleosts are thought to share, salmonids have had one additional genome duplication. We explored a methodology for untangling these duplications from each other to better understand them in Atlantic salmon. In this methodology, homeologous regions (paralogous/duplicated genomic regions originating from a whole genome duplication) from the most recent genome duplication were assumed to have duplicated genes at greater density and have greater sequence similarity. This assumption was used to differentiate duplicated gene pairs in Atlantic salmon that are either from the most recent genome duplication or from earlier duplications. From a comparison with multiple vertebrate species, it is clear that Atlantic salmon have retained more duplicated genes from ancient genome duplications than other vertebrates--often at higher density in the genome and containing fewer synonymous mutations. It may be that polysomic inheritance is the mechanism responsible for maintaining ancient gene duplicates in salmonids. Polysomic inheritance (when multiple chromosomes pair during meiosis) is thought to be relatively common in salmonids compared to other vertebrate species. These findings illuminate how genome duplications may not only increase the number of duplicated genes, but may also be involved in the maintenance of them from previous genome duplications as well.

Aquaculture genomics, genetics and breeding in the United States: current status, challenges, and priorities for future research

Advancing the production efficiency and profitability of aquaculture is dependent upon the ability to utilize a diverse array of genetic resources. The ultimate goals of aquaculture genomics, genetics and breeding research are to enhance aquaculture production efficiency, sustainability, product quality, and profitability in support of the commercial sector and for the benefit of consumers. In order to achieve these goals, it is important to understand the genomic structure and organization of aquaculture species, and their genomic and phenomic variations, as well as the genetic basis of traits and their interrelationships. In addition, it is also important to understand the mechanisms of regulation and evolutionary conservation at the levels of genome, transcriptome, proteome, epigenome, and systems biology. With genomic information and information between the genomes and phenomes, technologies for marker/causal mutation-assisted selection, genome selection, and genome editing can be developed for applications in aquaculture. A set of genomic tools and resources must be made available including reference genome sequences and their annotations (including coding and non-coding regulatory elements), genome-wide polymorphic markers, efficient genotyping platforms, high-density and high-resolution linkage maps, and transcriptome resources including non-coding transcripts. Genomic and genetic control of important performance and production traits, such as disease resistance, feed conversion efficiency, growth rate, processing yield, behaviour, reproductive characteristics, and tolerance to environmental stressors like low dissolved oxygen, high or low water temperature and salinity, must be understood. QTL need to be identified, validated across strains, lines and populations, and their mechanisms of control understood. Causal gene(s) need to be identified. Genetic and epigenetic regulation of important aquaculture traits need to be determined, and technologies for marker-assisted selection, causal gene/mutation-assisted selection, genome selection, and genome editing using CRISPR and other technologies must be developed, demonstrated with applicability, and application to aquaculture industries.Major progress has been made in aquaculture genomics for dozens of fish and shellfish species including the development of genetic linkage maps, physical maps, microarrays, single nucleotide polymorphism (SNP) arrays, transcriptome databases and various stages of genome reference sequences. This paper provides a general review of the current status, challenges and future research needs of aquaculture genomics, genetics, and breeding, with a focus on major aquaculture species in the United States: catfish, rainbow trout, Atlantic salmon, tilapia, striped bass, oysters, and shrimp. While the overall research priorities and the practical goals are similar across various aquaculture species, the current status in each species should dictate the next priority areas within the species. This paper is an output of the USDA Workshop for Aquaculture Genomics, Genetics, and Breeding held in late March 2016 in Auburn, Alabama, with participants from all parts of the United States.

Diversity as Opportunity: Insights from 600 Million Years of AHR Evolution

The aryl hydrocarbon receptor (AHR) was for many years of interest only to pharmacologists and toxicologists. However, this protein has fundamental roles in biology that are being revealed through studies in diverse animal species. The AHR is an ancient protein. AHR homologs exist in most major groups of modern bilaterian animals, including deuterostomes (chordates, hemichordates, echinoderms) and the two major clades of protostome invertebrates [ecdysozoans (e.g. arthropods and nematodes) and lophotrochozoans (e.g. molluscs and annelids)]. AHR homologs also have been identified in cnidarians such as the sea anemone Nematostella and in the genome of Trichoplax, a placozoan. Bilaterians, cnidarians, and placozoans form the clade Eumetazoa, whose last common ancestor lived approximately 600 million years ago (MYA). The presence of AHR homologs in modern representatives of all these groups indicates that the original eumetazoan animal possessed an AHR homolog. Studies in invertebrates and vertebrates reveal parallel functions of AHR in the development and function of sensory neural systems, suggesting that these may be ancestral roles. Vertebrate animals are characterized by the expansion and diversification of AHRs, via gene and genome duplications, from the ancestral protoAHR into at least five classes of AHR-like proteins: AHR, AHR1, AHR2, AHR3, and AHRR. The evolution of multiple AHRs in vertebrates coincided with the acquisition of high-affinity binding of halogenated and polynuclear aromatic hydrocarbons and the emergence of adaptive functions involving regulation of xenobiotic-metabolizing enzymes and roles in adaptive immunity. The existence of multiple AHRs may have facilitated subfunction partitioning and specialization of specific AHR types in some taxa. Additional research in diverse model and non-model species will continue to enrich our understanding of AHR and its pleiotropic roles in biology and toxicology.

A SNP Based Linkage Map of the Arctic Charr (Salvelinus alpinus) Genome Provides Insights into the Diploidization Process After Whole Genome Duplication

Diploidization, which follows whole genome duplication events, does not occur evenly across the genome. In salmonid fishes, certain pairs of homeologous chromosomes preserve tetraploid loci in higher frequencies toward the telomeres due to residual tetrasomic inheritance. Research suggests this occurs only in homeologous pairs where one chromosome arm has undergone a fusion event. We present a linkage map for Arctic charr (Salvelinus alpinus), a salmonid species with relatively fewer chromosome fusions. Genotype by sequencing identified 19,418 SNPs, and a linkage map consisting of 4508 markers was constructed from a subset of high quality SNPs and microsatellite markers that were used to anchor the new map to previous versions. Both male- and female-specific linkage maps contained the expected number of 39 linkage groups. The chromosome type associated with each linkage group was determined, and 10 stable metacentric chromosomes were identified, along with a chromosome polymorphism involving the sex chromosome AC04. Two instances of a weak form of pseudolinkage were detected in the telomeric regions of homeologous chromosome arms in both female and male linkage maps. Chromosome arm homologies within the Atlantic salmon (Salmo salar) and rainbow trout (Oncorhynchus mykiss) genomes were determined. Paralogous sequence variants (PSVs) were identified, and their comparative BLASTn hit locations showed that duplicate markers exist in higher numbers on seven pairs of homeologous arms, previously identified as preserving tetrasomy in salmonid species. Homeologous arm pairs where neither arm has been part of a fusion event in Arctic charr had fewer PSVs, suggesting faster diploidization rates in these regions.

A role for tight junction-associated MARVEL proteins in larval sea lamprey (Petromyzon marinus) osmoregulation

This study reports on tight junction-associated MARVEL proteins of larval sea lamprey (Petromyzon marinus) and their potential role in ammocoete osmoregulation. Two Occludin isoforms (designated Ocln and Ocln-a) and a tricellulin (Tric) were identified. Transcripts encoding ocln, ocln-a, and tric were broadly expressed in larval lamprey, with greatest abundance of ocln in gut, liver and kidney, ocln-a in the gill and skin, and tric in the kidney. Ocln and Ocln-a resolved as ∼63 kDa and ∼35 kDa MW proteins respectively while Tric resolved as a ∼50 kDa protein. Ocln immunolocalized to the gill vasculature and in gill mucous cells while Ocln-a localized to the gill pouch and gill epithelium. Both Ocln and Ocln-a localized in the nephron, the epidermis and the luminal side of the gut. In branchial tissue, Tric exhibited punctate localization, consistent with its presence at regions of tricellular contact. Following ion-poor water (IPW) acclimation of ammocoetes, serum [Na+] and [Cl−] reduced, but not [Ca++], and carcass moisture content increased. In association, Ocln abundance increased in skin and kidney, but reduced in gill of IPW-acclimated ammocoetes while Ocln-a abundance reduced in the kidney only. Tric abundance increased in the gill. Region-specific alterations in ocln, ocln-a and tric mRNA abundance was also observed in the gut. Data support a role for Ocln, Ocln-a and Tric in the osmoregulatory strategies of a basal vertebrate.

Molecular cloning and characterization of hagfish estrogen receptors

One or more distinct forms of the nuclear estrogen receptor (ER) have been isolated from many vertebrates to date. To better understand the molecular evolution of ERs, we cloned and characterized er cDNAs from the inshore hagfish, Eptatretus burgeri, a modern representative of the most primitive vertebrates, the agnathans. Two er cDNAs, er1 and er2, were isolated from the liver of a reproductive female hagfish. A phylogenetic analysis placed hagfish ER1 into a position prior to the divergence of vertebrate ERs. Conversely, hagfish ER2 was placed at the base of the vertebrate ERβ clade. The tissue distribution patterns of both ER subtype mRNAs appeared to be different, suggesting that each subtype has different physiological roles associated with estrogen actions. An estrogen responsive-luciferase reporter assay using mammalian HEK293 cells was used to functionally characterize these hagfish ERs. Both ER proteins displayed estrogen-dependent activation of transcription. These results clearly demonstrate that the hagfish has two functional ER subtypes.

Evolution and Diversity of Transposable Elements in Vertebrate Genomes

Transposable elements (TEs) are selfish genetic elements that mobilize in genomes via transposition or retrotransposition and often make up large fractions of vertebrate genomes. Here, we review the current understanding of vertebrate TE diversity and evolution in the context of recent advances in genome sequencing and assembly techniques. TEs make up 4-60% of assembled vertebrate genomes, and deeply branching lineages such as ray-finned fishes and amphibians generally exhibit a higher TE diversity than the more recent radiations of birds and mammals. Furthermore, the list of taxa with exceptional TE landscapes is growing. We emphasize that the current bottleneck in genome analyses lies in the proper annotation of TEs and provide examples where superficial analyses led to misleading conclusions about genome evolution. Finally, recent advances in long-read sequencing will soon permit access to TE-rich genomic regions that previously resisted assembly including the gigantic, TE-rich genomes of salamanders and lungfishes.

Evolution of the β-adrenoreceptors in vertebrates

The study of the evolutionary history of genes related to human disease lies at the interface of evolution and medicine. These studies provide the evolutionary context on which medical researchers should work, and are also useful in providing information to suggest further genetic experiments, especially in model species where genetic manipulations can be made. Here we studied the evolution of the β-adrenoreceptor gene family in vertebrates with the aim of adding an evolutionary framework to the already abundant physiological information. Our results show that in addition to the three already described vertebrate β-adrenoreceptor genes there is an additional group containing cyclostome sequences. We suggest that β-adrenoreceptors diversified as a product of the two whole genome duplications that occurred in the ancestor of vertebrates. Gene expression patterns are in general consistent across species, suggesting that expression dynamics were established early in the evolutionary history of vertebrates, and have been maintained since then. Finally, amino acid polymorphisms that are associated to pathological conditions in humans appear to be common in non-human mammals, suggesting that the phenotypic effects of these mutations depend on epistatic interaction with other positions. The evolutionary analysis of the β-adrenoreceptors delivers new insights about the diversity of these receptors in vertebrates, the evolution of the expression patterns and a comparative perspective regarding the polymorphisms that in humans are linked to pathological conditions.

Use of physiological knowledge to control the invasive sea lamprey (Petromyzon marinus) in the Laurentian Great Lakes

Sea lamprey (Petromyzon marinus) control in the Laurentian Great Lakes of North America is an example of using physiological knowledge to successfully control an invasive species and rehabilitate an ecosystem and valuable fishery. The parasitic sea lamprey contributed to the devastating collapse of native fish communities after invading the Great Lakes during the 1800s and early 1900s. Economic tragedy ensued with the loss of the fishery and severe impacts to property values and tourism resulting from sea lamprey-induced ecological changes. To control the sea lamprey and rehabilitate the once vibrant Great Lakes ecosystem and economy, the Great Lakes Fishery Commission (Commission) was formed by treaty between Canada and the United States in 1955. The Commission has developed a sea lamprey control programme based on their physiological vulnerabilities, which includes (i) the application of selective pesticides (lampricides), which successfully kill sedentary sea lamprey larvae in their natal streams; (ii) barriers to spawning migrations and associated traps to prevent infestations of upstream habitats and remove adult sea lamprey before they reproduce; and (iii) the release of sterilized males to reduce the reproductive potential of spawning populations in select streams. Since 1958, the application of the sea lamprey control programme has suppressed sea lamprey populations by ~90% from peak abundance. Great Lakes fish populations have rebounded and the economy is now thriving. In hopes of further enhancing the efficacy and selectivity of the sea lamprey control programme, the Commission is exploring the use of (i) sea lamprey chemosensory cues (pheromones and alarm cues) to manipulate behaviours and physiologies, and (ii) genetics to identify and manipulate genes associated with key physiological functions, for control purposes. Overall, the Commission capitalizes on the unique physiology of the sea lamprey and strives to develop a diverse integrated programme to successfully control a once devastating invasive species.

The presence of Anf/Hesx1 homeobox gene in lampreys suggests that it could play an important role in emergence of telencephalon

Accumulated evidence indicates that the core genetic mechanisms regulating early patterning of the brain rudiment in vertebrates are very similar to those operating during development of the anterior region of invertebrate embryos. However, the mechanisms underlying the morphological differences between the elaborate vertebrate brain and its simpler invertebrate counterpart remain poorly understood. Recently, we hypothesized that the emergence of the most anterior unit of the vertebrate brain, the telencephalon, could be related to the appearance in vertebrates’ ancestors of a unique homeobox gene, Anf/Hesx1(further Anf), which is absent from all invertebrates and regulates the earliest steps of telencephalon development in vertebrates. However, the failure of Anf to be detected in one of the most basal extant vertebrate species, the lamprey, seriously compromises this hypothesis. Here, we report the cloning of Anf in three lamprey species and demonstrate that this gene is indeed expressed in embryos in the same pattern as in other vertebrates and executes the same functions by inhibiting the expression of the anterior general regulator Otx2 in favour of the telencephalic regulator FoxG1. These results are consistent with the hypothesis that the Anf homeobox gene may have been important in the evolution of the telencephalon.

Evolutionary and Functional Diversification of the Vitamin D Receptor-Lithocholic Acid Partnership

The evolution, molecular behavior, and physiological function of nuclear receptors are of particular interest given their diverse roles in regulating essential biological processes. The vitamin D receptor (VDR) is well known for its canonical roles in calcium homeostasis and skeletal maintenance. Additionally, VDR has received an increased amount of attention due to the discovery of numerous non-calcemic functions, including the detoxification of lithocholic acid. Lithocholic acid is a toxic metabolite of chenodeoxycholic acid, a primary bile acid. The partnership between the VDR and lithocholic acid has been hypothesized to be a recent adaptation that evolved to mediate the detoxification and elimination of lithocholic acid from the gut. This partnership is speculated to be limited to higher vertebrates (birds and mammals), as lower vertebrates do not synthesize the parent compound of lithocholic acid. However, the molecular functions associated with the observed insensitivity of basal VDRs to lithocholic acid have not been explored. Here we characterize canonical nuclear receptor functions of VDRs from select species representing key nodes in vertebrate evolution and span a range of bile salt phenotypes. Competitive ligand binding assays revealed that the receptor's affinity for lithocholic acid is highly conserved across species, suggesting that lithocholic acid affinity is an ancient and non-adaptive trait. However, transient transactivation assays revealed that lithocholic acid-mediated VDR activation might have evolved more recently, as the non-mammalian receptors did not respond to lithocholic acid unless exogenous coactivator proteins were co-expressed. Subsequent functional assays indicated that differential lithocholic acid-mediated receptor activation is potentially driven by differential protein-protein interactions between VDR and nuclear receptor coregulator proteins. We hypothesize that the vitamin D receptor-lithocholic acid partnership evolved as a by-product of natural selection on the ligand-receptor partnership between the vitamin D receptor and the native VDR ligand: 1α,25-dihydroxyvitamin D3, the biologically active metabolite of vitamin D3.

Cloning of the GABAB Receptor Subunits B1 and B2 and their Expression in the Central Nervous System of the Adult Sea Lamprey

In vertebrates, γ-aminobutyric acid (GABA) is the main inhibitory transmitter in the central nervous system (CNS) acting through ionotropic (GABA_A) and metabotropic (GABA_B) receptors. The GABA_B receptor produces a slow inhibition since it activates second messenger systems through the binding and activation of guanine nucleotide-binding proteins [G-protein-coupled receptors (GPCRs)]. Lampreys are a key reference to understand molecular evolution in vertebrates. The importance of the GABA_B receptor for the modulation of the circuits controlling locomotion and other behaviors has been shown in pharmacological/physiological studies in lampreys. However, there is no data about the sequence of the GABA_B subunits or their expression in the CNS of lampreys. Our aim was to identify the sea lamprey GABA_B1 and GABA_B2 transcripts and study their expression in the CNS of adults. We cloned two partial sequences corresponding to the GABA_B1 and GABA_B2 cDNAs of the sea lamprey as confirmed by sequence analysis and comparison with known GABA_B sequences of other vertebrates. In phylogenetic analyses, the sea lamprey GABA_B sequences clustered together with GABA_Bs sequences of vertebrates and emerged as an outgroup to all gnathostome sequences. We observed a broad and overlapping expression of both transcripts in the entire CNS. Expression was mainly observed in neuronal somas of the periventricular regions including the identified reticulospinal cells. No expression was observed in identifiable fibers. Comparison of our results with those reported in other vertebrates indicates that a broad and overlapping expression of the GABA_B subunits in the CNS is a conserved character shared by agnathans and gnathostomes.

Identification and characteristic of three members of the C1q/TNF-related proteins (CTRPs) superfamily in Eudontomyzon morii

C1q is the target recognition protein of the classical complement pathway and a major connecting link between innate and acquired immunity. C1q and the multifunctional tumor necrosis factor (TNF) ligand family is of similar crystal structures, are designated the C1q/TNF-related proteins (CTRPs) superfamily. They are involved in processes as diverse as host defense, inflammation, apoptosis, autoimmunity, cell differentiation, organogenesis, hibernation and insulinresistant obesity. In this study, three members of the CTRPs superfamily were isolated and characterized in Yalu River lampreys (Eudontomyzon morii), and are respectively named LaC1qC, LaCTRP1, LaCTRP9. The full-length cDNAs of C1qC-like (LaC1qAL), CTRP1-like (LaCTRP1), and CTRP9-like (LaCTRP9) consist of 723, 762 and 825 bp of nucleotide sequence encoding polypeptides of 241, 254 and 275 amino acids, respectively. All-three proteins share three common domains: a signal peptide at the N terminus, a collagenous domain (characteristic Gly-X-Y repeats), and a C-terminal globular domain. In addition, the higher expression level of the three proteins in heart by RT-PCR and real-time PCR tissue profiling implied that they might involve in immune response or injury repair of the heart in lamprey.

EMBRYONIC DEVELOPMENT AND A QUANTITATIVE MODEL OF PROGRAMMED DNA ELIMINATION IN MESOCYCLOPS EDAX (S. A. FORBES, 1891) (COPEPODA: CYCLOPOIDA)

The highly programmed fragmentation of chromosomes and elimination of large amounts of nuclear DNA from the presomatic cell lineages (i.e., chromatin diminution), occurs in the embryos of the freshwater zooplankton Mesocyclops edax (S. A. Forbes, 1891) (Crustacea: Copepoda). The somatic genome is reorganized and reduced to a size five times smaller even though the germline genome remains intact. We present the first comprehensive, quantitative model of DNA content throughout embryogenesis in a copepod that possesses embryonic DNA elimination. We used densitometric image analysis to measure the DNA content of polar bodies, germline and somatic nuclei, and excised DNA "droplets." We report: 1) variable DNA contents of polar bodies, some of which do not contain the amount corresponding to the haploid germline genome size; 2) presence of pronuclei in newly laid embryo sacs; 3) gonomeric chromosomes in the second to fourth cleavage divisions and in the primordial germ cell and primordial endoderm cell during the fifth cleavage division; 4) timing of early embryonic cell stages, elimination of DNA, and divisions of the primordial germ cell and primordial endoderm cell at 22°C; and 5) persistence of a portion of the excised DNA "droplets" throughout embryogenesis. DNA elimination is a trait that spans multiple embryonic stages and a knowledge of the timing and variability of the associated cytological events with DNA elimination will promote the study of the molecular mechanisms involved in this trait. We propose the "genome yolk hypothesis" as a functional explanation for the persistence of the eliminated DNA that might serve as a resource during postdiminution cleavage divisions.

Evolution of the Rax family of developmental transcription factors in vertebrates

Rax proteins comprise a small family of paired-type, homeodomain-containing transcription factors with essential functions in eye and forebrain development. While invertebrates possess only one Rax gene, vertebrates can have several Rax paralogue genes, but the evolutionary history of the members of the family has not been studied in detail. Here, we present a thorough analysis of the evolutionary relationships between vertebrate Rax genes and proteins available in diverse genomic databases. Phylogenetic and synteny analyses indicate that Rax genes went through a duplication in an ancestor of all jawed vertebrates (Gnathostomata), giving rise to the ancestral vertebrate Rax1 and Rax2 genes. This duplication event is likely related to the proposed polyploidisations that occurred during early vertebrate evolution. Subsequent genome-wide duplications in the lineage of ray-finned fish (Actinopterygii) originated new Rax2 paralogues in the genomes of teleosts. In the lobe-finned fish lineage (Sarcopterygii), the N-terminal octapeptide domain of Rax2 was lost in a common ancestor of tetrapods, giving rise to a shorter version of Rax2 in this lineage. Within placental mammals, the Rax2 gene was lost altogether in an ancestor of rodents and lagomorphs (Glires). Finally, we discuss the scientific literature in the light of Rax gene evolution and propose new avenues of research on the function of this important family of transcriptional regulators.

Evolution of the RH gene family in vertebrates revealed by brown hagfish (Eptatretus atami) genome sequences

In vertebrates, there are four major genes in the RH (Rhesus) gene family, RH, RHAG, RHBG, and RHCG. These genes are thought to have been formed by the two rounds of whole-genome duplication (2R-WGD) in the common ancestor of all vertebrates. In our previous work, where we analyzed details of the gene duplications process of this gene family, three nucleotide sequences belonging to this family were identified in Far Eastern brook lamprey (Lethenteron reissneri), and the phylogenetic positions of the genes were determined. Lampreys, along with hagfishes, are cyclostomata (jawless fishes), which is a sister group of gnathostomata (jawed vertebrates). Although those results suggested that one gene was orthologous to the gnathostome RHCG genes, we did not identify clear orthologues for other genes. In this study, therefore, we identified three novel cDNA sequences that belong to the RH gene family using de novo transcriptome analysis of another cyclostome: the brown hagfish (Eptatretus atami). We also determined the nucleotide sequences for the RHBG and RHCG genes in a red stingray (Dasyatis akajei), which belongs to the cartilaginous fishes. The phylogenetic tree showed that two brown hagfish genes, which were probably duplicated in the cyclostome lineage, formed a cluster with the gnathostome RHAG genes, whereas another brown hagfish gene formed a cluster with the gnathostome RHCG genes. We estimated that the RH genes had a higher evolutionary rate than the RHAG, RHBG, and RHCG genes. Interestingly, in the RHBG genes, only the bird lineage showed a higher rate of nonsynonymous substitutions. It is likely that this higher rate was caused by a state of relaxed functional constraints rather than positive selection nor by pseudogenization.

Identification of three somatostatin genes in lampreys

Somatostatins (SSs) are a structurally diverse family of neuropeptides that play important roles in the regulation of growth, development and metabolism in vertebrates. It has been recently proposed that the common ancestor of gnathostomes possessed three SS genes, namely SS1, SS2 and SS5. SS1 and SS2 are still present in most extant gnathostome species investigated so far while SS5 primarily occurs in chondrichthyes, actinopterygians and actinistia but not in tetrapods. Very little is known about the repertoire of SSs in cyclostomes, which are extant jawless vertebrates. In the present study, we report the cloning of the cDNAs encoding three distinct lamprey SS variants that we call SSa, SSb and SSc. SSa and SSb correspond to the two SS variants previously characterized in lamprey, while SSc appears to be a totally novel one. SSa exhibits the same sequence as gnathostome SS1. SSb differs from SSa by only one substitution (Thr¹²→Ser). SSc exhibits a totally unique structure (ANCRMFYWKTMAAC) that shares only 50% identity with SSa and SSb. SSa, SSb and SSc precursors do not exhibit any appreciable sequence similarity outside the C-terminal region containing the SS sequence. Phylogenetic analyses failed to clearly assign orthology relationships between lamprey and gnathostome SS genes. Synteny analysis suggests that the SSc gene arose before the split of the three gnathostome genes SS1, SS2 and SS5.

A pre-metazoan origin of the CRK gene family and co-opted signaling network

CRK and CRKL adapter proteins play essential roles in development and cancer through their SRC homology 2 and 3 (SH2 and SH3) domains. To gain insight into the origin of their shared functions, we have investigated their evolutionary history. We propose a term, crk/crkl ancestral (crka), for orthologs in invertebrates before the divergence of CRK and CRKL in the vertebrate ancestor. We have isolated two orthologs expressed in the choanoflagellate Monosiga brevicollis, a unicellular relative to the metazoans. Consistent with its highly-conserved three-dimensional structure, the SH2 domain of M. brevicollis crka1 can bind to the mammalian CRK/CRKL SH2 binding consensus phospho-YxxP, and to the SRC substrate/focal adhesion protein BCAR1 (p130^CAS) in the presence of activated SRC. These results demonstrate an ancient origin of the CRK/CRKL SH2-target recognition specificity. Although BCAR1 orthologs exist only in metazoans as identified by an N-terminal SH3 domain, YxxP motifs, and a C-terminal FAT-like domain, some pre-metazoan transmembrane proteins include several YxxP repeats in their cytosolic region, suggesting that they are remotely related to the BCAR1 substrate domain. Since the tyrosine kinase SRC also has a pre-metazoan origin, co-option of BCAR1-related sequences may have rewired the crka-dependent network to mediate adhesion signals in the metazoan ancestor.

The Case of the Missing Ancient Fungal Polyploids

Polyploidy-the increase in the number of whole chromosome sets-is an important evolutionary force in eukaryotes. Polyploidy is well recognized throughout the evolutionary history of plants and animals, where several ancient events have been hypothesized to be drivers of major evolutionary radiations. However, fungi provide a striking contrast: while numerous recent polyploids have been documented, ancient fungal polyploidy is virtually unknown. We present a survey of known fungal polyploids that confirms the absence of ancient fungal polyploidy events. Three hypotheses may explain this finding. First, ancient fungal polyploids are indeed rare, with unique aspects of fungal biology providing similar benefits without genome duplication. Second, fungal polyploids are not successful in the long term, leading to few extant species derived from ancient polyploidy events. Third, ancient fungal polyploids are difficult to detect, causing the real contribution of polyploidy to fungal evolution to be underappreciated. We consider each of these hypotheses in turn and propose that failure to detect ancient events is the most likely reason for the lack of observed ancient fungal polyploids. We examine whether existing data can provide evidence for previously unrecognized ancient fungal polyploidy events but discover that current resources are too limited. We contend that establishing whether unrecognized ancient fungal polyploidy events exist is important to ascertain whether polyploidy has played a key role in the evolution of the extensive complexity and diversity observed in fungi today and, thus, whether polyploidy is a driver of evolutionary diversifications across eukaryotes. Therefore, we conclude by suggesting ways to test the hypothesis that there are unrecognized polyploidy events in the deep evolutionary history of the fungi.

A second estrogen receptor from Japanese lamprey (Lethenteron japonicum) does not have activities for estrogen binding and transcription

Estrogens regulate many physiological responses in vertebrates by binding to the estrogen receptor (ER), a ligand-activated transcription factor. To understand the evolution of vertebrate ERs and to investigate how estrogen acts in a jawless vertebrate, we used degenerate primer sets and PCR to isolate DNA fragments encoding two distinct ER subtypes, Esr1a and Esr1b from the Japanese lamprey, Lethenteron japonicum. Phylogenetic analysis indicates that these two ERs are the result of lineage-specific gene duplication within the jawless fishes, different from the previous duplication event of Esr1 (ERα) and Esr2 (ERβ) within the jawed vertebrates. Reporter gene assays show that lamprey Esr1a displays both constitutive and estrogen-dependent activation of gene transcription. Domain swapping experiments indicate that constitutive activity resides in the A/B domain of lamprey Esr1a. Unexpectedly, lamprey Esr1b does not bind estradiol and is not stimulated by other estrogens, androgens or corticosteroids. A 3D model of lamprey Esr1b suggests that although estradiol fits into the steroid binding site, some stabilizing contacts between the ligand and side chains that are found in human Esr1 and Esr2 are missing in lamprey Esr1b.

Neuronal fate specification by the Dbx1 transcription factor is linked to the evolutionary acquisition of a novel functional domain

Background

Dbx1 is a homeodomain transcription factor involved in neuronal fate specification belonging to a widely conserved family among bilaterians. In mammals, Dbx1 was proposed to act as a transcriptional repressor by interacting with the Groucho corepressors to allow the specification of neurons involved in essential biological functions such as locomotion or breathing.

Results

Sequence alignments of Dbx1 proteins from different species allowed us to identify two conserved domains related to the Groucho-dependent Engrailed repressor domain (RD), as well as a newly described domain composed of clusterized acidic residues at the C-terminus (Cter) which is present in tetrapods but also several invertebrates. Using a heterologous luciferase assay, we showed that the two putative repressor domains behave as such in a Groucho-dependent manner, whereas the Cter does not bear any intrinsic transcriptional activity. Consistently with in vitro data, we found that both RDs are involved in cell fate specification using in vivo electroporation experiments in the chick spinal cord. Surprisingly, we show that the Cter domain is required for Dbx1 function in vivo, acting as a modulator of its repressive activity and/or imparting specificity.

Conclusion

Our results strongly suggest that the presence of a Cter domain among tetrapods is essential for Dbx1 to regulate neuronal diversity and, in turn, nervous system complexity.

Electronic supplementary material

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

Efficient Gene Tree Correction Guided by Genome Evolution

MOTIVATIONS

Gene trees inferred solely from multiple alignments of homologous sequences often contain weakly supported and uncertain branches. Information for their full resolution may lie in the dependency between gene families and their genomic context. Integrative methods, using species tree information in addition to sequence information, often rely on a computationally intensive tree space search which forecloses an application to large genomic databases.

RESULTS

We propose a new method, called ProfileNJ, that takes a gene tree with statistical supports on its branches, and corrects its weakly supported parts by using a combination of information from a species tree and a distance matrix. Its low running time enabled us to use it on the whole Ensembl Compara database, for which we propose an alternative, arguably more plausible set of gene trees. This allowed us to perform a genome-wide analysis of duplication and loss patterns on the history of 63 eukaryote species, and predict ancestral gene content and order for all ancestors along the phylogeny.

AVAILABILITY

A web interface called RefineTree, including ProfileNJ as well as a other gene tree correction methods, which we also test on the Ensembl gene families, is available at: http://www-ens.iro.umontreal.ca/~adbit/polytomysolver.html. The code of ProfileNJ as well as the set of gene trees corrected by ProfileNJ from Ensembl Compara version 73 families are also made available.