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.

The hagfish genome and the evolution of vertebrates.

As the only surviving lineages of jawless fishes, hagfishes and lampreys provide a critical window into early vertebrate evolution. Here, we investigate the complex history, timing, and functional role of genome-wide duplications in vertebrates in the light of a chromosome-scale genome of the brown hagfish Eptatretus atami. Using robust chromosome-scale (paralogon-based) phylogenetic methods, we confirm the monophyly of cyclostomes, document an auto-tetraploidization (1RV) that predated the origin of crown group vertebrates ~517 Mya, and establish the timing of subsequent independent duplications in the gnathostome and cyclostome lineages. Some 1RV gene duplications can be linked to key vertebrate innovations, suggesting that this early genomewide event contributed to the emergence of pan-vertebrate features such as neural crest. The hagfish karyotype is derived by numerous fusions relative to the ancestral cyclostome arrangement preserved by lampreys. These genomic changes were accompanied by the loss of genes essential for organ systems (eyes, osteoclast) that are absent in hagfish, accounting in part for the simplification of the hagfish body plan; other gene family expansions account for hagfishes' capacity to produce slime. Finally, we characterise programmed DNA elimination in somatic cells of hagfish, identifying protein-coding and repetitive elements that are deleted during development. As in lampreys, the elimination of these genes provides a mechanism for resolving genetic conflict between soma and germline by repressing germline/pluripotency functions. Reconstruction of the early genomic history of vertebrates provides a framework for further exploration of vertebrate novelties.

Leukocyte Tyrosine Kinase (Ltk) Is the Mendelian Determinant of the Axolotl Melanoid Color Variant.

The great diversity of color patterns observed among amphibians is largely explained by the differentiation of relatively few pigment cell types during development. Mexican axolotls present a variety of color phenotypes that span the continuum from leucistic to highly melanistic. The melanoid axolotl is a Mendelian variant characterized by large numbers of melanophores, proportionally fewer xanthophores, and no iridophores. Early studies of melanoid were influential in developing the single-origin hypothesis of pigment cell development, wherein it has been proposed that all three pigment cell types derive from a common progenitor cell, with pigment metabolites playing potential roles in directing the development of organelles that define different pigment cell types. Specifically, these studies identified xanthine dehydrogenase (XDH) activity as a mechanism for the permissive differentiation of melanophores at the expense of xanthophores and iridophores. We used bulked segregant RNA-Seq to screen the axolotl genome for melanoid candidate genes and identify the associated locus. Dissimilar frequencies of single-nucleotide polymorphisms were identified between pooled RNA samples of wild-type and melanoid siblings for a region on chromosome 14q. This region contains gephyrin (Gphn), an enzyme that catalyzes the synthesis of the molybdenum cofactor that is required for XDH activity, and leukocyte tyrosine kinase (Ltk), a cell surface signaling receptor that is required for iridophore differentiation in zebrafish. Wild-type Ltk crispants present similar pigment phenotypes to melanoid, strongly implicating Ltk as the melanoid locus. In concert with recent findings in zebrafish, our results support the idea of direct fate specification of pigment cells and, more generally, the single-origin hypothesis of pigment cell development.

An improved germline genome assembly for the sea lamprey Petromyzon marinus illuminates the evolution of germline-specific chromosomes.

Programmed DNA loss is a gene silencing mechanism that is employed by several vertebrate and nonvertebrate lineages, including all living jawless vertebrates and songbirds. Reconstructing the evolution of somatically eliminated (germline-specific) sequences in these species has proven challenging due to a high content of repeats and gene duplications in eliminated sequences and a corresponding lack of highly accurate and contiguous assemblies for these regions. Here, we present an improved assembly of the sea lamprey (Petromyzon marinus) genome that was generated using recently standardized methods that increase the contiguity and accuracy of vertebrate genome assemblies. This assembly resolves highly contiguous, somatically retained chromosomes and at least one germline-specific chromosome, permitting new analyses that reconstruct the timing, mode, and repercussions of recruitment of genes to the germline-specific fraction. These analyses reveal major roles of interchromosomal segmental duplication, intrachromosomal duplication, and positive selection for germline functions in the long-term evolution of germline-specific chromosomes.

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.

