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.

Functional analyses of the polycomb-group genes in sea lamprey embryos undergoing programmed DNA loss.

During early development, sea lamprey embryos undergo programmatic elimination of DNA from somatic progenitor cells in a process termed programmed genome rearrangement (PGR). Eliminated DNA eventually becomes condensed into micronuclei, which are then physically degraded and permanently lost from the cell. Previous studies indicated that many of the genes eliminated during PGR have mammalian homologs that are bound by polycomb repressive complex (PRC) in embryonic stem cells. To test whether PRC components play a role in the faithful elimination of germline-specific sequences, we used a combination of CRISPR/Cas9 and lightsheet microscopy to investigate the impact of gene knockouts on early development and the progression through stages of DNA elimination. Analysis of knockout embryos for the core PRC2 subunits EZH, SUZ12, and EED show that disruption of all three genes results in an increase in micronucleus number, altered distribution of micronuclei within embryos, and an increase in micronucleus volume in mutant embryos. While the upstream events of DNA elimination are not strongly impacted by loss of PRC2 components, this study suggests that PRC2 plays a role in the later stages of elimination related to micronucleus condensation and degradation. These findings also suggest that other genes/epigenetic pathways may work in parallel during DNA elimination to mediate chromatin structure, accessibility, and the ultimate loss of germline-specific DNA.

Expression dynamics of dehydration tolerance in the tropical plant Marchantia inflexa.

Tolerance to prolonged water deficit occurs along a continuum in plants, with dehydration tolerance (DhT) and desiccation tolerance (DT) representing some of the most extreme adaptations to water scarcity. Although DhT and DT presumably vary among individuals of a single species, this variability remains largely unstudied. Here, we characterized expression dynamics throughout a dehydration-rehydration time-course in six diverse genotypes of the dioecious liverwort Marchantia inflexa. We identified classical signatures of stress response in M. inflexa, including major changes in transcripts related to metabolism, expression of LEA and ELIP genes, and evidence of cell wall remodeling. However, we detected very little temporal synchronization of these responses across different genotypes of M. inflexa, which may be related to genotypic variation among samples, constitutive expression of dehydration-associated transcripts, the sequestration of mRNAs in ribonucleoprotein partials prior to drying, or the lower tolerance of M. inflexa relative to most bryophytes studied to date. Our characterization of intraspecific variation in expression dynamics suggests that differences in the timing of transcriptional adjustments contribute to variation among genotypes, and that developmental differences impact the relative tolerance of meristematic and differentiated tissues. This work highlights the complexity and variability of water stress tolerance, and underscores the need for comparative studies that seek to characterize variation in DT and DhT.

A chromosome-scale assembly of the axolotl genome.

The axolotl (Ambystoma mexicanum) provides critical models for studying regeneration, evolution, and development. However, its large genome (∼32 Gb) presents a formidable barrier to genetic analyses. Recent efforts have yielded genome assemblies consisting of thousands of unordered scaffolds that resolve gene structures, but do not yet permit large-scale analyses of genome structure and function. We adapted an established mapping approach to leverage dense SNP typing information and for the first time assemble the axolotl genome into 14 chromosomes. Moreover, we used fluorescence in situ hybridization to verify the structure of these 14 scaffolds and assign each to its corresponding physical chromosome. This new assembly covers 27.3 Gb and encompasses 94% of annotated gene models on chromosomal scaffolds. We show the assembly's utility by resolving genome-wide orthologies between the axolotl and other vertebrates, identifying the footprints of historical introgression events that occurred during the development of axolotl genetic stocks, and precisely mapping several phenotypes including a large deletion underlying the cardiac mutant. This chromosome-scale assembly will greatly facilitate studies of the axolotl in biological research.

Sequencing of laser captured Z and W chromosomes of the tocantins paradoxical frog (Pseudis tocantins) provides insights on repeatome and chromosomal homology.

