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

Comparative analysis of amphibian genomes: An emerging resource for basic and applied research

Amphibians are the most threatened group of vertebrates and are in dire need of conservation intervention to ensure their continued survival. They exhibit unique features including a high diversity of reproductive strategies, permeable and specialized skin capable of producing toxins and antimicrobial compounds, multiple genetic mechanisms of sex determination and in some lineages, the ability to regenerate limbs and organs. Although genomic approaches would shed light on these unique traits and aid conservation, sequencing and assembly of amphibian genomes has lagged behind other taxa due to their comparatively large genome sizes. Fortunately, the development of long-read sequencing technologies and initiatives has led to a recent burst of new amphibian genome assemblies. Although growing, the field of amphibian genomics suffers from the lack of annotation resources, tools for working with challenging genomes and lack of high-quality assemblies in multiple clades of amphibians. Here, we analyse 51 publicly available amphibian genomes to evaluate their usefulness for functional genomics research. We report considerable variation in genome assembly quality and completeness and report some of the highest transposable element and repeat contents of any vertebrate. Additionally, we detected an association between transposable element content and climatic variables. Our analysis provides evidence of conserved genome synteny despite the long divergence times of this group, but we also highlight inconsistencies in chromosome naming and orientation across genome assemblies. We discuss sequencing gaps in the phylogeny and suggest key targets for future sequencing endeavours. Finally, we propose increased investment in amphibian genomics research to promote their conservation.

Key Proteins for Regeneration in A. mexicanum: Transcriptomic Insights From Aged and Juvenile Limbs

The axolotl, known for its remarkable regenerative abilities, is an excellent model for studying regenerative therapies. Nevertheless, the precise molecular mechanisms governing its regenerative potential remain uncertain. In this study, we collected samples from axolotls of different ages, including 8-year-old individuals and 8-month-old juveniles, obtaining their blastemas 10 days after amputation. Subsequently, we conducted a transcriptomic analysis comparing our samples to a set of previously published experiments. Our analysis unveiled a distinctive transcriptional response in the blastema, characterized by differential gene expression associated with processes such as bone and tissue remodeling, transcriptional regulation, angiogenesis, and intercellular communication. To gain deeper insights, we compared these findings with those from aged axolotls that showed no signs of regeneration 10 days after amputation. We identified four genes-FSTL1, ADAMTS17, GPX7, and CTHRC1-that showed higher expression in regenerating tissue compared to aged axolotls. Further scrutiny, including structural and homology analysis, revealed that these genes are conserved across vertebrate species. Our discoveries point to a group of proteins relevant to tissue regeneration, with their conservation in vertebrates suggesting critical roles in development. These findings also propose a novel gene set involved in axolotl regeneration, laying a promising foundation for future investigations across vertebrates.

The Molecular Basis of Differentiation Wave Activity in Embryogenesis

As development varies greatly across the tree of life, it may seem difficult to suggest a model that proposes a single mechanism for understanding collective cell behaviors and the coordination of tissue formation. Here we propose a mechanism called differentiation waves, which unify many disparate results involving developmental systems from across the tree of life. We demonstrate how a relatively simple model of differentiation proceeds not from function-related molecular mechanisms, but from so-called differentiation waves. A phenotypic model of differentiation waves is introduced, and its relation to molecular mechanisms is proposed. These waves contribute to a differentiation tree, which is an alternate way of viewing cell lineage and local action of the molecular factors. We construct a model of differentiation wave-related molecular mechanisms (genome, epigenome, and proteome) based on bioinformatic data from the nematode Caenorhabditis elegans. To validate this approach across different modes of development, we evaluate protein expression across different types of development by comparing Caenorhabditis elegans with several model organisms: fruit flies (Drosophila melanogaster), yeast (Saccharomyces cerevisiae), and mouse (Mus musculus). Inspired by gene regulatory networks, two Models of Interactive Contributions (fully-connected MICs and ordered MICs) are used to suggest potential genomic contributions to differentiation wave-related proteins. This, in turn, provides a framework for understanding differentiation and development.

The salamander blastema within the broader context of metazoan regeneration