The giant axolotl genome uncovers the evolution, scaling, and transcriptional control of complex gene loci.

Vertebrates harbor recognizably orthologous gene complements but vary 100-fold in genome size. How chromosomal organization scales with genome expansion is unclear, and how acute changes in gene regulation, as during axolotl limb regeneration, occur in the context of a vast genome has remained a riddle. Here, we describe the chromosome-scale assembly of the giant, 32 Gb axolotl genome. Hi-C contact data revealed the scaling properties of interphase and mitotic chromosome organization. Analysis of the assembly yielded understanding of the evolution of large, syntenic multigene clusters, including the Major Histocompatibility Complex (MHC) and the functional regulatory landscape of the Fibroblast Growth Factor 8 (Axfgf8) region. The axolotl serves as a primary model for studying successful regeneration.

Characterization of axolotl lampbrush chromosomes by fluorescence in situ hybridization and immunostaining.

The lampbrush chromosomes (LBCs) in oocytes of the Mexican axolotl (Ambystoma mexicanum) were identified some time ago by their relative lengths and predicted centromeres, but they have never been associated completely with the mitotic karyotype, linkage maps or genome assembly. We identified 9 of the axolotl LBCs using RNAseq to identify actively transcribed genes and 13 BAC (bacterial artificial clone) probes containing pieces of active genes. Using read coverage analysis to find candidate centromere sequences, we developed a centromere probe that localizes to all 14 centromeres. Measurements of relative LBC arm lengths and polymerase III localization patterns enabled us to identify all LBCs. This study presents a relatively simple and reliable way to identify each axolotl LBC cytologically and to anchor chromosome-length sequences (from the axolotl genome assembly) to the physical LBCs by immunostaining and fluorescence in situ hybridization. Our data will facilitate a more detailed transcription analysis of individual LBC loops.

HDAC Inhibitor Titration of Transcription and Axolotl Tail Regeneration.

New patterns of gene expression are enacted and regulated during tissue regeneration. Histone deacetylases (HDACs) regulate gene expression by removing acetylated lysine residues from histones and proteins that function directly or indirectly in transcriptional regulation. Previously we showed that romidepsin, an FDA-approved HDAC inhibitor, potently blocks axolotl embryo tail regeneration by altering initial transcriptional responses to injury. Here, we report on the concentration-dependent effect of romidepsin on transcription and regeneration outcome, introducing an experimental and conceptual framework for investigating small molecule mechanisms of action. A range of romidepsin concentrations (0-10 µM) were administered from 0 to 6 or 0 to 12 h post amputation (HPA) and distal tail tip tissue was collected for gene expression analysis. Above a threshold concentration, romidepsin potently inhibited regeneration. Sigmoidal and biphasic transcription response curve modeling identified genes with inflection points aligning to the threshold concentration defining regenerative failure verses success. Regeneration inhibitory concentrations of romidepsin increased and decreased the expression of key genes. Genes that associate with oxidative stress, negative regulation of cell signaling, negative regulation of cell cycle progression, and cellular differentiation were increased, while genes that are typically up-regulated during appendage regeneration were decreased, including genes expressed by fibroblast-like progenitor cells. Using single-nuclei RNA-Seq at 6 HPA, we found that key genes were altered by romidepin in the same direction across multiple cell types. Our results implicate HDAC activity as a transcriptional mechanism that operates across cell types to regulate the alternative expression of genes that associate with regenerative success versus failure outcomes.

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.

Molecular signatures of sexual communication in the phlebotomine sand flies.