Pseudis tocantins is the only frog species of the hylid genus Pseudis that possesses highly heteromorphic sex chromosomes. Z and W chromosomes of Ps. tocantins differ in size, morphology, position of the nucleolar organizer region (NOR) and the amount and distribution of heterochromatin. A chromosomal inversion and heterochromatin amplification on the W chromosome were previously inferred to be involved in the evolution of this sex chromosome pair. Despite these findings, knowledge related to the molecular composition of the large heterochromatic band of this W chromosome is restricted to the PcP190 satellite DNA, and no data are available regarding the gene content of either the W or the Z chromosome of Ps. tocantins. Here, we sequenced microdissected Z and W chromosomes of this species to further resolve their molecular composition. Comparative genomic analysis suggests that Ps. tocantins sex chromosomes are likely homologous to chromosomes 4 and 10 of Xenopus tropicalis. Analyses of the repetitive DNA landscape in the Z and W assemblies allowed for the identification of several transposable elements and putative satellite DNA sequences. Finally, some transposable elements from the W assembly were found to be highly diverse and divergent from elements found elsewhere in the genome, suggesting a rapid amplification of these elements on the W chromosome.

Population Genomics of New Zealand Pouched Lamprey (kanakana; piharau; Geotria australis).

Pouched lamprey (Geotria australis) or kanakana/piharau is a culturally and ecologically significant jawless fish that is distributed throughout Aotearoa New Zealand. Despite its importance, much remains unknown about historical relationships and gene flow between populations of this enigmatic species within New Zealand. To help inform management, we assembled a draft G. australis genome and completed the first comprehensive population genomics analysis of pouched lamprey within New Zealand using targeted gene sequencing (Cyt-b and COI) and restriction site-associated DNA sequencing (RADSeq) methods. Employing 16 000 genome-wide single nucleotide polymorphisms (SNPs) derived from RADSeq (n = 186) and sequence data from Cyt-b (766 bp, n = 94) and COI (589 bp, n = 20), we reveal low levels of structure across 10 sampling locations spanning the species range within New Zealand. F-statistics, outlier analyses, and STRUCTURE suggest a single panmictic population, and Mantel and EEMS tests reveal no significant isolation by distance. This implies either ongoing gene flow among populations or recent shared ancestry among New Zealand pouched lamprey. We can now use the information gained from these genetic tools to assist managers with monitoring effective population size, managing potential diseases, and conservation measures such as artificial propagation programs. We further demonstrate the general utility of these genetic tools for acquiring information about elusive species.

Large-scale variation in single nucleotide polymorphism density within the laboratory axolotl (Ambystoma mexicanum).

BACKGROUND: Recent efforts to assemble and analyze the Ambystoma mexicanum genome have dramatically improved the potential to develop molecular tools and pursue genome-wide analyses of genetic variation. RESULTS: To better resolve the distribution and origins of genetic variation with A mexicanum, we compared DNA sequence data for two laboratory A mexicanum and one A tigrinum to identify 702 million high confidence polymorphisms distributed across the 32 Gb genome. While the wild-caught A tigrinum was generally more polymorphic in a genome-wide sense, several multi-megabase regions were identified from A mexicanum genomes that were actually more polymorphic than A tigrinum. Analysis of polymorphism and repeat content reveals that these regions likely originated from the intentional hybridization of A mexicanum and A tigrinum that was used to introduce the albino mutation into laboratory stocks. CONCLUSIONS: Our findings show that axolotl genomes are variable with respect to introgressed DNA from a highly polymorphic species. It seems likely that other divergent regions will be discovered with additional sequencing of A mexicanum. This has practical implications for designing molecular probes and suggests a need to study A mexicanum phenotypic variation and genome evolution across the tiger salamander clade.

Genomic islands of divergence infer a phenotypic landscape in Pacific lamprey.

