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.

A lamprey neural cell type atlas illuminates the origins of the vertebrate brain.

The vertebrate brain emerged more than ~500 million years ago in common evolutionary ancestors. To systematically trace its cellular and molecular origins, we established a spatially resolved cell type atlas of the entire brain of the sea lamprey-a jawless species whose phylogenetic position affords the reconstruction of ancestral vertebrate traits-based on extensive single-cell RNA-seq and in situ sequencing data. Comparisons of this atlas to neural data from the mouse and other jawed vertebrates unveiled various shared features that enabled the reconstruction of cell types, tissue structures and gene expression programs of the ancestral vertebrate brain. However, our analyses also revealed key tissues and cell types that arose later in evolution. For example, the ancestral brain was probably devoid of cerebellar cell types and oligodendrocytes (myelinating cells); our data suggest that the latter emerged from astrocyte-like evolutionary precursors in the jawed vertebrate lineage. Altogether, our work illuminates the cellular and molecular architecture of the ancestral vertebrate brain and provides a foundation for exploring its diversification during evolution.

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.

Evolution of promoter-proximal pausing enabled a new layer of transcription control.

Promoter-proximal pausing of RNA polymerase II (Pol II) is a key regulatory step during transcription. Despite the central role of pausing in gene regulation, we do not understand the evolutionary processes that led to the emergence of Pol II pausing or its transition to a rate-limiting step actively controlled by transcription factors. Here we analyzed transcription in species across the tree of life. We found that unicellular eukaryotes display a slow acceleration of Pol II near transcription start sites. This proto-paused-like state transitioned to a longer, focused pause in derived metazoans which coincided with the evolution of new subunits in the NELF and 7SK complexes. Depletion of NELF reverts the mammalian focal pause to a proto-pause-like state and compromises transcriptional activation for a set of heat shock genes. Collectively, this work details the evolutionary history of Pol II pausing and sheds light on how new transcriptional regulatory mechanisms evolve.

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.

Early satellite cell communication creates a permissive environment for long-term muscle growth.

Using in vivo muscle stem cell (satellite cell)-specific extracellular vesicle (EV) tracking, satellite cell depletion, in vitro cell culture, and single-cell RNA sequencing, we show satellite cells communicate with other cells in skeletal muscle during mechanical overload. Early satellite cell EV communication primes the muscle milieu for proper long-term extracellular matrix (ECM) deposition and is sufficient to support sustained hypertrophy in adult mice, even in the absence of fusion to muscle fibers. Satellite cells modulate chemokine gene expression across cell types within the first few days of loading, and EV delivery of miR-206 to fibrogenic cells represses Wisp1 expression required for appropriate ECM remodeling. Late-stage communication from myogenic cells during loading is widespread but may be targeted toward endothelial cells. Satellite cells coordinate adaptation by influencing the phenotype of recipient cells, which extends our understanding of their role in muscle adaptation beyond regeneration and myonuclear donation.

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.

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.

Programmed DNA Elimination in Vertebrates.

Over the last few decades, an increasing number of vertebrate taxa have been identified that undergo programmed genome rearrangement, or programmed DNA loss, during development. In these organisms, the genome of germ cells is often reproducibly different from the genome of all other cells within the body. Although we clearly have not identified all vertebrate taxa that undergo programmed genome loss, the list of species known to undergo loss now represents ∼10% of vertebrate species, including several basally diverging lineages. Recent studies have shed new light on the targets and mechanisms of DNA loss and their association with canonical modes of DNA silencing. Ultimately, expansion of these studies into a larger collection of taxa will aid in reconstructing patterns of shared/independent ancestry of programmed DNA loss in the vertebrate lineage, as well as more recent evolutionary events that have shaped the structure and content of eliminated 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 tree frog (Boana pugnax) dataset of skin transcriptome for the identification of biomolecules with potential antimicrobial activities.