Phlebotomine sand flies employ an elaborate system of pheromone communication wherein males produce pheromones that attract other males to leks (thus acting as an aggregation pheromone) and females to the lekking males (sex pheromone). In addition, the type of pheromone produced varies among populations. Despite the numerous studies on sand fly chemical communication, little is known of their chemosensory genome. Chemoreceptors interact with chemicals in an organism's environment to elicit essential behaviors such as the identification of suitable mates and food sources. Thus, they play important roles during adaptation and speciation. Major chemoreceptor gene families, odorant receptors (ORs), gustatory receptors (GRs) and ionotropic receptors (IRs) together detect and discriminate the chemical landscape. Here, we annotated the chemoreceptor repertoire in the genomes of Lutzomyia longipalpis and Phlebotomus papatasi, major phlebotomine vectors in the New World and Old World, respectively. Comparison with other sequenced Diptera revealed a large and unique expansion where over 80% of the ~140 ORs belong to a single, taxonomically restricted clade. We next conducted a comprehensive analysis of the chemoreceptors in 63 L. longipalpis individuals from four different locations in Brazil representing allopatric and sympatric populations and three sex-aggregation pheromone types (chemotypes). Population structure based on single nucleotide polymorphisms (SNPs) and gene copy number in the chemoreceptors corresponded with their putative chemotypes, and corroborate previous studies that identified multiple populations. Our work provides genomic insights into the underlying behavioral evolution of sexual communication in the L. longipalpis species complex in Brazil, and highlights the importance of accounting for the ongoing speciation in central and South American Lutzomyia that could have important implications for vectorial capacity.

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.

Germline-Specific Repetitive Elements in Programmatically Eliminated Chromosomes of the Sea Lamprey (Petromyzon marinus).

The sea lamprey (Petromyzon marinus) is one of few vertebrate species known to reproducibly eliminate large fractions of its genome during normal embryonic development. This germline-specific DNA is lost in the form of large fragments, including entire chromosomes, and available evidence suggests that DNA elimination acts as a permanent silencing mechanism that prevents the somatic expression of a specific subset of "germline" genes. However, reconstruction of eliminated regions has proven to be challenging due to the complexity of the lamprey karyotype. We applied an integrative approach aimed at further characterization of the large-scale structure of eliminated segments, including: (1) in silico identification of germline-enriched repeats; (2) mapping the chromosomal location of specific repetitive sequences in germline metaphases; and (3) 3D DNA/DNA-hybridization to embryonic lagging anaphases, which permitted us to both verify the specificity of elements to physically eliminated chromosomes and characterize the subcellular organization of these elements during elimination. This approach resulted in the discovery of several repetitive elements that are found exclusively on the eliminated chromosomes, which subsequently permitted the identification of 12 individual chromosomes that are programmatically eliminated during early embryogenesis. The fidelity and specificity of these highly abundant sequences, their distinctive patterning in eliminated chromosomes, and subcellular localization in elimination anaphases suggest that these sequences might contribute to the specific targeting of chromosomes for elimination or possibly in molecular interactions that mediate their decelerated poleward movement in chromosome elimination anaphases, isolation into micronuclei and eventual degradation.

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.

Miniscule differences between sex chromosomes in the giant genome of a salamander.

In the Mexican axolotl (Ambystoma mexicanum), sex is determined by a single Mendelian factor, yet its sex chromosomes do not exhibit morphological differentiation typical of many vertebrate taxa that possess a single sex-determining locus. As sex chromosomes are theorized to differentiate rapidly, species with undifferentiated sex chromosomes provide the opportunity to reconstruct early events in sex chromosome evolution. Whole genome sequencing of 48 salamanders, targeted chromosome sequencing and in situ hybridization were used to identify the homomorphic sex chromosome that carries an A. mexicanum sex-determining factor and sequences that are present only on the W chromosome. Altogether, these sequences cover ~300 kb of validated female-specific (W chromosome) sequence, representing ~1/100,000th of the 32 Gb genome. Notably, a recent duplication of ATRX, a gene associated with mammalian sex-determining pathways, is one of few functional (non-repetitive) genes identified among these W-specific sequences. This duplicated gene (ATRW) was used to develop highly predictive markers for diagnosing sex and represents a strong candidate for a recently-acquired sex determining locus (or sexually antagonistic gene) in A. mexicanum.

Publisher Correction: The sea lamprey germline genome provides insights into programmed genome rearrangement and vertebrate evolution.

