Hellbender genome sequences shed light on genomic expansion at the base of crown salamanders

Spiralian genomics and the evolution of animal genome architecture

Recent developments in sequencing technologies have greatly improved our knowledge of phylogenetic relationships and genomic architectures throughout the tree of life. Spiralia, a diverse clade within Protostomia, is essential for understanding the evolutionary history of parasitism, gene conversion, nervous systems and animal body plans. In this review, we focus on the current hypotheses of spiralian phylogeny and investigate the impact of long-read sequencing on the quality of genome assemblies. We examine chromosome-level assemblies to highlight key genomic features that have driven spiralian evolution, including karyotype, synteny and the Hox gene organization. In addition, we show how chromosome rearrangement has influenced spiralian genomic structures. Although spiralian genomes have undergone substantial changes, they exhibit both conserved and lineage-specific features. We recommend increasing sequencing efforts and expanding functional genomics research to deepen insights into spiralian biology.

Transposon waves at the water-to-land transition

The major transitions in vertebrate evolution are associated with significant genomic reorganizations. In contrast to the evolutionary processes that occurred at the origin of vertebrates or prior to the radiation of teleost fishes, no whole-genome duplication events occurred during the water-to-land transition, and it remains an open question how did genome dynamics contribute to this prominent evolutionary event. Indeed, the recent sequencing of sarcopterygian and amphibian genomes has revealed that the extant lineages immediately preceding and succeeding this transition harbor an exceptional number of transposable elements and it is tempting to speculate that these sequences might have catalyzed the adaptations that enabled vertebrates to venture into land. Here, we review the genome dynamics associated with the major transitions in vertebrate evolution and discuss how the highly repetitive genomic landscapes revealed by recent efforts to characterize the genomes of amphibians and sarcopterygians argue for turbulent genome dynamics occurring before the water-to-land transition and possibly enabling it.

Metamorphosis Imposes Variable Constraints on Genome Expansion through Effects on Development

Genome size varies ∼100,000-fold across eukaryotes and has long been hypothesized to be influenced by metamorphosis in animals. Transposable element accumulation has been identified as a major driver of increase, but the nature of constraints limiting the size of genomes has remained unclear, even as traits such as cell size and rate of development co-vary strongly with genome size. Salamanders, which possess diverse metamorphic and non-metamorphic life histories, join the lungfish in having the largest vertebrate genomes-3 to 40 times that of humans-as well as the largest range of variation in genome size. We tested 13 biologically-inspired hypotheses exploring how the form of metamorphosis imposes varying constraints on genome expansion in a broadly representative phylogeny containing 118 species of salamanders. We show that metamorphosis during which animals undergo the most extensive and synchronous remodeling imposes the most severe constraint against genome expansion, with the severity of constraint decreasing with reduced extent and synchronicity of remodeling. More generally, our work demonstrates the potential for broader interpretation of phylogenetic comparative analysis in exploring the balance of multiple evolutionary pressures shaping phenotypic evolution.

Transposable element and host silencing activity in gigantic genomes

Transposable elements (TEs) and the silencing machinery of their hosts are engaged in a germline arms-race dynamic that shapes TE accumulation and, therefore, genome size. In animal species with extremely large genomes (>10 Gb), TE accumulation has been pushed to the extreme, prompting the question of whether TE silencing also deviates from typical conditions. To address this question, we characterize TE silencing via two pathways-the piRNA pathway and KRAB-ZFP transcriptional repression-in the male and female gonads of Ranodon sibiricus, a salamander species with a ∼21 Gb genome. We quantify 1) genomic TE diversity, 2) TE expression, and 3) small RNA expression and find a significant relationship between the expression of piRNAs and TEs they target for silencing in both ovaries and testes. We also quantified TE silencing pathway gene expression in R. sibiricus and 14 other vertebrates with genome sizes ranging from 1 to 130 Gb and find no association between pathway expression and genome size. Taken together, our results reveal that the gigantic R. sibiricus genome includes at least 19 putatively active TE superfamilies, all of which are targeted by the piRNA pathway in proportion to their expression levels, suggesting comprehensive piRNA-mediated silencing. Testes have higher TE expression than ovaries, suggesting that they may contribute more to the species' high genomic TE load. We posit that apparently conflicting interpretations of TE silencing and genomic gigantism in the literature, as well as the absence of a correlation between TE silencing pathway gene expression and genome size, can be reconciled by considering whether the TE community or the host is currently "on the attack" in the arms race dynamic.

Now that We Got There, What Next?

As seen in the protocols in this book, the opportunities to pursue work at the cellular and molecular work in salamanders have considerably broadened over the last years. The availability of genomic information and genome editing, and the possibility to image tissues live and other methods enhance the spectrum of biological questions accessible to all researchers. Here I provide a personal perspective on what I consider exciting future questions open for investigation.

