Genome duplication in early vertebrates: insights from agnathan cytogenetics

The hagfish genome and the evolution of vertebrates

As the only surviving lineages of jawless fishes, hagfishes and lampreys provide a crucial window into early vertebrate evolution1-3. Here we investigate the complex history, timing and functional role of genome-wide duplications4-7 and programmed DNA elimination8,9 in vertebrates in the light of a chromosome-scale genome sequence for the brown hagfish Eptatretus atami. Combining evidence from syntenic and phylogenetic analyses, we establish a comprehensive picture of vertebrate genome evolution, including an auto-tetraploidization (1RV) that predates the early Cambrian cyclostome-gnathostome split, followed by a mid-late Cambrian allo-tetraploidization (2RJV) in gnathostomes and a prolonged Cambrian-Ordovician hexaploidization (2RCY) in cyclostomes. Subsequently, hagfishes underwent extensive genomic changes, with chromosomal fusions accompanied by the loss of genes that are essential for organ systems (for example, genes involved in the development of eyes and in the proliferation of osteoclasts); these changes account, in part, for the simplification of the hagfish body plan1,2. Finally, we characterize programmed DNA elimination in hagfish, identifying protein-coding genes and repetitive elements that are deleted from somatic cell lineages during early development. The elimination of these germline-specific genes provides a mechanism for resolving genetic conflict between soma and germline by repressing germline and pluripotency functions, paralleling findings in lampreys10,11. Reconstruction of the early genomic history of vertebrates provides a framework for further investigations of the evolution of cyclostomes and jawed vertebrates.

Parallel evolution of fish bi-modal breathing and expansion of olfactory receptor (OR) genes: toward a universal ORs nomenclature

Olfactory receptors (ORs) play a key role in the prime sensorial perception, being highly relevant for intra/interspecific interactions. ORs are a subgroup of G-protein coupled receptors that exhibit highly complex subgenomes in vertebrates. However, OR repertoires remain poorly studied in fish lineages, precluding finely retracing their origin, evolution, and diversification, especially in the most basal groups. Here, we conduct an exhaustive gene screening upon 43 high-quality fish genomes exhibiting varied gene repertoires (2-583 genes). While the early vertebrates performed gas exchange through gills, we hypothesize that the emergence of new breathing structures (swim bladder and paired lungs) in early osteichthyans may be associated with expansions in the ORs gene families sensitive to airborne molecules. Additionally, we verify that the OR repertoire of moderns actinopterygians has not increased as expected following a whole genome duplication, likely due to regulatory mechanisms compensating the gene load excess. Finally, we identify 25 distinct OR families, allowing us to propose an updated universal nomenclature for the fish ORs.

An Evolutionary Conserved Signaling Network Between Mouse and Human Underlies the Differential Osteoskeletal Potential of Frontal and Parietal Calvarial Bones

The mammalian calvarial vault is an ancient and highly conserved structure among species, however, the mechanisms governing osteogenesis of the calvarial vault and how they might be conserved across mammalian species remain unclear. The aim of this study was to determine if regional differences in osteogenic potential of the calvarial vault, first described in mice, extend to humans. We derived human frontal and parietal osteoblasts from fetal calvarial tissue, demonstrating enhanced osteogenic potential both in vitro and in vivo of human frontal derived osteoblasts compared to parietal derived osteoblasts. Furthermore, we found shared differential signaling patterns in the canonical WNT, TGF-β, BMP, and FGF pathways previously described in the mouse to govern these regional differences in osteogenic potential. Taken together, our findings unveil evolutionary conserved similarities both at functional and molecular level between the mouse and human calvarial bones, providing further support that studies employing mouse models, are suitable for translational studies to human.

Improved Understanding of the Role of Gene and Genome Duplications in Chordate Evolution With New Genome and Transcriptome Sequences

Comparative approaches to understanding chordate genomes have uncovered a significant role for gene duplications, including whole genome duplications (WGDs), giving rise to and expanding gene families. In developmental biology, gene families created and expanded by both tandem and WGDs are paramount. These genes, often involved in transcription and signalling, are candidates for underpinning major evolutionary transitions because they are particularly prone to retention and subfunctionalisation, neofunctionalisation, or specialisation following duplication. Under the subfunctionalisation model, duplication lays the foundation for the diversification of paralogues, especially in the context of gene regulation. Tandemly duplicated paralogues reside in the same regulatory environment, which may constrain them and result in a gene cluster with closely linked but subtly different expression patterns and functions. Ohnologues (WGD paralogues) often diversify by partitioning their expression domains between retained paralogues, amidst the many changes in the genome during rediploidisation, including chromosomal rearrangements and extensive gene losses. The patterns of these retentions and losses are still not fully understood, nor is the full extent of the impact of gene duplication on chordate evolution. The growing number of sequencing projects, genomic resources, transcriptomics, and improvements to genome assemblies for diverse chordates from non-model and under-sampled lineages like the coelacanth, as well as key lineages, such as amphioxus and lamprey, has allowed more informative comparisons within developmental gene families as well as revealing the extent of conserved synteny across whole genomes. This influx of data provides the tools necessary for phylogenetically informed comparative genomics, which will bring us closer to understanding the evolution of chordate body plan diversity and the changes underpinning the origin and diversification of vertebrates.