Increases in the prevalence of multiply resistant microbes have necessitated the search for new molecules with antimicrobial properties. One noteworthy avenue in this search is inspired by the presence of native antimicrobial peptides in the skin of amphibians. Having the second highest diversity of frogs worldwide, Colombian anurans represent an extensive natural reservoir that could be tapped in this search. Among this diversity, species such as Boana pugnax (the Chirique-Flusse Treefrog) are particularly notable, in that they thrive in a diversity of marginal habitats, utilize both aquatic and arboreal habitats, and are members of one of few genera that are known to mount a robust immunological response against the fungus Batrachochytrium dendrobatidis, which has decimated the population of frogs worldwide. To search for molecules with potential antimicrobial activity, we have assembled and annotated a reference transcriptome from the skin of four wild captured B. pugnax from Antioquia, Colombia. Analysis of potential antimicrobial and immunological components was performed using ontology analyses, we identified several antimicrobial chemokines with particularly strong potential for exhibiting broadscale antimicrobial activities, as well as several genes related to rapid alteration of transcriptional (KRAB zinc finger protein) and phosphorylation (MAPK) responses to exogenous stressors. We also found eight families of transmembrane transport proteins, including sodium, potassium and voltage-dependent calcium channels, which will be invaluable in future studies aimed at more precisely defining the diversity and function of cationic antimicrobial peptides with alpha-helical structures. These data highlight the utility of frogs such as Boana pugnax in the search of new antimicrobial molecules. Moreover, the molecular datasets presented here allow us to expand our knowledge of this species and illustrate the importance of preserving the vast potential of Colombian biodiversity for the identification of useful biomolecules.

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.

A de novo reference transcriptome for Bolitoglossa vallecula, an Andean mountain salamander in Colombia.

The amphibian order Caudata, contains several important model species for biological research. However, there is need to generate transcriptome data from representative species of the primary salamander families. Here we describe a de novo reference transcriptome for a terrestrial salamander, Bolitoglossa vallecula (Caudata: Plethodontidae). We employed paired-end (PE) illumina RNA sequencing to assemble a de novo reference transcriptome for B. vallecula. Assembled transcripts were compared against sequences from other vertebrate taxa to identify orthologous genes, and compared to the transcriptome of a close plethodontid relative (Bolitoglossa ramosi) to identify commonly expressed genes in the skin. This dataset should be useful to future comparative studies aimed at understanding important biological process, such as immunity, wound healing, and the production of antimicrobial compounds.

Evolutionary Analysis of the Zinc Finger and Homeoboxes Family of Proteins Identifies Multiple Conserved Domains and a Common Early Chordate Ancestor.

The Zinc Fingers and Homeoboxes (Zhx) proteins, Zhx1, Zhx2, and Zhx3, comprise a small family of proteins containing two amino-terminal C2-H2 zinc fingers and four or five carboxy-terminal homeodomains. These multiple homeodomains make Zhx proteins unusual because the majority of homeodomain-containing proteins contain a single homeodomain. Studies in cultured cells and mice suggest that Zhx proteins can function as positive or negative transcriptional regulators. Zhx2 regulates numerous hepatic genes, and all three Zhx proteins have been implicated in different cancers. Because Zhx proteins contain multiple predicted homeodomains, are associated with interesting physiological traits, and seem to be only present in the vertebrate lineage, we investigated the evolutionary history of this small family by comparing Zhx homologs from a wide range of chordates. This analysis indicates that the zinc finger motifs and homeodomains are highly similar among all Zhx proteins and also identifies additional Zhx-specific conserved regions, including a 13 amino acid amino-terminal motif that is nearly identical among all gnathostome Zhx proteins. We found single Zhx proteins in the sea lamprey (Petromyzon marinus) and in the nonvertebrate chordates sea squirt (Ciona intestinalis) and lancelet (Branchiostoma floridae); these Zhx proteins are most similar to gnathostome Zhx3. Based on our analyses, we propose that a duplication of the primordial Zhx gene gave rise to Zhx3 and the precursor to Zhx1 and Zhx2. A subsequent tandem duplication of this precursor generated Zhx1 and Zhx2 found in gnathostomes.