Does genome size increase with water depth in marine fishes?

A growing body of research suggests that genome size in animals can be affected by ecological factors. Half a century ago, Ebeling et al. proposed that genome size increases with depth in some teleost fish groups and discussed a number of biological mechanisms that may explain this pattern (e.g., passive accumulation, adaptive acclimation). Using phylogenetic comparative approaches, we revisit this hypothesis based on genome size and ecological data from up to 708 marine fish species in combination with a set of large-scale phylogenies, including a newly inferred tree. We also conduct modeling approaches of trait evolution and implement a variety of regression analyses to assess the relationship between genome size and depth. Our reanalysis of Ebeling et al.'s dataset shows a weak association between these variables, but the overall pattern in their data is driven by a single clade. Although new analyses based on our "all-species" dataset resulted in positive correlations, providing some evidence that genome size evolves as a function of depth, only one subclade consistently yielded statistically significant correlations. By contrast, negative correlations are rare and nonsignificant. All in all, we find modest evidence for an increase in genome size along the depth axis in marine fishes. We discuss some mechanistic explanations for the observed trends.

Transposable Element Diversity Remains High in Gigantic Genomes

Transposable elements (TEs) are repetitive sequences of DNA that replicate and proliferate throughout genomes. Taken together, all the TEs in a genome form a diverse community of sequences, which can be studied to draw conclusions about genome evolution. TE diversity can be measured using models for ecological community diversity that consider species richness and evenness. Several models predict TE diversity decreasing as genomes expand because of selection against ectopic recombination and/or competition among TEs to garner host replicative machinery and evade host silencing mechanisms. Salamanders have some of the largest vertebrate genomes and highest TE loads. Salamanders of the genus Plethodon, in particular, have genomes that range in size from 20 to 70 Gb. Here, we use Oxford Nanopore sequencing to generate low-coverage genomic sequences for four species of Plethodon that encompass two independent genome expansion events, one in the eastern clade (Plethodon cinereus, 29.3 Gb vs. Plethodon glutinosus, 38.9 Gb) and one in the western clade (Plethodon vehiculum, 46.4 Gb vs Plethodon idahoensis, 67.0 Gb). We classified the TEs in these genomes and found > 40 TE superfamilies, accounting for 22-27% of the genomes. We calculated Simpson's and Shannon's diversity indices to quantify overall TE diversity. In both pairwise comparisons, the diversity index values for the smaller and larger genome were almost identical. This result indicates that, when genomes reach extremely large sizes, they maintain high levels of TE diversity at the superfamily level, in contrast to predictions made by previous studies on smaller genomes.

Gigantic genomes of salamanders indicate that body temperature, not genome size, is the driver of global methylation and 5-methylcytosine deamination in vertebrates

Transposable elements (TEs) are sequences that replicate and move throughout genomes, and they can be silenced through methylation of cytosines at CpG dinucleotides. TE abundance contributes to genome size, but TE silencing variation across genomes of different sizes remains underexplored. Salamanders include most of the largest C-values - 9 to 120 Gb. We measured CpG methylation levels in salamanders with genomes ranging from 2N = ∼58 Gb to 4N = ∼116 Gb. We compared these levels to results from endo- and ectothermic vertebrates with more typical genomes. Salamander methylation levels are approximately 90%, higher than all endotherms. However, salamander methylation does not differ from other ectotherms, despite an approximately 100-fold difference in nuclear DNA content. Because methylation affects the nucleotide compositional landscape through 5-methylcytosine deamination to thymine, we quantified salamander CpG dinucleotide levels and compared them to other vertebrates. Salamanders and other ectotherms have comparable CpG levels, and ectotherm levels are higher than endotherms. These data show no shift in global methylation at the base of salamanders, despite a dramatic increase in TE load and genome size. This result is reconcilable with previous studies that considered endothermy and ectothermy, which may be more important drivers of methylation in vertebrates than genome size.

Taming, Domestication and Exaptation: Trajectories of Transposable Elements in Genomes

During evolution, several types of sequences pass through genomes. Along with mutations and internal genetic tinkering, they are a useful source of genetic variability for adaptation and evolution. Most of these sequences are acquired by horizontal transfers (HT), but some of them may come from the genomes themselves. If they are not lost or eliminated quickly, they can be tamed, domesticated, or even exapted. Each of these processes results from a series of events, depending on the interactions between these sequences and the host genomes, but also on environmental constraints, through their impact on individuals or population fitness. After a brief reminder of the characteristics of each of these states (taming, domestication, exaptation), the evolutionary trajectories of these new or acquired sequences will be presented and discussed, emphasizing that they are not totally independent insofar as the first can constitute a step towards the second, and the second is another step towards the third.