Corticotropin-Releasing Hormone (CRH) Gene Family Duplications in Lampreys Correlate With Two Early Vertebrate Genome Doublings

The ancestor of gnathostomes (jawed vertebrates) is generally considered to have undergone two rounds of whole genome duplication (WGD). The timing of these WGD events relative to the divergence of the closest relatives of the gnathostomes, the cyclostomes, has remained contentious. Lampreys and hagfishes are extant cyclostomes whose gene families can shed light on the relationship between the WGDs and the cyclostome-gnathostome divergence. Previously, we have characterized in detail the evolution of the gnathostome corticotropin-releasing hormone (CRH) family and found that its five members arose from two ancestral genes that existed before the WGDs. The two WGDs resulted, after secondary losses, in one triplet consisting of CRH1, CRH2, and UCN1, and one pair consisting of UCN2 and UCN3. All five genes exist in representatives for cartilaginous fishes, ray-finned fishes, and lobe-finned fishes. Differential losses have occurred in some lineages. We present here analyses of CRH-family members in lamprey and hagfish by comparing sequences and gene synteny with gnathostomes. We found five CRH-family genes in each of two lamprey species (Petromyzon marinus and Lethenteron camtschaticum) and two genes in a hagfish (Eptatretus burgeri). Synteny analyses show that all five lamprey CRH-family genes have similar chromosomal neighbors as the gnathostome genes. The most parsimonious explanation is that the lamprey CRH-family genes are orthologs of the five gnathostome genes and thus arose in the same chromosome duplications. This suggests that lampreys and gnathostomes share the same two WGD events and that these took place before the lamprey-gnathostome divergence.

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

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

Persistence of the ABCC6 genes and the emergence of the bony skeleton in vertebrates

The ATP-binding cassette transporter 6 (ABCC6) gene encodes a cellular transmembrane protein transporter (MRP6) that is involved in the regulation of tissue calcification in mammals. Mutations in ABCC6 are associated with human ectopic calcification disorders. To gain insight into its evolution and involvement in tissue calcification we conducted a comparative analysis of the ABCC6 gene and the related gene ABCC1 from invertebrates to vertebrates where a bony endoskeleton first evolved. Taking into consideration the role of ABCC6 in ectopic calcification of human skin we analysed the involvement of both genes in the regeneration of scales, mineralized structures that develop in fish skin. The ABCC6 gene was only found in bony vertebrate genomes and was absent from Elasmobranchs, Agnatha and from invertebrates. In teleost fish the abcc6 gene duplicated but the two genes persisted only in some teleost genomes. Six disease causing amino acid mutations in human MRP6 are a normal feature of abcc6 in fish, suggesting they do not have a deleterious effect on the protein. After scale removal the abcc6 (5 and 10 days) and abcc1 (10 days) gene expression was up-regulated relative to the intact control skin and this coincided with a time of intense scale mineralization.

Vertebrate Genome Evolution in the Light of Fish Cytogenomics and rDNAomics

To understand the cytogenomic evolution of vertebrates, we must first unravel the complex genomes of fishes, which were the first vertebrates to evolve and were ancestors to all other vertebrates. We must not forget the immense time span during which the fish genomes had to evolve. Fish cytogenomics is endowed with unique features which offer irreplaceable insights into the evolution of the vertebrate genome. Due to the general DNA base compositional homogeneity of fish genomes, fish cytogenomics is largely based on mapping DNA repeats that still represent serious obstacles in genome sequencing and assembling, even in model species. Localization of repeats on chromosomes of hundreds of fish species and populations originating from diversified environments have revealed the biological importance of this genomic fraction. Ribosomal genes (rDNA) belong to the most informative repeats and in fish, they are subject to a more relaxed regulation than in higher vertebrates. This can result in formation of a literal 'rDNAome' consisting of more than 20,000 copies with their high proportion employed in extra-coding functions. Because rDNA has high rates of transcription and recombination, it contributes to genome diversification and can form reproductive barrier. Our overall knowledge of fish cytogenomics grows rapidly by a continuously increasing number of fish genomes sequenced and by use of novel sequencing methods improving genome assembly. The recently revealed exceptional compositional heterogeneity in an ancient fish lineage (gars) sheds new light on the compositional genome evolution in vertebrates generally. We highlight the power of synergy of cytogenetics and genomics in fish cytogenomics, its potential to understand the complexity of genome evolution in vertebrates, which is also linked to clinical applications and the chromosomal backgrounds of speciation. We also summarize the current knowledge on fish cytogenomics and outline its main future avenues.

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).