Throughout the animal kingdom regenerative ability varies greatly from species to species, and even tissue to tissue within the same organism. The sheer diversity of structures and mechanisms renders a thorough comparison of molecular processes truly daunting. Are "blastemas" found in organisms as distantly related as planarians and axolotls derived from the same ancestral process, or did they arise convergently and independently? Is a mouse digit tip blastema orthologous to a salamander limb blastema? In other fields, the thorough characterization of a reference model has greatly facilitated these comparisons. For example, the amphibian Spemann-Mangold organizer has served as an amazingly useful comparative template within the field of developmental biology, allowing researchers to draw analogies between distantly related species, and developmental processes which are superficially quite different. The salamander limb blastema may serve as the best starting point for a comparative analysis of regeneration, as it has been characterized by over 200 years of research and is supported by a growing arsenal of molecular tools. The anatomical and evolutionary closeness of the salamander and human limb also add value from a translational and therapeutic standpoint. Tracing the evolutionary origins of the salamander blastema, and its relatedness to other regenerative processes throughout the animal kingdom, will both enhance our basic biological understanding of regeneration and inform our selection of regenerative model systems.

Limb blastema formation: How much do we know at a genetic and epigenetic level?

Regeneration of missing body parts is an incredible ability which is present in a wide number of species. However, this regenerative capability varies among different organisms. Urodeles (salamanders) are able to completely regenerate limbs after amputation through the essential process of blastema formation. The blastema is a collection of relatively undifferentiated progenitor cells that proliferate and repattern to form the internal tissues of a regenerated limb. Understanding blastema formation in salamanders may enable comparative studies with other animals, including mammals, with more limited regenerative abilities and may inspire future therapeutic approaches in humans. This review focuses on the current state of knowledge about how limb blastemas form in salamanders, highlighting both the possible roles of epigenetic controls in this process as well as limitations to scientific understanding that present opportunities for research.

Somatic CNV Detection by Single-Cell Whole-Genome Sequencing in Postmortem Human Brain