From Cave Dragons to Genomics: Advancements in the Study of Subterranean Tetrapods

Throughout most of the kingdom Animalia, evolutionary transitions from surface life to a life permanently bound to caves and other subterranean habitats have occurred innumerous times. Not so in tetrapods, where a mere 14 cave-obligate species—all plethodontid and proteid salamanders—are known. We discuss why cave tetrapods are so exceptional and why only salamanders have made the transition. Their evolution follows predictable and convergent, albeit independent pathways. Among the many known changes associated with transitions to subterranean life, eye degeneration, starvation resistance, and longevity are especially relevant to human biomedical research. Recently, sequences of salamander genomes have become available opening up genomic research for cave tetrapods. We discuss new genomic methods that can spur our understanding of the evolutionary mechanisms behind convergent phenotypic change, the relative roles of selective and neutral evolution, cryptic species diversity, and data relevant for conservation such as effective population size and demography.

Investigation of the activity of transposable elements and genes involved in their silencing in the newt Cynops orientalis, a species with a giant genome

Caudata is an order of amphibians with great variation in genome size, which can reach enormous dimensions in salamanders. In this work, we analysed the activity of transposable elements (TEs) in the transcriptomes obtained from female and male gonads of the Chinese fire-bellied newt, Cynops orientalis, a species with a genome about 12-fold larger than the human genome. We also compared these data with genomes of two basal sarcopterygians, coelacanth and lungfish. In the newt our findings highlighted a major impact of non-LTR retroelements and a greater total TE activity compared to the lungfish Protopterus annectens, an organism also characterized by a giant genome. This difference in TE activity might be due to the presence of young copies in newt in agreement also with the increase in the genome size, an event that occurred independently and later than lungfish. Moreover, the activity of 33 target genes encoding proteins involved in the TE host silencing mechanisms, such as Ago/Piwi and NuRD complex, was evaluated and compared between the three species analysed. These data revealed high transcriptional levels of the target genes in both newt and lungfish and confirmed the activity of NuRD complex genes in adults. Finally, phylogenetic analyses performed on PRDM9 and TRIM28 allowed increasing knowledge about the evolution of these two key genes of the NuRD complex silencing mechanism in vertebrates. Our results confirmed that the gigantism of the newt genomes may be attributed to the activity and accumulation of TEs.

Cranial skeletogenesis of one of the largest amphibians, Andrias japonicus, provides insight into ontogenetic adaptations for feeding in salamanders

Skeletal development in salamanders is greatly influenced by their complex life history. Although the relationship between skull development and life history has been investigated in a variety of salamander species, few studies have compared skull development among lineages. Here, we describe the cranial skeletogenesis of the Japanese giant salamander, Andrias japonicus (family Cryptobranchidae), for the first time. To conduct interspecific comparisons, we also describe the cranial skeletogenesis of Hynobius nebulosus (Hynobiidae), Pleurodeles waltl (Salamandridae) and Ambystoma mexicanum (Ambystomatidae). Our comparative analyses reveal interspecific differences in cranial skeletogenesis, particularly in the timing of the onset of ossification, the direction of ossification (e.g. from anterior to posterior) and the ossification pattern of the hyobranchial skeleton. The early onset of jaw ossification, posterior extension of the maxilla and posterolateral inclination of the squamosal in An. japonicus might be adaptations that allow An. japonicus to gulp feed, primarily using the suction force generated by opening the mouth wide, while also supplementing that force by depressing the hyobranchial skeleton. Multivariate regression of skull shape on log10-transformed centroid size revealed that the unique skull shape of adult and subadult An. japonicus forms through hypermorphosis.

A bird-like genome from a frog: Mechanisms of genome size reduction in the ornate burrowing frog, Platyplectrum ornatum

The diversity of genome sizes across the tree of life is of key interest in evolutionary biology. Various correlates of variation in genome size, such as accumulation of transposable elements (TEs) or rate of DNA gain and loss, are well known, but the underlying molecular mechanisms driving or constraining genome size are poorly understood. Here, we study one of the smallest genomes among frogs characterized thus far, that of the ornate burrowing frog (Platyplectrum ornatum) from Australia, and compare it to other published frog and vertebrate genomes to examine the forces driving reduction in genome size. At ∼1.06 gigabases (Gb), the P. ornatum genome is like that of birds, revealing four major mechanisms underlying TE dynamics: reduced abundance of all major classes of TEs; increased net deletion bias in TEs; drastic reduction in intron lengths; and expansion via gene duplication of the repertoire of TE-suppressing Piwi genes, accompanied by increased expression of Piwi-interacting RNA (piRNA)-based TE-silencing pathway genes in germline cells. Transcriptomes from multiple tissues in both sexes corroborate these results and provide insight into sex-differentiation pathways in Platyplectrum Genome skimming of two closely related frog species (Lechriodus fletcheri and Limnodynastes fletcheri) confirms a reduction in TEs as a major driver of genome reduction in Platyplectrum and supports a macroevolutionary scenario of small genome size in frogs driven by convergence in life history, especially rapid tadpole development and tadpole diet. The P. ornatum genome offers a model for future comparative studies on mechanisms of genome size reduction in amphibians and vertebrates generally.