When published, this article did not initially appear open access. This error has been corrected, and the open access status of the paper is noted in all versions of the paper. Additionally, affiliation 16 denoting equal contribution was missing from author Robb Krumlauf in the PDF originally published. This error has also been corrected.

Publisher Correction: The sea lamprey germline genome provides insights into programmed genome rearrangement and vertebrate evolution.

In the version of this article initially published, the present addresses for authors Dorit Hockman and Chris Amemiya were switched. The error has been corrected in the HTML and PDF versions of the article.

The sea lamprey germline genome provides insights into programmed genome rearrangement and vertebrate evolution.

The sea lamprey (Petromyzon marinus) serves as a comparative model for reconstructing vertebrate evolution. To enable more informed analyses, we developed a new assembly of the lamprey germline genome that integrates several complementary data sets. Analysis of this highly contiguous (chromosome-scale) assembly shows that both chromosomal and whole-genome duplications have played significant roles in the evolution of ancestral vertebrate and lamprey genomes, including chromosomes that carry the six lamprey HOX clusters. The assembly also contains several hundred genes that are reproducibly eliminated from somatic cells during early development in lamprey. Comparative analyses show that gnathostome (mouse) homologs of these genes are frequently marked by polycomb repressive complexes (PRCs) in embryonic stem cells, suggesting overlaps in the regulatory logic of somatic DNA elimination and bivalent states that are regulated by early embryonic PRCs. This new assembly will enhance diverse studies that are informed by lampreys' unique biology and evolutionary/comparative perspective.

Identification of Mutant Genes and Introgressed Tiger Salamander DNA in the Laboratory Axolotl, Ambystoma mexicanum.

The molecular genetic toolkit of the Mexican axolotl, a classic model organism, has matured to the point where it is now possible to identify genes for mutant phenotypes. We used a positional cloning-candidate gene approach to identify molecular bases for two historic axolotl pigment phenotypes: white and albino. White (d/d) mutants have defects in pigment cell morphogenesis and differentiation, whereas albino (a/a) mutants lack melanin. We identified in white mutants a transcriptional defect in endothelin 3 (edn3), encoding a peptide factor that promotes pigment cell migration and differentiation in other vertebrates. Transgenic restoration of Edn3 expression rescued the homozygous white mutant phenotype. We mapped the albino locus to tyrosinase (tyr) and identified polymorphisms shared between the albino allele (tyr a ) and tyr alleles in a Minnesota population of tiger salamanders from which the albino trait was introgressed. tyr a has a 142 bp deletion and similar engineered alleles recapitulated the albino phenotype. Finally, we show that historical introgression of tyr a significantly altered genomic composition of the laboratory axolotl, yielding a distinct, hybrid strain of ambystomatid salamander. Our results demonstrate the feasibility of identifying genes for traits in the laboratory Mexican axolotl.

Initial characterization of the large genome of the salamander Ambystoma mexicanum using shotgun and laser capture chromosome sequencing.

Vertebrates exhibit substantial diversity in genome size, and some of the largest genomes exist in species that uniquely inform diverse areas of basic and biomedical research. For example, the salamander Ambystoma mexicanum (the Mexican axolotl) is a model organism for studies of regeneration, development and genome evolution, yet its genome is ~10× larger than the human genome. As part of a hierarchical approach toward improving genome resources for the species, we generated 600 Gb of shotgun sequence data and developed methods for sequencing individual laser-captured chromosomes. Based on these data, we estimate that the A. mexicanum genome is ~32 Gb. Notably, as much as 19 Gb of the A. mexicanum genome can potentially be considered single copy, which presumably reflects the evolutionary diversification of mobile elements that accumulated during an ancient episode of genome expansion. Chromosome-targeted sequencing permitted the development of assemblies within the constraints of modern computational platforms, allowed us to place 2062 genes on the two smallest A. mexicanum chromosomes and resolves key events in the history of vertebrate genome evolution. Our analyses show that the capture and sequencing of individual chromosomes is likely to provide valuable information for the systematic sequencing, assembly and scaffolding of large genomes.