The evidence for a role of somatic mutations, including copy-number variants (CNVs), in neurodegeneration has increased in the last decade. However, the understanding of the types and origins of these mutations, and their exact contributions to disease onset and progression, is still in its infancy. The use of single-cell (or nuclear) whole-genome sequencing (scWGS) has emerged as a powerful tool to answer these questions. In the present chapter, we provide laboratory and bioinformatic protocols used successfully in our lab to detect megabase-scale CNVs in single cells from multiple system atrophy (MSA) human postmortem brains, using immunolabeling prior to selection of nuclei for whole-genome amplification (WGA). We also present an unpublished comparison of scWGS generated from the same control substantia nigra (SN) sample, using the latest versions of popular WGA chemistries, MDA and PicoPLEX. We have used this protocol to focus on brain cell types most relevant to synucleinopathies (dopaminergic [DA] neurons in Parkinson's disease [PD] and oligodendrocytes in MSA), but it can be applied to any tissue and/or cell type with appropriate markers.

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.

Fluorescence In Situ Hybridization of DNA Probes on Mitotic Chromosomes of the Mexican Axolotl

Fluorescence in situ hybridization (FISH) is used extensively for visual localization of specific DNA fragments (and RNA fragments) in broad applications on chromosomes or nuclei at any stage of the cell cycle: metaphase, anaphase, or interphase. The cytogenetic slides that serve as a target for the labeled DNA probe might be prepared using any approach suitable for obtaining cells with appropriate morphology for imaging and analysis. In this chapter, we focus on the application of molecular cytogenetic methods such as DNA labeling, slide preparation, and in situ hybridization related to cells from Mexican axolotl.

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.

Evolutionary conservation of maternal RNA localization in fishes and amphibians revealed by TOMO-Seq

Asymmetrical localization of biomolecules inside the egg, results in uneven cell division and establishment of many biological processes, cell types and the body plan. However, our knowledge about evolutionary conservation of localized transcripts is still limited to a few models. Our goal was to compare localization profiles along the animal-vegetal axis of mature eggs from four vertebrate models, two amphibians (Xenopus laevis, Ambystoma mexicanum) and two fishes (Acipenser ruthenus, Danio rerio) using the spatial expression method called TOMO-Seq. We revealed that RNAs of many known important transcripts such as germ layer determinants, germ plasm factors and members of key signalling pathways, are localized in completely different profiles among the models. It was also observed that there was a poor correlation between the vegetally localized transcripts but a relatively good correlation between the animally localized transcripts. These findings indicate that the regulation of embryonic development within the animal kingdom is highly diverse and cannot be deduced based on a single model.

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.

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.

A genomic survey of LINE elements in Pipidae aquatic frogs shed light on Rex-elements evolution in these genomes

The transposable elements (TE) represent a large portion of anuran genomes that act as components of genetic diversification. The LINE order of retrotransposons is among the most representative and diverse TEs and is poorly investigated in anurans. Here we explored the LINE diversity with an emphasis on the elements generically called Rex in Pipidae species, more specifically, in the genomes ofXenopus tropicalis, used as a model genome in the study of anurans,the allotetraploid sister species Xenopus laevis and theAmerican species Pipa carvalhoi. We were able to identify a great diversity of LINEs from five clades, Rex1, L2, CR1, L1 and Tx1, in these three species, and the RTE clade was lost in X. tropicalis. It is clear that elements classified as Rex are distributed in distinct clades. The evolutionary pattern of Rex1 elements denote a complex evolution with independent losses of families and some horizontal transfer events between fishes and amphibians which were supported not only by the phylogenetic inconsistencies but also by the very low Ks values found for the TE sequences. The data obtained here update the knowledge of the LINEs diversity in X. laevis and represent the first study of TEs in P. carvalhoi.

The Axolotl's journey to the modern molecular era

The salamander Ambystoma mexicanum, commonly called "the axolotl" has a long, illustrious history as a model organism, perhaps with one of the longest track records as a laboratory-bred vertebrate, yet it also holds a prominent place among the emerging model organisms. Or rather it is by now an "emerged" model organism, boasting a full cohort molecular genetic tools that allows an expanding community of researchers in the field to explore the remarkable traits of this animal including regeneration, at cellular and molecular precision-which had been a dream for researchers over the years. This chapter describes the journey to this status, that could be helpful for those developing their respective animal or plant models.

Newt regeneration genes regulate Wingless signaling to restore patterning in Drosophila eye

Newts utilize their unique genes to restore missing parts by strategic regulation of conserved signaling pathways. Lack of genetic tools poses challenges to determine the function of such genes. Therefore, we used the Drosophila eye model to demonstrate the potential of 5 unique newt (Notophthalmus viridescens) gene(s), viropana1-viropana5 (vna1-vna5), which were ectopically expressed in L 2 mutant and GMR-hid, GMR-GAL4 eye. L 2 exhibits the loss of ventral half of early eye and head involution defective (hid) triggers cell-death during later eye development. Surprisingly, newt genes significantly restore missing photoreceptor cells both in L 2 and GMR>hid background by upregulating cell-proliferation and blocking cell-death, regulating evolutionarily conserved Wingless (Wg)/Wnt signaling pathway and exhibit non-cell-autonomous rescues. Further, Wg/Wnt signaling acts downstream of newt genes. Our data highlights that unique newt proteins can regulate conserved pathways to trigger a robust restoration of missing photoreceptor cells in Drosophila eye model with weak restoration capability.

Ambystoma mexicanum, a model organism in developmental biology and regeneration: a colombian experience

Ambystoma mexicanum is a urodele amphibian endemic to Xochimilco Lake in Mexico, it belongs to the salamander family Ambystomatidae. This species has frequently been used as model organism in developmental biology and regeneration laboratories around the world due to its broad regenerative capacities and adaptability to laboratory conditions. In this review we describe the establishment of the first colony of axolotls in Colombia to study tissue regeneration and our perspectives on the use A. mexicanum as a model organism in Colombia are discussed emphasizing its possible uses in regeneration and developmental biology

Salamanders: The molecular basis of tissue regeneration and its relevance to human disease

Salamanders are recognized for their ability to regenerate a broad range of tissues. They have also have been used for hundreds of years for classical developmental biology studies because of their large accessible embryos. The range of tissues these animals can regenerate is fascinating, from full limbs to parts of the brain or heart, a potential that is missing in humans. Many promising research efforts are working to decipher the molecular blueprints shared across the organisms that naturally have the capacity to regenerate different tissues and organs. Salamanders are an excellent example of a vertebrate that can functionally regenerate a wide range of tissue types. In this review, we outline some of the significant insights that have been made that are aiding in understanding the cellular and molecular mechanisms of tissue regeneration in salamanders and discuss why salamanders are a worthy model in which to study regenerative biology and how this may benefit research fields like regenerative medicine to develop therapies for humans in the future.

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.

Bridging the Gap between Vertebrate Cytogenetics and Genomics with Single-Chromosome Sequencing (ChromSeq)

The study of vertebrate genome evolution is currently facing a revolution, brought about by next generation sequencing technologies that allow researchers to produce nearly complete and error-free genome assemblies. Novel approaches however do not always provide a direct link with information on vertebrate genome evolution gained from cytogenetic approaches. It is useful to preserve and link cytogenetic data with novel genomic discoveries. Sequencing of DNA from single isolated chromosomes (ChromSeq) is an elegant approach to determine the chromosome content and assign genome assemblies to chromosomes, thus bridging the gap between cytogenetics and genomics. The aim of this paper is to describe how ChromSeq can support the study of vertebrate genome evolution and how it can help link cytogenetic and genomic data. We show key examples of ChromSeq application in the refinement of vertebrate genome assemblies and in the study of vertebrate chromosome and karyotype evolution. We also provide a general overview of the approach and a concrete example of genome refinement using this method in the species Anolis carolinensis.

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.

EXFI: Exon and splice graph prediction without a reference genome

For population genetic studies in nonmodel organisms, it is important to use every single source of genomic information. This paper presents EXFI, a Python pipeline that predicts the splice graph and exon sequences using an assembled transcriptome and raw whole-genome sequencing reads. The main algorithm uses Bloom filters to remove reads that are not part of the transcriptome, to predict the intron-exon boundaries, to then proceed to call exons from the assembly, and to generate the underlying splice graph. The results are returned in GFA1 format, which encodes both the predicted exon sequences and how they are connected to form transcripts. EXFI is written in Python, tested on Linux platforms, and the source code is available under the MIT License at https://github.com/jlanga/exfi.

Functional colour genes and signals of selection in colour-polymorphic salamanders

Coloration has been associated with multiple biologically relevant traits that drive adaptation and diversification in many taxa. However, despite the great diversity of colour patterns present in amphibians the underlying molecular basis is largely unknown. Here, we use insight from a highly colour-variable lineage of the European fire salamander (Salamandra salamandra bernardezi) to identify functional associations with striking variation in colour morph and pattern. The three focal colour morphs-ancestral black-yellow striped, fully yellow and fully brown-differed in pattern, visible coloration and cellular composition. From population genomic analyses of up to 4,702 loci, we found no correlations of neutral population genetic structure with colour morph. However, we identified 21 loci with genotype-phenotype associations, several of which relate to known colour genes. Furthermore, we inferred response to selection at up to 142 loci between the colour morphs, again including several that relate to coloration genes. By transcriptomic analysis across all different combinations, we found 196 differentially expressed genes between yellow, brown and black skin, 63 of which are candidate genes involved in animal coloration. The concordance across different statistical approaches and 'omic data sets provide several lines of evidence for loci linked to functional differences between colour morphs, including TYR, CAMK1 and PMEL. We found little association between colour morph and the metabolomic profile of its toxic compounds from the skin secretions. Our research suggests that current ecological and evolutionary hypotheses for the origins and maintenance of these striking colour morphs may need to be revisited.

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.

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.

HDAC Regulates Transcription at the Outset of Axolotl Tail Regeneration

Tissue regeneration is associated with complex changes in gene expression and post-translational modifications of proteins, including transcription factors and histones that comprise chromatin. We tested 172 compounds designed to target epigenetic mechanisms in an axolotl (Ambystoma mexicanum) embryo tail regeneration assay. A relatively large number of compounds (N = 55) inhibited tail regeneration, including 18 histone deacetylase inhibitors (HDACi). In particular, romidepsin, an FDA-approved anticancer drug, potently inhibited tail regeneration when embryos were treated continuously for 7 days. Additional experiments revealed that romidepsin acted within a very narrow, post-injury window. Romidepsin treatment for only 1-minute post amputation inhibited regeneration through the first 7 days, however after this time, regeneration commenced with variable outgrowth of tailfin tissue and abnormal patterning. Microarray analysis showed that romidepsin altered early, transcriptional responses at 3 and 6-hour post-amputation, especially targeting genes that are implicated in tumor cell death, as well as genes that function in the regulation of transcription, cell differentiation, cell proliferation, pattern specification, and tissue morphogenesis. Our results show that HDAC activity is required at the time of tail amputation to regulate the initial transcriptional response to injury and regeneration.

FGF, BMP, and RA signaling are sufficient for the induction of complete limb regeneration from non-regenerating wounds on Ambystoma mexicanum limbs

Some organisms, such as the Mexican axolotl, have the capacity to regenerate complicated biological structures throughout their lives. Which molecular pathways are sufficient to induce a complete endogenous regenerative response in injured tissue is an important question that remains unanswered. Using a gain-of-function regeneration assay, known as the Accessory Limb Model (ALM), we and others have begun to identify the molecular underpinnings of the three essential requirements for limb regeneration; wounding, neurotrophic signaling, and the induction of pattern from cells that retain positional memory. We have previously shown that treatment of Mexican axolotls with exogenous retinoic acid (RA) is sufficient to induce the formation of complete limb structures from blastemas that were generated by deviating a nerve bundle into an anterior-located wound site on the limb. Here we show that these ectopic structures are capable of regenerating and inducing new pattern to form when grafted into new anterior-located wounds. We additionally found that the expression of Alx4 decreases, and Shh expression increases in these anterior located blastemas, but not in the mature anterior tissues, supporting the hypothesis that RA treatment posteriorizes blastema tissue. Based on these and previous observations, we used the ALM assay to test the hypothesis that a complete regenerative response can be generated by treating anterior-located superficial limb wounds with a specific combination of growth factors at defined developmental stages. Our data shows that limb wounds that are first treated with a combination of FGF-2, FGF-8, and BMP-2, followed by RA treatment of the resultant mid-bud stage blastema, will result in the generation of limbs with complete proximal/distal and anterior/posterior limb axes. Thus, the minimal signaling requirements from the nerve and a positional disparity are achieved with the application of this specific combination of signaling molecules.

Taxon-specific expansion and loss of tektins inform metazoan ciliary diversity

Background

Cilia and flagella are complex cellular structures thought to have first evolved in a last ciliated eukaryotic ancestor due to the conserved 9 + 2 microtubule doublet structure of the axoneme and associated proteins. The Tektin family of coiled-coil domain containing proteins was previously identified in cilia of organisms as diverse as green algae and sea urchin. While studies have shown that some Tektins are necessary for ciliary function, there has been no comprehensive phylogenetic survey of tektin genes. To fill this gap, we sampled tektin sequences broadly among metazoan and unicellular lineages in order to determine how the tektin gene complements evolved in over 100 different extant species.

Results

Using Bayesian and Maximum Likelihood analyses, we have ascertained with high confidence that all metazoan tektins arose from a single ancestral tektin gene in the last common ancestor of metazoans and choanoflagellates. Gene duplications gave rise to two tektin genes in the metazoan ancestor, and a subsequent expansion to three and four tektin genes in early bilaterian ancestors. While all four tektin genes remained highly conserved in most deuterostome and spiralian species surveyed, most tektin genes in ecdysozoans are highly derived with extensive gene loss in several lineages including nematodes and some crustaceans. In addition, while tektin-1, − 2, and − 4 have remained as single copy genes in most lineages, tektin-3/5 has been duplicated independently several times, notably at the base of the spiralian, vertebrate and hymenopteran (Ecdysozoa) clades.

Conclusions

We provide a solid description of tektin evolution supporting one, two, three, and four ancestral tektin genes in a holozoan, metazoan, bilaterian, and nephrozoan ancestor, respectively. The isolated presence of tektin in a cryptophyte and a chlorophyte branch invokes events of horizontal gene transfer, and that the last common ciliated eukaryotic ancestor lacked a tektin gene. Reconstructing the evolutionary history of the tektin complement in each extant metazoan species enabled us to pinpoint lineage specific expansions and losses. Our analysis will help to direct future studies on Tektin function, and how gain and loss of tektin genes might have contributed to the evolution of various types of cilia and flagella.

Electronic supplementary material

The online version of this article (10.1186/s12862-019-1360-0) contains supplementary material, which is available to authorized users.

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,000^th 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.

Transcriptomic landscape of the blastema niche in regenerating adult axolotl limbs at single-cell resolution

Regeneration of complex multi-tissue structures, such as limbs, requires the coordinated effort of multiple cell types. In axolotl limb regeneration, the wound epidermis and blastema have been extensively studied via histology, grafting, and bulk-tissue RNA-sequencing. However, defining the contributions of these tissues is hindered due to limited information regarding the molecular identity of the cell types in regenerating limbs. Here we report unbiased single-cell RNA-sequencing on over 25,000 cells from axolotl limbs and identify a plethora of cellular diversity within epidermal, mesenchymal, and hematopoietic lineages in homeostatic and regenerating limbs. We identify regeneration-induced genes, develop putative trajectories for blastema cell differentiation, and propose the molecular identity of fibroblast-like blastema progenitor cells. This work will enable application of molecular techniques to assess the contribution of these populations to limb regeneration. Overall, these data allow for establishment of a putative framework for adult axolotl limb regeneration.

Limb regeneration requires a blastema with progenitor cells, immune cells, and an overlying wound epidermis, but molecular identities of these populations are unclear. Here, the authors use single-cell RNA-sequencing to identify transcriptionally distinct cell populations in adult axolotl limb blastemas.

Identification of critical sex-biased genes in Andrias davidianus by de novo transcriptome

The Chinese giant salamander Andrias davidianus is a protected amphibian with high nutritional and economic value. Understanding its sex determination mechanism is important for improving culture techniques and sex control in breeding. However, little information on the characterization of critical genes involved in sex is available. Herein, sequencing of ovary and test produced 40,783,222 and 46,128,902 raw reads, respectively, which were jointly assembled into 80,497 unigenes. Of these, 36,609 unigenes were annotated, of which 8907 were female-biased and 10,385 were male-biased. Several sex-related pathways were observed, including the Wnt signaling pathway. After elevated temperature and estrogen exposure, neomale and neofemale specimens were identified by a female-specific marker for the first time. RT-qPCR analysis showed the expression profile of ten selected sex-biased genes to be exhibited consistently in male and neomale and in female and neofemale, with the exception of the Amh and TfIIIa genes. Results suggested that these genes may play important roles in A. davidianus sex determination and gonad development. This provides a basis for further investigation of the molecular mechanisms of sex determination in amphibians.

Genomic data reject the hypothesis of sympatric ecological speciation in a clade of Desmognathus salamanders

Closely related taxa with dissimilar morphologies are often considered to have diverged via natural selection favoring different phenotypes. However, some studies have found these scenarios to be paired with limited or no genetic differentiation. Desmognathus quadramaculatus and D. marmoratus are sympatric salamander species thought to represent a case of ecological speciation based on distinct morphologies, but the results of previous studies have not resolved corresponding patterns of lineage divergence. Here, we use genome-wide data to test this hypothesis of ecological speciation. Population structure analyses partitioned individuals geographically, but not morphologically, into two adjacent regions of western North Carolina: Pisgah and Nantahala. Phylogenetic analyses confirmed the nominal species are nonmonophyletic and resolved deep divergence between the two geographic clusters. Model-testing overwhelmingly supported the hypothesis that lineage divergence followed geography. Finally, ecological niche modeling showed that Pisgah and Nantahala individuals occupy different climatic niches, and geographic boundaries for the two lineages correspond to differences in precipitation regimes across southern Appalachia. Overall, we reject the previous hypothesis of ecological speciation based on microhabitat partitioning. Instead, our results suggest that there are two cryptic lineages, each containing the same pair of morphotypes.

The Most Developmentally Truncated Fishes Show Extensive Hox Gene Loss and Miniaturized Genomes

The world’s smallest fishes belong to the genus Paedocypris. These miniature fishes are endemic to an extreme habitat: the peat swamp forests in Southeast Asia, characterized by highly acidic blackwater. This threatened habitat is home to a large array of fishes, including a number of miniaturized but also developmentally truncated species. Especially the genus Paedocypris is characterized by profound, organism-wide developmental truncation, resulting in sexually mature individuals of <8 mm in length with a larval phenotype. Here, we report on evolutionary simplification in the genomes of two species of the dwarf minnow genus Paedocypris using whole-genome sequencing. The two species feature unprecedented Hox gene loss and genome reduction in association with their massive developmental truncation. We also show how other genes involved in the development of musculature, nervous system, and skeleton have been lost in Paedocypris, mirroring its highly progenetic phenotype. Further, our analyses suggest two mechanisms responsible for the genome streamlining in Paedocypris in relation to other Cypriniformes: severe intron shortening and reduced repeat content. As the first report on the genomic sequence of a vertebrate species with organism-wide developmental truncation, the results of our work enhance our understanding of genome evolution and how genotypes are translated to phenotypes. In addition, as a naturally simplified system closely related to zebrafish, Paedocypris provides novel insights into vertebrate development.

Experimentally induced metamorphosis in highly regenerative axolotl (ambystoma mexicanum) under constant diet restructures microbiota

Axolotl (Ambystoma mexicanum) is a critically endangered salamander species and a model organism for regenerative and developmental biology. Despite life-long neoteny in nature and in captive-bred colonies, metamorphosis of these animals can be experimentally induced by administering Thyroid hormones (THs). However, microbiological consequences of this experimental procedure, such as host microbiota response, remain largely unknown. Here, we systematically compared host bacterial microbiota associated with skin, stomach, gut tissues and fecal samples, between neotenic and metamorphic axolotls based on 16S rRNA gene sequences. Our results show that distinct bacterial communities inhabit individual organs of axolotl and undergo substantial restructuring through metamorphosis. Skin microbiota among others, shifted sharply, as highlighted by a major transition from Firmicutes-enriched to Proteobacteria-enriched relative abundance and precipitously decreased diversity. Fecal microbiota of neotenic and metamorphic axolotl shared relatively higher similarity, suggesting that diet continues to shape microbiota despite fundamental transformations in the host digestive organs. We also reproduced the previous finding on reduction in regenerative capacity in limbs of axolotl following metamorphosis, highlighting the need to investigate whether shifts in microbiota is causally linked to regenerative capacity of axolotl. The initial results on axolotl microbiota provide novel insights into microbiological aspects of axolotl metamorphosis and will establish a baseline for future in-depth studies.

Genome-specific histories of divergence and introgression between an allopolyploid unisexual salamander lineage and two ancestral sexual species

Quantifying introgression between sexual species and polyploid lineages traditionally thought to be asexual is an important step in understanding what drives the longevity of putatively asexual groups. Here, we capitalize on three recent innovations-ultraconserved element (UCE) sequencing, bioinformatic techniques for identifying genome-specific variation in polyploids, and model-based methods for evaluating historical gene flow-to measure the extent and tempo of introgression over the evolutionary history of an allopolyploid lineage of all-female salamanders and two ancestral sexual species. Our analyses support a scenario in which the genomes sampled in unisexual salamanders last shared a common ancestor with genomes in their parental species ∼3.4 million years ago, followed by a period of divergence between homologous genomes. Recently, secondary introgression has occurred at different times with each sexual species during the last 500,000 years. Sustained introgression of sexual genomes into the unisexual lineage is the defining characteristic of their reproductive mode, but this study provides the first evidence that unisexual genomes have undergone long periods of divergence without introgression. Unlike other sperm-dependent taxa in which introgression is rare, the alternating periods of divergence and introgression between unisexual salamanders and their sexual relatives could explain why these salamanders are among the oldest described unisexual animals.

Polymorphic adaptations in metazoans to establish and maintain photosymbioses

Mutualistic symbioses are common throughout the animal kingdom. Rather unusual is a form of symbiosis, photosymbiosis, where animals are symbiotic with photoautotrophic organisms. Photosymbiosis is found among sponges, cnidarians, flatworms, molluscs, ascidians and even some amphibians. Generally the animal host harbours a phototrophic partner, usually a cyanobacteria or a unicellular alga. An exception to this rule is found in some sea slugs, which only retain the chloroplasts of the algal food source and maintain them photosynthetically active in their own cytosol - a phenomenon called 'functional kleptoplasty'. Research has focused largely on the biodiversity of photosymbiotic species across a range of taxa. However, many questions with regard to the evolution of the ability to establish and maintain a photosymbiosis are still unanswered. To date, attempts to understand genome adaptations which could potentially lead to the evolution of photosymbioses have only been performed in cnidarians. This knowledge gap for other systems is mainly due to a lack of genetic information, both for non-symbiotic and symbiotic species. Considering non-photosymbiotic species is, however, important to understand the factors that make symbiotic species so unique. Herein we provide an overview of the diversity of photosymbioses across the animal kingdom and discuss potential scenarios for the evolution of this association in different lineages. We stress that the evolution of photosymbiosis is probably based on genome adaptations, which (i) lead to recognition of the symbiont to establish the symbiosis, and (ii) are needed to maintain the symbiosis. We hope to stimulate research involving sequencing the genomes of various key taxa to increase the genomic resources needed to understand the most fundamental question: how have animals evolved the ability to establish and maintain a photosymbiosis?

Musashi and Plasticity of Xenopus and Axolotl Spinal Cord Ependymal Cells

The differentiated state of spinal cord ependymal cells in regeneration-competent amphibians varies between a constitutively active state in what is essentially a developing organism, the tadpole of the frog Xenopus laevis, and a quiescent, activatable state in a slowly growing adult salamander Ambystoma mexicanum, the Axolotl. Ependymal cells are epithelial in intact spinal cord of all vertebrates. After transection, body region ependymal epithelium in both Xenopus and the Axolotl disorganizes for regenerative outgrowth (gap replacement). Injury-reactive ependymal cells serve as a stem/progenitor cell population in regeneration and reconstruct the central canal. Expression patterns of mRNA and protein for the stem/progenitor cell-maintenance Notch signaling pathway mRNA-binding protein Musashi (msi) change with life stage and regeneration competence. Msi-1 is missing (immunohistochemistry), or at very low levels (polymerase chain reaction, PCR), in both intact regeneration-competent adult Axolotl cord and intact non-regeneration-competent Xenopus tadpole (Nieuwkoop and Faber stage 62+, NF 62+). The critical correlation for successful regeneration is msi-1 expression/upregulation after injury in the ependymal outgrowth and stump-region ependymal cells. msi-1 and msi-2 isoforms were cloned for the Axolotl as well as previously unknown isoforms of Xenopus msi-2. Intact Xenopus spinal cord ependymal cells show a loss of msi-1 expression between regeneration-competent (NF 50-53) and non-regenerating stages (NF 62+) and in post-metamorphosis froglets, while msi-2 displays a lower molecular weight isoform in non-regenerating cord. In the Axolotl, embryos and juveniles maintain Msi-1 expression in the intact cord. In the adult Axolotl, Msi-1 is absent, but upregulates after injury. Msi-2 levels are more variable among Axolotl life stages: rising between late tailbud embryos and juveniles and decreasing in adult cord. Cultures of regeneration-competent Xenopus tadpole cord and injury-responsive adult Axolotl cord ependymal cells showed an identical growth factor response. Epidermal growth factor (EGF) maintains mesenchymal outgrowth in vitro, the cells are proliferative and maintain msi-1 expression. Non-regeneration competent Xenopus ependymal cells, NF 62+, failed to attach or grow well in EGF+ medium. Ependymal Msi-1 expression in vivo and in vitro is a strong indicator of regeneration competence in the amphibian spinal cord.

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 .

Reading and editing the Pleurodeles waltl genome reveals novel features of tetrapod regeneration

Salamanders exhibit an extraordinary ability among vertebrates to regenerate complex body parts. However, scarce genomic resources have limited our understanding of regeneration in adult salamanders. Here, we present the ~20 Gb genome and transcriptome of the Iberian ribbed newt Pleurodeles waltl, a tractable species suitable for laboratory research. We find that embryonic stem cell-specific miRNAs mir-93b and mir-427/430/302, as well as Harbinger DNA transposons carrying the Myb-like proto-oncogene have expanded dramatically in the Pleurodeleswaltl genome and are co-expressed during limb regeneration. Moreover, we find that a family of salamander methyltransferases is expressed specifically in adult appendages. Using CRISPR/Cas9 technology to perturb transcription factors, we demonstrate that, unlike the axolotl, Pax3 is present and necessary for development and that contrary to mammals, muscle regeneration is normal without functional Pax7 gene. Our data provide a foundation for comparative genomic studies that generate models for the uneven distribution of regenerative capacities among vertebrates.

The Iberian ribbed newt Pleurodeles waltl has a wide spectrum of regeneration abilities. Here, Elewa et al. sequence its ~20 Gb genome and transcriptome to investigate the molecular features underlying its regenerative capacities.

A Tissue-Mapped Axolotl De Novo Transcriptome Enables Identification of Limb Regeneration Factors

Mammals have extremely limited regenerative capabilities; however, axolotls are profoundly regenerative and can replace entire limbs. The mechanisms underlying limb regeneration remain poorly understood, partly because the enormous and incompletely sequenced genomes of axolotls have hindered the study of genes facilitating regeneration. We assembled and annotated a de novo transcriptome using RNA-sequencing profiles for a broad spectrum of tissues that is estimated to have near-complete sequence information for 88% of axolotl genes. We devised expression analyses that identified the axolotl orthologs of cirbp and kazald1 as highly expressed and enriched in blastemas. Using morpholino anti-sense oligonucleotides, we find evidence that cirbp plays a cytoprotective role during limb regeneration whereas manipulation of kazald1 expression disrupts regeneration. Our transcriptome and annotation resources greatly complement previous transcriptomic studies and will be a valuable resource for future research in regenerative biology.

On the independent gene trees assumption in phylogenomic studies

Multilocus coalescent methods for inferring species trees or historical demographic parameters typically require the assumption that gene trees for sampled SNPs or DNA sequence loci are conditionally independent given their species tree. In practice, researchers have used different criteria to delimit "independent loci." One criterion identifies sampled loci as being independent of each other if they undergo Mendelian independent assortment (IA criterion). O'Neill et al. (2013, Molecular Ecology, 22, 111-129) used this approach in their phylogeographic study of North American tiger salamander species complex. In two other studies, researchers developed a pair of related methods that employ an independent genealogies criterion (IG criterion), which considers the effects of population-level recombination on correlations between the gene trees of intrachromosomal loci. Here, I explain these three methods, illustrate their use with example data, and evaluate their efficacies. I show that the IA approach is more conservative, is simpler to use and requires fewer assumptions than the IG approaches. However, IG approaches can identify much larger numbers of independent loci than the IA method, which, in turn, allows researchers to obtain more precise and accurate estimates of species trees and historical demographic parameters. A disadvantage of the IG methods is that they require an estimate of the population recombination rate. Despite their drawbacks, IA and IG approaches provide molecular ecologists with promising a priori methods for selecting SNPs or DNA sequence loci that likely meet the independence assumption in coalescent-based phylogenomic studies.

Rapid de novo assembly of the European eel genome from nanopore sequencing reads

We have sequenced the genome of the endangered European eel using the MinION by Oxford Nanopore, and assembled these data using a novel algorithm specifically designed for large eukaryotic genomes. For this 860 Mbp genome, the entire computational process takes two days on a single CPU. The resulting genome assembly significantly improves on a previous draft based on short reads only, both in terms of contiguity (N50 1.2 Mbp) and structural quality. This combination of affordable nanopore sequencing and light weight assembly promises to make high-quality genomic resources accessible for many non-model plants and animals.