Gigantic Genomes Provide Empirical Tests of Transposable Element Dynamics Models

Transposable elements (TEs) are a major determinant of eukaryotic genome size. The collective properties of a genomic TE community reveal the history of TE/host evolutionary dynamics and impact present-day host structure and function, from genome to organism levels. In rare cases, TE community/genome size has greatly expanded in animals, associated with increased cell size and changes to anatomy and physiology. Here, we characterize the TE landscape of the genome and transcriptome in an amphibian with a giant genome — the caecilianIchthyophis bannanicus, which we show has a genome size of 12.2 Gb. Amphibians are an important model system because the clade includes independent cases of genomic gigantism. The I. bannanicus genome differs compositionally from other giant amphibian genomes, but shares a low rate of ectopic recombination-mediated deletion. We examine TE activity using expression and divergence plots; TEs account for 15% of somatic transcription, and most superfamilies appear active. We quantify TE diversity in the caecilian, as well as other vertebrates with a range of genome sizes, using diversity indices commonly applied in community ecology. We synthesize previous models that integrate TE abundance, diversity, and activity, and test whether the caecilian meets model predictions for genomes with high TE abundance. We propose thorough, consistent characterization of TEs to strengthen future comparative analyses. Such analyses will ultimately be required to reveal whether the divergent TE assemblages found across convergent gigantic genomes reflect fundamental shared features of TE/host genome evolutionary dynamics.

A Field Guide to Eukaryotic Transposable Elements

Transposable elements (TEs) are mobile DNA sequences that propagate within genomes. Through diverse invasion strategies, TEs have come to occupy a substantial fraction of nearly all eukaryotic genomes, and they represent a major source of genetic variation and novelty. Here we review the defining features of each major group of eukaryotic TEs and explore their evolutionary origins and relationships. We discuss how the unique biology of different TEs influences their propagation and distribution within and across genomes. Environmental and genetic factors acting at the level of the host species further modulate the activity, diversification, and fate of TEs, producing the dramatic variation in TE content observed across eukaryotes. We argue that cataloging TE diversity and dissecting the idiosyncratic behavior of individual elements are crucial to expanding our comprehension of their impact on the biology of genomes and the evolution of species.

Cortisol is the predominant glucocorticoid in the giant paedomorphic hellbender salamander (Cryptobranchus alleganiensis)

Corticosterone is widely regarded to be the predominant glucocorticoid produced in amphibians. However, we recently described unusually low baseline and stress-induced corticosterone profiles in eastern hellbenders (Cryptobranchus alleganiensis alleganiensis), a giant, fully aquatic salamander. Here, we hypothesized that hellbenders might also produce cortisol, the predominant glucocorticoid used by fishes and non-rodent mammals. To test our hypothesis, we collected plasma samples in two field experiments and analyzed them using multiple analytical techniques to determine how plasma concentrations of cortisol and corticosterone co-varied after 1) physical restraint and 2) injection with adrenocorticotropic hormone (ACTH), the pituitary hormone responsible for triggering the release of glucocorticoids from amphibian interrenal glands. Using liquid chromatography-mass spectrometry, we found that baseline and restraint-induced plasma concentrations of cortisol were more than five times those of corticosterone. We then demonstrated that plasma concentrations of both glucocorticoids increased in response to ACTH in a dose-dependent manner, but cortisol concentrations were consistently higher (up to 10-fold) than corticosterone. Cortisol and corticosterone concentrations were not correlated with one another at basal or induced conditions. The extremely low plasma concentrations of corticosterone in hellbenders suggests that corticosterone could simply be a byproduct of cortisol production, and raises questions as to whether corticosterone has any distinct physiological function in hellbenders. Our results indicate that hellbenders produce cortisol as their predominant glucocorticoid, supporting a small and inconclusive body of literature indicating that some other amphibians may produce appreciable quantities of cortisol. We hypothesize that the use of cortisol by hellbenders could be an adaptation to their fully aquatic life history due to cortisol's ability to fulfill both mineralocorticoid and glucocorticoid functions, similar to its functions in fishes. Given the large number of amphibian species that are fully aquatic or have aquatic life stages, we suggest that the broadly held assumption that corticosterone is the predominant glucocorticoid in all amphibians requires further scrutiny. Ultimately, multi-species tests of this assumption will reveal the ecological factors that influenced the evolution of endocrine adaptations among amphibian lineages, and may provide insight into convergent evolution of endocrine traits in paedomorphic species.