High rates of dispersal can breakdown coadapted gene complexes. However, concentrated genomic architecture (i.e., genomic islands of divergence) can suppress recombination to allow evolution of local adaptations despite high gene flow. Pacific lamprey (Entosphenus tridentatus) is a highly dispersive anadromous fish. Observed trait diversity and evidence for genetic basis of traits suggests it may be locally adapted. We addressed whether concentrated genomic architecture could influence local adaptation for Pacific lamprey. Using two new whole genome assemblies and genotypes from 7,716 single nucleotide polymorphism (SNP) loci in 518 individuals from across the species range, we identified four genomic islands of divergence (on chromosomes 01, 02, 04, and 22). We determined robust phenotype-by-genotype relationships by testing multiple traits across geographic sites. These trait associations probably explain genomic divergence across the species' range. We genotyped a subset of 302 broadly distributed SNPs in 2,145 individuals for association testing for adult body size, sexual maturity, migration distance and timing, adult swimming ability, and larval growth. Body size traits were strongly associated with SNPs on chromosomes 02 and 04. Moderate associations also implicated SNPs on chromosome 01 as being associated with variation in female maturity. Finally, we used candidate SNPs to extrapolate a heterogeneous spatiotemporal distribution of these predicted phenotypes based on independent data sets of larval and adult collections. These maturity and body size results guide future elucidation of factors driving regional optimization of these traits for fitness. Pacific lamprey is culturally important and imperiled. This research addresses biological uncertainties that challenge restoration efforts.

Comparative transcriptomics of limb regeneration: Identification of conserved expression changes among three species of Ambystoma.

Transcriptome studies are revealing the complex gene expression basis of limb regeneration in the primary salamander model - Ambystoma mexicanum (axolotl). To better understand this complexity, there is need to extend analyses to additional salamander species. Using microarray and RNA-Seq, we performed a comparative transcriptomic study using A. mexicanum and two other ambystomatid salamanders: A. andersoni, and A. maculatum. Salamanders were administered forelimb amputations and RNA was isolated and analyzed to identify 405 non-redundant genes that were commonly, differentially expressed 24 h post amputation. Many of the upregulated genes are predicted to function in wound healing and developmental processes, while many of the downregulated genes are typically expressed in muscle. The conserved transcriptional changes identified in this study provide a high-confidence dataset for identifying factors that simultaneous orchestrate wound healing and regeneration processes in response to injury, and more generally for identifying genes that are essential for salamander limb regeneration.

The African coelacanth genome provides insights into tetrapod evolution.

The discovery of a living coelacanth specimen in 1938 was remarkable, as this lineage of lobe-finned fish was thought to have become extinct 70 million years ago. The modern coelacanth looks remarkably similar to many of its ancient relatives, and its evolutionary proximity to our own fish ancestors provides a glimpse of the fish that first walked on land. Here we report the genome sequence of the African coelacanth, Latimeria chalumnae. Through a phylogenomic analysis, we conclude that the lungfish, and not the coelacanth, is the closest living relative of tetrapods. Coelacanth protein-coding genes are significantly more slowly evolving than those of tetrapods, unlike other genomic features. Analyses of changes in genes and regulatory elements during the vertebrate adaptation to land highlight genes involved in immunity, nitrogen excretion and the development of fins, tail, ear, eye, brain and olfaction. Functional assays of enhancers involved in the fin-to-limb transition and in the emergence of extra-embryonic tissues show the importance of the coelacanth genome as a blueprint for understanding tetrapod evolution.

Cellular and Molecular Features of Developmentally Programmed Genome Rearrangement in a Vertebrate (Sea Lamprey: Petromyzon marinus).

The sea lamprey (Petromyzon marinus) represents one of the few vertebrate species known to undergo large-scale programmatic elimination of genomic DNA over the course of its normal development. Programmed genome rearrangements (PGRs) result in the reproducible loss of ~20% of the genome from somatic cell lineages during early embryogenesis. Studies of PGR hold the potential to provide novel insights related to the maintenance of genome stability during the cell cycle and coordination between mechanisms responsible for the accurate distribution of chromosomes into daughter cells, yet little is known regarding the mechanistic basis or cellular context of PGR in this or any other vertebrate lineage. Here we identify epigenetic silencing events that are associated with the programmed elimination of DNA and describe the spatiotemporal dynamics of PGR during lamprey embryogenesis. In situ analyses reveal that the earliest DNA methylation (and to some extent H3K9 trimethylation) events are limited to specific extranuclear structures (micronuclei) containing eliminated DNA. During early embryogenesis a majority of micronuclei (~60%) show strong enrichment for repressive chromatin modifications (H3K9me3 and 5meC). These analyses also led to the discovery that eliminated DNA is packaged into chromatin that does not migrate with somatically retained chromosomes during anaphase, a condition that is superficially similar to lagging chromosomes observed in some cancer subtypes. Closer examination of "lagging" chromatin revealed distributions of repetitive elements, cytoskeletal contacts and chromatin contacts that provide new insights into the cellular mechanisms underlying the programmed loss of these segments. Our analyses provide additional perspective on the cellular and molecular context of PGR, identify new structures associated with elimination of DNA and reveal that PGR is completed over the course of several successive cell divisions.

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.