An evaluation of the evolution of the gene structure of dystroglycan

Background

Dystroglycan (DG) is an adhesion receptor complex composed of two non-covalently associated subunits, transcribed from a single gene. The extracellular α-DG is highly and heterogeneously glycosylated and binds with high affinity to laminins, and the transmembrane β-DG binds intracellular dystrophin. Multiple cellular functions have been proposed for DG, notwithstanding that its role in skeletal muscle appears central as demonstrated by both primary and secondary severe muscular dystrophic phenotypes collectively known as dystroglycanopathies. We recently analysed the molecular phylogeny of the DG core protein and identified the α/β interface, transmembrane and cytoplasmic domains of β-DG as the most conserved region. It was also identified that the IG2_MAT_NU region has been independently duplicated in multiple lineages.

Results

To understand the evolution of dystroglycan in more depth, we investigated dystroglycan gene structure in 35 species representative of the phyla in which dystroglycan has been identified (i.e., all metazoan phyla except Ctenophora). The gene structure of three exons and two introns is remarkably conserved. However, additional lineage-specific introns were identified, which interrupt the coding sequence at distinct points, were identified in multiple metazoan groups, most prominently in ecdysozoans.

Conclusions

A coding DNA sequence (CDS) intron that interrupts the encoding of the IG1 domain is universally conserved and this intron is longer in gnathostomes (jawed vertebrates) than in other metazoans. Lineage-specific gain of additional introns has occurred notably in ecdysozoans, where multiple introns interrupt the large 3′ exon. More limited intron gain has also occurred in placozoa, cnidarians, urochordates and the DG paralogues of lamprey and teleost fish.

Electronic supplementary material

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

Evolutionary and Comparative Genomics to Drive Rational Drug Design, with Particular Focus on Neuropeptide Seven-Transmembrane Receptors

Seven transmembrane receptors (7TMRs), also known as G protein-coupled receptors, are popular targets of drug development, particularly 7TMR systems that are activated by peptide ligands. Although many pharmaceutical drugs have been discovered via conventional bulk analysis techniques the increasing availability of structural and evolutionary data are facilitating change to rational, targeted drug design. This article discusses the appeal of neuropeptide-7TMR systems as drug targets and provides an overview of concepts in the evolution of vertebrate genomes and gene families. Subsequently, methods that use evolutionary concepts and comparative analysis techniques to aid in gene discovery, gene function identification, and novel drug design are provided along with case study examples.

The Nuclear DNA Content and Genetic Diversity of Lampetra morii

We investigated the nuclear DNA content and genetic diversity of a river lamprey, the Korean lamprey Lampetra morii, which is distributed in the northeast of China. L. morii spends its whole life cycle in fresh water, and its adult size is relatively small (~160 mm long) compared with that of other lampreys. The haploid nuclear DNA content of L. morii is 1.618 pg (approximately 1.582 Gb) in germline cells, and there is ~15% germline DNA loss in somatic cells. These values are significantly smaller than those of Petromyzon marinus, a lamprey with a published draft genome. The chromosomes of L. morii are small and acrocentric, with a diploid modal number of 2n = 132, lower than some other lampreys. Sequence and AFLP analyses suggest that the allelic polymorphism rate (~0.14% based on examined nuclear and mitochondrial DNA sequences) of L. morii is much lower than that (~2%) of P. marinus. Phylogenetic analysis based on a mitochondrial DNA fragment confirms that L. morii belongs to the genus Lampetra, which, together with the genus Lethenteron, forms a sister group to P. marinus. These genetic background data are valuable for subsequent genetic and genomic research on L. morii.

Expression of a Novel D4 Dopamine Receptor in the Lamprey Brain. Evolutionary Considerations about Dopamine Receptors