Genomic Data Reveal Conserved Female Heterogamety in Giant Salamanders with Gigantic Nuclear Genomes

Systems of genetic sex determination and the homology of sex chromosomes in different taxa vary greatly across vertebrates. Much progress remains to be made in understanding systems of genetic sex determination in non-model organisms, especially those with homomorphic sex chromosomes and/or large genomes. We used reduced representation genome sequencing to investigate genetic sex determination systems in the salamander family Cryptobranchidae (genera Cryptobranchus and Andrias), which typifies both of these inherent difficulties. We tested hypotheses of male- or female-heterogamety by sequencing hundreds of thousands of anonymous genomic regions in a panel of known-sex cryptobranchids and characterized patterns of presence/absence, inferred zygosity, and depth of coverage to identify sex-linked regions of these 56 gigabase genomes. Our results strongly support the hypothesis that all cryptobranchid species possess homologous systems of female heterogamety, despite maintenance of homomorphic sex chromosomes over nearly 60 million years. Additionally, we report a robust, non-invasive genetic assay for sex diagnosis in Cryptobranchus and Andrias which may have great utility for conservation efforts with these endangered salamanders. Co-amplification of these W-linked markers in both cryptobranchid genera provides evidence for long-term sex chromosome stasis in one of the most divergent salamander lineages. These findings inform hypotheses about the ancestral mode of sex determination in salamanders, but suggest that comparative data from other salamander families are needed. Our results further demonstrate that massive genomes are not necessarily a barrier to effective genome-wide sequencing and that the resulting data can be highly informative about sex determination systems in taxa with homomorphic sex chromosomes.

Ecological constraints associated with genome size across salamander lineages

Salamanders have some of the largest, and most variable, genome sizes among the vertebrates. Larger genomes have been associated with larger cell sizes, lower metabolic rates, and longer embryonic and larval durations in many different taxonomic groups. These life-history traits are often important for dictating fitness under different environmental conditions, suggesting that a species' genome size may have the potential to constrain its ecological distribution. We test how genome size varies with the ephemerality of larval habitat across the salamanders, predicting that species with larger genomes will be constrained to more permanent habitats that permit slower development, while species with smaller genomes will be more broadly distributed across the gradient of habitat ephemerality. We found that salamanders with larger genomes are almost exclusively associated with permanent aquatic habitats. In addition, the evolutionary transition rate between permanent and ephemeral larval habitats is much higher in salamander lineages with smaller genome sizes. These patterns suggest that genome size may act as an evolutionary constraint on the ecological habitats of salamanders, restricting those species with large genomes and slower development to habitats with permanent sources of water.

Model systems for regeneration: salamanders

Salamanders have been hailed as champions of regeneration, exhibiting a remarkable ability to regrow tissues, organs and even whole body parts, e.g. their limbs. As such, salamanders have provided key insights into the mechanisms by which cells, tissues and organs sense and regenerate missing or damaged parts. In this Primer, we cover the evolutionary context in which salamanders emerged. We outline the varieties of mechanisms deployed during salamander regeneration, and discuss how these mechanisms are currently being explored and how they have advanced our understanding of animal regeneration. We also present arguments about why it is important to study closely related species in regeneration research.

The Evolution of Small-RNA-Mediated Silencing of an Invading Transposable Element

Transposable elements (TEs) are genomic parasites that impose fitness costs on their hosts by producing deleterious mutations and disrupting gametogenesis. Host genomes avoid these costs by regulating TE activity, particularly in germline cells where new insertions are heritable and TEs are exceptionally active. However, the capacity of different TE-associated fitness costs to select for repression in the host, and the role of selection in the evolution of TE regulation more generally remain controversial. In this study, we use forward, individual-based simulations to examine the evolution of small-RNA-mediated TE regulation, a conserved mechanism for TE repression that is employed by both prokaryotes and eukaryotes. To design and parameterize a biologically realistic model, we drew on an extensive survey of empirical studies of the transposition and regulation of P-element DNA transposons in Drosophila melanogaster. We observed that even under conservative assumptions, where small-RNA-mediated regulation reduces transposition only, repression evolves rapidly and adaptively after the genome is invaded by a new TE in simulated populations. We further show that the spread of repressor alleles through simulated populations is greatly enhanced by two additional TE-imposed fitness costs: dysgenic sterility and ectopic recombination. Finally, we demonstrate that the adaptive mutation rate to repression is a critical parameter that influences both the evolutionary trajectory of host repression and the associated proliferation of TEs after invasion in simulated populations. Our findings suggest that adaptive evolution of TE regulation may be stronger and more prevalent than previously appreciated, and provide a framework for interpreting empirical data.