A linkage map for the Newt Notophthalmus viridescens: Insights in vertebrate genome and chromosome evolution.

Genetic linkage maps are fundamental resources that enable diverse genetic and genomic approaches, including quantitative trait locus (QTL) analyses and comparative studies of genome evolution. It is straightforward to build linkage maps for species that are amenable to laboratory culture and genetic crossing designs, and that have relatively small genomes and few chromosomes. It is more difficult to generate linkage maps for species that do not meet these criteria. Here, we introduce a method to rapidly build linkage maps for salamanders, which are known for their enormous genome sizes. As proof of principle, we developed a linkage map with thousands of molecular markers (N=2349) for the Eastern newt (Notophthalmus viridescens). The map contains 12 linkage groups (152.3-934.7cM), only one more than the number of chromosome pairs. Importantly, this map was generated using RNA isolated from a single wild caught female and her 28 offspring. We used the map to reveal chromosome-scale conservation of synteny among N. viridescens, A. mexicanum (Urodela), and chicken (Amniota), and to identify large conserved segments between N. viridescens and Xenopus tropicalis (Anura). We also show that met1, a major effect QTL that regulates the expression of alternate metamorphic and paedomorphic modes of development in Ambystoma, associates with a chromosomal fusion that is not found in the N. viridescens map. Our results shed new light on the ancestral amphibian karyotype and reveal specific fusion and translocation events that shaped the genomes of three amphibian model taxa. The ability to rapidly build linkage maps for large salamander genomes will enable genetic and genomic analyses within this important vertebrate group, and more generally, empower comparative studies of vertebrate biology and evolution.

Using transcriptomics to enable a plethodontid salamander (Bolitoglossa ramosi) for limb regeneration research.

BACKGROUND: Tissue regeneration is widely distributed across the tree of life. Among vertebrates, salamanders possess an exceptional ability to regenerate amputated limbs and other complex structures. Thus far, molecular insights about limb regeneration have come from a relatively limited number of species from two closely related salamander families. To gain a broader perspective on the molecular basis of limb regeneration and enhance the molecular toolkit of an emerging plethodontid salamander (Bolitoglossa ramosi), we used RNA-Seq to generate a de novo reference transcriptome and identify differentially expressed genes during limb regeneration. RESULTS: Using paired-end Illumina sequencing technology and Trinity assembly, a total of 433,809 transcripts were recovered and we obtained functional annotation for 142,926 non-redundant transcripts of the B. ramosi de novo reference transcriptome. Among the annotated transcripts, 602 genes were identified as differentially expressed during limb regeneration. This list was further processed to identify a core set of genes that exhibit conserved expression changes between B. ramosi and the Mexican axolotl (Ambystoma mexicanum), and presumably their common ancestor from approximately 180 million years ago. CONCLUSIONS: We identified genes from B. ramosi that are differentially expressed during limb regeneration, including multiple conserved protein-coding genes and possible putative species-specific genes. Comparative analyses reveal a subset of genes that show similar patterns of expression with ambystomatid species, which highlights the importance of developing comparative gene expression data for studies of limb regeneration among salamanders.

The sea lamprey meiotic map improves resolution of ancient vertebrate genome duplications.

It is generally accepted that many genes present in vertebrate genomes owe their origin to two whole-genome duplications that occurred deep in the ancestry of the vertebrate lineage. However, details regarding the timing and outcome of these duplications are not well resolved. We present high-density meiotic and comparative genomic maps for the sea lamprey (Petromyzon marinus), a representative of an ancient lineage that diverged from all other vertebrates ∼550 million years ago. Linkage analyses yielded a total of 95 linkage groups, similar to the estimated number of germline chromosomes (1n ∼ 99), spanning a total of 5570.25 cM. Comparative mapping data yield strong support for the hypothesis that a single whole-genome duplication occurred in the basal vertebrate lineage, but do not strongly support a hypothetical second event. Rather, these comparative maps reveal several evolutionarily independent segmental duplications occurring over the last 600+ million years of chordate evolution. This refined history of vertebrate genome duplication should permit more precise investigations of vertebrate evolution.