Numerous data reported in lampreys, which belong to the phylogenetically oldest branch of vertebrates, show that the dopaminergic system was already well developed at the dawn of vertebrate evolution. The expression of dopamine in the lamprey brain is well conserved when compared to other vertebrates, and this is also true for the D2 receptor. Additionally, the key role of dopamine in the striatum, modulating the excitability in the direct and indirect pathways through the D1 and D2 receptors, has also been recently reported in these animals. The moment of divergence regarding the two whole genome duplications occurred in vertebrates suggests that additional receptors, apart from the D1 and D2 previously reported, could be present in lampreys. We used in situ hybridization to characterize the expression of a novel dopamine receptor, which we have identified as a D4 receptor according to the phylogenetic analysis. The D4 receptor shows in the sea lamprey a more restricted expression pattern than the D2 subtype, as reported in mammals. Its main expression areas are the striatum, lateral and ventral pallial sectors, several hypothalamic regions, habenula, and mesencephalic and rhombencephalic motoneurons. Some expression areas are well conserved through vertebrate evolution, as is the case of the striatum or the habenula, but the controversies regarding the D4 receptor expression in other vertebrates hampers for a complete comparison, especially in rhombencephalic regions. Our results further support that the dopaminergic system in vertebrates is well conserved and suggest that at least some functions of the D4 receptor were already present before the divergence of lampreys.

Identification of Thyroid Hormones and Functional Characterization of Thyroid Hormone Receptor in the Pacific Oyster Crassostrea gigas Provide Insight into Evolution of the Thyroid Hormone System

Thyroid hormones (THs) play important roles in development, metamorphosis, and metabolism in vertebrates. During the past century, TH functions were regarded as a synapomorphy of vertebrates. More recently, accumulating evidence has gradually convinced us that TH functions also occur in invertebrate chordates. To date, however, TH-related studies in non-chordate invertebrates have been limited. In this study, THs were qualitatively detected by two reliable methods (HPLC and LC/MS) in a well-studied molluscan species, the Pacific oyster Crassostrea gigas. Quantitative measurement of THs during the development of C. gigas showed high TH contents during embryogenesis and that oyster embryos may synthesize THs endogenously. As a first step in elucidating the TH signaling cascade, an ortholog of vertebrate TH receptor (TR), the most critical gene mediating TH effects, was cloned in C. gigas. The sequence of CgTR has conserved DNA-binding and ligand-binding domains that normally characterize these receptors. Experimental results demonstrated that CgTR can repress gene expression through binding to promoters of target genes and can interact with oyster retinoid X receptor. Moreover, CgTR mRNA expression was activated by T4 and the transcriptional activity of CgTR promoter was repressed by unliganded CgTR protein. An atypical thyroid hormone response element (CgDR5) was found in the promoter of CgTR, which was verified by electrophoretic mobility shift assay (EMSA). These results indicated that some of the CgTR function is conserved. However, the EMSA assay showed that DNA binding specificity of CgTR was different from that of the vertebrate TR and experiments with two dual-luciferase reporter systems indicated that l-thyroxine, 3,3′,5-triiodothyronine, and triiodothyroacetic acid failed to activate the transcriptional activity of CgTR. This is the first study to functionally characterize TR in mollusks. The presence of THs and the functions of CgTR in mollusks contribute to better understanding of the evolution of the TH system.

Lampreys as Diverse Model Organisms in the Genomics Era

Lampreys, one of the two surviving groups of ancient vertebrates, have become important models for study in diverse fields of biology. Lampreys (of which there are approximately 40 species) are being studied, for example, (a) to control pest sea lamprey in the North American Great Lakes and to restore declining populations of native species elsewhere; (b) in biomedical research, focusing particularly on the regenerative capability of lampreys; and (c) by developmental biologists studying the evolution of key vertebrate characters. Although a lack of genetic resources has hindered research on the mechanisms regulating many aspects of lamprey life history and development, formerly intractable questions are now amenable to investigation following the recent publication of the sea lamprey genome. Here, we provide an overview of the ways in which genomic tools are currently being deployed to tackle diverse research questions and suggest several areas that may benefit from the availability of the sea lamprey genome.

The CXC chemokine receptors of fish: Insights into CXCR evolution in the vertebrates

This article will review current knowledge on CXCR in fish, that represent three distinct vertebrate groups: Agnatha (jawless fishes), Chondrichthyes (cartilaginous fishes) and Osteichthyes (bony fishes). With the sequencing of many fish genomes, information on CXCR in these species in particular has expanded considerably. In mammals, 6 CXCRs have been described, and their homologues will be initially reviewed before considering a number of atypical CXCRs and a discussion of CXCR evolution.