Macroevolutionary shift in the size of amphibian genomes and the role of life history and climate

The evolution and great diversity of genome size has been of long-standing interest to biologists, but has seldom been investigated on a broad phylogenetic scale. Here we present a comparative quantitative analysis of factors shaping genome size evolution in amphibians, the extant class of vertebrates with the largest variation in genome size. We find that amphibian genomes have undergone saltations in size, although these are rare and the evolutionary history of genome size in amphibians has otherwise been one of gradual, time-dependent variation (that is, Brownian motion). This macroevolutionary homogeneity is remarkable given the evolutionary and ecological diversity of most other aspects of the natural history of amphibians. Contrary to previous claims, we find no evidence for associations between life cycle complexity and genome size despite the high diversity of reproductive modes and the multiple events of independent evolution of divergent life cycles in the group. Climate (temperature and humidity) affects genome size indirectly, at least in frogs, as a consequence of its effect on premetamorphic developmental period, although directionality of the relationship between developmental period and genome size is not unequivocal.

Influence of immunogenetics, sex and body condition on the cutaneous microbial communities of two giant salamanders

The complex association between hosts and microbial symbionts requires the implementation of multiple approaches to evaluate variation in host physiology. Within amphibians, heterogeneity in immunogenetic traits and cutaneous microbiota is associated with variation in disease resistance. Ozark (Cryptobranchus alleganiensis bishopi) and eastern hellbenders (C. a. alleganiensis) provide a model system to assess variation in host traits and microbial communities. Ozark hellbenders have experienced declines throughout their range, are federally endangered and experience wound retardation that is absent in the eastern subspecies. Previous microbial investigations indicate differentiation in the composition of the skin microbiota of both hellbender subspecies, but it is not clear whether these patterns are concurrent with diversity in the major histocompatibility complex (MHC) genes. We characterized the MHC IIB and the skin microbiota of hellbenders in Missouri, where both subspecies co-occur though not sympatric. We compared the microbiota composition and MHC diversity between both subspecies and investigated whether individual-level MHC diversity, sex and body condition were associated with microbiota composition. Overall, MHC IIB diversity was lower in Ozark hellbenders compared to the eastern subspecies. Multivariate statistical comparisons identified microbiota differentiation between Ozark and eastern hellbenders. MHC IIB allele presence/absence, allele divergence, body composition and sex defined grouping of hellbender microbiotas within populations. Differentiation of the cutaneous microbiotas and MHC IIB genes between eastern and Ozark hellbenders suggests that differences exist in immunity between the two subspecies. This study demonstrates how simultaneous assessments of host genetic traits and microbiotas can inform patterns of microbial community structure in natural systems.

The axolotl genome and the evolution of key tissue formation regulators

Salamanders serve as important tetrapod models for developmental, regeneration and evolutionary studies. An extensive molecular toolkit makes the Mexican axolotl (Ambystoma mexicanum) a key representative salamander for molecular investigations. Here we report the sequencing and assembly of the 32-gigabase-pair axolotl genome using an approach that combined long-read sequencing, optical mapping and development of a new genome assembler (MARVEL). We observed a size expansion of introns and intergenic regions, largely attributable to multiplication of long terminal repeat retroelements. We provide evidence that intron size in developmental genes is under constraint and that species-restricted genes may contribute to limb regeneration. The axolotl genome assembly does not contain the essential developmental gene Pax3. However, mutation of the axolotl Pax3 paralogue Pax7 resulted in an axolotl phenotype that was similar to those seen in Pax3^-/- and Pax7^-/- mutant mice. The axolotl genome provides a rich biological resource for developmental and evolutionary studies.

Virtual Genome Walking across the 32 Gb Ambystoma mexicanum genome; assembling gene models and intronic sequence

Large repeat rich genomes present challenges for assembly using short read technologies. The 32 Gb axolotl genome is estimated to contain ~19 Gb of repetitive DNA making an assembly from short reads alone effectively impossible. Indeed, this model species has been sequenced to 20× coverage but the reads could not be conventionally assembled. Using an alternative strategy, we have assembled subsets of these reads into scaffolds describing over 19,000 gene models. We call this method Virtual Genome Walking as it locally assembles whole genome reads based on a reference transcriptome, identifying exons and iteratively extending them into surrounding genomic sequence. These assemblies are then linked and refined to generate gene models including upstream and downstream genomic, and intronic, sequence. Our assemblies are validated by comparison with previously published axolotl bacterial artificial chromosome (BAC) sequences. Our analyses of axolotl intron length, intron-exon structure, repeat content and synteny provide novel insights into the genic structure of this model species. This resource will enable new experimental approaches in axolotl, such as ChIP-Seq and CRISPR and aid in future whole genome sequencing efforts. The assembled sequences and annotations presented here are freely available for download from https://tinyurl.com/y8gydc6n . The software pipeline is available from https://github.com/LooseLab/iterassemble .