Characterization of Somatically-Eliminated Genes During Development of the Sea Lamprey (Petromyzon marinus).

The sea lamprey (Petromyzon marinus) is a basal vertebrate that undergoes developmentally programmed genome rearrangements (PGRs) during early development. These events facilitate the elimination of ∼20% of the genome from the somatic cell lineage, resulting in distinct somatic and germline genomes. Thus far only a handful of germline-specific genes have been definitively identified within the estimated 500 Mb of DNA that is deleted during PGR, although a few thousand germline-specific genes are thought to exist. To improve our understanding of the evolutionary/developmental logic of PGR, we generated computational predictions to identify candidate germline-specific genes within a new transcriptomic dataset derived from adult germline and the early embryonic stages during which PGR occurs. Follow-up validation studies identified 44 germline-specific genes and further characterized patterns of transcription and DNA loss during early embryogenesis. Expression analyses reveal that many of these genes are differentially expressed during early embryogenesis and presumably function in the early development of the germline. Ontology analyses indicate that many of these germline-specific genes play known roles in germline development, pluripotency, and oncogenesis (when misexpressed). These studies provide support for the theory that PGR serves to segregate molecular functions related to germline development/pluripotency in order to prevent their potential misexpression in somatic cells. This larger set of eliminated genes also allows us to extend the evolutionary/developmental breadth of this theory, as some deleted genes (or their gnathostome homologs) appear to be associated with the early development of somatic lineages, perhaps through the evolution of novel functions within gnathostome lineages.

Origin of amphibian and avian chromosomes by fission, fusion, and retention of ancestral chromosomes.

Amphibian genomes differ greatly in DNA content and chromosome size, morphology, and number. Investigations of this diversity are needed to identify mechanisms that have shaped the evolution of vertebrate genomes. We used comparative mapping to investigate the organization of genes in the Mexican axolotl (Ambystoma mexicanum), a species that presents relatively few chromosomes (n = 14) and a gigantic genome (>20 pg/N). We show extensive conservation of synteny between Ambystoma, chicken, and human, and a positive correlation between the length of conserved segments and genome size. Ambystoma segments are estimated to be four to 51 times longer than homologous human and chicken segments. Strikingly, genes demarking the structures of 28 chicken chromosomes are ordered among linkage groups defining the Ambystoma genome, and we show that these same chromosomal segments are also conserved in a distantly related anuran amphibian (Xenopus tropicalis). Using linkage relationships from the amphibian maps, we predict that three chicken chromosomes originated by fusion, nine to 14 originated by fission, and 12-17 evolved directly from ancestral tetrapod chromosomes. We further show that some ancestral segments were fused prior to the divergence of salamanders and anurans, while others fused independently and randomly as chromosome numbers were reduced in lineages leading to Ambystoma and Xenopus. The maintenance of gene order relationships between chromosomal segments that have greatly expanded and contracted in salamander and chicken genomes, respectively, suggests selection to maintain synteny relationships and/or extremely low rates of chromosomal rearrangement. Overall, the results demonstrate the value of data from diverse, amphibian genomes in studies of vertebrate genome evolution.

Salamander Hox clusters contain repetitive DNA and expanded non-coding regions: a typical Hox structure for non-mammalian tetrapod vertebrates?