Targeted identification of TE insertions in a Drosophila genome through hemi-specific PCR

BACKGROUND

Transposable elements (TEs) are major components of eukaryotic genomes and drivers of genome evolution, producing intraspecific polymorphism and interspecific differences through mobilization and non-homologous recombination. TE insertion sites are often highly variable within species, creating a need for targeted genome re-sequencing (TGS) methods to identify TE insertion sites.

METHODS

We present a hemi-specific PCR approach for TGS of P-elements in Drosophila genomes on the Illumina platform. We also present a computational framework for identifying new insertions from TGS reads. Finally, we describe a new method for estimating the frequency of TE insertions from WGS data, which is based precise insertion sites provided by TGS annotations.

RESULTS

By comparing our results to TE annotations based on whole genome re-sequencing (WGS) data for the same Drosophilamelanogaster strain, we demonstrate that TGS is powerful for identifying true insertions, even in repeat-rich heterochromatic regions. We also demonstrate that TGS offers enhanced annotation of precise insertion sites, which facilitates estimation of TE insertion frequency.

CONCLUSIONS

TGS by hemi-specific PCR is a powerful approach for identifying TE insertions of particular TE families in species with a high-quality reference genome, at greatly reduced cost as compared to WGS. It may therefore be ideal for population genomic studies of particular TE families. Additionally, TGS and WGS can be used as complementary approaches, with TGS annotations identifying more annotated insertions with greater precision for a target TE family, and WGS data allowing for estimates of TE insertion frequencies, and a broader picture of the location of non-target TEs across the genome.

A reference gene set construction using RNA-seq of multiple tissues of Chinese giant salamander, Andrias davidianus

Background

Chinese giant salamander (CGS) is the largest extant amphibian species in the world. Owing to its evolutionary position and four peculiar phenomenon of life (longevity, starvation tolerance, regenerative ability, and hatch without sunshine), it is an invaluable model species for research. However, lack of genomic resources leads to fewer study progresses in these fields, due to its huge genome of ∼50 GB making it extremely difficult to be assembled.

Results

We reported the sequenced transcriptome of more than 20 tissues from adult CGS using Illumina Hiseq 2000 technology, and a total of 93 366 no-redundancy transcripts with a mean length of 1326 bp were obtained. We developed for the first time an efficient pipeline to construct a high-quality reference gene set of CGS and obtained 26 135 coding genes. BUSCO and homologous assessment showed that our assembly captured 70.6% of vertebrate universal single-copy orthologs, and this coding gene set had a higher proportion of completeness CDS with comparable quality of the protein sets of Tibetan frog.

Conclusions

These highest quality data will provide a valuable reference gene set to the subsequent research of CGS. In addition, our strategy of de novo transcriptome assembly and protein identification is applicable to similar studies.

Small RNAs from a Big Genome: The piRNA Pathway and Transposable Elements in the Salamander Species Desmognathus fuscus

Most of the largest vertebrate genomes are found in salamanders, a clade of amphibians that includes 686 species. Salamander genomes range in size from 14 to 120 Gb, reflecting the accumulation of large numbers of transposable element (TE) sequences from all three TE classes. Although DNA loss rates are slow in salamanders relative to other vertebrates, high levels of TE insertion are also likely required to explain such high TE loads. Across the Tree of Life, novel TE insertions are suppressed by several pathways involving small RNA molecules. In most known animals, TE activity in the germline is primarily regulated by the Piwi-interacting RNA (piRNA) pathway. In this study, we test the hypothesis that salamanders' unusually high TE loads reflect the loss of the ancestral piRNA-mediated TE-silencing machinery. We characterized the small RNA pool in the female and male adult gonads, testing for the presence of small RNA molecules that bear the characteristics of TE-targeting piRNAs. We also analyzed the amino acid sequences of piRNA pathway proteins from salamanders and other vertebrates, testing whether the overall patterns of sequence divergence are consistent with conserved pathway function across the vertebrate clade. Our results do not support the hypothesis of piRNA pathway loss; instead, they suggest that the piRNA pathway is expressed in salamanders. Given these results, we propose hypotheses to explain how the extraordinary TE loads in salamander genomes could have accumulated, despite the expression of TE-silencing machinery.

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.

Transposons, Genome Size, and Evolutionary Insights in Animals

The relationship between genome size and the percentage of transposons in 161 animal species evidenced that variations in genome size are linked to the amplification or the contraction of transposable elements. The activity of transposable elements could represent a response to environmental stressors. Indeed, although with different trends in protostomes and deuterostomes, comprehensive changes in genome size were recorded in concomitance with particular periods of evolutionary history or adaptations to specific environments. During evolution, genome size and the presence of transposable elements have influenced structural and functional parameters of genomes and cells. Changes of these parameters have had an impact on morphological and functional characteristics of the organism on which natural selection directly acts. Therefore, the current situation represents a balance between insertion and amplification of transposons and the mechanisms responsible for their deletion or for decreasing their activity. Among the latter, methylation and the silencing action of small RNAs likely represent the most frequent mechanisms.

Microsatellite landscape evolutionary dynamics across 450 million years of vertebrate genome evolution

The evolutionary dynamics of simple sequence repeats (SSRs or microsatellites) across the vertebrate tree of life remain largely undocumented and poorly understood. In this study, we analyzed patterns of genomic microsatellite abundance and evolution across 71 vertebrate genomes. The highest abundances of microsatellites exist in the genomes of ray-finned fishes, squamate reptiles, and mammals, while crocodilian, turtle, and avian genomes exhibit reduced microsatellite landscapes. We used comparative methods to infer evolutionary rates of change in microsatellite abundance across vertebrates and to highlight particular lineages that have experienced unusually high or low rates of change in genomic microsatellite abundance. Overall, most variation in microsatellite content, abundance, and evolutionary rate is observed among major lineages of reptiles, yet we found that several deeply divergent clades (i.e., squamate reptiles and mammals) contained relatively similar genomic microsatellite compositions. Archosauromorph reptiles (turtles, crocodilians, and birds) exhibit reduced genomic microsatellite content and the slowest rates of microsatellite evolution, in contrast to squamate reptile genomes that have among the highest rates of microsatellite evolution. Substantial branch-specific shifts in SSR content in primates, monotremes, rodents, snakes, and fish are also evident. Collectively, our results support multiple major shifts in microsatellite genomic landscapes among vertebrates.

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.

Low levels of LTR retrotransposon deletion by ectopic recombination in the gigantic genomes of salamanders

Across the tree of life, species vary dramatically in nuclear genome size. Mutations that add or remove sequences from genomes-insertions or deletions, or indels-are the ultimate source of this variation. Differences in the tempo and mode of insertion and deletion across taxa have been proposed to contribute to evolutionary diversity in genome size. Among vertebrates, most of the largest genomes are found within the salamanders, an amphibian clade with genome sizes ranging from ~14 to ~120 Gb. Salamander genomes have been shown to experience slower rates of DNA loss through small (i.e., <30 bp) deletions than do other vertebrate genomes. However, no studies have addressed DNA loss from salamander genomes resulting from larger deletions. Here, we focus on one type of large deletion-ectopic-recombination-mediated removal of LTR retrotransposon sequences. In ectopic recombination, double-strand breaks are repaired using a "wrong" (i.e., ectopic, or non-allelic) template sequence-typically another locus of similar sequence. When breaks occur within the LTR portions of LTR retrotransposons, ectopic-recombination-mediated repair can produce deletions that remove the internal transposon sequence and the equivalent of one of the two LTR sequences. These deletions leave a signature in the genome-a solo LTR sequence. We compared levels of solo LTRs in the genomes of four salamander species with levels present in five vertebrates with smaller genomes. Our results demonstrate that salamanders have low levels of solo LTRs, suggesting that ectopic-recombination-mediated deletion of LTR retrotransposons occurs more slowly than in other vertebrates with smaller genomes.

Variation in salamanders: an essay on genomes, development, and evolution

Regeneration is studied in a few model species of salamanders, but the ten families of salamanders show considerable variation, and this has implications for our understanding of salamander biology. The most recent classification of the families identifies the cryptobranchoidea as the basal group which diverged in the early Jurassic. Variation in the sizes of genomes is particularly obvious, and reflects a major contribution from transposable elements which is already present in the basal group.Limb development has been a focus for evodevo studies, in part because of the variable property of pre-axial dominance which distinguishes salamanders from other tetrapods. This is thought to reflect the selective pressures that operate on a free-living aquatic larva, and might also be relevant for the evolution of limb regeneration. Recent fossil evidence suggests that both pre-axial dominance and limb regeneration were present 300 million years ago in larval temnospondyl amphibians that lived in mountain lakes. A satisfying account of regeneration in salamanders may need to address all these different aspects in the future.