Hox genes encode transcription factors that regulate embryonic and post-embryonic developmental processes. The expression of Hox genes is regulated in part by the tight, spatial arrangement of conserved coding and non-coding sequences. The potential for evolutionary changes in Hox cluster structure is thought to be low among vertebrates; however, recent studies of a few non-mammalian taxa suggest greater variation than originally thought. Using next generation sequencing of large genomic fragments (>100 kb) from the red spotted newt (Notophthalamus viridescens), we found that the arrangement of Hox cluster genes was conserved relative to orthologous regions from other vertebrates, but the length of introns and intergenic regions varied. In particular, the distance between hoxd13 and hoxd11 is longer in newt than orthologous regions from vertebrate species with expanded Hox clusters and is predicted to exceed the length of the entire HoxD clusters (hoxd13-hoxd4) of humans, mice, and frogs. Many repetitive DNA sequences were identified for newt Hox clusters, including an enrichment of DNA transposon-like sequences relative to non-coding genomic fragments. Our results suggest that Hox cluster expansion and transposon accumulation are common features of non-mammalian tetrapod vertebrates.

A living fossil in the genome of a living fossil: Harbinger transposons in the coelacanth genome.

Emerging data from the coelacanth genome are beginning to shed light on the origin and evolution of tetrapod genes and noncoding elements. Of particular relevance is the realization that coelacanth retains active copies of transposable elements that once served as raw material for the evolution of new functional sequences in the vertebrate lineage. Recognizing the evolutionary significance of coelacanth genome in this regard, we employed an ab initio search strategy to further classify its repetitive complement. This analysis uncovered a class of interspersed elements (Latimeria Harbinger 1-LatiHarb1) that is a major contributor to coelacanth genome structure and gene content (∼1% to 4% or the genome). Sequence analyses indicate that 1) each ∼8.7 kb LatiHarb1 element contains two coding regions, a transposase gene and a gene whose function is as yet unknown (MYB-like) and 2) copies of LatiHarb1 retain biological activity in the coelacanth genome. Functional analyses verify transcriptional and enhancer activities of LatiHarb1 in vivo and reveal transcriptional decoupling that could permit MYB-like genes to play functional roles not directly linked to transposition. Thus, LatiHarb1 represents the first known instance of a harbinger-superfamily transposon with contemporary activity in a vertebrate genome. Analyses of LatiHarb1 further corroborate the notion that exaptation of anciently active harbinger elements gave rise to at least two vertebrate genes (harbi1 and naif1) and indicate that the vertebrate gene tsnare1 also traces its ancestry to this transposon superfamily. Based on our analyses of LatiHarb1, we speculate that several functional features of harbinger elements may predispose the transposon superfamily toward recurrent exaptive evolution of cellular coding genes. In addition, these analyses further reinforce the broad utility of the coelacanth genome and other "outgroup" genomes in understanding the ancestry and evolution of vertebrate genes and genomes.

Insight from the lamprey genome: glimpsing early vertebrate development via neuroendocrine-associated genes and shared synteny of gonadotropin-releasing hormone (GnRH).

Study of the ancient lineage of jawless vertebrates is key to understanding the origins of vertebrate biology. The establishment of the neuroendocrine system with the hypothalamic-pituitary axis at its crux is of particular interest. Key neuroendocrine hormones in this system include the pivotal gonadotropin-releasing hormones (GnRHs) responsible for controlling reproduction via the pituitary. Previous data incorporating several lines of evidence showed all known vertebrate GnRHs were grouped into four paralogous lineages: GnRH1, 2, 3 and 4; with proposed evolutionary paths. Using the currently available lamprey genome assembly, we searched genes of the neuroendocrine system and summarize here the details representing the state of the current lamprey genome assembly. Additionally, we have analyzed in greater detail the evolutionary history of the GnRHs based on the information of the genomic neighborhood of the paralogs in lamprey as compared to other gnathostomes. Significantly, the current evidence suggests that two genome duplication events (both 1R and 2R) that generated the different fish and tetrapod paralogs took place before the divergence of the ancestral agnathans and gnathostome lineages. Syntenic analysis supports this evidence in that the previously-classified type IV GnRHs in lamprey (lGnRH-I and -III) share a common ancestry with GnRH2 and 3, and thus are no longer considered type IV GnRHs. Given the single amino acid difference between lGnRH-II and GnRH2 we propose that a GnRH2-like gene existed before the lamprey/gnathostome split giving rise to lGnRH-II and GnRH2. Furthermore, paralogous type 3 genes (lGnRH-I/III and GnRH3) evolved divergent structure/function in lamprey and gnathostome lineages.